RNA-seq based whole transcriptome analysis of the cyclopoid copepod Paracyclopina nana focusing on xenobiotics metabolism

Transcriptomic and microbiome analyses of copepod Apocyclops royi in response to an AHPND-causing strain of Vibrio parahaemolyticus

Copepods are small crustaceans that live in microorganism-rich aquatic environments and provide a key supply of live food for fish and shellfish larviculture. To better understand the host-pathogen interaction between the copepod and Vibrio parahaemolyticus causing acute hepatopancreatic necrosis disease (VPAHPND), the comparative transcriptome and microbiome analyses were conducted in copepod Apocyclops royi-TH following VPAHPND infection. Transcriptome analysis identified a total of 836 differentially expressed genes, with 275 upregulated and 561 downregulated genes. Subsequent analysis showed that a total of 37 differentially expressed genes were associated with the innate immune system, including 16 upregulated genes related to Toll-like receptor signaling pathway, antimicrobial peptides, and stress response genes, and 21 downregulated genes associated with immunological modulators, signaling molecules, and apoptosis-related proteins. Analysis of the copepod microbiome following VPAHPND infection showed that the microbes changed significantly after bacterial infection, with a reduced alpha diversity accompanied by the increased level of Proteobacteria and decreased levels of Bdellovibrionota, Bacteroidota, and Verrucomicrobiota. The population of Vibrio genera were increased significantly, while several other genera, including Denitromonas, Nitrosomonas, Blastopirellula, Fusibacter, Alteromonas, KI89A_clade, and Ruegeria, were decreased significantly after infection. These findings suggest that VPAHPND infection has a significant impact on the immune defense and the composition of the copepod microbiota.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

RNA-seq transcriptome analysis and identification of the theromacin antimicrobial peptide of the copepod Apocyclops royi

Copepods, including Apocyclops royi, are small aquatic crustaceans and one of the important foods for fish and shellfish larvae. However, studies of the host-pathogen interactions and understanding of infectious disease in copepods are still very limited, yet they are likely to be a significant factor in the sustainable development of copepod aquaculture. In the present study, we performed de novo RNA sequence analysis of A. royi-TH (a Thai isolate of A. royi), which yielded 4.80 Gb bases of clean data and a total of 29,786 unigenes. Annotation was then performed by comparison against seven functional databases, yielding 17,617 (NR: 59.15%), 2,969 (NT: 9.97%), 15,023 (SwissProt: 50.44%), 14,543 (KOG: 48.82%), 15,077 (KEGG: 50.62%), 6,763(GO: 22.71%), and 15,841 (InterPro: 53.18%) unigenes. In comparison to the components of the shrimp Toll pathway, LGBP, Spätzle, Toll receptors, MyD88, Pelle, TRAF6, Dorsal, and Cactus homologs were successfully identified in A. royi-TH. Additionally, a novel antimicrobial peptide (Theromacin-like) was characterized in A. royi (ArTM-like). The ArTM-like ORF was 279 bp and predicted to encode for 92 amino acid residues, with a mature peptide of 75 amino acids and a molecular mass of 8.56 kDa. The genomic organization of the ArTM-like gene consisted of three exons and two introns. Expression analysis indicated that ArTM-like mRNA was abundantly expressed in copepodid and adult stages as an immune responsive gene after infection with the pathogenic Vibrio parahaemolyticus-(AHPND)-causing strain. Altogether, the knowledge obtained in this study will provide a basis for future functional studies of the molecular mechanisms in copepod immunity that may eventually be applied for disease prevention in copepod aquaculture.

First transcriptome of the copepod Gladioferens pectinatus subjected to chronic contaminant exposures

Contaminants are often at low concentrations in ecosystems and their effects on exposed organisms can occur over long periods of time and across multiple generations. Alterations to subcellular mechanistic pathways in response to exposure to contaminants can provide insights into mechanisms of toxicity that methods measuring higher levels of biological may miss. Analysis of the whole transcriptome can identify novel mechanisms of action leading to impacts in exposed biota. The aim of this study was to characterise how exposures to copper, benzophenone and diclofenac across multiple generations altered molecular expression pathways in the marine copepod Gladioferens pectinatus. Results of the study demonstrated differential gene expression was observed in cultures exposure to diclofenac (569), copper (449) and benzophenone (59). Pathways linked to stress, growth, cellular and metabolic processes were altered by exposure to all three contaminants with genes associated with oxidative stress and xenobiotic regulation also impacted. Protein kinase functioning, cytochrome P450, transcription, skeletal muscle contraction/relaxation, mitochondrial phosphate translocator, protein synthesis and mitochondrial methylation were all differentially expressed with all three chemicals. The results of the study also suggested that using dimethyl sulfoxide as a dispersant influenced the transcriptome and future research may want to investigate it's use in molecular studies. Data generated in this study provides a first look at transcriptomic response of G. pectinatus exposed to contaminants across multiple generations, future research is needed to validate the identified biomarkers and link these results to apical responses such as population growth to demonstrate the predictive capacity of molecular tools.

Transcriptome Analysis Reveals Olfactory System Expression Characteristics of Aquatic Snakes

Animal olfactory systems evolved with changes in habitat to detect odor cues from the environment. The aquatic environment, as a unique habitat, poses a formidable challenge for olfactory perception in animals, since the higher density and viscosity of water. The olfactory system in snakes is highly specialized, thus providing the opportunity to explore the adaptive evolution of such systems to unique habitats. To date, however, few studies have explored the changes in gene expression features in the olfactory systems of aquatic snakes. In this study, we carried out RNA sequencing of 26 olfactory tissue samples (vomeronasal organ and olfactory bulb) from two aquatic and two non-aquatic snake species to explore gene expression changes under the aquatic environment. Weighted gene co-expression network analysis showed significant differences in gene expression profiles between aquatic and non-aquatic habitats. The main olfactory systems of the aquatic and non-aquatic snakes were regulated by different genes. Among these genes, RELN may contribute to exploring gene expression changes under the aquatic environment by regulating the formation of inhibitory neurons in the granular cell layer and increasing the separation of neuronal patterns to correctly identify complex chemical information. The high expression of TRPC2 and V2R family genes in the accessory olfactory systems of aquatic snakes should enhance their ability to bind water-soluble odor molecules, and thus obtain more information in hydrophytic habitats. This work provides an important foundation for exploring the olfactory adaptation of snakes in special habitats.

Characteristics of RNA alternative splicing and its potential roles in ginsenoside biosynthesis in a single plant of ginseng, Panax ginseng C.A. Meyer

RNA alternative splicing (AS) is prevalent in higher organisms and plays a paramount role in biology; therefore, it is crucial to have comprehensive knowledge on AS to understand biology. However, knowledge is limited about how AS activates in a single plant and functions in a biological process. Ginseng is one of the most widely used medicinal herbs that is abundant in a number of medicinal bioactive components, especially ginsenosides. In this study, we sequenced the transcripts of 14 organs from a 4-year-old ginseng plant and quantified their ginsenoside contents. We identified AS genes by analyzing their transcripts with the ginseng genome and verified their AS events by PCR. The plant had a total of 13,863 AS genes subjected to 30,801 AS events with five mechanisms: skipped exon, retained intron, alternative 5'splice site, alternative 3' splice site, and mutually exclusive exon. The genes that were more conserved, had more exons, and/or expressed across organs were more likely to be subjected to AS. AS genes were enriched in over 500 GO terms in the plant even though the number of AS gene-enriched GO terms varied across organs. At least 24 AS genes were found to be involved in ginsenoside biosynthesis. These AS genes were significantly up-enriched and more likely to form a co-expression network, thus suggesting the functions of AS and correlations of the AS genes in the process. This study provides comprehensive insights into the molecular characteristics and biological functions of AS in a single plant; thus, helping better understand biology.

Four new Cyclopina (Copepoda, Cyclopinidae) from South Korea

Copepods are well studied in South Korea, with the exception of marine non-parasitic cyclopoids, and especially cyclopinids; only three species were found so far here, and only one of them is endemic. A survey of intertidal interstitial faunas from sandy beaches revealed four endemic members of the genus Cyclopina Claus, 1863, which represents the first record of the largest cyclopinid genus in South Korea. A detailed study of their morphology revealed numerous differences, including in rarely studied cuticular organs. Some of these micro-characters could easily be homologised and showed little intraspecific variability, which might prove invaluable for matching sexes and reconstructing phylogenetic relationships. Cyclopinabusanensissp. nov. is described from both sexes collected near Busan (South Coast of South Korea), and is most similar to the only congener from Japan: C.kiraensis Horomi, 1984. Cyclopinakoreanasp. nov. is described from both sexes collected near Gangneung (East Coast), and has no close relatives among currently known species. Cyclopinacurtijejusp. nov. is described from two females from Jeju (off South Coast); it is possibly closely related to C.smirnovi Herbst, 1982, but the latter is known from a single male from the Russian Far East. Cyclopinawidosp. nov. is described from both sexes from Wido (West Coast), and shows numerous reductions in segmentation and armature of appendages, most of them probably a consequence of its diminutive size. A table of 26 discrete and continuous characters commonly used in the taxonomy of this group is provided for 48 valid species and subspecies of Cyclopina.

Effects of different routes of exposure to metals on bioaccumulation and population growth of the cyclopoid copepod Paracyclopina nana

We examined effects of the three metals cadmium (Cd), copper (Cu), and nickel (Ni) on two subpopulations of the cyclopoid copepod Paracyclopina nana. We sought to investigate the effects of metal exposure on population growth and structure of P. nana and to understand the parameters affecting the metal bioaccumulation in copepods. A first experiment tested the hypothesis of competition between these metals in a mixture using a P. nana mass culture in 10 L beakers with the sublethal concentrations (1/3 of LC50) as determined for E. affinis. A second experiment pursued the same with a P. nana population which was adapted to a higher Cu concentration for several generations (226.9 ± 15.9 μg g^-1 dw Cu in copepods) and using the proper sublethal concentrations for P. nana. After 96 h of exposure, results from the first experiment showed a decreasing population growth and instead of an increasing metal accumulation in copepods. Cd also appeared to be more accumulated when it was alone, confirming the hypothesis of metal competition in mixture. Results from the second experiment revealed less marked effects. When metal concentrations increased in the treatment it decreased in copepods, indicating depuration activity in the population already adapted to metal exposure. This paper is the first one investigating the parameters affecting the bioaccumulation capacity of P. nana in response to metals. It offers a better understanding of copepod responses to metal contamination in a complex aquatic environment.

The transcriptome of the marine calanoid copepod Temora longicornis under heat stress and recovery

Understanding the impacts of global change in zooplankton communities is crucial, as alterations in the zooplankton communities can affect entire marine ecosystems. Despite the economic and ecological importance of the calanoid copepod Temora longicornis in the Belgian part of the North Sea, molecular data is still very limited for this species. Using HiSeq Illumina sequencing, we sequenced the whole transcriptome of T. longicornis, after being exposed to realistic temperatures of 14 and 17 °C. After both an acute (1 day) and a more sustained (5 days) thermal exposure to 17 °C, we investigated gene expression differences with animals exposed to 14 °C, which may be critical for the thermal acclimation and resilience of this copepod species. We also studied the possibility of a short term stress recovery of a heat shock. A total of 179,569 transcripts were yielded, of which 44,985 putative ORF transcripts were identified. These transcripts were subsequently annotated into roughly 22,000 genes based on known sequences using Gene Ontology (GO) and KEGG databases. Temora only showed a mild response to both the temperature and the duration of the exposure. We found that the expression of 27 transcripts varied significantly with an increase in temperature of 3 °C, of which eight transcripts were differentially expressed after acute exposure only. Gene set enrichment analysis revealed that, overall, T. longicornis was more impacted by a sustained thermal exposure, rather than an immediate (acute) exposure, with two times as many enriched GO terms in the sustained treatment. We also identified several general stress responses independent of exposure time, such as modified protein synthesis, energy mobilisation, cuticle and chaperone proteins. Finally, we highlighted candidate genes of a possible recovery from heat exposure, identifying similar terms as those enriched in the heat treatments, i.e. related to for example energy metabolism, cuticle genes and extracellular matrix. The data presented in this study provides the first transcriptome available for T. longicornis which can be used for future genomic studies.

De novo transcriptome assembly and differential gene expression analysis of the calanoid copepod Acartia tonsa exposed to nickel nanoparticles

The calanoid copepod Acartia tonsa is a reference species in standardized ecotoxicology bioassay. Despite this interest, there is a lack of knowledge on molecular responses of A. tonsa to contaminants. We generated a de novo assembled transcriptome of A. tonsa exposed 4 days to 8.5 and 17 mg/L nickel nanoparticles (NiNPs), which have been shown to reduce egg hatching success and larval survival but had no effects on the adults. Aims of our study were to 1) improve the knowledge on the molecular responses of A. tonsa copepod and 2) increase the genomic resources of this copepod for further identification of potential biomarkers of NP exposure. The de novo assembled transcriptome of A. tonsa consisted of 53,619 unigenes, which were further annotated to nr, GO, KOG and KEGG databases. In particular, most unigenes were assigned to Metabolic and Cellular processes (34-45%) GO terms, and to Human disease (28%) and Organismal systems (23%) KEGG categories. Comparison among treatments showed that 373 unigenes were differentially expressed in A. tonsa exposed to NiNPs at 8.5 and 17 mg/L, with respect to control. Most of these genes were downregulated and took part in ribosome biogenesis, translation and protein turnover, thus suggesting that NiNPs could affect the copepod ribosome synthesis machinery and functioning. Overall, our study highlights the potential of toxicogenomic approach in gaining more mechanistic and functional information about the mode of action of emerging compounds on marine organisms, for biomarker discovering in crustaceans.

Enrichment of genomic resources and identification of simple sequence repeats from medicinally important Clausena excavata

To broaden and delve into the genomic information of Clausena excavata, an important medicinal plant in many Asian countries, RNA sequencing (RNA-seq) analysis was performed and a total of 16,638 non-redundant unigenes (≥ 300 bp) with an average length of 755 bp were generated by de novo assembly from 17,580,456 trimmed clear reads. The functional categorization of the identified unigenes by a gene ontology (GO) term resulted in 2305 genes in the cellular component, 5577 in the biological processes, and 8056 in the molecular functions, respectively. The top sub-category in biological processes was the metabolic process with 4374 genes. Among annotated genes, 3006 were mapped to 123 metabolic pathways by the Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analysis tool. The search for simple sequence repeats (SSRs) resulted in 845 SSRs from 749 SSR-containing unigenes and the most abundant SSR motifs was AAG/CTT with 179 occurrences. Twelve SSR markers were tested for cross transferability among five Clausena species; eight of them exhibited polymorphism. Taken together, these data provide valuable resources for genomic or genetic studies of Clausena species and other relative studies. The transcriptome shotgun assembly data have been deposited at DDBJ/EMBL/GenBank under the accession GGEM00000000.

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.

Molecular Characterization of Copepod Photoreception

Copepod crustaceans are an abundant and ecologically significant group whose basic biology is guided by numerous visually guided behaviors. These behaviors are driven by copepod eyes, including naupliar eyes and Gicklhorn's organs, which vary widely in structure and function among species. Yet little is known about the molecular aspects of copepod vision. In this study we present a general overview of the molecular aspects of copepod vision by identifying phototransduction genes from newly generated and publicly available RNA-sequencing data and assemblies from 12 taxonomically diverse copepod species. We identify a set of 10 expressed transcripts that serve as a set of target genes for future studies of copepod phototransduction. Our more detailed evolutionary analyses of the opsin gene responsible for forming visual pigments found that all of the copepod species investigated express two main groups of opsins: middle-wavelength-sensitive (MWS) opsins and pteropsins. Additionally, there is evidence from a few species (e.g., Calanus finmarchicus, Eurytemora affinis, Paracyclopina nana, and Lernaea cyprinacea) for the expression of two additional groups of opsins-the peropsins and rhodopsin 7 (Rh7) opsins-at low levels or distinct developmental stages. An ontogenetic analysis of opsin expression in Calanus finmarchicus found the expression of a single dominant MWS opsin, as well as evidence for differences in expression across development in some MWS, pteropsin, and Rh7 opsins, with expression peaking in early naupliar through early copepodite stages.

Differential gene responses in the embryo of the green mussel Perna viridis exposed to dichlorodiphenyltrichloroethane (DDT)

The green-lipped mussel, Perna viridis, is considered to be an ideal indicator for marine environmental pollution. Dichlorodiphenyltrichloroethane (DDT), a typical persistent organic pollutant, is extensively distributed in marine environments. However, little is known about the toxic effects of DDT on the embryo of marine animals, particularly in marine bivalves. Using next-generation sequencing technology, we studied P. viridis embryo after DDT stress at the transcriptome level. A total of 99 202 unigenes were obtained based on the 2383 bp of unigene N50. These differentially expressed genes (DEGs) participated in the various molecular pathways of biological effects, including oxidative stress, detoxification, innate immunity and neurobehavioral disease. Quantitative real-time PCR was performed to verify the mRNA expression of several genes identified by differential gene expression (DGE) analysis. The results indicated that DDT was in induced a dose-dependent manner in the embryo of P. viridis, and most genes involved in oxidative stress and detoxification were up-regulated by DDT exposure; however, the immunity-related genes were down-regulated, except the genes involved in phagocytosis. Gene expression changes in embryo from P. viridis provide a preliminary basis to better understand the molecular toxic response mechanisms of embryo to DDT.

Identification of the Full 46 Cytochrome P450 (CYP) Complement and Modulation of CYP Expression in Response to Water-Accommodated Fractions of Crude Oil in the Cyclopoid Copepod Paracyclopina nana

The 46 cytochrome P450 (CYP) gene superfamily was identified in the marine copepod Paracyclopina nana after searching an RNA-seq database and comparing it with other copepod CYP gene families. To annotate the 46 Pn-CYP genes, a phylogenetic analysis of CYP genes was performed using a Bayesian method. Pn-CYP genes were separated into five different clans: CYP2, CYP3, CYP20, CYP26, and mitochondrial. Among these, the principal Pn-CYP genes involved in detoxification were identified by comparing them with those of the copepod Tigriopus japonicus and were examined with respect to their responses to exposure to a water-accommodated fraction (WAF) of crude oil and to the alkylated forms of two polycyclic aromatic hydrocarbons (PAHs; phenanthrene and fluorene). The expression of two Pn-CYP3027 genes (CYP3027F1 and CYP3027F2) was increased in response to WAF exposure and also was upregulated in response to the two alkylated PAHs. In particular, Pn-CYP3027F2 showed the most notable increase in response to 80% WAF exposure. These two responsive CYP genes (Pn-CYP3027F1 and CYP3027F2) were also phylogenetically clustered into the same clade of the WAF- and alkylated PAH-specific CYP genes of the copepod T. japonicus, suggesting that these CYP genes would be those chiefly involved in detoxification in response to WAF exposure in copepods. In this paper, we provide information on the copepod P. nana CYP gene superfamily and also speculate on its potential role in the detoxification of PAHs in marine copepods. Despite the nonlethality of WAF, Pn-CYP3027F2 was rapidly and significantly upregulated in response to WAF that may serve as a useful biomarker of 40% or higher concentration of WAF exposure. This paper will be helpful to better understand the molecular mechanistic events underlying the metabolism of environmental toxicants in copepods.