A sequential coalescent algorithm for chromosomal inversions.
Ontology highlight
ABSTRACT:
SUBMITTER: Peischl S
PROVIDER: S-EPMC3746820 | biostudies-other | 2013 Sep
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Bacteria can exchange and acquire new genetic material from other organisms directly and via the environment. This process, known as bacterial recombination, has a strong impact on the evolution of bacteria, for example, leading to the spread of antibiotic resistance across clades and species, and to the avoidance of clonal interference. Recombination hinders phylogenetic and transmission inference because it creates patterns of substitutions (homoplasies) inconsistent with the hypothesis of a single evolutionary tree. Bacterial recombination is typically modeled as statistically akin to gene conversion in eukaryotes, i.e., using the coalescent with gene conversion (CGC). However, this model can be very computationally demanding as it needs to account for the correlations of evolutionary histories of even distant loci. So, with the increasing popularity of whole genome sequencing, the need has emerged for a faster approach to model and simulate bacterial genome evolution. We present a new model that approximates the coalescent with gene conversion: the bacterial sequential Markov coalescent (BSMC). Our approach is based on a similar idea to the sequential Markov coalescent (SMC)-an approximation of the coalescent with crossover recombination. However, bacterial recombination poses hurdles to a sequential Markov approximation, as it leads to strong correlations and linkage disequilibrium across very distant sites in the genome. Our BSMC overcomes these difficulties, and shows a considerable reduction in computational demand compared to the exact CGC, and very similar patterns in simulated data. We implemented our BSMC model within new simulation software FastSimBac. In addition to the decreased computational demand compared to previous bacterial genome evolution simulators, FastSimBac provides more general options for evolutionary scenarios, allowing population structure with migration, speciation, population size changes, and recombination hotspots. FastSimBac is available from https://bitbucket.org/nicofmay/fastsimbac, and is distributed as open source under the terms of the GNU General Public License. Lastly, we use the BSMC within an Approximate Bayesian Computation (ABC) inference scheme, and suggest that parameters simulated under the exact CGC can correctly be recovered, further showcasing the accuracy of the BSMC. With this ABC we infer recombination rate, mutation rate, and recombination tract length of Bacillus cereus from a whole genome alignment.
| S-EPMC5419479 | biostudies-literature
Project description:Inversions are an important form of structural variation, but they are difficult to characterize, as their breakpoints often fall within inverted repeats. We have developed a method called 'haplotype fusion' in which an inversion breakpoint is genotyped by performing fusion PCR on single molecules of human genomic DNA. Fusing single-copy sequences bracketing an inversion breakpoint generates orientation-specific PCR products, exemplified by a genotyping assay for the int22 hemophilia A inversion on Xq28. Furthermore, we demonstrated that inversion events with breakpoints embedded within long (>100 kb) inverted repeats can be genotyped by haplotype-fusion PCR followed by bead-based single-molecule haplotyping on repeat-specific markers bracketing the inversion breakpoint. We illustrate this method by genotyping a Yp paracentric inversion sponsored by >300-kb-long inverted repeats. The generality of our methods to survey for, and genotype chromosomal inversions should help our understanding of the contribution of inversions to genomic variation, inherited diseases and cancer.
| S-EPMC2690135 | biostudies-literature
Project description:Chromosomal rearrangements occur readily in nature and are a major reshaping force during genome evolution. Such large scale modifications are usually deleterious causing several fitness defects, but sometimes can confer an advantage and become adaptive. For example the DAL metabolic cluster in yeast was assembled in recent evolutionary times in the Hemiascomycetes lineage, through a set of rearrangements that brought together the genes involved in the allantoin degradation pathway. In eukaryotes, the existence of physical clustering of genes with related functions supports the notion that neighbouring ORFs tend to be co-expressed and that the order of genes along the chromosomes may have biological significance, rather than being random as previously believed. In this study, we investigate the phenotypic effect that inversions have on the DAL gene cluster, expressed during nitrogen starvation. In all Saccharomyces "sensu stricto" species the order of the DAL cluster is conserved, while in the "sensu lato" species Naumovia castellii, which grows significantly worse than S. cerevisiae on allantoin, the cluster includes two nested inversions encompassing three DAL genes. We constructed several inverted and non-inverted S. cerevisiae strains possessing different inversions including those to mimic the configuration of the N. castellii DAL cluster. We showed that the inversion of DAL2 lower its own expression and reduces yeast fitness during nitrogen starvation. This rearrangement also altered the expression of the neighbouring genes DAL1 and DAL4. Moreover, we showed that the expression of the DAL4 anti-sense transcript (SUT614) does not change upon inversions of DAL2 and therefore is unlikely to be involved in its regulation. These results show that the order of the DAL cluster has an impact on the phenotype and gene expression, suggesting that these rearrangements may have been adaptive in the "sensu stricto" group in relation to the low availability of nitrogen in the environment.
| S-EPMC3419248 | biostudies-literature
Project description:The long interspersed element-1 (LINE-1 or L1) and Alu elements are the most abundant mobile elements comprising 21% and 11% of the human genome, respectively. Since the divergence of human and chimpanzee lineages, these elements have vigorously created chromosomal rearrangements causing genomic difference between humans and chimpanzees by either increasing or decreasing the size of genome. Here, we report an exotic mechanism, retrotransposon recombination-mediated inversion (RRMI), that usually does not alter the amount of genomic material present. Through the comparison of the human and chimpanzee draft genome sequences, we identified 252 inversions whose respective inversion junctions can clearly be characterized. Our results suggest that L1 and Alu elements cause chromosomal inversions by either forming a secondary structure or providing a fragile site for double-strand breaks. The detailed analysis of the inversion breakpoints showed that L1 and Alu elements are responsible for at least 44% of the 252 inversion loci between human and chimpanzee lineages, including 49 RRMI loci. Among them, three RRMI loci inverted exonic regions in known genes, which implicates this mechanism in generating the genomic and phenotypic differences between human and chimpanzee lineages. This study is the first comprehensive analysis of mobile element bases inversion breakpoints between human and chimpanzee lineages, and highlights their role in primate genome evolution.
| S-EPMC2603318 | biostudies-literature
Project description:Understanding the role of chromosomal inversions in speciation is a fundamental problem in evolutionary genetics. Here, we perform a comprehensive reconstruction of the evolutionary histories of the chromosomal inversions in Drosophila persimilis and D. pseudoobscura. We provide a solution to the puzzling origins of the selfish Sex-Ratio arrangement in D. persimilis and uncover surprising patterns of phylogenetic discordance on this chromosome. These patterns show that, contrary to widely held views, all fixed chromosomal inversions between D. persimilis and D. pseudoobscura were already present in their ancestral population long before the species split. Our results suggest that patterns of higher genomic divergence and an association of reproductive isolation genes with chromosomal inversions may be a direct consequence of incomplete lineage sorting of ancestral polymorphisms. These findings force a reconsideration of the role of chromosomal inversions in speciation, not as protectors of existing hybrid incompatibilities, but as fertile grounds for their formation.
| S-EPMC6085072 | biostudies-literature
Project description:We have used the inversion system of Drosophila pseudoobscura to investigate how genetic flux occurs among the gene arrangements. The patterns of nucleotide polymorphism at seven loci were used to infer gene conversion events between pairs of different gene arrangements. We estimate that the average gene conversion tract length is 205 bp and that the average conversion rate is 3.4 x 10(-6), which is 2 orders of magnitude greater than the mutation rate. We did not detect gene conversion events between all combinations of gene arrangements even though there was sufficient nucleotide variation for detection and sufficient opportunity for exchanges to occur. Genetic flux across the inverted chromosome resulted in higher levels of differentiation within 0.1 Mb of inversion breakpoints, but a slightly lower level of differentiation in central inverted regions. No gene conversion events were detected within 17 kb of an inversion breakpoint suggesting that the formation of double-strand breaks is reduced near rearrangement breakpoints in heterozygotes. At least one case where selection rather than proximity to an inversion breakpoint is responsible for reduction in polymorphism was identified.
| S-EPMC1456099 | biostudies-literature
Project description:There is increasing evidence regarding the role of chromosomal inversions in relevant biological processes such as local adaptation and speciation. A classic example of the adaptive role of chromosomal polymorphisms is given by the clines of inversion frequencies in Drosophila subobscura, repeatable across continents. Nevertheless, not much is known about the molecular variation associated with these polymorphisms. We characterized the genetic content of ca. 600 individuals from nine European populations following a latitudinal gradient by analysing 19 microsatellite loci from two autosomes (J and U) and the sex chromosome (A), taking into account their chromosomal inversions. Our results clearly demonstrate the molecular genetic uniformity within a given chromosomal inversion across a large latitudinal gradient, particularly from Groningen (Netherlands) in the north to Málaga (Spain) in the south, experiencing highly diverse environmental conditions. This low genetic differentiation within the same gene arrangement across the nine European populations is consistent with the local adaptation hypothesis for th evolutionof chromosomal polymorphisms. We also show the effective role of chromosomal inversions in maintaining different genetic pools within these inverted genomic regions even in the presence of high gene flow. Inversions represent thus an important barrier to gene flux and can help maintain specific allelic combinations with positive effects on fitness. Consistent patterns of microsatellite allele-inversion linkage disequilibrium particularly in loci within inversions were also observed. Finally, we identified areas within inversions presenting clinal variation that might be under selection.
| S-EPMC3525579 | biostudies-literature
Project description:Chromosomal inversions play a fundamental role in evolution and have been shown to be responsible for the epidemiologically important traits in malaria mosquitoes. However, they have never been characterized in the major vector of arboviruses Aedes aegypti because of the poor structure of its polytene chromosomes. In this study, we applied a Hi-C proximity ligation approach to identify chromosomal inversions in 23 recently collected strains of Ae. aegypti from its worldwide distribution, two old laboratory colonies, and Ae. mascarensis.
2024-09-06 | GSE243024 | GEO
Project description:Adaptation is defined as an evolutionary process allowing organisms to succeed in certain habitats or conditions. Chromosomal inversions have the potential to be key in the adaptation processes, since they can contribute to the maintenance of favoured combinations of adaptive alleles through reduced recombination between individuals carrying different inversions. We have analysed six genes (Pif1A, Abi, Sqd, Yrt, Atpα and Fmr1), located inside and outside three inversions of the O chromosome in European populations of Drosophila subobscura. Genetic differentiation was significant between inversions despite extensive recombination inside inverted regions, irrespective of gene distance to the inversion breakpoints. Surprisingly, the highest level of genetic differentiation between arrangements was found for the Atpα gene, which is located outside the O1 and O7 inversions. Two derived unrelated arrangements (O3+4+1 and O3+4+7) are nearly fixed for several amino acid substitutions at the Atpα gene that have been described to confer resistance in other species to the cardenolide ouabain, a plant toxin capable of blocking ATPases. Similarities in the Atpα variants, conferring ouabain resistance in both arrangements, may be the result of convergent substitution and be favoured in response to selective pressures presumably related to the presence of plants containing ouabain in the geographic locations where both inversions are present.
| S-EPMC4815013 | biostudies-literature
Project description:We propose a genealogy-sampling algorithm, Sequential Markov Ancestral Recombination Tree (SMARTree), that provides an approach to estimation from SNP haplotype data of the patterns of coancestry across a genome segment among a set of homologous chromosomes. To enable analysis across longer segments of genome, the sequence of coalescent trees is modeled via the modified sequential Markov coalescent (Marjoram and Wall, Genetics 7:16, 2006). To assess performance in estimating these local trees, our SMARTree implementation is tested on simulated data. Our base data set is of the SNPs in 10 DNA sequences over 50 kb. We examine the effects of longer sequences and of more sequences, and of a recombination and/or mutational hotspot. The model underlying SMARTree is an approximation to the full recombinant-coalescent distribution. However, in a small trial on simulated data, recovery of local trees was similar to that of LAMARC (Kuhner et al. Genetics 156:1393-1401, 2000a), a sampler which uses the full model.
| S-EPMC4104301 | biostudies-literature