LDJump: Estimating variable recombination rates from population genetic data.
Ontology highlight
ABSTRACT:
SUBMITTER: Hermann P
PROVIDER: S-EPMC6519033 | biostudies-literature | 2019 May
REPOSITORIES: biostudies-literature
Publications
LDJump: Estimating variable recombination rates from population genetic data.
Hermann Philipp P Heissl Angelika A Tiemann-Boege Irene I Futschik Andreas A
Molecular ecology resources 20190404 3
As recombination plays an important role in evolution, its estimation and the identification of hotspot positions is of considerable interest. We propose a novel approach for estimating population recombination rates based on genotyping or sequence data that involves a sequential multiscale change point estimator. Our method also permits demography to be taken into account. It uses several summary statistics within a regression model fitted on suitable scenarios. Our proposed method is accurate, ...[more]
Similar Datasets
Project description:BackgroundOur current understanding of archaic admixture in humans relies on statistical methods with large biases, whose magnitudes depend on the sizes and separation times of ancestral populations. To avoid these biases, it is necessary to estimate these parameters simultaneously with those describing admixture. Genetic estimates of population histories also confront problems of statistical identifiability: different models or different combinations of parameter values may fit the data equally well. To deal with this problem, we need methods of model selection and model averaging, which are lacking from most existing software.ResultsThe Legofit software package allows simultaneous estimation of parameters describing admixture, and the sizes and separation times of ancestral populations. It includes facilities for data manipulation, estimation, analysis of residuals, model selection, and model averaging.ConclusionsLegofit uses genetic data to study the history of a subdivided population. It is unaffected by recent history and can therefore focus on the deep history of population size, subdivision, and admixture. It outperforms several statistical methods that have been widely used to study population history and should be useful in any species for which DNA sequence data is available from several populations.
| S-EPMC6819480 | biostudies-literature
Project description:As the frequency of antifungal drug resistance continues to increase, understanding the genetic structure of fungal populations, where resistant isolates have emerged and spread, is of major importance. Aspergillus fumigatus is an ubiquitously distributed fungus and the primary causative agent of invasive aspergillosis (IA), a potentially lethal infection in immunocompromised individuals. In the last few years, an increasing number of A. fumigatus isolates has evolved resistance to triazoles, the primary drugs for treating IA infections. In most isolates, this multiple-triazole-resistance (MTR) phenotype is caused by mutations in the cyp51A gene, which encodes the protein targeted by the triazoles. We investigated the genetic differentiation and reproductive mode of A. fumigatus in the Netherlands, the country where the MTR phenotype probably originated, to determine their role in facilitating the emergence and distribution of resistance genotypes. Using 20 genome-wide neutral markers, we genotyped 255 Dutch isolates including 25 isolates with the MTR phenotype. In contrast to previous reports, our results show that Dutch A. fumigatus genotypes are genetically differentiated into five distinct populations. Four of the five populations show significant linkage disequilibrium, indicative of an asexual reproductive mode, whereas the fifth population is in linkage equilibrium, indicative of a sexual reproductive mode. Notably, the observed genetic differentiation among Dutch isolates does not correlate with geography, although all isolates with the MTR phenotype nest within a single, predominantly asexual, population. These results suggest that both reproductive mode and genetic differentiation contribute to the structure of Dutch A. fumigatus populations and are probably shaping the evolutionary dynamics of drug resistance in this potentially deadly pathogen.
| S-EPMC3258581 | biostudies-literature
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.
2011-03-16 | E-GEOD-25656 | biostudies-arrayexpress
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 × 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 ´ 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.
2011-03-16 | GSE25656 | GEO
Project description: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.
| S-EPMC2669376 | biostudies-literature
Project description:Copy-number-variable (CNV) loci differ from single nucleotide polymorphic (SNP) sites in size, mutation rate, and mechanisms of maintenance in natural populations. It is therefore hypothesized that population genetic divergence at CNV loci will differ from that found at SNP sites. Here, we test this hypothesis by analysing 856 CNV loci from the genomes of 1184 healthy individuals from 11 HapMap populations with a wide range of ancestry. The results show that population genetic divergence at the CNV loci is generally more than three times lower than at genome-wide SNP sites. Populations generally exhibit very small genetic divergence (Gst?=?0.05?±?0.049). The smallest divergence is among African populations (Gst?=?0.0081?±?0.0025), with increased divergence among non-African populations (Gst?=?0.0217?±?0.0109) and then among African and non-African populations (Gst?=?0.0324?±?0.0064). Genetic diversity is high in African populations (~0.13), low in Asian populations (~0.11), and intermediate in the remaining 11 populations. Few significant linkage disequilibria (LDs) occur between the genome-wide CNV loci. Patterns of gametic and zygotic LDs indicate the absence of epistasis among CNV loci. Mutation rate is about twice as large as the migration rate in the non-African populations, suggesting that the high mutation rates play dominant roles in producing the low population genetic divergence at CNV loci.
| S-EPMC5320550 | biostudies-literature
Project description:The knowledge of the frequency and relative weight of mutation and recombination events in evolution is essential for understanding how microorganisms reach fitted phenotypes. Traditionally, these evolutionary parameters have been inferred by using data from multilocus sequence typing (MLST), which is known to have yielded conflicting results. In the near future, these estimations will certainly be performed by computational analyses of full-genome sequences. However, it is not known whether this approach will yield accurate results as bacterial genomes exhibit heterogeneous representation of loci categories, and it is not clear how loci nature impacts such estimations. Therefore, we assessed how mutation and recombination inferences are shaped by loci with different genetic features, using the bacterium Chlamydia trachomatis as the study model. We found that loci assigning a high number of alleles and positively selected genes yielded nonconvergent estimates and incongruent phylogenies and thus are more prone to confound algorithms. Unexpectedly, for the model under evaluation, housekeeping genes and noncoding regions shaped estimations in a similar manner, which points to a nonrandom role of the latter in C. trachomatis evolution. Although the present results relate to a specific bacterium, we speculate that microbe-specific genomic architectures (such as coding capacity, polymorphism dispersion, and fraction of positively selected loci) may differentially buffer the effect of the confounding factors when estimating recombination and mutation rates and, thus, influence the accuracy of using full-genome sequences for such purpose. This putative bias associated with in silico inferences should be taken into account when discussing the results obtained by the analyses of full-genome sequences, in which the "one size fits all" approach may not be applicable.
| S-EPMC3385982 | biostudies-literature
Project description:Accurate estimation of recombination rates is critical for studying the origins and maintenance of genetic diversity. Because the inference of recombination rates under a full evolutionary model is computationally expensive, we developed an alternative approach using topological data analysis (TDA) on genome sequences. We find that this method can analyze datasets larger than what can be handled by any existing recombination inference software, and has accuracy comparable to commonly used model-based methods with significantly less processing time. Previous TDA methods used information contained solely in the first Betti number ([Formula: see text]) of a set of genomes, which aims to capture the number of loops that can be detected within a genealogy. These explorations have proven difficult to connect to the theory of the underlying biological process of recombination, and, consequently, have unpredictable behavior under perturbations of the data. We introduce a new topological feature, which we call ?, with a natural connection to coalescent models, and present novel arguments relating [Formula: see text] to population genetic models. Using simulations, we show that ? and [Formula: see text] are differentially affected by missing data, and package our approach as TREE (Topological Recombination Estimator). TREE's efficiency and accuracy make it well suited as a first-pass estimator of recombination rate heterogeneity or hotspots throughout the genome. Our work empirically and theoretically justifies the use of topological statistics as summaries of genome sequences and describes a new, unintuitive relationship between topological features of the distribution of sequence data and the footprint of recombination on genomes.
| S-EPMC6456321 | biostudies-literature
Project description:Despite the fundamental role of crossing-over in the pairing and segregation of chromosomes during human meiosis, the rates and placements of events vary markedly among individuals. Characterizing this variation and identifying its determinants are essential steps in our understanding of the human recombination process and its evolution.Using three large sets of European-American pedigrees, we examined variation in five recombination phenotypes that capture distinct aspects of crossing-over patterns. We found that the mean recombination rate in males and females and the historical hotspot usage are significantly heritable and are uncorrelated with one another. We then conducted a genome-wide association study in order to identify loci that influence them. We replicated associations of RNF212 with the mean rate in males and in females as well as the association of Inversion 17q21.31 with the female mean rate. We also replicated the association of PRDM9 with historical hotspot usage, finding that it explains most of the genetic variance in this phenotype. In addition, we identified a set of new candidate regions for further validation.These findings suggest that variation at broad and fine scales is largely separable and that, beyond three known loci, there is no evidence for common variation with large effects on recombination phenotypes.
| S-EPMC3117798 | biostudies-literature
Project description:Exome and whole-genome sequencing studies are becoming increasingly common, but little is known about the accuracy of the genotype calls made by the commonly used platforms. Here we use replicate high-coverage sequencing of blood and saliva DNA samples from four European-American individuals to estimate lower bounds on the error rates of Complete Genomics and Illumina HiSeq whole-genome and whole-exome sequencing. Error rates for nonreference genotype calls range from 0.1% to 0.6%, depending on the platform and the depth of coverage. Additionally, we found (1) no difference in the error profiles or rates between blood and saliva samples; (2) Complete Genomics sequences had substantially higher error rates than Illumina sequences had; (3) error rates were higher (up to 6%) for rare or unique variants; (4) error rates generally declined with genotype quality (GQ) score, but in a nonlinear fashion for the Illumina data, likely due to loss of specificity of GQ scores greater than 60; and (5) error rates increased with increasing depth of coverage for the Illumina data. These findings, especially (3)-(5), suggest that caution should be taken in interpreting the results of next-generation sequencing-based association studies, and even more so in clinical application of this technology in the absence of validation by other more robust sequencing or genotyping methods.
| S-EPMC4216915 | biostudies-literature