Genome-enabled hitchhiking mapping identifies QTLs for stress resistance in natural Drosophila

Natural and Artificial Selection for Parasitoid Resistance in Drosophila melanogaster Leave Different Genetic Signatures

Adaptation of complex traits depends on standing genetic variation at multiple loci. The allelic variants that have positive fitness effects, however, can differ depending on the genetic background and the selective pressure. Previously, we interrogated the Drosophila melanogaster genome at the population level for polymorphic positions and identified 215 single nucleotide polymorphisms (SNPs) that had significantly changed in frequency after experimental evolution for increased parasitoid resistance. In the current study, we follow up on 11 of these SNPs as putative targets of the experimental selection process (Jalvingh et al., 2014). We study the patterns of genetic variation for these SNPs in several European field populations. Furthermore, we associate the genetic variation of these SNPs to variation in resistance against the parasitoid Asobara tabida, by determining the individual phenotype and SNP genotype for 144 individuals from four Selection lines and four non-selected Control lines and for 400 individuals from 12 Field lines that differ in parasitoid resistance. For the Selection lines we additionally monitored the changes in allele frequencies throughout the five generations of experimental selection. For three genes, mbl (Zn-finger protein), mthl4 (G-protein coupled receptor) and CG17287 (protein-cysteine S-palmitoyltransferase) individual SNP genotypes were significantly associated with resistance level in the Selection and Control lines. Additionally, the minor allele in mbl and mthl4 were consistently and gradually favored throughout the five generations of experimental evolution. However, none of these alleles did appear to be associated to high resistance in the Field lines. We suggest that, within field populations, selection for parasitoid resistance is a gradual process that involves co-adapted gene complexes. Fast artificial selection, however, enforces the sudden cumulating of particular alleles that confer high resistance (genetic sweep). We discuss our findings in the context of local adaptation.

Dissection of Complex, Fitness-Related Traits in Multiple Drosophila Mapping Populations Offers Insight into the Genetic Control of Stress Resistance

We leverage two complementary Drosophila melanogaster mapping panels to genetically dissect starvation resistance-an important fitness trait. Using >1600 genotypes from the multiparental Drosophila Synthetic Population Resource (DSPR), we map numerous starvation stress QTL that collectively explain a substantial fraction of trait heritability. Mapped QTL effects allowed us to estimate DSPR founder phenotypes, predictions that were correlated with the actual phenotypes of these lines. We observe a modest phenotypic correlation between starvation resistance and triglyceride level, traits that have been linked in previous studies. However, overlap among QTL identified for each trait is low. Since we also show that DSPR strains with extreme starvation phenotypes differ in desiccation resistance and activity level, our data imply multiple physiological mechanisms contribute to starvation variability. We additionally exploited the Drosophila Genetic Reference Panel (DGRP) to identify sequence variants associated with starvation resistance. Consistent with prior work these sites rarely fall within QTL intervals mapped in the DSPR. We were offered a unique opportunity to directly compare association mapping results across laboratories since two other groups previously measured starvation resistance in the DGRP. We found strong phenotypic correlations among studies, but extremely low overlap in the sets of genomewide significant sites. Despite this, our analyses revealed that the most highly associated variants from each study typically showed the same additive effect sign in independent studies, in contrast to otherwise equivalent sets of random variants. This consistency provides evidence for reproducible trait-associated sites in a widely used mapping panel, and highlights the polygenic nature of starvation resistance.

Non-random transmission of parental alleles into crop-wild and crop-weed hybrid lineages separated by a transgene and neutral identifiers in rice

It is essential to assess environmental impact of transgene flow from genetically engineered crops to their wild or weedy relatives before commercialization. Measuring comparative trials of fitness in the transgene-flow-resulted hybrids plays the key role in the assessment, where the segregated isogenic hybrid lineages/subpopulations with or without a transgene of the same genomic background are involved. Here, we report substantial genomic differentiation between transgene-present and -absent lineages (F₂-F₃) divided by a glyphosate-resistance transgene from a crop-wild/weed hybrid population in rice. We further confirmed that such differentiation is attributed to increased frequencies of crop-parent alleles in transgenic hybrid lineages at multiple loci across the genome, as estimated by SSR (simple sequence repeat) markers. Such preferential transmission of parental alleles was also found in equally divided crop-wild/weed hybrid lineages with or without a particular neutral SSR identifier. We conclude that selecting either a transgene or neutral marker as an identifier to create hybrid lineages will result in different genomic background of the lineages due to non-random transmission of parental alleles. Non-random allele transmission may misrepresent the outcomes of fitness effects. We therefore propose seeking other means to evaluate fitness effects of transgenes for assessing environmental impact caused by crop-to-wild/weed gene flow.

Power analysis of artificial selection experiments using efficient whole genome simulation of quantitative traits

Evolve and resequence studies combine artificial selection experiments with massively parallel sequencing technology to study the genetic basis for complex traits. In these experiments, individuals are selected for extreme values of a trait, causing alleles at quantitative trait loci (QTL) to increase or decrease in frequency in the experimental population. We present a new analysis of the power of artificial selection experiments to detect and localize quantitative trait loci. This analysis uses a simulation framework that explicitly models whole genomes of individuals, quantitative traits, and selection based on individual trait values. We find that explicitly modeling QTL provides qualitatively different insights than considering independent loci with constant selection coefficients. Specifically, we observe how interference between QTL under selection affects the trajectories and lengthens the fixation times of selected alleles. We also show that a substantial portion of the genetic variance of the trait (50–100%) can be explained by detected QTL in as little as 20 generations of selection, depending on the trait architecture and experimental design. Furthermore, we show that power depends crucially on the opportunity for recombination during the experiment. Finally, we show that an increase in power is obtained by leveraging founder haplotype information to obtain allele frequency estimates.

Genomic changes under rapid evolution: selection for parasitoid resistance

In this study, we characterize changes in the genome during a swift evolutionary adaptation, by combining experimental selection with high-throughput sequencing. We imposed strong experimental selection on an ecologically relevant trait, parasitoid resistance in Drosophila melanogaster against Asobara tabida. Replicated selection lines rapidly evolved towards enhanced immunity. Larval survival after parasitization increased twofold after just five generations of selection. Whole-genome sequencing revealed that the fast and strong selection response in innate immunity produced multiple, highly localized genomic changes. We identified narrow genomic regions carrying a significant signature of selection, which were present across all chromosomes and covered in total less than 5% of the whole D. melanogaster genome. We identified segregating sites with highly significant changes in frequency between control and selection lines that fell within these narrow 'selected regions'. These segregating sites were associated with 42 genes that constitute possible targets of selection. A region on chromosome 2R was highly enriched in significant segregating sites and may be of major effect on parasitoid defence. The high genetic variability and small linkage blocks in our base population are likely responsible for allowing this complex trait to evolve without causing widespread erosive effects in the genome, even under such a fast and strong selective regime.

Promises and limitations of hitchhiking mapping

Building the connection between genetic and phenotypic variation is an important 'work in progress', and one that will enable proactive diagnosis and treatment in medicine, promote development of environment-targeted varieties in agriculture, and clarify the limits of species adaptation to changing environments in conservation. Quantitative trait loci (QTL) mapping and genome wide association (GWA) studies have recently been allied to an additional focus on 'hitchhiking' (HH) mapping--using changes in allele frequency due to artificial or natural selection. This older technique has been popularized by the falling costs of high throughput sequencing. Initial HH-resequensing experiments seem to have found many thousands of polymorphisms responding to selection. We argue that this interpretation appears too optimistic, and that the data might in fact be more consistent with dozens, rather than thousands, of loci under selection. We propose several developments required for sensible data analyses that will fully realize the great power of the HH technique, and outline ways of moving forward.

Maximum likelihood estimation of frequencies of known haplotypes from pooled sequence data

DNA samples are often pooled, either by experimental design or because the sample itself is a mixture. For example, when population allele frequencies are of primary interest, individual samples may be pooled together to lower the cost of sequencing. Alternatively, the sample itself may be a mixture of multiple species or strains (e.g., bacterial species comprising a microbiome or pathogen strains in a blood sample). We present an expectation–maximization algorithm for estimating haplotype frequencies in a pooled sample directly from mapped sequence reads, in the case where the possible haplotypes are known. This method is relevant to the analysis of pooled sequencing data from selection experiments, as well as the calculation of proportions of different species within a metagenomics sample. Our method outperforms existing methods based on single-site allele frequencies, as well as simple approaches using sequence read data. We have implemented the method in a freely available open-source software tool.

The opportunity for balancing selection in experimental populations of Caenorhabditis elegans

The role of balancing selection in maintaining diversity during the evolution of sexual populations to novel environments is poorly understood. To address this issue, we studied the impact of two mating systems, androdioecy and dioecy, on genotype distributions during the experimental evolution of Caenorhabditis elegans. We analyzed the temporal trajectories of 334 single nucleotide polymorphisms, covering 1/3 of the genome, and found extensive allele frequency changes and little loss of heterozygosities after 100 generations. As modeled with numerical simulations, SNP differentiation was consistent with genetic drift and average fitness effects of 2%, assuming that selection acted independently at each locus. Remarkably, inbreeding by self-fertilization was of little consequence to SNP differentiation. Modeling selection on deleterious recessive alleles suggests that the initial evolutionary dynamics can be explained by associative overdominance, but not the later stages because much lower heterozygosities would be maintained during experimental evolution. By contrast, models with selection on true overdominant loci can explain the heterozygote excess observed at all periods, particularly when negative epistasis or independent fitness effects were considered. Overall, these findings indicate that selection at single loci, including purging of recessive alleles, underlies most of the genetic differentiation accomplished during the experiment. Nonetheless, they also imply that maintenance of genetic diversity may in large part be due to balancing selection at multiple loci.

Genomic basis of aging and life-history evolution in Drosophila melanogaster

Natural diversity in aging and other life-history patterns is a hallmark of organismal variation. Related species, populations, and individuals within populations show genetically based variation in life span and other aspects of age-related performance. Population differences are especially informative because these differences can be large relative to within-population variation and because they occur in organisms with otherwise similar genomes. We used experimental evolution to produce populations divergent for life span and late-age fertility and then used deep genome sequencing to detect sequence variants with nucleotide-level resolution. Several genes and genome regions showed strong signatures of selection, and the same regions were implicated in independent comparisons, suggesting that the same alleles were selected in replicate lines. Genes related to oogenesis, immunity, and protein degradation were implicated as important modifiers of late-life performance. Expression profiling and functional annotation narrowed the list of strong candidate genes to 38, most of which are novel candidates for regulating aging. Life span and early age fecundity were negatively correlated among populations; therefore, the alleles we identified also are candidate regulators of a major life-history trade-off. More generally, we argue that hitchhiking mapping can be a powerful tool for uncovering the molecular bases of quantitative genetic variation.

Investigating natural variation in Drosophila courtship song by the evolve and resequence approach

A primary goal of population genetics is to determine the genetic basis of natural trait variation. We could significantly advance this goal by developing comprehensive genome-wide approaches to link genotype and phenotype in model organisms. Here we combine artificial selection with population-based resequencing to investigate the genetic basis of variation in the interpulse interval (IPI) of Drosophila melanogaster courtship song. We performed divergent selection on replicate populations for only 14 generations, but had considerable power to differentiate alleles that evolved due to selection from those that evolved stochastically. We identified a large number of variants that changed frequency in response to selection for this simple behavior, and they are highly underrepresented on the X chromosome. Though our power was adequate using this experimental technique, the ability to differentiate causal variants from those affected by linked selection requires further development.

High resolution mapping of candidate alleles for desiccation resistance in Drosophila melanogaster under selection

The ability to counter periods of low humidity is an important determinant of distribution range in Drosophila. Climate specialists with low physiological tolerance to desiccation stress are restricted to the tropics and may lack the ability to further increase resistance through evolution. Although the physiological adaptations to desiccation stress are well studied in Drosophila and other ectotherms, factors underlying evolutionary responses remain unknown because of a paucity of genetic data. We address this issue by mapping evolutionary shifts in D. melanogaster under selection for desiccation resistance. Genomic DNA from five independent replicate selected, and control lines were hybridized to high density Affymetrix Drosophila tiling arrays resulting in the detection of 691 single feature polymorphisms (SFPs) differing between the treatments. While randomly distributed throughout the genome, the SFPs formed specific clusters according to gene ontology. These included genes involved in ion transport and respiratory system development that provide candidates for evolutionary changes involving excretory and respiratory water balance. Changes to genes related to neuronal control of cell signaling, development, and gene regulation provide candidates to explore novel biological processes in stress resistance. Sequencing revealed the nucleotide shifts in a subset of the SFPs and highlighted larger regions of genomic diversity surrounding SFPs. The association between natural desiccation resistance and a 463-bp region of the 5' promoter region of the Dys gene undergoing allele frequency changes in response to selection in the experimental evolution lines was tested in an independent population from Coffs Harbour, Australia. The allele frequencies of 23 SNPs common to the two populations were inferred from the parents of the 10% most and 10% least resistant Coffs Harbour flies. The frequencies of the selected alleles were higher at all sites, with three sites significantly associated with the resistant Coffs Harbour flies. This study illustrates how rapid mapping can be used for discovering natural molecular variants associated with survival to low humidity and provides a wealth of candidate alleles to explore the genetic basis of physiological differences among resistant and susceptible Drosophila populations and species.

Mapping Mendelian traits in asexual progeny using changes in marker allele frequency

Linkage between markers and genes that affect a phenotype of interest may be determined by examining differences in marker allele frequency in the extreme progeny of a cross between two inbred lines. This strategy is usually employed when pooling is used to reduce genotyping costs. When the cross progeny are asexual, the extreme progeny may be selected by multiple generations of asexual reproduction and selection. We analyse this method of measuring phenotype in asexual progeny and examine the changes in marker allele frequency due to selection over many generations. Stochasticity in marker frequency in the selected population arises due to the finite initial population size. We derive the distribution of marker frequency as a result of selection at a single major locus, and show that in order to avoid spurious changes in marker allele frequency in the selected population, the initial population size should be in the low to mid hundreds.

Population-based resequencing of experimentally evolved populations reveals the genetic basis of body size variation in Drosophila melanogaster

Body size is a classic quantitative trait with evolutionarily significant variation within many species. Locating the alleles responsible for this variation would help understand the maintenance of variation in body size in particular, as well as quantitative traits in general. However, successful genome-wide association of genotype and phenotype may require very large sample sizes if alleles have low population frequencies or modest effects. As a complementary approach, we propose that population-based resequencing of experimentally evolved populations allows for considerable power to map functional variation. Here, we use this technique to investigate the genetic basis of natural variation in body size in Drosophila melanogaster. Significant differentiation of hundreds of loci in replicate selection populations supports the hypothesis that the genetic basis of body size variation is very polygenic in D. melanogaster. Significantly differentiated variants are limited to single genes at some loci, allowing precise hypotheses to be formed regarding causal polymorphisms, while other significant regions are large and contain many genes. By using significantly associated polymorphisms as a priori candidates in follow-up studies, these data are expected to provide considerable power to determine the genetic basis of natural variation in body size.

Understanding the causes and consequences of natural genetic variation is crucial to the characterization of biological evolution. Moreover, natural genetic variation is comprised of millions of perturbations, which are partially randomized across genotypes such that a small number of individuals can be used to combinatorially analyze a large number of differences, facilitating mechanistic understanding of biological systems. Here we demonstrate a powerful technique to parse genomic variation using artificial selection. By selecting replicate populations of Drosophila flies to become bigger and smaller, and then determining the evolutionary response at the genomic level, we have mapped hundreds of genes that respond to selection on body size. As our approach is powerful and cost-effective compared to existing approaches, we expect it to be a major component of diverse future efforts.

Information content in genome-wide scans: concordance between patterns of genetic differentiation and linkage mapping associations

BACKGROUND

Scanning the genome with high density SNP markers has become a standard approach for identifying regions of the genome showing substantial between-population genetic differentiation, and thus evidence of diversifying selection. Such regions may contain genes of large phenotypic effect. However, few studies have attempted to address the power or efficacy of such an approach.

RESULTS

In this study, the patterns of allele frequency differences between two cattle breeds based on the Bovine HapMap study were compared with statistical evidence for QTL based on a linkage mapping study of an experimental population formed by a cross between the same breeds. Concordance between the two datasets was seen for chromosomes carrying QTL with strong statistical support, such as BTA5 and BTA18, which carry genes associated with coat color. For these chromosomes, there was a correspondence between the strength of the QTL signal along the chromosome and the degree of genetic differentiation between breeds. However, such an association was not seen in a broader comparison that also included chromosomes carrying QTL with lower significance levels. In addition, other chromosomal regions with substantial QTL effects did not include markers showing extreme between-breed genetic differentiation. Furthermore, the overall consistency between the two studies was weak, with low genome-wide correlation between the statistical values obtained in the linkage mapping study and between-breed genetic differentiation from the HapMap study.

CONCLUSIONS

These results suggest that genomic diversity scans are capable of detecting regions associated with qualitative traits but may be limited in their power to detect regions associated with quantitative phenotypic differences between populations, which may depend on the marker resolution of the study and the level of LD in the populations under investigation.

Combined expression patterns of QTL-linked candidate genes best predict thermotolerance in Drosophila melanogaster

Knockdown resistance to high temperature (KRHT) is a thermal adaptation trait in Drosophila melanogaster. Here we used quantitative real-time PCR (qRT-PCR) to test for possible associations between KRHT and the expression of candidate genes within quantitative trait loci (QTL) in eight recombinant inbred lines (RIL). hsp60 and hsc70-3 map within an X-linked QTL, while CG10383, catsup, ddc, trap1, and cyp6a13 are linked in a KRHT-QTL on chromosome 2. hsc70-3 expression increased by heat-hardening. Principal Components analysis revealed that catsup, ddc and trap1 were either co-expressed or combined in their expression levels. This composite expression variable (e-PC1) was positively associated to KRHT in non-hardened RIL. In heat-hardened flies, hsp60 was negatively related to hsc70-3 on e-PC2, with effects on KRHT. These results are consistent with the notion that QTL can be shaped by expression variation in combined candidate loci. We found composite variables of gene expression (e-PCs) that best correlated to KRHT. Network effects with other untested linked loci are apparent because, in spite of their associations with KRHT phenotypes, e-PCs were sometimes uncorrelated with their QTL genotype.

Regulatory divergence in Drosophila melanogaster and D. simulans, a genomewide analysis of allele-specific expression

Species-specific regulation of gene expression contributes to the development and maintenance of reproductive isolation and to species differences in ecologically important traits. A better understanding of the evolutionary forces that shape regulatory variation and divergence can be developed by comparing expression differences among species and interspecific hybrids. Once expression differences are identified, the underlying genetics of regulatory variation or divergence can be explored. With the goal of associating cis and/or trans components of regulatory divergence with differences in gene expression, overall and allele-specific expression levels were assayed genomewide in female adult heads of Drosophila melanogaster, D. simulans, and their F1 hybrids. A greater proportion of cis differences than trans differences were identified for genes expressed in heads and, in accordance with previous studies, cis differences also explained a larger number of species differences in overall expression level. Regulatory divergence was found to be prevalent among genes associated with defense, olfaction, and among genes downstream of the Drosophila sex determination hierarchy. In addition, two genes, with critical roles in sex determination and micro RNA processing, Sxl and loqs, were identified as misexpressed in hybrid female heads, potentially contributing to hybrid incompatibility.

Experimental evolution withDrosophila

Experimental evolution is a powerful approach that can be used for the study of adaptation. Evolutionary biologists often use Drosophila as a model organism in experiments that test theories about the evolution of traits related to fitness. Such evolution experiments can take three forms: direct selection for a trait of interest; surveys of traits of interest in populations selected for other traits; and reverse selection. We review some of the Drosophila experiments that have provided insight into both the evolution of particular physiological traits and the correlations between physiological and life history traits, focusing on stress resistance. The most common artifacts that can obscure the results from evolution experiments are discussed. We also include a treatment of genomic technologies that are now available for the Drosophila model. The primary goal of this review is to introduce the kind of experimental evolution strategies and technologies that evolutionary physiologists might use in the future.

Genomic analysis of differentiation between soil types reveals candidate genes for local adaptation in Arabidopsis lyrata

Serpentine soil, which is naturally high in heavy metal content and has low calcium to magnesium ratios, comprises a difficult environment for most plants. An impressive number of species are endemic to serpentine, and a wide range of non-endemic plant taxa have been shown to be locally adapted to these soils. Locating genomic polymorphisms which are differentiated between serpentine and non-serpentine populations would provide candidate loci for serpentine adaptation. We have used the Arabidopsis thaliana tiling array, which has 2.85 million probes throughout the genome, to measure genetic differentiation between populations of Arabidopsis lyrata growing on granitic soils and those growing on serpentinic soils. The significant overrepresentation of genes involved in ion transport and other functions provides a starting point for investigating the molecular basis of adaptation to soil ion content, water retention, and other ecologically and economically important variables. One gene in particular, calcium-exchanger 7, appears to be an excellent candidate gene for adaptation to low Ca∶Mg ratio in A. lyrata.

Selection mapping of loci for quantitative disease resistance in a diverse maize population

The selection response of a complex maize population improved primarily for quantitative disease resistance to northern leaf blight (NLB) and secondarily for common rust resistance and agronomic phenotypes was investigated at the molecular genetic level. A tiered marker analysis with 151 simple sequence repeat (SSR) markers in 90 individuals of the population indicated that on average six alleles per locus were available for selection. An improved test statistic for selection mapping was developed, in which quantitative trait loci (QTL) are identified through the analysis of allele-frequency shifts at mapped multiallelic loci over generations of selection. After correcting for the multiple tests performed, 25 SSR loci showed evidence of selection. Many of the putatively selected loci were unlinked and dispersed across the genome, which was consistent with the diffuse distribution of previously published QTL for NLB resistance. Compelling evidence for selection was found on maize chromosome 8, where several putatively selected loci colocalized with published NLB QTL and a race-specific resistance gene. Analysis of F2 populations derived from the selection mapping population suggested that multiple linked loci in this chromosomal segment were, in part, responsible for the selection response for quantitative resistance to NLB.

Detecting genetic responses to environmental change

Changes in environmental conditions can rapidly shift allele frequencies in populations of species with relatively short generation times. Frequency shifts might be detectable in neutral genetic markers when stressful conditions cause a population decline. However, frequency shifts that are diagnostic of specific conditions depend on isolating sets of genes that are involved in adaptive responses. Shifts at candidate loci underlying adaptive responses and DNA regions that control their expression have now been linked to evolutionary responses to pollution, global warming and other changes. Conversely, adaptive constraints, particularly in physiological traits, are recognized through DNA decay in candidate genes. These approaches help researchers and conservation managers understand the power and constraints of evolution.