Transcriptional regulation of hypoxia inducible factors alpha (HIF-α) and their inhibiting factor (FIH-1) of channel catfish (Ictalurus punctatus) under hypoxia

Linking HIF oxygen-sensing system diversity to hypoxia fitness in Eleutheronema: Molecular characterization and transcriptional response to hypoxia exposure

Hypoxia is a mounting environmental problem affecting coastal waters globally, posing severe consequences for biodiversity and marine life. Metazoans respond to hypoxia stress via the hypoxia-inducible factor (HIF) pathway, but few studies have addressed the gene diversity of the functionally important HIF-pathway. Understanding whether functional diversity exists in the HIF-pathway is a key first step in identifying genes that may impact hypoxia fitness. Here, we leveraged whole-genome resequencing data and bioinformatics tools to identify the key members of the HIF-pathway (HIFα/β, EGLN, and VHL) and genetic diversity in the threatened Eleutheronema. Phylogenetic analysis revealed that teleost-specific duplicates of epas1 (epas1a/b) were followed by the loss of one of each hif1α and hif1αl in Eleutheronema species. Strong collinearity and similarity of gene characteristics suggested the functional conservation of the HIF-pathway during Eleutheronema evolution. Purifying selection was the major theme in HIF-pathway evolution, leading to a reduction in genetic diversity. Substantially low nucleotide diversity of the HIF-pathway was observed among populations, which might indicate the loss of hypoxia fitness in Eleutheronema. Additionally, the normoxic presence of the HIF-pathway differed among tissues and was species-dependent, indicating their diverse roles during development. Significant regulation of HIF-pathway expression levels was observed across tissues under hypoxic conditions, suggesting critical roles in the hypoxia stress response. Moreover, variant molecular characters suggested different roles in response to hypoxia of the HIF-pathway, which were reflected in the different expression patterns across tissues. Our present study provides novel insights into the interplay between gene diversity within the HIF-pathway and hypoxia fitness in threatened Eleutheronema. We highlighted the importance of HIF-pathway-mediated transcriptional regulation in response to hypoxia stress, which provided valuable information for the genetic mechanisms underlying hypoxia adaptation in fish. The bioinformatic methods developed here have broad applications for other species.

Hypoxia-inducible factor 1 alpha protein increases without changes in mRNA during acute hypoxic exposure of the Gulf killifish, Fundulus grandis

The hypoxia inducible factor 1 (HIF1) is a central regulator of the molecular responses of animals to low oxygen. While the hypoxia-responsiveness of HIF1 is generally attributed to the stabilization of the alpha protein subunit (HIF1α) at low oxygen, several studies on fish report increased tissue levels of HIF1A mRNA during hypoxia, suggesting transcriptional regulation. In the current study, HIF1α protein and HIF1A mRNA were determined in parallel in tissues of Gulf killifish, Fundulus grandis, exposed to short-term hypoxia (24 h at 1 mg O2 l-1). HIF1α protein was higher in brain, ovary, and skeletal muscle from fish exposed to hypoxia compared with normoxic controls by 6 h, and it remained elevated in brain and ovary at 24 h. In contrast, HIF1A mRNA levels were unaffected by hypoxia in any tissue. Moreover, HIF1α protein and HIF1A mRNA levels in the same tissues were not correlated with one another during either normoxia or hypoxia. Hence, an increase in HIF1α protein does not depend upon an increase in HIF1A mRNA during acute exposure to low oxygen in this species. The results support the widely accepted mechanism of post-translational protein stabilization, rather than new transcription, during the initial response of fish to hypoxia.

Role of HIF in fish inflammation

The hypoxia-inducing factor (HIF) is a central transcription factor in cellular oxygen sensing and regulation. It is common that the inflammation always appears in many diseases, like infectious diseases in fishes, and the inflammation is often accompanied by hypoxia, as a hallmark of inflammation. Besides coordinating cellular responses to low oxygen, HIF-mediated hypoxia signaling pathway is also crucial for immune responses such as the regulations of innate immune cell phenotype and function, as well as metabolic reprogramming under the inflammation. However, the understanding of the molecular mechanisms by which HIFs regulate the inflammatory response in fish is still very limited. Here, we review the characteristics of HIF as well as its roles in innate immune cells and the infections caused by bacteria and viruses. The regulatory effects of HIF on the metabolic reprogramming of innate immune cells are also discussed and the future research directions are outlooked. This paper will serve as a reference for elucidating the molecular mechanism of HIF regulating inflammation and identifying treatment strategies to target HIF for fish disease.

Transcriptomic Response of the Liver Tissue in Trachinotus ovatus to Acute Heat Stress

Trachinotus ovatus is a major economically important cultured marine fish in the South China Sea. However, extreme weather and increased culture density result in uncontrollable problems, such as increases in water temperature and a decline in dissolved oxygen (DO), hindering the high-quality development of aquaculture. In this study, liver transcriptional profiles of T. ovatus were investigated under acute high-temperature stress (31 °C and 34 °C) and normal water temperature (27 °C) using RNA sequencing (RNA-Seq) technology. Differential expression analysis and STEM analysis showed that 1347 differentially expressed genes (DEGs) and four significant profiles (profiles 0, 3, 4, and 7) were screened, respectively. Of these DEGs, some genes involved in heat shock protein (HSPs), hypoxic adaptation, and glycolysis were up-regulated, while some genes involved in the ubiquitin-proteasome system (UPS) and fatty acid metabolism were down-regulated. Our results suggest that protein dynamic balance and function, hypoxia adaptation, and energy metabolism transformation are crucial in response to acute high-temperature stress. Our findings contribute to understanding the molecular response mechanism of T. ovatus under acute heat stress, which may provide some reference for studying the molecular mechanisms of other fish in response to heat stress.

Alteration of antioxidant status, glucose metabolism, and hypoxia signal pathway in Eirocheir sinensis after acute hypoxic stress and reoxygenation

Dissolved oxygen (DO) is crucial for the survival of Chinese mitten crab (Eirocheir sinensis); low DO levels adversely affect the health of these crabs. In this study, we evaluated the underlying response mechanism of E. sinensis to acute hypoxic stress by analyzing antioxidant parameters, glycolytic indicators, and hypoxia signaling factors. The crabs were exposed to hypoxia for 0, 3, 6, 12, and 24 h and reoxygenated for 1, 3, 6, 12, and 24 h. The hepatopancreas, muscle, gill, and hemolymph were sampled at different exposure times to detect the biochemical parameters and gene expression. The results showed that the activity of catalase, antioxidants, and malondialdehyde in tissues significantly increased under acute hypoxia and gradually decreased during the reoxygenation phase. Under acute hypoxic stress, glycolysis indices, including hexokinase (HK), phosphofructokinase, pyruvate kinase (PK), pyruvic acid (PA), lactate dehydrogenase (LDH), lactic acid (LA), succinate dehydrogenase (SDH), glucose, and glycogen in the hepatopancreas, hemolymph, and gills increased to varying degrees but recovered to the control levels after reoxygenation. Gene expression data showed that hypoxia signaling pathway-related genes, including hypoxia-inducible factor-1α/β (HIF1α/β), prolyl hydroxylase (PHD), factor inhibiting hypoxia-inducible factor (FIH), and glycolysis-related factors (HK and PK) were upregulated, showing that the HIF signaling pathway was activated under hypoxic conditions. In conclusion, acute hypoxic exposure activated the antioxidant defense system, glycolysis, and HIF pathway to respond to adverse conditions. These data contribute to elucidating the defense and adaptive mechanisms of crustaceans to acute hypoxic stress and reoxygenation.

Transcriptomic Analysis of Takifugu obscurus Gills under Acute Hypoxic Stress

Takifugu obscurus has relatively small gills and gill pores, leading to a relatively low respiratory capacity and increased vulnerability to low dissolved oxygen (DO) levels compared to other fish. To investigate the responses of T. obscurus to acute hypoxic stress, high-throughput-sequencing-based transcriptomic analyses were conducted here to assess the responses of T. obscurus gills to acute hypoxic stress. Three environmental conditions were compared including normoxia (DO: 7.0 ± 0.2 mg/L), hypoxic stress (DO: 0.9 ± 0.2 mg/L), and reoxygenation (4, 8, 12, and 24 h after return to normoxia) conditions to identify differentially expressed genes (DEGs) responsive to hypoxia. A total of 992, 877, 1561, 1412, and 679 DEGs were identified in the normoxia and reoxygenation for 4, 8, 12, and 24 h groups in comparison to the hypoxia groups, respectively. The DEGs were primarily associated with oxidative stress, growth and development, and immune responses. Further functional annotation enrichment analysis of the DEGs revealed that they were primarily related to cytokine-cytokine interactions, transforming growth factor β receptor (TGF-β), cell adhesion molecules (CAMs), the vascular endothelial growth factor (VEGF) signaling pathway, and the mitogen-activated protein kinase (MAPK) signaling pathway. These results provide new insights into the physiological and biochemical mechanisms of T. obscurus adaptations to hypoxic stress. Furthermore, these results provide a framework for future studies into the molecular mechanisms of hypoxia tolerance and the healthy culture of T. obscurus and other fish.

Molecular Characterization and Functional Analysis of Hypoxia-Responsive Factor Prolyl Hydroxylase Domain 2 in Mandarin Fish (Siniperca chuatsi)

With increased breeding density, the phenomenon of hypoxia gradually increases in aquaculture. Hypoxia is primarily mediated by the hypoxia-inducible factor 1 (HIF-1) signaling pathway. Prolyl hydroxylase domain proteins (PHD) are cellular oxygen-sensing molecules that regulate the stability of HIF-1α through hydroxylation. In this study, the characterization of the PHD2 from mandarin fish Siniperca chuatsi (scPHD2) and its roles in the HIF-1 signaling pathway were investigated. Bioinformation analysis showed that scPHD2 had the conserved prolyl 4-hydroxylase alpha subunit homolog domains at its C-terminal and was more closely related to other Perciformes PHD2 than other PHD2. Tissue-distribution results revealed that scphd2 gene was expressed in all tissues tested and more highly expressed in blood and liver than in other tested tissues. Dual-luciferase reporter gene and RT-qPCR assays showed that scPHD2 overexpression could significantly inhibit the HIF-1 signaling pathway. Co-immunoprecipitation analysis showed that scPHD2 could interact with scHIF-1α. Protein degradation experiment results suggested that scPHD2 could promote scHIF-1α degradation through the proteasome degradation pathway. This study advances our understanding of how the HIF-1 signaling pathway is regulated by scPHD2 and will help in understanding the molecular mechanisms underlying hypoxia adaptation in teleost fish.

Genomic adaptation of Ethiopian indigenous cattle to high altitude

The mountainous areas of Ethiopia represent one of the most extreme environmental challenges in Africa faced by humans and other inhabitants. Selection for high-altitude adaptation is expected to have imprinted the genomes of livestock living in these areas. Here we assess the genomic signatures of positive selection for high altitude adaptation in three cattle populations from the Ethiopian mountainous areas (Semien, Choke, and Bale mountains) compared to three Ethiopian lowland cattle populations (Afar, Ogaden, and Boran), using whole-genome resequencing and three genome scan approaches for signature of selection (iHS, XP-CLR, and PBS). We identified several candidate selection signature regions and several high-altitude adaptation genes. These include genes such as ITPR2, MB, and ARNT previously reported in the human population inhabiting the Ethiopian highlands. Furthermore, we present evidence of strong selection and high divergence between Ethiopian high- and low-altitude cattle populations at three new candidate genes (CLCA2, SLC26A2, and CBFA2T3), putatively linked to high-altitude adaptation in cattle. Our findings provide possible examples of convergent selection between cattle and humans as well as unique African cattle signature to the challenges of living in the Ethiopian mountainous regions.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Impact of high temperature, CO2 and parasitic infection on inflammation, immunodepression and programmed cell death in Colossoma macropomum at the transcriptional level

The production of tambaqui Colossoma macropomum has recently reached a milestone, being considered the main native species produced in South American continental waters. Despite the importance of this fish, its immunity is poorly understood, and global warming could pose severe risks to its health as increasing water temperature leads to an increase in the incidence of parasitic diseases. In an experimental context based on the high-emission scenario of the 5th Intergovernmental Panel on Climate Change (IPCC) report, we evaluated the synergistic effect of exposure to the extreme climate change scenario (RCP8.5) during two exposure periods (7 and 30 days) and two levels of parasitism by monogeneans (low and high). The goal was to understand how the tambaqui immune system will react to this challenge. To achieve this goal, we analyzed the expression of nine immunity-related genes (jak3, stat3, il-10, socs1, casp1, il-1β, tp53, bcl2, and hif-1α) in the spleen. Our main findings showed downregulation in the jak3/stat3 pathway, genes related to the control of inflammation and apoptosis, in addition to upregulation of proinflammatory genes and those related to pyroptosis during the first 7 days of exposure to the extreme climate scenario, also indicating a stage of immunodepression in these animals. After 30 days of exposure, all genes tended to return to similar levels in the current scenario, possibly due to the decrease in parasite load caused by chronic exposure to the extreme scenario. Our data strongly suggest that the increase in parasitism intensity caused by the extreme climate change scenario is responsible for disturbances in the host's immune system. However, more studies are needed to clarify this poorly understood cascade of events.

Effect of secondary attack by Aeromonas hydrophila on the expression level of hif genes in common carp (Cyprinus carpio)

Hypoxia-inducible factors (hifs) are involved in infectious diseases inflammatory reactions, and immune regulation. Common carp, a representative allotetraploid species that has undergone genome-wide replication events, has important research value. In this study, common carp were infected twice with Aeromonas hydrophila. Liver tissues of common carp were collected at 4 h, 12 h, 24 h, 48 h, 3 days, 7 days post-first infection and 4 h, 12 h, 24 h post-second infection. The mRNA levels of hif genes were determined at different time points. The hif2a-2, hif3a-2, hif3b-1 and hif3b-2 expression levels in the infected group were upregulated when compared with those in the control group, whereas the expression levels of other genes were downregulated after the second infection. This indicates that the effect of A. hydrophila infection on gene expression pattern is dependent on the host, pathogen, infected tissue and gene. Pressure analysis of the hif gene family revealed that the non-synonymous substitution to synonymous substitution ratio of 12 hif genes was <1, which indicated that they were in a state of purification and selection. Combined with the differences between copy genes, the polyclonal antibodies against Hif1b-1 and Hif1b-2 were successfully prepared in this study. Western blot analysis showed that the protein expression of Hif1b-1 and Hif1b-2 reached to the highest level 48 h after the first infection. After the second A. hydrophila infection, the protein expression levels of Hif1b-1 and Hif1b-2 reached the highest levels at 4 and 48 h, respectively. This may indicate that the Hif1b-1 and Hif1b-2 genes in common carp play an important role in the immune mechanism at the protein level. The findings of this study will lay the foundation for future studies on the immune regulatory function of common carp hif genes, which may aid in devising novel therapeutic strategies for common carp diseases, such as A. hydrophila infection.

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.

Hypoxia- and hyperoxia-related gene expression dynamics during developmental critical windows of the tropical gar Atractosteus tropicus

Aquatic hypoxia is both a naturally-occurring and anthropogenically-generated event. Fish species have evolved different adaptations to cope with hypoxic environments, including gill modifications and air breathing. However, little is known about the molecular mechanisms involved in the respiration of embryonic and larval fishes during critical windows of development. We assessed expression of the genes hif-1α, fih-1, nhe1, epo, gr and il8 using the developing tropical gar as a piscine model during three developmental periods (fertilization to hatch, 1 to 6 days post hatch (dph) and 7 to 12 dph) when exposed to normoxia (~7.43 mg/L DO), hypoxia (~2.5 mg/L DO) or hyperoxia (~9.15 mg/L DO). All genes had higher expression when fish were exposed to either hypoxia or hyperoxia during the first two developmental periods. However, fish continuously exposed to hypoxia had increased expression of the six genes by hatching and 6 dph, and by 12 dph only hif-1α still had increased expression. The middle developmental period was the most hypoxia-sensitive, coinciding with several changes in physiology and morphology. The oldest larvae were the most resilient to gene expression change, with little variation in expression of the six genes compared. This study is the first to relate the molecular response of an air-breathing fish to oxygen availability to developmental critical windows and contributes to our understanding of some molecular responses of developing fish to changes in oxygen availability.

Effects of acute hypoxia and reoxygenation on oxygen sensors, respiratory metabolism, oxidative stress, and apoptosis in hybrid yellow catfish "Huangyou-1"

The regulation mechanism of the hybrid yellow catfish "Huangyou-1" was assessed under conditions of hypoxia and reoxygenation by examination of oxygen sensors and by monitoring respiratory metabolism, oxidative stress, and apoptosis. The expressions of genes related to oxygen sensors (HIF-1α, HIF-2α, VHL, HIF-1β, PHD2, and FIH-1) were upregulated in the brain and liver during hypoxia, and recovered compared with control upon reoxygenation. The expressions of genes related to glycolysis (HK1, PGK1, PGAM2, PFK, and LDH) were increased during hypoxia and then recovered compared with control upon reoxygenation. The mRNA levels of CS did not change during hypoxia in the brain and liver, but increased during reoxygenation. The mRNA levels of SDH decreased significantly only in the liver during hypoxia, but later increased compared with control upon reoxygenation in both tissues. Under hypoxic conditions, the expressions of genes related to oxidative stress (SOD1, SOD2, GSH-Px, and CAT) and the activity of antioxidant enzymes (SOD, CAT, and GSH-Px) and MDA were upregulated compared with control. The expressions of genes related to apoptosis (Apaf-1, Bax, Caspase 3, Caspase 9, and p53) were higher than those in control during hypoxic exposure, while the expressions of Bcl-2 and Cyt C were decreased. The findings of the transcriptional analyses will provide insights into the molecular mechanisms of hybrid yellow catfish "Huangyou-1" under conditions of hypoxia and reoxygenation. Overall, these findings showed that oxygen sensors of "Huangyou-1" are potentially useful biomarkers of environmental hypoxic exposure. Together with genes related to respiratory metabolism, oxidative stress and apoptosis occupy a quite high position in enhancing hypoxia tolerance. Our findings provided new insights into the molecular regulatory mechanism of hypoxia in "Huangyou-1."

The evolutionary and physiological significance of the Hif pathway in teleost fishes

The hypoxia-inducible factor (HIF) pathway is a key regulator of cellular O2 homeostasis and an important orchestrator of the physiological responses to hypoxia (low O2) in vertebrates. Fish can be exposed to significant and frequent changes in environmental O2, and increases in Hif-α (the hypoxia-sensitive subunit of the transcription factor Hif) have been documented in a number of species as a result of a decrease in O2. Here, we discuss the impact of the Hif pathway on the hypoxic response and the contribution to hypoxia tolerance, particularly in fishes of the cyprinid lineage, which includes the zebrafish (Danio rerio). The cyprinids are of specific interest because, unlike in most other fishes, duplicated paralogs of the Hif-α isoforms arising from a teleost-specific genome duplication event have been retained. Positive selection has acted on the duplicated paralogs of the Hif-α isoforms in some cyprinid sub-families, pointing to adaptive evolutionary change in the paralogs. Thus, cyprinids are valuable models for exploring the evolutionary significance and physiological impact of the Hif pathway on the hypoxic response. Knockout in zebrafish of either paralog of Hif-1α greatly reduces hypoxia tolerance, indicating the importance of both paralogs to the hypoxic response. Here, with an emphasis on the cardiorespiratory system, we focus on the role of Hif-1α in the hypoxic ventilatory response and the regulation of cardiac function. We explore the effects of the duration of the hypoxic exposure (acute, sustained or intermittent) on the impact of Hif-1α on cardiorespiratory function and compare relevant data with those from mammalian systems.

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.

Transcriptome analysis of hepatopancreas from different living states oriental river prawn (Macrobrachium nipponense) in response to hypoxia

As an important economical freshwater prawn, Macrobrachium nipponense has difficulty with adapting to hypoxia. In this study, comparative transcriptome analysis was used for the first time to explore the differences between different living states of Macrobrachium nipponense under hypoxia. A total of 94.22 Gb clean reads were obtained and assembled into 54,688 unigenes. A total of 224, 266, and 750 differently expressed genes were found in the comparison of the control and death groups, the control and moribund groups, and the control and survived groups, respectively. Three signal pathways closely related to hypoxia were found by enriching of the signal pathways in three comparison groups. In addition, much attention was focused on the differential genes in these pathways. Oxidative stress related genes, such as 70 kDa heat shock protein, phosphoenolpyruvate carboxykinase and cyclooxygenase were differentially expressed in different comparisons. After comparing with previous studies, cyclooxygenase was found to be an important hypoxia-related gene that is fully involved in the hypoxic response. Interestingly, two new genes with no Nr annotation were found in this manuscript. This manuscript will enrich our understanding of oxidative stress response to hypoxia and provide a theoretical basis for the subsequent solution of apoptosis caused by hypoxia.

Transcriptome sequencing provides insights into the mechanism of hypoxia adaption in bighead carp (Hypophthalmichthys nobilis)

Hypoxia negatively affects the behavior, immunology, physiology, and growth of fish. Therefore, uncovering the genetic mechanisms underlying hypoxia adaptation and tolerance in fish prior to any genetic improvement is essential. Bighead carp is one of the most important freshwater fish species in aquaculture worldwide; however, this species does not have a strong ability to tolerate hypoxia. In this study, the dissolved oxygen level (0.6 mg/L) was maintained above the asphyxiation point of bighead carp for a long time to simulate hypoxia stress. The liver, gills, and heart were sampled before (0 h) and after (1 h, 2 h, 4 h) the hypoxia tests. Glutathione peroxidase (GPx) and catalase (CAT) activities and malondialdehyde (MDA) levels in the liver were significantly (p < 0.05) elevated at 1 h after hypoxic stress. By observing tissue morphology, the cell structure of the liver and gill tissues was found to change to varying degrees before and after hypoxia stress. Transcriptome sequencing was performed on 36 samples of gill, liver, and heart at four time points, and a total of 293.55G of data was obtained. In the early phase (0-1 h), differentially expressed genes (DEGs, 807 genes upregulated, 654 genes downregulated) were mainly enriched in signal transduction, such as cytokine-cytokine receptor interactions and ECM-receptor interactions. In the middle phase (0-2 h), DEGs (1201 genes upregulated and 2036 genes downregulated) were mainly enriched in regulation and adaptation, such as the MAPK and FoxO signaling pathways. Finally, in the later phase (0-4 h), DEGs (3975 genes upregulated and 4412 genes downregulated) were mainly enriched in tolerance and apoptosis, such as the VEGF signaling pathway and apoptosis. The genes with the most remarkable upregulation at different time points in the three tissues had some similarities. Genetic differences in these genes may be responsible for the differences in hypoxia tolerance among individuals. Altogether, our study provides new insights into the molecular mechanisms of hypoxia adaptation in fish. Further, the key regulatory genes identified provide genetic resources for breeding hypoxia-tolerant bighead carp species.

Hypoxia inducible factor 1-α is minimally involved in determining the time domains of the hypoxic ventilatory response in adult zebrafish (Danio rerio)

In the current study, adult zebrafish (Danio rerio) were exposed to 72 h hypoxia (90 mmHg) to assess the time domains of the hypoxia ventilatory response (HVR) and the consequence on a subsequent more severe (40 mmHg) bout of acute hypoxia. Experiments were performed on wild-type fish and mutants in which one or both paralogs of hypoxia inducible factor-1α (hif-1α) were knocked out. Although there were subtle differences among the wild-type and knockout fish, resting f_V was reestablished after 2-8 h of continuous hypoxia in both groups, a striking example of hypoxic ventilatory decline (HVD). When fish were subsequently exposed to more severe hypoxia, a rapid increase in f_V was observed, the magnitude of which was independent of genotype or prior exposure history. During recovery, fish that had been exposed to 72 h of 90 mmHg hypoxia exhibited a pronounced undershoot in f_V, which was absent in the hif-1α double knockouts. Overall, the results revealed distinct time domains of the HVR in zebrafish that were largely Hif-1α-independent.

Contrasting strategies of hypoxic cardiac performance and metabolism in cichlids and armoured catfish

The heart of tropical fishes is a particularly useful model system in which to investigate mechanisms of hypoxic tolerance. Here we focus on insights gained from two groups of fishes, cichlids and armoured catfishes. Cichlids respond to hypoxia by entering a sustained hypometabolism with decreased heart performance to match whole animal circulatory needs. Heart rate is decreased along with protein turnover to reduce adenosine triphosphate demand. This occurs despite the inherent capacity for high levels of cardiac power development. Although highly hypoxic tolerant at the whole animal level, the heart of cichlids does not have high constitutive activities of glycolytic enzymes compared to other species. Information is conflicting with respect to changes in glycolytic gene expression and enzyme activity following hypoxic exposure with some studies showing increases and others decreases. In contrast to cichlids, species of armoured catfish, that are routinely exposed to water of low oxygen content, do not display hypoxic bradycardia. Under hypoxia there are early changes in glucose trafficking suggestive of activation of glycolysis before lactate accumulation. Thereafter, heart glycogen is mobilized and lactate accumulates in both heart and blood, in some species to very high levels. Heart performance under hypoxia is enhanced by defense of intracellular pH. A functional sarcoplasmic reticulum and binding of hexokinase to the outer mitochondrial membrane may also play a role in cardioprotection. Maintenance of heart performance under hypoxia may relate to a tradeoff between air breathing via a modified stomach and circulatory demands for digestion.

Molecular cloning and functional characterization of the hypoxia-inducible factor-1α in bighead carp (Aristichthys nobilis)

HIF-l is the earliest documented and most widely studied hypoxia-inducible factor (HIF) and plays a key role in the cell hypoxia signal transduction pathway. Particularly, the HIF-1α protein is sensitive to oxygen and plays a critical role in hypoxia regulation. This study is the first to report on the molecular cloning and characterization of HIF-1α in bighead carp (Aristichthys nobilis; anHIF-1α). The full-length cDNA of anHIF-1α was 2361 bp, and encodes an estimated 674 amino acids with a predicted molecular mass of 76.10 kDa and a theoretical isoelectric point of 7.72. Moreover, the conserved basic Helix-Loop-Helix domain along with two Per-ARNT-Sim domains (A/B), and C-TAD were identified in this protein. Interestingly, the tertiary structure of the anHIF-1α protein was found to be extremely similar to that of mice. Multiple comparison and phylogenetic tree results demonstrated that anHIF-1α was highly conserved. Under normoxic conditions, anHIF-1α mRNA transcripts could be detected in all tissues examined with the highest expression level in the heart. With gradually decreasing oxygen concentrations, anHIF-1α mRNA level was upregulated significantly in the gill, liver, kidney, spleen, intestine, brain, and muscle tissues (P < 0.05). Similarly, anHIF-1α was expressed in all examined bighead carp tissues, and the results suggested that the upregulation of anHIF-1α at the transcriptional level may be an important stress response adaptation to hypoxia in bighead carp. Finally, based on the tertiary structure comparative analyses between anHIF-1α with mouse HIF-1α, we think the physiological function, and protein structure of HIF-1α could be compared between fish and mammal in the future.

Insights into early ontogenesis: characterization of stress and development key genes of pikeperch (Sander lucioperca) in vivo and in vitro

There are still numerous difficulties in the successful farming of pikeperch in the anthropogenic environment of various aquaculture systems, especially during early developmental steps in the hatchery. To investigate the physiological processes involved on the molecular level, we determined the basal expression patterns of 21 genes involved in stress and immune responses and early ontogenesis of pikeperch between 0 and 175 days post hatch (dph). Their transcription patterns most likely reflect the challenges of growth and feed conversion. The gene coding for apolipoprotein A (APOE) was strongly expressed at 0 dph, indicating its importance for yolk sac utilization. Genes encoding bone morphogenetic proteins 4 and 7 (BMP4, BMP7), creatine kinase M (CKM), and SRY-box transcription factor 9 (SOX9) were highly abundant during the peak phases of morphological changes and acclimatization processes at 4-18 dph. The high expression of genes coding for peroxisome proliferator-activated receptors alpha and delta (PPARA, PPARD) at 121 and 175 dph, respectively, suggests their importance during this strong growth phase of juvenile stages. As an alternative experimental model to replace further in vivo investigations of ontogenetically important processes, we initiated the first approach towards a long-lasting primary cell culture from whole pikeperch embryos. The present study provides a set of possible biomarkers to support the monitoring of pikeperch farming and provides a first basis for the establishment of a suitable cell model of this emerging aquaculture species.

Comparative Transcriptome Analysis of Gill Tissue in Response to Hypoxia in Silver Sillago (Sillago sihama)

Silver sillago (Sillago sihama) is a commercially important marine fish species in East Asia. In this study, we compared the transcriptome response to hypoxia stress in the gill tissue of S. sihama. The fish were divided into four groups, such as 1 h of hypoxia (hypoxia1h, DO = 1.5 ± 0.1 mg/L), 4 h of hypoxia (hypoxia4h, DO = 1.5 ± 0.1 mg/L), 4 h of reoxygen (reoxygen4h, DO = 8.0 ± 0.2 mg/L) after 4 h of hypoxia (DO = 1.5 mg/L), and normoxia or control (DO = 8.0 ± 0.2 mg/L) groups. Compared to the normoxia group, a total of 3550 genes were identified as differentially expressed genes (DEGs) (log₂foldchange > 1 and padj < 0.05), including 1103, 1451 and 996 genes in hypoxia1h, hypoxia4h and reoxygen4h groups, respectively. Only 247 DEGs were differentially co-expressed in all treatment groups. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, DEGs were significantly enriched in steroid biosynthesis, biosynthesis of amino acids, glutathione metabolism and metabolism of xenobiotics by cytochrome P450, ferroptosis and drug metabolism-cytochrome P450 pathways. Of these, the cytochrome P450 (CYP) and glutathione S-transferase (GST) gene families were widely expressed. Our study represents the insights into the underlying molecular mechanisms of hypoxia stress.

The HIF1αn gene and its association with hypoxia tolerance in the Asian seabass

Hypoxia is one of the major challenges in aquaculture industry. Breeding of fish tolerant to hypoxia is an important task in genetic improvement of aquaculture species. The Asian seabass, Lates calcarifer, is an important foodfish species. We identified and characterized the hypoxia-inducible factor inhibitor (HIF1αn) gene in the Asian seabass. The full-length cDNA sequence of the HIF1αn was 3425 bp, with an ORF of 1065 bp, encoding 354 amino acids. The genomic sequence of the gene was 8667 bp in length, and contained eight exons and seven introns. Phylogenetic analysis of the gene in fish and tetrapods revealed that the HIF1αn in the Asian seabass was closely related to that of tilapia (Oreochromis niloticus). The HIF1αn was highly up-regulated in the gill, spleen and heart after 3.5-hours hypoxia treatment. We identified three SNPs in the third and fourth introns of the HIF1αn gene. The SNP (i.e. SNP 9332241 (C/T)) in the fourth intron was significantly (P < 0.01) associated with hypoxia tolerance. This SNP might be useful in selecting Asian seabass for improved tolerance to hypoxia. Our data also provide useful information for further detailed analysis of the function of the HIF1αn gene in hypoxia tolerance.

MicroRNA-124 regulates lactate transportation in the muscle of largemouth bass (micropterus salmoides) under hypoxia by targeting MCT1

Carbohydrate metabolism switches from aerobic to anaerobic (glycolysis) to supply energy in response to acute hypoxic stress. Acute hypoxic stress with dissolved oxygen (DO) levels of 1.2 ± 0.1 mg/L for 24 h and 12 h re-oxygenation was used to investigate the response of the anaerobic glycolytic pathway in Micropterus salmoides muscle. The results showed that the glucose concentration was significantly lower in muscle, while the lactic acid and pyruvic acid concentrations tended to increase during hypoxic stress. No significant difference was observed in muscle glycogen, and ATP content fluctuated significantly. The activities of gluconeogenesis-related enzymes were slightly elevated, such as phosphoenolpyruvate carboxykinase (PEPCK). The activities of the glycolytic enzymes increased after the induction of hypoxia, such as hexokinase (HK), pyruvate kinase (PK), and lactate dehydrogenase (LDH). Curiously, phosphofructokinase (PFK) activity was significantly down-regulated within 4 h during hypoxia, although these effects were transient, and most indices returned to control levels after 12 h of re-oxygenation. Upregulated hif-1α, ampkα, hk, glut1, and ldh mRNA expression suggested that carbohydrate metabolism was reprogrammed under hypoxia. Lactate transport was regulated by miR-124-5p according to quantitative polymerase chain reaction and dual luciferase reporter assays. Our findings provide new insight into the molecular regulatory mechanism of hypoxia in Micropterus salmoides muscle.

Effects of dietary supplementation of three strains of Lactococcus lactis on HIFs genes family expression of the common carp following Aeromonas hydrophila infection

HIFs (Hypoxia inducible factors) are the main regulators of the expression change of oxygen-dependent genes, in addition, they also play important roles in immune regulation. HIFs participate in infectious diseases and inflammatory responses, providing us a new therapeutic target for the treatment of diseases. In this study, 16 HIFs were identified in common carp genome database. Comparative genomics analysis showed large expansion of HIF gene family and approved the four round whole genome duplication (WGD) event in common carp. To further understand the function of HIFs, the domain architectures were predicted. All HIF proteins had the conserved HLH-PAS domain, which were essential for them to form dimer and bind to the downstream targets. The differences in domain of HIFα and HIFβ might result in their different functions. Phylogenetic analysis revealed that all HIFs were divided into two subfamilies and the HIFs in common carp were clustered with their teleost counterparts indicating they are highly conservative during evolution. In addition, the tissue distribution was examined by RT-PCR showed that most of HIF genes had a wide range of tissue distribution but exhibited tissue-specific expression patterns. The expression divergences were observed between the copy genes, for example, HIF1A-1, HIF2A-1, ARNT-2 had wide tissue distribution while their copies had limited tissue distribution, proving the function divergence of copies post the WGD event. In order to find an effective activation of HIFs and apply to treatment of aquatic diseases, we investigate the dietary supplementation effects of different strains of Lactococcus lactis on the expression of HIFα subfamily members in kidney of common carp infected with A. hydrophila. In addition, all of the HIF genes have a high expression in the early stages of infection, and decreased in the treatment time point of 48 h in common carp. This phenomenon confirms that as a switch, the main function of HIFs is to regulate the production of immune response factors in early infection. So activation of the switch may be an effective method for infectious disease treatment. As expected, the treatment groups improved the expression of HIFs compared with the control group, and the effects of the three strains are different. The strain1 of L. lactis had a stronger induction on HIF genes than strain2 and strain3, and it might be applied as a potential activation of HIF genes for disease treatment. So, adding befitting L. lactis maybe a well method to activate the HIF genes to protect them from mycobacterial infection.

Gene expression, genotoxicity, and physiological responses in an Amazonian fish, Colossoma macropomum (CUVIER 1818), exposed to Roundup® and subsequent acute hypoxia

Roundup® (RD) is a glyphosate-based herbicide used to control weeds in agriculture, and fishponds. In the Amazon, hypoxia is a natural phenomenon in flooded areas. Beyond the challenge of hypoxia, fish need to cope with the use of pesticides as RD that increases in the aquatic environment through the leaching of agricultural areas, and in aquaculture fish tanks. Thus, there is a need to better understand the combined effects of hypoxia and RD contamination for aquatic biota. The aim of this study was to investigate the effects of Roundup® (RD) and subsequent acute hypoxia in the gene expression, genotoxicity, histological and physiological responses of Colossoma macropomum. Fish were individually exposed to four different treatments during 96 h: normoxia (N), hypoxia (H), RD plus normoxia (NRD), and RD plus hypoxia (HRD) (RD concentration represents 75% of LC50 - nominal concentration 15 mg L-1 to C. macropomum). HRD fishes presented down-regulation of hif-1α gene and ras oncogene, while NRD fish presented overexpression of ras; no difference occurred in hif-1α gene expression in both normoxia treatments. The glutathione-S-transferase and catalase activities increased in HRD fish liver compared to NRD. Otherwise, there was no difference in lipoperoxidation (LPO) between all treatments. Genetic Damage Index, measured throughout comet assay in erythrocytes of all treatments, presented similar values, excepted by fish exposed to NRD. As regard as hypoxic exposure, hypoxic fish presented significantly lower values, compared to HRD fishes. An increase in liver histological injuries occurred in H and HRD fish groups. In conclusion, we may affirm that C. macropomum is sensitive concerning RD contamination and that this sensitivity increases when combined with hypoxia.

The conservation and functionality of the oxygen-sensing enzyme Factor Inhibiting HIF (FIH) in non-vertebrates

The asparaginyl hydroxylase, Factor Inhibiting HIF (FIH), is a cellular dioxygenase. Originally identified as oxygen sensor in the cellular response to hypoxia, where FIH acts as a repressor of the hypoxia inducible transcription factor alpha (HIF-α) proteins through asparaginyl hydroxylation, FIH also hydroxylates many proteins that contain ankyrin repeat domains (ARDs). Given FIH's promiscuity and the unclear functional effects of ARD hydroxylation, the biological relevance of HIF-α and ARD hydroxylation remains uncertain. Here, we have employed evolutionary and enzymatic analyses of FIH, and both HIF-α and ARD-containing substrates, in a broad range of metazoa to better understand their conservation and functional importance. Utilising Tribolium castaneum and Acropora millepora, we provide evidence that FIH from both species are able to hydroxylate HIF-α proteins, supporting conservation of this function beyond vertebrates. We further demonstrate that T. castaneum and A. millepora FIH homologs can also hydroxylate specific ARD proteins. Significantly, FIH is also conserved in several species with inefficiently-targeted or absent HIF, supporting the hypothesis of important HIF-independent functions for FIH. Overall, these data show that while oxygen-dependent HIF-α hydroxylation by FIH is highly conserved in many species, HIF-independent roles for FIH have evolved in others.

Molecular characterization and expression regulation of the factor-inhibiting HIF-1 (FIH-1) gene under hypoxic stress in bighead carp (Aristichthys nobilis)

Factor-inhibiting HIF-1 (FIH-1) is an asparagine hydroxylase that interacts with hypoxia-inducible factor 1α (HIF-1α) to regulate transcriptional activity of HIF-1. Few studies of fish FIH-1 have been reported to date. In this study, the cDNA of FIH-1 gene was cloned and characterized for bighead carp, Aristichthys nobilis (AnFIH-1). The AnFIH-1 cDNA is 2065 bp in length, encoding a protein of 357 amino acid (aa) residues, which contains a JmjC homology region of the jumonji transcription factors. AnFIH-1 shares high identities with other vertebrate FIH-1 (79.1-96.4%), especially in the JmjC homology region, suggesting its conserved function. During the embryonic stages of A. nobilis, AnFIH-1 had significantly high expression levels in unfertilized egg and blastula. In healthy tissues, its predominant mRNA expression was detected in muscle. The mRNA levels of AnFIH-1 were significantly upregulated in the liver, gill, hypothalamus, and spleen after hypoxic treatment, and then decreased to pretreatment levels after 6-h re-oxygenation. However, in the muscle, continual increasing of mRNA expression was observed after hypoxic shock and re-oxygenation. These results indicate that FIH-1 may play an important role in physiological regulation for adapting to hypoxia stress in A. nobilis.

The improved energy metabolism and blood oxygen-carrying capacity for pufferfish, Takifugu fasciatus, against acute hypoxia under the regulation of oxygen sensors

Hypoxia frequently occurs in aquatic ecosystem, which is influenced by salinity, water temperature, weather, and surface water runoff. In order to shed further light on the evolutionary and adaptive mechanisms in fish under hypoxic condition, the impact of acute hypoxia (1.63 ± 0.2 mg/L) and reoxygenation (7.0 ± 0.3 mg/L) on oxygen sensors, energy metabolism, and hematological indices was evaluated in Takifugu fasciatus. Data from transcriptional level analysis show that the expressions of genes related to oxygen sensors (HIF-1α, PHD2, and VHL) were upregulated in the brain and liver under hypoxia and recovered under reoxygenation. The upregulation of GLUT2, VEGF-A, and EPO in conjugation with VEGF-A protein and hematological indices conferred the rapid adjustments of cellular glucose uptake and blood oxygen-carrying capacities in pufferfish. Higher levels of glycolysis-related mRNAs (HK, PGK1, and PGAM2), HK activity, and proteins (PGK1 and PGAM2) were detected in the brain and liver under hypoxic condition compared with control. Interestingly, the expression of MDH1 at the mRNA, enzyme activity, and protein levels was significantly increased in the brain at 0 or 2 h and in the liver at 8 h under hypoxic condition. In addition, although the enzyme activity and mRNA expression of LDH in the brain were not significantly changed, a persistent upregulation was observed in the liver during hypoxia exposure. This study demonstrated that pufferfish could counterpoise the energetic demands and hematological functional properties evoked by oxygen sensors after hypoxia. Our findings provided new insights into the molecular regulatory mechanism of hypoxia in pufferfish.

Hypoxia alters the expression of hif-1a mRNA and downstream HIF-1 response genes in embryonic and larval lake whitefish (Coregonus clupeaformis)

Lake whitefish (Coregonus clupeaformis) embryos and larvae were exposed to hypoxia at different developmental ages to determine when the cellular response to hypoxia could be initiated. mRNA levels of hypoxia-inducible factor 1α (hif-1α), hsp70, and several HIF-1 target genes were quantified in embryos at 21, 38, 63, 83- and 103-days post fertilisation (dpf) and in larvae at 1, 2, 3- and 4-weeks post hatch (wph) following a 6-hour hypoxia exposure. hsp70 mRNA levels were increased in response to hypoxia at all embryonic ages. By comparison, the first observed change in hif-1α mRNA in response to hypoxia was at 38 dpf, where it was down-regulated from high basal levels, with this response persisting through to 83 dpf. Interestingly, this decrease in hif-1α mRNA coincided with increases in the mRNA levels of the HIF-1 target genes: vegfa (vascular endothelial growth factor A), igfbp1 (insulin-like growth factor binding protein 1), ldha (lactate dehydrogenase a), gapdh (glyceraldehyde-3-phosphate dehydrogenase) and epo (erythropoietin) at select ages. Collectively, this suggests a possible HIF-1-mediated response to hypoxia despite a decrease in hif-1α mRNA. Coinciding with a decrease in basal levels, increases in hif-1α were measured in response to hypoxia at 103 dpf and in larval fish at 1, 2 and 3 wph but there were no consistent increases in HIF-1 target genes at these ages. Overall, our findings indicate that lake whitefish can mount a response to hypoxia early in embryogenesis which may mitigate some of the damaging effects of exposure to low oxygen levels at these critical life history stages.

A Review of Molecular Responses of Catfish to Bacterial Diseases and Abiotic Stresses

Catfish is one of the major aquaculture species in the United States. However, the catfish industry is threatened by several bacterial diseases such as enteric septicemia of catfish (ESC), columnaris disease and Aeromonas disease, as well as by abiotic stresses such as high temperature and low oxygen. Research has been conducted for several decades to understand the host responses to these diseases and abiotic stresses. With the development of sequencing technologies, and the application of genome-wide association studies in aquaculture species, significant progress has been made. This review article summarizes recent progress in understanding the molecular responses of catfish after bacterial infection and stress challenges, and in understanding of genomic and genetic basis for disease resistance and stress tolerance.

Hypoxia-inducible transcription factors in fish: expression, function and interconnection with the circadian clock

The hypoxia-inducible transcription factors are key regulators for the physiological response to low oxygen availability. In vertebrates, typically three Hif-α isoforms, Hif-1α, Hif-2α and Hif-3α, are expressed, each of which, together with Hif-1β, may form a functional heterodimer under hypoxic conditions, controlling expression of hundreds of genes. A teleost-specific whole-genome duplication complicates the analysis of isoform-specific functions in fish, but recent studies suggest that the existence of paralogues of a specific isoform opens up the possibility for a subfunctionalization. In contrast to during development inside the uterus, fish eggs are freely accessible and studies analyzing Hif expression in fish embryos during development have revealed that Hif proteins are not only controlling the hypoxic response, but are also crucial for proper development and organ differentiation. Significant advances have been made in our knowledge about tissue-specific functions of Hif proteins, especially with respect to gill or gonadal tissue. The hypoxia signalling pathway is known to be tightly and mutually intertwined with the circadian clock in zebrafish and mammals. Recently, a mechanistic explanation for the hypoxia-induced dampening of the transcriptional clock was detected in zebrafish, including also metabolically induced alterations of cellular redox signalling. In turn, MAP kinase-mediated H2O2 signalling modulates the temporal expression of Hif-1α protein, similar to the redox regulation of the circadian clock itself. Once again, the zebrafish has emerged as an excellent model organism with which to explore these specific functional aspects of basic eukaryotic cell biology.

Comparative transcriptome analysis of the swimbladder reveals expression signatures in response to low oxygen stress in channel catfish, Ictalurus punctatus

Channel catfish is the leading aquaculture species in the US, and one of the reasons for its application in aquaculture is its relatively high tolerance against hypoxia. However, hypoxia can still cause huge economic losses to the catfish industry. Studies on hypoxia tolerance, therefore, are important for aquaculture. Fish swimbladder has been considered as an accessory respiration organ surrounded by a dense capillary countercurrent exchange system. In this regard, we conducted RNA-Seq analysis with swimbladder samples of catfish under hypoxic and normal conditions to determine if swimbladder was responsive to low oxygen treatment and to reveal genes, their expression patterns, and pathways involved in hypoxia responses in catfish. A total of 155 differentially expressed genes (DEGs) were identified from swimbladder of adult catfish, whereas a total of 2,127 DEGs were identified from swimbladder of fingerling catfish under hypoxic condition as compared with untreated controls. Subsequent pathway analysis revealed that many DEGs under hypoxia were involved in HIF signaling pathway ( nos2, eno2, camk2d2, prkcb, cdkn1a, eno1, and tfrc), MAPK signaling pathway (voltage-dependent calcium channel subunit genes), PI3K/Akt/mTOR signaling pathway ( itga6, g6pc, and cdkn1a), Ras signaling pathway ( efna3 and ksr2), and signaling by VEGF ( fn1, wasf3, and hspb1) in catfish swimbladder. This study provided insights into regulation of gene expression and their involved gene pathways in catfish swimbladder in response to low oxygen stresses.

The CC and CXC chemokine receptors in channel catfish (Ictalurus punctatus) and their involvement in disease and hypoxia responses

Chemokines are vital regulators of cell mobilization for immune surveillance, inflammation, and development. Chemokines signal through binding to their receptors that are a superfamily of seven-transmembrane domain G-coupled receptors. Recently, a complete repertoire of both CC and CXC chemokines have been identified in channel catfish, but nothing is known about their receptors. In this study, a set of 29 CC chemokine receptor (CCR) genes and 8 CXC chemokine receptor (CXCR) genes were identified and annotated from the channel catfish genome. Extensive phylogenetic and comparative genomic analyses were conducted to annotate these genes, revealing fish-specific CC chemokine receptors, and lineage-specific tandem duplications of chemokine receptors in the teleost genomes. With 29 genes, the channel catfish genome harbors the largest numbers of CC chemokine receptors among all the genomes characterized. Analysis of gene expression after bacterial infections indicated that the chemokine receptors were regulated in a gene-specific manner. Most differentially expressed chemokine receptors were up-regulated after Edwardsiella ictaluri and Flavobacterium columnare infection. Among which, CXCR3 and CXCR4 were observed to participate in immune responses to both bacterial infections, indicating their potential roles in catfish immune activities. In addition, CXCR3.2 was significantly up-regulated in ESC-susceptible fish, and CXCR4b was mildly induced in ESC-resistant fish, further supporting the significant roles of CXCR3 and CXCR4 in catfish immune responses. CXCR4b and CCR9a were both up-regulated not only after bacterial infection, but also after hypoxia stress, providing the linkage between bacterial infection and low oxygen stresses. These results should be valuable for comparative immunological studies and provide insights into their roles in disease and stress responses.

Dynamic mRNA and miRNA expression analysis in response to hypoxia and reoxygenation in the blunt snout bream (Megalobrama amblycephala)

Adaptation to hypoxia is a complex process involving various pathways and regulation mechanisms. A better understanding of the genetic influence on these mechanisms could permit selection for hypoxia-sensitive fish. To aid this understanding, an integrated analysis of miRNA and mRNA expression was performed in Megalobrama amblycephala under four acute hypoxia and reoxygenation stages. A number of significantly differentially-expressed miRNAs and genes associated with oxidative stress were identified, and their functional characteristics were revealed by GO function and KEGG pathway analysis. They were found to be involved in HIF-1 pathways known to affect energy metabolism and apoptosis. MiRNA-mRNA interaction pairs were detected from comparison of expression between the four different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating hypoxia stress. As a unique resource for gene expression and regulation during hypoxia and reoxygenation, this study could provide a starting point for further studies to better understand the genetic background of hypoxia stress.

Genome-Wide Association Study Reveals Multiple Novel QTL Associated with Low Oxygen Tolerance in Hybrid Catfish

Hypoxic condition is common in aquaculture, leading to major economic losses. Genetic analysis of hypoxia tolerance, therefore, is not only scientifically significant, but also economically important. Catfish is generally regarded as being highly tolerant to low dissolved oxygen, but variations exist among various populations, strains, and species. In this study, we conducted a genome-wide association study (GWAS) using the catfish 250 K SNP array to identify quantitative trait locus (QTL) associated with tolerance to low dissolved oxygen in the channel catfish × blue catfish interspecific system. Four linkage groups (LG2, LG4, LG23, and LG29) were found to be associated with low oxygen tolerance in hybrid catfish. Multiple significant SNPs were found to be physically linked in genomic regions containing significant QTL for low oxygen tolerance on LG2 and LG23, and in those regions containing suggestively significant QTL on LG2, LG4, LG23, and LG29, suggesting that the physically linked SNPs were genuinely segregating and related with low oxygen tolerance. Analysis of genes within the associated genomic regions suggested that many of these genes were involved in VEGF, MAPK, mTOR, PI3K-Akt, P53-mediated apoptosis, and DNA damage checkpoint pathways. Comparative analysis indicated that most of the QTL at the species level, as analyzed by using the interspecific system, did not overlap with those identified from six strains of channel catfish, confirming the complexity of the genetic architecture of hypoxia tolerance in catfish.

The chemokinome superfamily: II. The 64 CC chemokines in channel catfish and their involvement in disease and hypoxia responses

Chemokines are a superfamily of structurally related chemotactic cytokines exerting significant roles in regulating cell migration and activation. Based on the arrangement of the first four cysteine residues, they are classified into CC, CXC, C and CX3C subfamilies. In this study, a complete set of 64 CC chemokine ligand (CCL) genes was systematically identified, annotated, and characterized from the channel catfish genome. Extensive phylogenetic and comparative genomic analyses supported their annotations, allowing establishment of their orthologies, revealing fish-specific CC chemokines and the expansion of CC chemokines in the teleost genomes through lineage-specific tandem duplications. With 64 genes, the channel catfish genome harbors the largest numbers of CC chemokines among all the genomes characterized to date, however, they fall into 11 distinct CC chemokine groups. Analysis of gene expression after bacterial infections indicated that the CC chemokines were regulated in a gene-specific and time-dependent manner. While only one member of CCL19 (CCL19a.1) was significantly up-regulated after Edwardsiella ictaluri infection, all CCL19 members (CCL19a.1, CCL19a.2 and CCL19b) were significantly induced after Flavobacterium columnare infection. In addition, CCL19a.1, CCL19a.2 and CCL19b were also drastically up-regulated in ESC-susceptible fish, but not in resistant fish, suggesting potential significant roles of CCL19 in catfish immune responses. High expression levels of certain CC appeared to be correlated with susceptibility to diseases and intolerance to hypoxia.

Transcriptome comparison reveals insights into muscle response to hypoxia in blunt snout bream (Megalobrama amblycephala)

The economic and biological significance of blunt snout bream (Megalobrama amblycephala) makes this species important to explore the underlying molecular mechanism of hypoxia response. In the present study, we compared the transcriptional responses to serious hypoxia in skeletal muscle among hypoxia tolerant (MT), sensitive (MS) and control (without hypoxia treatment, MC) M. amblycephala obtained according to the time difference of losing balance after hypoxia treatment. A total of 88,200,889 clean reads were generated and assembled into 44,493 unigenes. Transcriptomic comparison revealed 463 genes differentially expressed among different groups. A similar hypoxia-induced transcription patterns suggested a common hypoxia response involved in cell cycle, p53 signaling pathway, apoptosis, heart contraction and blood circulation. Interesting, four genes, heat shock protein beta-8 (hspb8), cysteine/serine-rich nuclear protein 1 (csrnp1), salt-inducible kinase 1 (sik1), and visinin-like 1a (vsnl1a) were up-regulated in MT Vs MC but down-regulated in MS Vs MC. Additionally, FoxO signaling pathway was significantly enriched only in MT Vs MC. These results not only provided the first insights into the mechanism that muscle tissue coped with the hypoxia stress in cyprinid species, but offered a theory base for breeding of M. amblycephala with hypoxia-resistant traits.

Ras oncogene and Hypoxia-inducible factor-1 alpha (hif-1α) expression in the Amazon fish Colossoma macropomum (Cuvier, 1818) exposed to benzo[a]pyrene

Benzo[a]pyrene (B[a]P) is a petroleum derivative capable of inducing cancer in human and animals. In this work, under laboratory conditions, we analyzed the responses of Colossoma macropomum to B[a]P acute exposure through intraperitoneal injection of four different B[a]P concentrations (4, 8, 16 and 32 μmol/kg) or corn oil (control group). We analyzed expression of the ras oncogene and the Hypoxia-inducible factor-1 alpha (hif-1α) gene using quantitative real-time PCR. Additionally, liver histopathological changes and genotoxic effects were evaluated through the comet assay. Ras oncogene was overexpressed in fish exposed to 4, 8 of 16 μmol/kg B[a]P, showing 4.96, 7.10 and 6.78-fold increases, respectively. Overexpression also occurred in hif-1α in fish injected with 4 and 8 μmol/kg B[a]P, showing 8.82 and 4.64-fold increases, respectively. Histopathological damage in fish liver was classified as irreparable in fish exposed to 8, 16 and 32 μmol/kg μM B[a]P. The genotoxic damage increased in fish injected with 8 and 16 μmol/kg in comparison with the control group. Acute exposure of B[a]P was capable to interrupt the expression of ras oncogene and hif-1α, and increase DNA breaks due to tissue damage.

Identification of HIF-1 signaling pathway in Pelteobagrus vachelli using RNA-Seq: effects of acute hypoxia and reoxygenation on oxygen sensors, respiratory metabolism, and hematology indices

Oxygen is a vital element in aquatic environments. The concentration of oxygen to which aquatic organisms are exposed is influenced by salinity, water temperature, weather, and surface water runoff. Hypoxia has a serious effect on fish populations, and can lead to the loss of habitat and die-offs. Therefore, in the present study we used next-generation sequencing technology to characterize the transcriptomes of Pelteobagrus vachelli and identified 70 candidate genes in the HIF-1 signaling pathway that are important for the hypoxic response in all metazoan species. For the first time, the present study reported the effects of acute hypoxia and reoxygenation on oxygen sensors, respiratory metabolism, and hematology indices in P. vachelli. The predicted physiological adjustments show that P. vachelli's blood oxygen-carrying capacity was increased through increased RBC, HB, and SI after hypoxia exposure. Glycolysis-related enzyme activities (PFK, HK, and PK) and LDH in the brain and liver also increased, indicating a rise in anaerobic metabolism. The observed reduction in oxidative enzyme level (CS) in the liver during hypoxia suggests a concomitant depression in aerobic metabolism. There were significant increases in oxygen sensor mRNA expression and HIF-1α protein expression during hypoxia and reoxygenation exposure, suggesting that the HIF-1 signaling pathway was activated in the liver and brain of P. vachelli in response to acute hypoxia and reoxygenation. Our findings suggest that oxygen sensors (e.g., HIF-1α) of P. vachelli are potentially useful biomarkers of environmental hypoxic exposure. These data contribute to a better understanding of the molecular mechanisms of the hypoxia signaling pathway in fish under hypoxia and reoxygenation.

Characterization and functional analysis of hypoxia-inducible factor HIF1α and its inhibitor HIF1αn in tilapia

Hypoxia is a major cause of fish morbidity and mortality in the aquatic environment. Hypoxia-inducible factors are very important modulators in the transcriptional response to hypoxic stress. In this study, we characterized and conducted functional analysis of hypoxia-inducible factor HIF1α and its inhibitor HIF1αn in Nile tilapia (Oreochromis niloticus). By cloning and Sanger sequencing, we obtained the full length cDNA sequences for HIF1α (2686bp) and HIF1αn (1308bp), respectively. The CDS of HIF1α includes 15 exons encoding 768 amino acid residues and the CDS of HIF1αn contains 8 exons encoding 354 amino acid residues. The complete CDS sequences of HIF1α and HIF1αn cloned from tilapia shared very high homology with known genes from other fishes. HIF1α show differentiated expression in different tissues (brain, heart, gill, spleen, liver) and at different hypoxia exposure times (6h, 12h, 24h). HIF1αn expression level under hypoxia is generally increased (6h, 12h, 24h) and shows extremely highly upregulation in brain tissue under hypoxia. A functional determination site analysis in the protein sequences between fish and land animals identified 21 amino acid sites in HIF1α and 2 sites in HIF1αn as significantly associated sites (α = 0.05). Phylogenetic tree-based positive selection analysis suggested 22 sites in HIF1α as positively selected sites with a p-value of at least 95% for fish lineages compared to the land animals. Our study could be important for clarifying the mechanism of fish adaptation to aquatic hypoxia environment.

Differential expression of survival proteins during decreased intracellular oxygen tension in brain endothelial cells of grey mullets

The brain requires constant oxygen supply to perform its biological functions essential for survival. Because of low oxygen capacity and poor oxygen diffusibility of water, many fish species have evolved various adaptive mechanisms to cope with depleted oxygen. Endothelial cells (EC) are the primary components responsible for controlled environment of brain. Brain homeostasis largely depends on integrity of the EC. To elucidate their adaptive strategy, EC were isolated from the fish brain of Kovalam-control site and Ennore estuary-test/field hypoxic site and were subjected to low oxygen tension in laboratory. Cell viability, 4-hydroxynonenal (4HNE) and total antioxidant capacity (TAC) were analyzed to ascertain stress. Hypoxic insult, cytoprotective role of HSPs and apoptotic effect were analyzed by assessing hypoxia-inducible-factor-α (HIF1α), heat-shock-protein-70 (HSP70), heme-oxygenase 1 (HO-1), and apoptosis signal regulating kinase-1 (ASK1). This study evidenced that HSP70 and HO-1 are the key stress proteins, confer high tolerance to decreased oxygen tension mediated stress.

The chemokinome superfamily in channel catfish: I. CXC subfamily and their involvement in disease defense and hypoxia responses

Chemokines are a superfamily of structurally related chemotactic cytokines exerting significant roles in regulating cell migration and activation. They are defined by the presence of four conserved cysteine residues and are divided into four subfamilies depending on the arrangement of the first two conserved cysteines residues: CXC, CC, C and CX3C. In this study, a complete set of 17 CXC chemokine ligand (CXCL) genes was systematically identified and characterized from channel catfish genome through data mining of existing genomic resources. Phylogenetic analysis allowed annotation of the 17 CXC chemokines. Extensive comparative genomic analyses supported their annotations and orthologies, revealing the existence of fish-specific CXC chemokines and the expansion of CXC chemokines in the teleost genomes. The analysis of gene expression after bacterial infection indicated the CXC chemokines were expressed in a gene-specific manner. CXCL11.3 and CXCL20.3 were expressed significantly higher in resistant fish than in susceptible fish after ESC infection, while CXCL20.2 were expressed significantly higher in resistant fish than in susceptible fish after columnaris infection. The expression of those CXC chemokines, therefore can be a useful indicator of disease resistance. A similar pattern of expression was observed between resistant and susceptible fish with biotic and abiotic stresses, ESC, columnaris and hypoxia, suggesting that high levels of expression of the majority of CXC chemokines, with exception of CXC11 and CXC20, are detrimental to the host.

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.

Reactive oxygen species (ROS) and the heat stress response of Daphnia pulex: ROS-mediated activation of hypoxia-inducible factor 1 (HIF-1) and heat shock factor 1 (HSF-1) and the clustered expression of stress genes

BACKGROUND INFORMATION

Heat stress in ectotherms involves direct (e.g. protein damage) and/or indirect effects (temperature-induced hypoxia and ROS formation), which cause activation of the transcription factors (TF) heat shock factor 1 (HSF-1) and/or hypoxia-inducible factor 1 (HIF-1). The present study focused on the links between stress (ROS) signals, nuclear (n) and cytoplasmic (c) HSF-1/HIF-1 levels, and stress gene expression on mRNA and protein levels (e.g. heat-shock protein 90, HSP90) upon acute heat and ROS (H₂ O₂ ) stress.

RESULTS

Acute heat stress (30°C) evoked fluctuations in ROS level. Different feeding regimens, which affected the glutathione (GSH) level, allowed altering the frequency of ROS fluctuations. Other data showed fluctuation frequency to depend also on ROS production rate. The heat-induced slow or fast ROS fluctuations (at high or low GSH levels) evoked slow or fast fluctuations in the levels of nHIF-1α, nHSF-1 and gene products (mRNAs and protein), albeit after different time delays. Time delays to ROS fluctuations were, for example,shorter for nHIF-1α than for nHSF-1 fluctuations, and nHIF-1α fluctuations preceded and nHSF-1 fluctuations followed fluctuations in HSP90 mRNA level. Cytoplasmic TF levels either changed little (cHIF-1α) or showed a steady increase (cHSF-1). Applying acute H₂ O₂ stress (at 20°C) revealed effects on nHIF-1α and mRNA levels, but no significant effects on nHSF-1 level. Transcriptome data additionally showed coordinated fluctuations of mRNA levels upon acute heat stress, involving mRNAs for HSPs and other stress proteins, with all corresponding genes carrying DNA binding motifs for HIF-1 and HSF-1.

CONCLUSIONS

This study provided evidence for promoting effects of ROS and HIF-1 on early haemoglobin, HIF-1α and HSP90 mRNA expressions upon heat or ROS stress. The increasing cHSF-1 level likely affected nHSF-1 level and later HSP90 mRNA expression.

SIGNIFICANCE

Heat stress evoked ROS fluctuations, with this stress signal forwarded via nHIF-1 and nHSF-1 fluctuations to stress gene expression. The frequency of ROS fluctuations seemed to integrate information about ROS productionrate and GSH antioxidant buffer capacity, resulting in stress protein expression of different speed. Results of this study suggest ROS as early (pre-damage) and protein defects as later (post-damage) stress signals to trigger heat stress responses.

Comparative studies of hemolymph physiology response and HIF-1 expression in different strains of Litopenaeus vannamei under acute hypoxia

Litopenaeus vannamei has a high commercial value and is the primary cultured shellfish species globally. In this study, we have compared the hemolymph physiological responses between two L. vannamei strains under acute hypoxia. The results showed that hemocyanin concentration (HC) of strain A6410 was significantly higher than strain Zhengda; Total hemocyte counts (THC) decreased significantly in both strains under hypoxic stress (p < 0.05). We also investigated the temporal and spatial variations of hypoxia inducible factors 1 (HIF-1) by qRT-PCR. The results showed that hypoxia for 12 h increased the expression levels of HIF-1α in tissues of muscle and gill from the two strains (p < 0.05). In the hepatopancreas, the expression levels of HIF-1 increased significantly in strain Zhengda and decreased significantly in strain A6410 (p < 0.05). No significant changes of HIF-1 expression were detected in the same tissues between the two strains under hypoxia for 6 h (p > 0.05), but in the gills and hepatopancreas under hypoxia for 12 h (p < 0.05). Additionally, the expression level of HIF-1 was higher in the strain Zhengda than A6410 in the same tissue under hypoxia for 12 h. It was indicated that the hypoxic tolerance of Litopenaeus vannamei was closely correlated with the expression level of HIF-1, and the higher expression level of HIF-1 to hypoxia, the lower tolerance to hypoxia in the early stage of hypoxia. These results can help to better understand the molecular mechanisms of hypoxic tolerance and speed up the selective breeding process of hypoxia tolerance in L. vannamei.

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.