New insights into the all-testis differentiation in zebrafish with compromised endogenous androgen and estrogen synthesis.

The regulatory mechanism of gonadal sex differentiation, which is complex and regulated by multiple factors, remains poorly understood in teleosts. Recently, we have shown that compromised androgen and estrogen synthesis with increased progestin leads to all-male differentiation with proper testis development and spermatogenesis in cytochrome P450 17a1 (cyp17a1)-/- zebrafish. In the present study, the phenotypes of female-biased sex ratio were positively correlated with higher Fanconi anemia complementation group L (fancl) expression in the gonads of doublesex and mab-3 related transcription factor 1 (dmrt1)-/- and cyp17a1-/-;dmrt1-/- fish. The additional depletion of fancl in cyp17a1-/-;dmrt1-/- zebrafish reversed the gonadal sex differentiation from all-ovary to all-testis (in cyp17a1-/-;dmrt1-/-;fancl-/- fish). Luciferase assay revealed a synergistic inhibitory effect of Dmrt1 and androgen signaling on fancl transcription. Furthermore, an interaction between Fancl and the apoptotic factor Tumour protein p53 (Tp53) was found in vitro. The interaction between Fancl and Tp53 was observed via the WD repeat domain (WDR) and C-terminal domain (CTD) of Fancl and the DNA binding domain (DBD) of Tp53, leading to the K48-linked polyubiquitination degradation of Tp53 activated by the ubiquitin ligase, Fancl. Our results show that testis fate in cyp17a1-/- fish is determined by Dmrt1, which is thought to stabilize Tp53 by inhibiting fancl transcription during the critical stage of sexual fate determination in zebrafish.

Genomic Anatomy of Homozygous XX Females and YY Males Reveals Early Evolutionary Trajectory of Sex-determining Gene and Sex Chromosomes in Silurus Fishes.

Sex chromosomes display remarkable diversity and variability among vertebrates. Compared with research on the X/Y and Z/W chromosomes, which have long evolutionary histories in mammals and birds, studies on the sex chromosomes at early evolutionary stages are limited. Here, we precisely assembled the genomes of homozygous XX female and YY male Lanzhou catfish (Silurus lanzhouensis) derived from an artificial gynogenetic family and a self-fertilized family, respectively. Chromosome 24 (Chr24) was identified as the sex chromosome based on resequencing data. Comparative analysis of the X and Y chromosomes showed an approximate 320 kb Y-specific region with a Y-specific duplicate of anti-Mullerian hormone type II receptor (amhr2y), which is consistent with findings in 2 other Silurus species but on different chromosomes (Chr24 of Silurus meridionalis and Chr5 of Silurus asotus). Deficiency of amhr2y resulted in male-to-female sex reversal, indicating that amhr2y plays a male-determining role in S. lanzhouensis. Phylogenetic analysis and comparative genomics revealed that the common sex-determining gene amhr2y was initially translocated to Chr24 of the Silurus ancestor along with the expansion of transposable elements. Chr24 was maintained as the sex chromosome in S. meridionalis and S. lanzhouensis, whereas a sex-determining region transition triggered sex chromosome turnover from Chr24 to Chr5 in S. asotus. Additionally, gene duplication, translocation, and degeneration were observed in the Y-specific regions of Silurus species. These findings present a clear case for the early evolutionary trajectory of sex chromosomes, including sex-determining gene origin, repeat sequence expansion, gene gathering and degeneration in sex-determining region, and sex chromosome turnover.

A zebrafish NLRX1 isoform downregulates fish IFN responses by targeting the adaptor STING.

In mammals, NLRX1 is a unique member of the nucleotide-binding domain and leucine-rich repeat (NLR) family showing an ability to negatively regulate IFN antiviral immunity. Intron-containing genes, including NLRX1, have more than one transcript due to alternative splicing; however, little is known about the function of its splicing variants. Here, we identified a transcript variant of NLRX1 in zebrafish (Danio rerio), termed NLRX1-tv4, as a negative regulator of fish IFN response. Zebrafish NLRX1-tv4 was slightly induced by viral infection, with an expression pattern similar to the full-length NLRX1. Despite the lack of an N-terminal domain that exists in the full-length NLRX1, overexpression of NLRX1-tv4 still impaired fish IFN antiviral response and promoted viral replication in fish cells, similar to the full-length NLRX1. Mechanistically, NLRX1-tv4 targeted STING for proteasome-dependent protein degradation by recruiting an E3 ubiquitin ligase RNF5 to drive the K48-linked ubiquitination, eventually downregulating the IFN antiviral response. Mapping of NLRX1-tv4 domains showed that its N-terminal and C-terminal regions exhibited a similar potential to inhibit STING-mediated IFN antiviral response. Our findings reveal that like the full-length NLRX1, zebrafish NLRX-tv4 functions as an inhibitor to shape fish IFN antiviral response.IMPORTANCEIn this study, we demonstrate that a transcript variant of zebrafish NLRX1, termed NLRX1-tv4, downregulates fish IFN response and promotes virus replication by targeting STING for protein degradation and impairing the interaction of STING and TBK1 and that its N- and C-terminus exhibit a similar inhibitory potential. Our results are helpful in clarifying the current contradictory understanding of structure and function of vertebrate NLRX1s.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish.

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.

Fish female-biased gene cyp19a1a leads to female antiviral response attenuation between sexes by autophagic degradation of MITA.

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.

Establishment and characterization of a golden pompano (Trachinotus blochii) fin cell line for applications in marine fish pathogen immunology.

Pompano fishes have been widely farmed worldwide. As a representative commercial marine species of the Carangidae family, the golden pompano (Trachinotus blochii) has gained significant popularity in China and worldwide. However, because of rapid growth and high-density aquaculture, the golden pompano has become seriously threatened by various diseases. Cell lines are the most cost-effective resource for in vitro studies and are widely used for physiological and pathological research owing to their accessibility and convenience. In this study, we established a novel immortal cell line, GPF (Golden pompano fin cells). GPF has been passaged over 69 generations for 10 months. The morphology, adhesion and extension processes of GPF were evaluated using light and electron microscopy. GPF cells were passaged every 3 days with L-15 containing 20 % fetal bovine serum (FBS) at 1:3. The optimum conditions for GPF growth were 28 °C and a 20 % FBS concentration. DNA sequencing of 18S rRNA and mitochondrial 16S rRNA confirmed that GPF was derived from the golden pompano. Chromosomal analysis revealed that the number pattern of GPF was 48 chromosomes. Transfection experiments demonstrated that GPF could be utilized to express foreign genes. Furthermore, heavy metals (Cd, Cu, and Fe) exhibited dose-dependent cytotoxicity against GPF. After polyinosinic-polycytidylic acid (poly I:C) treatment, transcription of the retinoic acid-inducible gene I-like receptor (RLR) pathway genes, including mda5, mita, tbk1, irf3, and irf7 increased, inducing the expression of interferon (IFN) and anti-viral proteins in GPF cells. In addition, lipopolysaccharide (LPS) stimulation up-regulated the expression of inflammation-related factors, including myd88, irak1, nfκb, il1ß, il6, and cxcl10 expression. To the best of our knowledge, this is the first study on the immune response signaling pathways of the golden pompano using an established fin cell line. In this study, we describe a preliminary investigation of the GPF cell line immune response to poly I:C and LPS, and provide a more rapid and efficient experimental material for research on marine fish immunology.

A Novel Non-Mammalian-Specific HERC7 Negatively Regulates IFN Response through Degrading RLR Signaling Factors.

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.

A finTRIM Family Protein Acquires RNA-Binding Activity and E3 Ligase Activity to Shape the IFN Response in Fish.

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.

Genomic anatomy of male-specific microchromosomes in a gynogenetic fish.

Unisexual taxa are commonly considered short-lived as the absence of meiotic recombination is supposed to accumulate deleterious mutations and hinder the creation of genetic diversity. However, the gynogenetic gibel carp (Carassius gibelio) with high genetic diversity and wide ecological distribution has outlived its predicted extinction time of a strict unisexual reproduction population. Unlike other unisexual vertebrates, males associated with supernumerary microchromosomes have been observed in gibel carp, which provides a unique system to explore the rationales underlying male occurrence in unisexual lineage and evolution of unisexual reproduction. Here, we identified a massively expanded satellite DNA cluster on microchromosomes of hexaploid gibel carp via comparing with the ancestral tetraploid crucian carp (Carassius auratus). Based on the satellite cluster, we developed a method for single chromosomal fluorescence microdissection and isolated three male-specific microchromosomes in a male metaphase cell. Genomic anatomy revealed that these male-specific microchromosomes contained homologous sequences of autosomes and abundant repetitive elements. Significantly, several potential male-specific genes with transcriptional activity were identified, among which four and five genes displayed male-specific and male-biased expression in gonads, respectively, during the developmental period of sex determination. Therefore, the male-specific microchromosomes resembling common features of sex chromosomes may be the main driving force for male occurrence in gynogenetic gibel carp, which sheds new light on the evolution of unisexual reproduction.

Gonadal transcriptomes reveal sex-biased expression genes associated with sex determination and differentiation in red-tail catfish (Hemibagrus wyckioides).

BACKGROUND: Red-tail catfish (Hemibagrus wyckioides) is an important commercially farmed catfish in southern China. Males of red-tail catfish grow faster than females, suggesting that all-male catfish will produce more significant economic benefits in aquaculture practice. However, little research has been reported on sex determination and gonadal development in red-tail catfish. RESULTS: In this study, we performed the first transcriptomic analysis of male and female gonads at four developmental stages at 10, 18, 30, and 48 days post hatching (dph) using RNA-seq technology. A total of 23,588 genes were screened in 24 sequenced samples, of which 28, 213, 636, and 1381 differentially expressed genes (DEGs) were detected at four developmental stages, respectively. Seven candidate genes of sex determination and differentiation were further identified. Real-time quantitative PCR (RT-qPCR) further confirmed that anti-Mullerian hormone (amh), growth differentiation factor 6a (gdf6a), testis-specific gene antigen 10 (tsga10), and cytochrome P450 family 17 subfamily A (cyp17a) were highly expressed mainly in the male, while cytochrome P450 family 19 subfamily A polypeptide 1b (cyp19a1b), forkhead box L2 (foxl2), and hydroxysteroid 17-beta dehydrogenase 1 (hsd17b1) were highly expressed in the female. The KEGG pathway enrichment data showed that these identified DEGs were mainly involved in steroid hormone biosynthesis and TGF-ß signaling pathways. CONCLUSIONS: Based on RNA-seq data of gonads at the early developmental stages, seven DEGs shared by the four developmental stages were identified, among which amh and gdf6a may be the male-biased expression genes, while foxl2, cyp19a1b and hsd17b1 may be the female-biased expression genes in red-tail catfish. Our study will provide crucial genetic information for the research on sex control in red-tail catfish, as well as for exploring the evolutionary processes of sex determination mechanisms in fish.

Changes in Ploidy Drive Reproduction Transition and Genomic Diversity in a Polyploid Fish Complex.

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.

Interferon regulatory factor 1 (IRF1) and anti-pathogen innate immune responses.

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.

Molecular and functional characterization of golden pompano (Trachinotus blochii) TBK1 on IFN regulation.

The golden pompano (Trachinotus blochii), a pivotal commercial marine species in China, has gained significant popularity worldwide. However, accompanied with rapid growth and high density aquaculture, golden pompano has been seriously threatened by Nervous necrosis virus (NNV), while its molecular biology research regarding the innate immune system remains unexplored, which is crucial for understanding the activation of interferon (IFN) production and antiviral responses. In this study, we aimed to identify the characterization and function of golden pompano TANK-binding kinase 1 (gpTBK1), thereby providing evidence of the conservation of this classical factor in the RLR pathway among marine fish. Initially, we found the expression of gpTBK1 upregulation in diseased golden pompano with NNV infection and we successfully cloned the full-length open reading frame (ORF) of gpTBK1, consisting of 2172 nucleotides encoding 723 amino acids, from the head kidney. Subsequent analysis of the amino acid sequence revealed homology between gpTBK1 and other fish TBK1 proteins, with conserved N-terminal Serine/Threonine protein kinases catalytic domain (S_TKc) and C-terminal coiled coil domain (CCD). Moreover, the expression pattern showed that gpTBK1 exhibited ubiquitous expression across all evaluated tissues. Furthermore, functional identification experiments indicated that gpTBK1 activated interferon promoters' activity in golden pompano and induced the expression of downstream IFN-stimulated genes (ISGs). Notably, gpTBK1 was found to co-localize and interact with gpIRF3 in the cytoplasm. Collectively, these data provide a comprehensive analysis of the characterization and functional role of gpTBK1 in promoting interferon production. This research may facilitate the further study of the innate antiviral response, particularly the anti-NNV mechanisms, in golden pompano.

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism.

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex.

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Zebrafish HERC7c Acts as an Inhibitor of Fish IFN Response.

In humans, four small HERCs (HERC3-6) exhibit differential degrees of antiviral activity toward HIV-1. Recently we revealed a novel member HERC7 of small HERCs exclusively in non-mammalian vertebrates and varied copies of herc7 genes in distinct fish species, raising a question of what is the exact role for a certain fish herc7 gene. Here, a total of four herc7 genes (named HERC7a-d sequentially) are identified in the zebrafish genome. They are transcriptionally induced by a viral infection, and detailed promoter analyses indicate that zebrafish herc7c is a typical interferon (IFN)-stimulated gene. Overexpression of zebrafish HERC7c promotes SVCV (spring viremia of carp virus) replication in fish cells and concomitantly downregulates cellular IFN response. Mechanistically, zebrafish HERC7c targets STING, MAVS, and IRF7 for protein degradation, thus impairing cellular IFN response. Whereas the recently-identified crucian carp HERC7 has an E3 ligase activity for the conjugation of both ubiquitin and ISG15, zebrafish HERC7c only displays the potential to transfer ubiquitin. Considering the necessity for timely regulation of IFN expression during viral infection, these results together suggest that zebrafish HERC7c is a negative regulator of fish IFN antiviral response.

Functional Divergence of Multiple Duplicated Foxl2 Homeologs and Alleles in a Recurrent Polyploid Fish.

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.

Characterization of DNA Binding and Nuclear Retention Identifies Zebrafish IRF11 as a Positive Regulator of IFN Antiviral Response.

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.

Comparative mitogenome analyses uncover mitogenome features and phylogenetic implications of the subfamily Cobitinae.

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.

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

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