Search

Goldfish (Carassius auratus) have long fascinated evolutionary biologists and geneticists because of their diverse morphological and color variations. Recent genome-wide association studies have provided a clue to uncover genomic basis underlying these phenotypic variations, but the causality between phenotypic and genotypic variations have not yet been confirmed. Here, we edited proposed candidate genes to recreate phenotypic traits and developed a rapid biotechnology approach which combines gene editing with high-efficiency breeding, artificial gynogenesis, and temperature-induced sex reversal to establish homozygous mutants within two generations (approximately eight months). We first verified that low-density lipoprotein receptor-related protein 2B (lrp2aB) is the causal gene for the dragon-eye variation and recreated the dragon-eye phenotype in side-view Pleated-skirt Lion-head goldfish. Subsequently, we demonstrated that the albino phenotype was determined by both homeologs of oculocutaneous albinism type II (oca2), which has subfunctionalized to differentially govern melanogenesis in the goldfish body surface and pupils. Overall, we determined two causal genes for dragon-eye and albino phenotypes, and created four stable homozygous strains and more appealing goldfish with desirable traits. The developed biotechnology approach facilitates precise genetic breeding, which will accelerate re-domestication and recreation of phenotypically desirable goldfish.

Tripartite motif (TRIM) family proteins have come forth as important modulators of innate signaling dependent on of E3 ligase activity. Recently, several human TRIM proteins have been identified as unorthodox RNA-binding proteins by RNA interactome analyses; however, their targets and functions remain largely unknown. FTRCA1 is a crucian carp (Carassius auratus)–specific finTRIM (fish novel TRIM) member and negatively regulates the IFN antiviral response by targeting two retinoic acid–inducible gene-I (RIG-I)–like receptor (RLR) pathway molecules, that is, TANK-binding kinase 1 (TBK1) and IFN regulatory factor 7 (IRF7). In this study, we identify FTRCA1 as an RNA-binding E3 ligase and characterize the contribution of its RNA-binding activity and E3 ligase activity to fish IFN response. Besides targeting TBK1 and IRF7, FTRCA1 downregulates fish IFN response also by targeting stimulator of IFN response cGAMP interactor 1 (STING1). E3 ligase activity is required for full inhibition on the TBK1- and IRF7-mediated IFN response, but partial inhibition on the STING1-mediated IFN response. However, FTRCA1 has a general binding potential to mRNAs in vitro, it selectively binds STING1 and IRF7 mRNAs in vivo to attenuate mRNA levels, and it directly interacts with TBK1 protein to target protein degradation for downregulating the IFN response. Our results present an interesting example of a fish species–specific finTRIM protein that has acquired RNA-binding activity and E3 ligase activity to fine-tune fish IFN response.

Triploids are rare in nature because of difficulties in meiotic and gametogenic processes, especially in vertebrates. The Carassius complex of cyprinid teleosts contains sexual tetraploid crucian carp/goldfish (C. auratus) and unisexual hexaploid gibel carp/Prussian carp (C. gibelio) lineages, providing a valuable model for studying the evolution and maintenance mechanism of unisexual polyploids in vertebrates. Here we sequence the genomes of the two species and assemble their haplotypes, which contain two subgenomes (A and B), to the chromosome level. Sequencing coverage analysis reveals that C. gibelio is an amphitriploid (AAABBB) with two triploid sets of chromosomes; each set is derived from a different ancestor. Resequencing data from different strains of C. gibelio show that unisexual reproduction has been maintained for over 0.82 million years. Comparative genomics show intensive expansion and alterations of meiotic cell cycle-related genes and an oocyte-specific histone variant. Cytological assays indicate that C. gibelio produces unreduced oocytes by an alternative ameiotic pathway; however, sporadic homologous recombination and a high rate of gene conversion also exist in C. gibelio. These genomic changes might have facilitated purging deleterious mutations and maintaining genome stability in this unisexual amphitriploid fish. Overall, the current results provide novel insights into the evolutionary mechanisms of the reproductive success in unisexual polyploid vertebrates.

The small HERC family currently comprises four members (HERC3–6) involved in the regulation of various physiological activities. Little is known about the role of HERCs in IFN response. In this study, we identify a novel fish HERC member, named crucian carp HERC7, as a negative regulator of fish IFN response. Genome-wide search of homologs and comprehensive phylogenetic analyses reveal that the small HERC family, apart from HERC3–6 that have been well-characterized in mammals, contains a novel HERC7 subfamily exclusively in nonmammalian vertebrates. Lineage-specific and even species-specific expansion of HERC7 subfamily in fish indicates that crucian carp HERC7 might be species-specific. In virally infected fish cells, HERC7 is induced by IFN and selectively targets three retinoic acid–inducible gene-I–like receptor signaling factors for degradation to attenuate IFN response by two distinct strategies. Mechanistically, HERC7 delivers mediator of IFN regulatory factor 3 activator and mitochondrial antiviral signaling protein for proteasome-dependent degradation at the protein level and facilitates IFN regulatory factor 7 transcript decay at the mRNA level, thus abrogating cellular IFN induction to promote virus replication. Whereas HERC7 is a putative E3 ligase, the E3 ligase activity is not required for its negative regulatory function. These results demonstrate that the ongoing expansion of the small HERC family generates a novel HERC7 to fine-tune fish IFN antiviral response.

Aquatic animal viruses infect and transmit in aquatic environments, causing serious harm to the aquaculture industry and a variety of wild aquatic animals. How are they affected by environmental factors and do they represent potential threat to mammalian heath or not? Here, the effects of environmental factors (ultraviolet radiation (UV), temperature, pH, and drying) and their threshold on five epidemic aquatic animal viruses infecting amphibians and bony fish, including Rana grylio virus (RGV), Andrias davidianus ranavirus (ADRV), Grass carp reovirus (GCRV), Paralichthys olivaceus rhabdovirus (PORV), and Scophthalmus maximus rhabdovirus (SMRV), were measured and compared in a fish cell line. The examination of virus titers after different treatment in fish cells showed that the two iridoviruses, RGV and ADRV, had a higher tolerance to all of the environmental factors, such as they only had a decay rate of 22–36% when incubated at 37 °C for 7 days. However, the rhabdovirus SMRV was sensitive to all of the factors, with a decay rate of more than 80% in most of the treatments; even a complete inactivation (100%) can be observed after drying treatment. To address the potential threat to mammals, infectivity and limitation factors of the five viruses in Baby hamster kidney fibroblast cells (BHK-21) were tested, which showed that three of the five viruses can replicate at a low temperature, but a high temperature strongly inhibited their infection and none of them could replicate at 37 °C. This study clarified the sensitivity or tolerance of several different types of aquatic animal viruses to the main environmental factors in the aquatic environment and proved that the viruses cannot replicate in mammalian cells at normal physiological temperature.

Unisexual animals are commonly found in some polyploid species complexes, and most of these species have had a long evolutionary history. However, their method for avoiding genomic decay remains unclear. The polyploid Carassius complex naturally comprises the sexual amphidiploid C. auratus (crucian carp or goldfish) (AABB) and the gynogenetic amphitriploid C. gibelio (gibel carp) (AAABBB). Recently, we developed a fertile synthetic amphitetraploid (AAAABBBB) male from C. gibelio by incorporating a C. auratus genome. In this study, we generated novel amphitriploids (AAABBB) by backcrossing the amphitetraploid male with the amphidiploid C. auratus. Whole-genome resequencing revealed the genomic changes, including recombination and independent assortment between homologs of C. gibelio and C. auratus. The fertility, sex determination system, oocyte development, and fertilization behaviors of the novel amphitriploids were investigated. Approximately 80% of the novel amphitriploid females recovered the unisexual gynogenesis ability. Intriguingly, two types of primary oocyte (with and without homolog synapsis) were discovered, and their distinct development fates were observed. Type I oocytes entered apoptosis due to improper synaptonemal complex assembly and incomplete double-strand break repair, whereas subsequent type II oocytes bypassed meiosis through an alternative ameiotic pathway to develop into mature eggs. Moreover, gynogenesis was stabilized in their offspring, and a new array of diverse gynogenetic amphitriploid clones was produced. These revealed genomic changes and detailed cytological data provide comprehensive evidence that changes in ploidy drive unisexual and sexual reproduction transition, thereby resulting in genomic diversity and allowing C. gibelio avoid genomic decay.

Chedrinae fish, which belong to Danionidae, have important ornamental, economic, and scientific value. At present, however, their mitogenomic features are unclear and their phylogenetic relationships remain controversial. In this study, we presented five new Chedrinae mitochondrial genomes (mitogenomes) and analyzed the conserved and variable mitogenomic characteristics of 17 Chedrinae fish. The gene composition and arrangement and secondary structure of transfer RNAs (tRNAs) were highly conserved among the Chedrinae mitogenomes. However, the length of the control region and base composition were variable. Interestingly, the mitogenome of Barilius barila was unusual, with lower A+T content in the first codon of protein-coding genes (PCGs) (47.32% versus average of 54.47%) and distinct pattern of codons per thousand codons (CDspT). Three Chedrinae fish had a long tandem repeat (>291 bp) in the 5'-end of the control region, which may increase their adaptability. In addition, tRNALys had notably larger DHU and TΨC loops than other tRNAs. The phylogenetic trees of the Chedrinae fish suggested that the Barilius genus was not a monophyletic group but could be divided into two main groups based on significant differences in A+T content. This study provides insights into the mitogenomic features and phylogenetic implications of Chedrinae fish, which should benefit their systematics and conservation.

Retinoic acid inducible gene-I (RIG-I)-like receptors (RLRs) are viral RNA sensors that regulate host interferon (IFN)-mediated antiviral signaling. LGP2 (laboratory genetics and physiology 2) lacks the N-terminal caspase activation and recruitment domains (CARDs) responsible for signaling transduction in the other two RLR proteins, RIG-I and melanoma differentiation associated gene-5 (MDA5). How LGP2 regulates IFN signaling is controversial, and inconsistent results have often been obtained in overexpression assays when performed in fish cells and mammalian cells. Here we report that the differential sensitivity of fish cells and mammalian cells to poly(I:C) transfection conceals the function conservation of zebrafish and human LGP2. In fish cells, overexpression of zebrafish or human LGP2 initially activates IFN signaling in a dose-dependent manner, followed by inhibition at a critical threshold of LGP2 expression. A similar trend exists for LGP2-dependent IFN induction in response to stimulation by low and high concentrations of poly(I:C). In contrast, overexpression of zebrafish or human LGP2 alone in mammalian cells does not activate IFN signaling, but co-stimulation with very low or very high concentrations of poly(I:C) shows LGP2-dependent enhancement or inhibition of IFN signaling, respectively. Titration assays show that LGP2 promotes MDA5 signaling in mammalian cells mainly under low concentration of poly(I:C) and inhibits RIG-I/MDA5 signaling mainly under high concentration of poly(I:C). Our results suggest that fish and human LGP2s switch regulatory roles from a positive one to a negative one in increasing concentrations of poly(I:C)-triggered IFN response.

Foxl2 plays conserved central function in ovarian differentiation and maintenance in several fish species. However, its expression pattern and function in fish embryogenesis are still largely unknown. In this study, we first presented a sequential expression pattern of zebrafish foxl2a and foxl2b during embryo development. They were predominantly expressed in the cranial paraxial mesoderm (CPM) and cranial venous vasculature (CVV) during somitogenesis and subsequently expressed in the pharyngeal arches after 48 h post-fertilization (hpf). Then, we compared the brain structures among zebrafish wildtype (WT) and three homozygous foxl2 mutants (foxl2a

Improvement in fish feed conversion efficiency (FCE) is beneficial for sustaining global food fish supplies. Here, we show that a set of polymorphisms at locus of the corticotropin releasing hormone receptor 2 (crhr2), which is involved in hypothalamus-pituitary-interrenal (HPI) axis signaling, is associated with improved FCE in farmed allogynogenetic gibel carp strain CAS III compared with that in the wild gibel carp strain Dongting (DT). This set of polymorphisms downregulates the expression levels of crhr2 mRNA in the brain and pituitary tissues in gibel carp strain CAS III compared with those in strain DT. Furthermore, compromised HPI axis signaling is observed in gibel carp strain CAS III, such as decreased α-melanocyte stimulating hormone protein levels, plasma cortisol content, and stress responses. Moreover, enhanced activation of protein kinase B/mammalian target of rapamycin complex 1 signaling observed in the muscle tissue of strain CAS III in comparison to that in strain DT indicated elevated anabolic metabolism in strain CAS III. Thus, these studies demonstrate that the genetic markers associated with compromised HPI axis signaling, such as crhr2, are potentially useful for genetic selection toward improvement in farmed fish growth and FCE, which would reduce fishmeal consumption and thereby indirectly facilitate sustainable fisheries.

Goldfish comprise around 300 different strains with drastically altered and aesthetical morphologies making them suitable models for evolutionary developmental biology. The dragon-eye strain is characterized by protruding eyes (analogous to those of Chinese dragons). Although the strain has been selected for about 400 years, the mechanism of its eye development remains unclear. In this study, a stable dragon-eye goldfish strain with a clear genetic background was rapidly established and studied. We found that upregulation of the PPAR signaling pathway accompanied by an increase in lipid accumulation might trigger the morphological and structural transformation of the eye in dragon-eye goldfish. At the developmental stage of proptosis (eye protrusion), downregulation of the phototransduction pathway was consistent with the structural defects and myopia of the dragon-eye strain. With the impairment of retinal development, cytokine-induced inflammation was activated, especially after proptosis, similar to the pathologic symptoms of many human ocular diseases. In addition, differentially expressed transcription factors were significantly enriched in the PAX and homeobox families, two well-known transcription factor families involved in eye development. Therefore, our findings reveal the dynamic changes in key pathways during eye development in dragon-eye goldfish, and provide insights into the molecular mechanisms underlying drastically altered eyes in goldfish and human ocular disease.

Interferon regulatory factors (IRFs) constitute a family of transcription factors that synchronize interferon (IFN) antiviral response through translocating to nucleus and binding to the promoters of IFN and IFN-stimulated genes (ISGs). Fish contain 11 IRF members; however, whether or how fish IRF family genes function in IFN response remains limited. Herein, we determine the regulatory roles of 11 zebrafish IRF family members in IFN response relevant to their subcellular localization and promoter binding. Zebrafish IRF family members display three patterns of constitutive localization, only in nucleus (IRF1/2/9/11), only in cytoplasm (IRF3/5/7), and largely in nucleus with small amounts in cytoplasm (IRF4b/6/8/10). DNA pull-down assays confirm that all zebrafish IRF proteins are capable to bind fish IFN promoters, albeit to various degrees, thus regulating IFN gene transcription as activators (IRF1/3/5/6/7/8/9/11) or repressors (IRF2/4b/10). Further characterization of distinct IFN gene activation reveals that IRF1/3/5/6/7/8/9/11 efficiently stimulate zebrafish IFNφ1 expression, and IRF1/7/11 are responsible for zebrafish IFNφ3 expression. Two conserved basic residues within the helix α3 of DNA binding domains (DBDs) contribute to constitutive or inducible nuclear import for all zebrafish IRF family members and DNA binding for most members, thereby enabling them to function as transcription factors. Our results reveal a conserved and general mechanism that specifies zebrafish IRF family proteins to nuclear import and DNA binding, thereby regulating fish IFN response.

Early embryogenesis relies on maternally inherited mRNAs. Although the mechanism of maternal mRNA degradation during maternal-to-zygotic transition (MZT) has been extensively studied in vertebrates, how the embryos maintain maternal mRNA stability remains unclear. Here, we identify Igf2bp3 as an important regulator of maternal mRNA stability in zebrafish. Depletion of maternal igf2bp3 destabilizes maternal mRNAs prior to MZT and leads to severe developmental defects, including abnormal cytoskeleton organization and cell division. However, the process of oogenesis and the expression levels of maternal mRNAs in unfertilized eggs are normal in maternal igf2bp3 mutants. Gene ontology analysis revealed that these functions are largely mediated by Igf2bp3-bound mRNAs. Indeed, Igf2bp3 depletion destabilizes while its overexpression enhances its targeting maternal mRNAs. Interestingly, igf2bp3 overexpression in wild-type embryos also causes a developmental delay. Altogether, these findings highlight an important function of Igf2bp3 in controlling early zebrafish embryogenesis by binding and regulating the stability of maternal mRNAs., In early embryogenesis, maternal mRNAs are degraded at the maternal-to-zygotic transition (MZT). Mei et al find that the stability of maternal mRNAs before MZT is maintained by the mRNA binding protein Igf2bp3, which is also required for normal zebrafish development.

Fish interferon (IFN)-mediated antiviral innate immunity is the first line of defense against virus invasion. In the present study, we identify two fish IFN genes (here tentatively named IFNa and IFNc) with different-sized 3' UTRs from clone F strain of gibel carp Carassius auratus gibelio. Carp IFNa has a relatively short 3'UTR without AU-rich elements (AREs) but IFNc has a long one with 9 AREs. Functionally, carp IFNa and IFNc display significantly antiviral potential to viral infection, likely through induction of downstream IFN-stimulated genes (ISGs). Both carp IFN genes are induced by viral infection, poly(I:C) treatment and IRF3/7, which are ascribed to the IFN-sensitive response elements (ISRE) within their promoters. Carp IFN genes are also induced by each other and by themselves, indicating existence of a positive feedback loop in fish IFN-mediated antiviral immune response. Comparative analyses of 3'UTR-mediated expression regulation at mRNA and protein levels show that the ARE-containing 3'UTR of carp IFNc rather than the short 3'UTR of carp IFNa promotes mRNA decay but instead results in high-level protein expression, indicating that 3'UTR of fish IFN mRNAs might be a potential factor for regulation of IFN-mediated antiviral immune response. Considering a fact that a given protein function is largely related to its protein level, these results suggest that both promoter and 3'UTR contribute to the transcription and translation of fish IFN genes, thus shaping their eventually antiviral potential.

How unisexual animals eliminate deleterious mutations to avoid dead ends is one of the most interesting puzzles in evolutionary genetics. Incorporation of microchromosomes derived from exogenous sperm had been observed in gynogenetic animals, but little is known about their detailed process and hereditary fate. Here, we show a stable genome incorporation case in an artificial clone F of gynogenetic gibel carp (Carassius gibelio). A total of 12 exogenous DNA fragments were screened through a read depth-dependent comparison strategy and confirmed to be specific to the clone F and the paternal blunt snout bream (Megalobrama amblycephala Yin) by SCAR (sequence characterized amplified regions) marker detection. Moreover, these sperm-derived DNA fragments were not detected in some samples in early gynogenetic generations, but they were found to exist in all examined individuals through artificial gynogenetic selections of 13 generations, implying that they might have stably incorporated into the genome of clone F. Furthermore, chromosome localization and sequence characterization indicate that the largest fragment CgA22_34 is derived from blunt snout bream non-LTR retrotransposon and durably incorporated into only one of three homologous chromosomes of gibel carp clone F. Our results suggest that the incorporated sperm-derived DNA fragments by allogynogenesis should increase genetic diversity and introduce new traits into unisexual animals which will benefit genetic breeding of gibel carp. During the process, transposable elements (TEs) may play significant roles in shaping the genome structures. Simultaneously, the incorporated DNA fragments are able to be used as genetic markers to perform selective breeding programs in aquaculture practices of gibel carp.

From insects to mammals, both innate and adaptive immune response are usually higher in females than in males, with the sex chromosome and hormonal differences considered the main reasons. Here, we report that zebrafish cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a), an autosomal gene with female-biased expression, causes female fish to exhibit a lower antiviral response. First, we successfully constructed an infection model by intraperitoneal injection of spring viremia of carp virus (SVCV) into zebrafish (Danio rerio) and Carassius auratus herpesvirus (CaHV) in gibel carp (Carassius gibelio). Specifically, female fish were more vulnerable to viral infection than males, accompanied by a significantly weaker interferon (IFN) expression. After screening several candidates, cyp19a1a, which was highly expressed in female fish tissues, was selected for further analysis. The IFN expression and antiviral response were significantly higher in cyp19a1a-/- than in cyp19a1a+/+. Further investigation of the molecular mechanism revealed that Cyp19a1a targets mediator of IRF3 activation (MITA) for autophagic degradation. Interestingly, in the absence of MITA, Cyp19a1a alone could not elicit an autophagic response. Furthermore, the autophagy factor ATG14 (autophagy-related 14) was found interacted with Cyp19a1a to either promote or attenuate Cyp19a1a-mediated MITA degradation by either being overexpressed or knocked down, respectively. At the cellular level, both the normal and MITA-enhanced cellular antiviral responses were diminished by Cyp19a1a. These findings demonstrated a sex difference in the antiviral response based on a regulation mechanism controlled by a female-biased gene besides sex chromosome and hormonal differences, supplying the current understanding of sex differences in fish.

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Src homology region 2 domain-containing phosphatase 1 (SHP1), encoded by the protein tyrosine phosphatase nonreceptor type 6 (ptpn6) gene, belongs to the family of protein tyrosine phosphatases (PTPs) and participates in multiple signaling pathways of immune cells. However, the mechanism of SHP1 in regulating fish immunity is largely unknown. In this study, we first identified two gibel carp (Carassius gibelio) ptpn6 homeologs (Cgptpn6-A and Cgptpn6-B), each of which had three alleles with high identities. Then, relative to Cgptpn6-B, dominant expression in adult tissues and higher upregulated expression of Cgptpn6-A induced by polyinosinic-polycytidylic acid (poly I:C), poly deoxyadenylic-deoxythymidylic (dA:dT) acid and spring viremia of carp virus (SVCV) were uncovered. Finally, we demonstrated that CgSHP1-A (encoded by the Cgptpn6-A gene) and CgSHP1-B (encoded by the Cgptpn6-B gene) act as negative regulators of the RIG-I-like receptor (RLR)-mediated interferon (IFN) response via two mechanisms: the inhibition of CaTBK1-induced phosphorylation of CaMITA shared by CgSHP1-A and CgSHP1-B, and the autophagic degradation of CaMITA exclusively by CgSHP1-A. Meanwhile, the data support that CgSHP1-A and CgSHP1-B have sub-functionalized and that CgSHP1-A overwhelmingly dominates CgSHP1-B in the process of RLR-mediated IFN response. The current study not only sheds light on the regulative mechanism of SHP1 in fish immunity, but also provides a typical case of duplicated gene evolutionary fates.

Recent studies have related the membrane-associated RING-CH-type finger (MARCH) family proteins to host innate immune response. Zebrafish (Danio rerio) MARCH8 is reported to target SVCV glycoprotein for degradation; however, little is known about whether fish MARCH8 is involved in innate interferon (IFN) response. In this study, zebrafish march8 was significantly induced by SVCV infection. Overexpression of MARCH8 diminished fish IFN-mediated antiviral response, thus promoting the replication of SVCV and GCRV in fish cells. Mechanistically, MARCH8 interacts with and degrades MITA and TBK1 proteins to inhibit IFN response. Moreover, MARCH8 has an E3 ligase activity and enhances MITA and TBK1 polyubiquitination. Our findings reveal a mechanism whereby zebrafish MARCH8 downregulates fish IFN response and facilitates viral replication by targeting MITA and TBK1 for protein degradation.

In mammals and fish, emerging evidence highlights that TRIM family members play important roles in the interferon (IFN) antiviral immune response. Fish TRIM family has undergone an unprecedented expansion leading to generation of finTRIM subfamily, which is exclusively specific to fish. Our recent results have shown that FTRCA1 (finTRIM C. auratus 1) is likely a fish species-specific finTRIM member in crucian carp C. auratus and acts as a negative modulator to downregulate fish IFN response by autophage-lysosomal degradation of protein kinase TBK1. In the present study, we found that FTRCA1 also impedes the activation of crucian carp IFN promoter by IRF7 but not by IRF3. Mechanistically, FTRCA1 attenuates IRF7 transcription levels likely due to enhanced decay of IRF7 mRNA, leading to reduced IRF7 protein levels and subsequently reduced fish IFN expression. E3 ligase activity is required for FTRCA1 to negatively regulate IRF7-mediated IFN response, because ligase-inactive mutants and the RING-deleted mutant of FTRCA1 lose the ability to block the activation of crucian carp IFN promoter by IRF7. These results together indicate that FTRCA1 is a multifaceted modulator to target different signaling factors for shaping fish IFN response in crucian carp.

Background Accompanied with rapid growth and high density aquaculture, gibel carp has been seriously threatened by Carassius auratus herpesvirus (CaHV) since 2012. In previous study, distinct CaHV resistances and immune responses were revealed in the diseased individuals of three gibel carp gynogenetic clones (A+, F and H). However, little is known about the gene expression changes in the survivors after CaHV challenge, particularly their differences of innate and adaptive immune system between susceptible clone and resistant clone. Results We firstly confirmed the CaHV carrier state in the survivors of three gibel carp clones after CaHV challenge by evaluating the abundances of five CaHV genes. The assay of viral loads indicated the resistant clone H possessed not only stronger resistance but also higher tolerance to CaHV. Then, 2818, 4047 and 3323 differentially expressed unigenes (DEUs) were screened from the head-kidney transcriptome profiles of survivors compared with controls from clone A+, F and H. GO and KEGG analysis suggested that a persistent immune response might sustain in resistant clone H and F, while susceptible clone A+ had a long-term impact on the circulatory system which was consistent with the major symptoms of bleeding caused by CaHV. Among the top 30 enriched pathways of specifically up-regulated DEUs in respective clones, 26, 7 and 15 pathways in clone H, F and A+ were associated with infections, diseases, or immune-related pathways respectively. In addition, 20 pathways in clone F belonged to “metabolism” or “biogenesis”, and 7 pathways involved in “circulatory system” were enriched in clone A+. Significantly, we revealed the differential expression changes of IFN system genes and immunoglobulin (Ig) genes among the survivors of three clones. Finally, myosins and Igs were identified as co-expression modules which were positively or negatively correlated to CaHV viral loads respectively. Conclusions Our results revealed the common and distinct gene expression changes in immune and circulatory system in the survivors of three gibel carp gynogenetic clones with different CaHV resistances. The current study represents a paradigm of differential innate and adaptive immune reactions in teleost, and will be beneficial to the disease-resistance breeding of gibel carp. Electronic supplementary material The online version of this article (10.1186/s12864-019-5777-z) contains supplementary material, which is available to authorized users.

Sinorhodeus microlepis (S. microlepis) is recently described as a new species and represents a new genus Sinorhodeu of the subfamily Acheilognathinae. In this study, we first sequenced the complete mitogenome of S. microlepis and compared with the other 29 bitterling mitogenomes. The S. microlepis mitogenome is 16,591 bp in length and contains 37 genes. Gene distribution pattern is identical among 30 bitterling mitogenomes. A significant linear correlation between A+T% and AT-skew were found among 29 bitterling mitogenomes, except S. microlepis shows unusual AT-skew with slightly negative in tRNAs and PCGs. Bitterling mitogenomes exhibit highly conserved usage bias of start codon, relative synonymous codons and amino acids, overlaps and non-coding intergenic spacers. Phylogenetic trees constructed by 13 PCGs strongly support the polyphyly of the genus Acheilognathus and the paraphyly of Rhodeus and Tanakia. Together with the unusual characters of S. microlepis mitogenomes and phylogenetic trees, S. microlepis should be a sister species to the genus Rhodeu that might diverge about 13.69 Ma (95% HPD: 12.96–14.48 Ma).

In mammals, tripartite motif (TRIM) proteins have emerged as pivotal players endowed with, directly, antiviral effects and, indirectly, modulatory capacity of the innate immune response. An unprecedented expansion of TRIM family has occurred in fish; however, the functional role of fish TRIM family members remains largely unknown. In this study, we identify a species-specific TRIM gene from crucian carp Carassius auratus, named FTRCA1, phylogenetically similar to the members of finTRIM, a subfamily of TRIM exclusively in teleost fish. FTRCA1 is induced by IFN and IFN stimuli as a typical IFN-stimulated gene. Overexpression of FTRCA1 negatively regulates IFN antiviral response by inhibition of IRF3 phosphorylation; consistently, knockdown of FTRCA1 results in enhanced levels of IRF3 phosphorylation and also IFN expression following poly(I:C) transfection. Whereas FTRCA1 is associated with several pivotal signaling molecules of RIG-I–like receptor pathway, its association with TBK1 results in autophage-lysosomal degradation of TBK1, thus abrogating the downstream IFN induction. Interestingly, FTRCA1 is phosphorylated by TBK1, but this phosphorylation is not required for downregulation of TBK1 protein. Transfection assays indicate that FTRCA1 is likely an E3 ligase with the requirement of RING finger domain, and deletion of N-terminal RING domain or mutation of seven conservative sites abolishes the negative regulatory function of FTRCA1. Collectively, these results illuminate a novel finTRIM-mediated innate immune modulatory pathway, thus providing insights into species-specific regulation of fish IFN response.

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

Background Fatty liver has become a main problem that causes huge economic losses in many aquaculture modes. It is a common physiological or pathological phenomenon in aquaculture, but the causes and occurring mechanism are remaining enigmatic. Methods Each three liver samples from the control group of allogynogenetic gibel carp with normal liver and the overfeeding group with fatty liver were collected randomly for the detailed comparison of histological structure, lipid accumulation, transcriptomic profile, latent pathway identification analysis (LPIA), marker gene expression, and hepatocyte mitochondria analyses. Results Compared to normal liver, larger hepatocytes and more lipid accumulation were observed in fatty liver. Transcriptomic analysis between fatty liver and normal liver showed a totally different transcriptional trajectory. GO terms and KEGG pathways analyses revealed several enriched pathways in fatty liver, such as lipid biosynthesis, degradation accumulation, peroxidation, or metabolism and redox balance activities. LPIA identified an activated ferroptosis pathway in the fatty liver. qPCR analysis confirmed that gpx4, a negative regulator of ferroptosis, was significantly downregulated while the other three positively regulated marker genes, such as acsl4, tfr1 and gcl, were upregulated in fatty liver. Moreover, the hepatocytes of fatty liver had more condensed mitochondria and some of their outer membranes were almost ruptured. Conclusions We reveal an association between ferroptosis and fish fatty liver for the first time, suggesting that ferroptosis might be activated in liver fatty. Therefore, the current study provides a clue for future studies on fish fatty liver problems.

Polyploidy and subsequent diploidization provide genomic opportunities for evolutionary innovations and adaptation. The researches on duplicated gene evolutionary fates in recurrent polyploids have seriously lagged behind that in paleopolyploids with diploidized genomes. Moreover, the antiviral mechanisms of Viperin remain largely unclear in fish. Here, we elaborate the distinct antiviral mechanisms of two viperin homeologs (Cgviperin-A and Cgviperin-B) in auto-allo-hexaploid gibel carp (Carassius gibelio). First, Cgviperin-A and Cgviperin-B showed differential and biased expression patterns in gibel carp adult tissues. Subsequently, using co-immunoprecipitation (Co-IP) screening analysis, both CgViperin-A and CgViperin-B were found to interact with crucian carp (C. auratus) herpesvirus (CaHV) open reading frame 46 right (ORF46R) protein, a negative herpesvirus regulator of host interferon (IFN) production, and to promote the proteasomal degradation of ORF46R via decreasing K63-linked ubiquitination. Additionally, CgViperin-B also mediated ORF46R degradation through autophagosome pathway, which was absent in CgViperin-A. Moreover, we found that the N-terminal α-helix domain was necessary for the localization of CgViperin-A and CgViperin-B at the endoplasmic reticulum (ER), and the C-terminal domain of CgViperin-A and CgViperin-B was indispensable for the interaction with degradation of ORF46R. Therefore, the current findings clarify the divergent antiviral mechanisms of the duplicated viperin homeologs in a recurrent polyploid fish, which will shed light on the evolution of teleost duplicated genes.

Anti-Mullerian hormone receptor type II (Amhr2) is a key receptor of Amh signaling in regulating gonad development. The amhr2 gene has been identified in numerous species, including a few teleost fishes. However, the roles of Amhr2 in Amh signaling in fish are poorly studied.In this study, an amhr2 homolog from obscure puffer (Takifugu obscurus) was identified, and its molecular characteristics were systematically analyzed. Expression analysis revealed that amhr2 was highly expressed in the gonads of adult pufferfish and significantly upregulated during sex differentiation. Significantly, a sex-linked SNP site was verified in obscure puffer amhr2. Females exhibited a homozygous genotype (C/C), while males possessed a heterozygous genotype (C/G), resulting in an amino acid variation (His/Asp384) in the kinase domain of Amhr2. Then, the functions of the different Amhr2 genotypes were further investigated. The male genotype protein (Amhr2This study provides direct experimental evidence for the roles of different Amhr2 genotypes in pufferfish and suggests that amhr2 is responsible for male sex differentiation in obscure puffer.

Carassius auratus herpesvirus (CaHV) has been identified as a high-virulence pathogenic virus that infects aquatic animals, but the key factor for virus–host interaction is still unclear. Five Really interesting new genes (RING) finger proteins (39L, 52L, 131R, 136L, and 143R) of CaHV were screened to determine structural diversity. RING finger proteins were also predicted in other known fish herpesviruses, with an arrangement and number similar to CaHV. We performed multifaceted analyses of the proteins, including protein sizes, skeleton structures, subcellular localizations, and ubiquitination activities, to determine their precise roles in virus–host interactions. The five proteins were overexpressed and detected different levels of ubiquitination activities, and 143R showed the highest activity. Then, the prokaryotic expressed and purified full-length proteins (131R and 136L), RING domain isolates (131R12–43 and 136L45–87), and RING domain-deleted mutants (131RΔ12–43 and 136LΔ45–87) were prepared to detect their activities through ubiquitination assays. The results indicate that both full-length proteins and their isolates have activities that catalyze ubiquitination, and the full-length proteins possess higher activity than the isolates, but RING domain-deleted mutants lose their activities. Furthermore, the activities of the five proteins were verified as E3 ubiquitin ligase activity, showing that the RING domains determine the ubiquitination activity. These proteins present different subcellular localization. RING domain-deleted mutants showed similar subcellular localization with their full-length proteins, and all the isolates diffused in whole cells. The current results indicate that the sequence outside the RING domain determines subcellular localization and the level of ubiquitination activity, suggesting that the RING finger proteins of fish herpesviruses might have diverse functions in virus–host interaction.

Background Loaches of Cobitinae, widely distributed in Eurasian continent, have high economic, ornamental and scientific value. However, the phylogeny of Cobitinae fishes within genera or family level remains complex and controversial. Up to now, about 60 Cobitinae mitogenomes had been deposited in GenBank, but their integrated characteristics were not elaborated. Results In this study, we sequenced and analyzed the complete mitogenomes of a female Cobits macrostigma. Then we conducted a comparative mitogenome analysis and revealed the conserved and unique characteristics of 58 Cobitinae mitogenomes, including C. macrostigma. Cobitinae mitogenomes display highly conserved tRNA secondary structure, overlaps and non-coding intergenic spacers. In addition, distinct base compositions were observed among different genus and significantly negative linear correlation between AT% and AT-skew were found among Cobitinae, genus Cobitis and Pangio mitogenomes, respectively. A specific 3 bp insertion (GCA) in the atp8-atp6 overlap was identified as a unique feature of loaches, compared to other Cypriniformes fish. Additionally, all protein coding genes underwent a strong purifying selection. Phylogenetic analysis strongly supported the paraphyly of Cobitis and polyphyly of Misgurnus. The strict molecular clock predicted that Cobitinae might have split into northern and southern lineages in the late Eocene (42.11 Ma), furthermore, mtDNA introgression might occur (14.40 Ma) between ancestral species of Cobitis and ancestral species of Misgurnus. Conclusions The current study represents the first comparative mitogenomic and phylogenetic analyses within Cobitinae and provides new insights into the mitogenome features and evolution of fishes belonging to the cobitinae family.

Evolutionary fates of duplicated genes have been widely investigated in many polyploid plants and animals, but research is scarce in recurrent polyploids. In this study, we focused on foxl2, a central player in ovary, and elaborated the functional divergence in gibel carp (Carassius gibelio), a recurrent auto-allo-hexaploid fish. First, we identified three divergent foxl2 homeologs (Cgfoxl2a-B, Cgfoxl2b-A, and Cgfoxl2b-B), each of them possessing three highly conserved alleles and revealed their biased retention/loss. Then, their abundant sexual dimorphism and biased expression were uncovered in hypothalamic–pituitary–gonadal axis. Significantly, granulosa cells and three subpopulations of thecal cells were distinguished by cellular localization of CgFoxl2a and CgFoxl2b, and the functional roles and the involved process were traced in folliculogenesis. Finally, we successfully edited multiple foxl2 homeologs and/or alleles by using CRISPR/Cas9. Cgfoxl2a-B deficiency led to ovary development arrest or complete sex reversal, whereas complete disruption of Cgfoxl2b-A and Cgfoxl2b-B resulted in the depletion of germ cells. Taken together, the detailed cellular localization and functional differences indicate that Cgfoxl2a and Cgfoxl2b have subfunctionalized and cooperated to regulate folliculogenesis and gonad differentiation, and Cgfoxl2b has evolved a new function in oogenesis. Therefore, the current study provides a typical case of homeolog/allele diversification, retention/loss, biased expression, and sub-/neofunctionalization in the evolution of duplicated genes driven by polyploidy and subsequent diploidization from the recurrent polyploid fish.

The eponymous member of the interferon regulatory factor (IRF) family, IRF1, was originally identified as a nuclear factor that binds and activates the promoters of type I interferon genes. However, subsequent studies using genetic knockouts or RNAi-mediated depletion of IRF1 provide a much broader view, linking IRF1 to a wide range of functions in protection against invading pathogens. Conserved throughout vertebrate evolution, IRF1 has been shown in recent years to mediate constitutive as well as inducible host defenses against a variety of viruses. Fine-tuning of these ancient IRF1-mediated host defenses, and countering strategies by pathogens to disarm IRF1, play crucial roles in pathogenesis and determining the outcome of infection.

Environmental bacteria have a great impact on fish gut microbiota, yet little is known as to where fish acquire their gut symbionts, and how gut microbiota response to the disturbance from environmental bacteria. Through the integrative analysis by community profiling and source tracking, we show that feed-associated bacteria can impose a strong disturbance upon the hindgut microbiota of cultured fugu. Consequently, marked alterations in the composition and function of gut microbiota in slow growth fugu were observed, implying a reduced stability upon bacterial disturbance from feed. Moreover, quantitative ecological analyses indicated that homogeneous selection and dispersal limitation largely contribute to the community stability and partial variations among hosts in the context of lower degree of disturbance. While the disturbance peaked, variable selection leads to an augmented interaction within gut microbiota, entailing community unstability and shift. Our findings emphasized the intricate linkage between feed and gut microbiota and highlighted the importance of resolving the feed source signal before the conclusion of comparative analysis of microbiota can be drawn. Our results provide a deeper insight into aquaculture of fugu and other economically important fishes and have further implications for an improved understanding of host-microbe interactions in the vertebrate gastrointestinal tract.

Hepatocellular carcinoma (HCC) is one of the deadliest cancer types worldwide. The centromere proteins (CENPs) are critical for the mitosis-related protein complex and are involved in kinetochore assembly and spindle checkpoint signaling during mitosis. However, the clinical significance of CENPs in the recurrence and progression of HCC remains poorly understood. Here, we examined the expression of all CENPs and their association with recurrence and survival of HCC patients using the global gene expression profile dataset established in our laboratory. The effect of silencing CENPF on cell viability, migration, and epithelial-mesenchymal transition (EMT) were detected using CCK-8, transwell, and western blot, respectively. RT-qPCR and western blot were performed to confirm the silencing of CENPF and the relationship between STAT5A and CENPF, while tumorigenesis was tested using the HCC Huh7 xenograft mouse model. Most of the CENPs is overexpressed in HCC, and overexpression of CENPF was significantly associated with the poor survival of HCC patients. CENPF promoted HCC cell lines migration and EMT progression. Knockdown CENPF inhibited cell growth activity against human HCC cells in vitro and xenograft tumors in vivo. Bioinformatics analysis revealed that CENPF genes are enriched in the cell cycle. Silencing CENPF arrested cell cycle at the G2/M phase and inhibited Cyclin B1 and Cyclin E1 expressions. Meanwhile, silencing CENPF prohibited phosphorylation of ERK and the expression of NEK2. Additionally, we found that STAT5A down-regulated CENPF expression and inhibited cancer cell growth viability. In conclusion, our data suggested that CENPF could be potentially developed into a theranostic biomarker to tackle HCC progression.

Chemokine receptor cxcr4 and its ligand cxcl12 have evolved two paralogs in the teleost lineage. In this study, we have identified four duplicated cxcr4 and cxcl12 genes from hexaploid gibel carp, Carassius gibelio, respectively. Cgcxcr4bs and Cgcxcl12as were dynamically and differentially expressed in immune-related tissues, and significantly up-regulated in head kidney and spleen after crucian carp herpesvirus (CaHV) infection. Blocking Cxcr4/Cxcl12 axis by injecting AMD3100 brought more severe bleeding symptom and lower survival rate in CaHV-infected fish. AMD3100 treatment also suppressed the up-regulation of key antiviral genes in head kidney and spleen, and resulted in more acute replication of CaHV in vivo. Consistently, the similar suppression of up-regulated expression of key antiviral genes were also observed in CAB cells treated by AMD3100 after poly(I:C) stimulation. Finally, MAPK3 and JAK/STAT were identified as the possible pathways that CgCxcr4s and CgCxcl12s participate in to promote the antiviral response in vitro.

In mammals, transcription factors of IFN-regulatory factors (IRFs) family translate viral recognition into IFN antiviral responses through translocating to nucleus and subsequently binding to the promoters of IFN and IFN-stimulated genes (ISGs). In addition to IRF1-9 conserved across vertebrates and IRF10 in teleost fish and bird, teleost fish has another novel member, IRF11; however, little is known about its role in IFN response. In this study, we provide evidence that IRF11 is present only in Osteichthyes (bony fish) but lost in tetrapods and subsequently characterize the stimulatory potential of zebrafish IRF11 to IFN antiviral response relevant to its subcellular localization and promoter binding. Overexpression of zebrafish IRF11 restricts virus replication through induction of IFN and ISGs. Zebrafish IRF11 is constitutively localized to nucleus, which is driven by a tripartite NLS motif, consisting of three interdependent basic clusters, two in DNA binding domain (DBD) and one in the region immediately C-terminal to DBD. Nuclear IRF11 binds to the IRF-binding element/IFN-stimulated response element motifs of zebrafish IFN promoters depending on the two conserved amino acids (K78, R82) within DBD helix α3. K78 and R82 also benefit zebrafish IRF11 nuclear import as two key residues positioned at the first basic cluster of the tripartite NLS motif. Such features enable zebrafish IRF11 to function as a positive transcription factor for fish IFN antiviral response. Our results identify a unique tripartite NLS motif that integrates DNA-binding activity and nuclear import ability, allowing zebrafish IRF11 to initiate IFN and ISG expression.

The reproductive process is usually controlled by the hypothalamic-pituitary-gonad axis in vertebrates, while Kiss/gonadotropin-releasing hormone (GnRH) system in the hypothalamus is required for mammalian reproduction but dispensable for fish reproduction. The regulation of follicle stimulating hormone/luteinizing hormone (LH) expression in fish species is still unknown. Here, we identified miR-200s on chromosome 23 (chr23-miR-200s) as important regulators for female zebrafish reproduction. Knockout of chr23-miR-200s (chr23-miR-200s-KO) resulted in dysregulated expression of luteinizing hormone beta lhb (luteinizing hormone beta) and some hormone genes in the pituitary as revealed by comparative transcriptome profiling, leading to failure of oocyte maturation and ovulation as well as defects in reproductive duct development. Chr23-miR-200s mainly expressed in the pituitary and regulated lhb expression by targeting the transcription repressor wt1a. Injection of human chorionic gonadotropin (hCG) could rescue the defects of oocyte maturation in chr23-miR-200s-KO zebrafish, whereas GnRH or LHRH-A2 could not, suggesting that Chr23-miR-200s regulated lhb expression in a GnRH-independent pathway. It was remarkable that either injection of carp pituitary extraction, or co-injection of hCG with synthetic oxytocin and vasotocin could greatly rescue the defects of both oocyte maturation and ovulation in chr23-miR-200s-KO zebrafish. Altogether, our findings highlight an important function of chr23-miR-200s in controlling oocyte maturation by regulation LH expression, and oxytocin and vasotocin are potentially responsible for the ovulation in fish species.

Polyploids in vertebrates are generally associated with unisexual reproduction, but the direct consequences of polyploidy on sex determination system and reproduction mode remain unknown. Here, we synthesized a group of artificial octoploids between unisexual gynogenetic hexaploid Carassius gibelio and sexual tetraploid Carassius auratus. The synthetic octoploids were revealed to have more than 200 chromosomes, in which 50 chromosomes including the X/Y sex determination system were identified to transfer from sexual tetraploid C. auratus into the unisexual gynogenetic hexaploid C. gibelio. Significantly, a few synthetic octoploid males were found to be fertile, and one octoploid male was confirmed to regain sexual reproduction ability, which exhibits characteristics that are the same to sexual reproduction tetraploid males, such as 1:1 sex ratio occurrence, meiosis completion and euploid sperm formation in spermatogenesis, as well as normal embryo development and gene expression pattern during embryogenesis. Therefore, the current finding provides a unique case to explore the effect of sex determination system incorporation on reproduction mode transition from unisexual gynogenesis to sexual reproduction along with genome synthesis of recurrent polyploidy in vertebrates.

Sex is prevalent in nature and sex determination is one of the most fundamental biological processes, while the way of initiating female and male development exhibits remarkable diversity and variability across vertebrates. The knowledge on why and how sex determination mechanisms evolve unusual plasticity remains limited. Here, we summarize sex determination systems, master sex-determining genes and gene-regulatory networks among vertebrates. Recent research advancements on sex determination system transition are also introduced and discussed in some non-model animals with multiple sex determination mechanisms. This review will provide insights into the origin, transition and evolutionary adaption of different sex determination strategies in vertebrates, as well as clues for future perspectives in this field.

Aquatic viruses include infected viruses in aquatic animals, plants and microorganisms, and free-floating viruses (virioplankton) in water environments. In the last three decades, a huge number of aquatic viruses, especially diverse free-floating viruses, including cyanophages, phycoviruses, archaea viruses, giant viruses, and even virophages, have been identified by virological experiments and metagenomic analyses. Based on a comprehensive introduction of aquatic virus classification and their morphological and genetic diversity, here, we summarize and outline main virus species, their evolutionary contribution to aquatic communities through horizontal gene transfer, and their ecological roles for cyanobacterial bloom termination and global biogeochemical cycling in freshwater and marine ecosystems. Thereby, some novel insights of aquatic viruses and virus-host interactions, especially their evolutionary contribution and ecological rolesin diverse aquatic communities and ecosystems, are highlighted in this review.

Grass carp is one of the most important freshwater aquaculture species in China. However, the mechanisms underlying the growth of muscle tissue in the fish are unclear. High-throughput RNA-Seq was used to analyze the transcriptome of grass carp muscle tissue between fast- and slow-growing fish family groups. Twenty-four individuals each from 4 fast-growing families and 4 slow-growing families were used to reduce background noise. 71 up-regulated and 35 down-regulated genes were identified in the differentially expressed genes (DEGs). GO and KEGG enrichment analyses revealed the DEGs were involved in the GH/IGF axis, calcium metabolism, protein and glycogen synthesis, oxygen transport, cytoskeletal and myofibrillar components. IGFBP1 was up-regulated in big fish while GHR2 was down-regulated. Glutamic pyruvate transaminase 2, an indicator of liver tissue damage, was down-regulated in big grass carp, which indicates that the fish was better adapted to an artificially formulated diet. GAPDH, the rate-limiting enzyme in glycolytic flux was highly expressed in fast-growing grass carp, reflecting enhanced carbohydrate metabolism. Higher expression of ALAS2 and myoglobin 1 in big grass carp, related to oxygen transport might promote aerobic exercise along with food intake and muscle growth. Genes for cytoskeletal and myofibrillar components such as tropomyosin, meromyosin, and troponin I were also up-regulated in big grass carp. These results provide valuable information about the key genes for use as biomarkers of growth in selective breeding programs for grass carp and contribute to our understanding of the molecular mechanisms and regulative pathways regulating growth in fish.

Diverse immunoglobulin (Ig) domain-containing protein (DICP) family is a novel bony fish-specific multi-gene family encoding diversified immune receptors. However, their function and the implication of binding partners remain unknown. In this study, we first identified 28 DICPs from three gibel carp gynogenetic clones and revealed their high variability and clone-specific feature. After crucian carp herpesvirus (CaHV) infection, these DICPs were significantly upregulated in head kidney, kidney and spleen. The up-regulation folds in clone A+, F and H were related to the susceptibility to CaHV, progressively increasing from resistant clone to susceptible clone. Overexpression of gibel carp DICPs inhibited interferon (IFN) and viperin promoter-driven luciferase activity. The additions of E. coli extracts and lipid A significantly enhanced the inhibition effect. In addition, gibel carp DICPs can interact with SHP-1 and SHP-2. These findings suggest that gible carp DICPs, as inhibitory receptors, might specifically recognize lipid A, and then interact with SHP-1 and SHP-2 to inhibit the induction of IFN and ISGs.

Long non-coding RNAs (lncRNAs) perform distinct functions in various biological processes in mammals, including sex differentiation. However, the roles of lncRNAs in other vertebrates, especially in the Chinese soft-shell turtle (Pelodiscus sinensis), remain to be clarified. In this study, we performed genome-wide analysis of the lncRNA expression profiles in gonad tissues and screened numerous sex-specific lncRNAs in the Chinese soft-shell turtle. Of the 363,310,650 clean reads obtained, 5,994 sequences were typed as lncRNAs, of which 4,463 were novel. A selection of sex-specific lncRNAs (♀ 932, ♂ 449) from female ovaries and male testis were shown to act on target genes in cis and in trans, and most were involved in gonad differentiation based on Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Furthermore, interactions among the differentially expressed lncRNA-mRNAs and protein coding genes were identified by construction of correlation networks. Overall, our systematic analysis of lncRNA expression profiles in gonad tissues revealed numerous sex-specific lncRNAs in P. sinensis. Thereby, these findings provide new insights into the function of lncRNAs in sex differentiation and highlight a group of candidate lncRNAs for future studies.

Background Ranaviruses (family Iridoviridae, nucleocytoplasmic large DNA viruses) have been reported as promiscuous pathogens of cold-blooded vertebrates. Rana grylio virus (RGV, a ranavirus), from diseased frog Rana grylio with a genome of 105.79 kb and Andrias davidianus ranavirus (ADRV), from diseased Chinese giant salamander (CGS) with a genome of 106.73 kb, contains 99% homologous genes. Results To uncover the differences in virus replication and host responses under interspecies infection, we analyzed transcriptomes of CGS challenged with RGV and ADRV in different time points (1d, 7d) for the first time. A total of 128,533 unigenes were obtained from 820,858,128 clean reads. Transcriptome analysis revealed stronger gene expression of RGV than ADRV at 1 d post infection (dpi), which was supported by infection in vitro. RGV replicated faster and had higher titers than ADRV in cultured CGS cell line. RT-qPCR revealed the RGV genes including the immediate early gene (RGV-89R) had higher expression level than that of ADRV at 1 dpi. It further verified the acute infection of RGV in interspecies infection. The number of differentially expressed genes and enriched pathways from RGV were lower than that from ADRV, which reflected the variant host responses at transcriptional level. No obvious changes of key components in pathway “Antigen processing and presentation” were detected for RGV at 1 dpi. Contrarily, ADRV infection down-regulated the expression levels of MHC I and CD8. The divergent host immune responses revealed the differences between interspecies and natural infection, which may resulted in different fates of the two viruses. Altogether, these results revealed the differences in transcriptome responses among ranavirus interspecies infection of amphibian and new insights in DNA virus-host interactions in interspecies infection. Conclusion The DNA virus (RGV) not only expressed self-genes and replicated quickly after entry into host under interspecies infection, but also avoided the over-activation of host responses. The strategy could gain time for the survival of interspecies pathogen, and may provide opportunity for its adaptive evolution and interspecies transmission. Electronic supplementary material The online version of this article (10.1186/s12864-018-4596-y) contains supplementary material, which is available to authorized users.

In mammals, RIG-I like receptors (RLRs) RIG-I and melanoma differentiation–associated gene 5 (MDA5) sense cytosolic viral RNA, leading to IFN antiviral response; however, LGP2 exhibits controversial functions. The same happens to fish LGP2. In this study we report that three zebrafish LGP2 splicing transcripts, a full-length LGP2, and two truncating variants, LGP2v1 and LGP2v2, play distinct roles during IFN antiviral response. Overexpression of the full-length LGP2 not only potentiates IFN response through the RLR pathway, in the absence or presence of poly(I:C) at limited concentrations, but also inhibits IFN response by relative high concentrations of poly(I:C) through functional attenuation of signaling factors involved in the RLR pathway; however, LGP2v1 and LGP2v2 only retain the inhibitory role. Consistently, LGP2 but not LGP2v1 and LGP2v2 confers protection on fish cells against spring viremia of carp virus (SVCV) infection and at limited expression levels, LGP2 exerts more significant protection than either RIG-I or MDA5. Further data suggest that in the early phase of SVCV infection, LGP2 functions as a positive regulator but along with SVCV replicating in cells up to a certain titer, which leads to a far more robust expression of IFN, LGP2 switches to a negative role. These in vitro results suggest an ingenious mechanism where the three zebrafish LGP2 splicing transcripts work cooperatively to shape IFN antiviral responses.

Hybridization is an efficient method to breed new strains of aquatic animals. In the present study, we produced a hybrid puffer by crossing female obscure puffer with male tiger puffer. The hybrid puffer could live in fresh water like obscure puffer and exhibited growth superiority. The averaged body weight of 4- and 6-month-old hybrid puffer were respectively 38.06% and 38.93% higher than that of obscure puffer. Then, we analyzed the underlying genetic basis for the growth advantage of hybrid puffer by comparative transcriptome analysis. A total number of 4264 and 1285 differentially expressed genes (DEGs) were respectively identified from pituitary and liver transcriptome profiles between hybrid puffer and obscure puffer. Comprehensive analysis showed that the DEGs related with cell proliferation and differentiation, and protein synthesis and export, specifically showed higher expression levels in hybrid puffer, such as "ECM-receptor interaction", "focal adhesion", "protein export" and "protein processing in endoplasmic reticulum". While the DEGs involved in gametogenesis and carbohydrate and energy metabolism highly expressed in obscure puffer, such as "oxidative phosphorylation", "citrate cycle", "progesterone-mediated oocyte maturation" and "oocyte meiosis". Furthermore, a series of candidate genes related to the growth superiority of hybrid puffer were identified, such as fn1a, ptprc, plcg2, igf1, tgfβ1, bmp4, abl1, col1a2, col1a1a, and myl9a. These results will be beneficial to understand the molecular basis of growth superiority and helpful to the hybrid breeding of pufferfish.

Determining the sex and controlling the sex ratio are essential aspects of fish genetics that can assist in developing successful fish breeding programs. High quality genome assembly and annotations are prerequisites to determine sex-specific genes and their expression. In addition, analysis of resequencing data can identify genomic difference between male and female fishes. In this study, we performed chromosome-level genome assembly in female Ancherythroculter nigrocauda fish having low heterozygosity using PacBio reads. High-throughput chromatin conformation capture (HiC) yielded a genome of size 1054.05 Mb, with a contig N50 length of 3.40 Mb and a scaffold N50 length of 42.68 Mb. In addition, we sequenced 5 female and 5 male A. nigrocauda samples and identified sex-specific regions on LG20 Furthermore, diet-specific amino acid mutation were found on 582 genes between herbivorous and carnivorous fishes, with 26 of them exhibiting significantly different expression patterns in the liver tissue of these two types of fishes. The chromosome-level genome assembly of A. nigrocauda provides valuable resources for conducting in-depth comparative genomic studies with immense applications in fish genetic breeding and farming. Similarly, the diet-specific amino acid mutations are useful in the breeding of new strains of carnivorous fishes with an herbivorous diet.

Signal transducer and activator of transcription 5b (STAT5b) has been identified as a key downstream mediator of growth hormone (GH) signaling in somatic growth of mammalian. However, the corresponding homologue gene of Stat5b is unknown in fish species. In this study, we generated loss-of-function mutants in stat5.1 and stat5.2, two stat5 homologues existing in zebrafish. In stat5.1-deficient zebrafish, a significant reduction of body length and body weight was detected in the embryos/larvae and adults compared with the wild-type control fish, and sexual size dimorphism in adult zebrafish was also eliminated. However, the stat5.2-deficient zebrafish displayed a normal developmental phenotype during all lifespan. Chromatin immunoprecipitation combined with deep sequencing (ChIP-seq) method was adopted to further investigate the potential transcriptional targets of Stat5 protein and cast much light upon the biological function of Stat5. We identified more than 800 genes as transcriptional targets of Stat5 during zebrafish embryogenesis. KEGG analysis indicated that the Stat5 target gene network is predominantly linked to the metabolic pathways, neuroactive ligand-receptor interaction and JAK-STAT signaling pathways. Further validation studies suggested that Stat5.1 protein could directly regulate the expression of gh1, and stat5.1-mutated zebrafish showed a reduction of gh1 mRNA level. In the present study, stat5.1 was revealed as the corresponding homologue gene of Stat5b in fish species. Additionally, we found a novel molecular interaction between Stat5.1/Stat5b and GH, and unraveled a positive feedback loop Stat5.1-GH-Stat5.1 which is necessary for somatic growth and body development in zebrafish.

The turtle carapace is composed of severely deformed fused dorsal vertebrae, ribs, and bone plates. In particular, the lateral growth in the superficial layer of turtle ribs in the dorsal trunk causes an encapsulation of the scapula and pelvis. The recent study suggested that the carapacial ridge (CR) is a new model of epithelial–mesenchymal transition which is essential for the arrangement of the ribs. Therefore, it is necessary to explore the regulatory mechanism of carapacial ridge development to analyze the formation of the turtle shell. However, the current understanding of the regulatory network underlying turtle carapacial ridge development is poor due to the lack of both systematic gene screening at different carapacial ridge development stages and gene function verification studies. In this study, we obtained genome-wide gene transcription and gene translation profiles using RNA sequencing and ribosome nascent-chain complex mRNA sequencing from carapacial ridge tissues of Chinese soft-shell turtle at different development stages. A correlation analysis of the transcriptome and translatome revealed that there were 129, 670, and 135 codifferentially expressed genes, including homodirection and opposite-direction differentially expressed genes, among three comparison groups, respectively. The pathway enrichment analysis of codifferentially expressed genes from the Kyoto Encyclopedia of Genes and Genomes showed dynamic changes in signaling pathways involved in carapacial ridge development. Especially, the results revealed that the Wnt signaling pathway and MAPK signaling pathway may play important roles in turtle carapacial ridge development. In addition, Wnt and Fgf were expressed during the carapacial ridge development. Furthermore, we discovered that Wnt5a regulated carapacial ridge development through the Wnt5a/JNK pathway. Therefore, our studies uncover that the morphogenesis of the turtle carapace might function through the co-operation between conserved WNT and FGF signaling pathways. Consequently, our findings revealed the dynamic signaling pathways acting on the carapacial ridge development of Chinese soft-shell turtle and provided new insights into uncover the molecular mechanism underlying turtle shell morphogenesis.

Interferon regulatory factor 7 (IRF7) is a key mediator in regulating the type Ι IFN response. Although irf7 has been identified in more than twenty fish species, alternative splicing has not been found in teleost irf7. Alternative splicing is an important mechanism expanding the transcriptomic and proteomic diversity, and has been found in several IRF family members. Here, three alternative splicing variants of irf7 were identified and characterized in obscure puffer. The first splicing transcript (Toirf7v1) was predicted to encode 428 amino acids with a DNA-binding domain (DBD), an interaction-associated domain (IAD) and a serine-rich domain (SRD). Toirf7v2 encoded 430 amino acids caused by the intron retention, and contained the whole conserved domains. Toirf7v3 encoded a truncated protein with 337 amino acids resulting from the alternative 5' splice-site selection, and lacked part of IAD domain and the entire SRD domain. Functional studies demonstrated that all of the three isoforms could activate the expression of type I IFN and IFN-stimulated genes (ISGs). Nevertheless, the two variants (Toirf7v2 and Toirf7v3) exhibited much less ability to induce transcription of IFN and ISGs compared to the Toirf7v1. Our findings suggest that these splicing variants may have distinct roles in the regulation of immune response. These results will be beneficial to understand the functional characteristics of irf7 variants in fish.