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

Transcriptomic atlas for hypoxia and following re-oxygenation in Ancherythroculter nigrocauda heart and brain tissues: insights into gene expression, alternative splicing, and signaling pathways

Hypoxia is a mounting problem that affects the world's freshwaters, with severe consequence for many species, including death and large economical loss. The hypoxia problem has increased recently due to the combined effects of water eutrophication and global warming. In this study, we investigated the transcriptome atlas for the bony fish Ancherythroculter nigrocauda under hypoxia for 1.5, 3, and 4.5 h and its recovery to normal oxygen levels in heart and brain tissues. We sequenced 21 samples for brain and heart tissues (a total of 42 samples) plus three control samples and obtained an average of 32.40 million raw reads per sample, and 95.24% mapping rate of the filtered clean reads. This robust transcriptome dataset facilitated the discovery of 52,428 new transcripts and 6,609 novel genes. In the heart tissue, the KEGG enrichment analysis showed that genes linked to the Vascular smooth muscle contraction and MAPK and VEGF signaling pathways were notably altered under hypoxia. Re-oxygenation introduced changes in genes associated with abiotic stimulus response and stress regulation. In the heart tissue, weighted gene co-expression network analysis pinpointed a module enriched in insulin receptor pathways that was correlated with hypoxia. Conversely, in the brain tissue, the response to hypoxia was characterized by alterations in the PPAR signaling pathway, and re-oxygenation influenced the mTOR and FoxO signaling pathways. Alternative splicing analysis identified an average of 27,226 and 28,290 events in the heart and brain tissues, respectively, with differential events between control and hypoxia-stressed groups. This study offers a holistic view of transcriptomic adaptations in A. nigrocauda heart and brain tissues under oxygen stress and emphasizes the role of gene expression and alternative splicing in the response mechanisms.

Effect of Agaricus bisporus Polysaccharides (ABPs) on anti-CCV immune response of channel catfish

Herein, the effects of Agaricus bisporus Polysaccharides (ABPs) on anti-channel catfish virus (CCV) infections to promote their application in channel catfish culture were explored. Transcriptome and metabolome analyses were conducted on the spleen of a CCV-infected channel catfish model fed with or without ABPs. CCV infections upregulated many immune and apoptosis-related genes, such as IL-6, IFN-α3, IFN-γ1, IL-26, Casp3, Casp8, and IL-10, and activated specific immunity mediated by B cells. However, after adding ABPs, the expression of inflammation-related genes decreased in CCV-infected channel catfish, and the inflammatory inhibitors NLRC3 were upregulated. Meanwhile, the expression of apoptosis-related genes was reduced, indicating that ABPs can more rapidly and strongly enhance the immunity of channel catfish to resist viral infection. Moreover, the metabonomic analysis showed that channel catfish had a high energy requirement during CCV infection, and ABPs could enhance the immune function of channel catfish. In conclusion, ABPs can enhance the antiviral ability of channel catfish by enhancing immune response and regulating inflammation. Thus, these findings provided new insights into the antiviral response effects of ABPs, which might support their application in aquaculture.

Living in a hypoxic world: A review of the impacts of hypoxia on aquaculture

Hypoxia is a harmful result of anthropogenic climate change. With the expansion of global low-oxygen zones (LOZs), many organisms have faced unprecedented challenges affecting their survival and reproduction. Extensive research has indicated that oxygen limitation has drastic effects on aquatic animals, including on their development, morphology, behavior, reproduction, and physiological metabolism. In this review, the global distribution and formation of LOZs were analyzed, and the impacts of hypoxia on aquatic animals and the molecular responses of aquatic animals to hypoxia were then summarized. The commonalities and specificities of the response to hypoxia in aquatic animals in different LOZs were discussed lastly. In general, this review will deepen the knowledge of the impacts of hypoxia on aquaculture and provide more information and research directions for the development of fishery resource protection strategies.

Effects of diets rich in Agaricus bisporus polysaccharides on the growth, antioxidant, immunity, and resistance to Yersinia ruckeri in channel catfish

To promote the application of Agaricus bisporus polysaccharides (ABPs) in channel catfish (Ictalurus punctatus) culture, we evaluated the effects of ABPs on the growth, immunity, antioxidant, and antibacterial activity of channel catfish. When the amount of ABPs was 250 mg/kg, channel catfish's weight gain and specific growth rates increased significantly while the feed coefficient decreased. We also found that adding ABPs in the feed effectively increased the activities of ACP, MDA, T-SOD, AKP, T-AOC, GSH, and CAT enzymes and immune-related genes such as IL-1β, Hsp70, and IgM in the head kidney of channel catfish. Besides, long-term addition will not cause pathological damage to the head kidney. When the amount of ABPs was over 125 mg/kg, the protection rate of channel catfish was more than 60%. According to the intestinal transcriptome analysis, the addition of ABPs promoted the expression of intestinal immunity genes and growth metabolism-related genes and enriched multiple related KEEG pathways. When challenged by Yersinia ruckeri infection, the immune response of channel catfish fed with ABPs was intenser and quicker. Additionally, the 16S rRNA gene sequencing analysis showed that the composition of the intestinal microbial community of channel catfish treated with ABPs significantly changed, and the abundance of microorganisms beneficial to channel catfish growth, such as Firmicutes and Bacteroidota increased. In conclusion, feeding channel catfish with ABPs promoted growth, enhanced immunity and antioxidant, and improved resistance to bacterial infections. Our current results might promote the use of ABPs in channel catfish and even other aquacultured fish species.

Transcriptome and 16S rRNA Analyses Reveal That Hypoxic Stress Affects the Antioxidant Capacity of Largemouth Bass (Micropterus salmoides), Resulting in Intestinal Tissue Damage and Structural Changes in Microflora

Dissolved oxygen (DO) is a key factor affecting the health of aquatic organisms in an intensive aquaculture environment. In this study, largemouth bass (Micropterus salmoides) were subjected to acute hypoxic stress for 96 h (DO: 1.00 mg/L) followed by recovery under sufficient DO conditions (DO: 7.50 mg/L) for 96 h. Serum biochemical indices, intestinal histomorphology, the transcriptome, and intestinal microbiota were compared between hypoxia-treated fish and those in a control group. The results showed that hypoxia caused oxidative stress, exfoliation of the intestinal villus epithelium and villus rupture, and increased cell apoptosis. Transcriptome analyses revealed that antioxidant-, inflammation-, and apoptosis-related pathways were activated, and that the MAPK signaling pathway played an important role under hypoxic stress. In addition, 16S rRNA sequencing analyses revealed that hypoxic stress significantly decreased bacterial richness and identified the dominant phyla (Proteobacteria, Firmicutes) and genera (Mycoplasma, unclassified Enterobacterales, Cetobacterium) involved in the intestinal inflammatory response of largemouth bass. Pearson's correlation analyses showed that differentially expressed genes in the MAPK signaling pathway were significantly correlated with some microflora. The results of this study will help to develop strategies to reduce damage caused by hypoxic stress in aquacultured fish.

Comparative Transcriptome Analysis of Head Kidney of Aeromonas hydrophila-infected Hypoxia-tolerant and Normal Large Yellow Croaker

The large yellow croaker (Larimichthys crocea) is one of the most economically important marine fish on the southeast coast of China and much of its yield is usually lost by hypoxia. To address this problem and lay a foundation for culturing a new strain of large yellow croaker with hypoxia tolerance, our research group screened a hypoxia-tolerant population of L. crocea. Surprisingly, we also found that hypoxia-tolerant population exhibited higher survival when infected with pathogens compared to the normal population during the farming operation. In order to understand the mechanism underlying the higher survival rate of the hypoxia-tolerant population and enrich the head kidney immune mechanism of L. crocea infected with pathogens, we compared and analyzed the head kidney transcriptome of the hypoxia-tolerant and normal individuals under Aeromonas hydrophila infection. We obtained 159.68 GB high-quality reads, of which more than 87.61% were successfully localized to the reference genome of L. crocea. KEGG analysis revealed differentially expressed genes in the signaling pathways involving immunity, cell growth and death, transport and catabolism, and metabolism. Among these, the toll-like receptor signaling pathway, Nod-like receptor signaling pathway, cytokine-cytokine receptor interaction, phagosome, apoptosis, and OXPHOS pathways were enriched in both groups after infection compared to before, and were enriched in infected tolerant individuals compared to normal individuals. In addition, we found that the expression of hif1α and its downstream genes were higher in the hypoxia-sensitive group of fish than in the normal group. In conclusion, our results showed some signaling pathways and hub genes, which may participate in A. hydrophila defense in the head kidney of two populations, and may contribute to the higher survival rate in the hypoxia-tolerant population. Overall, these findings increase our understanding of the defense mechanism within the head kidney of L. crocea under A. hydrophila infection, and suggest a preliminary hypothesis for why hypoxia-tolerant individuals may exhibit a higher survival rates after infection. Our study provides scientific evidence for the breeding of a new hypoxia-tolerant strain of L. crocea for aquaculture.

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.

Genome-wide identification and expression analysis of mitogen-activated protein kinase (MAPK) genes in response to salinity stress in channel catfish (Ictalurus punctatus)

The mitogen-activated protein kinase (MAPK) gene family has been systematically described in several fish species, but less so in channel catfish (Ictalurus punctatus), which is an important global aquaculture species. In this study, 16 MAPK genes were identified in the channel catfish genome and classified into three subfamilies based on phylogenetic analysis, including six extracellular signal regulated kinase (ERK) genes, six p38-MAPK genes and four C-Jun N-terminal kinase (JNK) genes. All MAPK genes were distributed unevenly across 10 chromosomes, of which three (IpMAPK8, IpMAPK12 and IpMAPK14) underwent teleost-specific whole genome duplication during evolution. Gene expression profiles in channel catfish during salinity stress were analysed using transcriptome sequencing and qRT-PCR (quantitative reverse transcription PCR). Results from reads per kilobase million (RPKM) analysis showed IpMAPK13, IpMAPK14a and IpMAPK14b genes were differentially expressed when compared with other genes between treatment and control groups. Furthermore, three of these genes were validated by qRT-PCR, of which IpMAPK14a expression levels were significantly upregulated in treatment groups (high and low salinity) when compared with the control group, with the highest expression levels in the low salinity group (P < 0.05). Therefore, IpMAPK14a may have important response roles to salinity stress in channel catfish.

Comparative genomics reveals putative evidence for high-elevation adaptation in the American pika (Ochotona princeps)

High-elevation environments have lower atmospheric oxygen content, reduced temperatures, and higher levels of UV radiation than found at lower elevations. As such, species living at high elevations must overcome these challenges to survive, grow, and reproduce. American pikas (Ochotona princeps) are alpine lagomorphs that are habitat specialists typically found at elevations >2,000 m. Previous research has shown putative evidence for high-elevation adaptation; however, investigations to date have been limited to a fraction of the genome. Here, we took a comparative genomics approach to identify putative regions under selection using a chromosomal reference genome assembly for the American pika relative to 8 other mammalian species targeted based on phylogenetic relatedness and (dis)similarity in ecology. We first identified orthologous gene groups across species and then extracted groups containing only American pika genes as well as unclustered pika genes to inform functional enrichment analyses; among these, we found 141 enriched terms with many related to hypoxia, metabolism, mitochondrial function/development, and DNA repair. We identified 15 significantly expanded gene families within the American pika across all orthologous gene groups that displayed functionally enriched terms associated with hypoxia adaptation. We further detected 196 positively selected genes, 41 of which have been associated with putative adaptation to hypoxia, cold tolerance, and response to UV following a literature review. In particular, OXNAD1, NRDC, and those genes critical in DNA repair represent important targets for future research to examine their functional implications in the American pika, especially as they may relate to adaptation to rapidly changing environments.

The transcriptomic responses of blunt snout bream (Megalobrama amblycephala) to acute hypoxia stress alone, and in combination with bortezomib

BACKGROUND

Blunt snout bream (Megalobrama amblycephala) is sensitive to hypoxia. A new blunt snout bream strain, "Pujiang No.2", was developed to overcome this shortcoming. As a proteasome inhibitor, bortezomib (PS-341) has been shown to affect the adaptation of cells to a hypoxic environment. In the present study, bortezomib was used to explore the hypoxia adaptation mechanism of "Pujiang No.2". We examined how acute hypoxia alone (hypoxia-treated, HN: 1.0 mg·L^- 1), and in combination with bortezomib (hypoxia-bortezomib-treated, HB: Use 1 mg bortezomib for 1 kg fish), impacted the hepatic ultrastructure and transcriptome expression compared to control fish (normoxia-treated, NN).

RESULTS

Hypoxia tolerance was significantly decreased in the bortezomib-treated group (LOE_crit, loss of equilibrium, 1.11 mg·L^- 1 and 1.32 mg·L^- 1) compared to the control group (LOE_crit, 0.73 mg·L^- 1 and 0.85 mg·L^- 1). The HB group had more severe liver injury than the HN group. Specifically, the activities of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) in the HB group (52.16 U/gprot, 32 U/gprot) were significantly (p < 0.01) higher than those in the HN group (32.85 U/gprot, 21. 68 U/gprot). In addition, more severe liver damage such as vacuoles, nuclear atrophy, and nuclear lysis were observed in the HB group. RNA-seq was performed on livers from the HN, HB and NN groups. KEGG pathway analysis disclosed that many DEGs (differently expressed genes) were enriched in the HIF-1, FOXO, MAPK, PI3K-Akt and AMPK signaling pathway and their downstream.

CONCLUSION

We explored the adaptation mechanism of "Pujiang No.2" to hypoxia stress by using bortezomib, and combined with transcriptome analysis, accurately captured the genes related to hypoxia tolerance advantage.

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.

Dual RNA-Seq of Trunk Kidneys Extracted From Channel Catfish Infected With Yersinia ruckeri Reveals Novel Insights Into Host-Pathogen Interactions

Host-pathogen intectarions are complex, involving large dynamic changes in gene expression through the process of infection. These interactions are essential for understanding anti-infective immunity as well as pathogenesis. In this study, the host-pathogen interaction was analyzed using a model of acute infection where channel catfish were infected with Yersinia ruckeri. The infected fish showed signs of body surface hyperemia as well as hyperemia and swelling in the trunk kidney. Double RNA sequencing was performed on trunk kidneys extracted from infected channel catfish and transcriptome data was compared with data from uninfected trunk kidneys. Results revealed that the host-pathogen interaction was dynamically regulated and that the host-pathogen transcriptome fluctuated during infection. More specifically, these data revealed that the expression levels of immune genes involved in Cytokine-cytokine receptor interactions, the NF-kappa B signaling pathway, the JAK-STAT signaling pathway, Toll-like receptor signaling and other immune-related pathways were significantly upregulated. Y. ruckeri mainly promote pathogenesis through the flagellum gene fliC in channel catfish. The weighted gene co-expression network analysis (WGCNA) R package was used to reveal that the infection of catfish is closely related to metabolic pathways. This study contributes to the understanding of the host-pathogen interaction between channel catfish and Y. ruckeri, more specifically how catfish respond to infection through a transcriptional perspective and how this infection leads to enteric red mouth disease (ERM) in these fish.

Genetics Responses to Hypoxia and Reoxygenation Stress in Larimichthys crocea Revealed via Transcriptome Analysis and Weighted Gene Co-Expression Network

Simple Summary

Hypoxia, which occurs frequently in aquaculture, can cause serious harm to all aspects of the growth, reproduction and metabolism of cultured fish. Due to the intolerance of Larimichthys crocea to hypoxia, Larimichthys crocea often floats head or even dies under hypoxic environment. However, the molecular mechanism of hypoxia tolerance in Larimichthys crocea has not been fully described. Therefore, the aim of this study was to explore the hub regulatory genes under hypoxic stress environment by transcriptome analysis of three key tissues (liver, blood and gill) in Larimichthys crocea. We identified a number of important genes that exercise different regulatory functions. Overall, this study will provide important clues to the molecular mechanisms of hypoxia tolerance in Larimichthys crocea.

Abstract

The large yellow croaker (Larimichthys crocea) is an important marine economic fish in China; however, its intolerance to hypoxia causes widespread mortality. To understand the molecular mechanisms underlying hypoxia tolerance in L. crocea, the transcriptome gene expression profiling of three different tissues (blood, gills, and liver) of L. crocea exposed to hypoxia and reoxygenation stress were performed. In parallel, the gene relationships were investigated based on weighted gene co-expression network analysis (WGCNA). Accordingly, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that several pathways (e.g., energy metabolism, signal transduction, oxygen transport, and osmotic regulation) may be involved in the response of L. crocea to hypoxia and reoxygenation stress. In addition, also, four key modules (darkorange, magenta, saddlebrown, and darkolivegreen) that were highly relevant to the samples were identified by WGCNA. Furthermore, some hub genes within the association module, including RPS16, EDRF1, KCNK5, SNAT2, PFKL, GSK-3β, and PIK3CD, were found. This is the first study to report the co-expression patterns of a gene network after hypoxia stress in marine fish. The results provide new clues for further research on the molecular mechanisms underlying hypoxia tolerance in L. crocea.

Selection and introgression facilitated the adaptation of Chinese native endangered cattle in extreme environments

Although persistent efforts have identified and characterized a few candidate genes and related biological processes with potential functions in the adaptation of many species to extreme environments, few works have been conducted to determine the genomic basis of adaptation in endangered livestock breeds that have been living in extreme conditions for more than thousands of years. To fill this gap, we sequenced the whole genomes of nine individuals from three Chinese native endangered cattle breeds that are living in high-altitude or arid environments. Phylogenetic and evolutionary history analyses of these three and other six breeds showed that the genetic structure of the cattle populations is primarily related to geographic location. Interestingly, we identified pervasive introgression from the yak to Zhangmu cattle (ZMC) that cover several genes (e.g., NOS2, EGLN1 and EPAS1) involved in the hypoxia response and previously identified as positive selection genes in other species, which suggested that the adaptive introgression from yak may have contributed to the adaptation of ZMC to high-altitude environments. In addition, by contrasting the breeds in opposite living conditions, we revealed a set of candidate genes with various functions from hypoxia response, water metabolism, immune response and body shape change to embryo development and skeletal system development, etc., that may be related to high-altitude or arid adaptation. Our research provides new insights into the recovery and adaptation of endangered native cattle and other species in extreme environments and valuable resources for future research on animal husbandry to cope with climate change.

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

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

NOD-like receptor signaling pathway activation: A potential mechanism underlying negative effects of benzo(α)pyrene on zebrafish

Benzo(α)pyrene (BaP) is one of typical polycyclic aromatic hydrocarbons (PAHs) in aquatic environments and has been shown to cause toxic effects to aquatic animals. Although the negative effects of BaP have been investigated, the potential toxic mechanisms remain uncharacterized. To explore the potential mechanisms mediating the toxic effects of BaP, zebrafish (Danio rerio) were exposed to BaP for 15 days and the toxic effects of BaP in zebrafish liver were investigated using physiological and transcriptomic analyses. After 15-day BaP exposure, zebrafish liver exhibited abnormalities including increased cytoplasmic vacuolation, inflammatory cell infiltration, swelled nuclei and irregular pigmentation. BaP exposure also induced oxidative stress to the liver of zebrafish. Transcriptomic profiles revealed 5129 differentially expressed genes (DEGs) after 15-days of BaP exposure, and the vast majority of DEGs were up-regulated under BaP treatment. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggest that genes related to immune response were significantly dysregulated. Furthermore, the nucleotide-binding, oligomerization domain (NOD)-like receptor signaling pathway was significantly enriched and most of the genes in this pathway exhibited enhanced expression after BaP exposure. These results partially explained the mechanisms underlying the toxic effects of BaP on zebrafish liver. In conclusion, BaP has the potential to induce physiological responses in zebrafish liver through altering associated genes.

Knockdown of α-enolase (Eno1) genes by RNAi does not increase the sensitivity of Propsilocerus akamusi to cadmium stress

α-enolase (Eno1) is a multifunctional enzyme which can as a stress protein under various environmental stresses. Recent researches also reported that Eno1 appears to have Cd²⁺ stress-related functions in cadmium tolerant plants. Our previous study inferred that the Eno1 gene might play an important role in the response of Propsilocerus akamusi to exogenous Cd²⁺. However, reports on the role of the Eno1 gene in coping with cadmium stress are still limited. In this study, we evaluated the roles of PaEno1 in the tolerance of P. akamusi to Cd²⁺ using RNAi technology and the response of recombinant proteins of PaEno1 in an E. coli expression system under Cd²⁺ stress. Our results showed that knockdown of PaEno1 did not increase but reduce the sensitivity of P. akamusi larvae to Cd²⁺ stress. However, bioassays showed the expression of recombinant PaEno1 protein in Rosetta cells enhanced the growth ability of E. coli under Cd²⁺ stress. These results suggested that overexpression of PaEno1 can significantly enhance the tolerance to heavy metal cadmium stresses in E. coli cells. However, knockdown of PaEno1 genes by RNAi does not increase the sensitivity of P. akamusi to cadmium stress.

Identification of potential markers for type 2 diabetes mellitus via bioinformatics analysis

Type 2 diabetes mellitus (T2DM) is a multifactorial and multigenetic disease, and its pathogenesis is complex and largely unknown. In the present study, microarray data (GSE201966) of β-cell enriched tissue obtained by laser capture microdissection were downloaded, including 10 control and 10 type 2 diabetic subjects. A comprehensive bioinformatics analysis of microarray data in the context of protein-protein interaction (PPI) networks was employed, combined with subcellular location information to mine the potential candidate genes for T2DM and provide further insight on the possible mechanisms involved. First, differential analysis screened 108 differentially expressed genes. Then, 83 candidate genes were identified in the layered network in the context of PPI via network analysis, which were either directly or indirectly linked to T2DM. Of those genes obtained through literature retrieval analysis, 27 of 83 were involved with the development of T2DM; however, the rest of the 56 genes need to be verified by experiments. The functional analysis of candidate genes involved in a number of biological activities, demonstrated that 46 upregulated candidate genes were involved in ‘inflammatory response’ and ‘lipid metabolic process’, and 37 downregulated candidate genes were involved in ‘positive regulation of cell death’ and ‘positive regulation of cell proliferation’. These candidate genes were also involved in different signaling pathways associated with ‘PI3K/Akt signaling pathway’, ‘Rap1 signaling pathway’, ‘Ras signaling pathway’ and ‘MAPK signaling pathway’, which are highly associated with the development of T2DM. Furthermore, a microRNA (miR)-target gene regulatory network and a transcription factor-target gene regulatory network were constructed based on miRNet and NetworkAnalyst databases, respectively. Notably, hsa-miR-192-5p, hsa-miR-124-5p and hsa-miR-335-5p appeared to be involved in T2DM by potentially regulating the expression of various candidate genes, including procollagen C-endopeptidase enhancer 2, connective tissue growth factor and family with sequence similarity 105, member A, protein phosphatase 1 regulatory inhibitor subunit 1 A and C-C motif chemokine receptor 4. Smad5 and Bcl6, as transcription factors, are regulated by ankyrin repeat domain 23 and transmembrane protein 37, respectively, which might also be used in the molecular diagnosis and targeted therapy of T2DM. Taken together, the results of the present study may offer insight for future genomic-based individualized treatment of T2DM and help determine the underlying molecular mechanisms that lead to T2DM.

De novo transcriptome assembly and sex-biased gene expression in the gonads of Amur catfish (Silurus asotus)

Amur catfish is extensively distributed and cultured in Asian countries. Despite of economic importance, the genomic information of this species remains limited. A reference transcriptome of Amur catfish was assembled and the sex-biased gene expression in the gonads was characterized using RNA-sequencing. The assembled transcriptome of Amur catfish consisted of 74,840 transcripts. The N50, mean length and max length of transcripts are 1970, 1235 and 16,748 bp. Putative sex-specific transcripts were identified and sex-specific expression of the representative genes was verified by RT-PCR. Differential expression analysis identified 5401 ovary-biased and 5618 testis-biased genes. The ovary-biased genes were mainly enriched in pathways such as RNA transport and ribosome biogenesis in eukaryotes. The testis-biased genes were enriched in calcium signaling and cytokine-cytokine receptor interaction, etc. Our data provide a valuable genomic resource for further investigating the genetic basis of sex determination, sex differentiation and sexual dimorphism of catfish.

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.