Dataset: Gene expression of Eurytemora affinis (saline and freshwater)

Name: Evolutionary and acclimatory shifts in gene expression of Eurytemora affinis copepods reared in saline and freshwater conditions during laboratory experiments from 2011-2014
Published: 2022-11-10
Keywords: gene expression, Invasive species, osmoregulation, evolutionary adaptation, acclimation, biota, oceans

Cite This Dataset

DOI:10.26008/1912/bco-dmo.883426.1

Spatial Extent: N:48.003889 E:-69.425278 S:42.729444 W:-87.778889

Temporal Extent: 2011-09-28 - 2014-03-18

Project:

Evolutionary Responses to Global Changes in Salinity and Temperature (Evolutionary genomics of a copepod)

Principal Investigator:

Carol E. Lee (University of Wisconsin, UW-Madison)

Student:

Marijan Posavi (University of Wisconsin, UW-Madison)

BCO-DMO Data Manager:

Dana Stuart Gerlach (Woods Hole Oceanographic Institution, WHOI BCO-DMO)

Version:

Version Date:

2022-11-10

Restricted:

Validated:

Yes

Current State:

Final no updates expected

Evolutionary and acclimatory shifts in gene expression of Eurytemora affinis copepods reared in saline and freshwater conditions during laboratory experiments from 2011-2014

Abstract:

To explore mechanisms of freshwater adaptation and distinguish between adaptive (evolutionary) and acclimatory (plastic) responses to salinity change, we examined genome‐wide patterns of gene expression between ancestral saline and derived freshwater populations of the Eurytemora affinis species complex, reared under two different common‐garden conditions (0 vs. 15 PSU). These data include the RNA-seq Illumina short paired end reads (101 base pairs) of 10 freshwater and 12 saline copepods Eurytemora affinis samples. The freshwater copepod samples were collected in Lake Michigan (Racine Harbor), while the saline copepods were collected in Baie de L'Isle Verte, St. Lawrence marsh, Quebec, Canada. These data have important implications for our understanding of the evolutionary and physiological mechanisms of range expansions by some of the most widespread invaders in aquatic habitats.

Expand/Collapse All

Archival Copy

Version Date	Archive	Persistent Identifier (PID)	Date PID Issued
2022-11-10	Marine Biological Laboratory/Woods Hole Oceanographic Institution Library (MBLWHOI DLA)	10.26008/1912/bco-dmo.883426.1	2022-11-12

Methods & Sampling

Study objectives
The goal of this study was to explore evolutionary shifts in gene expression between ancestral saline and freshwater invading populations of the Eurytemora affinis (copepod) species complex on a genome-wide scale.

To explore mechanisms of freshwater adaptation and distinguish between adaptive (evolutionary) and acclimatory (plastic) responses to salinity change, laboratory experiments were conducted using both ancestral saline and derived freshwater populations of Eurytemora affinis. Then RNA-seq data -- Illumina short paired-end (PE) reads (101 base pairs) of 10 freshwater and 12 saline E.affinis samples -- were used to answer the following questions:

(1) What are the patterns of evolutionary shifts in gene expression between the ancestral saline and freshwater invading populations?
(2) What are the plastic (acclimatory) changes in gene expression between salinities (0 PSU vs. 15 PSU conditions) within each of the saline and freshwater populations?
(3) Are the magnitude and direction of plasticity in gene expression correlated with evolutionary responses?
(4) Has plasticity in gene expression evolved following freshwater invasions?

Collection of ancestral populations
The copepods were collected using a plankton net mesh size of 50 μm in diameter, from a depth of 1-4 meters from near the shore. The freshwater copepods were collected in April-May 2006 by throwing the plankton net off the dock in Racine Harbor, Lake Michigan in Wisconsin, USA (42.729444 N, 87.778889 W). The saline copepods were collected by small boats near the shore in Baie de L'Isle Verte, St. Lawrence marsh, Quebec, Canada (48.003889 N, 69.425278 W) in May-June 2006. Collected samples were transported to the laboratory where Eurytemora affinis individuals were identified and sampled under the microscope.

Laboratory cultures and experiments
Four inbred lines of Eurytemora affinis (two each from the two populations) were generated through full-sibling mating for 30 generations. Two independent saltwater inbred lines (SW1 and SW2) were derived from the ancestral saline population in Baie de L’Isle Verte (Canada) and reared at their native salinity of 15 PSU. The two freshwater inbred lines (FW1 and FW2) were derived from the freshwater invading population in Racine Harbor (USA) and reared in Lake Michigan water (0 PSU, conductivity 300 μS/cm). In addition, reciprocal F1 crosses between freshwater and saline inbred lines were created and reared to test for allele-specific expression by comparing gene expression in parental lines and their F1 crosses.

Two replicate common-garden reaction norm experiments, each consisting of a 2 × 2 factorial design, were performed to compare patterns of the gene expression of the FW and SW inbred lines (see Materials and Methods in Posavi et al. 2020). Total RNA from whole bodies of 50 copepods (25 females and 25 males) per sample was extracted using Trizol reagent (Ambion RNA) and Qiagen RNeasy Mini Kit for purification (Qiagen cat. no. 74104). Extracted and purified RNA samples were stored at -80 degrees Celsius until sequencing. The strand-specific Illumina RNA-seq libraries (Parkhomchuk et al., 2009) of polyA purified mRNA were constructed using the TruSeq RNA Sample Prep kit (Illumina). Three biological replicates per inbred line were sequenced using the Illumina HiSeq 2000 platform in the Institute for Genome Sciences at the University of Maryland School of Medicine and generated 101-bp-long paired-end read data.

These data have important implications for understanding the evolutionary and physiological mechanisms of range expansions by some of the most widespread invaders in aquatic habitats.

Problem report
One replicate of the FW1 inbred line was excluded because of bacterial infection

Additional information
~ Detailed methods, results, and figures can be found in Posavi et al. (2020) (see Related Publications section).
~ The sequence data can be viewed under NCBI BioProject PRJNA278152 (see Related Datasets).

Data Processing Description

To assess the taxonomic composition and ensure the provenance of the sample, the RNA-seq reads were screened by the Institute for Genome Sciences' QC pipeline against a local installation of the NCBI nucleotide database. After that, each sample was run through the data processing pipeline to detect the presence of adaptor sequences or low read quality using FastQC (Andrews, 2010). Reads were trimmed with Trimmomatic version 032 (Bolger et al., 2014). On average, 3.5 × 107 paired-end (101 bp) reads per sample passed these filtering steps.

To quantify transcript (gene) expression levels, an expectation maximization approach employing the RSEM (RNA-seq by Expectation Maximization) package (Li & Dewey, 2011) was used. The Eurytemora affinis complex (Atlantic clade, aka E. carolleeae) draft genome served as the reference genome. Automated gene annotation of this genome was conducted at the Baylor College of Medicine Human Genome Sequencing Center within the i5K pilot project (using Maker2.2 following methods of Holt & Yandell, 2011) resulting in 29,783 gene models. For additional details on gene annotation, see Eyun et al. (2017).

To improve the automated gene annotation, the Cufflinks Tuxedo protocol (Trapnell et al., 2012) was used. The merged gene annotation file was used as input into RSEM to (a) build reference transcript sequences using the prepare-reference and (b) align RNA-seq reads to the reference transcripts and estimate gene and transcript abundances (using rsem-calculate-expression).

To map RNA-seq reads to the E. affinis complex genome, Bowtie2 was employed resulting in 16–30 million mapped paired-end reads. To verify the annotation of DE genes, the manual annotation of the E. affinis complex genome (using the Web Apollo platform on the i5k Workspace, https://i5k.nal.usda.gov/Eurytemora_affinis) was performed.

Structural gene annotation, generated by merging transcriptomes of 22 RNA samples, resulted in 37,827 putative genes. To increase the power to detect differential expression, genes with less than one count-per million (CPM) in at least two samples were filtered out, as were transcripts with best BLASTx matches to bacteria (n = 3), and transcripts without BLASTx hits. These filtering steps left 14,082 putative genes remaining for the differential gene expression analysis, all of which mapped to the E. affinis genome assembly. Subsequently, the normalization on 14,082 genes was performed, using the Trimmed Mean of M values method (Robinson & Oshlack, 2010), available in Bioconductor's edgeR package for R software (Chen et al., 2018; Robinson et al., 2009).

To identify significant differences in gene expression between saline and freshwater inbred lines (Goal 1) and between salinities (0 and 15 PSU) (Goal 2), statistical analyses using a generalized linear model (GLM) were performed. To detect DE genes, a negative binomial generalized linear model was used that accommodated the complex designs of the common-garden experiments (Bioconductor’s edgeR package with function glmQLFit with option robust= TRUE). To conduct the test for each genotype (inbred line) and salinity combination, the read counts were modeled as the result of the fixed effects of genotype (inbred line effect), salinity (0 and15 PSU), batch, and genotype-by-salinity interactions. For multiple hypothesis testing, we adjusted p-values using the Benjamini and Hochberg (1995) method with a false discovery rate (FDR) threshold of 0.05.

Software
(See also Related Publications section below)

Bowtie2 (Langmead & Salzberg, 2012)
BLASTx (Mukhopadhyay et al., 2017)
Cufflinks Tuxedo protocol (Trapnell et al., 2012)
edgeR (Chen et al., 2018; Robinson et al, 2009)
FastQC (Andrews, 2010)
Maker2.2 (Holt & Yandell, 2011)
RSEM (Li & Dewey, 2011)
Trimmomatic (Bolger et al., 2014)

BCO-DMO Processing
- Converted date to Y-M-D format
- Added information about inbred lines for better readability

Data Files

File
gene_expression.csv (Comma Separated Values (.csv), 8.64 KB) MD5:d28111ecea446ecf40dec8018e1dc61d Primary data file for dataset ID 883426

More Information about this dataset

Funding Sources

Funding Source	Award Number
NSF Division of Ocean Sciences	OCE-1658517

Deployments

None available yet

Instruments

Automated DNA Sequencer

Supplied Name: Illumina HiSeq 2000

Supplied Description:

Libraries were sequenced on an Illumina HiSeq platform at the University of Maryland, School of Medicine, Institute for Genome Sciences.

Instrument Type

Generic Name: Automated DNA Sequencer

Acronym: Automated Sequencer

Generic Description:

General term for a laboratory instrument used for deciphering the order of bases in a strand of DNA. Sanger sequencers detect fluorescence from different dyes that are used to identify the A, C, G, and T extension reactions. Contemporary or Pyrosequencer methods are based on detecting the activity of DNA polymerase (a DNA synthesizing enzyme) with another chemoluminescent enzyme. Essentially, the method allows sequencing of a single strand of DNA by synthesizing the complementary strand along it, one base pair at a time, and detecting which base was actually added at each step.

Plankton Net

Supplied Name:

Supplied Description:

Instrument Type

Generic Name: Plankton Net

Community Identifier: http://vocab.nerc.ac.uk/collection/L05/current/22/

Generic Description:

A Plankton Net is a generic term for a sampling net that is used to collect plankton. It is used only when detailed instrument documentation is not available.

Parameters

Date and time formats are described in python strftime() syntax.

Supplied Name	Supplied description	Supplied Units	Standard Name
Record_num	Record number	unitless	sampling_sequence
Culture_sample_date	Date of sampling for the inbred line and F1 crosses Format: %Y-%m-%d	unitless	date
Sample_name	Sample name indicating the E.affinis line, salinity PSU, and replicate number	unitless	sample
Organism	Organism that was raised and studied (copepod Eurytemora affinis)	unitless	sample_descrip
Sex	Sex of copepods comprising the sample; male, female, or pooled	unitless	sample_descrip
Female_inbred_line	Identifier of inbred line for the females in the sample; VA = SW1 (saline inbred line 1), VE = SW2 (saline inbred line 2) , RA = FW1 (freshwater inbred line 1), RB = FW2 (freshwater inbred line 2)	unitless	sample_descrip
female_ancestor_location	Home location of the female ancestor	unitless	site_descrip
Male_inbred_line	Identifier of inbred line for the males in the sample; VA = SW1 (saline inbred line 1), VE = SW2 (saline inbred line 2) , RA = FW1 (freshwater inbred line 1), RB = FW2 (freshwater inbred line 2)	unitless	sample_descrip
male_ancestor_location	Home location of the male ancestor	unitless	site_descrip
Experimental_culture	Culture identification of inbred lines	unitless	exp_id
Description	Experiment description	unitless	exp_type
Salinity_conditions_rearing	Salinity conditions in which the copepods were reared	unitless	treatment
PSU_culture	Practical salinity unit value for the culture water (or water in which the copepods were reared)	PSU	sal
BioProject	NCBI BioProject identifier	unitless	BioProject
SRA	NCBI Sequence Read Archive identifier for sample accession	unitless	SRA_sample
BioSample	NCBI BioSample number	unitless	BioSample

Status message

Database

Contribute Data

Dataset: Gene expression of Eurytemora affinis (saline and freshwater)