Balancing selection and the functional effects of shared polymorphism in cryptic Daphnia species.

The patterns of genetic variation within and between related taxa represent the genetic history of a species. Shared polymorphisms, loci with identical alleles across species, are of unique interest as they may represent cases of ancient selection maintaining functional variation post-speciation. In this study, we investigate the abundance of shared polymorphism in the Daphnia pulex species complex. We test whether shared mutations are consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting, or convergent evolution. We analyzed over 2,000 genomes from North American and European D. pulex and several outgroup species to examine the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that while North American and European D. pulex diverged over ten million years ago, they retained tens of thousands of shared alleles. We found that the number of shared polymorphisms between North American and European D. pulex cannot be explained by hybridization or incomplete lineage sorting alone. Instead, we show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting young and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild. Ultimately, our findings provide insights into the genetic basis of adaptation and the maintenance of genetic diversity between species.

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.

A small survey of introduced Zaprionus indianus (Diptera: Drosophilidae) in orchards of the eastern United States.

The African fig fly, Zaprionus indianus (Gupta), is a generalist fruit fly that typically breeds in decaying fruits from over 70 plant species. The species has spread globally from its native range in tropical Africa, becoming an invasive pest on ripening figs in Brazil. First reported in the United States in 2005 in Florida, Z. indianus has since been documented as far north as Canada and is hypothesized to recolonize northwards from southern refugia each year. We sampled drosophilid communities over the growing season at 2 orchards in Virginia from 2020 to 2022 and 11 orchards along the East Coast during the fall of 2022 to quantify the abundance of Z. indianus relative to other drosophilids across locations, seasons, and fruit crops. Massachusetts had the northernmost population, with no Z. indianus detected in Maine and no correlation between latitude and relative abundance. Variation in Z. indianus relative abundance was high between nearby orchards and abundance was higher on peaches relative to apples within orchards. Comparisons of seasonal abundance curves between 2 Virginia orchards showed similar dynamics across years with individuals first detected around July and becoming absent around December, with peaks in late summer and mid-fall. The variation in seasonal and latitudinal abundance shown here highlights a need for broader sampling to accurately characterize the range, spread, and environmental tolerances of Z. indianus in North America.

Spatial and temporal variation in abundance of introduced African fig fly (Zaprionus indianus) (Diptera: Drosophilidae) in the eastern United States.

The African fig fly, Zaprionus indianus (Gupta), has spread globally from its native range in tropical Africa, becoming an invasive crop pest in select areas such as Brazil. Z. indianus was first reported in the United States in 2005 and has since been documented as far north as Canada. As a tropical species, Z. indianus is expected to have low cold tolerance, likely limiting its ability to persist at northern latitudes. In North America, the geographic regions where Z. indianus can thrive and seasonal fluctuations in its abundance are not well understood. The purpose of this study was to characterize the temporal and spatial variation in Z. indianus abundance to better understand its invasion of the eastern United States. We sampled drosophilid communities over the growing season at two orchards in Virginia from 2020-2022 and several locations along the East Coast during the fall of 2022. Virginia abundance curves showed similar seasonal dynamics across years with individuals first detected around July and becoming absent around December. Massachusetts was the northernmost population and no Z. indianus were detected in Maine. Variation in Z. indianus relative abundance was high between nearby orchards and across different fruits within orchards but was not correlated with latitude. Fitness of wild-caught females decreased later in the season and at higher latitudes. The patterns of Z. indianus abundance shown here demonstrate an apparent susceptibility to cold and highlight a need for systematic sampling to accurately characterize the range and spread of Z. indianus.

Distinct signals of clinal and seasonal allele frequency change at eQTLs in Drosophila melanogaster.

Populations of short-lived organisms can respond to spatial and temporal environmental heterogeneity through local adaptation. Local adaptation can be reflected on both phenotypic and genetic levels, and it has been documented in many organisms. Although complex fitness-related phenotypes have been shown to vary across latitudinal clines and seasons in similar ways in Drosophila melanogaster populations, the comparative signals of local adaptation across space and time remain poorly understood. Here, we examined patterns of allele frequency change across a latitudinal cline and between seasons at previously reported expression quantitative trait loci (eQTLs). We divided eQTLs into groups by using differential expression profiles of fly populations collected across latitudinal clines or exposed to different environmental conditions. In general, we find that eQTLs are enriched for clinally varying polymorphisms, and that these eQTLs change in frequency in concordant ways across the cline and in response to starvation and chill-coma. The enrichment of eQTLs among seasonally varying polymorphisms is more subtle, and the direction of allele frequency change at eQTLs appears to be somewhat idiosyncratic. Taken together, we suggest that clinal adaptation at eQTLs is at least partially distinct from seasonal adaptation.

Adaptive phenotypic plasticity is under stabilizing selection in Daphnia.

The adaptive nature of phenotypic plasticity is widely documented. However, little is known about the evolutionary forces that shape genetic variation of plasticity within populations. Whether genetic variation in plasticity is driven by stabilizing or diversifying selection and whether the strength of such forces remains constant through time, remain open questions. Here, we address this issue by assessing the evolutionary forces that shape genetic variation in antipredator developmental plasticity of Daphnia pulex. Antipredator plasticity in D. pulex is characterized by the growth of a pedestal and spikes in the dorsal head region upon exposure to predator cue. We characterized genetic variation in plasticity using a method that describes the entire dorsal shape amongst >100 D. pulex strains recently derived from the wild. We observed the strongest reduction in genetic variation in dorsal areas where plastic responses were greatest, consistent with stabilizing selection. We compared mutational variation (Vm) to standing variation (Vg) and found that Vg/Vm is lowest in areas of greatest plasticity, again consistent with stabilizing selection. Our results suggest that stabilizing selection operates directly on phenotypic plasticity in Daphnia and provide a rare glimpse into the evolution of fitness-related traits in natural populations.

Genetic Variation in Reproductive Investment Across an Ephemerality Gradient in Daphnia pulex.

Species across the tree of life can switch between asexual and sexual reproduction. In facultatively sexual species, the ability to switch between reproductive modes is often environmentally dependent and subject to local adaptation. However, the ecological and evolutionary factors that influence the maintenance and turnover of polymorphism associated with facultative sex remain unclear. We studied the ecological and evolutionary dynamics of reproductive investment in the facultatively sexual model species, Daphnia pulex. We found that patterns of clonal diversity, but not genetic diversity varied among ponds consistent with the predicted relationship between ephemerality and clonal structure. Reconstruction of a multi-year pedigree demonstrated the coexistence of clones that differ in their investment into male production. Mapping of quantitative variation in male production using lab-generated and field-collected individuals identified multiple putative quantitative trait loci (QTL) underlying this trait, and we identified a plausible candidate gene. The evolutionary history of these QTL suggests that they are relatively young, and male limitation in this system is a rapidly evolving trait. Our work highlights the dynamic nature of the genetic structure and composition of facultative sex across space and time and suggests that quantitative genetic variation in reproductive strategy can undergo rapid evolutionary turnover.

Broad geographic sampling reveals the shared basis and environmental correlates of seasonal adaptation in Drosophila.

To advance our understanding of adaptation to temporally varying selection pressures, we identified signatures of seasonal adaptation occurring in parallel among Drosophila melanogaster populations. Specifically, we estimated allele frequencies genome-wide from flies sampled early and late in the growing season from 20 widely dispersed populations. We identified parallel seasonal allele frequency shifts across North America and Europe, demonstrating that seasonal adaptation is a general phenomenon of temperate fly populations. Seasonally fluctuating polymorphisms are enriched in large chromosomal inversions, and we find a broad concordance between seasonal and spatial allele frequency change. The direction of allele frequency change at seasonally variable polymorphisms can be predicted by weather conditions in the weeks prior to sampling, linking the environment and the genomic response to selection. Our results suggest that fluctuating selection is an important evolutionary force affecting patterns of genetic variation in Drosophila.

Accurate, ultra-low coverage genome reconstruction and association studies in Hybrid Swarm mapping populations.

Genetic association studies seek to uncover the link between genotype and phenotype, and often utilize inbred reference panels as a replicable source of genetic variation. However, inbred reference panels can differ substantially from wild populations in their genotypic distribution, patterns of linkage-disequilibrium, and nucleotide diversity. As a result, associations discovered using inbred reference panels may not reflect the genetic basis of phenotypic variation in natural populations. To address this problem, we evaluated a mapping population design where dozens to hundreds of inbred lines are outbred for few generations, which we call the Hybrid Swarm. The Hybrid Swarm approach has likely remained underutilized relative to pre-sequenced inbred lines due to the costs of genome-wide genotyping. To reduce sequencing costs and make the Hybrid Swarm approach feasible, we developed a computational pipeline that reconstructs accurate whole genomes from ultra-low-coverage (0.05X) sequence data in Hybrid Swarm populations derived from ancestors with phased haplotypes. We evaluate reconstructions using genetic variation from the Drosophila Genetic Reference Panel as well as variation from neutral simulations. We compared the power and precision of Genome-Wide Association Studies using the Hybrid Swarm, inbred lines, recombinant inbred lines (RILs), and highly outbred populations across a range of allele frequencies, effect sizes, and genetic architectures. Our simulations show that these different mapping panels vary in their power and precision, largely depending on the architecture of the trait. The Hybrid Swam and RILs outperform inbred lines for quantitative traits, but not for monogenic ones. Taken together, our results demonstrate the feasibility of the Hybrid Swarm as a cost-effective method of fine-scale genetic mapping.

Host plants and Wolbachia shape the population genetics of sympatric herbivore populations.

Changing climate and land-use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host-plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Sulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome-wide SNP data revealed two major genetic clusters-psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non-crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co-occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host-plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non-crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

Unique genetic signatures of local adaptation over space and time for diapause, an ecologically relevant complex trait, in Drosophila melanogaster.

Organisms living in seasonally variable environments utilize cues such as light and temperature to induce plastic responses, enabling them to exploit favorable seasons and avoid unfavorable ones. Local adapation can result in variation in seasonal responses, but the genetic basis and evolutionary history of this variation remains elusive. Many insects, including Drosophila melanogaster, are able to undergo an arrest of reproductive development (diapause) in response to unfavorable conditions. In D. melanogaster, the ability to diapause is more common in high latitude populations, where flies endure harsher winters, and in the spring, reflecting differential survivorship of overwintering populations. Using a novel hybrid swarm-based genome wide association study, we examined the genetic basis and evolutionary history of ovarian diapause. We exposed outbred females to different temperatures and day lengths, characterized ovarian development for over 2800 flies, and reconstructed their complete, phased genomes. We found that diapause, scored at two different developmental cutoffs, has modest heritability, and we identified hundreds of SNPs associated with each of the two phenotypes. Alleles associated with one of the diapause phenotypes tend to be more common at higher latitudes, but these alleles do not show predictable seasonal variation. The collective signal of many small-effect, clinally varying SNPs can plausibly explain latitudinal variation in diapause seen in North America. Alleles associated with diapause are segregating in Zambia, suggesting that variation in diapause relies on ancestral polymorphisms, and both pro- and anti-diapause alleles have experienced selection in North America. Finally, we utilized outdoor mesocosms to track diapause under natural conditions. We found that hybrid swarms reared outdoors evolved increased propensity for diapause in late fall, whereas indoor control populations experienced no such change. Our results indicate that diapause is a complex, quantitative trait with different evolutionary patterns across time and space.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses.

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Common structuring principles of the Drosophila melanogaster microbiome on a continental scale and between host and substrate.

The relative importance of host control, environmental effects and stochasticity in the assemblage of host-associated microbiomes is being debated. We analysed the microbiome among fly populations that were sampled across Europe by the European Drosophila Population Genomics Consortium (DrosEU). In order to better understand the structuring principles of the natural D. melanogaster microbiome, we combined environmental data on climate and food-substrate with dense genomic data on host populations and microbiome profiling. Food-substrate, temperature, and host population structure correlated with microbiome structure. Microbes, whose abundance was co-structured with host populations, also differed in abundance between flies and their substrate in an independent survey. This finding suggests common, host-related structuring principles of the microbiome on different spatial scales.

Phenotypic plasticity, but not adaptive tracking, underlies seasonal variation in post-cold hardening freeze tolerance of Drosophila melanogaster.

In temperate regions, an organism's ability to rapidly adapt to seasonally varying environments is essential for its survival. In response to seasonal changes in selection pressure caused by variation in temperature, humidity, and food availability, some organisms exhibit plastic changes in phenotype. In other cases, seasonal variation in selection pressure can rapidly increase the frequency of genotypes that offer survival or reproductive advantages under the current conditions. Little is known about the relative influences of plastic and genetic changes in short-lived organisms experiencing seasonal environmental fluctuations. Cold hardening is a seasonally relevant plastic response in which exposure to cool, but nonlethal, temperatures significantly increases the organism's ability to later survive at freezing temperatures. In the present study, we demonstrate seasonal variation in cold hardening in Drosophila melanogaster and test the extent to which plasticity and adaptive tracking underlie that seasonal variation. We measured the post-cold hardening freeze tolerance of flies from outdoor mesocosms over the summer, fall, and winter. We bred outdoor mesocosm-caught flies for two generations in the laboratory and matched each outdoor cohort to an indoor control cohort of similar genetic background. We cold hardened all flies under controlled laboratory conditions and then measured their post-cold hardening freeze tolerance. Comparing indoor and field-caught flies and their laboratory-reared G1 and G2 progeny allowed us to determine the roles of seasonal environmental plasticity, parental effects, and genetic changes on cold hardening. We also tested the relationship between cold hardening and other factors, including age, developmental density, food substrate, presence of antimicrobials, and supplementation with live yeast. We found strong plastic responses to a variety of field- and laboratory-based environmental effects, but no evidence of seasonally varying parental or genetic effects on cold hardening. We therefore conclude that seasonal variation in post-cold hardening freeze tolerance results from environmental influences and not genetic changes.

Spatiotemporal dynamics and genome-wide association genome-wide association analysis of desiccation tolerance in Drosophila melanogaster.

Water availability is a major environmental challenge to a variety of terrestrial organisms. In insects, desiccation tolerance varies predictably over spatial and temporal scales and is an important physiological determinant of fitness in natural populations. Here, we examine the dynamics of desiccation tolerance in North American populations of Drosophila melanogaster using: (a) natural populations sampled across latitudes and seasons; (b) experimental evolution in field mesocosms over seasonal time; (c) genome-wide associations to identify SNPs/genes associated with variation for desiccation tolerance; and (d) subsequent analysis of patterns of clinal/seasonal enrichment in existing pooled sequencing data of populations sampled in both North America and Australia. A cline in desiccation tolerance was observed, for which tolerance exhibited a positive association with latitude; tolerance also varied predictably with culture temperature, demonstrating a significant degree of thermal plasticity. Desiccation tolerance evolved rapidly in field mesocosms, although only males showed differences in desiccation tolerance between spring and autumn collections from natural populations. Water loss rates did not vary significantly among latitudinal or seasonal populations; however, changes in metabolic rates during prolonged exposure to dry conditions are consistent with increased tolerance in higher latitude populations. Genome-wide associations in a panel of inbred lines identified twenty-five SNPs in twenty-one loci associated with sex-averaged desiccation tolerance, but there is no robust signal of spatially varying selection on genes associated with desiccation tolerance. Together, our results suggest that desiccation tolerance is a complex and important fitness component that evolves rapidly and predictably in natural populations.

Rapid seasonal evolution in innate immunity of wild Drosophila melanogaster.

Understanding the rate of evolutionary change and the genetic architecture that facilitates rapid adaptation is a current challenge in evolutionary biology. Comparative studies show that genes with immune function are among the most rapidly evolving genes across a range of taxa. Here, we use immune defence in natural populations of Drosophila melanogaster to understand the rate of evolution in natural populations and the genetics underlying rapid change. We probed the immune system using the natural pathogens Enterococcus faecalis and Providencia rettgeri to measure post-infection survival and bacterial load of wild D. melanogaster populations collected across seasonal time along a latitudinal transect along eastern North America (Massachusetts, Pennsylvania and Virginia). There are pronounced and repeatable changes in the immune response over the approximately 10 generations between spring and autumn collections, with a significant but less distinct difference observed among geographical locations. Genes with known immune function are not enriched among alleles that cycle with seasonal time, but the immune function of a subset of seasonally cycling alleles in immune genes was tested using reconstructed outbred populations. We find that flies containing seasonal alleles in Thioester-containing protein 3 (Tep3) have different functional responses to infection and that epistatic interactions among seasonal Tep3 and Drosomycin-like 6 (Dro6) alleles underlie the immune phenotypes observed in natural populations. This rapid, cyclic response to seasonal environmental pressure broadens our understanding of the complex ecological and genetic interactions determining the evolution of immune defence in natural populations.

Seasonally fluctuating selection can maintain polymorphism at many loci via segregation lift.

Most natural populations are affected by seasonal changes in temperature, rainfall, or resource availability. Seasonally fluctuating selection could potentially make a large contribution to maintaining genetic polymorphism in populations. However, previous theory suggests that the conditions for multilocus polymorphism are restrictive. Here, we explore a more general class of models with multilocus seasonally fluctuating selection in diploids. In these models, the multilocus genotype is mapped to fitness in two steps. The first mapping is additive across loci and accounts for the relative contributions of heterozygous and homozygous loci-that is, dominance. The second step uses a nonlinear fitness function to account for the strength of selection and epistasis. Using mathematical analysis and individual-based simulations, we show that stable polymorphism at many loci is possible if currently favored alleles are sufficiently dominant. This general mechanism, which we call "segregation lift," requires seasonal changes in dominance, a phenomenon that may arise naturally in situations with antagonistic pleiotropy and seasonal changes in the relative importance of traits for fitness. Segregation lift works best under diminishing-returns epistasis, is not affected by problems of genetic load, and is robust to differences in parameters across loci and seasons. Under segregation lift, loci can exhibit conspicuous seasonal allele-frequency fluctuations, but often fluctuations may be small and hard to detect. An important direction for future work is to formally test for segregation lift in empirical data and to quantify its contribution to maintaining genetic variation in natural populations.

Using NextRAD sequencing to infer movement of herbivores among host plants.

Herbivores often move among spatially interspersed host plants, tracking high-quality resources through space and time. This dispersal is of particular interest for vectors of plant pathogens. Existing molecular tools to track such movement have yielded important insights, but often provide insufficient genetic resolution to infer spread at finer spatiotemporal scales. Here, we explore the use of Nextera-tagmented reductively-amplified DNA (NextRAD) sequencing to infer movement of a highly-mobile winged insect, the potato psyllid (Bactericera cockerelli), among host plants. The psyllid vectors the pathogen that causes zebra chip disease in potato (Solanum tuberosum), but understanding and managing the spread of this pathogen is limited by uncertainty about the insect's host plant(s) outside of the growing season. We identified 1,978 polymorphic loci among psyllids separated spatiotemporally on potato or in patches of bittersweet nightshade (S. dulcumara), a weedy plant proposed to be the source of potato-colonizing psyllids. A subset of the psyllids on potato exhibited genetic similarity to insects on nightshade, consistent with regular movement between these two host plants. However, a second subset of potato-collected psyllids was genetically distinct from those collected on bittersweet nightshade; this suggests that a currently unrecognized source, i.e., other nightshade patches or a third host-plant species, could be contributing to psyllid populations in potato. Oftentimes, dispersal of vectors of pathogens must be tracked at a fine scale in order to understand, predict, and manage disease spread. We demonstrate that emerging sequencing technologies that detect genome-wide SNPs of a vector can be used to infer such localized movement.

The Genetic Basis of Natural Variation in Drosophila (Diptera: Drosophilidae) Virgin Egg Retention.

Drosophila melanogaster is able to thrive in harsh northern climates through adaptations in life-history traits and physiological mechanisms that allow for survival through the winter. We examined the genetic basis of natural variation in one such trait, female virgin egg retention, which was previously shown to vary clinally and seasonally. To further our understanding of the genetic basis and evolution of virgin egg retention, we performed a genome-wide association study (GWAS) using the previously sequenced Drosophila Genetic Reference Panel (DGRP) mapping population. We found 29 single nucleotide polymorphisms (SNPs) associated with virgin egg retention and assayed 6 available mutant lines, each harboring a mutation in a candidate gene, for effects on egg retention time. We found that four out of the six mutant lines had defects in egg retention time as compared with the respective controls: mun, T48, Mes-4, and Klp67A Surprisingly, none of these genes has a recognized role in ovulation control, but three of the four genes have known effects on fertility or have high expression in the ovaries. We also found that the SNP set associated with egg retention time was enriched for clinal SNPs. The majority of clinal SNPs had alleles associated with longer egg retention present at higher frequencies in higher latitudes. Our results support previous studies that show higher frequency of long retention times at higher latitude, providing evidence for the adaptive value of virgin egg-retention.

Eco-evolutionary Biology: Feeding and Feedback Loops.

A dogma in ecology and evolution holds that the environment is an extrinsic force that is not, in turn, shaped by the adaptive evolution of species. Recent work on stickleback life history, community ecology and speciation challenges this dogma.