Genome sequencing of the perciform fish Larimichthys crocea provides insights into molecular and genetic mechanisms of stress adaptation

Effects of Gene Alternative Splicing Events on Resistance to Cryptocaryonosis of Large Yellow Croaker (Larimichthys crocea)

Large yellow croaker (L. crocea) is a productive species in marine aquaculture with great economic value in China. However, the sustainable development of large yellow croaker is hampered by various diseases including cryptocaryonosis caused by Cryptocaryon irritans. The genetic regulation processes for cryptocaryonosis in large yellow croaker are still unclear. In this present study, we analyzed differential alternative splicing events between a C. irritans resistance strain (RS) and a commercial strain (CS). We identified 678 differential alternative splicing (DAS) events from 453 genes in RS and 719 DAS events from 500 genes in CS. A set of genes that are specifically alternatively spliced in RS was identified including mfap5, emp1, and trim33. Further pathway analysis revealed that the specifically alternative spliced genes in RS were involved in innate immune responses through the PRR pathway and the Toll and Imd pathway, suggesting their important roles in the genetic regulation processes for cryptocaryonosis in large yellow croaker. This study would be helpful for the studies of the pathogenesis of cryptocaryonosis and dissection of C. irritans resistance for L. crocea.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

Aggregation and proliferation of B cells and T cells in MALTs upon Cryptocaryon irritans infection in large yellow croaker Larimichthys crocea

Mucosal immunity in mucosa-associated lymphoid tissues (MALTs) plays crucial roles in resisting infection by pathogens, including parasites, bacteria and viruses. However, the mucosal immune response in the MALTs of large yellow croaker (Larimichthys crocea) upon parasitic infection remains largely unknown. In this study, we investigated the role of B cells and T cells in the MALTs of large yellow croaker following Cryptocaryon irritans infection. Upon C. irritans infection, the total IgM and IgT antibody levels were significantly increased in the skin mucus and gill mucus. Notably, parasite-specific IgM antibody level was increased in the serum, skin and gill mucus following parasitic infection, while the level of parasite-specific IgT antibody was exclusively increased in MALTs. Moreover, parasitic infection induced both local and systemic aggregation and proliferation of IgM+ B cells, suggesting that the increased levels of IgM in mucus may be derived from both systemic and mucosal immune tissues. In addition, we observed significant aggregation and proliferation of T cells in the gill, head kidney and spleen, suggesting that T cells may also be involved in the systemic and mucosal immune responses upon parasitic infection. Overall, our findings provided further insights into the role of immunoglobulins against pathogenic infection, and the simultaneous aggregation and proliferation of both B cells and T cells at mucosal surfaces suggested potential interactions between these two major lymphocyte populations during parasitic infection.

Structural and Functional Characterization of a Fish Type I Subgroup d IFN Reveals Its Binding to Receptors

Teleost fish type I IFNs and the associated receptors from the cytokine receptor family B (CRFB) are characterized by remarkable diversity and complexity. How the fish type I IFNs bind to their receptors is still not fully understood. In this study, we demonstrate that CRFB1 and CRFB5 constitute the receptor pair through which type I subgroup d IFN (IFNd) from large yellow croaker, Larimichthys crocea, activates the conserved JAK-STAT signaling pathway as a part of the antiviral response. Our data suggest that L. crocea IFNd (LcIFNd) has a higher binding affinity with L. crocea CRFB5 (LcCRFB5) than with LcCRFB1. Furthermore, we report the crystal structure of LcIFNd at a 1.49-Å resolution and construct structural models of LcIFNd in binary complexes with predicted structures of extracellular regions of LcCRFB1 and LcCRFB5, respectively. Despite striking similarities in overall architectures of LcIFNd and its ortholog human IFN-ω, the receptor binding patterns between LcIFNd and its receptors show that teleost and mammalian type I IFNs may have differentially selected helices that bind to their homologous receptors. Correspondingly, key residues mediating binding of LcIFNd to LcCRFB1 and LcCRFB5 are largely distinct from the receptor-interacting residues in other fish and mammalian type I IFNs. Our findings reveal a ligand/receptor complex binding mechanism of IFNd in teleost fish, thus providing new insights into the function and evolution of type I IFNs.

Biodirected Screening and Preparation of Larimichthys crocea Angiotensin-I-Converting Enzyme-Inhibitory Peptides by a Combined In Vitro and In Silico Approach

Food-derived angiotensin-I-converting enzyme (ACE)-inhibitory peptides have gained attention for their potent and safe treatment of hypertensive disorders. However, there are some limitations of conventional methods for preparing ACE-inhibitory peptides. In this study, in silico hydrolysis, the quantitative structure-activity relationship (QSAR) model, LC-MS/MS, inhibition kinetics, and molecular docking were used to investigate the stability, hydrolyzability, in vitro activity, and inhibition mechanism of bioactive peptides during the actual hydrolysis process. Six novel ACE-inhibitory peptides were screened from the Larimichthys crocea protein (LCP) and had low IC50 values (from 0.63 ± 0.09 µM to 10.26 ± 0.21 µM), which were close to the results of the QSAR model. After in vitro gastrointestinal simulated digestion activity of IPYADFK, FYEPFM and NWPWMK were found to remain almost unchanged, whereas LYDHLGK, INEMLDTK, and IHFGTTGK were affected by gastrointestinal digestion. Meanwhile, the inhibition kinetics and molecular docking results were consistent in that ACE-inhibitory peptides of different inhibition forms could effectively bind to the active or non-central active centers of ACE through hydrogen bonding. Our proposed method has better reproducibility, accuracy, and higher directivity than previous methods. This study can provide new approaches for the deep processing, identification, and preparation of Larimichthys crocea.

DNA methylation regulates B cell activation via repressing Pax5 expression in teleost

In mammals, the transcription factor Pax5 is a key regulator of B cell development and maturation and specifically expressed in naive/mature B cells but repressed upon B cell activation. Despite the long-standing proposal that Pax5 repression is essential for proper B cell activation, the underlying mechanisms remain largely elusive. In this study, we used a teleost model to elucidate the mechanisms governing Pax5 repression during B cell activation. Treatment with lipopolysaccharide (LPS) and chitosan oligosaccharide (COS) significantly enhanced the antibody secreting ability and phagocytic capacity of IgM+ B cells in large yellow croaker (Larimichthys crocea), coinciding with upregulated expression of activation-related genes, such as Bcl6, Blimp1, and sIgM, and downregulated expression of Pax5. Intriguingly, two CpG islands were identified within the promoter region of Pax5. Both CpG islands exhibited hypomethylation in naive/mature B cells, while CpG island1 was specifically transited into hypermethylation upon B cell activation. Furthermore, treatment with DNA methylation inhibitor 5-aza-2'-deoxycytidine (AZA) prevented the hypermethylation of CpG island1, and concomitantly impaired the downregulation of Pax5 and activation of B cells. Finally, through in vitro methylation experiments, we demonstrated that DNA methylation exerts an inhibitory effect on promoter activities of Pax5. Taken together, our findings unveil a novel mechanism underlying Pax5 repression during B cell activation, thus promoting the understanding of B cell activation process.

The identification and expression of an interleukin-21 receptor in large yellow croaker (Larimichthys crocea)

BACKGROUND

We identified a homologue of IL-21R (LcIL-21R) in large yellow croaker (Larimichthys crocea, Lc). Our investigation focused on understanding the molecular structural features and immune function of LcIL-21R.

METHODS

We cloned the LcIL-21R gene from the genome of Larimichthys crocea by RT‒PCR, and the molecular and structural characteristics of LcIL-21R were analyzed by a series of protein analysis tools. We used real-time PCR to investigate the tissue distribution of LcIL-21R, and LcIL-21R gene expression regulation was also measured in head kidney leukocytes under trivalent bacterial vaccine or poly (I:C) stimulation.

RESULTS

The open reading frame (ORF) of the LcIL-21R gene is 1629 bp long and encodes a precursor protein of 542 amino acids (aa), with a 23-aa signal peptide and a 519-aa mature peptide containing four putative N-glycosylation sites. LcIL-21R has two fibronectin type III (FNIII)-like domains (D1 and D2), a transmembrane domain, and a cytoplasmic region. A conserved WSXWS motif was also found in the D2 domain. The predicted structure of the extracellular region of LcIL-21R (LcIL-21R-Ex) is highly similar to that of human IL-21R. LcIL-21R was constitutively expressed in all tissues examined, and LcIL-21R mRNA levels were increased in the head kidney and spleen upon inactivated trivalent bacterial vaccine or poly(I:C) stimulation.

CONCLUSIONS

Our results suggest that LcIL-21R shares structural and functional properties with IL-21Rs found in other vertebrates, indicating its potential involvement in the IL-21-mediated immune response to pathogenic infections. These findings contribute to our understanding of the evolutionary conservation of IL-21 signaling and its role in the immune system.

Improved high-quality reference genome of red drum facilitates the processes of resistance-related gene exploration

Sciaenops ocellatus is among the most important artificially introduced farmed fish across 11 countries and regions. However, the frequent occurrence of extreme weather events and breeding escapes have placed great pressure on local marine biodiversity and ecosystems. We reported the de novo assembly and annotation with a contig N50 of 28.30 Mb using PacBio HiFi sequencing and Hi-C technologies, which resulted in a 283-fold increase in contig N50 length and improvement in continuity and quality in complex repetitive region for S. ocellatus compared to the previous version. In total, 257.36 Mb of repetitive sequences accounted for 35.48% of the genome, and 22,845 protein-coding genes associated with a BUSCO value of 98.32%, were identified by genome annotation. Moreover, 54 hub genes rapidly responding to hypoosmotic stress were identified by WGCNA. The high-quality chromosome-scale S. ocellatus genome and candidate resistance-related gene sets will not only provide a genomic basis for genetic improvement via molecular breeding, but will also lay an important foundation for investigating the molecular regulation of rapid responses to stress.

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

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

Large yellow croaker (Larimichthys crocea) mitofusin 2 inhibits type I IFN responses by degrading MAVS via enhanced K48-linked ubiquitination

In mammals, mitofusin 2 (MFN2) is involved in mitochondrial fusion, and suppresses the virus-induced RIG-I-like receptor (RLR) signaling pathway. However, little is known about the function of MFN2 in non-mammalian species. In the present study, we cloned an MFN2 ortholog (LcMFN2) in large yellow croaker (Larimichthys crocea). Phylogenetic analysis showed that MFN2 emerged after the divergence of amphioxus and vertebrates. The protein sequences of MFN2 were well conserved from fish to mammals. LcMFN2 was expressed in all the tissues/organs examined at different levels, and its expression was upregulated in response to poly(I:C) stimulation. Overexpression of LcMFN2 inhibited MAVS-induced type I interferon (IFN) promoter activation and antiviral gene expression. In contrast, knockdown of endogenous LcMFN2 enhanced poly(I:C) induced production of type I IFNs. Additionally, LcMFN2 enhanced K48-linked polyubiquitination of MAVS, promoting its degradation. Also, overexpression of LcMFN2 impaired the cellular antiviral response, as evidenced by the increased expression of viral genes and more severe cytopathic effects (CPE) in cells infected with spring viremia of carp virus (SVCV). These results indicated that LcMFN2 inhibited type I IFN response by degrading MAVS, suggesting its negative regulatory role in cellular antiviral response. Therefore, our study sheds a new light on the regulatory mechanisms of the cellular antiviral response in teleosts.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00189-8.

Genome-Wide RADseq Reveals Genetic Differentiation of Wild and Cultured Populations of Large Yellow Croaker

Larimichthys crocea (also known as the large yellow croaker) is one of the most economically important marine fishes in China, and research on the ecology and genetics of this species is of immense significance. In this study, we performed restriction site-associated DNA sequencing (RAD-seq) of 54 individuals collected from four sites in China to analyze the genetic structure and diversity of large yellow croaker at the genome level. It revealed that the large yellow croaker populations in the Ningde and Zhoushan coastal waters can be clearly distinguished. Different genetic diversity indices were used to analyze the genetic diversity of the large yellow croaker, which showed that there was a differentiation trend between the wild and farmed populations in Ningde. Moreover, we identified genetically differentiated genomic regions between the populations. GO gene enrichment analysis identified genes that are related to fatty acid metabolism and growth. These findings enhance our understanding of genetic differentiation and adaptation to different living environments, providing a theoretical basis for the preservation and restoration of the genetic resources of the large yellow croaker.

Fish Skin Mucus Extracts: An Underexplored Source of Antimicrobial Agents

The slow discovery of new antibiotics combined with the alarming emergence of antibiotic-resistant bacteria underscores the need for alternative treatments. In this regard, fish skin mucus has been demonstrated to contain a diverse array of bioactive molecules with antimicrobial properties, including peptides, proteins, and other metabolites. This review aims to provide an overview of the antimicrobial molecules found in fish skin mucus and its reported in vitro antimicrobial capacity against bacteria, fungi, and viruses. Additionally, the different methods of mucus extraction, which can be grouped as aqueous, organic, and acidic extractions, are presented. Finally, omic techniques (genomics, transcriptomics, proteomics, metabolomics, and multiomics) are described as key tools for the identification and isolation of new antimicrobial compounds. Overall, this study provides valuable insight into the potential of fish skin mucus as a promising source for the discovery of new antimicrobial agents.

Molecular and functional characterization of two IgM subclasses in large yellow croaker (Larimichthys crocea)

As the predominant immunoglobulin (Ig) isotype, IgM plays a crucial role in the acquired immunity of vertebrates. There is only one Igμ gene in mammals, except cattle, while the number of Igμ gene varies among teleost fish. In the current study, we found two functional Igμ genes (Igμ1 and Igμ2) and a pseudo Cμ gene (ψIgμ) in large yellow croaker (Larimichthys crocea). Both Igμ1 and Igμ2 genes possessed two transcript variants, which encoded the heavy chains of secreted (sIgM1 and sIgM2) and membrane-bound IgM1 and IgM2 (mIgM1 and mIgM2), respectively. Both the heavy chains of sIgM1 and sIgM2 consisted of a variable Ig domain, four constant Ig domains (CH1, CH2, CH3 and CH4) and a secretory tail, while those of mIgM1 and mIgM2 consisted of a variable Ig domain, three constant Ig domains (CH1, CH2 and CH3), a transmembrane domain and a short cytoplasmic tail. Cysteine residues that are necessary for the formation of intrachain and interchain disulfide bonds and tryptophan residues that are important for the folding of the Ig superfamily domain were well conserved in large yellow croaker IgM1 and IgM2. Interestingly, large yellow croaker IgM2 had an extra cysteine (C94) in the CH1 domain compared with IgM1, which may cause the structural difference between IgM1 and IgM2. A liquid chromatography-tandem mass spectrometry analysis revealed that both IgM1 and IgM2 were present at the protein level in large yellow croaker serum. Both the Igμ1 and Igμ2 genes were mainly expressed in systemic immune tissues, such as head kidney and spleen, but the expression level of Igμ2 was much lower than that of Igμ1. After Pseudomonas plecoglossicida infection, the expression levels of Igμ1 and Igμ2 in both the spleen and head kidney were significantly upregulated, with a higher upregulation of Igμ2 than that of Igμ1. These results suggested that Igμ1 and Igμ2 may play a differential role in the immune response of large yellow croaker against bacterial infection.

Chromosome genome assembly for the meagre, Argyrosomus regius, reveals species adaptations and sciaenid sex-related locus evolution

The meagre, Argyrosomus regius, has recently become a species of increasing economic interest for the Mediterranean aquaculture and there is ongoing work to boost production efficiency through selective breeding. Access to the complete genomic sequence will provide an essential resource for studying quantitative trait-associated loci and exploring the genetic diversity of different wild populations and aquaculture stocks in more detail. Here, we present the first complete genome for A. regius, produced through a combination of long and short read technologies and an efficient in-house developed pipeline for assembly and polishing. Scaffolding using previous linkage map data allowed us to reconstruct a chromosome level assembly with high completeness, complemented with gene annotation and repeat masking. The 696 Mb long assembly has an N50 = 27.87 Mb and an L50 = 12, with 92.85% of its length placed in 24 chromosomes. We use this new resource to study the evolution of the meagre genome and other Sciaenids, via a comparative analysis of 25 high-quality teleost genomes. Combining a rigorous investigation of gene duplications with base-wise conservation analysis, we identify candidate loci related to immune, fat metabolism and growth adaptations in the meagre. Following phylogenomic reconstruction, we show highly conserved synteny within Sciaenidae. In contrast, we report rapidly evolving syntenic rearrangements and gene copy changes in the sex-related dmrt1 neighbourhood in meagre and other members of the family. These novel genomic datasets and findings will add important new tools for aquaculture studies and greatly facilitate husbandry and breeding work in the species.

Effect of dietary phenylalanine on growth performance and intestinal health of triploid rainbow trout (Oncorhynchus mykiss) in low fishmeal diets

This study aimed to investigate the effects of phenylalanine on the growth, digestive capacity, antioxidant capability, and intestinal health of triploid rainbow trout (Oncorhynchus mykiss) fed a low fish meal diet (15%). Five isonitrogenous and isoenergetic diets with different dietary phenylalanine levels (1.82, 2.03, 2.29, 2.64, and 3.01%) were fed to triplicate groups of 20 fish (initial mean body weight of 36.76 ± 3.13 g). The weight gain rate and specific growth rate were significantly lower (p < 0.05) in the 3.01% group. The trypsin activity in the 2.03% group was significantly higher than that in the control group (p < 0.05). Amylase activity peaked in the 2.64% treatment group. Serum superoxide dismutase, catalase, and lysozyme had the highest values in the 2.03% treatment group. Liver superoxide dismutase and catalase reached their maximum values in the 2.03% treatment group, and lysozyme had the highest value in the 2.29% treatment group. Malondialdehyde levels in both the liver and serum were at their lowest in the 2.29% treatment group. Interleukin factors IL-1β and IL-6 both reached a minimum in the 2.03% group and were significantly lower than in the control group, while IL-10 reached a maximum in the 2.03% group (p < 0.05). The tight junction protein-related genes occludin, claudin-1, and ZO-1 all attained their highest levels in the 2.03% treatment group and were significantly higher compared to the control group (p < 0.05). The intestinal villi length and muscle layer thickness were also improved in the 2.03% group (p < 0.05). In conclusion, dietary phenylalanine effectively improved the growth, digestion, absorption capacity, antioxidant capacity, and intestinal health of O. mykiss. Using a quadratic curve model analysis based on WGR, the dietary phenylalanine requirement of triploid O. mykiss fed a low fish meal diet (15%) was 2.13%.

The genome of a hadal sea cucumber reveals novel adaptive strategies to deep-sea environments

How organisms cope with coldness and high pressure in the hadal zone remains poorly understood. Here, we sequenced and assembled the genome of hadal sea cucumber Paelopatides sp. Yap with high quality and explored its potential mechanisms for deep-sea adaptation. First, the expansion of ACOX1 for rate-limiting enzyme in the DHA synthesis pathway, increased DHA content in the phospholipid bilayer, and positive selection of EPT1 may maintain cell membrane fluidity. Second, three genes for translation initiation factors and two for ribosomal proteins underwent expansion, and three ribosomal protein genes were positively selected, which may ameliorate the protein synthesis inhibition or ribosome dissociation in the hadal zone. Third, expansion and positive selection of genes associated with stalled replication fork recovery and DNA repair suggest improvements in DNA protection. This is the first genome sequence of a hadal invertebrate. Our results provide insights into the genetic adaptations used by invertebrate in deep oceans.

Enhance energy supply of largemouth bass (Micropterus salmoides) in gills during acute hypoxia exposure

Gills are the location of gas exchange and also the first target organ of fish response for environmental stress. As a multifunctional organ, its energy supply, when faced with insufficient dissolved oxygen in the water, remains unclear. In this study, largemouth bass was subjected to hypoxia stress (1.2 mg/L) for 24 h and 12 h reoxygenation (R12) to evaluate energy supply strategy of gills. Under hypoxia exposure, the respiratory rate of largemouth bass increased by an average of 20 breaths per minute. A total of 2026, 1744, 1003, 579, 485, and 265 differentially expressed genes (DGEs) were identified at 0 h, 4 h, 8 h, 12 h, 24 h, and R12h in gills after hypoxia exposure. KEGG functional analysis of DEGs revealed that the glycolysis/gluconeogenesis pathway was enriched across all the sampling points (0, 4, 8, 12, 24 h, R12). The gene expression and enzyme activity of three rate-limiting enzymes (hexokinase, phosphofructokinase-6, pyruvate kinase) in glycolysis pathway were significantly increased. Increased levels of glycolysis products pyruvate and lactic acid, as well as the number of mitochondria (1.8-fold), suggesting an enhancement of aerobic and anaerobic metabolism of glucose in gills. These results suggest that the gill of largemouth bass enhanced the energy supply during acute exposure to hypoxia stress.

Transcriptomic Analysis of Large Yellow Croaker (Larimichthys crocea) during Early Development under Hypoxia and Acidification Stress

Simple Summary

The large yellow croaker is one of the most economically important fish in China. In recent years, the deterioration of the water environment and unregulated aquaculture have caused great economic losses to the large yellow croaker breeding industry. The aim of this study was to analyze the effects of hypoxia and acidification stress on large yellow croaker. This study revealed that hypoxia and acidification stress suppressed the growth of the large yellow croaker. Transcriptome analysis revealed that genes of the collagen family play an important role in the response of large yellow croaker to hypoxia and acidification stress. The study elucidates the mechanism underlying the response of large yellow croaker to hypoxia–acidification stress during early development and provides a basic understanding of the potential combined effects of reduced pH and dissolved oxygen on Sciaenidae fishes.

Abstract

Fishes live in aquatic environments and several aquatic environmental factors have undergone recent alterations. The molecular mechanisms underlying fish responses to hypoxia and acidification stress have become a serious concern in recent years. This study revealed that hypoxia and acidification stress suppressed the growth of body length and height of the large yellow croaker (Larimichthys crocea). Subsequent transcriptome analyses of L. crocea juveniles under hypoxia, acidification, and hypoxia–acidification stress led to the identification of 5897 differentially expressed genes (DEGs) in the five groups. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that several DEGs were enriched in the ‘protein digestion and absorption’ pathway. Enrichment analysis revealed that this pathway was closely related to hypoxia and acidification stress in the five groups, and we found that genes of the collagen family may play a key role in this pathway. The zf-C2H2 transcription factor may play an important role in the hypoxia and acidification stress response, and novel genes were additionally identified. The results provide new clues for further research on the molecular mechanisms underlying hypoxia–acidification tolerance in L. crocea and provides a basic understanding of the potential combined effects of reduced pH and dissolved oxygen on Sciaenidae fishes.

Linked selection, differential introgression and recombination rate variation promote heterogeneous divergence in a pair of yellow croakers

Understanding the mechanisms underlying heterogeneous genomic divergence is of particular interest in evolutionary biology. Highly differentiated genomic regions, known as genomic islands, often evolve between diverging lineages. These genomic islands may be related to selection promoting adaptation or reproductive isolation. Based on whole genome assembly and genome-wide RAD sequencing in a pair of yellow croakers (genus: Larimichthys), we investigated the evolutionary processes shaping genomic landscapes of divergence. Demographic modelling indicated that the two species diverged following a secondary contact scenario, where differential introgression and linked selection were suggested to be involved in heterogeneous genomic divergence. We identified reduced recombination rate in genomic islands and a relatively good conservation of both genetic diversity and recombination landscapes between species, which highlight the roles of linked selection and recombination rate variation in promoting heterogeneous divergence in the common ancestral lineage of the two species. In addition, we found a positive correlation between differentiation (F_ST ) and absolute sequence divergence (D_xy ), and elevated D_xy in genomic islands, indicating that the genomic landscape of divergence was not shaped by linked selection alone. Restricted gene flow in highly differentiated regions has probably remodelled the landscape of heterogeneous genomic divergence. This study highlights that highly differentiated genomic regions can also arise from a combination of linked selection and differential gene flow in interaction with varying recombination rates.

Extensive involvement of CD40 and CD154 costimulators in multiple T cell-mediated responses in a perciform fish Larimichthys crocea

CD40 and CD154 are well-characterized costimulatory molecules involved in adaptive humoral immunity in humans and other mammals. These two costimulatory molecules were found to be originated from teleost fish during vertebrate evolution. However, the functionality of fish CD40 and CD154 remains to be explored. In this study, we identified the CD40 and CD154 homologs (LcCD40 and LcCD154) from large yellow croaker (Larimichthys crocea), a marine species of the perciform fish family. The LcCD40 and LcCD154 share conserved structural features to their mammalian counterparts, and are widely expressed in immune-relevant tissues and leukocytes at different transcriptional levels. Immunofluorescence staining and FCM analysis showed that LcCD40 and LcCD154 proteins are distributed on MHC-II⁺ APCs and CD4-2⁺ T cells, and are significantly upregulated in response to antigen stimulation. Co-IP assay exhibited strong association between LcCD40 and LcCD154 proteins. Blockade of LcCD154 with anti-LcCD154 antibody (Ab) or recombinant soluble LcCD40-Ig fusion protein remarkably decreased the MHC-II⁺ APC-initiated CD4⁺ T cell response upon Aeromonas hydrophila stimulation, and alloreactive T cell activation as examined by mixed lymphocyte reaction (MLR). These findings highlight the costimulatory role of LcCD40 and LcCD154 in T cell activities in Larimichthys crocea. Thus, the CD40 and CD154 costimulators may extensively participate in the regulation of multiple T cell-mediated immune responses in teleost fish. It is anticipated that this study would provide a cross-species understanding of the evolutionary history of CD40 and CD154 costimulatory signals from fish to mammals.

Molecular and Structural Basis of Receptor Binding and Signaling of a Fish Type I IFN with Three Disulfide Bonds

In mammals, type I IFNs, which commonly contain one or two disulfide bonds, activate the JAK-STAT signaling pathway through binding to the common cell surface receptor formed by IFN-α/β receptor (IFNAR)1 and IFNAR2 subunits. Although type I IFNs are also known to be essential for antiviral defense in teleost fish, very little is known about mechanisms underlying the recognition of fish type I IFNs by associated receptors. In this study, we demonstrate that a type I IFN of large yellow croaker Larimichthys crocea (LcIFNi), belonging to a new subgroup of fish type I IFNs, triggers antiviral response via the conserved JAK-STAT pathway through stable binding with a heterodimeric receptor comprising subunits LcCRFB5 and LcCRFB2. LcIFNi binds to LcCRFB5 with a much higher affinity than to LcCRFB2. Furthermore, we determined the crystal structure of LcIFNi at a 1.39 Å resolution. The high-resolution structure is, to our knowledge, the first reported structure of a type I IFN with three disulfide bonds, all of which were found to be indispensable for folding and stability of LcIFNi. Using structural analysis, mutagenesis, and biochemical assays, we identified key LcIFNi residues involved in receptor interaction and proposed a structural model of LcIFNi bound to the LcCRFB2–LcCRFB5 receptor. The results show that LcIFNi–LcCRFB2 exhibits a similar binding pattern to human IFN-ω–IFNAR2, whereas the binding pattern of LcIFNi–LcCRFB5 is quite different from that of IFN-ω–IFNAR1. Altogether, our findings reveal the structural basis for receptor interaction and signaling of a type I IFN with three disulfide bonds and provide new insights into the mechanisms underlying type I IFN recognition in teleosts.

Disruption of mstn Gene by CRISPR/Cas9 in Large Yellow Croaker (Larimichthys crocea)

The large yellow croaker (Larimichthys crocea) plays an economically vital role in the marine aquaculture in China. Suffering from infection of bacteria and protozoon, effect of extreme weather and stress from high-density farming, genome editing is thought to be an important tool applied to L. croea for enhancing commercial traits such as growth rate, disease resistance, and nutrition component. In this study, we identified two mstn genes in L. croea and investigated the different phylogenetic clades, gene structures, and conserved syntenic relationships. To obtain fast-growing large yellow croaker, we specially selected two validated targets for mstnb knockout, which was homologous to mammalian myostatin gene (MSTN) and downregulated skeletal muscle growth and development. Five significant mutation types were generated in two mosaic mutants by transferring specific CRISPR/Cas9 RNPs (ribonucleoprotein) into the one-cell fertilized embryos based on CRISPR/Cas9 technology. Subsequently, we also elucidated the obstacles and possible measures to improve the success rate of inducing modified large yellow croaker. Our results would provide valuable method and reference for facilitating genome editing programs of the large yellow croaker in the future.

High Genetic Connectivity Inferred from Whole-Genome Resequencing Provides Insight into the Phylogeographic Pattern of Larimichthys polyactis

Compared with terrestrial biota, marine fishes usually present lower genetic differentiation among different geographical populations because of high-level gene flow and lack of physical barriers. Understanding the genetic structure of marine fishes is essential for dividing management unit and making reasonable protection measures. The small yellow croaker (Larimichthys polyactis) belongs to the family Sciaenidae, which is an economic fish and widely distributed in the Western Pacific. To delineate genetic diversity and phylogeographic pattern, whole-genome resequencing was used to evaluate genetic connectivity, genetic diversity, and spatial pattern of L. polyactis for the first time. We obtained 6,645,711 high-quality single nucleotide polymorphisms (SNPs) markers from 40 L. polyactis individuals. The phylogenetic analysis, STRUCTURE, principal component analysis, and F_st results all indicated that no genetic structure consistent with the distribution pattern was found. This result revealed high genetic connectivity of L. polyactis in different sampling sites. High genetic diversity was also detected, indicating that there was sufficient evolutionary potential to maintain its effective population size. Besides, a similar result of high genetic connectivity and genetic diversity was also detected by mitochondria DNA marker. Our study demonstrated the persistence of high levels of genetic connectivity and a lack of population structure across L. polyactis in different sea areas. This study aimed to analyze the division of population structure and the reason for the decline and not exhaustion of L. polyactis resource on a genetic level.

A chromosome-level genome assembly of the jade perch (Scortum barcoo)

Endemic to Australia, jade perch (Scortum barcoo) is a highly profitable freshwater bass species. It has extraordinarily high levels of omega-3 polyunsaturated fatty acids (PUFAs), which detailed genes involved in are largely unclear. Meanwhile, there were four chromosome-level bass species have been previous sequenced, while the bass ancestor genome karyotypes have not been estimated. Therefore, we sequenced, assembled and annotated a genome of jade perch to characterize the detailed genes for biosynthesis of omega-3 PUFAs and to deduce the bass ancestor genome karyotypes. We constructed a chromosome-level genome assembly with 24 pairs of chromosomes, 657.7 Mb in total length, and the contig and the scaffold N50 of 4.8 Mb and 28.6 Mb respectively. We also identified repetitive elements (accounting for 19.7% of the genome assembly) and predicted 26,905 protein-coding genes. Meanwhile, we performed genome-wide localization and characterization of several important genes encoding some key enzymes in the biosynthesis pathway of PUFAs. These genes may contribute to the high concentration of omega-3 in jade perch. Moreover, we conducted a series of comparative genomic analyses among four representative bass species at a chromosome level, resulting in a series of sequences of a deductive bass ancestor genome.

Measurement(s)	whole genome sequencing
Technology Type(s)	Illumina Sequencing • Oxford Nanopore Sequencing

The chromosome-level genome assembly of the Japanese yellowtail jack Seriola aureovittata provides insights into genome evolution and efficient oxygen transport

Fishes of the genus Seriola are widely farmed and highly valued in global aquaculture production. To further understand their economically important traits and help improve aquaculture product quality and sustainability, we performed a chromosome-level genome construction for Seriola aureovittata. Combining two technologies, PacBio and BGISEQ-500, we assembled 649.86 Mb S. aureovittata genome sequences with a contig N50 of 22.21 Mb, and 98% of BUSCO genes were detected in total. The initial assembly was then further scaffolded into 24 pseudochromosomes using Hi-C data, indicating the high quality of the genome. Genome evolution analysis showed that many genes related to fatty acid metabolism and oxygen binding, or transport were expanded, which provided insights into the metabolic characteristics of fatty acids and efficient oxygen transport. Based on the genome data, we confirmed the evolutionary relationship of S. aureovittata, S. dorsalis and S. lalandi and identified chr12 as the putative sex chromosome of S. aureovittata. Our chromosome-level genome assembly provides a genetic foundation for the phylogenetic and taxonomic investigation of different Seriola species. Moreover, the genome will provide an important genomic resource for further biological and aquaculture studies of S. aureovittata.

IL-2 Signaling Couples the MAPK and mTORC1 Axes to Promote T Cell Proliferation and Differentiation in Teleosts

IL-2 is a pleiotropic cytokine that is critical for T cell immunity. Although the IL-2-mediated regulation of T cell immunity in mammals is relatively well understood, it remains largely unknown whether and how IL-2 regulates T cell immunity in lower vertebrates. To address this knowledge gap, we investigated the role played by IL-2 in the regulation of T cell response, as well as the associated underlying mechanisms in a teleost fish, large yellow croaker (Larimichthys crocea). We found that large yellow croaker (L. crocea) IL-2 (LcIL-2) significantly promoted T cell proliferation both in vivo and in vitro; significantly induced the differentiation of Th1, Th2, regulatory T, and cytotoxic T cells while inhibiting Th17 differentiation; and participated in the elimination of invading pathogenic bacteria. Mechanistically, the binding of LcIL-2 to its heterotrimer receptor complex (LcIL-15Rα/LcIL-2Rβ/Lcγc) triggered the conserved JAK-STAT5 pathway, which in turn regulated the expression of genes involved in T cell expansion, differentiation, and biological function. The MAPK and mammalian target of rapamycin complex 1 (mTORC1) axes, which are involved in TCR-mediated signaling, were also required for LcIL-2-mediated T cell response. Collectively, our results demonstrated that fish IL-2 plays a comprehensive regulatory role in T cell response and highlighted the complex and delicate network regulating T cell-driven immune response. We propose that T cell immunity is regulated by the interplay between TCR signaling and cytokine signaling, and that this basic strategy evolved before the emergence of the tetrapod lineage. Our findings provide valuable insights into the regulatory mechanisms underlying T cell response in teleosts.

PU.1 Regulates Cathepsin S Expression in Large Yellow Croaker (Larimichthys crocea) Macrophages

Different morphologies have been detected in teleost macrophages. In this study, two macrophage cell lines were sub-cloned from a large yellow croaker head kidney cell line, LYCK. One type of sub-cloned cells was fusiform but the other was round, named LYC-FM and LYC-RM cells respectively, based on their morphologies. Both types showed the characteristics of macrophages, including expression of macrophage-specific marker genes, possession of phagocytic and bactericidal activities, and production of reactive oxygen species (ROS) and nitric oxide (NO). The transcription factor PU.1, crucial for the development of macrophages in mammals, was found to exist in two transcripts, PU.1a and PU.1b, in large yellow croaker, and constitutively expressed in LYC-FM and LYC-RM cells. The expression levels of PU.1a and PU.1b could be upregulated by recombinant large yellow croaker IFN-γ protein (rLcIFN-γ). Further studies showed that both PU.1a and PU.1b increased the expression of cathepsin S (CTSS) by binding to different E26−transformation−specific (Ets) motifs of the CTSS promoter. Additionally, we demonstrated that all three domains of PU.1a and PU.1b were essential for initiating CTSS expression by truncated mutation experiments. Our results therefore provide the first evidence that teleost PU.1 has a role in regulating the expression of CTSS.

Genomic deciphering of sex determination and unique immune system of a potential model species rare minnow (Gobiocypris rarus)

Gobiocypris rarus is sensitive to environmental pollution, especially to heavy metal and grass carp reovirus (GCRV). Hence, it has potential utility as a biological monitor. Genetic deciphering of its unique immune system will advance our understanding of its unique adaptive strategies, which provide cues for its better application. A de novo genome of rare minnow was obtained, and its sex determination mechanism is ZZ/ZW. We identified several specific mutation genes and specific lost genes of rare minnow, and these might be related to the sensitivity of rare minnow to environmental stimuli. We also analyzed the gene expression level of different organs/tissues and found that several IFIT genes may play key roles in GCRV resistance. In addition, knockout of the gene PCDH10L indicates that PCDH10L affects Pb²⁺-induced mortality in rare minnow. Rare minnow is ready for genetic manipulation and shows potential as an emerging experimental model.

Identification and characterization of mkk genes and their expression profiles in rainbow trout (Oncorhynchus mykiss) symptomatically or asymptomatically infected with Vibrio anguillarum

Mitogen-activated protein kinase kinases (MKKs) are intermediate kinases of mitogen-activated protein kinases (MAPKs) signaling pathways. MKKs are activated by mitogen-activated protein kinase kinase kinase (MKKK) and then the activated MKKs trigger the activation of downstream MAPKs. MAPK signaling pathways play an important role in regulating immune functions including apoptosis and inflammation. However, studies on identification and characterization of mkk repertoire in rainbow trout (Oncorhynchus mykiss) are still limited. Trout experienced 4 rounds (4R) of whole genome duplication (WGD), thus exhibiting increased paralogs of mkks with potentially functional diversity. In this study, we identified 17 mkk genes in trout and the following bacterial challenge (Vibrio anguillarum) studies showed functional diversity of different mkk subtypes. Vibrio anguillarum infection resulted in significantly up-regulated mkk2 subtypes in spleen and liver, and mkk4b3 in spleen, suggesting immunomodulation was regulated by activation of ERK, p38 and JNK pathways. Compared to other mkk subtypes, mkk6s were down-regulated in symptomatic group, rather than asymptomatic group. The organisms present negative feedback on MAPK activation, thus reducing extra damage to cells. We observed down-regulated mkk6s with up-regulated genes (dusp1 & dusp2) involved in negative feedback of MAPK activation. Based on these results, we might propose the distinct expression patterns of genes associated with MAPK pathways resulted in different phenotypes and symptoms of trout in response to bacterial challenge.

Identification and Characterization of Immunoglobulin T Heavy Chain in Large Yellow Croaker (Larimichthys crocea)

Three immunoglobulin (Ig) isotypes have been identified in teleosts, IgM, IgD, and IgT or IgZ. IgT, a new teleost Ig isotype, plays a vital role in mucosal immunity. However, information on molecular and functional characteristics of fish IgT is still limited. In this study, an IgT heavy chain (LcIgT) gene was cloned and characterized in large yellow croaker (Larimichthys crocea). Complete cDNA of LcIgT was 1930 bp in length, encoding a protein of 554 amino acids. The deduced LcIgT contains a VH region and only three CH regions (CH1, CH2, CH4), but no transmembrane region was predicted. Phylogenetic analysis showed that IgT heavy chain sequences from all fish species are grouped together. Homology comparison showed that LcIgT shares the highest amino acid identity of 58.73% with IgT heavy chain in Scophthalmus maximus. The VH domain of LcIgT has the highest identity of 72.50% with that of Scophthalmus maximus IgT. Relatively, each constant domain of LcIgT exhibits the highest amino acid identity with that of IgT in Oreochromis niloticus (67.61% identity for CH1, 61.11% identity for CH2, and 63.74% identity for CH4). LcIgT was constitutively expressed in various tissues tested, with the highest levels in mucosa-associated tissues such as gills and skin. After Cryptocaryon irritans infection, the mRNA levels of LcIgT were significantly up-regulated in the spleen (3.27-fold) at 4 d, in the head kidney (3.98-fold) and skin (2.11-fold) at 7 d, and in gills (4.45-fold) at 14 d. The protein levels in these detected tissues were all significantly up-regulated; the peak of its up-regulation was 6.33-fold at 28d in gills, 3.44-fold at 7d in skin, and 3.72-fold at 14d in spleen. These results showed that IgT response could be simultaneously induced in both systemic and mucosal tissues after parasitic infection and that IgT may be involved in systemic immunity and mucosal immunity against parasitic infection.

Development of Recombinant Dihydrolipoamide Dehydrogenase Subunit Vaccine against Vibrio Infection in Large Yellow Croaker

Large yellow croaker (Larimichthys crocea), an economically important marine fish in China, has suffered from serious vibriosis, which has resulted in great economic losses for the large yellow croaker industry. Vaccination has been considered to be a safe and effective method to prevent and control vibriosis. However, due to the complex diversity and serotypes of the Vibrio genus, the progress of Vibrio vaccine development is still slow. In this study, we prepared recombinant Vibrio dihydrolipoamide dehydrogenase (rDLD) protein and investigated its potential as a candidate to be a subunit vaccine against Vibrio. The lysozyme activity and the rDLD-specific antibody level in sera of large yellow croakers immunized with rDLD were significantly higher than those in the control group, and the transcript levels of proinflammatory cytokines (IL-6, IL-8, IL-1β), MHC IIα/β, CD40, CD8α, IL-4/13A, and IL-4/13B were significantly up-regulated in the spleen and head kidney of large yellow croakers immunized with rDLD, suggesting that rDLD could induce both specific and nonspecific immune responses in this species. In addition, rDLD protein increased the survival rate of large yellow croakers against Vibrio alginolyticus and Vibrio parahaemolyticus, with the relative percent of survival (RPS) being 74.5% and 66.9%, respectively. These results will facilitate the development of a potential subunit vaccine against Vibrio in large yellow croaker aquaculture.

A high-quality chromosome-level genome of the endangered roughskin sculpin provides insights into its evolution and adaptation

The cottids (Cottidae) are a taxonomically diverse and ecologically important component of many marine and freshwater ecosystems. Despite recent breakthroughs in long-read sequencing, high quality genomic resources are still limited for studies of ecological and evolutionary processes in cottids. Here we generated a high-quality, chromosome-scale genome assembly (521.26 Mb) of the catadromous roughskin sculpin (Trachidermus fasciatus Heckel) with a contig N50 of 2.93 Mb and a scaffold N50 of 24.06 Mb. Approximately 21.97% of the genome was composed of repetitive elements. A total of 21,872 protein-coding genes were predicted, of which 19,900 genes (90.98%) were functionally annotated. Phylogenetic analysis supported the validity of Scorpaenoidei and Cottioidei as two suborders of the Perciformes. Chromosome-scale collinearity analyses identified four chromosome fusions leading to the reduction of chromosome number in T. fasciatus. Gene families related to cell apoptosis and cell death were expanded and those related to immune system were contracted, suggesting that these gene families might be relevant to a host of phenotypic differences between T. fasciatus and other teleosts. Gene families associated with osmoregulation were also expanded, which might be associated with its catadromous life history. A total of 50 aging-associated genes were found to be under positive selection, which might be associated with the short lifespan of T. fasciatus. The high-quality genome assembly and annotation will promote researches into the evolution of catadromous life history and short lifespan for T. fasciatus and facilitate comparative genomic studies of cottids in the near future.

Molecular characterization, expression analysis and cellular location of IL-4/13 receptors in large yellow croaker (Larimichthys crocea)

Interleukin (IL)-4 and IL-13 are closely related class I cytokines that play key roles in the T helper (Th)-2 immune response via heterodimeric receptors. IL-4 signals via both the type I (IL-4Rα/γc) and type II (IL-4Rα/IL-13Rα1) receptor complexes, while IL-13 signals only via the type II receptor complex. IL-13Rα2 is traditionally considered a "decoy" receptor for IL-13. However, the IL-4/13 system and its response to pathogenic infection are still not fully understood in fish. In this study, we identified four IL-4/13 receptor subunit genes in the large yellow croaker (Larimichthys crocea): LcIL-4Rα1, LcIL-4Rα2, LcIL-13Rα1, and LcIL-13Rα2. Sequence analysis showed that these receptors possessed typical characteristic domains, including a signal peptide, two fibronectin type III (FN III)-like domains, and a transmembrane domain, but their cytoplasmic regions were not well conserved. The mRNA and protein of the four IL-4/13 receptors were constitutively expressed in all examined tissues of large yellow croaker. Their mRNAs were also detected in primary head kidney macrophages (PKMs), primary head kidney granulocytes (PKGs), and primary head kidney lymphocytes (PKLs). Immunofluorescence assay further showed that LcIL-4Rα and LcIL-13Rα1 were expressed on the membrane of IgM ⁺ B cells. After stimulation by Vibrio alginolyticus and poly (I:C) (a viral dsRNA mimic), the mRNA levels of LcIL-4/13 receptors were significantly upregulated in the head kidney and spleen. Their mRNA levels were also upregulated in head kidney leukocytes in response to poly (I:C) and lipopolysaccharide (LPS) treatment. Moreover, both recombinant LcIL-4/13A and LcIL-4/13B upregulated LcIL-4Rα1 and LcIL-4Rα2 in primary leukocytes, but only recombinant LcIL-4/13A upregulated LcIL-13Rα1 and LcIL-13Rα2. These results indicated that LcIL-4/13 receptors, containing conserved functional domains, may be involved in the IL-4/13-mediated immune response to pathogenic infections in the large yellow croaker.

HSP90 and HSP70 Families in Lateolabrax maculatus: Genome-Wide Identification, Molecular Characterization, and Expression Profiles in Response to Various Environmental Stressors

Heat shock proteins (HSPs) are a large class of highly conserved chaperons, which play important roles in response to elevated temperature and other environmental stressors. In the present study, 5 HSP90 genes and 17 HSP70 genes were systematically characterized in spotted seabass (Lateolabrax maculatus). The evolutionary footprint of HSP genes was revealed via the analysis of phylogeny, chromosome location, and gene copy numbers. In addition, the gene structure features and the putative distribution of heat shock elements (HSEs) and hypoxia response elements (HREs) in the promoter regions were analyzed. The protein-protein interaction (PPI) network analyses results indicated the potential transcriptional regulation between the heat shock factor 1 (HSF1) and HSPs and a wide range of interactions among HSPs. Furthermore, quantitative (q)PCR was performed to detect the expression profiles of HSP90 and HSP70 genes in gill, liver, and muscle tissues after heat stress, meanwhile, the expression patterns in gills under alkalinity and hypoxia stresses were determined by analyzing RNA-Seq datasets. Results showed that after heat stress, most of the examined HSP genes were significantly upregulated in a tissue-specific and time-dependent manners, and hsp90aa1.1, hsp90aa1.2, hsp70.1, and hsp70.2 were the most intense responsive genes in all three tissues. In response to alkalinity stress, 11 out of 13 significantly regulated HSP genes exhibited suppressed expression patterns. Alternatively, among the 12 hypoxia-responsive-expressed HSP genes, 7 genes showed induced expressions, while hsp90aa1.2, hsp70.1, and hsp70.2 had more significant upregulated changes after hypoxic challenge. Our findings provide the essential basis for further functional studies of HSP genes in response to abiotic stresses in spotted seabass.

Recurrent expansions of B30.2-associated immune receptor families in fish

B30.2 domains, also known as PRY/SPRY, are key components of specific subsets of two large families of proteins involved in innate immunity: the tripartite motif proteins (TRIMs) and the Nod-like receptors (NLRs). TRIM proteins are important, often inducible factors of antiviral innate immunity, targeting multiple steps of viral cycles through a variety of mechanisms. NLRs prime and regulate systemic innate defenses, especially against bacteria, and control inflammation. Large TRIM and NLR subsets characterized by the presence of a B30.2 domain have been reported from a few fish species including zebrafish and seem to be strongly prone to gene duplication/expansion. Here, we performed a large-scale survey of these receptors across about 150 fish genomes, focusing on ray-finned fishes. We assessed the number and genomic distribution of domains and domain combinations associated with TRIMs, NLRs, and other genes containing B30.2 domains and looked for gene expansion patterns across fish groups. We then used a model to test the impact of taxonomy, genome size, and environmental variables on the copy numbers of these genes. Our findings reveal novel domain structures, clade-specific gains and losses. They also assist with the timing of the gene expansions, reveal patterns associated with the MHC, and lay the groundwork for further studies delving deeper into the forces that drive the copy number variation of immune genes on a species level.

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.

Whole-Genome Sequencing and Genome-Wide Studies of Spiny Head Croaker (Collichthys lucidus) Reveals Potential Insights for Well-Developed Otoliths in the Family Sciaenidae

Spiny head croaker (Collichthys lucidus), belonging to the family Sciaenidae, is a small economic fish with a main distribution in the coastal waters of Northwestern Pacific. Here, we constructed a nonredundant chromosome-level genome assembly of spiny head croaker and also made genome-wide investigations on genome evolution and gene families related to otolith development. A primary genome assembly of 811.23 Mb, with a contig N50 of 74.92 kb, was generated by a combination of 49.12-Gb Illumina clean reads and 35.24 Gb of PacBio long reads. Contigs of this draft assembly were further anchored into chromosomes by integration with additional 185.33-Gb Hi-C data, resulting in a high-quality chromosome-level genome assembly of 817.24 Mb, with an improved scaffold N50 of 26.58 Mb. Based on our phylogenetic analysis, we observed that C. lucidus is much closer to Larimichthys crocea than Miichthys miiuy. We also predicted that many gene families were significantly expanded (p-value <0.05) in spiny head croaker; among them, some are associated with "calcium signaling pathway" and potential "inner ear functions." In addition, we identified some otolith-related genes (such as otol1a that encodes Otolin-1a) with critical deletions or mutations, suggesting possible molecular mechanisms for well-developed otoliths in the family Sciaenidae.

Chromosome-scale assembly and high-density genetic map of the yellow drum, Nibea albiflora

The yellow drum (Nibea albiflora) is an economically important sciaenid fish in East Asian countries. In this study, we sequenced and assembled a near-complete gynogenetic yellow drum genome. We generated 45.63 Gb of Illumina short-reads and 80.27 Gb of PacBio long-reads and assembled them into a 628.01-Mb genome with a contig N50 of 4.42 Mb. Twenty-four chromosomes with a scaffold N50 of 26.73 Mb were obtained using the Hi-C analysis. We predicted a set of 27,069 protein-coding genes, of which 1,581 and 2,583 were expanded and contracted gene families, respectively. The most expanded genes were categorised into the protein binding, zinc-ion binding and ATP binding functional pathways. We built a high-density genetic linkage map that spanned 4,300.2 cM with 24 linkage groups and a resolution of 0.69 cM. The high-quality reference genome and annotated profiles that we produced will not only increase our understanding of the genetic architecture of economic traits in the yellow drum, but also help us explore the evolution and unique biological characteristics of sciaenid fishes.

Molecular characterization and role in virus infection of Beclin-1 in large yellow croaker (Larimichthys crocea)

Beclin-1, the ortholog of yeast autophagy-related gene 6 (Atg6), has a central role in autophagy, which has been linked to diverse biological processes including immunity, development, tumor suppression, and lifespan extension. However, understanding of function of fish Beclin-1 is limited now. In this study, the complete Beclin-1 cDNA of large yellow croaker Larimichthys crocea (LcBeclin-1) was cloned, whose open reading frame (ORF) is 1344 bp long and encodes a protein of 447 amino acids (aa). The deduced LcBeclin-1 possesses a typical Bcl-2 homology domain 3(BH3) and an APG6 domain that contains a central coiled-coil domain (CCD, residues 174 to 231) and a C-terminal evolutionarily conserved domain (ECD, residues 241 to 334). LcBeclin-1 shared a high amino acid identity of 81.66-98.66% with reported Beclin-1 molecules from other vertebrate species. LcBeclin-1 gene was constitutively expressed in all tissues tested, with the highest levels in heart. LcBeclin-1 transcripts were also detected in primary head kidney granulocytes (PKGs), primary head kidney macrophages (PKMs), primary head kidney leukocytes (PKLs), and large yellow croaker head kidney cell line (LYCK), and were significantly upregulated by poly (I:C) in PKMs and LYCK cells. Subcellular localization showed that LcBeclin-1 was evenly distributed in the cytoplasm and nucleus of LYCK cells. Overexpression of LcBeclin-1 significantly increased the replication of SVCV, as evidenced by increased severity of the cytopathic effects, enhanced viral titre, and upregulated transcriptional levels of viral genes. Further studies showed that LcBeclin-1 induced the occurrence of autophagy in LYCK cells. Additionally, LcBeclin-1 also decreased the expression levels of large yellow croaker interferons (IFNs; IFNc, IFNd, and IFNh), interferon regulatory factor 3 (IRF3) and IRF7, IFN-stimulated genes (ISGs; Mx, PKR, and Viperin) in LYCK cells. All these data suggest that LcBeclin-1 promoted the viral replication possibly by inducing autophagy or negatively modulating IFN response, which will help us to further understand the function of fish Beclin-1.

Hypoxic acclimation improves mitochondrial bioenergetic function in large yellow croaker Larimichthys crocea under Cu stress

The purpose of this study was to investigate how pre-hypoxia exposure affected the mitochondrial structure and bioenergetic function of large yellow croaker in responding to Cu stress. Fish were acclimated to normoxia and 3.0 mg DO L-1 for 48 h, then subjected to 0 and 120 μg Cu L-1 for another 48 h. Hypoxic acclimation did not affect mitochondrial ultrastructure and reactive oxygen species (ROS), but reduced oxidative phosphorylation (OXPHOS) efficiency. Cu exposure impaired mitochondrial ultrastructure, increased ROS generation and inhibited OXPHOS efficiency. Compared with Cu exposure alone, hypoxic acclimation plus Cu exposure reduced ROS production and improved OXPHOS efficiency by enhancing mitochondrial respiratory control ratio, mitochondrial membrane potential, and activities and gene expressions of electron transport chain enzymes. In conclusion, hypoxic acclimation improved the mitochondrial energy metabolism of large yellow croaker under Cu stress, facilitating our understanding of the molecular mechanisms regarding adaptive responses of hypoxia-acclimated fish under Cu stress.

The application of genome editing technology in fish

The advent and development of genome editing technology has opened up the possibility of directly targeting and modifying genomic sequences in the field of life sciences with rapid developments occurring in the last decade. As a powerful tool to decipher genome data at the molecular biology level, genome editing technology has made important contributions to elucidating many biological problems. Currently, the three most widely used genome editing technologies include: zinc finger nucleases (ZFN), transcription activator like effector nucleases (TALEN), and clustered regularly interspaced short palindromic repeats (CRISPR). Researchers are still striving to create simpler, more efficient, and accurate techniques, such as engineered base editors and new CRISPR/Cas systems, to improve editing efficiency and reduce off-target rate, as well as a near-PAMless SpCas9 variants to expand the scope of genome editing. As one of the important animal protein sources, fish has significant economic value in aquaculture. In addition, fish is indispensable for research as it serves as the evolutionary link between invertebrates and higher vertebrates. Consequently, genome editing technologies were applied extensively in various fish species for basic functional studies as well as applied research in aquaculture. In this review, we focus on the application of genome editing technologies in fish species detailing growth, gender, and pigmentation traits. In addition, we have focused on the construction of a zebrafish (Danio rerio) disease model and high-throughput screening of functional genes. Finally, we provide some of the future perspectives of this technology.

Cross-species Exon Capture and Whole Exome Sequencing: Application, Utility and Challenges for Genomic Resource Development in Non-model Species

Comprehending the genetic architecture of complex traits has many applications in evolution, ecology, conservation biology and plant and animal production systems. Underlying research questions in these fields are diverse species that often have limited genetic information available. In aquaculture, for example, genetic progress has been slow in many species due to a lack in such genetic information. In this study, zebrafish (as a well-studied model species) was used in cross-species transfer to develop genomic resources and identify candidate genes underling growth differentials in dusky kob. Dusky kob is a Sciaenid finfish and an emerging aquaculture species. The zebrafish All Exon Predesigned Probe-set capture protocol was used to enrich fractionated DNA samples from kob, classified as either large or small, before massive parallel sequencing on the Ion Torrent platform. Although vast quantities of sequence data were generated, only about 30% of contigs could be identified as zebrafish homologues. There were numerous species-specific sequences and inconsistent coverage of sequencing products across samples, likely due to non-specific binding of the probe-set as a result of the evolutionary divergence between zebrafish and kob. Nonetheless, more than 55,000 SNPs could be reliably identified and genotyped to the individual level. Using SNP genotypic divergence estimates, between large and small cohorts, a number of candidate genes associated with growth was also identified for future investigation. These findings contribute to the growing body of evidence demonstrating the utility of a cross-species capture approach in the development of important genomic resources for understanding traits of interest in species without reference genomes.

Hypoxic responses in Oncorhynchus mykiss involve angiogenesis, lipid, and lactate metabolism, which may be triggered by the cortisol stress response and epigenetic methylation

The incidence of hypoxia in water bodies is increasing more rapidly than aquatic life can adapt. This study aimed to determine the effects of hypoxia on fish physiology, as well as protein expression through proteomics. To do this, 40 rainbow trout were divided into normoxic control (11.5 mg/L dissolved oxygen) and hypoxic treatment (5 mg/L dissolved oxygen) tanks for a period of 7 days. Fish were then anesthetized and blood was sampled. Fish were then euthanized and heart and liver samples were taken. Blood glucose, cortisol and lipid, body and liver mass, fork length, hematocrit and, blood cell counts and global heart methylation were measured. Red blood cell counts were significantly lower, while hematocrit and mean corpuscular volume were significantly higher in the hypoxic treatment. Global DNA methylation was significantly decreased in hypoxic heart tissue. Plasma cortisol and 18:1 monoacylglyerol increased, while 15:0-18:1 phosphatidylethanolamine, and 18:1 lysophosphatidylethanolamine decreased in plasma of rainbow trout under hypoxic conditions. Plasma proteomics revealed 70 significantly altered proteins (p < 0.05) in the hypoxia treatment (Data are available via ProteomeXchange with identifier PXD026589). Many of these molecular changes appear to be related to the observed increase in red blood cell volume and epigenetic modifications, as well as to angiogenesis, lipid, and glucose metabolism. This study highlights a range of cellular and molecular responses in the blood and plasma of freshwater fish that may be phenotypic adaptions to hypoxia, and that could aid in diagnosing the health status of wild fish populations using several, potential, discovered biomarkers.

Development of monoclonal antibody against IgT of a perciform fish, large yellow croaker (Larimichthys crocea) and characterization of IgT+ B cells

Teleost immunoglobulin T (IgT) is considered to be a primitive immunoglobulin class specialized in mucosal immunity. In the present study, a recombinant protein containing the CH2 region of large yellow croaker (Larimichthys crocea) IgT heavy chain was expressed, purified, and used as an immunogen to produce a monoclonal antibody (mAb) against large yellow croaker IgT. Western blotting results indicated that the obtained mouse anti-IgT mAb could specifically recognize a 45 kDa protein in the skin mucus of large yellow croaker, which was identified as the IgT heavy chain by mass spectrometric analysis. Immunofluorescence assay (IFA) analysis further demonstrated that this mouse anti-IgT mAb could recognize membrane-bound IgT (mIgT) molecules on large yellow croaker IgT⁺ leukocytes. This mAb also could be used for sorting of large yellow croaker IgT⁺ B cells by flow cytometry sorting technology. Then, flow cytometric immunofluorescence analysis (FCIA) results showed that the percentages of IgT⁺ B cells in skin, gills, gut, spleen, head kidney and peripheral blood lymphocytes were 27.553% ± 3.312%, 12.588% ± 3.538%, 12.355% ± 3.352%, 13.075 ± 2.258%, 5.552 ± 3.275%, and 2.600 ± 0.521%, respectively, indicating that mucosal tissues (skin, gills, and gut) contained a high ratio of IgT⁺ B cells. Accordingly, the high protein levels of IgT were also detected in these mucosal tissues, suggesting that IgT may play a role in mucosal immunity in large yellow croaker. Taken together, our data demonstrated that the mouse anti-IgT mAb developed in this study could be used for characterizing IgT⁺ B cells and studying the functions of IgT in large yellow croaker.

Molecular cloning, inducible expression with SGIV and Vibrio alginolyticus challenge, and function analysis of Epinephelus coioides PDCD4

Programmed cell death 4 (PDCD4) in mammals, a gene closely associated with apoptosis, is involved in many biological processes, such as cell aging, differentiation, regulation of cell cycle, and inflammatory response. In this study, grouper Epinephelus coioides PDCD4, EcPDCD4-1 and EcPDCD4-2, were obtained. The open reading frame (ORF) of EcPDCD4-1 is 1413 bp encoding 470 amino acids with a molecular mass of 52.39 kDa and a theoretical pI of 5.33. The ORF of EcPDCD4-2 is 1410 bp encoding 469 amino acids with a molecular mass of 52.29 kDa and a theoretical pI of 5.29. Both EcPDCD4-1 and EcPDCD4-2 proteins contain two conserved MA3 domains, and their mRNA were detected in all eight tissues of E. coioides by quantitative real-time PCR (qRT-PCR) with the highest expression in liver. The expressions of two EcPDCD4s were significantly up-regulated after Singapore grouper iridovirus (SGIV) or Vibrio alginolyticus infection. In addition, over-expression of EcPDCD4-1 or EcPDCD4-2 can inhibit the activity of the nuclear factor-κB (NF-κB) and activator protein-1 (AP-1), and regulate SGIV-induced apoptosis. The results demonstrated that EcPDCD4s might play important roles in E. coioides tissues during pathogen-caused inflammation.

Whole-genome resequencing of large yellow croaker (Larimichthys crocea) reveals the population structure and signatures of environmental adaptation

Large yellow croaker is an economically important fish in China and East Asia. Despite its economic importance, genome-wide adaptions of domesticated large yellow croaker are largely unknown. Here, we performed whole-genome resequencing of 198 individuals of large yellow croaker obtained in the sea or from farmers in Zhoushan or Ningde. Population genomics analyses revealed the genetic population structure of our samples, reflecting the living environment. Each effective population size is estimated to be declining over generations. Moreover, we identified genetically differentiated genomic regions between the sea-captured population in the Zhoushan Sea area and that of the Ningde Sea area or between the sea-captured population and the farmed population in either area. Gene ontology analyses revealed the gene groups under selective sweep for the adaptation to the domesticated environment. All these results suggest that individuals of the large yellow croaker populations show genomic signatures of adaptation to different living environments.

A chromosome-level assembly of blunt snout bream (Megalobrama amblycephala) reveals an expansion of olfactory receptor genes in freshwater fish

The number of olfactory receptor genes (ORs), which are responsible for detecting diverse odor molecules varies extensively among mammals as a result of frequent gene gains and losses that contribute to olfactory specialization. However, how OR expansions/contractions in fish are influenced by habitat and feeding habit and which OR subfamilies are important in each ecological niche is unknown. Here, we report a major OR expansion in a freshwater herbivorous fish, Megalobrama amblycephala, using a highly contiguous, chromosome-level assembly. We evaluate the possible contribution of OR expansion to habitat and feeding specialization by comparing the OR repertoire in 28 phylogenetically and ecologically diverse teleosts. In total, we analyzed > 4,000 ORs including 3,253 intact, 122 truncated, and 913 pseudogenes. The number of intact ORs is highly variable ranging from 20 to 279. We estimate that the most recent common ancestor of Osteichthyes had 62 intact ORs, which declined in most lineages except the freshwater Otophysa clade that has a substantial expansion in subfamily β and ε ORs. Across teleosts, we found a strong association between duplications of β and ε ORs and freshwater habitat. Nearly, all ORs were expressed in the olfactory epithelium (OE) in three tested fish species. Specifically, all the expanded β and ε ORs were highly expressed in OE of M. amblycephala. Together, we provide molecular and functional evidence for how OR repertoires in fish have undergone gain and loss with respect to ecological factors and highlight the role of β and ε OR in freshwater adaptation.

Whole genome sequencing of a snailfish from the Yap Trench (~7,000 m) clarifies the molecular mechanisms underlying adaptation to the deep sea

Hadal environments (depths below 6,000 m) are characterized by extremely high hydrostatic pressures, low temperatures, a scarce food supply, and little light. The evolutionary adaptations that allow vertebrates to survive in this extreme environment are poorly understood. Here, we constructed a high-quality reference genome for Yap hadal snailfish (YHS), which was captured at a depth of ~7,000 m in the Yap Trench. The final YHS genome assembly was 731.75 Mb, with a contig N50 of 0.75 Mb and a scaffold N50 of 1.26 Mb. We predicted 24,329 protein-coding genes in the YHS genome, and 24,265 of these genes were successfully functionally annotated. Phylogenetic analyses suggested that YHS diverged from a Mariana Trench snailfish approximately 0.92 million years ago. Many genes associated with DNA repair show evidence of positive selection and have expanded copy numbers in the YHS genome, possibly helping to maintain the integrity of DNA under increased hydrostatic pressure. The levels of trimethylamine N-oxide (TMAO), a potent protein stabilizer, are much higher in the muscles of YHS than in those of shallow-water fish. This difference is perhaps due to the five copies of the TMAO-generating enzyme flavin-containing monooxygenase-3 gene (fmo3) in the YHS genome and the abundance of trimethylamine (TMA)-generating bacteria in the YHS gut. Thus, the high TMAO content might help YHS adapt to high hydrostatic pressure by improving protein stability. Additionally, the evolutionary features of the YHS genes encoding sensory-related proteins are consistent with the scarce food supply and darkness in the hadal environments. These results clarify the molecular mechanisms underlying the adaptation of hadal organisms to the deep-sea environment and provide valuable genomic resources for in-depth investigations of hadal biology.

Hadal environments (depths below 6,000 m) are characterized by extremely high hydrostatic pressures, low temperatures, a scarce food supply, and little light. Fish are the only vertebrates inhabiting the hadal zone, and hadal snailfishes have been found in at least five geographically separated marine trenches. However, the genetic mechanisms that allow vertebrates to live in such extreme conditions are not well understood. Here, we constructed a high-quality reference genome for Yap hadal snailfish (YHS) captured at a depth of ~7,000 m in the Yap Trench, using long reads obtained by Pacific Biosciences Sequel sequencing. Comparative genomic analyses revealed that many genes associated with DNA repair show evidence of positive selection and have expanded copy numbers in the YHS genome, which potentially reflect the difficulty of maintaining DNA integrity under high hydrostatic pressure. Moreover, the five copies of the trimethylamine N-oxide (TMAO)-generating enzyme flavin-containing monooxygenase-3 gene (fmo3) and the abundance of trimethylamine (TMA)-generating bacteria in the YHS gut could provide enough TMAO to improve protein stability under hadal conditions. In addition, characteristics of the YHS sensory system genes were consistent with the scarce food supply and darkness in the hadal zone. Our results provide new insights into the molecular mechanisms underlying the adaptation of hadal organisms to the deep-sea environment and valuable genomic resources that will help further clarify hadal adaptations.