Transcriptomic analysis reveals growth-related genes in juvenile grass carp, Ctenopharyngodon idella

Identification of Genes Related to Resistance to Ichthyophthirius multifiliis Based on Co-expression Network Analysis in Grass Carp

The ciliate protozoan Ichthyophthirius multifiliis is an essential parasite causing white spot disease in grass carp, leading to significant economic losses. Understanding the molecular basis of grass carp's response to I. multifiliis has important scientific and environmental values. The transcriptional network analysis offers a valuable strategy to decipher the changes in gene expression in grass carp infected with I. multifiliis. Our goal was to screen the genes and pathways involved in resistance to I. multifiliis in grass carp. The different traits exhibited by grass carp infected with I. multifiliis may be caused by the differences in gene expression among grass carp individuals. Herein, to reveal those resistance-associated genes against I. multifiliis infection, we performed RNA sequencing using weighted gene co-expression network analysis (WGCNA). The biological function analysis and hub gene annotation for highly relevant modules revealed that different pathogen recognition and clearance responses resulted in different resistance to I. multifiliis infection. Furthermore, gene enrichment analysis revealed that I. multifiliis invasion in the disease-resistant group mainly activated immune pathways, including scavenger receptor activity and kappa B kinase/NF-kappa B signaling. By the annotation of the highly correlated module of the hub gene, we revealed that the apoptosis and ribosome biogenesis-related genes were enriched in the disease-resistant grass carp. The results of the dark grey module showed that several genes were mainly enriched in the two-component system (ko02020) and steroid biosynthesis (ko00100), suggesting that they are resistance-associated and energy metabolism-associated genes. In the disease resistance group, hub genes mainly included Nlrc3, fos, AAP8, HAP2, HAX, cho2, and zgc:113,036. This study revealed the gene network associated with disease resistance after I. multifiliis infection. The disease resistance-related pathways and central genes identified in this study are candidate references for breeders breeding disease-resistant. The results of this study may also provide some references for the development of drugs to antagonize I. multifiliis infection.

Integrated analysis of mRNAs and lncRNAs reveals candidate marker genes and potential hub lncRNAs associated with growth regulation of the Pacific Oyster, Crassostrea gigas

BACKGROUND

The Pacific oyster, Crassostrea gigas, is an economically important shellfish around the world. Great efforts have been made to improve its growth rate through genetic breeding. However, the candidate marker genes, pathways, and potential lncRNAs involved in oyster growth regulation remain largely unknown. To identify genes, lncRNAs, and pathways involved in growth regulation, C. gigas spat was cultured at a low temperature (15 ℃) to yield a growth-inhibited model, which was used to conduct comparative transcriptome analysis with spat cultured at normal temperature (25 ℃).

RESULTS

In total, 8627 differentially expressed genes (DEGs) and 1072 differentially expressed lncRNAs (DELs) were identified between the normal-growth oysters (cultured at 25 ℃, hereinafter referred to as NG) and slow-growth oysters (cultured at 15 ℃, hereinafter referred to as SG). Functional enrichment analysis showed that these DEGs were mostly enriched in the AMPK signaling pathway, MAPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process. LncRNAs analysis identified 265 cis-acting pairs and 618 trans-acting pairs that might participate in oyster growth regulation. The expression levels of LNC_001270, LNC_003322, LNC_011563, LNC_006260, and LNC_012905 were inducible to the culture temperature and food abundance. These lncRNAs were located at the antisense, upstream, or downstream of the SREBP1/p62, CDC42, CaM, FAS, and PIK3CA genes, respectively. Furthermore, the expression of the trans-acting lncRNAs, including XR_9000022.2, LNC_008019, LNC_015817, LNC_000838, LNC_00839, LNC_011859, LNC_007294, LNC_006429, XR_002198885.1, and XR_902224.2 was also significantly associated with the expression of genes enriched in AMPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process.

CONCLUSIONS

In this study, we identified the critical growth-related genes and lncRNAs that could be utilized as candidate markers to illustrate the molecular mechanisms underlying the growth regulation of Pacific oysters.

Identification of differentially expressed genes and SNPs linked to harvest body weight of genetically improved rohu carp, Labeo rohita

In most of the aquaculture selection programs, harvest body weight has been a preferred performance trait for improvement. Molecular interplay of genes linked to higher body weight is not elucidated in major carp species. The genetically improved rohu carp with 18% average genetic gain per generation with respect to harvest body weight is a promising candidate for studying genes' underlying performance traits. In the present study, muscle transcriptome sequencing of two groups of individuals, with significant difference in breeding value, belonging to the tenth generation of rohu carp was performed using the Illumina HiSeq 2000 platform. A total of 178 million paired-end raw reads were generated to give rise to 173 million reads after quality control and trimming. The genome-guided transcriptome assembly and differential gene expression produced 11,86,119 transcripts and 451 upregulated and 181 downregulated differentially expressed genes (DEGs) between high-breeding value and low-breeding value (HB & LB) groups, respectively. Similarly, 39,158 high-quality coding SNPs were identified with the Ts/Tv ratio of 1.23. Out of a total of 17 qPCR-validated transcripts, eight were associated with cellular growth and proliferation and harbored 13 SNPs. The gene expression pattern was observed to be positively correlated with RNA-seq data for genes such as myogenic factor 6, titin isoform X11, IGF-1 like, acetyl-CoA, and thyroid receptor hormone beta. A total of 26 miRNA target interactions were also identified to be associated with significant DETs (p-value < 0.05). Genes such as Myo6, IGF-1-like, and acetyl-CoA linked to higher harvest body weight may serve as candidate genes in marker-assisted breeding and SNP array construction for genome-wide association studies and genomic selection.

RNA-Seq Transcriptome Analysis of the Liver and Brain of the Black Carp (Mylopharyngodon piceus) During Fasting

The black carp (Mylopharyngodon piceus) is an important carnivorous freshwater-cultured species. To understand the molecular basis underlying the response of black carp to fasting, we used RNA-Seq to analyze the liver and brain transcriptome of fasting fish. Annotation to the NCBI database identified 66,609 unigenes, of which 22,841 were classified into the Gene Ontology database and 15,925 were identified in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Comparative analysis of the expression profile between fasting and normal feeding fish revealed 13,737 differentially expressed genes (P < 0.05), of which 12,480 were found in liver tissue and 1257 were found in brain tissue. The KEGG pathway analysis showed significant differences in expression of genes involved in metabolic and immune pathways, such as the insulin signaling pathway, PI3K-Akt signaling pathway, cAMP signaling pathway, FoxO signaling pathway, AMPK signaling pathway, endocytosis, and apoptosis. Quantitative real-time PCR analysis confirmed that expression of the genes encoding the factors involved in those pathways differed between fasting and feeding fish. These results provide valuable information about the molecular response mechanism of black carp under fasting conditions.