Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history.

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future.

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

Polymorphism-aware models in RevBayes: Species trees, disentangling Balancing Selection and CG-biased gene conversion.

The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the polymorphism-aware phylogenetic models (PoMos) framework rooted in the Moran model, we introduce PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analysing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in PoMos' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited colour dimorphism in Drosophila erecta.

Increased Positive Selection in Highly Recombining Genes Does not Necessarily Reflect an Evolutionary Advantage of Recombination.

It is commonly thought that the long-term advantage of meiotic recombination is to dissipate genetic linkage, allowing natural selection to act independently on different loci. It is thus theoretically expected that genes with higher recombination rates evolve under more effective selection. On the other hand, recombination is often associated with GC-biased gene conversion (gBGC), which theoretically interferes with selection by promoting the fixation of deleterious GC alleles. To test these predictions, several studies assessed whether selection was more effective in highly recombining genes (due to dissipation of genetic linkage) or less effective (due to gBGC), assuming a fixed distribution of fitness effects (DFE) for all genes. In this study, I directly derive the DFE from a gene's evolutionary history (shaped by mutation, selection, drift, and gBGC) under empirical fitness landscapes. I show that genes that have experienced high levels of gBGC are less fit and thus have more opportunities for beneficial mutations. Only a small decrease in the genome-wide intensity of gBGC leads to the fixation of these beneficial mutations, particularly in highly recombining genes. This results in increased positive selection in highly recombining genes that is not caused by more effective selection. Additionally, I show that the death of a recombination hotspot can lead to a higher dN/dS than its birth, but with substitution patterns biased towards AT, and only at selected positions. This shows that controlling for a substitution bias towards GC is therefore not sufficient to rule out the contribution of gBGC to signatures of accelerated evolution. Finally, although gBGC does not affect the fixation probability of GC-conservative mutations, I show that by altering the DFE, gBGC can also significantly affect nonsynonymous GC-conservative substitution patterns.

High-throughput classification of S. cerevisiae tetrads using deep learning.

Meiotic crossovers play a vital role in proper chromosome segregation and evolution of most sexually reproducing organisms. Meiotic recombination can be visually observed in Saccharomyces cerevisiae tetrads using linked spore-autonomous fluorescent markers placed at defined intervals within the genome, which allows for analysis of meiotic segregation without the need for tetrad dissection. To automate the analysis, we developed a deep learning-based image recognition and classification pipeline for high-throughput tetrad detection and meiotic crossover classification. As a proof of concept, we analyzed a large image data set from wild-type and selected gene knock-out mutants to quantify crossover frequency, interference, chromosome missegregation, and gene conversion events. The deep learning-based method has the potential to accelerate the discovery of new genes involved in meiotic recombination in S. cerevisiae such as the underlying factors controlling crossover frequency and interference.

Spontaneous and double-strand break repair-associated quasipalindrome and frameshift mutagenesis in budding yeast: role of mismatch repair.

Although gene conversion (GC) in Saccharomyces cerevisiae is the most error-free way to repair double-strand breaks (DSBs), the mutation rate during homologous recombination is 1,000 times greater than during replication. Many mutations involve dissociating a partially copied strand from its repair template and re-aligning with the same or another template, leading to -1 frameshifts in homonucleotide runs, quasipalindrome (QP)-associated mutations and microhomology-mediated interchromosomal template switches. We studied GC induced by HO endonuclease cleavage at MATα, repaired by an HMR::KI-URA3 donor. We inserted into HMR::KI-URA3 an 18-bp inverted repeat where one arm had a 4-bp insertion. Most GCs yield MAT::KI-ura3::QP + 4 (Ura-) outcomes, but template-switching produces Ura+ colonies, losing the 4-bp insertion. If the QP arm without the insertion is first encountered by repair DNA polymerase and is then (mis)used as a template, the palindrome is perfected. When the QP + 4 arm is encountered first, Ura+ derivatives only occur after second-end capture and second-strand synthesis. QP + 4 mutations are suppressed by mismatch repair (MMR) proteins Msh2, Msh3, and Mlh1, but not Msh6. Deleting Rdh54 significantly reduces QP mutations only when events creating Ura+ occur in the context of a D-loop but not during second-strand synthesis. A similar bias is found with a proofreading-defective DNA polymerase mutation (poI3-01). DSB-induced mutations differed in several genetic requirements from spontaneous events. We also created a + 1 frameshift in the donor, expanding a run of 4 Cs to 5 Cs. Again, Ura+ recombinants markedly increased by disabling MMR, suggesting that MMR acts during GC but favors the unbroken, template strand.

On the estimation of genome-average recombination rates.

The rate at which recombination events occur in a population is an indicator of its effective population size and the organism's reproduction mode. It determines the extent of linkage disequilibrium along the genome and, thereby, the efficacy of both purifying and positive selection. The population recombination rate can be inferred using models of genome evolution in populations. Classic methods based on the patterns of linkage disequilibrium provide the most accurate estimates, providing large sample sizes are used and the demography of the population is properly accounted for. Here, the capacity of approaches based on the sequentially Markov coalescent (SMC) to infer the genome-average recombination rate from as little as a single diploid genome is examined. SMC approaches provide highly accurate estimates even in the presence of changing population sizes, providing that (1) within genome heterogeneity is accounted for and (2) classic maximum-likelihood optimization algorithms are employed to fit the model. SMC-based estimates proved sensitive to gene conversion, leading to an overestimation of the recombination rate if conversion events are frequent. Conversely, methods based on the correlation of heterozygosity succeed in disentangling the rate of crossing over from that of gene conversion events, but only when the population size is constant and the recombination landscape homogeneous. These results call for a convergence of these two methods to obtain accurate and comparable estimates of recombination rates between populations.

Enigmatic Nodal and Lefty gene repertoire discrepancy: Latent evolutionary history revealed by vertebrate-wide phylogeny.

Homology in vertebrate body plans is traditionally ascribed to the high-level conservation of regulatory components within the genetic programs governing them, particularly during the "phylotypic stage." However, advancements in embryology and molecular phylogeny have unveiled the dynamic nature of gene repertoires responsible for early development. Notably, the Nodal and Lefty genes, members of the transforming growth factor-beta superfamily producing intercellular signaling molecules and crucial for left-right (L-R) symmetry breaking, exhibit distinctive features within their gene repertoires. These features encompass among-species gene repertoire variations resulting from gene gain and loss, as well as gene conversion. Despite their significance, these features have been largely unexplored in a phylogenetic context, but accumulating genome-wide sequence information is allowing the scrutiny of these features. It has exposed hidden paralogy between Nodal1 and Nodal2 genes resulting from differential gene loss in amniotes. In parallel, the tandem cluster of Lefty1 and Lefty2 genes, which was thought to be confined to mammals, is observed in sharks and rays, with an unexpected phylogenetic pattern. This article provides a comprehensive review of the current understanding of the origins of these vertebrate gene repertoires and proposes a revised nomenclature based on the elucidated history of vertebrate genome evolution.

The Patterns of Codon Usage between Chordates and Arthropods are Different but Co-evolving with Mutational Biases.

Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic.

GC Content Across Insect Genomes: Phylogenetic Patterns, Causes and Consequences.

The proportions of A:T and G:C nucleotide pairs are often unequal and can vary greatly between animal species and along chromosomes. The causes and consequences of this variation are incompletely understood. The recent release of high-quality genome sequences from the Darwin Tree of Life and other large-scale genome projects provides an opportunity for GC heterogeneity to be compared across a large number of insect species. Here we analyse GC content along chromosomes, and within protein-coding genes and codons, of 150 insect species from four holometabolous orders: Coleoptera, Diptera, Hymenoptera, and Lepidoptera. We find that protein-coding sequences have higher GC content than the genome average, and that Lepidoptera generally have higher GC content than the other three insect orders examined. GC content is higher in small chromosomes in most Lepidoptera species, but this pattern is less consistent in other orders. GC content also increases towards subtelomeric regions within protein-coding genes in Diptera, Coleoptera and Lepidoptera. Two species of Diptera, Bombylius major and B. discolor, have very atypical genomes with ubiquitous increase in AT content, especially at third codon positions. Despite dramatic AT-biased codon usage, we find no evidence that this has driven divergent protein evolution. We argue that the GC landscape of Lepidoptera, Diptera and Coleoptera genomes is influenced by GC-biased gene conversion, strongest in Lepidoptera, with some outlier taxa affected drastically by counteracting processes.

Prevalent role of homologous recombination in the repair of specific double-strand breaks in Rhizobium etli.

Double-strand breaks (DSBs) are the most dangerous injuries for a genome. When unrepaired, death quickly ensues. In most bacterial systems, DSBs are repaired through homologous recombination. Nearly one-quarter of bacterial species harbor a second system, allowing direct ligation of broken ends, known as Non-Homologous End Joining (NHEJ). The relative role of both systems in DSBs repair in bacteria has been explored only in a few cases. To evaluate this in the bacterium Rhizobium etli, we used a modified version of the symbiotic plasmid (264 kb), containing a single copy of the nifH gene. In this plasmid, we inserted an integrative plasmid harboring a modified nifH gene fragment containing an I-SceI site. DSBs were easily inflicted in vivo by conjugating a small, replicative plasmid that expresses the I-SceI nuclease into the appropriate strains. Repair of a DSB may be achieved through homologous recombination (either between adjacent or distant repeats) or NHEJ. Characterization of the derivatives that repaired DSB in different configurations, revealed that in most cases (74%), homologous recombination was the prevalent mechanism responsible for repair, with a relatively minor contribution of NHEJ (23%). Inactivation of the I-SceI gene was detected in 3% of the cases. Sequence analysis of repaired derivatives showed the operation of NHEJ. To enhance the number of derivatives repaired through NHEJ, we repeated these experiments in a recA mutant background. Derivatives showing NHEJ were readily obtained when the DSB occurred on a small, artificial plasmid in a recA mutant. However, attempts to deliver a DSB on the symbiotic plasmid in a recA background failed, due to the accumulation of mutations that inactivated the I-SceI gene. This result, coupled with the absence of derivatives that lost the nonessential symbiotic plasmid, may be due to an unusual stability of the symbiotic plasmid, possibly caused by the presence of multiple toxin-antitoxin modules.

Half a Century of Controversy: The Neutralist/Selectionist Debate in Molecular Evolution.

The neutral and nearly neutral theories, introduced more than 50 yr ago, have raised and still raise passionate discussion regarding the forces governing molecular evolution and their relative importance. The debate, initially focused on the amount of within-species polymorphism and constancy of the substitution rate, has spread, matured, and now underlies a wide range of topics and questions. The neutralist/selectionist controversy has structured the field and influences the way molecular evolutionary scientists conceive their research.

Molecular analysis of the role of polymerase theta in gene targeting in Arabidopsis thaliana.

Precise genetic modification can be achieved via a sequence homology-mediated process known as gene targeting (GT). Whilst established for genome engineering purposes, the application of GT in plants still suffers from a low efficiency for which an explanation is currently lacking. Recently reported reduced rates of GT in A. thaliana deficient in polymerase theta (Polθ), a core component of theta-mediated end joining (TMEJ) of DNA breaks, have led to the suggestion of a direct involvement of this enzyme in the homology-directed process. Here, by monitoring homology-driven gene conversion in plants with CRISPR reagent and donor sequences pre-integrated at random sites in the genome (in planta GT), we demonstrate that Polθ action is not required for GT, but instead suppresses the process, likely by promoting the repair of the DNA break by end-joining. This finding indicates that lack of donor integration explains the previously established reduced GT rates seen upon transformation of Polθ-deficient plants. Our study additionally provides insight into ectopic gene targeting (EGT), recombination events between donor and target that do not map to the target locus. EGT, which occurs at similar frequencies as "true" GT during transformation, was rare in our in planta GT experiments arguing that EGT predominantly results from target locus recombination with nonintegrated T-DNA molecules. By describing mechanistic features of GT our study provides directions for the improvement of precise genetic modification of plants.

The influence of GC-biased gene conversion on non-adaptive sequence evolution in short introns of Drosophila melanogaster.

Population genetic inference of selection on the nucleotide sequence level often proceeds by comparison to a reference sequence evolving only under mutation and population demography. Among the few candidates for such a reference sequence is the 5' part of short introns (5SI) in Drosophila. In addition to mutation and population demography, however, there is evidence for a weak force favouring GC bases, likely due to GC-biased gene conversion (gBGC), and for the effect of linked selection. Here, we use polymorphism and divergence data of Drosophila melanogaster to detect and describe the forces affecting the evolution of the 5SI. We separately analyse mutation classes, compare them between chromosomes, and relate them to recombination rate frequencies. GC-conservative mutations seem to be mainly influenced by mutation and drift, with linked selection mostly causing differences between the central and the peripheral (i.e., telomeric and centromeric) regions of the chromosome arms. Comparing GC-conservative mutation patterns between autosomes and the X chromosome showed differences in mutation rates, rather than linked selection, in the central chromosomal regions after accounting for differences in effective population sizes. On the other hand, GC-changing mutations show asymmetric site frequency spectra, indicating the presence of gBGC, varying among mutation classes and in intensity along chromosomes, but approximately equal in strength in autosomes and the X chromosome.

Shared evolutionary processes shape landscapes of genomic variation in the great apes.

For at least the past 5 decades, population genetics, as a field, has worked to describe the precise balance of forces that shape patterns of variation in genomes. The problem is challenging because modeling the interactions between evolutionary processes is difficult, and different processes can impact genetic variation in similar ways. In this paper, we describe how diversity and divergence between closely related species change with time, using correlations between landscapes of genetic variation as a tool to understand the interplay between evolutionary processes. We find strong correlations between landscapes of diversity and divergence in a well-sampled set of great ape genomes, and explore how various processes such as incomplete lineage sorting, mutation rate variation, GC-biased gene conversion and selection contribute to these correlations. Through highly realistic, chromosome-scale, forward-in-time simulations, we show that the landscapes of diversity and divergence in the great apes are too well correlated to be explained via strictly neutral processes alone. Our best fitting simulation includes both deleterious and beneficial mutations in functional portions of the genome, in which 9% of fixations within those regions is driven by positive selection. This study provides a framework for modeling genetic variation in closely related species, an approach which can shed light on the complex balance of forces that have shaped genetic variation.

Novel Insights into the Landscape of Crossover and Noncrossover Events in Rhesus Macaques (Macaca mulatta).

Meiotic recombination landscapes differ greatly between distantly and closely related taxa, populations, individuals, sexes, and even within genomes; however, the factors driving this variation are yet to be well elucidated. Here, we directly estimate contemporary crossover rates and, for the first time, noncrossover rates in rhesus macaques (Macaca mulatta) from four three-generation pedigrees comprising 32 individuals. We further compare these results with historical, demography-aware, linkage disequilibrium-based recombination rate estimates. From paternal meioses in the pedigrees, 165 crossover events with a median resolution of 22.3 kb were observed, corresponding to a male autosomal map length of 2,357 cM-approximately 15% longer than an existing linkage map based on human microsatellite loci. In addition, 85 noncrossover events with a mean tract length of 155 bp were identified-similar to the tract lengths observed in the only other two primates in which noncrossovers have been studied to date, humans and baboons. Consistent with observations in other placental mammals with PRDM9-directed recombination, crossover (and to a lesser extent noncrossover) events in rhesus macaques clustered in intergenic regions and toward the chromosomal ends in males-a pattern in broad agreement with the historical, sex-averaged recombination rate estimates-and evidence of GC-biased gene conversion was observed at noncrossover sites.

No Evidence that Selection on Synonymous Codon Usage Affects Patterns of Protein Evolution in Bacteria.

Bias in synonymous codon usage has been reported across all kingdoms of life. Evidence suggests that codon usage bias is often driven by selective pressures, typically for translational efficiency. These selective pressures have been shown to depress the rate at which synonymous sites evolve. We hypothesize that selection on synonymous codon use could also slow the rate of protein evolution if a non-synonymous mutation changes the codon from being preferred to unpreferred. We test this hypothesis by looking at patterns of protein evolution using polymorphism and substitution data in two bacterial species, Escherichia coli and Streptococcus pneumoniae. We find no evidence that non-synonymous mutations that change a codon from being unpreferred to preferred are more common than the opposite. Overall, selection on codon bias seems to have little influence over non-synonymous polymorphism or substitution patterns.

Haplotype-based phylogenetic analysis and population genomics uncover the origin and domestication of sweetpotato.

The hexaploid sweetpotato (Ipomoea batatas) is one of the most important root crops worldwide. However, its genetic origin remains controversial, and its domestication history remains unknown. In this study, we used a range of genetic evidence and a newly developed haplotype-based phylogenetic analysis to identify two probable progenitors of sweetpotato. The diploid progenitor was likely closely related to Ipomoea aequatoriensis and contributed the B1 subgenome, IbT-DNA2, and the lineage 1 type of chloroplast genome to sweetpotato. The tetraploid progenitor of sweetpotato was most likely I. batatas 4x, which donated the B2 subgenome, IbT-DNA1, and the lineage 2 type of chloroplast genome. Sweetpotato most likely originated from reciprocal crosses between the diploid and tetraploid progenitors, followed by a subsequent whole-genome duplication. In addition, we detected biased gene exchanges between the subgenomes; the rate of B1 to B2 subgenome conversions was nearly three times higher than that of B2 to B1 subgenome conversions. Our analyses revealed that genes involved in storage root formation, maintenance of genome stability, biotic resistance, sugar transport, and potassium uptake were selected during the speciation and domestication of sweetpotato. This study sheds light on the evolution of sweetpotato and paves the way for improvement of this crop.

GC heterogeneity reveals sequence-structures evolution of angiosperm ITS2.

BACKGROUND: Despite GC variation constitutes a fundamental element of genome and species diversity, the precise mechanisms driving it remain unclear. The abundant sequence data available for the ITS2, a commonly employed phylogenetic marker in plants, offers an exceptional resource for exploring the GC variation across angiosperms. RESULTS: A comprehensive selection of 8666 species, comprising 165 genera, 63 families, and 30 orders were used for the analyses. The alignment of ITS2 sequence-structures and partitioning of secondary structures into paired and unpaired regions were performed using 4SALE. Substitution rates and frequencies among GC base-pairs in the paired regions of ITS2 were calculated using RNA-specific models in the PHASE package. The results showed that the distribution of ITS2 GC contents on the angiosperm phylogeny was heterogeneous, but their increase was generally associated with ITS2 sequence homogenization, thereby supporting the occurrence of GC-biased gene conversion (gBGC) during the concerted evolution of ITS2. Additionally, the GC content in the paired regions of the ITS2 secondary structure was significantly higher than that of the unpaired regions, indicating the selection of GC for thermodynamic stability. Furthermore, the RNA substitution models demonstrated that base-pair transformations favored both the elevation and fixation of GC in the paired regions, providing further support for gBGC. CONCLUSIONS: Our findings highlight the significance of secondary structure in GC investigation, which demonstrate that both gBGC and structure-based selection are influential factors driving angiosperm ITS2 GC content.

Phenotype reversion as "natural gene therapy" in Fanconi anemia by a gene conversion event.

Somatic mosaicism appears as a recurrent phenomenon among patients suffering from Fanconi anemia (FA), but its direct prognostic significance mostly remains an open question. The clinical picture of FA mosaic subjects could indeed vary from just mild features to severe hematologic failure. Here, we illustrate the case of a proband whose FA familiarity, modest signs (absence of hematological anomalies and fertility issues), and chromosome fragility test transition to negative overtime were suggestive of somatic mosaicism. In line with this hypothesis, genetic testing on patient's peripheral blood and buccal swab reported the presence of the only FANCA paternal variant (FANCA:c.2638C>T, p. Arg880*) and of both parental alleles (the additional FANCA:c.3164G>A, p. Arg1055Gln), respectively. Moreover, the SNP analysis performed on the same biological specimens allowed us to attribute the proband's mosaicism status to a possible gene conversion mechanism. Our case clearly depicts the positive association between somatic mosaicism and the proband's favorable clinical course due to the occurrence of the reversion event at the hematopoietic stem cell level. Since this condition concerns only a limited subgroup of FA individuals, the accurate evaluation of the origin and extent of clonality would be key to steer clinicians toward the most appropriate therapeutic decision for their FA mosaic patients.