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

Genome-wide identification and expression profiling of the hemoglobin gene family under hypoxia stress in Sillago sihama

Hemoglobin (Hb) plays a pivotal role in oxygen transport and is essential for the adaptive response to hypoxic stress. Studying the molecular evolution characteristics and expression pattern of S. sihama Hb gene family can provide theoretical basis for further revealing the roles of the Hb genes against hypoxia stress. This study aimed to identify the molecular evolution characteristics and expression patterns of Hb genes in S. sihama. Utilizing bioinformatics methods, we identified ten members of the Hb gene family within the whole genome of S. sihama, located on chromosomes LG16 and LG21 in a tandem repeat pattern. Phylogenetic analysis classified the Hb gene family into two subfamilies (alpha and beta), showing close relationships with the Hb gene family members of the large yellow croaker (Larimichthys crocea). Both phylogenetic and collinearity analysis indicated the conservation of Hb genes throughout evolution. The nonsynonymous substitution rate / synonymous substitution rate (Ka/Ks) results showed that its evolution was purifying selection. RNA-seq and qRT-PCR analyses revealed that five Hb genes were highlighted high expression levels in the heart and liver of S. sihama, with sshba1.4, sshbad, and sshba1 showing significant up-regulation in response to hypoxic stress. In summary, this study highlights the vital role of the Hb gene family in the hypoxic stress response of S. sihama, suggesting their potential as crucial candidate genes for regulating environmental stress in this species.

Transcriptomic analysis of immune-related genes in Pacific white snook (Centropomus viridis) gills infected with the monogenean parasite Rhabdosynochus viridisi

The parasite Rhabdosynochus viridisi (Platyhelminthes: Monogenea) infects the Pacific white snook Centropomus viridis gills and can cause adverse effects in the aquaculture industry. The immune responses of Pacific white snook to monogenean infections are poorly understood. Thus, this study aimed to identify differentially expressed genes (DEGs) in the gills of Pacific white snook juveniles experimentally infected with R. viridisi, emphasizing immune-related genes and pathways activated or suppressed during the infection. RNA sequencing was performed on the gills of uninfected (control) and infected fish. The algorithm Seq2Fun was selected without a reference transcriptome to map the reads to transcripts of fishes available from a database for gene orthologs (EcoOmics) and obtain the counting table. The ExpressAnalyst software was used for differential expression and functional analyses. A total of 20,106 transcripts were found, and 1430 (7 %) were differentially expressed genes (DEGs) between infected and control groups. We identified 860 (60 %) downregulated and 570 (40 %) upregulated genes. Thirteen canonical pathways after the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were overrepresented, and most of the DEGs were downregulated, suggesting the inactivation of these pathways. The functions of most of the DEGs with higher fold change found in this study are poorly understood in fish. Even though the well-known pro-inflammatory cytokines remained unchanged in infected gills of C. viridis, and transforming growth factor β (tgfβ) was downregulated, interleukin-17 ligands il17d and il17a/f1, as well as C-X-C motif chemokine receptor 2 (cxcr2) genes were upregulated, indicating that the infection with R. viridisi promotes Th17-like immunity. Overexpression of plasma B cell activity markers such as immunoglobulin light chain-like genes and the v-set pre-B cell surrogate light chain 3 (vpreb3) was also detected in this study. The possible implications of DEGs related to calcium imbalance, hypoxia adaptation, hemostasis, and immunity are discussed. These results will support future studies to improve the prevention and treatment of monogenean infections in finfish aquaculture.

Multi-omics analysis of the toxic effects on gill tissues of crucian carp (Carassius auratus) from chronic exposure to environmentally relevant concentrations of Di(2-ethylhexyl) phthalate (DEHP)

The pervasive use of the plasticizer di(2-ethylhexyl) phthalate (DEHP) poses potential risks to global aquatic ecosystems. This study systematically evaluated the adverse effects of chronic exposure to environmentally relevant concentrations of DEHP on gill tissues of crucian carp, utilizing histological examination, metabolomic, and transcriptomic analysis. The results demonstrated that DEHP induced significant histopathological alterations in gill tissues, with significant enrichment observed in multiple pathways associated with amino acid, hormone, lipid, and xenobiotic metabolism. Metabonomics-transcriptomics analyses indicated that DEHP-induced significantly over-activation of cytochrome P450 1B1-like (p < 0.001) and cytochrome P450 3A30-like (p < 0.05) via the nuclear xenobiotic receptors pathway was a key factor contributing to the disruption of tryptophan metabolism and steroid hormone biosynthesis, as well as inducing circadian rhythm disruption. Moreover, circadian rhythm disruption further exacerbated the imbalance of cytochrome P450 (CYP450) enzyme system as well as linoleic acid, arachidonic acid, sphingolipid, and glycerophospholipid metabolism. Overall, the feedback regulation between the CYP450 enzyme system and circadian rhythms emerged as the primary mechanism underlying DEHP-induced metabolic and transcriptional disruptions, ultimately resulting in gill toxicity. This study not only enriched the toxic effects on aquatic organisms of chronic exposure to DEHP, but provided potential biomarkers for the environmental risk assessment of DEHP.

Comparative analysis of intestinal structure, enzyme activity, intestinal microbiota and gene expression in different segments of pufferfish (Takifugu Obscurus)

The structure of fish intestines does not have a clear regional division, while the function of the intestines may be related to their structure. Therefore, in this study, the delimitation of intestinal segments in pufferfish (Takifugu obscurus) was achieved by morphological analysis. Subsequently, enzyme activity, intestinal microbiota, and gene expression were examined to compare the differences among the pufferfish various segments. According to four morphological parameters: height of mucosa folds (HF), width of mucosa folds (WF), thickness of muscularis (TM), and cross-sectional area (CSA), the pufferfish's intestine was divided into anterior intestine (AI), middle intestine (MI), and posterior intestine (PI). The activity levels of amylase, lipase, and trypsin in the AI and MI were significantly higher than these in the PI. According to the analysis of 16S rDNA, the dominant microbiota at the phylum level in the different segments were Epsilonbacteraeota, Spirochaetes, and Proteobacteria. At the genus level, there were variations observed in the relative abundance of Brevinema, Mycobacterium, Bradyrhizobium, and Microvirga. α diversity analysis revealed that the richness indexes (Ace and Chao1) were the lowest in the MI, while β diversity analysis revealed significant difference in intestinal microbial community composition among the three intestinal segments. Furthermore, RNA-Seq was used to identify differential expression genes (DEGs) and biological pathways among the different intestinal segments. The DEGs between the AI and MI were enriched in pancreatic secretion and protein digestion and absorption, those between AI and PI were involved in ascorbate and aldarate metabolism and glutathione metabolism, and those between MI and PI were involved in steroid biosynthesis, fat digestion and absorption, vitamin digestion and absorption, and glycine, serine and threonine metabolism. In conclusion, the presented results compare and analyze the differences in various intestinal segments of pufferfish, which will be conductive to future exploration of the functions of these different segments.

Antioxidant Capacity, Enzyme Activities Related to Energy Metabolism, and Transcriptome Analysis of Crassostrea hongkongensis Exposed to Hypoxia

Crassostrea hongkongensis (C. hongkongensis) is one of the three most commonly cultivated oyster species in China. Seasonal hypoxia is one of the most serious threats to its metabolism, reproductive behavior, and survival. To investigate the effects of hypoxia stress on the antioxidant capacity and energy metabolism of C. hongkongensis, the total antioxidant capacity (T-AOC), glycogen content, and enzyme activities (phosphofructokinase, PFK; pyruvate kinase, PK; phosphoenolpyruvate carboxykinase, PEPCK) of oysters were determined under normoxic (DO 6 ± 0.2 mg/L) and hypoxic (DO 1.5 mg/L) conditions at 0 h, 6 h, 48 h, and 72 h. We also determined the T-AOC, glycogen content, and enzyme activities of oysters under reoxygenation (recovered to normoxia for 24 h). To further examine the potential molecular regulatory mechanism of hypoxic adaptation, a transcriptome analysis was conducted on the gill of C. hongkongensis under normoxia (N, 72 h), hypoxia (H, 72 h), and reoxygenation (R). After being exposed to hypoxia for 6 h, the T-AOC, glycogen content, and enzyme activities of PK, PFK, and PEPCK in C. hongkongensis were significantly decreased. However, after prolonging the duration of hypoxia exposure for 72 h, the T-AOC, glycogen content, and enzyme activities increased compared to that of 48 h. After 24 h reoxygenation, the T-AOC, glycogen content, and enzyme activity of PK and PFK returned to close to initial levels. In addition, a transcriptome analysis discovered 6097 novel genes by mapping the C. hongkongensis genome with the clean reads. In total, 352 differentially expressed genes (DEGs) were identified in the H vs. N comparison group (235 upregulated and 117 downregulated genes). After recovery to normoxia, 292 DEGs (134 upregulated and 158 downregulated genes) were identified in the R vs. N comparison group, and 632 DEGs were identified (253 upregulated and 379 downregulated genes) in the R vs. H comparison group. The DEGs included some hypoxia-tolerant genes, such as phosphoenolpyruvate carboxykinase (PEPCK), mitochondrial (AOX), tyramine beta-hydroxylase (TBH), superoxide dismutase (SOD), glutathione S-transferase (GST), and egl nine homolog 1 isoform X2 (EGLN1). Additionally, DEGs were significantly enriched in the KEGG pathways that are involved in hypoxia tolerance, including the metabolism of xenobiotics by cytochrome P450 pathways and the HIF-1 signaling pathway. Then, we selected the five hypoxic-tolerant candidate DEGs for real-time quantitative polymerase chain reaction (RT-qPCR) validation, and the results were consistent with the transcriptome sequencing data. These discoveries have increased our understanding of hypoxia tolerance, recovery ability after reoxygenation, and molecular mechanisms governing the responses to hypoxia in C. hongkongensis.

Integrated analysis of miRNA profiles and gut bacterial changes in Altica viridicyanea following antibiotic treatment

The gut bacteria involves in insect homeostasis by playing essential roles in host physiology, metabolism, innate immunity, and so forth. microRNAs (miRNAs) are endogenous small noncoding RNAs that posttranscriptionally regulate gene expression to affect immune or metabolic processes in insects. For several non-model insects, the available knowledge on the relationship between changes in the gut bacteria and miRNA profiles is limited. In this study, we investigated the gut bacterial diversity, composition, and function from Altica viridicyanea feeding on normal- and antibiotic-treated host plants using 16S rRNA amplicon sequencing; antibiotics have been shown to affect the body weight and development time in A. viridicyanea, suggesting that the gut bacteria of the normal sample were more diverse and abundant than those of the antibiotic-fed group, and most of them were involved in various physical functions by enrichment analysis. Furthermore, we executed small RNA transcriptome sequencing using the two experimental groups to obtain numerous sRNAs, such as piRNAs, siRNAs, and known and novel miRNAs, by data mapping and quality control, and furthermore, a total of 224 miRNAs were identified as significantly differentially expressed miRNAs, of which some DEMs and their target genes participated in immune- and metabolism-related pathways based on GO and KEGG annotation. Besides, regarding the regulatory roles of miRNA and target genes, a interaction network of DEM-target gene pairs from eight immune- or metabolism-related signaling pathways were constructed. Finally, we discovered that DEMs from above pathways were significantly positively or negatively correlated with gut bacterial alterations following antibiotic treatment. Collectively, the observations of this study expand our understanding of how the disturbance of gut bacteria affects miRNA profiles in A. viridicyanea and provide new valuable resources from extreme ranges for future studies on the adaptive evolution in insects.

Hypoxia-induced physiological responses in fish: From organism to tissue to molecular levels

Dissolved oxygen (DO) in water bodies is a prerequisite for fish survival and plays a crucial role in fish growth, development, and physiological processes. However, with increasing eutrophication, greenhouse effects, and extreme weather conditions, DO levels in aquatic environments often become lower than normal. This leads to stress in fish, causing them to exhibit escape behavior, inhibits their growth and development, and causes tissue damage. Moreover, oxidative stress, decreased immune function, and altered metabolism have been observed. Severe hypoxia can cause massive fish mortality, resulting in significant economic losses to the aquaculture industry. In response to hypoxia, fish exhibit a series of behavioral and physiological changes that are self-protective mechanisms formed through long-term evolution. This review summarizes the effects of hypoxic stress on fish, including the asphyxiation point, behavior, growth and reproduction, tissue structure, physiological and biochemical processes, and regulation of gene expression. Furthermore, future research directions are discussed to provide new insights and references.

Gut flora alterations among aquatic firefly Aquatica leii inhabiting various dissolved oxygen in fresh water

Knowledge about the impact of different dissolved oxygen (DO) on the composition and function of gut bacteria of aquatic insects is largely unknown. Herein, we constructed freshwater environments with different DOs (hypoxia: 2.50 ± 0.50, normoxia: 7.00 ± 0.50, and hyperoxia: 13.00 ± 0.50 mg/L) where aquatic firefly Aquatica leii larvae lived for three months. Their gut flora was analyzed using the combination of 16S rRNA amplicon sequencing and metagenomics. The results showed no difference in alpha diversity of the gut flora between A. leii inhabiting various DOs. However, the relative abundance of several bacterial lineages presented significant changes, such as Pseudomonas. In addition, bacterial genes with an altered relative abundance in response to various DOs were primarily related to metabolism. The alteration of these functions correlated with the DO change. This is the first to uncover structure of gut flora under various DOs in aquatic insect larvae.

High temperature influences DNA methylation and transcriptional profiles in sea urchins (Strongylocentrotus intermedius)

BACKGROUND

DNA methylation plays an important role in life processes by affecting gene expression, but it is still unclear how DNA methylation is controlled and how it regulates gene transcription under high temperature stress conditions in Strongylocentrotus intermedius. The potential link between DNA methylation variation and gene expression changes in response to heat stress in S. intermedius was investigated by MethylRAD-seq and RNA-seq analysis. We screened DNA methylation driver genes in order to comprehensively elucidate the regulatory mechanism of its high temperature adaptation at the DNA/RNA level.

RESULTS

The results revealed that high temperature stress significantly affected not only the DNA methylation and transcriptome levels of S. intermedius (P < 0.05), but also growth. MethylRAD-seq analysis revealed 12,129 CG differential methylation sites and 966 CWG differential methylation sites, and identified a total of 189 differentially CG methylated genes and 148 differentially CWG methylated genes. Based on KEGG enrichment analysis, differentially expressed genes (DEGs) are mostly enriched in energy and cell division, immune, and neurological damage pathways. Further RNA-seq analysis identified a total of 1968 DEGs, of which 813 genes were upregulated and 1155 genes were downregulated. Based on the joint MethylRAD-seq and RNA-seq analysis, metabolic processes such as glycosaminoglycan degradation, oxidative phosphorylation, apoptosis, glutathione metabolism, thermogenesis, and lysosomes are regulated by DNA methylation.

CONCLUSIONS

High temperature affected the DNA methylation and expression levels of genes such as MOAP-1, GGT1 and RDH8, which in turn affects the metabolism of HPSE, Cox, glutathione, and retinol, thereby suppressing the immune, energy metabolism, and antioxidant functions of the organism and finally manifesting as stunted growth. In summary, the observations in the present study improve our understanding of the molecular mechanism of the response to high temperature stress in sea urchin.

Study on the Adaptive Regulation of Light on the Stress Response of Mandarin Fish (Siniperca chuatsi) with Re-Feeding after Starvation

Light influences the stress response to environmental stimuli and feeding behaviors of Siniperca chuatsi and, thus, is an important regulator of normal growth and development. In this study, we first explored the important role of light on the digestive and stress capacity of S. chuatsi by studying the changes in physiological and biochemical indicators of S. chuatsi, taking the re-feeding after starvation as the constant environmental stimulus and the light intensity as the adjustable environmental stimulus. The activity of protease and lipase was generally higher in the stomach tissues than in the intestinal tissues, especially lipase, which was higher in stomach tissues under all light conditions, and the protease and lipase activity peaked in the stomach tissues of S. chuatsi at a light intensity of 18.44 ± 3.00 lx and in intestinal tissues at 11.15 ± 2.01 lx, respectively, indicating that greater light intensity increased the digestive capacity of stomach tissues, whereas lower light intensity facilitated the digestive capacity of intestinal tissues. The tissues of the gill, stomach, and intestine had relatively high activity of stress-related enzymes, whereas the tissues of the brain, kidney, liver, and plasma samples had relatively low activity of enzymes. Collectively, the results show that light intensity at 11.15 ± 2.01 lx promoted digestive capacity in the intestine and enhanced the anti-stress ability of S. chuatsi in response to stress induced by re-feeding after starvation. These findings should prove useful for artificial breeding of S. chuatsi.

Cytochrome P450 superfamily in spotted sea bass: Genome-wide identification and expression profiles under trichlorfon and environmental stresses

Cytochrome P450s (CYPs), as one of the most diverse enzyme superfamilies in nature, play critical functions in antioxidant reactions against endogenous and exogenous compounds. In this study, we performed genome-wide characterization of CYP superfamily members and analyzed their expression patterns under several abiotic stresses in spotted sea bass, which is known as an economically important fish species in the Chinese aquaculture industry. A total of 55 CYP genes were identified and divided into 17 families within 10 clans. The analysis of phylogeny, gene structure, and syntenic relationships provided evidence for the evolution of CYP genes and confirmed their annotation and orthology. The expression of CYP genes was examined in the liver during trichlorfon stress using quantitative real-time PCR. The results showed that 20 tested CYP genes displayed significant mRNA expression changes, indicating that they may play crucial roles in the metabolism of trichlorfon and can be potential biomarkers for trichlorfon pollution. Moreover, by screening transcriptomic databases, 10, 3 and 19 CYP genes exhibited differential expression patterns in response to hypoxia, alkalinity and heat stress, respectively. Taken together, this study provided insights into the regulation of CYP genes by toxicological and environmental stresses, laid basis for extensive functional studies of the CYP superfamily in spotted sea bass and other teleost species.

Genome-Wide Identification and Characterization of Olfactory Receptor Genes in Silver Sillago (Sillago sihama)

Olfactory receptor (OR) genes are essential in the specific recognition of diverse stimuli in fish. In this study, a total of 141 OR genes were identified in silver sillago (Sillago sihama), a marine fish sensitive to environmental stimuli, including 112 intact genes, 26 truncated genes, and three pseudogenes. A phylogenetic tree analysis elucidated that the OR genes of S. sihama were classified into six groups, of which β, γ, δ, ε, and ζ groups belonged to type I, and the η group belonged to type II. The type I OR genes contained almost all conserved motifs (n = 62), while type II OR genes mainly retained conserved motifs 7(3), 1, 10, 4, and 2 (n = 39). OR genes were mainly distributed on LG1, LG9, LG11, and LG12. Of all OR genes, 36.23% (50 genes) showed significant expansion in S. sihama. Ka/Ks analysis demonstrated that 227 sites were under purifying selection, while 12 sites were under positive selection, including eight genes in the OR2A12 gene subfamily. Sixty-one genes (44.20%) displayed differential expression under hypoxic stress. The identified OR genes explored the mechanism of environmental stress and ecological adaptation of S. sihama, and provided valuable genomic resources for further research on the olfaction of teleosts.

RNA-seq and qRT-PCR analyses reveal the physiological response to acute hypoxia and reoxygenation in Epinephelus coioides

Hypoxia is a critical problem in intensive Epinephelus coioides aquaculture systems. In the present study, the physiological responses of E. coioides muscle to acute hypoxic stress (DO = 0.6 ± 0.1 mg/L) and reoxygenation (DO = 6.0 ± 0.1 mg/L) were analyzed by transcriptome sequencing (RNA-seq) and quantitative real-time PCR (qRT–PCR). RNA-seq was conducted on the muscle tissues of E. coioides in the hypoxia-tolerant (EMS), hypoxia-sensitive (EMW), and normoxic (CM) groups. Among the three groups, a total of 277 differentially expressed genes (DEGs) were identified. KEGG analysis revealed that the pathways significantly enriched after hypoxic stress are involved in the immune response, glycolysis/gluconeogenesis, energy metabolism, vasodilation and proliferation, cell proliferation, and apoptosis. qRT‒PCR verified that the differentially expressed genes FIH-1, PHD-2, PPARα, BCL-XL, LDH-A, and Flt-1 were significantly upregulated after hypoxic stress and returned to normal levels after reoxygenation, suggesting that these DEGs play important roles in responding to hypoxia treatment. In addition, the HIF-1 signaling pathway was also activated under hypoxic stress, and qRT‒PCR confirmed that the expression level of HIF-1α was significantly elevated under acute hypoxic stress, indicating that the HIF-1 signaling pathway is the central pathway in the E. coioides hypoxic response mechanism and activates other related pathways to adapt to hypoxic stress. These pathways jointly regulate energy metabolism, substance synthesis, blood vessel proliferation, cell proliferation, and differentiation and prolong survival time. These results provide ideas for understanding physiological regulation after hypoxic stress and reoxygenation and provide basic insights for the future breeding of hypoxia-tolerant E. coioides.

PacBio Full-Length and Illumina Transcriptomes of the Gill Reveal the Molecular Response of Corbicula fluminea under Aerial Exposure

Air exposure is a common stress for Corbicula fluminea, an economically important freshwater shellfish consumed in China, during aquaculture and transportation. However, little is known about its molecular responses to air exposure. Therefore, this study used a combination of PacBio full-length and Illumina transcriptomes to investigate its molecular responses to air exposure. A total of 36,772 transcripts were obtained using PacBio sequencing. Structural analysis identified 32,069 coding sequences, 1906 transcription factors, 8873 simple sequence repeats, and 17,815 long non-coding RNAs. Subcellular localization analysis showed that most transcripts were located in the cytoplasm and nucleus. After 96-h of air exposure, 210 differentially expressed genes (DEGs) in the gill were obtained via Illumina sequencing. Among these DEGs, most of the genes related to glycolysis, tricarboxylic acid cycle, lipid metabolism, and amino acid metabolism were upregulated. Additionally, many DEGs associated with immunity, cytoskeleton reorganization, autophagy, and ferroptosis were identified. These findings indicated that metabolic strategy change, immune response, cytoskeleton reconstruction, autophagy, and ferroptosis might be the important mechanisms that C. fluminea use to cope with air exposure. This study will enrich the gene resources of C. fluminea and provide valuable data for studying the molecular mechanisms coping with air exposure in C. fluminea and other freshwater mollusks.

Comparative transcriptome analysis provides novel insights into the molecular mechanism of the silver carp (Hypophthalmichthys molitrix) brain in response to hypoxia stress

The brain of fish plays an important role in regulating growth and adapting to environmental changes. However, few studies have been performed to address the changes in gene expression profiles in fish brains under hypoxic stress. In the present study, silver carp (Hypophthalmichthys molitrix) were kept under hypoxic experimental conditions by using the method of natural oxygen consumption, which resulted in a significant decrease in malondialdehyde (MDA) and glutathione (GSH) content and superoxide dismutase (SOD) activity in the brain. In addition, RNA sequencing (RNA-Seq) was performed to analyze transcriptional regulation in the brains of silver carp under normoxia (control group), hypoxia, semi-asphyxia, and asphyxia conditions. The results of KEGG enrichment pathway analysis showed that the immune system, such as antigen processing and presentation, natural killer cell-mediated cytotoxicity, was enriched in the hypoxia group; the nervous system (e.g., "glutamatergic synapse"), signal transduction (e.g., "calcium signaling pathway"; "foxo signaling pathway"), and signaling molecules and interactions (e.g., "neuroactive ligand-receptor interaction") were enriched in the semi-asphyxia group; and signaling molecules and interactions (e.g., "neuroactive ligand-receptor interaction") were enriched in the asphyxia group. These results provide novel insights into the molecular regulatory mechanism of the fish brain coping with hypoxia stress.

Gene expression profiles in Malpighian tubules of the vector leafhopper Psammotettix striatus (L.) revealed regional functional diversity and heterogeneity

Background

Many leafhoppers are known as pests and disease vectors of economically important plants. Previous studies of the physiological functions of vector leafhoppers have mainly focused on the salivary glands and the alimentary tract that are deemed to be associated with digestion, host defense and phytoplasma and/or virus transmission. By contrast, the significance of Malpighian tubules (MTs) is less studied. To clarify the physiological function of MTs of the vector leafhopper Psammotettix striatus that transmits phytoplasma triggering the wheat blue dwarf disease, we performed a transcriptome study on P. striatus MTs and compared gene expression profiles among different anatomical regions in the tubules (i.e., MT1+2, the anterior segment together with the sub-anterior segment; MT3, the inflated segment; and MT4, the distal segment).

Results

Transcriptome of P. striatus MTs generate a total of 42,815 high-quality unigenes, among which highly expressed unigenes are mainly involved in organic solute transport, detoxification and immunity in addition to osmoregulation. Region-specific comparative analyses reveal that all these MT regions have functions in osmoregulation, organic solute transport and detoxification, but each region targets different substrates. Differential expression and regional enrichment of immunity-related effector activities and molecules involved in phagocytosis and the biosynthesis of antimicrobial peptides among different regions indicate that MT1+2 and MT4 have the ability to eliminate the invading pathogens. However, in MT3 which secrets brochosomes to the integument and eggs as physical barriers, disulfide-isomerase, acidic ribosomal protein P and many other unigenes were highly expressed, which can be attractive candidate genes for future studies of the biosynthesis and the origin of brochosomes.

Conclusions

Psammotettix striatus MTs perform multiple physiological functions as versatile organs than just excretory organs with osmoregulatory function. Heterogeneity of physiological functions among different MT regions is related to organic solute transport, detoxification, immunity and brochosome biosynthesis in addition to osmoregulation, and each region targets different substrates. These functions may be helpful for P. striatus to resist pathogens from habitats and to utilize a wider range of host plants, which may assist the transmission and spread of phytoplasmas. The results provide potential molecular targets for the exploit of chemical and/or gene-silencing insecticides.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08300-6.

Physiological Effects and Transcriptomic Analysis of sbGnRH on the Liver in Pompano (Trachinotus ovatus)

Pompano (Trachinotus ovatus) is one of the important economic marine fishes in the south coast of China. At present, the research on the basic biology of pompano is relatively weak, which has seriously affected the development of this economic important fish. The liver is an important digestive and metabolic organ of fish which plays an important regulatory role in its growth and development. It is necessary to clarify the effects of sea bream gonadotropin releasing hormone (sbGnRH) on liver physiology and metabolic enzyme activity. The effects of sbGnRH peptides (10 ng/gbw) on the physiological and biochemical indices and metabolic enzyme activities of pompano liver were studied. It was found that after injection of 10 ng/gbw sbGnRH peptides, the contents of albumin, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, glucose, creatine kinase, iron, magnesium, aspartate aminotransferase, alanine aminotransferase and creatinine increased, while of cholesterol and calcium contents decreased. The activities of amylase, lipase, pyruvate kinase, acyl CoA oxidase, superoxide dismutase, phospholipid hydroperoxide glutathione peroxidase, catalase, glucose-6-phosphate dehydrogenase, fatty acid synthase and lipoprotein lipase increased, while the activities of malic enzyme, carnitine acyl, carnitine translocation, acetyl CoA carboxylase and malondialdehyde decreased. Three hours after the injection of different concentrations of sbGnRH peptides (0 and 10 ng/gbw), the transcriptome sequences of the two groups of livers were sequenced. After quality control and removal of some low-quality data, clean reads of 21,283,647、19,427,359、21,873,990、21,732,174、23,660,062 and 21,592,338 were obtained respectively. In this study, 99 genes were screened and identified as differentially expressed genes, including 77 up-regulated genes and 22 down-regulated genes. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analyses, these pathways and the typical genes involved can be divided into cellular processes, environmental information processing, genetic information processing, diseases, metabolism and organismal systems. The results from this study provide a the oretical basis for studying the effects of sbGnRH on the physiology, biochemistry and metabolic enzyme activities of liver in pompano.

Transcriptomic analysis of pituitary in female and male spotted scat (Scatophagus argus) after 17β-estradiol injection

Spotted scat (Scatophagus argus) is a popular species of marine fish cultured in China. It shows normal sexual growth dimorphism. Female spotted scat grows quicker and bigger than males. Growth and reproduction are the most important traits in aquaculture. In vertebrates, the pituitary gland occupies an important position in the growth and reproduction axis. Estrogen is involved in regulating growth and reproduction in the pituitary gland in an endocrine fashion. Transcriptome sequencing of the pituitary was performed in female and male fish at 6 h after 17β-estradiol injection (4.0 μg E₂/g body weight, BW). Compared with the pituitary of female and male groups, 144 and 64 genes [|log₂^{(fold change)}| ≥ 1.0 and false discovery rate (FDR) < 0.05] were significantly differentially expressed in E₂-injected females and males, respectively (p < 0.05). Of these, 59 and 48 were up-regulated, and 85 and 16 were down-regulated. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analyses, DEGs were involved in signal pathways, such as growth, reproduction, oocyte meiosis and steroid biosynthesis. Of these, estrogen affected the expression of some sex steroid synthesis and receptor genes in the pituitary gland through feedback, such as hsd17b7, pgr and cyp19a1b, regulating the reproductive activities. Besides, some growth-related genes, such as gap43, junbb, mstn2 and insm1a responded to estrogen. E₂ might affect the expression level of gh mRNA by regulating the expression levels of growth-related genes. Our results provide a theoretical basis for studying the molecular mechanism of growth and reproduction regulation at the pituitary level of spotted scat responded to E₂.

Transcriptome profiling reveals differential expression of immune-related genes in gills of hybrid yellow catfish (Tachysurus fulvidraco ♀ × Pseudobagrus vachellii ♂) under hypoxic stress: Potential NLR-mediated immune response

Fish gills are the primary organ that respond to sudden changes in the dissolved oxygen (DO) level in the aquatic environment. Hypoxic stress impairs the normal function of gill tissues. However, little is known about the mechanisms of the response of yellow catfish gills to hypoxic stress. In this study, we compared transcriptomic and physiological changes in gill tissues of hybrid yellow catfish (Tachysurus fulvidraco ♀ × Pseudobagrus vachellii ♂) between a hypoxia-treated group (DO: 1.5 mg/L) and a control group (DO: 6.5 mg/L). In fish in the hypoxia-treated group, gill filaments underwent adaptive changes, and the number of vacuoles in gill tissues increased. Exposure to hypoxic conditions for 96 h resulted in increased anaerobic metabolism and decreased antioxidant and immune capacity in gill tissues. Transcriptome analyses revealed 1556 differentially expressed genes, including 316 up-regulated and 1240 down-regulated genes, between fish in the hypoxia-treated and control groups. Functional analyses indicated that the main pathway enriched with differentially expressed genes was immune response, followed by energy metabolism and signal transduction. Under hypoxic stress, the transcript levels of genes involved in the NOD-like receptor signaling pathway initially increased rapidly but then decreased over time, suggesting that the NOD-like receptor-mediated immune response plays an essential role in hypoxia tolerance and resistance in hybrid yellow catfish. Our results provide novel insights into which immune-related genes and pathways are activated under hypoxic stress, and reveal details of early adaptation of the immune response and defense mechanisms under hypoxic stress.

Genome-wide identification and expression profiling of glutathione S-transferase family under hypoxia stress in silver sillago (Sillago sihama)

Glutathione S-transferase (GST) is an important detoxification enzyme in organisms. GSTs play an important role in responding to environmental stresses. This study aimed to identify the GST gene superfamily in silver sillago (Sillago sihama) and analyze its expression pattern under hypoxia stress. A total of 17 GST genes were identified in silver sillago. Phylogenetic analysis showed that the GST gene family contained two subgroups (cytosolic and MAPEGs), and lacked three subgroups (i.e. Pi, Kappa, and MGST2). Phylogenetic and syntenic analysis revealed that GST genes were conserved in evolution. Eight SsGSTs were significantly differentially expressed under hypoxia stress in silver sillago by RNA-seq and qRT-PCR analysis. The expression levels of SsMGST3b, SsGSTO1, SsGSTT1b and SsGSTR2 genes were significantly up-regulated after 4 h of reoxygenation in the gill tissue. In the heart tissue, the expression of SsGSTR3 was significantly up-regulated after 1 h of hypoxia while the expression levels of SsGSTT1b and SsFLAP genes were significantly down-regulated after 4 h of hypoxia. In summary, this study provides for the first time a comprehensive analysis of the GST gene superfamily of silver sillago.

Comparative transcriptomic analysis reveals an association of gibel carp fatty liver with ferroptosis pathway

BACKGROUND

Fatty liver has become a main problem that causes huge economic losses in many aquaculture modes. It is a common physiological or pathological phenomenon in aquaculture, but the causes and occurring mechanism are remaining enigmatic.

METHODS

Each three liver samples from the control group of allogynogenetic gibel carp with normal liver and the overfeeding group with fatty liver were collected randomly for the detailed comparison of histological structure, lipid accumulation, transcriptomic profile, latent pathway identification analysis (LPIA), marker gene expression, and hepatocyte mitochondria analyses.

RESULTS

Compared to normal liver, larger hepatocytes and more lipid accumulation were observed in fatty liver. Transcriptomic analysis between fatty liver and normal liver showed a totally different transcriptional trajectory. GO terms and KEGG pathways analyses revealed several enriched pathways in fatty liver, such as lipid biosynthesis, degradation accumulation, peroxidation, or metabolism and redox balance activities. LPIA identified an activated ferroptosis pathway in the fatty liver. qPCR analysis confirmed that gpx4, a negative regulator of ferroptosis, was significantly downregulated while the other three positively regulated marker genes, such as acsl4, tfr1 and gcl, were upregulated in fatty liver. Moreover, the hepatocytes of fatty liver had more condensed mitochondria and some of their outer membranes were almost ruptured.

CONCLUSIONS

We reveal an association between ferroptosis and fish fatty liver for the first time, suggesting that ferroptosis might be activated in liver fatty. Therefore, the current study provides a clue for future studies on fish fatty liver problems.

Chromosomal-Level Genome Assembly of Silver Sillago (Sillago sihama)

Silver sillago, Sillago sihama is a member of the family Sillaginidae and found in all Chinese inshore waters. It is an emerging commercial marine aquaculture species in China. In this study, high-quality chromosome-level reference genome of S. sihama was first constructed using PacBio Sequel sequencing and high-throughput chromosome conformation capture (Hi-C) technique. A total of 66.16 Gb clean reads were generated by PacBio sequencing platforms. The genome-scale was 521.63 Mb with 556 contigs, and 13.54 Mb of contig N50 length. Additionally, Hi-C scaffolding of the genome resulted in 24 chromosomes containing 96.93% of the total assembled sequences. A total of 23,959 protein-coding genes were predicted in the genome, and 96.51% of the genes were functionally annotated in public databases. A total of 71.86 Mb repetitive elements were detected, accounting for 13.78% of the genome. The phylogenetic relationships of silver sillago with other teleosts showed that silver sillago was separated from the common ancestor of Sillago sinica ∼7.92 Ma. Comparative genomic analysis of silver sillago with other teleosts showed that 45 unique and 100 expansion gene families were identified in silver sillago. In this study, the genomic resources provide valuable reference genomes for functional genomics research of silver sillago.