A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

Status of single-cell RNA sequencing for reproductive toxicology in zebrafish and the transcriptomic trade-off

The utilization of transcriptomic studies identifying profiles of gene expression, especially in toxicogenomics, has catapulted next-generation sequencing to the forefront of reproductive toxicology. An innovative yet underutilized RNA sequencing technique emerging into this field is single-cell RNA sequencing (scRNA-seq), which provides sequencing at the individual cellular level of gonad tissue. ScRNA-seq provides a novel and unique perspective for identifying distinct cellular profiles, including identification of rare cell subtypes. The specificity of scRNA-seq is a powerful tool for reproductive toxicity research, especially for translational animal models including zebrafish. Studies to date not only have focused on 'tissue atlassing' or characterizing what cell types make up different tissues but have also begun to include toxicant exposure as a factor that this review aims to explore. Future scRNA-seq studies will contribute to understanding exposure-induced outcomes; however, the trade-offs with traditional methods need to be considered.

Estrogen Signaling Inhibits the Expression of anti-Müllerian hormone (amh) and gonadal-soma-derived factor (gsdf) during the Critical Time of Sexual Fate Determination in Zebrafish

The mechanism of fish gonadal sex differentiation is complex and regulated by multiple factors. It has been widely known that proper steroidogenesis in Leydig cells and sex-related genes in Sertoli cells play important roles in gonadal sex differentiation. In teleosts, the precise interaction of these signals during the sexual fate determination remains elusive, especially their effect on the bi-potential gonad during the critical stage of sexual fate determination. Recently, all-testis phenotypes have been observed in the cyp17a1-deficient zebrafish and common carp, as well as in cyp19a1a-deficient zebrafish. By mating cyp17a1-deficient fish with transgenic zebrafish Tg(piwil1:EGFP-nanos3UTR), germ cells in the gonads were labelled with enhanced green fluorescent protein (EGFP). We classified the cyp17a1-deficient zebrafish and their control siblings into primordial germ cell (PGC)-rich and -less groups according to the fluorescence area of the EGFP labelling. Intriguingly, the EGFP-labelled bi-potential gonads in cyp17a1+/+ fish from the PGC-rich group were significantly larger than those of the cyp17a1-/- fish at 23 days post-fertilization (dpf). Based on the transcriptome analysis, we observed that the cyp17a1-deficient fish of the PGC-rich group displayed a significantly upregulated expression of amh and gsdf compared to that of control fish. Likewise, the upregulated expressions of amh and gsdf were observed in cyp19a1a-deficient fish as examined at 23 dpf. This upregulation of amh and gsdf could be repressed by treatment with an exogenous supplement of estradiol. Moreover, tamoxifen, an effective antagonist of both estrogen receptor α and β (ERα and Erβ), upregulates the expression of amh and gsdf in wild-type (WT) fish. Using the cyp17a1- and cyp19a1a-deficient zebrafish, we provide evidence to show that the upregulated expression of amh and gsdf due to the compromised estrogen signaling probably determines their sexual fate towards testis differentiation. Collectively, our data suggest that estrogen signaling inhibits the expression of amh and gsdf during the critical time of sexual fate determination, which may broaden the scope of sex steroid hormones in regulating gonadal sex differentiation in fish.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Significance of Premature Vertebral Mineralization in Zebrafish Models in Mechanistic and Pharmaceutical Research on Hereditary Multisystem Diseases

Zebrafish are increasingly becoming an important model organism for studying the pathophysiological mechanisms of human diseases and investigating how these mechanisms can be effectively targeted using compounds that may open avenues to novel treatments for patients. The zebrafish skeleton has been particularly instrumental in modeling bone diseases as-contrary to other model organisms-the lower load on the skeleton of an aquatic animal enables mutants to survive to early adulthood. In this respect, the axial skeletons of zebrafish have been a good read-out for congenital spinal deformities such as scoliosis and degenerative disorders such as osteoporosis and osteoarthritis, in which aberrant mineralization in humans is reflected in the respective zebrafish models. Interestingly, there have been several reports of hereditary multisystemic diseases that do not affect the vertebral column in human patients, while the corresponding zebrafish models systematically show anomalies in mineralization and morphology of the spine as their leading or, in some cases, only phenotype. In this review, we describe such examples, highlighting the underlying mechanisms, the already-used or potential power of these models to help us understand and amend the mineralization process, and the outstanding questions on how and why this specific axial type of aberrant mineralization occurs in these disease models.

PUBERTAL VIRILIZATION IN AN ADOLESCENT WITH 46, XY DISORDER OF SEXUAL DEVELOPMENT: A NOVEL MUTATION IN NR5A1 GENE

Background

NR5A1 [Steroidogenic factor 1 (SF1)] is a nuclear receptor that is essential for the development of gonads and adrenal glands as well as the establishment of steroidogenesis in these organs. The clinical findings of the mutations of NR5A1 gene in 46, XY individuals are variable. Virilization at puberty can be seen in some of the 46, XY children who have a female phenotype and are raised as female.A girl aged 13 years and 10 months old was brought by the family for deepening of her voice. On physical examination, her breast development was Tanner stage 2, axillary hair (+) and pubic hair was Tanner stage 4. She had labioscrotal fusion and 4.4 cm phallus (External Masculinisation Score was 6). Hypergonadotropic hypogonadism, low AMH and high testosterone levels were detected in laboratory tests. Uterus was not visualized in pelvic ultrasonography. Karyotype analysis was reported as 46, XY. Sequence analysis of the NR5A1 gene revealed a novel heterozygote c.1075_1089del (p.Leu359_Leu363del) variant. The patient was raised as a female and oestrogen replacement was started following gonadectomy.

Conclusion

It should be kept in mind that virilization may develop at puberty in individuals with 46, XY disorder of sexual development due to NR5A1 mutation.

Tandem Mass Tag-Based Quantitative Proteomics Analysis of Gonads Reveals New Insight into Sexual Reversal Mechanism in Chinese Soft-Shelled Turtles

Chinese soft-shelled turtles display obvious sex dimorphism. The exogenous application of hormones (estradiol and methyltestosterone) can change the direction of gonadal differentiation of P. sinensis to produce sex reversed individuals. However, the molecular mechanism remains unclear. In this study, TMT-based quantitative proteomics analysis of four types of P. sinensis (female, male, pseudo-female, and pseudo-male) gonads were compared. Quantitative analysis of 6107 labeled proteins in the four types of P. sinensis gonads was performed. We identified 440 downregulated and 423 upregulated proteins between pseudo-females and males, as well as 394 downregulated and 959 upregulated proteins between pseudo-males and females. In the two comparisons, the differentially expressed proteins, including K7FKG1, K7GIQ2, COL4A6, K7F2U2, and K7FF80, were enriched in some important pathways, such as focal adhesion, endocytosis, apoptosis, extracellular matrix-receptor interaction, and the regulation of actin cytoskeleton, which were upregulated in pseudo-female vs. male and downregulated in pseudo-male vs. female. In pathways such as ribosome and spliceosome, the levels of RPL28, SRSF3, SNRNP40, and HNRNPK were increased from male to pseudo-female, while they decreased from female to pseudo-male. All differentially expressed proteins after sexual reversal were divided into six clusters, according to their altered levels in the four types of P. sinensis, and associated with cellular processes, such as embryonic development and catabolic process, that were closely related to sexual reversal. These data will provide clues for the sexual reversal mechanism in P. sinensis.

The SARS-CoV-2 receptor and other key components of the Renin-Angiotensin-Aldosterone System related to COVID-19 are expressed in enterocytes in larval zebrafish

People with underlying conditions, including hypertension, obesity, and diabetes, are especially susceptible to negative outcomes after infection with coronavirus SARS-CoV-2, which causes COVID-19. Hypertension and respiratory inflammation are exacerbated by the Renin-Angiotensin-Aldosterone System (RAAS), which normally protects from rapidly dropping blood pressure via Angiotensin II (Ang II) produced by the enzyme Ace. The Ace paralog Ace2 degrades Ang II, counteracting its chronic effects, and serves as the SARS-CoV-2 receptor. Ace, the coronavirus, and COVID-19 comorbidities all regulate Ace2, but we do not yet understand how. To exploit zebrafish (Danio rerio) to help understand the relationship of the RAAS to COVID-19, we must identify zebrafish orthologs and co-orthologs of human RAAS genes and understand their expression patterns. To achieve these goals, we conducted genomic and phylogenetic analyses and investigated single cell transcriptomes. Results showed that most human RAAS genes have one or more zebrafish orthologs or co-orthologs. Results identified a specific type of enterocyte as the specific site of expression of zebrafish orthologs of key RAAS components, including Ace, Ace2, Slc6a19 (SARS-CoV-2 co-receptor), and the Angiotensin-related peptide cleaving enzymes Anpep (receptor for the common cold coronavirus HCoV-229E), and Dpp4 (receptor for the Middle East Respiratory Syndrome virus, MERS-CoV). Results identified specific vascular cell subtypes expressing Ang II receptors, apelin, and apelin receptor genes. These results identify genes and cell types to exploit zebrafish as a disease model for understanding mechanisms of COVID-19.