Gene conversion generates evolutionary novelty that fuels genetic conflicts

WHO elements - A new category of selfish genetic elements at the borderline between homing elements and transposable elements

Homing genetic elements are a form of selfish DNA that inserts into a specific target site in the genome and spreads through the population by a process of biased inheritance. Two well-known types of homing element, called inteins and homing introns, were discovered decades ago. In this review we describe WHO elements, a newly discovered type of homing element that constitutes a distinct third category but is rare, having been found only in a few yeast species so far. WHO elements are inferred to spread using the same molecular homing mechanism as inteins and introns: they encode a site-specific endonuclease that cleaves the genome at the target site, making a DNA break that is subsequently repaired by copying the element. For most WHO elements, the target site is in the glycolytic gene FBA1. WHO elements differ from inteins and homing introns in two fundamental ways: they do not interrupt their host gene (FBA1), and they occur in clusters. The clusters were formed by successive integrations of different WHO elements into the FBA1 locus, the result of an 'arms race' between the endonuclease and its target site. We also describe one family of WHO elements (WHO10) that is no longer specifically associated with the FBA1 locus and instead appears to have become transposable, inserting at random genomic sites in Torulaspora globosa with up to 26 copies per strain. The WHO family of elements is therefore at the borderline between homing genetic elements and transposable elements.

Evolutionary dynamics of pro-inflammatory caspases in primates and rodents

Caspase-1 and related proteases are key players in inflammation and innate immunity. Here, we characterize the evolutionary history of caspase-1 and its close relatives across 19 primates and 21 rodents, focusing on differences that may cause discrepancies between humans and animal studies. While caspase-1 has been retained in all these taxa, other members of the caspase-1 subfamily (caspase-4, -5, -11, -12, and CARD16, 17, and 18) each have unique evolutionary trajectories. Caspase-4 is found across simian primates, whereas we identified multiple pseudogenization and gene loss events in caspase-5, caspase-11, and the CARDs. Because caspases-4 and -11 are both key players in the non-canonical inflammasome pathway, we expected that these proteins would be likely to evolve rapidly. Instead, we found that these two proteins are largely conserved, whereas caspase-4's close paralog, caspase-5, showed significant indications of positive selection, as did primate caspase-1. Caspase-12 is a non-functional pseudogene in humans. We find this extends across most primates, although many rodents and some primates retain an intact, and likely functional, caspase-12. In mouse laboratory lines, we found that 50% of common strains carry non-synonymous variants that may impact the functions of caspase-11 and -12, and therefore recommend specific strains to be used (and avoided). Finally, unlike rodents, primate caspases have undergone repeated rounds of gene conversion, duplication, and loss leading to a highly dynamic pro-inflammatory caspase repertoire. Thus we uncovered many differences in the evolution of primate and rodent pro-inflammatory caspases, and discuss the potential implications of this history for caspase gene functions.

Intermolecular Gene Conversion for the Equalization of Genome Copies in the Polyploid Haloarchaeon Haloferax volcanii: Identification of Important Proteins

The model haloarchaeon Haloferax volcanii is polyploid with about 20 copies of its major chromosome. Recently it has been described that highly efficient intermolecular gene conversion operates in H. volcanii to equalize the chromosomal copies. In the current study, 24 genes were selected that encode proteins with orthologs involved in gene conversion or homologous recombination in archaea, bacteria, or eukaryotes. Single gene deletion strains of 22 genes and a control gene were constructed in two parent strains for a gene conversion assay; only radA and radB were shown to be essential. Protoplast fusions were used to generate strains that were heterozygous for the gene HVO_2528, encoding an enzyme for carotinoid biosynthesis. It was revealed that a lack of six of the proteins did not influence the efficiency of gene conversion, while sixteen mutants had severe gene conversion defects. Notably, lack of paralogous proteins of gene families had very different effects, e.g., mutant Δrad25b had no phenotype, while mutants Δrad25a, Δrad25c, and Δrad25d were highly compromised. Generation of a quadruple rad25 and a triple sph deletion strain also indicated that the paralogs have different functions, in contrast to sph2 and sph4, which cannot be deleted simultaneously. There was no correlation between the severity of the phenotypes and the respective transcript levels under non-stressed conditions, indicating that gene expression has to be induced at the onset of gene conversion. Phylogenetic trees of the protein families Rad3/25, MutL/S, and Sph/SMC/Rad50 were generated to unravel the history of the paralogous proteins of H. volcanii. Taken together, unselected intermolecular gene conversion in H. volcanii involves at least 16 different proteins, the molecular roles of which can be studied in detail in future projects.

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.

Repeat-Induced Point Mutation and Gene Conversion Coinciding with Heterochromatin Shape the Genome of a Plant-Pathogenic Fungus

Meiosis is associated with genetic changes in the genome-via recombination, gene conversion, and mutations. The occurrence of gene conversion and mutations during meiosis may further be influenced by the chromatin conformation, similar to the effect of the chromatin conformation on the mitotic mutation rate. To date, however, the exact distribution and type of meiosis-associated changes and the role of the chromatin conformation in this context are largely unexplored. Here, we determine recombination, gene conversion, and de novo mutations using whole-genome sequencing of all meiotic products of 23 individual meioses in Zymoseptoria tritici, an important pathogen of wheat. We confirm a high genome-wide recombination rate of 65 centimorgan (cM)/Mb and see higher recombination rates on the accessory compared to core chromosomes. A substantial fraction of 0.16% of all polymorphic markers was affected by gene conversions, showing a weak GC-bias and occurring at higher frequency in regions of constitutive heterochromatin, indicated by the histone modification H3K9me3. The de novo mutation rate associated with meiosis was approximately three orders of magnitude higher than the corresponding mitotic mutation rate. Importantly, repeat-induced point mutation (RIP), a fungal defense mechanism against duplicated sequences, is active in Z. tritici and responsible for the majority of these de novo meiotic mutations. Our results indicate that the genetic changes associated with meiosis are a major source of variability in the genome of an important plant pathogen and shape its evolutionary trajectory. IMPORTANCE The impact of meiosis on the genome composition via gene conversion and mutations is mostly poorly understood, in particular, for non-model species. Here, we sequenced all four meiotic products for 23 individual meioses and determined the genetic changes caused by meiosis for the important fungal wheat pathogen Zymoseptoria tritici. We found a high rate of gene conversions and an effect of the chromatin conformation on gene conversion rates. Higher conversion rates were found in regions enriched with the H3K9me3-a mark for constitutive heterochromatin. Most importantly, meiosis was associated with a much higher frequency of de novo mutations than mitosis; 78% of the meiotic mutations were caused by repeat-induced point mutations-a fungal defense mechanism against duplicated sequences. In conclusion, the genetic changes associated with meiosis are therefore a major factor shaping the genome of this fungal pathogen.

Egg Yolk Protein Homologs Identified in Live-Bearing Sharks: Co-Opted in the Lecithotrophy-to-Matrotrophy Shift?

Reproductive modes of vertebrates are classified into two major embryonic nutritional types: yolk deposits (i.e., lecithotrophy) and maternal investment (i.e., matrotrophy). Vitellogenin (VTG), a major egg yolk protein synthesized in the female liver, is one of the molecules relevant to the lecithotrophy-to-matrotrophy shift in bony vertebrates. In mammals, all VTG genes are lost following the lecithotrophy-to-matrotrophy shift, and it remains to be elucidated whether the lecithotrophy-to-matrotrophy shift in nonmammalians is also associated with VTG repertoire modification. In this study, we focused on chondrichthyans (cartilaginous fishes)-a vertebrate clade that underwent multiple lecithotrophy-to-matrotrophy shifts. For an exhaustive search of homologs, we performed tissue-by-tissue transcriptome sequencing for two viviparous chondrichthyans, the frilled shark Chlamydoselachus anguineus and the spotless smooth-hound Mustelus griseus, and inferred the molecular phylogeny of VTG and its receptor very low-density lipoprotein receptor (VLDLR), across diverse vertebrates. As a result, we identified either three or four VTG orthologs in chondrichthyans including viviparous species. We also showed that chondrichthyans had two additional VLDLR orthologs previously unrecognized in their unique lineage (designated as VLDLRc2 and VLDLRc3). Notably, VTG gene expression patterns differed in the species studied depending on their reproductive mode; VTGs are broadly expressed in multiple tissues, including the uterus, in the two viviparous sharks, and in addition to the liver. This finding suggests that the chondrichthyans VTGs do not only function as the yolk nutrient but also as the matrotrophic factor. Altogether, our study indicates that the lecithotrophy-to-matrotrophy shift in chondrichthyans was achieved through a distinct evolutionary process from mammals.

The wtf meiotic driver gene family has unexpectedly persisted for over 100 million years

Meiotic drivers are selfish elements that bias their own transmission into more than half of the viable progeny produced by a driver+/driver- heterozygote. Meiotic drivers are thought to exist for relatively short evolutionary timespans because a driver gene or gene family is often found in a single species or in a group of very closely related species. Additionally, drivers are generally considered doomed to extinction when they spread to fixation or when suppressors arise. In this study, we examine the evolutionary history of the wtf meiotic drivers first discovered in the fission yeast Schizosaccharomyces pombe. We identify homologous genes in three other fission yeast species, S. octosporus, S. osmophilus, and S. cryophilus, which are estimated to have diverged over 100 million years ago from the S. pombe lineage. Synteny evidence supports that wtf genes were present in the common ancestor of these four species. Moreover, the ancestral genes were likely drivers as wtf genes in S. octosporus cause meiotic drive. Our findings indicate that meiotic drive systems can be maintained for long evolutionary timespans.

Evolution of host-microbe cell adherence by receptor domain shuffling

Stable adherence to epithelial surfaces is required for colonization by diverse host-associated microbes. Successful attachment of pathogenic microbes to host cells via adhesin molecules is also the first step in many devastating infections. Despite the primacy of epithelial adherence in establishing host-microbe associations, the evolutionary processes that shape this crucial interface remain enigmatic. Carcinoembryonic antigen-related cell adhesion molecules (CEACAMs) encompass a multifunctional family of vertebrate cell surface proteins which are recurrent targets of bacterial adhesins at epithelial barriers. Here, we show that multiple members of the primate CEACAM family exhibit evidence of repeated natural selection at protein surfaces targeted by bacteria, consistent with pathogen-driven evolution. Divergence of CEACAM proteins between even closely related great apes is sufficient to control molecular interactions with a range of bacterial adhesins. Phylogenetic analyses further reveal that repeated gene conversion of CEACAM extracellular domains during primate divergence plays a key role in limiting bacterial adhesin host tropism. Moreover, we demonstrate that gene conversion has continued to shape CEACAM diversity within human populations, with abundant human CEACAM1 variants mediating evasion of adhesins from pathogenic Neisseria. Together this work reveals a mechanism by which gene conversion shapes first contact between microbes and animal hosts.

Trillions of bacteria live in and on the human body. Most of them are harmless but some can cause serious infections. To grow in or on the body, bacteria often attach to proteins on the surface of cells that make up the lining of tissues like the gut or the throat. In some cases, bacteria use these proteins to invade the cells causing an infection. Genetic mutations in the genes encoding these proteins that protect against infection are more likely to be passed on to future generations. This may lead to rapid spread of these beneficial genes in a population.

A family of proteins called CEACAMs are frequent targets of infection-causing bacteria. These proteins have been shown to play a role in cancer progression. But they also play many helpful roles in the body, including helping transmit messages between cells, aiding cell growth, and helping the immune system recognize pathogens. Scientists are not sure if these multi-tasking CEACAM proteins can evolve to evade bacteria without affecting their other roles.

Baker et al. show that CEACAM proteins targeted by bacteria have undergone rapid evolution in primates. In the experiments, human genes encoding CEACAMs were compared with equivalent genes from 19 different primates. Baker et al. found the changes in human and primate CEACAMs often occur through a process called gene conversion. Gene conversion occurs when DNA sections are copied and pasted from one gene to another. Using laboratory experiments, they showed that some of these changes enabled CEACAM proteins to prevent certain harmful bacteria from binding.

The experiments suggest that some versions of CEACAM genes may protect humans or other primates against bacterial infections. Studies in natural populations are needed to test if this is the case. Learning more about how CEACAM proteins evolve and what they do may help scientists better understand the role they play in cancer and help improve cancer care. Studying CEACAM evolution may also help scientists understand how bacteria and other pathogens drive protein evolution in the body.

Satellite DNA-mediated diversification of a sex-ratio meiotic drive gene family in Drosophila

Sex chromosomes are susceptible to the evolution of selfish meiotic drive elements that bias transmission and distort progeny sex ratios. Conflict between such sex-ratio drivers and the rest of the genome can trigger evolutionary arms races resulting in genetically suppressed 'cryptic' drive systems. The Winters cryptic sex-ratio drive system of Drosophila simulans comprises a driver, Distorter on the X (Dox) and an autosomal suppressor, Not much yang, a retroduplicate of Dox that suppresses via production of endogenous small interfering RNAs (esiRNAs). Here we report that over 22 Dox-like (Dxl) sequences originated, amplified and diversified over the ~250,000-year history of the three closely related species, D. simulans, D. mauritiana and D. sechellia. The Dxl sequences encode a rapidly evolving family of protamines. Dxl copy numbers amplified by ectopic exchange among euchromatic islands of satellite DNAs on the X chromosome and separately spawned four esiRNA-producing suppressors on the autosomes. Our results reveal the genomic consequences of evolutionary arms races and highlight complex interactions among different classes of selfish DNAs.

Illegitimate Recombination between Duplicated Genes Generated from Recursive Polyploidizations Accelerated the Divergence of the Genus Arachis

The peanut (Arachis hypogaea L.) is the leading oil and food crop among the legume family. Extensive duplicate gene pairs generated from recursive polyploidizations with high sequence similarity could result from gene conversion, caused by illegitimate DNA recombination. Here, through synteny-based comparisons of two diploid and three tetraploid peanut genomes, we identified the duplicated genes generated from legume common tetraploidy (LCT) and peanut recent allo-tetraploidy (PRT) within genomes. In each peanut genome (or subgenomes), we inferred that 6.8–13.1% of LCT-related and 11.3–16.5% of PRT-related duplicates were affected by gene conversion, in which the LCT-related duplicates were the most affected by partial gene conversion, whereas the PRT-related duplicates were the most affected by whole gene conversion. Notably, we observed the conversion between duplicates as the long-lasting contribution of polyploidizations accelerated the divergence of different Arachis genomes. Moreover, we found that the converted duplicates are unevenly distributed across the chromosomes and are more often near the ends of the chromosomes in each genome. We also confirmed that well-preserved homoeologous chromosome regions may facilitate duplicates’ conversion. In addition, we found that these biological functions contain a higher number of preferentially converted genes, such as catalytic activity-related genes. We identified specific domains that are involved in converted genes, implying that conversions are associated with important traits of peanut growth and development.

Conversion between 100-million-year-old duplicated genes contributes to rice subspecies divergence

BACKGROUND

Duplicated gene pairs produced by ancient polyploidy maintain high sequence similarity over a long period of time and may result from illegitimate recombination between homeologous chromosomes. The genomes of Asian cultivated rice Oryza sativa ssp. indica (XI) and Oryza sativa ssp. japonica (GJ) have recently been updated, providing new opportunities for investigating ongoing gene conversion events and their impact on genome evolution.

RESULTS

Using comparative genomics and phylogenetic analyses, we evaluated gene conversion rates between duplicated genes produced by polyploidization 100 million years ago (mya) in GJ and XI. At least 5.19-5.77% of genes duplicated across the three rice genomes were affected by whole-gene conversion after the divergence of GJ and XI at ~ 0.4 mya, with more (7.77-9.53%) showing conversion of only portions of genes. Independently converted duplicates surviving in the genomes of different subspecies often use the same donor genes. The ongoing gene conversion frequency was higher near chromosome termini, with a single pair of homoeologous chromosomes, 11 and 12, in each rice genome being most affected. Notably, ongoing gene conversion has maintained similarity between very ancient duplicates, provided opportunities for further gene conversion, and accelerated rice divergence. Chromosome rearrangements after polyploidization are associated with ongoing gene conversion events, and they directly restrict recombination and inhibit duplicated gene conversion between homeologous regions. Furthermore, we found that the converted genes tended to have more similar expression patterns than nonconverted duplicates. Gene conversion affects biological functions associated with multiple genes, such as catalytic activity, implying opportunities for interaction among members of large gene families, such as NBS-LRR disease-resistance genes, contributing to the occurrence of the gene conversion.

CONCLUSION

Duplicated genes in rice subspecies generated by grass polyploidization ~ 100 mya remain affected by gene conversion at high frequency, with important implications for the divergence of rice subspecies.

High and Highly Variable Spontaneous Mutation Rates in Daphnia

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⁻⁹/bp/generation [95% CI: 6.66–11.97 × 10⁻⁹/bp/generation] in the nuclear genome and µ-m,bs = 8.7 × 10⁻⁷/bp/generation [95% CI: 4.40–15.12 × 10⁻⁷/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.

Dynamic Evolution of De Novo DNA Methyltransferases in Rodent and Primate Genomes

Transcriptional silencing of retrotransposons via DNA methylation is paramount for mammalian fertility and reproductive fitness. During germ cell development, most mammalian species utilize the de novo DNA methyltransferases DNMT3A and DNMT3B to establish DNA methylation patterns. However, many rodent species deploy a third enzyme, DNMT3C, to selectively methylate the promoters of young retrotransposon insertions in their germline. The evolutionary forces that shaped DNMT3C’s unique function are unknown. Using a phylogenomic approach, we confirm here that Dnmt3C arose through a single duplication of Dnmt3B that occurred ∼60 Ma in the last common ancestor of muroid rodents. Importantly, we reveal that DNMT3C is composed of two independently evolving segments: the latter two-thirds have undergone recurrent gene conversion with Dnmt3B, whereas the N-terminus has instead evolved under strong diversifying selection. We hypothesize that positive selection of Dnmt3C is the result of an ongoing evolutionary arms race with young retrotransposon lineages in muroid genomes. Interestingly, although primates lack DNMT3C, we find that the N-terminus of DNMT3A has also evolved under diversifying selection. Thus, the N-termini of two independent de novo methylation enzymes have evolved under diversifying selection in rodents and primates. We hypothesize that repression of young retrotransposons might be driving the recurrent innovation of a functional domain in the N-termini on germline DNMT3s in mammals.

Rapid Evolution of HERC6 and Duplication of a Chimeric HERC5/6 Gene in Rodents and Bats Suggest an Overlooked Role of HERCs in Mammalian Immunity

Studying the evolutionary diversification of mammalian antiviral defenses is of main importance to better understand our innate immune repertoire. The small HERC proteins are part of a multigene family, including HERC5 and HERC6, which have probably diversified through complex evolutionary history in mammals. Here, we performed mammalian-wide phylogenetic and genomic analyses of HERC5 and HERC6, using 83 orthologous sequences from bats, rodents, primates, artiodactyls, and carnivores—the top five representative groups of mammalian evolution. We found that HERC5 has been under weak and differential positive selection in mammals, with only primate HERC5 showing evidences of pathogen-driven selection. In contrast, HERC6 has been under strong and recurrent adaptive evolution in mammals, suggesting past and widespread genetic arms-races with viral pathogens. Importantly, the rapid evolution of mammalian HERC6 spacer domain suggests that it might be a host-pathogen interface, targeting viral proteins and/or being the target of virus antagonists. Finally, we identified a HERC5/6 chimeric gene that arose from independent duplication in rodent and bat lineages and encodes for a conserved HERC5 N-terminal domain and divergent HERC6 spacer and HECT domains. This duplicated chimeric gene highlights adaptations that potentially contribute to rodent and bat immunity. Our findings open new research avenues on the functions of HERC6 and HERC5/6 in mammals, and on their implication in antiviral innate immunity.

DGINN, an automated and highly-flexible pipeline for the detection of genetic innovations on protein-coding genes

Adaptive evolution has shaped major biological processes. Finding the protein-coding genes and the sites that have been subjected to adaptation during evolutionary time is a major endeavor. However, very few methods fully automate the identification of positively selected genes, and widespread sources of genetic innovations such as gene duplication and recombination are absent from most pipelines. Here, we developed DGINN, a highly-flexible and public pipeline to Detect Genetic INNovations and adaptive evolution in protein-coding genes. DGINN automates, from a gene's sequence, all steps of the evolutionary analyses necessary to detect the aforementioned innovations, including the search for homologs in databases, assignation of orthology groups, identification of duplication and recombination events, as well as detection of positive selection using five methods to increase precision and ranking of genes when a large panel is analyzed. DGINN was validated on nineteen genes with previously-characterized evolutionary histories in primates, including some engaged in host-pathogen arms-races. Our results confirm and also expand results from the literature, including novel findings on the Guanylate-binding protein family, GBPs. This establishes DGINN as an efficient tool to automatically detect genetic innovations and adaptive evolution in diverse datasets, from the user's gene of interest to a large gene list in any species range.