Comprehensive transcriptome study to develop molecular resources of the copepod Calanus sinicus for their potential ecological applications

Insights into the species evolution of Calanus copepods in the northern seas revealed by de novo transcriptome sequencing

Copepods of the zooplankton genus Calanus play a key role in marine ecosystems in the northern seas. Although being among the most studied organisms on Earth, due to their ecological importance, genomic resources for Calanus spp. remain scarce, mostly due to their large genome size (from 6 to 12 Gbps). As an alternative to whole-genome sequencing in Calanus spp., we sequenced and de novo assembled transcriptomes of five Calanus species: Calanus glacialis, C. hyperboreus, C. marshallae, C. pacificus, and C. helgolandicus. Functional assignment of protein families based on clusters of orthologous genes (COG) and gene ontology (GO) annotations showed analogous patterns of protein functions across species. Phylogenetic analyses using maximum likelihood (ML) of 191 protein-coding genes mined from RNA-seq data fully resolved evolutionary relationships among seven Calanus species investigated (five species sequenced for this study and two species with published datasets), with gene and site concordance factors showing that 109 out of 191 protein-coding genes support a separation between three groups: the C. finmarchicus group (including C. finmarchicus, C. glacialis, and C. marshallae), the C. helgolandicus group (including C. helgolandicus, C. sinicus, and C. pacificus) and the monophyletic C. hyperboreus group. The tree topology obtained in ML analyses was similar to a previously proposed phylogeny based on morphological criteria and cleared certain ambiguities from past studies on evolutionary relationships among Calanus species.

De Novo Transcriptome Assembly and Gene Expression Profiling of the Copepod Calanus helgolandicus Feeding on the PUA-Producing Diatom Skeletonema marinoi

Diatoms are the dominant component of the marine phytoplankton. Several diatoms produce secondary metabolites, namely oxylipins, with teratogenic effects on their main predators, crustacean copepods. Our study reports the de novo assembled transcriptome of the calanoid copepod Calanus helgolandicus feeding on the oxylipin-producing diatom Skeletonema marinoi. Differential expression analysis was also performed between copepod females exposed to the diatom and the control flagellate Prorocentrum minimum, which does not produce oxylipins. Our results showed that transcripts involved in carbohydrate, amino acid, folate and methionine metabolism, embryogenesis, and response to stimulus were differentially expressed in the two conditions. Expression of 27 selected genes belonging to these functional categories was also analyzed by RT-qPCR in C. helgolandicus females exposed to a mixed solution of the oxylipins heptadienal and octadienal at the concentration of 10 µM, 15 µM, and 20 µM. The results confirmed differential expression analysis, with up-regulation of genes involved in stress response and down-regulation of genes associated with folate and methionine metabolism, embryogenesis, and signaling. Overall, we offer new insights on the mechanism of action of oxylipins on maternally-induced embryo abnormality. Our results may also help identify biomarker genes associated with diatom-related reproductive failure in the natural copepod population at sea.

Prevalence and proliferation of antibiotic resistance genes in the subtropical mangrove wetland ecosystem of South China Sea

The emerging pollutants antibiotic resistance genes (ARGs) are prevalent in aquatic environments such as estuary. Coastal mangrove ecosystems always serve as natural wetlands for receiving sewage which always carry ARGs. Currently, the research considering ARG distribution in mangrove ecosystems gains more interest. In this work, we investigated the diversity of ARGs in an urban estuary containing mangrove and nonmangrove areas of the South China Sea. A total of 163 ARGs that classified into 22 resistance types and six resistance mechanisms were found. ARG abundance of the samples in the estuary is between 0.144 and 0.203. This is within the general range of Chinese estuaries. The difference analysis showed that abundances of total ARGs, six most abundant ARGs (mtrA, rpoB, rpoC, rpsL, ef‐Tu, and parY), the most abundant resistance types (elfamycin, multidrug, and peptide), and the most abundant resistance mechanism (target alteration) were significantly lower in mangrove sediment than that in nonmangrove sediment (p < 0.05). Network and partial redundancy analysis showed that sediment properties and mobile genetic elements were the most influential factors impacting ARG distribution rather than microbial community. The two factors collectively explain 51.22% of the differences of ARG distribution. Our study indicated that mangrove sediments have the capacity to remove ARGs. This work provides a research paradigm for analysis of ARG prevalence and proliferation in the subtropical marine coastal mangrove ecosystem.

A deep transcriptomic resource for the copepod crustacean Labidocera madurae: A potential indicator species for assessing near shore ecosystem health

Coral reef ecosystems of many sub-tropical and tropical marine coastal environments have suffered significant degradation from anthropogenic sources. Research to inform management strategies that mitigate stressors and promote a healthy ecosystem has focused on the ecology and physiology of coral reefs and associated organisms. Few studies focus on the surrounding pelagic communities, which are equally important to ecosystem function. Zooplankton, often dominated by small crustaceans such as copepods, is an important food source for invertebrates and fishes, especially larval fishes. The reef-associated zooplankton includes a sub-neustonic copepod family that could serve as an indicator species for the community. Here, we describe the generation of a de novo transcriptome for one such copepod, Labidocera madurae, a pontellid from an intensively-studied coral reef ecosystem, Kāne'ohe Bay, Oahu, Hawai'i. The transcriptome was assembled using high-throughput sequence data obtained from whole organisms. It comprised 211,002 unique transcripts, including 72,391 with coding regions. It was assessed for quality and completeness using multiple workflows. Bench-marking-universal-single-copy-orthologs (BUSCO) analysis identified transcripts for 88% of expected eukaryotic core proteins. Targeted gene-discovery analyses included searches for transcripts coding full-length "giant" proteins (>4,000 amino acids), proteins and splice variants of voltage-gated sodium channels, and proteins involved in the circadian signaling pathway. Four different reference transcriptomes were generated and compared for the detection of differential gene expression between copepodites and adult females; 6,229 genes were consistently identified as differentially expressed between the two regardless of reference. Automated bioinformatics analyses and targeted manual gene curation suggest that the de novo assembled L. madurae transcriptome is of high quality and completeness. This transcriptome provides a new resource for assessing the global physiological status of a planktonic species inhabiting a coral reef ecosystem that is subjected to multiple anthropogenic stressors. The workflows provide a template for generating and assessing transcriptomes in other non-model species.

Mortality and transcriptional effects of inorganic mercury in the marine copepod Calanus finmarchicus

Inorganic mercury (Hg) is highly toxic to organisms including crustaceans and displays multiple toxic modes of action (MoA). The main aim of this investigation was to assess the acute and sublethal toxicity mediated by mercury chloride (HgCl₂) in the marine copepod Calanus finmarchicus. A combination of short-term static studies to determine acute toxicity and a transcriptional investigation to characterize the sublethal MoA of HgCl₂ were conducted with an in-house continuous culture of C. finmarchicus. Transcriptional changes were determined by a custom 6.6 k C. finmarchicus Agilent oligonucleotide microarray and quantitative RT-PCR analysis. Data demonstrate that HgCl₂ produced a concentration- and time-dependent reduction in survival (NOEC_{48 h} = 6.9 μg/L [Hg²⁺] and LC₅₀ of 279, 73, 48, and 34 µg/L [Hg²⁺] after 24, 48, 72, and 96 h, respectively) and that exposure to sublethal concentrations of HgCl₂ (5 μg/L [Hg2⁺]) induced differential expression of 98 features (probes) on the microarray. Gene ontology (GO) and toxicological pathway analyses suggested that the main MOA were (1) uncoupling of mitochondrial oxidative phosphorylation (OXPHOS) and ATP production, (2) oxidative stress and macromolecular damage, (3) inactivation of cellular enzymes, (4) induction of cellular apoptosis and autophagocytosis, (5) over-excitation of glutamate receptors (neurotoxicity), (6) disruption of calcium homeostasis and signaling, and (7) modulation of nuclear receptor activity involved in vitamin D receptor signaling. Quantitative RT-PCR analysis verified that oligoarray performed reliably in terms of specificity and response, thus demonstrating that Hg²⁺ exerts multiple potential MoA in C. finmarchicus.

Exploring developmental gene toolkit and associated pathways in a potential new model crustacean using transcriptomic analysis

The crustaceans are one of the largest, most diverse, and most successful groups of invertebrates. The diversity among the crustaceans is also reflected in embryonic development models. However, the molecular genetics that regulates embryonic development is not known in those crustaceans that have a short germ-band development with superficial cleavage, such as Macrobrachium olfersi. This species is a freshwater decapod and has great potential to become a model for developmental biology, as well as for evolutionary and environmental studies. To obtain sequence data of M. olfersi from an embryonic developmental perspective, we performed de novo assembly and annotation of the embryonic transcriptome. Using a pooling strategy of total RNA, paired-end Illumina sequencing, and assembly with multiple k-mers, a total of 25,636,097 pair reads were generated. In total, 99,751 unigenes were identified, and 20,893 of these returned a Blastx hit. KEGG pathway analysis mapped a total of 6866 unigenes related to 129 metabolic pathways. In general, 21,845 unigenes were assigned to gene ontology (GO) categories: molecular function (19,604), cellular components (10,254), and biological processes (13,841). Of these, 2142 unigenes were assigned to the developmental process category. More specifically, a total of 35 homologs of embryonic development toolkit genes were identified, which included maternal effect (one gene), gap (six), pair-rule (six), segment polarity (seven), Hox (four), Wnt (eight), and dorsoventral patterning genes (three). In addition, genes of developmental pathways were found, including TGF-β, Wnt, Notch, MAPK, Hedgehog, Jak-STAT, VEGF, and ecdysteroid-inducible nuclear receptors. RT-PCR analysis of eight genes related to embryonic development from gastrulation to late morphogenesis/organogenesis confirmed the applicability of the transcriptome analysis.

Identification of functional genes involved in Cd(2+) response of Chinese surf clam (Mactra chinensis) through transcriptome sequencing

The Chinese surf clam Mactra chinensis is an economically important bivalve species in the coastal waters of Liaoning and Shandong Province, China. In this study, we carried out transcriptome sequencing to develop molecular resources for M. chinensis and conducted an acute test of Cd(2+) stimulation through quantitative real-time PCR (qRT-PCR) to analyze the relative expression of six functional genes. A total of 100,839 transcripts and 56,712 unigenes were obtained from 39.9 million filtered reads and 21,305 unigenes were annotated by hitting against NCBI database. According to the results of qRT-PCR, heat shock protein 22 (Hsp22) and cytochrome P450 (CYP450(2C31)) were inhibited in the low concentration, and induced in the high concentration of Cd(2+); thioredoxin peroxidase (TPx-A) was at normal level in low concentration, but induced in high concentration of Cd(2+); glutathione peroxidase A (GPA), glutathione peroxidase 1 (GPA1) and Mn superoxide dismutase gene (MnSOD) were down-regulated when exposed to any treatment groups. Expression levels of the six functional genes following Cd(2+) exposure indicated that these genes were linked to environmental stress. Moreover, the present work enriched the molecule genetic data of M. chinensis.

Glutathione S-Transferase (GST) Gene Diversity in the Crustacean Calanus finmarchicus--Contributors to Cellular Detoxification

Detoxification is a fundamental cellular stress defense mechanism, which allows an organism to survive or even thrive in the presence of environmental toxins and/or pollutants. The glutathione S-transferase (GST) superfamily is a set of enzymes involved in the detoxification process. This highly diverse protein superfamily is characterized by multiple gene duplications, with over 40 GST genes reported in some insects. However, less is known about the GST superfamily in marine organisms, including crustaceans. The availability of two de novo transcriptomes for the copepod, Calanus finmarchicus, provided an opportunity for an in depth study of the GST superfamily in a marine crustacean. The transcriptomes were searched for putative GST-encoding transcripts using known GST proteins from three arthropods as queries. The identified transcripts were then translated into proteins, analyzed for structural domains, and annotated using reciprocal BLAST analysis. Mining the two transcriptomes yielded a total of 41 predicted GST proteins belonging to the cytosolic, mitochondrial or microsomal classes. Phylogenetic analysis of the cytosolic GSTs validated their annotation into six different subclasses. The predicted proteins are likely to represent the products of distinct genes, suggesting that the diversity of GSTs in C. finmarchicus exceeds or rivals that described for insects. Analysis of relative gene expression in different developmental stages indicated low levels of GST expression in embryos, and relatively high expression in late copepodites and adult females for several cytosolic GSTs. A diverse diet and complex life history are factors that might be driving the multiplicity of GSTs in C. finmarchicus, as this copepod is commonly exposed to a variety of natural toxins. Hence, diversity in detoxification pathway proteins may well be key to their survival.

Effects of the oxylipin-producing diatom Skeletonema marinoi on gene expression levels of the calanoid copepod Calanus sinicus

Diatoms are eukaryotic unicellular plants that constitute one of the major components of marine phytoplankton, comprising up to 40% of annual productivity at sea and representing 25% of global carbon-fixation. Diatoms have traditionally been considered a preferential food for zooplankton grazers such as copepods, but, in the last two decades, this beneficial role has been challenged after the discovery that many species of diatoms produce toxic metabolites, collectively termed oxylipins, that induce reproductive failure in zooplankton grazers. Diatoms are the dominant natural diet of Calanus sinicus, a cold-temperate calanoid copepod that supports secondary production of important fisheries in the shelf ecosystems of the Northwest Pacific Ocean, Yellow Sea, Sea of Japan and South China Sea. In this study, the effect of the oxylipin-producing diatom Skeletonema marinoi on C. sinicus has been evaluated by analyzing expression level changes of genes involved in defense and detoxification systems. Results show that C. sinicus is more resistant to a diet of this diatom species in terms of gene expression patterns, compared to the congeneric species Calanus helgolandicus which is an important constituent of the temperate waters of the Atlantic Ocean and northern Mediterranean Sea. These findings contribute to the better understanding of genetic and/or phenotypic flexibility of copepod species and their capabilities to cope with stress by identifying molecular markers (such as stress and detoxification genes) as biosensors for environmental perturbations (e.g. toxins and contaminants) affecting marine copepods.