Search

International audience; Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation(1,2). Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox(3,4), it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

Gene duplication and loss are major sources of genetic polymorphism in populations, and are important forces shaping the evolution of genome content and organization. We have reconstructed the origin and history of a 127 kbp segmental duplication,R2d, in the house mouse (Mus musculus). R2dcontains a single protein-coding gene,Cwc22. De novoassembly of both the ancestral (R2d1) and the derived (R2d2) copies reveals that they have been subject to non-allelic gene conversion events spanning tens of kilobases.R2d2is also a hotspot for structural variation: its diploid copy number ranges from zero in the mouse reference genome to more than 80 in wild mice sampled from around the globe. Hemizgyosity for high-copy-number alleles ofR2d2is associated inciswith meiotic drive, suppression of meiotic crossovers, and copy-number instability, with a mutation rate in excess of 1 per 100 transmissions in laboratory populations. We identify an additional 57 loci covering 0.8% of the mouse genome with patterns of sequence variation similar to those atR2d1andR2d2. Our results provide a striking example of allelic diversity generated by duplication and demonstrate the value ofde novoassembly in a phylogenetic context for understanding the mutational processes affecting duplicate genes.

Genetic elements that are inherited at super-Mendelian frequencies could be used in a 'gene drive' to spread an allele to high prevalence in a population with the goal of eliminating invasive species or disease vectors. We recently demonstrated that the gene conversion mechanism underlying a CRISPR-Cas9-mediated gene drive is feasible in mice. Although substantial technical hurdles remain, overcoming these could lead to strategies that might decrease the spread of rodent-borne Lyme disease or eliminate invasive populations of mice and rats that devastate island ecology. Perhaps more immediately achievable at moderate gene conversion efficiency, applications in a laboratory setting could produce complex genotypes that reduce the time and cost in both dollars and animal lives compared with Mendelian inheritance strategies. Here, we discuss what we have learned from early efforts to achieve CRISPR-Cas9-mediated gene conversion, potential for broader applications in the laboratory, current limitations, and plans for optimizing this potentially powerful technology.

Cytonuclear coordination between biparental-nuclear genomes and uniparental-cytoplasmic organellar genomes in plants is often resolved by genetic and transcriptional cytonuclear responses. Whether this mechanism also acts in allopolyploid members of other kingdoms is not clear. Additionally, cytonuclear coordination of interleaved allopolyploid cells/individuals within the same population is underexplored. The yeast Saccharomyces pastorianus provides the opportunity to explore cytonuclear coevolution during different growth stages and from novel dimensions. Using S. pastorianus cells from multiple growth stages in the same environment, we show that nuclear mitochondria-targeted genes have undergone both asymmetric gene conversion and growth stage-specific biased expression favoring genes from the mitochondrial genome donor (Saccharomyces eubayanus). Our results suggest that cytonuclear coordination in allopolyploid lager yeast species entails an orchestrated and compensatory genetic and transcriptional evolutionary regulatory shift. The common as well as unique properties of cytonuclear coordination underlying allopolyploidy between unicellular yeasts and higher plants offers novel insights into mechanisms of cytonuclear evolution associated with allopolyploid speciation.

Wrasses are dominant components of major coral reef systems. Among wrasses, Cheilinus undulatus is an endangered species with high economic and ecological value that exhibits sex reversal of females to males, while sexual selection occurs in breeding aggregations. However, the molecular-associated mechanism underlying this remains unclear. Opsin gene diversification is regarded as a potent force in sexual selection. Here we present a genome assembly of C. undulatus, using Illumina, Nanopore and Hi-C sequencing. The 1.17 Gb genome was generated from 328 contigs with an N50 length of 16.5 Mb and anchored to 24 chromosomes. In total, 22,218 genes were functionally annotated, and 96.36% of BUSCO genes were fully represented. Transcriptomic analyses showed that 96.79% of the predicted genes were expressed. Transposons were most abundant, accounting for 39.88% of the genome, with low divergence, owing to their evolution with close species ~60.53 million years ago. In total, 567/1,826 gene families were expanded and contracted in the reconstructed phylogeny, respectively. Forty-six genes were under positive selection. Comparative genomic analyses with other fish revealed expansion of opsin SWS2B, LWS1 and Rh2. The elevated duplicates of SWS2B were generated by gene conversions via transposition of transposons followed by nonallelic homologous recombination. Amino acid substitutions of opsin paralogues occurred at key tuning sites, causing a spectral shift in maximal absorbance of visual pigment to capture functional changes. Among these opsin genes, SWS2B-3 and 4 and Rh1 are expressed in the retina. The genome sequence of C. undulatus provides valuable resources for future investigation of the conservation, evolution and behaviour of fishes.

Meiotic recombination is essential for chromosome segregation at meiosis and fertility. It is initiated by programmed DNA double-strand breaks (DSBs) introduced by Spo11, a eukaryotic homologue of an archaeal topoisomerase (Topo VIA)1. Here we describe previously uncharacterized Spo11-induced lesions, 34 to several hundred base pair-long gaps, which are generated by coordinated pairs of DSBs termed double DSBs. Isolation and genome-wide mapping of the resulting fragments with single base-pair precision revealed enrichment at DSB hotspots but also a widely dispersed distribution across the genome. Spo11 prefers to cut sequences with similarity to a DNA-bending motif2, which indicates that bendability contributes to the choice of cleavage site. Moreover, fragment lengths have a periodicity of approximately (10.4n + 3) base pairs, which indicates that Spo11 favours cleavage on the same face of underwound DNA. Consistently, double DSB signals overlap and correlate with topoisomerase II-binding sites, which points to a role for topological stress and DNA crossings in break formation, and suggests a model for the formation of DSBs and double DSBs in which Spo11 traps two DNA strands. Double DSB gaps, which make up an estimated 20% of all initiation events, can account for full gene conversion events that are independent of both Msh2-dependent heteroduplex repair3,4 and the MutLγ endonuclease4. Because non-homologous gap repair results in deletions, and ectopically re-integrated double DSB fragments result in insertions, the formation of double DSBs is a potential source of evolutionary diversity and pathogenic germline aberrations. Meiotic recombination in yeast is not only initiated by single break sites, but also caused by closely spaced Spo11-dependent double-stranded DNA breaks that create chromosomal gaps.

The majority of homeobox genes are highly conserved across animals, but the eutherian-specific ETCHbox genes, embryonically expressed and highly divergent duplicates of CRX, are a notable exception. Here we compare the ETCHbox genes of 34 mammalian species, uncovering dynamic patterns of gene loss and tandem duplication, including the presence of a large tandem array of LEUTX loci in the genome of the European rabbit (Oryctolagus cuniculus). Despite extensive gene gain and loss, all sampled species possess at least two ETCHbox genes, suggesting their collective role is indispensable. We find evidence for positive selection and show that TPRX1 and TPRX2 have been the subject of repeated gene conversion across the Boreoeutheria, homogenising their sequences and preventing divergence, especially in the homeobox region. Together, these results are consistent with a model where mammalian ETCHbox genes are dynamic in evolution due to functional overlap, yet have collective indispensable roles.

BACKGROUND: Plasma tumor DNA fraction is prognostic in metastatic cancers. This could improve risk stratification before commencing a new treatment. We hypothesized that a second sample collected after one cycle of treatment could refine outcome prediction of patients identified as poor prognosis based on plasma DNA collected pre-treatment. PATIENTS AND METHODS: Plasma DNA [128 pre-treatment, 134 cycle 2 day 1 (C2D1), and 49 progression] from 151 chemotherapy-naive metastatic castration-resistant prostate cancer (mCRPC) patients in a phase II study of abiraterone acetate (NCT01867710) were subjected to custom targeted next-generation sequencing covering exons of these genes: TP53, AR, RB1, PTEN, PIK3CA, BRCA1, BRCA2, ATM, CDK12, CHEK2, FANCA HDAC2 and PALB2. We also captured 1500 pan-genome regions enriched for single nucleotide polymorphisms to allow detection of tumor DNA using the rolling B-allele method. We tested associations with overall survival (OS) and progression-free survival (PFS). RESULTS: Plasma tumor DNA detection was associated with shorter OS [hazard ratio (HR): 2.89, 95% confidence intervals (CI): 1.77-4.73, P ≤ 0.0001] and PFS (HR: 2.05; 95% CI: 1.36-3.11, P < 0.001). Using a multivariable model including plasma tumor DNA, patients who had a TP53, RB1 or PTEN gene alteration pre-treatment and at C2D1 had a significantly shorter OS than patients with no alteration at either time point (TP53: HR 7.13, 95% CI 2.37-21.47, P < 0.001; RB1: HR 6.24, 95% CI 1.97-19.73, P = 0.002; PTEN: HR 11.9, 95% CI 3.6-39.34, P < 0.001). Patients who were positive pre-treatment and converted to undetectable had no evidence of a difference in survival compared with those who were undetectable pre-treatment (P = 0.48, P = 0.43, P = 0.5, respectively). Progression samples harbored AR gain in all patients who had gain pre-treatment (9/49) and de novo AR somatic point mutations were detected in 8/49 patients. CONCLUSIONS: Plasma gene testing after one cycle treatment refines prognostication and could provide an early indication of treatment benefit.

The present study aims to analyze the effect of apricot kernels' extract (AKE) and amygdalin (AMY) on bleomycin-induced genetic alternations. Five endpoints were analyzed: cell survival, Ty1 retrotransposition, mitotic gene conversion in the trp-5 locus, reverse point mutations in ilv1-92 allele, and mitotic crossing-over in the ade2 locus. The present work provides the first experimental evidence that bleomycin induces Ty1 retrotransposition in Saccharomyces cerevisiae. New data is obtained that the degree of DNA protection of AMY and AKE depends on the studied genetic event. AKE has been found to provide significant protection against bleomycin-induced Ty1 retrotransposition due to better-expressed antioxidant potential. On the other side, AMY better-expressed protection against bleomycin-induced mitotic gene conversion and reverse mutations may be attributed to the activation of the repair enzymes.

Analysis of human cancer genome sequences has revealed specific mutational signatures associated with BRCA1-deficient tumors, but the underlying mechanisms remain poorly understood. Here, we show that one-ended DNA double strand breaks (DSBs) converted from CRISPR/Cas9-induced nicks by DNA replication, not two-ended DSBs, cause more characteristic chromosomal aberrations and micronuclei in Brca1-deficient cells than in wild-type cells. BRCA1 is required for efficient homologous recombination of these nick-converted DSBs and suppresses bias towards long tract gene conversion and tandem duplication (TD) mediated by two-round strand invasion in a replication strand asymmetry. However, aberrant repair of these nick-converted one-ended DSBs, not that of two-ended DSBs in Brca1-deficient cells, generates mutational signatures such as small indels with microhomology (MH) at the junctions, translocations and small MH-mediated TDs, resembling those in BRCA1-deficient tumors. These results suggest a major contribution of DNA nicks to mutational signatures associated with BRCA1 deficiency in cancer and the underlying mechanisms.

Protein coding genes terminate with one of three stop codons (TAA, TGA, or TAG) that, like synonymous codons, are not employed equally. With TGA and TAG having identical nucleotide content, analysis of their differential usage provides an unusual window into the forces operating on what are ostensibly functionally identical residues. Across genomes and between isochores within the human genome, TGA usage increases with G + C content but, with a common G + C → A + T mutation bias, this cannot be explained by mutation bias-drift equilibrium. Increased usage of TGA in G + C-rich genomes or genomic regions is also unlikely to reflect selection for the optimal stop codon, as TAA appears to be universally optimal, probably because it has the lowest read-through rate. Despite TAA being favored by selection and mutation bias, as with codon usage bias G + C pressure is the prime determinant of between-species TGA usage trends. In species with strong G + C-biased gene conversion (gBGC), such as mammals and birds, the high usage and conservation of TGA is best explained by an A + T → G + C repair bias. How to explain TGA enrichment in other G + C-rich genomes is less clear. Enigmatically, across bacterial and archaeal species and between human isochores TAG usage is mostly unresponsive to G + C pressure. This unresponsiveness we dub the TAG paradox as currently no mutational, selective, or gBGC model provides a well-supported explanation. That TAG does increase with G + C usage across eukaryotes makes the usage elsewhere yet more enigmatic. We suggest resolution of the TAG paradox may provide insights into either an unknown but common selective preference (probably at the DNA/RNA level) or an unrecognized complexity to the action of gBGC.

Alternative synonymous codons are often used at unequal frequencies. Classically, studies of such codon usage bias (CUB) attempted to separate the impact of neutral from selective forces by assuming that deviations from a predicted neutral equilibrium capture selection. However, GC-biased gene conversion (gBGC) can also cause deviation from a neutral null. Alternatively, selection has been inferred from CUB in highly expressed genes, but the accuracy of this approach has not been extensively tested, and gBGC can interfere with such extrapolations (e.g., if expression and gene conversion rates covary). It is therefore critical to examine deviations from a mutational null in a species with no gBGC. To achieve this goal, we implement such an analysis in the highly AT rich genome of Dictyostelium discoideum, where we find no evidence of gBGC. We infer neutral CUB under mutational equilibrium to quantify “adaptive codon preference,” a nontautologous genome wide quantitative measure of the relative selection strength driving CUB. We observe signatures of purifying selection consistent with selection favoring adaptive codon preference. Preferred codons are not GC rich, underscoring the independence from gBGC. Expression-associated “preference” largely matches adaptive codon preference but does not wholly capture the influence of selection shaping patterns across all genes, suggesting selective constraints associated specifically with high expression. We observe patterns consistent with effects on mRNA translation and stability shaping adaptive codon preference. Thus, our approach to quantifying adaptive codon preference provides a framework for inferring the sources of selection that shape CUB across different contexts within the genome.

Gerbils are a subfamily of rodents living in arid regions of Asia and Africa. Recent studies have shown that several gerbil species have unusual amino acid changes in the PDX1 protein, a homeodomain transcription factor essential for pancreatic development and β-cell function. These changes were linked to strong GC-bias in the genome that may be caused by GC-biased gene conversion, and it has been hypothesized that this caused accumulation of deleterious changes. Here we use two approaches to examine if the unusual changes are adaptive or deleterious. First, we compare PDX1 protein sequences between 38 rodents to test for association with habitat. We show the PDX1 homeodomain is almost totally conserved in rodents, apart from gerbils, regardless of habitat. Second, we use ectopic gene overexpression and gene editing in cell culture to compare functional properties of PDX1 proteins. We show that the divergent gerbil PDX1 protein inefficiently binds an insulin gene promoter and ineffectively regulates insulin expression in response to high glucose in rat cells. The protein has, however, retained the ability to regulate some other β-cell genes. We suggest that during the evolution of gerbils, the selection-blind process of biased gene conversion pushed fixation of mutations adversely affecting function of a normally conserved homeodomain protein. We argue these changes were not entirely adaptive and may be associated with metabolic disorders in gerbil species on high carbohydrate diets. This unusual pattern of molecular evolution could have had a constraining effect on habitat and diet choice in the gerbil lineage.

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.

The evolution of tRNA multigene families remains poorly understood, exhibiting unusual phenomena such as functional conversions of tRNA genes through anticodon shift substitutions. We improved FlyBase tRNA gene annotations from twelve Drosophila species, incorporating previously identified ortholog sets to compare substitution rates across tRNA bodies at single-site and base-pair resolution. All rapidly evolving sites fell within the same metal ion-binding pocket that lies at the interface of the two major stacked helical domains. We applied our tRNA Structure–Function Mapper (tSFM) method independently to each Drosophila species and one outgroup species Musca domestica and found that, although predicted tRNA structure–function maps are generally highly conserved in flies, one tRNA Class-Informative Feature (CIF) within the rapidly evolving ion-binding pocket—Cytosine 17 (C17), ancestrally informative for lysylation identity—independently gained asparaginylation identity and substituted in parallel across tRNAAsn paralogs at least once, possibly multiple times, during evolution of the genus. In D. melanogaster, most tRNALys and tRNAAsn genes are co-arrayed in one large heterologous gene cluster, suggesting that heterologous gene conversion as well as structural similarities of tRNA-binding interfaces in the closely related asparaginyl-tRNA synthetase (AsnRS) and lysyl-tRNA synthetase (LysRS) proteins may have played a role in these changes. A previously identified Asn-to-Lys anticodon shift substitution in D. ananassae may have arisen to compensate for the convergent and parallel gains of C17 in tRNAAsn paralogs in that lineage. Our results underscore the functional and evolutionary relevance of our tRNA structure–function map predictions and illuminate multiple genomic and structural factors contributing to rapid, parallel and compensatory evolution of tRNA multigene families.

Diversification of the avian primary immunoglobulin (Ig) repertoire is achieved in developing B cells by somatic hypermutation (SHM) and gene conversion (GCV). GCV is a type of homologous recombination that unidirectionally transfers segments of Ig pseudogenes to Ig variable domains. It is regulated by epigenetic mechanisms like histone modifications, but the role of DNA methylation remains unclear. Here, we demonstrate that the chicken B‐cell line DT40 lacking TET3, a member of the TET (Ten‐eleven translocation) family dioxygenases that facilitate DNA demethylation, exhibited a marked reduction in GCV activity in Ig variable regions. This was accompanied by a drop in the bulk levels of 5‐hydroxymethylcytosine, an oxidized derivative of 5‐methylcytosine, whereas TET1‐deficient or TET2‐deficient DT40 strains did not exhibit such effects. Deletion of TET3 caused little effects on the expression of proteins required for SHM and GCV, but induced hypermethylation in some Ig pseudogene templates. Notably, the enhanced methylation occurred preferably on non‐CpG cytosines. Disruption of both TET1 and TET3 significantly inhibited the expression of activation‐induced cytidine deaminase (AID), an essential player in Ig diversification. These results uncover unique roles of TET proteins in avian Ig diversification, highlighting the potential importance of TET3 in maintaining hypomethylation In Ig pseudogenes., We demonstrated the chicken B‐cell line DT40 lacking TET3, a member of the TET family proteins that facilitate DNA demethylation, exhibited a marked reduction in immunoglobulin (Ig) gene diversification. This was accompanied by enhanced methylation on non‐CpG cytosines in Ig pseudogene templates. These results uncover unique roles of TET3 in avian Ig diversification likely via maintaining hypomethylation in Ig pseudogenes.

Proximal spinal muscular atrophy (SMA), a leading genetic cause of infant death worldwide, is an early-onset motor neuron disease characterized by loss of α-motor neurons and associated muscle atrophy. SMA is caused by deletion or other disabling mutations of survival motor neuron 1 (SMN1) but retention of one or more copies of the paralog SMN2. Within the SMA population, there is substantial variation in SMN2 copy number (CN); in general, those individuals with SMA who have a high SMN2 CN have a milder disease. Because SMN2 functions as a disease modifier, its accurate CN determination may have clinical relevance. In this study, we describe the development of array digital PCR (dPCR) to quantify SMN1 and SMN2 CNs in DNA samples using probes that can distinguish the single nucleotide difference between SMN1 and SMN2 in exon 8. This set of dPCR assays can accurately and reliably measure the number of SMN1 and SMN2 copies in DNA samples. In a cohort of SMA patient-derived cell lines, the assay confirmed a strong inverse correlation between SMN2 CN and disease severity. We can detect SMN1-SMN2 gene conversion events in DNA samples by comparing CNs at exon 7 and exon 8. Partial deletions of SMN1 can also be detected with dPCR by comparing CNs at exon 7 or exon 8 with those at intron 1. Array dPCR is a practical technique to determine, accurately and reliably, SMN1 and SMN2 CNs from SMA samples as well as identify gene conversion events and partial deletions of SMN1.

Leucine-rich repeat proteins and antimicrobial peptides are the key components of the innate immune response to

The single-cell parasite Trypanosoma brucei causes the fatal disease human African trypanosomiasis and is able to colonize the blood, fat, skin, and central nervous system. Trypanosomes survive in the mammalian host owing to a dense protective protein coat that consists of a single-variant surface glycoprotein species.

The human MSY ampliconic region is mainly composed of large duplicated sequences that are organized in eight palindromes (termed P1–P8), and may undergo arm-to-arm gene conversion. Although the importance of these elements is widely recognized, their evolutionary dynamics are still nuanced. Here, we focused on the P8 palindrome, which shows a complex evolutionary history, being involved in intra- and inter-chromosomal gene conversion. To disclose its evolutionary complexity, we performed a high-depth (50×) targeted next-generation sequencing of this element in 157 subjects belonging to the most divergent lineages of the Y chromosome tree. We found a total of 72 polymorphic paralogous sequence variants that have been exploited to identify 41 Y-Y gene conversion events that occurred during recent human history. Through our analysis, we were able to categorize P8 arms into three portions, whose molecular diversity was modelled by different evolutionary forces. Notably, the outer region of the palindrome is not involved in any gene conversion event and evolves exclusively through the action of mutational pressure. The inner region is affected by Y-Y gene conversion occurring at a rate of 1.52 × 10−5 conversions/base/year, with no bias towards the retention of the ancestral state of the sequence. In this portion, GC-biased gene conversion is counterbalanced by a mutational bias towards AT bases. Finally, the middle region of the arms, in addition to intra-chromosomal gene conversion, is involved in X-to-Y gene conversion (at a rate of 6.013 × 10−8 conversions/base/year) thus being a major force in the evolution of the VCY/VCX gene family.

The strength of selection varies among populations and across the genome, but the determinants of efficacy of selection remain unclear. In this study, we used whole-genome sequencing data from 467 Boechera stricta accessions to quantify the strength of selection and characterize the pattern of local adaptation. We found low genetic diversity on 0-fold degenerate sites and conserved non-coding sites, indicating functional constraints on these regions. The estimated distribution of fitness effects and the proportion of fixed substitutions suggest relaxed negative and positive selection in B. stricta. Among the four population groups, the NOR and WES groups have smaller effective population size (Ne), higher proportions of effectively neutral sites, and lower rates of adaptive evolution compared with UTA and COL groups, reflecting the effect of Ne on the efficacy of natural selection. We also found weaker selection on GC-biased sites compared with GC-conservative (unbiased) sites, suggested that GC-biased gene conversion has affected the strength of selection in B. stricta. We found mixed evidence for the role of the recombination rate on the efficacy of selection. The positive and negative selection was stronger in high-recombination regions compared with low-recombination regions in COL but not in other groups. By scanning the genome, we found different subsets of selected genes suggesting differential adaptation among B. stricta groups. These results show that differences in effective population size, nucleotide composition, and recombination rate are important determinants of the efficacy of selection. This study enriches our understanding of the roles of natural selection and local adaptation in shaping genomic variation.

Quantifying and comparing the amount of adaptive evolution among different species is key to understanding how evolution works. Previous studies have shown differences in adaptive evolution across species; however, their specific causes remain elusive. Here, we use improved modeling of weakly deleterious mutations and the demographic history of the outgroup species and ancestral population and estimate that at least 20% of nonsynonymous substitutions between humans and an outgroup species were fixed by positive selection. This estimate is much higher than previous estimates, which did not correct for the sizes of the outgroup species and ancestral population. Next, we jointly estimate the proportion and selection coefficient (p+ and s+, respectively) of newly arising beneficial nonsynonymous mutations in humans, mice, and Drosophila melanogaster by examining patterns of polymorphism and divergence. We develop a novel composite likelihood framework to test whether these parameters differ across species. Overall, we reject a model with the same p+ and s+ of beneficial mutations across species and estimate that humans have a higher p+s+ compared with that of D. melanogaster and mice. We show that this result cannot be caused by biased gene conversion or hypermutable CpG sites. We discuss possible biological explanations that could generate the observed differences in the amount of adaptive evolution across species.

In the last 20 years, severe wheat stem rust outbreaks have been recorded in Africa, Europe, and Central Asia. This previously well controlled disease, caused by the fungus Puccinia graminis f. sp. tritici, has reemerged as a major threat to wheat cultivation. The stem rust (Sr) resistance gene Sr22 encodes a nucleotide-binding and leucine-rich repeat receptor which confers resistance to the highly virulent African stem rust isolate Ug99. Here, we show that the Sr22 gene is conserved among grasses in the Triticeae and Poeae lineages. Triticeae species contain syntenic loci with single-copy orthologs of Sr22 on chromosome 7, except Hordeum vulgare, which has experienced major expansions and rearrangements at the locus. We also describe 14 Sr22 sequence variants obtained from both Triticum boeoticum and the domesticated form of this species, T. monococcum, which have been postulated to encode both functional and nonfunctional Sr22 alleles. The nucleotide sequence analysis of these alleles identified historical sequence exchange resulting from recombination or gene conversion, including breakpoints within codons, which expanded the coding potential at these positions by introduction of nonsynonymous substitutions. Three Sr22 alleles were transformed into wheat cultivar Fielder and two postulated resistant alleles from Schomburgk (hexaploid wheat introgressed with T. boeoticum segment carrying Sr22) and T. monococcum accession PI190945, respectively, conferred resistance to P. graminis f. sp. tritici race TTKSK, thereby unequivocally confirming Sr22 effectiveness against Ug99. The third allele from accession PI573523, previously believed to confer susceptibility, was confirmed as nonfunctional against Australian P. graminis f. sp. tritici race 98-1,2,3,5,6. [Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.

Significance The male-specific Y chromosome harbors genes important for sperm production. Because Y is repetitive, its DNA sequence was deciphered for only a few species, and its evolution remains elusive. Here we compared the Y chromosomes of great apes (human, chimpanzee, bonobo, gorilla, and orangutan) and found that many of their repetitive sequences and multicopy genes were likely already present in their common ancestor. Y repeats had increased intrachromosomal contacts, which might facilitate preservation of genes and gene regulatory elements. Chimpanzee and bonobo, experiencing high sperm competition, underwent many DNA changes and gene losses on the Y. Our research is significant for understanding the role of the Y chromosome in reproduction of nonhuman great apes, all of which are endangered., The mammalian male-specific Y chromosome plays a critical role in sex determination and male fertility. However, because of its repetitive and haploid nature, it is frequently absent from genome assemblies and remains enigmatic. The Y chromosomes of great apes represent a particular puzzle: their gene content is more similar between human and gorilla than between human and chimpanzee, even though human and chimpanzee share a more recent common ancestor. To solve this puzzle, here we constructed a dataset including Ys from all extant great ape genera. We generated assemblies of bonobo and orangutan Ys from short and long sequencing reads and aligned them with the publicly available human, chimpanzee, and gorilla Y assemblies. Analyzing this dataset, we found that the genus Pan, which includes chimpanzee and bonobo, experienced accelerated substitution rates. Pan also exhibited elevated gene death rates. These observations are consistent with high levels of sperm competition in Pan. Furthermore, we inferred that the great ape common ancestor already possessed multicopy sequences homologous to most human and chimpanzee palindromes. Nonetheless, each species also acquired distinct ampliconic sequences. We also detected increased chromatin contacts between and within palindromes (from Hi-C data), likely facilitating gene conversion and structural rearrangements. Our results highlight the dynamic mode of Y chromosome evolution and open avenues for studies of male-specific dispersal in endangered great ape species.

The aim of this study was to investigate a family with nonhomologous sequence recombination of HBA1 and HBA2 genes and provide a favorable basis for genetic counseling and eugenics. Peripheral blood of family members was collected. Hematological parameters were determined by an automated cell counter and hemoglobin (Hb) analysis was performed using high performance liquid chromatography (HPLC). Villus samples were taken for prenatal diagnosis (PND). Gap-polymerase chain reaction (gap-PCR) and reverse dot-blot were used for thalassemia genotyping. DNA sequencing was used to analyze the gene sequence of HBA1 (α1-globin) and HBA2 (α2-globin). The nonhomologous sequence recombination allelic variant of HBA1 and HBA2 genes were identified, namely, a gene conversion on the HBA2 gene called α12 (HBA12). The α12 allele consists primarily of the HBA2 gene sequence except for a segment of the IVS-II in which HBA2-specific sequences have been replaced by HBA1-specific sequences. The following genotypes were observed: - -SEA/αα12 (Southeast Asian deletion), αα/αα12 and αQSα/αα12 (Hb Quong Sze or Hb QS; HBA2: c.377T>C), and all manifested as small cell hypochromic anemia. To find the α12 allele in the Chinese population and clarify the influence of the α12 allele and its common inheritance with abnormal Hb and α-thalassemia (α-thal) on α-globin gene expression can help guide clinical diagnosis and genetic counseling.

Synonymous substitution rates in plant mitochondrial genomes vary by orders of magnitude among species, whereas synonymous rates among genes within a genome are generally consistent. Exceptionally, genes within the Ajuga reptans (Lamiaceae) mitochondrial genome exhibit unprecedented intragenomic heterogeneity in synonymous sequence divergence, but the biological mechanisms underlying this rate variation remain unclear. We tracked the origin and evolutionary trajectory of mitochondrial rate variations by dense sampling in Ajugoideae and found differences in the timing and magnitude of rate acceleration for particular genes. The most divergent genes accelerated earlier, retained a high rate across Ajugoideae, and are generally devoid of RNA editing, whereas moderately diverged genes accelerated later and retained relatively higher RNA editing frequency. The acceleration of mutation rates correlates with increased guanine-cytosine (GC) content, suggesting a key role for GC-biased gene conversion and/or repair after the breakage of ancestral gene clusters.

Reduction of fitness due to deleterious mutations imposes a limit to adaptive evolution. By characterizing features that influence this genetic load we may better understand constraints on responses to both natural and human-mediated selection. Here, using whole-genome, transcriptome, and methylome data from >600 Arabidopsis thaliana individuals, we set out to identify important features influencing selective constraint. Our analyses reveal that multiple factors underlie the accumulation of maladaptive mutations, including gene expression level, gene network connectivity, and gene-body methylation. We then focus on a feature with major effect, nucleotide composition. The ancestral vs. derived status of segregating alleles suggests that GC-biased gene conversion, a recombination-associated process that increases the frequency of G and C nucleotides regardless of their fitness effects, shapes sequence patterns in A. thaliana. Through estimation of mutational effects, we present evidence that biased gene conversion hinders the purging of deleterious mutations and contributes to a genome-wide signal of decreased efficacy of selection. By comparing these results to two outcrossing relatives, Arabidopsis lyrata and Capsella grandiflora, we find that protein evolution in A. thaliana is as strongly affected by biased gene conversion as in the outcrossing species. Last, we perform simulations to show that natural levels of outcrossing in A. thaliana are sufficient to facilitate biased gene conversion despite increased homozygosity due to selfing. Together, our results show that even predominantly selfing taxa are susceptible to biased gene conversion, suggesting that it may constitute an important constraint to adaptation among plant species.

The rate and spectrum of spontaneous mutations are critical parameters in basic and applied biology because they dictate the pace and character of genetic variation introduced into populations, which is a prerequisite for evolution. We use a mutation–accumulation approach to estimate mutation parameters from whole-genome sequence data from multiple genotypes from multiple populations of Daphnia magna, an ecological and evolutionary model system. We report extremely high base substitution mutation rates (µ-n,bs = 8.96 × 10−9/bp/generation [95% CI: 6.66–11.97 × 10−9/bp/generation] in the nuclear genome and µ-m,bs = 8.7 × 10−7/bp/generation [95% CI: 4.40–15.12 × 10−7/bp/generation] in the mtDNA), the highest of any eukaryote examined using this approach. Levels of intraspecific variation based on the range of estimates from the nine genotypes collected from three populations (Finland, Germany, and Israel) span 1 and 3 orders of magnitude, respectively, resulting in up to a ∼300-fold difference in rates among genomic partitions within the same lineage. In contrast, mutation spectra exhibit very consistent patterns across genotypes and populations, suggesting the mechanisms underlying the mutational process may be similar, even when the rates at which they occur differ. We discuss the implications of high levels of intraspecific variation in rates, the importance of estimating gene conversion rates using a mutation–accumulation approach, and the interacting factors influencing the evolution of mutation parameters. Our findings deepen our knowledge about mutation and provide both challenges to and support for current theories aimed at explaining the evolution of the mutation rate, as a trait, across taxa.

Vision represents an excellent model for studying adaptation, given the genotype-to-phenotype-map that has been characterized in a number of taxa. Fish possess a diverse range of visual sensitivities and adaptations to underwater light making them an excellent group to study visual system evolution. In particular, some speciose but understudied lineages can provide a unique opportunity to better understand aspects of visual system evolution such as opsin gene duplication and neofunctionalization. In this study, we characterized the visual system of Neotropical Characiformes, which is the result of several spectral tuning mechanisms acting in concert including gene duplications and losses, gene conversion, opsin amino acid sequence and expression variation, and A1/A2-chromophore shifts. The Characiforms we studied utilize three cone opsin classes (SWS2, RH2, LWS) and a rod opsin (RH1). However, the characiform’s entire opsin gene repertoire is a product of dynamic evolution by opsin gene loss (SWS1, RH2) and duplication (LWS, RH1). The LWS- and RH1-duplicates originated from a teleost specific whole-genome duplication as well as characiform-specific duplication events. Both LWS-opsins exhibit gene conversion and, through substitutions in key tuning sites, one of the LWS-paralogs has acquired spectral sensitivity to green light. These sequence changes suggest reversion and parallel evolution of key tuning sites. In addition, characiforms exhibited species-specific differences in opsin expression. Finally, we found interspecific and intraspecific variation in the use of A1/A2-chromophores correlating with the light environment. These multiple mechanisms may be a result of the highly diverse visual environments where Characiformes have evolved.

Monoclonal antibodies (mAbs) are widely utilized as therapeutic drugs for various diseases, such as cancer, autoimmune diseases, and infectious diseases. Using the avian-derived B cell line DT40, we previously developed an antibody display technology, namely, the ADLib system, which rapidly generates antigen-specific mAbs. Here, we report the development of a human version of the ADLib system and showcase the streamlined generation and optimization of functional human mAbs. Tailored libraries were first constructed by replacing endogenous immunoglobulin genes with designed human counterparts. From these libraries, clones producing full-length human IgGs against distinct antigens can be isolated, as exemplified by the selection of antagonistic mAbs. Taking advantage of avian biology, effective affinity maturation was achieved in a straightforward manner by seamless diversification of the parental clones into secondary libraries followed by single-cell sorting, quickly affording mAbs with improved affinities and functionalities. Collectively, we demonstrate that the human ADLib system could serve as an integrative platform with unique diversity for rapid de novo generation and optimization of therapeutic or diagnostic antibody leads. Furthermore, our results suggest that libraries can be constructed by introducing exogenous genes into DT40 cells, indicating that the ADLib system has the potential to be applied for the rapid and effective directed evolution and optimization of proteins in various fields beyond biomedicine.

Shwachman-Diamond syndrome characterized by metaphyseal dysplasia, pancreatic insufficiency, and pancytopenia is caused by biallelic mutations in SBDS. Gene conversion between SBDS and its pseudogene SBDSP1 is the major cause. Here, we report two unrelated patients with Shwachman-Diamond syndrome who were shown to be compound heterozygotes for relatively frequent pathogenic alleles (the 258+2T>C allele and another allele composed of 183-184TA>CT and 201A>G) using an established polymerase chain reaction sequencing assay with SBDS-specific primers. Exome analysis of the patients showed discrepant results: 258+2T>C with variant allele frequency around 0.85, and no variants detected for the 183-184TA>CT allele. Parental exome analysis of the two families further supported this notion. Confronted with two patients with an unexpected segregation pattern, we performed a transcriptome analysis of peripheral blood-derived mRNA to demonstrate that the results were compatible with those obtained using SBDS-specific PCR primers. Both alleles could be accounted for by gene conversion events. The diagnostic discrepancy can be accounted for by a decreased efficiency in the computational mapping of the reads with 183-184TA>CT and 201A>G to the reference sequence of the SBDS locus during exome analysis. This report highlights the pitfall of exome analysis for genes with pseudogenes, such as SBDS and the alternative use of RNA-seq is recommended to circumvent this problem.

Recombination increases the local GC-content in genomic regions through GC-biased gene conversion (gBGC). The recent discovery of a large genomic region with extreme GC-content in the fat sand rat Psammomys obesus provides a model to study the effects of gBGC on chromosome evolution. Here, we compare the GC-content and GC-to-AT substitution patterns across protein-coding genes of four gerbil species and two murine rodents (mouse and rat). We find that the known high-GC region is present in all the gerbils, and is characterized by high substitution rates for all mutational categories (AT-to-GC, GC-to-AT, and GC-conservative) both at synonymous and nonsynonymous sites. A higher AT-to-GC than GC-to-AT rate is consistent with the high GC-content. Additionally, we find more than 300 genes outside the known region with outlying values of AT-to-GC synonymous substitution rates in gerbils. Of these, over 30% are organized into at least 17 large clusters observable at the megabase-scale. The unusual GC-skewed substitution pattern suggests the evolution of genomic regions with very high recombination rates in the gerbil lineage, which can lead to a runaway increase in GC-content. Our results imply that rapid evolution of GC-content is possible in mammals, with gerbil species providing a powerful model to study the mechanisms of gBGC.

Glutamine synthetase is an essential enzyme in ammonium assimilation and glutamine biosynthesis. The Haloferax mediterranei genome has two other glnA-type genes (glnA2 and glnA3) in addition to the glutamine synthetase gene glnA. To determine whether the glnA2 and glnA3 genes can replace glnA in nitrogen metabolism, we generated deletion mutants of glnA. The glnA deletion mutants could not be generated in a medium without glutamine, and thus, glnA is an essential gene in H. mediterranei. The glnA deletion mutant was achieved by adding 40 mM glutamine to the selective medium. This conditional HM26-ΔglnA mutant was characterised with different approaches in the presence of distinct nitrogen sources and nitrogen starvation. Transcriptomic analysis was performed to compare the expression profiles of the strains HM26-ΔglnA and HM26 under different growth conditions. The glnA deletion did not affect the expression of glnA2, glnA3 and nitrogen assimilation genes under nitrogen starvation. Moreover, the results showed that glnA, glnA2 and glnA3 were not expressed under the same conditions. These results indicated that glnA is an essential gene for H. mediterranei and, therefore, glnA2 and glnA3 cannot replace glnA in the conditions analysed. This work was funded by MICINN Grant Number BIO2013-42921P (to MJB), Generalitat Valenciana Grant Number ACIF/2018/200 (to GP) and Universidad de Alicante (VIGROB-016).

Gene conversion is a process of transferring genetic material from one homologous sequence to another. Most reported gene conversions are meiotic although mitotic gene conversion is also described. When using CRISPR/Cas9 to target the human hemoglobin subunit beta (HBB) gene, hemoglobin subunit delta (HBD) gene footprints were observed in HBB gene. However, it is unclear whether these were the results of gene conversion or PCR-mediated sequence shuffling between highly homologous sequences. Here we provide evidence that the HBD footprints in HBB were indeed results of gene conversion. We demonstrated that the CRISPR/Cas9 facilitated unidirectional sequence transfer from the homologous gene without double-strand breaks (DSB) to the one with DSBs, and showed that the rates of HBD footprint in HBB were positively correlated to the HBB insertion and deletion rates. We further showed that when targeting HBD gene, HBB footprints could also be observed in HBD gene. The mitotic gene conversion was observed not only in immortalized HEK293T cells, but also in human primary cells. Our work reveals mitotic gene conversion as an often overlooked effect of CRISPR/Cas9-mediated genome editing.

Lyme disease, also known as Lyme borreliosis, is the most common tick-transmitted disease in the Northern Hemisphere. The disease is caused by the bacterial spirochete Borrelia burgdorferi and other related Borrelia species. One of the many fascinating features of this unique pathogen is an elaborate system for antigenic variation, whereby the sequence of the surface-bound lipoprotein VlsE is continually modified through segmental gene conversion events. This perpetual changing of the guard allows the pathogen to remain one step ahead of the acquired immune response, enabling persistent infection. Accordingly, the vls locus is the most evolutionarily diverse genetic element in Lyme disease–causing borreliae. Small stretches of information are transferred from a series of silent cassettes in the vls locus to generate an expressed mosaic vlsE gene version that contains genetic information from several different silent cassettes, resulting in ∼10(40) possible vlsE sequences. Yet, despite its extreme evolutionary flexibility, the locus has rigidly conserved structural features. These include a telomeric location of the vlsE gene, an inverse orientation of vlsE and the silent cassettes, the presence of nearly perfect inverted repeats of ∼100 bp near the 5′ end of vlsE, and an exceedingly high concentration of G runs in vlsE and the silent cassettes. We discuss the possible roles of these evolutionarily conserved features, highlight recent findings from several studies that have used next-generation DNA sequencing to unravel the switching process, and review advances in the development of a mini-vls system for genetic manipulation of the locus.

TaON and Ta

The major processes in carcinogenesis include the inactivation of tumor-suppressor genes (TSGs). Although Knudson's two-hit model requires two independent inactivating mutations, perhaps more frequently, a TSG inactivation can occur through a loss of heterozygosity (LOH) of an inactivating mutation. Deletion and uniparental disomy (UPD) have been well documented as LOH mechanisms, but the role of gene conversion is poorly understood. Here, we developed a simple algorithm to detect somatic gene conversion from short-read sequencing data. We applied it to 6285 cancer patient samples, from which 4978 somatic mutations that underwent gene conversion to achieve LOH were found. This number accounted for 14.8% of the total LOH mutations. We further showed that LOH by gene conversion was enriched in TSGs compared with non-TSG genes, showing a significant contribution of gene conversion to carcinogenesis.

Genotyping by sequencing (GBS) is an emerging technology to rapidly call an abundance of Single Nucleotide Polymorphisms (SNPs) using genome sequencing technology. Several different methodologies and approaches have recently been established, most of these relying on a specific preparation of data. Here we describe our GBS-pipeline, which uses high coverage reads from two parents and low coverage reads from their double haploid offspring to call SNPs on a large scale. The upside of this approach is the high resolution and scalability of the method.

Defects in BRCA1, BRCA2 and other genes of the homology-dependent DNA repair (HR) pathway cause an elevated rate of mutagenesis, eliciting specific mutation patterns including COSMIC signature SBS3. Using genome sequencing of knock-out cell lines we show that Y family translesion synthesis (TLS) polymerases contribute to the spontaneous generation of base substitution and short insertion/deletion mutations in BRCA1 deficient cells, and that TLS on DNA adducts is increased in BRCA1 and BRCA2 mutants. The inactivation of 53BP1 in BRCA1 mutant cells markedly reduces TLS-specific mutagenesis, and rescues the deficiency of template switch–mediated gene conversions in the immunoglobulin V locus of BRCA1 mutant chicken DT40 cells. 53BP1 also promotes TLS in human cellular extracts in vitro. Our results show that HR deficiency–specific mutagenesis is largely caused by TLS, and suggest a function for 53BP1 in regulating the choice between TLS and error-free template switching in replicative DNA damage bypass., Loss of BRCA1 or BRCA2 results in genomic instability; however most studies have focused on the role of these proteins in double-strand break repair. Here the authors coupled cell line genetics and whole genome sequencing to investigate the formation of base substitutions and short indels in BRCA1-deficient cells, revealing a role for translesion DNA synthesis regulated by 53BP1.

Y chromosomes across diverse species convergently evolve a gene-poor, heterochromatic organization enriched for duplicated genes, LTR retrotransposable elements, and satellite DNA. Sexual antagonism and a loss of recombination play major roles in the degeneration of young Y chromosomes. However, the processes shaping the evolution of mature, already degenerated Y chromosomes are less well-understood. Because Y chromosomes evolve rapidly, comparisons between closely related species are particularly useful. We generated de novo long read assemblies complemented with cytological validation to reveal Y chromosome organization in three closely related species of the Drosophila simulans complex, which diverged only 250,000 years ago and share >98% sequence identity. We find these Y chromosomes are divergent in their organization and repetitive DNA composition and discover new Y-linked gene families whose evolution is driven by both positive selection and gene conversion. These Y chromosomes are also enriched for large deletions, suggesting that the repair of double-strand breaks on Y chromosomes may be biased toward microhomology-mediated end joining over canonical non-homologous end-joining. We propose that this repair mechanism generally contributes to the convergent evolution of Y chromosome organization.

Blue cone monochromacy (BCM) is an X-linked retinal disorder characterized by low vision, photoaversion, and poor color discrimination. BCM is due to the lack of long-wavelength-sensitive and middle-wavelength-sensitive cone photoreceptor function and caused by mutations in the OPN1LW/OPN1MW gene cluster on Xq28. Here, we investigated the prevalence and the landscape of submicroscopic structural variants (SVs) at single-base resolution in BCM patients. We found that about one-third (n = 73) of the 213 molecularly confirmed BCM families carry an SV, most commonly deletions restricted to the OPN1LW/OPN1MW gene cluster. The structure and precise breakpoints of the SVs were resolved in all but one of the 73 families. Twenty-two families-all from the United States-showed the same SV, and we confirmed a common ancestry of this mutation. In total, 42 distinct SVs were identified, including 40 previously unreported SVs, thereby quadrupling the number of precisely mapped SVs underlying BCM. Notably, there was no "region of overlap" among these SVs. However, 90% of SVs encompass the upstream locus control region, an essential enhancer element. Its minimal functional extent based on deletion mapping in patients was refined to 358 bp. Breakpoint analyses suggest diverse mechanisms underlying SV formation as well as in one case the gene conversion-based exchange of a 142-bp deletion between opsin genes. Using parsimonious assumptions, we reconstructed the composition and copy number of the OPN1LW/OPN1MW gene cluster prior to the mutation event and found evidence that large gene arrays may be predisposed to the occurrence of SVs at this locus.

Ehrlichia ruminantium is a tick-borne intracellular pathogen of ruminants that causes heartwater, a disease present in Sub-saharan Africa, islands in the Indian Ocean and the Caribbean, inducing significant economic losses. At present, three avirulent strains of E. ruminantium (Gardel, Welgevonden and Senegal isolates) have been produced by a process of serial passaging in mammalian cells in vitro, but unfortunately their use as vaccines do not offer a large range of protection against other strains, possibly due to the genetic diversity present within the species. So far no genetic basis for virulence attenuation has been identified in any E. ruminantium strain that could offer targets to facilitate vaccine production. Virulence attenuated Senegal strains have been produced twice independently, and require many fewer passages to attenuate than the other strains. We compared the genomes of a virulent and attenuated Senegal strain and identified a likely attenuator gene, ntrX, a global transcription regulator and member of a two-component system that is linked to environmental sensing. This gene has an inverted partial duplicate close to the parental gene that shows evidence of gene conversion in different E. ruminantium strains. The pseudogenisation of the gene in the avirulent Senegal strain occurred by gene conversion from the duplicate to the parent, transferring a 4 bp deletion which is unique to the Senegal strain partial duplicate amongst the wild isolates. We confirmed that the ntrX gene is not expressed in the avirulent Senegal strain by RT-PCR. The inverted duplicate structure combined with the 4 bp deletion in the Senegal strain can explain both the attenuation and the faster speed of attenuation in the Senegal strain relative to other strains of E. ruminantium. Our results identify nrtX as a promising target for the generation of attenuated strains of E. ruminantium by random or directed mutagenesis that could be used for vaccine production.

Meiosis is undoubtedly the mechanism that underpins Mendelian genetics. Meiosis is a specialised, reductional cell division which generates haploid gametes (reproductive cells) carrying a single chromosome complement from diploid progenitor cells harbouring two chromosome sets. Through this process, the hereditary material is shuffled and distributed into haploid gametes such that upon fertilisation, when two haploid gametes fuse, diploidy is restored in the zygote. During meiosis the transient physical connection of two homologous chromosomes (one originally inherited from each parent) each consisting of two sister chromatids and their subsequent segregation into four meiotic products (gametes), is what enables genetic marker assortment forming the core of Mendelian laws. The initiating events of meiotic recombination are DNA double-strand breaks (DSBs) which need to be repaired in a certain way to enable the homologous chromosomes to find each other. This is achieved by DSB ends searching for homologous repair templates and invading them. Ultimately, the repair of meiotic DSBs by homologous recombination physically connects homologous chromosomes through crossovers. These physical connections provided by crossovers enable faithful chromosome segregation. That being said, the DSB repair mechanism integral to meiotic recombination also produces genetic transmission distortions which manifest as postmeiotic segregation events and gene conversions. These processes are non-reciprocal genetic exchanges and thus non-Mendelian.

Flowering pattern is a major taxonomic characteristic differentiating the two main subspecies of cultivated peanut (Arachis hypogaea L.). subsp. fastigiata possessing flowers on the mainstem (MSF) and a sequential flowering pattern, whereas subsp. hypogaea lacks flowers on the mainstem and exhibits an alternate flowering pattern. This character is considered the main contributor to plant architecture and the adaptability of each subgroup to specific growing conditions. Evidence indicates that flowering pattern differentiation occurred during the several thousand years of domestication and diversification in South America. However, the exact genetic mechanism that controls flowering pattern and the molecular changes that led to its historical diversification in peanut are unknown. We investigated the genetics of the flowering pattern in a recombinant inbred population of 259 lines (RILs), derivatives of an A. hypogaea and A. fastigiata cross. RILs segregated 1:1 in both the sequential/alternative and the MSF-plus/MSF-minus traits, indicating a single gene effect. Using the Axiom_Arachis2 SNP-array, MSF was mapped to a 1.7 Mbp segment on chromosome B02 of the cultivated A. hypogaea. Significant haplotype conservation was found for this locus in the USA peanut mini core collection, suggesting a possible selection upon hypogaea/fastigiata speciation. Furthermore, a candidate Terminal Flowering 1-like (AhTFL1) gene was identified within the MSF region, in which a 1492 bp deletion occurred in the fastigiata line that leads to a truncated protein product. Remapping MSF in the RIL population with the AhTFL1 deletion as a marker increased the LOD score from 53.3 to 158.8 with no recombination. The same deletion was also found to co-segregate with the phenotype in two EMS-mutagenized M2 families, suggesting a hotspot for large mutational deletion or gene conversion that may play a role in evolution. BLASTX analysis showed that the most similar homologous gene for TFL1-like in soybean is Det1, which previously was shown to control shoot determination. Sequence analysis of the TFL-1 in a series of domesticated lines showed that TFL1 was subjected to gain/loss events of the deletion, partly explaining the evolution of MSF in Arachis. Altogether, these results support the role of AhTFL-1 in peanut speciation during domestication and modern cultivation.

Muller’s ratchet is a process in which deleterious mutations are fixed irreversibly in the absence of recombination. The degeneration of the Y chromosome, and the gradual loss of its genes, can be explained by Muller’s ratchet. However, most theories consider single-copy genes, and may not be applicable to Y chromosomes, which have a number of duplicated genes in many species, which are probably undergoing concerted evolution by gene conversion. We developed a model of Muller’s ratchet to explore the evolution of the Y chromosome. The model assumes a nonrecombining chromosome with both single-copy and duplicated genes. We used analytical and simulation approaches to obtain the rate of gene loss in this model, with special attention to the role of gene conversion. Homogenization by gene conversion makes both duplicated copies either mutated or intact. The former promotes the ratchet, and the latter retards, and we ask which of these counteracting forces dominates under which conditions. We found that the effect of gene conversion is complex, and depends upon the fitness effect of gene duplication. When duplication has no effect on fitness, gene conversion accelerates the ratchet of both single-copy and duplicated genes. If duplication has an additive fitness effect, the ratchet of single-copy genes is accelerated by gene duplication, regardless of the gene conversion rate, whereas gene conversion slows the degeneration of duplicated genes. Our results suggest that the evolution of the Y chromosome involves several parameters, including the fitness effect of gene duplication by increasing dosage and gene conversion rate.

GBA variants cause the autosomal recessive Gaucher disease, and carriers are at increased risk of Parkinson’s disease (PD) and Lewy body dementia (LBD). The presence of a highly homologous nearby pseudogene (GBAP1) predisposes to a range of structural variants arising from either gene conversion or reciprocal recombination, the latter resulting in copy number gains or losses, complicating genetic testing and analysis. To date, short-read sequencing has not been able to fully resolve these or other variants in the key homology region, and targeted long-read sequencing has not previously resolved reciprocal recombinants. We present and validate two independent methods to resolve recombinant alleles and other variants in GBA: Gauchian, a novel bioinformatics tool for short-read, whole-genome sequencing data analysis, and Oxford Nanopore long-read sequencing after enrichment with appropriate PCR. The methods were concordant for 42 samples including 30 with a range of recombinants and GBAP1-related mutations, and Gauchian outperforms the GATK Best Practices pipeline. Applying Gauchian to Illumina sequencing of over 10,000 individuals from publicly available cohorts shows that copy number variants (CNVs) spanning GBAP1 are relatively common in Africans. CNV frequencies in PD and LBD are similar to controls, but gains may coexist with other mutations in patients, and a modifying effect cannot be excluded. Gauchian detects a higher frequency of GBA variants in LBD than PD, especially severe ones. These findings highlight the importance of accurate GBA mutation detection in these patients, which is possible by either Gauchian analysis of short-read whole genome sequencing, or targeted long-read sequencing.

Recombination can occur either as a result of crossover or gene conversion events. Population genetic methods for inferring the rate of recombination from patterns of linkage disequilibrium generally assume a simple model of recombination that only involves crossover events and ignore gene conversion. However, distinguishing the two processes is not only necessary for a complete description of recombination, but also essential for understanding the evolutionary consequences of inversions and other genomic partitions in which crossover (but not gene conversion) is reduced. We present heRho, a simple composite likelihood scheme for co-estimating the rate of crossover and gene conversion from individual diploid genomes. The method is based on analytic results for the distance-dependent probability of heterozygous and homozygous states at two loci. We apply heRho to simulations and data from the house mouse Mus musculus castaneus, a well studied model. Our analyses show i) that the rates of crossover and gene conversion can be accurately co-estimated at the level of individual chromosomes and ii) that previous estimates of the population scaled rate of recombination under a pure crossover model are likely biased.