Revealing the profound influence of diapause on gene expression: Insights from the annual transcriptome of the copepod Calanus finmarchicus

Annual rhythms are observed in living organisms with numerous ecological implications. In the zooplanktonic copepod Calanus finmarchicus, such rhythms are crucial regarding its phenology, body lipid accumulation, and global carbon storage. Climate change drives annual biological rhythms out of phase with the prevailing environmental conditions with yet unknown but potentially catastrophic consequences. However, the molecular dynamics underlying phenology are still poorly described. In a rhythmic analysis of C. finmarchicus annual gene expression, results reveal that more than 90% of the transcriptome shows significant annual rhythms, with abrupt and dramatic upheaval between the active and diapause life cycle states. This work explores the implication of the circadian clock in the annual timing, which may control epigenetic mechanisms to profoundly modulate gene expression in response to calendar time. Results also suggest an increased light sensitivity during diapause that would ensure the photoperiodic entrainment of the endogenous annual clock.

Annual transcriptome of a key zooplankton species, the copepod Calanus finmarchicus

The copepod Calanus finmarchicus (Crustacea, Copepoda) is a key zooplanktonic species with a crucial position in the North Atlantic food web and significant contributor to ocean carbon flux. Like many other high latitude animals, it has evolved a programmed arrested development called diapause to cope with long periods of limited food supply, while growth and reproduction are timed to take advantage of seasonal peaks in primary production. However, anthropogenic warming is inducing changes in the expected timing of phytoplankton blooms, suggesting phenological mismatches with negative consequences for the N. Atlantic ecosystem. While diapause mechanisms are mainly studied in terrestrial arthropods, specifically on laboratory model species, such as the fruit fly Drosophila, the molecular investigations of annual rhythms in wild marine species remain fragmentary. Here we performed a rigorous year-long monthly sampling campaign of C. finmarchicus in a Scottish Loch (UK; 56.45°N, 5.18°W) to generate an annual transcriptome. The mRNA of 36 samples (monthly triplicate of 25 individuals) have been deeply sequenced with an average depth of 137 ± 4 million reads (mean ± SE) per sample, aligned to the reference transcriptome, and filtered. We detail the quality assessment of the datasets and provide a high-quality resource for the investigation of wild annual transcriptomic rhythms (35,357 components) in a key diapausing zooplanktonic species.

Diapause vs. reproductive programs: transcriptional phenotypes in a keystone copepod

Many arthropods undergo a seasonal dormancy termed "diapause" to optimize timing of reproduction in highly seasonal environments. In the North Atlantic, the copepod Calanus finmarchicus completes one to three generations annually with some individuals maturing into adults, while others interrupt their development to enter diapause. It is unknown which, why and when individuals enter the diapause program. Transcriptomic data from copepods on known programs were analyzed using dimensionality reduction of gene expression and functional analyses to identify program-specific genes and biological processes. These analyses elucidated physiological differences and established protocols that distinguish between programs. Differences in gene expression were associated with maturation of individuals on the reproductive program, while those on the diapause program showed little change over time. Only two of six filters effectively separated copepods by developmental program. The first one included all genes annotated to RNA metabolism and this was confirmed using differential gene expression analysis. The second filter identified 54 differentially expressed genes that were consistently up-regulated in individuals on the diapause program in comparison with those on the reproductive program. Annotated to oogenesis, RNA metabolism and fatty acid biosynthesis, these genes are both indicators for diapause preparation and good candidates for functional studies.

De novo Transcriptome Assembly of Myllocerinus aurolineatus Voss in Tea Plants

Myllocerinus aurolineatus Voss is a species of the insecta class in the arthropod. In this study, we first observed and identified M. aurolineatus Voss in tea plants in Guizhou, China, where it caused severe quantity and quality losses in tea plants. Knowledge on M. aurolineatus Voss genome is inadequate, especially for biological or functional research. We performed the first transcriptome sequencing by using the Illumina Hiseq™ technique on M. aurolineatus Voss. Over 55.9 million high-quality paired-end reads were generated and assembled into 69,439 unigenes using the Trinity short read software, resulting in a cluster of 1,207 bp of the N50 length. A total of 69,439 genes were predicted by BLAST to known proteins in the NCBI database and were distributed into Gene Ontology (20,190), eukaryotic complete genomes (12,488), and the Kyoto Encyclopedia of Genes and Genomes (3,170). We also identified 96,790 single-nucleotide polymorphisms and 13,121 simple sequence repeats in these unigenes. Our transcriptome data provide a useful resource for future functional studies of M. aurolineatus Voss for dispersal control in tea plants.

Transcriptomic analysis elucidates the molecular processes associated with hydrogen peroxide-induced diapause termination in Artemia-encysted embryos

Treatment with hydrogen peroxide (H₂O₂) raises the hatching rate through the development and diapause termination of Artemia cysts. To comprehend the upstream genetic regulation of diapause termination activated by exterior H₂O₂ elements, an Illumina RNA-seq analysis was performed to recognize and assess comparative transcript amounts to explore the genetic regulation of H₂O₂ in starting the diapause termination of cysts in Artemia salina. We examined three groupings treated with no H₂O₂ (control), 180 μM H₂O₂ (low) and 1800 μM H₂O₂ (high). The results showed a total of 114,057 unigenes were identified, 41.22% of which were functionally annotated in at least one particular database. When compared to control group, 34 and 98 differentially expressed genes (DEGs) were upregulated in 180 μM and 1800 μM H₂O₂ treatments, respectively. On the other hand, 162 and 30 DEGs were downregulated in the 180 μM and 1800 μM H₂O₂ treatments, respectively. Cluster analysis of DEGs demonstrated significant patterns among these types of 3 groups. GO and KEGG enrichment analysis showed the DEGs involved in the regulation of blood coagulation (GO: 0030193; GO: 0050818), regulation of wound healing (GO:0061041), regulation of hemostasis (GO: 1900046), antigen processing and presentation (KO04612), the Hippo signaling pathway (KO04391), as well as the MAPK signaling pathway (KO04010). This research helped to define the diapause-related transcriptomes of Artemia cysts using RNA-seq technology, which might fill up a gap in the prevailing body of knowledge.

De novo transcriptome assembly of the Southern Ocean copepod Rhincalanus gigas sheds light on developmental changes in gene expression

Copepods are small crustaceans that dominate most zooplankton communities in terms of both abundance and biomass. In the polar oceans, a subset of large lipid-storing copepods occupy central positions in the food web because of their important role in linking phytoplankton and microzooplankton with higher trophic levels. In this paper, we generated a high-quality de novo transcriptome for Rhincalanus gigas, the largest-and among the most abundant-of the Southern Ocean copepods. We then conducted transcriptional profiling to characterize the developmental transition between late-stage juveniles and adult females. We found that juvenile R. gigas substantially upregulate lipid synthesis and glycolysis pathways relative to females, as part of a developmental gene expression program that also implicates processes such as muscle growth, chitin formation, and ion transport. This study provides the first transcriptional profile of a developmental transition within Rhincalanus gigas or any endemic Southern Ocean copepod, thereby extending our understanding of copepod molecular physiology.

Lipid metabolism in Calanus finmarchicus is sensitive to variations in predation risk and food availability

Late developmental stages of the marine copepods in the genus Calanus can spend extended periods in a dormant stage (diapause) that is preceded by the accumulation of large lipid stores. We assessed how lipid metabolism during development from the C4 stage to adult is altered in response to predation risk and varying food availability, to ultimately understand more of the metabolic processes during development in Calanus copepods. We used RNA sequencing to assess if perceived predation risk in combination with varied food availability affects expression of genes associated with lipid metabolism and diapause preparation in C. finmarchicus. The lipid metabolism response to predation risk differed depending on food availability, time and life stage. Predation risk caused upregulation of lipid catabolism with high food, and downregulation with low food. Under low food conditions, predation risk disrupted lipid accumulation. The copepods showed no clear signs of diapause preparation, supporting earlier observations of the importance of multiple environmental cues in inducing diapause in C. finmarchicus. This study demonstrates that lipid metabolism is a sensitive endpoint for the interacting environmental effects of predation pressure and food availability. As diapause may be controlled by lipid accumulation, our findings may contribute towards understanding processes that can ultimately influence diapause timing.

The β-oxidation pathway is downregulated during diapause termination in Calanus copepods

Calanus copepods are keystone species in marine ecosystems, mainly due to their high lipid content, which is a nutritious food source for e.g. juvenile fish. Accumulated lipids are catabolized to meet energy requirements during dormancy (diapause), which occurs during the last copepodite stage (C5). The current knowledge of lipid degradation pathways during diapause termination is limited. We characterized changes in lipid fullness and generated transcriptional profiles in C5s during termination of diapause and progression towards adulthood. Lipid fullness of C5s declined linearly during developmental progression, but more β-oxidation genes were upregulated in early C5s compared to late C5s and adults. We identified four possible master regulators of energy metabolism, which all were generally upregulated in early C5s, compared to late C5s and adults. We discovered that one of two enzymes in the carnitine shuttle is absent from the calanoid copepod lineage. Based on the geographical location of the sampling site, the field-samples were initially presumed to consist of C. finmarchicus. However, the identification of C. glacialis in some samples underlines the need for performing molecular analyses to reliably identify Calanus species. Our findings contributes to a better understanding of molecular events occurring during diapause and diapause termination in calanoid copepods.

Diapause within the Context of Life-History Strategies in Calanid Copepods (Calanoida: Crustacea)

Post-embryonic diapause in copepods is an adaptation that allows species in the copepod family Calanidae to thrive in high-latitude environments by transforming a short spring phytoplankton bloom into large numbers of lipid-rich individuals capable of surviving a long period of starvation. The copepods, with their high-energy lipid reservoirs, are a critical food source for higher trophic levels, making the Calanidae a key component of high-latitude marine ecosystems. The physiological ecology of the developmental program remains poorly understood. However, new studies using high-throughput RNA sequencing approaches are giving detailed access to physiological status by generating gene expression profiles for both field-collected and laboratory-incubated individuals. These are beginning to characterize the diapause phenotype, elucidate the transcriptional and physiological progression through the diapause program, and illustrate the effects of organism-environment interactions. This paper reviews gene expression profiling studies on the life cycle and diapause program of Neocalanus flemingeri Miller (1988) that were conducted as part of a long-term observation program in the northern Gulf of Alaska. It summarizes recent findings and relates them to the ecology of this species and to that of other calanids.

Noncoding RNA Regulation of Dormant States in Evolutionarily Diverse Animals

Dormancy is evolutionarily widespread and can take many forms, including diapause, dauer formation, estivation, and hibernation. Each type of dormancy is characterized by distinct features; but accumulating evidence suggests that each is regulated by some common processes, often referred to as a common "toolkit" of regulatory mechanisms, that likely include noncoding RNAs that regulate gene expression. Noncoding RNAs, especially microRNAs, are well-known regulators of biological processes associated with numerous dormancy-related processes, including cell cycle progression, cell growth and proliferation, developmental timing, metabolism, and environmental stress tolerance. This review provides a summary of our current understanding of noncoding RNAs and their involvement in regulating dormancy.

A Crude Awakening: Effects of Crude Oil on Lipid Metabolism in Calanoid Copepods Terminating Diapause

Calanus finmarchicus and Calanus glacialis are keystone zooplankton species in North Atlantic and Arctic marine ecosystems because they form a link in the trophic transfer of nutritious lipids from phytoplankton to predators on higher trophic levels. These calanoid copepods spend several months of the year in deep waters in a dormant state called diapause, after which they emerge in surface waters to feed and reproduce during the spring phytoplankton bloom. Disruption of diapause timing could have dramatic consequences for marine ecosystems. In the present study, Calanus C5 copepodites were collected in a Norwegian fjord during diapause and were subsequently experimentally exposed to the water-soluble fraction of a naphthenic North Sea crude oil during diapause termination. The copepods were sampled repeatedly while progressing toward adulthood and were analyzed for utilization of lipid stores and for differential expression of genes involved in lipid metabolism. Our results indicate that water-soluble fraction exposure led to a temporary pause in lipid catabolism, suggested by (i) slower utilization of lipid stores in water-soluble fraction-exposed C5 copepodites and (ii) more genes in the β-oxidation pathway being downregulated in water-soluble fraction-exposed C5 copepodites than in the control C5 copepodites. Because lipid content and/or composition may be an important trigger for termination of diapause, our results imply that the timing of diapause termination and subsequent migration to the surface may be delayed if copepods are exposed to oil pollution during diapause or diapause termination. This delay could have detrimental effects on ecosystem dynamics.

ArHsp40 and ArHsp40-2 contribute to stress tolerance and longevity in Artemia franciscana, but only ArHsp40 influences diapause entry

Embryos of the crustacean Artemia franciscana develop either ovoviviparously or oviparously, yielding swimming larvae (nauplii) or encysted gastrulae (cysts), respectively. Nauplii moult several times and become adults whereas cysts enter diapause, a state of dormancy characterized by exceptionally low metabolism and high stress tolerance. Synthesis of molecular chaperones such as the J-domain proteins ArHsp40 and ArHsp40-2 occurs during embryo development and post-diapause growth of A. franciscana and they influence development and stress tolerance. To further investigate J-domain protein function, ArHsp40 and ArHsp40-2 were each knocked down by RNA interference. Reductions in ArHsp40 and ArHsp40-2 had no effect on adult survival, time to release of cysts and nauplii from females and first-brood size. However, knockdown of both A. franciscana J-domain proteins reduced the longevity and heat tolerance of nauplii, with the loss of ArHsp40 having a greater effect. The knockdown of ArHsp40, but not of ArHsp40-2, caused approximately 50% of cysts to abort diapause entry and hatch without exposure to an exogenous signal such as low temperature and/or desiccation. Cysts lacking ArHsp40 that entered diapause exhibited decreased stress tolerance as did cysts with reduced ArHsp40-2, the latter to a lesser degree. The longevity of nauplii hatching prematurely from cysts was less than for nauplii arising by other means. The results expand our understanding of Hsp40 function in A. franciscana stress tolerance and development, especially during diapause, and they provide the first example of a molecular chaperone that influences diapause entry.

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.

The Physiology and Ecology of Diapause in Marine Copepods

Diapause is a type of dormancy that requires preparation, typically precedes the onset of unfavorable conditions, and necessitates a period of arrest before development can proceed. Two ecologically important groups of copepods have incorporated diapausing stages into their life histories. In freshwater, estuarine, and coastal environments, species within the Centropagoidea superfamily can produce resting eggs containing embryos that may be quiescent, diapausing, or in some intermediate state. Resting eggs sink into the sediments, remain viable over months to years, and form a reservoir from which the planktonic population is reestablished. In coastal and oceanic environments, copepods within the Calanidae and Eucalanidae families can enter diapause during late juvenile (copepodid) or adult stages. These copepods accumulate large amounts of lipids before they migrate into deep water and diapause for several months. Through respiration, diapausing copepods may sequester more carbon in the deep ocean than any other biogeochemical process, and changes in diapause phenology associated with climate change (particularly reduction in diapause duration) could have a significant impact not only on regional ecosystems, but on global climate as well.

Resources and Recommendations for Using Transcriptomics to Address Grand Challenges in Comparative Biology

High-throughput RNA sequencing (RNA-seq) technology has become an important tool for studying physiological responses of organisms to changes in their environment. De novo assembly of RNA-seq data has allowed researchers to create a comprehensive catalog of genes expressed in a tissue and to quantify their expression without a complete genome sequence. The contributions from the "Tapping the Power of Crustacean Transcriptomics to Address Grand Challenges in Comparative Biology" symposium in this issue show the successes and limitations of using RNA-seq in the study of crustaceans. In conjunction with the symposium, the Animal Genome to Phenome Research Coordination Network collated comments from participants at the meeting regarding the challenges encountered when using transcriptomics in their research. Input came from novices and experts ranging from graduate students to principal investigators. Many were unaware of the bioinformatics analysis resources currently available on the CyVerse platform. Our analysis of community responses led to three recommendations for advancing the field: (1) integration of genomic and RNA-seq sequence assemblies for crustacean gene annotation and comparative expression; (2) development of methodologies for the functional analysis of genes; and (3) information and training exchange among laboratories for transmission of best practices. The field lacks the methods for manipulating tissue-specific gene expression. The decapod crustacean research community should consider the cherry shrimp, Neocaridina denticulata, as a decapod model for the application of transgenic tools for functional genomics. This would require a multi-investigator effort.

Tapping the Power of Crustacean Transcriptomics to Address Grand Challenges in Comparative Biology: An Introduction to the Symposium

Crustaceans, and decapods in particular (i.e., crabs, shrimp, and lobsters), are a diverse and ecologically and commercially important group of organisms. Understanding responses to abiotic and biotic factors is critical for developing best practices in aquaculture and assessing the effects of changing environments on the biology of these important animals. A relatively small number of decapod crustacean species have been intensively studied at the molecular level; the availability, experimental tractability, and economic relevance factor into the selection of a particular species as a model. Transcriptomics, using high-throughput next generation sequencing (NGS, coupled with RNA sequencing or RNA-seq) is revolutionizing crustacean biology. The 11 symposium papers in this volume illustrate how RNA-seq is being used to study stress response, molting and limb regeneration, immunity and disease, reproduction and development, neurobiology, and ecology and evolution. This symposium occurred on the 10th anniversary of the symposium, "Genomic and Proteomic Approaches to Crustacean Biology", held at the Society for Integrative and Comparative Biology 2006 meeting. Two participants in the 2006 symposium, the late Paul Gross and David Towle, were recognized as leaders who pioneered the use of molecular techniques that would ultimately foster the transcriptomics research reviewed in this volume. RNA-seq is a powerful tool for hypothesis-driven research, as well as an engine for discovery. It has eclipsed the technologies available in 2006, such as microarrays, expressed sequence tags, and subtractive hybridization screening, as the millions of "reads" from NGS enable researchers to de novo assemble a comprehensive transcriptome without a complete genome sequence. The symposium series concludes with a policy paper that gives an overview of the resources available and makes recommendations for developing better tools for functional annotation and pathway and network analysis in organisms in which the genome is not available or is incomplete.