Transcriptomic characterization of cold acclimation in larval zebrafish

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.

The effect of astaxanthin on the alkalinity stress resistance of Exopalaemon carinicauda

Astaxanthin (Axn), a feed additive, can improve growth performance and enhance the environmental stress tolerance of shrimp at all growth stages. High carbonate alkalinity is considered a major stressor that affects the survival, growth, and reproduction of aquatic animals in saline-alkaline waters. In this study, a combined analysis of physiology, transcriptomics, and metabolomics was performed to explore the effected mechanism of Axn on Exopalaemon carinicauda (E. carinicauda) under alkalinity stress. The results revealed that dietary Axn can inhibit oxidative stress damage caused by alkalinity stress and maintain the normal cell structure and mitochondrial membrane potential. Transcriptomic data indicated that differentially expressed genes (DEGs) under alkalinity stress and those under alkalinity stress after Axn feeding were associated with apoptosis. The metabolic data suggested that alkalinity stress has adverse effects on ammonia metabolism, unsaturated fatty acid metabolism, and TCA cycle, and dietary Axn can improve the metabolic processes in E. carinicauda. In addition, transcriptomics and metabolomics analyses showed that Axn could help maintain the cytoskeletal structure and inhibit apoptosis under alkalinity stress; a TUNEL assay further confirmed these effects. Lastly, metabolic responses to alkalinity stress included changes in multiple amino acids and unsaturated fatty acids, and pathways related to energy metabolism were downregulated in the hepatopancreas of E. carinicauda under alkalinity stress. Collectively, all these results provide new insights into the molecular mechanisms underlying alkalinity stress tolerance in E. carinicauda after Axn feeding.

A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation

Identifying the genetic basis of local adaptation and fitness trade-offs across environments is a central goal of evolutionary biology. Cold acclimation is an adaptive plastic response for surviving seasonal freezing, and costs of acclimation may be a general mechanism for fitness trade-offs across environments in temperate zone species. Starting with locally adapted ecotypes of Arabidopsis thaliana from Italy and Sweden, we examined the fitness consequences of a naturally occurring functional polymorphism in CBF2. This gene encodes a transcription factor that is a major regulator of cold-acclimated freezing tolerance and resides within a locus responsible for a genetic trade-off for long-term mean fitness. We estimated the consequences of alternate genotypes of CBF2 on 5-y mean fitness and fitness components at the native field sites by comparing near-isogenic lines with alternate genotypes of CBF2 to their genetic background ecotypes. The effects of CBF2 were validated at the nucleotide level using gene-edited lines in the native genetic backgrounds grown in simulated parental environments. The foreign CBF2 genotype in the local genetic background reduced long-term mean fitness in Sweden by more than 10%, primarily via effects on survival. In Italy, fitness was reduced by more than 20%, primarily via effects on fecundity. At both sites, the effects were temporally variable and much stronger in some years. The gene-edited lines confirmed that CBF2 encodes the causal variant underlying this genetic trade-off. Additionally, we demonstrated a substantial fitness cost of cold acclimation, which has broad implications for potential maladaptive responses to climate change.

2Gs and plant architecture: breaking grain yield ceiling through breeding approaches for next wave of revolution in rice (Oryza sativa L.)

Rice is a principal food crop for more than half of the global population. Grain number and grain weight (2Gs) are the two complex traits controlled by several quantitative trait loci (QTLs) and are considered the most critical components for yield enhancement in rice. Novel molecular biology and QTL mapping strategies can be utilized in dissecting the complex genetic architecture of these traits. Discovering the valuable genes/QTLs associated with 2Gs traits hidden in the rice genome and utilizing them in breeding programs may bring a revolution in rice production. Furthermore, the positional cloning and functional characterization of identified genes and QTLs may aid in understanding the molecular mechanisms underlying the 2Gs traits. In addition, knowledge of modern genomic tools aids the understanding of the nature of plant and panicle architecture, which enhances their photosynthetic activity. Rice researchers continue to combine important yield component traits (including 2Gs for the yield ceiling) by utilizing modern breeding tools, such as marker-assisted selection (MAS), haplotype-based breeding, and allele mining. Physical co-localization of GW7 (for grain weight) and DEP2 (for grain number) genes present on chromosome 7 revealed the possibility of simultaneous introgression of these two genes, if desirable allelic variants were found in the single donor parent. This review article will reveal the genetic nature of 2Gs traits and use this knowledge to break the yield ceiling by using different breeding and biotechnological tools, which will sustain the world's food requirements.

Regulation of PGC-1α of the Mitochondrial Energy Metabolism Pathway in the Gills of Indian Medaka (Oryzias dancena) under Hypothermal Stress

Ectothermic fish exposure to hypothermal stress requires adjusting their metabolic molecular machinery, which was investigated using Indian medaka (Oryzias dancena; 10 weeks old, 2.5 ± 0.5 cm) cultured in fresh water (FW) and seawater (SW; 35‱) at room temperature (28 ± 1 °C). The fish were fed twice a day, once in the morning and once in the evening, and the photoperiod was 12 h:12 h light: dark. In this study, we applied two hypothermal treatments to reveal the mechanisms of energy metabolism via pgc-1α regulation in the gills of Indian medaka; cold-stress (18 °C) and cold-tolerance (extreme cold; 15 °C). The branchial ATP content was significantly higher in the cold-stress group, but not in the cold-tolerance group. In FW- and SW-acclimated medaka, the expression of genes related to mitochondrial energy metabolism, including pgc-1α, prc, Nrf2, tfam, and nd5, was analyzed to illustrate differential responses of mitochondrial energy metabolism to cold-stress and cold-tolerance environments. When exposed to cold-stress, the relative mRNA expression of pgc-1α, prc, and Nrf2 increased from 2 h, whereas that of tfam and nd5 increased significantly from 168 h. When exposed to a cold-tolerant environment, prc was significantly upregulated at 2 h post-cooling in the FW and SW groups, and pgc-1α was significantly upregulated at 2 and 12 h post-cooling in the FW group, while tfam and nd5 were downregulated in both FW and SW fish. Hierarchical clustering revealed gene interactions in the cold-stress group, which promoted diverse mitochondrial energy adaptations, causing an increase in ATP production. However, the cold-tolerant group demonstrated limitations in enhancing ATP levels through mitochondrial regulation via the PGC-1α energy metabolism pathway. These findings suggest that ectothermic fish may develop varying degrees of thermal tolerance over time in response to climate change. This study provides insights into the complex ways in which fish adjust their metabolism when exposed to cold stress, contributing to our knowledge of how they adapt.

Physiological, metabolomic, and transcriptomic reveal metabolic pathway alterations in Gymnocypris przewalskii due to cold exposure

Teleost fish have evolved various adaptations that allow them to tolerate cold water conditions. However, the underlying mechanism of this adaptation is poorly understood in Tibetan Plateau fish. RNA-seq combined with liquid chromatography‒mass spectrometry (LC‒MS/MS) metabolomics was used to investigate the physiological responses of a Tibetan Plateau-specific teleost, Gymnocypris przewalskii, under cold conditions. The 8-month G. przewalskii juvenile fish were exposed to cold (4 ℃, cold acclimation, CA) and warm (17 ℃, normal temperature, NT) temperature water for 15 days. Then, the transcript profiles of eight tissues, including the brain, gill, heart, intestine, hepatopancreas, kidney, muscle, and skin, were evaluated by transcriptome sequencing. The metabolites of the intestine, hepatopancreas, and muscle were identified by LC‒MS/MS. A total of 5,745 differentially expressed genes (DEGs) were obtained in the CA group. The key DEGs were annotated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. The DEGs from the eight tissues were significantly enriched in spliceosome pathways, indicating that activated alternative splicing is a critical biological process that occurs in the tissues to help fish cope with cold stress. Additionally, 82, 97, and 66 differentially expressed metabolites were identified in the intestine, hepatopancreas, and muscle, respectively. Glutathione metabolism was the only overlapping significant pathway between the transcriptome and metabolome analyses in these three tissues, indicating that an activated antioxidative process was triggered during cold stress. In combination with the multitissue transcriptome and metabolome, we established a physiology-gene‒metabolite interaction network related to energy metabolism during cold stress and found that gluconeogenesis and long-chain fatty acid metabolism played critical roles in glucose homeostasis and energy supply.

The role of antifreeze genes in the tolerance of cold stress in the Nile tilapia (Oreochromis niloticus)

BACKGROUND

Tilapia is one of the most essential farmed fishes in the world. It is a tropical and subtropical freshwater fish well adapted to warm water but sensitive to cold weather. Extreme cold weather could cause severe stress and mass mortalities in tilapia. The present study was carried out to investigate the effects of cold stress on the up-regulation of antifreeze protein (AFP) genes in Nile tilapia (Oreochromis niloticus). Two treatment groups of fish were investigated (5 replicates of 15 fish for each group in fibreglass tanks/70 L each): 1) a control group; the fish were acclimated to lab conditions for two weeks and the water temperature was maintained at 25 °C during the whole experimental period with feeding on a commercial diet (30% crude protein). 2) Cold stress group; the same conditions as the control group except for the temperature. Initially, the temperature was decreased by one degree every 12 h. The fish started showing death symptoms when the water temperature reached 6-8 °C. In this stage the tissue (muscle) samples were taken from both groups. The immune response of fish exposed to cold stress was detected and characterized using Differential Display-PCR (DD-PCR).

RESULTS

The results indicated that nine different up-regulation genes were detected in the cold-stressed fish compared to the control group. These genes are Integrin-alpha-2 (ITGA-2), Gap junction gamma-1 protein-like (GJC1), WD repeat-containing protein 59 isoform X2 (WDRP59), NUAK family SNF1-like kinase, G-protein coupled receptor-176 (GPR-176), Actin cytoskeleton-regulatory complex protein pan1-like (PAN-1), Whirlin protein (WHRN), Suppressor of tumorigenicity 7 protein isoform X2 (ST7P) and ATP-binding cassette sub-family A member 1-like isoform X2 (ABCA1). The antifreeze gene type-II amplification using a specific PCR product of 600 bp, followed by cloning and sequencing analysis revealed that the identified gene is antifreeze type-II, with similarity ranging from 70 to 95%. The in-vitro transcribed gene induced an antifreeze protein with a molecular size of 22 kDa. The antifreeze gene, ITGA-2 and the WD repeat protein belong to the lectin family (sugar-protein).

CONCLUSIONS

In conclusion, under cold stress, Nile tilapia express many defence genes, an antifreeze gene consisting of one open reading frame of approximately 0.6 kbp.

Delayed First Feeding Chronically Impairs Larval Fish Growth Performance, Hepatic Lipid Metabolism, and Visceral Lipid Deposition at the Mouth-Opening Stage

During the mouth-opening stage, fish larvae are susceptible to delayed first feeding (DFF). In this study, we explored the effects of DFF for two days on later growth and energy metabolism in larval fish. Results showed that DFF chronically impaired larval growth performance, thereby reducing the efficiency of feed utilization by larvae. In DFF larvae, the mRNA levels of growth inhibitors (i.e., igfbp1a and igfbp1b) were significantly upregulated and consistently maintained at high expression levels, which may be an important attribution of larval growth retardation. Concomitantly, DFF retarded the growth of adipose tissue and reduced lipid deposition in larval viscera, suggesting lipid metabolism is disordered in DFF larvae and generates inefficient lipid reserves. In the liver, we observed that DFF resulted in a significant accumulation of neutral lipids, and this phenotype did not disappear rapidly after DFF larvae received exogenous nutrition. As to the transcript analyses, we found that the expression of genes related to hepatic lipid synthesis (e.g., srebf1, srebf2, dgat1a, dgat1b, fasn, and scdb) in DFF larvae was consistently upregulated, while the expression of genes involved in lipid transport (e.g., apoa2, apoa4b.1, and apoa4b.3) was downregulated. Therefore, it appears that the inefficient lipid reserves in DFF larvae are associated with their hepatic lipid transport dysfunction. Taken together, our findings contribute to understanding the impairments to fish larvae caused by delayed first feeding during the mouth-opening stage and to aiding larval management in the aquaculture industry.

Flower transcriptional response to long term hot and cold environments in Antirrhinum majus

Short term experiments have identified heat shock and cold response elements in many biological systems. However, the effect of long-term low or high temperatures is not well documented. To address this gap, we grew Antirrhinum majus plants from two-weeks old until maturity under control (normal) (22/16°C), cold (15/5°C), and hot (30/23°C) conditions for a period of two years. Flower size, petal anthocyanin content and pollen viability obtained higher values in cold conditions, decreasing in middle and high temperatures. Leaf chlorophyll content was higher in cold conditions and stable in control and hot temperatures, while pedicel length increased under hot conditions. The control conditions were optimal for scent emission and seed production. Scent complexity was low in cold temperatures. The transcriptomic analysis of mature flowers, followed by gene enrichment analysis and CNET plot visualization, showed two groups of genes. One group comprised genes controlling the affected traits, and a second group appeared as long-term adaptation to non-optimal temperatures. These included hypoxia, unsaturated fatty acid metabolism, ribosomal proteins, carboxylic acid, sugar and organic ion transport, or protein folding. We found a differential expression of floral organ identity functions, supporting the flower size data. Pollinator-related traits such as scent and color followed opposite trends, indicating an equilibrium for rendering the organs for pollination attractive under changing climate conditions. Prolonged heat or cold cause structural adaptations in protein synthesis and folding, membrane composition, and transport. Thus, adaptations to cope with non-optimal temperatures occur in basic cellular processes.

Rising floor and dropping ceiling: organ heterogeneity in response to cold acclimation of the largest extant amphibian

Low temperature imposes strong selective pressure on ectotherms. To maximize their overall fitness under cold conditions, the ectotherms may either try to maintain their physiological activities through metabolic compensation or enter into metabolic depression; however, some species adopt both strategies to cope with different degrees of cold. Nevertheless, how these two seemingly opposite strategies are coordinated has rarely been elucidated. Here, we investigated the molecular strategy underlying the cold acclimation of Andrias davidianus, the largest extant amphibian, using multi-organ metabolomics and transcriptomics. The results showed remarkable organ heterogeneity in response to cold. While most organs showed transcriptional upregulation of metabolic processes, the heart exhibited downregulation. This heterogeneity explained the adaptive reorganization in resource allocation, which compensates for metabolic maintenance by compromising growth. Importantly, the cardiac function might constitute a ‘ceiling’ to constrain the space for compensation, especially under colder conditions. Additionally, the opposite transcriptional regulation of oxidative phosphorylation and other pathways might also shape the overall metabolic capacity under cold conditions. The heterogeneity in cold responses may have directed a shift in cold adaptive strategy from compensation to depression with a drop in temperature. These results provide a novel insight into the regulatory mechanisms underlying cold survival strategies of ectotherms.

Understanding the Function and Mechanism of Zebrafish Tmem39b in Regulating Cold Resistance

Autophagy and endoplasmic reticulum (ER) stress response are among the key pathways regulating cold resistance of fish through eliminating damaged cellular components and facilitating the restoration of cell homeostasis upon exposure to acute cold stress. The transmembrane protein 39A (TMEM39A) was reported to regulate both autophagy and ER stress response, but its vertebrate-specific paralog, the transmembrane protein 39B (TMEM39B), has not been characterized. In the current study, we generate tmem39b-knockout zebrafish lines and characterize their survival ability under acute cold stress. We observed that the dysfunction of Tmem39b remarkably decreased the cold resilience of both the larval and adult zebrafish. Gene transcription in the larvae exposed to cold stress and rewarming were characterized by RNA sequencing (RNA-seq) to explore the mechanisms underlying functions of Tmem39b in regulating cold resistance. The results indicate that the deficiency of Tmem39b attenuates the up-regulation of both cold- and rewarming-induced genes. The cold-induced transcription factor genes bif1.2, fosab, and egr1, and the rewarming-activated immune genes c3a.3, il11a, and sting1 are the representatives influenced by Tmem39b dysfunction. However, the loss of tmem39b has little effect on the transcription of the ER stress response- and autophagy-related genes. The measurements of the phosphorylated H2A histone family member X (at Ser 139, abbreviated as γH2AX) demonstrate that zebrafish Tmem39b protects the cells against DNA damage caused by exposure to the cold-warming stress and facilitates tissue damage repair during the recovery phase. The gene modules underlying the functions of Tmem39b in zebrafish are highly enriched in biological processes associated with immune response. The dysfunction of Tmem39b also attenuates the up-regulation of tissue C-reactive protein (CRP) content upon rewarming. Together, our data shed new light on the function and mechanism of Tmem39b in regulating the cold resistance of fish.

Mild cold stress specifically disturbs clustering movement of DFCs and sequential organ left-right patterning in zebrafish

In poikilothermic animals, the distinct acclimatization ability of different organs has been previously addressed, while the tissue-specific role of cold stress in early development is largely unknown. In this study, we discovered that despite its role in delaying embryonic development, mild cold stress (22°C) does not disturb multiple-organ progenitor specification, but does give rise to organ left-right (LR) patterning defects. Regarding the mechanism, the data showed that mild cold stress downregulated the expression of cell-adhesion genes cdh1 and cdh2 during gastrulation, especially in dorsal forerunner cells (DFCs), which partially disturbed the clustering movement of DFCs, Kupffer’s vesicle (KV) morphogenesis, and ciliogenesis. As a result, the defects of KV/cilia disrupted asymmetric nodal signaling and subsequent heart and liver LR patterning. In conclusion, our data novelly identified that, in early development, DFCs are more sensitive to mild cold stress, and mild cold stress repressed the expression of cell adhesion-related gene cdh1 and cdh2. This role partially disturbed the clustering movement of DFCs, which resulted in defective KV/cilia development and sequential organ LR patterning defects.

Warming Accelerates the Onset of the Molecular Stress Response and Increases Mortality of Larval Atlantic Cod

Temperature profoundly affects ectotherm physiology. Although differential thermal responses influence fitness, thus driving population dynamics and species distributions, our understanding of the molecular architecture underlying these responses is limited, especially during the critical larval stage. Here, using RNA-sequencing of laboratory-reared Atlantic cod (Gadus morhua) larvae of wild origin, we find changes in gene expression in thousands of transcripts consistent with a severe cellular stress response at both ambient and projected (+2°C and +4°C) temperatures. In addition, specific responses to stress, heat, and hypoxia were commonly identified in gene ontology enrichment analyses and 33 of the 44 genes comprising the minimum stress proteome of all organisms were upregulated. Earlier onset of the stress response was evident at higher temperatures; concomitant increased growth and mortality suggests a reduction in fitness. Temporal differences in gene expression levels do not correspond to differences in growing degree days, suggesting negative physiological consequences of warming beyond accelerated development. Because gene expression is costly, we infer that the upregulation of thousands of transcripts in response to warming in larval cod might act as an energetic drain. We hypothesize that the energetically costly stress response, coupled with increased growth rate at warmer temperatures, leads to faster depletion of energy reserves and increased risk of mortality in larval cod. As sea surface temperatures continue to rise over the next century, reduced fitness of Atlantic cod larvae might lead to population declines in this ecologically and socioeconomically important species. Further, our findings expand our understanding of transcriptomic responses to temperature by ectothermic vertebrate larvae beyond the critical first-feeding stage, a time when organisms begin balancing the energetic demands of growth, foraging, development, and maintenance. Linking the molecular basis of a thermal response to key fitness-related traits is fundamentally important to predicting how global warming will affect ectotherms.

Transcriptome analysis reveals the early resistance of zebrafish larvae to oxidative stress

Oxidative stress is one of most common environmental stresses encountered by fish, especially during their fragile larval stage. More and more studies are aimed at understanding the antioxidant defense mechanism of fish larvae. Herein we characterized the early resistance of zebrafish larvae to oxidative stress and investigated the underlying transcriptional regulations using RNA-seq. We found that pre-exposure of zebrafish larvae to 2 mM H₂O₂ for 1 or 3 h significantly improved their survival under higher doses of H₂O₂ (3 mM), suggesting the antioxidant defenses of zebrafish larvae were rapidly built under pre-exposure of H₂O₂. Comparative transcriptome analysis showed that 310 (185 up and 125 down) and 512 (331 up and 181 down) differentially expressed genes were generated after 1 and 3 h of pre-exposure, respectively. KEGG enrichment analysis revealed that protein processing in endoplasmic reticulum is a highly enriched pathway; multiple genes (e.g., hsp70.1, hsp70.2, and hsp90aa1.2) encoding heat shock proteins in this pathway were sharply upregulated presumably to correct protein misfolding and maintaining the cellular normal functions during oxidative stress. More importantly, the Keap1/Nrf2 system-mediated detoxification enzyme system was significantly activated, which regulates the upregulation of target genes (e.g., gstp1, gsr, and prdx1) to scavenger reactive oxygen species, thereby defending against apoptosis. In addition, the MAPK, as a transmitter of stress signals, was activated, which may play an important role in activating antioxidant system in the early stages of oxidative stress. Altogether, these findings demonstrate that zebrafish larvae rapidly establish resistance to oxidative stress, and this involves changes in protein processing, stress signal transmission, and the activation of detoxification pathways.

Tissue-Specific and Differential Cold Responses in the Domesticated Cold Tolerant Fugu

Domestication can be defined as the artificial selection in animals to achieve morphological, physiological, and developmental conformity to human needs, with the aim of improving various limitations in species under a human feeding environment. The future sustainability of aquaculture may rely partly on the availability of numerous domesticated fish species. However, the underlying adaptive mechanisms that result in the domestication of fish are still unclear. Because they are poikilothermic, temperature is a key environmental element that affects the entire life of fish, so studying the association between physiological and behavioral changes in low-temperature domesticated fish can provide a model for understanding the response mechanisms of fish under cold stress. Through 5 generations and 10 years of artificial selection at low temperatures, we used cold-tolerant fugu as a biological model to compare transcriptome changes in brain and liver tissues to study the effects of cold stress on fish. It was found that the expression of genes such as apoptosis, p53, oxidative phosphorylation, and mitochondrial β-oxidation in the brain of cold-tolerant fugu was significantly lower than the wild type due to cold stress, while excessive energy metabolism would lead to the production of reactive oxygen species (ROS) and exacerbate the brain damage, thus causing rollover and coma. Meanwhile, under cold stress, the signaling pathways involved in glycogenolysis and lipid metabolism, such as insulin signaling, adipocytokines, and mTOR signaling pathways, were significantly up-regulated in the liver of cold-tolerant fugu. Although the mitochondrial β-oxidation pathway was increased in cold-tolerant fugu liver tissues, the transcriptome was not enriched in apoptotic. These phenomena predict that in response to low-temperature conditions, cold-tolerant fugu employs a dynamic inter-organ metabolic regulation strategy to cope with cold stress and reduce damage to brain tissues.

A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies

Organisms residing in regions with alternating seasons often develop different phenotypes, or forms, in each season. These forms are often adaptations to each season and result from an altered developmental response to specific environmental cues such as temperature. While multiple studies have examined form-specific gene expression profiles in a diversity of species, little is known about how environments and developmental transitions, cued by hormone pulses, alter post-transcriptional patterns. In this study, we examine how gene expression, alternative splicing, and miRNA-mediated gene silencing in Bicyclus anynana butterfly hindwing tissue, varies across two rearing temperatures at four developmental timepoints. These timepoints flank two temperature-sensitive periods that coincide with two pulses of the insect hormone 20E. Our results suggest that developmental transitions, coincident with 20E pulses, elicit a greater impact on all these transcriptomic patterns than rearing temperatures per se. More similar transcriptomic patterns are observed pre-20E pulses than those observed post-20E pulses. We also found functionally distinct sets of differentially expressed and differentially spliced genes in the seasonal forms. Furthermore, around 10% of differentially expressed genes are predicted to be direct targets of, and regulated by, differentially expressed miRNAs between the seasonal forms. Many differentially expressed genes, miRNAs, or differentially spliced genes potentially regulate eyespot size plasticity, and we validated the differential splicing pattern of one such gene, daughterless. We present a comprehensive and interactive transcriptomic atlas of the hindwing tissue of both seasonal forms of B. anynana throughout development, a model organism of seasonal plasticity.

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.

An updated review of cold shock and cold stress in fish

Temperature is critical in regulating virtually all biological functions in fish. Low temperature stress (cold shock/stress) is an often-overlooked challenge that many fish face as a result of both natural events and anthropogenic activities. In this study, we present an updated review of the cold shock literature based on a comprehensive literature search, following an initial review on the subject by M.R. Donaldson and colleagues, published in a 2008 volume of this journal. We focus on how knowledge on cold shock and fish has evolved over the past decade, describing advances in the understanding of the generalized stress response in fish under cold stress, what metrics may be used to quantify cold stress and what knowledge gaps remain to be addressed in future research. We also describe the relevance of cold shock as it pertains to environmental managers, policymakers and industry professionals, including practical applications of cold shock. Although substantial progress has been made in addressing some of the knowledge gaps identified a decade ago, other topics (e.g., population-level effects and interactions between primary, secondary and tertiary stress responses) have received little or no attention despite their significance to fish biology and thermal stress. Approaches using combinations of primary, secondary and tertiary stress responses are crucial as a research priority to better understand the mechanisms underlying cold shock responses, from short-term physiological changes to individual- and population-level effects, thereby providing researchers with better means of quantifying cold shock in laboratory and field settings.

Antioxidant, Transcriptome and the Metabolome Response to Dietary Astaxanthin in Exopalaemon carinicauda

Astaxanthin (Axn), a feed additive, is becoming increasingly important for modulating the metabolism, growth, development, and reproduction of aquatic organisms in aquaculture. In this study, Exopalaemon carinicauda (E. carinicauda) is an economically important fishery species in China that has been found to exhibit increased body weight following Axn feeding as compared to a standard diet. The antioxidant, transcriptomic, and metabolomic analyses of the response of E. carinicauda after Axn feeding were investigated. Axn could reduce the content of malondialdehyde and increase the activities of various antioxidant enzymes, which also proved that axn can improve the antioxidant capacity Transcriptomic analysis suggested that synthesis and secretion of immune proteins, cytoskeleton structure, and apoptosis signaling were altered after Axn feeding. The metabolic response to axn mainly includes the up regulation of different amino acids and the change of unsaturated fatty acids. Combined transcriptomic and metabolomic data indicated that amino acid metabolic pathways were upregulated in the muscles after Axn feeding. For good measure, energy metabolism pathways were upregulated in the muscles to improve ATP and unsaturated fatty acid production. This study provides key information to increase our understanding of the effects of Axn in shrimp.

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.

Low Temperature Effect on the Endocrine and Circadian Systems of Adult Danio rerio

The control of the biological rhythms begins with the activation of photo- and thermosensitive cells located in various organs of the fish such as brain, eye, and skin, but a central clock is still to be identified in teleosts. Thermal changes are stressors which increase cortisol and affect the rhythm of other hormones such as melatonin and growth hormone (GH), in both endo- and ectothermic organisms. Our aim was to investigate how temperature (23°C for 6 days) lower than the optimal (28°C) modulates expression of several gene pathways including growth hormone (gh1) and its receptors (ghra, ghrb), insulin-like growth factor1 (igf1a, igf1b) and its receptors (igf1ra, igf1rb), cortisol and its receptor (gr), the limiting enzyme of melatonin synthesis (arylalkylamine N-acetyltransferase, aanat) and melatonin receptors (mtnr1aa, mtnr1bb), as well as their relationship with clock genes in Danio rerio in early light and early dark phases of the day. Lower temperature reduced the expression of the hormone gene gh1, and of the related receptors ghra, ghrb, igf1ra, and igf1rb. Cortisol levels were higher at the lower temperature, with a decrease of its receptor (gr) transcripts in the liver. Interestingly, we found higher levels of aanat transcripts in the brain at 23°C. Overall, lower temperature downregulated the transcription of hormone related genes and clock genes. The results suggest a strong correlation of temperature challenge with the clock molecular mechanism and the endocrine systems analyzed, especially the growth hormone and melatonin axes, in D. rerio tissues.

Identification of Genes Related to Cold Tolerance and Novel Genetic Markers for Molecular Breeding in Taiwan Tilapia (Oreochromis spp.) via Transcriptome Analysis

Taiwan tilapia is one of the primary species used in aquaculture practices in Taiwan. However, as a tropical fish, it is sensitive to cold temperatures that can lead to high mortality rates during winter months. Genetic and broodstock management strategies using marker-assisted selection and breeding are the best tools currently available to improve seed varieties for tilapia species. The purpose of this study was to develop molecular markers for cold stress-related genes using digital gene expression analysis of next-generation transcriptome sequencing in Taiwan tilapia (Oreochromis spp.). We constructed and sequenced cDNA libraries from the brain, gill, liver, and muscle tissues of cold-tolerance (CT) and cold-sensitivity (CS) strains. Approximately 35,214,833,100 nucleotides of raw sequencing reads were generated, and these were assembled into 128,147 unigenes possessing a total length of 185,382,926 bp and an average length of 1446 bp. A total of 25,844 unigenes were annotated using five protein databases and Venny analysis, and 38,377 simple sequence repeats (SSRs) and 65,527 single nucleotide polymorphisms (SNPs) were identified. Furthermore, from the 38-cold tolerance-related genes that were identified using differential gene expression analysis in the four tissues, 13 microsatellites and 37 single nucleotide polymorphism markers were identified. The results of the genotype analysis revealed that the selected markers could be used for population genetics. In addition to the diversity assessment, one of the SNP markers was determined to be significantly related to cold-tolerance traits and could be used as a molecular marker to assist in the selection and verification of cold-tolerant populations. The specific genetic markers explored in this study can be used for the identification of genetic polymorphisms and cold tolerance traits in Taiwan tilapia, and they can also be used to further explore the physiological and biochemical molecular regulation pathways of fish that are involved in their tolerance to environmental temperature stress.

Cloning and characterisation of a Δ9 fatty acyl desaturase-like gene from the red claw crayfish (Cherax quadricarinatus) and its expression analysis under cold stress

Desaturase is one of the key enzymes in the unsaturated fatty acid synthesis pathway. Δ9 desaturase catalyzes the synthesis of oleic acid from stearic acid by introducing double bonds in the 9th and 10th carbon chains, thereby increasing the content of MUFAs in the body. In order to explore the main function of the Δ9 desaturase gene under low temperature stress, RACE-PCR technology was used in this study to clone the full-length sequence of the CqFAD9-like from the hepatopancreas of red claw crayfish, Cherax quadricarinatus. The full length of the sequence is 1236 bp, and the open reading frame is 1041 bp, encoding 346 amino acid residues. The 5 'UTR is 116 bp, the 3' UTR is 79 bp, and the 3 'UTR contains a PloyA tail. The predicted theoretical isoelectric point and molecular weight are 8.68 and 40.28 kDa, respectively. Homology analysis showed that the sequence had the highest similarity with FAD9 from crustaceans. The results of real-time PCR showed that the expression level of this gene was highest in the hepatopancreas, which was significantly higher than other tissues, followed by the ovaries, brain ganglion and stomach. At the same time, the expression of the CqFAD9-like in hepatopancreas of crayfish cultured at 25, 20, 15 and 9 °C for four weeks was detected. The results showed that expression of the FAD9 gene increased gradually with decreasing temperature, indicating that metabolic desaturation might play a regulatory role during cold stress.

Transcriptional response to cold and fasting acclimation in Onychostoma macrolepis during the overwintering stage

In this study, we investigated the transcriptome responses of the liver of Onychostoma macrolepis in by RNA sequencing. The sampling process involved three groups: 1G (0 week, 10 °C), 2G (12 weeks, 0 °C) and 3G (24 weeks, 10 °C). The body weight, viscera index, hepatopancreas index and intraperitoneal fat index of O. macrolepis showed a decreasing trend with the prolonging of overwintering time. The crude fat contents of whole fish, muscle and liver in O. macrolepis after overwintering were significantly lower than those of the fish before overwintering (p < 0.05). In 1G versus 2G group, 2G versus 3G group and 1G versus 3G group, the differently expressed genes (DEGs) were 4630, 3976 and 2311, respectively. These results indicated that different stages of overwintering period had significant effects on gene expression of O. macrolepis, and the influence degree gradually decreased with the extension of overwintering period. The results of Gene ontology (GO) enrichment showed that these DEGs were mainly related to metabolism and immunity, and most of them were down-regulated. In this study, the KEGG pathway classification results showed that signal transduction was the most representative. In addition, KOG enrichment results showed that many DEGs associated with lipid transport and metabolism were down-regulated during the overwintering period. These observations suggested that slowing metabolism and delaying immunity may be the strategies for overwintering adaptation of O. macrolepis.

Different ways to play it cool: Transcriptomic analysis sheds light on different activity patterns of three amphipod species under long-term cold exposure

Species of littoral freshwater environments in regions with continental climate experience pronounced seasonal temperature changes. Coping with long cold winters and hot summers requires specific physiological and behavioural adaptations. Endemic amphipods of Lake Baikal, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, show high metabolic activity throughout the year; E. verrucosus even reproduces in winter. In contrast, the widespread Holarctic amphipod Gammarus lacustris overwinters in torpor. This study investigated the transcriptomic hallmarks of E. verrucosus, E. cyaneus and G. lacustris exposed to low water temperatures. Amphipods were exposed to 1.5°C and 12°C (corresponding to the mean winter and summer water temperatures, respectively, in the Baikal littoral) for one month. At 1.5°C, G. lacustris showed upregulation of ribosome biogenesis and mRNA processing genes, as well as downregulation of genes related to growth, reproduction and locomotor activity, indicating enhanced energy allocation to somatic maintenance. Our results suggest that the mitogen-activated protein kinase (MAPK) signalling pathway is involved in the preparation for hibernation; downregulation of the actin cytoskeleton pathway genes could relate to the observed low locomotor activity of G. lacustris at 1.5°C. The differences between the transcriptomes of E. verrucosus and E. cyaneus from the 1.5°C and 12°C exposures were considerably smaller than for G. lacustris. In E. verrucosus, cold-exposure triggered reproductive activity was indicated by upregulation of respective genes, whereas in E. cyaneus, genes related to mitochondria functioning were upregulated, indicating cold compensation in this species. Our data elucidate the molecular characteristics behind the different adaptations of amphipod species from the Lake Baikal area to winter conditions.

Improving water quality does not guarantee fish health: Effects of ammonia pollution on the behaviour of wild-caught pre-exposed fish

Ammonia is a pollutant frequently found in aquatic ecosystems. In fish, ammonia can cause physical damage, alter its behaviour, and even cause death. Exposure to ammonia also increases fish physiological stress, which can be measured through biomarkers. In this study, we analysed the effect of sublethal ammonia concentrations on the behaviour and the oxidative stress of Barbus meridionalis that had been pre-exposed to this compound in the wild. Wild-caught fish from a polluted site (pre-exposed fish) and from an unpolluted site (non-pre-exposed fish) were exposed, under experimental conditions, to total ammonia concentrations (TAN) of 0, 1, 5, and 8 mg/L. Swimming activity, feeding behaviour, and oxidative stress response based on biomarkers were analysed. Pre-exposed fish showed both an altered behaviour and an altered oxidative stress response in the control treatment (0 mg/L). Differences in swimming activity were also found as pre-exposed fish swam less. Lower feeding activity (voracity and satiety) and altered response to oxidative stress were also observed at ≥ 1 mg/L TAN. Biomarker results confirmed pre-exposed fish suffer from a reduction in their antioxidant defences and, hence, showed increased oxidative tissue damage. In summary, pre-exposed fish showed more sensitivity to ammonia exposure than fish from a pristine site.

Gene Expression Profile and Co-Expression Network of Pearl Gentian Grouper under Cold Stress by Integrating Illumina and PacBio Sequences

Simple Summary

In this study, we investigated the liver transcriptomic responses of pearl gentian grouper towards cold stress. Some cold-related key genes and biological pathways were screened, of which energy-related metabolic pathways and genes had higher expression levels under cold stress. This suggested that energy homeostasis plays a crucial role in the physiological adjustments of pearl gentian grouper when exposed to the cold stress environment. Our results will expedite the understanding of different fishes adaptive mechanisms to profound environmental temperature changes and provide insights into the molecular breeding of cold-tolerant pearl gentian grouper varieties.

Abstract

Pearl gentian grouper (Epinephelus fuscoguttatus ♀ × Epinephelus lanceolatus ♂) is a fish of high commercial value in the aquaculture industry in Asia. However, this hybrid fish is not cold-tolerant, and its molecular regulation mechanism underlying cold stress remains largely elusive. This study thus investigated the liver transcriptomic responses of pearl gentian grouper by comparing the gene expression of cold stress groups (20, 15, 12, and 12 °C for 6 h) with that of control group (25 °C) using PacBio SMRT-Seq and Illumina RNA-Seq technologies. In SMRT-Seq analysis, a total of 11,033 full-length transcripts were generated and used as reference sequences for further RNA-Seq analysis. In RNA-Seq analysis, 3271 differentially expressed genes (DEGs), two low-temperature specific modules (tan and blue modules), and two significantly expressed gene sets (profiles 0 and 19) were screened by differential expression analysis, weighted gene co-expression networks analysis (WGCNA), and short time-series expression miner (STEM), respectively. The intersection of the above analyses further revealed some key genes, such as PCK, ALDOB, FBP, G6pC, CPT1A, PPARα, SOCS3, PPP1CC, CYP2J, HMGCR, CDKN1B, and GADD45Bc. These genes were significantly enriched in carbohydrate metabolism, lipid metabolism, signal transduction, and endocrine system pathways. All these pathways were linked to biological functions relevant to cold adaptation, such as energy metabolism, stress-induced cell membrane changes, and transduction of stress signals. Taken together, our study explores an overall and complex regulation network of the functional genes in the liver of pearl gentian grouper, which could benefit the species in preventing damage caused by cold stress.

Generation of GCaMP6s-Expressing Zebrafish to Monitor Spatiotemporal Dynamics of Calcium Signaling Elicited by Heat Stress

The ability of organisms to quickly sense and transduce signals of environmental stresses is critical for their survival. Ca²⁺ is a versatile intracellular messenger involved in sensing a wide variety of stresses and regulating the subsequent cellular responses. So far, our understanding for calcium signaling was mostly obtained from ex vivo tissues and cultured cell lines, and the in vivo spatiotemporal dynamics of stress-triggered calcium signaling in a vertebrate remains to be characterized. Here, we describe the generation and characterization of a transgenic zebrafish line with ubiquitous expression of GCaMP6s, a genetically encoded calcium indicator (GECI). We developed a method to investigate the spatiotemporal patterns of Ca²⁺ events induced by heat stress. Exposure to heat stress elicited immediate and transient calcium signaling in developing zebrafish. Cells extensively distributed in the integument of the head and body trunk were the first batch of responders and different cell populations demonstrated distinct response patterns upon heat stress. Activity of the heat stress-induced calcium signaling peaked at 30 s and swiftly decreased to near the basal level at 120 s after the beginning of exposure. Inhibition of the heat-induced calcium signaling by LaCl₃ and capsazepine and treatment with the inhibitors for CaMKII (Ca²²/calmodulin-dependent protein kinase II) and HSF1 (Heat shock factor 1) all significantly depressed the enhanced heat shock response (HSR). Together, we delineated the spatiotemporal dynamics of heat-induced calcium signaling and confirmed functions of the Ca²⁺-CaMKII-HSF1 pathway in regulating the HSR in zebrafish.

Characterization of Biological Pathways Regulating Acute Cold Resistance of Zebrafish

Low temperature stress represents a major threat to the lives of both farmed and wild fish species. However, biological pathways determining the development of cold resistance in fish remain largely unknown. Zebrafish larvae at 96 hpf were exposed to lethal cold stress (10 °C) for different time periods to evaluate the adverse effects at organism, tissue and cell levels. Time series RNA sequencing (RNA-seq) experiments were performed to delineate the transcriptomic landscape of zebrafish larvae under cold stress and during the subsequent rewarming phase. The genes regulated by cold stress were characterized by progressively enhanced or decreased expression, whereas the genes associated with rewarming were characterized by rapid upregulation upon return to normal temperature (28 °C). Genes such as trib3, dusp5 and otud1 were identified as the representative molecular markers of cold-induced damages through network analysis. Biological pathways involved in cold stress responses were mined from the transcriptomic data and their functions in regulating cold resistance were validated using specific inhibitors. The autophagy, FoxO and MAPK (mitogen-activated protein kinase) signaling pathways were revealed to be survival pathways for enhancing cold resistance, while apoptosis and necroptosis were the death pathways responsible for cold-induced mortality. Functional mechanisms of the survival-enhancing factors Foxo1, ERK (extracellular signal-regulated kinase) and p38 MAPK were further characterized by inhibiting their activities upon cold stress and analyzing gene expression though RNA-seq. These factors were demonstrated to determine the cold resistance of zebrafish through regulating apoptosis and p53 signaling pathway. These findings have provided novel insights into the stress responses elicited by lethal cold and shed new light on the molecular mechanisms underlying cold resistance of fish.

Comparative transcriptomics of spotted seatrout (Cynoscion nebulosus) populations to cold and heat stress

Resilience to climate change depends on a species' adaptive potential and phenotypic plasticity. The latter can enhance survival of individual organisms during short periods of extreme environmental perturbations, allowing genetic adaptation to take place over generations. Along the U.S. East Coast, estuarine-dependent spotted seatrout (Cynoscion nebulosus) populations span a steep temperature gradient that provides an ideal opportunity to explore the molecular basis of phenotypic plasticity. Genetically distinct spotted seatrout sampled from a northern and a southern population were exposed to acute cold and heat stress (5 biological replicates in each treatment and control group), and their transcriptomic responses were compared using RNA-sequencing (RNA-seq). The southern population showed a larger transcriptomic response to acute cold stress, whereas the northern population showed a larger transcriptomic response to acute heat stress compared with their respective population controls. Shared transcripts showing significant differences in expression levels were predominantly enriched in pathways that included metabolism, transcriptional regulation, and immune response. In response to heat stress, only the northern population significantly upregulated genes in the apoptosis pathway, which could suggest greater vulnerability to future heat waves in this population as compared to the southern population. Genes showing population-specific patterns of expression, including hpt, acot, hspa5, and hsc71, are candidates for future studies aiming to monitor intraspecific differences in temperature stress responses in spotted seatrout. Our findings contribute to the current understanding of phenotypic plasticity and provide a basis for predicting the response of a eurythermal fish species to future extreme temperatures.

Acute environmental temperature variation affects brain protein expression, anxiety and explorative behaviour in adult zebrafish

This study investigated the effect of 4-d acute thermal treatments at 18 °C, 26 °C (control) and 34 °C on the nervous system of adult zebrafish (Danio rerio) using a multidisciplinary approach based on behavioural tests and brain proteomic analysis. The behavioural variations induced by thermal treatment were investigated using five different tests, the novel tank diving, light and dark preference, social preference, mirror biting, and Y-Maze tests, which are standard paradigms specifically tailored for zebrafish to assess their anxiety-like behaviour, boldness, social preference, aggressiveness, and explorative behaviour, respectively. Proteomic data revealed that several proteins involved in energy metabolism, messenger RNA translation, protein synthesis, folding and degradation, cytoskeleton organisation and synaptic vesiculation are regulated differently at extreme temperatures. The results showed that anxiety-like behaviours increase in zebrafish at 18 °C compared to those at 26 °C or 34 °C, whereas anxiety-related protein signalling pathways are downregulated. Moreover, treatments at both 18 °C and 34 °C affect the exploratory behaviour that appears not to be modulated by past experiences, suggesting the impairment of fish cognitive abilities. This study is the continuation of our previous work on the effect of 21-d chronic treatment at the same constant temperature level and will enable the comparison of acute and chronic treatment effects on the nervous system function in adult zebrafish.

Footprints of global change in marine life: Inferring past environment based on DNA methylation and gene expression marks

Ocean global warming affects the distribution, life history and physiology of marine life. Extreme events, like marine heatwaves, are increasing in frequency and intensity. During sensitive stages of early fish development, the consequences may be long-lasting and mediated by epigenetic mechanisms. Here, we used European sea bass as a model to study the possible long-lasting effects of a marine heatwave during early development. We measured DNA methylation and gene expression in four tissues (brain, muscle, liver and testis) and detected differentially methylated regions (DMRs). Six genes were differentially expressed and contained DMRs three years after exposure to increased temperature, indicating direct phenotypic consequences and representing persistent changes. Interestingly, nine genes contained DMRs around the same genomic regions across tissues, therefore consisting of common footprints of developmental temperature in environmentally responsive loci. These loci are, to our knowledge, the first metastable epialleles (MEs) described in fish. MEs may serve as biomarkers to infer past life history events linked with persistent consequences. These results highlight the importance of subtle phenotypic changes mediated by epigenetics to extreme weather events during sensitive life stages. Also, to our knowledge, it is the first time the molecular effects of a marine heatwave during the lifetime of individuals are assessed. MEs could be used in surveillance programs aimed at determining the footprints of climate change on marine life. Our study paves the way for the identification of conserved MEs that respond equally to environmental perturbations across species. Conserved MEs would constitute a tool of assessment of global change effects in marine life at a large scale.

Inhibition of Embryonic HSP 90 Function Promotes Variation of Cold Tolerance in Zebrafish

Accumulating evidence indicates that heat shock protein 90 (HSP90) plays essential roles in modulation of phenotypic plasticity in vertebrate development, however, the roles of HSP90 in modulation of cold tolerance capacity in fish are still unclear. In the present study, we showed that transient inhibition of embryonic HSP90 function by a chemical inhibitor or low conductivity stress promoted variation of cold tolerance capacity in adult zebrafish. Further work showed that embryonic HSP90 inhibition enhanced cold tolerance in adult zebrafish could be transmitted to their offspring. RNA-seq data showed that embryonic HSP90 inhibition enhanced cold tolerance involves variation of gene expression related to proteasome, lysosome, autophagy, and ribosome. Experiments with zebrafish ZF4 cells showed that two differentially expressed genes atg9b and psmd12 were up-regulated by radicicol treatment and provided protective roles for cells under cold stress, indicating that up-regulation of autophagy and proteasome function contributes to enhanced cold tolerance. The present work sheds a light on the roles of HSP90 in regulation of phenotypic plasticity associated with thermal adaptation in fish.

Integrated transcriptomic and metabolomic responses in the hepatopancreas of kuruma shrimp (Marsupenaeus japonicus) under cold stress

In aquatic ecosystems, the temperature of the water is an important ecological factor that modulates aquatic organisms' metabolism, growth, development, and reproduction. In this study, the morphological, transcriptomic, and metabolomic analyses of response of Marsupeneus japonicus to acute cold stress was investigated. The results revealed that low temperature caused profound morphological damage to the hepatopancreas. Transcriptomic responses suggested that energy and primary metabolism, cytoskeleton structure, and apoptosis signaling were altered. The metabolic responses to cold stress included changes of multiple amino acids and unsaturated fatty acids. Combined transcriptomic and metabolomic data indicated that energy metabolism pathways were downregulated in the hepatopancreas under cold stress. However, M. japonicus increased ATP and unsaturated fatty acids production to ameliorate. Moreover, cold stress caused significant attenuation of macrophage apoptosis. This study provides key information to increase our understanding of low-temperature tolerance in shrimp.

Gene expression changes in response to low temperatures in embryos of the kelp grouper, Epinephelus moara

The kelp grouper Epinephelus moara has high economic value and is popular in fisheries and aquaculture in China. In the previous study, we treated the embryos at 16-22 somite stage at 4 °C, -25.7 °C, -140 °C and -196 °C, and successfully obtained surviving embryos in each group. To better understand the molecular changes affected by the low temperatures, we conducted a comparative transcriptome analysis among embryos exposed at 4 °C for 30 min, embryos exposed at -25.7 °C for 30 min and the control group. qPCR assays were conducted for the validation. Signal transduction pathways were highly enriched for the differentially expressed genes. c-Fos, c-Jun, JunD, GADD45, involved in MAPK signaling pathway, were upregulated when embryos were treated at low temperatures. As immediate early genes, Egr-1a and b, and IER2, that respond quickly to the environment stress, their expression increased as well. Hsp70 showed similar expression pattern as immediate early genes. Meanwhile, transcription factors Sox, HES, TFIID, muscle movement and protein synthesis-related genes were downregulated. Taken together, our findings suggest that cooling disrupts gene expression patterns in E. moara embryos. The differentially expressed genes may be involved in cellular resistance against low temperatures, possibly through neural activation, apoptosis, proliferation, differentiation, cellular recovery and heat shock regulation. This study also provides transcriptome dataset of E. moara embryos exposed to cold temperatures for future studies focusing on the molecular effects of cryopreservation.

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

BACKGROUND

Tropical stenothermal fish exhibit special tolerance and response to cold stress. However current knowledge of the molecular mechanisms response to cold stress in aquatic ectotherms is largely drawn from eurythermal or extreme stenothermal species. The tiger barb Puntius tetrazona is a tropical stenothermal fish, with great popularity in aquarium trade and research.

RESULTS

To investigate the response mechanism of P. tetrazona to low temperature, fish were exposed to increasing levels of acute cold stress. Histopathological analysis showed that the brain, gill, liver and muscle tissues appeared serious damage after cold stress (13 °C). Brain, gill, liver and muscle tissues from control (CTRL) groups (27 °C) and COLD stress groups (13 °C) of eight-month fish (gender-neutral) were sampled and assessed for transcriptomic profiling by high-throughput sequencing. 83.0 Gb of raw data were generated, filtered and assembled for de novo transcriptome assembly. According to the transcriptome reference, we obtained 392,878 transcripts and 238,878 unigenes, of which 89.29% of the latter were annotated. There were 23,743 differently expressed genes (DEGs) been filtered from four pairs of tissues (brain, gill, liver and muscle) between these cold stress and control groups. These DEGs were mainly involved in circadian entrainment, circadian rhythm, biosynthesis of steroid and fatty acid. There were 64 shared DEGs between the four pairs of groups, and five were related to ubiquitylation/deubiquitylation. Our results suggested that ubiquitin-mediated protein degradation might be necessary for tropical stenothermal fish coping with acute cold stress. Also, the significant cold-induced expression of heat shock 70 kDa protein (HSP70) and cold-induced RNA-binding protein (CIRBP) was verified. These results suggested that the expression of the molecular chaperones HSP70 and CIRBP in P. tetrazona might play a critical role in coping with acute cold stress.

CONCLUSIONS

This is the first transcriptome analysis of P. tetrazona using RNA-Seq technology. Novel findings about tropical stenothermal fish under cold stress (such as HSP70 and CIRBP genes) are presented here. This study contributes new insights into the molecular mechanisms of tropical stenothermal species response to acute cold stress.

Transcriptional Programs Underlying Cold Acclimation of Common Carp (Cyprinus carpio L.)

Properly regulated transcriptional responses to environmental perturbations are critical for the fitness of fish. Although gene expression profiles in the tissues of common carp upon cold stress were previously characterized, the transcriptional programs underlying cold acclimation are still not well known. In this study, the ability of three common carp strains including Hebao red carp (HB), Songpu mirror carp (SPM) and Yellow river carp (YR) to establish cold resistance after acclimation to a mild hypothermia stress at 18°C for 24 h was confirmed by measurements of the critical thermal minimums (CTMin). The gene expression profiles of the brain and the heart from these strains under both control and cold-acclimated conditions were characterized with RNA-sequencing. The data of the three common carp strains with different genetic background were combined in the differential gene expression analyses to balance the effects of genetic diversity on gene expression. Marked effects of tissue origins on the cold-induced transcriptional responses were revealed by comparing the differentially expressed gene (DEG) lists of the two tissues. Functional categories including spliceosome and RNA splicing were highly enriched in the DEGs of both tissues. However, steroid biosynthesis was specifically enriched in DEGs of the brain and response to unfolded protein was solely enriched in DEGs of the heart. Consistent with the up-regulation of the genes involved in cholesterol biosynthesis, total cholesterol content of the brain was significantly increased upon cold stress. Moreover, cold-induced alternative splicing (AS) events were explored and AS of the rbmx (RNA-binding motif protein, X chromosome) gene was confirmed by real-time quantitative PCR. Finally, a core set of cold responsive genes (CRGs) were defined by comparative transcriptomic analyses. Our data provide insights into the transcriptional programs underlying cold acclimation of common carp and offer valuable clues for further investigating the genetic determinants for cold resistance of farmed fish.

Modeling Cancer Using Zebrafish Xenografts: Drawbacks for Mimicking the Human Microenvironment

The first steps towards establishing xenografts in zebrafish embryos were performed by Lee et al., 2005 and Haldi et al., 2006, paving the way for studying human cancers using this animal species. Since then, the xenograft technique has been improved in different ways, ranging from optimizing the best temperature for xenografted embryo incubation, testing different sites for injection of human tumor cells, and even developing tools to study how the host interacts with the injected cells. Nonetheless, a standard protocol for performing xenografts has not been adopted across laboratories, and further research on the temperature, microenvironment of the tumor or the cell-host interactions inside of the embryo during xenografting is still needed. As a consequence, current non-uniform conditions could be affecting experimental results in terms of cell proliferation, invasion, or metastasis; or even overestimating the effects of some chemotherapeutic drugs on xenografted cells. In this review, we highlight and raise awareness regarding the different aspects of xenografting that need to be improved in order to mimic, in a more efficient way, the human tumor microenvironment, resulting in more robust and accurate in vivo results.

Transcriptome analysis of the Larimichthys polyactis under heat and cold stress

The small yellow croaker (Larimichthys polyactis) is an economically important marine fish that is widely distributed in the East Sea, Yellow Sea and Bohai of China. However, the wild populations of L. polyactis are severely depleted, and there is currently a developing large-scale artificial propagation of this fish for aquaculture. However, the current variety of L. polyactis that is cultivated is not capable to coping with large fluctuations in temperature. Therefore, it is important to understand the molecular mechanisms that are activated in response to temperature stress in the L. polyactis. Here, we conducted transcriptomic analysis of the liver of L. polyactis under heat and cold stress. A total of 270,844,888, 265,727,006 and 259,666,218 clean reads were generated from high temperature group, low temperature group and control group, respectively, and comparing expression of genes in these transcriptomes, 10,878 unigenes that were differential expressed were identified. Seventeen of the differentially expressed unigenes were validated by qRT-PCR. Pathway enrichment analysis identified that the ER pathway, immune signaling pathway and metabolic response pathway were affected by temperature stress. The results of this study provide a comprehensive overview of temperature stress-induced transcriptional patterns in liver tissues of the L. polyactis. In addition, these results can guide future molecular studies of heat and cold stress response in this species for improving the stock used for aquaculture.

Transcriptomic profiling revealed key signaling pathways for cold tolerance and acclimation of two carp species

Background

Closely related species of the carp family (Cyprinidae) have evolved distinctive abilities to survive under cold stress, but molecular mechanisms underlying the generation of cold resistance remain largely unknown. In this study, we compared transcriptomic profiles of two carp species to identify key factors and pathways for cold tolerance and acclimation.

Results

Larvae of Songpu mirror carp and Barbless carp that were pretreated at 18 °C for 24 h significantly improved their survival rates under lethal cold temperature at 8 °C or 10 °C, indicating that two carp species possess the ability of cold acclimation. However, Songpu mirror carp exhibited stronger abilities of cold tolerance and acclimation than Barbless carp. Transcriptomic profiles of Songpu mirror carp and Barbless carp larvae at 28 °C and 18 °C were compared during cold acclimation through RNA-seq. Differentially expressed genes that are closely associated with the differences in cold acclimation between two carp species were identified through bioinformatics and Venn’s diagram analysis. GO enrichment analysis of these genes indicated that cellular component assembly involved in morphogenesis, secondary alcohol metabolism and drug transport were the most up-regulated biological processes during cold acclimation of Songpu mirror carp. Conversely, positive regulation of macroautophagy, intracellular protein transport, and organonitrogen compound catabolism were the most down-regulated biological processes during cold acclimation of Barbless carp. KEGG enrichment analysis revealed that factors in the FoxO-related signaling pathways are mainly responsible for the development of differences in cold tolerance and acclimation between two carp species since altering the phosphorylation of key proteins in the FoxO-related signaling pathways with inhibitors or an activator significantly decreased the cold tolerance and acclimation of Songpu mirror carp. These data provided key clues for dissection of molecular mechanisms underlying the development of cold tolerance and acclimation in carps.

Conclusions

These findings indicate that larvae of two carp species possess different abilities of cold tolerance and can build cold acclimation under mild low temperature. Multiple biological processes and FoxO-related signaling pathways are closely associated with the development of differences in cold tolerance and acclimation between two carp species.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

Transcriptomic responses to pollution in natural riverine environment in Rita rita

Aquatic pollution is one of the most common threats to the ecological health of aquatic ecosystems and its biota. Fish as lower vertebrates are excellent model to study the impact and responses of aquatic pollution. In fish, gill is the main organ indicator of whole animal health as it comes in contact with the surrounding water and absorbs many pollutants and contaminants; therefore, investigations on alterations in fish gill at transcriptome level could provide newer insights to the stress response mechanism(s) and pathways. For comprehensive evaluation of the impacts of pollutants (joint toxicity) prevalent in the riverine environment, comparative transcriptome analysis, by Next Generation Sequencing under Illumina HiSeq 2500 platform, was carried out in gill tissues of Rita rita collected from two stretches of river Ganga (Kanpur and Farakka) and results were validated by RT-qPCR. Out of 154,077 unigenes (Accession SRR548008), a total of 2024 differentially expressed genes (DEGs) including 942 up-regulated and 1082 down-regulated genes were identified by DESeq program. Further, Gene Ontology (GO) of DEGs showed that ribosomal large subunit biogenesis, mitochondrial ribosome and box H/ACA SnoRNA binding categories are highly affected by pollution. Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis showed the involvement of the DEGs in energy metabolism, translational and transcriptional machinery, protein folding and degradation suggesting that these signalling pathways are highly affected by aquatic pollution. Among the DEGs, up-regulation of cytochrome c oxidase subunit (cox) 7a2 (69.47 fold), hsp70 subunit 14 (hsp70-14, 5.27 fold), muscle related coiled-coil protein (MURC, 21.55 fold), lysozyme G (40.14 fold), cox17 (29.36 fold) were the conspicuous ones which showed similar trends in expression when analysed by RT-qPCR. Based on fold change, perturbation values, correlation analysis by PCA and RT-qPCR validation, up-regulation of cox7a2, MURC and hsp70-14 appeared to be the most promising biomarker responses and could be useful in the evaluation of gill health and possibly be extended towards aquatic ecosystem health assessment.

Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism for cells to generate transcript variability and proteome diversity. No systematic investigation of AS events among different tissues in response to stressors is available for tilapia currently. In this study, AS among different tissues was identified and the cold stress-related AS events were explored in a Nile tilapia (Oreochromis niloticus) line based on 42 RNA-seq datasets using a bioinformatics pipeline. 14,796 (82.76%; SD = 2,840) of the expression genes showed AS events. The two most abundant AS types were alternative transcription start site (TSS) and terminal site (TTS) in tilapia. Testis, brain and kidney possess the most abundant AS gene number, while the blood, muscle and liver possess the least number in each tissue. Furthermore, 208 differentially alternative splicing (DAS) genes in heart and 483 DAS in brain in response to cold stress. The number of AS types for alternative exon end, exon skipping and retention of single intron increased significantly under cold stress. GO enrichment and pathway overrepresentation analysis indicated that many DAS genes, e.g., genes in circadian clock pathway, may influence expression of downstream genes under cold stress. Our study revealed that AS exists extensively in tilapia and plays an important role in cold adaption.

Identification and characterization of miRNAs involved in cold acclimation of zebrafish ZF4 cells

MicroRNAs (miRNAs) play vital roles in various biological processes under multiple stress conditions by leading to mRNA cleavage or translational repression. However, the detailed roles of miRNAs in cold acclimation in fish are still unclear. In the present study, high-throughput sequencing was performed to identify miRNAs from 6 small RNA libraries from the zebrafish embryonic fibroblast ZF4 cells under control (28°C, 30 days) and cold-acclimation (18°C, 30 days) conditions. A total of 414 miRNAs, 349 known and 65 novel, were identified. Among those miRNAs, 24 (19 known and 5 novel) were up-regulated, and 23 (9 known and 14 novel) were down-regulated in cold acclimated cells. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses indicated that the target genes of known differentially expressed miRNAs (DE-miRNA) are involved in cold acclimation by regulation of phosphorylation, cell junction, intracellular signal transduction, ECM-receptor interaction and so on. Moreover, both miR-100-3p inhibitor and miR-16b mimics could protect ZF4 cells under cold stress, indicating the involvement of miRNA in cold acclimation. Further study showed that miR-100-3p and miR-16b could regulate inversely the expression of their target gene (atad5a, cyp2ae1, lamp1, rilp, atxn7, tnika, btbd9), and that overexpression of miR-100-3p disturbed the early embryonic development of zebrafish. In summary, the present data show that miRNAs are closely involved in cold acclimation in zebrafish ZF4 cells and provide information for further understanding of the roles of miRNAs in cold acclimation in fish.

Dynamic Phosphoproteome Profiling of Zebrafish Embryonic Fibroblasts during Cold Acclimation

Temperature affects almost all aspects of the fish life. To cope with low temperature, fish have evolved the ability of cold acclimation for survival. However, intracellular signaling events underlying cold acclimation in fish remain largely unknown. Here, the formation of cold acclimation in zebrafish embryonic fibroblasts (ZF4) is monitored and the phosphorylation events during the process are investigated through a large-scale quantitative phosphoproteomic approach. In total, 11 474 phosphorylation sites are identified on 4066 proteins and quantified 5772 phosphosites on 2519 proteins. Serine, threonine, and tyrosine (Ser/Thr/Tyr) phosphorylation accounted for 85.5%, 13.3%, and 1.2% of total phosphosites, respectively. Among all phosphosites, 702 phosphosites on 510 proteins show differential regulation during cold acclimation of ZF4 cells. These phosphosites are divided into six clusters according to their dynamic changes during cold exposure. Kinase-substrate prediction reveals that mitogen-activated protein kinase (MAPK) among the kinase groups is predominantly responsible for phosphorylation of these phosphosites. The differentially regulated phosphoproteins are functionally associated with various cellular processes such as regulation of actin cytoskeleton and MAPK signaling pathway. These data enrich the database of protein phosphorylation sites in zebrafish and provide key clues for the elucidation of intracellular signaling networks during cold acclimation of fish.

The STK16 inhibitor STK16-IN-1 inhibits non-adrenergic and non-neurogenic smooth muscle contractions in the human prostate and the human male detrusor

Mixed lower urinary tract symptoms (LUTS) (voiding symptoms suggestive of benign prostatic hyperplasia plus storage symptoms, which can be caused by overactive bladder) are common in men. Unwanted contraction of prostate and/or bladder smooth muscle has been implied in the pathophysiology of male LUTS. Here, we examined effects of the serine/threonine kinase 16 (STK16) inhibitor STK16-IN-1 on contraction of human tissues from the prostate and male detrusor. Tissues were obtained from radical prostatectomy and radical cystectomy. Contractions were studied in an organ bath and STK16 expressions by Western blot analyses and fluorescence staining. In prostate tissues, STK16-IN-1 (1 μM) inhibited contractions induced by endothelin-1 and the thromboxane A₂ analog U46619. Contractions of prostate tissues induced by noradrenaline, the α₁-agonists phenylephrine and methoxamine, or electric field stimulation (EFS) were not changed by STK16-IN-1. In male detrusor tissues, STK16-IN-1 inhibited contractions induced by the cholinergic agonists carbachol and metacholine, and contractions induced by U46619. EFS-induced contractions of detrusor tissues were not changed by STK16-IN-1. Western blot analyses of prostate and detrusor tissues revealed bands matching the molecular weight of STK16. Fluorescence staining of prostate tissues using STK16 antibodies resulted in immunoreactivity in smooth muscle cells. STK16-IN-1 selectively inhibits non-adrenergic/non-neurogenic smooth muscle contractions in the male prostate and to limited extent in the bladder. Because non-adrenergic contractions in the male LUTS may account for limited efficacy of α₁-blockers and for α₁-blocker-resistant symptoms, studies assessing add-on of STK16-IN-1 to α₁-blockers in mixed LUTS appear feasible.

Transcriptome analyses provide the first insight into the molecular basis of cold tolerance in Larimichthys polyactis

Larimichthys polyactis is one of the most economically important marine fish species that have become newly cultured in China in recent years. The gene expression changes that L. polyactis experiences in cold toleranceis still unknown, limiting the expansion of its cultivation, fast growth, and high yield. To investigate the molecular mechanism behind L. polyactis's cold tolerance and to provide a resource for conducting genetic research on L. polyactis, transcriptome sequencing using RNA-seq was performed on individuals that survived cold stress at 4 °C (cold tolerant, CT), and individuals that barely survived 4 °C (cold sensitive, CS), which was considered as the control. A number of 387,607,550 clean reads were obtained from the transcriptomes, and comparative transcriptomic analysis identified 141 differently expressed genes (DEGs), of which 67 were up-regulated and 74 were down-regulated in CT compared to CS under cold stress. Furthermore, ten differently expressed genes were selected from the RNA-Seq analysis to be further validated by real-time PCR. Functional network analysis indicated that L. polyactis adapted to cold stress by employing a series of mechanisms to minimize damages caused by exposure to cold temperatures. The molecular mechanisms identified through RNA-Seq included Extracellular matrix (ECM) receptor interaction, glycerolipid metabolism, regulation of autophagy and focal adhesion pathway as playing vital roles in cold tolerance in L. polyactis. This study may help elucidate how L. polyactis tolerates cold, which is of value for breeding cold-tolerant L. polyactis stocks for cultivation.

Transcriptome analysis and weighted gene co-expression network reveals potential genes responses to heat stress in turbot Scophthalmus maximus

Turbot (Scophthalmus maximus) is an economically important marine fish cultured in China. In this study, we performed transcriptome gene expression profiling of kidney tissue in turbot exposed to heat stress (20, 23, 25 and 28 °C); control fish were maintained at 14 °C. We investigated gene relationships based on weighted gene co-expression network analysis (WGCNA). Accordingly, enrichment analyses of GO terms and KEGG pathways showed that several pathways (e.g., fat metabolism, cell apoptosis, immune system, and insulin signaling) may be involved in the response of turbot to heat stress. Moreover, via WGCNA, we identified 19 modules: the dark grey module was mainly enriched in pathways associated with fat metabolism and the FOXO and Jak-STAT signaling pathways. The ivory module was significantly enriched in the P53 signaling pathway. Furthermore, the key hub genes CBP, AKT3, CCND2, PIK3r2, SCOS3, mdm2, cyc-B, and p48 were enriched in the FOXO, Jak-STAT and P53 signaling pathways. This is the first study reporting co-expression patterns of a gene network after heat stress in marine fish. Our results may contribute to our understanding of the underlying molecular mechanism of thermal tolerance.