A sequential coalescent algorithm for chromosomal inversions

Chromosomal inversions are common in natural populations and are believed to be involved in many important evolutionary phenomena, including speciation, the evolution of sex chromosomes and local adaptation. While recent advances in sequencing and genotyping methods are leading to rapidly increasing amounts of genome-wide sequence data that reveal interesting patterns of genetic variation within inverted regions, efficient simulation methods to study these patterns are largely missing. In this work, we extend the sequential Markovian coalescent, an approximation to the coalescent with recombination, to include the effects of polymorphic inversions on patterns of recombination. Results show that our algorithm is fast, memory-efficient and accurate, making it feasible to simulate large inversions in large populations for the first time. The SMC algorithm enables studies of patterns of genetic variation (for example, linkage disequilibria) and tests of hypotheses (using simulation-based approaches) that were previously intractable.     

Contrasting patterns of phenotypic variation linked to chromosomal inversions in native and colonizing populations of Drosophila subobscura

In fewer than two decades after invading the Americas, the fly Drosophila subobscura evolved latitudinal clines for chromosomal inversion frequencies and wing size that are parallel to the long‐standing ones in native Palearctic populations. By sharp contrast, wing shape clines also evolved in the New World, but the relationship with latitude was opposite to that in the Old World. Previous work has suggested that wing trait differences among individuals are partially due to the association between chromosomal inversions and particular alleles which influence the trait under consideration. Furthermore, it is well documented that a few number of effective individuals founded the New World populations, which might have modified the biometrical effect of inversions on quantitative traits. Here we evaluate the relative contribution of chromosomal inversion clines in shaping the parallel clines in wing size and contrasting clines in wing shape in native and colonizing populations of the species. Our results reveal that inversion‐size and inversion‐shape associations in native and colonizing (South America) populations are generally different, probably due to the bottleneck effect. Contingent, unpredictable evolution was suggested as an explanation for the different details involved in the otherwise parallel wing size clines between Old and New World populations of D. subobscura. We challenge this assertion and conclude that contrasting wing shape clines came out as a correlated response of inversion clines that might have been predicted considering the genetic background of colonizers.     

A Coalescent Model for a Sweep of a Unique Standing Variant

The use of genetic polymorphism data to understand the dynamics of adaptation and identify the loci that are involved has become a major pursuit of modern evolutionary genetics. In addition to the classical “hard sweep” hitchhiking model, recent research has drawn attention to the fact that the dynamics of adaptation can play out in a variety of different ways and that the specific signatures left behind in population genetic data may depend somewhat strongly on these dynamics. One particular model for which a large number of empirical examples are already known is that in which a single derived mutation arises and drifts to some low frequency before an environmental change causes the allele to become beneficial and sweeps to fixation. Here, we pursue an analytical investigation of this model, bolstered and extended via simulation study. We use coalescent theory to develop an analytical approximation for the effect of a sweep from standing variation on the genealogy at the locus of the selected allele and sites tightly linked to it. We show that the distribution of haplotypes that the selected allele is present on at the time of the environmental change can be approximated by considering recombinant haplotypes as alleles in the infinite-alleles model. We show that this approximation can be leveraged to make accurate predictions regarding patterns of genetic polymorphism following such a sweep. We then use simulations to highlight which sources of haplotypic information are likely to be most useful in distinguishing this model from neutrality, as well as from other sweep models, such as the classic hard sweep and multiple-mutation soft sweeps. We find that in general, adaptation from a unique standing variant will likely be difficult to detect on the basis of genetic polymorphism data from a single population time point alone, and when it can be detected, it will be difficult to distinguish from other varieties of selective sweeps. Samples from multiple populations and/or time points have the potential to ease this difficulty.     

The foldback-like transposon Galileo is involved in the generation of two different natural chromosomal inversions of Drosophila buzzatii

Chromosomal inversions are the most common type of genome rearrangement in the genus Drosophila. Although the potential of transposable elements (TEs) for generating inversions has been repeatedly demonstrated in the laboratory, little is known on their role in the generation of natural inversions, which are those effectively contributing to the adaptation and/or evolution of species. We have cloned and sequenced the two breakpoints of the polymorphic inversion 2q7 of D. buzzatii. The sequence analysis of the breakpoint regions revealed the presence in the inverted chromosomes of large insertions, formed by complex assemblies of transposons, that are absent from the chromosomes without the inversion. Among the transposons inserted, the Foldback-like element Galileo, that was previously found responsible of the generation of the widespread inversion 2j of D. buzzatii, is present at both 2q7 breakpoints and is the most likely inducer of the inversion. A detailed study of the nucleotide and structural variation in the breakpoint regions of six chromosomal lines with the 2q7 inversion detected no nucleotide differences between them, which suggests a monophyletic and recent origin. In contrast, a remarkable degree of structural variation was observed in the same six chromosomal lines. It thus appears that the two breakpoints of the inverted chromosomes have become genetically unstable hotspots, as was previously found for the 2j inversion breakpoints. The possibility that this instability is caused by structural properties of Foldback elements is discussed.     

The DAIBAM MITE element is involved in the origin of one fixed and two polymorphic Drosophila virilis phylad inversions

Chromosomal inversions can originate from breakage and repair by non-homologous end-joining. Nevertheless, they can also originate from ectopic recombination between transposable elements located on the same chromosome inserted in opposite orientations. Here, we show that a MITE element (DAIBAM), previously involved in the origin of one Drosophila americana polymorphic inversion, is also involved in the origin of one fixed inversion between D. virilis and D. americana and another D. americana polymorphic inversion. Therefore, DAIBAM is responsible for at least 20% of the chromosomal rearrangements that are observed within and between species of the virilis phylad (D. virilis, D. lummei, D. novamexicana and D. americana), having thus played a significant role in the chromosomal evolution of this group of closely related species.     

Genome-wide evolutionary response to a heat wave in Drosophila

Extreme climatic events can substantially affect organismal performance and Darwinian fitness. In April 2011, a strong heat wave struck extensive geographical areas of the world, including Western Europe. At that time, we happened to resume and extend a long-term time series of seasonal genetic data in the widespread fly Drosophila subobscura, which provided a unique opportunity to quantify the intensity of the genetic perturbation caused by the heat wave. We show that the spring 2011 genetic constitution of the populations transiently shifted to summer-like frequencies, and that the magnitude of the genetic anomaly quantitatively matched the temperature anomaly. The results provide compelling evidence that direct effects of rising temperature are driving adaptive evolutionary shifts, and also suggest a strong genetic resilience in this species.     

Coalescent approximation for structured populations in a stationary random environment

We establish convergence to the Kingman coalescent for the genealogy of a geographically-or otherwise-structured version of the Wright-Fisher population model with fast migration. The new feature is that migration probabilities may change in a random fashion. This brings a novel formula for the coalescent effective population size (EPS). We call it a quenched EPS to emphasize the key feature of our model — random environment. The quenched EPS is compared with an annealed (mean-field) EPS which describes the case of constant migration probabilities obtained by averaging the random migration probabilities over possible environments.     

Tracking the origin of the American colonization by Drosophila subobscura: genetic comparison between Eastern and Western Mediterranean populations

Several pieces of evidence indicate a Mediterranean origin of the colonization of America by Drosophila subobscura . To ascertain whether the origin was from the Eastern or the Western Mediterranean region, samples from Barcelona (Spain) and Mt Parnes (Greece) were collected and O chromosomal inversion polymorphism and lethal genes were analysed. The frequencies of lethal chromosomes were 0.244 ± 0.039 in Barcelona and 0.336 ± 0.043 in Mt Parnes, consistent with the expectations for large populations located in the central area of the species distribution. Lethal genes seem to be distributed at random along the O chromosome in both populations. The intra-populational allelism frequencies of Barcelona and Mt Parnes were 0.016 ± 0.007 and 0.012 ± 0.005 respectively. Thus, the estimates of the effective population size were high in both populations (between 6964 and 13 004 in Barcelona and 11 874 to 26 828 in Mt Parnes). The cases of allelism in Mt Parnes were observed only between individual lethal genes, but in Barcelona some concatenated clusters of allelism were detected. This pattern of allelism can be explained by synthetic lethality, hybrid dysgenesis, the induction of recurrent lethal mutations by different factors or an effect of microdifferentiation in subpopulations. In both populations, a reduction in fitness in the heterozygotes for lethal genes has been detected. Furthermore, the estimates of the migration coefficient (between 0.0085 and 0.0120 in Barcelona, and 0.0057 and 0.0087 in Mt Parnes) confirm the existence of gene flow between Palearctic populations of D. subobscura . Our lethal genes and chromosomal inversion results are consistent with a Mediterranean origin of the colonization of America by D. subobscura , but are inconclusive with regard to the identification of the population from which these colonizers came.     

Genetic loads and coadaptivity of chromosomal inversions I. Inversion polymorphism and genetic load on chromosome O in a Drosophila subobscura population from Petnica

Inversion polymorphism on chromosome O and polymorphism for the viability of determining genes have been studied in a natural population of Drosophila subobscura from Petnica (Serbia). The range of inversion polymorphism and the abundance of particular gene arrangements in the study population agree with a general pattern of inversion polymorphism of D. subobscura in Europe. The data obtained on the amount of genetic loads show that the D. subobscura population from Petnica displays a moderate degree of that polymorphism, compared to the other studied populations of these species. Therefore, the D. subobscura population from Petnica could be tentatively classified as an ecologically central population. Examination association of chromosomal, thus, inversion polymorphism with gene polymorphism, in the form of genetic loads show that differences exist in the mean viability among certain gene arrangements. The distribution of deleterious genes among chromosome O gene arrangements were non-random.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

A microsatellite linkage map for Drosophila montana shows large variation in recombination rates,and a courtship song trait maps to an area of low recombination

Current advances in genetic analysis are opening up our knowledge of the genetics of species differences, but challenges remain, particularly for out‐bred natural populations. We constructed a microsatellite‐based linkage map for two out‐bred lines of Drosophila montana derived from divergent populations by taking advantage of the Drosophila virilis genome and available cytological maps of both species. Although the placement of markers was quite consistent with cytological predictions, the map indicated large heterogeneity in recombination rates along chromosomes. We also performed a quantitative trait locus (QTL) analysis on a courtship song character (carrier frequency), which differs between populations and is subject to strong sexual selection. Linkage mapping yielded two significant QTLs, which explained 3% and 14% of the variation in carrier frequency, respectively. Interestingly, as in other recent studies of traits which can influence speciation, the strongest QTL mapped to a genomic region partly covered by an inversion polymorphism.     

Adaptation genomics: next generation sequencing reveals a shared haplotype for rapid early development in geographically and genetically distant populations of rainbow trout

Local adaptation occurs when a population evolves a phenotype that confers a selective advantage in its local environment, but which may not be advantageous in other habitats. Restricted gene flow and strong selection pressures are prerequisites for local adaptation. Fishes in the family Salmonidae are predicted to provide numerous examples of local adaptation because of the high fidelity of returning to spawn in their natal streams, which results in highly structured populations, and the wide diversity of environments that salmonids have colonized. These conditions are ideally suited for producing a set of specialist phenotypes, whose fitness is maximized for one specific habitat, rather than a generalist phenotype similarly viable in several environments. Understanding patterns and processes leading to local adaptations has long been a goal of evolutionary biology, but it is only recently that identifying the molecular basis for local adaptation has become feasible because of advances in genomic technologies. The study of shared adaptive phenotypes in populations that are both geographically distant and genetically distinct should reveal some of the fundamental molecular mechanisms associated with local adaptation. In this issue of Molecular Ecology, Miller et al. (2012) make a significant contribution to the development of adaptation genomics. This study suggests that salmonids use standing genetic variation to select beneficial alleles for local adaptations rather than de novo mutations in the same gene or alternative physiological pathways. Identifying the genetic basis for local adaptation has major implications for the management, conservation and potential restoration of salmonid populations.     

Genetic diversity of photosynthetic characters in native populations of Triticum dicoccoides

Wild emmer wheat (Triticum dicoccoides Korn) has shown wide genetic diversity for disease resistance and morpho-physiological traits of economic importance. Our objectives were to test for genetic variation (V_G) in photosynthetic characteristics residing within and between native populations sampled from three ecogeographical regions of Israel, and to identify potential sources of high photosynthetic efficiency for future wheat improvement. Accessions sampled in the center of wild emmer distribution (upper Jordan Valley) in a relatively narrow geographical range showed the greatest diversity in CO₂-assimilation rate per unit leaf area (A) or per unit chlorophyll (A/Chl). Genetic variation was absent for internal CO₂ concentration (Ci) and water-use efficiency (WUE) and generally lacking for stomatal conductance (gs). Leaf area, although quite variable, was not a significant cofactor in assessing genetic potential for photosynthesis. Accessions within a given population showed 10-times more variation in A and A/Chl than populations sampled from different locations in a region. Accessions with the highest photosynthetic efficiency were derived from upland steppic populations located in marginal habitats extending southward into Israel. Some accessions having high photosynthetic capacity (A=32 mol m^-2 s^-1) with no significant reduction in leaf size constitute a potentially valuable genetic resource yet untapped for genetic improvement of hexaploid (T. aestivum L.) wheat.Abbreviations LA total leaf area - V_G total genetic variance - PPFD photosynthetic photon flux density - WUE water-use efficiency     

Altitudinal and seasonal variation in microsatellite allele frequencies of Drosophila buzzatii

Variation in climate, particularly temperature, is known to affect the genetic composition of populations. Although there have been many studies of latitudinal variation, comparisons of populations across altitudes or seasons, particularly for animal species, are less common. Here, we study genetic variation (microsatellite markers) in populations of Drosophila buzzatii collected along altitudinal gradients and in different seasons. We found no differences in genetic variation between 2 years or between seasons within years. However, there were numerous cases of significant associations between allele frequencies or expected heterozygosities and altitude, with more than half showing nonlinear relationships. While these associations indicate possible selection and local altitudinal adaptation, direct tests gave strong evidence for selection affecting two loci and weaker evidence for five other loci. Two loci that are located within an inversion (including the one with strongest evidence for selection) show a linear increase in genetic diversity with altitude, likely due to thermal selection. Parallel associations with altitude here and with latitude in Australian populations indicate that selection is operating on chromosomal regions marked by some of the loci.     

Paracentric inversions in human chromosome 7 as a graphic example of reverse chromosomal mutation

A famillialy inherited paracentric inversion of human chromosome 7 is described. The breakpoints are localized in the bands q1 1.2 and q2 2.1. The similarity between the inverted chromosome 7 and the chromosome 6 of the gorilla characterizes this inversion as a reverse chromosomal mutation. The parallels between human chromosome pathology and the chromosomal evolution of Pongidae and man are discussed.     

Evidence for local adaptation in neighbouring Daphnia populations: a laboratory transplant experiment

1. We tested whether two neighbouring Daphnia galeata populations (from Lake Blankaart and Fish Pond), separated only by 5 m of land and with the occasional exchange of water were genetically differentiated in allozyme markers and life history traits. Allozyme electrophoresis revealed that the populations differed in allelic as well as in genotypic composition. 2. In a laboratory transplant experiment, in which animals of four clones of each of the populations were raised in the water of both ponds, survival in Blankaart water was high for both the Blankaart and Fish Pond clones, whereas survival in Fish Pond water was high for the Fish Pond clones, but low for the Blankaart clones. 3. Fish Pond clones produced fewer neonates than Blankaart clones when cultured in Blankaart water. High egg mortality was observed for animals that were raised in Blankaart water, and this egg mortality was higher for Fish Pond clones than for Blankaart clones. 4. Our results provide evidence for genetic differentiation between Daphnia populations inhabiting neighbouring water bodies and suggest local adaptation to environmental conditions other than direct predation.