Identification and expression pattern of the sex determination gene fruitless-like in Cherax quadricarinatus.

The fruitless (fru) gene has an important function in the courtship behavior and sex determination pathway of Drosophila melanogaster; however, the fru gene has never been reported in shrimps. In this study, the fruitless-like gene was identified in Cherax quadricarinatus (Cqfru) and is reported here for the first time. A sequence analysis revealed a conserved BTB domain in Cqfru which is the same as fru in D. melanogaster. An analysis of the expression level of Cqfru showed that it was highly expressed in the gastrula stage during embryonic development. Furthermore, in situ hybridization and expression distribution in tissues showed that its sexually dimorphic expression may be focused on the hepatopancreas, brains, and gonads. The gonads, brains, and hepatopancreas of males had a higher expression level of Cqfru than those of females; however, the expression level of the abdominal ganglion was found to be higher in females than in males in this study. The results of an RNA interference treatment showed that a knockdown of Cqfru reduced the expression of the insulin-like androgenic gland hormone (IAG) and tumor necrosis factor (TNF). The characteristic fru gene in shrimps is reported here for the first time, with the results providing basic information for research into the sex-determination mechanism in C. quadricarinatus.

Effects of gonadotropin-releasing hormone analog (GnRHa) immunization on the gonadal transcriptome and proteome of tilapia (Oreochromis niloticus).

Gonadotropin releasing hormone (GnRH) plays an important role in the regulation of vertebrate reproduction. Studies have shown that immunization against GnRHa can induce sexually sterile tilapia. To explore the mechanism behind this, in this study, RNA-seq and data-independent acquisition (DIA) techniques were used to study the transcriptome and proteome of the gonad of tilapia immunized with GnRHa. 644 differentially expressed genes (80 upregulated and 564 downregulated) and 1150 differentially expressed proteins (351 upregulated and 799 downregulated) were identified. There were 209 genes with consistent differential expression patterns in the transcriptomic and proteomic analyses, of which 9 were upregulated and 200 downregulated, indicating that the gonad gene expression was inhibited by GnRHa immunization. The downregulated genes were particularly involved in the functions of single-organism process, binding, cellular process, metabolic process and catalytic activity, and associated with the pathways including ECM-receptor interaction, focal adhesion, cardiac muscle contraction and oxidative phosphorylation. The expression of six differentially expressed genes involved in the GnRH signaling pathway was all downregulated. In addition, several important functional genes related to gonadal development after GnRHa immunization were screened. This study confirmed the expression of corresponding genes was affected by GnRHa on the gonad development in tilapia at the molecular level, and laid a foundation for elucidating the mechanism of GnRHa immunization.

Single molecule real-time sequencing and RNA-seq unravel the role of long non-coding and circular RNA in the regulatory network during Nile tilapia (Oreochromis niloticus) infection with Streptococcus agalactiae.

BACKGROUND: The tilapia aquaculture industry is facing heavy economic losses due to Streptococcus agalactiae (S. agalactiae) infections. While progress has been made in past years, the lack of a high-quality tilapia genome and transcript annotations makes systematic and comprehensive exploration for a non-coding RNA regulatory network associated with the infection process unfeasible, and it stunts further research focused on disease defense and treatment. Herein, single molecular real time sequencing (SMRT-Seq) and RNA-seq data were utilized to generate a high-quality transcript annotation. In addition, Changes in mRNA and non-coding RNA expression were also analyzed during a S. agalactiae infection in tilapia. FINDINGS: In total, 16.79 Gb of clean data were obtained by sequencing on six SMRT cells, with 712,294 inserts (326,645 full-length non-chimeric reads and 354,188 non-full-length reads). A total of 197,952 consensus transcripts were obtained. Additionally, 55,857 transcript sequences were acquired, with 12,297 previously annotated and 43,560 newly identified transcripts. To further examine the immune response in Oreochromis niloticus following a S. agalactiae infection, a total of 470.62 Gb of clean data was generated by sequencing a library containing 18 S. agalactiae infected tilapia samples. Of the identified genes, 9911 were newly exploited, of which 7102 were functional annotated. Furthermore, 7874 mRNAs, 1281 long non-coding RNAs (out of 21,860 long non-coding RNAs), and 61 circular RNAs (out of 1026 circular RNAs) were found to be differentially expressed during infection, with the 1026 circRNAs not previously identified in tilapia. Moreover, k-means clustering and WGCNA analyses revealed that the immune response of tilapia to a S. agalactiae infection can be divided into three stages: cytokines driven rapid immune response, energy metabolism promotion, and the production of lysosomes and phagosomes. During this response, the head kidney and spleen have synergistic effects, while maintaining independent characteristics. Finally, lncRNA-mRNA (trans and cis), lncRNA-miRNA-mRNA, and circRNA-miRNA-mRNA regulatory networks were constructed and revealed that non-coding RNA is involved in the regulation of immune-related genes. CONCLUSIONS: This study generated a greatly-improved transcript annotation for tilapia using long-read PacBio sequencing technology, and revealed the presence of a regulatory network comprised of non-coding RNAs in Nile tilapia infected with S. agalactiae.

Transcriptomics of Cherax quadricarinatus hepatopancreas during infection with Decapod iridescent virus 1 (DIV1).

Cherax quadricarinatus is a large-sized, highly fecund, and fast-growing species of freshwater crayfish, and has become one of the world's most intensely studied crustaceans. Decapod iridescent virus 1 (DIV1), a newly described species in the family Iridoviridae, is known to infect various crustaceans, including C. quadricarinatus, and may pose a new threat in the shrimp-farming industry. The present study performed de novo transcriptome sequencing of C. quadricarinatus hepatopancreas during DIV1 infection. A total of 114,784 transcripts and 56,418 genes were obtained; 1070 genes were upregulated and 775 genes were downregulated when compared with the uninfected samples (controls). Three pattern recognition receptor genes (fibrinogen-related protein, C-type lectin, and beta-1,3-glucan-binding protein) were upregulated during DIV1 infection. Among the top-30 upregulated unigenes, 9 unigenes were identified as vitellogenin (Vg) genes, and the top-3 upregulated unigenes were identified as involved in Vg lipid transport, lipid localization, and lipid transporter activity, which were all significantly over-representative GO terms in the GO enrichment analysis of total and upregulated differentially expressed genes (DEGs). Many genes associated with Jak-STAT signaling pathway, Endocytosis, Phagosome, MAPK signaling pathway, Apoptosis and Lysosome were positively modified after DIV1 infection. The predicted protein-protein interaction (PPI) analysis showed NF1 and TUBA, CRM1 and TUBB were involved in protein interactions. This research showed that DIV1 infection has a significant impact on the transcriptome profile of C. quadricarinatus hepatopancreas, and the results enhance our understanding of virus-host interactions. Furthermore, the high number of transcripts generated in the present study will provide information for identifying novel genes in the absence of a full C. quadricarinatus genome sequence.

Establishment of pyrosequencing technology to detect White Spot Syndrome Virus (WSSV) in cultured aquatic animals.

This study aimed to establish pyrosequencing methods to detect white spot syndrome virus (WSSV). One pair of polymerase chain reaction (PCR) primers, and one pyrosequencing primer, were designed for WSSV. The pyrosequencing reaction system and conditions were optimized and a pyrosequencing method for detecting WSSV was successfully established. This method was able to specifically detect WSSV in eight viruses, with high sensitivity. The minimum detectable limit for nucleic acid was 23 copies/µL. The method was verified by detecting WSSV in 1881 batches of samples collected from domestic and imported shrimps. The detection results were more sensitive than conventional PCR. This research has therefore provided a new detection method for monitoring, and controlling aquatic animal virus diseases.

LncRNA and mRNA profiling during activation of tilapia macrophages by HSP70 and Streptococcus agalactiae antigen.

OBJECTIVES: To investigate the lncRNA profiling during tilapia peritoneal macrophages (TPMs) activation and discuss the relationship between lncRNA and mRNA. MATERIALS AND METHODS: RNA sequencing was used to investigate the lncRNA and mRNA profiles of TPMs activation following stimulation with Streptococcus agalactiae (Sa) antigen, heat shock protein 70 (HSP70) and HSP70+Sa. The expressions of lncRNA and mRNA were confirmed by qPCR. 356 lncRNA, 10173 mRNA and 1782 transcripts of uncertain coding potential (TUCP) were differentially expressed by pairwise comparison. These lncRNAs were shorter in length, fewer in exon number and higher in expression levels as compared with mRNAs. 683 lncRNAs and 4320 mRNAs were co-located, while 316 lncRNAs and 9997 mRNAs were in co-expression networks. Seven mRNAs (ANKRD34A, FMODA, GJA3, CNTN5, BMP10, BAI2 and HS3ST6) were involved in both networks of LNC_00035 and LNC_000466. Differentially expressed genes were involved in signaling pathways, such as "phosphorylation", "cytokine-cytokine receptor interaction", "endocytosis" and "MHC protein complex". LNC_000792, LNC_000215, LNC_000035 and LNC_000310, with cis and/or trans relationships with mRNAs, were also involved in ceRNA network. CONCLUSIONS: These results might represent the first identified expression profile of lncRNAs and mRNAs in tilapia macrophages activated by HSP70 and Sa.

A physiologically based toxicokinetic and toxicodynamic model links the tissue distribution of benzo[a]pyrene and toxic effects in the scallop Chlamys farreri.

A physiologically based toxicokinetic and toxicodynamic (PBTK-TD) model was developed for benzo[a]pyrene (B[a]P) in scallop Chlamys farreri. The PBTK model structure consisted of gill, digestive gland, adductor muscle, hemolymph and other tissues. In TD modeling, aryl hydrocarbon hydroxylase (AHH) activity assay, comet assay, protein carbonyl measurement and lipid peroxidation level determination in digestive gland were used as biomarkers to reflect toxic effects. We integrated B[a]P concentration and biomarkers by using sigmoid Emax model in digestive gland. The PBTK-TD model predicted the B[a]P concentrations within each organ compartment and the toxic effects in digestive gland. The results showed that the predicted and measured data in different organs were in good agreement and comet assay was considered as the best biomarker. This model would serve as a useful tool for pollution monitoring and food security.

Application of the biomarker responses in scallop (Chlamys farreri) to assess metals and PAHs pollution in Jiaozhou Bay, China.

Scallops (Chlamys farreri) were collected at three sites in the Jiaozhou Bay, China in April and July 2007. The responses of biomarkers were evaluated for assessment of physiological status of scallops and pollution in the areas. Compared with site S1, the activity of aryl hydrocarbon hydroxylase was induced notably (P < 0.05) and glutathione S-transferase activity was inhibited significantly (P < 0.05) in S2 and S3. The levels of DNA alkaline unwinding, DNA protein crosslinks, protein carbonyl and lipid peroxidation were tested and indicated the oxidative stress situation of scallops. There was significant difference of biomarker levels between sampling seasons (P < 0.05), and between digestive gland and gill (P < 0.05). The results provided the reference data for multiple-pollution assessment in the marine environment and indicated that tissue type and seasons affect biomarkers, therefore these factors should be taken into consideration when biomarkers are used for environmental assessment purposes.