Population-specific responses in eastern oysters exposed to low salinity in the northern Gulf of Mexico.

Eastern oysters, Crassostrea virginica, are facing rapid environmental changes in the northern Gulf of Mexico and can respond to these changes via plasticity or evolution. Plastic responses can immediately buffer against environmental changes, although this buffering may impact the organism's ability to evolve in subsequent generations. While plasticity and evolution are not mutually exclusive, the relative contribution and interaction between them remains unclear. In this study, we investigated the roles of plastic and evolved responses of C. virginica acclimated to low salinity using a common garden experiment with four populations exposed to two salinities. We used three transcriptomic analyses (edgeR, PERMANOVA and WGCNA) combined with physiology data to identify the effect of genotype (population), environment (salinity) and the genotype-environment interaction on both whole-organism and molecular phenotypes. We demonstrate that variation in gene expression is mainly driven by population, with relatively small changes in response to salinity. In contrast, the morphology and physiology data reveal that salinity has a larger influence on oyster performance than the population of origin. All analyses lacked signatures of the genotype×environment interaction and, in contrast to previous studies, we found no evidence for population-specific responses to low salinity. However, individuals from the highest salinity estuary displayed highly divergent gene expression from that of other populations, which could potentially drive population-specific responses to other stressors. Our findings suggest that C. virginica largely rely on plasticity in physiology to buffer the effects of low salinity, but that these changes in physiology do not rely on large persistent changes in gene expression.

Constitutive gene expression differs in three brain regions important for cognition in neophobic and non-neophobic house sparrows (Passer domesticus).

Neophobia (aversion to new objects, food, and environments) is a personality trait that affects the ability of wildlife to adapt to new challenges and opportunities. Despite the ubiquity and importance of this trait, the molecular mechanisms underlying repeatable individual differences in neophobia in wild animals are poorly understood. We evaluated wild-caught house sparrows (Passer domesticus) for neophobia in the lab using novel object tests. We then selected a subset of neophobic and non-neophobic individuals (n = 3 of each, all females) and extracted RNA from four brain regions involved in learning, memory, threat perception, and executive function: striatum, caudal dorsomedial hippocampus, medial ventral arcopallium, and caudolateral nidopallium (NCL). Our analysis of differentially expressed genes (DEGs) used 11,889 gene regions annotated in the house sparrow reference genome for which we had an average of 25.7 million mapped reads/sample. PERMANOVA identified significant effects of brain region, phenotype (neophobic vs. non-neophobic), and a brain region by phenotype interaction. Comparing neophobic and non-neophobic birds revealed constitutive differences in DEGs in three of the four brain regions examined: hippocampus (12% of the transcriptome significantly differentially expressed), striatum (4%) and NCL (3%). DEGs included important known neuroendocrine mediators of learning, memory, executive function, and anxiety behavior, including serotonin receptor 5A, dopamine receptors 1, 2 and 5 (downregulated in neophobic birds), and estrogen receptor beta (upregulated in neophobic birds). These results suggest that some of the behavioral differences between phenotypes may be due to underlying gene expression differences in the brain. The large number of DEGs in neophobic and non-neophobic birds also implies that there are major differences in neural function between the two phenotypes that could affect a wide variety of behavioral traits beyond neophobia.

Comparative transcriptomics reveals altered species interaction between the bioeroding sponge Cliona varians and the coral Porites furcata under ocean acidification.

Bioeroding sponges interact and compete with corals on tropical reefs. Experimental studies have shown global change alters this biotic interaction, often in favour of the sponge. Ocean acidification in particular increases sponge bioerosion and reduces coral calcification, yet little is known about the molecular basis of these changes. We used RNA-Seq data to understand how acidification impacts the interaction between the bioeroding sponge, Cliona varians, and the coral, Porites furcata, at the transcriptomic level. Replicate sponge and coral genets were exposed to ambient (8.1 pH) and acidified (7.6 pH) conditions in isolation and in treatments where they were joined for 48 h. The coral had a small gene expression response (tens of transcripts) to the sponge, suggesting it does little at the transcriptomic level to deter sponge overgrowth. By contrast, the sponge differentially expressed 7320 transcripts in response to the coral under ambient conditions and 3707 transcripts in response to acidification. Overlap in the responses to acidification and the coral, 2500 transcripts expressed under both treatments, suggests a similar physiological response to both cues. The sponge expressed 50× fewer transcripts in response to the coral under acidification, suggesting energetic costs of bioerosion, and other cellular processes, are lower for sponges under acidification. Our results suggest how acidification drives ecosystem-level changes in the accretion/bioerosion balance on coral reefs. This shift is not only the result of changes to the thermodynamic balance of these chemical reactions but also the result of active physiological responses of organisms to each other and their abiotic environment.

Differential DNA methylation across environments has no effect on gene expression in the eastern oyster.

It has been hypothesized that environmentally induced changes to gene body methylation could facilitate adaptive transgenerational responses to changing environments. We compared patterns of global gene expression (Tag-seq) and gene body methylation (reduced representation bisulfite sequencing) in 80 eastern oysters Crassostrea virginica from six full-sib families, common gardened for 14 months at two sites in the northern Gulf of Mexico that differed in mean salinity. At the time of sampling, oysters from the two sites differed in mass by 60% and in parasite loads by nearly two orders of magnitude. They also differentially expressed 35% of measured transcripts. However, we observed differential methylation at only 1.4% of potentially methylated loci in comparisons between individuals from these different environments, and little correspondence between differential methylation and differential gene expression. Instead, methylation patterns were largely driven by genetic differences among families, with a PERMANOVA analysis indicating nearly a two orders of magnitude greater number of genes differentially methylated between families than between environments. An analysis of CpG observed/expected values (CpG O/E) across the C. virginica genome showed a distinct bimodal distribution, with genes from the first cluster showing the lower CpG O/E values, greater methylation and higher and more stable gene expression, while genes from the second cluster showed lower methylation, and lower and more variable gene expression. Taken together, the differential methylation results suggest that only a small portion of the C. virginica genome is affected by environmentally induced changes in methylation. At this point, there is little evidence to suggest that environmentally induced methylation states would play a leading role in regulating gene expression responses to new environments.

Tolerance of northern Gulf of Mexico eastern oysters to chronic warming at extreme salinities.

The eastern oyster, Crassostrea virginica, provides critical ecosystem services and supports valuable fishery and aquaculture industries in northern Gulf of Mexico (nGoM) subtropical estuaries where it is grown subtidally. Its upper critical thermal limit is not well defined, especially when combined with extreme salinities. The cumulative mortalities of the progenies of wild C. virginica from four nGoM estuaries differing in mean annual salinity, acclimated to low (4.0), moderate (20.0), and high (36.0) salinities at 28.9 °C (84 °F) and exposed to increasing target temperatures of 33.3 °C (92 °F), 35.6 °C (96 °F) or 37.8 °C (100 °F), were measured over a three-week period. Oysters of all stocks were the most sensitive to increasing temperatures at low salinity, dying quicker (i.e., lower median lethal time, LT50) than at the moderate and high salinities and resulting in high cumulative mortalities at all target temperatures. Oysters of all stocks at moderate salinity died the slowest with high cumulative mortalities only at the two highest temperatures. The F1 oysters from the more southern and hypersaline Upper Laguna Madre estuary were generally more tolerant to prolonged higher temperatures (higher LT50) than stocks originating from lower salinity estuaries, most notably at the highest salinity. Using the measured temperatures oysters were exposed to, 3-day median lethal Celsius degrees (LD50) were estimated for each stock at each salinity. The lowest 3-day LD50 (35.1-36.0 °C) for all stocks was calculated at a salinity of 4.0, while the highest 3-day LD50 (40.1-44.0 °C) was calculated at a salinity of 20.0.

Lack of genotype-by-environment interaction suggests limited potential for evolutionary changes in plasticity in the eastern oyster, Crassostrea virginica.

Eastern oysters in the northern Gulf of Mexico are facing rapid environmental changes and can respond to this change via plasticity or evolution. Plasticity can act as an immediate buffer against environmental change, but this buffering could impact the organism's ability to evolve in subsequent generations. While plasticity and evolution are not mutually exclusive, the relative contribution and interaction between them remains unclear. In this study, we investigate the roles of plastic and evolved responses to environmental variation and Perkinsus marinus infection in Crassostrea virginica by using a common garden experiment with 80 oysters from six families outplanted at two field sites naturally differing in salinity. We use growth data, P. marinus infection intensities, 3' RNA sequencing (TagSeq) and low-coverage whole-genome sequencing to identify the effect of genotype, environment and genotype-by-environment interaction on the oyster's response to site. As one of first studies to characterize the joint effects of genotype and environment on transcriptomic and morphological profiles in a natural setting, we demonstrate that C. virginica has a highly plastic response to environment and that this response is parallel among genotypes. We also find that genes responding to genotype have distinct and opposing profiles compared to genes responding to environment with regard to expression levels, Ka/Ks ratios and nucleotide diversity. Our findings suggest that C. virginica may be able to buffer the immediate impacts of future environmental changes by altering gene expression and physiology, but the lack of genetic variation in plasticity suggests limited capacity for evolved responses.

Evolutionary Change in the Eastern Oyster, Crassostrea Virginica, Following Low Salinity Exposure.

The presence of standing genetic variation will play a role in determining a population's capacity to adapt to environmentally relevant stressors. In the Gulf of Mexico, extreme climatic events and anthropogenic changes to local hydrology will expose productive oyster breeding grounds to stressful low salinity conditions. We identified genetic variation for performance under low salinity (due to the combined effects of low salinity and genetic load) using a single-generation selection experiment on larvae from two populations of the eastern oyster, Crassostrea virginica. We used pool-sequencing to test for allele frequency differences at 152 salinity-associated genes for larval families pre- and post-low salinity exposure. Our results have implications for how evolutionary change occurs during early life history stages at environmentally relevant salinities. Consistent with observations of high genetic load observed in oysters, we demonstrate evidence for purging of deleterious alleles at the larval stage in C. virginica. In addition, we observe increases in allele frequencies at multiple loci, suggesting that natural selection for low salinity performance at the larval stage can act as a filter for genotypes found in adult populations.

Selection Experiments in the Sea: What Can Experimental Evolution Tell Us About How Marine Life Will Respond to Climate Change?

AbstractRapid evolution may provide a buffer against extinction risk for some species threatened by climate change; however, the capacity to evolve rapidly enough to keep pace with changing environments is unknown for most taxa. The ecosystem-level consequences of climate adaptation are likely to be the largest in marine ecosystems, where short-lived phytoplankton with large effective population sizes make up the bulk of primary production. However, there are substantial challenges to predicting climate-driven evolution in marine systems, including multiple simultaneous axes of change and considerable heterogeneity in rates of change, as well as the biphasic life cycles of many marine metazoans, which expose different life stages to disparate sources of selection. A critical tool for addressing these challenges is experimental evolution, where populations of organisms are directly exposed to controlled sources of selection to test evolutionary responses. We review the use of experimental evolution to test the capacity to adapt to climate change stressors in marine species. The application of experimental evolution in this context has grown dramatically in the past decade, shedding light on the capacity for evolution, associated trade-offs, and the genetic architecture of stress-tolerance traits. Our goal is to highlight the utility of this approach for investigating potential responses to climate change and point a way forward for future studies.

Transgenerational plasticity and the capacity to adapt to low salinity in the eastern oyster, Crassostrea virginica.

Salinity conditions in oyster breeding grounds in the Gulf of Mexico are expected to drastically change due to increased precipitation from climate change and anthropogenic changes to local hydrology. We determined the capacity of the eastern oyster, Crassostrea virginica, to adapt via standing genetic variation or acclimate through transgenerational plasticity (TGP). We outplanted oysters to either a low- or medium-salinity site in Louisiana for 2 years. We then crossed adult parents using a North Carolina II breeding design, and measured body size and survival of larvae 5 dpf raised under low or ambient salinity. We found that TGP is unlikely to significantly contribute to low-salinity tolerance since we did not observe increased growth or survival in offspring reared in low salinity when their parents were also acclimated at a low-salinity site. However, we detected genetic variation for body size, with an estimated heritability of 0.68 ± 0.25 (95% CI). This suggests there is ample genetic variation for this trait to evolve, and that evolutionary adaptation is a possible mechanism through which oysters will persist with future declines in salinity. The results of this experiment provide valuable insights into successfully breeding low-salinity tolerance in this commercially important species.

Transcriptomic signatures of temperature adaptation in the eastern oyster Crassostrea virginica.

The large geographic distribution of the eastern oyster, Crassostrea virginica, makes it an ideal species to test how populations have adapted to latitudinal gradients in temperature. Despite inhabiting distinct thermal regimes, populations of C. virginica near the species' southern and northern geographic range show no population differences in their physiological response to temperature. In this study, we used comparative transcriptomics to understand how oysters from either end of the species' range maintain enantiostasis across three acclimation temperatures (10, 20, and 30°C). With this approach, we identified genes that were differentially expressed in response to temperature between individuals of C. virginica collected from New Brunswick, Canada and Louisiana, USA. We observed a core set of genes whose expression responded to temperature in both populations, but also an even larger set of genes with expression patterns that were unique to each population. Intriguingly, the genes with population-specific responses to temperature had elevated FST and Ka/Ks ratios compared to the genome-wide average. In contrast, genes showing only a response to temperature were found to only have elevated FST values suggesting that divergent FST may be due to selection on linked regulatory regions rather than positive selection on protein coding regions. Taken together, our results suggest that, despite coarse-scale physiological similarities, natural selection has shaped divergent gene expression responses to temperature in geographically separated populations of this broadly eurythermal marine invertebrate.

An Experimental Test of Adaptive Introgression in Locally Adapted Populations of Splash Pool Copepods.

As species struggle to keep pace with the rapidly warming climate, adaptive introgression of beneficial alleles from closely related species or populations provides a possible avenue for rapid adaptation. We investigate the potential for adaptive introgression in the copepod, Tigriopus californicus, by hybridizing two populations with divergent heat tolerance limits. We subjected hybrids to strong heat selection for 15 generations followed by whole-genome resequencing. Utilizing a hybridize evolve and resequence (HER) technique, we can identify loci responding to heat selection via a change in allele frequency. We successfully increased the heat tolerance (measured as LT50) in selected lines, which was coupled with higher frequencies of alleles from the southern (heat tolerant) population. These repeatable changes in allele frequencies occurred on all 12 chromosomes across all independent selected lines, providing evidence that heat tolerance is polygenic. These loci contained genes with lower protein-coding sequence divergence than the genome-wide average, indicating that these loci are highly conserved between the two populations. In addition, these loci were enriched in genes that changed expression patterns between selected and control lines in response to a nonlethal heat shock. Therefore, we hypothesize that the mechanism of heat tolerance divergence is explained by differential gene expression of highly conserved genes. The HER approach offers a unique solution to identifying genetic variants contributing to polygenic traits, especially variants that might be missed through other population genomic approaches.

Population epigenetic divergence exceeds genetic divergence in the Eastern oyster Crassostrea virginica in the Northern Gulf of Mexico.

Populations may respond to environmental heterogeneity via evolutionary divergence or phenotypic plasticity. While evolutionary divergence occurs through DNA sequence differences among populations, plastic divergence among populations may be generated by changes in the epigenome. Here, we present the results of a genome-wide comparison of DNA methylation patterns and genetic structure among four populations of Eastern oyster (Crassostrea virginica) in the northern Gulf of Mexico. We used a combination of restriction site-associated DNA sequencing (RADseq) and reduced representation bisulfite sequencing (RRBS) to explore population structure, gene-wide averages of F ST, and DNA methylation differences between oysters inhabiting four estuaries with unique salinity profiles. This approach identified significant population structure despite a moderately low F ST (0.02) across the freshwater boundary of the Mississippi river, a finding that may reflect recent efforts to restore oyster stock populations. Divergence between populations in CpG methylation was greater than for divergence in F ST, likely reflecting environmental effects on DNA methylation patterns. Assessment of CpG methylation patterns across all populations identified that only 26% of methylated DNA was intergenic; and, only 17% of all differentially methylated regions (DMRs) were within these same regions. DMRs within gene bodies between sites were associated with genes known to be involved in DNA damage repair, ion transport, and reproductive timing. Finally, when assessing the correlation between genomic variation and DNA methylation between these populations, we observed population-specific DNA methylation profiles that were not directly associated with single nucleotide polymorphisms or broader gene-body mean F ST trends. Our results suggest that C. virginica may use DNA methylation to generate environmentally responsive plastic phenotypes and that there is more divergence in methylation than divergence in allele frequencies.

Differential responses to ocean acidification between populations of Balanophyllia elegans corals from high and low upwelling environments.

Ocean acidification (OA), the global decrease in surface water pH from absorption of anthropogenic CO2 , may put many marine taxa at risk. However, populations that experience extreme localized conditions, and are adapted to these conditions predicted in the global ocean in 2,100, may be more tolerant to future OA. By identifying locally adapted populations, researchers can examine the mechanisms used to cope with decreasing pH. One oceanographic process that influences pH is wind-driven upwelling. Here we compare two Californian populations of the coral Balanophyllia elegans from distinct upwelling regimes, and test their physiological and transcriptomic responses to experimental seawater acidification. We measured respiration rates, protein and lipid content, and gene expression in corals from both populations exposed to pH levels of 7.8 and 7.4 for 29 days. Corals from the population that experiences lower pH due to high upwelling maintained the same respiration rate throughout the exposure. In contrast, corals from the low upwelling site had reduced respiration rates, protein content and lipid-class content at low pH exposure, suggesting they have depleted their energy reserves. Using RNA-Seq, we found that corals from the high upwelling site upregulated genes involved in calcium ion binding and ion transport, most likely related to pH homeostasis and calcification. In contrast, corals from the low upwelling site downregulated stress response genes at low pH exposure. Divergent population responses to low pH observed in B. elegans highlight the importance of multi-population studies for predicting a species' response to future OA.

Phenotypic and transcriptomic responses to salinity stress across genetically and geographically divergent Tigriopus californicus populations.

Species inhabiting the North American west coast intertidal must tolerate an extremely variable environment, with large fluctuations in both temperature and salinity. Uncovering the mechanisms for this tolerance is key to understanding species' persistence. We tested for differences in salinity tolerance between populations of Tigriopus californicus copepods from locations in northern (Bodega Reserve) and southern (San Diego) California known to differ in temperature, precipitation and humidity. We also tested for differences between populations in their transcriptomic responses to salinity. Although these two populations have ~20% mtDNA sequence divergence and differ strongly in other phenotypic traits, we observed similarities in their phenotypic and transcriptomic responses to low and high salinity stress. Salinity significantly affected respiration rate (increased under low salinity and reduced under high salinity), but we found no significant effect of population on respiration or a population by salinity interaction. Under high salinity, there was no population difference in knock-down response, but northern copepods had a smaller knock-down under low salinity stress, corroborating previous results for T. californicus. Northern and southern populations had a similar transcriptomic response to salinity based on a principle components analysis, although differential gene expression under high salinity stress was three times lower in the northern population compared to the southern population. Transcripts differentially regulated under salinity stress were enriched for "amino acid transport" and "ion transport" annotation categories, supporting previous work demonstrating that the accumulation of free amino acids is important for osmotic regulation in T. californicus.

Transcriptomics reveal transgenerational effects in purple sea urchin embryos: Adult acclimation to upwelling conditions alters the response of their progeny to differential pCO2 levels.

Understanding the mechanisms with which organisms can respond to a rapidly changing ocean is an important research priority in marine sciences, especially in the light of recent predictions regarding the pace of ocean change in the coming decades. Transgenerational effects, in which the experience of the parental generation can shape the phenotype of their offspring, may serve as such a mechanism. In this study, adult purple sea urchins, Strongylocentrotus purpuratus, were conditioned to regionally and ecologically relevant pCO2 levels and temperatures representative of upwelling (colder temperature and high pCO2 ) and nonupwelling (average temperature and low pCO2 ) conditions typical of coastal upwelling regions in the California Current System. Following 4.5 months of conditioning, adults were spawned and offspring were raised under either high or low pCO2 levels, to examine the role of maternal effects. Using RNA-seq and comparative transcriptomics, our results indicate that differential conditioning of the adults had an effect on the gene expression patterns of the progeny during the gastrula stage of early development. For example, maternal conditioning under upwelling conditions intensified the transcriptomic response of the progeny when they were raised under high versus low pCO2 conditions. Additionally, mothers that experienced upwelling conditions produced larger progeny. The overall findings of this study are complex, but do suggest that transgenerational plasticity in situ could act as an important mechanism by which populations might keep pace with rapid environmental change.

Trait Correlations in the Genomics Era.

Thinking about the evolutionary causes and consequences of trait correlations has been dominated by quantitative genetics theory that is focused on hypothetical loci. Since this theory was initially developed, technology has enabled the identification of specific genetic variants that contribute to trait correlations. Here, we review studies of the genetic basis of trait correlations to ask: What has this new information taught us? We find that causal variants can be pleiotropic and/or linked in different ways, indicating that pleiotropy and linkage are not alternative genetic mechanisms. Further, many trait correlations have a polygenic basis, suggesting that both pleiotropy and linkage likely contribute. We discuss implications of these findings for the evolutionary causes and consequences of trait correlations.

Physiological plasticity and local adaptation to elevated pCO2 in calcareous algae: an ontogenetic and geographic approach.

To project how ocean acidification will impact biological communities in the future, it is critical to understand the potential for local adaptation and the physiological plasticity of marine organisms throughout their entire life cycle, as some stages may be more vulnerable than others. Coralline algae are ecosystem engineers that play significant functional roles in oceans worldwide and are considered vulnerable to ocean acidification. Using different stages of coralline algae, we tested the hypothesis that populations living in environments with higher environmental variability and exposed to higher levels of pCO 2 would be less affected by high pCO 2 than populations from a more stable environment experiencing lower levels of pCO 2. Our results show that spores are less sensitive to elevated pCO 2 than adults. Spore growth and mortality were not affected by pCO 2 level; however, elevated pCO 2 negatively impacted the physiology and growth rates of adults, with stronger effects in populations that experienced both lower levels of pCO 2 and lower variability in carbonate chemistry, suggesting local adaptation. Differences in physiological plasticity and the potential for adaptation could have important implications for the ecological and evolutionary responses of coralline algae to future environmental changes.

Adaptation to climate change: trade-offs among responses to multiple stressors in an intertidal crustacean.

Trade-offs may influence both physiological and evolutionary responses to co-occurring stressors, but their effects on both plastic and adaptive responses to climate change are poorly understood. To test for genetic and physiological trade-offs incurred in tolerating multiple stressors, we hybridized two populations of the intertidal copepod Tigriopus californicus that were divergent for both heat and salinity tolerance. Starting in the F2 generation, we selected for increased tolerance of heat, low salinity, and high salinity in replicate lines. After five generations of selection, heat-selected lines had greater heat tolerance but lower fecundity, indicating an energetic cost to tolerance. Lines selected for increased salinity tolerance did not show evidence of adaptation to their respective environments; however, hypo-osmotic selection lines showed substantial loss of tolerance to hyperosmotic stress. Neither of the salinity selection regimes resulted in diminished heat tolerance at ambient salinity; however, simultaneous exposure to heat and hypo-osmotic stress led to decreased heat tolerance, implying a physiological trade-off in tolerance to the two stressors. When we quantified the transcriptomic response to heat and salinity stress via RNA sequencing, we observed little overlap in the stress responses, suggesting the observed synergistic effects of heat and salinity stress were driven by competing energetic demands, rather than shared stress response pathways.

Plastic and Evolved Responses to Global Change: What Can We Learn from Comparative Transcriptomics?

Physiological plasticity and adaptive evolution may facilitate persistence in a changing environment. As a result, there is an interest in understanding species' capacities for plastic and evolved responses, and the mechanisms by which these responses occur. Transcriptome sequencing has become a powerful tool for addressing these questions, providing insight into otherwise unobserved effects of changing conditions on organismal physiology and variation in these effects among individuals and populations. Here, we review recent studies using comparative transcriptomics to understand plastic and evolutionary responses to changing environments. We focus on 2 areas where transcriptomics has played an important role: first, in understanding the genetic basis for local adaptation to current gradients as a proxy for future adaptation, and second, in understanding organismal responses to multiple stressors. We find most studies examining multiple stressors have tested the effects of each stressor individually; the few studies testing multiple stressors simultaneously have found synergistic effects on gene expression that would not have been predicted from single stressor studies. We discuss the importance of robust experimental design to allow for a more sophisticated characterization of transcriptomic responses and conclude by offering recommendations for future research, including integrating genomics with transcriptomics, testing gene regulatory networks, and comparing the equivalence of transcription to translation and the effects of environmental stress on the proteome.