Global identification of the genetic networks and cis-regulatory elements of the cold response in zebrafish

Metabolomic-based analysis reveals bile acid-mediated ovarian failure induced by low temperature in zebrafish

As ectotherms, fish are highly sensitive to temperature fluctuations, which can profoundly impact their reproductive cycles. In this study, we investigated the fertility and histological characteristics of zebrafish ( Danio rerio) ovaries exposed to a temperature gradient ranging from the thermopreferendum temperature of the species, 27°C, to lower temperatures of 22°C, 20°C, and 13°C over a period of two weeks. Comparative metabolomic (six biological replicates for each temperature) and transcriptomic (four biological replicates for each temperature) analyses were conducted under the four temperature conditions. Results indicated that lower temperatures inhibited oocyte development and differential metabolites were involved in steroid hormone production, antioxidant function, and lipid and protein catabolism. Disrupted reproductive hormones, increased proteolysis, and lipid degradation significantly impeded oocyte development and egg maturation. Notably, a significant increase in bile acid content was noted in the ovaries of the cold-treated fish, indicating that bile acids play a critical role in ovarian failure. Overall, these findings provide valuable insights into the mechanisms governing the reproductive response of fish to cold stress.

Non-concordant epigenetic and transcriptional responses to acute thermal stress in western mosquitofish (Gambusia affinis)

Climate change is intensifying the frequency and severity of extreme temperatures. Understanding the molecular mechanisms underlying the ability to cope with acute thermal stress is key for predicting species' responses to extreme temperature events. While many studies have focused on the individual roles of gene expression, post-transcriptional processes and epigenetic modifications in response to acute thermal stress, the relative contribution of these molecular mechanisms remains unclear. The wide range of thermal limits of western mosquitofish (Gambusia affinis) provides an opportunity to explore this interplay. Here, we quantified changes in gene expression, alternative splicing, DNA methylation and microRNA (miRNA) expression in muscle tissue dissected from mosquitofish immediately after reaching high (CTmax) or low thermal limit (CTmin). Although the numbers of genes showing expression and splicing changes in response to acute temperature stress were small, we found a possibly larger and non-redundant role of splicing compared to gene expression, with more genes being differentially spliced (DSGs) than differentially expressed (DEGs), and little overlap between DSGs and DEGs. We also identified a small proportion of CpGs showing significant methylation change (i.e. differentially methylated cytosines, DMCs) in fish at thermal limits; however, there was no overlap between DEGs and genes annotated with DMCs in both CTmax and CTmin experiments. The weak interplay between epigenetic modifications and gene expression was further supported by our discoveries of no differentially expressed miRNAs. These findings provide novel insights into the relative role of different molecular mechanisms underlying immediate responses to extreme temperatures and demonstrate non-concordant responses of epigenetic and transcriptional mechanisms to acute temperature stress.

Tissue-specific temperature dependence of RNA editing levels in zebrafish

BACKGROUND

RNA editing by adenosine deaminase acting on RNA (ADAR) occurs in all metazoans and fulfils several functions. Here, we examined effects of acclimation temperature (27 °C, 18 °C,13 °C) on editing patterns in six tissues of zebrafish (Danio rerio).

RESULTS

Sites and total amounts of editing differed among tissues. Brain showed the highest levels, followed by gill and skin. In these highly edited tissues, decreases in temperatures led to large increases in total amounts of editing and changes in specific edited sites. Gene ontology analysis showed both similarities (e.g., endoplasmic reticulum stress response) and differences in editing among tissues. The majority of edited sites were in transcripts of transposable elements and the 3'UTR regions of protein coding genes. By experimental validation, translation efficiency was directly related to extent of editing of the 3'UTR region of an mRNA.

CONCLUSIONS

RNA editing increases 3'UTR polymorphism and affects efficiency of translation. Such editing may lead to temperature-adaptive changes in the proteome through altering relative amounts of synthesis of different proteins.

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.

Mitophagy: A Bridge Linking HMGB1 and Parkinson's Disease Using Adult Zebrafish as a Model Organism

High-mobility group box 1 (HMGB1) has been implicated as a key player in two critical factors of Parkinson's disease (PD): mitochondrial dysfunction and neuroinflammation. However, the specific role of HMGB1 in PD remains elusive. We investigated the effect of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) administration on mitochondrial dysfunction and HMGB1-associated inflammatory genes as well as locomotor activity in zebrafish, aiming to elucidate the role of HMGB1 in PD. Adult zebrafish received MPTP injections, and locomotor activity was measured at 24- and 48-h post-administration. Gene expression levels related to mitophagy (fis1, pink1, and park2) and HMGB1-mediated inflammation (hmgb1, tlr4, and nfkb) were quantified through RT-qPCR analysis. Following MPTP injection, the significant increase in transcript levels of fis1, pink1, and park2 indicated notable changes in PINK1/Parkin mitophagy, while the upregulation of hmgb1, tlr4, and nfkb genes pointed to the activation of the HMGB1/TLR4/NFκB inflammatory pathway. Furthermore, MPTP-injected zebrafish exhibited decreased locomotor activity, evident through reduced distance travelled, mean speed, and increased freezing durations. HMGB1 plays a major role in cellular processes as it is involved in both the mitophagy process and functions as a pro-inflammatory protein. MPTP administration in adult zebrafish activated mitophagy and inflammatory signaling, highlighting the significant role of HMGB1 as a mediator in both processes and further emphasizing its significant contribution to PD pathogenesis.

Transcriptomic Networks Reveal the Tissue-Specific Cold Shock Responses in Japanese Flounder (Paralichthys olivaceus)

Low temperature is among the important factors affecting the distribution, survival, growth, and physiology of aquatic animals. In this study, coordinated transcriptomic responses to 10 °C acute cold stress were investigated in the gills, hearts, livers, and spleens of Japanese flounder (Paralichthys olivaceus), an important aquaculture species in east Asia. Histological examination suggested different levels of injury among P. olivaceus tissues after cold shock, mainly in the gills and livers. Based on transcriptome and weighted gene coexpression network analysis, 10 tissue-specific cold responsive modules (CRMs) were identified, revealing a cascade of cellular responses to cold stress. Specifically, five upregulated CRMs were enriched with induced differentially expressed genes (DEGs), mainly corresponding to the functions of "extracellular matrix", "cytoskeleton", and "oxidoreductase activity", indicating the induced cellular response to cold shock. The "cell cycle/division" and "DNA complex" functions were enriched in the downregulated CRMs for all four tissues, which comprised inhibited DEGs, suggesting that even with tissue-specific responses, cold shock may induce severely disrupted cellular functions in all tissues, reducing aquaculture productivity. Therefore, our results revealed the tissue-specific regulation of the cellular response to low-temperature stress, which warrants further investigation and provides more comprehensive insights for the conservation and cultivation of P. olivaceus in cold water.

Dusp1 regulates thermal tolerance limits in zebrafish by maintaining mitochondrial integrity

Temperature tolerance restricts the distribution of a species. However, the molecular and cellular mechanisms that set the thermal tolerance limits of an organism are poorly understood. Here, we report on the function of dual-specificity phosphatase 1 (DUSP1) in thermal tolerance regulation. Notably, we found that dusp1 ^-/- zebrafish grew normally but survived within a narrowed temperature range. The higher susceptibility of these mutant fish to both cold and heat challenges was attributed to accelerated cell death caused by aggravated mitochondrial dysfunction and over-production of reactive oxygen species in the gills. The DUSP1-MAPK-DRP1 axis was identified as a key pathway regulating these processes in both fish and human cells. These observations suggest that DUSP1 may play a role in maintaining mitochondrial integrity and redox homeostasis. We therefore propose that maintenance of cellular redox homeostasis may be a key mechanism for coping with cellular thermal stress and that the interplay between signaling pathways regulating redox homeostasis in the most thermosensitive tissue (i.e., gills) may play an important role in setting the thermal tolerance limit of zebrafish.

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.

Prediction and Experimental Validation of a New Salinity-Responsive Cis-Regulatory Element (CRE) in a Tilapia Cell Line

Transcriptional regulation is a major mechanism by which organisms integrate gene x environment interactions. It can be achieved by coordinated interplay between cis-regulatory elements (CREs) and transcription factors (TFs). Euryhaline tilapia (Oreochromis mossambicus) tolerate a wide range of salinity and thus are an appropriate model to examine transcriptional regulatory mechanisms during salinity stress in fish. Quantitative proteomics in combination with the transcription inhibitor actinomycin D revealed 19 proteins that are transcriptionally upregulated by hyperosmolality in tilapia brain (OmB) cells. We searched the extended proximal promoter up to intron1 of each corresponding gene for common motifs using motif discovery tools. The top-ranked motif identified (STREME1) represents a binding site for the Forkhead box TF L1 (FoxL1). STREME1 function during hyperosmolality was experimentally validated by choosing two of the 19 genes, chloride intracellular channel 2 (clic2) and uridine phosphorylase 1 (upp1), that are enriched in STREME1 in their extended promoters. Transcriptional induction of these genes during hyperosmolality requires STREME1, as evidenced by motif mutagenesis. We conclude that STREME1 represents a new functional CRE that contributes to gene x environment interactions during salinity stress in tilapia. Moreover, our results indicate that FoxL1 family TFs are contribute to hyperosmotic induction of genes in euryhaline fish.

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.

Scale Development-Related Genes Identified by Transcriptome Analysis

Scales, as key structures of fish skin, play an important role in physiological function. The study of fish scale development mechanisms provides a basis for exploring the molecular-level developmental differences between scaled and non-scaled fishes. In this study, alizarin red staining was used to divide the different stages of zebrafish (Danio rerio) scale development. Four developmental stages, namely stage I (~17 dpf, scales have not started to grow), stage II (~33 dpf, the point at which scales start to grow), stage III (~41 dpf, the period in which the scales almost cover the whole body), and stage IV (~3 mpf, scales cover the whole body), were determined and used for subsequent transcriptome analysis. WGCNA (weighted correlation network analysis) and DEG (differentially expressed gene) analysis were used for screening the key genes. Based on the comparison between stage II and stage I, 54 hub-genes were identified by WGCNA analysis. Key genes including the Scpp family (Scpp7, Scpp6, Scpp5, and Scpp8), the Fgf family (Fgfr1b and Fgfr3), Tcf7, Wnt10b, Runx2b, and Il2rb were identified by DEG analysis, which indicated that these genes played important roles in the key nodes of scale development signal pathways. Combined with this analysis, the TGF-β, Wnt/β-catenin, and FGF signaling pathways were suggested to be the most important signal pathways for scales starting to grow. This study laid a foundation for exploring the scale development mechanism of other fishes. The scale development candidate genes identified in the current study will facilitate functional gene identifications in the future.

Alternative splicing of medaka bcl6aa and its repression by Prdm1a and Prdm1b

Bcl6 and Prdm1 (Blimp1) are a pair of transcriptional factors that repressing each other in mammals. Prdm1 represses the expression of bcl6 by binding a cis-element of the bcl6 gene in mammals. The homologs of Bcl6 and Prdm1 have been identified in teleost fish. However, whether these two factors regulate each other in the same way in fish like that in mammals is not clear. In this study, the regulation of bcl6aa by Prdm1 was investigated in medaka. The mRNA of bcl6aa has three variants (bcl6aaX1-X3) at the 5'-end by alternative splicing detected by RT-PCR. The three variants can be detected in adult tissues and developing embryos of medaka. Prdm1a and prdm1b are expressed in the tissues and embryos where and when bcl6aa is expressed. The expression of prdm1a was high while the expression of bcl6aa was low, and vice versa, detected in the spleen after stimulation with LPS or polyI:C. In vitro reporter assay indicated that bcl6aa could be directly repressed by both Prdm1a and Prdm1b in a dosage-dependent manner. After mutation of the key base, G, of all predicted binding sites in the core promoter region of bcl6aa, the repression by Prdm1a and/or Prdm1b disappeared. The binding site of Prdm1 in the bcl6aa gene is GAAAA(T/G). These results indicate that both Prdm1a and Prdm1b directly repress the expression of bcl6aa by binding their binding sites where the 5'-G is critical in medaka fish.

Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9

Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.

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.

metaRE R Package for Meta-Analysis of Transcriptome Data to Identify the cis-Regulatory Code behind the Transcriptional Reprogramming

At the molecular level, response to an external factor or an internal condition causes reprogramming of temporal and spatial transcription. When an organism undergoes physiological and/or morphological changes, several signaling pathways are activated simultaneously. Examples of such complex reactions are the response to temperature changes, dehydration, various biologically active substances, and others. A significant part of the regulatory ensemble in such complex reactions remains unidentified. We developed metaRE, an R package for the systematic search for cis-regulatory elements enriched in the promoters of the genes significantly changed their transcription in a complex reaction. metaRE mines multiple expression profiling datasets generated to test the same organism’s response and identifies simple and composite cis-regulatory elements systematically associated with differential expression of genes. Here, we showed metaRE performance for the identification of low-temperature-responsive cis-regulatory code in Arabidopsis thaliana and Danio rerio. MetaRE identified potential binding sites for known as well as unknown cold response regulators. A notable part of cis-elements was found in both searches discovering great conservation in low-temperature responses between plants and animals.

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.

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.

Deficiency in the membrane protein Tmbim3a/Grinaa initiates cold-induced ER stress and cell death by activating an intrinsic apoptotic pathway in zebrafish

Most members of the family of proteins containing a transmembrane BAX inhibitor motif (TMBIM) have anti-apoptotic activity, but their in vivo functions and intracellular mechanisms remain obscure. Here, we report that zebrafish Tmbim3a/Grinaa functions in the prevention of cold-induced endoplasmic reticulum (ER) stress and apoptosis. Using a gene-trapping approach, we obtained a mutant zebrafish line in which the expression of the tmbim3a/grinaa gene is disrupted by a Tol2 transposon insertion. Homozygous tmbim3a/grinaa mutant larvae exhibited time-dependently increased mortality and apoptosis under cold exposure (at 16 °C). Mechanistically, using immunofluorescence, fluorescence-based assessments of intracellular/mitochondrial Ca²⁺ levels, mitochondrial membrane potential measurements, and Ca²⁺-ATPase assays, we found that cold exposure suppresses sarcoplasmic/ER Ca²⁺-ATPase (SERCA) activity and induces the unfolded protein response (UPR) and ER stress. We also found that the cold-induced ER stress is increased in homozygous tmbim3a/grinaa mutant embryos. The cold-stress hypersensitivity of the tmbim3a/grinaa mutants was tightly associated with disrupted intracellular Ca²⁺ homeostasis, followed by mitochondrial Ca²⁺ overload and cytochrome c release, leading to the activation of caspase 9- and caspase-3-mediated intrinsic apoptotic pathways. Treatment of zebrafish larvae with the intracellular Ca²⁺ chelator 1,2-bis(2-aminophenoxy)ethane-N,N,N',N'-tetraacetate-acetoxymethyl ester (BAPTA-AM) or with 2-aminoethoxydiphenyl borate (2-APB), an inhibitor of the calcium-releasing protein inositol 1,4,5-trisphosphate receptor (IP₃R), alleviated cold-induced cell death. Together, these findings unveil a key role of Tmbim3a/Grinaa in relieving cold-induced ER stress and in protecting cells against caspase 9- and caspase 3-mediated apoptosis during zebrafish development.

Bcl6aa and bcl6ab are ubiquitously expressed and are inducible by lipopolysaccharide and polyI:C in adult tissues of medaka Oryzias latipes

B-cell lymphoma-6 (Bcl6) is a transcriptional repressor that plays important roles in various physiological activities such as innate and adaptive immune response, lymphocyte differentiation, and cell cycle regulation in mammals. Two homologs of Bcl6a, namely Bcl6aa and Bcl6ab, are identified in teleost fish including medaka Oryzias latipes. The expression profiles of bcl6aa and bcl6ab in medaka were studied using reverse-transcription polymerase chain reaction and in situ hybridization. The transcripts of bcl6aa and bcl6ab were detected from very early embryos such as the four-cell stage until hatching. Bcl6aa and bcl6ab were clearly detected in the embryonic body from 5 days postfertilization onward by in situ hybridization. Bcl6aa was specifically expressed in the retina, whereas bcl6ab was expressed in entire embryonic body. The results referred to that both bcl6aa and bcl6ab originate maternally in the zygotes and may play major roles in embryogenesis of medaka. The transcripts of bcl6aa and bcl6ab were detected in all examined adult tissues, including immune organs such as the gill, spleen, kidney, liver, and intestine. The expression of bcl6aa and bcl6ab in the liver, spleen, head-kidney, and intestine could be upregulated or downregulated by lipopolysaccharide and polyriboinosinic-polyribocytidylic acid. These results indicate that both bcl6aa and bcl6ab may be involved in immune response in medaka.

Transcriptomic responses to low temperature stress in the Nile tilapia, Oreochromis niloticus

The Nile tilapia, Oreochromis niloticus, is a species of high economic value and extensively cultured. The limited stress tolerance of this species to a low temperature usually leads to mass mortality and great loss. Nevertheless, there is limited information on the molecular mechanisms underlying the susceptibility to low temperature in the tilapia. In this study, tilapia was treated at 28 °C to a lethal temperature of 8 °C by a gradual decrement. Transcriptomic response of the immune organ, kidney, in tilapia was characterized using RNA-seq. In total, 2191 genes were annotated for significant expression, which were mainly associated with metabolism and immunity. Pathway analysis showed that immune-related pathways of phagosome and cell adhesion molecules (CAMs) pathway were significantly down-regulated under low temperature. Moreover, ferroptosis, a significantly changed pathway involved in tissue damage and acute renal failure, is reported here for the first time. The levels of serum parameters associated with kidney damage such as urea and uric acid (UA) increased significantly under low temperature. The immunofluorescence staining of the kidney showed that cell apoptosis occurred at low temperature. The results of the present study indicate that exposure to low temperature can cause kidney disfunction and down-regulate the immune-related pathway in the kidney of tilapia. This study provides new insight into the mechanism of kidney damage in fish under low temperature.

Physiological responses to cold stress in the gills of discus fish (Symphysodon aequifasciatus) revealed by conventional biochemical assays and GC-TOF-MS metabolomics

Discus fish (Symphysodon aequifasciatus) is a cichlid that is among the most popular fish for warm-water aquaria and also frequently used as the model animal for environmental science. However, little is known about the responses of S. aequifasciatus to low temperatures caused by environmental variation. Here, by using conventional biochemical assays and gas chromatography time-of-flight mass spectrometry metabolomics, we investigated the physiological responses of S. aequifasciatus gills exposed for 30 days to two temperature regimes: 28 °C and 20 °C. Low temperature resulted in elevated production of reactive oxygen species but not increased malondialdehyde. This might be partially related to protective responses in the antioxidant system, revealed by increased activities of superoxide dismutase and glutathione peroxidase, and level of reduced glutathione (GSH), compensating for the depletion of catalase activity. A total of 35 metabolites were identified as potential biomarkers of cold stress, showing the most influenced pathways including starch and sucrose metabolism, pentose phosphate pathway, glycerolipid metabolism, sphingolipid metabolism, glutathione metabolism, and arginine and proline metabolism. Moreover, the activation of glutathione metabolism agreed with the increased GSH level detected by biochemical assays. Overall, the results of this study suggest that low temperature can activate a protective antioxidant defence response and modify the metabolic pathways in gills of S. aequifasciatus, providing insights into the physiological regulation in response to cold stress in this tropical fish.

Spatiotemporal control of zebrafish (Danio rerio) gene expression using a light-activated CRISPR activation system

CRISPR activation (CRISPRa) system is the convenient tool for targeted-gene activation, it has been developed and combined with a lighting-based system that can control transcription initiation spatially and temporally by utilizing photoreceptor derived from plant Arabidopsis thaliana. A blue light photoreceptor the Cryptochrome 2 (CRY2), and its binding partner CIB1 will dimerize by exposure to the blue light and it has been applied to human cells. However, the application of a combination of these two systems to zebrafish cell is still not explored. We performed zebrafish gene activation using p65 and VP64 activators in the zebrafish cells (ZF4). Our study demonstrated that we have successfully controlled the transcription level of ASCL1a, BCL6a, and HSP70 genes using blue light-activated CRISPR activation system. The result showed that using this system, mRNA level expression of ASCL1a, BCL6a, and HSP70 genes increased after irradiated under blue light for several hours and significantly different to those which treated in the dark.

Under control: The innate immunity of fish from the inhibitors' perspective

The innate immune response involves a concerted network of induced gene products, preformed immune effectors, biochemical signalling cascades and specialised cells. However, the multifaceted activation of these defensive measures can derail or overshoot and, if left unchecked, overwhelm the host. A plenty of regulatory devices therefore mediate the fragile equilibrium between pathogen defence and pathophysiological manifestations. Over the past decade in particular, an almost complete set of teleostean sequences orthologous to mammalian immunoregulatory factors has been identified in various fish species, which prove the remarkable conservation of innate immune-control concepts among vertebrates. This review will present the current knowledge on more than 50 teleostean regulatory factors (plus additional fish-specific paralogs) that are of paramount importance for controlling the clotting cascade, the complement system, pattern-recognition pathways and cytokine-signalling networks. A special focus lies on those immunoregulatory features that have emerged as potential biomarker genes in transcriptome-wide research studies. Moreover, we report on the latest progress in elucidating control elements that act directly with immune-gene-encoding nucleic acids, such as transcription factors, hormone receptors and micro- and long noncoding RNAs. Investigations into the function of teleostean inhibitory factors are still mainly based on gene-expression profiling or overexpression studies. However, in support of structural and in-vitro analyses, evidence from in-vivo trials is also available and revealed many biochemical details on piscine immune regulation. The presence of multiple gene copies in fish adds a degree of complexity, as it is so far hardly understood if they might play distinct roles during inflammation. The present review addresses this and other open questions that should be tackled by fish immunologists in future.

Characterization of lncRNAs involved in cold acclimation of zebrafish ZF4 cells

Long non-coding RNAs (lncRNAs) are increasingly regarded as a key role in regulating diverse biological processes in various tissues and species. Although the cold responsive lncRNAs have been reported in plants, no data is available on screening and functional prediction of lncRNAs in cold acclimation in fish so far. Here we compared the expression profile of lncRNAs in cold acclimated zebrafish embryonic fibroblast cells (ZF4) cultured at 18°C for 30 days with that of cells cultured at 28°C as control by high-throughput sequencing. Totally 8,363 novel lncRNAs were identified. Including known and novel lncRNAs, there are 347 lncRNAs up-regulated and 342 lncRNAs down-regulated in cold acclimated cells. Among the differentially expressed lncRNAs, 74 and 61 were detected only in control cells or cold-acclimated cells, respectively. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses of adjacent genes to the differentially expressed lncRNAs showed that the enriched genes are involved in electron transport, cell adhesion, oxidation-reduction process, and so on. We also predicted the target genes of the differentially expressed lncRNAs by looking for interactions between lncRNAs and mRNAs, and constructed an interaction network. In summary, our genome-wide systematic identification and functional prediction of cold responsive lncRNAs in zebrafish cells suggests a crucial role of lincRNAs in cold acclimation in fish.

Temporally distinct transcriptional regulation of myocyte dedifferentiation and Myofiber growth during muscle regeneration

BACKGROUND

Tissue regeneration requires a series of steps, beginning with generation of the necessary cell mass, followed by cell migration into damaged area, and ending with differentiation and integration with surrounding tissues. Temporal regulation of these steps lies at the heart of the regenerative process, yet its basis is not well understood. The ability of zebrafish to dedifferentiate mature "post-mitotic" myocytes into proliferating myoblasts that in turn regenerate lost muscle tissue provides an opportunity to probe the molecular mechanisms of regeneration.

RESULTS

Following subtotal excision of adult zebrafish lateral rectus muscle, dedifferentiating residual myocytes were collected at two time points prior to cell cycle reentry and compared to uninjured muscles using RNA-seq. Functional annotation (GAGE or K-means clustering followed by GO enrichment) revealed a coordinated response encompassing epigenetic regulation of transcription, RNA processing, and DNA replication and repair, along with protein degradation and translation that would rewire the cellular proteome and metabolome. Selected candidate genes were phenotypically validated in vivo by morpholino knockdown. Rapidly induced gene products, such as the Polycomb group factors Ezh2 and Suz12a, were necessary for both efficient dedifferentiation (i.e. cell reprogramming leading to cell cycle reentry) and complete anatomic regeneration. In contrast, the late activated gene fibronectin was important for efficient anatomic muscle regeneration but not for the early step of myocyte cell cycle reentry.

CONCLUSIONS

Reprogramming of a "post-mitotic" myocyte into a dedifferentiated myoblast requires a complex coordinated effort that reshapes the cellular proteome and rewires metabolic pathways mediated by heritable yet nuanced epigenetic alterations and molecular switches, including transcription factors and non-coding RNAs. Our studies show that temporal regulation of gene expression is programmatically linked to distinct steps in the regeneration process, with immediate early expression driving dedifferentiation and reprogramming, and later expression facilitating anatomical regeneration.

Transcriptomic responses to environmental change in fishes: Insights from RNA sequencing

The need to better understand how plasticity and evolution affect organismal responses to environmental variability is paramount in the face of global climate change. The potential for using RNA sequencing (RNA-seq) to study complex responses by non-model organisms to the environment is evident in a rapidly growing body of literature. This is particularly true of fishes for which research has been motivated by their ecological importance, socioeconomic value, and increased use as model species for medical and genetic research. Here, we review studies that have used RNA-seq to study transcriptomic responses to continuous abiotic variables to which fishes have likely evolved a response and that are predicted to be affected by climate change (e.g., salinity, temperature, dissolved oxygen concentration, and pH). Field and laboratory experiments demonstrate the potential for individuals to respond plastically to short- and long-term environmental stress and reveal molecular mechanisms underlying developmental and transgenerational plasticity, as well as adaptation to different environmental regimes. We discuss experimental, analytical, and conceptual issues that have arisen from this work and suggest avenues for future study.

Genetic Adaptation of Schizothoracine Fish to the Phased Uplifting of the Qinghai-Tibetan Plateau

Many species of Schizothoracine, a subfamily of Cyprinidae, are highly endemic to the Qinghai–Tibetan Plateau (QTP). To characterize the adaptive changes associated with the Schizothoracine expansion at high altitudes, we sequenced tissue transcriptomes of two highland and two subhighland Schizothoracines and analyzed gene evolution patterns by comparing with lowland cyprinids. Phylogenetic tree reconstruction and divergence time estimation indicated that the common ancestor of Schizothoracine fish lived ∼32.7 million years ago (MYA), coinciding with the timing of the first phase of QTP uplifting. Both high- and subhigh-Schizothoracines demonstrated elevated dN/dS ratios in the protein-coding genes compared to lowland cyprinids, from which some biological processes implicated in altitude adaptation were commonly identified. On the other hand, the highland and subhighland lineages presented drastically divergent landscapes of positively selected genes (PSGs), enriched with very different gene ontology (GO) profiles, including those in “sensory organ morphogenesis,” “regulation of protein ubiquitination,” “blood circulation,” and “blood vessel development.” These results indicated different selection pressures imposed on the highland and subhighland lineages of the Schizothoracine subfamily, with a higher number of genes in the high-altitude species involved in adaptations such as sensory perception, blood circulation, and protein metabolism. Our study indicated divergent genetic adaptations in the aquatic species facing the phased uplifting of QTP.

Transcriptome comparison reveals a genetic network regulating the lower temperature limit in fish

Transcriptional plasticity is a major driver of phenotypic differences between species. The lower temperature limit (LTL), namely the lower end of survival temperature, is an important trait delimiting the geographical distribution of a species, however, the genetic mechanisms are poorly understood. We investigated the inter-species transcriptional diversification in cold responses between zebrafish Danio rerio and tilapia Oreochromis niloticus, which were reared at a common temperature (28 °C) but have distinct LTLs. We identified significant expressional divergence between the two species in the orthologous genes from gills when the temperature cooled to the LTL of tilapia (8 °C). Five KEGG pathways were found sequentially over-represented in the zebrafish/tilapia divergently expressed genes in the duration (12 hour) of 8 °C exposure, forming a signaling cascade from metabolic regulation to apoptosis via FoxO signaling. Consistently, we found differential progression of apoptosis in the gills of the two species in which zebrafish manifested a delayed and milder apoptotic phenotype than tilapia, corresponding with a lower LTL of zebrafish. We identified diverged expression in 25 apoptosis-related transcription factors between the two species which forms an interacting network with diverged factors involving the FoxO signaling and metabolic regulation. We propose a genetic network which regulates LTL in fishes.

A σE-Mediated Temperature Gauge Controls a Switch from LuxR-Mediated Virulence Gene Expression to Thermal Stress Adaptation in Vibrio alginolyticus

In vibrios, the expression of virulence factors is often controlled by LuxR, the master quorum-sensing regulator. Here, we investigate the interplay between LuxR and σ^E, an alternative sigma factor, during the control of virulence-related gene expression and adaptations to temperature elevations in the zoonotic pathogen Vibrio alginolyticus. An rpoE null V. alginolyticus mutant was unable to adapt to various stresses and was survival-deficient in fish. In wild type V. alginolyticus, the expression of LuxR-regulated virulence factors increased as the temperature was increased from 22°C to 37°C, but mutants lacking σ^E did not respond to temperature, indicating that σ^E is critical for the temperature-dependent upregulation of virulence genes. Further analyses revealed that σ^E binds directly to -10 and -35 elements in the luxR promoter that drive its transcription. ChIP assays showed that σ^E binds to the promoter regions of luxR, rpoH and rpoE at high temperatures (e.g., 30°C and 37°C). However, at higher temperatures (42°C) that induce thermal stress, σ^E binding to the luxR promoter decreased, while its binding to the rpoH and rpoE promoters was unchanged. Thus, the temperature-dependent binding of σ^E to distinct promoters appears to underlie a σ^E-controlled switch between the expression of virulence genes and adaptation to thermal stress. This study illustrates how a conserved temperature response mechanism integrates into quorum-sensing circuits to regulate both virulence and stress adaptation.

Zoonotic Vibrio outbreaks are believed to be closely associated with increases in environmental temperature. The mechanisms underlying this phenomenon have not been defined. Here, we show that the expression of the V. alginolyticus exotoxin Asp and other quorum-sensing (QS)-regulated virulence factors are induced by increasing temperatures, with the maximum expression observed at approximately 37°C. σ^E plays an essential role in regulating the QS master regulator LuxR in response to temperature shifts by binding directly to the -10 and -35 regions of the luxR promoter to drive its transcription. However, at higher thermal stress temperatures, σ^E binding to the luxR promoter decreased, resulting in a reduction in luxR transcription. This change underlies a binomial switch mechanism that regulates σ^E-controlled virulence gene expression patterns. Furthermore, we found that anti-σ^E signaling was involved in this stress and virulence reciprocal switch. This study suggests that a common temperature response mechanism is integrated into QS circuits to regulate both virulence and adaptation in related Vibrio taxa.