Cosmopolitan Inversions Research Articles

A cosmopolitan inversion facilitates seasonal adaptation in overwintering Drosophila.

Fluctuations in the strength and direction of natural selection through time are a ubiquitous feature of life on Earth. One evolutionary outcome of such fluctuations is adaptive tracking, wherein populations rapidly adapt from standing genetic variation. In certain circumstances, adaptive tracking can lead to the long-term maintenance of functional polymorphism despite allele frequency change due to selection. Although adaptive tracking is likely a common process, we still have a limited understanding of aspects of its genetic architecture and its strength relative to other evolutionary forces such as drift. Drosophila melanogaster living in temperate regions evolve to track seasonal fluctuations and are an excellent system to tackle these gaps in knowledge. By sequencing orchard populations collected across multiple years, we characterized the genomic signal of seasonal demography and identified that the cosmopolitan inversion In(2L)t facilitates seasonal adaptive tracking and shows molecular footprints of selection. A meta-analysis of phenotypic studies shows that seasonal loci within In(2L)t are associated with behavior, life history, physiology, and morphological traits. We identify candidate loci and experimentally link them to phenotype. Our work contributes to our general understanding of fluctuating selection and highlights the evolutionary outcome and dynamics of contemporary selection on inversions.

Pool-GWAS on reproductive dormancy in Drosophila simulans suggests a polygenic architecture.

The genetic basis of adaptation to different environments has been of long-standing interest to evolutionary biologists. Dormancy is a well-studied adaptation to facilitate overwintering. In Drosophila melanogaster, a moderate number of genes with large effects have been described, which suggests a simple genetic basis of dormancy. On the other hand, genome-wide scans for dormancy suggest a polygenic architecture in insects. In D. melanogaster, the analysis of the genetic architecture of dormancy is complicated by the presence of cosmopolitan inversions. Here, we performed a genome-wide scan to characterize the genetic basis of this ecologically extremely important trait in the sibling species of D. melanogaster, D. simulans that lacks cosmopolitan inversions. We performed Pool-GWAS in a South African D. simulans population for dormancy incidence at 2 temperature regimes (10 and 12°C, LD 10:14). We identified several genes with SNPs that showed a significant association with dormancy (P-value < 1e-13), but the overall modest response suggests that dormancy is a polygenic trait with many loci of small effect. Our results shed light on controversies on reproductive dormancy in Drosophila and have important implications for the characterization of the genetic basis of this trait.

Genetic structuring of Drosophila bipectinata in Indian natural populations based on the distribution of cosmopolitan inversions

Genetic differentiation among different natural populations of a species depends upon the environmental factors and the evolutionary forces that operate on them. In this study, seven Indian natural populations of D. bipectinata, two from north and five from south India, have been studied for their chromosomal inversion polymorphism. A total of nine paracentric autosomal inversions were recorded from these seven places but only three of them, present on the 2L, 2R and 3L were found to be cosmopolitan in distribution. In all the populations, the frequency of standard gene arrangement was found to be high than their respective cosmopolitan inversion gene arrangement. The average heterozygosity (Ho) of cosmopolitan inversions increases from north to south. There is a latitudinal cline in the distribution of three cosmopolitan inversion arrangements because their frequency increases with the decreasing latitude, i.e. from north to south India. A comparison of the genetic profile of two north Indian and five south Indian natural populations of D. bipectinata reveals the role of natural selection as well as bottleneck effect in the genetic structuring of these populations which may be due to their varying ecological conditions to which they are constantly encountered. Further, the presence of all kinds of paracentric inversions in individual populations was analysed following Poisson distribution to see whether these inversions occur randomly in natural populations or not and the results indicate that north Indian populations show the random occurrence of these inversions than the populations derived from the south.

Dobzhansky’s concept of genetic coadaptation: Drosophila ananassae is an exception to this concept

Dobzhansky was the first to show that the inversion polymorphism in Drosophila pseudoobscura is subject to natural selection and is a device to cope with the diversity of environments. His extensive work on D. pseudoobscura has revealed interesting phenomena of population genetics. In continuation of his work on this species, he constructed a number of homozygous lines for different gene arrangements in the third chromosome, and while employing these lines in intrapopulation and interpopulation crosses, he quantified the fitness of inversion homokaryotypes and heterokaryotypes. Interestingly, his results showed that heterokaryotypes formed by chromosomes originating from the same geographic area exhibited superiority over the corresponding homokaryotypes. However, superiority of heterokaryotypes was lost in the crosses when chromosomes were derived from different localities. Based on these results, Dobzhansky suggested the concept of genetic coadaptation. According to this concept, 'in each locality, the chromosomes with different gene arrangements aremutually adjusted or coadapted to yield highly fit inversion heterozygotes through long continued natural selection. However, this adaptive superiority of inversion heterozygotes breaks down in interracial hybridization experiments when two gene arrangements are derived from different localities'. This concept has received experimental evidence in its favour on the basis of work done in other species of Drosophila, such as D. willistoni, D. paulistorum, D. pavani and D. bipectinata. In all these species, interracial hybridization led to the loss of superiority of inversion heterozygotes. Further, it has been suggested that coadapted polygenic complexes contained in the chromosomes are disrupted as a result of recombination in interstrain crosses. This concept was also tested in D. ananassae, a cosmopolitan and domestic species commonly found in India, while employing three cosmopolitan inversions exhibiting heterotic buffering. In interstrain crosses involving monomorphic and polymorphic strains due to three cosmopolitan inversions, the persistence of heterosis was observed, which does not support the above-mentioned hypothesis of Dobzhansky. Thus, evidence for coadaptation is lacking in natural populations of D. ananassae, which is considered as an exception to the Dobzhansky's concept of genetic coadaptation. Thus, heterotic buffering associated with the three cosmopolitan inversions in D. ananassae is not populational heterosis; rather, it appears to be simple luxuriance.

Pleiotropic effects of regulatory variation in tan result in correlation of two pigmentation traits in Drosophila melanogaster.

Traits with a common genetic basis frequently display correlated phenotypic responses to selection or environmental conditions. In Drosophila melanogaster, pigmentation of the abdomen and a trident‐shaped region on the thorax are genetically correlated. Here, we used a pooled replicated genomewide association approach (Pool‐GWAS) to identify the genetic basis of variation in thoracic trident pigmentation in two Drosophila melanogaster populations. We confirmed the previously reported large effect of ebony and the association of the cosmopolitan inversion In(3R)Payne. For the first time, we identified tan as another major locus contributing to variation in trident pigmentation. Intriguingly, the regulatory variants of tan that were most strongly associated with female abdominal pigmentation also showed a strong association with trident pigmentation. We validated this common genetic basis in transgenic assays and found qualitatively similar effects on trident and abdominal pigmentation. Further work is required to determine whether this genetic correlation is favoured by natural selection or reflects a neutral by‐product of a shared regulatory architecture.

Frequencies of chromosomal inversions in Drosophila melanogaster in Fukushima after the nuclear power plant accident.

In order to investigate genetic impact of a large amount of radionuclides released by the Fukushima Dai-ichi Nuclear Power Plant accident in 2011, we surveyed 2,304 haploid genomes of Drosophila melanogaster collected in three localities in Fukushima in 2012 and 2013 for chromosomal inversions. No unique inversion was found in 298 genomes in 2012 and only two in 2,006 genomes in 2013. The observed frequencies were even lower than the long-term average frequency of unique inversions in Japan. The common cosmopolitan inversions were also examined in Fukushima, Kyoto, and Iriomote (Okinawa) in 2012. Among three samples in Fukushima, the flies in Iizaka, where environmental radiation level was the highest, showed the lowest frequency of In(2L)t, but the highest frequency of In(3R)P, contrary to the expectation of decreasing of their frequencies in higher polluted areas. These results suggest that, at this level of genetic analysis, Fukushima populations of D. melanogaster would not have been negatively impacted following the release of radionuclides. Transposable P-element mobility was not likely to induce DNA damage solely or synergistically with radioactivity, because their transposition activity was totally repressed in the Fukushima strains. However, it should be noted that, because of limitations in access to the exclusion zone, we could only sample the populations in areas of relatively low radioactive contamination (0.39–0.63 μSv/h). Therefore, the present study is likely to be underpowered to detect any effects that might be expected in heavily contaminated areas.

The Effect of Common Inversion Polymorphisms In(2L)t and In(3R)Mo on Patterns of Transcriptional Variation in Drosophila melanogaster.

Chromosomal inversions are a ubiquitous feature of genetic variation. Theoretical models describe several mechanisms by which inversions can drive adaptation and be maintained as polymorphisms. While inversions have been shown previously to be under selection, or contain genetic variation under selection, the specific phenotypic consequences of inversions leading to their maintenance remain unclear. Here we use genomic sequence and expression data from the Drosophila Genetic Reference Panel (DGRP) to explore the effects of two cosmopolitan inversions, In(2L)t and In(3R)Mo, on patterns of transcriptional variation. We demonstrate that each inversion has a significant effect on transcript abundance for hundreds of genes across the genome. Inversion-affected loci (IAL) appear both within inversions as well as on unlinked chromosomes. Importantly, IAL do not appear to be influenced by the previously reported genome-wide expression correlation structure. We found that five genes involved with sterol uptake, four of which are Niemann-Pick Type 2 orthologs, are upregulated in flies with In(3R)Mo but do not have SNPs in linkage disequilibrium (LD) with the inversion. We speculate that this upregulation is driven by genetic variation in mod(mdg4) that is in LD with In(3R)Mo. We find that there is little evidence for a regional or position effect of inversions on gene expression at the chromosomal level, but do find evidence for the distal breakpoint of In(3R)Mo interrupting one gene and possibly disassociating the two flanking genes from regulatory elements.

Parallel effects of the inversion In(3R)Payne on body size across the North American and Australian clines in Drosophila melanogaster.

Chromosomal inversions are thought to play a major role in climatic adaptation. In D. melanogaster, the cosmopolitan inversion In(3R)Payne exhibits latitudinal clines on multiple continents. As many fitness traits show similar clines, it is tempting to hypothesize that In(3R)P underlies observed clinal patterns for some of these traits. In support of this idea, previous work in Australian populations has demonstrated that In(3R)P affects body size but not development time or cold resistance. However, similar data from other clines of this inversion are largely lacking; finding parallel effects of In(3R)P across multiple clines would considerably strengthen the case for clinal selection. Here, we have analysed the phenotypic effects of In(3R)P in populations originating from the endpoints of the latitudinal cline along the North American east coast. We measured development time, egg-to-adult survival, several size-related traits (femur and tibia length, wing area and shape), chill coma recovery, oxidative stress resistance and triglyceride content in homokaryon lines carrying In(3R)P or the standard arrangement. Our central finding is that the effects of In(3R)P along the North American cline match those observed in Australia: standard arrangement lines were larger than inverted lines, but the inversion did not influence development time or cold resistance. Similarly, In(3R)P did not affect egg-to-adult survival, oxidative stress resistance and lipid content. In(3R)P thus seems to specifically affect size traits in populations from both continents. This parallelism strongly suggests an adaptive pattern, whereby the inversion has captured alleles associated with growth regulation and clinal selection acts on size across both continents.

Population genetic evidence for cold adaptation in European Drosophila melanogaster populations

We studied Drosophila melanogaster populations from Europe (the Netherlands and France) and Africa (Rwanda and Zambia) to uncover genetic evidence of adaptation to cold. We present here four lines of evidence for genes involved in cold adaptation from four perspectives: (i) the frequency of SNPs at genes previously known to be associated with chill-coma recovery time (CCRT), startle reflex (SR) and resistance to starvation stress (RSS) vary along environmental gradients and therefore among populations; (ii) SNPs of genes that correlate significantly with latitude and altitude in African and European populations overlap with SNPs that correlate with a latitudinal cline from North America; (iii) at the genomewide level, the top candidate genes are enriched in gene ontology (GO) terms that are related to cold tolerance; (iv) GO enriched terms from North American clinal genes overlap significantly with those from Africa and Europe. Each SNP was tested in 10 independent runs of Bayenv2, using the median Bayes factors to ascertain candidate genes. None of the candidate genes were found close to the breakpoints of cosmopolitan inversions, and only four candidate genes were linked to QTLs related to CCRT. To overcome the limitation that we used only four populations to test correlations with environmental gradients, we performed simulations to estimate the power of our approach for detecting selection. Based on our results, we propose a novel network of genes that is involved in cold adaptation.

Genetic differentiation in natural populations and their mass culture stocks of Drosophila ananassae

Natural populations and their mass culture stocks of Drosophila ananassae were analyzed for the frequencies of three cosmopolitan inversions. The frequency of ST gene rearrangement was recorded to be high in both natural populations as well as in the mass culture stocks. Furthermore, out of three cosmopolitan inversions, the AL and DE had higher occurrence showing their adaptive significance to this species. ET inversion that persists in all populations, was thinly represented than the other two. These natural populations, coming from different geographical localities, showed very little genetic variation. The mass culture stocks also showed almost similar genetic composition as their respective natural populations. Various parameters of population genetics were computed in order to find out genetic structure of the natural populations. The results indicate that genetic changes in the natural populations have scarcely undergone, and no significant genetic differentiation could be observed when natural populations were compared with the respective mass culture stocks.

Inference of chromosomal inversion dynamics fromPool‐Seq data in natural and laboratory populations ofDrosophila melanogaster

Sequencing of pools of individuals (Pool-Seq) represents a reliable and cost-effective approach for estimating genome-wide SNP and transposable element insertion frequencies. However, Pool-Seq does not provide direct information on haplotypes so that, for example, obtaining inversion frequencies has not been possible until now. Here, we have developed a new set of diagnostic marker SNPs for seven cosmopolitan inversions in Drosophila melanogaster that can be used to infer inversion frequencies from Pool-Seq data. We applied our novel marker set to Pool-Seq data from an experimental evolution study and from North American and Australian latitudinal clines. In the experimental evolution data, we find evidence that positive selection has driven the frequencies of In(3R)C and In(3R)Mo to increase over time. In the clinal data, we confirm the existence of frequency clines for In(2L)t, In(3L)P and In(3R)Payne in both North America and Australia and detect a previously unknown latitudinal cline for In(3R)Mo in North America. The inversion markers developed here provide a versatile and robust tool for characterizing inversion frequencies and their dynamics in Pool-Seq data from diverse D. melanogaster populations.

Genome‐wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome-wide next-generation sequencing of DNA pools ('pool-seq') from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome-wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well-known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.

Complex patterns of local adaptation in heat tolerance in Drosophila simulans from eastern Australia

Latitudinal clines are considered a powerful means of investigating evolutionary responses to climatic selection in nature. However, most clinal studies of climatic adaptation in Drosophila have involved species that contain cosmopolitan inversion polymorphisms that show clinal patterns themselves, making it difficult to determine whether the traits or inversions are under selection. Further, although climatic selection is unlikely to act on only one life stage in metamorphic organisms, a few studies have examined clinal patterns across life stages. Finally, clinal patterns of heat tolerance may also depend on the assay used. To unravel these potentially confounding effects on clinal patterns of thermal tolerance, we examined adult and larval heat tolerance traits in populations of Drosophila simulans from eastern Australia using static and dynamic (ramping 0.06 °C min(-1)) assays. We also used microsatellites markers to clarify whether demographic factors or selection are responsible for population differentiation along clines. Significant cubic clinal patterns were observed for adult static basal, hardened and dynamic heat knockdown time and static basal heat survival in larvae. In contrast, static, hardened larval heat survival increased linearly with latitude whereas no clinal association was found for larval ramping survival. Significant associations between adult and larval traits and climatic variables, and low population differentiation at microsatellite loci, suggest a role for climatic selection, rather than demographic processes, in generating these clinal patterns. Our results suggest that adaptation to thermal stress may be species and life-stage specific, complicating our efforts to understand the evolutionary responses to selection for increasing thermotolerance.

Divergent enhancer haplotype of ebony on inversion In(3R)Payne associated with pigmentation variation in a tropical population of Drosophila melanogaster

The pattern and intensity of pigmentation have direct impact on individual fitness through various ecological factors. In a Drosophila melanogaster population from southern Japan, thoracic trident pigmentation intensity of most of the strains could be classified into Dark or Light-type. The expression level variation of the ebony gene correlated well with this phenotype and the allelic differences in expression indicated that the variation is partly due to cis-regulatory changes. In the ∼13 kb gene region, we identified 17 nucleotide sites and 2 indels that were in complete association with the thoracic trident pigmentation intensity. Interestingly, 11 out of 19 sites located within ∼0.5 kb of the core epidermis enhancer. These sites had no obvious association with the abdominal pigmentation intensity in the previously analysed African populations from Uganda and Kenya, which suggested that multiple potential mutational pathways in the cis-regulatory control region of a single gene could lead to similar phenotypic variation within this species. We also found that the Light-type enhancer haplotype is strongly linked to a cosmopolitan inversion, In(3R)Payne, which is predominant in warmer climatic regions in both hemispheres. The sequence pattern suggested that the strong linkage may be due to selective forces related to thermal adaptation. The inferred selection for lighter pigmentation in the Japanese population is in the opposite direction of the previously reported case of selection for darker individuals in African populations. Nevertheless, both adaptive changes involved cis-regulatory changes of ebony, which shows that this gene is likely to be a common target of natural selection.

Genomic differentiation between temperate and tropical Australian populations of Drosophila melanogaster.

Determining the genetic basis of environmental adaptation is a central problem of evolutionary biology. This issue has been fruitfully addressed by examining genetic differentiation between populations that are recently separated and/or experience high rates of gene flow. A good example of this approach is the decades-long investigation of selection acting along latitudinal clines in Drosophila melanogaster. Here we use next-generation genome sequencing to reexamine the well-studied Australian D. melanogaster cline. We find evidence for extensive differentiation between temperate and tropical populations, with regulatory regions and unannotated regions showing particularly high levels of differentiation. Although the physical genomic scale of geographic differentiation is small--on the order of gene sized--we observed several larger highly differentiated regions. The region spanned by the cosmopolitan inversion polymorphism In(3R)P shows higher levels of differentiation, consistent with the major difference in allele frequencies of Standard and In(3R)P karyotypes in temperate vs. tropical Australian populations. Our analysis reveals evidence for spatially varying selection on a number of key biological processes, suggesting fundamental biological differences between flies from these two geographic regions.

Population Genetics of Drosophila ananassae: Evidence for Population Sub-Structuring at the Level of Inversion Polymorphism in Indian Natural Populations

Drosophila ananassae is a cosmopolitan and domestic species distributed in the tropical, subtropical and mildlytemperate regions. Population structure analysis in forty-five Indian natural populations of D. ananassae was performedemploying three cosmopolitan inversions as markers. Pairwise FST analysis and genetic distance (D) values showed stronggenetic differentiation. Though, lowermost values correspond to geographically closest populations, we did not find anysignificant ‘isolation by distance’ effect. Values of gene flow based on FST estimates are very low (Nm < 5). All thesefindings, viz. strong genetic differentiation and minimal gene flow indicate strong sub-structuring in Indian naturalpopulations of D. ananassae at the level of inversion polymorphism. This finding is particularly intriguing in case of D.ananassae as it is frequently transported via human traffic. Given limited gene flow, populations are expected to divergegenetically due to drift. Low level of gene flow coupled with high degree of genetic differentiation might have occurredhistorically and is maintained currently. Demographic properties, historical and contemporary events and other factors aremore important in shaping the patterns of population sub-structuring, genetic differentiation and gene flow than mereterrestrial habitat characteristics (un) favorable for migration.

Investigating latitudinal clines for life history and stress resistance traits in Drosophila simulans from eastern Australia

Latitudinal clines have been demonstrated for many quantitative traits in Drosophila and are assumed to be due to climatic selection. However, clinal studies are often performed in species of Drosophila that contain common cosmopolitan inversion polymorphisms that also show clinal patterns. These inversion polymorphisms may be responsible for much of the observed clinal variation. Here, we consider latitudinal clines for quantitative traits in Drosophila simulans from eastern Australia. Drosophila simulans does not contain cosmopolitan inversion polymorphisms, so allows the study of clinal selection on quantitative traits that are not confounded by associations with inversions. Body size showed a strong linear cline for both females and males. Starvation resistance exhibited a weak linear cline in females, whereas chill-coma recovery exhibited a significant nonlinear cline in females only. No clinal pattern was evident for development time, male chill-coma recovery, desiccation or heat resistance. We discuss these results with reference to the role inversion polymorphisms play in generating clines in quantitative traits of Drosophila.

Population genetics ofDrosophila ananassae

Drosophila ananassae Doleschall is a cosmopolitan and domestic species. It occupies a unique status among Drosophila species due to certain peculiarities in its genetic behaviour and is of common occurrence in India. Quantitative genetics of sexual and non-sexual traits provided evidence for genetic control of these traits. D. ananassae exhibits high level of chromosomal polymorphism in its natural populations. Indian natural populations of D. ananassae show geographic differentiation of inversion polymorphism due to their adaptation to varying environments and natural selection operates to maintain three cosmopolitan inversions. Populations do not show divergence on temporal scale, an evidence for rigid polymorphism. D. ananassae populations show substantial degree of sub-structuring and exist as semi-isolated populations. Gene flow is low despite co-transportation with human goods. There is persistence of cosmopolitan inversions when populations are transferred to laboratory conditions, which suggests that heterotic buffering is associated with these inversions in D. ananassae. Populations collected from similar environmental conditions that initially show high degree of genetic similarity have diverged to different degrees in laboratory environment. This randomness could be due to genetic drift. Interracial hybridization does not lead to breakdown of heterosis associated with cosmopolitan inversions, which shows that there is lack of genetic co-adaptation in D. ananassae. Linkage disequilibrium between independent inversions in laboratory populations has often been observed, which is likely to be due to suppression of crossing-over and random genetic drift. No evidence for chromosomal interactions has been found in natural and laboratory populations of D. ananassae. This strengthens the previous suggestion that there is lack of genetic co-adaptation in D. ananassae.

Inversion polymorphism in Indian populations of Drosophila ananassae

Twelve natural populations of Drosophila ananassae from India were analysed for chromosome inversions. The chromosomal analysis revealed the presence of all three previously described cosmopolitan inversions in the populations studied. The frequency of inversions and the level of inversion heterozygosity vary significantly among populations. The x2 values calculated to measure the difference in karyotype frequencies between populations demonstrate highly significant variation in most of the cases. To quantify the amount of genetic differentiation at the level of chromosomal polymorphism in these populations, the average genetic identity (I) has been calculated on the basis of frequencies of different gene sequences. The values of I suggest that Indian populations of D. ananassae have undergone considerable amount of genetic divergence at the level of chromosomal polymorphism.

Mapping Regions Within Cosmopolitan Inversion In(3R)Payne Associated With Natural Variation in Body Size in Drosophila melanogaster

Associations between genotypes for inversions and quantitative traits have been reported in several organisms, but little has been done to localize regions within inversions controlling variation in these traits. Here, we use an association mapping technique to identify genomic regions controlling variation in wing size within the cosmopolitan inversion In(3R)Payne in Drosophila melanogaster. Previous studies have shown that this inversion strongly influences variation in wing size, a trait highly correlated with body size. We found three alleles from two separate regions within In(3R)Payne with significant additive effects on wing size after the additional effect of the inversion itself had been taken into account. There were also several alleles with significant genotype-by-inversion interaction effects on wing size. None of the alleles tested had a significant additive effect on development time, suggesting different genes control these traits and that clinal patterns in them have therefore arisen independently. The presence of multiple regions within In(3R)Payne controlling size is consistent with the idea that inversions persist in populations because they contain multiple sets of locally adapted alleles, but more work needs to be done to test if they are indeed coadapted.