Search

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.

About one-quarter of the euchromatic portion of the male-specific region of the human Y chromosome consists of large duplicated sequences that are organized in eight palindromes (termed P1–P8), which undergo arm-to arm gene conversion, a proposed mechanism for maintaining their sequence integrity. Although the relevance of gene conversion in the evolution of palindromic sequences has been profoundly recognized, the dynamic of this mechanism is still nuanced. To shed light into the evolution of these genomic elements, we performed a high-depth (50×) targeted next-generation sequencing of the palindrome P6 in 157 subjects belonging to the most divergent evolutionary lineages of the Y chromosome. We found 118 new paralogous sequence variants, which were placed into the context of a robust Y chromosome phylogeny based on 7240 SNPs of the X-degenerate region. We mapped along the phylogeny 80 gene conversion events that shaped the diversity of P6 arms during recent human history. In contrast to previous studies, we demonstrated that arm-to-arm gene conversion, which occurs at a rate of 6.01 × 10 −6 conversions/base/year, is not biased toward the retention of the ancestral state of sequences. We also found a significantly lower mutation rate of the arms (6.18 × 10−10 mutations/base/year) compared with the spacer (9.16 × 10−10 mutations/base/year), a finding that may explain the observed higher inter-species conservation of arms, without invoking any bias of conversion. Finally, by formally testing the mutation/conversion balance in P6, we found that the arms of this palindrome reached a steady-state equilibrium between mutation and gene conversion.

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.

Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Nav) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Nav gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Nav paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 million years. TTX-bearing newts possess additional unique substitutions across the entire Nav gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Nav1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Nav1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.

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.

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.

Sex chromosomes often carry large nonrecombining regions that can extend progressively over time, generating evolutionary strata of sequence divergence. However, some sex chromosomes display an incomplete suppression of recombination. Large genomic regions without recombination and evolutionary strata have also been documented around fungal mating-type loci, but have been studied in only a few fungal systems. In the model fungus Podospora anserina (Ascomycota, Sordariomycetes), the reference S strain lacks recombination across a 0.8-Mb region around the mating-type locus. The lack of recombination in this region ensures that nuclei of opposite mating types are packaged into a single ascospore (pseudohomothallic lifecycle). We found evidence for a lack of recombination around the mating-type locus in the genomes of ten P. anserina strains and six closely related pseudohomothallic Podospora species. Importantly, the size of the nonrecombining region differed between strains and species, as indicated by the heterozygosity levels around the mating-type locus and experimental selfing. The nonrecombining region is probably labile and polymorphic, differing in size and precise location within and between species, resulting in occasional, but infrequent, recombination at a given base pair. This view is also supported by the low divergence between mating types, and the lack of strong linkage disequilibrium, chromosomal rearrangements, transspecific polymorphism and genomic degeneration. We found a pattern suggestive of evolutionary strata in P. pseudocomata. The observed heterozygosity levels indicate low but nonnull outcrossing rates in nature in these pseudohomothallic fungi. This study adds to our understanding of mating-type chromosome evolution and its relationship to mating systems.

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.

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.

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 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.

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.

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.

Genetic stability depends on the maintenance of a variety of chromosome structures and the precise repair of DNA breaks. During meiosis, programmed double-strand breaks (DSBs) made in prophase I are normally repaired as gene conversions or crossovers. DSBs can also be made by other mechanisms, such as the movement of transposable elements (TEs), which must also be resolved. Incorrect repair of these DNA lesions can lead to mutations, copy-number changes, translocations, and/or aneuploid gametes. In Drosophila melanogaster, as in most organisms, meiotic DSB repair occurs in the presence of a rapidly evolving multiprotein structure called the synaptonemal complex (SC). Here, whole-genome sequencing is used to investigate the fate of meiotic DSBs in D. melanogaster mutant females lacking functional SC, to assay for de novo CNV formation, and to examine the role of the SC in transposable element movement in flies. The data indicate that, in the absence of SC, copy-number variation still occurs and meiotic DSB repair by gene conversion occurs infrequently. Remarkably, an 856-kilobase de novo CNV was observed in two unrelated individuals of different genetic backgrounds and was identical to a CNV recovered in a previous wild-type study, suggesting that recurrent formation of large CNVs occurs in Drosophila. In addition, the rate of novel TE insertion was markedly higher than wild type in one of two SC mutants tested, suggesting that SC proteins may contribute to the regulation of TE movement and insertion in the genome. Overall, this study provides novel insight into the role that the SC plays in genome stability and provides clues as to why the sequence, but not structure, of SC proteins is rapidly evolving.

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.

Since the inception of the molecular clock model for sequence evolution, the investigation of protein divergence has revolved around the question of a more or less constant change of amino acid sequences, with specific overall rates for each family. Although anomalies in clock-like divergence are well known, the assumption of a constant decay rate for a given protein family is usually taken as the null model for protein evolution. However, systematic tests of this null model at a genome-wide scale have lagged behind, despite the databases’ enormous growth. We focus here on divergence rate comparisons between very closely related lineages since this allows clear orthology assignments by synteny and reliable alignments, which are crucial for determining substitution rate changes. We generated a high-confidence dataset of syntenic orthologs from four ape species, including humans. We find that despite the appearance of an overall clock-like substitution pattern, several hundred protein families show lineage-specific acceleration and deceleration in divergence rates, or combinations of both in different lineages. Hence, our analysis uncovers a rather dynamic history of substitution rate changes, even between these closely related lineages, implying that one should expect that a large fraction of proteins will have had a history of episodic rate changes in deeper phylogenies. Furthermore, each of the lineages has a separate set of particularly fast diverging proteins. The genes with the highest percentage of branch-specific substitutions are ADCYAP1 in the human lineage (9.7%), CALU in chimpanzees (7.1%), SLC39A14 in the internal branch leading to humans and chimpanzees (4.1%), RNF128 in gorillas (9%), and S100Z in gibbons (15.2%). The mutational pattern in ADCYAP1 suggests a biased mutation process, possibly through asymmetric gene conversion effects. We conclude that a null model of constant change can be problematic for predicting the evolutionary trajectories of individual proteins.

Highly efficient gene conversion systems have the potential to facilitate the study of complex genetic traits using laboratory mice and, if implemented as a “gene drive,” to limit loss of biodiversity and disease transmission caused by wild rodent populations. We previously showed that such a system of gene conversion from heterozygous to homozygous after a sequence targeted CRISPR/Cas9 double-strand DNA break (DSB) is feasible in the female mouse germline. In the male germline, however, all DSBs were instead repaired by end joining (EJ) mechanisms to form an “insertion/deletion” (indel) mutation. These observations suggested that timing Cas9 expression to coincide with meiosis I is critical to favor conditions when homologous chromosomes are aligned and interchromosomal homology-directed repair (HDR) mechanisms predominate. Here, using a Cas9 knock-in allele at the Spo11 locus, we show that meiotic expression of Cas9 does indeed mediate gene conversion in the male as well as in the female germline. However, the low frequency of both HDR and indel mutation in both male and female germlines suggests that Cas9 may be expressed from the Spo11 locus at levels too low for efficient DSB formation. We suggest that more robust Cas9 expression initiated during early meiosis I may improve the efficiency of gene conversion and further increase the rate of “super-mendelian” inheritance from both male and female mice., This study shows that while Cas9 expression during meiosis I promotes genotype conversion - the mechanism underlying CRISPR ’gene drive’ - in both male and female mice, timing and high levels of Cas9 protein are critical to achieve robust efficiency.

Antigenic variation is an immune evasion strategy used by Trypanosoma brucei that results in the periodic exchange of the surface protein coat. This process is facilitated by the movement of variant surface glycoprotein genes in or out of a specialized locus known as bloodstream form expression site by homologous recombination, facilitated by blocks of repetitive sequence known as the 70-bp repeats, that provide homology for gene conversion events. DNA double strand breaks are potent drivers of antigenic variation, however where these breaks must fall to elicit a switch is not well understood. To understand how the position of a break influences antigenic variation we established a series of cell lines to study the effect of an I-SceI meganuclease break in the active expression site. We found that a DNA break within repetitive regions is not productive for VSG switching, and show that the break position leads to a distinct gene expression profile and DNA repair response which dictates how antigenic variation proceeds in African trypanosomes., Author summary Crucial to triggering antigenic variation is the formation of DNA double strand breaks (DSB). These lesions have been shown to be potent drivers of variant surface glycoprotein (VSG) switching, albeit highly toxic. Trypanosomes immune evasion strategy relies on their ability to rapidly exchange the singly expressed VSG for one that is antigenically distinct. It has been previously shown that the subtelomeric ends, here the locus from which the VSG is expressed, accumulate DSBs. Using the I-SceI meganuclease system we established a series of cell lines to assess how the position of a DSB influences antigenic variation and the cellular response to a break. We show that a DSB in highly repetitive regions are poor triggers for antigenic variation. Contrastingly, a DSB that does lead to VSG switching via recombination results in the upregulation of DNA damage linked genes. Our results provide new insights into how the position of a DSB influences repair pathway choice and the subsequent gene expression changes. We propose that where repair is not dominated by recombination, but rather by an error prone mechanism, silent BES promoters are partially activated to facilitate rapid transcriptional switching should repair be deleterious to the cell.

Homologous recombination (HR) typically occurs during meiosis between homologs, at a few unplanned locations along the chromosomes. In this study, we tested whether targeted recombination between homologous chromosomes can be achieved via Clustered Regulatory Interspaced Short Palindromic Repeat associated protein Cas9 (CRISPR-Cas9)-induced DNA double-strand break (DSB) repair in Arabidopsis thaliana. Our experimental system includes targets for DSB induction in euchromatic and heterochromatic genomic regions of hybrid F1 plants, in one or both parental chromosomes, using phenotypic and molecular markers to measure Non-Homologous End Joining and HR repair. We present a series of evidence showing that targeted DSBs can be repaired via HR using a homologous chromosome as the template in various chromatin contexts including in pericentric regions. Targeted crossover was rare, but gene conversion events were the most frequent outcome of HR and were found in both “hot and cold” regions. The length of the conversion tracts was variable, ranging from 5 to 7505 bp. In addition, a typical feature of these tracks was that they often were interrupted. Our findings pave the way for the use of targeted gene-conversion for precise breeding.

Immune evasion by Treponema pallidum subspecies pallidum (T. pallidum) has been attributed to antigenic variation of its putative outer-membrane protein TprK. In TprK, amino acid diversity is confined to seven variable (V) regions, and generation of sequence diversity within the V regions occurs via a non-reciprocal segmental gene conversion mechanism where donor cassettes recombine into the tprK expression site. Although previous studies have shown the significant role of immune selection in driving accumulation of TprK variants, the contribution of baseline gene conversion activity to variant diversity is less clear. Here, combining longitudinal tprK deep sequencing of near clonal Chicago C from immunocompetent and immunosuppressed rabbits along with the newly developed in vitro cultivation system for T. pallidum, we directly characterized TprK alleles in the presence and absence of immune selection. Our data confirm significantly greater sequence diversity over time within the V6 region during syphilis infection in immunocompetent rabbits compared to immunosuppressed rabbits, consistent with previous studies on the role of TprK in evasion of the host immune response. Compared to strains grown in immunocompetent rabbits, strains passaged in vitro displayed low level changes in allele frequencies of TprK variable region sequences similar to that of strains passaged in immunosuppressed rabbits. Notably, we found significantly increased rates of V6 allele generation relative to other variable regions in in vitro cultivated T, pallidum strains, illustrating that the diversity within these hypervariable regions occurs in the complete absence of immune selection. Together, our results demonstrate antigenic variation in T. pallidum can be studied in vitro and occurs even in the complete absence of immune pressure, allowing the T. pallidum population to continuously evade the immune system of the infected host., Author summary Syphilis continues to be a disease of global and public health concern, even though the infection can be easily diagnosed and effectively treated with penicillin. Although infected individuals often develop a strong immunity to the pathogen, repeated infection with the syphilis agent, Treponema pallidum subspecies pallidum (T. pallidum), is possible. Several studies point at antigenic variation of the T. pallidum TprK protein as the mechanism responsible for evasion of the immunity that develops during infection, pathogen persistence, and re-infection. Past studies have highlighted the importance of immune clearance of dominant variants that, in turn, allows less represented variants to emerge. The contribution of an immunity-independent baseline generation of variability in the tprK gene is less clear. Here, we used deep sequencing to profile tprK variants using a laboratory-isolated T. pallidum strain nearly isogenic for tprK that was propagated over time in vitro, where no immune pressure is exerted on the pathogen, as well as in samples obtained from immunosuppressed and immunocompetent rabbits infected with the same strain. We confirmed that tprK accumulates significantly more diversity under immune pressure, and demonstrated a low but discernible basal rate of gene conversion in complete absence of immune pressure.

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 that involve branch migration as an integral feature and that can be separated from break repair 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.

Population genetics studies often make use of a class of nucleotide site free from selective pressures, in order to make inferences about population size changes or natural selection at other sites. If such neutral sites can be identified, they offer the opportunity to avoid any confounding effects of selection. Here, we investigate evolution at putatively neutrally evolving short intronic sites in natural populations of Drosophila melanogaster and Drosophila simulans, in order to understand the properties of spontaneous mutations and the extent of GC-biased gene conversion in these species. Use of data on the genetics of natural populations is advantageous because it integrates information from large numbers of individuals over long timescales. In agreement with direct evidence from observations of spontaneous mutations in Drosophila, we find a bias in the spectrum of mutations toward AT basepairs. In addition, we find that this bias is stronger in the D. melanogaster lineage than in the D. simulans lineage. The evidence for GC-biased gene conversion in Drosophila has been equivocal. Here, we provide evidence for a weak force favoring GC in both species, which is correlated with the GC content of introns and is stronger in D. simulans than in D. melanogaster.

Phylogenetic codon models are routinely used to characterize selective regimes in coding sequences. Their parametric design, however, is still a matter of debate, in particular concerning the question of how to account for differing nucleotide frequencies and substitution rates. This problem relates to the fact that nucleotide composition in protein-coding sequences is the result of the interactions between mutation and selection. In particular, because of the structure of the genetic code, the nucleotide composition differs between the three coding positions, with the third position showing a more extreme composition. Yet, phylogenetic codon models do not correctly capture this phenomenon and instead predict that the nucleotide composition should be the same for all three positions. Alternatively, some models allow for different nucleotide rates at the three positions, an approach conflating the effects of mutation and selection on nucleotide composition. In practice, it results in inaccurate estimation of the strength of selection. Conceptually, the problem comes from the fact that phylogenetic codon models do not correctly capture the fixation bias acting against the mutational pressure at the mutation–selection equilibrium. To address this problem and to more accurately identify mutation rates and selection strength, we present an improved codon modeling approach where the fixation rate is not seen as a scalar, but as a tensor. This approach gives an accurate representation of how mutation and selection oppose each other at equilibrium and yields a reliable estimate of the mutational process, while disentangling the mean fixation probabilities prevailing in different mutational directions.

Among extant cetaceans, mysticetes are filter feeders that do not possess teeth and use their baleen for feeding, while most odontocetes are considered suction feeders, which capture prey by suction without biting or chewing with teeth. In the present study, we address the functionality of amelogenin (AMEL) genes in cetaceans. AMEL encodes a protein that is specifically involved in dental enamel formation and is located on the sex chromosomes in eutherians. The X-copy AMELX is functional in enamel-bearing eutherians, whereas the Y-copy AMELY appears to have undergone decay and was completely lost in some species. Consistent with these premises, we detected various deleterious mutations and/or non-canonical splice junctions in AMELX of mysticetes and four suction feeding odontocetes, Delphinapterus leucas, Monodon monoceros, Kogia breviceps, and Physeter macrocephalus, and in AMELY of mysticetes and odontocetes. Regardless of the functionality, both AMELX and AMELY are equally and unusually small in cetaceans, and even their functional AMELX genes presumably encode a degenerate core region, which is thought to be essential for enamel matrix assembly and enamel crystal growth. Furthermore, our results suggest that the most recent common ancestors of extant cetaceans had functional AMELX and AMELY, both of which are similar to AMELX of Platanista minor. Similar small AMELX and AMELY in archaic cetaceans can be explained by gene conversion between AMELX and AMELY. We speculate that common ancestors of modern cetaceans employed a degenerate AMELX, transferred from a decaying AMELY by gene conversion, at an early stage of their transition to suction feeders.

Meiotic recombination is a fundamental biological process that leads to crossover and gene conversion. High resolution maps of meiotic recombination have been reported in several model organisms. However, few studies have studied how rapidly selection affects the products of meiotic recombination. Here we constructed and sequenced a yeast population of 38 haploid strains derived from hybridizations of two common used strains S288c and YJM789. We identified 20 regions with strong biased allele frequency across the genome, revealing signatures of selection in a rather short period. These regions harbor ample crossovers and gene conversions, which enable us to trace how selection works on the genomic fragments after meiosis. The total length of such regions under selection accounts for 5% of the entire genome, and those regions contain many functional-related genes. In addition, recombination breaks down linkage disequilibrium to half of its maximum within 42 kb and reduces nucleotide diversity significantly in selected regions. Our study thus provides details of directional selection on the outcomes of meiotic recombination using experimental approaches, and will shed light on our understanding of the fast reshaping of population structure by selection, as well as the important roles of recombination.

Ribosomal protein (RP) genes encode structural components of ribosomes, the cellular machinery for protein synthesis. A single functional copy has been maintained in most of 78–80 RP families in animals due to evolutionary constraints imposed by gene dosage balance. Some fungal species have maintained duplicate copies in most RP families. The mechanisms by which the RP genes were duplicated and maintained and their functional significance are poorly understood. To address these questions, we identified all RP genes from 295 fungi and inferred the timing and nature of gene duplication events for all RP families. We found that massive duplications of RP genes have independently occurred by different mechanisms in three distantly related lineages: budding yeasts, fission yeasts, and Mucoromycota. The RP gene duplicates in budding yeasts and Mucoromycota were mainly created by whole genome duplication events. However, duplicate RP genes in fission yeasts were likely generated by retroposition, which is unexpected considering their dosage sensitivity. The sequences of most RP paralogs have been homogenized by repeated gene conversion in each species, demonstrating parallel concerted evolution, which might have facilitated the retention of their duplicates. Transcriptomic data suggest that the duplication and retention of RP genes increased their transcript abundance. Physiological data indicate that increased ribosome biogenesis allowed these organisms to rapidly consume sugars through fermentation while maintaining high growth rates, providing selective advantages to these species in sugar-rich environments.

Copy number variants (CNVs), deletions and duplications of segments of DNA, account for at least five times more variable base pairs in humans than single-nucleotide variants. Several common CNVs were shown to change coding and regulatory sequences and thus dramatically affect adaptive phenotypes involving immunity, perception, metabolism, skin structure, among others. Some of these CNVs were also associated with susceptibility to cancer, infection, and metabolic disorders. These observations raise the possibility that CNVs are a primary contributor to human phenotypic variation and consequently evolve under selective pressures. Indeed, locus-specific haplotype-level analyses revealed signatures of natural selection on several CNVs. However, more traditional tests of selection which are often applied to single-nucleotide variation often have diminished statistical power when applied to CNVs because they often do not show strong linkage disequilibrium with nearby variants. Recombination-based formation mechanisms of CNVs lead to frequent recurrence and gene conversion events, breaking the linkage disequilibrium involving CNVs. Similar methodological challenges also prevent routine genome-wide association studies to adequately investigate the impact of CNVs on heritable human disease. Thus, we argue that the full relevance of CNVs to human health and evolution is yet to be elucidated. We further argue that a holistic investigation of formation mechanisms within an evolutionary framework would provide a powerful framework to understand the functional and biomedical impact of CNVs. In this paper, we review several cases where studies reveal diverse evolutionary histories and unexpected functional consequences of CNVs. We hope that this review will encourage further work on CNVs by both evolutionary and medical geneticists.

Published genomes frequently contain erroneous gene models that represent issues associated with identification of open reading frames, start sites, splice sites, and related structural features. The source of these inconsistencies is often traced back to integration across text file formats designed to describe long read alignments and predicted gene structures. In addition, the majority of gene prediction frameworks do not provide robust downstream filtering to remove problematic gene annotations, nor do they represent these annotations in a format consistent with current file standards. These frameworks also lack consideration for functional attributes, such as the presence or absence of protein domains that can be used for gene model validation. To provide oversight to the increasing number of published genome annotations, we present a software package, the Gene Filtering, Analysis, and Conversion (gFACs), to filter, analyze, and convert predicted gene models and alignments. The software operates across a wide range of alignment, analysis, and gene prediction files with a flexible framework for defining gene models with reliable structural and functional attributes. gFACs supports common downstream applications, including genome browsers, and generates extensive details on the filtering process, including distributions that can be visualized to further assess the proposed gene space. gFACs is freely available and implemented in Perl with support from BioPerl libraries at https://gitlab.com/PlantGenomicsLab/gFACs. Keywords: Genome annotation, Bioinformatics, Protein annotation, Gene prediction, Alignment

Natural selection works best when the two alleles in a diploid organism are transmitted to offspring at equal frequencies. Despite this, selfish loci known as meiotic drivers that bias their own transmission into gametes are found throughout eukaryotes. Drive is thought to be a powerful evolutionary force, but empirical evolutionary analyses of drive systems are limited by low numbers of identified meiotic drive genes. Here, we analyze the evolution of the wtf gene family of Schizosaccharomyces pombe that contains both killer meiotic drive genes and suppressors of drive. We completed assemblies of all wtf genes for two S. pombe isolates, as well as a subset of wtf genes from over 50 isolates. We find that wtf copy number can vary greatly between isolates and that amino acid substitutions, expansions and contractions of DNA sequence repeats, and nonallelic gene conversion between family members all contribute to dynamic wtf gene evolution. This work demonstrates the power of meiotic drive to foster rapid evolution and identifies a recombination mechanism through which transposons can indirectly mobilize meiotic drivers.

As limits on O2 availability during submergence impose severe constraints on aerobic respiration, the oxygen binding globin proteins of marine mammals are expected to have evolved under strong evolutionary pressures during their land-to-sea transition. Here, we address this question for the order Sirenia by retrieving, annotating, and performing detailed selection analyses on the globin repertoire of the extinct Steller’s sea cow (Hydrodamalis gigas), dugong (Dugong dugon), and Florida manatee (Trichechus manatus latirostris) in relation to their closest living terrestrial relatives (elephants and hyraxes). These analyses indicate most loci experienced elevated nucleotide substitution rates during their transition to a fully aquatic lifestyle. While most of these genes evolved under neutrality or strong purifying selection, the rate of nonsynonymous/synonymous replacements increased in two genes (Hbz-T1 and Hba-T1) that encode the α-type chains of hemoglobin (Hb) during each stage of life. Notably, the relaxed evolution of Hba-T1 is temporally coupled with the emergence of a chimeric pseudogene (Hba-T2/Hbq-ps) that contributed to the tandemly linked Hba-T1 of stem sirenians via interparalog gene conversion. Functional tests on recombinant Hb proteins from extant and ancestral sirenians further revealed that the molecular remodeling of Hba-T1 coincided with increased Hb–O2 affinity in early sirenians. Available evidence suggests that this trait evolved to maximize O2 extraction from finite lung stores and suppress tissue O2 offloading, thereby facilitating the low metabolic intensities of extant sirenians. In contrast, the derived reduction in Hb–O2 affinity in (sub)Arctic Steller’s sea cows is consistent with fueling increased thermogenesis by these once colossal marine herbivores.

Agrocybe aegerita is a cultivated edible mushroom in numerous countries, which also serves as a model basidiomycete to study fruiting body formation. Aiming to create an easily expandable customised molecular toolset for transformation and constitutive gene of interest expression, we first created a homologous dominant marker for transformant selection. Progeny monokaryons of the genome-sequenced dikaryon A. aegerita AAE-3 used here were identified as sensitive to the systemic fungicide carboxin. We cloned the wild-type gene encoding the iron-sulphur protein subunit of succinate dehydrogenase AaeSdi1 including its up- and downstream regions, and introduced a single-point mutation (His237 to Leu) to make it confer carboxin resistance. PEG-mediated transformation of protoplasts derived from either oidia or vegetative monokaryotic mycelium with the resulting carboxin resistance marker (CbxR) plasmid pSDI1E3 yielded carboxin-resistant transformants in both cases. Plasmid DNA linearised within the selection marker resulted in transformants with ectopic multiple insertions of plasmid DNA in a head-to-tail repeat-like fashion. When circular plasmid was used, ectopic single integration into the fungal genome was favoured, but also gene conversion at the homologous locus was seen in 1 out of 11 analysed transformants. Employing CbxR as selection marker, two versions of a reporter gene construct were assembled via Golden Gate cloning which allows easy recombination of its modules. These consisted of an eGFP expression cassette controlled by the native promoter PAaeGPDII and the heterologous terminator Tnos, once with and once without an intron in front of the eGFP start codon. After protoplast transformation with either construct as circular plasmid DNA, GFP fluorescence was detected with either transformants, indicating that expression of eGFP is intron-independent in A. aegerita. This paves the way for functional genetics approaches to A. aegerita, e.g., via constitutive expression of fruiting-related genes.

The sex chromosomes often follow unusual evolutionary trajectories. In particular, the sex-limited chromosomes frequently exhibit a small but unusual gene content in numerous species, where many genes have undergone massive gene amplification. The reasons for this remain elusive with a number of recent studies implicating meiotic drive, sperm competition, genetic drift, and gene conversion in the expansion of gene families. However, our understanding is primarily based on Y chromosome studies as few studies have systematically tested for copy number variation on W chromosomes. Here, we conduct a comprehensive investigation into the abundance, variability, and evolution of ampliconic genes on the avian W. First, we quantified gene copy number and variability across the duck W chromosome. We find a limited number of gene families as well as conservation in W-linked gene copy number across duck breeds, indicating that gene amplification may not be such a general feature of sex chromosome evolution as Y studies would initially suggest. Next, we investigated the evolution of HINTW, a prominent ampliconic gene family hypothesized to play a role in female reproduction and oogenesis. In particular, we investigated the factors driving the expansion of HINTW using contrasts between modern chicken and duck breeds selected for different female-specific selection regimes and their wild ancestors. Although we find the potential for selection related to fecundity in explaining small-scale gene amplification of HINTW in the chicken, purifying selection seems to be the dominant mode of evolution in the duck. Together, this challenges the assumption that HINTW is key for female fecundity across the avian phylogeny.

Gene conversion is GC-biased across a wide range of taxa. Large palindromes on mammalian sex chromosomes undergo frequent gene conversion that maintains arm-to-arm sequence identity greater than 99%, which may increase their susceptibility to the effects of GC-biased gene conversion. Here, we demonstrate a striking history of GC-biased gene conversion in 12 palindromes conserved on the X chromosomes of human, chimpanzee, and rhesus macaque. Primate X-chromosome palindrome arms have significantly higher GC content than flanking single-copy sequences. Nucleotide replacements that occurred in human and chimpanzee palindrome arms over the past 7 million years are one-and-a-half times as GC-rich as the ancestral bases they replaced. Using simulations, we show that our observed pattern of nucleotide replacements is consistent with GC-biased gene conversion with a magnitude of 70%, similar to previously reported values based on analyses of human meioses. However, GC-biased gene conversion since the divergence of human and rhesus macaque explains only a fraction of the observed difference in GC content between palindrome arms and flanking sequence, suggesting that palindromes are older than 29 million years and/or had elevated GC content at the time of their formation. This work supports a greater than 2:1 preference for GC bases over AT bases during gene conversion and demonstrates that the evolution and composition of mammalian sex chromosome palindromes is strongly influenced by GC-biased gene conversion.

DNA synthesis during homologous recombination is highly mutagenic and prone to template switches. Two‐ended DNA double‐strand breaks (DSBs) are usually repaired by gene conversion with a short patch of DNA synthesis, thus limiting the mutation load to the vicinity of the DSB. Single‐ended DSBs are repaired by break‐induced replication (BIR), which involves extensive and mutagenic DNA synthesis spanning up to hundreds of kilobases. It remains unknown how mutagenic BIR is suppressed at two‐ended DSBs. Here, we demonstrate that BIR is suppressed at two‐ended DSBs by proteins coordinating the usage of two ends of a DSB: (i) ssDNA annealing proteins Rad52 and Rad59 that promote second end capture, (ii) D‐loop unwinding helicase Mph1, and (iii) Mre11‐Rad50‐Xrs2 complex that promotes synchronous resection of two ends of a DSB. Finally, BIR is also suppressed when Sir2 silences a normally heterochromatic repair template. All of these proteins are particularly important for limiting BIR when recombination occurs between short repetitive sequences, emphasizing the significance of these mechanisms for species carrying many repetitive elements such as humans.

Next-generation sequencing studies are dependent on a high-quality reference genome for single nucleotide variant (SNV) calling. Although the two most recent builds of the human genome are widely used, position information is typically not directly comparable between them. Re-alignment gives the most accurate position information, but this procedure is often computationally expensive, and therefore, tools such as liftOver and CrossMap are used to convert data from one build to another. However, the positions of converted SNVs do not always match SNVs derived from aligned data, and in some instances, SNVs are known to change chromosome when converted. This is a significant problem when compiling sequencing resources or comparing results across studies. Here, we describe a novel algorithm to identify positions that are unstable when converting between human genome reference builds. These positions are detected independent of the conversion tools and are determined by the chain files, which provide a mapping of contiguous positions from one build to another. We also provide the list of unstable positions for converting between the two most commonly used builds GRCh37 and GRCh38. Pre-excluding SNVs at these positions, prior to conversion, results in SNVs that are stable to conversion. This simple procedure gives the same final list of stable SNVs as applying the algorithm and subsequently removing variants at unstable positions. This work highlights the care that must be taken when converting SNVs between genome builds and provides a simple method for ensuring higher confidence converted data. Unstable positions and algorithm code, available at https://github.com/cathaloruaidh/genomeBuildConversion