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.     

Critical assessment of coalescent simulators in modeling recombination hotspots in genomic sequences

Background

Coalescent simulation is pivotal for understanding population evolutionary models and demographic histories, as well as for developing novel analytical methods for genetic association studies for DNA sequence data. A plethora of coalescent simulators are developed, but selecting the most appropriate program remains challenging.

Results

We extensively compared performances of five widely used coalescent simulators – Hudson’s ms, msHOT, MaCS, Simcoal2, and fastsimcoal, to provide a practical guide considering three crucial factors, 1) speed, 2) scalability and 3) recombination hotspot position and intensity accuracy. Although ms represents a popular standard coalescent simulator, it lacks the ability to simulate sequences with recombination hotspots. An extended program msHOT has compensated for the deficiency of ms by incorporating recombination hotspots and gene conversion events at arbitrarily chosen locations and intensities, but remains limited in simulating long stretches of DNA sequences. Simcoal2, based on a discrete generation-by-generation approach, could simulate more complex demographic scenarios, but runs comparatively slow. MaCS and fastsimcoal, both built on fast, modified sequential Markov coalescent algorithms to approximate standard coalescent, are much more efficient whilst keeping salient features of msHOT and Simcoal2, respectively. Our simulations demonstrate that they are more advantageous over other programs for a spectrum of evolutionary models. To validate recombination hotspots, LDhat 2.2 rhomap package, sequenceLDhot and Haploview were compared for hotspot detection, and sequenceLDhot exhibited the best performance based on both real and simulated data.

Conclusions

While ms remains an excellent choice for general coalescent simulations of DNA sequences, MaCS and fastsimcoal are much more scalable and flexible in simulating a variety of demographic events under different recombination hotspot models. Furthermore, sequenceLDhot appears to give the most optimal performance in detecting and validating cross-over hotspots.     

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.     

Gene conversion homogenizes the CMT1A paralogous repeats

Background

Non-allelic homologous recombination between paralogous repeats is increasingly being recognized as a major mechanism causing both pathogenic microdeletions and duplications, and structural polymorphism in the human genome. It has recently been shown empirically that gene conversion can homogenize such repeats, resulting in longer stretches of absolute identity that may increase the rate of non-allelic homologous recombination.

Results

Here, a statistical test to detect gene conversion between pairs of non-coding sequences is presented. It is shown that the 24 kb Charcot-Marie-Tooth type 1A paralogous repeats (CMT1A-REPs) exhibit the imprint of gene conversion processes whilst control orthologous sequences do not. In addition, Monte Carlo simulations of the evolutionary divergence of the CMT1A-REPs, incorporating two alternative models for gene conversion, generate repeats that are statistically indistinguishable from the observed repeats. Bounds are placed on the rate of these conversion processes, with central values of 1.3 × 10^-4 and 5.1 × 10^-5 per generation for the alternative models.

Conclusions

This evidence presented here suggests that gene conversion may have played an important role in the evolution of the CMT1A-REP paralogous repeats. The rates of these processes are such that it is probable that homogenized CMT1A-REPs are polymorphic within modern populations. Gene conversion processes are similarly likely to play an important role in the evolution of other segmental duplications and may influence the rate of non-allelic homologous recombination between them.

     

The Red Queen theory of recombination hotspots

Recombination hotspots are small chromosomal regions, where meiotic crossover events happen with high frequency. Recombination is initiated by a double‐strand break (DSB) that requires the intervention of the molecular repair mechanism. The DSB repair mechanism may result in the exchange of homologous chromosomes (crossover) and the conversion of the allelic sequence that breaks into the one that does not break (biased gene conversion). Biased gene conversion results in a transmission advantage for the allele that does not break, thus preventing recombination and rendering recombination hotspots transient. How is it possible that recombination hotspots persist over evolutionary time (maintaining the average chromosomal crossover rate) when they are self‐destructive? This fundamental question is known as the recombination hotspot paradox and has attracted much attention in recent years. Yet, that attention has not translated into a fully satisfactory answer. No existing model adequately explains all aspects of the recombination hotspot paradox. Here, we formulate an intragenomic conflict model resulting in Red Queen dynamics that fully accounts for all empirical observations regarding the molecular mechanisms of recombination hotspots, the nonrandom targeting of the recombination machinery to hotspots and the evolutionary dynamics of hotspot turnover.     

Support vector machine for classification of meiotic recombination hotspots and coldspots in Saccharomyces cerevisiae based on codon composition

Background  

Meiotic double-strand breaks occur at relatively high frequencies in some genomic regions (hotspots) and relatively low frequencies in others (coldspots). Hotspots and coldspots are receiving increasing attention in research into the mechanism of meiotic recombination. However, predicting hotspots and coldspots from DNA sequence information is still a challenging task.     

Drosophila duplication hotspots are associated with late-replicating regions of the genome

Duplications play a significant role in both extremes of the phenotypic spectrum of newly arising mutations: they can have severe deleterious effects (e.g. duplications underlie a variety of diseases) but can also be highly advantageous. The phenotypic potential of newly arisen duplications has stimulated wide interest in both the mutational and selective processes shaping these variants in the genome. Here we take advantage of the Drosophila simulans-Drosophila melanogaster genetic system to further our understanding of both processes. Regarding mutational processes, the study of two closely related species allows investigation of the potential existence of shared duplication hotspots, and the similarities and differences between the two genomes can be used to dissect its underlying causes. Regarding selection, the difference in the effective population size between the two species can be leveraged to ask questions about the strength of selection acting on different classes of duplications. In this study, we conducted a survey of duplication polymorphisms in 14 different lines of D. simulans using tiling microarrays and combined it with an analogous survey for the D. melanogaster genome. By integrating the two datasets, we identified duplication hotspots conserved between the two species. However, unlike the duplication hotspots identified in mammalian genomes, Drosophila duplication hotspots are not associated with sequences of high sequence identity capable of mediating non-allelic homologous recombination. Instead, Drosophila duplication hotspots are associated with late-replicating regions of the genome, suggesting a link between DNA replication and duplication rates. We also found evidence supporting a higher effectiveness of selection on duplications in D. simulans than in D. melanogaster. This is also true for duplications segregating at high frequency, where we find evidence in D. simulans that a sizeable fraction of these mutations is being driven to fixation by positive selection.     

Nucleosome occupancy landscape and dynamics at mouse recombination hotspots

During meiosis, paternal and maternal homologous chromosomes recombine at specific recombination sites named hotspots. What renders 2% of the mammalian genomes permissive to meiotic recombination by allowing Spo11 endonuclease to initiate double‐strand breaks is largely unknown. Work in yeast has shown that chromatin accessibility seems to be important for this activity. Here, we define nucleosome profiles and dynamics at four mouse recombination hotspots by purifying highly enriched fractions of meiotic cells. We found that nucleosome occupancy is generally stable during meiosis progression. Interestingly, the cores of recombination hotspots have largely open chromatin structure, and the localization of the few nucleosomes present in these cores correlates precisely with the crossover‐free zones in recombinogenic domains. Collectively, these high‐resolution studies suggest that nucleosome occupancy seems to direct, at least in part, how meiotic recombination events are processed.     

Transmission of infectious diseases en route to habitat hotspots

Background

The spread of infectious diseases in wildlife populations is influenced by patterns of between-host contacts. Habitat “hotspots” - places attracting a large numbers of individuals or social groups - can significantly alter contact patterns and, hence, disease propagation. Research on the importance of habitat hotspots in wildlife epidemiology has primarily focused on how inter-individual contacts occurring at the hotspot itself increase disease transmission. However, in territorial animals, epidemiologically important contacts may primarily occur as animals cross through territories of conspecifics en route to habitat hotspots. So far, the phenomenon has received little attention. Here, we investigate the importance of these contacts in the case where infectious individuals keep visiting the hotspots and in the case where these individuals are not able to travel to the hotspot any more.

Methodology and Principal Findings

We developed a simulation epidemiological model to investigate both cases in a scenario when transmission at the hotspot does not occur. We find that (i) hotspots still exacerbate epidemics, (ii) when infectious individuals do not travel to the hotspot, the most vulnerable individuals are those residing at intermediate distances from the hotspot rather than nearby, and (iii) the epidemiological vulnerability of a population is the highest when the number of hotspots is intermediate.

Conclusions and Significance

By altering animal movements in their vicinity, habitat hotspots can thus strongly increase the spread of infectious diseases, even when disease transmission does not occur at the hotspot itself. Interestingly, when animals only visit the nearest hotspot, creating additional artificial hotspots, rather than reducing their number, may be an efficient disease control measure.     

Transboundary conservation hotspots in China and potential impacts of the belt and road initiative

Aim

Biodiversity hotspots often span international borders, thus conservation efforts must as well. China is one of the most biodiverse countries and the length of its international land borders is the longest in the world; thus, there is a strong need for transboundary conservation. We identify China's transboundary conservation hotspots and analyse the potential effects of the Belt and Road Initiative (BRI) on them to provide recommendations for conservation actions.

Location

China, Asia.

Methods

We compiled a species list of terrestrial vertebrates that span China's borders. Using their distribution, we extracted the top 30% of the area with the highest richness value weighted by Red List category and considered these transboundary hotspots for conservation priority. Then we analysed protected area (PA) coverage and connectivity to identify conservation gaps. To measure potential impact of the BRI, we counted the species whose distribution range is traversed by the BRI, and calculated the aggregation index, proportion of natural land and night light index along its routes.

Results

We identified 1964 terrestrial vertebrate species living in the border region. We identified four transboundary hotspots and found insufficient PA coverage and low connectivity in three of them. The BRI routes intersected all four hotspots and traversed 82.4% (1619/1964) of the transboundary species, half of which (918) are sensitive to the potential risks brought by the BRI. Night light index increased generally along the BRI. However, the proportion of natural land and the aggregation index near the BRI showed different trends in hotspots.

Main Conclusions

There is an urgent need for conservation action in China's transboundary region. The BRI should put biodiversity conservation at the core of its development strategy. Furthermore, we suggest using the planned BRI as a platform for dialogue and consultation, knowledge and data sharing, and joint planning to promote transboundary conservation.     

Meiotic homoeologous recombination-based mapping of wheat chromosome 2B and its homoeologues in <Emphasis Type="Italic">Aegilops speltoides</Emphasis> and <Emphasis Type="Italic">Thinopyrum elongatum</Emphasis>

Key message

We physically dissected and mapped wheat chromosome 2B and its homoeologues in Aegilops speltoides and Thinopyrum elongatum based on meiotic homoeologous recombination, providing a unique physical framework for genome studies.

Abstract

Common wheat has a large and complex genome with narrow genetic diversity and various degrees of recombination between the A, B, and D subgenomes. This has limited the homologous recombination-based genome studies in wheat. Here, we exploited meiotic homoeologous recombination for molecular mapping of wheat chromosome 2B and its homoeologue 2S from Aegilops speltoides and 2E from Thinopyrum elongatum. The 2B–2S and 2B–2E recombination was induced by the ph1b mutant, and recovered using molecular markers and fluorescent genomic in situ hybridization (FGISH). A total of 112 2B–2S and 87 2B–2E recombinants involving different chromosome regions were developed and physically delineated by FGISH. The 2B–2S and 2B–2E recombination hotspots mapped to the subterminal regions on both arms. Recombination hotspots with the highest recombination rates mapped to the short arms. Eighty-three 2B–2S and 67 2B–2E recombinants were genotyped using the wheat 90 K SNP arrays. Based on the genotyping results and FGISH patterns of the recombinants, chromosomes 2B, 2S, and 2E were partitioned into 93, 66, and 46 bins, respectively. In total, 1037 SNPs physically mapped onto distinct bins of these three homoeologous chromosomes. A homoeologous recombination-based bin map was constructed for chromosome 2B, providing a unique physical framework for genome studies in wheat and its relatives. Meiotic homoeologous recombination also facilitates gene introgression to diversify the wheat genome for germplasm development. Therefore, homoeologous recombination-based studies enhance understanding of the wheat genome and its homoeologous counterparts from wild grasses, and expand the genetic variability of the wheat genome.

     

Dengue virus serotype infection specifies the activation of the unfolded protein response

Background

Tomato-infecting begomoviruses are widely distributed across the world and cause diseases of high economic impact on wide range of agriculturally important crops. Though recombination plays a pivotal role in diversification and evolution of these viruses, it is currently unknown whether there are differences in the number and quality of recombination events amongst different tomato-infecting begomovirus species. To examine this we sought to characterize the recombination events, estimate the frequency of recombination, and map recombination hotspots in tomato-infecting begomoviruses of South and Southeast Asia.

Results

Different methods used for recombination breakpoint analysis provided strong evidence for presence of recombination events in majority of the sequences analyzed. However, there was a clear evidence for absence or low Recombination events in viruses reported from North India. In addition, we provide evidence for non-random distribution of recombination events with the highest frequency of recombination being mapped in the portion of the N-terminal portion of Rep.

Conclusion

The variable recombination observed in these viruses signified that all begomoviruses are not equally prone to recombination. Distribution of recombination hotspots was found to be reliant on the relatedness of the genomic region involved in the exchange. Overall the frequency of phylogenetic violations and number of recombination events decreased with increasing parental sequence diversity. These findings provide valuable new information for understanding the diversity and evolution of tomato-infecting begomoviruses in Asia.     

Patterns of species richness hotspots and estimates of their protection are sensitive to spatial resolution

Aim

Species richness is a measure of biodiversity often used in spatial conservation assessments and mapped by summing species distribution maps. Commission errors inherent those maps influence richness patterns and conservation assessments. We sought to further the understanding of the sensitivity of hotspot delineation methods and conservation assessments to commission errors, and choice of threshold for hotspot delineation.

Location

United States.

Methods

We created range maps and 30‐m and 1‐km resolution habitat maps for terrestrial vertebrates in the United States and generated species richness maps with each dataset. With the richness maps and the GAP Protected Areas Dataset, we created species richness hotspot maps and calculated the proportion of hotspots within protected areas; calculating protection under a range of thresholds for defining hotspots. Our method allowed us to identify the influence of commission errors by comparing hotspot maps.

Results

Commission errors from coarse spatial grain data and lack of porosity in the range data inflated richness estimates and altered their spatial patterns. Coincidence of hotspots from different data types was low. The 30‐m hotspots were spatially dispersed, and some were very long distances from the hotspots mapped with coarser data. Estimates of protection were low for each of the taxa. The relationship between estimates of hotspot protection and threshold choice was nonlinear and inconsistent among data types (habitat and range) and grain size (30‐m and 1‐km).

Main conclusions

Coarse mapping methods and grain sizes can introduce commission errors into species distribution data that could result in misidentifications of the regions where hotspots occur and affect estimates of hotspot protection. Hotspot conservation assessments are also sensitive to choice of threshold for hotspot delineation. There is value in developing species distribution maps with high resolution and low rates of commission error for conservation assessments.     

SequenceLDhot: detecting recombination hotspots

MOTIVATION: There is much local variation in recombination rates across the human genome--with the majority of recombination occurring in recombination hotspots--short regions of around approximately 2 kb in length that have much higher recombination rates than neighbouring regions. Knowledge of this local variation is important, e.g. in the design and analysis of association studies for disease genes. Population genetic data, such as that generated by the HapMap project, can be used to infer the location of these hotspots. We present a new, efficient and powerful method for detecting recombination hotspots from population data. RESULTS: We compare our method with four current methods for detecting hotspots. It is orders of magnitude quicker, and has greater power, than two related approaches. It appears to be more powerful than HotspotFisher, though less accurate at inferring the precise positions of the hotspot. It was also more powerful than LDhot in some situations: particularly for weaker hotspots (10-40 times the background rate) when SNP density is lower (< 1/kb). AVAILABILITY: Program, data sets, and full details of results are available at: http://www.maths.lancs.ac.uk/~fearnhea/Hotspot.     

No dosage effect of recombinational hotspots in the mouse major histocompatibility complex

The sites of meiotic recombination in the proximal region of the mouse major histocompatibility complex (MHC) are clustered at hotspots. Some MHC haplotypes derived from Asian wild mice increase the frequency of recombination at such hotspots when heterozygous with standard laboratory haplotypes. The wm7 and cas3 haplotypes, have a hotspot close to the Lmp-2 gene (Lmp-2 hotspot), and the cas4 haplotype has a hotspot about 100 kilobase (kb) proximal, close to the Pb gene (Pb hotspot). To examine the effect of a double dose of hotspots, we estimated the rate of recombination and determined the location of the breakpoints in crosses of wm7/cas3 and wm7/cas4. In 3570 backcross progeny we identified 29 new recombinants in the H-2K to Ab interval, at a frequency of 0.81%. This frequency is 40-fold higher than in crosses between laboratory haplotypes and very similar to those previously obtained in crosses between these wild and standard laboratory haplotypes. Thus, a double dose of hotspots has no additive effect on the frequency of meiotic recombination. The site-specificity of recombination was also conserved. Twenty-three breakpoints were confined within 5.4 kb in the Lmp-2 hotspot, and six breakpoints from the cas4 cross were located in the Pb hotspot, which we have now confined to a 15 kb segment. Correspondence to: T. Shiroishi.     

A coalescent model of recombination hotspots

Recent experimental findings suggest that the assumption of a homogeneous recombination rate along the human genome is too naive. These findings point to block-structured recombination rates; certain regions (called hotspots) are more prone than other regions to recombination. In this report a coalescent model incorporating hotspot or block-structured recombination is developed and investigated analytically as well as by simulation. Our main results can be summarized as follows: (1) The expected number of recombination events is much lower in a model with pure hotspot recombination than in a model with pure homogeneous recombination, (2) hotspots give rise to large variation in recombination rates along the genome as well as in the number of historical recombination events, and (3) the size of a (nonrecombining) block in the hotspot model is likely to be overestimated grossly when estimated from SNP data. The results are discussed with reference to the current debate about block-structured recombination and, in addition, the results are compared to genome-wide variation in recombination rates. A number of new analytical results about the model are derived.     

Population demographic history can cause the appearance of recombination hotspots

Although the prevailing view among geneticists suggests that recombination hotspots exist ubiquitously across the human genome, there is only limited experimental evidence from a few genomic regions to support the generality of this claim. A small number of true recombination hotspots are well supported experimentally, but the vast majority of hotspots have been identified on the basis of population genetic inferences from the patterns of linkage disequilibrium (LD) seen in the human population. These inferences are made assuming a particular model of human history, and one of the assumptions of that model is that the effective population size of humans has remained constant throughout our history. Our results show that relaxation of the constant population size assumption can create LD and variation patterns that are qualitatively and quantitatively similar to human populations without any need to invoke localized hotspots of recombination. In other words, apparent recombination hotspots could be an artifact of variable population size over time. Several lines of evidence suggest that the vast majority of hotspots identified on the basis of LD information are unlikely to have elevated recombination rates.     

Mutational hotspots in the TP53 gene and, possibly, other tumor suppressors evolve by positive selection

Background  

The mutation spectra of the TP53 gene and other tumor suppressors contain multiple hotspots, i.e., sites of non-random, frequent mutation in tumors and/or the germline. The origin of the hotspots remains unclear, the general view being that they represent highly mutable nucleotide contexts which likely reflect effects of different endogenous and exogenous factors shaping the mutation process in specific tissues. The origin of hotspots is of major importance because it has been suggested that mutable contexts could be used to infer mechanisms of mutagenesis contributing to tumorigenesis.