Analysis of the black-chinned tilapia Sarotherodon melanotheron heudelotii reproducing under a wide range of salinities: from RNA-seq to candidate genes

Integrative transcriptome analysis and discovery of signaling pathways involved in the protective effects of curcumin against oxidative stress in tilapia hepatocytes

Summer outbreaks of the hepatobiliary syndrome in fish impose a heavy burden on aquaculture in China. Curcumin is a polyphenol with antioxidant activity that has been used to protect the health of fish livers, but the mechanism underlying its protective effect is unclear. In this study, an in vitro model of hepatocyte oxidative damage in Oreochromis niloticus was established using H₂O₂. Treatment with 5 mM H₂O₂ for 2.5 h markedly reduced cell viability and antioxidant activity and elevated lactate dehydrogenase (LDH) activity, indicating conditions that can be used to establish an oxidative stress model. Under H₂O₂ stress, curcumin pretreatment significantly maintained cell viability, reduced malondialdehyde (MDA) levels, and increased superoxide dismutase (SOD) activity. RNA-seq results showed that acute H₂O₂ treatment resulted in minor changes in gene expression, whereas curcumin changed the expression profile and affected cytochrome P450 (Cyp 450), glutathione (GSH) metabolism, and the peroxisome proliferator-activated receptor (PPAR) signaling pathway. Several critical antioxidant defense signaling pathways were identified, and altered expression was confirmed by q-PCR. These results indicate that curcumin might upregulate PPAR expression by increasing Cyp2J2 expression. Further experiments showed that curcumin can upregulate the Nrf2-Keap1 signaling pathway at the transcriptional level, and this upregulation can induce downstream defense genes, including glutamate cysteine ligase catalytic subunit(GCLC) and glutamate cysteine ligase modifier subunit (GCLM), and thereby promote GSH synthesis and the expression of related antioxidases. This study might shed light on the effects of curcumin on the prevention and alleviation of liver diseases in fish.

Phenotypic plasticity in gene expression and physiological response in red drum Sciaenops ocellatus exposed to a long-term freshwater environment

Red drum (Sciaenops ocellatus) is a euryhaline fish commonly found in the Gulf of Mexico and along the Atlantic coast of North America. Because of high commercial demand and its euryhaline characteristics, aquaculture of this species has diversified from marine to low-salinity aquaculture systems. In recent years, interest in the feasibility of producing red drum in inland freshwater systems has grown and this prompted us to investigate its osmoregulatory capacity after rearing for 8 months in a freshwater aquaculture system. We compared the activities of several genes and enzymes involved in the osmoregulatory process in freshwater-acclimatized (FW) and seawater (SW) red drum. The gene expression profiles were variable: the expression of genes encoding Na⁺/K⁺-ATPase (NKA) and the cystic fibrosis transmembrane regulator (CFTR) was slightly higher in SW than FW fish, while phosphoenolpyruvate carboxykinase (PEPCK) and the glucocorticoid receptor messenger RNA (mRNA) levels were higher in FW red drum. The total plasma K concentration was 60.3% lower, and gill NKA activity was 63.5% lower in FW than in SW fish. PEPCK activity was twofold higher in FW than in SW red drum. Similarly, liver glycogen was 60% higher in FW fish. In summary, both gene expression and the enzyme activity data support the phenotypic plasticity of red drum and suggest that the limited capacity for ion homeostasis observed, in particular the low plasma K concentration, was due to the composition of freshwater and does not necessarily reflect a physiological inability to osmoregulate.

Identifying a Major QTL Associated with Salinity Tolerance in Nile Tilapia Using QTL-Seq

Selection of new lines with high salinity tolerance allows for economically feasible production of tilapias in brackish water areas. Mapping QTLs and identifying the markers linked to salinity-tolerant traits are the first steps in the improvement of the tolerance in tilapia through marker-assisted selection techniques. By using QTL-seq strategy and linkage-based analysis, two significant QTL intervals (chrLG4 and chrLG18) on salinity-tolerant traits were firstly identified in the Nile tilapia. Fine mapping with microsatellite and SNP markers suggested a major QTL region that located at 23.0 Mb of chrLG18 and explained 79% of phenotypic variation with a LOD value of 95. Expression analysis indicated that at least 10 genes (e.g., LACTB2, KINH, NCOA2, DIP2C, LARP4B, PEX5R, and KCNJ9) near or within the QTL interval were significantly differentially expressed in intestines, brains, or gills under 10, 15, or 20 ppt challenges. Our findings suggest that QTL-seq can be effectively utilized in QTL mapping of salinity-tolerant traits in fish. The identified major QTL is a promising locus to improve our knowledge on the genetic mechanism of salinity tolerance in tilapia.

Genome-Wide QTL Analysis Identified Significant Associations Between Hypoxia Tolerance and Mutations in the GPR132 and ABCG4 Genes in Nile Tilapia

Exposure to hypoxia induces both acute and chronic stress responses, which plays an important role in health of cultured organisms including growth, reproduction, immunity, and other energy demanding activities. Application of advanced genomic technologies allows rapid identification of hypoxia trait-associated genes and precise selection of superior brood stocks with high tolerance in tilapia. By applying QTL-seq and double-digest restriction-site associated DNA sequencing (ddRAD-seq) techniques, we identified four genome-wide significant quantitative trait loci (QTLs) for hypoxia tolerance and many suggestive QTLs in Nile tilapia. These QTLs explained 6.6-14.7% of the phenotypic variance. Further analysis revealed that single nucleotide polymorphisms (SNPs) in exons of both GPR132 and ABCG4 genes located in genome-wide QTL intervals were significantly associated with hypoxia-tolerant traits. Expression analysis of both genes suggested that they were strong candidate genes involved into hypoxia tolerance in tilapia. Our findings suggest that both QTL-seq and ddRAD-seq techniques can be effectively utilized in QTL mapping of hypoxia traits in fish. Our data supply a basis for further marker-assisted selection of super lines with a high level of tolerance against low oxygen stress in the tilapia.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.

Plasticity of gene expression according to salinity in the testis of broodstock and F1 black-chinned tilapia, Sarotherodon melanotheron heudelotii

The black-chinned tilapia Sarotherodon melanotheron heudelotii Rüppell 1852 (Teleostei, Cichlidae) displays remarkable acclimation capacities. When exposed to drastic changes of salinity, which can be the case in its natural habitat, it develops quick physiological responses and keeps reproducing. The present study focused on the physiological impact of salinity on male reproductive capacities, using gene expression as a proxy of acclimation process. Two series of experimental fish were investigated: the first one was composed of fish maintained in freshwater for several generations and newly acclimated to salinities of 35 and 70, whereas the second one consisted of the descendants of the latter born and were raised under their native salinity. Expression patterns of 43 candidate genes previously identified from the testes of wild males was investigated in the three salinities and two generations. Twenty of them showed significant expression differences between salinities, and their predicted function revealed that most of them are involved in the osmotic tolerance of sperm cells and/or in the maintenance of sperm motility. A high level of expression variation was evidenced, especially for fish maintained in freshwater. In spite of this, gene expression patterns allowed the differentiation between fish raised in freshwater and those maintained in hypersaline water in both generations. Altogether, the results presented here suggest that this high variability of expression is likely to ensure the reproductive success of this species under varying salinities.

Comparative transcriptome analysis of two oysters, Crassostrea gigas and Crassostrea hongkongensis provides insights into adaptation to hypo-osmotic conditions

Environmental salinity creates a key barrier to limit the distribution of most aquatic organisms. Adaptation to osmotic fluctuation is believed to be a factor facilitating species diversification. Adaptive evolution often involves beneficial mutations at more than one locus. Bivalves hold great interest, with numerous species living in waters, as osmoconformers, who maintain the osmotic pressure balance mostly by free amino acids. In this study, 107,076,589 reads from two groups of Crassostrea hongkongensis were produced and the assembled into 130,629 contigs. Transcripts putatively involved in stress-response, innate immunity and cell processes were identified according to Gene ontology and KEGG pathway analyses. Comparing with the transcriptome of C. gigas to characterize the diversity of transcripts between species with osmotic divergence, we identified 182,806 high-quality single nucleotide polymorphisms (SNPs) for C. hongkongensis, and 196,779 SNPs for C. gigas. Comparison of 11,602 pairs of putative orthologs allowed for identification of 14 protein-coding genes that experienced strong positive selection (Ka/Ks>1). In addition, 45 genes that may show signs of moderate positive selection (1 ≥ Ka/Ks>0.5) were also identified. Based on Ks ratios and divergence time between the two species published previously, we estimated a neutral transcriptome-wide substitution mutation rate of 1.39 × 10(-9) per site per year. Several genes were differentially expressed across the control and treated groups of each species. This is the first time to sequence the transcriptome of C. hongkongensis and provide the most comprehensive transcriptomic resource available for it. The increasing amount of transcriptome data on Crassostrea provides an excellent resource for phylogenetic analysis. A large number of SNPs identified in this work are expected to provide valuable resources for future marker and genotyping assay development. The analysis of natural selection provides an innovative view on the adaptation within species and sets the basis for future genetic and evolutionary studies.