Dataset Information

Population genomic inference of recombination rates and hotspots.

ABSTRACT: As more human genomic data become available, fine-scale recombination rate variation can be inferred on a genome-wide scale. Current statistical methods to infer recombination rates that can be applied to moderate, or large, genomic regions are limited to approximated likelihoods. Here, we develop a Bayesian full-likelihood method using Markov Chain Monte Carlo (MCMC) to estimate background recombination rates and hotspots. The probability model is inspired by the observed patterns of recombination at several genomic regions analyzed in sperm-typing studies. Posterior probabilities and Bayes factors of recombination hotspots along chromosomes are inferred. For moderate-size genomic regions (e.g., with <100 SNPs), the full-likelihood method is used. Larger regions are split into subintervals (typically each having between 20 and 50 markers). The likelihood is approximated based on the genealogies for each subinterval. The background recombination rates, hotspots, and parameters are evaluated by using a parallel computing approach and assuming shared parameters across the subintervals. Simulation analyses show that our method can accurately estimate the variation in recombination rates across genomic regions. In particular, clusters of hotspots can be distinguished even though weaker hotspots are present. The method is applied to SNP data from the HLA region, the MS32, and chromosome 19.

SUBMITTER: Wang Y

PROVIDER: S-EPMC2669376 | biostudies-literature | 2009 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Population genomic inference of recombination rates and hotspots.

Wang Ying Y Rannala Bruce B

Proceedings of the National Academy of Sciences of the United States of America 20090402 15

As more human genomic data become available, fine-scale recombination rate variation can be inferred on a genome-wide scale. Current statistical methods to infer recombination rates that can be applied to moderate, or large, genomic regions are limited to approximated likelihoods. Here, we develop a Bayesian full-likelihood method using Markov Chain Monte Carlo (MCMC) to estimate background recombination rates and hotspots. The probability model is inspired by the observed patterns of recombinat ...[more]

PMID: 19342488

Similar Datasets

Project description:We used a high-density tiling array to estimate genetic recombination rate among 32 independent recombinant progeny of a P. falciparum genetic cross (7G8 M-CM-^W GB4). We detected 3184 segregating multi-probe single-feature polymorphisms (mSFPs) and 638 recombination events (496 excluding those from subtelomeric regions). These data, in combination with results from 254 previously reported microsatellites, enabled us to construct a high-resolution genetic map. Comparing genetic and physical maps, we obtained an overall recombination rate of 9.6 kb/cM (12.8 kb/cM excluding subtelomeric regions) and identified 54 hotspots, some of which occurred in genes encoding surface antigens or proteins with repetitive motifs that might play a role in genetic recombination in the parasite. Motifs enriched in hotspots were also identified. In agreement with results from a previous cross (HB3 M-BM-4 Dd2), there was positive correlation between sizes of individual chromosomes and their recombination events. These results show that the P. falciparum genome is highly recombinogenic, providing an important genetic basis for parasite survival under various selection pressures. GC-rich repetitive motifs identified in the hotspot sequences may play a role in the high recombination frequency observed. Ten microgram of genomic DNA, extracted and purified from 3D7 (reference), thirty-two P. falciparum independent recombinant progeny of the 7G8 x GB4 cross, and the two parental lines (Hayton, 2008), were hybridized to the PFSANGER GenechipM-BM-. (Affymetrix, Inc., Santa Clara, CA, USA). The scanned image CEL files were first processed using the RMA method, then averaged and compared with reference genome 3D7, and lastly assigned either 7G8 or GB4 alleles based on similarities to the two parental lines. Total of 35 genomic DNA samples (biological replicates: 6 for 3D7, 4 for 7G8, 4 for GB4, and 2 for Pf_WE2). The supplementary file 'GSE25656_QuantNormData_Log2_AllSamples.txt' contains the RMA-normalized data for all of the samples. The supplementary files 'GSE25656_chr*' contain the parental allele assignment of each chromosome and include probe-level annotation.

Project description:Genetic recombination is a very important evolutionary mechanism that mixes parental haplotypes and produces new raw material for organismal evolution. As a result, information on recombination rates is critical for biological research. In this paper, we introduce a new extremely fast open-source software package (FastEPRR) that uses machine learning to estimate recombination rate [Formula: see text] (=[Formula: see text]) from intraspecific DNA polymorphism data. When [Formula: see text] and the number of sampled diploid individuals is large enough ([Formula: see text]), the variance of [Formula: see text] remains slightly smaller than that of [Formula: see text] The new estimate [Formula: see text] (calculated by averaging [Formula: see text] and [Formula: see text]) has the smallest variance of all cases. When estimating [Formula: see text], the finite-site model was employed to analyze cases with a high rate of recurrent mutations, and an additional method is proposed to consider the effect of variable recombination rates within windows. Simulations encompassing a wide range of parameters demonstrate that different evolutionary factors, such as demography and selection, may not increase the false positive rate of recombination hotspots. Overall, accuracy of FastEPRR is similar to the well-known method, LDhat, but requires far less computation time. Genetic maps for each human population (YRI, CEU, and CHB) extracted from the 1000 Genomes OMNI data set were obtained in less than 3 d using just a single CPU core. The Pearson Pairwise correlation coefficient between the [Formula: see text] and [Formula: see text] maps is very high, ranging between 0.929 and 0.987 at a 5-Mb scale. Considering that sample sizes for these kinds of data are increasing dramatically with advances in next-generation sequencing technologies, FastEPRR (freely available at http://www.picb.ac.cn/evolgen/) is expected to become a widely used tool for establishing genetic maps and studying recombination hotspots in the population genomic era.

Dataset Information

Population genomic inference of recombination rates and hotspots.

Publications

Population genomic inference of recombination rates and hotspots.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets