Alan O. Bergland

Empirical Validation of Pooled Whole Genome Population Re-Sequencing in Drosophila melanogaster

The sequencing of pooled non-barcoded individuals is an inexpensive and efficient means of assessing genome-wide population allele frequencies, yet its accuracy has not been thoroughly tested. We assessed the accuracy of this approach on whole, complex eukaryotic genomes by resequencing pools of largely isogenic, individually sequenced Drosophila melanogaster strains. We called SNPs in the pooled data and estimated false positive and false negative rates using the SNPs called in individual strain as a reference. We also estimated allele frequency of the SNPs using “pooled” data and compared them with “true” frequencies taken from the estimates in the individual strains. We demonstrate that pooled sequencing provides a faithful estimate of population allele frequency with the error well approximated by binomial sampling, and is a reliable means of novel SNP discovery with low false positive rates. However, a sufficient number of strains should be used in the pooling because variation in the amount of DNA derived from individual strains is a substantial source of noise when the number of pooled strains is low. Our results and analysis confirm that pooled sequencing is a very powerful and cost-effective technique for assessing of patterns of sequence variation in populations on genome-wide scales, and is applicable to any dataset where sequencing individuals or individual cells is impossible, difficult, time consuming, or expensive.

Quantitative trait loci affecting phenotypic plasticity and the allometric relationship of ovariole number and thorax length in Drosophila melanogaster

Environmental factors during juvenile growth such as temperature and nutrition have major effects on adult morphology and life-history traits. In Drosophila melanogaster, ovary size, measured as ovariole number, and body size, measured as thorax length, are developmentally plastic traits with respect to larval nutrition. Herein we investigated the genetic basis for plasticity of ovariole number and body size, as well the genetic basis for their allometric relationship using recombinant inbred lines (RILs) derived from a natural population in Winters, California. We reared 196 RILs in four yeast concentrations and measured ovariole number and body size. The genetic correlation between ovariole number and thorax length was positive, but the strength of this correlation decreased with increasing yeast concentration. Genetic variation and genotype-by-environment (G × E) interactions were observed for both traits. We identified quantitative trait loci (QTL), epistatic, QTL-by-environment, and epistatic-by-environment interactions for both traits and their scaling relationships. The results are discussed in the context of multivariate trait evolution.

Secondary contact and local adaptation contribute to genome‐wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome‐wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome‐wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an FST outlier approach. Our results provide novel insight into the well‐studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.

Drosophila suzukii: The Genetic Footprint of a Recent, Worldwide Invasion

Native to Asia, the soft-skinned fruit pest Drosophila suzukii has recently invaded the United States and Europe. The eastern United States represents the most recent expansion of their range, and presents an opportunity to test alternative models of colonization history. Here, we investigate the genetic population structure of this invasive fruit fly, with a focus on the eastern United States. We sequenced six X-linked gene fragments from 246 individuals collected from a total of 12 populations. We examine patterns of genetic diversity within and between populations and explore alternative colonization scenarios using approximate Bayesian computation. Our results indicate high levels of nucleotide diversity in this species and suggest that the recent invasions of Europe and the continental United States are independent demographic events. More broadly speaking, our results highlight the importance of integrating population structure into demographic models, particularly when attempting to reconstruct invasion histories. Finally, our simulation results illustrate the general challenge in reconstructing invasion histories using genetic data and suggest that genome-level data are often required to distinguish among alternative demographic scenarios.

A highly pleiotropic amino acid polymorphism in the Drosophila insulin receptor contributes to life-history adaptation

Finding the specific nucleotides that underlie adaptive variation is a major goal in evolutionary biology, but polygenic traits pose a challenge because the complex genotype–phenotype relationship can obscure the effects of individual alleles. However, natural selection working in large wild populations can shift allele frequencies and indicate functional regions of the genome. Previously, we showed that the two most common alleles of a complex amino acid insertion–deletion polymorphism in the Drosophila insulin receptor show independent, parallel clines in frequency across the North American and Australian continents. Here, we report that the cline is stable over at least a five‐year period and that the polymorphism also demonstrates temporal shifts in allele frequency concurrent with seasonal change. We tested the alleles for effects on levels of insulin signaling, fecundity, development time, body size, stress tolerance, and life span. We find that the alleles are associated with predictable differences in these traits, consistent with patterns of Drosophila life‐history variation across geography that likely reflect adaptation to the heterogeneous climatic environment. These results implicate insulin signaling as a major mediator of life‐history adaptation in Drosophila, and suggest that life‐history trade‐offs can be explained by extensive pleiotropy at a single locus.

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster.

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.

LDx: Estimation of Linkage Disequilibrium from High-Throughput Pooled Resequencing Data

High-throughput pooled resequencing offers significant potential for whole genome population sequencing. However, its main drawback is the loss of haplotype information. In order to regain some of this information, we present LDx, a computational tool for estimating linkage disequilibrium (LD) from pooled resequencing data. LDx uses an approximate maximum likelihood approach to estimate LD (r2) between pairs of SNPs that can be observed within and among single reads. LDx also reports r2 estimates derived solely from observed genotype counts. We demonstrate that the LDx estimates are highly correlated with r2 estimated from individually resequenced strains. We discuss the performance of LDx using more stringent quality conditions and infer via simulation the degree to which performance can improve based on read depth. Finally we demonstrate two possible uses of LDx with real and simulated pooled resequencing data. First, we use LDx to infer genomewide patterns of decay of LD with physical distance in D. melanogaster population resequencing data. Second, we demonstrate that r2 estimates from LDx are capable of distinguishing alternative demographic models representing plausible demographic histories of D. melanogaster.

Fine-scale mapping of natural variation in fly fecundity identifies neuronal domain of expression and function of an aquaporin.

To gain insight into the molecular genetic basis of standing variation in fitness related traits, we identify a novel factor that regulates the molecular and physiological basis of natural variation in female Drosophila melanogaster fecundity. Genetic variation in female fecundity in flies derived from a wild orchard population is heritable and largely independent of other measured life history traits. We map a portion of this variation to a single QTL and then use deficiency mapping to further refine this QTL to 5 candidate genes. Ubiquitous expression of RNAi against only one of these genes, an aquaporin encoded by Drip, reduces fecundity. Within our mapping population Drip mRNA level in the head, but not other tissues, is positively correlated with fecundity. We localize Drip expression to a small population of corazonin producing neurons located in the dorsolateral posterior compartments of the protocerebrum. Expression of Drip–RNAi using both the pan-neuronal ELAV-Gal4 and the Crz-Gal4 drivers reduces fecundity. Low-fecundity RILs have decreased Crz expression and increased expression of pale, the enzyme encoding the rate-limiting step in the production of dopamine, a modulator of insect life histories. Taken together these data suggest that natural variation in Drip expression in the corazonin producing neurons contributes to standing variation in fitness by altering the concentration of two neurohormones.

Global Transcriptional Profiling of Diapause and Climatic Adaptation in Drosophila melanogaster

Abstract Wild populations of the model organism Drosophila melanogaster experience highly heterogeneous environments over broad geographical ranges as well as over seasonal and annual timescales. Diapause is a primary adaptation to environmental heterogeneity, and in D. melanogaster the propensity to enter diapause varies predictably with latitude and season. Here we performed global transcriptomic profiling of naturally occurring variation in diapause expression elicited by short day photoperiod and moderately low temperature in two tissue types associated with neuroendocrine and endocrine signaling, heads, and ovaries. We show that diapause in D. melanogaster is an actively regulated phenotype at the transcriptional level, suggesting that diapause is not a simple physiological or reproductive quiescence. Differentially expressed genes and pathways are highly distinct in heads and ovaries, demonstrating that the diapause response is not uniform throughout the soma and suggesting that it may be comprised of functional modules associated with specific tissues. Genes downregulated in heads of diapausing flies are significantly enriched for clinally varying single nucleotide polymorphism (SNPs) and seasonally oscillating SNPs, consistent with the hypothesis that diapause is a driving phenotype of climatic adaptation. We also show that chromosome location-based coregulation of gene expression is present in the transcriptional regulation of diapause. Taken together, these results demonstrate that diapause is a complex phenotype actively regulated in multiple tissues, and support the hypothesis that natural variation in diapause propensity underlies adaptation to spatially and temporally varying selective pressures.

Adaptive dynamics of cuticular hydrocarbons in Drosophila

Cuticular hydrocarbons (CHCs) are hydrophobic compounds deposited on the arthropod cuticle that are of functional significance with respect to stress tolerance, social interactions and mating dynamics. We characterized CHC profiles in natural populations of Drosophila melanogaster at five levels: across a latitudinal transect in the eastern United States, as a function of developmental temperature during culture, across seasonal time in replicate years, and as a function of rapid evolution in experimental mesocosms in the field. Furthermore, we also characterized spatial and temporal changes in allele frequencies for SNPs in genes that are associated with the production and chemical profile of CHCs. Our data demonstrate a striking degree of parallelism for clinal and seasonal variation in CHCs in this taxon; CHC profiles also demonstrate significant plasticity in response to rearing temperature, and the observed patterns of plasticity parallel the spatiotemporal patterns observed in nature. We find that these congruent shifts in CHC profiles across time and space are also mirrored by predictable shifts in allele frequencies at SNPs associated with CHC chain length. Finally, we observed rapid and predictable evolution of CHC profiles in experimental mesocosms in the field. Together, these data strongly suggest that CHC profiles respond rapidly and adaptively to environmental parameters that covary with latitude and season, and that this response reflects the process of local adaptation in natural populations of D. melanogaster.