Transcriptome profiling of laser-captured germ cells and functional characterization of zbtb40 during 17alpha-methyltestosterone-induced spermatogenesis in orange-spotted grouper (Epinephelus coioides)

Synchronously sexual maturity in hermaphrodite fish as revealed by transcriptome analysis in Plectropomus leopardus

Leopard coral grouper (Plectropomus leopardus) is a type of hermaphrodite fish, but the mechanisms of gonadal development and gametogenesis remain unclear. In the present study, we performed histological observation and transcriptomic analysis during the process of sexual differentiation in P. leopardus. According to the histological results, sexual differentiation was completed at 15 months old, developed synchronously in male and female individuals at 2 years old, and matured synchronously at 3 years old. Comparative transcriptomic analyses showed that the gonadal had differentiated by 15 months old, with enrichment of pathways associated with cell proliferation, transcriptional metabolism, and germline stem cell differentiation. Furthermore, cilium movement and fatty acid anabolism, which are associated with spermatogenesis and oocyte growth, were significantly enriched at 3 years old. In addition, key genes associated with male and female sex differentiation, such as amh, dmrt1, dmrt2a, zp4, sox3, gdf9, and gsdf, were identified by weighted gene co-expression network analysis (WGCNA). Finally, the localization and expression of the key genes amh and sox3 were observed in different cell types within the testes and ovaries, reflecting the development of the testes and ovaries, respectively. All the evidence indicates that P. leopardus is a hermaphrodite and synchronously sexually mature fish. Our study complements the gonadal development patterns of hermaphroditic fish by providing new insights into the molecular mechanisms underlying sexual differentiation and sex change in hermaphroditic groupers.

Comparative genomic survey and functional analysis of DKKL1 during spermatogenesis in the Chinese soft-shelled turtle (Pelodiscus sinensis)

A feature of the Chinese soft-shelled turtle (Pelodiscus sinensis) is seasonal spermatogenesis; however, the underlying molecular mechanism is not well clarified. Here, we firstly cloned and characterized P. sinensis DKKL1, and then performed comparative genomic studies, expression analysis, and functional validation. P. sinensis DKKL1 had 2 putative N-glycosylation sites and 16 phosphorylation sites. DKKL1 also had classic transmembrane structures that were extracellularly localized. DKKL1's genetic distance was close to turtles, followed by amphibians and mammals, but its genetic distance was far from fishes. DKKL1 genes from different species shared distinct genomic characteristics. Meanwhile, they were also relatively conserved among themselves, at least from the perspective of classes. Notably, the transcription factors associated with spermatogenesis were also identified, containing CTCF, EWSR1, and FOXL2. DKKL1 exhibited sexually dimorphic expression only in adult gonads, which was significantly higher than that in other somatic tissues (P < 0.001), and was barely expressed in embryonic gonads. DKKL1 transcripts showed a strong signal in sperm, while faint signals were detected in other male germ cells. DKKL1 in adult testes progressively increased per month (P < 0.05), displaying a seasonal expression trait. DKKL1 was significantly downregulated in testes cells after the sex hormones (17β-estradiol and 17α-methyltestosterone) and Wnt/β-catenin inhibitor treatment (P < 0.05). Likewise, the Wnt/β-catenin inhibitor treatment dramatically repressed CTCF, EWSR1, and FOXL2 expression. Conversely, they were markedly upregulated after the 17β-estradiol and 17α-methyltestosterone treatment, suggesting that the three transcription factors might bind to different promoter regions, thereby negatively regulating DKKL1 transcription in response to the changes in the estrogen and androgen pathways, and positively controlling DKKL1 transcription in answer to the alterations in the Wnt/β-catenin pathway. Knockdown of DKKL1 significantly reduced the relative expression of HMGB2 and SPATS1 (P < 0.01), suggesting that it may be involved in seasonal spermatogenesis of P. sinensis through a positive regulatory interaction with these two genes. Overall, our findings provide novel insights into the genome evolution and potential functions of seasonal spermatogenesis of P. sinensis DKKL1.

ZBTB40 is a telomere-associated protein and protects telomeres in human ALT cells

Alternative lengthening of telomeres (ALTs) mechanism is activated in some somatic, germ cells, and human cancer cells. However, the key regulators and mechanisms of the ALT pathway remain elusive. Here we demonstrated that ZBTB40 is a novel telomere-associated protein and binds to telomeric dsDNA through its N-terminal BTB (BR-C, ttk and bab) or POZ (Pox virus and Zinc finger) domain in ALT cells. Notably, the knockout or knockdown of ZBTB40 resulted in the telomere dysfunction-induced foci and telomere lengthening in the ALT cells. The results also show that ZBTB40 is associated with ALT-associated promyelocytic leukemia nuclear bodies, and the loss of ZBTB40 induces the accumulation of the ALT-associated promyelocytic leukemia nuclear bodies in U2OS cells. Taken together, our results implicate that ZBTB40 is a key player of telomere protection and telomere lengthening regulation in human ALT cells.

Zbtb40 Deficiency Leads to Morphological and Phenotypic Abnormalities of Spermatocytes and Spermatozoa and Causes Male Infertility

Studies on the gene regulation of spermatogenesis are of unusual significance for maintaining male reproduction and treating male infertility. Here, we have demonstrated, for the first time, that a loss of ZBTB40 function leads to abnormalities in the morphological and phenotypic characteristics of mouse spermatocytes and spermatids as well as male infertility. We revealed that Zbtb40 was expressed in spermatocytes of mouse testes, and it was co-localized with γH2AX in mouse secondary spermatocytes. Interestingly, spermatocytes of Zbtb40 knockout mice had longer telomeres, compromised double-strand break (DSB) repair in the sex chromosome, and a higher apoptosis ratio compared to wild-type (WT) mice. The testis weight, testicular volume, and cauda epididymis body weight of the Zbtb40^+/- male mice were significantly lower than in WT mice. Mating tests indicated that Zbtb40^+/- male mice were able to mate normally, but they failed to produce any pups. Notably, sperm of Zbtb40^+/- mice showed flagellum deformities and abnormal acrosome biogenesis. Furthermore, a ZBTB40 mutation was associated with non-obstructive azoospermia. Our results implicate that ZBTB40 deficiency leads to morphological and phenotypic abnormalities of spermatocytes and spermatids and causes male infertility. This study thus offers a new genetic mechanism regulating mammalian spermatogenesis and provides a novel target for gene therapy in male infertility.

Effects of 17α-Methyltestosterone on the Transcriptome and Sex Hormones in the Brain of Gobiocypris rarus

17α-Methyltestosterone (MT), a synthetic environmental endocrine disruptor with androgenic effects, has been shown to disrupt the reproductive system and inhibit germ cell maturation in Gobiocypris rarus. To further investigate the regulation of gonadal development by MT through the hypothalamic-pituitary-gonadal (HPG) axis, G. rarus were exposed to 0, 25, 50, and 100 ng/L of MT for 7, 14, and 21 days. We analyzed its biological indicators, gonadotropin-releasing hormone (GnRH), gonadotropins, reproduction-related gene expression, and brain tissue transcriptome profiles. We found a significant decrease in the gonadosomatic index (GSI) in G. rarus males exposed to MT for 21 days compared to the control group. GnRH, follicle-stimulating hormone (FSH), and luteinizing hormone (LH) levels, as well as the expressions of the gnrh3, gnrhr1, gnrhr3, fshβ, and cyp19a1b genes, were significantly reduced in the brains of both male and female fish when exposed to 100 ng/L MT for 14 days compared to the controls. Therefore, we further constructed four RNA-seq libraries from 100 ng/L MT-treated groups of male and female fish, obtaining 2412 and 2509 DEGs in male and female brain tissue, respectively. Three common pathways were observed to be affected in both sexes after exposure to MT, namely, nicotinate and nicotinamide metabolism, focal adhesion, and cell adhesion molecules. Furthermore, we found that MT affected the PI3K/Akt/FoxO3a signaling pathway through the upregulation of foxo3 and ccnd2, and the downregulation of pik3c3 and ccnd1. Therefore, we hypothesize that MT interferes with the levels of gonadotropin-releasing hormone (GnRH, FSH, and LH) in G. rarus brains through the PI3K/Akt/FoxO3a signaling pathway, and affects the expression of key genes in the hormone production pathway (gnrh3, gnrhr1 and cyp19a1b) to interfere with the stability of the HPG axis, thus leading to abnormal gonadal development. This study provides a multidimensional perspective on the damaging effects of MT on fish and confirms that G. rarus is a suitable model animal for aquatic toxicology.

Recombinant Vibrio harveyi flagellin A protein and partial deletions of middle variable region and D0 domain induce immune related genes in Epinephelus coioides and Cyprinus carpio

Vibrio harveyi is a Gram-negative bacterium that causes vibriosis in various aquaculture species, including the orange-spotted grouper (Epinephelus coioides). Bacterial flagellin is a potent pathogen-associated molecule that stimulates the innate and adaptive immune systems through toll-like receptor 5 (TLR5) signaling. In this study, we isolated V. harveyi flagellin A (VhFliA) gene from V. harveyi (originated from orange-spotted grouper) and investigated the in vivo activities of recombinant VhFliA protein. Multiple sequence alignment showed that the amino acid sequence of VhFliA has conserved domains of N- and C-terminals (D0 and D1) and a middle variable (MV) region. We produced the VhFliA recombinant protein (wild type (WT)-VhFliA) by Escherichia coli and investigated its in vivo biological activity. Additionally, we prepared the VhFliA recombinant proteins with deletion of domains (ΔMV-VhFliA and ΔD0MV-VhFliA) to identify the domain for biological activity in the orange-spotted grouper. WT and ΔMV-VhFliA induced the expression of inflammatory cytokines (IFNγ, IL-1β, and IL-8) in groupers. However, ΔD0MV-VhFliA did not induce the expression of inflammatory cytokines. Additionally, to demonstrate the applicability of recombinant VhFliA to teleost species, we performed an in vivo assay of the recombinant proteins in koi carp (Cyprinus carpio). WT-VhFliA stimulates the expression of inflammatory cytokines (IL-1β, IL-6, and IL-8) in carp. ΔMV-VhFliA did not upregulate IL-1β and IL-6, whereas ΔD0MV-VhFliA induced expression in carp. These findings showed the potential of VhFliA as an effective immune stimulant adjuvant and comparative studies of flagellin - TLR5 signaling in teleosts.

Transcriptome sequencing reveals the effects of circRNA on testicular development and spermatogenesis in Qianbei Ma goats

Circular RNAs (circRNAs) play an important role in regulating the mammalian reproductive system, especially testicular development and spermatogenesis. However, their functions in testicular development and spermatogenesis in the Qianbei Ma goat, the Guizhou endemic breed are still unclear. In this study, tissue sectioning and circRNAs transcriptome analysis were conducted to compare the changes of morphology and circular RNAs gene expression profile at four different developmental stages (0Y, 0-month-old; 6Y, 6-month-old; 12Y, 12-month-old; 18Y, 18-month-old). The results showed that the circumferences and area of the seminiferous tubule gradually increased with age, and the lumen of the seminiferous tubule in the testis differentiated significantly. 12,784 circRNAs were detected from testicular tissues at four different developmental stages by RNA sequencing, and 8,140 DEcircRNAs (differentially expressed circRNAs) were found in 0Y vs. 6Y, 6Y vs. 12Y, 12Y vs. 18Y and 0Y vs. 18Y, 0Y vs. 12Y, 6Y vs. 18Y Functional enrichment analysis of the source genes showed that they were mainly enriched in testicular development and spermatogenesis. In addition, the miRNAs and mRNAs associated with DECircRNAs in 6 control groups were predicted by bioinformatics, and 81 highly expressed DECircRNAs and their associated miRNAs and mRNAs were selected to construct the ceRNA network. Through functional enrichment analysis of the target genes of circRNAs in the network, some candidate circRNAs related to testicular development and spermatogenesis were obtained. Such as circRNA_07172, circRNA_04859, circRNA_07832, circRNA_00032 and circRNA_07510. These results will help to reveal the mechanism of circRNAs in testicular development and spermatogenesis, and also provide some guidance for goat reproduction.

Establishment of a Spermatogonial Stem Cell Line with Potential of Meiosis in a Hermaphroditic Fish, Epinephelus coioides

Spermatogonial stem cells (SSCs) are unique adult stem cells capable of self-renewal and differentiation into sperm. Grouper is a protogynous hermaphroditic fish farmed widely in the tropical and subtropical seas. In this study, we established an SSC line derived from adult testis of orange-spotted grouper, Epinephelus coioides. In the presence of basic fibroblast growth factor (bFGF) and leukemia inhibitory factor (LIF), the cells could be maintained with proliferation and self-renewal over 20 months and 120 passages under in vitro culture conditions. The cells exhibited strong alkaline phosphatase activity and the characteristics of SSCs with the expression of germ cell markers, including Vasa, Dazl, and Plzf, as well as the stem cell markers Nanog, Oct4, and Ssea1. Furthermore, the cultured cells could be induced by 11-ketotestosterone treatment to highly express the meiotic markers Rec8, Sycp3, and Dmc1, and produce some spherical cells, and even sperm-like cells with a tail. The findings of this study suggested that the cultured grouper SSC line would serve as an excellent tool to study the molecular mechanisms behind SSCs self-renewal and differentiation, meiosis during spermatogenesis, and sex reversal in hermaphroditic vertebrates. Moreover, this SSC line has great application value in grouper fish aquaculture, such as germ cell transplantation, genetic manipulation, and disease research.

Screening and Validation of p38 MAPK Involved in Ovarian Development of Brachymystax lenok

Brachymystax lenok (lenok) is a rare cold-water fish native to China that is of high meat quality. Its wild population has declined sharply in recent years, and therefore, exploring the molecular mechanisms underlying the development and reproduction of lenoks for the purposes of artificial breeding and genetic improvement is necessary. The lenok comparative transcriptome was analyzed by combining single molecule, real-time, and next generation sequencing (NGS) technology. Differentially expressed genes (DEGs) were identified in five tissues (head kidney, spleen, liver, muscle, and gonad) between immature [300 days post-hatching (dph)] and mature [three years post-hatching (ph)] lenoks. In total, 234,124 and 229,008 full-length non-chimeric reads were obtained from the immature and mature sequencing data, respectively. After NGS correction, 61,405 and 59,372 non-redundant transcripts were obtained for the expression level and pathway enrichment analyses, respectively. Compared with the mature group, 719 genes with significantly increased expression and 1,727 genes with significantly decreased expression in all five tissues were found in the immature group. Furthermore, DEGs and pathways involved in the endocrine system and gonadal development were identified, and p38 mitogen-activated protein kinases (MAPKs) were identified as potentially regulating gonadal development in lenok. Inhibiting the activity of p38 MAPKs resulted in abnormal levels of gonadotropin-releasing hormone, follicle-stimulating hormone, and estradiol, and affected follicular development. The full-length transcriptome data obtained in this study may provide a valuable reference for the study of gene function, gene expression, and evolutionary relationships in B. lenok and may illustrate the basic regulatory mechanism of ovarian development in teleosts.

Pou5f1 and Nanog Are Reliable Germ Cell-Specific Genes in Gonad of a Protogynous Hermaphroditic Fish, Orange-Spotted Grouper (Epinephelus coioides)

Pluripotency markers Pou5f1 and Nanog are core transcription factors regulating early embryonic development and maintaining the pluripotency and self-renewal of stem cells. Pou5f1 and Nanog also play important roles in germ cell development and gametogenesis. In this study, Pou5f1 (EcPou5f1) and Nanog (EcNanog) were cloned from orange-spotted grouper, Epinephelus coioides. The full-length cDNAs of EcPou5f1 and EcNanog were 2790 and 1820 bp, and encoded 475 and 432 amino acids, respectively. EcPou5f1 exhibited a specific expression in gonads, whereas EcNanog was expressed highly in gonads and weakly in some somatic tissues. In situ hybridization analyses showed that the mRNA signals of EcNanog and EcPou5f1 were exclusively restricted to germ cells in gonads. Likewise, immunohistofluorescence staining revealed that EcNanog protein was limited to germ cells. Moreover, both EcPou5f1 and EcNanog mRNAs were discovered to be co-localized with Vasa mRNA, a well-known germ cell maker, in male and female germ cells. These results implied that EcPou5f1 and EcNanog could be also regarded as reliable germ cell marker genes. Therefore, the findings of this study would pave the way for elucidating the mechanism whereby EcPou5f1 and EcNanog regulate germ cell development and gametogenesis in grouper fish, and even in other protogynous hermaphroditic species.

Single-Cell Atlas of Adult Testis in Protogynous Hermaphroditic Orange-Spotted Grouper, Epinephelus coioides

Spermatogenesis is a process of self-renewal and differentiation in spermatogonial stem cells. During this process, germ cells and somatic cells interact intricately to ensure long-term fertility and accurate genome propagation. Spermatogenesis has been intensely investigated in mammals but remains poorly understood with regard to teleosts. Here, we performed single-cell RNA sequencing of ~9500 testicular cells from the male, orange-spotted grouper. In the adult testis, we divided the cells into nine clusters and defined ten cell types, as compared with human testis data, including cell populations with characteristics of male germ cells and somatic cells, each of which expressed specific marker genes. We also identified and profiled the expression patterns of four marker genes (calr, eef1a, s100a1, vasa) in both the ovary and adult testis. Our data provide a blueprint of male germ cells and supporting somatic cells. Moreover, the cell markers are candidates that could be used for further cell identification.

Integration of ATAC-seq and RNA-seq Unravels Chromatin Accessibility during Sex Reversal in Orange-Spotted Grouper (Epinephelus coioides)

Chromatin structure plays a pivotal role in maintaining the precise regulation of gene expression. Accessible chromatin regions act as the binding sites of transcription factors (TFs) and cis-elements. Therefore, information from these open regions will enhance our understanding of the relationship between TF binding, chromatin status and the regulation of gene expression. We employed an assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and RNA-seq analyses in the gonads of protogynous hermaphroditic orange-spotted groupers during sex reversal to profile open chromatin regions and TF binding sites. We focused on several crucial TFs, including ZNF263, SPIB, and KLF9, and analyzed the networks of TF-target genes. We identified numerous transcripts exhibiting sex-preferred expression among their target genes, along with their associated open chromatin regions. We then investigated the expression patterns of sex-related genes as well as the mRNA localization of certain genes during sex reversal. We found a set of sex-related genes that—upon further study—might be identified as the sex-specific or cell-specific marker genes that trigger sex reversal. Moreover, we discovered the core genes (gnas, ccnb2, and cyp21a) of several pathways related to sex reversal that provide the guideposts for future study.