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

Hypothalamic transcriptome response to simulated diel earthen pond hypoxia cycles in channel catfish (Ictalurus punctatus)

Commercial culture of channel catfish (Ictalurus punctatus) occurs in earthen ponds that are characterized by diel swings in dissolved oxygen concentration that can fall to severe levels of hypoxia, which can suppress appetite and lead to suboptimal growth. Given the significance of the hypothalamus in regulating these processes in other fishes, an investigation into the hypothalamus transcriptome was conducted to identify specific genes and expression patterns responding to hypoxia. Channel catfish in normoxic water were compared with catfish subjected to 12 h of hypoxia (20% oxygen saturation; 1.8 mg O2/L; 27°C) followed by 12 h of recovery in normoxia to mimic 24 h in a catfish aquaculture pond. Fish were sampled at 0-, 6-, 12-, 18-, and 24-h timepoints, with the 6- and 12-h samplings occurring during hypoxia. A total of 190 genes were differentially expressed during the experiment, with most occurring during hypoxia and returning to baseline values within 6 h of normoxia. Differentially expressed genes were sorted by function into Gene Ontology biological processes and revealed that most were categorized as "response to hypoxia," "sprouting angiogenesis," and "cellular response to xenobiotic stimulus." The patterns of gene expression reported here suggest that transcriptome responses to hypoxia are broad and quickly reversibly with the onset of normoxia. Although no genes commonly reported to modulate appetite were found to be differentially expressed in this experiment, several candidates were identified for future studies investigating the interplay between hypoxia and appetite in channel catfish, including adm, igfbp1a, igfbp7, and stc2b.NEW & NOTEWORTHY Channel catfish are an economically important species that experience diel episodic periods of hypoxia that can reduce appetite. This is the first study to investigate their transcriptome from the hypothalamus in a simulated 24-h span in a commercial catfish pond, with 12 h of hypoxia and 12 h of normoxia. The research revealed functional groups of genes relating to hypoxia, angiogenesis, and glycolysis as well as individual target genes possibly involved in appetite regulation.

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.

Whole Transcriptome Analysis of Differentially Expressed Genes in Cultured Nile Tilapia (O. niloticus) Subjected to Chronic Stress Reveals Signaling Pathways Associated with Depressed Growth

Chronic stress is a serious threat to aquaculture as it lowers fish growth performance and compromises fish welfare. The exact mechanism by which growth is retarded is, however, not clearly understood. This study sought to elucidate the gene expression profiles associated with chronic stress in cultured Nile tilapia (Oreochromis niloticus) reared for 70 days at different ammonia concentrations and stocking densities. Fish in the treatment groups showed negative growth, while the controls showed positive allometric growth. The specific condition factor (Kn) ranged from 1.17 for the controls to 0.93 for the ammonia and 0.91 for the stocking density treatments. RNA was extracted from muscle tissue using TRIzol followed by library construction and Illumina sequencing. Comparative transcriptome analysis revealed 209 differentially expressed genes (DEGs) (156 up- and 53 down-regulated) in the ammonia and 252 DEGs (175 up- and 77 down-regulated) in the stocking density treatment. In both treatments, 24 and 17 common DEGs were up- and down-regulated, respectively. DEGs were significantly enriched in six pathways associated with muscle activity, energy mobilization and immunity. The heightened muscular activity consumes energy which would otherwise have been utilized for growth. These results bring to fore the molecular mechanisms underlying chronic stress’ suppression of growth in cultured Nile tilapia.

Fipronil impairs the GABAergic brain responses of Nile Tilapia during the transition from normoxia to acute hypoxia

γ-aminobutyric acid (GABA) is one of the main neurotransmitters involved in the adaptation processes against the damage that hypoxia can cause to the brain. Due to its antagonist action on GABA receptors, the insecticide fipronil can turn the fish more susceptible to the negative effects of hypoxia. This study aimed to understand better if fipronil affects these GABAergic responses of Tilapia ahead to hypoxia. Oreochromis Niloticus (Nile Tilapia) were exposed for 3 and 8 h to fipronil (0.0, 0.1, and 0.5 µg.L^-1 ) under normoxia (dissolved O₂  > 6 mg.L^-1 ) and moderate hypoxia (dissolved O₂  < 2 mg.L^-1 ) conditions. Briefly, hypoxia caused opposite effects on the gene transcription of the evaluated ionotropic and metabotropic GABA receptors. Unexpectedly, we obtained reduced HIF1A mRNA and brain GABA levels, mostly in the first 3 h of the experiment, for the hypoxic group compared with the normoxia one. Besides that, we also demonstrated that the insecticide fipronil impairs the brain GABAergic signaling of a hypoxia-tolerant fish during the transition from a normoxic to an acute hypoxic state. Thus, these results predict the relevant impact on the brain metabolic adaptations of fishes exposed to such stressful conditions in an aquatic environment, as well as the effects of fipronil in the GABAergic responses to hypoxia, which in turn may have ecological and physiological significance to hypoxia-tolerant fishes exposed to this insecticide.

The Asparagine Hydroxylase FIH: A Unique Oxygen Sensor

Significance: Limited oxygen availability (hypoxia) commonly occurs in a range of physiological and pathophysiological conditions, including embryonic development, physical exercise, inflammation, and ischemia. It is thus vital for cells and tissues to monitor their local oxygen availability to be able to adjust in case the oxygen supply is decreased. The cellular oxygen sensor factor inhibiting hypoxia-inducible factor (FIH) is the only known asparagine hydroxylase with hypoxia sensitivity. FIH uniquely combines oxygen and peroxide sensitivity, serving as an oxygen and oxidant sensor. Recent Advances: FIH was first discovered in the hypoxia-inducible factor (HIF) pathway as a modulator of HIF transactivation activity. Several other FIH substrates have now been identified outside the HIF pathway. Moreover, FIH enzymatic activity is highly promiscuous and not limited to asparagine hydroxylation. This includes the FIH-mediated catalysis of an oxygen-dependent stable (likely covalent) bond formation between FIH and selected substrate proteins (called oxomers [oxygen-dependent stable protein oligomers]). Critical Issues: The (patho-)physiological function of FIH is only beginning to be understood and appears to be complex. Selective pharmacologic inhibition of FIH over other oxygen sensors is possible, opening new avenues for therapeutic targeting of hypoxia-associated diseases, increasing the interest in its (patho-)physiological relevance. Future Directions: The contribution of FIH enzymatic activity to disease development and progression should be analyzed in more detail, including the assessment of underlying molecular mechanisms and relevant FIH substrate proteins. Also, the molecular mechanism(s) involved in the physiological functions of FIH remain(s) to be determined. Furthermore, the therapeutic potential of recently developed FIH-selective pharmacologic inhibitors will need detailed assessment. Antioxid. Redox Signal. 37, 913-935.

Influence of hypoxia on biochemical aspects and on expression of genes related to oxygen-homeostasis of the Amazonian cichlid Astronotus ocellatus (Agassiz, 1831)

Variations in dissolved oxygen levels are common in the Amazonian aquatic environments and the aquatic organisms that inhabit these environments developed a variety of adaptive responses to deal with such conditions. Some Amazonian fish species are tolerant to low oxygen levels and the cichlid Astronotus ocellatus is one of the most hypoxia-tolerant species. Herein, we aimed to unveil the biochemical and molecular responses that A. ocellatus presents when submitted to hypoxia. Hypoxia indicators were measured, such as plasma glucose, plasma lactate, hepatic glycogen and relative transcript levels of prolyl hydroxylase 2 (phd2) and hypoxia-inducible factor-1α (hif-1α) in juveniles of approximately 50 g exposed to 1, 3, and 5 hours of hypoxia (0.7 mg O₂.L^-1), followed by 3 hours of recovery in normoxia (6 mg O₂.L^-1). Fish exposed to hypoxia reduced liver glycogen levels within 3 hours of hypoxia, when comparing with 1 hour, and increased plasma glucose and lactate. Under the same condition, phd2 transcripts levels increased in gills, but decreased in liver. In contrast, hypoxia did not affect relative gene expression of hif-1α in both tissues. Based on the transcription pattern of phd2, these results showed that liver and gills of A. ocellatus have different molecular strategies to cope with environmental hypoxia.

Biomarkers of Exposure and Effect in Nile Tilapia (Oreochromis niloticus) Environmentally Exposed to Multiple Stressors in Egypt

Chemical pollution poses adverse effects on aquatic organisms, including altered gene expression. This study, therefore, investigated the hepatic expression of biomarkers of exposure and effect in Oreochromis niloticus inhabiting the drainage canal of a sugarbeet factory (an industrially polluted habitat) and the Nile River (a reference habitat). Compared to the Nile River, the drainage canal contained significantly elevated levels of heavy metals (Cd, Fe, Mn, Co, and Pb), nitrate, ammonia, nitrite, chemical oxygen demand, and biological oxygen demand. The liver of O. niloticus from the drainage canal accumulated significantly higher concentrations of the heavy metals. O. niloticus from the drainage canal had significantly higher hepatic expressions of genes related to oxidative stress (GSTa, GPx, and NRF2), metal toxicity (MT), endocrine disruption (Vtg), and hypoxia (HIF1a) and a significantly lower hepatic expression of the apoptosis-associated CAS gene. CYP1A (a detoxification biomarker) expression did not exhibit significant differences. These results provide insights into complex adaptive responses of O. niloticus to multiple chemicals and support the utilization of molecular biomarkers as warning signals for water pollution.

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.

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.

Replacement of dietary fish oil with plant oils improves the immunological responses and the antioxidant status in Oreochromis niloticus exposed to suboptimal temperature

Here, we investigated the effects of total dietary fish oil (FO) substitution with plant oil (PO) on hematological indices, immune status, antioxidant activity, IL1β and TNF-α gene expression, and hypoxia stress resistance in Oreochromis niloticus at suboptimal temperatures. Fish (n = 360) were randomly divided into 12 circular fiberglass tanks (500 L; 3 replicates for each dietary group, 30 fish/replicate, 90 fish/group). The control group was fed a basal diet with FO as the lipid source. The CO, SFO, and LSO groups were fed a basal diet with complete replacement of FO with corn, sunflower, and linseed oils, respectively. After 4 weeks, no effects on hematological indices were observed in fish fed the experimental diets. The LSO and CO groups showed a significant increase in γ globulin levels. The highest levels of non-specific immune parameters, antioxidant activity, and IL-1β and TNF-α gene expressions were recorded in the LSO group. The LSO group also exhibited good resistance to hypoxia stress. Therefore, the total dietary substitution of FO with PO (especially LSO and CO) is recommended as a valuable strategy to ameliorate the immunosuppressive effects of suboptimal temperatures and enhance the resistance of O. niloticus to hypoxia stress.

Low oxygen: A (tough) way of life for Okavango fishes

Botswana's Okavango Delta is a World Heritage Site and biodiverse wilderness. In 2016-2018, following arrival of the annual flood of rainwater from Angola's highlands, and using continuous oxygen logging, we documented profound aquatic hypoxia that persisted for 3.5 to 5 months in the river channel. Within these periods, dissolved oxygen rarely exceeded 3 mg/L and dropped below 0.5 mg/L for up to two weeks at a time. Although these dissolved oxygen levels are low enough to qualify parts of the Delta as a dead zone, the region is a biodiversity hotspot, raising the question of how fish survive. In association with the hypoxia, histological samples, collected from native Oreochromis andersonii (threespot tilapia), Coptodon rendalli (redbreast tilapia), and Oreochromis macrochir (greenhead tilapia), exhibited widespread hepatic and splenic inflammation with marked granulocyte infiltration, melanomacrophage aggregates, and ceroid and hemosiderin accumulations. It is likely that direct tissue hypoxia and polycythemia-related iron deposition caused this pathology. We propose that Okavango cichlids respond to extended natural hypoxia by increasing erythrocyte production, but with significant health costs. Our findings highlight seasonal hypoxia as an important recurring stressor, which may limit fishery resilience in the Okavango as concurrent human impacts rise. Moreover, they illustrate how fish might respond to hypoxia elsewhere in the world, where dead zones are becoming more common.

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

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

Baicalin relieves hypoxia-aroused H9c2 cell apoptosis by activating Nrf2/HO-1-mediated HIF1α/BNIP3 pathway

Background: Myocardial ischemia is the main reason for ischemic heart disease. Baicalin is a plant-derived flavonoid with cardio-protective activity. Herein, we tested the influences of baicalin on cardiomyocytes H9c2 apoptosis aroused by hypoxia stimulation. Methods: Firstly, H9c2 cells were subjected to hypoxia and/or baicalin exposure. Cell viability and apoptosis, along with hypoxia-inducible factor 1α (HIF1α) and Bcl-2/adenovirus E1B 19-KDa interacting protein 3 (BNIP3) expressions were tested respectively. Then, si-HIF1α was transfected into H9c2 cells to probe whether up-regulation of HIF1α attended to the influences of baicalin on hypoxia-stimulated H9c2 cells. Finally, the regulatory effect of nuclear factor E2-related factor 2 (Nrf2)/heme oxygenase 1 (HO-1) pathway on HIF1α expression was analyzed. Results: Hypoxia exposure aroused H9c2 cell viability reduction and apoptosis. Baicalin mitigated H9c2 cell viability reduction and apoptosis aroused by hypoxia. Moreover, HIF1α/BNIP3 pathway was further activated by baicalin in hypoxia-exposed H9c2 cells. Silencing HIF1α lowered the functions of baicalin on hypoxia-exposed H9c2 cells. Besides, baicalin enhanced hypoxia-caused activation of Nrf2/HO-1 pathway. Activation of Nrf2/HO-1 pathway was associated with the up-regulation of HIF1α and protective functions of baicalin on hypoxia-exposed H9c2 cells. Conclusion: Baicalin relieved cardiomyocytes H9c2 apoptosis aroused by hypoxia might be achieved through activating Nrf2/HO-1-mediated HIF1α/BNIP3 pathway. Highlights Baicalin mitigates H9c2 cell viability loss and apoptosis aroused by hypoxia; Baicalin activates HIF1a/BNIP3 pathway in hypoxia-exposed H9c2 cells; Silencing HIF1α weakens the influences of baicalin on hypoxia-exposed H9c2 cells; Baicalin promotes Nrf2/HO-1 pathway in hypoxia-exposed H9c2 cells; Promotion of Nrf2/HO-1 pathway is related to the up-regulation of HIF1α.

Transcriptome signatures associated with meningioma progression

Meningiomas are the most common primary brain tumor of adults. The majority are benign (WHO grade I), with a mostly indolent course; 20% of them (WHO grade II and III) are, however, considered aggressive and require a more complex management. WHO grade II and III tumors are heterogeneous and, in some cases, can develop from a prior lower grade meningioma, although most arise de novo. Mechanisms leading to progression or implicated in de novo grade II and III tumorigenesis are poorly understood. RNA-seq was used to profile the transcriptome of grade I, II, and III meningiomas and to identify genes that may be involved in progression. Bioinformatic analyses showed that grade I meningiomas that progress to a higher grade are molecularly different from those that do not. As such, we identify GREM2, a regulator of the BMP pathway, and the snoRNAs SNORA46 and SNORA48, as being significantly reduced in meningioma progression. Additionally, our study has identified several novel fusion transcripts that are differentially present in meningiomas, with grade I tumors that did not progress presenting more fusion transcripts than all other tumors. Interestingly, our study also points to a difference in the tumor immune microenvironment that correlates with histopathological grade.

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.

Acting locally - affecting globally: RNA sequencing of gilthead sea bream with a mild Sparicotyle chrysophrii infection reveals effects on apoptosis, immune and hypoxia related genes

BACKGROUND

Monogenean flatworms are the main fish ectoparasites inflicting serious economic losses in aquaculture. The polyopisthocotylean Sparicotyle chrysophrii parasitizes the gills of gilthead sea bream (GSB, Sparus aurata) causing anaemia, lamellae fusion and sloughing of epithelial cells, with the consequent hypoxia, emaciation, lethargy and mortality. Currently no preventive or curative measures against this disease exist and therefore information on the host-parasite interaction is crucial to find mitigation solutions for sparicotylosis. The knowledge about gene regulation in monogenean-host models mostly comes from freshwater monopysthocotyleans and almost nothing is known about polyopisthocotyleans. The current study aims to decipher the host response at local (gills) and systemic (spleen, liver) levels in farmed GSB with a mild natural S. chrysophrii infection by transcriptomic analysis.

RESULTS

Using Illumina RNA sequencing and transcriptomic analysis, a total of 2581 differentially expressed transcripts were identified in infected fish when compared to uninfected controls. Gill tissues in contact with the parasite (P gills) displayed regulation of fewer genes (700) than gill portions not in contact with the parasite (NP gills) (1235), most likely due to a local silencing effect of the parasite. The systemic reaction in the spleen was much higher than that at the parasite attachment site (local) (1240), and higher than in liver (334). NP gills displayed a strong enrichment of genes mainly related to immune response and apoptosis. Processes such as apoptosis, inflammation and cell proliferation dominated gills, whereas inhibition of apoptosis, autophagy, platelet activation, signalling and aggregation, and inflammasome were observed in spleen. Proteasome markers were increased in all tissues, whereas hypoxia-related genes were down-regulated in gills and spleen.

CONCLUSIONS

Contrasting forces seem to be acting at local and systemic levels. The splenic down-regulation could be part of a hypometabolic response, to counteract the hypoxia induced by the parasite damage to the gills and to concentrate the energy on defence and repair responses. Alternatively, it can be also interpreted as the often observed action of helminths to modify host immunity in its own interest. These results provide the first toolkit for future studies towards understanding and management of this parasitosis.

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.

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.