Chromosome-level genome assembly of the yellow-cheek carp Elopichthys bambusa.

Yellow-cheek carp (Elopichthys bambusa) is a typical large and ferocious carnivorous fish endemic to East Asia, with high growth rate, nutritional value and economic value. In this study, a chromosome-level genome of yellow-cheek carp was generated by combining PacBio reads, Illumina reads and Hi-C data. The genome size is 827.63 Mb with a scaffold N50 size of 33.65 Mb, and 99.51% (823.61 Mb) of the assembled sequences were anchored to 24 pseudo-chromosomes. The genome is predicted to contain 24,153 protein-coding genes, with 95.54% having functional annotations. Repeat elements account for approximately 55.17% of the genomic landscape. The completeness of yellow-cheek carp genome assembly is highlighted by a BUSCO score of 98.4%. This genome will help us understand the genetic diversity of yellow-cheek carp and facilitate its conservation planning.

Combining genome-wide association study and transcriptome analysis to identify molecular markers and genetic basis of population-asynchronous ovarian development in Coilia nasus.

Coilia nasus, a migratory fish species found in the middle and lower reaches of the Yangtze River and along offshore areas of China, possesses considerable aquacultural and economic potential. However, the species faces challenges due to significant variation in the gonadal development rate among females, resulting in inconsistent ovarian maturation times at the population level, an extended reproductive period, and limitations on fish growth rate due to ovarian prematurity. In the present study, we combined genome-wide association study (GWAS) and comparative transcriptome analysis to investigate the potential single nucleotide polymorphisms (SNPs) and candidate genes associated with population-asynchronous ovarian development in C. nasus. Genotyping of the female population based on whole-genome resequencing yielded 2 120 695 high-quality SNPs, 39 of which were suggestively associated with ovarian development. Of note, a significant SNP peak on LG21 containing 30 suggestively associated SNPs was identified, with cpne5a determined as the causal gene of the peak. Therefore, single-marker and haplotype association analyses were performed on cpne5a, revealing four genetic markers ( P<0.05) and seven haplotypes (r 2>0.9) significantly associated with the phenotype. Comparative transcriptome analysis of precociously and normally maturing individuals screened out 29 and 426 overlapping differentially expressed genes in the brain and ovary, respectively, between individuals of different body sizes. Integrating the GWAS and transcriptome analysis results, this study identified genes and pathways related to hypothalamic-pituitary-gonadal axis hormone secretion, extracellular matrix, angiogenesis, and gap junctions involved in population-asynchronous ovarian development. The insights gained from this study provide a basis for a deeper understanding of the molecular mechanisms underlying ovarian development in fish and may facilitate the genetic breeding of C. nasus strains exhibiting population-synchronous ovarian development in the future.

Molecular Characterization and Expression Changes of the bcl2l13 Gene in Response to Hypoxia in Megalobrama amblycephala.

Hypoxia is a unique environmental stress, which not only reflects the insufficient oxygen supply of cells and tissues, but also occurs in various physiological and pathological environments. Mitophagy as a selective autophagy can recover and utilize damaged organelles and misfolded proteins to ensure normal cell functions and promote cell survival. Bcl2l13 (B-cell lymphoma-2 like 13) is reported to induce mitophagy as a functional mammalian homolog of Atg32. However, the function of the bcl2l13 gene is still unclear in fish. Here the sequence and structure of the bcl2l13 gene in Megalobrama amblycephala were identified and showed that bcl2l13 contained an open reading frame (ORF) of 1458 bp for encoding 485 aa. Amino acid sequence analysis indicated that Bcl2l13, as a typical anti-apoptotic protein of the Bcl2 family, contained four BH domains, one BHNo domain, and one TM domain. Further study showed that Bcl2l13 was mainly located in the mitochondria, while its localization was changed within the whole cell after the TM domain was deleted. Real-time PCR analysis revealed that bcl2l13 showed higher expression levels in early embryos. After hypoxia treatment, the mRNA levels of the bcl2l13 and autophagy-related genes were significantly up-regulated in most detected tissues, and the bcl2l13 transcription was regulated by Hif-1α mediated pathway. Additionally, the transcription activity of the bcl2l13 promoter was further analyzed using luciferase reporter assays and showed the highest activity in the promoter region from -475 to +111. These results indicated that bcl2l13 may play important roles in embryogenesis and hypoxia mediated autophagy in fish.

Muscle Fiber Characteristics and Transcriptome Analysis in Slow- and Fast-Growing Megalobrama amblycephala.

Growth is an important trait in aquaculture that is influenced by various factors, among which genetic regulation plays a crucial role. Megalobrama amblycephala, one of the most important freshwater species in China, exhibits wide variations in body mass among individuals of the same age within the same pool. But the molecular mechanisms underlying wide variation in body mass remain unclear. Here, we performed muscle histological and transcriptome analysis of muscle tissues from Fast-Growing (FG) and Slow-Growing (SG) M. amblycephala at the age of 4 months old (4 mo) and 10 months old (10 mo) to elucidate its muscle development and growth mechanism. The muscle histological analysis showed smaller diameter and higher total number of muscle fibers in FG compared to SG at 4 mo, while larger diameter and total number of muscle fibers were detected in FG at 10 mo. The transcriptome analysis of muscle tissue detected 1171 differentially expressed genes (DEGs) between FG and SG at 4 mo, and 718 DEGs between FG and SG at 10 mo. Furthermore, 44 DEGs were consistently up-regulated in FG at both 4 mo and 10 mo. Up-regulated DEGs in FG at 4 mo were mainly enriched in the pathways related to cell proliferation, while down-regulated DEGs were significantly enriched in cell fusion and muscle contraction. Up-regulated DEGs in FG at 10 mo were mainly enriched in the pathways related to cell proliferation and protein synthesis. Therefore, these results provide novel insights into the molecular mechanism of M. amblycephala muscle growth at different stages, and will be of great guiding significance to promote the fast growth of M. amblycephala.

Effect of runx2b deficiency in intermuscular bones on the regulatory network of lncRNA-miRNA-mRNA.

Intermuscular bones (IBs) are mineralized spicules that negatively impact the quality and value of fish products. Runx2b is a crucial modulator in promoting bone formation through regulating osteoblast differentiation. Previous studies suggested that loss of runx2b gene completely inhibited IBs formation in zebrafish. However, how the whole transcriptome, including mRNA and non-coding RNA (ncRNA), affects the IBs development in runx2b-/- zebrafish are not known. The aim of this study was to identify the regulatory networks of differentially expressed (DE) lncRNAs, miRNAs, and mRNAs in zebrafish with and without IBs (runx2b+/+ fish and runx2b-/- fish) utilizing high-throughput sequencing techniques. All together there are 1051 mRNAs, 456 lncRNAs, and 18 miRNAs differentially expressed were found between these two strains. The analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) has highlighted significant pathways linked to the development of IBs, specifically the TGF-beta and Wnt signaling pathways, and a number of genes concentrated on these two signaling pathways related to the formation of IBs. Further, 1989 competing endogenous RNA (ceRNA) networks were created according to the correlation among mRNAs, miRNAs and lncRNAs. The ceRNA networks results revealed 52 ceRNA pairs related to the IBs formation, consisting of 52 mRNAs, 37 lncRNAs, and 6 miRNAs. Of these, we found that dre-miR-2189 was the key element of ceRNA pairs, interacting with 19 mRNAs and 11 lncRNAs, and MSTRG.13175.1 could regulate sp7 expression by interacting with dre-miR-2189 to function in osteogenic differentiation. Subsequent experiments at the cellular level also revealed the interaction mechanism. The outcomes indicated a crucial role of miRNAs and lncRNAs in the development of fish IBs, which offer new views into the functions of ncRNAs involved in IBs formation.

The genetic basis of the leafy seadragon's unique camouflage morphology and avenues for its efficient conservation derived from habitat modeling.

The leafy seadragon certainly is among evolution's most "beautiful and wonderful" species aptly named for its extraordinary camouflage mimicking its coastal seaweed habitat. However, limited information is known about the genetic basis of its phenotypes and conspicuous camouflage. Here, we revealed genomic signatures of rapid evolution and positive selection in core genes related to its camouflage, which allowed us to predict population dynamics for this species. Comparative genomic analysis revealed that seadragons have the smallest olfactory repertoires among all ray-finned fishes, suggesting adaptations to the highly specialized habitat. Other positively selected and rapidly evolving genes that serve in bone development and coloration are highly expressed in the leaf-like appendages, supporting a recent adaptive shift in camouflage appendage formation. Knock-out of bmp6 results in dysplastic intermuscular bones with a significantly reduced number in zebrafish, implying its important function in bone formation. Global climate change-induced loss of seagrass beds now severely threatens the continued existence of this enigmatic species. The leafy seadragon has a historically small population size likely due to its specific habitat requirements that further exacerbate its vulnerability to climate change. Therefore, taking climate change-induced range shifts into account while developing future protection strategies.

Genetic Diversity and Population Structure of Coilia nasus Revealed by 2b-RAD Sequencing.

Coilia nasus is a threatened migratory species in the Yangtze River Basin. To reveal the genetic diversity of natural and farmed populations of C. nasus and the status of germplasm resources in the Yangtze River, the genetic diversity and structure of two wild populations (Yezhi Lake: YZ; Poyang Lake: PY) and two farmed populations (Zhenjiang: ZJ; Wuhan: WH) of C. nasus were analyzed using 44,718 SNPs obtained via 2b-RAD sequencing. The results indicate that both the wild and farmed populations had low genetic diversity, and germplasm resources have undergone varying degrees of degradation. Population genetic structure analyses indicated that the four populations may have come from two ancestral groups. Different amounts of gene flow were identified among WH, ZJ, and PY populations, but gene flow among YZ and other populations was low. It is speculated that the river-lake isolation of Yezhi Lake is the main cause of this phenomenon. In conclusion, this study revealed that genetic diversity reduction and germplasm resource degradation had occurred in both wild and farmed C. nasus, suggesting that conservation of its resources is of great urgency. This study provides a theoretical basis for the conservation and rational exploitation of germplasm resources for C. nasus.

Identification of thyroid hormone response genes in the remodeling of dorsal muscle during Microhyla fissipes metamorphosis.

Introduction: Extensive morphological, biochemical, and cellular changes occur during anuran metamorphosis, which is triggered by a single hormone, thyroid hormone (TH). The function of TH is mainly mediated through thyroid receptor (TR) by binding to the specific thyroid response elements (TREs) of direct response genes, in turn regulating the downstream genes in the cascade. The remodeling of dorsal skeletal muscle during anuran metamorphosis provides the perfect model to identify the immediate early and direct response genes that are important during apoptosis, proliferation, and differentiation of the muscle. Methods: In our current study, we performed Illumina sequencing combined with single-molecule real-time (SMRT) sequencing in the dorsal muscle of Microhyla fissipes after TH, cycloheximide (CHX), and TH_CHX treatment. Results and Discussion: We first identified 1,245 differentially expressed transcripts (DETs) after TH exposure, many of which were involved in DNA replication, protein processing in the endoplasmic reticulum, cell cycle, apoptosis, p53 signaling pathway, and protein digestion and absorption. In the comparison of the TH group vs. control group and TH_CHX group vs. CHX group overlapping gene, 39 upregulated and 6 downregulated genes were identified as the TH directly induced genes. Further analysis indicated that AGGTCAnnTnAGGTCA is the optimal target sequence of target genes for TR/RXR heterodimers in M. fissipes. Future investigations on the function and regulation of these genes and pathways should help to reveal the mechanisms governing amphibian dorsal muscle remodeling. These full-length and high-quality transcriptomes in this study also provide an important foundation for future studies in M. fissipes metamorphosis.

DNA methylation mediates gonadal development via regulating the expression levels of cyp19a1a in loach Misgurnus anguillicaudatus.

DNA methylation-mediated transcriptional regulation has been considered to significantly impact some steroidogenic enzyme genes expression. To uncover the roles of DNA methylation on the regulation of aromatase gene expression during gametogenesis in Misgurnus anguillicaudatus, the expression profiles and cellular localization of cyp19a1a and cyp19a1b were analyzed, and the landscape of DNA methylation dynamics was investigated. We found that cyp19a1a was predominantly expressed in granulosa cells of oocytes, while cyp19a1b expression was enriched in radial glial cells of the forebrain. In ovary, cyp19a1a was highly expressed until the vitellogenesis stage. The average methylation levels, especially for two CpG sites within the cAMP response element, were negatively correlated with cyp19a1a expression levels, indicating that methylation could regulate cyp19a1a transcriptional activity by modulating the binding efficiency of cAMP to its response elements. Compared with in ovary, cyp19a1a showed lower expression in testis but was hypermethylated. Cyp19a1b in female brain weakly expressed before the vitellogenesis stage, but significantly elevated at the maturation stage. In both sexes, it maintained high methylation levels in brain despite the obvious fluctuation of the cyp19a1b expression. This study revealed that DNA methylation plays a key role in establishing cyp19a1a spatiotemporal expression patterns and thus mediates gonadal development in teleosts.

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.

Comparative transcriptome profiles of four sexually size dimorphic fish.

Sexual size dimorphism is widespread in fish species. Although sex growth differences in multiple species have been studied successively, the commonalities of regulatory mechanisms across sexually dimorphic species are unknown. In this study, we performed RNA-seq analysis of four representative fish (loach, half-smooth tongue sole, yellow catfish, and Nile tilapia) with significant growth differences between females and males. Clean reads were identified from four fish species, ranging from 45,718,052 to 57,733,120. Following comparison transcriptome analysis, there were 1,132 and 1,108, 1,290 and 1,102, 4,732 and 4,266, 748 and 192 differentially expressed genes (DEGs) in the brain and muscle of loach, half-smooth tongue sole, yellow catfish, and Nile tilapia, respectively. Furthermore, the expression levels were validated by quantitative real-time PCR (qRT-PCR). Comparative transcriptome profiles of four fish described here will provide fundamental information for further studies on the commonalities of sexually size dimorphic fish in regulating growth differences between females and males.

Single-cell transcriptomes and runx2b-/- mutants reveal the genetic signatures of intermuscular bone formation in zebrafish.

Intermuscular bones (IBs) are mineralized spicules, present in the myosepta of many, but not all, teleost species. IBs are often small and sharp, and they consequently limit how the fish can be processed; the IBs may cause injury or trauma if lodged in consumers' throats or mouths, and therefore affect the appeal of the fish to many consumers. The development of IBs in teleosts is still not fully understood and the molecular basis of IB development remains to be established. Here, the characteristics of IB tissue are evaluated based on single-cell transcriptomics in wild-type zebrafish. The analysis defined 18 distinct cell types. Differentiation trajectories showed that IBs are derived from tendons and that a core tendon-osteoblast cell lineage is related to IB formation. In particular, the functions of 10 candidate genes were evaluated via CRISPR-Cas9 mutants. Among those, runx2b-/- mutants completely lost IBs, while swimming performance, growth and bone mineral density were not significantly different from runx2b+/+ zebrafish. Comparative single-cell RNA sequencing (scRNA-seq) analysis in runx2b-/- and runx2b+/+ zebrafish revealed the role of osteoblasts in IB formation. In addition, differentially expressed genes were enriched in the transforming growth factor ß/bone morphogenetic protein (TGF-ß/BMP) pathway after runx2b deletion. This study provides evidence for the crucial role of runx2b regulation in IB formation. Genetic breeding can target runx2b regulation and generate strains of commercial fish species without IBs, which can improve the safe consumption and economic value of many farmed fish species.

Functional differentiation of bmp2a and bmp2b genes in zebrafish.

Bone morphogenetic protein 2 plays an important role in the regulation of osteoblast proliferation and differentiation. Phylogenetic analysis showed that the bmp2 ortholog evolved from the same ancestral gene family in vertebrates and was duplicated in teleost, which were named bmp2a and bmp2b. The results of whole-mount in situ hybridization showed that the expression locations of bmp2a and bmp2b in zebrafish were different in different periods (24 hpf, 48 hpf, 72 hpf), which revealed potential functional differentiation between bmp2a and bmp2b. Phenotypic analysis showed that bmp2a mutations caused partial rib and vertebral deformities in zebrafish, while bmp2b-/- embryos died massively after 12 hpf due to abnormal somite formation. We further explored the expression pattern changes of genes (bmp2a, bmp2b, smad1, fgf4, runx2b, alp) related to skeletal development at different developmental stages (20 dpf, 60 dpf, 90 dpf) in wild-type and bmp2a-/- zebrafish. The results showed that the expression of runx2b in bmp2a-/- was significantly downregulated at three stages and the expression of other genes were significantly downregulated at 90 dpf compared with wild-type zebrafish. The study revealed functional differentiation of bmp2a and bmp2b in zebrafish embryonic and skeletal development.

High-temperature stress will put the thermo-sensitive teleost yellow catfish (Tachysurus fulvidraco) in danger through reducing reproductivity.

Recently, concerns for species that sex differentiation is influenced by temperature in the context of global warming have increased because disrupted operational sex ratios could threaten population maintenance. In contrast, little attention has been given to the reproductive ability of populations that experienced elevated temperatures. In this study, we demonstrated that high temperature (HT) would decrease population size via three different aspects of reproductive ability for the first time. We show that, in a thermo-sensitive teleost yellow catfish, a short period of HT (+3 °C) exposure during the critical period of sex differentiation leads to a different percentage of masculinization of XX genotypic females (1-23%) in wet-lab and natural water bodies. Combining the results of gonadal appearance, histology, sperm parameters, and fertilization rate, we found that XX pseudo-males induced by HT display significantly discounted fertility and reproductive performance compared to XY normal males. We demonstrate that the survival of the XY genotype is lower than XX genotype under environmental stress, including HT, hypoxia, and parasite infection, and the differential survival seems unrelated to male-biased sexual size dimorphism. The mathematical model predicts that the phenotypic female percent will be stabilized at 50% and the population will be sustainably maintained when masculinizing force is less than 0.5, while HT will put the population in danger when the masculinizing force exceeds 0.5. However, when we combine the real-world data of reproductive ability and mathematic model, our results suggest the population size decreases and the long-term survival of the studied species are threatened under the projected pace of increasing temperature. These findings will be useful for understanding the long-term effects of increasing temperature on sex ratio, reproduction and population maintenance in teleost.

Functional Differentiation of BMP7 Genes in Zebrafish: bmp7a for Dorsal-Ventral Pattern and bmp7b for Melanin Synthesis and Eye Development.

Bone morphogenetic protein 7 (BMP7) belongs to the transforming growth factor ß (TGF-ß) family, which not only induces cartilage and bone formation, but also regulates eye development and melanoma tumorigenesis in mammals. In teleosts, BMP7 differentiates into two subtypes, bmp7a and bmp7b, which have clearly differentiated structures. To fully understand the functional differentiation of bmp7a and bmp7b in fish species, we successfully constructed bmp7a and bmp7b gene deletion mutants in zebrafish using CRISPR/Cas9-mediated gene editing technology. Our results showed that bmp7a mutation caused abnormal development of the embryo's dorsal-ventral pattern that led to death; bmp7b mutation induced growth inhibition and increased melanin production in the skin and eye of mutants. Histological analysis revealed that melanin in the retina of the eyes in bmp7b mutants increased, and behavioral observation showed that the vision and sensitivity to food of the mutants were reduced. Transcriptome analysis of the skin and eye tissues showed that the expression changes of wnt7ba and gna14 in bmp7b mutants might promote the increase of melanin. Additionally, the eye transcriptome analysis indicated that changes in the structure of the eyes in bmp7b mutants led to defects in phototransduction, and seven DEGs (rgs9a, rgs9b, rcvrn2, guca1d, grk1b, opn1mw4, and gc2) were identified as key candidate genes that affected the photonic response of the eyes. The study revealed the functional differentiation of bmp7a and bmp7b in teleosts and the first report about the inhibitory effect of bmp7b on melanogenesis may provide useful information for the future research on human melanoma-related diseases.

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

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.

Comparative Transcriptome Analysis Revealed Genes Involved in Sexual and Polyploid Growth Dimorphisms in Loach (Misgurnus anguillicaudatus).

Sexual and polyploidy size dimorphisms are widespread phenomena in fish, but the molecular mechanisms remain unclear. Loach (Misgurnus anguillicaudatus) displays both sexual and polyploid growth dimorphism phenomena, and are therefore ideal models to study these two phenomena. In this study, RNA-seq was used for the first time to explore the differentially expressed genes (DEGs) between both sexes of diploid and tetraploid loaches in four tissues (brain, gonad, liver, and muscle). Results showed that 21,003, 17, and 1 DEGs were identified in gonad, liver, and muscle tissues, respectively, between females and males in both diploids and tetraploids. Regarding the ploidy levels, 4956, 1496, 2187, and 1726 DEGs were identified in the brain, gonad, liver, and muscle tissues, respectively, between tetraploids and diploids of the same sex. When both sexual and polyploid size dimorphisms were considered simultaneously in the four tissues, only 424 DEGs were found in the gonads, indicating that these gonadal DEGs may play an important regulatory role in regulating sexual and polyploid size dimorphisms. Regardless of the sex or ploidy comparison, the significant DEGs involved in glycolysis/gluconeogenesis and oxidative phosphorylation pathways were upregulated in faster-growing individuals, while steroid hormone biosynthesis-related genes and fatty acid degradation and elongation-related genes were downregulated. This suggests that fast-growing loaches (tetraploids, females) have higher energy metabolism levels and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing loaches (diploids, males). Our findings provide an archive for future systematic research on fish sexual and polyploid dimorphisms.

Identification and Characterization of Novel circRNAs Involved in Muscle Growth of Blunt Snout Bream (Megalobrama amblycephala).

Circular RNAs (circRNAs), a novel class of endogenous RNAs, have been recognized to play important roles in the growth of animals. However, the regulatory mechanism of circRNAs on fish muscle growth is still unclear. In this study, we performed whole transcriptome analysis of skeletal muscles from two populations with different growth rates (fast-growing and slow-growing) of blunt snout bream (Megalobrama amblycephala), an important fish species for aquaculture. The selected circRNAs were validated by qPCR and Sanger sequencing. Pairs of circRNA-miRNA-mRNA networks were constructed with the predicted differentially expressed (DE) pairs, which revealed regulatory roles in muscle myogenesis and hypertrophy. As a result, a total of 445 circRNAs were identified, including 42 DE circRNAs between fast-growing (FG) and slow-growing (SG) groups. Many of these DE circRNAs were related with aminoglycan biosynthetic and metabolic processes, cytokinetic processes, and the adherens junction pathway. The functional prediction results showed that novel_circ_0001608 and novel_circ_0002886, competing to bind with dre-miR-153b-5p and dre-miR-124-6-5p, might act as competing endogenous RNAs (ceRNAs) to control MamblycephalaGene14755 (pik3r1) and MamblycephalaGene10444 (apip) level, respectively, thus playing an important regulatory role in muscle growth. Overall, these results will not only help us to further understand the novel RNA transcripts in M. amblycephala, but also provide new clues to investigate the potential mechanism of circRNAs regulating fish growth and muscle development.

Seadragon genome analysis provides insights into its phenotype and sex determination locus.

The iconic phenotype of seadragons includes leaf-like appendages, a toothless tubular mouth, and male pregnancy involving incubation of fertilized eggs on an open "brood patch." We de novo-sequenced male and female genomes of the common seadragon (Phyllopteryx taeniolatus) and its closely related species, the alligator pipefish (Syngnathoides biaculeatus). Transcription profiles from an evolutionary novelty, the leaf-like appendages, show that a set of genes typically involved in fin development have been co-opted as well as an enrichment of transcripts for potential tissue repair and immune defense genes. The zebrafish mutants for scpp5, which is lost in all syngnathids, were found to lack or have deformed pharyngeal teeth, supporting the hypothesis that the loss of scpp5 has contributed to the loss of teeth in syngnathids. A putative sex-determining locus encoding a male-specific amhr2y gene shared by common seadragon and alligator pipefish was identified.