The draft genome of blunt snout bream (Megalobrama amblycephala) reveals the development of intermuscular bone and adaptation to herbivorous diet

A novel lncRNA MSTRG.59348.1 regulates muscle cells proliferation and innate immunity of Megalobrama amblycephala

In mammals, long non-coding RNAs (lncRNAs) play a regulatory role in gene expression, contribute to immune responses, and aid in pathogen elimination, primarily through interactions with RNA-binding proteins (RBPs). However, the role of lncRNAs in fish innate immunity and their interaction with RBPs remains uncertain. To investigate the immunomodulatory role of lncRNAs in Megalobrama amblycephala, we identified the novel lncRNA MSTRG.59348.1 and examined its function in the innate immune response to Aeromonas hydrophila infection. Localization studies in hepatocytes revealed that MSTRG.59348.1 is primarily located in the nucleus, suggesting its potential involvement in gene regulation, possibly through chromatin modification or other nuclear processes. The expression of MSTRG.59348.1 was significantly up-regulated after lipopolysaccharide (LPS) stimulation in liver cells. RNA-seq analysis of muscle cells revealed that genes differentially expressed following MSTRG.59348.1 overexpression were enriched in immune pathways. MSTRG.59348.1 overexpression significantly inhibited the expression of sting and ifn, and significantly up-regulated muscle cell viability and promoted cell proliferation by targeting sting, ifn, nf-κb1, and bcl2. Screening by RNA pull-down and mass spectrometry identified 57 RBPs interacting with MSTRG.59348.1, with functions enriched in immune pathways. Our results suggest that MSTRG.59348.1 plays a crucial regulatory role in fish antibacterial response, marking it as a significant subject for future research in innate immunity.

Integrated transcriptomic analysis reveals evolutionary and developmental characteristics of tendon ossification in teleost

BACKGROUND

Intermuscular bones (IBs) are segmental intramembranous ossifications located within myosepta. They share similarities with tendon ossification, a form of heterotopic ossification (HO). The mechanisms underlying IB formation remain incompletely understood.

RESULTS

In this study, we systematically analyzed transcriptome data across multiple tissues, species, time points, and resolutions in teleosts. First, we identified IB-specific expression genes using the tau index method. Through cross-species comparisons of the tendon development process, we discovered that candidate genes were primarily enriched in extracellular matrix organization, ossification, regulation of angiogenesis, and other related processes. We also revealed that some of these candidate genes are abnormally expressed in runx2b-/- zebrafish, which lack IBs. To clarify the trajectory of cell differentiation during IB formation, we demonstrated that myoseptal stem cells differentiate into osteoblasts, fibroblasts, and tenocytes in runx2b+/+ zebrafish. However, in runx2b-/- zebrafish, the differentiation of myoseptal stem cell into osteoblast was inhibited, while differentiation into clec3bb + tenocyte and fibroblast was enhanced. Additionally, runx2b deficiency led to the upregulation of clec3bb expression in the clec3bb + tenocyte cluster. Notably, a compensatory effect was observed in cell differentiation and gene expression in runx2b-/- zebrafish, suggesting that runx2b and the candidate genes, such as clec3bb, were involved in the gene network of IB development.

CONCLUSIONS

The results elucidate cell differentiation process during tendon ossification in teleosts and identify the key factor clec3bb involved in this process. These findings provide a foundation for understanding tendon ossification in teleosts and for further research on tendon ossification in mammals.

Whole genome resequencing reveals the correlation between selection signatures and adaptability of Micropterus salmoides to artificial fed

Largemouth bass (Micropterus salmoides, LMB) is an important aquaculture species due to its excellent flesh quality and environmental adaptability. It has been continuously introduced to many countries and cultured for decades. Here, an LMB population was used for selective breeding to improve growth rate and feed adaptability. After five generations of breeding, the growth rate improved by 38%, and feed adaptability improved by 22% compared to the non-breeding population. To study the underlying genetic mechanism, 100 LMB from the breeding population and 100 from the non-breeding population were sampled for whole-genome resequencing. The population genetics analysis shows that the breeding population has a higher inbreeding coefficient and linkage disequilibrium (LD) level, a lower nucleic acid diversity and effective population size (Ne). Using [Formula: see text] (fixation index), we found that the average [Formula: see text] value between the two populations was 0.07, with the highest [Formula: see text] value reaching 0.38, which overlaps with the trypsin gene. Additionally, other genes exhibiting high [Formula: see text] values are associated with functions such as neural development, glucose metabolism, and growth. Using [Formula: see text] and nucleic acid diversity as criteria, we identified 698 genes that are positively selected in the breeding population, and gene functional enrichment analysis shows that 36 genes are related to the olfactory receptor pathway. Overall, our study found that multiple genes were selected in the LMB breeding population. These genes may be associated with adaptation and digestion of artificial feed in fish.

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

Betaine Supplementation Into High-Carbohydrate Diets Improves Feed Efficiency and Liver Health of Megalobrama amblycephala by Increasing Taurine Synthesis

Dietary betaine supplementation has been reported to alleviate the adverse effects of high-carbohydrate diets on Megalobrama amblycephala, while the regulatory mechanism remains largely unknown. In the present study, a 79-day feeding trial was conducted with 450 juvenile Megalobrama amblycephala (average weight 6.75 ± 0.10 g), which were fed with five high-carbohydrate diets (43%) supplementing betaine at 0% (CD group), 0.2% (0.2Bet group), 0.4% (0.4Bet group), 0.8% (0.8Bet group), and 1.6% (1.6Bet group), respectively. Results showed M. amblycephala in 0.8Bet group exhibited the best growth performance, indicated by the largest weight gain ratio (142.88%) and least feed conversion ratio (1.63). Moreover, liver health was promoted in 0.8Bet group, with decreased number of non-nucleated cells and less lipid accumulation, which was accompanied by the lowest hepatosomatic index (1.38%). In order to further illustrate the regulatory mechanism, metabolites assay indicated that dietary betaine supplementation significantly increased plasma contents of methionine, serine, hypotaurine, and taurine, but did not affect plasma contents of cystathionine, cystine, or cysteic acid. Accordingly, the mRNA expressions of cysteine sulfinate decarboxylase in cysteine sulfinic acid pathway and cysteamine dioxygenase (ADO) in sulfinic acid (CS) pathway, which were both involved in taurine synthesis, were also upregulated in the liver. Meanwhile, the microbial communities in M. amblycephala intestine were more stable and uniform with betaine supplementation. Therefore, dietary betaine supplementation may exert its protective roles in improving feed efficiency and liver health of M. amblycephala via promoting de novo taurine synthesis and stabilizing intestinal microbial communities.

The Chinese longsnout catfish genome provides novel insights into the feeding preference and corresponding metabolic strategy of carnivores

Fish show variation in feeding habits to adapt to complex environments. However, the genetic basis of feeding preference and the corresponding metabolic strategies that differentiate feeding habits remain elusive. Here, by comparing the whole genome of a typical carnivorous fish (Leiocassis longirostris Günther) with that of herbivorous fish, we identify 250 genes through both positive selection and rapid evolution, including taste receptor taste receptor type 1 member 3 (tas1r3) and trypsin We demonstrate that tas1r3 is required for carnivore preference in tas1r3-deficient zebrafish and in a diet-shifted grass carp model. We confirm that trypsin correlates with the metabolic strategies of fish with distinct feeding habits. Furthermore, marked alterations in trypsin activity and metabolic profiles are accompanied by a transition of feeding preference in tas1r3-deficient zebrafish and diet-shifted grass carp. Our results reveal a conserved adaptation between feeding preference and corresponding metabolic strategies in fish, and provide novel insights into the adaptation of feeding habits over the evolution course.

Genomic insights into the seawater adaptation in Cyprinidae

BACKGROUND

Cyprinidae, the largest fish family, encompasses approximately 367 genera and 3006 species. While they exhibit remarkable adaptability to diverse aquatic environments, it is exceptionally rare to find them in seawater, with the Far Eastern daces being of few exceptions. Therefore, the Far Eastern daces serve as a valuable model for studying the genetic mechanisms underlying seawater adaptation in Cyprinidae.

RESULTS

Here, we sequenced the chromosome-level genomes of two Far Eastern daces (Pseudaspius brandtii and P. hakonensis), the two known cyprinid fishes found in seawater, and performed comparative genomic analyses to investigate their genetic mechanism of seawater adaptation. Demographic history reconstruction of the two species reveals that their population dynamics are correlated with the glacial-interglacial cycles and sea level changes. Genomic analyses identified Pseudaspius-specific genetic innovations related to seawater adaptation, including positively selected genes, rapidly evolving genes, and conserved non-coding elements (CNEs). Functional assays of Pseudaspius-specific variants of the prolactin (prl) gene showed enhanced cell adaptation to greater osmolarity. Functional assays of Pseudaspius specific CNEs near atg7 and usp45 genes suggest that they exhibit higher promoter activity and significantly induced at high osmolarity.

CONCLUSIONS

Our results reveal the genome-wide evidence for the evolutionary adaptation of cyprinid fishes to seawater, offering valuable insights into the molecular mechanisms supporting the survival of migratory fish in marine environments. These findings are significant as they contribute to our understanding of how cyprinid fishes navigate and thrive in diverse aquatic habitats, providing useful implications for the conservation and management of marine ecosystems.

Characterization of type II IFNs and their receptors in a cyprinid fish, the blunt snout bream Megalobrama amblycephala

Type II interferons (IFNs) are a key class of molecules regulating innate and adaptive immunity in vertebrates. In the present study, two members of the type II IFNs, IFN-γ and IFNγ-rel, were identified in the blunt snout bream (Megalobrama amblycephala). The open reading frame (ORF) of IFN-γ and IFNγ-rel was found to have 564 bp and 492 bp, encoding 187 and 163 amino acids, with the first 26 and 24 amino acids being the signal peptide, respectively. IFN-γ and IFNγ-rel genes showed a high degree of similarity to their zebrafish homologues, being 76.9 % and 58.9 %, respectively. In the phylogenetic tree, IFN-γ and IFNγ-rel were clustered with homologous genes in cyprinids. In blunt snout bream, IFN-γ and IFNγ-rel were constitutively expressed in trunk kidney, head kidney, spleen, liver, heart, muscle, gill, intestine and brain and were significantly up-regulated by poly (I:C) induction in head kidney, spleen, liver, gill and intestine. Using recombinant proteins of IFN-γ and IFNγ-rel, the surface plasmon resonance (SPR) results showed that IFN-γ was bound to CRFB6, CRFB13 and CRFB17, but mainly to CRFB6 and CRFB13, whereas IFN-γrel bound mainly to CRFB17 and had no affinity with CRFB6. These results contribute to a better understanding on type II IFNs and their receptor usage in teleost fish.

PacBio Full-Length Transcriptome of a Tetraploid Sinocyclocheilus multipunctatus Provides Insights into the Evolution of Cavefish

Sinocyclocheilus multipunctatus is a second-class nationally protected wild animal in China. As one of the cavefish, S. multipunctatus has strong adaptability to harsh subterranean environments. In this study, we used PacBio SMRT sequencing technology to generate a first representative full-length transcriptome for S. multipunctatus. Sequence clustering analysis obtained 232,126 full-length transcripts. Among all transcripts, 40,487 were annotated in public databases, while 70,300 microsatellites, 2384 transcription factors, and 16,321 long non-coding RNAs were identified. The phylogenetic tree showed that S. multipunctatus shows a closer relationship to Carassius auratus and Cyprinus carpio, phylogenetically diverging from the common ancestor ~14.74 million years ago (Mya). We also found that between 15.6 and 17.5 Mya, S. multipunctatus also experienced an additional whole-genome duplication (WGD) event, which may have promoted the species evolution of S. multipunctatus. Meanwhile, the overall rates of evolutionary of polyploid S. multipunctatus were significantly higher than those of the other cyprinids, and 220 positively selected genes (PSGs) were identified in two sub-genomes of S. multipunctatus. These PSGs are likely to fulfill critical roles in the process of adapting to diverse cave environments. This study has the potential to facilitate future investigations into the genomic characteristics of S. multipunctatus and provide valuable insights into revealing the evolutionary history of polyploid S. multipunctatus.

Characterization and Expression Analysis of Genes from Megalobrama amblycephala Encoding Hemoglobins with Extracellular Microbicidal Activity

Hemoglobin (Hb) usually comprises two α and two β subunits, forming a tetramer responsible for oxygen transportation and storage. Few studies have elucidated fish hemoglobin immune functions. Megalobrama amblycephala is a freshwater-cultured fish prevalent in China. We identified two M. amblycephala hemoglobin subunits and analyzed their expression patterns and antibacterial activities. The respective full-length cDNA sequences of the M. amblycephala Hb α (MaHbα) and β (MaHbβ) subunits were 588 and 603 bp, encoding 143 and 148 amino acids. MaHbα and MaHbβ were highly homologous to hemoglobins from other fish, displaying typical globin-like domains, most heme-binding sites, and tetramer interface regions highly conserved in teleosts. In phylogenetic analyses, the hemoglobin genes from M. amblycephala and other cypriniformes clustered into one branch, and those from other fishes and mammals clustered into other branches, revealing fish hemoglobin conservation. These M. amblycephala Hb subunits exhibit different expression patterns in various tissues and during development. MaHbα is mainly expressed in the blood and brain, while MaHbβ gene expression is highest in the muscle. MaHbα expression was detectable and abundant post-fertilization, with levels fluctuating during the developmental stages. MaHbβ expression began at 3 dph and gradually increased. Expression of both M. amblycephala Hb subunits was down-regulated in most examined tissues and time points post-Aeromonas hydrophila infection, which might be due to red blood cell (RBC) and hematopoietic organ damage. Synthetic MaHbα and MaHbβ peptides showed excellent antimicrobial activities, which could inhibit survival and growth in five aquatic pathogens. Two M. amblycephala hemoglobin subunits were identified, and their expression patterns and antibacterial activities were analyzed, thereby providing a basis for the understanding of evolution and functions of fish hemoglobins.

Molecular characterization of nineteen cytokine receptor family B (CRFB) members, CRFB1, CRFB2, CRFB4-17, with three CRFB9 and two CRFB14 in a cyprinid fish, the blunt snout bream Megalobrama amblycephala

The class II cytokine receptor family members are receptors of class 2 helical cytokines in mammals, and are named cytokine receptor family B (CRFB) in fish. In zebrafish, sixteen members, including CRFB1, CRFB2 and CRFB4-17 were reported. With the availability of genome sequence, a total of nineteen CRFBs was identified in the blunt snout bream (Megalobrama amblycephala), including CRFB1, CRFB2, CRFB4-17 with the presence of three CRFB9 isoforms, and two CRFB14 isoforms. These CRFB molecules contain well conserved features, such as fibronectin type III (FNIII) domain, transmembrane and intracellular domains as other class II cytokine receptors, and are phylogenetically grouped into thirteen clades with their homologues from other species of fish. The CRFB genes were constitutively expressed in organs/tissues examined in the fish. The finding of more CRFB members in the bream may provide clues to understand possible receptor-ligand interaction and their diversity from an evolutionary point of view.

Metabolites, gene expression, and gut microbiota profiles suggest the putative mechanisms via which dietary creatine increases the serum taurine and g-ABA contents in Megalobrama amblycephala

A 90-day experiment was conducted to explore the effects of creatine on growth performance, liver health status, metabolites, and gut microbiota in Megalobrama amblycephala. There were 6 treatments as follows: control (CD, 29.41% carbohydrates), high carbohydrate (HCD, 38.14% carbohydrates), betaine (BET, 1.2% betaine + 39.76% carbohydrates), creatine 1 (CRE1, 0.5% creatine + 1.2% betaine + 39.29% carbohydrates), creatine 2 (CRE2, 1% creatine + 1.2% betaine + 39.50% carbohydrates), and creatine 3 (CRE3, 2% creatine + 1.2% betaine + 39.44% carbohydrates). The results showed that supplementing creatine and betaine together reduced the feed conversion ratio significantly (P < 0.05, compared to CD and HCD) and improved liver health (compared to HCD). Compared with the BET group, dietary creatine significantly increased the abundances of Firmicutes, Bacteroidota, ZOR0006, and Bacteroides and decreased the abundances of Proteobacteria, Fusobacteriota, Vibrio, Crenobacter, and Shewanella in the CRE1 group. Dietary creatine increased the content of taurine, arginine, ornithine, γ-aminobutyric acid (g-ABA), and creatine (CRE1 vs. BET group) and the expression of creatine kinase (ck), sulfinoalanine decarboxylase (csad), guanidinoacetate N-methyltransferase (gamt), glycine amidinotransferase (gatm), agmatinase (agmat), diamine oxidase1 (aoc1), and glutamate decarboxylase (gad) in the CRE1 group. Overall, these results suggested that dietary supplementation of creatine (0.5-2%) did not affect the growth performance, but it altered the gut microbial composition at the phylum and genus levels, which might be beneficial to the gut health of M. amblycephala; dietary creatine also increased the serum content of taurine by enhancing the expressions of ck and csad and increased the serum content of g-ABA by enhancing the arginine content and the expressions of gatm, agmat, gad, and aoc1.

Janus kinase 1 in Megalobrama amblycephala: Identification, phylogenetic analysis and expression profiling after Aeromonas hydrophila infection

Janus kinase 1 (JAK1), a member of the JAK family, plays an essential and non-redundant role in the mammalian immune system. However, the potential role of JAK1 in fish immune response remains largely unclear. In the present study, the JAK1 gene of Megalobrama amblycephala (MamJAK1) was identified and characterized. The open reading frame (ORF) of MamJAK1 was 3462 bp, encoding 1153 amino acids. MamJAK1 consists of four common domains of the JAK family, including B41, SH2, STyrKc (a pseudo kinase domain), and TyrKc (a kinase domain). Phylogenetic analysis showed that JAK1s are divided into two evolutionary clades, one containing fish JAK1s, and the other containing JAK1s from other vertebrates. The results of quantitative real-time PCR (qPCR) showed that in healthy M. amblycephala, MamJAK1 mRNA was highest expressed in blood, followed by spleen, intestine and mid-kidney, and lowly expressed in other tissues including gill, liver, head kidney, muscle, brain and heart. After Aeromonas hydrophila infection, the expression of MamJAK1 mRNA was significantly induced in four selected tissues including spleen, mid-kidney, liver and intestine, reaching a peak at 24 hpi (hour post infection) in spleen and mid-kidney, at 12 hpi in liver and at 4 hpi in intestine, and then the expression level was restricted to control levels at 72 or 120 hpi. In addition, the results of Western blot showed that the phosphorylation level of MamJAK1 protein in spleen and mid-kidney increased significantly after A. hydrophila infection, although MamJAK1 protein did not change obviously. Further, the JAK1 phosphorylation in Ctenopharyngodon idellus kidney (CIK) cells was found to be significantly induced by LPS stimulation and IL-6R over-expression. The results above suggest that MamJAK1 may play an essential role in the immune response against bacterial infection through the IL-6R mediated JAK1/STAT signaling pathway, which further deepen our understanding of JAK1 and provides a potential target for the treatment and prevention of bacterial diseases in teleost.

Functional exploration of SNP mutations in HIF2αb gene correlated with hypoxia tolerance in blunt snout bream (Megalobrama amblycephala)

Blunt snout bream (Megalobrama amblycephala) is sensitive to hypoxia environment. Hypoxia-inducible factor (HIF) is the most critical factor in the HIF pathway, which strictly regulates the hypoxia stress process of fish. In this study, we found six hifα genes in blunt snout bream that demonstrated different expressions under hypoxia conditions. In HEK293T cells, all six hifαs were detected to activate the HRE region by luciferase reporter assay. More importantly, we identified two linkage-disequilibrium SNP sites at exon 203 and 752 of the hif2αb gene in blunt snout bream. Haplotype II (A²⁰³A⁷⁵²) and its homozygous diplotype II (A²⁰³A²⁰³A⁷⁵²A⁷⁵²) appeared frequently in a selected strain of blunt snout bream with hypoxia tolerance. Diplotype II has a lower oxygen tension threshold for loss of equilibrium (LOE_crit) over a similar range of temperatures. Moreover, its erythrocyte number increased significantly (p < 0.05) than those in diplotype I and diplotype III strains at 48 h of hypoxia. The enzymes related with hypoxia tolerant traits, i.e., reduced glutathione, superoxide dismutase, and catalase, were also significantly (p < 0.05) induced in diplotype II than in diplotype I or III. In addition, the expression of epo in the liver of diplotype II was significantly (p < 0.01) higher than that in the diplotype I or III strains at 48 h of hypoxia. Taken together, our results found that the hypoxia-tolerant-related diplotype II of hif2αb has the potential to be used as a molecular marker in future genetic breeding of hypoxia-tolerant strain.

Intelectin enhances the phagocytosis of macrophages via CDC42-WASF2-ARPC2 signaling axis in Megalobrama amblycephala

Intelectin has been identified in various vertebrates and plays an important role in the host immune system. In our previous studies, recombinant Megalobrama amblycephala intelectin (rMaINTL) protein with excellent bacterial binding and agglutination activities enhances the phagocytic and killing activities of macrophages in M. amblycephala; however, the underlying regulatory mechanisms remain unclear. The present study showed that treatment with Aeromonas hydrophila and LPS induced the expression of rMaINTL in macrophages, and its level and distribution in macrophages or kidney tissue markedly increased after incubation or injection with rMaINTL. The cellular structure of macrophages was significantly affected after incubation with rMaINTL, resulting in an increased surface area and pseudopodia extension, which might contribute to enhancing the phagocytic ability of macrophages. Then, digital gene expression profiling analysis of the kidneys from rMaINTL-treated juvenile M. amblycephala identified some phagocytosis-related signaling factors that were enriched in pathways involved in the regulation of the actin cytoskeleton. In addition, qRT-PCR and western blotting verified that rMaINTL upregulated the expression of CDC42, WASF2, and ARPC2 in vitro and in vivo; however, the expression of these proteins was inhibited by a CDC42 inhibitor in macrophages. Moreover, CDC42 mediated the promotion of rMaINTL on actin polymerization by increasing the F-actin/G-actin ratio, which led to the extension of pseudopodia and remodeling of the macrophage cytoskeleton. Furthermore, the enhancement of macrophage phagocytosis by rMaINTL was blocked by the CDC42 inhibitor. These results suggested that rMaINTL induced the expression of CDC42 as well as the downstream signaling molecules WASF2 and ARPC2, thereby facilitating actin polymerization to promote cytoskeletal remodeling and phagocytosis. Overall, MaINTL enhanced the phagocytosis activity of macrophages in M. amblycephala via activation of the CDC42-WASF2-ARPC2 signaling axis.

Comparative Expression Profiling Reveals the Regulatory Effects of Dietary Mannan Oligosaccharides on the Intestinal Immune Response of Juvenile Megalobrama amblycephala against Aeromonas hydrophila Infection

Mannan oligosaccharides (MOS) are functional oligosaccharides with beneficial effects on the non-specific immunity of Megalobrama amblycephala, but systematic studies on the immunomodulatory mechanisms of MOS are still lacking. To investigate the protective mechanisms of three different levels of dietary MOS supplementation on the intestinal immunity of juvenile M. amblycephala, comparative digital gene expression (DGE) profiling was performed. In this study, 622 differentially expressed genes (DEGs) were identified, while the similar expression tendency of 34 genes by qRT-PCR validated the accuracy of the DGE analyses. Gene Ontology (GO) enrichment revealed that the DEGs were mainly enriched in two functional categories of biological process and molecular function. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the DEGs were mainly related to complement and coagulation cascades, coagulation cascades, platelet activation, natural killer cell mediated cytotoxicity, Fc gamma R-mediated phagocytosis and antigen processing and presentation. In addition, the pro-inflammatory, apoptosis and tight junction-related genes were more significantly up-regulated upon infection in the dietary MOS groups to enhance host immune functions and maintain the stability of the intestinal barrier. These results will be helpful to clarify the regulatory mechanism of MOS on the intestinal immunity of M. amblycephala and lay the theoretical foundation for the prevention and protection of fish bacterial diseases.

Dietary Leucine Supplementation Improves Muscle Fiber Growth and Development by Activating AMPK/Sirt1 Pathway in Blunt Snout Bream (Megalobrama amblycephala)

This research is aimed at evaluating the effects of leucine supplementation on muscle fibers growth and development of blunt snout bream through a feeding trial and a primary muscle cells treatment. An 8-week trial with diets containing 1.61% leucine (LL) or 2.15% leucine (HL) was conducted in blunt snout bream (mean initial weight = 56.56 ± 0.83 g). Results demonstrated that the specific gain rate and the condition factor of fish in the HL group were the highest. The essential amino acids content of fish fed HL diets was significantly higher than that fed LL diets. The texture (hardness, springiness, resilience, and chewiness), the small-sized fiber ratio, fibers density, and sarcomere lengths in fish all obtained the highest in the HL group. Additionally, the proteins expression related with the activation of the AMPK pathway (p-Ampk, Ampk, p-Ampk/Ampk, and Sirt1) and the expression of genes (myogenin (myog), myogenic regulatory factor 4 (mrf4) and myoblast determination protein (myod), and protein (Pax7) related to muscle fiber formation were significantly upregulated with increasing level of dietary leucine. In vitro, the muscle cells were treated with 0, 40 and 160 mg/L leucine for 24 h. The results showed that treated with 40 mg/L leucine significantly raised the protein expressions of BCKDHA, Ampk, p-Ampk, p-Ampk/Ampk, Sirt1, and Pax7 and the gene expressions of myog, mrf4, and myogenic factor 5 (myf5) in muscle cells. In summary, leucine supplementation promoted muscle fibers growth and development, which may be related to the activation of BCKDH and AMPK.

Transcriptomic analysis of intermuscular bone development in barbel steed (Hemibarbus labeo)

Intermuscular bones (IBs), which are little, bony spicules in muscle, are embedded in lower teleosts' myosepta. Despite the importance of studying IB development in freshwater aquaculture species, the genes associated with IB development need to be further explored. In the present study, we identified four stages of IB development in barbel steed (Hemibarbus labeo), namely stage 1: IBs have not emerged, stage 2: a few small IBs have emerged in the tail, stage 3: longer IBs gradually emerged in the tail and stage 4: all of the IBs in the tail are mature and long, via Alizarin red staining. Subsequently, we used the HiseqXTen platform to sequence and de novo assemble the transcriptome of epaxial muscle (between 35th and 40th myomere) of barbel steed at 29 days (stage 1) and 42 days (stage 3) after hatching. A total of 190,814 unigenes were obtained with an average length and N50 of 648 bp and 1027 bp, respectively. We found 2174 differentially expressed genes (DEGs) between stages 1 and 3, of which 378 and 1796 were up- and down-regulated, respectively. Functional enrichment analysis showed that several DEGs functioned in ossification, positive regulation of osteoblast differentiation, osteoblast differentiation, and BMP signaling pathway, and were further enriched in signal pathway, including osteoclast differentiation, TGF-β signaling pathway, cytokine-cytokine receptor interaction, Jak-STAT signaling pathway, and other KEEG pathways. In conclusion, we identified genes that may be related to IB development, such as kazal type serine peptidase inhibitor domain 1 (KAZALD1), extracellular matrix protein 1 (ECM1), tetranectin, bone morphogenetic protein 1 (bmp1), acid phosphatase 5 (ACP5), collagen type XI alpha 1 chain (COL11A1), matrix metallopeptidase 9 (MMP9), pannexin-3 (PANX3), sp7 transcription factor (Sp7), and c-x-c motif chemokine ligand 8 (CXCL8), by comparing the transcriptomes of epaxial muscle before and after IB ossification. This study provided a theoretical basis for identifying the molecular mechanisms underlying IB development in fish.

Role of feeding specialization in taste receptor loss: insights from sweet and umami receptor evolution in Carnivora

Controversy and misunderstanding surround the role of feeding specialization in taste receptor loss in vertebrates. We refined and tested the hypothesis that this loss is caused by feeding specializations. Specifically, feeding specializations were proposed to trigger time-dependent process of taste receptor loss through deprivation of benefit of using the receptor's gustatory function. We propose that this process may be accelerated by abiotic environmental conditions or decelerated/stopped because of extragustatory functions of the receptor's protein(s). As test case we used evolution of the sweet (TAS1R2+TAS1R3) and umami (TAS1R1+TAS1R3) receptors in Carnivora (dogs, cats, and kin). We predicted these receptors' absence/presence using data on presence/absence of inactivating mutations in these receptors' genes and data from behavioral sweet/umami preference tests. We identified 20 evolutionary events of sweet (11) or umami (9) receptor loss. These events affected species with feeding specializations predicted to favor sweet/umami receptor loss (27 and 22 species, respectively). All species with feeding habits predicted to favor sweet/umami receptor retention (11 and 24, respectively) were found to retain that receptor. Six species retained the sweet (5) or umami (1) receptor despite feeding specialization predicted to favor loss of that receptor, which can be explained by the time dependence of sweet/umami receptor loss process and the possible decelerating effect of TAS1R extragustatory functions so that the sweet/umami receptor process is ongoing in these species. Our findings support the idea that feeding specialization leads to taste receptor loss and is the main if not only triggering factor for evolutionary loss of taste receptors.

Molecular Characterization, Expression, Evolutionary Selection, and Biological Activity Analysis of CD68 Gene from Megalobrama amblycephala

CD68 is a highly glycosylated transmembrane glycoprotein that belongs to the lysosome-associated membrane glycoprotein family and is involved in various immune processes. In this study, Megalobrama amblycephala CD68 (MaCD68) was cloned and characterized, and its expression patterns and evolutionary characteristics were analyzed. The coding region of MaCD68 was 987 bp, encoding 328 amino acids, and the predicted protein molecular weight was 34.9 kDa. MaCD68 contained two transmembrane helical structures and 18 predicted N-glycosylation sites. Multiple sequence alignments showed that the MaCD68 protein had high homology with other fish, and their functional sites were also highly conserved. Phylogenetic analysis revealed that MaCD68 and other cypriniformes fish clustered into one branch. Adaptive evolution analysis identified several positively selected sites of teleost CD68 using site and branch-site models, indicating that it was under positive selection pressure during evolution. Quantitative real-time reverse transcription polymerase chain reaction analysis showed that MaCD68 was highly expressed in the head kidney, spleen, and heart. After Aeromonas hydrophila infection, MaCD68 was significantly upregulated in all tested tissues, peaking at 12 h post-infection (hpi) in the kidney and head kidney and at 120 hpi in the liver and spleen, suggesting that MaCD68 participated in the innate immune response of the host against bacterial infection. Immunohistochemical and immunofluorescence analyses also showed that positive signals derived from the MaCD68 protein were further enhanced after bacterial and lipopolysaccharide treatment, which suggested that MaCD68 is involved in the immune response and could be used as a macrophage marker. Biological activity analysis indicated that recombinant MaCD68 (rMaCD68) protein had no agglutination or bactericidal effects on A. hydrophila but did have these effects on Escherichia coli. In conclusion, these results suggest that MaCD68 plays a vital role in the immune response against pathogens, which is helpful in understanding the immune responses and mechanisms of M. amblycephala.

Breaking the reproductive barrier of divergent species to explore the genomic landscape

Background: Climate change will have significant consequences for species. Species range shifts induce the emergence of new hybrid zones or the spatial displacement of pre-existing ones. These hybrid zones may become more porous as alleles are passed from one species to another. Currently, hybridization between highly divergent species living in sympatry seems extremely limited. Indeed, this phenomenon involves breaking two barriers. The first is the pre-mating barrier, related to the reproductive phenology of the two species. The second is the post-zygotic barrier, related to the genetic divergence between these species. Here, we were interested in identifying new hybridization patterns and potential implications, especially in the context of environmental modifications.

Methods: We sampled Telestes souffia and Parachondrostoma toxostoma wild specimens from different locations across France and genotyped them for SNP markers. We identified discriminant loci using F1-hybrid specimens and parental species and performed principal component analysis and Bayesian model-based clustering to analyze phylogenetic information. Furthermore, we assessed deviation in allele frequency from F1 to F2 and for Hardy–Weinberg equilibrium for F2 and assessed gene function associated with two F2 cohorts.

Results: We demonstrate that by breaking the ecological barrier, massive introgressive hybridization is possible between two endemic lineages of Cyprinidae belonging to two distinct genera. For both cohorts studied (=2 cm and >2 cm), a large majority of loci (>88%) presented no deviation in allele frequency and no departure from the Hardy–Weinberg equilibrium. For individuals beyond the 2 cm stage, two phenomena were observed. The first was an allelic imbalance in favor of P. toxostoma, for some genomic regions, with genes involved in developmental regulatory processes, cytoskeletal organization, and chromosome organization. The second was an excess of heterozygous loci coupled with an equilibrium of allelic frequencies for genes involved in immune response and kidney/liver development. Moreover, the 2 cm-sized specimens with high mortality yielded a particular genomic signature.

Conclusion: Our study displayed important results for understanding the early stages of hybridization between divergent lineages and predicting the emergence of future hybrid zones in the wild. Moreover, this hybridization generates a wide spectrum of hybrids that are a potential source of important evolutionary novelties.

Chromosome-Level Genome Assembly and Transcriptome Comparison Analysis of Cephalopholis sonnerati and Its Related Grouper Species

The tomato hind, Cephalopholis sonnerati, is a bottom-dwelling coral reef fish, which is widely distributed in the Indo-Pacific and Red Sea. C. sonnerati also features complex social structures and behaviour mechanisms. Here, we present a high-quality, chromosome-level genome assembly for C. sonnerati that was derived using PacBio sequencing and Hi-C technologies. A 1043.66 Mb genome with an N50 length of 2.49 Mb was assembled, produced containing 795 contigs assembled into 24 chromosomes. Overall, 97.2% of the complete BUSCOs were identified in the genome. A total of 26,130 protein-coding genes were predicted, of which 94.26% were functionally annotated. Evolutionary analysis revealed that C. sonnerati diverged from its common ancestor with E. lanceolatus and E. akaara approximately 41.7 million years ago. In addition, comparative genome analyses indicated that the expanded gene families were highly enriched in the sensory system. Finally, we found the tissue-specific expression of 8108 genes. We found that these tissue-specific genes were highly enriched in the brain. In brief, the high-quality, chromosome-level reference genome will provide a valuable genome resource for studies of the genetic conservation, resistance breeding, and evolution of C. sonnerati.

Identification, expression patterns, evolutionary characteristics and recombinant protein activities analysis of CD209 gene from Megalobrama amblycephala

CD209 is a type II transmembrane protein in the C-type lectin family, which is involved in the regulation of innate and adaptive immune system. Although it has been widely studied in mammals, but little has been reported about fish CD209 genes. In the present study, Megalobrama amblycephala CD209 (MaCD209) gene was cloned and characterized, its expression patterns, evolutionary characteristics, agglutinative and bacteriostatic activities were analyzed. These results showed that the open reading frame (ORF) of MaCD209 gene was 795 bp, encoding 264 aa, and the calculated molecular weight of the encoded protein was 29.7 kDa. MaCD209 was predicted to contain 2 N-glycosylation sites, 1 functional domain (C-LECT-DC-SIGN-like) and 1 transmembrane domain. Multiple sequence alignment showed that the amino acid sequence of MaCD209 was highly homologous with that of partial fishes, especially the highly conserved C-LECT-DC-SIGN-like domain and functional sites of CD209. Phylogenetic analysis showed that the CD209 genes from M. amblycephala and other cypriniformes fishes were clustered into one group, which was reliable and could be used for evolutionary analysis. Then, adaptive evolutionary analysis of teleost CD209 was conducted, and several positively selected sites were identified using site and branch-site models. Quantitative real-time PCR analysis showed that MaCD209 gene was highly expressed in the liver and heart. Moreover, the expression of MaCD209 was significantly induced upon Aeromonas hydrophila infection, with the peak levels at 4 h or 12 h post infection. The immunohistochemical analysis also revealed increased distribution of MaCD209 protein post bacterial infection. In addition, recombinant MaCD209 (rMaCD209) protein was prepared using a pET32a expression system, which showed excellent bacterial binding and agglutinative activities in a Ca²⁺-independent manner. However, rMaCD209 could only inhibit the proliferation of Escherichia coli rather than A. hydrophila. In conclusion, this study identified the MaCD209 gene, detected its expression and evolutionary characteristics, and evaluated the biological activities of rMaCD209 protein, which would provide a theoretical basis for understanding the evolution and functions of fish CD209 genes.

Evolution and expression analysis of STAT family members in blunt snout bream (Megalobrama amblycephala)

The Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway is involved in regulating the body's immunity, cell proliferation, differentiation, and apoptosis. Members of the STAT family have been extensively studied in different mammalian species. However, there are few studies on the STAT family genes in farmed economic fish. In this study, eight STAT genes including STAT1a, STAT1b, STAT2, STAT3, STAT4, STAT5a, STAT5b and STAT6, in blunt snout bream (Megalobrama amblycephala), an economically important fish in China, were identified and characterized. Analyses of gene location, phylogeny and conserved synteny were conducted to infer the evolutionary origin of these STAT family genes. Furthermore, the evolutionary origin model of STATs was constructed based on the 2R hypothesis and teleost genome duplication (TGD) hypothesis, which clarified the evolutionary origin of the eight STATs in blunt snout bream. Besides, expression of the eight STATs was detected in 10 tissues of healthy blunt snout bream, which showed different expression patterns, and all had the highest level in the blood. In addition, expression of the STATs was significantly induced in the spleen, liver, and kidney after infection of Aeromonas hydrophila, suggesting that they play an important role in protecting the host from pathogens. In general, the evolution of cytokine-related genes parallels that of the immune system, which has likely been a main evolutionary driver. Therefore, the evolutionary model of STAT genes, constructed in the non-model organism pioneeringly, may provide some enlightenment for the evolution of the fish STAT family genes and their involvement in the immune function.

Population Genomics of Megalobrama Provides Insights into Evolutionary History and Dietary Adaptation

Simple Summary

Megalobrama is the economically most important freshwater fish genus in China. In recent years, germplasm resources of Megalobrama have been depleting as a result of environmental degradation and artificial factors. In this study, we established the whole genome database of Megalobrama populations using the whole genome re-sequencing technology, explored population genetic structure, and inferred comprehensive evolutionary relationships using principal component analysis and population structure analysis. Furthermore, employing selective sweep analysis, we identified candidate genes related to variations in feeding habits, revealing the molecular mechanisms of environmental adaptability in Megalobrama populations. Taken together, this study describes the population history and genetic diversity of Megalobrama populations and also the molecular mechanisms likely involved their environmental adaptability. These findings will make a substantial contribution to conservation and utilization of Megalobrama germplasm resources.

Abstract

Megalobrama, a genus of cyprinid fish, is an economically important freshwater fish widely distributed in major waters of China. Here, we report the genome resequencing of 180 Megalobrama fish including M. amblycephala, M. skolkovii, M. hoffmanni, and M. pellegrini. Population structure indicated that geographically divergent Megalobrama populations were separated into six subgroups. A phylogenetic tree showed that M. skolkovii was more closely related to M. pellegrini than other species and M. hoffmanni was clustered apart from other Megalobrama species, showing a high nucleotide diversity in geographic groups. Treemix validated gene flow from M. amblycephala to M. skolkovii, suggesting that introgression may provide an important source of genetic variation in the M. skolkovii populations. According to the demographic history analysis, it is speculated that Megalobrama might have been originally distributed in the Pearl River with some spread to Hainan Island and northern China due to lower sea levels during the glacial period. Whole-genome selective sweeps analysis demonstrated that M. amblycephala likely developed an enhanced energy metabolism mostly through fatty acid degradation pathways whereas M. hoffmanni possibly regulate lipid absorption via the cholesterol metabolism pathway. Taken together, this study provides a valuable genomic resource for future genetic investigations aiming to improve genome-assisted breeding of Megalobrama species.

Comparison of Myosepta Development and Transcriptome Profiling between Blunt Snout Bream with and Tilapia without Intermuscular Bones

Simple Summary

The presence or absence of intermuscular bones (IBs) is directly related to the economic and edible value of fish. The specific regulatory mechanism of IB formation is not completely known yet. Here, we explored the molecular mechanisms that regulate the formation of IBs based on histological analysis, transcriptome profiling, and gene expression quantification using M. amblycephala (with IBs) and O. niloticus (without IBs) as models. As a result, we identified several bone-related genes and elucidated their regulatory roles in the development of IBs.

Abstract

Intermuscular bones (IBs) are small spicule-like bones located in the myosepta of basal teleosts, which negatively affect the edibleness and economic value of fish. Blunt snout bream (Megalobrama amblycephala, with epineural and epipleural IBs) and tilapia (Oreochromis niloticus, without epineural and epipleural IBs) are two major aquaculture species and ideal models for studying the formation mechanisms of fish IBs. Here, we compared myosepta development between M. amblycephala and O. niloticus, based on histological analysis, transcriptome profiling, and expression analysis of bone-related genes. The histological results showed that dye condensation began to appear in the myosepta 20 days post hatching (dph) in M. amblycephala, and IBs could be clearly observed 50 dph in the myosepta, based on different staining methods. However, in O. niloticus, dye condensation was not observed in the myosepta from 10 to 60 dph. Differentially expressed genes (DEGs) at different developmental stages were screened by comparing the transcriptomes of M. amblycephala and O. niloticus, and KEGG analysis demonstrated that these DEGs were enriched in many bone-related pathways, such as focal adhesion, calcium, and Wnt signaling pathways. Quantitative PCR was performed to further compare the expression levels of some bone-related genes (scxa, scxb, runx2a, runx2b, bgp, sp7, col1a2, entpd5a, entpd5b, phex, alpl, and fgf23). All the tested genes (except for alpl) exhibited higher expression levels in M. amblycephala than in O. niloticus. A comparison of the dorsal and abdominal muscle tissues between the two species also revealed significant expression differences for most of the tested genes. The results suggest that scxa, scxb, runx2a, runx2b, entpd5a, col1a2, and bgp may play important roles in IB development. Our findings provide some insights into the molecular mechanisms of IB formation, as well as clues for further functional analysis of the identified genes to better understand the development of IBs.

Transcriptome Analysis Reveals Sexual Disparities between Olfactory and Immune Gene Expression in the Olfactory Epithelium of Megalobrama amblycephala

The olfactory organ is an important chemoreceptor in vertebrates. However, the sexual disparities in gene expression patterns in the olfactory organ in fish remain unstudied. Here, we conducted a transcriptome analysis of the olfactory epithelium (OE) of male and female blunt snout bream (Megalobrama amblycephala) to identify the differences. The histological analysis showed that there were 22 leaf-like olfactory lamellaes on one side of the OE of the adult blunt snout bream. The sensory area of OE is enriched with ciliated receptor cells and microvilli receptor cells. The transcriptome analysis showed that only 10 out of 336 olfactory receptor genes (224 ORs, 5 V1Rs, 55 V2Rs, and 52 TAARs) exhibited significant expression differences between males and females, and most of the differentially expressed genes were related to the immune system. We also validated these results using qPCR: 10 OR genes and 6 immunity-related genes significantly differed between males and females. The FISH analysis results indicated that the ORs were mainly expressed at the edge of the olfactory lamellae. Collectively, our study reveals that gender is not an important factor influencing the expression of olfactory receptors, but the expression of immune genes varies greatly between the genders in blunt snout bream.

Identification and characterization of circRNAs in the liver of blunt snout bream (Megalobrama amblycephala) infected with Aeromonas hydrophila

Circular RNAs (circRNAs), a class of non-coding RNAs, play an important role in regulating various biological processes. In the present study, circRNAs from the Megalobrama amblycephala liver were identified at five different time points post Aeromonas hydrophila using RNA-seq technology. A total of 250 circRNAs were identified, of which 106 were differentially expressed (DE) in ten pairwise comparisons. GO and KEGG analyses showed that the parental genes of DE circRNAs were enriched in phagocytosis, complement and coagulation cascades, and Fc gamma R-mediated phagocytosis pathways. According to ceRNA hypothesis, the interaction network of circRNAs, miRNAs and mRNAs was constructed. Moreover, WGCNA was conducted, and five specific modules significantly related to bacterial infection were identified. All the above results reveal the important role of circRNAs in immune response, which enriches the information of circRNAs in teleost, and helps to understand the immune response mechanism of M. amblycephala to A. hydrophila.

Integrated analysis of lncRNA and mRNA in liver of Megalobrama amblycephala post Aeromonas hydrophila infection

Background

As non-coding RNA molecules of more than 200 bp in length, long non-coding RNAs (lncRNAs) play a variety of roles in biological processes, including regulating the immune responses to bacterial infections. In recent years, there have been many in-depth studies on mammalian lncRNAs, but the relevant studies in fish are very limited. Meanwhile, since lncRNAs are not conserved among species, it is difficult to apply the existing results directly to unstudied species.

Results

To obtain the information of lncRNAs in Megalobrama amblycephala, one of the most economically important freshwater fish in China, also to better understand the biological significance of lncRNAs in the immunity system, the fish liver at 0, 4, 12, 24, and 72 h post Aeromonas hydrophila infection (hpi) were obtained for lncRNA-sequencing (lncRNA-seq). A total of 14,849 lncRNAs were identified, and 2196 lncRNAs showed significant differences at different time points post A. hydrophila infection. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that the target genes of the differentially expressed lncRNAs were enriched in several pathways related to immune such as apoptosis, inflammation, and immune response. Time-specific modules were then identified, using weighted correlation network analysis (WGCNA), and 28 modules significantly correlated with different time point after infection were found. Furthermore, four immune-related genes and six lncRNAs in the time-specific modules were subsequently verified by RT-qPCR.

Conclusions

The above findings reveal the discovery of widespread differentially expressed lncRNAs in the M. amblycephala liver post A. hydrophila infection, suggesting that lncRNAs might participate in the regulation of host response to bacterial infection, enriching the information of lncRNAs in teleost and providing a resources basis for further studies on the immune function of lncRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07969-5.

Chromosome-level assembly of the Hypophthalmichthys molitrix (Cypriniformes: Cyprinidae) genome provides insights into its ecological adaptation

Hypophthalmichthys molitrix (silver carp) is phytoplanktivorous and is an economically and ecologically important fish species. As a well-known invasive species, a number of factors associated with the ecological adaptations of this species are largely unknown. Here, we present a chromosomal-level assembly of the species based on the PacBio Sequel II platform and Hi-C scaffolding technology. Based on the high-quality genome sequences and previous genome sequencing projects, a number of genes that were probably subject to positive selection reside in the genome of H. molitrix, and the last common ancestors of H. molitrix and H. nobilis were identified. Some of these genes may partially explain the mechanisms of H. molitrix for surviving damage due to toxic algae. Demographic history estimation suggests that the effective population size (EPS) of the species may have constantly increased along with the uplift of the Qinghai-Tibet Plateau, started to decline when quaternary glaciation started, and further declined during the Younger Dryas Period. Moreover, the introgression from H. nobilis to H. molitrix in North America was corroborated based on the whole-genome sequencing data, and the proportion of introgressed regions was estimated to be approximately 5.8%. Based on the high-quality assembly, the possible mechanisms by which H. molitrix adapts to its endemic and invaded locations were profiled.

Molecular mechanisms of intermuscular bone development in fish: a review

Intermuscular bones (IBs) are slender linear bones embedded in muscle, which ossify from tendons through a process of intramembranous ossification, and only exist in basal teleosts. IBs are essential for fish swimming, but they present a choking risk during human consumption, especially in children, which can lead to commercial risks that have a negative impact on the aquaculture of these fish. In this review, we discuss the morphogenesis and functions of IBs, including their underlying molecular mechanisms, as well as the advantages and disadvantages of different methods for IB studies and techniques for breeding and generating IB-free fish lines. This review reveals that the many key genes involved in tendon development, osteoblast differentiation, and bone formation, e.g., scxa, msxC, sost, twist, bmps, and osterix, also play roles in IB development. Thus, this paper provides useful information for the breeding of new fish strains without IBs via genome editing and artificial selection.

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.

Genome of the butterfly hillstream loach provides insights into adaptations to torrential mountain stream life

Butterfly hillstream loach (Beaufortia kweichowensis), a benthic fish in the torrential mountain streams, possesses a totally flat ventrum, flattened craniofacial and body skeletons, and enlarged paired fins covered by substantially small keratinous structures. However, little is known about the genetic basis of these specialized morphological adaptations. Here we present a 448.52-Mb genome assembly with contig N50 length of 5.53 Mb by integrating Illumina short-read sequencing, Nanopore long-read sequencing and HiC-based chromatin map. Demographic history reconstruction of the butterfly hillstream loach reveals that the population dynamics is correlated with the different stages of uplifting of the Tibetan Plateau. Comparative genomic analysis finds evidence of six keratin genes in butterfly hillstream loach evolving under positive selection. Within these genes, two keratin genes exhibit species-specific and divergent amino acid changes, suggesting a role in the formation of the unculi. Additionally, a series of positively selected genes, rapid evolving genes, specific variant genes and expanded gene families are found, including genes related to Hedgehog, Notch and BMP pathways, which may be involved in craniofacial development. These findings may have important implications for understanding the genetic basis of phenotypic adaptation to torrential mountain stream life.

A chromosome-level genome assembly of the yellowfin seabream (Acanthopagrus latus; Hottuyn, 1782) provides insights into its osmoregulation and sex reversal

The yellowfin seabream Acanthopagrus latus is the economically most important Sparidae fish in the northern South China Sea. As euryhaline fish, they are perfect model for investigating osmoregulatory mechanisms in teleosts. Moreover, the reproductive biology of hermaphrodites has long been intriguing; however, little information is known about the molecular pathways underlying their sex change. Here, we report a chromosome level reference genome of A. latus generated by employing the PacBio single molecule sequencing technique (SMRT) and high-throughput chromosome conformation capture (Hi-C) technologies. The draft genome of yellowfin seabream was 806 Mb, with 732 Mb scaffolds anchored on 24 chromosomes. The contig N50 and scaffold N50 were 2.6 Mb and 30.17 Mb, respectively. The assembly is of high integrity and includes 92.23% universal single-copy orthologues based on benchmarking universal single-copy orthologs (BUSCO) analysis. A total of 19,631 protein-coding genes were functionally annotated in the reference genome. Moreover, ARRDC3 and GSTA gene families which related to osmoregulation underwent an extensive expansion in two euryhaline sparids fish genomes compared to other teleost genomes. Moreover, integrating sex-specific transcriptome analyses, several genes related to the transforming growth factor beta (TGF-β) signalling pathway involved in sex differentiation and development. This genomic resource will not only be valuable for studying the osmoregulatory mechanisms in estuarine fish and sex determination in hermaphrodite vertebrate species, but also provide useful genomic tools for facilitating breeding of the yellowfin seabream.

Transcriptome profiling towards understanding of the morphogenesis in the scale development of blunt snout bream (Megalobrama amblycephala)

Skin appendages in vertebrates have individual morphological differences, but share the same evolutionary origin. In this study, we used Megalobrama amblycephala as a fish model to study the developmental regulation mechanism of a common skin appendage in fish: scales. By combining in-toto live imaging method and transcriptomic analysis during the scale development, we elucidated core features of scale patterning containing three distinct regions and experiencing four stages. Differentially expressed genes in skin tissues at the initial site before and after scale development were analyzed and some key regulatory genes (Wnt3, Wnt6, Fgf8, Fgf10, Fgf16, Fgfr1a, Ihhb and BMP6) which are crucial for scale morphogenesis were selected. This study provides a strong reference for further exploration of the function of genes related to the molecular regulation mechanism of scale development in M. amblycephala, as well as in other fishes.

Genome-Wide Integrated Analysis Revealed Functions of lncRNA-miRNA-mRNA Interaction in Growth of Intermuscular Bones in Megalobrama amblycephala

Intermuscular bone (IB) occurs in the myosepta of teleosts. Its existence has an adverse influence on the edible and economic value of fish, especially for aquaculture species belonging to Cypriniformes. The growth mechanism of IBs is quite lacking. In this study, we firstly used single molecular real-time sequencing (SMRT) technology to improve the draft genome annotation and full characterization of the transcriptome for one typical aquaculture species, blunt snout bream (Megalobrama amblycephala). The long non-coding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) expression profiles in two IB growth stages (1 and 3 years old) were compared through transcriptome and degradome analyses. A total of 126 miRNAs, 403 mRNAs, and 353 lncRNAs were found to be differentially expressed between the two stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the significantly upregulated map2k6 and cytc in the MAPK/p53 signaling pathway and the significantly downregulated lama3 and thbs4b in the extracellular matrix (ECM)–receptor pathway may play a key regulatory role in IB growth. Bioinformatics analysis subsequently revealed 14 competing endogenous RNA (ceRNA) pairs related to the growth of IBs, consisting of 10 lncRNAs, 7 miRNAs, and 10 mRNAs. Of these, dre-miR-24b-3p and dre-miR-193b-3p are core regulatory factors interacting with four lncRNAs and three mRNAs, the interaction mechanism of which was also revealed by subsequent experiments at the cellular level. In conclusion, our data showed that IBs had higher activity of cell apoptosis and lower mineralization activity in IB_III compared to IB_I via interaction of MAPK/p53 and ECM–receptor signaling pathways. The downregulated zip1 interacted with miR-24a-3p and lnc017705, decreased osteoblast differentiation and Ca²⁺ deposition in the IB_III stage. Our identified functional mRNAs, lncRNAs, and miRNAs provide a data basis for in-depth elucidation of the growth mechanism of teleost IB.

Whole genome sequencing of silver carp (Hypophthalmichthys molitrix) and bighead carp (Hypophthalmichthys nobilis) provide novel insights into their evolution and speciation

The edible silver carp (Hypophthalmichthys molitrix) and bighead carp (H. nobilis), which are two of the "Four Domesticated Fish" of China, are cultivated intensively worldwide. Here, we constructed 837- and 845-Mb draft genome assemblies for the silver carp and the bighead carp, respectively, including 24,571 and 24,229 annotated protein-coding genes. Genetic maps, anchoring 71.7% and 83.8% of all scaffolds, were obtained for the silver and bighead carp, respectively. Phylogenetic analysis showed that the bighead carp formed a clade with the silver carp, with an estimated divergence time of 3.6 million years ago; the time of divergence between the silver carp and zebrafish was 50.7 million years ago. An East Asian cyprinid genome-specific chromosome fusion took place ~9.2 million years after this clade diverged from the clade containing the common carp and Sinocyclocheilus. KEGG and GO analyses indicated that the expanded gene families in the silver and bighead carp were associated with diseases, the immune system and environmental adaptations. Genomic regions differentiating the silver and bighead carp populations were detected based on the whole-genome sequences of 42 individuals. Genes associated with the divergent regions were associated with reproductive system development and the development of primary female sexual characteristics. Thus, our results provided a novel systematic genomic analysis of the East Asian cyprinids, as well as the evolution and speciation of the silver carp and bighead carp.

Genome-wide analysis of intermuscular bone development reveals changes of key genes expression and signaling pathways in blunt snout bream (Megalobrama amblycephala)

Intermuscular bone (IB) is a hard-boned spicule exist in lower teleost, which brings a lot of detrimental effects on palatability and economic value of blunt snout bream (Megalobrama amblycephala). Masson trichrome staining for ossific IB indicated that some osteoblasts appeared at the edge of the bone matrix and a few osteocytes are present in the center of the mineralized bone matrix. By comparing the orthologous gene families of fish with IBs and without IBs, we screened the key signaling pathways associated with IB formation. Furthermore, the transcriptomic data demonstrated the functional importance of these gene families. The candidate genes involved in chondrocyte development were highly expressed in stage 1 compared with stage 2 and stage 3, suggesting that the development process of IB might mainly involve in intramembranous ossification. Our research reveals the molecular mechanism of IBs formation, and provides molecular evidence for the further study on intermuscular boneless stains.

TRAF3 of blunt snout bream participates in host innate immune response to pathogenic bacteria via NF-κB signaling pathway

Tumor necrosis factor receptor-associated factor 3 (TRAF3) is a multifunctional adaptor protein primarily involved in both bacterial defense and antiviral immunity in living organisms. However, the knowledge on TRAF3 in blunt snout bream (Megalobrama amblycephala), a freshwater fish with economic values, remained unclear. In the present study, we identified and characterized successfully Traf3 gene from M. amblycephala (maTraf3). The maTraf3 cDNA contained a 1722 bp open reading frame that encoded a protein of 573 amino acid residues. The deduced amino acid sequence comprised of a RING finger domain, two zinc finger motifs, a coiled-coil region and a MATH domain. Analysis of the transcriptional patterns of maTraf3 revealed that it was ubiquitously distributed in various tissues tested from M. amblycephala, with the abundance of expression in spleen and muscle. Following a challenge with Aeromonas hydrophila and lipopolysaccharide stimulation, the expression of maTraf3 was strongly enhanced at different time points in vitro and in vivo. MaTRAF3 was identified as a cytosolic protein and suggested to form aggregates or be associated with vesicles scattering in the cytoplasm. NF-κB transcription was activated by maTraf3 in reporter assay. The overexpression of maTraf3 produced high levels of pro-inflammatory cytokines such as IL-1β, IL-6, IL-8 and TNF-α, implying its immune-regulatory role in M. amblycephala. Taken together, our results obtained in this study demonstrated the crucial role of maTraf3 in mediating host innate immune response to pathogen invasion via NF-κB signaling pathway, which might indicate a novel therapeutic approach to combat bacterial infection in fish.

Insights into the Synthesis, Secretion and Curing of Barnacle Cyprid Adhesive via Transcriptomic and Proteomic Analyses of the Cement Gland

Barnacles represent one of the model organisms used for antifouling research, however, knowledge regarding the molecular mechanisms underlying barnacle cyprid cementation is relatively scarce. Here, RNA-seq was used to obtain the transcriptomes of the cement glands where adhesive is generated and the remaining carcasses of Megabalanus volcano cyprids. Comparative transcriptomic analysis identified 9060 differentially expressed genes, with 4383 upregulated in the cement glands. Four cement proteins, named Mvcp113k, Mvcp130k, Mvcp52k and Mvlcp1-122k, were detected in the cement glands. The salivary secretion pathway was significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the differentially expressed genes, implying that the secretion of cyprid adhesive might be analogous to that of saliva. Lysyl oxidase had a higher expression level in the cement glands and was speculated to function in the curing of cyprid adhesive. Furthermore, the KEGG enrichment analysis of the 352 proteins identified in the cement gland proteome partially confirmed the comparative transcriptomic results. These results present insights into the molecular mechanisms underlying the synthesis, secretion and curing of barnacle cyprid adhesive and provide potential molecular targets for the development of environmentally friendly antifouling compounds.

The Hybrid Genome of a New Goldfish-Like Fish Lineage Provides Insights Into the Origin of the Goldfish

Distant hybridization leads to obvious changes in genotypes and phenotypes, giving rise to species with novel capabilities. However, the fusion of distinct genomes also polymerizes the DNA or gene variations that occur during the course of evolution. Knowledge of the early stages of post-hybridization evolution is particularly important. Here, we investigated the full-length (FL) transcriptomes and the sequences resulting from the genome resequencing of the red crucian carp-like homodiploid fish (RCC-L) and goldfish-like homodiploid fish (GF-L) derived from the interspecific hybridization of koi carp (KOC) and blunt snout bream (BSB) to provide molecular evidence for the hybrid origin of the goldfish (GF). We compared the orthologous genes in the transcriptomes of RCC-L and GF-L with those of KOC and BSB. We also mapped the orthologous genes to the common carp (CC) and BSB genomes and classified them into eight gene patterns in three categories (chimaera, mutant, and biparental origin genes). The results showed that 48.20% and 46.50% of the genes were chimaera and that 3.70% and 8.30% of the genes were mutations of orthologous genes in RCC-L and GF-L, respectively. In RCC-L and GF-L, 63.70% and 68.20% of the genetic materials were from KOC, and 12.30% and 11.90% of the genetic materials were from BSB. The sequences from the genome resequencing of RCC-L and GF-L were mapped to the genome sequences of CC and BSB, revealing that the similarities of both RCC-L and GF-L to the CC genome (92.52%, 90.18%) were obviously higher than to the BSB genome (50.33%, 49.18%), supporting the suggestion that the genomes of both RCC-L and GF-L were mainly inherited from KOC but had some DNA fragments from BSB. Overall, our results provide molecular biological evidence for the hybrid origin of red crucian carp (RCC) and GF.

Expansion of sweet taste receptor genes in grass carp (Ctenopharyngodon idellus) coincided with vegetarian adaptation

BACKGROUND

Taste is fundamental to diet selection in vertebrates. Genetic basis of sweet taste receptor in the shaping of food habits has been extensively studied in mammals and birds, but scarcely studied in fishes. Grass carp is an excellent model for studying vegetarian adaptation, as it exhibits food habit transition from carnivory to herbivory.

RESULTS

We identified six sweet taste receptors (gcT1R2A-F) in grass carp. The four gcT1R2s (gcT1R2C-F) have been suggested to be evolved from and paralogous to the two original gcT1R2s (gcT1R2A and gcT1R2B). All gcT1R2s were expressed in taste organs and mediated glucose-, fructose- or arginine-induced intracellular calcium signaling, revealing they were functional. In addition, grass carp was performed to prefer fructose to glucose under a behavioral experiment. Parallelly, compared with gcT1R2A-F/gcT1R3 co-transfected cells, gcT1R2C-F/gcT1R3 co-transfected cells showed a higher response to plant-specific fructose. Moreover, food habit transition from carnivory to herbivory in grass carp was accompanied by increased gene expression of certain gcT1R2s.

CONCLUSIONS

We suggested that the gene expansion of T1R2s in grass carp was an adaptive strategy to accommodate the change in food environment. Moreover, the selected gene expression of gcT1R2s might drive the food habit transition from carnivory to herbivory in grass carp. This study provided some evolutional and physiological clues for the formation of herbivory in grass carp.

Functional conservation and divergence of duplicated the suppressor of cytokine signaling 1 in blunt snout bream (Megalobrama amblycephala)

The suppressor of cytokine signaling 1 (SOCS1) is an essential feedback regulator extensively involved in many different cytokine signaling pathways, such as regulation of the immune system and growth of organism. However, the molecular and functional information on socs1 genes in freshwater fish is unclear. In the present paper, we identified and characterized the full-length closely related but distinct socs1 genes (socs 1a and -1b) in blunt snout bream (Megalobrama amblycephala). The bioinformatic analysis results showed that duplicated socs1s shared majority conserved motifs with other vertebrates. Both socs1a and -1b mRNAs were detected throughout embryogenesis, and gradually increase and then constantly expressed after 16 hpf. Whole-mount in situ hybridization demonstrated that socs1a and socs1b mRNAs were detected in the brain at 12hpf and 24hpf, and in the notochord and brain at 36hpf. In adult fish, the socs1a mRNA were strongly expressed in the heart, eye, kidney, spleen and gonad, but were found to be relatively low in the intestine and liver. On the other hand, the expression of socs1b mRNA was significantly high in the muscle, eye and spleen, and relatively low in the intestine, liver, skin and heart. The results of hGH treatment experiment showed that socs1a and 1b mRNAs were upregulated markedly in the kidney, muscle and liver. Overexpression of socs1s significantly inhibit the GH and JAK/STAT factor stat3 and the inhibitory effect of SOCS1s on GH may be involved in JAK-STAT signaling pathway. These results indicate that SOCS1 plays an important role in regulating growth and development.

Isolation and expression of four Megalobrama amblycephala toll-like receptor genes in response to a bacterial infection

Toll-like receptors (TLRs) are a category of pattern recognition receptors (PRRs), which recognize pathogen associated molecular patterns (PAMPs) and participate in the immune responses. We identified tlr5a, tlr5b, tlr9 and tlr21 from the genome of blunt snout bream (Megalobrama amblycephala). All four tlrs were constitutively expressed in all examined tissues. After an immune bacterial challenge with Aeromonas hydrophila, their expressionwas up-regulated in lymphoid organs and tissues. Recombinant eukaryotic plasmid pEGFP-N1 was transfected into the common carp (Cyprinus carpio) EPC (epithelioma papulosum cyprini) cells for the purpose of subcellular localization. pcDNA3.1(+) recombinant eukaryotic plasmid was used to investigate the effects of overexpression of tlrs on the expression of downstream interferon-associated immune factors. The four Tlrs were distributed in the cytoplasm of transfected cells and appeared as filamentous or reticular. The expression of irf3, irf7, isg15, mx1, pkr and viperin at 0, 6, 12, 18, 24, 36, 48 and 72 h post-transfection in transfected EPC cells was quantified by qPCR. Overexpression of tlrs upregulated the expression of viperin, isg15, irf3, irf7, mx1 and pkr (in that order of magnitude). We also cloned the following promoters of irfs: Irf1-p, irf2-p, irf6-p, irf7-p, irf8-p and irf9-p. Results of the dual luciferase reporter assay suggested that tlr5a, tlr5b and tlr9 enhanced the activities of irf7-p, while tlr5b enhanced the activities of irf1-p and irf7-p. This suggests that they all play a role in the innate immunity. The experiments also indicated that TLRs activate irf3 or irf7 signaling to induce IFN secretion and subsequent upregulation of IFN-stimulated genes. These results indicate that tlrs and irfs play an important immune role in response to A. hydrophila infection in blunt snout bream, and pave the way for further studies of immune mechanisms mediated by TLRs in fish.

Expansion of vomeronasal receptor genes (OlfC) in the evolution of fright reaction in Ostariophysan fishes

Ostariophysans are the most diverse group of freshwater fishes and feature a pheromone-elicited fright reaction. However, the genetic basis of fright reaction is unclear. Here, we compared vomeronasal type 2 receptor-like (OlfC) genes from fishes having and lacking fright reaction, to provide insight into evolution of pheromonal olfaction in fishes. We found OlfC genes expanded remarkably in ostariophysans having fright reaction compared with fishes lacking fright reaction. Phylogenetic analysis indicates OlfC subfamily 9 expanded specifically in ostariophysans having fright reaction. Principle component and phylogenetic logistic regression analysis partitioned fishes by ecotype (having or lacking fright reaction) and identified OlfC subfamily 9 as being an important factor for fright reaction. Expression levels of expanded OlfC subfamily genes after fright reaction in zebrafish changed more than did genes that had not expanded. Furthermore, evidence of positive selection was found in the expanded OlfC proteins in ostariophysan fishes having fright reaction. These results provide new insight into the genetic basis of fright reaction in ostariophysan fish and will enable future research into the mechanism of action of OlfC proteins.

Molecular characterization, expression and antimicrobial activity of complement factor D in Megalobrama amblycephala

Complement factor D (Df) is a serine protease, which can activate the alternative pathway by cleaving complement factor B, and involves in the innate defense against pathogens infection in teleost. In this study, we cloned, characterized the Df gene from blunt snout bream (Megalobrama amblycephala) (Mamdf), and examined its expression pattern and antimicrobial activity. The open reading frame (ORF) of Mamdf was 753 bp, encoding 250 amino acids with a molecular mass of 27.2 kDa. Mamdf consisted of a single serine protease trypsin superfamily domain, 3 substrate binding sites and 3 active sites, but no potential N-glycosylation site. Pairwise alignment showed that Mamdf shared the highest identity (94%) with grass carp (Ctenopharyngodon idellus). Phylogenetic analysis indicated that Mamdf and other vertebrate Df had a common ancestral origin. Mamdf structured with 4 introns and 5 exons. The Mamdf mRNA expressed relatively high at the intestine appearance stage during early development and constitutively expressed in various tissues with the highest expression in the kidney in healthy adults. After challenged with Aeromonas hydrophila, significant changes of Mamdf at both mRNA and protein levels in the kidney, spleen, liver and head-kidney were observed. The recombinant Mamdf protein showed antimicrobial activity against both gram-positive bacteria and gram-negative bacteria. The above results suggested the immune function of Mamdf, and would benefit further detailed Df function research in the immune process in teleost.