Transcriptomic Analysis of Differentially Expressed Genes During Larval Development of Rapana venosa by Digital Gene Expression Profiling

Identification and characterization of olfactory receptor genes and olfactory perception in rapa whelk Rapana venosa (Valenciennes, 1846) during larval settlement and metamorphosis

The rapa whelk Rapana venosa, an economically important marine fishery resource in China but a major invader all over the world, changes from a phytophagous to a carnivorous form following settlement and metamorphosis. However, the low settlement and metamorphosis rates (<1%) of larvae limit the abundance of R. venosa. This critical step (settlement and metamorphosis) remains poorly characterized but may be related to how larvae perceive the presence of shellfish, their new source of food. Here, we report that larvae may use olfactory perception to sense shellfish. Olfactory receptor (OR) genes are involved in odor sensing in animals. We identified a total of 463 OR genes, which could be grouped into nine clades based on phylogenetic analysis. When assessing the attraction of larvae at different developmental stages to oyster odor, R. venosa showed active settlement and metamorphosis behavior only at the J4 stage (competent larva, 1000-1500 μm shell length) and in the presence of shellfish odor at the same time. Expression of OR gene family members differed between stage 2 (four-spiral whorl stage) and stage 1 (single- to three-spiral whorl stage), indicating significant changes in the olfactory system during larval development. These findings broaden our understanding of olfactory perception, settlement, and metamorphosis in gastropods and can be used to improve R. venosa harvesting, as well as the sustainable development and utilization of this resource.

Temporally balanced selection during development of larval Pacific oysters (Crassostrea gigas) inherently preserves genetic diversity within offspring

Balancing selection is one of the mechanisms which has been proposed to explain the maintenance of genetic diversity in species across generations. For species with large populations and complex life histories, however, heterogeneous selection pressures may create a scenario in which the net effects of selection are balanced across developmental stages. With replicated cultures and a pooled sequencing approach, we show that genotype-dependent mortality in larvae of the Pacific oyster (Crassostrea gigas) is largely temporally dynamic and inconsistently in favour of a single genotype or allelic variant at each locus. Overall, the patterns of genetic change we observe to be taking place are more complex than what would be expected under classical examples of additive or dominant genetic interactions. They are also not easily explained by our current understanding of the effects of genetic load. Collectively, temporally heterogeneous selection pressures across different larval developmental stages may act to maintain genetic diversity, while also inherently sheltering genetic load within oyster populations.

Development and Interrogation of a Transcriptomic Resource for the Giant Triton Snail (Charonia tritonis)

Gastropod molluscs are among the most abundant species that inhabit coral reef ecosystems. Many are specialist predators, along with the giant triton snail Charonia tritonis (Linnaeus, 1758) whose diet consists of Acanthaster planci (crown-of-thorns starfish), a corallivore known to consume enormous quantities of reef-building coral. C. tritonis are considered vulnerable due to overexploitation, and a decline in their populations is believed to have contributed to recurring A. planci population outbreaks. Aquaculture is considered one approach that could help restore natural populations of C. tritonis and mitigate coral loss; however, numerous questions remain unanswered regarding their life cycle, including the molecular factors that regulate their reproduction and development. In this study, we have established a reference C. tritonis transcriptome derived from developmental stages (embryo and veliger) and adult tissues. This was used to identify genes associated with cell signalling, such as neuropeptides and G protein-coupled receptors (GPCRs), involved in endocrine and olfactory signalling. A comparison of developmental stages showed that several neuropeptide precursors are exclusively expressed in post-hatch veligers and functional analysis found that FFamide stimulated a significant (20.3%) increase in larval heart rate. GPCRs unique to veligers, and a diversity of rhodopsin-like GPCRs located within adult cephalic tentacles, all represent candidate olfactory receptors. In addition, the cytochrome P450 superfamily, which participates in the biosynthesis and degradation of steroid hormones and lipids, was also found to be expanded with at least 91 genes annotated, mostly in gill tissue. These findings further progress our understanding of C. tritonis with possible application in developing aquaculture methods.

Changes in Symbiotic Microbiota and Immune Responses in Early Development Stages of Rapana venosa (Valenciennes, 1846) Provide Insights Into Immune System Development in Gastropods

The symbiotic microbiota can stimulate modulation of immune system, which also can promote immune system mature in critical developmental periods. In this study, we have investigated the symbiotic microbiota in Rapana venosa at five early development stages using Illumina high-throughput sequencing, and detected immune responses in larvae. Analysis of the symbiotic microbiota sequences identified that the most abundant phylum was Proteobacteria. Beta diversity analysis indicated that the structure of the symbiotic microbiota dramatically shifted in early development stages. The abundance of immune-related KEGG Orthologs (KOs) also increased in competent larval (J4, 30-day post-hatching) and postlarval after 3 days of metamorphosis (Y5, 33-day post-hatching) stages. Acid phosphatase activity decreased significantly in the Y5 stage, and alkaline phosphatase activity also at a lower level in Y5 stage, whereas lysozyme activities exhibited no remarkable change. Also, the activities of catalase and superoxide dismutase activities decreased dramatically during early development stages of R. venosa. Dramatic changes in the symbiotic microbiota and the immune response mainly occurred in the initially hatched veliger (C1), competent larval (J4) and postlarval (Y5) stages, during which the hosts might experience substantial environmental changes or changes in physiological structure and function. These findings expand our understanding of the stage-specific symbiotic microbiota in R. venosa and the close association between immune system and symbiotic microbiota in mollusks, however, the specific relationship may need more researches are needed to investigated in the future.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Transcriptome analysis of air-breathing land slug, Incilaria fruhstorferi reveals functional insights into growth, immunity, and reproduction

Background

Incilaria (= Meghimatium) fruhstorferi is an air-breathing land slug found in restricted habitats of Japan, Taiwan and selected provinces of South Korea (Jeju, Chuncheon, Busan, and Deokjeokdo). The species is on a decline due to depletion of forest cover, predation by natural enemies, and collection. To facilitate the conservation of the species, it is important to decide on a number of traits related to growth, immunity and reproduction addressing fitness advantage of the species.

Results

The visceral mass transcriptome of I. fruhstorferi was enabled using the Illumina HiSeq 4000 sequencing platform. According to BUSCO (Benchmarking Universal Single-Copy Orthologs) method, the transcriptome was considered complete with 91.8% of ortholog genes present (Single: 70.7%; Duplicated: 21.1%). A total of 96.79% of the raw read sequences were processed as clean reads. TransDecoder identified 197,271 contigs that contained candidate-coding regions. Of a total of 50,230 unigenes, 34,470 (68.62% of the total unigenes) annotated to homologous proteins in the Protostome database (PANM-DB). The GO term and KEGG pathway analysis indicated genes involved in metabolism, phosphatidylinositol signalling system, aminobenzoate degradation, and T-cell receptor signalling pathway. Many genes associated with molluscan innate immunity were categorized under pathogen recognition receptor, TLR signalling pathway, MyD88 dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related. The reproduction-associated unigenes showed homology to protein fem-1, spermatogenesis-associated protein, sperm associated antigen, and testis expressed sequences, among others. In addition, we identified key growth-related genes categorized under somatotrophic axis, muscle growth, chitinases and collagens. A total of 4822 Simple Sequence Repeats (SSRs) were also identified from the unigene sequences of I. fruhstorferi.

Conclusions

This is the first available genomic information for non-model land slug, I. fruhstorferi focusing on genes related to growth, immunity, and reproduction, with additional focus on microsatellites and repeating elements. The transcriptome provides access to greater number of traits of unknown relevance in the species that could be exploited for in-depth analyses of evolutionary plasticity and making informed choices during conservation planning. This would be appropriate for understanding the dynamics of the species on a priority basis considering the ecological, health, and social benefits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5526-3) contains supplementary material, which is available to authorized users.

Identification of genes associated with ricinoleic acid accumulation in Hiptage benghalensis via transcriptome analysis

BACKGROUND

Ricinoleic acid is a high-value hydroxy fatty acid with broad industrial applications. Hiptage benghalensis seed oil contains a high amount of ricinoleic acid (~ 80%) and represents an emerging source of this unusual fatty acid. However, the mechanism of ricinoleic acid accumulation in H. benghalensis is yet to be explored at the molecular level, which hampers the exploration of its potential in ricinoleic acid production.

RESULTS

To explore the molecular mechanism of ricinoleic acid biosynthesis and regulation, H. benghalensis seeds were harvested at five developing stages (13, 16, 19, 22, and 25 days after pollination) for lipid analysis. The results revealed that the rapid accumulation of ricinoleic acid occurred at the early-mid-seed development stages (16-22 days after pollination). Subsequently, the gene transcription profiles of the developing seeds were characterized via a comprehensive transcriptome analysis with second-generation sequencing and single-molecule real-time sequencing. Differential expression patterns were identified in 12,555 transcripts, including 71 enzymes in lipid metabolic pathways, 246 putative transcription factors (TFs) and 124 long noncoding RNAs (lncRNAs). Twelve genes involved in diverse lipid metabolism pathways, including fatty acid biosynthesis and modification (hydroxylation), lipid traffic, triacylglycerol assembly, acyl editing and oil-body formation, displayed high expression levels and consistent expression patterns with ricinoleic acid accumulation in the developing seeds, suggesting their primary roles in ricinoleic acid production. Subsequent co-expression network analysis identified 57 TFs and 35 lncRNAs, which are putatively involved in the regulation of ricinoleic acid biosynthesis. The transcriptome data were further validated by analyzing the expression profiles of key enzyme-encoding genes, TFs and lncRNAs with quantitative real-time PCR. Finally, a network of genes associated with ricinoleic acid accumulation in H. benghalensis was established.

CONCLUSIONS

This study was the first step toward the understating of the molecular mechanisms of ricinoleic acid biosynthesis and oil accumulation in H. benghalensis seeds and identified a pool of novel genes regulating ricinoleic acid accumulation. The results set a foundation for developing H. benghalensis into a novel ricinoleic acid feedstock at the transcriptomic level and provided valuable candidate genes for improving ricinoleic acid production in other plants.

Identification of genes associated with ricinoleic acid accumulation in Hiptage benghalensis via transcriptome analysis

BACKGROUND

Ricinoleic acid is a high-value hydroxy fatty acid with broad industrial applications. Hiptage benghalensis seed oil contains a high amount of ricinoleic acid (~ 80%) and represents an emerging source of this unusual fatty acid. However, the mechanism of ricinoleic acid accumulation in H. benghalensis is yet to be explored at the molecular level, which hampers the exploration of its potential in ricinoleic acid production.

RESULTS

To explore the molecular mechanism of ricinoleic acid biosynthesis and regulation, H. benghalensis seeds were harvested at five developing stages (13, 16, 19, 22, and 25 days after pollination) for lipid analysis. The results revealed that the rapid accumulation of ricinoleic acid occurred at the early-mid-seed development stages (16-22 days after pollination). Subsequently, the gene transcription profiles of the developing seeds were characterized via a comprehensive transcriptome analysis with second-generation sequencing and single-molecule real-time sequencing. Differential expression patterns were identified in 12,555 transcripts, including 71 enzymes in lipid metabolic pathways, 246 putative transcription factors (TFs) and 124 long noncoding RNAs (lncRNAs). Twelve genes involved in diverse lipid metabolism pathways, including fatty acid biosynthesis and modification (hydroxylation), lipid traffic, triacylglycerol assembly, acyl editing and oil-body formation, displayed high expression levels and consistent expression patterns with ricinoleic acid accumulation in the developing seeds, suggesting their primary roles in ricinoleic acid production. Subsequent co-expression network analysis identified 57 TFs and 35 lncRNAs, which are putatively involved in the regulation of ricinoleic acid biosynthesis. The transcriptome data were further validated by analyzing the expression profiles of key enzyme-encoding genes, TFs and lncRNAs with quantitative real-time PCR. Finally, a network of genes associated with ricinoleic acid accumulation in H. benghalensis was established.

CONCLUSIONS

This study was the first step toward the understating of the molecular mechanisms of ricinoleic acid biosynthesis and oil accumulation in H. benghalensis seeds and identified a pool of novel genes regulating ricinoleic acid accumulation. The results set a foundation for developing H. benghalensis into a novel ricinoleic acid feedstock at the transcriptomic level and provided valuable candidate genes for improving ricinoleic acid production in other plants.

Understanding microRNA Regulation Involved in the Metamorphosis of the Veined Rapa Whelk (Rapana venosa)

The veined rapa whelk (Rapana venosa) is widely consumed in China. Nevertheless, it preys on oceanic bivalves, thereby reducing this resource worldwide. Its larval metamorphosis comprises a transition from pelagic to benthic form, which involves considerable physiological and structural changes and has vital roles in its natural populations and commercial breeding. Thus, understanding the endogenous microRNAs (miRNAs) that drive metamorphosis is of great interest. This is the first study to use high-throughput sequencing to examine the alterations in miRNA expression that occur during metamorphosis in a marine gastropod. A total of 195 differentially expressed miRNAs were obtained. Sixty-five of these were expressed during the transition from precompetent to competent larvae. Thirty-three of these were upregulated and the others were downregulated. Another 123 miRNAs were expressed during the transition from competent to postlarvae. Ninety-six of these were upregulated and the remaining 27 were downregulated. The expression of miR-276-y, miR-100-x, miR-183-x, and miR-263-x showed a >100-fold change during development, while the miR-242-x and novel-m0052-3p expression levels changed over 3000-fold. Putative target gene coexpression, gene ontology, and pathway analyses suggest that these miRNAs play important parts in cell proliferation, migration, apoptosis, metabolic regulation, and energy absorption. Twenty miRNAs and their target genes involved in ingestion, digestion, cytoskeleton, cell adhesion, and apoptosis were identified. Nine of them were analyzed with real-time polymerase chain reaction (PCR), which showed an inverse correlation between the miRNAs and their relative expression levels. Our data elucidate the role of miRNAs in R. venosa metamorphic transition and serve as a solid basis for further investigations into regulatory mechanisms of gastropod metamorphosis.

Selection of housekeeping genes as internal controls for quantitative RT-PCR analysis of the veined rapa whelk (Rapana venosa)

Background

The veined rapa whelk Rapana venosa is an important commercial shellfish in China and quantitative real-time PCR (qRT-PCR) has become the standard method to study gene expression in R. venosa. For accurate and reliable gene expression results, qRT-PCR assays require housekeeping genes as internal controls, which display highly uniform expression in different tissues or stages of development. However, to date no studies have validated housekeeping genes in R. venosa for use as internal controls for qRT-PCR.

Methods

In this study, we selected the following 13 candidate genes for suitability as internal controls: elongation factor-1α (EF-1α), α-actin (ACT), cytochrome c oxidase subunit 1 (COX1), nicotinamide adenine dinucleotide dehydrogenase (ubiquinone) 1α subcomplex subunit 7 (NDUFA7), 60S ribosomal protein L5 (RL5), 60S ribosomal protein L28 (RL28), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), β-tubulin (TUBB), 40S ribosomal protein S25 (RS25), 40S ribosomal protein S8 (RS8), ubiquitin-conjugating enzyme E2 (UBE2), histone H3 (HH3), and peptidyl-prolyl cis-trans isomerase A (PPIA). We measured the expression levels of these 13 candidate internal controls in eight different tissues and twelve larvae developmental stages by qRT-PCR. Further analysis of the expression stability of the tested genes was performed using GeNorm and RefFinder algorithms.

Results

Of the 13 candidate genes tested, we found that EF-1α was the most stable internal control gene in almost all adult tissue samples investigated with RL5 and RL28 as secondary choices. For the normalization of a single specific tissue, we suggested that EF-1α and NDUFA7 are the best combination in gonad, as well as COX1 and RL28 for intestine, EF-1α and RL5 for kidney, EF-1α and COX1 for gill, EF-1α and RL28 for Leiblein and mantle, EF-1α, RL5, and NDUFA7 for liver, GAPDH, PPIA, and RL28 for hemocyte. From a developmental perspective, we found that RL28 was the most stable gene in all developmental stages measured, and COX1 and RL5 were appropriate secondary choices. For the specific developmental stage, we recommended the following combination for normalization, PPIA, RS25, and RL28 for stage 1, RL5 and RL28 for stage 2 and 5, RL28 and NDUFA7 for stage 3, and PPIA and TUBB for stage 4.

Discussion

Our results are instrumental for the selection of appropriately validated housekeeping genes for use as internal controls for gene expression studies in adult tissues or larval development of R. venosa in the future.