Validation of reference genes for real-time PCR of reproductive system in the black tiger shrimp

Reference genes for qPCR expression in black tiger shrimp, Penaeus monodon

BACKGROUND

Gene expression profiling via qPCR is an essential tool for unraveling the intricate molecular mechanisms underlying growth and development. Identifying and validating the most appropriate reference genes is essential for qPCR experiments. Nevertheless, there exists a deficiency in a thorough assessment of reference genes concerning the expression of the genes in the research in the context of the growth and development of the Black Tiger Shrimp, P. monodon. This popular marine crustacean is extensively raised for human consumption. In this study, we assessed the expression stability of seven reference genes (ACTB, 18S, EF-1α, AK, PK, cox1, and CLTC) in adult tissues (hepatopancreas, gills, and stomach) of small and large polymorphs of P. monodon.

METHODS AND RESULTS

The stability of gene expressions was assessed utilizing NormFinder, BestKeeper, and geNorm, and a comprehensive ranking of these genes was conducted through the online tool RefFinder. In the overall ranking, 18S and CLTC emerged as the most stable genes in the hepatopancreas and stomach, while CLTC and AK exhibited significant statistical reliability in the gills of adult P. monodon. The validation of these identified stable genes was carried out using a growth-associated gene, insr-1.

CONCLUSION

The results indicated that 18S and CLTC stand out as the most versatile reference genes for conducting qPCR analysis focused on the growth of P. monodon. This study represents the first comprehensive exploration that identifies and assesses reference genes for qPCR analysis in P. monodon, providing valuable tools for research involving similar crustaceans.

Screening of Suitable Reference Genes for Immune Gene Expression Analysis Stimulated by Vibrio anguillarum and Copper Ions in Chinese Mitten Crab (Eriocheir sinensis)

The reference gene expression is not always stable under different experimental conditions, and screening of suitable reference genes is a prerequisite in quantitative real-time polymerase chain reaction (qRT-PCR). In this study, we investigated gene selection, and the most stable reference gene for the Chinese mitten crab (Eriocheir sinensis) was screened under the stimulation of Vibrio anguillarum and copper ions, respectively. Ten candidate reference genes were selected, including arginine kinase (AK), ubiquitin-conjugating enzyme E2b (UBE), glutathione S-transferase (GST), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), elongation factor 1α (EF-1α), α-tubulin (α-TUB), heat shock protein 90 (HSP90), β-actin (β-ACTIN), elongation factor 2 (EF-2) and phosphoglucomutase 2 (PGM2). Expression levels of these reference genes were detected under the stimulation of V. anguillarum at different times (0 h, 6 h, 12 h, 24 h, 48 h and 72 h) and copper ions in different concentrations (11.08 mg/L, 2.77 mg/L, 0.69 mg/L and 0.17 mg/L). Four types of analytical software, namely geNorm, BestKeeper, NormFinder and Ref-Finder, were applied to evaluate the reference gene stability. The results showed that the stability of the 10 candidate reference genes was in the following order: AK > EF-1α > α-TUB > GAPDH > UBE > β-ACTIN > EF-2 > PGM2 > GST > HSP90 under V. anguillarum stimulation. It was GAPDH > β-ACTIN > α-TUB > PGM2 > EF-1α > EF-2 > AK > GST > UBE > HSP90 under copper ion stimulation. The expression of E. sinensis Peroxiredoxin4 (EsPrx4) was detected when the most stable and least stable internal reference genes were selected, respectively. The results showed that reference genes with different stability had great influence on the accurate results of the target gene expression. In the Chinese mitten crab (E. sinensis), AK and EF-1α were the most suitable reference genes under the stimulation of V. anguillarum. Under the stimulation of copper ions, GAPDH and β-ACTIN were the most suitable reference genes. This study provided important information for further research on immune genes in V. anguillarum or copper ion stimulation.

Overexpression of LAS21 in Cellulase-Displaying Saccharomyces cerevisiae for High-Yield Ethanol Production from Pretreated Sugarcane Bagasse

The valorization of lignocellulosic feedstocks into biofuels and biochemicals has received much attention due to its environmental friendliness and sustainability. However, engineering an ideal microorganism that can both produce sufficient cellulases and ferment ethanol is highly challenging. In this study, we have tested seven different genes that are involved in glycosylphosphatidylinositol (GPI) biosynthesis and remodeling for the improvement of cellulase activity tethered on the S. cerevisiae cell surface. It was found that the overexpression of LAS21 can improve β-glucosidase activity by 48.8% compared to the original strain. Then, the three cellulase genes (cellobiohydrolase, endoglucanase, and β-glucosidase) and the LAS21 gene were co-introduced into a diploid thermotolerant S. cerevisiae strain by a multiple-round transformation approach, resulting in the cellulolytic ECBLCCE5 strain. Further optimization of the bioprocess parameters was found to enhance the ethanol yield of the ECBLCCE5 strain. Scaling up the valorization of pretreated sugarcane bagasses in a 1 L bioreactor resulted in a maximum ethanol concentration of 28.0 g/L (86.5% of theoretical yield). Our study provides a promising way to improve the economic viability of second-generation ethanol production. Moreover, the engineering of genes involved in GPI biosynthesis and remodeling can be applied to other yeast cell surface display applications.

Molecular characterization and expression profiling of transformer 2 and fruitless-like homologs in the black tiger shrimp, Penaeus monodon

Transformer 2 (tra 2) and fruitless (fru) genes have been proven to play a key role in sex determination pathways in many Arthropods, including insects and crustaceans. In this study, a paralog of P. monodon tra 2 (Pmtra 2), P. monodon ovarian associated transformer 2 (PmOvtra 2) and 2 isoforms of P. monodon fruitless-like gene (Pmfru-1 and Pmfru-2) were identified and characterized. The full cDNA sequence of PmOvtra 2 consisted of 1,774 bp with the longest open reading frame (ORF) of 744 bp encoding for 247 amino acids. The PmOvtra 2 exhibited a predicted RNA-recognition motif (RRM) domain and two arginine-serine (RS) regions, suggesting its function in RNA splicing. The full cDNA sequence of Pmfru-1 consisted of 1,306 bp with 1,182 bp ORF encoding for 393 amino acids, whereas the full cDNA sequence of Pmfru-2 consisted of 1,858 bp with 1,437 bp ORF encoding 478 amino acids. The deduced amino acid sequences of Pmfru-1 and Pmfru-2 exhibited highly conserved domains of Fru proteins, including Broad-complex, Tramtrack and Bric-a-brac (BTB), and zinc finger (ZF) domains. In addition, Pmfru-1 and Pmfru-2 were suggestively originated from the same single genomic locus by genomic sequence analysis. Specifically, Pmfru pre-mRNA was alternatively spliced for Pmfru-1 and Pmfru-2 to include mutually exclusive exon 7 and exon 6, respectively. Temporal and spatial expression of PmOvtra 2, Pmfru-1, and Pmfru-2 were also investigated by qPCR. The results showed that all were expressed in early developmental stages with undifferentiated gonads starting from nauplius until postlarvae. The expression of PmOvtra 2 started at nauplius stage and gradually increased from mysis to postlarvae (PL) 1. However, the expression of Pmfru-1 was low at the nauplii stage and slightly increased from protozoea to PL5, whereas the expression of Pmfru-2 maintained a low level from nauplius to mysis and then gradually increased at the PL stages. Expressions of PmOvtra 2, Pmfru-1, and Pmfru-2 were detected in various tissues including nervous tissue, gill, heart, hepatopancreas, gut, and gonads. Interestingly, the sexually dimorphic expression of PmOvtra 2, Pmfru-1, and Pmfru-2 was demonstrated in fully developed gonads in which the ovary showed significantly higher expressions than the testis. The great difference in the expression pattern of PmOvtra 2, Pmfru-1, and Pmfru-2 in the ovary and testis suggested their roles in the female sex determination in P. monodon.

Validation of reference genes for quantitative real-time PCR in chemical exposed and at different age's brackish water flea Diaphanosoma celebensis

Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), a primary approach for evaluating gene expression, requires an appropriate normalization strategy to confirm relative gene expression levels by comparison, and rule out variations that might occur in analytical procedures. The best option is to use a reference gene whose expression level is stable across various experimental conditions to compare the mRNA levels of a target gene. However, there is limited information on how the reference gene is differentially expressed at different ages (growth) in small invertebrates with notable changes such as molting. In this study, expression profiles of nine candidate reference genes from the brackish water flea, Diaphanosoma celebensis, were evaluated under diverse exposure to toxicants and according to growth. As a result, four different algorithms showed similar stabilities of genes for chemical exposures in the case of limited conditions using the same developmental stage (H2A was stable, whereas Act was fairly unstable in adults), while the results according to age showed a significantly different pattern in suite of candidate reference genes. This affected the results of genes EcRA and GST, which are involved in development and detoxification mechanisms, respectively. Our finding is the first step towards establishing a standardized real-time qRT-PCR analysis of this environmentally important invertebrate that has potential for aquatic ecotoxicology, particularly in estuarine environments.

Long non-coding RNA profile in banana shrimp, Fenneropenaeus merguiensis and the potential role of lncPV13 in vitellogenesis

The long non-coding RNAs (lncRNAs) have been known to play important roles in several biological processes as well as in reproduction. This study aimed to identify lncRNA in ovary female banana shrimp, Fenneropenaeus merguiensis, and investigate the potential role of lncPV13 in the vitellogenesis. After the in silico identification of the ovarian transcriptome, a total of 24,733 putative lncRNAs were obtained, and only 147 putative lncRNAs were significantly differentially expressed among the ovarian development stages. To validate the in silico identification of lncRNAs, the 16 lncRNAs with the highest differential expression in the transcriptome analysis were evaluated by RT-qPCR. The 6 lncRNAs showed higher expression levels in the mature stage than in the previtellogenic stage and were found in several tissues such as in eyestalks, brains, thoracic ganglia, gills, and muscle. Furthermore, most candidate lncRNAs were amplifiable in Litopenaeus vannamei's and Penaeus monodon's DNA but not in Macrobrachium rosenbergii's DNA, suggesting some lncRNAs are expressed in a species-specific manner among penaeid shrimp. In this study, the lncPV13 was investigated for its vitellogenin regulating function by RNA interference. The result indicates that the lncPV13 expression was suppressed in the ovary on day 7 after the injection of double-stranded RNA specific to lncPV13 (dslncPV13), while vitellogenin (Vg) expression was significantly decreased. In contrast, the gonad inhibiting hormone (GIH) expression was significantly increased in the lncPV13 knockdown shrimp. However, the oocyte proliferation was not significantly different between control and lncPV13 knockdown shrimp. This suggests that lncPV13 regulate Vg synthesis through GIH inhibition. Finally, our findings provide lncRNA information and potential lncRNAs involved in the vitellogenesis of female banana shrimp.

Reference gene and tropomyosin expression in mud crab Scylla olivacea, Scylla paramamosain and Scylla tranquebarica

Tropomyosin, a muscle tissue protein is a major allergen in most of shellfish including mud crab. Quantitative real time-PCR (qRT-PCR) using a stable reference gene is the most sensitive approach to produce accurate relative gene expression that has yet to be demonstrated for allergenic tropomyosin in mud crab species. This study was conducted to identify the suitable reference gene and tropomyosin expression in different body parts of local mud crabs, Scylla olivacea, Scylla paramamosain and Scylla tranquebarica. Myosin, 18S rRNA, GADPH and EF1α were selected as candidate reference genes and their expression was measured in the abdomen, walking leg and cheliped tissues of local Scylla spp. The expression stability was analyzed using the comparative delta-Ct method, BestKeeper, NormFinder and geNorm then comprehensively ranked by RefFinder algorithm. Findings showed that EF1α was the most suitable reference gene across three mud crab species. Meanwhile, the abdomen, walking leg and cheliped selected their own suitable reference gene either Myosin, 18S rRNA, EF1α or GADPH. Overall, tropomyosin was the highest in S. tranquebarica, whereas the least was in S. paramamosain. Interestingly, tropomyosin was the highest in the abdomen of all mud crab species. This is the first analysis on reference genes selection for qRT-PCR data normalization of tropomyosin expression in mud crab. These results will provide more accurate findings for further gene expression and allergen analysis in Scylla spp.

Validation of Real-time PCR Reference Genes of Muscle Metabolism in Harvested Spiny-Cheek Crayfish (Faxonius limosus) Exposed to Seasonal Variation

Real-time quantitative reverse transcription PCR (RT-qPCR) is a sensitive and broadly used technique of assessing gene activity. To obtain a reliable result, stably expressed reference genes are essential for normalization of transcripts in various samples. To our knowledge, this is the first systematic analysis of reference genes for normalization of RT-qPCR data in spiny-cheek crayfish (Faxonius limosus). In this study, expression of five candidate reference genes (actb, β-actin; gapdh, glyceraldehyde-3-phosphate dehydrogenase; eif, eukaryotic translation initiation factor 5a; ef-1α, elongation factor-1α; and tub, α-tubulin) in muscle samples from male and female F. limosus in spring and autumn was analyzed. Additionally, the most stable reference genes were used for accurate normalization of five target genes, i.e., tnnc, troponin c; ak, arginine kinase; fr, ferritin; ccbp-23, crustacean calcium-binding protein 23; and actinsk8, skeletal muscle actin 8. Results obtained using the geNorm and NormFinder algorithms showed high consistency, and differences in the activity of the selected actb with eif genes were successfully identified. The spring and autumn activities of the target genes (except ak) in the muscle tissue of males and females differed significantly, showing that both sexes are immensely involved in an array of breeding behaviors in spring, and females intensively recover in the autumn season. Characterization of first reference genes in spiny-cheek crayfish will facilitate more accurate and reliable expression studies in this key species.

Transcriptome analyses reveal the synergistic effects of feeding and eyestalk ablation on ovarian maturation in black tiger shrimp

Unilateral eyestalk ablation in the female black tiger shrimp Penaeus monodon is commonly employed to induce ovarian maturation. However, the importance of complementing this practice with the provision of live feed supplement (such as polychaetes) has not been emphasized in previous studies. Indeed, it has been less emphasized that female broodstock must be fed with live feeds such as polychaetes for this practice to be effective. While the effects of eyestalk ablation have been thoroughly studied in various aspects, the synergistic effects of feeding with live feeds and the ablation have never been elucidated at a transcriptome-wide level. With recent advances in the next-generation sequencing platforms, it is now possible to investigate the effects of eyestalk ablation and live feeds at the transcriptomic levels. This study employed both short-read Illumina RNA sequencing and long-read Pacific Biosciences (PacBio) isoform sequencing (Iso-seq) to generate the first high-quality ovarian reference transcriptome in P. monodon. This novel assembly allowed us to dissect the effects of feeds and eyestalk ablation and reveal their synergistic effects at the transcriptomic level through the regulation of important genes involved in fatty acid regulation, energy production, and hormone-mediated oocyte maturation pathways. The synergistic effects between the polychaete feeding and the eyestalk ablation in the process of ovarian maturation in black tiger shrimp suggest that without having proper nutrients from the polychaetes, female broodstock might not be ready to develop its ovary. However, even with proper nutrients, the eyestalk ablation is still necessary to perhaps manipulate the female endocrine of the black tiger shrimp. These findings shed the light on molecular mechanisms and key molecular pathways that lead to successful ovarian maturation.

Variation of protein kinase C-α expression in eyestalk removal-activated ovaries in whiteleg shrimp, Litopenaeus vannamei

The regulation of female reproduction in crustaceans is controlled by a variety of hormones. Previous studies of hormone-initiated cellular mechanisms controlling ovarian maturation focused mainly on those initiated by inhibitory hormones. In order to facilitate research on ovarian development, it is necessary to gain a better understanding of factors that promote ovarian growth on the cellular level. Here, we used eyestalk ablation to firstly induce a state of ovarian maturation in Litopenaeus vannamei. Gonadosomatic index, hemolymph vitellogenin (Vg) concentrations, and Vg gene transcript levels in the ovaries were significantly elevated at 1 and 2 weeks after eyestalk ablation (P < 0.05). Correspondingly, immunoblot analysis revealed a remarkable decrease in anti-PKC-α immunoreactivity in both cytosol and membrane fractions of ovarian tissue homogenates: it was strongly apparent in intact animals, but decreased with time after eyestalk ablation, showing a stronger tendency to do so in the membrane fraction than in the cytosol fraction. Considered overall, the data presented strongly suggest that PKC-α isoform plays a role in the regulation of ovarian growth in L. vannamei through a negative-based regulating mechanism.

Effects of 5-Aminolevulinic Acid on Gene Expression, Immunity, and ATP Levels in Pacific White Shrimp, Litopenaeus vannamei

With the emergence of several infectious diseases in shrimp aquaculture, there is a growing interest in the use of feed additives to enhance shrimp immunity. Recently, the use of 5-aminolevulinic acid (5-ALA), a non-protein amino acid that plays a rate-limiting role in heme biosynthesis, has received attention for its positive effect on immunity in livestock animals. To evaluate the effect of 5-ALA in the Pacific white shrimp, Litopenaeus vannamei, we conducted microarray analysis, a Vibrio parahaemolyticus immersion challenge test, an ATP level assay, and gene expression analysis of some hemoproteins and genes associated with heme synthesis and degradation. Out of 15,745 L. vannamei putative genes on the microarray, 101 genes were differentially expressed by more than fourfold (p < 0.05) between 5-ALA-supplemented and control shrimp hepatopancreas. 5-ALA upregulated 99 of the 101 genes, 41 of which were immune- and defense-related genes based on sequence homology. Compared to the control, the 5-ALA-supplemented group had a higher survival rate in the challenge test, higher transcript levels of porphobilinogen synthase, ferrochelatase, catalase, nuclear receptor E75, and heme oxygenase-1 and higher levels of ATP. These findings suggest that dietary 5-ALA enhanced the immune response of L. vannamei to V. parahaemolyticus, upregulated immune- and defense-related genes, and enhanced aerobic energy metabolism, respectively. Further studies are needed to elucidate the extent of 5-ALA use in shrimp culture.

Validation and Evaluation of Reference Genes for Quantitative Real-Time PCR in Macrobrachium Nipponense

Quantitative real-time PCR (qPCR) is widely used in molecular biology, although the accuracy of the quantitative results is determined by the stability of the reference genes used. Recent studies have investigated suitable reference genes for some crustaceans under various conditions, but studies in Macrobrachium nipponense are currently lacking. In this study, we selected the following seven genes from among 35 commonly used housekeeping genes as candidate qPCR reference genes for temporal and spatial expression: EIF (eukaryotic translation initiation factor 5A), 18S (18S ribosomal RNA), EF-1α (elongation factor-1α), GAPDH (glyceraldehyde-3-phosphate dehydrogenase), TUB (α-tubulin), β-act (β-actin), and RPL18 (Ribosomal protein L18). The stability of each reference gene was evaluated by GeNorm, NormFinder, BestKeeper, and comparative ∆C _t methods, and was comprehensively ranked using RefFinder. RPL18 was shown to be the most suitable reference gene for adult M. nipponense tissues, while EIF was the most stable in different ovarian and embryo stages and in white spot syndrome virus infection, and β-act was the most stable reference gene under hypoxia stress. The reliability of the rankings was confirmed by RNA interference experiments. To the best of our knowledge, this represents the first systematic analysis of reference genes for qPCR experiments in M. nipponense, and the results will provide invaluable information for future research in closely related crustaceans.

Gene expression of benthic amphipods (genus: Diporeia) in relation to a circular ssDNA virus across two Laurentian Great Lakes

Circular rep-encoding ssDNA (CRESS-DNA) viruses are common constituents of invertebrate viral consortia. Despite their ubiquity and sequence diversity, the effects of CRESS-DNA viruses on invertebrate biology and ecology remain largely unknown. This study assessed the relationship between the transcriptional profile of benthic amphipods of genus Diporeia and the presence of the CRESS-DNA virus, LM29173, in the Laurentian Great Lakes to provide potential insight into the influence of these viruses on invertebrate gene expression. Twelve transcriptomes derived from Diporeia were compared, representing organisms from two amphipod haplotype clades (Great Lakes Michigan and Superior, defined by COI barcode sequencing) with varying viral loads (up to 3 × 10⁶ genome copies organism⁻¹). Read recruitment to de novo assembled transcripts revealed 2,208 significantly over or underexpressed contigs in transcriptomes with above average LM29173 load. Of these contigs, 31.5% were assigned a putative function. The greatest proportion of annotated, differentially expressed transcripts were associated with functions including: (1) replication, recombination, and repair, (2) cell structure/biogenesis, and (3) post-translational modification, protein turnover, and chaperones. Contigs putatively associated with innate immunity displayed no consistent pattern of expression, though several transcripts were significantly overexpressed in amphipods with high viral load. Quantitation (RT-qPCR) of target transcripts, non-muscular myosin heavy chain, β-actin, and ubiquitin-conjugating enzyme E2, corroborated transcriptome analysis and indicated that Lake Michigan and Lake Superior amphipods with high LM29173 load exhibit lake-specific trends in gene expression. While this investigation provides the first comparative survey of the transcriptional profile of invertebrates of variable CRESS-DNA viral load, additional inquiry is required to define the scope of host-specific responses to potential infection.

Sex and tissue specific gene expression patterns identified following de novo transcriptomic analysis of the Norway lobster, Nephrops norvegicus

BACKGROUND

The Norway lobster, Nephrops norvegicus, is economically important in European fisheries and is a key organism in local marine ecosystems. Despite multi-faceted scientific interest in this species, our current knowledge of genetic resources in this species remains very limited. Here, we generated a reference de novo transcriptome for N. norvegicus from multiple tissues in both sexes. Bioinformatic analyses were conducted to detect transcripts that were expressed exclusively in either males or females. Patterns were validated via RT-PCR.

RESULTS

Sixteen N. norvegicus libraries were sequenced from immature and mature ovary, testis and vas deferens (including the masculinizing androgenic gland). In addition, eyestalk, brain, thoracic ganglia and hepatopancreas tissues were screened in males and both immature and mature females. RNA-Sequencing resulted in >600 million reads. De novo assembly that combined the current dataset with two previously published libraries from eyestalk tissue, yielded a reference transcriptome of 333,225 transcripts with an average size of 708 base pairs (bp), with an N50 of 1272 bp. Sex-specific transcripts were detected primarily in gonads followed by hepatopancreas, brain, thoracic ganglia, and eyestalk, respectively. Candidate transcripts that were expressed exclusively either in males or females were highlighted and the 10 most abundant ones were validated via RT-PCR. Among the most highly expressed genes were Serine threonine protein kinase in testis and Vitellogenin in female hepatopancreas. These results align closely with gene annotation results. Moreover, a differential expression heatmap showed that the majority of differentially expressed transcripts were identified in gonad and eyestalk tissues. Results indicate that sex-specific gene expression patterns in Norway lobster are controlled by differences in gene regulation pattern between males and females in somatic tissues.

CONCLUSIONS

The current study presents the first multi-tissue reference transcriptome for the Norway lobster that can be applied to future biological, wild restocking and fisheries studies. Sex-specific markers were mainly expressed in males implying that males may experience stronger selection than females. It is apparent that differential expression is due to sex-specific gene regulatory pathways that are present in somatic tissues and not from effects of genes located on heterogametic sex chromosomes. The N. norvegicus data provide a foundation for future gene-based reproductive studies.

Differential regulation of the lipoxygenase pathway in shrimp hepatopancreases and ovaries during ovarian development in the black tiger shrimp Penaeus monodon

Dietary polyunsaturated fatty acids (PUFAs) are critical to the success of ovarian development in marine crustaceans, especially for domesticated species such as the black tiger shrimp Penaeus monodon. These fatty acids are stored in a midgut gland called the hepatopancreas and subsequently serve as an energy source or are incorporated in yolk during ovarian development. PUFAs are known precursors of hydroxy fatty acids, including hydroxyeicosatetraenoic acid and hydroxyeicosapentaenoic acid (HEPE), which are catalyzed by lipoxygenases (LOX). In previous studies, 8-HEPE has been shown to regulate female reproduction and adipogenesis in marine crustaceans. However, whether the biosynthesis of 8-HEPE in these species is the result of LOX activity has yet to be investigated. In this study, 8-HEPE was identified exclusively in P. monodon hepatopancreases using liquid chromatography-mass spectrometry. Treatment with nordihydroguaiaretic acid resulted in the reduction of 8-HEPE, suggesting the enzyme-dependent catalysis of 8-HEPE in hepatopancreases. Additionally, a full-length P. monodon LOX (PmLOX) was amplified from shrimp ovary cDNA. Sequence analysis revealed that the putative PmLOX contains domains and catalytic residues required for LOX catalytic function. Furthermore, PmLOX expression increased steadily as shrimp ovary maturation progressed, while PmLOX expression and the amount of 8-HEPE decreased in shrimp hepatopancreases. These findings not only suggest differential requirements for hydroxy fatty acid biosynthesis in shrimp ovaries and hepatopancreases during the P. monodon ovarian development, but also provide insights into the LOX pathway in marine crustaceans.

Transcriptomic Analysis of Male Black Tiger Shrimp (Penaeus monodon) After Polychaete Feeding to Enhance Testicular Maturation

To reveal molecular mechanism of how polychaetes enhanced reproductive maturation in the male black tiger shrimp (Penaeus monodon), transcriptomic profiles of male reproductive organs (testes and vas deferens) between polychaete-fed and commercial pellet-fed male brooders were compared using cDNA microarray. The overall profiles were distinguishingly different between the two feed groups as well as between testes and vas deferens. Additionally, six of 11 differentially expressed gene identified by the microarray (HNRPUL1 and GCP4 in testes, MAT2B, CDC16, and CSN5 in vas deferens, and SLD5 in both organs) were validated by quantitative real-time PCR (qPCR) and found to exhibit significantly higher expression levels in polychaete-fed shrimp than those in commercial pellet-fed shrimp. From microarray and qPCR results, the differentially expressed transcripts in both testes and vas deferens between different feeds belonged to DNA replication and microtubule nucleation pathways. Interestingly, while the transcripts involved in nutrient uptake and nucleotide biosynthesis were increased only in testes, those involved in protein refolding and apoptosis were increased only in vas deferens. These findings suggest that polychaetes may enhance spermatogenesis by increasing spermatogonia proliferation in testes and by regulating mature spermatozoa in vas deferens.

Analysis of variance, normal quantile-quantile correlation and effective expression support of pooled expression ratio of reference genes for defining expression stability

Identification of a reference gene unaffected by the experimental conditions is obligatory for accurate measurement of gene expression through relative quantification. Most existing methods directly analyze variability in crossing point (Cp) values of reference genes and fail to account for template-independent factors that affect Cp values in their estimates. We describe the use of three simple statistical methods namely analysis of variance (ANOVA), normal quantile-quantile correlation (NQQC) and effective expression support (EES), on pooled expression ratios of reference genes in a panel to overcome this issue. The pooling of expression ratios across the genes in the panel nullify the sample specific effects uniformly affecting all genes that are falsely reflected as instability. Our methods also offer the flexibility to include sample specific PCR efficiencies in estimations, when available, for improved accuracy. Additionally, we describe a correction factor from the ANOVA method to correct the relative fold change of a target gene if no truly stable reference gene could be found in the analyzed panel. The analysis is described on a synthetic data set to simplify the explanation of the statistical treatment of data.

Identification and evaluation of reference genes for expression studies by RT-qPCR during embryonic development of the emerging model organism, Macrobrachium olfersii

RT-qPCR is a sensitive and highly efficient technique that is widely used in gene expression analysis and to provide insight into the molecular mechanisms underlying embryonic development. The freshwater prawn, Macrobrachium olfersii is an emerging model organism, but, the stable reference genes of this species need to be identified and validated for RT-qPCR analysis. Thus, the aim of this study was to evaluate the expression stability of six genes (β-act, GAPDH, EF-1α, RpL8, RpS6, AK) in embryos and in adult tissues (cerebral ganglia, muscle and hepatopancreas) of M. olfersii. The expression stabilities of these genes were evaluated using geNorm, NormFinder, BestKeeper, ΔCt method and integrated tool RefFinder. In the general ranking, RpL8 and RpS6 were the most stable genes in embryos, while RpS6 and RpL8 were the most stable in a combined adult tissue analysis. Analysis of the adult tissues revealed that β-act and AK were the most stable genes in cerebral ganglia, RpL8 and AK in muscle, and RpS6 and β-act in hepatopancreas. EF-1α and GAPDH were the least stable genes and as normalizer genes in RT-qPCR affected expression of the Distal-less gene during M. olfersii development. This study provides suitable reference genes for RT-qPCR analysis and allows future studies of the gene expression in M. olfersii for understanding the molecular mechanisms of their development. To our knowledge, this is the first published study that identifies and evaluates reference genes for RT-qPCR analysis in M. olfersii and could be useful as basis for evaluations of reference genes in other prawns.

Neuroendocrine disruption in the shore crab Carcinus maenas: Effects of serotonin and fluoxetine on chh- and mih-gene expression, glycaemia and ecdysteroid levels

Serotonin, a highly conserved neurotransmitter, controls many biological functions in vertebrates, but also in invertebrates. Selective serotonin reuptake inhibitors (SSRIs), such as fluoxetine, are commonly used in human medication to ease depression by affecting serotonin levels. Their residues and metabolites can be detected in the aquatic environment and its biota. They may also alter serotonin levels in aquatic invertebrates, thereby perturbing physiological functions. To investigate whether such perturbations can indeed be expected, shore crabs (Carcinus maenas) were injected either with serotonin, fluoxetine or a combination of both. Dose-dependent effects of fluoxetine ranging from 250 to 750nM were investigated. Gene expression of crustacean hyperglycemic hormone (chh) as well as moult inhibiting hormone (mih) was assessed by RT-qPCR at 2h and 12h after injection. Glucose and ecdysteroid levels in the haemolymph were monitored in regular intervals until 12h. Serotonin led to a rapid increase of chh and mih expression. On the contrary, fluoxetine only affected chh and mih expression after several hours, but kept expression levels significantly elevated. Correspondingly, serotonin rapidly increased glycaemia, which returned to normal or below normal levels after 12h. Fluoxetine, however, resulted in a persistent low-level increase of glycaemia, notably during the period when negative feedback regulation reduced glycaemia in the serotonin treated animals. Ecdysteroid levels were significantly decreased by serotonin and fluoxetine, with the latter showing less pronounced and less rapid, but longer lasting effects. Impacts of fluoxetine on glycaemia and ecdysteroids were mostly observed at higher doses (500 and 750nM) and affected principally the response dynamics, but not the amplitude of glycaemia and ecdysteroid-levels. These results suggest that psychoactive drugs are able to disrupt neuroendocrine control in decapod crustaceans, as they interfere with the normal regulation of the serotonergic system.

Silencing of Gonad-Inhibiting Hormone Transcripts in Litopenaeus vannamei Females by use of the RNA Interference Technology

The method usually employed to stimulate gonadal maturation and spawning of captive shrimp involves unilateral eyestalk ablation, which results in the removal of the endocrine complex responsible for gonad-inhibiting hormone (GIH) synthesis and release. In the present study, RNAi technology was used to inhibit transcripts of GIH in Litopenaeus vannamei females. The effect of gene silencing on gonad development was assessed by analyzing the expression of GIH and vitellogenin, respectively, in the eyestalk and ovaries of L. vannamei females, following ablation or injection with dsRNA-GIH, dsRNA-IGSF4D (non-related dsRNA), or saline solution. Histological analyses were performed to determine the stage of gonadal development and to assess the diameter of oocytes throughout the experimental procedure. Only oocytes at pre-vitellogenesis and primary vitellogenesis stages were identified in females injected with dsRNA-GIH, dsRNA-IGSF4D, or saline solution. Oocytes at all developmental stages were observed in eyestalk-ablated females, with predominance of later stages, such as secondary vitellogenesis and mature oocytes. Despite achieving 64, 73, and 71% knockdown of eyestalk GIH mRNA levels by 15, 30, and 37 days post-injection (dpi), respectively, in dsRNA-GIH-injected females, the expected increase in ovary vitellogenin mRNA expression was only observed on the 37th dpi. This is the first report of the use of RNAi technology to develop an alternative method to eyestalk ablation in captive L. vannamei shrimps.

Gonadal transcriptomic analysis and differentially expressed genes in the testis and ovary of the Pacific white shrimp (Litopenaeus vannamei)

Background

The Pacific white shrimp (Litopenaeus vannamei) is the world’s most prevalent cultured crustacean species. However, the supply of high-quality broodstocks is limited and baseline information related to its reproductive activity and molecular issues related to gonad development are scarce. In this study, we performed transcriptome sequencing on the gonads of adult male and female L. vannamei to identify sex-related genes.

Results

A total of 25.16 gigabases (Gb) of sequences were generated from four L. vannamei gonadal tissue libraries. After quality control, 24.11 Gb of clean reads were selected from the gonadal libraries. De-novo assembly of all the clean reads generated a total of 65,218 unigenes with a mean size of 1021 bp and a N50 of 2000 bp. A search of all-unigene against Nr, SwissProt, KEGG, COG and NT databases resulted in 26,482, 23,062, 20,659, 11,935 and 14,626 annotations, respectively, providing a total of 30,304 annotated unigenes. Among annotated unigenes, 12,320 unigenes were assigned to gene ontology categories and 20,659 unigenes were mapped to 258 KEGG pathways. By comparing the ovary and testis libraries, 19,279 testicular up-regulated and 3,529 ovarian up-regulated unigenes were identified. Enrichment analysis of differentially expressed unigenes resulted in 1060 significantly enriched GO terms and 34 significantly enriched KEGG pathways. Nine ovary-specific, 6 testis-specific, 45 testicular up-regulated and 39 ovarian up-regulated unigenes were then confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, using all-unigenes as a reference, a total of 13,233 simple sequence repeats (SSRs) were identified in 10,411 unigene sequences.

Conclusions

The present study depicts the first large-scale RNA sequencing of shrimp gonads. We have identified many important sex-related functional genes, GO terms and pathways, all of which will facilitate future research into the reproductive biology of shrimp. We expect that the SSRs detected in this study can then be used as genetic markers for germplasm evaluation of breeding and imported populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2219-4) contains supplementary material, which is available to authorized users.

Identification and Characterization of Reference Genes for Normalizing Expression Data from Red Swamp Crawfish Procambarus clarkii

qRT-PCR is a widely used technique for rapid and accurate quantification of gene expression data. The use of reference genes for normalization of the expression levels is crucial for accuracy. Several studies have shown that there is no perfect reference gene that is appropriate for use in all experimental conditions, and research on suitable reference genes in red swamp crawfish (Procambarus clarkii) is particularly scarce. In this study, eight commonly used crustacean reference genes were chosen from P. clarkii transcriptome data and investigated as potential candidates for normalization of qRT-PCR data. Expression of these genes under different experimental conditions was examined by qRT-PCR, and the stability of their expression was evaluated using three commonly used statistical algorithms, geNorm, NormFinder and BestKeeper. A final comprehensive ranking determined that EIF and 18S were the optimal reference genes for expression data from different tissues, while TBP and EIF were optimal for expression data from different ovarian developmental stages. To our knowledge, this is the first systematic analysis of reference genes for normalization of qRT-PCR data in P. clarkii. These results will facilitate more accurate and reliable expression studies of this and other crustacean species.

Behavioural and transcriptional changes in the amphipod Echinogammarus marinus exposed to two antidepressants, fluoxetine and sertraline

In the past decade, there have been increasing concerns over the effects of pharmaceutical compounds in the aquatic environment, however very little is known about the effects of antidepressants such as the selective serotonin re-uptake inhibitors (SSRIs). Many biological functions within invertebrates are under the control of serotonin, such as reproduction, metabolism, moulting and behaviour. The effects of serotonin and fluoxetine have recently been shown to alter the behaviour of the marine amphipod, Echinogammarus marinus (Leach, 1815). The purpose of this study was to observe behavioural and transcriptional modifications in this crustacean exposed to the two most prescribed SSRIs (fluoxetine and sertraline) and to develop biomarkers of neurological endocrine disruption. The animals were exposed to both drugs at environmentally relevant concentrations from 0.001 to 1μg/L during short-term (1h and 1day) and medium-term (8 days) experiments. The movement of the amphipods was tracked using the behavioural analysis software during 12min alternating dark/light conditions. The behavioural analysis revealed a significant effect on velocity which was observed after 1h exposure to sertraline at 0.01μg/L and after 1 day exposure to fluoxetine as low as 0.001μg/L. The most predominant effect of drugs on velocity was recorded after 1 day exposure for the 0.1 and 0.01μg/L concentrations of fluoxetine and sertraline, respectively. Subsequently, the expression (in this article gene expression is taken to represent only transcription, although it is acknowledged that gene expression can also be regulated at translation, mRNA and protein stability levels) of several E. marinus neurological genes, potentially involved in the serotonin metabolic pathway or behaviour regulation, were analysed in animals exposed to various SSRIs concentrations using RT-qPCR. The expression of a tryptophan hydroxylase (Ph), a neurocan core protein (Neuc), a Rhodopsin (Rhod1) and an Arrestin (Arr) were measured following exposure to fluoxetine or sertraline for 8 days. The levels of Neuc, Rhod1 and Arr were significantly down-regulated to approximately 0.5-, 0.29- and 0.46-fold, respectively, for the lower concentrations of fluoxetine suggesting potential changes in the phototransduction pathway. The expression of Rhod1 tended to be up-regulated for the lower concentration of sertraline but not significantly. In summary, fluoxetine and sertraline have a significant impact on the behaviour and neurophysiology of this amphipod at environmentally relevant concentrations with effects observed after relatively short periods of time.

The intestinal microbiome of fish under starvation

Background

Starvation not only affects the nutritional and health status of the animals, but also the microbial composition in the host’s intestine. Next-generation sequencing provides a unique opportunity to explore gut microbial communities and their interactions with hosts. However, studies on gut microbiomes have been conducted predominantly in humans and land animals. Not much is known on gut microbiomes of aquatic animals and their changes under changing environmental conditions. To address this shortcoming, we determined the microbial gene catalogue, and investigated changes in the microbial composition and host-microbe interactions in the intestine of Asian seabass in response to starvation.

Results

We found 33 phyla, 66 classes, 130 orders and 278 families in the intestinal microbiome. Proteobacteria (48.8%), Firmicutes (15.3%) and Bacteroidetes (8.2%) were the three most abundant bacteria taxa. Comparative analyses of the microbiome revealed shifts in bacteria communities, with dramatic enrichment of Bacteroidetes, but significant depletion of Betaproteobacteria in starved intestines. In addition, significant differences in clusters of orthologous groups (COG) functional categories and orthologous groups were observed. Genes related to antibiotic activity in the microbiome were significantly enriched in response to starvation, and host genes related to the immune response were generally up-regulated.

Conclusions

This study provides the first insights into the fish intestinal microbiome and its changes under starvation. Further detailed study on interactions between intestinal microbiomes and hosts under dynamic conditions will shed new light on how the hosts and microbes respond to the changing environment.

Insights into the prostanoid pathway in the ovary development of the penaeid shrimp Penaeus monodon

The prostanoid pathway converts polyunsaturated fatty acids (PUFAs) into bioactive lipid mediators, including prostaglandins, thromboxanes and prostacyclins, all of which play vital roles in the immune and reproductive systems in most animal phyla. In crustaceans, PUFAs and prostaglandins have been detected and often associated with female reproductive maturation. However, the presence of prostanoid biosynthesis genes remained in question in these species. In this study, we outlined the prostanoid pathway in the black tiger shrimp Penaeus monodon based on the amplification of nine prostanoid biosynthesis genes: cytosolic phospholipase A2, hematopoietic prostaglandin D synthase, glutathione-dependent prostaglandin D synthase, prostaglandin E synthase 1, prostaglandin E synthase 2, prostaglandin E synthase 3, prostaglandin F synthase, thromboxane A synthase and cyclooxygenase. TBLASTX analysis confirmed the identities of these genes with 51-99% sequence identities to their closest homologs. In addition, prostaglandin F2α (PGF2α), which is a product of the prostaglandin F synthase enzyme, was detected for the first time in P. monodon ovaries along with the previously identified PUFAs and prostaglandin E2 (PGE2) using RP-HPLC and mass-spectrometry. The prostaglandin synthase activity was also observed in shrimp ovary homogenates using in vitro activity assay. When prostaglandin biosynthesis was examined in different stages of shrimp ovaries, we found that the amounts of prostaglandin F synthase gene transcripts and PGF2α decreased as the ovaries matured. These findings not only indicate the presence of a functional prostanoid pathway in penaeid shrimp, but also suggest a possible role of the PGF2α biosynthesis in shrimp ovarian development.