Molecular features of meiotic recombination hot spots

Chromatin-mediated regulators of meiotic recombination revealed by proteomics of a recombination hotspot

Background

Meiotic recombination hotspots control the frequency and distribution of Spo11 (Rec12)-initiated recombination in the genome. Recombination occurs within and is regulated in part by chromatin structure, but relatively few of the many chromatin remodeling factors and histone posttranslational modifications (PTMs) have been interrogated for a role in the process.

Results

We developed a chromatin affinity purification and mass spectrometry-based approach to identify proteins and histone PTMs that regulate recombination hotspots. Small (4.2 kbp) minichromosomes (MiniCs) bearing the fission yeast ade6-M26 hotspot or a basal recombination control were purified approximately 100,000-fold under native conditions from meiosis; then, associated proteins and histone PTMs were identified by mass spectrometry. Proteins and PTMs enriched at the hotspot included known regulators (Atf1, Pcr1, Mst2, Snf22, H3K14ac), validating the approach. The abundance of individual histones varied dynamically during meiotic progression in hotspot versus basal control MiniCs, as did a subset of 34 different histone PTMs, implicating these as potential regulators. Measurements of basal and hotspot recombination in null mutants confirmed that additional, hotspot-enriched proteins are bona fide regulators of hotspot activation within the genome. These chromatin-mediated regulators include histone H2A-H2B and H3-H4 chaperones (Nap1, Hip1/Hir1), subunits of the Ino80 complex (Arp5, Arp8), a DNA helicase/E3 ubiquitin ligase (Rrp2), components of a Swi2/Snf2 family remodeling complex (Swr1, Swc2), and a nucleosome evictor (Fft3/Fun30).

Conclusions

Overall, our findings indicate that a remarkably diverse collection of chromatin remodeling factors and histone PTMs participate in designating where meiotic recombination occurs in the genome, and they provide new insight into molecular mechanisms of the process.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0233-x) contains supplementary material, which is available to authorized users.

High-Resolution Mapping of Crossover and Non-crossover Recombination Events by Whole-Genome Re-sequencing of an Avian Pedigree

Recombination is an engine of genetic diversity and therefore constitutes a key process in evolutionary biology and genetics. While the outcome of crossover recombination can readily be detected as shuffled alleles by following the inheritance of markers in pedigreed families, the more precise location of both crossover and non-crossover recombination events has been difficult to pinpoint. As a consequence, we lack a detailed portrait of the recombination landscape for most organisms and knowledge on how this landscape impacts on sequence evolution at a local scale. To localize recombination events with high resolution in an avian system, we performed whole-genome re-sequencing at high coverage of a complete three-generation collared flycatcher pedigree. We identified 325 crossovers at a median resolution of 1.4 kb, with 86% of the events localized to <10 kb intervals. Observed crossover rates were in excellent agreement with data from linkage mapping, were 52% higher in male (3.56 cM/Mb) than in female meiosis (2.28 cM/Mb), and increased towards chromosome ends in male but not female meiosis. Crossover events were non-randomly distributed in the genome with several distinct hot-spots and a concentration to genic regions, with the highest density in promoters and CpG islands. We further identified 267 non-crossovers, whose location was significantly associated with crossover locations. We detected a significant transmission bias (0.18) in favour of 'strong' (G, C) over 'weak' (A, T) alleles at non-crossover events, providing direct evidence for the process of GC-biased gene conversion in an avian system. The approach taken in this study should be applicable to any species and would thereby help to provide a more comprehensive portray of the recombination landscape across organism groups.

Recombination hotspots: Models and tools for detection

Recombination hotspots are the regions within the genome where the rate, and the frequency of recombination are optimum with a size varying from 1 to 2kb. The recombination event is mediated by the double-stranded break formation, guided by the combined enzymatic action of DNA topoisomerase and Spo 11 endonuclease. These regions are distributed non-uniformly throughout the human genome and cause distortions in the genetic map. Numerous lines of evidence suggest that the number of hotspots known in humans has increased manifold in recent years. A few facts about the hotspot evolutions were also put forward, indicating the differences in the hotspot position between chimpanzees and humans. In mice, recombination hot spots were found to be clustered within the major histocompatibility complex (MHC) region. Several models, that help explain meiotic recombination has been proposed. Moreover, scientists also developed some computational tools to locate the hotspot position and estimate their recombination rate in humans is of great interest to population and medical geneticists. Here we reviewed the molecular mechanisms, models and in silico prediction techniques of hot spot residues.

Meiotic DSB patterning: A multifaceted process

Meiosis is a specialized two-step cell division responsible for genome haploidization and the generation of genetic diversity during gametogenesis. An integral and distinctive feature of the meiotic program is the evolutionarily conserved initiation of homologous recombination (HR) by the developmentally programmed induction of DNA double-strand breaks (DSBs). The inherently dangerous but essential act of DSB formation is subject to multiple forms of stringent and self-corrective regulation that collectively ensure fruitful and appropriate levels of genetic exchange without risk to cellular survival. Within this article we focus upon an emerging element of this control--spatial regulation--detailing recent advances made in understanding how DSBs are evenly distributed across the genome, and present a unified view of the underlying patterning mechanisms employed.

Organization and roles of nucleosomes at mouse meiotic recombination hotspots

Meiotic double strand breaks (DSBs) occur at discrete regions in the genome coined hotspots. Precisely what directs site selection of these DSBs is hotly debated and in particular it is unclear which chromatin features, and regulatory factors are necessary for a genomic region to initiate and resolve DSBs as a crossover (CO) event. In human and mouse, one layer of hotspot selection control is a recognition sequence element present at these sites that is bound by the Prdm9 zinc-finger protein. Furthermore, an overall open chromatin structure is thought to be required to allow access of the recombination machinery, and this is often dictated by the packaging of DNA around nucleosomes. We recently defined the nucleosome occupancy maps of four mouse recombination hotspots throughout meiosis. These analyses revealed no obvious dynamic changes in nucleosome occupancy, suggesting an intrinsic nature of recombinogenic sites, yet they also revealed that nucleosomes define zones of exclusion for CO resolution. Here, we discuss new evidence implicating nucleosome occupancy in recombinogenic repair and its potential roles in controlling chromatin structure at mouse meiotic hotspots.

Prdm9, a major determinant of meiotic recombination hotspots, is not functional in dogs and their wild relatives, wolves and coyotes

Meiotic recombination is a fundamental process needed for the correct segregation of chromosomes during meiosis in sexually reproducing organisms. In humans, 80% of crossovers are estimated to occur at specific areas of the genome called recombination hotspots. Recently, a protein called PRDM9 was identified as a major player in determining the location of genome-wide meiotic recombination hotspots in humans and mice. The origin of this protein seems to be ancient in evolutionary time, as reflected by its fairly conserved structure in lineages that diverged over 700 million years ago. Despite its important role, there are many animal groups in which Prdm9 is absent (e.g. birds, reptiles, amphibians, diptera) and it has been suggested to have disruptive mutations and thus to be a pseudogene in dogs. Because of the dog's history through domestication and artificial selection, we wanted to confirm the presence of a disrupted Prdm9 gene in dogs and determine whether this was exclusive of this species or whether it also occurred in its wild ancestor, the wolf, and in a close relative, the coyote. We sequenced the region in the dog genome that aligned to the last exon of the human Prdm9, containing the entire zinc finger domain, in 4 dogs, 17 wolves and 2 coyotes. Our results show that the three canid species possess mutations that likely make this gene non functional. Because these mutations are shared across the three species, they must have appeared prior to the split of the wolf and the coyote, millions of years ago, and are not related to domestication. In addition, our results suggest that in these three canid species recombination does not occur at hotspots or hotspot location is controlled through a mechanism yet to be determined.

Nucleosomal organization of replication origins and meiotic recombination hotspots in fission yeast

In Schizosaccharomyces pombe, DNA replication origins (ORIs) and meiotic recombination hotspots lack consensus sequences and show a bias towards mapping to large intergenic regions (IGRs). To explore whether this preference depended on underlying chromatin features, we have generated genome-wide nucleosome profiles during mitosis and meiosis. We have found that meiotic double-strand break sites (DSBs) colocalize with nucleosome-depleted regions (NDRs) and that large IGRs include clusters of NDRs that overlap with almost half of all DSBs. By contrast, ORIs do not colocalize with NDRs and they are regulated independently of DSBs. Physical relocation of NDRs at ectopic loci or modification of their genomic distribution during meiosis was paralleled by the generation of new DSB sites. Over 80% of all meiotic DSBs colocalize with NDRs that are also present during mitosis, indicating that the recombination pattern is largely dependent on constitutive properties of the genome and, to a lesser extent, on the transcriptional profile during meiosis. The organization of ORIs and of DSBs regions in S. pombe reveals similarities and differences relative to Saccharomyces cerevisiae.

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.

Integrative mapping analysis of chicken microchromosome 16 organization

BACKGROUND

The chicken karyotype is composed of 39 chromosome pairs, of which 9 still remain totally absent from the current genome sequence assembly, despite international efforts towards complete coverage. Some others are only very partially sequenced, amongst which microchromosome 16 (GGA16), particularly under-represented, with only 433 kb assembled for a full estimated size of 9 to 11 Mb. Besides the obvious need of full genome coverage with genetic markers for QTL (Quantitative Trait Loci) mapping and major genes identification studies, there is a major interest in the detailed study of this chromosome because it carries the two genetically independent MHC complexes B and Y. In addition, GGA16 carries the ribosomal RNA (rRNA) genes cluster, also known as the NOR (nucleolus organizer region). The purpose of the present study is to construct and present high resolution integrated maps of GGA16 to refine its organization and improve its coverage with genetic markers.

RESULTS

We developed 79 STS (Sequence Tagged Site) markers to build a physical RH (radiation hybrid) map and 34 genetic markers to extend the genetic map of GGA16. We screened a BAC (Bacterial Artificial Chromosome) library with markers for the MHC-B, MHC-Y and rRNA complexes. Selected clones were used to perform high resolution FISH (Fluorescent In Situ Hybridization) mapping on giant meiotic lampbrush chromosomes, allowing meiotic mapping in addition to the confirmation of the order of the three clusters along the chromosome. A region with high recombination rates and containing PO41 repeated elements separates the two MHC complexes.

CONCLUSIONS

The three complementary mapping strategies used refine greatly our knowledge of chicken microchromosome 16 organisation. The characterisation of the recombination hotspots separating the two MHC complexes demonstrates the presence of PO41 repetitive sequences both in tandem and inverted orientation. However, this region still needs to be studied in more detail.

Detecting sequence polymorphisms associated with meiotic recombination hotspots in the human genome

Background

Meiotic recombination events tend to cluster into narrow spans of a few kilobases long, called recombination hotspots. Such hotspots are not conserved between human and chimpanzee and vary between different human ethnic groups. At the same time, recombination hotspots are heritable. Previous studies showed instances where differences in recombination rate could be associated with sequence polymorphisms.

Results

In this work we developed a novel computational approach, LDsplit, to perform a large-scale association study of recombination hotspots with genetic polymorphisms. LDsplit was able to correctly predict the association between the FG11 SNP and the DNA2 hotspot observed by sperm typing. Extensive simulation demonstrated the accuracy of LDsplit under various conditions. Applying LDsplit to human chromosome 6, we found that for a significant fraction of hotspots, there is an association between variations in intensity of historical recombination and sequence polymorphisms. From flanking regions of the SNPs output by LDsplit we identified a conserved 11-mer motif GGNGGNAGGGG, whose complement partially matches 13-mer CCNCCNTNNCCNC, a critical motif for the regulation of recombination hotspots.

Conclusions

Our result suggests that computational approaches based on historical recombination events are likely to be more powerful than previously anticipated. The putative associations we identified may be a promising step toward uncovering the mechanisms of recombination hotspots.

Surprising fitness consequences of GC-biased gene conversion: I. Mutation load and inbreeding depression

GC-biased gene conversion (gBGC) is a recombination-associated process mimicking selection in favor of G and C alleles. It is increasingly recognized as a widespread force in shaping the genomic nucleotide landscape. In recombination hotspots, gBGC can lead to bursts of fixation of GC nucleotides and to accelerated nucleotide substitution rates. It was recently shown that these episodes of strong gBGC could give spurious signatures of adaptation and/or relaxed selection. There is also evidence that gBGC could drive the fixation of deleterious amino acid mutations in some primate genes. This raises the question of the potential fitness effects of gBGC. While gBGC has been metaphorically termed the "Achilles' heel" of our genome, we do not know whether interference between gBGC and selection merely has practical consequences for the analysis of sequence data or whether it has broader fundamental implications for individuals and populations. I developed a population genetics model to predict the consequences of gBGC on the mutation load and inbreeding depression. I also used estimates available for humans to quantitatively evaluate the fitness impact of gBGC. Surprising features emerged from this model: (i) Contrary to classical mutation load models, gBGC generates a fixation load independent of population size and could contribute to a significant part of the load; (ii) gBGC can maintain recessive deleterious mutations for a long time at intermediate frequency, in a similar way to overdominance, and these mutations generate high inbreeding depression, even if they are slightly deleterious; (iii) since mating systems affect both the selection efficacy and gBGC intensity, gBGC challenges classical predictions concerning the interaction between mating systems and deleterious mutations, and gBGC could constitute an additional cost of outcrossing; and (iv) if mutations are biased toward A and T alleles, very low gBGC levels can reduce the load. A robust prediction is that the gBGC level minimizing the load depends only on the mutational bias and population size. These surprising results suggest that gBGC may have nonnegligible fitness consequences and could play a significant role in the evolution of genetic systems. They also shed light on the evolution of gBGC itself.

Discrete DNA sites regulate global distribution of meiotic recombination

Homologous recombination is induced to high levels in meiosis, is initiated by Spo11-catalyzed DNA double-strand breaks (DSBs) and is clustered at hotspots that regulate its positioning in the genome. Recombination is required for proper chromosome segregation in meiosis and defects in its frequency or positioning cause chromosome mis-segregation and, consequently, congenital birth defects such as Down's syndrome. Therefore, elucidating how meiotic recombination is positioned is of fundamental and biomedical interest. Our integration of historical and contemporary advances in the field, plus the re-analysis of published microarray data on the genome-wide distribution of recombination supports a unifying model for such regulation. We posit that discrete DNA sequence motifs position and regulate essentially all recombination across the genome, in much the same way that DNA sites position and regulate transcription. Moreover, we illustrate the use of overlapping mechanisms for the regulation of transcription and meiotic recombination. Bound transcription factors induce histone modifications that position recombination at hotspots.

Gene conversion causing human inherited disease: evidence for involvement of non-B-DNA-forming sequences and recombination-promoting motifs in DNA breakage and repair

A variety of DNA sequence motifs including inverted repeats, minisatellites, and the chi recombination hotspot, have been reported in association with gene conversion in human genes causing inherited disease. However, no methodical statistically based analysis has been performed to formalize these observations. We have performed an in silico analysis of the DNA sequence tracts involved in 27 nonoverlapping gene conversion events in 19 different genes reported in the context of inherited disease. We found that gene conversion events tend to occur within (C+G)- and CpG-rich regions and that sequences with the potential to form non-B-DNA structures, and which may be involved in the generation of double-strand breaks that could, in turn, serve to promote gene conversion, occur disproportionately within maximal converted tracts and/or short flanking regions. Maximal converted tracts were also found to be enriched (P<0.01) in a truncated version of the chi-element (a TGGTGG motif), immunoglobulin heavy chain class switch repeats, translin target sites and several novel motifs including (or overlapping) the classical meiotic recombination hotspot, CCTCCCCT. Finally, gene conversions tend to occur in genomic regions that have the potential to fold into stable hairpin conformations. These findings support the concept that recombination-inducing motifs, in association with alternative DNA conformations, can promote recombination in the human genome.

Evolution of mutational robustness in the yeast genome: a link to essential genes and meiotic recombination hotspots

Deleterious mutations inevitably emerge in any evolutionary process and are speculated to decisively influence the structure of the genome. Meiosis, which is thought to play a major role in handling mutations on the population level, recombines chromosomes via non-randomly distributed hot spots for meiotic recombination. In many genomes, various types of genetic elements are distributed in patterns that are currently not well understood. In particular, important (essential) genes are arranged in clusters, which often cannot be explained by a functional relationship of the involved genes. Here we show by computer simulation that essential gene (EG) clustering provides a fitness benefit in handling deleterious mutations in sexual populations with variable levels of inbreeding and outbreeding. We find that recessive lethal mutations enforce a selective pressure towards clustered genome architectures. Our simulations correctly predict (i) the evolution of non-random distributions of meiotic crossovers, (ii) the genome-wide anti-correlation of meiotic crossovers and EG clustering, (iii) the evolution of EG enrichment in pericentromeric regions and (iv) the associated absence of meiotic crossovers (cold centromeres). Our results furthermore predict optimal crossover rates for yeast chromosomes, which match the experimentally determined rates. Using a Saccharomyces cerevisiae conditional mutator strain, we show that haploid lethal phenotypes result predominantly from mutation of single loci and generally do not impair mating, which leads to an accumulation of mutational load following meiosis and mating. We hypothesize that purging of deleterious mutations in essential genes constitutes an important factor driving meiotic crossover. Therefore, the increased robustness of populations to deleterious mutations, which arises from clustered genome architectures, may provide a significant selective force shaping crossover distribution. Our analysis reveals a new aspect of the evolution of genome architectures that complements insights about molecular constraints, such as the interference of pericentromeric crossovers with chromosome segregation.

Phosphorylation-independent regulation of Atf1-promoted meiotic recombination by stress-activated, p38 kinase Spc1 of fission yeast

Background

Stress-activated protein kinases regulate multiple cellular responses to a wide variety of intracellular and extracellular conditions. The conserved, multifunctional, ATF/CREB protein Atf1 (Mts1, Gad7) of fission yeast binds to CRE-like (M26) DNA sites. Atf1 is phosphorylated by the conserved, p38-family kinase Spc1 (Sty1, Phh1) and is required for many Spc1-dependent stress responses, efficient sexual differentiation, and activation of Rec12 (Spo11)-dependent meiotic recombination hotspots like ade6-M26.

Methodology/Principal Findings

We sought to define mechanisms by which Spc1 regulates Atf1 function at the ade6-M26 hotspot. The Spc1 kinase was essential for hotspot activity, but dispensable for basal recombination. Unexpectedly, a protein lacking all eleven MAPK phospho-acceptor sites and detectable phosphorylation (Atf1-11M) was fully proficient for hotspot recombination. Furthermore, tethering of Atf1 to ade6 in the chromosome by a heterologous DNA binding domain bypassed the requirement for Spc1 in promoting recombination.

Conclusions/Significance

The Spc1 protein kinase regulates the pathway of Atf1-promoted recombination at or before the point where Atf1 binds to chromosomes, and this pathway regulation is independent of the phosphorylation status of Atf1. Since basal recombination is Spc1-independent, the principal function of the Spc1 kinase in meiotic recombination is to correctly position Atf1-promoted recombination at hotspots along chromosomes. We also propose new hypotheses on regulatory mechanisms for shared (e.g., DNA binding) and distinct (e.g., osmoregulatory vs. recombinogenic) activities of multifunctional, stress-activated protein Atf1.

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.

A bird's-eye view of sex chromosome dosage compensation

Intensive study of a few genetically tractable species with XX/XY sex chromosomes has produced generalizations about the process of sex chromosome dosage compensation that do not fare well when applied to ZZ/ZW sex chromosome systems, such as those in birds. The inherent sexual imbalance in dose of sex chromosome genes has led to the evolution of sex-chromosome-wide mechanisms for balancing gene dosage between the sexes and relative to autosomal genes. Recent advances in our knowledge of avian genomes have led to a reexamination of sex-specific dosage compensation (SSDC) in birds, which is less effective than in known XX/XY systems. Insights about the mechanisms of SSDC in birds also suggest similarities to and differences from those in XX/XY species. Birds are thus offering new opportunities for studying dosage compensation in a ZZ/ZW system, which should shed light on the evolution of SSDC more broadly.

Cut thy neighbor: cyclic birth and death of recombination hotspots via genetic conflict

Most recombination takes place in numerous, localized regions called hotspots. However, empirical evidence indicates that nascent hotspots are susceptible to removal due to biased gene conversion, so it is paradoxical that they should be so widespread. Previous modeling work has shown that hotspots can evolve due to genetic drift overpowering their intrinsic disadvantage. Here we synthesize recent theoretical and empirical results to show how natural selection can favor hotspots. We propose that hotspots are part of a cycle of antagonistic coevolution between two tightly linked chromosomal regions: an inducer region that initiates recombination during meiosis by cutting within a nearby region of DNA and the cut region itself, which can evolve to be resistant to cutting. Antagonistic coevolution between inducers and their cut sites is driven by recurrent episodes of Hill-Robertson interference, genetic hitchhiking, and biased gene conversion.

Meiotic recombination hotspots of fission yeast are directed to loci that express non-coding RNA

Background

Polyadenylated, mRNA-like transcripts with no coding potential are abundant in eukaryotes, but the functions of these long non-coding RNAs (ncRNAs) are enigmatic. In meiosis, Rec12 (Spo11) catalyzes the formation of dsDNA breaks (DSBs) that initiate homologous recombination. Most meiotic recombination is positioned at hotspots, but knowledge of the mechanisms is nebulous. In the fission yeast genome DSBs are located within 194 prominent peaks separated on average by 65-kbp intervals of DNA that are largely free of DSBs.

Methodology/Principal Findings

We compared the genome-wide distribution of DSB peaks to that of polyadenylated ncRNA molecules of the prl class. DSB peaks map to ncRNA loci that may be situated within ORFs, near the boundaries of ORFs and intergenic regions, or most often within intergenic regions. Unconditional statistical tests revealed that this colocalization is non-random and robust (P≤5.5×10⁻⁸). Furthermore, we tested and rejected the hypothesis that the ncRNA loci and DSB peaks localize preferentially, but independently, to a third entity on the chromosomes.

Conclusions/Significance

Meiotic DSB hotspots are directed to loci that express polyadenylated ncRNAs. This reveals an unexpected, possibly unitary mechanism for what directs meiotic recombination to hotspots. It also reveals a likely biological function for enigmatic ncRNAs. We propose specific mechanisms by which ncRNA molecules, or some aspect of RNA metabolism associated with ncRNA loci, help to position recombination protein complexes at DSB hotspots within chromosomes.

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.

Variation in patterns of human meiotic recombination

In the last 30 years it has become evident that patterns of meiotic recombination can be highly variable among individuals. The evidence comes from both low and high resolution analyses of hotspots of recombination in human and other species. In addition, a comparison of the recombination profiles in closely related species such as human and chimpanzee reveals essentially no correlation in the position of hotspots. Although the variation in hotspots of meiotic recombination is clearly documented, the mechanisms responsible for such variation are far from being understood. Here we will review the available evidence of natural variation in meiotic recombination and will discuss potential implications of this variation on the functional mechanisms of crossover formation and control.

Gene-interleaving patterns of synteny in the Saccharomyces cerevisiae genome: are they proof of an ancient genome duplication event?

Background

Recent comparative genomic studies claim local syntenic gene-interleaving relationships in Ashbya gossypii and Kluyveromyces waltii are compelling evidence for an ancient whole-genome duplication event in Saccharomyces cerevisiae. We here test, using Hannenhalli-Pevzner rearrangement algorithms that address the multiple genome rearrangement problem, whether syntenic patterns are proof of paleopolyploidization.

Results

We focus on (1) pairwise comparison of gene arrangement sequences in A. gossypii and S. cerevisiae, (2) reconstruction of gene arrangements ancestral to A. gossypii, S. cerevisiae, and K. waltii, (3) synteny patterns arising within and between lineages, and (4) expected gene orientation of duplicate gene sets. The existence of syntenic patterns between ancestral gene sets and A. gossypii, S. cerevisiae, and K. waltii, and other evidence, suggests that gene-interleaving relationships are the natural consequence of topological rearrangements in chromosomes and that a more gradual scenario of genome evolution involving segmental duplication and recombination constitutes a more parsimonious explanation. Furthermore, phylogenetic trees reconstructed under alternative hypotheses placed the putative whole-genome duplication event after the divergence of the S. cerevisiae and K. waltii lineages, but in the lineage leading to K. waltii. This is clearly incompatible with an ancient genome duplication event in S. cerevisiae.

Conclusion

Because the presence of syntenic patterns appears to be a condition that is necessary, but not sufficient, to support the existence of the whole-genome duplication event, our results prompt careful re-evaluation of paleopolyploidization in the yeast lineage and the evolutionary meaning of syntenic patterns.

Reviewers

This article was reviewed by Kenneth H. Wolfe (nominated by Nicolas Galtier), Austin L. Hughes (nominated by Eugene Koonin), Mikhail S. Gelfand, and Mark Gerstein.

Association of Spot14α Gene Polymorphisms with Body Weight in the Chicken

In mammals, thyroid hormone responsive Spot14 (THRSP) is a small acidic protein that responds to thyroid hormone stimulation and, therefore, is thought to play a role in growth. The current study was designed to investigate the associations of Spot14alpha gene polymorphisms on chicken growth and body composition traits. The Northeast Agricultural University Resource Population (NEAURP) was used in the present study. The NEAURP was established by crossing broiler sires, derived from Northeast Agricultural University broiler lines divergently selected for abdominal fat content, with Baier layer dams, a local Chinese breed. The F(1) birds were intercrossed to produce the F(2) population. Body weight and body composition traits were measured in the F(2) population. Polymorphisms of the gene were detected between parental lines by DNA sequencing. Primers were designed according to the chicken Spot14alpha gene (AY568628). The PCR-RFLP and PCR-length polymorphisms methods were then developed to genotype polymorphisms in the NEAURP. The A213C and 9 bp insertion-deletion of the Spot14alpha gene in the F(2) population was found to be associated with BW, which implied that Spot14alpha gene or a tightly linked gene had an important effect on growth in the chicken.

A high-resolution linkage map for the Z chromosome in chicken reveals hot spots for recombination

A comprehensive linkage map for chicken chromosome Z was constructed as the result of a large-scale screening of single nucleotide polymorphisms (SNPs). A total of 308 SNPs were assigned to Z based on the genotype distribution among 182 birds representing several populations. A linkage map comprising 210 markers and spanning 200.9 cM was established by analyzing a small Red junglefowl/White Leghorn intercross. There was excellent agreement between the linkage map for Z and a recently released assembly of the chicken genome (May 2006). Almost all SNPs assigned to chromosome Z in the present study are on Z in the new genome assembly. The remaining 12 loci are all found on unassigned contigs that can now be assigned to Z. The average recombination rate was estimated at 2.7 cM/Mb but there was a very uneven distribution of recombination events with both cold and hot spots of recombination. The existence of one of the major hot spots of recombination, located around position 39.4 Mb, was supported by the observed pattern of linkage disequilibrium. Thirteen markers from unassigned contigs were shown to be located on chromosome W. Three of these contigs included genes that have homologues on chromosome Z. The preliminary assignment of three more genes to the gene-poor W chromosome may be important for studies on the mechanism of sex determination and dosage compensation in birds.

Identification of a mitotic recombination hotspot on chromosome III of the asexual fungus Aspergillus niger and its possible correlation with [corrected] elevated basal transcription

Genetic recombination is an important tool in strain breeding in many organisms. We studied the possibilities of mitotic recombination in strain breeding of the asexual fungus Aspergillus niger. By identifying genes that complemented mapped auxotrophic mutations, the physical map was compared to the genetic map of chromosome III using the genome sequence. In a program to construct a chromosome III-specific marker strain by selecting mitotic crossing-over in diploids, a mitotic recombination hotspot was identified. Analysis of the mitotic recombination hotspot revealed some physical features, elevated basal transcription and a possible correlation with purine stretches.

Understanding the recent evolution of the human genome: insights from human-chimpanzee genome comparisons

The sequencing of the chimpanzee genome and the comparison with its human counterpart have begun to reveal the spectrum of genetic changes that has accompanied human evolution. In addition to gross karyotypic rearrangements such as the fusion that formed human chromosome 2 and the human-specific pericentric inversions of chromosomes 1 and 18, there is considerable submicroscopic structural variation involving deletions, duplications, and inversions. Lineage-specific segmental duplications, detected by array comparative genomic hybridization and direct sequence comparison, have made a very significant contribution to this structural divergence, which is at least three-fold greater than that due to nucleotide substitutions. Since structural genomic changes may have given rise to irreversible functional differences between the diverging species, their detailed analysis could help to identify the biological processes that have accompanied speciation. To this end, interspecies comparisons have revealed numerous human-specific gains and losses of genes as well as changes in gene expression. The very considerable structural diversity (polymorphism) evident within both lineages has, however, hampered the analysis of the structural divergence between the human and chimpanzee genomes. The concomitant evaluation of genetic divergence and diversity at the nucleotide level has nevertheless served to identify many genes that have evolved under positive selection and may thus have been involved in the development of human lineage-specific traits. Genes that display signs of weak negative selection have also been identified and could represent candidate loci for complex genomic disorders. Here, we review recent progress in comparing the human and chimpanzee genomes and discuss how the differences detected have improved our understanding of the evolution of the human genome.

Recombination hotspots flank the Cryptococcus mating-type locus: implications for the evolution of a fungal sex chromosome

Recombination increases dramatically during meiosis to promote genetic exchange and generate recombinant progeny. Interestingly, meiotic recombination is unevenly distributed throughout genomes, and, as a consequence, genetic and physical map distances do not have a simple linear relationship. Recombination hotspots and coldspots have been described in many organisms and often reflect global features of chromosome structure. In particular, recombination frequencies are often distorted within or outside sex-determining regions of the genome. Here, we report that recombination is elevated adjacent to the mating-type locus (MAT) in the pathogenic basidiomycete Cryptococcus neoformans. Among fungi, C. neoformans has an unusually large MAT locus, and recombination is suppressed between the two >100-kilobase mating-type specific alleles. When genetic markers were introduced at defined physical distances from MAT, we found the meiotic recombination frequency to be ~20% between MAT and a flanking marker at 5, 10, 50, or 100 kilobases from the right border. As a result, the physical/genetic map ratio in the regions adjacent to MAT is distorted ~10- to 50-fold compared to the genome-wide average. Moreover, recombination frequently occurred on both sides of MAT and negative interference between crossovers was observed. MAT heterozygosity was not required for enhanced recombination, implying that this process is not due to a physical distortion from the two non-paired alleles and could also occur during same-sex mating. Sequence analysis revealed a correlation between high G + C content and these hotspot regions. We hypothesize that the presence of recombinational activators may have driven several key events during the assembly and reshaping of the MAT locus and may have played similar roles in the origins of both metabolic and biosynthetic gene clusters. Our findings suggest that during meiosis the MAT locus may be exchanged onto different genetic backgrounds and therefore have broad evolutionary implications with respect to mating-type switching in both model and pathogenic yeasts.

It is hypothesized that sexual reproduction enables genetic recombination between individuals to generate diversity, thus increasing population fitness. However, a well-known meiotic feature is that recombination is not randomly distributed across the genome: “coldspots” and “hotspots” exist, implying some regions undergo exchange more frequently. In this paper, the authors report the discovery of recombination hotspots linked to the sex-determining region mating-type locus (MAT) in the human pathogenic fungus Cryptococcus neoformans. Through genetic analysis, hotspots were found to reside on both sides of MAT and are associated with DNA regions marked by high G + C base pair composition. Moreover, recombination on one side of MAT is associated with a recombination on the other. As a result, the MAT locus can be replaced onto the homologous chromosome as a unit—an effective switch of MAT. Based on these findings, they propose that the MAT-linked recombination hotspots impacted key steps during MAT evolution. This study has broader implications on how gene clusters (including those involved in metabolism or secondary metabolite production) are assembled and maintained and explains how recombination is distorted in sex-determining regions in eukaryotes.

A new method for detecting human recombination hotspots and its applications to the HapMap ENCODE data

Computational detection of recombination hotspots from population polymorphism data is important both for understanding the nature of recombination and for applications such as association studies. We propose a new method for this task based on a multiple-hotspot model and an (approximate) log-likelihood ratio test. A truncated, weighted pairwise log-likelihood is introduced and applied to the calculation of the log-likelihood ratio, and a forward-selection procedure is adopted to search for the optimal hotspot predictions. The method shows a relatively high power with a low false-positive rate in detecting multiple hotspots in simulation data and has a performance comparable to the best results of leading computational methods in experimental data for which recombination hotspots have been characterized by sperm-typing experiments. The method can be applied to both phased and unphased data directly, with a very fast computational speed. We applied the method to the 10 500-kb regions of the HapMap ENCODE data and found 172 hotspots among the three populations, with average hotspot width of 2.4 kb. By comparisons with the simulation data, we found some evidence that hotspots are not all identical across populations. The correlations between detected hotspots and several genomic characteristics were examined. In particular, we observed that DNaseI-hypersensitive sites are enriched in hotspots, suggesting the existence of human beta hotspots similar to those found in yeast.

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

Background

Meiotic recombination events have been found to concentrate in 1–2.5 kilo base regions, but these recombination hot spots do not share a consensus sequence and why they occur at specific sites is not fully understood. Some previous evidence suggests that poly-purine/poly-pyrimidine (poly-pu/py) tracts (PPTs), a class of sequence with distinctive biochemical properties, could be involved in recombination, but no general association of PPTs with meiotic recombination hot spots has previously been reported.

Results

We used computational methods to investigate in detail the relationship between PPTs and hot spots. We show statistical associations of PPT frequency with hot spots of meiotic recombination initiating lesions, double-strand breaks, in the genome of the yeast S. cerevisiae and with experimentally well characterized human meiotic recombination hot spots. Supporting a possible role of poly-pu/py-rich sequences in hot spot recombination, we also found that all three single nucleotide polymorphisms previously shown to be associated with human hot spot activity changes occur within sequence contexts of 14 bp or longer that are 85% or more poly-pu/py and at least 70% G/C. These polymorphisms are all close to the hot spot mid points. Comparing the sequences of experimentally characterized human hot spots with the orthologous regions of the chimpanzee genome previously shown not to contain hot spots, we found that in all five cases in which comparisons for the hot spot central regions are possible with publicly available sequence data, there are differences near the human hot spot mid points within sequences 14 bp or longer consisting of more than 80% poly-pu/py and at least 50% G/C.

Conclusion

Our results, along with previous evidence for the unique biochemical properties and recombination-stimulating potential of poly-pu/py-rich sequences, suggest that the possible functional involvement of this type of sequence in meiotic recombination hot spots deserves further experimental exploration.