Transcriptome analysis reveals molecular mechanisms responsive to acute cold stress in the tropical stenothermal fish tiger barb (Puntius tetrazona)

Molecular underpinnings underlying behaviors changes in the brain of juvenile common carp (Cyrinus carpio) in response to warming

INTRODUCTION

Global warming is increasing interest in how aquatic animals can adjust their physiological performance and cope with temperature changes. Therefore, understanding the behavioral changes and molecular underpinnings in fish under warming is crucial for both the individual and groups survival. This could provide experimental evidence and resource for evaluating the impact of global warming.

OBJECTIVE

Three genetic families of common carp (Cyprinus carpio) were generated. These juveniles were constructed short-term (4 days) and long-term (30 days) warming groups to investigate the effects of warming on behavioral responses and to elucidate the potential underlying mechanisms of warming-driven behavior.

METHODS

Behavioral tests were used to explore the effects of short- and long-term exposure to warming on the swimming behavior of C. carpio. Brain transcriptome combined with measurement of nervous system activity was used to further investigated the comprehensive neuromolecular mechanisms under warming.

RESULTS

Long-term warming groups had a more significant impact on the decline of swimming behavior in juvenile C. carpio. Furthermore, brain comparative transcriptomic analysis combined with measurement of nervous system activity revealed that genes involved in cytoskeletal organization, mitochondrial regulation, and energy metabolism are major regulators of behavior in the juvenile under warming. Importantly, especially in the long-term warming groups, enrichment analysis of associated gene expression suggested functional alterations of synaptic transmission and signal transduction leading to swimming function impairment in the central nervous system, as revealed by behavioral tests.

CONCLUSIONS

Our study provides evidence of the neurogenomic mechanism underlying the decreased swimming activity in juvenile C. carpio under warming. These findings have important implications for understanding the impacts of climate change on aquatic ecosystems and the organisms that inhabit them.

Transcriptomics reveals that NAD(P)H affects the development of the Zig-zag eel (Mastacembelus armatus ♀) × Spiny eel (Sinobdella sinensis ♂) hybrid offspring leading to low hatching rates

Zig-zag eel (Mastacembelus armatus (2 n = 48)) and Spiny eel (Sinobdella sinensis (2 n = 48)) are two species of the Mastacembelidae family commonly found in southern China. Hybridization between the two has a very high deformity rate and a very low hatching rate. In order to investigate the reasons for this, the first hybridization between M. armatus and S. sinensis was carried out using artificial insemination, and the embryonic development of the hybrid offspring was examined using microphotography, and the malformations of the hybrid offspring were investigated by transcriptomics. The experiments showed that the average egg production was 4265.7 ± 322.94 (Mean ± SD), the average fertilization rate of hybrid offspring was 98.67 ± 0.58 % (Mean ± SD), the hatching rate was 12.06 ± 3.44 % (Mean ± SD), the deformity rate was 98.15 ± 3.21 % (Mean ± SD), and the embryonic development successively went through the five main stages of fertilized egg, egg cleavage, embryo formation, organogenesis, and exertion of membranes. Transcriptomics showed that the expression of NAD(P)H-related enzyme activity DEGs was increased, and many DEGs related to cell signaling molecule transmission and metabolic regulation are enriched in KEGG pathways, such as IL-17 signaling pathway, Osteoclast differentiation, TNF signaling pathway and MAPK signaling pathway. etc. The major types of DEGs corresponded to those coding for proteins. This study suggests that the high malformation rate in hybrid offspring may be caused by impaired synthesis of proteins during embryonic development.

Histology, fatty acid composition, antioxidant and glycolipid metabolism, and transcriptome analyses of the acute cold stress response in Phoxinus lagowskii

Water temperature is a crucial environmental factor that significantly affects the physiological and biochemical processes of fish. Due to the occurrence of cold events in aquaculture, it is imperative to investigate how fish respond to cold stress. This study aims to uncover the mechanisms responds to acute cold stress by conducting a comprehensive analysis of the histomorphology, glycolipid metabolic and antioxidant enzymes, fatty acid composition and transcriptome at three temperatures (16 °C, 10 °C and 4 °C) in Phoxinus lagowskii. Our results showed that cold stress not damaged muscle microstructure but caused autophagy (at 10 °C). In addition, serum glucose (Glu) and triglycerides (TG) increased during cold stress. The activities of reactive oxygen species (ROS), superoxide dismutase (SOD), catalase (CAT), fructose phosphokinase (PFK), hexokinase (HK), pyruvate kinase (PK), and malondialdehyde (MDA) content in muscle were measured and analyzed. During cold stress, superoxide dismutase and catalase activities increased, reactive oxygen species content decreased. No significant difference in Glutathione peroxidase (GPx) activity, malondialdehyde and total cholesterol (T-CHO) contents among groups. Phosphokinase and pyruvate kinase activities decreased, and HK activity increased during cold stress. Our study resulted in the identification of a total of 25,400 genes, with 2524 genes showing differential expression across different temperature treatments. Furthermore, KEGG pathway indicated that some pathways upregulated during light cold stress (at 10 °C, including autophagy, and AMP-activated protein kinase (AMPK) signaling pathway. Additionally, circadian rhythm is among the most enriched pathways in genes up-regulated during severe cold stress (at 4 °C). Our findings offer valuable insights into how cold-water fish respond to cold stress.

First neurotranscriptome of adults Tambaquis (Colossoma macropomum) with characterization and differential expression between males and females

The Tambaqui is one of the most representative Amazon fish species, being highly exploited in fisheries, aquaculture and as a research model. Nonetheless, data about functional genome are still required to evaluate reproductive and nutrition parameters as well as resistance to pathogens. The of next-generation sequencing has allows assessing the transcriptional processes in non-model species by providing comprehensive gene collections to be used as a database in further genomic applications and increased performance of captive populations. In this study, we relied on RNAseq approach to generate the first transcriptome of the telencephalon from adult males and females of Colossoma macropomum, resulting in a reference dataset for future functional studies. We retrieved 896,238 transcripts, including the identification of 267,785 contigs and 203,790 genes. From this total, 91 transcripts were differentially expressed, being 63 and 28 of them positively regulated for females and males, respectively. The functional annotation resulted in a library of 40 candidate genes for females and 20 for males. The functional enrichment classes comprised reproductive processes (GO:0,048,609; GO:0,003,006; GO:0,044,703; GO:0,032,504; GO:0,019,953) being related to sex differentiation (e.g., SAFB) and immune response (e.g., SLC2A6, AHNAK, NLRC3, NLRP3 and IgC MHC I alpha3), thus indicating that the genes in the neurotranscriptome of Tambaqui participate in sex differentiation and homeostasis of captive specimens. These data are useful to design the selection of genes related to sex determination and animal welfare in raising systems of Tambaqui.

Comparative Transcriptome Analyses Provide New Insights into the Evolution of Divergent Thermal Resistance in Two Eel Gobies

Adaptation to thermal conditions in tidal mudflats always involves tolerating frequent fluctuations and often extreme environmental temperatures. Regulation of gene expression plays a fundamental role in the evolution of these thermal adaptations. To identify the key gene regulatory networks associated with the thermal adaptation, we investigated the capability of cold tolerance, as well as the transcriptomic changes under cold stress in two mudflat inhabitants (Odontamblyopus lacepedii and O. rebecca) with contrasting latitude affinity. Our results revealed a remarkable divergent capacity of cold tolerance (CTmin: 0.61 °C vs. 9.57 °C) between the two gobies. Analysis of transcriptomic changes under cold stress unveiled 193 differentially expressed genes exhibiting similar expression profiles across all tissues and species, including several classic metabolic and circadian rhythm molecules such as ACOD and CIART that may represent the core cold response machinery in eel gobies. Meanwhile, some genes show a unique expression spectrum in the more cold-tolerant O. lacepedii suggesting their roles in the enhanced cold tolerance and hence the extreme thermal adaptations. In addition, a weighted gene co-expression network analysis (WGCNA) revealed a subset of metabolic hub genes including MYH11 and LIPT2 showing distinct down-regulation in O. lacepedii when exposed to cold stress which highlights the role of reduced energy consumption in the enhanced cold tolerance of eel gobies. These findings not only provide new insights into how mudflat teleosts could cope with cold stress and their potential evolutionary strategies for adapting to their thermal environment, but also have important implications for sound management and conservation of their fishery resources in a scenario of global climate warming in the marine realm.

Acute cold stress leads to zebrafish ovarian dysfunction by regulating miRNA and mRNA

Temperature is a critical factor that regulates the reproduction processes in teleost. However, the gonadal response mechanism to cold stress in fish remains largely unknown. In the present study, female zebrafish were exposed to different extents of low temperatures at 18 °C and 10 °C for 48 h. The ovarian histology was remarkably damaged after cold stress exposure. Integrated analysis of miRNA-mRNA was used to investigate the ovarian response to acute cold stress. A large number of mRNAs and miRNAs were altered by cold stress, which are involved in extensive biological processes. It is indicated that the signal transduction of MAPK and Calcium signaling pathway is highly engaged in zebrafish ovary to adapt to cold stress. The immune system was dysregulated by cold stress while the ovarian autophagy was activated. Remarkably increased gene number related to reproductive functions was identified in the cold stress at 10 °C compared to the control. The cold stress-induced dysregulated reproductive genes include star, hsd3b1, hsd17b1, inha, insl3, amh, nanos1 and foxl2. Combined with the dysregulated insulin, IGF and progesterone signaling, it is suggested that cold stress affects ovarian function in multiple aspects, including oocyte meiosis, folliculogenesis, final maturation and ovarian maintenance. On the other hand, the ovarian miRNA-mRNA regulatory network response to cold stress was also constructed. Overall, our result revealed the ovarian response to cold stress in zebrafish and provided insight into the fish adaptation mechanism to acute temperature change.

Construction of a high-density genetic linkage map and QTL mapping of growth and cold tolerance traits in Takifugu fasciatus

Takifugu fasciatus is an aquaculture species with high economic value. In recent years, problems such as environmental pollution and inbreeding have caused a serious decline in T. fasciatus germplasm resources. In this study, a high-density genetic linkage map was constructed by whole-genome resequencing. The map consists of 4891 bin markers distributed across 22 linkage groups (LGs), with a total genetic coverage of 2381.353 cM and a mean density of 0.535 cM. Quantitative trait locus (QTL) localization analysis showed that a total of 19 QTLs associated with growth traits of T. fasciatus in the genome-wide significance threshold range, distributed on 11 LGs. In addition, 11 QTLs associated with cold tolerance traits were identified, each scattered on a different LG. Furthermore, we used QTL localization analysis to screen out three candidate genes (IGF1, IGF2, ADGRB) related to growth in T. fasciatus. Meanwhile, we screened three candidate genes (HSP90, HSP70, and HMGB1) related to T. fasciatus cold tolerance. Our study can provide a theoretical basis for the selection and breeding of cold-tolerant or fast-growing T. fasciatus.

Chromosome-Level Genome Assembly of Protosalanx chinensis and Response to Air Exposure Stress

Protosalanx chinensis is a suitable particular species for genetic studies on nearly scaleless skin, transparency and high sensitivity to hypoxia stress. Here, we generated a high-quality chromosome-level de novo assembly of P. chinensis. The final de novo assembly yielded 379.47 Mb with 28 pseudo-chromosomes and a scaffold N50 length of 14.52 Mb. In total, 21,074 protein-coding genes were predicted. P. chinensis, Esox lucius and Hypomesus transpacificus had formed a clade, which diverged about 115.5 million years ago. In the air exposure stress experiment, we found that some genes play an essential role during P. chinensis hypoxia, such as bhlh, Cry1, Clock, Arntl and Rorb in the circadian rhythm pathway. These genomic data offer a crucial foundation for P. chinensis ecology and adaptation studies, as well as a deeper understanding of the response to air exposure stress.

Transcriptome analysis of Chinese sucker (Myxocyprinus asiaticus) head kidney and discovery of key immune-related genes to cold stress after swimming fatigue

For Chinese sucker (Myxocyprinus asiaticus), passing through a dam with fast flow and cold water are always unavoidable, and this process can cause stress, disease or even death. In this study, comparative transcriptome analysis was conducted to investigate the potential immune mechanism in head kidney of M. asiaticus with swimming fatigue stress and cold stress after fatigue. In general, a total of 181,781 unigenes were generated, and 38,545 differentially expressed genes (DEGs) were identified. In these DEGs, 22,593, 7286 and 8666 DEGs were identified among groups of fatigue vs. cold, control vs. cold, and control vs. fatigue, respectively. Enrichment analysis revealed these DEGs were involved in coagulation cascades and complement, natural killer cell mediated cytotoxicity, antigen processing and presentation, Toll-like receptor signaling pathways, and chemokine signaling pathway. Notably, immune genes including heat shock protein 4a (HSP4a), HSP70 and HSP90α genes were significantly up-regulated in fishes with cold stress after fatigue. Differently, more immune genes in control vs. cold compared with that in control vs. fatigue were significantly down-regulated expression, such as claudin-15-like, Toll-like receptor 13, antimicrobial peptide (hepcidin), immunoglobulin, CXCR4 chemokine receptor, T-cell receptor, complement factor B/C2-A3, and interleukin 8. In this study, the number of DEGs in the head kidney was less than that our previous study in the spleen, which we speculated was more sensitive to changes in water temperature than the head kidney. In summary, lots of immune-related genes in the head kidney were down-regulated under cold stress after fatigue, suggesting that M. asiaticus might have experienced severe immunosuppression in the process of passing through the dam.

Comparative analysis of quantitative phosphoproteomics between two tilapias (Oreochromis niloticus and Oreochromis aureus) under low-temperature stress

As an important farmed fish, tilapia has poor tolerance to low-temperatures. At the same time, different tilapia strains have apparent differences in low-temperature tolerance. In this study, using the iTRAQ method, the phosphorylated proteomics of two tilapia strains (Oreochromis niloticus and Oreochromis aureus) with different tolerances to low-temperature stress were quantitatively and comparatively analyzed, to clarify the physiological mechanism of tilapia's response to low-temperature stress. Through the GO and IPR analyses of differentially phosphorylated proteins, a number of similarities in physiological activities and regulatory effects were found between the two tilapias in response to low-temperature stress. Many differentially phosphorylated proteins are mainly involved in lipid metabolism, cell proliferation and apoptosis. However, the difference in endurance of low temperature of these two tilapias might be related to the differences in categories, expression and modification level of genetic products which were involved in the aforementioned physiological processes. And meanwhile, the enrichment results of KEGG showed the changes of multiple immune-related and growth-related phosphorylated proteins in the cytokine-cytokine receptor interaction pathway in O. aureus are more prominent. Furthermore, the significantly enriched pathway of carbohydrate digestion and absorption in O. niloticus may indicate that low-temperature stress exerts a more severe impact on energy metabolism. The relative results would help elucidating the molecular mechanism by which tilapia responds to low-temperature stress, and developing culture of tilapia species.

Genome-wide identification of TNF receptor-associated factors in Japanese flounder (Paralichthys olivaceus) and functional analysis of resistance to temperature and Edwardsiella tarda stress

Tumor necrosis factor receptor-associated factors (TRAFs), as the signaling mediators of the tumor necrosis factor (TNFR) superfamily, toll-like receptors (TLR) and interleukin-1 receptor (IL-1R) superfamily, can activate downstream signal transduction pathways and play an important role in the body's immune process. In this study, six TRAF genes, namely PoTRAF2a, PoTRAF2b, PoTRAF3, PoTRAF4, PoTRAF6 and PoTRAF7, were identified and annotated in Japanese flounder by using bioinformatics methods. Phylogenetic analysis confirmed that TRAF genes can be divided into seven groups. Analysis of motif composition and gene structure demonstrated that all PoTRAF members were evolutionarily conserved. The expression patterns of PoTRAF genes were then further investigated in six different developmental stages and eleven tissues of healthy fish, and it was found that there were spatial and tissue specificities among the members. To investigate the immune response of Japanese flounder to abiotic and biotic stresses, we further analyzed the expression profile of PoTRAFs after temperature stress and pathogen challenge. The result showed that PoTRAF3 and PoTRAF4 were observably differentially expressed under temperature stress, indicating that they were involved in the immune response after temperature stress. The expression of PoTRAF2a, PoTRAF2b and PoTRAF4 was significantly different after E. tarda infection, suggesting that they might have antibacterial effects. These results would help to clarify the molecular roles of PoTRAF genes in the regulation of immune and inflammatory responses in Japanese flounder.

Genome-wide identification of olive flounder (Paralichthys olivaceus) SOCS genes: Involvement in immune response regulation to temperature stress and Edwardsiella tarda infection

The suppressors of cytokine signaling (SOCS) gene family participates in development and immunity through negative regulation of cytokine signaling pathways. Although the immune response of SOCS gene family members has been extensively characterized in teleost, no similar study has been reported in olive flounder yet. In our present study, a total of 13 SOCSs in olive flounder were identified and characterized systematically. By querying the SOCS sequences of ten teleost fish species, we found there were exactly more members of SOCSs in fish than mammals, which indicated that there were more duplication events occurred in fish than in higher vertebrates. Phylogenetic analysis clearly illuminated that SOCS genes were highly conserved. The analysis of gene structure and motif showed SOCS proteins of olive flounder shared a high level of sequence similarity strikingly. The expression profiles of tissues and developmental stages indicated that SOCS members had a kind of specificity in temporality and spatiality. RNA-Seq analysis of temperature stress and E. Tarda infection demonstrated SOCS members were involved in inflammatory response. In a word, our results would provide a further reference for understanding the mechanism of SOCS genes in olive flounder.

Transposable Element Tissue-Specific Response to Temperature Stress in the Stenothermal Fish Puntius tetrazona

Ray-finned fish represent a very interesting group of vertebrates comprising a variety of organisms living in different aquatic environments worldwide. In the case of stenothermal fish, thermal fluctuations are poorly tolerated, thus ambient temperature represents a critical factor. In this paper, we considered the tiger barb Puntius tetrazona, a freshwater fish belonging to the family Cyprinidae, living at 21-28 °C. We analyzed the available RNA-Seq data obtained from specimens exposed at 27 °C and 13 °C to investigate the transcriptional activity of transposable elements (TEs) and genes encoding for proteins involved in their silencing in the brain, gill, and liver. TEs are one of the tools generating genetic variability that underlies biological evolution, useful for organisms to adapt to environmental changes. Our findings highlighted a different response of TEs in the three analyzed tissues. While in the brain and gill, no variation in TE transcriptional activity was observed, a remarkable increase at 13 °C was recorded in the liver. Moreover, the transcriptional analysis of genes encoding proteins involved in TE silencing such as heterochromatin formation, the NuRD complex, and the RISC complex (e.g., AGO and GW182 proteins) highlighted their activity in the hepatic tissue. Overall, our findings suggested that this tissue is a target organ for this kind of stress, since TE activation might regulate the expression of stress-induced genes, leading to a better response of the organism to temperature changes. Therefore, this view corroborates once again the idea of a potential role of TEs in organism rapid adaptation, hence representing a promising molecular tool for species resilience.

Effects of temperature and combinational exposures on lipid metabolism in aquatic invertebrates

Studies of changes in fatty acids in response to environmental temperature changes have been conducted in many species, particularly mammals. However, few studies have considered aquatic invertebrates, even though they are particularly vulnerable to changes in environmental temperature. In this review, we summarize the process by which animals synthesize common fatty acids and point out differences between the fatty acid profiles of vertebrates and those of aquatic invertebrates. Unlike vertebrates, some aquatic invertebrates can directly synthesize polyunsaturated fatty acids (PUFAs), which can be used to respond to temperature changes. Various studies have shown that aquatic invertebrates increase the degree of saturation in their fatty acids through an increase in saturated fatty acid production or a decrease in PUFAs as the temperature increases. In addition, we summarize recent studies that have examined the complex effects of temperature and combinational stressors to determine whether the degree of saturation in aquatic invertebrates is influenced by other factors. The combined effects of carbon dioxide partial pressure, food quality, starvation, salinity, and chemical exposures have been confirmed, and fatty acid profile changes in response to high temperature were greater than those from combinational stressors.

Identification and functional analysis of dual-specificity phosphatases (DUSP) genes in Japanese flounder (Paralichthys olivaceus) against temperature and Edwardsiella tarda stress

Dual-specificity Phosphatases (DUSPs) are not only the key regulators of dephosphorylating and inactivating mitogen-activated protein kinases (MAPKs), but play a crucial role in the immune response. However, the role of DUSP genes in Japanese flounder (PoDUSPs) is still unclear. In this study, 28 DUSP genes in Japanese flounder were identified and classified based on the whole genome database. Phylogenetic analysis and protein structure analysis revealed that DUSPs had highly conserved domains in teleosts. Molecular evolution analysis indicated that the PoDUSP genes were conservative during evolution and were functional-constrained. Meanwhile, PoDUSP genes were found to express in different embryonic and larval stages which might play the role of sentinel in healthy organisms. Furthermore, PoDUSP genes' expression profiles after temperature stress and Edwardsiella tarda (E. tarda) infection were determined in Japanese flounder without precedent, and the results demonstrated that Podusp1, Podusp2 and Podusp16 were more respective to temperature variation whereas Podusp1 and Podusp6 were more respective to E. tarda infection. In summary, our results provide useful resources for understanding the immune responsibilities of DUSP genes in flatfish.

Transcriptome analysis revealed potential genes involved in thermogenesis in muscle tissue in cold-exposed lambs

Cold tolerance is an important trait for sheep raised at high altitudes. Muscle tissue, comprising 30–40% of the total body mass, produces heat during cold exposure. However, little is known about the genetic mechanisms of this tissue and its role in thermogenesis in lambs. We examined genes in skeletal muscle tissue in a cold-adapted sheep breed, Altay, and a cold-intolerant sheep breed, Hu, when exposed to low air temperature. Three ewe-lambs of each breed were maintained at −5°C and three ewe-lambs of each breed were maintained at 20°C. After cold exposure for 25 days, the longissimus dorsi of each lamb was collected, and transcriptome profiles were sequenced and analyzed. The results of RNA-seq showed that the average reads among the four groups were 11.0 Gbase. The genome mapping rate averaged 88.1% and the gene mapping rate averaged 82.5%. The analysis of differentially expressed genes (DEGs) indicated that the peroxisome proliferator-activated receptors (PPAR), cAMP, and calcium signaling pathways and muscle contraction in muscle tissue were linked to thermogenesis in cold-exposed lambs. Furthermore, PCK1 (phosphoenolpyruvate carboxykinase1) increased glyceroneogenesis in cold-exposed Altay lambs, and APOC3 (apolipoprotein C3), LPL (lipoprotein lipase), and FABP4 (fatty acid binding protein 4, adipocyte) were involved in the intake and transport of free fatty acids. In Hu sheep, cAMP biosynthesis from ATP hydrolysis was regulated by ADCY10 (adenylate cyclase) and ADORA2a (adenosine A2a receptor). Skeletal muscle contraction was regulated by MYL2 (myosin light chain 2). In conclusion, cold exposure altered the expression level of genes involved in heat production in muscle tissue. Some potential mechanisms were revealed, including calcium ion transport in the calcium signaling pathway, fatty acid metabolism in the PPAR signaling pathway, and cAMP biosynthesis in the cAMP signaling pathway. This study implied that skeletal muscle plays an important role in thermoregulation in lambs.

Cell Junction and Vesicle Trafficking-Mediated Melanosome/Melanin Transfer Are Involved in the Dynamic Transformation of Goldfish Carassius auratus Skin Color

Goldfish are one of the most popular models for studying the genetic diversity of skin color. Transcriptome sequencing (RNA-seq) and whole genome bisulfate sequencing (WGBS) of skin tissues from the third filial (F3) cyan (CN), black (BK), and white (WH) goldfish were conducted to analyze the molecular mechanism of color transformation in fish. The RNA-seq yielded 56 Gb of clean data and 56,627 transcripts from nine skin samples. The DEGs (differentially expressed genes) were enriched in cell junction cellular components and the tight junction pathway. Ninety-five homologs of the claudin family were predicted and 16 claudins were identified in correlation with skin color transformation. WGBS yielded 1079 Gb of clean data from 15 samples. Both the DEGs and the DMRs (differentially methylated regions) in the BK_CN group were found to be enriched in cytoskeleton reorganization and vesicle trafficking. Masson staining and TEM (transmission electron microscopy) confirmed the varied distribution and processes of melanosome/melanin in skin tissues. Our results suggested that cytoskeleton reorganization, cell junction, and the vesicle trafficking system played key roles in the transfer of the melanosome/melanin, and it was the extracellular translocation rather than the biosynthesis or metabolism of the melanin process that resulted in the color transformation of cyan goldfish. The data will facilitate the understanding of the molecular mechanisms underlying dynamic skin color transformation in goldfish.

Min pig skeletal muscle response to cold stress

The increased sensitivity of pigs to ambient temperature is due to today's intensive farming. Frequent climate disasters increase the pressure on healthy pig farming. Min pigs are an indigenous pig breed in China with desirable cold resistance characteristics, and hence are ideal for obtaining cold-resistant pig breeds. Therefore, it is important to discover the molecular mechanisms that are activated in response to cold stress in the Min pig. Here, we conducted a transcriptomic analysis of the skeletal muscle of Min pigs under chronic low-temperature acclimation (group A) and acute short cold stress (group B). Cold exposure caused more genes to be upregulated. Totals of 125 and 96 differentially expressed genes (DEGs) were generated from groups A and B. Sixteen common upregulated DEGs were screened; these were concentrated in oxidative stress (SRXN1, MAFF), immune and inflammatory responses (ITPKC, AREG, MMP25, FOSL1), the nervous system (RETREG1, GADD45A, RCAN1), lipid metabolism (LRP11, LIPG, ITGA5, AMPD2), solute transport (SLC19A2, SLC28A1, SLCO4A1), and fertility (HBEGF). There were 102 and 73 genes that were specifically differentially expressed in groups A and B, respectively. The altered mRNAs were enriched in immune, endocrine, and cancer pathways. There were 186 and 91 differentially expressed lncRNAs generated from groups A and B. Analysis of the target genes suggested that they may be involved in regulating the MAPK signaling pathway for resistance to cold. The results of this study provide a comprehensive overview of cold exposure-induced transcriptional patterns in skeletal muscle of the Min pig. These results can guide future molecular studies of cold stress response in pigs for improving cold tolerance as a goal in breeding programs.

Regulating Strategies of Transcription and Alternative Splicing for Cold Tolerance Harpadon nehereus Fish

In recent years, Harpadon nehereus gradually become a dominant species with great potential for exploitation in the East China Sea, and it is worth investigating whether H. nehereus would tolerate cold stress to continue to expand into the colder northern waters. The molecular regulation level is favorable evidence to explore the cold tolerance of H. nehereus, a total of 6,650, 1,936, and 2,772 differentially expressed genes (DEGs) in transcription regulation, and 4,409, 1,250, and 2,303 differential alternative splicing genes (DASGs) in alternative splicing regulation were identified in H. nehereus at 13, 15, and 17°C, respectively, importantly, 47 genes were identified as the key candidate genes for cold tolerance in H. nehereus. In transcription regulation, up-regulated DEGs were enriched in metabolic process terms and ribosome, spliceosome pathway, etc., while down-regulated DEGs were enriched in signal transduction terms, focal adhesion, proteoglycans in cancer pathway, etc., at 13, 15, and 17°C, respectively. In alternative splicing regulation, spliceosome, mRNA surveillance pathway, etc., were significantly enriched in DASGs. In a word, H. nehereus adapts to cold environments mainly through transcription and translation, transmembrane transport, protein modification, etc., while cold stress may also induce some diseases in H. nehereus.

Assembly-free rapid differential gene expression analysis in non-model organisms using DNA-protein alignment

BACKGROUND

RNA-seq is being increasingly adopted for gene expression studies in a panoply of non-model organisms, with applications spanning the fields of agriculture, aquaculture, ecology, and environment. For organisms that lack a well-annotated reference genome or transcriptome, a conventional RNA-seq data analysis workflow requires constructing a de-novo transcriptome assembly and annotating it against a high-confidence protein database. The assembly serves as a reference for read mapping, and the annotation is necessary for functional analysis of genes found to be differentially expressed. However, assembly is computationally expensive. It is also prone to errors that impact expression analysis, especially since sequencing depth is typically much lower for expression studies than for transcript discovery.

RESULTS

We propose a shortcut, in which we obtain counts for differential expression analysis by directly aligning RNA-seq reads to the high-confidence proteome that would have been otherwise used for annotation. By avoiding assembly, we drastically cut down computational costs - the running time on a typical dataset improves from the order of tens of hours to under half an hour, and the memory requirement is reduced from the order of tens of Gbytes to tens of Mbytes. We show through experiments on simulated and real data that our pipeline not only reduces computational costs, but has higher sensitivity and precision than a typical assembly-based pipeline. A Snakemake implementation of our workflow is available at: https://bitbucket.org/project_samar/samar .

CONCLUSIONS

The flip side of RNA-seq becoming accessible to even modestly resourced labs has been that the time, labor, and infrastructure cost of bioinformatics analysis has become a bottleneck. Assembly is one such resource-hungry process, and we show here that it can be avoided for quick and easy, yet more sensitive and precise, differential gene expression analysis in non-model organisms.

Transcriptome, antioxidant enzymes and histological analysis reveal molecular mechanisms responsive to long-term cold stress in silver pomfret (Pampus argenteus)

Temperature is a major environmental factor influence fish growth, development, metabolism and physiological performance. Silver pomfret (Pampus argenteus) is an economically important fishery species, however, the molecular mechanisms responsive to long-term cold stress are still unclear. Hence, we altered water temperature from 13 °C to 8 °C, a logistic fit curve for the survival rate of P. argenteus under a gradient cold stress were thus achieved, 50% survival rate at a measured temperature of 7 °C-7.5 °C. After stimulation, the gill, liver and muscle tissues were investigated through transcriptome, antioxidant enzymes and histological observation. The results showed that antioxidant enzyme and Na⁺-k⁺ ATPase activity in gill tissue was significantly increased, tissue damage and apoptosis were observed in multi-tissues. By high-throughput sequencing, a total of 618,097,404 reads of raw data and 598,855,490 reads of clean data were obtained, containing 12,489 differently expressed genes (DEGs). KEGG pathway enrichment analysis showed that DNA replication, protein digestion and absorption, cardiac muscle contraction, adrenergic signaling in cardiomyocytes, and metabolic pathways were significantly enriched in multi-tissues. Fifteen DEGs were selected for real-time PCR (RT-qPCR) analysis, and the results were consistent with transcriptome profiling. Based on the results, we inferred that P. argenteus survived at low temperatures may be achieved by improving the ability to scavenge oxyradical substance and enhancing cell fluidity. This present study indicated that the effects of long-term cold stress on P. argenteus, which is valuable for breeding cold-tolerant P. argenteus stocks for cultivation.

Characterization and analysis of transcriptome complexity using SMRT-Seq combined with RNA-Seq for a better understanding of Acanthogobius ommaturus in response to temperature stress

Acanthogobius ommaturus is a eurythermic fish, which is widely distributed in coastal, estuarine and bay waters of China, Japan and Korea. Due to the lack of whole genomic information, full-length transcriptome of A. ommaturus was firstly generated by single molecule real-time sequencing (SMRT-seq) in this study. A total of 49,833 full-length non-redundant transcripts (FLNRTs), 2255 alternative splices, 46,856 simple sequence repeats, 5094 long non-coding RNAs and 2708 transcription factors were obtained. Additionally, FLNRTs were used as reference sequences for the following transcriptome analysis of the temperature stress (7 °C, 14 °C, 21 °C (control), 28 °C and 35 °C). GO and KEGG enrichment analysis using GSEA were performed on all genes in 10 response modules which were screened out by WGCNA. Enrichment analysis showed that protein degradation, immune response and energy metabolism play an active role in the temperature stress of A. ommaturus. The differentially expressed hub genes (DEHGs) in response modules were closely related to adhesion, vascular remodeling and disease. The results of this study provided the first systematical full-length transcriptome profile of A. ommaturus and characterized its temperature stress responses, which will serve as the foundation for further exploring the molecular mechanism of the temperature stress in fish.

The effects of population and thermal acclimation on the growth, condition and cold responsive mRNA expression of age-0 lake sturgeon (Acipenser fulvescens)

In Manitoba, Canada, wild lake sturgeon (Acipenser fulvescens) populations exist along a latitudinal gradient and are reared in hatcheries to bolster threatened populations. We reared two populations of lake sturgeon, one from each of the northern and southern ends of Manitoba and examined the effects of typical hatchery temperatures (16°C) as well as 60-day acclimation to elevated rearing temperatures (20°C) on mortality, growth and condition throughout early development. Additionally, we examined the cold shock response, which may be induced during stocking, through the hepatic mRNA expression of genes involved in the response to cold stress and homeoviscous adaptation (HSP70, HSP90a, HSP90b, CIRP and SCD). Sturgeon were sampled after 1 day and 1 week following stocking into temperatures of 8, 6 and 4°C in a controlled laboratory environment. The southern population showed lower condition and higher mortality during early life than the northern population while increased rearing temperature impacted the growth and condition of developing northern sturgeon. During the cold shock, HSP70 and HSP90a mRNA expression increased in all sturgeon treatments as stocking temperature decreased, with higher expression observed in the southern population. Expression of HSP90b, CIRP and SCD increased as stocking temperature decreased in northern sturgeon with early acclimation to 20°C. Correlation analyses indicated the strongest molecular relationships were in the expression of HSP90b, CIRP and SCD, across all treatments, with a correlation between HSP90b and body condition in northern sturgeon with early acclimation to 20°C. Together, these observations highlight the importance of population and rearing environment throughout early development and on later cellular responses induced by cold stocking temperatures.

Transcriptome analysis provides the first insight into the molecular basis of temperature plasticity in Banggai cardinalfish, Pterapogon kauderni

Banggai cardinalfish, Pterapogon kauderni, is a tropical fish listed as an endangered species by IUCN. Its distribution and survival condition are extremely limited, and the changes of living environment caused by global warming may seriously threaten its geographical distribution. In order to understand the survival temperature range and the potential mechanism of temperature plasticity of P. kauderni, transcriptome analysis was performed under five temperature conditions (18 °C, 22 °C, 26 °C, 30 °C and 34 °C). A total of 432,444,497 clean reads were obtained from the mix tissues of whole head, viscera (except intestine), and muscle. All clean data were spliced into 194,832 unigenes. Compared with 26 °C, 57, 107, 187 and 174 differentially expressed genes (DEGs) were obtained at 18 °C, 22 °C, 30 °C and 34 °C, respectively. Gene Ontology (GO) analysis showed the most highly enriched in the DEGs were cellular processes, binding, metabolic processes and biological regulation. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated circadian rhythm, protein processing in endoplasmic reticulum, influenza A and prion disease were significantly enriched. 47 genes that may be related to temperature stress were identified, such as Per1, MLP, IGFBP1, HSP70, HSP90α, HSPA4, DNAJB1, CALR. This is the first RNA-Seq study of P. kauderni. This information should be valuable for further targeted studies on temperature tolerance, thereby assisting the protection and development of P. kauderni.

Fibronectin 1B Gene Plays an Important Role in Loach Barbel Air-Breathing

Loach (Misgurnus anguillicaudatus) is well known to perform air-breathing through the posterior intestine and skin. However, we find here for the first time a unique central vascular structure in the loach barbel, with a blood–gas diffusion distance as short as that of the posterior intestine. Under acute hypoxia, the distance of loach barbels became significantly shorter. Moreover, barbel removal significantly decreased air-breathing frequency of the loach. These findings imply that the barbel is another air-breathing organ of the loach. For further investigation of loach barbel air-breathing, a transcriptome analysis of barbels with air exposure treatment was performed. A total of 2546 differentially expressed genes (DEGs) between the T-XU (air exposure) and C-XU (control) group were identified, and 13 key DEGs related to barbel air-breathing were screened out. On this foundation, sequence, expression, and location analysis results indicated an important positive role of fibronectin 1b (fn1b) in loach barbel air-breathing. We further generated an fn1b-depletion loach (MT for short) using the CRISPR/Cas9 technique. It was indicated that depletion of fn1b could weaker barbel air-breathing ability. In conclusion, due to nonlethal and regenerative characteristics, the loach barbel, a newly discovered and fn1b-related fish air-breathing organ, can be a good model for fish air-breathing research.

Genome-wide association analyses identify genotype-by-environment interactions of growth traits in Simmental cattle

Understanding genotype-by-environment interactions (G × E) is crucial to understand environmental adaptation in mammals and improve the sustainability of agricultural production. Here, we present an extensive study investigating the interaction of genome-wide SNP markers with a vast assortment of environmental variables and searching for SNPs controlling phenotypic variance (vQTL) using a large beef cattle dataset. We showed that G × E contribute 10.1%, 3.8%, and 2.8% of the phenotypic variance of birth weight, weaning weight, and yearling weight, respectively. G × E genome-wide association analysis (GWAA) detected a large number of G × E loci affecting growth traits, which the traditional GWAA did not detect, showing that functional loci may have non-additive genetic effects regardless of differences in genotypic means. Further, variance-heterogeneity GWAA detected loci enriched with G × E effects without requiring prior knowledge of the interacting environmental factors. Functional annotation and pathway analysis of G × E genes revealed biological mechanisms by which cattle respond to changes in their environment, such as neurotransmitter activity, hypoxia-induced processes, keratinization, hormone, thermogenic and immune pathways. We unraveled the relevance and complexity of the genetic basis of G × E underlying growth traits, providing new insights into how different environmental conditions interact with specific genes influencing adaptation and productivity in beef cattle and potentially across mammals.