The effect of short-term hypoxic exposure on metabolic gene expression

Differential expression analyses and detection of SNP loci associated with environmental variables: Are salinity and temperature factors involved in population differentiation and speciation in Odontesthes?

Environmental factors play a key role in individual adaptation to different local conditions. Because of this, studies about the physiological and genetic responses of individuals exposed to different natural environments offer clues about mechanisms involved in population differentiation, and as a subsequent result, speciation. Marine environments are especially suited to survey this kind of phenomena because they commonly harbor species adapted to different local conditions along a geographic continuum. Silversides belonging to Odontesthes are commonly distributed in tropical and temperate regions of South America and exhibit noticeable phenotypic plasticity, which allows them to adapt to contrasting environments. In this study, the genetic expression of O. argentinensis sampled along the Uruguayan Atlantic coast and estuarine adjacent areas was investigated. In addition, the correlation between individual genotypes and environmental variables was also analysed in O. argentinensis and O. bonariensis. Results obtained suggest a differential expression pattern of low magnitude among individuals from the different areas sampled and a correlation between several SNP loci and environmental variables. The analyses carried out did not show a clear differentiation among individuals sampled along different salinity regimens, but enriched GOTerms seem to be driven by water oxygen content. On the other hand, a total of 46 SNPs analysed in O. argentinensis and O. bonariensis showed a correlation with salinity and temperature. Although none of the correlated SNPs and corresponding genes from our both analyses were directly associated with hypoxia, genes related to the cardiovascular system and muscle cell differentiation were found. All these genes are interesting candidates for future studies since they are closely related to the differentially expressed genes. Although salinity was also mentioned as an important parameter limiting introgression between O. argentinensis and O. bonariensis, it was found that salinity does not drive differential expression in O. argentinensis, but rather oxygen levels. Moreover, salinity does not directly affect the structure and genetic divergence of the populations, they appear to be structured based on their degree of isolation and geographical distance between them. Further studies, like genome-wide analyses, could help to elucidate additional genes adapted to the different environments in these silverside species.

The heat shock protein 20 gene family in large yellow croaker (Larimichthys crocea): Identification, phylogenetic relationships, expression analyses

Large yellow croaker (Larimichthys crocea) is an economically important fish in China, but its aquaculture industry has been threatened by both biotic and abiotic stressors such as hypoxia and pathogens. In the current study, hsp20 genes were identified and analyzed systematically for the first time from the genome of large yellow croaker, and their roles in hypoxia response and Aeromonas hydrophila, Pseudomonas plecoglossicida infection were investigated. Herein, 11 hsp20 genes were identified and annotated, phylogenetic analysis and selection pressure analysis showed that the hsp20 genes were evolutionarily-constrained and their function was conserved among fishes. Besides, we observed the expression patterns of the hsp20 genes under hypoxia and two pathogens' stress. In brief, seven, four, seven genes responded to hypoxia stress, A. hydrophila infection and P. plecoglossicida challenge, respectively, which indicated that they were involved in hypoxia and disease responses. Furthermore, pathogen- and time-specific pattern was observed after A. hydrophila and P. plecoglossicida infection whereas tissue-specific pattern was observed after hypoxia exposure, revealing that hsp20 genes showed differential functions in response to hypoxia and immune stress. Taken together, these results provided preliminary information for future analysis of the roles of hsp20 genes in both biotic and abiotic stress response in fish.

Are the effects of catch-and-release angling evident in changes to mRNA abundances related to metabolism, acid-base regulation and stress in lake trout (Salvelinus namaycush) gills?

Catch-and-release (C&R) angling is a conservation-oriented practice intended to reduce the impact recreational angling has on fish populations. Even though most recreationally angled fish are released, little is known about how C&R angling impacts fish at the cellular or tissue level. As the first to explore the impacts of C&R angling on mRNA abundances, our study aimed to identify how the stress of angling influenced metabolism, acid-base regulation and cellular stress in the gills of lake trout (Salvelinus namaycush). Because gills are responsible for metabolic gas exchange, are crucial sites of acid-base homeostasis and respond to stressors quickly, we hypothesized that the relative mRNA abundance of genes related to these three physiological processes would be altered after angling. We took gill samples of live lake trout at 0, 2 or 48 h after fish were angled by rod and reel, and then used quantitative PCR (qPCR) to measure the relative abundance of nine candidate mRNA transcripts. Heat shock protein 70 (hsp70) mRNA levels significantly increased over 5-fold 2 h after angling, indicating a potential activation of a cytoprotective response. However, contrary to our hypothesis, we observed no change in the relative mRNA abundance of genes related to metabolism or acid-base regulation in response to C&R angling within a 48-h period. As C&R angling can negatively impact fish populations, further use of transcript-level studies will allow us to understand the impact C&R has on specific tissues and improve our knowledge of how C&R influences overall fish health.

Whole-Genome Methylation Sequencing of Large Yellow Croaker (Larimichthys crocea) Liver Under Hypoxia and Acidification Stress

Large yellow croaker (Larimichthys crocea) is an important aquaculture species in China. This study analysed whole-genome methylation differences in liver tissues of young fish under different hypoxic and acidification conditions. Differentially methylated regions (DMRs) and differentially methylated genes (DMGs) were identified. Gene ontology (GO) and Kyoto encyclopaedia of genes and genomes (KEGG) enrichment analyses of DMGs were conducted to explore the mechanism of coping with hypoxic acidification. The main methylation type was CG, accounting for > 70% of total methylation, significantly higher than CHG and CHH methylation types. GO enrichment analysis of DMGs revealed strong enrichment of nervous system development, cell periphery, plasma membrane, cell junction organisation, cell junction, signalling receptor activity, molecular sensor activity, cell-linked tissue junction organisation, cell-cell adhesion and nervous system development. KEGG enrichment analysis of DMR-related genes identified cell adhesion molecules, cortisol synthesis and secretion and aldosterone synthesis and secretion as the three key pathways regulating the physiological responses to hypoxia and acidification. Long-term hypoxic and acidification stress affected the immune system, nervous system and stress responses of large yellow croaker. Whole-genome sequencing analysis of exposed tissues was used to investigate changes that occur in L. crocea in response to hypoxic and acidic conditions at the DNA methylation level. The findings contribute to our comprehensive understanding of functional methylation in large yellow croaker and will support future research on the response mechanisms of this species under different environmental pressures.

Physiological and behavioural strategies of aquatic animals living in fluctuating environments

Shallow or near-shore environments, such as ponds, estuaries and intertidal zones, are among the most physiologically challenging of all aquatic settings. Animals inhabiting these environments experience conditions that fluctuate markedly over relatively short temporal and spatial scales. Living in these habitats requires the ability to tolerate the physiological disturbances incurred by these environmental fluctuations. This tolerance is achieved through a suite of physiological and behavioural responses that allow animals to maintain homeostasis, including the ability to dynamically modulate their physiology through reversible phenotypic plasticity. However, maintaining the plasticity to adjust to some stresses in a dynamic environment may trade off with the capacity to deal with other stressors. This paper will explore studies on select fishes and invertebrates exposed to fluctuations in dissolved oxygen, salinity and pH. We assess the physiological mechanisms these species employ to achieve homeostasis, with a focus on the plasticity of their responses, and consider the resulting physiological trade-offs in function. Finally, we discuss additional factors that may influence organismal responses to fluctuating environments, such as the presence of multiple stressors, including parasites. We echo recent calls from experimental biologists to consider physiological responses to life in naturally fluctuating environments, not only because they are interesting in their own right but also because they can reveal mechanisms that may be crucial for living with increasing environmental instability as a consequence of climate change.

Integrated application of multi-omics strategies provides insights into the environmental hypoxia response in Pelteobagrus vachelli muscle

Increasing pressures on aquatic ecosystems because of pollutants, nutrient enrichment, and global warming have severely depleted oxygen concentrations. This sudden and significant lack of oxygen has resulted in persistent increases in fish mortality rates. Revealing the molecular mechanism of fish hypoxia adaptation will help researchers to find markers for hypoxia induced by environmental stress. Here, we used a multiomics approach to identify several hypoxia-associated miRNAs, mRNAs, proteins, and metabolites involved in diverse biological pathways in the muscles of Pelteobagrus vachelli. Our findings revealed significant hypoxia-associated changes in muscles over 4 h of hypoxia exposure and discrete tissue-specific patterns. We have previously reported that P. vachelli livers exhibit increased anaerobic glycolysis, heme synthesis, erythropoiesis, and inhibit apoptosis when exposed to hypoxia for 4 h. However, the opposite was observed in muscles. According to our comprehensive analysis, fishes show an acute response to hypoxia, including activation of catabolic pathways to generate more energy, reduction of biosynthesis to decrease energy consumption, and shifting from aerobic to anaerobic metabolic contributions. Also, we found that hypoxia induced muscle dysfunction by impairing mitochondrial function, activating inflammasomes, and apoptosis. The hypoxia-induced mitochondrial dysfunction enhanced oxidative stress, apoptosis, and further triggered interleukin-1β production via inflammasome activation. In turn, interleukin-1β further impaired mitochondrial function or apoptosis by suppressing downstream mitochondrial biosynthesis–related proteins, thus resulting in a vicious cycle of inflammasome activation and mitochondrial dysfunction. Our findings contribute meaningful insights into the molecular mechanisms of hypoxia, and the methods and study design can be utilized across different fish species.

•

First multiomics analysis of mRNA, miRNA, protein, and metabolite in fishes.

•

Liver and muscle were tissue-specific induced by hypoxia.

•

About 70 genes and 16 miRNAs related to hypoxia adaptation were detected.

•

Hypoxia affects muscle function by mediating energy metabolism via HIF pathway.

Identification of Candidate Genes Associated With Hypoxia Tolerance in Trachinotus blochii Using Bulked Segregant Analysis and RNA-Seq

Golden Pompano (Trachinotus blochii) has rapidly developed into the one of the main valuable fish species in Chinese marine aquaculture. Due to its rapid growth, active metabolism, and high oxygen consumption, hypoxia will increase its mortality and cause serious economic losses. We constructed two experimental groups of fish with different degrees of tolerance to hypoxia, used BSR-Seq analysis based on genome and genetic linkage groups to locate SNPs and genes that were related to the differences in hypoxia tolerance. The results showed that hypoxia tolerance SNPs of golden pompano may be jointly determined by multiple linkage groups, especially linkage groups 18 and 22. There were 768 and 348 candidate genes located in the candidate regions of the brain and liver, respectively. These genes were mainly involved in anaerobic energy metabolism, stress response, immune response, waste discharge, and cell death. The prostaglandin-endoperoxide synthase 2 (PTGS2) on LG8, which is involved in the metabolism of arachidonic acid, has a G/A nonsynonymous mutation at position 20641628, and the encoded amino acid was changed from hydrophobic aspartic acid to asparaginate. The specific pathway of the RIG-I-like receptor signaling pathway in the liver may mediate the metabolic system and the immune system, linking glucose metabolism with immune regulation. The death of the hypoxia-intolerant group may be due to the accumulation of lactic acid caused by the activation of anaerobic glycolysis during the early stage of hypoxia stress, and the activation of type I interferon was inhibited, which resulted in decreased immunity. Among the genes involved in the RIG-I-like receptor signaling pathway, the CYLD Lysine 63 Deubiquitinase (CYLD) located on LG16 had a G/T nonsynonymous mutation at position 13629651, and the encoded amino acid was changed from alanine acid to valine. The interferon induced with helicase C domain 1 (Ifih1) located on LG18 has a G/C nonsynonymous mutation at position 16153700, and the encoded hydrophilic glycine was changed to hydrophobic alanine. Our findings suggest these SNPs may assist in the molecular breeding of hypoxia-tolerant golden pompano, and speculate that the balance of glucose and lipid metabolism plays a key role in Trachinotus blochii under acute hypoxia.

Physiological responses and molecular strategies in heart of silver carp (Hypophthalmichthys molitrix) under hypoxia and reoxygenation

A sufficient oxygen concentration is essential for fish growth, reproduction, and metabolism. Silver carp (Hypophthalmichthys molitrix) is sometimes challenged by hypoxia during intensive aquaculture or because of environmental changes. However, the response to hypoxic stress in the heart of silver carp remains relatively unknown. In the present study, we reported the effects of hypoxia on histological structures, enzyme activities, and gene expression in the heart of silver carp. Hematoxylin and eosin (H&E) staining of heart sections showed that the myocardial fibers gradually became disordered, swollen, and even ruptured during hypoxic treatment. These phenotypes were also supported by increased activities of injury-related enzymes. Moreover, the transcriptome was analyzed to determine the molecular strategies of hypoxia adaptation in the heart. PI3K-Akt signaling pathway, FoxO signaling pathway, and JAK-STAT signaling pathway were the most prominent pathways activated by hypoxia. Twenty significantly differentially expressed genes were selected to create a network diagram related to cell proliferation, carbohydrate metabolism, oxidative stress, and angiogenesis. Additionally, reoxygenation could ameliorate cardiac injury and eliminate the effects of hypoxia on gene expression. This was the first comparative transcriptomic study to explore the molecular mechanism of the response to hypoxia and reoxygenation in the heart of silver carp. Our results provide a theoretical basis for cultivating hypoxia-tolerant carp varieties in the future.

Effects of hypoxia and hyperoxia on growth parameters and transcription levels of growth, immune system and stress related genes in rainbow trout

Hypoxia and hyperoxia are disparate stressors which can have destructive influences on fish growth and physiology. It is yet to be determined if hypoxia and hyperoxia have a cumulative effect in aquatic ecosystems that affect biological parameters in fish, and to understand if this is associated with gene expression. Here we address whether growth performance and expressions of growth, immune system and stress related genes were affected by hypoxia and hyperoxia in fish. Rainbow trout was chosen as the study organism due to its excellent service as biomonitor. After an acclimatization period, fish were exposed to hypoxia (4.0 ± 0.5 ppm O₂), normoxia (7.5 ± 0.5 ppm O₂) and hyperoxia (12 ± 1.2 ppm O₂) for 28 days. At 6 h, 12 h, 24 h, 48 h, 72 h and 28 days, samples were collected. Hypoxia and hyperoxia negatively affected weight gain (WG), specific growth rate (SGR), survival rate (SR) and feed conversion ratio (FCR). The best WG, SGR, SR and FCR values occurred in fish exposed to normoxia, whereas hypoxia was most suppressive on growth and hyperoxia showed intermediate suppression of these parameters. Gene expression analyses were performed in liver and results revealed that long term exposure caused reduced growth hormone-I (GH-I) and insulin like growth factor I-II (IGF I-II) levels in both hypoxia and hyperoxia-treated fish. Heat shock protein (HSP70) levels increased in both hypoxia and hyperoxia treatment, and both exposures caused elevation of leptin (LEP) expression in long-term exposure. Overall data indicate that both hypoxia and hyperoxia cause stress in rainbow trout and negatively affects growth parameters.

The influence of hypoxia on the cardiac transcriptomes of two estuarine species - C. variegatus and F. grandis

Increased nutrient loading has led to eutrophication of coastal shelf waters which has resulted in increased prevalence of persistent hypoxic zones - areas in which the dissolved oxygen content of the water drops below 2 mg/L. The northern Gulf of Mexico, fed primarily by the Mississippi River watershed, undergoes annual establishment of one of the largest hypoxic zones in the world. Exposure to hypoxia can induce physiological impacts in fish cardiac systems that include bradycardia, changes in stroke volume, and altered cardiovascular vessel development. While these impacts have been addressed at the functional level, there is little information regarding the molecular basis for these changes. This study used transcriptomic analysis techniques to interrogate the effects of hypoxia exposure on the developing cardiovascular system in newly hatched larvae of two estuarine species that occupy the same ecological niche - the sheepshead minnow (Cyprinodon variegatus) and the Gulf killifish (Fundulus grandis). Results suggest that while differential gene expression is largely distinct between the two species, downstream impacts on pathways and functional responses such as reduced cardiac hypertrophy, modulation of blood pressure, and increased incidence of apoptosis appear to be conserved. Further, differences in the magnitude of these conserved responses may suggest that the length of embryonic development could impart a level of resiliency to hypoxic perturbation in early life stage fish.

Transcriptome analysis reveals molecular strategies in gills and heart of large yellow croaker (Larimichthys crocea) under hypoxia stress

The gills and heart are two major targets of hypoxia in fish. However, the molecular responses in fish gills and heart to hypoxia challenge remain unclear. Here, RNA-Seq technology was used to study the gene expression profiles in gills and heart of large yellow croaker (Larimichthys crocea) at 6, 24, and 48 h after hypoxia stress. A total of 1,546 and 2,746 differentially expressed genes (DEGs) were identified in gills and heart, respectively. Expression changes of nine genes in each tissue were further validated by the qPCR. Based on KEGG and Gene ontology enrichments, we found that various innate immunity-related genes, such as complement components (C1qs, C2, C3, C6, and C7), chemokines (CCL3, CCL17, CCL19, CCL25, and CXCL8_L3), chemokine receptors (CCR9, CXCR1, and CXCR3), and nitric oxide synthase (NOS), were significantly down-regulated in gills and/or heart, suggesting that innate immune processes mediated by these genes may be inhibited by hypoxia. The genes involved in both glycolysis pathway (LDHA) and tricarboxylic acid cycle (IDH2 and OGDH) were up-regulated in gills and heart of hypoxic large yellow croakers, possibly because gill and heart tissues need enough energy to accelerate gas exchange and blood circulation. Hypoxia also affected the ion transport in gills of large yellow croaker, through down-regulating the expression levels of numerous classical ion transporters, including HVCN1, SLC20A2, SLC4A4, RHBG, RHCG, and SCN4A, suggesting an energy conservation strategy to hypoxia stress. All these results indicate that the immune processes, glycolytic pathways, and ion transport were significantly altered in gills and/or heart of large yellow croaker under hypoxia, possibly contributing to maintain cellular energy balance during hypoxia. Our data, therefore, afford new information to understand the tissue-specific molecular responses of bony fish to hypoxia stress.

Acute hypoxia changes the mode of glucose and lipid utilization in the liver of the largemouth bass (Micropterus salmoides)

Dissolved oxygen (DO) undountedly affects fish distribution, metabolism, and evern survival. Intensive aquaculture and environmental changes will inevitably lead to hypoxic stress for largemouth bass (Micropterus salmoides). The different metabolic responses and mechanism still remains relatively unknown during acute hypoxia exposure. In this study, largemouth bass were subjected to hypoxic stress (3.0 ± 0.2 mg/L and 1.2 ± 0.2 mg/L) for 24 h and 12 h reoxygenation to systemically evaluate indicators of glucose and lipid metabolism. A regulatory network was constructed using RNA-seq to further elucidate the transcriptional regulation of glucose and lipid metabolism. During hypoxia for 4 h, the liver glycogen, glucose and pyruvic acid contents significantly decreased, whereas plasma glucose content and liver lactic acid content increased significantly. The accumulation of liver triglycerides and non-esterified fatty acids was enhanced during hypoxia for 8 h. The activity of key enzymes revealed the different metabolic responses to hypoxia exposure for 4 h, including the enhancement of glycolysis, and inhibition of gluconeogenesis. Furthermore, hypoxia exposure for 8 h increased lipid mobilization, and inhibited the β-oxidation. In addition, an integrated regulatory network of 9 major pathways involved in the response to hypoxia exposure was constructed, including HIF signaling pathway, VEGF signaling pathway, AMPK signaling pathway, insulin signaling pathway and PPAR signaling pathway; glycolysis/gluconeogenesis, pyruvate metabolism, fatty acid degradation and fatty acid biosynthesis. Additionally, reoxygenation inhibited glycolysis, and promoted gluconeogenesis and lipid oxidation, but energy deficits persisted. In short, although the mobilization and activation of fatty acid in liver were enhanced in the early stage of hypoxia, glycolysis was the main energy source under acute hypoxia. The extent and duration of hypoxia determine the degree of change in energy metabolism.

Tissue-Specific Orchestration of Gilthead Sea Bream Resilience to Hypoxia and High Stocking Density

Two different O₂ levels (normoxia: 75–85% O₂ saturation; moderate hypoxia: 42–43% O₂ saturation) and stocking densities (LD: 9.5, and HD: 19 kg/m³) were assessed on gilthead sea bream (Sparus aurata) in a 3-week feeding trial. Reduced O₂ availability had a negative impact on feed intake and growth rates, which was exacerbated by HD despite of the improvement in feed efficiency. Blood physiological hallmarks disclosed the enhancement in O₂-carrying capacity in fish maintained under moderate hypoxia. This feature was related to a hypo-metabolic state to cope with a chronic and widespread environmental O₂ reduction, which was accompanied by a differential regulation of circulating cortisol and growth hormone levels. Customized PCR-arrays were used for the simultaneous gene expression profiling of 34–44 selected stress and metabolic markers in liver, white skeletal muscle, heart, and blood cells. The number of differentially expressed genes ranged between 22 and 19 in liver, heart, and white skeletal muscle to 5 in total blood cells. Partial Least-Squares Discriminant Analysis (PLS-DA) explained [R2Y(cum)] and predicted [Q2Y(cum)] up to 95 and 65% of total variance, respectively. The first component (R2Y = 0.2889) gathered fish on the basis of O₂ availability, and liver and cardiac genes on the category of energy sensing and oxidative metabolism (cs, hif-1α, pgc1α, pgc1β, sirts 1-2-4-5-6-7), antioxidant defense and tissue repair (prdx5, sod2, mortalin, gpx4, gr, grp-170, and prdx3) and oxidative phosphorylation (nd2, nd5, and coxi) highly contributed to this separation. The second component (R2Y = 0.2927) differentiated normoxic fish at different stocking densities, and the white muscle clearly promoted this separation by a high over-representation of genes related to GH/IGF system (ghr-i, igfbp6b, igfbp5b, insr, igfbp3, and igf-i). The third component (R2Y = 0.2542) discriminated the effect of stocking density in fish exposed to moderate hypoxia by means of hepatic fatty acid desaturases (fads2, scd1a, and scd1b) and muscle markers of fatty acid oxidation (cpt1a). All these findings disclose the different contribution of analyzed tissues (liver ≥ heart > muscle > blood) and specific genes to the hypoxic- and crowding stress-mediated responses. This study will contribute to better explain and understand the different stress resilience of farmed fish across individuals and species.

Genome-Wide Microarray Analysis Suggests Transcriptomic Response May Not Play a Major Role in High- to Low-Altitude Acclimation in Harvest Mouse (Micromys minutus)

Simple Summary

Micromys minutus is a small rodent species that has a wide range of vertical distribution in Taiwan. By comparing the gene expression profile of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site, the Tnfrsf12a gene was demonstrated to have a differential expression pattern. Although this finding may be correlated with the altitude acclimation, the observation of only one gene transcript with significant alteration leads us to suggest that genetic response may not play a major role in altitude acclimation in M. minutus. Future comparative functional genomics studies involving other organ systems (in addition to skeletal muscles) and alarger sample size are warranted for better insight into the altitude acclimation of this small rodent species.

Abstract

The harvest mouse (Micromys minutus) is a small rodent species with a wide range of vertical distribution in Taiwan, extending from the sea level to 3100 m altitude. This species has recently suffered from habitat loss in high-altitude areas due to orchard cultivation, which may have resulted in mouse migration from high to low altitude. To investigate whether there is any physiological mechanism involved in altitude acclimation, rat cDNA microarray was used to compare transcriptomic patterns of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site. Of the 23,188 genes being analyzed, 47 (33 up-regulated and 14 down-regulated) were found to have differential expression (fold change > 4 or < −4, ANOVA p < 0.05). However, after multiple testing correction with a false discovery rate (FDR), only the result for Tnfrsf12a was found to be statistically significant (fold change = 13, FDR p < 0.05). The result was confirmed by quantitative polymerase chain reaction (q-PCR). The expression of Tnfrsf12a possibly relates to the skeletal muscle biology and thus can be correlated with altitude acclimation. However, finding only one gene transcript with significant alteration suggests that transcriptomic response may not play a major role in high- to low-altitude acclimation in harvest mouse.

Can variation among hypoxic environments explain why different fish species use different hypoxic survival strategies?

In aquatic environments, hypoxia is a multi-dimensional stressor that can vary in O₂ level (partial pressure of O₂ in water, Pw_O2 ), rate of induction and duration. Natural hypoxic environments can therefore be very different from one another. For the many fish species that have evolved to cope with these different hypoxic environments, survival requires adjusting energy supply and demand pathways to maintain energy balance. The literature describes innumerable ways that fishes combine aerobic metabolism, anaerobic metabolism and metabolic rate depression (MRD) to accomplish this, but it is unknown whether the evolutionary paths leading to these different strategies are determined primarily by species' phylogenetic histories, genetic constraint or their native hypoxic environments. We explored this idea by devising a four-quadrant matrix that bins different aquatic hypoxic environments according to their duration and Pw_O2  characteristics. We then systematically mined the literature for well-studied species native to environments within each quadrant, and, for each of 10 case studies, described the species' total hypoxic response (THR), defined as its hypoxia-induced combination of sustained aerobic metabolism, enhanced anaerobic metabolism and MRD, encompassing also the mechanisms underlying these metabolic modes. Our analysis revealed that fishes use a wide range of THRs, but that distantly related species from environments within the same matrix quadrant have converged on similar THRs. For example, environments of moderately hypoxic Pw_O2  favoured predominantly aerobic THRs, whereas environments of severely hypoxic Pw_O2  favoured MRD. Capacity for aerial emergence as well as predation pressure (aquatic and aerial) also contributed to these responses, in addition to other biotic and abiotic factors. Generally, it appears that the particular type of hypoxia experienced by a fish plays a major role in shaping its particular THR.

Genome-wide identification and differentially expression analysis of lncRNAs in tilapia

Background

Long noncoding RNAs (LncRNAs) play important roles in fundamental biological processes. However, knowledge about the genome-wide distribution and stress-related expression of lncRNAs in tilapia is still limited.

Results

Genome-wide identification of lncRNAs in the tilapia genome was carried out in this study using bioinformatics tools. 103 RNAseq datasets that generated in our laboratory or collected from NCBI database were analyzed. In total, 72,276 high-confidence lncRNAs were identified. The averaged positive correlation coefficient (r_mean = 0.286) between overlapped lncRNA and mRNA pairs showed significant differences with the values for all lncRNA-mRNA pairs (r_mean = 0.176, z statistics = − 2.45, p value = 0.00071) and mRNA-mRNA pairs (r_mean = 0.186, z statistics = − 2.23, p value = 0.0129). Weighted correlation network analysis of the lncRNA and mRNA datasets from 12 tissues identified 21 modules and many interesting mRNA genes that clustered with lncRNAs. Overrepresentation test indicated that these mRNAs enriched in many biological processes, such as meiosis (p = 0.00164), DNA replication (p = 0.00246), metabolic process (p = 0.000838) and in molecular function, e.g., helicase activity (p = 0.000102) and catalytic activity (p = 0.0000612). Differential expression (DE) analysis identified 99 stress-related lncRNA genes and 1955 tissue-specific DE lncRNA genes. MiRNA-lncRNA interaction analysis detected 72,267 lncRNAs containing motifs with sequence complementary to 458 miRNAs.

Conclusions

This study provides an invaluable resource for further studies on molecular bases of lncRNAs in tilapia genomes. Further function analysis of the lncRNAs will help to elucidate their roles in regulating stress-related adaptation in tilapia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5115-x) contains supplementary material, which is available to authorized users.

Acute hypoxia stress induced abundant differential expression genes and alternative splicing events in heart of tilapia

Hypoxia is one of the critical environmental stressors for fish in aquatic environments. Although accumulating evidences indicate that gene expression is regulated by hypoxia stress in fish, how genes undergoing differential gene expression and/or alternative splicing (AS) in response to hypoxia stress in heart are not well understood. Using RNA-seq, we surveyed and detected 289 differential expressed genes (DEG) and 103 genes that undergo differential usage of exons and splice junctions events (DUES) in heart of a hypoxia tolerant fish, Nile tilapia, Oreochromis niloticus following 12h hypoxic treatment. The spatio-temporal expression analysis validated the significant association of differential exon usages in two randomly selected DUES genes (fam162a and ndrg2) in 5 tissues (heart, liver, brain, gill and spleen) sampled at three time points (6h, 12h, and 24h) under acute hypoxia treatment. Functional analysis significantly associated the differential expressed genes with the categories related to energy conservation, protein synthesis and immune response. Different enrichment categories were found between the DEG and DUES dataset. The Isomerase activity, Oxidoreductase activity, Glycolysis and Oxidative stress process were significantly enriched for the DEG gene dataset, but the Structural constituent of ribosome and Structural molecule activity, Ribosomal protein and RNA binding protein were significantly enriched only for the DUES genes. Our comparative transcriptomic analysis reveals abundant stress responsive genes and their differential regulation function in the heart tissues of Nile tilapia under acute hypoxia stress. Our findings will facilitate future investigation on transcriptome complexity and AS regulation during hypoxia stress in fish.

Gene expression profiling of whole blood cells supports a more efficient mitochondrial respiration in hypoxia-challenged gilthead sea bream (Sparus aurata)

BACKGROUND

Acclimation to abiotic challenges, including decreases in O2 availability, requires physiological and anatomical phenotyping to accommodate the organism to the environmental conditions. The retention of a nucleus and functional mitochondria in mature fish red blood cells makes blood a promising tissue to analyse the transcriptome and metabolic responses of hypoxia-challenged fish in an integrative and non-invasive manner.

METHODS

Juvenile gilthead sea bream (Sparus aurata) were reared at 20-21 °C under normoxic conditions (> 85% O2 saturation) followed by exposure to a gradual decrease in water O2 concentration to 3.0 ppm (41-42% O2 saturation) for 24 h or 1.3 ppm (18-19% O2 saturation) for up to 4 h. Blood samples were collected at three different sampling points for haematological, biochemical and transcriptomic analysis.

RESULTS

Blood physiological hallmarks remained almost unaltered at 3.0 ppm, but the haematocrit and circulating levels of haemoglobin, glucose and lactate were consistently increased when fish were maintained below the limiting oxygen saturation at 1.3 ppm. These findings were concurrent with an increase in total plasma antioxidant activity and plasma cortisol levels, whereas the opposite trend was observed for growth-promoting factors, such as insulin-like growth factor I. Additionally, gene expression profiling of whole blood cells revealed changes in upstream master regulators of mitochondria (pgcβ and nrf1), antioxidant enzymes (gpx1, gst3, and sod2), outer and inner membrane translocases (tom70, tom22, tim44, tim10, and tim9), components of the mitochondrial dynamics system (mfn2, miffb, miro1a, and miro2), apoptotic factors (aifm1), uncoupling proteins (ucp2) and oxidative enzymes of fatty acid β-oxidation (acca2, ech, and hadh), the tricarboxylic acid cycle (cs) and the oxidative phosphorylation pathway. The overall response is an extensive reduction in gene expression of almost all respiratory chain enzyme subunits of the five complexes, although mitochondrial-encoded catalytic subunits and nuclear-encoded regulatory subunits of Complex IV were primarily increased in hypoxic fish.

CONCLUSIONS

Our results demonstrate the re-adjustment of mitochondrial machinery at transcriptional level to cope with a decreased basal metabolic rate, consistent with a low risk of oxidative stress, diminished aerobic ATP production and higher O2-carrying capacity. Taken together, these results suggest that whole blood cells can be used as a highly informative target tissue of metabolic condition.

Dietary supplementation of heat-treated Gracilaria and Ulva seaweeds enhanced acute hypoxia tolerance in gilthead sea bream (Sparus aurata)

Intensive aquaculture practices involve rearing fish at high densities. In these conditions, fish may be exposed to suboptimal dissolved O2 levels with an increased formation of reactive O2 species (ROS) in tissues. Seaweeds (SW) contain biologically active substances with efficient antioxidant capacities. This study evaluated the effects of dietary supplementation of heat-treated SW (5% Gracilaria vermiculophylla or 5% Ulva lactuca) on stress bioindicators in sea bream subjected to a hypoxic challenge. 168 fish (104.5 g average weight) were distributed in 24 tanks, in which eight tanks were fed one of three experimental diets for 34 days: (i) a control diet without SW supplementation, (ii) a control diet supplemented with Ulva, or (iii) a control diet with Gracilaria Thereafter, fish from 12 tanks (n=4 tanks/dietary treatment) were subjected to 24 h hypoxia (1.3 mg O2 l-1) and subsequent recovery normoxia (8.6 mg O2 l-1). Hypoxic fish showed an increase in hematocrit values regardless of dietary treatment. Dietary modulation of the O2-carrying capacity was conspicuous during recovery, as fish fed SW supplemented diets displayed significantly higher haemoglobin concentration than fish fed the control diet. After the challenge, survival rates in both groups of fish fed SW were higher, which was consistent with a decrease in hepatic lipid peroxidation in these groups. Furthermore, the hepatic antioxidant enzyme activities were modulated differently by changes in environmental O2 condition, particularly in sea bream fed the Gracilaria diet. After being subjected to hypoxia, the gene expression of antioxidant enzymes and molecular chaperones in liver and heart were down regulated in sea bream fed SW diets. This study suggests that the antioxidant properties of heat-treated SW may have a protective role against oxidative stress. The nature of these compounds and possible mechanisms implied are currently being investigated.

Sequence and functional characterization of hypoxia-inducible factors, HIF1α, HIF2αa, and HIF3α, from the estuarine fish, Fundulus heteroclitus

The hypoxia-inducible factor (HIF) family of transcription factors plays central roles in the development, physiology, pathology, and environmental adaptation of animals. Because many aquatic habitats are characterized by episodes of low dissolved oxygen, fish represent ideal models to study the roles of HIF in the response to aquatic hypoxia. The estuarine fish Fundulus heteroclitus is found in habitats prone to hypoxia. It responds to low oxygen via behavioral, physiological, and molecular changes, and one member of the HIF family, HIF2α, has been previously described. Herein, cDNA sequencing, phylogenetic analyses, and genomic approaches were used to determine other members of the HIFα family from F. heteroclitus and their relationships to HIFα subunits from other vertebrates. In vitro and cellular approaches demonstrated that full-length forms of HIF1α, HIF2α, and HIF3α independently formed complexes with the β-subunit, aryl hydrocarbon receptor nuclear translocator, to bind to hypoxia response elements and activate reporter gene expression. Quantitative PCR showed that HIFα mRNA abundance varied among organs of normoxic fish in an isoform-specific fashion. Analysis of the F. heteroclitus genome revealed a locus encoding a second HIF2α-HIF2αb-a predicted protein lacking oxygen sensing and transactivation domains. Finally, sequence analyses demonstrated polymorphism in the coding sequence of each F. heteroclitus HIFα subunit, suggesting that genetic variation in these transcription factors may play a role in the variation in hypoxia responses among individuals or populations.

Multiple across-strain and within-strain QTLs suggest highly complex genetic architecture for hypoxia tolerance in channel catfish

The ability to survive hypoxic conditions is important for various organisms, especially for aquatic animals. Teleost fish, representing more than 50 % of vertebrate species, are extremely efficient in utilizing low levels of dissolved oxygen in water. However, huge variations exist among various taxa of fish in their ability to tolerate hypoxia. In aquaculture, hypoxia tolerance is among the most important traits because hypoxia can cause major economic losses. Genetic enhancement for hypoxia tolerance in catfish is of great interest, but little was done with analysis of the genetic architecture of hypoxia tolerance. The objective of this study was to conduct a genome-wide association study to identify QTLs for hypoxia tolerance using the catfish 250K SNP array with channel catfish families from six strains. Multiple significant and suggestive QTLs were identified across and within strains. One significant QTL and four suggestive QTLs were identified across strains. Six significant QTLs and many suggestive QTLs were identified within strains. There were rare overlaps among the QTLs identified within the six strains, suggesting a complex genetic architecture of hypoxia tolerance. Overall, within-strain QTLs explained larger proportion of phenotypic variation than across-strain QTLs. Many of genes within these identified QTLs have known functions for regulation of oxygen metabolism and involvement in hypoxia responses. Pathway analysis indicated that most of these genes were involved in MAPK or PI3K/AKT/mTOR signaling pathways that were known to be important for hypoxia-mediated angiogenesis, cell proliferation, apoptosis and survival.

The Effect of Hypoxia and Hyperoxia on Growth and Expression of Hypoxia-Related Genes and Proteins in Spotted Gar Lepisosteus oculatus Larvae and Juveniles

We studied the molecular responses to different water oxygen levels in gills and swim bladder of spotted gar (Lepisosteus oculatus), a bimodal breather. Fish at swim-up stage were exposed for 71 days to normoxic, hypoxic, and hyperoxic water conditions. Then, all aquaria were switched to normoxic conditions for recovery until the end of the experiment (120 days). Fish were sampled at the beginning of the experiment, and then at 71 days of exposure and at 8 days of recovery. We first cloned three hypoxia-related genes, hypoxia-inducible factor 2α (HIF-2α), Na(+) /H(+) exchanger 1 (NHE-1), and NHE-3, and uploaded their cDNA sequences in the GeneBank database. We then used One Step Taqman® real-time PCR to quantify the mRNA copies of target genes in gills and swim bladder of fish exposed to different water O2 levels. We also determined the protein expression of HIF-2α and neuronal nitric oxide synthase (nNOS) in the swim bladder by using confocal immunofluorescence. Hypoxic stress for 71 days significantly increased the mRNA copies of HIF-2α and NHE-1 in gills and swim bladder, whereas normoxic recovery for 8 days decreased the HIF-2α mRNA copies to control values in both tissues. We did not found significant changes in mRNA copies of the NHE-3 gene in either gills or swim bladder in response to hypoxia and hyperoxia. Unlike in normoxic swim bladder, double immunohistochemical staining in hypoxic and hyperoxic swim bladder using nNOS/HIF-2α showed extensive bundles of HIF-2α-positive nerve fibers in the trabecular musculature associated with a few varicose nNOS immunoreactive nerve terminals.

Integrated analysis of mRNA-seq and miRNA-seq in the liver of Pelteobagrus vachelli in response to hypoxia

Pelteobagrus vachelli is a well-known commercial species in Asia. However, a sudden lack of oxygen will result in mortality and eventually to pond turnover. Studying the molecular mechanisms of hypoxia adaptation in fishes will not only help us to understand fish speciation and the evolution of the hypoxia-signaling pathway, but will also guide us in the breeding of hypoxia-tolerant fish strains. Despite this, the genetic regulatory network for miRNA-mRNA and the signaling pathways involved in hypoxia responses in fish have remained unexamined. In the present study, we used next-generation sequencing technology to characterise mRNA-seq and miRNA-seq of control- and hypoxia-treated P. vachelli livers to elucidate the molecular mechanisms of hypoxia adaptation. We were able to find miRNA-mRNA pairs using bioinformatics analysis and miRNA prediction algorithms. Furthermore, we compared several key pathways which were identified as involved in the hypoxia response of P. vachelli. Our study is the first report on integrated analysis of mRNA-seq and miRNA-seq in fishes and offers a deeper insight into the molecular mechanisms of hypoxia adaptation. qRT-PCR analysis further confirmed the results of mRNA-Seq and miRNA-Seq analysis. We provide a good case study for analyzing mRNA/miRNA expression and profiling a non-model fish species using next-generation sequencing technology.

Proteomic analysis of the European flounder Platichthys flesus response to experimental PAH-PCB contamination

Platichthys flesus is often used as a sentinel species to monitor the estuarine water quality. In this study, we carried out an experimental contamination of fish using a PAHs/PCBs mixture, which was designed to mimic the concentrations found in the Seine estuary (C1) and 10 times these concentrations (C2). We used a proteomic approach to understand the molecular mechanisms implied in the response of P. flesus to these xenobiotics. We showed that 54 proteins were differentially accumulated in one or several conditions, which 34 displayed accumulation factors higher than two. 18 of these proteins were identified by MALDI TOF-TOF mass spectrometry. The results indicated the deregulation of oxidative stress- and glutathione metabolism-(GST, GPx) proteins as well as of several proteins belonging to the betaine demethylation pathway and the methionine cycle (BHMT, SHMT, SAHH), suggesting a role for these different pathways in the P. flesus response to chemical contamination.

Unraveling the molecular signatures of oxidative phosphorylation to cope with the nutritionally changing metabolic capabilities of liver and muscle tissues in farmed fish

Mitochondrial oxidative phosphorylation provides over 90% of the energy produced by aerobic organisms, therefore the regulation of mitochondrial activity is a major issue for coping with the changing environment and energy needs. In fish, there is a large body of evidence of adaptive changes in enzymatic activities of the OXPHOS pathway, but less is known at the transcriptional level and the first aim of the present study was to define the molecular identity of the actively transcribed subunits of the mitochondrial respiratory chain of a livestock animal, using gilthead sea bream as a model of farmed fish with a high added value for European aquaculture. Extensive BLAST searches in our transcriptomic database (www.nutrigroup-iats.org/seabreamdb) yielded 97 new sequences with a high coverage of catalytic, regulatory and assembly factors of Complex I to V. This was the basis for the development of a PCR array for the simultaneous profiling of 88 selected genes. This new genomic resource allowed the differential gene expression of liver and muscle tissues in a model of 10 fasting days. A consistent down-regulated response involving 72 genes was made by the liver, whereas an up-regulated response with 29 and 10 differentially expressed genes was found in white skeletal muscle and heart, respectively. This differential regulation was mostly mediated by nuclear-encoded genes (skeletal muscle) or both mitochondrial- and nuclear-encoded genes (liver, heart), which is indicative of a complex and differential regulation of mitochondrial and nuclear genomes, according to the changes in the lipogenic activity of liver and the oxidative capacity of glycolytic and highly oxidative muscle tissues. These insights contribute to the identification of the most responsive elements of OXPHOS in each tissue, which is of relevance for the appropriate gene targeting of nutritional and/or environmental metabolic disturbances in livestock animals.

Identification of hypoxia-regulated genes in the liver of common sole (Solea solea) fed different dietary lipid contents

Coastal systems could be affected by hypoxic events brought about by global change. These areas are essential nursery habitats for several fish species including the common sole (Solea solea L.). Tolerance of fish to hypoxia depends on species and also on their physiological condition and nutritional status. Indeed, high dietary lipid content has been recently shown to negatively impact the resistance of sole to a severe hypoxic challenge. In order to study the molecular mechanisms involved in the early response to hypoxic stress, the present work examined the hepatic transcriptome in common sole fed diets with low and high lipid content, exposed to severe hypoxia. The activity of AMP-activated protein kinase (AMPK) was also investigated through the quantification of threonine-172 phosphorylation in the alpha subunit. The results show that hypoxia consistently regulates several actors involved in energy metabolism pathways and particularly AMPKα, as well as some involved in cell growth and maintenance or unfolded protein response. Our findings reveal that (1) the expression of genes involved in biological processes with high energy cost or implicated in aerobic ATP synthesis was down-regulated by hypoxia, contrary to genes involved in neoglucogenesis or in angiogenesis, (2) the consumption of high lipid induced regulation of metabolic pathways going against this energy saving, and (3) this control was fine-tuned by the regulation of several transcriptomic factors. These results provide insight into the biological processes involved in the hepatic response to hypoxic stress and underline the negative impact of high lipid consumption on the tolerance of common sole to hypoxia.

The relationship between phenotypic and environmental variation: do physiological responses reduce interindividual differences?

What is the effect of a variable environment on phenotypic variation? Does the physiological response to a new environment increase or decrease the differences among individuals? We provide a speculative hypothesis suggesting that the induction of a physiological response to environmental change minimizes phenotypic differences among individuals in outbred genetically variable populations. Although this suggestion runs counter to the general idea that environmental variation induces phenotypic variation, we provide evidence that this is not always the case. One explanation for this counterintuitive hypothesis is that in a variable environment, the physiological mechanism that maintains homeostasis changes the concentrations of active transcription factors (TFs). This change in TFs reduces the effectiveness of nucleotide polymorphisms in TF binding sites and thus reduces the variation among individuals in mRNA expression and in the phenotypes affected by these mRNA transcripts. Thus, there are fewer differences among individuals in a variable environment compared with the variation observed in a constant environment. Our conjecture is that the physiological mechanisms that maintain homeostasis in response to environmental variation canalize phenotypic variation. If our hypothesis is correct, then the physiological canalization of gene expression in a variable environment hides genetic variation and thereby reduces the evolutionary costs of polymorphism. This hypothesis provides a new perspective on the mechanisms by which high levels of genetic variation can persist in real-world populations.