Evolutionary Insights from a Large-Scale Survey of Population-Genomic Variation.

The field of genomics has ushered in new methods for studying molecular-genetic variation in natural populations. However, most population-genomic studies still rely on small sample sizes (typically, <100 individuals) from single time points, leaving considerable uncertainties with respect to the behavior of relatively young (and rare) alleles and, owing to the large sampling variance of measures of variation, to the specific gene targets of unusually strong selection. Genomic sequences of ∼1,700 haplotypes distributed over a 10-year period from a natural population of the microcrustacean Daphnia pulex reveal evolutionary-genomic features at a refined scale, including previously hidden information on the behavior of rare alleles predicted by recent theory. Background selection, resulting from the recurrent introduction of deleterious alleles, appears to strongly influence the dynamics of neutral alleles, inducing indirect negative selection on rare variants and positive selection on common variants. Temporally fluctuating selection increases the persistence of nonsynonymous alleles with intermediate frequencies, while reducing standing levels of variation at linked silent sites. Combined with the results from an equally large metapopulation survey of the study species, classes of genes that are under strong positive selection can now be confidently identified in this key model organism. Most notable among rapidly evolving Daphnia genes are those associated with ribosomes, mitochondrial functions, sensory systems, and lifespan determination.

Environment and shipping drive environmental DNA beta-diversity among commercial ports.

The spread of nonindigenous species by shipping is a large and growing global problem that harms coastal ecosystems and economies and may blur coastal biogeographical patterns. This study coupled eukaryotic environmental DNA (eDNA) metabarcoding with dissimilarity regression to test the hypothesis that ship-borne species spread homogenizes port communities. We first collected and metabarcoded water samples from ports in Europe, Asia, Australia and the Americas. We then calculated community dissimilarities between port pairs and tested for effects of environmental dissimilarity, biogeographical region and four alternative measures of ship-borne species transport risk. We predicted that higher shipping between ports would decrease community dissimilarity, that the effect of shipping would be small compared to that of environment dissimilarity and shared biogeography, and that more complex shipping risk metrics (which account for ballast water and stepping-stone spread) would perform better. Consistent with our hypotheses, community dissimilarities increased significantly with environmental dissimilarity and, to a lesser extent, decreased with ship-borne species transport risks, particularly if the ports had similar environments and stepping-stone risks were considered. Unexpectedly, we found no clear effect of shared biogeography, and that risk metrics incorporating estimates of ballast discharge did not offer more explanatory power than simpler traffic-based risks. Overall, we found that shipping homogenizes eukaryotic communities between ports in predictable ways, which could inform improvements in invasive species policy and management. We demonstrated the usefulness of eDNA metabarcoding and dissimilarity regression for disentangling the drivers of large-scale biodiversity patterns. We conclude by outlining logistical considerations and recommendations for future studies using this approach.

The road not taken: Evolution of tetrodotoxin resistance in the Sierra garter snake (Thamnophis couchii) by a path less travelled.

The repeated evolution of tetrodotoxin (TTX) resistance provides a model for testing hypotheses about the mechanisms of convergent evolution. This poison is broadly employed as a potent antipredator defence, blocking voltage-gated sodium channels (Nav ) in muscles and nerves, paralysing and sometimes killing predators. Resistance in taxa bearing this neurotoxin and a few predators appears to come from convergent replacements in specific Nav residues that interact with TTX. This stereotyped genetic response suggests molecular and phenotypic evolution may be constrained and predictable. Here, we investigate the extent of mechanistic convergence in garter snakes (Thamnophis) that prey on TTX-bearing newts (Taricha) by examining the physiological and genetic basis of TTX resistance in the Sierra garter snake (Th. couchii). We characterize variation in this predatory adaptation across populations at several biological scales: whole-animal TTX resistance; skeletal muscle resistance; functional genetic variation in three Nav encoding loci; and levels of gene expression for one of these loci. We found Th. couchii possess extensive geographical variation in resistance at the whole-animal and skeletal muscle levels. As in other Thamnophis, resistance at both levels is highly correlated, suggesting convergence across the biological levels linking organism to organ. However, Th. couchii shows no functional variation in Nav loci among populations or difference in candidate gene expression. Local variation in TTX resistance in Th. couchii cannot be explained by the same relationship between genotype and phenotype seen in other taxa. Thus, historical contingencies may lead different species of Thamnophis down alternative routes to local adaptation.

Global invasion history of the agricultural pest butterfly Pieris rapae revealed with genomics and citizen science.

The small cabbage white butterfly, Pieris rapae, is a major agricultural pest of cruciferous crops and has been introduced to every continent except South America and Antarctica as a result of human activities. In an effort to reconstruct the near-global invasion history of P. rapae, we developed a citizen science project, the "Pieris Project," and successfully amassed thousands of specimens from 32 countries worldwide. We then generated and analyzed nuclear (double-digest restriction site-associated DNA fragment procedure [ddRAD]) and mitochondrial DNA sequence data for these samples to reconstruct and compare different global invasion history scenarios. Our results bolster historical accounts of the global spread and timing of P. rapae introductions. We provide molecular evidence supporting the hypothesis that the ongoing divergence of the European and Asian subspecies of P. rapae (∼1,200 y B.P.) coincides with the diversification of brassicaceous crops and the development of human trade routes such as the Silk Route (Silk Road). The further spread of P. rapae over the last ∼160 y was facilitated by human movement and trade, resulting in an almost linear series of at least 4 founding events, with each introduced population going through a severe bottleneck and serving as the source for the next introduction. Management efforts of this agricultural pest may need to consider the current existence of multiple genetically distinct populations. Finally, the international success of the Pieris Project demonstrates the power of the public to aid scientists in collections-based research addressing important questions in invasion biology, and in ecology and evolutionary biology more broadly.

Climate-mediated hybrid zone movement revealed with genomics, museum collection, and simulation modeling.

Climate-mediated changes in hybridization will dramatically alter the genetic diversity, adaptive capacity, and evolutionary trajectory of interbreeding species. Our ability to predict the consequences of such changes will be key to future conservation and management decisions. Here we tested through simulations how recent warming (over the course of a 32-y period) is affecting the geographic extent of a climate-mediated developmental threshold implicated in maintaining a butterfly hybrid zone (Papilio glaucus and Papilio canadensis; Lepidoptera: Papilionidae). These simulations predict a 68-km shift of this hybrid zone. To empirically test this prediction, we assessed genetic and phenotypic changes using contemporary and museum collections and document a 40-km northward shift of this hybrid zone. Interactions between the two species appear relatively unchanged during hybrid zone movement. We found no change in the frequency of hybridization, and regions of the genome that experience little to no introgression moved largely in concert with the shifting hybrid zone. Model predictions based on climate scenarios predict this hybrid zone will continue to move northward, but with substantial spatial heterogeneity in the velocity (55-144 km/1 °C), shape, and contiguity of movement. Our findings suggest that the presence of nonclimatic barriers (e.g., genetic incompatibilities) and/or nonlinear responses to climatic gradients may preserve species boundaries as the species shift. Further, we show that variation in the geography of hybrid zone movement could result in evolutionary responses that differ for geographically distinct populations spanning hybrid zones, and thus have implications for the conservation and management of genetic diversity.

Evolutionary transition from blood feeding to obligate nonbiting in a mosquito.

The spread of blood-borne pathogens by mosquitoes relies on their taking a blood meal; if there is no bite, there is no disease transmission. Although many species of mosquitoes never take a blood meal, identifying genes that distinguish blood feeding from obligate nonbiting is hampered by the fact that these different lifestyles occur in separate, genetically incompatible species. There is, however, one unique extant species with populations that share a common genetic background but blood feed in one region and are obligate nonbiters in the rest of their range: Wyeomyia smithii Contemporary blood-feeding and obligate nonbiting populations represent end points of divergence between fully interfertile southern and northern populations. This divergence has undoubtedly resulted in genetic changes that are unrelated to blood feeding, and the challenge is to winnow out the unrelated genetic factors to identify those related specifically to the evolutionary transition from blood feeding to obligate nonbiting. Herein, we determine differential gene expression resulting from directional selection on blood feeding within a polymorphic population to isolate genetic differences between blood feeding and obligate nonbiting. We show that the evolution of nonbiting has resulted in a greatly reduced metabolic investment compared with biting populations, a greater reliance on opportunistic metabolic pathways, and greater reliance on visual rather than olfactory sensory input. W. smithii provides a unique starting point to determine if there are universal nonbiting genes in mosquitoes that could be manipulated as a means to control vector-borne disease.

The roles of vicariance and isolation by distance in shaping biotic diversification across an ancient archipelago: evidence from a Seychelles caecilian amphibian.

BACKGROUND: Island systems offer excellent opportunities for studying the evolutionary histories of species by virtue of their restricted size and easily identifiable barriers to gene flow. However, most studies investigating evolutionary patterns and processes shaping biotic diversification have focused on more recent (emergent) rather than ancient oceanic archipelagos. Here, we focus on the granitic islands of the Seychelles, which are unusual among island systems because they have been isolated for a long time and are home to a monophyletic radiation of caecilian amphibians that has been separated from its extant sister lineage for ca. 65-62 Ma. We selected the most widespread Seychelles caecilian species, Hypogeophis rostratus, to investigate intraspecific morphological and genetic (mitochondrial and nuclear) variation across the archipelago (782 samples from nine islands) to identify patterns and test processes that shaped their evolutionary history within the Seychelles. RESULTS: Overall a signal of strong geographic structuring with distinct northern- and southern-island clusters were identified across all datasets. We suggest that these distinct groups have been isolated for ca. 1.26 Ma years without subsequent migration between them. Populations from the somewhat geographically isolated island of Frégate showed contrasting relationships to other islands based on genetic and morphological data, clustering alternatively with northern-island (genetic) and southern-island (morphological) populations. CONCLUSIONS: Although variation in H. rostratus across the Seychelles is explained more by isolation-by-distance than by adaptation, the genetic-morphological incongruence for affinities of Frégate H. rostratus might be caused by local adaptation over-riding the signal from their vicariant history. Our findings highlight the need of integrative approaches to investigate fine-scale geographic structuring to uncover underlying diversity and to better understand evolutionary processes on ancient, continental islands.

A comprehensive epigenomic analysis of phenotypically distinguishable, genetically identical female and male Daphnia pulex.

BACKGROUND: Daphnia species reproduce by cyclic parthenogenesis involving both sexual and asexual reproduction. The sex of the offspring is environmentally determined and mediated via endocrine signalling by the mother. Interestingly, male and female Daphnia can be genetically identical, yet display large differences in behaviour, morphology, lifespan and metabolic activity. Our goal was to integrate multiple omics datasets, including gene expression, splicing, histone modification and DNA methylation data generated from genetically identical female and male Daphnia pulex under controlled laboratory settings with the aim of achieving a better understanding of the underlying epigenetic factors that may contribute to the phenotypic differences observed between the two genders. RESULTS: In this study we demonstrate that gene expression level is positively correlated with increased DNA methylation, and histone H3 trimethylation at lysine 4 (H3K4me3) at predicted promoter regions. Conversely, elevated histone H3 trimethylation at lysine 27 (H3K27me3), distributed across the entire transcript length, is negatively correlated with gene expression level. Interestingly, male Daphnia are dominated with epigenetic modifications that globally promote elevated gene expression, while female Daphnia are dominated with epigenetic modifications that reduce gene expression globally. For examples, CpG methylation (positively correlated with gene expression level) is significantly higher in almost all differentially methylated sites in male compared to female Daphnia. Furthermore, H3K4me3 modifications are higher in male compared to female Daphnia in more than 3/4 of the differentially regulated promoters. On the other hand, H3K27me3 is higher in female compared to male Daphnia in more than 5/6 of differentially modified sites. However, both sexes demonstrate roughly equal number of genes that are up-regulated in one gender compared to the other sex. Since, gene expression analyses typically assume that most genes are expressed at equal level among samples and different conditions, and thus cannot detect global changes affecting most genes. CONCLUSIONS: The epigenetic differences between male and female in Daphnia pulex are vast and dominated by changes that promote elevated gene expression in male Daphnia. Furthermore, the differences observed in both gene expression changes and epigenetic modifications between the genders relate to pathways that are physiologically relevant to the observed phenotypic differences.

The geographic mosaic in parallel: Matching patterns of newt tetrodotoxin levels and snake resistance in multiple predator-prey pairs.

The Geographic Mosaic Theory of Coevolution predicts that coevolutionary arms races will vary over time and space because of the diverse ecological settings and population histories of interacting species across the landscape. Thus, understanding coevolution may require investigating broad sets of populations sampled across the range of the interaction. In addition, comparing coevolutionary dynamics between similar systems may reveal the importance of specific factors that structure coevolution. Here, we examine geographic patterns of prey traits and predator traits in the relatively unstudied interaction between the Sierra garter snake (Thamnophis couchii) and sympatric prey, the rough-skinned newt (Taricha granulosa), Sierra newt (Ta. sierrae) and California newt (Ta. torosa). This system parallels, in space and phenotypes, a classic example of coevolution between predatory common garter snakes (Th. sirtalis) and their toxic newt prey exhibiting hotspots of newt tetrodotoxin (TTX) levels and matching snake TTX resistance. We quantified prey and predator traits from hundreds of individuals across their distributions, and functional trait matching at sympatric sites. We show strong regional patterns of trait covariation across the shared ranges of Th. couchii and newt prey. Traits differ significantly among localities, with lower newt TTX levels and snake TTX resistance at the northern latitudes, and higher TTX levels and snake resistance at southern latitudes. Newts and snakes in northern populations show the highest degree of functional trait matching despite possessing the least extreme traits. Conversely, newts and snakes in southern populations show the greatest mismatch despite possessing exaggerated traits, with some snakes so resistant to TTX they would be unaffected by any sympatric newt. Nevertheless, individual variation was substantial, and appears to offer the opportunity for continued reciprocal selection in most populations. Overall, the three species of newts appear to be engaged in a TTX-mediated arms race with Th. couchii. These patterns are congruent with those seen between newts and Th. sirtalis, including the same latitudinal gradient in trait covariation, and the potential 'escape' from the arms race by snake predators. Such concordance in broad scale patterns across two distinct systems suggests common phenomena might structure geographic mosaics in similar ways.

Early transcriptional response pathways in Daphnia magna are coordinated in networks of crustacean-specific genes.

Natural habitats are exposed to an increasing number of environmental stressors that cause important ecological consequences. However, the multifarious nature of environmental change, the strength and the relative timing of each stressor largely limit our understanding of biological responses to environmental change. In particular, early response to unpredictable environmental change, critical to survival and fitness in later life stages, is largely uncharacterized. Here, we characterize the early transcriptional response of the keystone species Daphnia magna to twelve environmental perturbations, including biotic and abiotic stressors. We first perform a differential expression analysis aimed at identifying differential regulation of individual genes in response to stress. This preliminary analysis revealed that a few individual genes were responsive to environmental perturbations and they were modulated in a stressor and genotype-specific manner. Given the limited number of differentially regulated genes, we were unable to identify pathways involved in stress response. Hence, to gain a better understanding of the genetic and functional foundation of tolerance to multiple environmental stressors, we leveraged the correlative nature of networks and performed a weighted gene co-expression network analysis. We discovered that approximately one-third of the Daphnia genes, enriched for metabolism, cell signalling and general stress response, drives transcriptional early response to environmental stress and it is shared among genetic backgrounds. This initial response is followed by a genotype- and/or condition-specific transcriptional response with a strong genotype-by-environment interaction. Intriguingly, genotype- and condition-specific transcriptional response is found in genes not conserved beyond crustaceans, suggesting niche-specific adaptation.

The molecular evolutionary dynamics of oxidative phosphorylation (OXPHOS) genes in Hymenoptera.

BACKGROUND: The primary energy-producing pathway in eukaryotic cells, the oxidative phosphorylation (OXPHOS) system, comprises proteins encoded by both mitochondrial and nuclear genes. To maintain the function of the OXPHOS system, the pattern of substitutions in mitochondrial and nuclear genes may not be completely independent. It has been suggested that slightly deleterious substitutions in mitochondrial genes are compensated by substitutions in the interacting nuclear genes due to positive selection. Among the four largest insect orders, Coleoptera (beetles), Hymenoptera (sawflies, wasps, ants, and bees), Diptera (midges, mosquitoes, and flies) and Lepidoptera (moths and butterflies), the mitochondrial genes of Hymenoptera exhibit an exceptionally high amino acid substitution rate while the evolution of nuclear OXPHOS genes is largely unknown. Therefore, Hymenoptera is an excellent model group for testing the hypothesis of positive selection driving the substitution rate of nuclear OXPHOS genes. In this study, we report the evolutionary rate of OXPHOS genes in Hymenoptera and test for evidence of positive selection in nuclear OXPHOS genes of Hymenoptera. RESULTS: Our analyses revealed that the amino acid substitution rate of mitochondrial and nuclear OXPHOS genes in Hymenoptera is higher than that in other studied insect orders. In contrast, the amino acid substitution rate of non-OXPHOS genes in Hymenoptera is lower than the rate in other insect orders. Overall, we found the dN/dS ratio of the nuclear OXPHOS genes to be higher in Hymenoptera than in other insect orders. However, nuclear OXPHOS genes with high dN/dS ratio did not always exhibit a high amino acid substitution rate. Using branch-site and site model tests, we identified various codon sites that evolved under positive selection in nuclear OXPHOS genes. CONCLUSIONS: Our results showed that nuclear OXPHOS genes in Hymenoptera are evolving faster than the genes in other three insect orders. The branch test suggested that while some nuclear OXPHOS genes in Hymenoptera show a signature of positive selection, the pattern is not consistent across all nuclear OXPHOS genes. As only few codon sites were under positive selection, we suggested that positive selection might not be the only factor contributing to the rapid evolution of nuclear OXPHOS genes in Hymenoptera.

Reproductive isolation and environmental adaptation shape the phylogeography of mountain pine beetle (Dendroctonus ponderosae).

Chromosomal rearrangement can be an important mechanism driving population differentiation and incipient speciation. In the mountain pine beetle (MPB, Dendroctonus ponderosae), deletions on the Y chromosome that are polymorphic among populations are associated with reproductive incompatibility. Here, we used RAD sequencing across the entire MPB range in western North America to reveal the extent of the phylogeographic differences between Y haplotypes compared to autosomal and X-linked loci. Clustering and geneflow analyses revealed three distinct Y haplogroups geographically positioned within and on either side of the Great Basin Desert. Despite close geographic proximity between populations on the boundaries of each Y haplogroup, there was extremely low Y haplogroup mixing among populations, and gene flow on the autosomes was reduced across Y haplogroup boundaries. These results are consistent with a previous study suggesting that independent degradation of a recently evolved neo-Y chromosome in previously isolated populations causes male sterility or inviability among Y haplotype lineages. Phylogeographic results supported historic contraction of MPB into three separate Pleistocene glacial refugia followed by postglacial range expansion and secondary contact. Distinct sets of SNPs were statistically associated with environmental data among the most genetically distinct sets of geographic populations. This finding suggests that the process of adaptation to local climatic conditions is influenced by population genetic structure, with evidence for largely independent evolution in the most genetically isolated Y haplogroup.

Patterns of divergence across the geographic and genomic landscape of a butterfly hybrid zone associated with a climatic gradient.

Hybrid zones are a valuable tool for studying the process of speciation and for identifying the genomic regions undergoing divergence and the ecological (extrinsic) and nonecological (intrinsic) factors involved. Here, we explored the genomic and geographic landscape of divergence in a hybrid zone between Papilio glaucus and Papilio canadensis. Using a genome scan of 28,417 ddRAD SNPs, we identified genomic regions under possible selection and examined their distribution in the context of previously identified candidate genes for ecological adaptations. We showed that differentiation was genomewide, including multiple candidate genes for ecological adaptations, particularly those involved in seasonal adaptation and host plant detoxification. The Z chromosome and four autosomes showed a disproportionate amount of differentiation, suggesting genes on these chromosomes play a potential role in reproductive isolation. Cline analyses of significantly differentiated genomic SNPs, and of species-diagnostic genetic markers, showed a high degree of geographic coincidence (81%) and concordance (80%) and were associated with the geographic distribution of a climate-mediated developmental threshold (length of the growing season). A relatively large proportion (1.3%) of the outliers for divergent selection were not associated with candidate genes for ecological adaptations and may reflect the presence of previously unrecognized intrinsic barriers between these species. These results suggest that exogenous (climate-mediated) and endogenous (unknown) clines may have become coupled and act together to reinforce reproductive isolation. This approach of assessing divergence across both the genomic and geographic landscape can provide insight about the interplay between the genetic architecture of reproductive isolation and endogenous and exogenous selection.

Environmental DNA metabarcoding: Transforming how we survey animal and plant communities.

The genomic revolution has fundamentally changed how we survey biodiversity on earth. High-throughput sequencing ("HTS") platforms now enable the rapid sequencing of DNA from diverse kinds of environmental samples (termed "environmental DNA" or "eDNA"). Coupling HTS with our ability to associate sequences from eDNA with a taxonomic name is called "eDNA metabarcoding" and offers a powerful molecular tool capable of noninvasively surveying species richness from many ecosystems. Here, we review the use of eDNA metabarcoding for surveying animal and plant richness, and the challenges in using eDNA approaches to estimate relative abundance. We highlight eDNA applications in freshwater, marine and terrestrial environments, and in this broad context, we distill what is known about the ability of different eDNA sample types to approximate richness in space and across time. We provide guiding questions for study design and discuss the eDNA metabarcoding workflow with a focus on primers and library preparation methods. We additionally discuss important criteria for consideration of bioinformatic filtering of data sets, with recommendations for increasing transparency. Finally, looking to the future, we discuss emerging applications of eDNA metabarcoding in ecology, conservation, invasion biology, biomonitoring, and how eDNA metabarcoding can empower citizen science and biodiversity education.

Bisulfite Sequencing with Daphnia Highlights a Role for Epigenetics in Regulating Stress Response to Microcystis through Preferential Differential Methylation of Serine and Threonine Amino Acids.

Little is known about the influence that environmental stressors may have on genome-wide methylation patterns, and to what extent epigenetics may be involved in environmental stress response. Yet, studies of methylation patterns under stress could provide crucial insights on stress response and toxicity pathways. Here, we focus on genome-wide methylation patterns in the microcrustacean Daphnia magna, a model organism in ecotoxicology and risk assessment, exposed to the toxic cyanobacterium Microcystis aeruginosa. Bisulfite sequencing of exposed and control animals highlighted differential methylation patterns in Daphnia upon exposure to Microcystis primarily in exonic regions. These patterns are enriched for serine/threonine amino acid codons and genes related to protein synthesis, transport and degradation. Furthermore, we observed that genes with differential methylation corresponded well with genes susceptible to alternative splicing in response to Microcystis stress. Overall, our results suggest a complex mechanistic response in Daphnia characterized by interactions between DNA methylation and gene regulation mechanisms. These results underscore that DNA methylation is modulated by environmental stress and can also be an integral part of the toxicity response in our study species.

Genome-wide profiling of 24 hr diel rhythmicity in the water flea, Daphnia pulex: network analysis reveals rhythmic gene expression and enhances functional gene annotation.

BACKGROUND: Marine and freshwater zooplankton exhibit daily rhythmic patterns of behavior and physiology which may be regulated directly by the light:dark (LD) cycle and/or a molecular circadian clock. One of the best-studied zooplankton taxa, the freshwater crustacean Daphnia, has a 24 h diel vertical migration (DVM) behavior whereby the organism travels up and down through the water column daily. DVM plays a critical role in resource tracking and the behavioral avoidance of predators and damaging ultraviolet radiation. However, there is little information at the transcriptional level linking the expression patterns of genes to the rhythmic physiology/behavior of Daphnia. RESULTS: Here we analyzed genome-wide temporal transcriptional patterns from Daphnia pulex collected over a 44 h time period under a 12:12 LD cycle (diel) conditions using a cosine-fitting algorithm. We used a comprehensive network modeling and analysis approach to identify novel co-regulated rhythmic genes that have similar network topological properties and functional annotations as rhythmic genes identified by the cosine-fitting analyses. Furthermore, we used the network approach to predict with high accuracy novel gene-function associations, thus enhancing current functional annotations available for genes in this ecologically relevant model species. Our results reveal that genes in many functional groupings exhibit 24 h rhythms in their expression patterns under diel conditions. We highlight the rhythmic expression of immunity, oxidative detoxification, and sensory process genes. We discuss differences in the chronobiology of D. pulex from other well-characterized terrestrial arthropods. CONCLUSIONS: This research adds to a growing body of literature suggesting the genetic mechanisms governing rhythmicity in crustaceans may be divergent from other arthropod lineages including insects. Lastly, these results highlight the power of using a network analysis approach to identify differential gene expression and provide novel functional annotation.

Temporal genetic stability in natural populations of the waterflea Daphnia magna in response to strong selection pressure.

Studies monitoring changes in genetic diversity and composition through time allow a unique understanding of evolutionary dynamics and persistence of natural populations. However, such studies are often limited to species with short generation times that can be propagated in the laboratory or few exceptional cases in the wild. Species that produce dormant stages provide powerful models for the reconstruction of evolutionary dynamics in the natural environment. A remaining open question is to what extent dormant egg banks are an unbiased representation of populations and hence of the species' evolutionary potential, especially in the presence of strong environmental selection. We address this key question using the water flea Daphnia magna, which produces dormant stages that accumulate in biological archives over time. We assess temporal genetic stability in three biological archives, previously used in resurrection ecology studies showing adaptive evolutionary responses to rapid environmental change. We show that neutral genetic diversity does not decline with the age of the population and it is maintained in the presence of strong selection. In addition, by comparing temporal genetic stability in hatched and unhatched populations from the same biological archive, we show that dormant egg banks can be consulted to obtain a reliable measure of genetic diversity over time, at least in the multidecadal time frame studied here. The stability of neutral genetic diversity through time is likely mediated by the buffering effect of the resting egg bank.

Parallel evolution of tetrodotoxin resistance in three voltage-gated sodium channel genes in the garter snake Thamnophis sirtalis.

Members of a gene family expressed in a single species often experience common selection pressures. Consequently, the molecular basis of complex adaptations may be expected to involve parallel evolutionary changes in multiple paralogs. Here, we use bacterial artificial chromosome library scans to investigate the evolution of the voltage-gated sodium channel (Nav) family in the garter snake Thamnophis sirtalis, a predator of highly toxic Taricha newts. Newts possess tetrodotoxin (TTX), which blocks Nav's, arresting action potentials in nerves and muscle. Some Thamnophis populations have evolved resistance to extremely high levels of TTX. Previous work has identified amino acid sites in the skeletal muscle sodium channel Nav1.4 that confer resistance to TTX and vary across populations. We identify parallel evolution of TTX resistance in two additional Nav paralogs, Nav1.6 and 1.7, which are known to be expressed in the peripheral nervous system and should thus be exposed to ingested TTX. Each paralog contains at least one TTX-resistant substitution identical to a substitution previously identified in Nav1.4. These sites are fixed across populations, suggesting that the resistant peripheral nerves antedate resistant muscle. In contrast, three sodium channels expressed solely in the central nervous system (Nav1.1-1.3) showed no evidence of TTX resistance, consistent with protection from toxins by the blood-brain barrier. We also report the exon-intron structure of six Nav paralogs, the first such analysis for snake genes. Our results demonstrate that the molecular basis of adaptation may be both repeatable across members of a gene family and predictable based on functional considerations.

Conserved transcriptional responses to cyanobacterial stressors are mediated by alternate regulation of paralogous genes in Daphnia.

Despite a significant increase in genomic data, our knowledge of gene functions and their transcriptional responses to environmental stimuli remains limited. Here, we use the model keystone species Daphnia pulex to study environmental responses of genes in the context of their gene family history to better understand the relationship between genome structure and gene function in response to environmental stimuli. Daphnia were exposed to five different treatments, each consisting of a diet supplemented with one of five cyanobacterial species, and a control treatment consisting of a diet of only green algae. Differential gene expression profiles of Daphnia exposed to each of these five cyanobacterial species showed that genes with known functions are more likely to be shared by different expression profiles, whereas genes specific to the lineage of Daphnia are more likely to be unique to a given expression profile. Furthermore, while only a small number of nonlineage-specific genes were conserved across treatment type, there was a high degree of overlap in expression profiles at the functional level. The conservation of functional responses across the different cyanobacterial treatments can be attributed to the treatment-specific expression of different paralogous genes within the same gene family. Comparison with available gene expression data in the literature suggests differences in nutritional composition in diets with cyanobacterial species compared to diets of green algae as a primary driver for cyanobacterial effects on Daphnia. We conclude that conserved functional responses in Daphnia across different cyanobacterial treatments are mediated through alternate regulation of paralogous gene families.

Constraint shapes convergence in tetrodotoxin-resistant sodium channels of snakes.

Natural selection often produces convergent changes in unrelated lineages, but the degree to which such adaptations occur via predictable genetic paths is unknown. If only a limited subset of possible mutations is fixed in independent lineages, then it is clear that constraint in the production or function of molecular variants is an important determinant of adaptation. We demonstrate remarkably constrained convergence during the evolution of resistance to the lethal poison, tetrodotoxin, in six snake species representing three distinct lineages from around the globe. Resistance-conferring amino acid substitutions in a voltage-gated sodium channel, Na(v)1.4, are clustered in only two regions of the protein, and a majority of the replacements are confined to the same three positions. The observed changes represent only a small fraction of the experimentally validated mutations known to increase Na(v)1.4 resistance to tetrodotoxin. These results suggest that constraints resulting from functional tradeoffs between ion channel function and toxin resistance led to predictable patterns of evolutionary convergence at the molecular level. Our data are consistent with theoretical predictions and recent microcosm work that suggest a predictable path is followed during an adaptive walk along a mutational landscape, and that natural selection may be frequently constrained to produce similar genetic outcomes even when operating on independent lineages.