Methylation and gene expression data from: Differential DNA methylation across environments has no effect on gene expression in the eastern oyster
Data files
Dec 07, 2021 version files 3.30 GB
-
All_nonsig_methylation.matrix
103.24 MB
-
C_virginica-3.0_TE-all.gff
66.06 MB
-
CV_Gene_Annotations_Final.txt
9.50 MB
-
CV_GO_annotations_FINAL.tab
10.78 MB
-
Cvirginica_gene_CpG_OE_ratio.tab
720.93 KB
-
expression_data.matrix
31.08 MB
-
GCF_002022765.2_C_virginica-3.0_genomic.gff
431.80 MB
-
GL01_bismark.cov
8.80 MB
-
GL02_bismark.cov
26.13 MB
-
GL03_bismark.cov
32.39 MB
-
GL05_bismark.cov
39.73 MB
-
GL06_bismark.cov
23.51 MB
-
GL07_bismark.cov
20.47 MB
-
GL08_bismark.cov
56.02 MB
-
GL09_bismark.cov
20.81 MB
-
GL10_bismark.cov
20.63 MB
-
GL11_bismark.cov
30.79 MB
-
GL12_bismark.cov
4.90 MB
-
GL13_bismark.cov
48.44 MB
-
GL14_bismark.cov
21.53 MB
-
GL15_bismark.cov
32.47 MB
-
GL16_bismark.cov
30.65 MB
-
GL17_bismark.cov
45.94 MB
-
GL18_bismark.cov
37.01 MB
-
GL19_bismark.cov
37.32 MB
-
GL21_bismark.cov
21.11 MB
-
GL23_bismark.cov
11.40 MB
-
GV01_bismark.cov
43.52 MB
-
GV02_bismark.cov
37.25 MB
-
GV03_bismark.cov
72.39 MB
-
GV04_bismark.cov
30.68 MB
-
GV05_bismark.cov
11.43 MB
-
GV07_bismark.cov
51.03 MB
-
GV08_bismark.cov
17.92 MB
-
GV09_bismark.cov
63.78 MB
-
GV10_bismark.cov
36.39 MB
-
GV11_bismark.cov
52.27 MB
-
GV12_bismark.cov
26.10 MB
-
GV13_bismark.cov
32.82 MB
-
GV14_bismark.cov
38.08 MB
-
GV15_bismark.cov
168.74 MB
-
GV16_bismark.cov
64.67 MB
-
GV17_bismark.cov
59.80 MB
-
GV18_bismark.cov
45.09 MB
-
GV21_bismark.cov
74.92 MB
-
GV22_bismark.cov
56.53 MB
-
GV24_bismark.cov
71.71 MB
-
GV25_bismark.cov
34.44 MB
-
LL01_bismark.cov
19.90 MB
-
LL02_bismark.cov
15.95 MB
-
LL03_bismark.cov
28.56 MB
-
LL04_bismark.cov
17.39 MB
-
LL05_bismark.cov
25.59 MB
-
LL06_bismark.cov
22.25 MB
-
LL07_bismark.cov
31.33 MB
-
LL08_bismark.cov
24.55 MB
-
LL09_bismark.cov
22.67 MB
-
LL10_bismark.cov
35.77 MB
-
LL11_bismark.cov
24.85 MB
-
LL12_bismark.cov
22.46 MB
-
LL13_bismark.cov
22.26 MB
-
LL14_bismark.cov
3.49 MB
-
LL15_bismark.cov
19.17 MB
-
LL16_bismark.cov
21.92 MB
-
LL17_bismark.cov
30.50 MB
-
LL18_bismark.cov
29.88 MB
-
LL19_bismark.cov
32.49 MB
-
LL20_bismark.cov
14.96 MB
-
LV01_bismark.cov
18.44 MB
-
LV05_bismark.cov
33.34 MB
-
LV06_bismark.cov
30.32 MB
-
LV07_bismark.cov
39.13 MB
-
LV09_bismark.cov
30.01 MB
-
LV10_bismark.cov
28.37 MB
-
LV11_bismark.cov
29.52 MB
-
LV12_bismark.cov
20.24 MB
-
LV13_bismark.cov
38.50 MB
-
LV14_bismark.cov
28.79 MB
-
LV15_bismark.cov
17.84 MB
-
LV16_bismark.cov
20.64 MB
-
LV17_bismark.cov
16.68 MB
-
LV18_bismark.cov
35.56 MB
-
LV19_bismark.cov
27.13 MB
-
LV20_bismark.cov
1.29 MB
-
LV21_bismark.cov
28.37 MB
-
LV22_bismark.cov
22.34 MB
-
LV24_bismark.cov
26.41 MB
-
LV25_bismark.cov
31.58 MB
-
Mapping_Stats.xlsx
14.67 KB
-
README.txt
725 B
Abstract
1. It has been hypothesized that environmentally induced changes to gene body methylation could facilitate adaptive transgenerational responses to changing environments.
2. 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.
3. 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.
4. 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.
5. 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.
6. 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.
In May 2016, adult oysters (C. virginica) were collected by dredging from Sister Lake, LA (29˚14’57” N, 90˚55’16” W). These oysters were transported to the Louisiana Department of Wildlife and Fisheries Michael C. Voisin Oyster Hatchery in Grand Isle, LA (29°14'20.3" N, 90°00'11.2" W) and placed into off-bottom mesh cages for long-term acclimation. In October 2016, after six months of acclimation, the oysters were spawned at the MCV oyster hatchery using 3 males and 2 females. Oyster spat were reared in an upwelling system, individually tagged, and outplanted in one of three adjustable long-line mesh bags at both the Grand Isle Hatchery farm and near the Louisiana Universities Marine Consortium (LUMCON) (29°15'12.6" N, 90°39'45.9" W) on February 20th, 2017. Because larvae from all six families were combined for culturing, parentage assignments were unknown at the time of outplant, and as such families were unequally outplanted between sites. Oysters within each bag were monitored for mortality and cleaned of epibionts approximately every 3 months over a 14-month period. On April 24th, 2018, after 14-months at the two outplant sites, 40 individuals were haphazardly chosen from each site. Shell height of each individual was measured from shell umbo to distal edge using a digital caliper (ABS Coolant Proof Calipers, Mituyoto Corporation, Japan). Approximately 1 cm2 piece of gill tissue was sampled in the field from each individual and preserved with either Invitrogen RNAlater (gene expression) or 95% ethanol (DNA methylation). The remaining whole animal was placed in a pre-weighed 50 ml test tube and used to measure wet meat weight and Perkinsus marinus infection intensities.
Gene expression
Total RNA was extracted using a E.Z.N.A.® Total RNA Kit I (Omega BIO-TEK Inc., Norcross, GA, USA) following the manufacturer's instructions. The yield and quantity were initially assessed using a NanoDrop 2000 spectrophotometer. Total RNA extracted from the 80 individuals was sent to the University of Texas at Austin’s Genomic Sequencing and Analysis Facility where RNA quality control was confirmed using a 2100 Agilent Bioanalyzer on a Eukaryote Total RNA Nano chip and libraries were produced using the Tag-Sequencing approach (Meyer, Aglyamova, & Matz, 2011). The resulting 80 libraries were sequenced on two lanes of an Illumina HiSeq 2500 platform, with 100 base pair single-end reads.
Sequencing reads were trimmed of adapter sequences using Trimmomatic (version 0.39) (Bolger, Lohse, & Usadel, 2014) and base pairs with quality scores below 30 were removed (Table S2). The trimmed reads were mapped to the C. virginica reference genome (Gómez-Chiarri, Warren, Guo, & Proestou, 2015) with known haplotigs removed (https://github.com/jpuritz/OysterGenomeProject/tree/master/Haplotig_Masked_Genome) using the single pass option for STAR RNA-seq aligner (version 2.6.0a) (Dobin et al., 2013). Reads were mapped to gene features with the options (--quantMode GeneCounts --outFilterScoreMinOverLread 0.50 --outFilterMatchNminOverLread 0.50) specified to adjust for poly-A tail contamination. A count matrix was generated from the ReadsPerGene.out.tab output from STAR.
The R scripts files will provide data analysis steps needed for reproducing the results.