Association of poly-purine/poly-pyrimidine sequences with meiotic recombination hot spots

DNA polymerase ε and δ variants drive mutagenesis in polypurine tracts in human tumors

Alterations in the exonuclease domain of DNA polymerase ε cause ultramutated cancers. These cancers accumulate AGA>ATA transversions; however, their genomic features beyond the trinucleotide motifs are obscure. We analyze the extended DNA context of ultramutation using whole-exome sequencing data from 524 endometrial and 395 colorectal tumors. We find that G>T transversions in POLE-mutant tumors predominantly affect sequences containing at least six consecutive purines, with a striking preference for certain positions within polypurine tracts. Using this signature, we develop a machine-learning classifier to identify tumors with hitherto unknown POLE drivers and validate two drivers, POLE-E978G and POLE-S461L, by functional assays in yeast. Unlike other pathogenic variants, the E978G substitution affects the polymerase domain of Pol ε. We further show that tumors with POLD1 drivers share the extended signature of POLE ultramutation. These findings expand the understanding of ultramutation mechanisms and highlight peculiar mutagenic properties of polypurine tracts in the human genome.

Extreme clustering of type-1 NF1 deletion breakpoints co-locating with G-quadruplex forming sequences

The breakpoints of type-1 NF1 deletions encompassing 1.4-Mb are located within NF1-REPa and NF1-REPc, which exhibit a complex structure comprising different segmental duplications in direct and inverted orientation. Here, we systematically assessed the proportion of type-1 NF1 deletions caused by nonallelic homologous recombination (NAHR) and those mediated by other mutational mechanisms. To this end, we analyzed 236 unselected type-1 deletions and observed that 179 of them (75.8%) had breakpoints located within the NAHR hotspot PRS2, whereas 39 deletions (16.5%) had breakpoints located within PRS1. Sixteen deletions exhibited breakpoints located outside of these NAHR hotspots but were also mediated by NAHR. Taken together, the breakpoints of 234 (99.2%) of the 236 type-1 NF1 deletions were mediated by NAHR. Thus, NF1-REPa and NF1-REPc are strongly predisposed to recurrent NAHR, the main mechanism underlying type-1 NF1 deletions. We also observed a non-random overlap between type-1 NF1-deletion breakpoints and G-quadruplex forming sequences (GQs) as well as regions flanking PRDM9_A binding-sites. These findings imply that GQs and PRDM9_A binding-sites contribute to the clustering of type-1 deletion breakpoints. The co-location of both types of sequence was at its highest within PRS2, indicative of their synergistic contribution to the greatly increased NAHR activity within this hotspot.

The consequences of sequence erosion in the evolution of recombination hotspots

Meiosis is initiated by a double-strand break (DSB) introduced in the DNA by a highly controlled process that is repaired by recombination. In many organisms, recombination occurs at specific and narrow regions of the genome, known as recombination hotspots, which overlap with regions enriched for DSBs. In recent years, it has been demonstrated that conversions and mutations resulting from the repair of DSBs lead to a rapid sequence evolution at recombination hotspots eroding target sites for DSBs. We still do not fully understand the effect of this erosion in the recombination activity, but evidence has shown that the binding of trans-acting factors like PRDM9 is affected. PRDM9 is a meiosis-specific, multi-domain protein that recognizes DNA target motifs by its zinc finger domain and directs DSBs to these target sites. Here we discuss the changes in affinity of PRDM9 to eroded recognition sequences, and explain how these changes in affinity of PRDM9 can affect recombination, leading sometimes to sterility in the context of hybrid crosses. We also present experimental data showing that DNA methylation reduces PRDM9 binding in vitro. Finally, we discuss PRDM9-independent hotspots, posing the question how these hotspots evolve and change with sequence erosion.

This article is part of the themed issue ‘Evolutionary causes and consequences of recombination rate variation in sexual organisms’.

Probing the Potential Role of Non-B DNA Structures at Yeast Meiosis-Specific DNA Double-Strand Breaks

A plethora of evidence suggests that different types of DNA quadruplexes are widely present in the genome of all organisms. The existence of a growing number of proteins that selectively bind and/or process these structures underscores their biological relevance. Moreover, G-quadruplex DNA has been implicated in the alignment of four sister chromatids by forming parallel guanine quadruplexes during meiosis; however, the underlying mechanism is not well defined. Here we show that a G/C-rich motif associated with a meiosis-specific DNA double-strand break (DSB) in Saccharomyces cerevisiae folds into G-quadruplex, and the C-rich sequence complementary to the G-rich sequence forms an i-motif. The presence of G-quadruplex or i-motif structures upstream of the green fluorescent protein-coding sequence markedly reduces the levels of gfp mRNA expression in S. cerevisiae cells, with a concomitant decrease in green fluorescent protein abundance, and blocks primer extension by DNA polymerase, thereby demonstrating the functional significance of these structures. Surprisingly, although S. cerevisiae Hop1, a component of synaptonemal complex axial/lateral elements, exhibits strong affinity to G-quadruplex DNA, it displays a much weaker affinity for the i-motif structure. However, the Hop1 C-terminal but not the N-terminal domain possesses strong i-motif binding activity, implying that the C-terminal domain has a distinct substrate specificity. Additionally, we found that Hop1 promotes intermolecular pairing between G/C-rich DNA segments associated with a meiosis-specific DSB site. Our results support the idea that the G/C-rich motifs associated with meiosis-specific DSBs fold into intramolecular G-quadruplex and i-motif structures, both in vitro and in vivo, thus revealing an important link between non-B form DNA structures and Hop1 in meiotic chromosome synapsis and recombination.

Identification of genes containing expanded purine repeats in the human genome and their apparent protective role against cancer

Purine repeat sequences present in a gene are unique as they have high propensity to form unusual DNA-triple helix structures. Friedreich's ataxia is the only human disease that is well known to be associated with DNA-triplexes formed by purine repeats. The purpose of this study was to recognize the expanded purine repeats (EPRs) in human genome and find their correlation with cancer pathogenesis. We developed "PuRepeatFinder.pl" algorithm to identify non-overlapping EPRs without pyrimidine interruptions in the human genome and customized for searching repeat lengths, n ≥ 200. A total of 1158 EPRs were identified in the genome which followed Wakeby distribution. Two hundred and ninety-six EPRs were found in geneic regions of 282 genes (EPR-genes). Gene clustering of EPR-genes was done based on their cellular function and a large number of EPR-genes were found to be enzymes/enzyme modulators. Meta-analysis of 282 EPR-genes identified only 63 EPR-genes in association with cancer, mostly in breast, lung, and blood cancers. Protein-protein interaction network analysis of all 282 EPR-genes identified proteins including those in cadherins and VEGF. The two observations, that EPRs can induce mutations under malignant conditions and that identification of some EPR-gene products in vital cell signaling-mediated pathways, together suggest the crucial role of EPRs in carcinogenesis. The new link between EPR-genes and their functionally interacting proteins throws a new dimension in the present understanding of cancer pathogenesis and can help in planning therapeutic strategies. Validation of present results using techniques like NGS is required to establish the role of the EPR genes in cancer pathology.

Specific modifications of histone tails, but not DNA methylation, mirror the temporal variation of mammalian recombination hotspots

Recombination clusters nonuniformly across mammalian genomes at discrete genomic loci referred to as recombination hotspots. Despite their ubiquitous presence, individual hotspots rapidly lose their activities, and the molecular and evolutionary mechanisms underlying such frequent hotspot turnovers (the so-called “recombination hotspot paradox”) remain unresolved. Even though some sequence motifs are significantly associated with hotspots, multiple lines of evidence indicate that factors other than underlying sequences, such as epigenetic modifications, may affect the evolution of recombination hotspots. Thus, identifying epigenetic factors that covary with recombination at fine-scale is a promising step for this important research area. It was previously reported that recombination rates correlate with indirect measures of DNA methylation in the human genome. Here, we analyze experimentally determined DNA methylation and histone modification of human sperms, and show that the correlation between DNA methylation and recombination in long-range windows does not hold with respect to the spatial and temporal variation of recombination at hotspots. On the other hand, two histone modifications (H3K4me3 and H3K27me3) overlap extensively with recombination hotspots. Similar trends were observed in mice. These results indicate that specific histone modifications rather than DNA methylation are associated with the rapid evolution of recombination hotspots. Furthermore, many human recombination hotspots occupy “bivalent” chromatin regions that harbor both active (H3K4me3) and repressive (H3K27me3) marks. This may explain why human recombination hotspots tend to occur in nongenic regions, in contrast to yeast and Arabidopsis hotspots that are characterized by generally active chromatins. Our results highlight the dynamic epigenetic underpinnings of recombination hotspot evolution.

Small rare recurrent deletions and reciprocal duplications in 2q21.1, including brain-specific ARHGEF4 and GPR148

We have identified a rare small (~450 kb unique sequence) recurrent deletion in a previously linked attention-deficit hyperactivity disorder (ADHD) locus at 2q21.1 in five unrelated families with developmental delay (DD)/intellectual disability (ID), ADHD, epilepsy and other neurobehavioral abnormalities from 17 035 samples referred for clinical chromosomal microarray analysis. Additionally, a DECIPHER (http://decipher.sanger.ac.uk) patient 2311 was found to have the same deletion and presented with aggressive behavior. The deletion was not found in either six control groups consisting of 13 999 healthy individuals or in the DGV database. We have also identified reciprocal duplications in five unrelated families with autism, developmental delay (DD), seizures and ADHD. This genomic region is flanked by large, complex low-copy repeats (LCRs) with directly oriented subunits of ~109 kb in size that have 97.7% DNA sequence identity. We sequenced the deletion breakpoints within the directly oriented paralogous subunits of the flanking LCR clusters, demonstrating non-allelic homologous recombination as a mechanism of formation. The rearranged segment harbors five genes: GPR148, FAM123C, ARHGEF4, FAM168B and PLEKHB2. Expression of ARHGEF4 (Rho guanine nucleotide exchange factor 4) is restricted to the brain and may regulate the actin cytoskeletal network, cell morphology and migration, and neuronal function. GPR148 encodes a G-protein-coupled receptor protein expressed in the brain and testes. We suggest that small rare recurrent deletion of 2q21.1 is pathogenic for DD/ID, ADHD, epilepsy and other neurobehavioral abnormalities and, because of its small size, low frequency and more severe phenotype might have been missed in other previous genome-wide screening studies using single-nucleotide polymorphism analyses.

Searching microsatellites in DNA sequences: approaches used and tools developed

Microsatellite instability associated genomic activities and evolutionary changes have led to a renewed focus on microsatellite research. In last decade, a number of microsatellite mining tools have been introduced based on different computational approaches. The choice is generally made between slow but exhaustive dynamic programming based approaches, or fast and incomplete heuristic methods. Tools based on stochastic approaches are more popular due to their simplicity and added ornamental features. We have performed a comparative evaluation of the relative efficiency of some microsatellite search tools with their default settings. The graphical user interface, the statistical analysis of the output and ability to mine imperfect repeats are the most important criteria in selecting a tool for a particular investigation. However, none of the available tools alone provides complete and accurate information about microsatellites, and a lot depends on the discretion of the user.

High recombination rates and hotspots in a Plasmodium falciparum genetic cross

Background

The human malaria parasite Plasmodium falciparum survives pressures from the host immune system and antimalarial drugs by modifying its genome. Genetic recombination and nucleotide substitution are the two major mechanisms that the parasite employs to generate genome diversity. A better understanding of these mechanisms may provide important information for studying parasite evolution, immune evasion and drug resistance.

Results

Here, we used a high-density tiling array to estimate the genetic recombination rate among 32 progeny of a P. falciparum genetic cross (7G8 × GB4). We detected 638 recombination events and constructed a high-resolution genetic map. Comparing genetic and physical maps, we obtained an overall recombination rate of 9.6 kb per centimorgan and identified 54 candidate recombination hotspots. Similar to centromeres in other organisms, the sequences of P. falciparum centromeres are found in chromosome regions largely devoid of recombination activity. Motifs enriched in hotspots were also identified, including a 12-bp G/C-rich motif with 3-bp periodicity that may interact with a protein containing 11 predicted zinc finger arrays.

Conclusions

These results show that the P. falciparum genome has a high recombination rate, although it also follows the overall rule of meiosis in eukaryotes with an average of approximately one crossover per chromosome per meiosis. GC-rich repetitive motifs identified in the hotspot sequences may play a role in the high recombination rate observed. The lack of recombination activity in centromeric regions is consistent with the observations of reduced recombination near the centromeres of other organisms.

Recurrent distal 7q11.23 deletion including HIP1 and YWHAG identified in patients with intellectual disabilities, epilepsy, and neurobehavioral problems

We report 26 individuals from ten unrelated families who exhibit variable expression and/or incomplete penetrance of epilepsy, learning difficulties, intellectual disabilities, and/or neurobehavioral abnormalities as a result of a heterozygous microdeletion distally adjacent to the Williams-Beuren syndrome region on chromosome 7q11.23. In six families with a common recurrent ∼1.2 Mb deletion that includes the Huntingtin-interacting protein 1 (HIP1) and tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein gamma (YWHAG) genes and that is flanked by large complex low-copy repeats, we identified sites for nonallelic homologous recombination in two patients. There were no cases of this ∼1.2 Mb distal 7q11.23 deletion copy number variant identified in over 20,000 control samples surveyed. Three individuals with smaller, nonrecurrent deletions (∼180-500 kb) that include HIP1 but not YWHAG suggest that deletion of HIP1 is sufficient to cause neurological disease. Mice with targeted mutation in the Hip1 gene (Hip1⁻(/)⁻) develop a neurological phenotype characterized by failure to thrive, tremor, and gait ataxia. Overall, our data characterize a neurodevelopmental and epilepsy syndrome that is likely caused by recurrent and nonrecurrent deletions, including HIP1. These data do not exclude the possibility that YWHAG loss of function is also sufficient to cause neurological phenotypes. Based on the current knowledge of Hip1 protein function and its proposed role in AMPA and NMDA ionotropic glutamate receptor trafficking, we believe that HIP1 haploinsufficiency in humans will be amenable to rational drug design for improved seizure control and cognitive and behavioral function.

Important characteristics of sequence-specific recombination hotspots in Schizosaccharomyces pombe

In many organisms, meiotic recombination occurs preferentially at a limited number of sites in the genome known as hotspots. In the fission yeast Schizosaccharomyces pombe, simple sequence motifs determine the location of at least some, and possibly most or all, hotspots. Recently, we showed that a large number of different sequences can create hotspots. Among those sequences we identified some recurring motifs that fell into at least five distinct families, including the well-characterized CRE family of hotspots. Here we report the essential sequence for activity of two of the novel hotspots, the oligo-C and CCAAT hotspots, and identify associated trans-acting factors required for hotspot activity. The oligo-C hotspot requires a unique 8-bp sequence, CCCCGCAC, though hotspot activity is also significantly affected by adjacent nucleotides. The CCAAT hotspot requires a more complex and degenerate sequence, including the originally identified seven nucleotide CCAATCA sequence at its core. We identified transcription factors, the CCAAT-binding factor (CBF) and Rst2, which are required specifically for activity of the CCAAT hotspots and oligo-C hotspots, respectively. Each of these factors binds to its respective motifs in vitro. However, unlike CRE, the sequence required for hotspot activity is larger than the sequence required for binding, suggesting the involvement of additional factors.

Structures and molecular mechanisms for common 15q13.3 microduplications involving CHRNA7: benign or pathological?

We have investigated four approximately 1.6-Mb microduplications and 55 smaller 350-680-kb microduplications at 15q13.2-q13.3 involving the CHRNA7 gene that were detected by clinical microarray analysis. Applying high-resolution array-CGH, we mapped all 118 chromosomal breakpoints of these microduplications. We also sequenced 26 small microduplication breakpoints that were clustering at hotspots of nonallelic homologous recombination (NAHR). All four large microduplications likely arose by NAHR between BP4 and BP5 LCRs, and 54 small microduplications arose by NAHR between two CHRNA7-LCR copies. We identified two classes of approximately 1.6-Mb microduplications and five classes of small microduplications differing in duplication size, and show that they duplicate the entire CHRNA7. We propose that size differences among small microduplications result from preexisting heterogeneity of the common BP4-BP5 inversion. Clinical data and family histories of 11 patients with small microduplications involving CHRNA7 suggest that these microduplications might be associated with developmental delay/mental retardation, muscular hypotonia, and a variety of neuropsychiatric disorders. However, we conclude that these microduplications and their associated potential for increased dosage of the CHRNA7-encoded alpha 7 subunit of nicotinic acetylcholine receptors are of uncertain clinical significance at present. Nevertheless, if they prove to have a pathological effects, their high frequency could make them a common risk factor for many neurobehavioral disorders.

Novel nucleotide sequence motifs that produce hotspots of meiotic recombination in Schizosaccharomyces pombe

In many organisms, including yeasts and humans, meiotic recombination is initiated preferentially at a limited number of sites in the genome referred to as recombination hotspots. Predicting precisely the location of most hotspots has remained elusive. In this study, we tested the hypothesis that hotspots can result from multiple different sequence motifs. We devised a method to rapidly screen many short random oligonucleotide sequences for hotspot activity in the fission yeast Schizosaccharomyces pombe and produced a library of approximately 500 unique 15- and 30-bp sequences containing hotspots. The frequency of hotspots found suggests that there may be a relatively large number of different sequence motifs that produce hotspots. Within our sequence library, we found many shorter 6- to 10-bp motifs that occurred multiple times, many of which produced hotspots when reconstructed in vivo. On the basis of sequence similarity, we were able to group those hotspots into five different sequence families. At least one of the novel hotspots we found appears to be a target for a transcription factor, as it requires that factor for its hotspot activity. We propose that many hotspots in S. pombe, and perhaps other organisms, result from simple sequence motifs, some of which are identified here.

High frequency of microsatellites in S. cerevisiae meiotic recombination hotspots

Background

Microsatellites are highly abundant in eukaryotic genomes but their function and evolution are not yet well understood. Their elevated mutation rate makes them ideal markers of genetic difference, but high levels of unexplained heterogeneity in mutation rates among microsatellites at different genomic locations need to be elucidated in order to improve the power and accuracy of the many types of study that use them as genetic markers. Recombination could contribute to this heterogeneity, since while replication errors are thought to be the predominant mechanism for microsatellite mutation, meiotic recombination is involved in some mutation events. There is also evidence suggesting that microsatellites could function as recombination signals. The yeast S. cerevisiae is a useful model organism with which to further explore the link between microsatellites and recombination, since it is very amenable to genetic study, and meiotic recombination hotspots have been mapped throughout its entire genome.

Results

We examined in detail the relationship between microsatellites and hotspots of meiotic double-strand breaks, the precursors of meiotic recombination, throughout the S. cerevisiae genome. We included all tandem repeats with motif length (repeat period) between one and six base pairs. Long, short and two-copy arrays were considered separately. We found that long, mono-, di- and trinucleotide microsatellites are around twice as frequent in hot than non-hot intergenic regions. The associations are weak or absent for repeats with less than six copies, and also for microsatellites with 4–6 base pair motifs, but high-copy arrays with motif length greater than three are relatively very rare throughout the genome. We present evidence that the association between high-copy, short-motif microsatellites and recombination hotspots is not driven by effects on microsatellite distribution of other factors previously linked to both recombination and microsatellites, including transcription, GC-content and transposable elements.

Conclusion

Our findings suggest that a mutation bias relating to recombination hotspots causing repeats to form and grow, and/or regulation of a subset of hotspots by simple sequences, may be significant processes in yeast. Some previous evidence has cast doubt on both of these possibilities, and as a result they have not been explored on a large scale, but the strength of the association we report suggests that they deserve further experimental testing.

The rise, fall and renaissance of microsatellites in eukaryotic genomes

Microsatellites are among the most versatile of genetic markers, being used in an impressive number of biological applications. However, the evolutionary dynamics of these markers remain a source of contention. Almost 20 years after the discovery of these ubiquitous simple sequences, new genomic data are clarifying our understanding of the structure, distribution and variability of microsatellites in genomes, especially for the eukaryotes. While these new data provide a great deal of descriptive information about the nature and abundance of microsatellite sequences within eukaryotic genomes, there have been few attempts to synthesise this information to develop a global concept of evolution. This review provides an up-to-date account of the mutational processes, biases and constraints believed to be involved in the evolution of microsatellites, particularly with respect to the creation and degeneration of microsatellites, which we assert may be broadly viewed as a life cycle. In addition, we identify areas of contention that require further research and propose some possible directions for future investigation.