Mapping the relationship of white matter lesions to depression in multiple sclerosis.

Importance: Multiple sclerosis (MS) is an immune-mediated neurological disorder that affects nearly one million people in the United States. Up to 50% of patients with MS experience depression. Objective: To investigate how white matter network disruption is related to depression in MS. Design: Retrospective case-control study of participants who received research-quality 3-tesla neuroimaging as part of MS clinical care from 2010-2018. Analyses were performed from May 1 to September 30, 2022. Setting: Single-center academic medical specialty MS clinic. Participants: Participants with MS were identified via the electronic health record (EHR). All participants were diagnosed by an MS specialist and completed research-quality MRI at 3T. After excluding participants with poor image quality, 783 were included. Inclusion in the depression group (MS+Depression) required either: 1) ICD-10 depression diagnosis (F32-F34.*); 2) prescription of antidepressant medication; or 3) screening positive via Patient Health Questionnaire-2 (PHQ-2) or -9 (PHQ-9). Age- and sex-matched nondepressed comparators (MS-Depression) included persons with no depression diagnosis, no psychiatric medications, and were asymptomatic on PHQ-2/9. Exposure: Depression diagnosis. Main Outcomes and Measures: We first evaluated if lesions were preferentially located within the depression network compared to other brain regions. Next, we examined if MS+Depression patients had greater lesion burden, and if this was driven by lesions specifically in the depression network. Outcome measures were the burden of lesions (e.g., impacted fascicles) within a network and across the brain. Secondary measures included between-diagnosis lesion burden, stratified by brain network. Linear mixed-effects models were employed. Results: Three hundred-eighty participants met inclusion criteria, (232 MS+Depression: age[SD]=49[12], %females=86; 148 MS-Depression: age[SD]=47[13], %females=79). MS lesions preferentially affected fascicles within versus outside the depression network (ß=0.09, 95% CI=0.08-0.10, P<0.001). MS+Depression had more white matter lesion burden (ß=0.06, 95% CI=0.01-0.10, P=0.015); this was driven by lesions within the depression network (ß=0.02, 95% CI 0.003-0.040, P=0.020). Conclusions and Relevance: We provide new evidence supporting a relationship between white matter lesions and depression in MS. MS lesions disproportionately impacted fascicles in the depression network. MS+Depression had more disease than MS-Depression, which was driven by disease within the depression network. Future studies relating lesion location to personalized depression interventions are warranted.

De novo transcriptomes of six calanoid copepods (Crustacea): a resource for the discovery of novel genes.

This study presents eight new high-quality de novo transcriptomes from six co-occurring species of calanoid copepods, the first published for Neocalanus plumchrus, N. cristatus, Eucalanus bungii and Metridia pacifica and additional ones for N. flemingeri and Calanus marshallae. They are ecologically-important members of sub-arctic North Pacific marine zooplankton communities. 'Omics data for this diverse and numerous taxonomic group are sparse and difficult to obtain. Total RNA from single individuals was used to construct gene libraries that were sequenced on an Illumina Next-Seq platform. Quality filtered reads were assembled with Trinity software and validated using multiple criteria. The study's primary purpose is to provide a resource for gene expression studies. The integrated database can be used for quantitative inter- and intra-species comparisons of gene expression patterns across biological processes. An example of an additional use is provided for discovering novel and evolutionarily-significant proteins within the Calanoida. A workflow was designed to find and characterize unannotated transcripts with homologies across de novo assemblies that have also been shown to be eco-responsive.

Physiological acclimatization in high-latitude zooplankton.

How individual organisms adapt to nonoptimal conditions through physiological acclimatization is central to predicting the consequences of unusual abiotic and biotic conditions such as those produced by marine heat waves. The Northeast Pacific, including the Gulf of Alaska, experienced an extreme warming event (2014-2016, "The Blob") that affected all trophic levels and led to large-scale changes in the community. The marine copepod Neocalanus flemingeri is a key member of the subarctic Pacific pelagic ecosystem. During the spring phytoplankton bloom this copepod builds substantial lipid stores as it prepares for its nonfeeding adult phase. A 3-year comparison of gene expression profiles of copepods collected in Prince William Sound in the Gulf of Alaska between 2015 and 2017 included two high-temperature years (2015 and 2016) and one year with very low phytoplankton abundances (2016). The largest differences in gene expression were between high and low chlorophyll years, and not between warm and cool years. The observed gene expression patterns were indicative of physiological acclimatization. The predominant signal in 2016 was the down-regulation of genes involved in glycolysis and its incoming pathways, consistent with the modulation of metabolic rates in response to prolonged low food conditions. Despite the down-regulation of genes involved in metabolism, there was no evidence of suppression of protein synthesis based on gene expression or behavioural activity. Genes involved in muscle function were up-regulated, and the copepods were actively swimming and responsive to stimuli at collection. However, genes involved in fatty acid metabolism were down-regulated in 2016, suggesting reduced lipid accumulation.

Post-diapause transcriptomic restarts: insight from a high-latitude copepod.

BACKGROUND: Diapause is a seasonal dormancy that allows organisms to survive unfavorable conditions and optimizes the timing of reproduction and growth. Emergence from diapause reverses the state of arrested development and metabolic suppression returning the organism to an active state. The physiological mechanisms that regulate the transition from diapause to post-diapause are still unknown. In this study, this transition has been characterized for the sub-arctic calanoid copepod Neocalanus flemingeri, a key crustacean zooplankter that supports the highly productive North Pacific fisheries. Transcriptional profiling of females, determined over a two-week time series starting with diapausing females collected from > 400 m depth, characterized the molecular mechanisms that regulate the post-diapause trajectory. RESULTS: A complex set of transitions in relative gene expression defined the transcriptomic changes from diapause to post-diapause. Despite low temperatures (5-6 °C), the switch from a "diapause" to a "post-diapause" transcriptional profile occurred within 12 h of the termination stimulus. Transcriptional changes signaling the end of diapause were activated within one-hour post collection and included the up-regulation of genes involved in the 20E cascade pathway, the TCA cycle and RNA metabolism in combination with the down-regulation of genes associated with chromatin silencing. By 12 h, females exhibited a post-diapause phenotype characterized by the up-regulation of genes involved in cell division, cell differentiation and multiple developmental processes. By seven days post collection, the reproductive program was fully activated as indicated by up-regulation of genes involved in oogenesis and energy metabolism, processes that were enriched among the differentially expressed genes. CONCLUSIONS: The analysis revealed a finely structured, precisely orchestrated sequence of transcriptional changes that led to rapid changes in the activation of biological processes paving the way to the successful completion of the reproductive program. Our findings lead to new hypotheses related to potentially universal mechanisms that terminate diapause before an organism can resume its developmental program.

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.

First De Novo Transcriptome of the Copepod Rhincalanus gigas from Antarctic Waters.

Antarctic waters are the largest almost untapped diversified resource of our planet. Molecular resources for Antarctic organisms are very limited and mostly represented by sequences used for species genotyping. In this study, we present the first transcriptome for the copepod Rhincalanus gigas, one of the predominant zooplankton species of Antarctic waters. This transcriptome represents also the first molecular resource for an eucalanoid copepod. The transcriptome is of high quality and completeness. The presence of three predicted genes encoding antifreeze proteins and gene duplication within the glutathione metabolism pathway are suggested as possible adaptations to cope with this harsh environment. The R. gigas transcriptome represents a powerful new resource for investigating the molecular basis associated with polar biological processes and ecology.

t-Distributed Stochastic Neighbor Embedding (t-SNE): A tool for eco-physiological transcriptomic analysis.

High-throughput RNA sequencing (RNA-Seq) has transformed the ecophysiological assessment of individual plankton species and communities. However, the technology generates complex data consisting of millions of short-read sequences that can be difficult to analyze and interpret. New bioinformatics workflows are needed to guide experimentation, environmental sampling, and to develop and test hypotheses. One complexity-reducing tool that has been used successfully in other fields is "t-distributed Stochastic Neighbor Embedding" (t-SNE). Its application to transcriptomic data from marine pelagic and benthic systems has yet to be explored. The present study demonstrates an application for evaluating RNA-Seq data using previously published, conventionally analyzed studies on the copepods Calanus finmarchicus and Neocalanus flemingeri. In one application, gene expression profiles were compared among different developmental stages. In another, they were compared among experimental conditions. In a third, they were compared among environmental samples from different locations. The profile categories identified by t-SNE were validated by reference to published results using differential gene expression and Gene Ontology (GO) analyses. The analyses demonstrate how individual samples can be evaluated for differences in global gene expression, as well as differences in expression related to specific biological processes, such as lipid metabolism and responses to stress. As RNA-Seq data from plankton species and communities become more common, t-SNE analysis should provide a powerful tool for determining trends and classifying samples into groups with similar transcriptional physiology, independent of collection site or time.

Regional heterogeneity impacts gene expression in the subarctic zooplankter Neocalanus flemingeri in the northern Gulf of Alaska.

Marine pelagic species are being increasingly challenged by environmental change. Their ability to persist will depend on their capacity for physiological acclimatization. Little is known about limits of physiological plasticity in key species at the base of the food web. Here we investigate the capacity for acclimatization in the copepod Neocalanus flemingeri, which inhabits the Gulf of Alaska, a heterogeneous and highly seasonal environment. RNA-Seq analysis of field-collected pre-adults identified large regional differences in expression of genes involved in metabolic and developmental processes and response to stressors. We found that lipid synthesis genes were up-regulated in individuals from Prince William Sound and down-regulated in the Gulf of Alaska. Up-regulation of lipid catabolic genes in offshore individuals suggests they are experiencing nutritional deficits. The expression differences demonstrate physiological plasticity in response to a steep gradient in food availability. Our transcriptional analysis reveals mechanisms of acclimatization that likely contribute to the observed resilience of this population.

Physiological characterization of the emergence from diapause: A transcriptomics approach.

Organisms inhabiting high-latitude environments have evolved adaptations, such as diapause to time reproduction and growth to optimize their survival. However, the physiological regulation of the timing of complex life histories is poorly understood, particularly for marine copepods, that diapause at depth. A member of the pelagic community of the sub-Arctic Pacific Ocean, Neocalanus flemingeri enters diapause in June. Egg production occurs in winter/spring. In order to characterize the transition from diapause to egg release, females were collected in late September from 400-700 m depth, incubated in the dark at 4-5 °C and sampled for RNASeq at weekly intervals. The diapause phenotype showed down-regulation of protein turnover and up-regulation of stress genes. Activation of the reproductive program was marked by the up-regulation of genes involved in germline development. Thereafter, progress through phases of oocyte development could be linked to changes in gene expression. At 5 weeks, females showed up-regulation of spermatogenesis, indicating that stored sperm had been in a quiescent stage and completed their maturation inside the female. Gene expression profiles provide a framework to stage field-collected females. The 7-week progression from diapause to late oogenesis suggests that females typically spawning in January initiated the reproductive program in November.

Molecular evidence for an intrinsic circadian pacemaker in the cardiac ganglion of the American lobster, Homarus americanus - Is diel cycling of heartbeat frequency controlled by a peripheral clock system?

Whether cardiac output in decapod crustaceans is under circadian control has long been debated, with mixed evidence for and against the hypothesis. Moreover, the locus of the clock system controlling cardiac activity, if it is under circadian control, is unknown. However, a report that the crayfish heart in organ culture maintains a circadian oscillation in heartbeat frequency suggests the presence of a peripheral pacemaker within the cardiac neuromuscular system itself. Because the decapod heart is neurogenic, with contractions controlled by the five motor and four premotor neurons that make up the cardiac ganglion (CG), a likely locus for a circadian clock is the CG itself. Here, a CG-specific transcriptome was generated for the lobster, Homarus americanus, and was used to assess the presence/absence of transcripts encoding putative clock-related proteins in the ganglion. Using known Homarus brain/eyestalk ganglia clock-related proteins as queries, BLAST searches of the CG transcriptome were conducted for the five proteins that form the core clock, i.e., clock, cryptochrome 2, cycle, period and timeless, as well as for a variety of clock-associated, clock input pathway and clock output pathway proteins. With the exception of pigment dispersing hormone receptor [PDHR], a putative clock output pathway protein, one or more transcripts encoding each of the proteins searched for were identified from the CG assembly; no PDHR-encoding transcripts were found. RT-PCR confirmed the expression of all core clock transcripts in multiple independent CG cDNAs; RNA-Seq data suggest that both the motor and premotor neurons could contribute to the cellular locus of a pacemaker. These data provide support for the possible existence of an intrinsic circadian clock in the H. americanus CG, and form a foundation for guiding future anatomical, molecular and physiological investigations of circadian signaling in the lobster cardiac neuromuscular system.

Circadian signaling in Homarus americanus: Region-specific de novo assembled transcriptomes show that both the brain and eyestalk ganglia possess the molecular components of a putative clock system.

Essentially all organisms exhibit recurring patterns of physiology/behavior that oscillate with a period of ~24-h and are synchronized to the solar day. Crustaceans are no exception, with robust circadian rhythms having been documented in many members of this arthropod subphylum. However, little is known about the molecular underpinnings of their circadian rhythmicity. Moreover, the location of the crustacean central clock has not been firmly established, although both the brain and eyestalk ganglia have been hypothesized as loci. The American lobster, Homarus americanus, is known to exhibit multiple circadian rhythms, and immunodetection data suggest that its central clock is located within the eyestalk ganglia rather than in the brain. Here, brain- and eyestalk ganglia-specific transcriptomes were generated and used to assess the presence/absence of transcripts encoding the commonly recognized protein components of arthropod circadian signaling systems in these two regions of the lobster central nervous system. Transcripts encoding putative homologs of the core clock proteins clock, cryptochrome 2, cycle, period and timeless were found in both the brain and eyestalk ganglia assemblies, as were transcripts encoding similar complements of putative clock-associated, clock input pathway and clock output pathway proteins. The presence and identity of transcripts encoding core clock proteins in both regions were confirmed using PCR. These findings suggest that both the brain and eyestalk ganglia possess all of the molecular components needed for the establishment of a circadian signaling system. Whether the brain and eyestalk clocks are independent of one another or represent a single timekeeping system remains to be determined. Interestingly, while most of the proteins deduced from the identified transcripts are shared by both the brain and eyestalk ganglia, assembly-specific isoforms were also identified, e.g., several period variants, suggesting the possibility of region-specific variation in clock function, especially if the brain and eyestalk clocks represent independent oscillators.

Prediction of a peptidome for the ecotoxicological model Hyalella azteca (Crustacea; Amphipoda) using a de novo assembled transcriptome.

Due to its sensitivity to many environmental and anthropogenic stressors, including a wide range of chemical compounds, Hyalella azteca, a freshwater amphipod, has emerged as one of the most commonly used invertebrates for ecotoxicological assessment.Peptidergic signaling systems are key components in the control of organism-environment interactions, and there is a growing literature suggesting that they are targets of a number of aquatic toxicants.Interestingly, and despite its model species status in the field of ecotoxicology, little is known about the peptide hormones of H. azteca.Here, a transcriptome was produced for this species using the de novo assembler Trinity and mined for sequences encoding putative peptide precursors; the transcriptome was assembled from 460,291,636 raw reads and consists of 133,486 unique transcripts.Seventy-six sequences encoding peptide pre/preprohormones were identified from this transcriptome, allowing for the prediction of 202 distinct peptides, which included members of the allatostatin A, allatostatin B, allatostatin C, allatotropin, bursicon, CCHamide, corazonin, crustacean cardioactive peptide, crustacean hyperglycemic hormone/molt-inhibiting hormone, ecdysis-triggering hormone, eclosion hormone, elevenin, FMRFamide-like peptide, glycoprotein hormone, GSEFLamide, inotocin, leucokinin, myosuppressin, neuropeptide F, orcokinin, orcomyotropin, pigment dispersing hormone, proctolin, pyrokinin, red pigment concentrating hormone, RYamide, short neuropeptide F, SIFamide, sulfakinin, tachykinin-related peptide and trissin families.These peptides expand the known peptidome for H. azteca approximately nine-fold, forming a strong foundation for future studies of peptidergic control, including disruption by aquatic toxicants, in this important ecotoxicological model.

De novo transcriptome assembly of the calanoid copepod Neocalanus flemingeri: A new resource for emergence from diapause.

Copepods, small planktonic crustaceans, are key links between primary producers and upper trophic levels, including many economically important fishes. In the subarctic North Pacific, the life cycle of copepods like Neocalanus flemingeri includes an ontogenetic migration to depth followed by a period of diapause (a type of dormancy) characterized by arrested development and low metabolic activity. The end of diapause is marked by the production of the first brood of eggs. Recent temperature anomalies in the North Pacific have raised concerns about potential negative effects on N. flemingeri. Since diapause is a developmental program, its progress can be tracked using through global gene expression. Thus, a reference transcriptome was developed as a first step towards physiological profiling of diapausing females using high-throughput Illumina sequencing. The de novo transcriptome, the first for this species was designed to investigate the diapause period. RNA-Seq reads were obtained for dormant to reproductive N. flemingeri females. A high quality de novo transcriptome was obtained by first assembling reads from each individual using Trinity software followed by clustering with CAP3 Assembly Program. This assembly consisted of 140,841transcripts (contigs). Bench-marking universal single-copy orthologs analysis identified 85% of core eukaryotic genes, with 79% predicted to be complete. Comparison with other calanoid transcriptomes confirmed its quality and degree of completeness. Trinity assembly of reads originating from multiple individuals led to fragmentation. Thus, the workflow applied here differed from the one recommended by Trinity, but was required to obtain a good assembly.

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, Kane'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.

Complementary mechanisms for neurotoxin resistance in a copepod.

Toxin resistance is a recurring evolutionary response by predators feeding on toxic prey. These adaptations impact physiological interaction and community ecology. Mechanisms for resistance vary depending on the predator and the nature of the toxin. Potent neurotoxins like tetrodotoxin (TTX) and saxitoxin (STX) that are highly toxic to humans and other vertebrates, target conserved voltage-gated sodium channels (NaV) of nerve and muscle, causing paralysis. The copepod Calanus finmarchicus consumes the STX-producing dinoflagellate, Alexandrium fundyense with no effect on survival. Using transcriptomic approaches to search for the mechanism that confers resistance in C. finmarchicus, we identified splice variants of NaVs that were predicted to be toxin resistant. These were co-expressed with putatively non-resistant form in all developmental stages. However its expression was unresponsive to toxin challenge nor was there any up-regulation of genes involved in multi-xenobiotic resistance (MXR) or detoxification (phases I or II). Instead, adults consistently regulated genes encoding digestive enzymes, possibly to complement channel resistance by limiting toxin assimilation via the digestive process. The nauplii, which were more susceptible to STX, did not regulate these enzymes. This study demonstrates how deep-sequencing technology can elucidate multiple mechanisms of toxin resistance concurrently, revealing the linkages between molecular/cellular adaptations and the ecology of an organism.

Prediction of a neuropeptidome for the eyestalk ganglia of the lobster Homarus americanus using a tissue-specific de novo assembled transcriptome.

In silico transcriptome mining is a powerful tool for crustacean peptidome prediction. Using homology-based BLAST searches and a simple bioinformatics workflow, large peptidomes have recently been predicted for a variety of crustaceans, including the lobster, Homarus americanus. Interestingly, no in silico studies have been conducted on the eyestalk ganglia (lamina ganglionaris, medulla externa, medulla interna and medulla terminalis) of the lobster, although the eyestalk is the location of a major neuroendocrine complex, i.e., the X-organ-sinus gland system. Here, an H. americanus eyestalk ganglia-specific transcriptome was produced using the de novo assembler Trinity. This transcriptome was generated from 130,973,220 Illumina reads and consists of 147,542 unique contigs. Eighty-nine neuropeptide-encoding transcripts were identified from this dataset, allowing for the deduction of 62 distinct pre/preprohormones. Two hundred sixty-two neuropeptides were predicted from this set of precursors; the peptides include members of the adipokinetic hormone-corazonin-like peptide, allatostatin A, allatostatin B, allatostatin C, bursicon α, CCHamide, corazonin, crustacean cardioactive peptide, crustacean hyperglycemic hormone (CHH), CHH precursor-related peptide, diuretic hormone 31, diuretic hormone 44, eclosion hormone, elevenin, FMRFamide-like peptide, glycoprotein hormone α2, glycoprotein hormone ß5, GSEFLamide, intocin, leucokinin, molt-inhibiting hormone, myosuppressin, neuroparsin, neuropeptide F, orcokinin, orcomyotropin, pigment dispersing hormone, proctolin, pyrokinin, red pigment concentrating hormone, RYamide, short neuropeptide F, SIFamide, sulfakinin, tachykinin-related peptide and trissin families. The predicted peptides expand the H. americanus eyestalk ganglia neuropeptidome approximately 7-fold, and include 78 peptides new to the lobster. The transcriptome and predicted neuropeptidome described here provide new resources for investigating peptidergic signaling within/from the lobster eyestalk ganglia.

Transcriptomic responses of the calanoid copepod Calanus finmarchicus to the saxitoxin producing dinoflagellate Alexandrium fundyense.

In the Gulf of Maine, the copepod Calanus finmarchicus co-occurs with the neurotoxin-producing dinoflagellate, Alexandrium fundyense. The copepod is resistant to this toxic alga, but little is known about other effects. Gene expression profiles were used to investigate the physiological response of females feeding for two and five days on a control diet or a diet containing either a low or a high dose of A. fundyense. The physiological responses to the two experimental diets were similar, but changed between the time points. At 5-days the response was characterized by down-regulated genes involved in energy metabolism. Detoxification was not a major component of the response. Instead, genes involved in digestion were consistently regulated, suggesting that food assimilation may have been affected. Thus, predicted increases in the frequency of blooms of A. fundyense could affect C. finmarchicus populations by changing the individuals' energy budget and reducing their ability to build lipid reserves.

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.

Diffusible gas transmitter signaling in the copepod crustacean Calanus finmarchicus: identification of the biosynthetic enzymes of nitric oxide (NO), carbon monoxide (CO) and hydrogen sulfide (H2S) using a de novo assembled transcriptome.

Neurochemical signaling is a major component of physiological/behavioral control throughout the animal kingdom. Gas transmitters are perhaps the most ancient class of molecules used by nervous systems for chemical communication. Three gases are generally recognized as being produced by neurons: nitric oxide (NO), carbon monoxide (CO) and hydrogen sulfide (H2S). As part of an ongoing effort to identify and characterize the neurochemical signaling systems of the copepod Calanus finmarchicus, the biomass dominant zooplankton in much of the North Atlantic Ocean, we have mined a de novo assembled transcriptome for sequences encoding the neuronal biosynthetic enzymes of these gases, i.e. nitric oxide synthase (NOS), heme oxygenase (HO) and cystathionine ß-synthase (CBS), respectively. Using Drosophila proteins as queries, two NOS-, one HO-, and one CBS-encoding transcripts were identified. Reverse BLAST and structural analyses of the deduced proteins suggest that each is a true member of its respective enzyme family. RNA-Seq data collected from embryos, early nauplii, late nauplii, early copepodites, late copepodites and adults revealed the expression of each transcript to be stage specific: one NOS restricted primarily to the embryo and the other was absent in the embryo but expressed in all other stages, no CBS expression in the embryo, but present in all other stages, and HO expressed across all developmental stages. Given the importance of gas transmitters in the regulatory control of a number of physiological processes, these data open opportunities for investigating the roles these proteins play under different life-stage and environmental conditions in this ecologically important species.