Increased mutation and gene conversion within human segmental duplications.

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.

Disclosing complex mutational dynamics at a Y chromosome palindrome evolving through intra- and inter-chromosomal gene conversion.

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.

Applications of and considerations for using CRISPR-Cas9-mediated gene conversion systems in rodents.

Genetic elements that are inherited at super-Mendelian frequencies could be used in a 'gene drive' to spread an allele to high prevalence in a population with the goal of eliminating invasive species or disease vectors. We recently demonstrated that the gene conversion mechanism underlying a CRISPR-Cas9-mediated gene drive is feasible in mice. Although substantial technical hurdles remain, overcoming these could lead to strategies that might decrease the spread of rodent-borne Lyme disease or eliminate invasive populations of mice and rats that devastate island ecology. Perhaps more immediately achievable at moderate gene conversion efficiency, applications in a laboratory setting could produce complex genotypes that reduce the time and cost in both dollars and animal lives compared with Mendelian inheritance strategies. Here, we discuss what we have learned from early efforts to achieve CRISPR-Cas9-mediated gene conversion, potential for broader applications in the laboratory, current limitations, and plans for optimizing this potentially powerful technology.

Meiotic Cas9 expression mediates gene conversion in the male and female mouse germline.

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.

Genomic GC content drifts downward in most bacterial genomes.

In every kingdom of life, GC->AT transitions occur more frequently than any other type of mutation due to the spontaneous deamination of cytidine. In eukaryotic genomes, this slow loss of GC base pairs is counteracted by biased gene conversion which increases genomic GC content as part of the recombination process. However, this type of biased gene conversion has not been observed in bacterial genomes, so we hypothesized that GC->AT transitions cause a reduction of genomic GC content in prokaryotic genomes on an evolutionary time scale. To test this hypothesis, we used a phylogenetic approach to analyze triplets of closely related genomes representing a wide range of the bacterial kingdom. The resulting data indicate that genomic GC content is drifting downward in bacterial genomes where GC base pairs comprise 40% or more of the total genome. In contrast, genomes containing less than 40% GC base pairs have fewer opportunities for GC->AT transitions to occur so genomic GC content is relatively stable or actually increasing. It should be noted that this observed change in genomic GC content is the net change in shared parts of the genome and does not apply to parts of the genome that have been lost or acquired since the genomes being compared shared common ancestor. However, a more detailed analysis of two Caulobacter genomes revealed that the acquisition of mobile elements by the two genomes actually reduced the total genomic GC content as well.

Mitotic gene conversion can be as important as meiotic conversion in driving genetic variability in plants and other species without early germline segregation.

In contrast to common meiotic gene conversion, mitotic gene conversion, because it is so rare, is often ignored as a process influencing allelic diversity. We show that if there is a large enough number of premeiotic cell divisions, as seen in many organisms without early germline sequestration, such as plants, this is an unsafe position. From examination of 1.1 million rice plants, we determined that the rate of mitotic gene conversion events, per mitosis, is 2 orders of magnitude lower than the meiotic rate. However, owing to the large number of mitoses between zygote and gamete and because of long mitotic tract lengths, meiotic and mitotic gene conversion can be of approximately equivalent importance in terms of numbers of markers converted from zygote to gamete. This holds even if we assume a low number of premeiotic cell divisions (approximately 40) as witnessed in Arabidopsis. A low mitotic rate associated with long tracts is also seen in yeast, suggesting generality of results. For species with many mitoses between each meiotic event, mitotic gene conversion should not be overlooked.

Detection of SMN1 to SMN2 gene conversion events and partial SMN1 gene deletions using array digital PCR.

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.

Genome-wide mapping of spontaneous genetic alterations in diploid yeast cells.

Genomic alterations including single-base mutations, deletions and duplications, translocations, mitotic recombination events, and chromosome aneuploidy generate genetic diversity. We examined the rates of all of these genetic changes in a diploid strain of Saccharomyces cerevisiae by whole-genome sequencing of many independent isolates (n = 93) subcloned about 100 times in unstressed growth conditions. The most common alterations were point mutations and small (<100 bp) insertion/deletions (n = 1,337) and mitotic recombination events (n = 1,215). The diploid cells of most eukaryotes are heterozygous for many single-nucleotide polymorphisms (SNPs). During mitotic cell divisions, recombination can produce derivatives of these cells that have become homozygous for the polymorphisms, termed loss-of-heterozygosity (LOH) events. LOH events can change the phenotype of the cells and contribute to tumor formation in humans. We observed two types of LOH events: interstitial events (conversions) resulting in a short LOH tract (usually less than 15 kb) and terminal events (mostly cross-overs) in which the LOH tract extends to the end of the chromosome. These two types of LOH events had different distributions, suggesting that they may have initiated by different mechanisms. Based on our results, we present a method of calculating the probability of an LOH event for individual SNPs located throughout the genome. We also identified several hotspots for chromosomal rearrangements (large deletions and duplications). Our results provide insights into the relative importance of different types of genetic alterations produced during vegetative growth.

Resection and repair of a Cas9 double-strand break at CTG trinucleotide repeats induces local and extensive chromosomal deletions.

Microsatellites are short tandem repeats, ubiquitous in all eukaryotes and represent ~2% of the human genome. Among them, trinucleotide repeats are responsible for more than two dozen neurological and developmental disorders. Targeting microsatellites with dedicated DNA endonucleases could become a viable option for patients affected with dramatic neurodegenerative disorders. Here, we used the Streptococcus pyogenes Cas9 to induce a double-strand break within the expanded CTG repeat involved in myotonic dystrophy type 1, integrated in a yeast chromosome. Repair of this double-strand break generated unexpected large chromosomal deletions around the repeat tract. These deletions depended on RAD50, RAD52, DNL4 and SAE2, and both non-homologous end-joining and single-strand annealing pathways were involved. Resection and repair of the double-strand break (DSB) were totally abolished in a rad50Δ strain, whereas they were impaired in a sae2Δ mutant, only on the DSB end containing most of the repeat tract. This observation demonstrates that Sae2 plays significant different roles in resecting a DSB end containing a repeated and structured sequence as compared to a non-repeated DSB end. In addition, we also discovered that gene conversion was less efficient when the DSB could be repaired using a homologous template, suggesting that the trinucleotide repeat may interfere with gene conversion too. Altogether, these data show that SpCas9 may not be the best choice when inducing a double-strand break at or near a microsatellite, especially in mammalian genomes that contain many more dispersed repeated elements than the yeast genome.

Adapting Biased Gene Conversion theory to account for intensive GC-content deterioration in the human genome by novel mutations.

We examined seventy million well-characterized human mutations, and their impact on G+C-compositional dynamics, in order to understand the formation and maintenance of major genomic nucleotide sequence patterns. Among novel mutations, those that change a strong (S) base pair G:C/C:G to a weak (W) pair A:T/T:A occur at nearly twice the frequency of the opposite mutations. Such imbalance puts strong downward pressure on overall GC-content. However, along protracted paths to fixation, S→W mutations are much less likely to propagate than W→S mutations. The magnitude of relative propagation disadvantages for S→W mutations is inexplicable by any currently-accepted model. This fact forced us to re-examine the quantitative features of Biased Gene Conversion (BGC) theory. Revised parameters of BGC that, per average individual, convert 7-14 W base pairs into S pairs, would account for the S-content turnover differences between new and old mutations, and make BGC an instrumental force for nucleotide dynamics and evolution. BGC should thus be considered seriously in both theories and biomedical practice. In particular, BGC should be taken into account during allele imputations, where missing SNP alleles are computationally predicted based on the information about several neighboring alleles. Finally, we analyzed the effect of neighboring nucleotide context on the mutation frequencies, dynamics, and GC-composition turnover. For this purpose, we examined genomic regions having extremely biased nucleotide compositions (enriched for S-, W-, purine/pyrimidine strand asymmetry, or AC/GT-strand asymmetry). It was found that point mutations in these regions preferentially degrade the nucleotide inhomogeneities, decreasing the sequence biases. Degradation of sequence bias is highest for novel mutations, and considerably lower for older mutations (those widespread across populations). Besides BGC, there may be additional, still uncharacterized molecular mechanisms that either preserve genomic regions with biased nucleotide compositions from mutational degradation or fail to degrade such inhomogeneities in specific chromosomal regions.

Strong suppression of gene conversion with increasing DNA double-strand break load delimited by 53BP1 and RAD52.

In vertebrates, genomic DNA double-strand breaks (DSBs) are removed by non-homologous end-joining processes: classical non-homologous end-joining (c-NHEJ) and alternative end-joining (alt-EJ); or by homology-dependent processes: gene-conversion (GC) and single-strand annealing (SSA). Surprisingly, these repair pathways are not real alternative options restoring genome integrity with equal efficiency, but show instead striking differences in speed, accuracy and cell-cycle-phase dependence. As a consequence, engagement of one pathway may be associated with processing-risks for the genome absent from another pathway. Characterization of engagement-parameters and their consequences is, therefore, essential for understanding effects on the genome of DSB-inducing agents, such as ionizing-radiation (IR). Here, by addressing pathway selection in G2-phase, we discover regulatory confinements in GC with consequences for SSA- and c-NHEJ-engagement. We show pronounced suppression of GC with increasing DSB-load that is not due to RAD51 availability and which is delimited but not defined by 53BP1 and RAD52. Strikingly, at low DSB-loads, GC repairs ∼50% of DSBs, whereas at high DSB-loads its contribution is undetectable. Notably, with increasing DSB-load and the associated suppression of GC, SSA gains ground, while alt-EJ is suppressed. These observations explain earlier, apparently contradictory results and advance our understanding of logic and mechanisms underpinning the wiring between DSB repair pathways.

Linking high GC content to the repair of double strand breaks in prokaryotic genomes.

Genomic GC content varies widely among microbes for reasons unknown. While mutation bias partially explains this variation, prokaryotes near-universally have a higher GC content than predicted solely by this bias. Debate surrounds the relative importance of the remaining explanations of selection versus biased gene conversion favoring GC alleles. Some environments (e.g. soils) are associated with a high genomic GC content of their inhabitants, which implies that either high GC content is a selective adaptation to particular habitats, or that certain habitats favor increased rates of gene conversion. Here, we report a novel association between the presence of the non-homologous end joining DNA double-strand break repair pathway and GC content; this observation suggests that DNA damage may be a fundamental driver of GC content, leading in part to the many environmental patterns observed to-date. We discuss potential mechanisms accounting for the observed association, and provide preliminary evidence that sites experiencing higher rates of double-strand breaks are under selection for increased GC content relative to the genomic background.

Adaptation by Loss of Heterozygosity in Saccharomyces cerevisiae Clones Under Divergent Selection.

Loss of heterozygosity (LOH) is observed during vegetative growth and reproduction of diploid genotypes through mitotic crossovers, aneuploidy caused by nondisjunction, and gene conversion. We aimed to test the role that LOH plays during adaptation of two highly heterozygous Saccharomyces cerevisiae genotypes to multiple environments over a short time span in the laboratory. We hypothesized that adaptation would be observed through parallel LOH events across replicate populations. Using genome resequencing of 70 clones, we found that LOH was widespread with 5.2 LOH events per clone after ∼500 generations. The most common mode of LOH was gene conversion (51%) followed by crossing over consistent with either break-induced replication or double Holliday junction resolution. There was no evidence that LOH involved nondisjunction of whole chromosomes. We observed parallel LOH in both an environment-specific and environment-independent manner. LOH largely involved recombining existing variation between the parental genotypes, but also was observed after de novo, presumably beneficial, mutations occurred in the presence of canavanine, a toxic analog of arginine. One highly parallel LOH event involved the ENA salt efflux pump locus on chromosome IV, which showed repeated LOH to the allele from the European parent, an allele originally derived by introgression from S. paradoxus Using CRISPR-engineered LOH we showed that the fitness advantage provided by this single LOH event was 27%. Overall, we found extensive evidence that LOH could be adaptive and is likely to be a greater source of initial variation than de novo mutation for rapid evolution of diploid genotypes.

Transcriptional initiation of a small RNA, not R-loop stability, dictates the frequency of pilin antigenic variation in Neisseria gonorrhoeae.

Neisseria gonorrhoeae, the sole causative agent of gonorrhea, constitutively undergoes diversification of the Type IV pilus. Gene conversion occurs between one of the several donor silent copies located in distinct loci and the recipient pilE gene, encoding the major pilin subunit of the pilus. A guanine quadruplex (G4) DNA structure and a cis-acting sRNA (G4-sRNA) are located upstream of the pilE gene and both are required for pilin antigenic variation (Av). We show that the reduced sRNA transcription lowers pilin Av frequencies. Extended transcriptional elongation is not required for Av, since limiting the transcript to 32 nt allows for normal Av frequencies. Using chromatin immunoprecipitation (ChIP) assays, we show that cellular G4s are less abundant when sRNA transcription is lower. In addition, using ChIP, we demonstrate that the G4-sRNA forms a stable RNA:DNA hybrid (R-loop) with its template strand. However, modulating R-loop levels by controlling RNase HI expression does not alter G4 abundance quantified through ChIP. Since pilin Av frequencies were not altered when modulating R-loop levels by controlling RNase HI expression, we conclude that transcription of the sRNA is necessary, but stable R-loops are not required to promote pilin Av.

The Role of Gene Conversion between Transposable Elements in Rewiring Regulatory Networks.

Nature has found many ways to utilize transposable elements (TEs) throughout evolution. Many molecular and cellular processes depend on DNA-binding proteins recognizing hundreds or thousands of similar DNA motifs dispersed throughout the genome that are often provided by TEs. It has been suggested that TEs play an important role in the evolution of such systems, in particular, the rewiring of gene regulatory networks. One mechanism that can further enhance the rewiring of regulatory networks is nonallelic gene conversion between copies of TEs. Here, we will first review evidence for nonallelic gene conversion in TEs. Then, we will illustrate the benefits nonallelic gene conversion provides in rewiring regulatory networks. For instance, nonallelic gene conversion between TE copies offers an alternative mechanism to spread beneficial mutations that improve the network, it allows multiple mutations to be combined and transferred together, and it allows natural selection to work efficiently in spreading beneficial mutations and removing disadvantageous mutations. Future studies examining the role of nonallelic gene conversion in the evolution of TEs should help us to better understand how TEs have contributed to evolution.

Topoisomerase II-Induced Chromosome Breakage and Translocation Is Determined by Chromosome Architecture and Transcriptional Activity.

Topoisomerase II (TOP2) relieves torsional stress by forming transient cleavage complex intermediates (TOP2ccs) that contain TOP2-linked DNA breaks (DSBs). While TOP2ccs are normally reversible, they can be "trapped" by chemotherapeutic drugs such as etoposide and subsequently converted into irreversible TOP2-linked DSBs. Here, we have quantified etoposide-induced trapping of TOP2ccs, their conversion into irreversible TOP2-linked DSBs, and their processing during DNA repair genome-wide, as a function of time. We find that while TOP2 chromatin localization and trapping is independent of transcription, it requires pre-existing binding of cohesin to DNA. In contrast, the conversion of trapped TOP2ccs to irreversible DSBs during DNA repair is accelerated 2-fold at transcribed loci relative to non-transcribed loci. This conversion is dependent on proteasomal degradation and TDP2 phosphodiesterase activity. Quantitative modeling shows that only two features of pre-existing chromatin structure-namely, cohesin binding and transcriptional activity-can be used to predict the kinetics of TOP2-induced DSBs.

The impact of poly-A microsatellite heterologies in meiotic recombination.

Meiotic recombination has strong, but poorly understood effects on short tandem repeat (STR) instability. Here, we screened thousands of single recombinant products with sperm typing to characterize the role of polymorphic poly-A repeats at a human recombination hotspot in terms of hotspot activity and STR evolution. We show that the length asymmetry between heterozygous poly-A's strongly influences the recombination outcome: a heterology of 10 A's (9A/19A) reduces the number of crossovers and elevates the frequency of non-crossovers, complex recombination products, and long conversion tracts. Moreover, the length of the heterology also influences the STR transmission during meiotic repair with a strong and significant insertion bias for the short heterology (6A/7A) and a deletion bias for the long heterology (9A/19A). In spite of this opposing insertion-/deletion-biased gene conversion, we find that poly-A's are enriched at human recombination hotspots that could have important consequences in hotspot activation.

Genome-wide analysis of genomic alterations induced by oxidative DNA damage in yeast.

Oxidative DNA damage is a threat to genome stability. Using a genetic system in yeast that allows detection of mitotic recombination, we found that the frequency of crossovers is greatly elevated when cells are treated with hydrogen peroxide (H2O2). Using a combination of microarray analysis and genomic sequencing, we mapped the breakpoints of mitotic recombination events and other chromosome rearrangements at a resolution of about 1 kb. Gene conversions and crossovers were the two most common types of events, but we also observed deletions, duplications, and chromosome aneuploidy. In addition, H2O2-treated cells had elevated rates of point mutations (particularly A to T/T to A and C to G/G to C transversions) and small insertions/deletions (in/dels). In cells that underwent multiple rounds of H2O2 treatments, we identified a genetic alteration that resulted in improved H2O2 tolerance by amplification of the CTT1 gene that encodes cytosolic catalase T. Lastly, we showed that cells grown in the absence of oxygen have reduced levels of recombination. This study provided multiple novel insights into how oxidative stress affects genomic instability and phenotypic evolution in aerobic cells.

Pervasive gene conversion in chromosomal inversion heterozygotes.

Chromosomal inversions shape recombination landscapes, and species differing by inversions may exhibit reduced gene flow in these regions of the genome. Though single crossovers within inversions are not usually recovered from inversion heterozygotes, the recombination barrier imposed by inversions is nuanced by noncrossover gene conversion. Here, we provide a genomewide empirical analysis of gene conversion rates both within species and in species hybrids. We estimate that gene conversion occurs at a rate of 1 × 10-5 to 2.5 × 10-5 converted sites per bp per generation in experimental crosses within Drosophila pseudoobscura and between D. pseudoobscura and its naturally hybridizing sister species D. persimilis. This analysis is the first direct empirical assessment of gene conversion rates within inversions of a species hybrid. Our data show that gene conversion rates in interspecies hybrids are at least as high as within-species estimates of gene conversion rates, and gene conversion occurs regularly within and around inverted regions of species hybrids, even near inversion breakpoints. We also found that several gene conversion events appeared to be mitotic rather than meiotic in origin. Finally, we observed that gene conversion rates are higher in regions of lower local sequence divergence, yet our observed gene conversion rates in more divergent inverted regions were at least as high as in less divergent collinear regions. Given our observed high rates of gene conversion despite the sequence differentiation between species, especially in inverted regions, gene conversion has the potential to reduce the efficacy of inversions as barriers to recombination over evolutionary time.