Evolution of Alternative Adaptive Immune Systems in Vertebrates.

Adaptive immunity in jawless fishes is based on antigen recognition by three types of variable lymphocyte receptors (VLRs) composed of variable leucine-rich repeats, which are differentially expressed by two T-like lymphocyte lineages and one B-like lymphocyte lineage. The T-like cells express either VLRAs or VLRCs of yet undefined antigen specificity, whereas the VLRB antibodies secreted by B-like cells bind proteinaceous and carbohydrate antigens. The incomplete VLR germline genes are assembled into functional units by a gene conversion-like mechanism that employs flanking variable leucine-rich repeat sequences as templates in association with lineage-specific expression of cytidine deaminases. B-like cells develop in the hematopoietic typhlosole and kidneys, whereas T-like cells develop in the thymoid, a thymus-equivalent region at the gill fold tips. Thus, the dichotomy between T-like and B-like cells and the presence of dedicated lymphopoietic tissues emerge as ancestral vertebrate features, whereas the somatic diversification of structurally distinct antigen receptor genes evolved independently in jawless and jawed vertebrates.

Efficient embryonic homozygous gene conversion via RAD51-enhanced interhomolog repair.

Searching for factors to improve knockin efficiency for therapeutic applications, biotechnology, and generation of non-human primate models of disease, we found that the strand exchange protein RAD51 can significantly increase Cas9-mediated homozygous knockin in mouse embryos through an interhomolog repair (IHR) mechanism. IHR is a hallmark of meiosis but only occurs at low frequencies in somatic cells, and its occurrence in zygotes is controversial. Using multiple approaches, we provide evidence for an endogenous IHR mechanism in the early embryo that can be enhanced by RAD51. This process can be harnessed to generate homozygotes from wild-type zygotes using exogenous donors and to convert heterozygous alleles into homozygous alleles without exogenous templates. Furthermore, we identify additional IHR-promoting factors and describe features of IHR events. Together, our findings show conclusive evidence for IHR in mouse embryos and describe an efficient method for enhanced gene conversion.

Age-Dependent Alterations in Meiotic Recombination Cause Chromosome Segregation Errors in Spermatocytes.

Faithful chromosome segregation in meiosis requires crossover (CO) recombination, which is regulated to ensure at least one CO per homolog pair. We investigate the failure to ensure COs in juvenile male mice. By monitoring recombination genome-wide using cytological assays and at hotspots using molecular assays, we show that juvenile mouse spermatocytes have fewer COs relative to adults. Analysis of recombination in the absence of MLH3 provides evidence for greater utilization in juveniles of pathways involving structure-selective nucleases and alternative complexes, which can act upon precursors to generate noncrossovers (NCOs) at the expense of COs. We propose that some designated CO sites fail to mature efficiently in juveniles owing to inappropriate activity of these alternative repair pathways, leading to chromosome mis-segregation. We also find lower MutLγ focus density in juvenile human spermatocytes, suggesting that weaker CO maturation efficiency may explain why younger men have a higher risk of fathering children with Down syndrome.

Repeated strand invasion and extensive branch migration are hallmarks of meiotic recombination.

Currently favored models for meiotic recombination posit that both noncrossover and crossover recombination are initiated by DNA double-strand breaks but form by different mechanisms: noncrossovers by synthesis-dependent strand annealing and crossovers by formation and resolution of double Holliday junctions centered around the break. This dual mechanism hypothesis predicts different hybrid DNA patterns in noncrossover and crossover recombinants. We show that these predictions are not upheld, by mapping with unprecedented resolution parental strand contributions to recombinants at a model locus. Instead, break repair in both noncrossovers and crossovers involves synthesis-dependent strand annealing, often with multiple rounds of strand invasion. Crossover-specific double Holliday junction formation occurs via processes involving branch migration as an integral feature, one that can be separated from repair of the break itself. These findings reveal meiotic recombination to be a highly dynamic process and prompt a new view of the relationship between crossover and noncrossover recombination.

Identification of a robust DNA methylation signature for Fanconi anemia.

Fanconi anemia (FA) is a clinically variable and genetically heterogeneous cancer-predisposing disorder representing the most common bone marrow failure syndrome. It is caused by inactivating predominantly biallelic mutations involving >20 genes encoding proteins with roles in the FA/BRCA DNA repair pathway. Molecular diagnosis of FA is challenging due to the wide spectrum of the contributing gene mutations and structural rearrangements. The assessment of chromosomal fragility after exposure to DNA cross-linking agents is generally required to definitively confirm diagnosis. We assessed peripheral blood genome-wide DNA methylation (DNAm) profiles in 25 subjects with molecularly confirmed clinical diagnosis of FA (FANCA complementation group) using Illumina's Infinium EPIC array. We identified 82 differentially methylated CpG sites that allow to distinguish subjects with FA from healthy individuals and subjects with other genetic disorders, defining an FA-specific DNAm signature. The episignature was validated using a second cohort of subjects with FA involving different complementation groups, documenting broader genetic sensitivity and demonstrating its specificity using the EpiSign Knowledge Database. The episignature properly classified DNA samples obtained from bone marrow aspirates, demonstrating robustness. Using the selected probes, we trained a machine-learning model able to classify EPIC DNAm profiles in molecularly unsolved cases. Finally, we show that the generated episignature includes CpG sites that do not undergo functional selective pressure, allowing diagnosis of FA in individuals with reverted phenotype due to gene conversion. These findings provide a tool to accelerate diagnostic testing in FA and broaden the clinical utility of DNAm profiling in the diagnostic setting.

Mechanistic View and Genetic Control of DNA Recombination during Meiosis.

Meiotic recombination is essential for fertility and allelic shuffling. Canonical recombination models fail to capture the observed complexity of meiotic recombinants. Here, by combining genome-wide meiotic heteroduplex DNA patterns with meiotic DNA double-strand break (DSB) sites, we show that part of this complexity results from frequent template switching during synthesis-dependent strand annealing that yields noncrossovers and from branch migration of double Holliday junction (dHJ)-containing intermediates that mainly yield crossovers. This complexity also results from asymmetric positioning of crossover intermediates relative to the initiating DSB and Msh2-independent conversions promoted by the suspected dHJ resolvase Mlh1-3 as well as Exo1 and Sgs1. Finally, we show that dHJ resolution is biased toward cleavage of the pair of strands containing newly synthesized DNA near the junctions and that this bias can be decoupled from the crossover-biased dHJ resolution. These properties are likely conserved in eukaryotes containing ZMM proteins, which includes mammals.

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.

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.

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.

Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-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 a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing 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 polymorphism-aware phylogenetic models' 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 color 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.

PIF1 helicase promotes break-induced replication in mammalian cells.

Break-induced replication (BIR) is a specialized homologous-recombination pathway for DNA double-strand break (DSB) repair, which often induces genome instability. In this study, we establish EGFP-based recombination reporters to systematically study BIR in mammalian cells and demonstrate an important role of human PIF1 helicase in promoting BIR. We show that at endonuclease cleavage sites, PIF1-dependent BIR is used for homology-initiated recombination requiring long track DNA synthesis, but not short track gene conversion (STGC). We also show that structure formation-prone AT-rich DNA sequences derived from common fragile sites (CFS-ATs) induce BIR upon replication stress and oncogenic stress, and PCNA-dependent loading of PIF1 onto collapsed/broken forks is critical for BIR activation. At broken replication forks, even STGC-mediated repair of double-ended DSBs depends on POLD3 and PIF1, revealing an unexpected mechanism of BIR activation upon replication stress that differs from the conventional BIR activation model requiring DSB end sensing at endonuclease-generated breaks. Furthermore, loss of PIF1 is synthetically lethal with loss of FANCM, which is involved in protecting CFS-ATs. The breast cancer-associated PIF1 mutant L319P is defective in BIR, suggesting a direct link of BIR to oncogenic processes.

Regulation of hetDNA Length during Mitotic Double-Strand Break Repair in Yeast.

Heteroduplex DNA (hetDNA) is a key molecular intermediate during the repair of mitotic double-strand breaks by homologous recombination, but its relationship to 5' end resection and/or 3' end extension is poorly understood. In the current study, we examined how perturbations in these processes affect the hetDNA profile associated with repair of a defined double-strand break (DSB) by the synthesis-dependent strand-annealing (SDSA) pathway. Loss of either the Exo1 or Sgs1 long-range resection pathway significantly shortened hetDNA, suggesting that these pathways normally collaborate during DSB repair. In addition, altering the processivity or proofreading activity of DNA polymerase δ shortened hetDNA length or reduced break-adjacent mismatch removal, respectively, demonstrating that this is the primary polymerase that extends both 3' ends. Data are most consistent with the extent of DNA synthesis from the invading end being the primary determinant of hetDNA length during SDSA.

Ancient variation of the AvrPm17 gene in powdery mildew limits the effectiveness of the introgressed rye Pm17 resistance gene in wheat.

Introgressions of chromosomal segments from related species into wheat are important sources of resistance against fungal diseases. The durability and effectiveness of introgressed resistance genes upon agricultural deployment is highly variable-a phenomenon that remains poorly understood, as the corresponding fungal avirulence genes are largely unknown. Until its breakdown, the Pm17 resistance gene introgressed from rye to wheat provided broad resistance against powdery mildew (Blumeria graminis). Here, we used quantitative trait locus (QTL) mapping to identify the corresponding wheat mildew avirulence effector AvrPm17. It is encoded by two paralogous genes that exhibit signatures of reoccurring gene conversion events and are members of a mildew sublineage specific effector cluster. Extensive haplovariant mining in wheat mildew and related sublineages identified several ancient virulent AvrPm17 variants that were present as standing genetic variation in wheat powdery mildew prior to the Pm17 introgression, thereby paving the way for the rapid breakdown of the Pm17 resistance. QTL mapping in mildew identified a second genetic component likely corresponding to an additional resistance gene present on the 1AL.1RS translocation carrying Pm17. This gene remained previously undetected due to suppressed recombination within the introgressed rye chromosomal segment. We conclude that the initial effectiveness of 1AL.1RS was based on simultaneous introgression of two genetically linked resistance genes. Our results demonstrate the relevance of pathogen-based genetic approaches to disentangling complex resistance loci in wheat. We propose that identification and monitoring of avirulence gene diversity in pathogen populations become an integral part of introgression breeding to ensure effective and durable resistance in wheat.

The landscape of submicroscopic structural variants at the OPN1LW/OPN1MW gene cluster on Xq28 underlying blue cone monochromacy.

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.

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.

Interlocus Gene Conversion, Natural Selection, and Paralog Homogenization.

Following a duplication, the resulting paralogs tend to diverge. While mutation and natural selection can accelerate this process, they can also slow it. Here, we quantify the paralog homogenization that is caused by point mutations and interlocus gene conversion (IGC). Among 164 duplicated teleost genes, the median percentage of postduplication codon substitutions that arise from IGC rather than point mutation is estimated to be between 7% and 8%. By differentiating between the nonsynonymous codon substitutions that homogenize the protein sequences of paralogs and the nonhomogenizing nonsynonymous substitutions, we estimate the homogenizing nonsynonymous rates to be higher for 163 of the 164 teleost data sets as well as for all 14 data sets of duplicated yeast ribosomal protein-coding genes that we consider. For all 14 yeast data sets, the estimated homogenizing nonsynonymous rates exceed the synonymous rates.

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.

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.

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.