A sequential coalescent algorithm for chromosomal inversions

How chromosomal inversions reorient the evolutionary process

Inversions are structural mutations that reverse the sequence of a chromosome segment and reduce the effective rate of recombination in the heterozygous state. They play a major role in adaptation, as well as in other evolutionary processes such as speciation. Although inversions have been studied since the 1920s, they remain difficult to investigate because the reduced recombination conferred by them strengthens the effects of drift and hitchhiking, which in turn can obscure signatures of selection. Nonetheless, numerous inversions have been found to be under selection. Given recent advances in population genetic theory and empirical study, here we review how different mechanisms of selection affect the evolution of inversions. A key difference between inversions and other mutations, such as single nucleotide variants, is that the fitness of an inversion may be affected by a larger number of frequently interacting processes. This considerably complicates the analysis of the causes underlying the evolution of inversions. We discuss the extent to which these mechanisms can be disentangled, and by which approach.

Efficient ancestry and mutation simulation with msprime 1.0

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this, a large number of specialized simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and the tskit library. We summarize msprime’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialized alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

Chromosomal Speciation in the Genomics Era: Disentangling Phylogenetic Evolution of Rock-wallabies

The association of chromosome rearrangements (CRs) with speciation is well established, and there is a long history of theory and evidence relating to "chromosomal speciation." Genomic sequencing has the potential to provide new insights into how reorganization of genome structure promotes divergence, and in model systems has demonstrated reduced gene flow in rearranged segments. However, there are limits to what we can understand from a small number of model systems, which each only tell us about one episode of chromosomal speciation. Progressing from patterns of association between chromosome (and genic) change, to understanding processes of speciation requires both comparative studies across diverse systems and integration of genome-scale sequence comparisons with other lines of evidence. Here, we showcase a promising example of chromosomal speciation in a non-model organism, the endemic Australian marsupial genus Petrogale. We present initial phylogenetic results from exon-capture that resolve a history of divergence associated with extensive and repeated CRs. Yet it remains challenging to disentangle gene tree heterogeneity caused by recent divergence and gene flow in this and other such recent radiations. We outline a way forward for better integration of comparative genomic sequence data with evidence from molecular cytogenetics, and analyses of shifts in the recombination landscape and potential disruption of meiotic segregation and epigenetic programming. In all likelihood, CRs impact multiple cellular processes and these effects need to be considered together, along with effects of genic divergence. Understanding the effects of CRs together with genic divergence will require development of more integrative theory and inference methods. Together, new data and analysis tools will combine to shed light on long standing questions of how chromosome and genic divergence promote speciation.

The Evolution of Genome Structure by Natural and Sexual Selection

Progress on understanding how genome structure evolves is accelerating with the arrival of new genomic, comparative, and theoretical approaches. This article reviews progress in understanding how chromosome inversions and sex chromosomes evolve, and how their evolution affects species' ecology. Analyses of clines in inversion frequencies in flies and mosquitoes imply strong local adaptation, and roles for both over- and under dominant selection. Those results are consistent with the hypothesis that inversions become established when they capture locally adapted alleles. Inversions can carry alleles that are beneficial to closely related species, causing them to introgress following hybridization. Models show that this "adaptive cassette" scenario can trigger large range expansions, as recently happened in malaria mosquitoes. Sex chromosomes are the most rapidly evolving genome regions of some taxa. Sexually antagonistic selection may be the key force driving transitions of sex determination between different pairs of chromosomes and between XY and ZW systems. Fusions between sex-chromosomes and autosomes most often involve the Y chromosome, a pattern that can be explained if fusions are mildly deleterious and fix by drift. Sexually antagonistic selection is one of several hypotheses to explain the recent discovery that the sex determination system has strong effects on the adult sex ratios of tetrapods. The emerging view of how genome structure evolves invokes a much richer constellation of forces than was envisioned during the Golden Age of research on Drosophila karyotypes.

Efficient Coalescent Simulation and Genealogical Analysis for Large Sample Sizes

A central challenge in the analysis of genetic variation is to provide realistic genome simulation across millions of samples. Present day coalescent simulations do not scale well, or use approximations that fail to capture important long-range linkage properties. Analysing the results of simulations also presents a substantial challenge, as current methods to store genealogies consume a great deal of space, are slow to parse and do not take advantage of shared structure in correlated trees. We solve these problems by introducing sparse trees and coalescence records as the key units of genealogical analysis. Using these tools, exact simulation of the coalescent with recombination for chromosome-sized regions over hundreds of thousands of samples is possible, and substantially faster than present-day approximate methods. We can also analyse the results orders of magnitude more quickly than with existing methods.

Our understanding of the distribution of genetic variation in natural populations has been driven by mathematical models of the underlying biological and demographic processes. A key strength of such coalescent models is that they enable efficient simulation of data we might see under a variety of evolutionary scenarios. However, current methods are not well suited to simulating genome-scale data sets on hundreds of thousands of samples, which is essential if we are to understand the data generated by population-scale sequencing projects. Similarly, processing the results of large simulations also presents researchers with a major challenge, as it can take many days just to read the data files. In this paper we solve these problems by introducing a new way to represent information about the ancestral process. This new representation leads to huge gains in simulation speed and storage efficiency so that large simulations complete in minutes and the output files can be processed in seconds.

Genomic evidence for role of inversion 3RP of Drosophila melanogaster in facilitating climate change adaptation

Chromosomal inversion polymorphisms are common in animals and plants, and recent models suggest that alternative arrangements spread by capturing different combinations of alleles acting additively or epistatically to favour local adaptation. It is also thought that inversions typically maintain favoured combinations for a long time by suppressing recombination between alternative chromosomal arrangements. Here, we consider patterns of linkage disequilibrium and genetic divergence in an old inversion polymorphism in Drosophila melanogaster (In(3R)Payne) known to be associated with climate change adaptation and a recent invasion event into Australia. We extracted, karyotyped and sequenced whole chromosomes from two Australian populations, so that changes in the arrangement of the alleles between geographically separated tropical and temperate areas could be compared. Chromosome-wide linkage disequilibrium (LD) analysis revealed strong LD within the region spanned by In(3R)Payne. This genomic region also showed strong differentiation between the tropical and the temperate populations, but no differentiation between different karyotypes from the same population, after controlling for chromosomal arrangement. Patterns of differentiation across the chromosome arm and in gene ontologies were enhanced by the presence of the inversion. These data support the notion that inversions are strongly selected by bringing together combinations of genes, but it is still not clear if such combinations act additively or epistatically. Our data suggest that climatic adaptation through inversions can be dynamic, reflecting changes in the relative abundance of different forms of an inversion and ongoing evolution of allelic content within an inversion.

Advances in computer simulation of genome evolution: toward more realistic evolutionary genomics analysis by approximate bayesian computation

NGS technologies present a fast and cheap generation of genomic data. Nevertheless, ancestral genome inference is not so straightforward due to complex evolutionary processes acting on this material such as inversions, translocations, and other genome rearrangements that, in addition to their implicit complexity, can co-occur and confound ancestral inferences. Recently, models of genome evolution that accommodate such complex genomic events are emerging. This letter explores these novel evolutionary models and proposes their incorporation into robust statistical approaches based on computer simulations, such as approximate Bayesian computation, that may produce a more realistic evolutionary analysis of genomic data. Advantages and pitfalls in using these analytical methods are discussed. Potential applications of these ancestral genomic inferences are also pointed out.

Matrix inversions for chromosomal inversions: a method to construct summary statistics in complex coalescent models

Chromosomal inversions allow genetic divergence of locally adapted populations by reducing recombination between chromosomes with different arrangements. While patterns of genetic variation within inverted regions are increasingly documented, inferential methods are largely missing to analyze such data. Previous work has provided expectations for coalescence patterns of neutral sites linked to an inversion polymorphism in two locally adapted populations. Here, we define a method to construct summary statistics in such complex population structure models. Under a scenario of selection on the inversion breakpoints, we first construct estimators of the migration rate between the two habitats, and of the recombination rate of a nucleotide site between the two inversion backgrounds. Next, we analyze the disequilibrium between two sites within an inversion and provide an estimator of the distinct recombination rate between these two sites in homokaryotypes and heterokaryotypes. These estimators should be suitable summary statistics for simulation-based methods that can handle the complex dependences in the data.

Long range linkage disequilibrium across the human genome

Long-range linkage disequilibria (LRLD) between sites that are widely separated on chromosomes may suggest that population admixture, epistatic selection, or other evolutionary forces are at work. We quantified patterns of LRLD on a chromosome-wide level in the YRI population of the HapMap dataset of single nucleotide polymorphisms (SNPs). We calculated the disequilibrium between all pairs of SNPs on each chromosome (a total of >2×10(11) values) and evaluated significance of overall disequilibrium using randomization. The results show an excess of associations between pairs of distant sites (separated by >0.25 cM) on all of the 22 autosomes. We discuss possible explanations for this observation.