Epigenetic-associated phenotypic plasticity of the ocean acidification-acclimated edible oyster in the mariculture environment.

For marine invertebrates with a pelagic-benthic life cycle, larval exposure to ocean acidification (OA) can affect adult performance in response to another environmental stressor. This carry-over effect has the potential to alter phenotypic traits. However, the molecular mechanisms that mediate "OA"-triggered carry-over effects have not been explored despite such information being key to improving species fitness and management strategies for aquafarming. This study integrated the genome-wide DNA methylome and transcriptome to examine epigenetic modification-mediated carry-over OA impacts on phenotypic traits of the ecologically and commercially important oyster species Crassostrea hongkongensis under field conditions. Larvae of C. hongkongensis were exposed to control pH 8.0 and low pH 7.4 conditions, mimicking near future OA scenario in their habitat, before being outplanted as post-metamorphic juveniles at two mariculture field sites with contrasting environmental stressors for 9 months. The larval carry-over OA effect was found to have persistent impacts on the growth and survival trade-off traits on the outplanted juveniles, although the beneficial or adverse effect depended on the environmental conditions at the outplanted sites. Site-specific plasticity was demonstrated with a diverse DNA methylation-associated gene expression profile, with signal transduction and the endocrine system being the most common and highly enriched functions. Highly methylated exons prevailed in the key genes related to general metabolic and endocytic responses and these genes are evolutionarily conserved in various marine invertebrates in response to OA. These results suggest that oysters with prior larval exposure history to OA had the ability to trigger rapid local adaptive responses via epigenetic modification to cope with multiple stressors in the field.

Triploid Pacific oysters exhibit stress response dysregulation and elevated mortality following heatwaves.

Polyploidy has been suggested to negatively impact environmental stress tolerance, resulting in increased susceptibility to extreme climate events. In this study, we compared the genomic and physiological response of diploid (2n) and triploid (3n) Pacific oysters (Crassostrea gigas) to conditions present during an atmospheric heatwave that impacted the Pacific Northwestern region of the United States in the summer of 2021. Climate stressors were applied either singly (single stressor; elevated seawater temperature, 30°C) or in succession (multiple stressor; elevated seawater temperature followed by aerial emersion at 44°C), replicating conditions present within the intertidal over a tidal cycle during the event. Oyster mortality rate was elevated within stress treatments with respect to the control and was significantly higher in triploids than diploids following multiple stress exposure (36.4% vs. 14.8%). Triploids within the multiple stressor treatment exhibited signs of energetic limitation, including metabolic depression, a significant reduction in ctenidium Na+ /K+ ATPase activity, and the dysregulated expression of genes associated with stress response, innate immunity, glucose metabolism, and mitochondrial function. Functional enrichment analysis of ploidy-specific gene sets identified that biological processes associated with metabolism, stress tolerance, and immune function were overrepresented within triploids across stress treatments. Our results suggest that triploidy impacts the transcriptional regulation of key processes that underly the stress response of Pacific oysters, resulting in downstream shifts in physiological tolerance limits that may increase susceptibility to extreme climate events that present multiple environmental stressors. The impact of chromosome set manipulation on the climate resilience of marine organisms has important implications for domestic food security within future climate scenarios, especially as triploidy induction becomes an increasingly popular tool to elicit reproductive control across a wide range of species used within marine aquaculture.

Differential DNA methylation in Pacific oyster reproductive tissue in response to ocean acidification.

BACKGROUND: There is a need to investigate mechanisms of phenotypic plasticity in marine invertebrates as negative effects of climate change, like ocean acidification, are experienced by coastal ecosystems. Environmentally-induced changes to the methylome may regulate gene expression, but methylome responses can be species- and tissue-specific. Tissue-specificity has implications for gonad tissue, as gonad-specific methylation patterns may be inherited by offspring. We used the Pacific oyster (Crassostrea gigas) - a model for understanding pH impacts on bivalve molecular physiology due to its genomic resources and importance in global aquaculture- to assess how low pH could impact the gonad methylome. Oysters were exposed to either low pH (7.31 ± 0.02) or ambient pH (7.82 ± 0.02) conditions for 7 weeks. Whole genome bisulfite sequencing was used to identify methylated regions in female oyster gonad samples. C- > T single nucleotide polymorphisms were identified and removed to ensure accurate methylation characterization. RESULTS: Analysis of gonad methylomes revealed a total of 1284 differentially methylated loci (DML) found primarily in genes, with several genes containing multiple DML. Gene ontologies for genes containing DML were involved in development and stress response, suggesting methylation may promote gonad growth homeostasis in low pH conditions. Additionally, several of these genes were associated with cytoskeletal structure regulation, metabolism, and protein ubiquitination - commonly-observed responses to ocean acidification. Comparison of these DML with other Crassostrea spp. exposed to ocean acidification demonstrates that similar pathways, but not identical genes, are impacted by methylation. CONCLUSIONS: Our work suggests DNA methylation may have a regulatory role in gonad and larval development, which would shape adult and offspring responses to low pH stress. Combined with existing molluscan methylome research, our work further supports the need for tissue- and species-specific studies to understand the potential regulatory role of DNA methylation.

Acclimatory gene expression of primed clams enhances robustness to elevated pCO2.

Sublethal exposure to environmental challenges may enhance ability to cope with chronic or repeated change, a process known as priming. In a previous study, pre-exposure to seawater enriched with pCO2 improved growth and reduced antioxidant capacity of juvenile Pacific geoduck Panopea generosa clams, suggesting that transcriptional shifts may drive phenotypic modifications post-priming. To this end, juvenile clams were sampled and TagSeq gene expression data were analysed after (i) a 110-day acclimation under ambient (921 µatm, naïve) and moderately elevated pCO2 (2870 µatm, pre-exposed); then following (ii) a second 7-day exposure to three pCO2 treatments (ambient: 754 µatm; moderately elevated: 2750 µatm; severely elevated: 4940 µatm), a 7-day return to ambient pCO2 and a third 7-day exposure to two pCO2 treatments (ambient: 967 µatm; moderately elevated: 3030 µatm). Pre-exposed geoducks frontloaded genes for stress and apoptosis/innate immune response, homeostatic processes, protein degradation and transcriptional modifiers. Pre-exposed geoducks were also responsive to subsequent encounters, with gene sets enriched for mitochondrial recycling and immune defence under elevated pCO2 and energy metabolism and biosynthesis under ambient recovery. In contrast, gene sets with higher expression in naïve clams were enriched for fatty-acid degradation and glutathione components, suggesting naïve clams could be depleting endogenous fuels, with unsustainable energetic requirements if changes in carbonate chemistry persist. Collectively, our transcriptomic data indicate that pCO2 priming during post-larval periods could, via gene expression regulation, enhance robustness in bivalves to environmental change. Such priming approaches may be beneficial for aquaculture, as seafood demand intensifies concurrent with increasing climate change in marine systems.

Oyster biomineralization under ocean acidification: From genes to shell.

Biomineralization is one of the key processes that is notably affected in marine calcifiers such as oysters under ocean acidification (OA). Understanding molecular changes in the biomineralization process under OA and its heritability, therefore, is key to developing conservation strategies for protecting ecologically and economically important oyster species. To do this, in this study, we have explicitly chosen the tissue involved in biomineralization (mantle) of an estuarine commercial oyster species, Crassostrea hongkongensis. The primary aim of this study is to understand the influence of DNA methylation over gene expression of mantle tissue under decreased ~pH 7.4, a proxy of OA, and to extrapolate if these molecular changes can be observed in the product of biomineralization-the shell. We grew early juvenile C. hongkongensis, under decreased ~pH 7.4 and control ~pH 8.0 over 4.5 months and studied OA-induced DNA methylation and gene expression patterns along with shell properties such as microstructure, crystal orientation and hardness. The population of oysters used in this study was found to be moderately resilient to OA at the end of the experiment. The expression of key biomineralization-related genes such as carbonic anhydrase and alkaline phosphatase remained unaffected; thus, the mechanical properties of the shell (shell growth rate, hardness and crystal orientation) were also maintained without any significant difference between control and OA conditions with signs of severe dissolution. In addition, this study makes three major conclusions: (1) higher expression of Ca2+ binding/signalling-related genes in the mantle plays a key role in maintaining biomineralization under OA; (2) DNA methylation changes occur in response to OA; however, these methylation changes do not directly control gene expression; and (3) OA would be more of a 'dissolution problem' rather than a 'biomineralization problem' for resilient species that maintain calcification rate with normal shell growth and mechanical properties.

Repeat exposure to hypercapnic seawater modifies growth and oxidative status in a tolerant burrowing clam.

Although low levels of thermal stress, irradiance and dietary restriction can have beneficial effects for many taxa, stress acclimation remains little studied in marine invertebrates, even though they are threatened by climate change stressors such as ocean acidification. To test the role of life-stage and stress-intensity dependence in eliciting enhanced tolerance under subsequent stress encounters, we initially conditioned pediveliger Pacific geoduck (Panopea generosa) larvae to ambient and moderately elevated PCO2 (920 µatm and 2800 µatm, respectively) for 110 days. Then, clams were exposed to ambient, moderate or severely elevated PCO2 (750, 2800 or 4900 µatm, respectively) for 7 days and, following 7 days in ambient conditions, a 7-day third exposure to ambient (970 µatm) or moderate PCO2 (3000 µatm). Initial conditioning to moderate PCO2 stress followed by second and third exposure to severe and moderate PCO2 stress increased respiration rate, organic biomass and shell size, suggesting a stress-intensity-dependent effect on energetics. Additionally, stress-acclimated clams had lower antioxidant capacity compared with clams under ambient conditions, supporting the hypothesis that stress over postlarval-to-juvenile development affects oxidative status later in life. Time series and stress intensity-specific approaches can reveal life-stages and magnitudes of exposure, respectively, that may elicit beneficial phenotypic variation.

Temporal proteomic profiling reveals insight into critical developmental processes and temperature-influenced physiological response differences in a bivalve mollusc.

BACKGROUND: Protein expression patterns underlie physiological processes and phenotypic differences including those occurring during early development. The Pacific oyster (Crassostrea gigas) undergoes a major phenotypic change in early development from free-swimming larval form to sessile benthic dweller while proliferating in environments with broad temperature ranges. Despite the economic and ecological importance of the species, physiological processes occurring throughout metamorphosis and the impact of temperature on these processes have not yet been mapped out. RESULTS: Towards this, we comprehensively characterized protein abundance patterns for 7978 proteins throughout metamorphosis in the Pacific oyster at different temperature regimes. We used a multi-statistical approach including principal component analysis, ANOVA-simultaneous component analysis, and hierarchical clustering coupled with functional enrichment analysis to characterize these data. We identified distinct sets of proteins with time-dependent abundances generally not affected by temperature. Over 12 days, adhesion and calcification related proteins acutely decreased, organogenesis and extracellular matrix related proteins gradually decreased, proteins related to signaling showed sinusoidal abundance patterns, and proteins related to metabolic and growth processes gradually increased. Contrastingly, different sets of proteins showed temperature-dependent abundance patterns with proteins related to immune response showing lower abundance and catabolic pro-growth processes showing higher abundance in animals reared at 29 °C relative to 23 °C. CONCLUSION: Although time was a stronger driver than temperature of metamorphic proteome changes, temperature-induced proteome differences led to pro-growth physiology corresponding to larger oyster size at 29 °C, and to altered specific metamorphic processes and possible pathogen presence at 23 °C. These findings offer high resolution insight into why oysters may experience high mortality rates during this life transition in both field and culture settings. The proteome resource generated by this study provides data-driven guidance for future work on developmental changes in molluscs. Furthermore, the analytical approach taken here provides a foundation for effective shotgun proteomic analyses across a variety of taxa.

Carryover effects of temperature and pCO2 across multiple Olympia oyster populations.

Predicting how populations will respond to ocean change across generations is critical to effective conservation of marine species. One emerging factor is the influence of parental exposures on offspring phenotype, known as intergenerational carryover effects. Parental exposure may deliver beneficial or detrimental characteristics to offspring that can influence larval recruitment patterns, thus shaping how populations and community structure respond to ocean change. Impacts of adult exposure to elevated winter temperature and pCO2 on reproduction and offspring viability were examined in the Olympia oyster (Ostrea lurida) using three populations of adult, hatchery-reared O. lurida, plus an additional cohort spawned from one of the populations. Oysters were sequentially exposed to elevated temperature (+4°C, at 10°C), followed by elevated pCO2 (+2,204 µatm, at 3,045 µatm) during winter months. Male gametes were more developed after elevated temperature exposure and less developed after high pCO2 exposure, but there was no impact on female gametes or sex ratios. Oysters previously exposed to elevated winter temperature released larvae earlier, regardless of pCO2 exposure. Those exposed to elevated winter temperature as a sole treatment released more larvae on a daily basis but, when also exposed to high pCO2 , there was no effect. These combined results indicate that elevated winter temperature accelerates O. lurida spermatogenesis, resulting in earlier larval release and increased production, with elevated pCO2 exposure negating effects of elevated temperature. Altered recruitment patterns may therefore follow warmer winters due to precocious spawning, but these effects may be masked by coincidental high pCO2 . Offspring were reared in common conditions for 1 yr, then deployed for 3 months in four estuarine bays with distinct environmental conditions. Offspring of parents exposed to elevated pCO2 had higher survival rates in two of the four bays. This carryover effect demonstrates that parental conditions can have substantial ecologically relevant impacts that should be considered when predicting impacts of environmental change. Furthermore, Olympia oysters may be more resilient in certain environments when progenitors are pre-conditioned in stressful conditions. Combined with other recent studies, our work suggests that the Olympia may be more equipped than other oysters for the challenge of a changing ocean.

Integrating Discovery-driven Proteomics and Selected Reaction Monitoring To Develop a Noninvasive Assay for Geoduck Reproductive Maturation.

Geoduck clams (Panopea generosa) are an increasingly important fishery and aquaculture product along the eastern Pacific coast from Baja California, Mexico, to Alaska. These long-lived clams are highly fecund, although sustainable hatchery production of genetically diverse larvae is hindered by the lack of sexual dimorphism, resulting in asynchronous spawning of broodstock, unequal sex ratios, and low numbers of breeders. The development of assays of gonad physiology could indicate sex and maturation stage as well as be used to assess the status of natural populations. Proteomic profiles were determined for three reproductive maturation stages in both male and female clams using data-dependent acquisition (DDA) of gonad proteins. Gonad proteomes became increasingly divergent between males and females as maturation progressed. The DDA data were used to develop targets analyzed with selected reaction monitoring (SRM) in gonad tissue as well as hemolymph. The SRM assay yielded a suite of indicator peptides that can be used as an efficient assay to determine geoduck gonad maturation status. Application of SRM in hemolymph samples demonstrates that this procedure could effectively be used to assess reproductive status in marine mollusks in a nonlethal manner.

Genetic and epigenetic insight into morphospecies in a reef coral.

Incongruence between conventional and molecular systematics has left the delineation of many species unresolved. Reef-building corals are no exception, with phenotypic plasticity among the most plausible explanations for alternative morphospecies. As potential molecular signatures of phenotypic plasticity, epigenetic processes may contribute to our understanding of morphospecies. We compared genetic and epigenetic variation in Caribbean branching Porites spp., testing the hypothesis that epigenetics-specifically, differential patterns of DNA methylation-play a role in alternative morphotypes of a group whose taxonomic status has been questioned. We used reduced representation genome sequencing to analyse over 1,000 single nucleotide polymorphisms and CpG sites in 27 samples of Porites spp. exhibiting a range of morphotypes from a variety of habitats in Belize. We found stronger evidence for genetic rather than epigenetic structuring, identifying three well-defined genetic groups. One of these groups exhibited significantly thicker branches, and branch thickness was a better predictor of genetic groups than depth, habitat or symbiont type. In contrast, no clear epigenetic patterns emerged with respect to phenotypic or habitat variables. While there was a weak positive correlation between pairwise genetic and epigenetic distance, two pairs of putative clones exhibited substantial epigenetic differences, suggesting a strong environmental effect. We speculate that epigenetic patterns are a complex mosaic reflecting diverse environmental histories superimposed over a relatively small heritable component. Given the role of genetics in branching Porites spp. morphospecies we were able to detect with genomewide sequencing, use of such techniques throughout the geographic range of these corals may help settle their phylogeny.

Germline DNA methylation in reef corals: patterns and potential roles in response to environmental change.

DNA methylation is an epigenetic mark that plays an inadequately understood role in gene regulation, particularly in nonmodel species. Because it can be influenced by the environment, DNA methylation may contribute to the ability of organisms to acclimatize and adapt to environmental change. We evaluated the distribution of gene body methylation in reef-building corals, a group of organisms facing significant environmental threats. Gene body methylation in six species of corals was inferred from in silico transcriptome analysis of CpG O/E, an estimate of germline DNA methylation that is highly correlated with patterns of methylation enrichment. Consistent with what has been documented in most other invertebrates, all corals exhibited bimodal distributions of germline methylation suggestive of distinct fractions of genes with high and low levels of methylation. The hypermethylated fractions were enriched with genes with housekeeping functions, while genes with inducible functions were highly represented in the hypomethylated fractions. High transcript abundance was associated with intermediate levels of methylation. In three of the coral species, we found that genes differentially expressed in response to thermal stress and ocean acidification exhibited significantly lower levels of methylation. These results support a link between gene body hypomethylation and transcriptional plasticity that may point to a role of DNA methylation in the response of corals to environmental change.

Evaluating hypoxia-inducible factor-1α mRNA expression in a pelagic fish, Pacific herring Clupea pallasii, as a biomarker for hypoxia exposure.

Hypoxia [dissolved oxygen (DO)<2 mg L(-1)] is a major environmental perturbation for many aquatic ecosystems, particularly highly productive estuaries. Most research attention and understanding about the impacts of hypoxia on estuarine species has focused on the benthos, where hypoxia is most common. Although the pelagic zone is also susceptible to the effects of hypoxia, the biological interactions and consequences are not as well understood in marine environments because documenting exposure or avoidance of hypoxia is often difficult. Physiological biomarkers may provide a way to gain more detailed spatiotemporal information regarding species' exposure to hypoxia. Here, we identified and tested a hypoxia-specific responsive gene, hypoxia-inducible factor-1α (hif-1α), to evaluate its potential as a biomarker for hypoxia exposure in Pacific herring (Clupea pallasii). We conducted controlled laboratory experiments to establish the level of hepatic hif-1α elevated gene expression (>1 sd normoxic mean), exposure amplification (≥2 hours), reduction rate (ca. 24 hours), and some evidence of a lethal hypoxic limit (ca. 2 mg L(-1), ≥4 hours). We then used these findings to evaluate the spatiotemporal patterns of hif-1α for Pacific herring in a seasonally hypoxia estuary, Hood Canal, Washington, USA. Although expression did not parallel the local hypoxic conditions in the estuary, herring from the more severe hypoxic year (2013) had a higher probability of having elevated mRNA levels. These patterns indicate that hepatic hif-1α levels may not be directly indicative of local DO levels for pelagic marine fish, but rather provide insight into hypoxia exposure over broader scales.

Shotgun proteomics reveals physiological response to ocean acidification in Crassostrea gigas.

BACKGROUND: Ocean acidification as a result of increased anthropogenic CO2 emissions is occurring in marine and estuarine environments worldwide. The coastal ocean experiences additional daily and seasonal fluctuations in pH that can be lower than projected end-of-century open ocean pH reductions. In order to assess the impact of ocean acidification on marine invertebrates, Pacific oysters (Crassostrea gigas) were exposed to one of four different p CO2 levels for four weeks: 400 µatm (pH 8.0), 800 µatm (pH 7.7), 1000 µatm (pH 7.6), or 2800 µatm (pH 7.3). RESULTS: At the end of the four week exposure period, oysters in all four p CO2 environments deposited new shell, but growth rate was not different among the treatments. However, micromechanical properties of the new shell were compromised by elevated p CO2. Elevated p CO2 affected neither whole body fatty acid composition, nor glycogen content, nor mortality rate associated with acute heat shock. Shotgun proteomics revealed that several physiological pathways were significantly affected by ocean acidification, including antioxidant response, carbohydrate metabolism, and transcription and translation. Additionally, the proteomic response to a second stress differed with p CO2, with numerous processes significantly affected by mechanical stimulation at high versus low p CO2 (all proteomics data are available in the ProteomeXchange under the identifier PXD000835). CONCLUSIONS: Oyster physiology is significantly altered by exposure to elevated p CO2, indicating changes in energy resource use. This is especially apparent in the assessment of the effects of p CO2 on the proteomic response to a second stress. The altered stress response illustrates that ocean acidification may impact how oysters respond to other changes in their environment. These data contribute to an integrative view of the effects of ocean acidification on oysters as well as physiological trade-offs during environmental stress.

Epigenetic and Genetic Population Structure is Coupled in a Marine Invertebrate.

Delineating the relative influence of genotype and the environment on DNA methylation is critical for characterizing the spectrum of organism fitness as driven by adaptation and phenotypic plasticity. In this study, we integrated genomic and DNA methylation data for two distinct Olympia oyster (Ostrea lurida) populations while controlling for within-generation environmental influences. In addition to providing the first characterization of genome-wide DNA methylation patterns in the oyster genus Ostrea, we identified 3,963 differentially methylated loci between populations. Our results show a clear coupling between genetic and epigenetic patterns of variation, with 27% of variation in interindividual methylation differences explained by genotype. Underlying this association are both direct genetic changes in CpGs (CpG-SNPs) and genetic variation with indirect influence on methylation (mQTLs). When comparing measures of genetic and epigenetic population divergence at specific genomic regions this relationship surprisingly breaks down, which has implications for the methods commonly used to study epigenetic and genetic coupling in marine invertebrates.

Transcriptome profile in heat resilient Pacific oyster Crassostrea gigas families under thermal challenge.

Since the introduction of the Pacific oyster Crassostrea gigas in Baja California Sur, Mexico, its culture has faced environmental challenges, specifically increasing temperatures that result in high mortalities. The inter-tidal zone seawater temperature during a year at the Baja California Peninsula broadly ranges from 7 °C to 39 °C. Therefore, to understand how oysters respond to heat stress during daily temperature oscillations, heat-resistant (RR, father, and mother resistant) and heat-susceptible (SS, both parents susceptible) phenotypes families from a C. gigas breeding program were exposed to a thermal challenge. Based on a laboratory-simulated daily oscillatory thermal challenge (26 to 34 °C) for 30 days, RR phenotype presented differences compared to SS phenotype since the beginning (day 0) of the thermal challenge. Gene expression analyses revealed 1822 differentially expressed up-regulated transcripts in RR, related to functions of metabolic processes, biological regulation, and response to stimulus and signaling. At the end of the experiment (day 30), 2660 differentially expressed up-regulated transcripts were identified in RR. Functional analysis of the genes expressed indicates responses of regulation of biological processes and response to a stimulus. Additionally, 340 genes were differentially expressed among RR vs. SS from the beginning to the end of the thermal challenge, where 170 genes were up-regulated, and 170 were down-regulated. These transcriptomic profiles represent the first report to identify gene expression markers associated with RR phenotypes for the Pacific oyster to the future broodstock selection.

Invertebrate methylomes provide insight into mechanisms of environmental tolerance and reveal methodological biases.

There is a growing focus on the role of DNA methylation in the ability of marine invertebrates to rapidly respond to changing environmental factors and anthropogenic impacts. However, genome-wide DNA methylation studies in nonmodel organisms are currently hampered by a limited understanding of methodological biases. Here, we compare three methods for quantifying DNA methylation at single base-pair resolution-whole genome bisulfite sequencing (WGBS), reduced representation bisulfite sequencing (RRBS), and methyl-CpG binding domain bisulfite sequencing (MBDBS)-using multiple individuals from two reef-building coral species with contrasting environmental sensitivity. All methods reveal substantially greater methylation in Montipora capitata (11.4%) than the more sensitive Pocillopora acuta (2.9%). The majority of CpG methylation in both species occurs in gene bodies and flanking regions. In both species, MBDBS has the greatest capacity for detecting CpGs in coding regions at our sequencing depth, but MBDBS may be influenced by intrasample methylation heterogeneity. RRBS yields robust information for specific loci albeit without enrichment of any particular genome feature and with significantly reduced genome coverage. Relative genome size strongly influences the number and location of CpGs detected by each method when sequencing depth is limited, illuminating nuances in cross-species comparisons. As genome-wide methylation differences, supported by data across bisulfite sequencing methods, may contribute to environmental sensitivity phenotypes in critical marine invertebrate taxa, these data provide a genomic resource for investigating the functional role of DNA methylation in environmental tolerance.

DNA methylation profiling of a cnidarian-algal symbiosis using nanopore sequencing.

Symbiosis with protists is common among cnidarians such as corals and sea anemones and is associated with homeostatic and phenotypic changes in the host that could have epigenetic underpinnings, such as methylation of CpG dinucleotides. We leveraged the sensitivity to base modifications of nanopore sequencing to probe the effect of symbiosis with the chlorophyte Elliptochloris marina on methylation in the sea anemone Anthopleura elegantissima. We first validated the approach by comparison of nanopore-derived methylation levels with CpG depletion analysis of a published transcriptome, finding that high methylation levels are associated with CpG depletion as expected. Next, using reads generated exclusively from aposymbiotic anemones, a largely complete draft genome comprising 243 Mb was assembled. Reads from aposymbiotic and symbiotic sea anemones were then mapped to this genome and assessed for methylation using the program Nanopolish, which detects signal disruptions from base modifications as they pass through the nanopore. Based on assessment of 452,841 CpGs for which there was adequate read coverage (approximately 8% of the CpGs in the genome), symbiosis with E. marina was, surprisingly, associated with only subtle changes in the host methylome. However, we did identify one extended genomic region with consistently higher methylation among symbiotic individuals. The region was associated with a DNA polymerase zeta that is noted for its role in translesion synthesis, which opens interesting questions about the biology of this symbiosis. Our study highlights the power and relative simplicity of nanopore sequencing for studies of nucleic acid base modifications in non-model species.

DNA methylation changes in response to ocean acidification at the time of larval metamorphosis in the edible oyster, Crassostrea hongkongensis.

Unprecedented rate of increased CO2 level in the ocean and the subsequent changes in carbonate system including decreased pH, known as ocean acidification (OA), is predicted to disrupt not only the calcification process but also several other physiological and developmental processes in a variety of marine organisms, including edible oysters. Nonetheless, not all species are vulnerable to those OA threats, e.g. some species may be able to cope with OA stress using environmentally induced modifications on gene and protein expressions. For example, external environmental stressors including OA can influence the addition and removal of methyl groups through epigenetic modification (e.g. DNA methylation) process to turn gene expression "on or off" as part of a rapid adaptive mechanism to cope with OA. In this study, we tested the above hypothesis through testing the effect of OA, using decreased pH 7.4 as proxy, on DNA methylation pattern of an endemic and a commercially important estuary oyster species, Crassostrea hongkongensis at the time of larval habitat selection and metamorphosis. Larval growth rate did not differ between control pH 8.1 and treatment pH 7.4. The metamorphosis rate of the pediveliger larvae was higher at pH 7.4 than those in control pH 8.1, however over one-third of the larvae raised at pH 7.4 failed to attach on optimal substrate as defined by biofilm presence. During larval development, a total of 130 genes were differentially methylated across the two treatments. The differential methylation in the larval genes may have partially accounted for the higher metamorphosis success rate under decreased pH 7.4 but with poor substratum selection ability. Differentially methylated loci were concentrated in the exon regions and appear to be associated with cytoskeletal and signal transduction, oxidative stress, metabolic processes, and larval metamorphosis, which implies the high potential of C. hongkongensis larvae to acclimate and adapt through non-genetic ways to OA threats within a single generation.

DNA methylation changes in response to ocean acidification at the time of larval metamorphosis in the edible oyster, Crassostrea hongkongensis.

Unprecedented rate of increased CO2 level in the ocean and the subsequent changes in carbonate system including decreased pH, known as ocean acidification (OA), is predicted to disrupt not only the calcification process but also several other physiological and developmental processes in a variety of marine organisms, including edible oysters. Nonetheless, not all species are vulnerable to those OA threats, e.g. some species may be able to cope with OA stress using environmentally induced modifications on gene and protein expressions. For example, external environmental stressors including OA can influence the addition and removal of methyl groups through epigenetic modification (e.g. DNA methylation) process to turn gene expression "on or off" as part of a rapid adaptive mechanism to cope with OA. In this study, we tested the above hypothesis through testing the effect of OA, using decreased pH 7.4 as proxy, on DNA methylation pattern of an endemic and a commercially important estuary oyster species, Crassostrea hongkongensis at the time of larval habitat selection and metamorphosis. Larval growth rate did not differ between control pH 8.1 and treatment pH 7.4. The metamorphosis rate of the pediveliger larvae was higher at pH 7.4 than those in control pH 8.1, however over one-third of the larvae raised at pH 7.4 failed to attach on optimal substrate as defined by biofilm presence. During larval development, a total of 130 genes were differentially methylated across the two treatments. The differential methylation in the larval genes may have partially accounted for the higher metamorphosis success rate under decreased pH 7.4 but with poor substratum selection ability. Differentially methylated loci were concentrated in the exon regions and appear to be associated with cytoskeletal and signal transduction, oxidative stress, metabolic processes, and larval metamorphosis, which implies the high potential of C. hongkongensis larvae to acclimate and adapt through non-genetic ways to OA threats within a single generation.