Search

Mendel proposed that the heritable material is particulate and that transmission of alleles is unbiased. An assumption of unbiased transmission was necessary to show how variation can be preserved in the absence of selection, so overturning an early objection to Darwinism. In the second half of the twentieth century, it was widely recognised that even strongly deleterious alleles can invade if they have strongly biased transmission (i.e. strong segregation distortion). The spread of alleles with distorted segregation can explain many curiosities. More recently, the selectionist-neutralist duopoly was broken by the realisation that biased gene conversion can explain phenomena such as mammalian isochore structures. An initial focus on unbiased transmission in 1919, has thus given way to an interest in biased transmission in 2019. A focus on very weak bias is now possible owing to technological advances, although technical biases may put a limit on resolving power. To understand the relevance of weak bias we could profit from having the concept of the effectively Mendelian allele, a companion to the effectively neutral allele. Understanding the implications of unbiased and biased transmission may, I suggest, be a good way to teach evolution so as to avoid psychological biases.

Azolla is the only plant with a co-evolving nitrogen-fixing (diazotrophic) cyanobacterial symbiont (cyanobiont), Nostoc azollae, resulting from whole-genome duplication (WGD) 80 million years ago in Azolla's ancestor. Additional genes from the WGD resulted in genetic, biochemical, and morphological changes in the plant that enabled the transmission of the cyanobiont to successive generations via its megaspores. The resulting permanent symbiosis and co-evolution led to the loss, downregulation, or conversion of non-essential genes to pseudogenes in the cyanobiont, changing it from a free-living organism to an obligate symbiont. The upregulation of other genes in the cyanobiont increased its atmospheric dinitrogen fixation and the provision of nitrogen-based products to the plant. As a result, Azolla can double its biomass in less than two days free-floating on fresh water and sequester large amounts of atmospheric CO2, giving it the potential to mitigate anthropogenic climate change through carbon capture and storage. Azolla's biomass can also provide local, low-cost food, biofertiliser, feed, and biofuel that are urgently needed as our population increases by a billion every twelve years. This paper integrates data from biology, genetics, geology, and palaeontology to identify the location, timing and mechanism for the acquisition of a co-evolving diazotrophic cyanobiont by Azolla's ancestor in the Late Cretaceous (Campanian) of North America. [ABSTRACT FROM AUTHOR]

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

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

Single-nucleotide variants (SNVs) in segmental duplications (SDs) have not been systematically assessed because of the limitations of mapping short-read sequencing data1,2. Here we constructed 1:1 unambiguous alignments spanning high-identity SDs across 102 human haplotypes and compared the pattern of SNVs between unique and duplicated regions3,4. We find that human SNVs are elevated 60% in SDs compared to unique regions and estimate that at least 23% of this increase is due to interlocus gene conversion (IGC) with up to 4.3 megabase pairs of SD sequence converted on average per human haplotype. We develop a genome-wide map of IGC donors and acceptors, including 498 acceptor and 454 donor hotspots affecting the exons of about 800 protein-coding genes. These include 171 genes that have 'relocated' on average 1.61 megabase pairs in a subset of human haplotypes. Using a coalescent framework, we show that SD regions are slightly evolutionarily older when compared to unique sequences, probably owing to IGC. SNVs in SDs, however, show a distinct mutational spectrum: a 27.1% increase in transversions that convert cytosine to guanine or the reverse across all triplet contexts and a 7.6% reduction in the frequency of CpG-associated mutations when compared to unique DNA. We reason that these distinct mutational properties help to maintain an overall higher GC content of SD DNA compared to that of unique DNA, probably driven by GC-biased conversion between paralogous sequences5,6.

CRISPR-based active genetic elements, or gene-drives, copied via homology-directed repair (HDR) in the germline, are transmitted to progeny at super-Mendelian frequencies. Active genetic elements also can generate widespread somatic mutations, but the genetic basis for such phenotypes remains uncertain. It is generally assumed that such somatic mutations are generated by non-homologous end-joining (NHEJ), the predominant double stranded break repair pathway active in somatic cells. Here, we develop CopyCatcher systems in Drosophila to detect and quantify somatic gene conversion (SGC) events. CopyCatchers inserted into two independent genetic loci reveal unexpectedly high rates of SGC in the Drosophila eye and thoracic epidermis. Focused RNAi-based genetic screens identify several unanticipated loci altering SGC efficiency, one of which (c-MYC), when downregulated, promotes SGC mediated by both plasmid and homologous chromosome-templates in human HEK293T cells. Collectively, these studies suggest that CopyCatchers can serve as effective discovery platforms to inform potential gene therapy strategies.

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.

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

Gene families underlie genetic innovation and phenotypic diversification. However, our understanding of the early genomic and functional evolution of tandemly arranged gene families remains incomplete as paralog sequence similarity hinders their accurate characterization. The Drosophila melanogaster-specific gene family Sdic is tandemly repeated and impacts sperm competition. We scrutinized Sdic in 20 geographically diverse populations using reference-quality genome assemblies, read-depth methodologies, and qPCR, finding that ∼90% of the individuals harbor 3-7 copies as well as evidence of population differentiation. In strains with reliable gene annotations, copy number variation (CNV) and differential transposable element insertions distinguish one structurally distinct version of the Sdic region per strain. All 31 annotated copies featured protein-coding potential and, based on the protein variant encoded, were categorized into 13 paratypes differing in their 3' ends, with 3-5 paratypes coexisting in any strain examined. Despite widespread gene conversion, the only copy present in all strains has functionally diverged at both coding and regulatory levels under positive selection. Contrary to artificial tandem duplications of the Sdic region that resulted in increased male expression, CNV in cosmopolitan strains did not correlate with expression levels, likely as a result of differential genome modifier composition. Duplicating the region did not enhance sperm competitiveness, suggesting a fitness cost at high expression levels or a plateau effect. Beyond facilitating a minimally optimal expression level, Sdic CNV acts as a catalyst of protein and regulatory diversity, showcasing a possible evolutionary path recently formed tandem multigene families can follow toward long-term consolidation in eukaryotic genomes.

Genetic conflicts arise when the evolutionary interests of two genetic elements are not aligned. Conflicts between genomes (e.g. pathogen versus host) or within the same genome (e.g. internal parasitic DNA sequences versus the rest of the host genome) can both foster 'molecular arms races', in which genes on both sides of the conflict rapidly evolve due to bouts of adaptation and counter-adaptation. Importantly, a source of genetic novelty is needed to fuel these arms races. In this review, we highlight gene conversion as a major force in generating the novel alleles on which selection can act. Using examples from both intergenomic and intragenomic conflicts, we feature the mechanisms by which gene conversion facilitates the rapid evolution of genes in conflict.

A gene drive biases the transmission of one of the two copies of a gene such that it is inherited more frequently than by random segregation. Highly efficient gene drive systems have recently been developed in insects, which leverage the sequence-targeted DNA cleavage activity of CRISPR-Cas9 and endogenous homology-directed repair mechanisms to convert heterozygous genotypes to homozygosity1-4. If implemented in laboratory rodents, similar systems would enable the rapid assembly of currently impractical genotypes that involve multiple homozygous genes (for example, to model multigenic human diseases). To our knowledge, however, such a system has not yet been demonstrated in mammals. Here we use an active genetic element that encodes a guide RNA, which is embedded in the mouse tyrosinase (Tyr) gene, to evaluate whether targeted gene conversion can occur when CRISPR-Cas9 is active in the early embryo or in the developing germline. Although Cas9 efficiently induces double-stranded DNA breaks in the early embryo and male germline, these breaks are not corrected by homology-directed repair. By contrast, Cas9 expression limited to the female germline induces double-stranded breaks that are corrected by homology-directed repair, which copies the active genetic element from the donor to the receiver chromosome and increases its rate of inheritance in the next generation. These results demonstrate the feasibility of CRISPR-Cas9-mediated systems that bias inheritance of desired alleles in mice and that have the potential to transform the use of rodent models in basic and biomedical research.

During meiotic recombination, homologue-templated repair of programmed DNA double-strand breaks (DSBs) produces relatively few crossovers and many difficult-to-detect non-crossovers. By intercrossing two diverged mouse subspecies over five generations and deep-sequencing 119 offspring, we detect thousands of crossover and non-crossover events genome-wide with unprecedented power and spatial resolution. We find that both crossovers and non-crossovers are strongly depleted at DSB hotspots where the DSB-positioning protein PRDM9 fails to bind to the unbroken homologous chromosome, revealing that PRDM9 also functions to promote homologue-templated repair. Our results show that complex non-crossovers are much rarer in mice than humans, consistent with complex events arising from accumulated non-programmed DNA damage. Unexpectedly, we also find that GC-biased gene conversion is restricted to non-crossover tracts containing only one mismatch. These results demonstrate that local genetic diversity profoundly alters meiotic repair pathway decisions via at least two distinct mechanisms, impacting genome evolution and Prdm9-related hybrid infertility.

Understanding the plasticity of genomes has been greatly aided by assays for recombination, repair and mutagenesis. These assays have been developed in microbial systems that provide the advantages of genetic and molecular reporters that can readily be manipulated. Cellular assays comprise genetic, molecular, and cytological reporters. The assays are powerful tools but each comes with its particular advantages and limitations. Here the most commonly used assays are reviewed, discussed, and presented as the guidelines for future studies.

Double strand DNA breaks (DSBs) are dangerous events that can result from various causes including environmental assaults or the collapse of DNA replication. While the efficient and precise repair of DSBs is essential for cell survival, faulty repair can lead to genetic instability, making the choice of DSB repair an important step. Here we report that inverted DNA repeats (IRs) placed near a DSB can channel its repair from an accurate pathway that leads to gene conversion to instead a break-induced replication (BIR) pathway that leads to genetic instabilities. The effect of IRs is explained by their ability to form unusual DNA structures when present in ssDNA that is formed by DSB resection. We demonstrate that IRs can form two types of unusual DNA structures, and the choice between these structures depends on the length of the spacer separating IRs. In particular, IRs separated by a long (1-kb) spacer are predominantly involved in inter-molecular single-strand annealing (SSA) leading to the formation of inverted dimers; IRs separated by a short (12-bp) spacer participate in intra-molecular SSA, leading to the formation of fold-back (FB) structures. Both of these structures interfere with an accurate DSB repair by gene conversion and channel DSB repair into BIR, which promotes genomic destabilization. We also report that different protein complexes participate in the processing of FBs containing short (12-bp) versus long (1-kb) ssDNA loops. Specifically, FBs with short loops are processed by the MRX-Sae2 complex, whereas the Rad1-Rad10 complex is responsible for the processing of long loops. Overall, our studies uncover the mechanisms of genomic destabilization resulting from re-routing DSB repair into unusual pathways by IRs. Given the high abundance of IRs in the human genome, our findings may contribute to the understanding of IR-mediated genomic destabilization associated with human disease.

Some of the fastest evolving regions of the human genome are conserved noncoding elements with many human-specific DNA substitutions. These human accelerated regions (HARs) are enriched nearby regulatory genes, and several HARs function as developmental enhancers. To investigate if this evolutionary signature is unique to humans, we quantified evidence of accelerated substitutions in conserved genomic elements across multiple lineages and applied this approach simultaneously to the genomes of five apes: human, chimpanzee, gorilla, orangutan, and gibbon. We find roughly similar numbers and genomic distributions of lineage-specific accelerated regions (linARs) in all five apes. In particular, apes share an enrichment of linARs in regulatory DNA nearby genes involved in development, especially transcription factors and other regulators. Many developmental loci harbor clusters of nonoverlapping linARs from multiple apes, suggesting that accelerated evolution in each species affected distinct regulatory elements that control a shared set of developmental pathways. Our statistical tests distinguish between GC-biased and unbiased accelerated substitution rates, allowing us to quantify the roles of different evolutionary forces in creating linARs. We find evidence of GC-biased gene conversion in each ape, but unbiased acceleration consistent with positive selection or loss of constraint is more common in all five lineages. It therefore appears that similar evolutionary processes created independent accelerated regions in the genomes of different apes, and that these lineage-specific changes to conserved noncoding sequences may have differentially altered expression of a core set of developmental genes across ape evolution.

The use of engineered nucleases combined with a homologous DNA donor template can result in targeted gene correction of the sickle cell disease mutation in hematopoietic stem and progenitor cells. However, because of the high homology between the adjacent human β- and δ-globin genes, off-target cleavage is observed at δ-globin when using some endonucleases targeted to the sickle mutation in β-globin. Introduction of multiple double-stranded breaks by endonucleases has the potential to induce intergenic alterations. Using a novel droplet digital PCR assay and high-throughput sequencing, we characterized the frequency of rearrangements between the β- and δ-globin paralogs when delivering these nucleases. Pooled CD34+ cells and colony-forming units from sickle bone marrow were treated with nuclease only or including a donor template and then analyzed for potential gene rearrangements. It was observed that, in pooled CD34+ cells and colony-forming units, the intergenic β-δ-globin deletion was the most frequent rearrangement, followed by inversion of the intergenic fragment, with the inter-chromosomal translocation as the least frequent. No rearrangements were observed when endonuclease activity was restricted to on-target β-globin cleavage. These findings demonstrate the need to develop site-specific endonucleases with high specificity to avoid unwanted gene alterations.

BackgroundGene innovation by duplication is a fundamental evolutionary process but is difficult to study in humans due to the large size, high sequence identity, and mosaic nature of segmental duplication blocks. The human-specific gene hydrocephalus-inducing 2, HYDIN2, was generated by a 364 kbp duplication of 79 internal exons of the large ciliary gene HYDIN from chromosome 16q22.2 to chromosome 1q21.1. Because the HYDIN2 locus lacks the ancestral promoter and seven terminal exons of the progenitor gene, we sought to characterize transcription at this locus by coupling reverse transcription polymerase chain reaction and long-read sequencing.Results5' RACE indicates a transcription start site for HYDIN2 outside of the duplication and we observe fusion transcripts spanning both the 5' and 3' breakpoints. We observe extensive splicing diversity leading to the formation of altered open reading frames (ORFs) that appear to be under relaxed selection. We show that HYDIN2 adopted a new promoter that drives an altered pattern of expression, with highest levels in neural tissues. We estimate that the HYDIN duplication occurred ~3.2 million years ago and find that it is nearly fixed (99.9%) for diploid copy number in contemporary humans. Examination of 73 chromosome 1q21 rearrangement patients reveals that HYDIN2 is deleted or duplicated in most cases.ConclusionsTogether, these data support a model of rapid gene innovation by fusion of incomplete segmental duplications, altered tissue expression, and potential subfunctionalization or neofunctionalization of HYDIN2 early in the evolution of the Homo lineage.

BackgroundSiglecs-11 and -16 are members of the sialic acid recognizing Ig-like lectin family, and expressed in same cells. Siglec-11 functions as an inhibitory receptor, whereas Siglec-16 exhibits activating properties. In humans, SIGLEC11 and SIGLEC16 gene sequences are extremely similar in the region encoding the extracellular domain due to gene conversions. Human SIGLEC11 was converted by the nonfunctional SIGLEC16P allele, and the converted SIGLEC11 allele became fixed in humans, possibly because it provides novel neuroprotective functions in brain microglia. However, the detailed evolutionary history of SIGLEC11 and SIGLEC16 in other primates remains unclear.ResultsWe analyzed SIGLEC11 and SIGLEC16 gene sequences of multiple primate species, and examined glycan binding profiles of these Siglecs. The phylogenetic tree demonstrated that gene conversions between SIGLEC11 and SIGLEC16 occurred in the region including the exon encoding the sialic acid binding domain in every primate examined. Functional assays showed that glycan binding preference is similar between Siglec-11 and Siglec-16 in all analyzed hominid species. Taken together with the fact that Siglec-11 and Siglec-16 are expressed in the same cells, Siglec-11 and Siglec-16 are regarded as paired receptors that have maintained similar ligand binding preferences via gene conversions. Relaxed functional constraints were detected on the SIGLEC11 and SIGLEC16 exons that underwent gene conversions, possibly contributing to the evolutionary acceptance of repeated gene conversions. The frequency of nonfunctional SIGLEC16P alleles is much higher than that of SIGLEC16 alleles in every human population.ConclusionsOur findings indicate that Siglec-11 and Siglec-16 have been maintained as paired receptors by repeated gene conversions under relaxed functional constraints in the primate lineage. The high prevalence of the nonfunctional SIGLEC16P allele and the fixation of the converted SIGLEC11 imply that the loss of Siglec-16 and the gain of Siglec-11 in microglia might have been favored during the evolution of human lineage.

Gene clusters of recently duplicated genes are hotbeds for evolutionary change. However, our understanding of how mutational mechanisms and evolutionary forces shape the structural and functional evolution of these clusters is hindered by the high sequence identity among the copies, which typically results in their inaccurate representation in genome assemblies. The presumed testis-specific, chimeric gene Sdic originated, and tandemly expanded in Drosophila melanogaster, contributing to increased male-male competition. Using various types of massively parallel sequencing data, we studied the organization, sequence evolution, and functional attributes of the different Sdic copies. By leveraging long-read sequencing data, we uncovered both copy number and order differences from the currently accepted annotation for the Sdic region. Despite evidence for pervasive gene conversion affecting the Sdic copies, we also detected signatures of two episodes of diversifying selection, which have contributed to the evolution of a variety of C-termini and miRNA binding site compositions. Expression analyses involving RNA-seq datasets from 59 different biological conditions revealed distinctive expression breadths among the copies, with three copies being transcribed in females, opening the possibility to a sexually antagonistic effect. Phenotypic assays using Sdic knock-out strains indicated that should this antagonistic effect exist, it does not compromise female fertility. Our results strongly suggest that the genome consolidation of the Sdic gene cluster is more the result of a quick exploration of different paths of molecular tinkering by different copies than a mere dosage increase, which could be a recurrent evolutionary outcome in the presence of persistent sexual selection.

In yeast, DNA breaks are usually repaired by homologous recombination (HR). An early step for HR pathways is formation of a heteroduplex, in which a single-strand from the broken DNA molecule pairs with a strand derived from an intact DNA molecule. If the two strands of DNA are not identical, there will be mismatches within the heteroduplex DNA (hetDNA). In wild-type strains, these mismatches are repaired by the mismatch repair (MMR) system, producing a gene conversion event. In strains lacking MMR, the mismatches persist. Most previous studies involving hetDNA formed during mitotic recombination were restricted to one locus. Below, we present a global mapping of hetDNA formed in the MMR-defective mlh1 strain. We find that many recombination events are associated with repair of double-stranded DNA gaps and/or involve Mlh1-independent mismatch repair. Many of our events are not explicable by the simplest form of the double-strand break repair model of recombination.

Events in primate evolution are often dated by assuming a constant rate of substitution per unit time, but the validity of this assumption remains unclear. Among mammals, it is well known that there exists substantial variation in yearly substitution rates. Such variation is to be expected from differences in life history traits, suggesting it should also be found among primates. Motivated by these considerations, we analyze whole genomes from 10 primate species, including Old World Monkeys (OWMs), New World Monkeys (NWMs), and apes, focusing on putatively neutral autosomal sites and controlling for possible effects of biased gene conversion and methylation at CpG sites. We find that substitution rates are up to 64% higher in lineages leading from the hominoid-NWM ancestor to NWMs than to apes. Within apes, rates are ∼2% higher in chimpanzees and ∼7% higher in the gorilla than in humans. Substitution types subject to biased gene conversion show no more variation among species than those not subject to it. Not all mutation types behave similarly, however; in particular, transitions at CpG sites exhibit a more clocklike behavior than do other types, presumably because of their nonreplicative origin. Thus, not only the total rate, but also the mutational spectrum, varies among primates. This finding suggests that events in primate evolution are most reliably dated using CpG transitions. Taking this approach, we estimate the human and chimpanzee divergence time is 12.1 million years,​ and the human and gorilla divergence time is 15.1 million years​.

Plant chloroplast genomes (plastomes) are characterized by an inverted repeat (IR) region and two larger single copy (SC) regions. Patterns of molecular evolution in the IR and SC regions differ, most notably by a reduced rate of nucleotide substitution in the IR compared to the SC region. In addition, the organization and structure of plastomes is fluid, and rearrangements through time have repeatedly shuffled genes into and out of the IR, providing recurrent natural experiments on how chloroplast genome structure can impact rates and patterns of molecular evolution. Here we examine four loci (psbA, ycf2, rps7, and rps12 exon 2-3) that were translocated from the SC into the IR during fern evolution. We use a model-based method, within a phylogenetic context, to test for substitution rate shifts. All four loci show a significant, 2- to 3-fold deceleration in their substitution rate following translocation into the IR, a phenomenon not observed in any other, nontranslocated plastid genes. Also, we show that after translocation, the GC content of the third codon position and of the noncoding regions is significantly increased, implying that gene conversion within the IR is GC-biased. Taken together, our results suggest that the IR region not only reduces substitution rates, but also impacts nucleotide composition. This finding highlights a potential vulnerability of correlating substitution rate heterogeneity with organismal life history traits without knowledge of the underlying genome structure.

Our understanding of the chronology of human evolution relies on the "molecular clock" provided by the steady accumulation of substitutions on an evolutionary lineage. Recent analyses of human pedigrees have called this understanding into question by revealing unexpectedly low germline mutation rates, which imply that substitutions accrue more slowly than previously believed. Translating mutation rates estimated from pedigrees into substitution rates is not as straightforward as it may seem, however. We dissect the steps involved, emphasizing that dating evolutionary events requires not "a mutation rate" but a precise characterization of how mutations accumulate in development in males and females-knowledge that remains elusive.

IFIT (interferon-induced with tetratricopeptide repeats) proteins are critical mediators of mammalian innate antiviral immunity. Mouse IFIT1 selectively inhibits viruses that lack 2'O-methylation of their mRNA 5' caps. Surprisingly, human IFIT1 does not share this antiviral specificity. Here, we resolve this discrepancy by demonstrating that human and mouse IFIT1 have evolved distinct functions using a combination of evolutionary, genetic and virological analyses. First, we show that human IFIT1 and mouse IFIT1 (renamed IFIT1B) are not orthologs, but are paralogs that diverged >100 mya. Second, using a yeast genetic assay, we show that IFIT1 and IFIT1B proteins differ in their ability to be suppressed by a cap 2'O-methyltransferase. Finally, we demonstrate that IFIT1 and IFIT1B have divergent antiviral specificities, including the discovery that only IFIT1 proteins inhibit a virus encoding a cap 2'O-methyltransferase. These functional data, combined with widespread turnover of mammalian IFIT genes, reveal dramatic species-specific differences in IFIT-mediated antiviral repertoires.

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

Abstract Background Sex chromosomes have arisen independently in a wide variety of species, yet they share common characteristics, including the presence of suppressed recombination surrounding sex determination loci. Mammalian sex chromosomes contain multiple palindromic repeats across the non-recombining region that show sequence conservation through gene conversion and contain genes that are crucial for sexual reproduction. In plants, it is not clear if palindromic repeats play a role in maintaining sequence conservation in the absence of homologous recombination. Results Here we present the first evidence of large palindromic structures in a plant sex chromosome, based on a highly contiguous assembly of the W chromosome of the dioecious shrub Salix purpurea. The W chromosome has an expanded number of genes due to transpositions from autosomes. It also contains two consecutive palindromes that span a region of 200 kb, with conspicuous 20-kb stretches of highly conserved sequences among the four arms that show evidence of gene conversion. Four genes in the palindrome are homologous to genes in the sex determination regions of the closely related genus Populus, which is located on a different chromosome. These genes show distinct, floral-biased expression patterns compared to paralogous copies on autosomes. Conclusion The presence of palindromes in sex chromosomes of mammals and plants highlights the intrinsic importance of these features in adaptive evolution in the absence of recombination. Convergent evolution is driving both the independent establishment of sex chromosomes as well as their fine-scale sequence structure.

As in other organisms, CRISPR/Cas9 methods provide a powerful approach for genome editing in the nematode Caenorhabditis elegans. Oligonucleotides are excellent repair templates for introducing substitutions and short insertions, as they are cost effective, require no cloning, and appear in other organisms to target changes by homologous recombination at DNA double-strand breaks (DSBs). Here, I describe a methodology in C. elegans to efficiently knock in epitope tags in 8-9 days, using a temperature-sensitive lethal mutation in the pha-1 gene as a co-conversion marker. I demonstrate that 60mer oligos with 29 bp of homology drive efficient knock-in of point mutations, and that disabling nonhomologous end joining by RNAi inactivation of the cku-80 gene significantly improves knock-in efficiency. Homology arms of 35-80 bp are sufficient for efficient editing and DSBs up to 54 bp away from the insertion site produced knock-ins. These findings will likely be applicable for a range of genome editing approaches in C. elegans, which will improve editing efficiency and minimize screening efforts.

Facilitated by recent advances using CRISPR/Cas9, genome editing technologies now permit custom genetic modifications in a wide variety of organisms. Ideally, modified animals could be both efficiently made and easily identified with minimal initial screening and without introducing exogenous sequence at the locus of interest or marker mutations elsewhere. To this end, we describe a coconversion strategy, using CRISPR/Cas9 in which screening for a dominant phenotypic oligonucleotide-templated conversion event at one locus can be used to enrich for custom modifications at another unlinked locus. After the desired mutation is identified among the F1 progeny heterozygous for the dominant marker mutation, F2 animals that have lost the marker mutation are picked to obtain the desired mutation in an unmarked genetic background. We have developed such a coconversion strategy for Caenorhabditis elegans, using a number of dominant phenotypic markers. Examining the coconversion at a second (unselected) locus of interest in the marked F1 animals, we observed that 14-84% of screened animals showed homologous recombination. By reconstituting the unmarked background through segregation of the dominant marker mutation at each step, we show that custom modification events can be carried out recursively, enabling multiple mutant animals to be made. While our initial choice of a coconversion marker [rol-6(su1006)] was readily applicable in a single round of coconversion, the genetic properties of this locus were not optimal in that CRISPR-mediated deletion mutations at the unselected rol-6 locus can render a fraction of coconverted strains recalcitrant to further rounds of similar mutagenesis. An optimal marker in this sense would provide phenotypic distinctions between the desired mutant/+ class and alternative +/+, mutant/null, null/null, and null/+ genotypes. Reviewing dominant alleles from classical C. elegans genetics, we identified one mutation in dpy-10 and one mutation in sqt-1 that meet these criteria and demonstrate that these too can be used as effective conversion markers. Coconversion was observed using a variety of donor molecules at the second (unselected) locus, including oligonucleotides, PCR products, and plasmids. We note that the coconversion approach described here could be applied in any of the variety of systems where suitable coconversion markers can be identified from previous intensive genetic analyses of gain-of-function alleles.

The yeast Exo1p nuclease functions in multiple cellular roles: resection of DNA ends generated during recombination, telomere stability, DNA mismatch repair, and expansion of gaps formed during the repair of UV-induced DNA damage. In this study, we performed high-resolution mapping of spontaneous and UV-induced recombination events between homologs in exo1 strains, comparing the results with spontaneous and UV-induced recombination events in wild-type strains. One important comparison was the lengths of gene conversion tracts. Gene conversion events are usually interpreted as reflecting heteroduplex formation between interacting DNA molecules, followed by repair of mismatches within the heteroduplex. In most models of recombination, the length of the gene conversion tract is a function of the length of single-stranded DNA generated by end resection. Since the Exo1p has an important role in end resection, a reduction in the lengths of gene conversion tracts in exo1 strains was expected. In accordance with this expectation, gene conversion tract lengths associated with spontaneous crossovers in exo1 strains were reduced about twofold relative to wild type. For UV-induced events, conversion tract lengths associated with crossovers were also shorter for the exo1 strain than for the wild-type strain (3.2 and 7.6 kb, respectively). Unexpectedly, however, the lengths of conversion tracts that were unassociated with crossovers were longer in the exo1 strain than in the wild-type strain (6.2 and 4.8 kb, respectively). Alternative models of recombination in which the lengths of conversion tracts are determined by break-induced replication or oversynthesis during strand invasion are proposed to account for these observations.

Abstract Background Sequence exchange between homologous chromosomes through crossing over and gene conversion is highly conserved among eukaryotes, contributing to genome stability and genetic diversity. A lack of recombination limits breeding efforts in crops; therefore, increasing recombination rates can reduce linkage drag and generate new genetic combinations. Results We use computational analysis of 13 recombinant inbred mapping populations to assess crossover and gene conversion frequency in the hexaploid genome of wheat (Triticum aestivum). We observe that high-frequency crossover sites are shared between populations and that closely related parents lead to populations with more similar crossover patterns. We demonstrate that gene conversion is more prevalent and covers more of the genome in wheat than in other plants, making it a critical process in the generation of new haplotypes, particularly in centromeric regions where crossovers are rare. We identify quantitative trait loci for altered gene conversion and crossover frequency and confirm functionality for a novel RecQ helicase gene that belongs to an ancient clade that is missing in some plant lineages including Arabidopsis. Conclusions This is the first gene to be demonstrated to be involved in gene conversion in wheat. Harnessing the RecQ helicase has the potential to break linkage drag utilizing widespread gene conversions.

Abstract Background The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention. Results We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events—transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database. Conclusions We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

In the yeast Saccharomyces cerevisiae and most other eukaryotes, mitotic recombination is important for the repair of double-stranded DNA breaks (DSBs). Mitotic recombination between homologous chromosomes can result in loss of heterozygosity (LOH). In this study, LOH events induced by ultraviolet (UV) light are mapped throughout the genome to a resolution of about 1 kb using single-nucleotide polymorphism (SNP) microarrays. UV doses that have little effect on the viability of diploid cells stimulate crossovers more than 1000-fold in wild-type cells. In addition, UV stimulates recombination in G1-synchronized cells about 10-fold more efficiently than in G2-synchronized cells. Importantly, at high doses of UV, most conversion events reflect the repair of two sister chromatids that are broken at approximately the same position whereas at low doses, most conversion events reflect the repair of a single broken chromatid. Genome-wide mapping of about 380 unselected crossovers, break-induced replication (BIR) events, and gene conversions shows that UV-induced recombination events occur throughout the genome without pronounced hotspots, although the ribosomal RNA gene cluster has a significantly lower frequency of crossovers.

GC-biased gene conversion (gBGC) is a recombination-associated process that favors the fixation of G/C alleles over A/T alleles. In mammals, gBGC is hypothesized to contribute to variation in GC content, rapidly evolving sequences, and the fixation of deleterious mutations, but its prevalence and general functional consequences remain poorly understood. gBGC is difficult to incorporate into models of molecular evolution and so far has primarily been studied using summary statistics from genomic comparisons. Here, we introduce a new probabilistic model that captures the joint effects of natural selection and gBGC on nucleotide substitution patterns, while allowing for correlations along the genome in these effects. We implemented our model in a computer program, called phastBias, that can accurately detect gBGC tracts about 1 kilobase or longer in simulated sequence alignments. When applied to real primate genome sequences, phastBias predicts gBGC tracts that cover roughly 0.3% of the human and chimpanzee genomes and account for 1.2% of human-chimpanzee nucleotide differences. These tracts fall in clusters, particularly in subtelomeric regions; they are enriched for recombination hotspots and fast-evolving sequences; and they display an ongoing fixation preference for G and C alleles. They are also significantly enriched for disease-associated polymorphisms, suggesting that they contribute to the fixation of deleterious alleles. The gBGC tracts provide a unique window into historical recombination processes along the human and chimpanzee lineages. They supply additional evidence of long-term conservation of megabase-scale recombination rates accompanied by rapid turnover of hotspots. Together, these findings shed new light on the evolutionary, functional, and disease implications of gBGC. The phastBias program and our predicted tracts are freely available.

During meiosis, the recombination-initiating DNA double-strand breaks (DSBs) are repaired by crossovers or noncrossovers (gene conversions). While crossovers are easily detectable, noncrossover identification is hampered by the small size of their converted tracts and the necessity of sequence polymorphism. We report identification and characterization of a mouse chromosome-wide set of noncrossovers by next-generation sequencing of 10 mouse intersubspecific chromosome substitution strains. Based on 94 identified noncrossovers, we determined the mean length of a conversion tract to be 32 bp. The spatial chromosome-wide distribution of noncrossovers and crossovers significantly differed, although both sets overlapped the known hotspots of PRDM9-directed histone methylation and DNA DSBs, thus supporting their origin in the standard DSB repair pathway. A significant deficit of noncrossovers descending from asymmetric DSBs proved their proposed adverse effect on meiotic recombination and pointed to sister chromatids as an alternative template for their repair. The finding has implications for the molecular mechanism of hybrid sterility in mice from crosses between closely related Mus musculus musculus and Mus musculus domesticus subspecies. [ABSTRACT FROM AUTHOR]

The fastest-evolving regions in the human and chimpanzee genomes show a remarkable excess of weak (A,T) to strong (G,C) nucleotide substitutions since divergence from their common ancestor. We investigated the phylogenetic extent and possible causes of this weak to strong (W → S) bias in divergent sequences (BDS) using recently sequenced genomes and recombination maps from eight trios of eukaryotic species. To quantify evidence for BDS, we inferred substitution histories using an efficient maximum likelihood approach with a context-dependent evolutionary model. We then annotated all lineage-specific substitutions in terms of W → S bias and density on the chromosomes. Finally, we used the inferred substitutions to calculate a BDS score-a log odds ratio between substitution type and density-and assessed its statistical significance with Fisher's exact test. Applying this approach, we found significant BDS in the coding and noncoding sequence of human, mouse, dog, stickleback, fruit fly, and worm. We also observed a significant lack of W → S BDS in chicken and yeast. The BDS score varies between species and across the chromosomes within each species. It is most strongly correlated with different genomic features in different species, but a strong correlation with recombination rates is found in several species. Our results demonstrate that a W → S substitution bias in fast-evolving sequences is a widespread phenomenon. The patterns of BDS observed suggest that a recombination-associated process, such as GC-biased gene conversion, is involved in the production of the bias in many species, but the strength of the BDS likely depends on many factors, including genome stability, variability in recombination rate over time and across the genome, the frequency of meiosis, and the amount of outcrossing in each species.

Abstract Background Meiotic recombination events include crossovers and non-crossovers or gene conversions. Although the rate of crossovers is often used for genetic mapping, the gene conversion events are not well studied especially in outbred species, which could produce distorted markers and thus affect the precision of genetic maps. Results We proposed a strategy for identifying gene conversion events in Populus with the next-generation sequencing (NGS) data from the two parents and their progeny in an F1 hybrid population. The strategy first involved phasing the heterozygous SNPs of the parents to obtain the parental haplotype blocks by NGS analytical tools, permitting to identify the parental gene conversion events with progeny genotypes. By incorporating available genetic linkage maps, longer haplotype blocks each corresponding to a chromosome can be created, not only allowing to detect crossover events but also possibly to locate a crossover in a small region. Our analysis revealed that gene conversions are more abundant than crossovers in Populus, with a higher probability to generate distorted markers in the regions involved than in the other regions on genome. The analytical procedures were implemented with Perl scripts as a freely available package, findGCO at https://github.com/tongchf/findGCO. Conclusions The novel strategy and the new developed Perl package permit to identify gene conversion events with the next-generation sequencing technology in a hybrid population of outbred species. The new method revealed that in a genetic mapping population some distorted genetic markers are possibly due to the gene conversion events.

Centromeres are the most dynamic regions of the genome, yet they are typified by little or no crossing over, making it difficult to explain the origin of this diversity. To address this question, we developed a novel CENH3 ChIP display method that maps kinetochore footprints over transposon-rich areas of centromere cores. A high level of polymorphism made it possible to map a total of 238 within-centromere markers using maize recombinant inbred lines. Over half of the markers were shown to interact directly with kinetochores (CENH3) by chromatin immunoprecipitation. Although classical crossing over is fully suppressed across CENH3 domains, two gene conversion events (i.e., non-crossover marker exchanges) were identified in a mapping population. A population genetic analysis of 53 diverse inbreds suggests that historical gene conversion is widespread in maize centromeres, occurring at a rate >1x10(-5)/marker/generation. We conclude that gene conversion accelerates centromere evolution by facilitating sequence exchange among chromosomes.