Single-cell transcriptome reveals insights into the development and function of the zebrafish ovary

Reproductive toxicity of DDT in the Japanese medaka fish model: Revisiting the impacts of DDT+ on female reproductive health

The organochlorine pesticide dichlorodiphenyltrichloroethane (DDT) is an endocrine-disrupting compound (EDC) that has been banned by most countries for decades. However, it continues to be detected in nearly all humans and wildlife due to its biological and environmental persistence. The ovarian dysgenesis syndrome hypothesis speculates that exposure to EDCs during sensitive developmental windows such as early gonadal differentiation lead to reproductive disorders later in life. Yet, mechanisms by which DDT affects developing gonads remain unclear due to the inherent challenge of getting developmental exposure data from adults presenting with reproductive disease. The Japanese medaka (Oryzias latipes) is a valuable fish model for sex-specific toxicological studies due to its chromosomal sex determination, external embryonic development, short generation time, and extensively mapped genome. It is well documented that medaka exposed to DDT and its metabolites and byproducts (herein referred to as DDT+) at different developmental time points experience permanent alterations in gonadal morphology, reproductive success, and molecular and hormonal signaling. However, the overwhelming majority of studies focus primarily on functional and morphological outcomes in males and females and have rarely investigated long-term transcriptional or molecular effects. This review summarizes previous experimental findings and the state of our knowledge concerning toxic effects DDT + on reproductive development, fertility, and health in the valuable medaka model. It also identifies gaps in knowledge, emphasizing a need for more focus on molecular mechanisms of ovarian endocrine disruption using enhanced molecular tools that have become increasingly available over the past few decades. Furthermore, DDT forms a myriad of over 45 metabolites and transformation products in biota and the environment, very few of which have been evaluated for environmental abundance or health effects. This reinforces the demand for high throughput and economical in vivo models for predictive toxicology screening, and the Japanese medaka is uniquely positioned to meet this need.

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.

The killifish germline regulates longevity and somatic repair in a sex-specific manner

Classical evolutionary theories propose tradeoffs among reproduction, damage repair and lifespan. However, the specific role of the germline in shaping vertebrate aging remains largely unknown. In this study, we used the turquoise killifish (Nothobranchius furzeri) to genetically arrest germline development at discrete stages and examine how different modes of infertility impact life history. We first constructed a comprehensive single-cell gonadal atlas, providing cell-type-specific markers for downstream phenotypic analysis. We show here that germline depletion-but not arresting germline differentiation-enhances damage repair in female killifish. Conversely, germline-depleted males instead showed an extension in lifespan and rejuvenated metabolic functions. Through further transcriptomic analysis, we highlight enrichment of pro-longevity pathways and genes in germline-depleted male killifish and demonstrate functional conservation of how these factors may regulate longevity in germline-depleted Caenorhabditis elegans. Our results, therefore, demonstrate that different germline manipulation paradigms can yield pronounced sexually dimorphic phenotypes, implying alternative responses to classical evolutionary tradeoffs.

The killifish germline regulates longevity and somatic repair in a sex-specific manner

Classical evolutionary theories propose tradeoffs between reproduction, damage repair, and lifespan. However, the specific role of the germline in shaping vertebrate aging remains largely unknown. Here, we use the turquoise killifish ( N. furzeri ) to genetically arrest germline development at discrete stages, and examine how different modes of infertility impact life-history. We first construct a comprehensive single-cell gonadal atlas, providing cell-type-specific markers for downstream phenotypic analysis. Next, we show that germline depletion - but not arresting germline differentiation - enhances damage repair in female killifish. Conversely, germline-depleted males instead showed an extension in lifespan and rejuvenated metabolic functions. Through further transcriptomic analysis, we highlight enrichment of pro-longevity pathways and genes in germline-depleted male killifish and demonstrate functional conservation of how these factors may regulate longevity in germline-depleted C. elegans . Our results therefore demonstrate that different germline manipulation paradigms can yield pronounced sexually dimorphic phenotypes, implying alternative responses to classical evolutionary tradeoffs.

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 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.

Dmrt1 is the only male pathway gene tested indispensable for sex determination and functional testis development in tilapia

Sex is determined by multiple factors derived from somatic and germ cells in vertebrates. We have identified amhy, dmrt1, gsdf as male and foxl2, foxl3, cyp19a1a as female sex determination pathway genes in Nile tilapia. However, the relationship among these genes is largely unclear. Here, we found that the gonads of dmrt1;cyp19a1a double mutants developed as ovaries or underdeveloped testes with no germ cells irrespective of their genetic sex. In addition, the gonads of dmrt1;cyp19a1a;cyp19a1b triple mutants still developed as ovaries. The gonads of foxl3;cyp19a1a double mutants developed as testes, while the gonads of dmrt1;cyp19a1a;foxl3 triple mutants eventually developed as ovaries. In contrast, the gonads of amhy;cyp19a1a, gsdf;cyp19a1a, amhy;foxl2, gsdf;foxl2 double and amhy;cyp19a1a;cyp19a1b, gsdf;cyp19a1a;cyp19a1b triple mutants developed as testes with spermatogenesis via up-regulation of dmrt1 in both somatic and germ cells. The gonads of amhy;foxl3 and gsdf;foxl3 double mutants developed as ovaries but with germ cells in spermatogenesis due to up-regulation of dmrt1. Taking the respective ovary and underdeveloped testis of dmrt1;foxl3 and dmrt1;foxl2 double mutants reported previously into consideration, we demonstrated that once dmrt1 mutated, the gonad could not be rescued to functional testis by mutating any female pathway gene. The sex reversal caused by mutation of male pathway genes other than dmrt1, including its upstream amhy and downstream gsdf, could be rescued by mutating female pathway gene. Overall, our data suggested that dmrt1 is the only male pathway gene tested indispensable for sex determination and functional testis development in tilapia.

miR-210 promotes immune- and suppresses oocyte meiosis-related genes in the zebrafish ovarian cells

microRNA-210 (miRNA), a well-documented miRNA, has been implicated in a myriad of biological processes, including responses to hypoxia, angiogenesis, cell proliferation, and male infertility in humans. However, a comprehensive understanding of its functions in fish requires further investigation. This study pursued to elucidate the downstream effect of dre-miR-210-5p on primary ovarian cell culture in zebrafish (Danio rerio), an animal model. A protocol was settled down by incubations with either an miR-210 mimic or a scrambled miRNA in the isolated ovaries. RNA-sequencing analysis identified ∼6000 differentially expressed target genes revealing that downregulated genes were associated with reproduction-related pathways while immune-related pathways displayed an upregulated pattern. To identify molecular markers, predicted target genes were classified into reproduction and immune cell types. These findings underscore the existence of a profound interplay between the reproductive and immune systems, with miR-210 emerging as a pivotal player in orchestrating transcriptomic alterations within fish ovaries.

Rbpms2 promotes female fate upstream of the nutrient sensing Gator2 complex component, Mios

Reproductive success relies on proper establishment and maintenance of biological sex. In many animals, including mammals, the primary gonad is initially ovary in character. We previously showed the RNA binding protein (RNAbp), Rbpms2, is required for ovary fate in zebrafish. Here, we identified Rbpms2 targets in oocytes (Rbpms2-bound oocyte RNAs; rboRNAs). We identify Rbpms2 as a translational regulator of rboRNAs, which include testis factors and ribosome biogenesis factors. Further, genetic analyses indicate that Rbpms2 promotes nucleolar amplification via the mTorc1 signaling pathway, specifically through the mTorc1-activating Gap activity towards Rags 2 (Gator2) component, Missing oocyte (Mios). Cumulatively, our findings indicate that early gonocytes are in a dual poised, bipotential state in which Rbpms2 acts as a binary fate-switch. Specifically, Rbpms2 represses testis factors and promotes oocyte factors to promote oocyte progression through an essential Gator2-mediated checkpoint, thereby integrating regulation of sexual differentiation factors and nutritional availability pathways in zebrafish oogenesis.

The Ancient Origin and Function of Germline Cysts

Gamete production in most animal species is initiated within an evolutionarily ancient multicellular germline structure, the germline cyst, whose interconnected premeiotic cells synchronously develop from a single progenitor arising just downstream from a stem cell. Cysts in mice, Drosophila, and many other animals protect developing sperm, while in females, cysts generate nurse cells that guard sister oocytes from transposons (TEs) and help them grow and build a Balbiani body. However, the origin and extreme evolutionary conservation of germline cysts remains a mystery. We suggest that cysts arose in ancestral animals like Hydra and Planaria whose multipotent somatic and germline stem cells (neoblasts) express genes conserved in all animal germ cells and frequently begin differentiation in cysts. A syncytial state is proposed to help multipotent stem cell chromatin transition to an epigenetic state with heterochromatic domains suitable for TE repression and specialized function. Most modern animals now lack neoblasts but have retained stem cells and cysts in their early germlines, which continue to function using this ancient epigenetic strategy.

Tilapia, a good model for studying reproductive endocrinology

The Nile tilapia (Oreochromis niloticus), with a system of XX/XY sex determination, is a worldwide farmed fish with a shorter sexual maturation time than that of most cultured fish. Tilapia show a spawning cycle of approximately 14 days and can be artificially propagated in the laboratory all year round to obtain genetically all female (XX) and all male (XY) fry. Its genome sequence has been opened, and a perfect gene editing platform has been established. With a moderate body size, it is convenient for taking enough blood to measure hormone level. In recent years, using tilapia as animal model, we have confirmed that estrogen is crucial for female development because 1) mutation of star2, cyp17a1 or cyp19a1a (encoding aromatase, the key enzyme for estrogen synthesis) results in sex reversal (SR) due to estrogen deficiency in XX tilapia, while mutation of star1, cyp11a1, cyp17a2, cyp19a1b or cyp11c1 affects fertility due to abnormal androgen, cortisol and DHP levels in XY tilapia; 2) when the estrogen receptors (esr2a/esr2b) are mutated, the sex is reversed from female to male, while when the androgen receptors are mutated, the sex cannot be reversed; 3) the differentiated ovary can be transdifferentiated into functional testis by inhibition of estrogen synthesis, and the differentiated testis can be transdifferentiated into ovary by simultaneous addition of exogenous estrogen and androgen synthase inhibitor; 4) loss of male pathway genes amhy, dmrt1, gsdf causes SR with upregulation of cyp19a1a in XY tilapia. Disruption of estrogen synthesis rescues the male to female SR of amhy and gsdf but not dmrt1 mutants; 5) mutation of female pathway genes foxl2 and sf-1 causes SR with downregulation of cyp19a1a in XX tilapia; 6) the germ cell SR of foxl3 mutants fails to be rescued by estrogen treatment, indicating that estrogen determines female germ cell fate through foxl3. This review also summarized the effects of deficiency of other steroid hormones, such as androgen, DHP and cortisol, on fish reproduction. Overall, these studies demonstrate that tilapia is an excellent animal model for studying reproductive endocrinology of 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.

Acute cold stress leads to zebrafish ovarian dysfunction by regulating miRNA and mRNA

Temperature is a critical factor that regulates the reproduction processes in teleost. However, the gonadal response mechanism to cold stress in fish remains largely unknown. In the present study, female zebrafish were exposed to different extents of low temperatures at 18 °C and 10 °C for 48 h. The ovarian histology was remarkably damaged after cold stress exposure. Integrated analysis of miRNA-mRNA was used to investigate the ovarian response to acute cold stress. A large number of mRNAs and miRNAs were altered by cold stress, which are involved in extensive biological processes. It is indicated that the signal transduction of MAPK and Calcium signaling pathway is highly engaged in zebrafish ovary to adapt to cold stress. The immune system was dysregulated by cold stress while the ovarian autophagy was activated. Remarkably increased gene number related to reproductive functions was identified in the cold stress at 10 °C compared to the control. The cold stress-induced dysregulated reproductive genes include star, hsd3b1, hsd17b1, inha, insl3, amh, nanos1 and foxl2. Combined with the dysregulated insulin, IGF and progesterone signaling, it is suggested that cold stress affects ovarian function in multiple aspects, including oocyte meiosis, folliculogenesis, final maturation and ovarian maintenance. On the other hand, the ovarian miRNA-mRNA regulatory network response to cold stress was also constructed. Overall, our result revealed the ovarian response to cold stress in zebrafish and provided insight into the fish adaptation mechanism to acute temperature change.

Macrophage activation drives ovarian failure and masculinization in zebrafish

BMP15 is a conserved regulator of ovarian development and maintenance in vertebrates. In humans, premature ovarian insufficiency is caused by autoimmunity and genetic factors, including mutation of BMP15. The cellular mechanisms underlying ovarian failure caused by BMP15 mutation and immune contributions are not understood. Using zebrafish, we established a causal link between macrophage activation and ovarian failure, which, in zebrafish, causes sex reversal. We define a germline-soma signaling axis that activates macrophages and drives ovarian failure and female-to-male sex reversal. Germline loss of zebrafish Bmp15 impairs oogenesis and initiates this cascade. Single-cell RNA sequencing and genetic analyses implicate ovarian somatic cells that express conserved macrophage-activating ligands as mediators of ovarian failure and sex reversal. Genetic ablation of macrophages or elimination of Csf1Rb ligands, Il34 or Csf1a, delays or blocks premature oocyte loss and sex reversal. The axis identified here provides insight into the cells and pathways governing oocyte and ovary maintenance and potential therapeutic targets to preserve female fertility.

Genetic evidence for differential functions of figla and nobox in zebrafish ovarian differentiation and folliculogenesis

FIGLA and NOBOX are important oocyte-specific transcription factors. Both figla-/- and nobox-/- mutants showed all-male phenotype in zebrafish due to increased dominance of the male-promoting pathway. The early diversion towards males in these mutants has precluded analysis of their roles in folliculogenesis. In this study, we attenuated the male-promoting pathway by deleting dmrt1, a key male-promoting gene, in figla-/- and nobox-/- fish, which allows a sufficient display of defects in folliculogenesis. Germ cells in figla-/-;dmrt1-/- double mutant remained in cysts without forming follicles. In contrast, follicles could form well but exhibited deficient growth in nobox-/-;dmrt1-/- double mutants. Follicles in nobox-/-;dmrt1-/- ovary could progress to previtellogenic (PV) stage but failed to enter vitellogenic growth. Such arrest at PV stage suggested a possible deficiency in estrogen signaling. This was supported by lines of evidence in nobox-/-;dmrt1-/-, including reduced expression of ovarian aromatase (cyp19a1a) and level of serum estradiol (E2), regressed genital papilla (female secondary sex characteristics), and more importantly the resumption of vitellogenic growth by E2 treatment. Expression analysis suggested Nobox might regulate cyp19a1a by controlling Gdf9 and/or Bmp15. Our discoveries indicate that Figla is essential for ovarian differentiation and follicle formation whereas Nobox is important for driving subsequent follicle development.

Defining the cellular complexity of the zebrafish bipotential gonad

Zebrafish are routinely used to model reproductive development, function, and disease, yet we still lack a clear understanding of the fundamental steps that occur during early bipotential gonad development, including when endothelial cells, pericytes, and macrophage arrive at the bipotential gonad to support gonad growth and differentiation. Here, we use a combination of transgenic reporters and single-cell sequencing analyses to define the arrival of different critical cell types to the larval zebrafish gonad. We determined that blood initially reaches the gonad via a vessel formed from the swim bladder artery, which we have termed the gonadal artery. We find that vascular and lymphatic development occurs concurrently in the bipotential zebrafish gonad and our data suggest that similar to what has been observed in developing zebrafish embryos, lymphatic endothelial cells in the gonad may be derived from vascular endothelial cells. We mined preexisting sequencing datasets to determine whether ovarian pericytes had unique gene expression signatures. We identified 215 genes that were uniquely expressed in ovarian pericytes, but not expressed in larval pericytes. Similar to what has been shown in the mouse ovary, our data suggest that pdgfrb+ pericytes may support the migration of endothelial tip cells during ovarian angiogenesis. Using a macrophage-driven photoconvertible protein, we found that macrophage established a nascent resident population as early as 12 dpf and can be observed removing cellular material during gonadal differentiation. This foundational information demonstrates that the early bipotential gonad contains complex cellular interactions, which likely shape the health and function of the mature gonad.

Gonacin: A germ cell-derived hormone with glucogenic, orexigenic, and gonadal activities

Fish require abundant nutrients to generate a large number of eggs for spawning. Based on the evolutionary conservation of human FBN2 and its C-terminal placensin-like sequences in fish, we identified a peptide hormone gonacin (GONAdal Cell placensIN) and found its high expression in early-stage germ cells in the ovary and testis of zebrafish. We demonstrated that gonacin is essential for food intake, glucose release, and ovarian development in zebrafish. Similar expression patterns and functions of gonacin were also demonstrated in rainbow trout. Gonacin represents the first hormone secreted by germ cells with endocrine functions in vertebrates, bridging the energy homeostasis and reproduction.

Germ cells do not progress through spermatogenesis in the infertile zebrafish testis

Vertebrate spermatogonial stem cells maintain sperm production over the lifetime of an animal but fertility declines with age. While morphological studies have greatly informed our understanding of typical spermatogenesis, the molecular and cellular mechanisms underlying spermatogenesis are not yet understood, particularly with respect to the onset of fertility. We used single-cell RNA sequencing to generate a developmental atlas of the zebrafish testis. Using 5 timepoints across the adult life of a zebrafish, we described cellular profiles in the testis during and after fertility. While all germ cell stages of spermatogenesis are detected in testes from fertile adult zebrafish, testes from older infertile males only contained spermatogonia and a reduced population of spermatocytes. These remaining germ cells are transcriptionally distinct from fertile spermatogonia. Immune cells including macrophages and lymphocytes drastically increase in abundance in infertile testes. Our developmental atlas reveals the cellular changes as the testis ages and defines a molecular roadmap for the regulation of male fertility.

Rescue of bmp15 deficiency in zebrafish by mutation of inha reveals mechanisms of BMP15 regulation of folliculogenesis

As an oocyte-specific growth factor, bone morphogenetic protein 15 (BMP15) plays a critical role in controlling folliculogenesis. However, the mechanism of BMP15 action remains elusive. Using zebrafish as the model, we created a bmp15 mutant using CRISPR/Cas9 and demonstrated that bmp15 deficiency caused a significant delay in follicle activation and puberty onset followed by a complete arrest of follicle development at previtellogenic (PV) stage without yolk accumulation. The mutant females eventually underwent female-to-male sex reversal to become functional males, which was accompanied by a series of changes in secondary sexual characteristics. Interestingly, the blockade of folliculogenesis and sex reversal in bmp15 mutant could be partially rescued by the loss of inhibin (inha-/-). The follicles of double mutant (bmp15-/-;inha-/-) could progress to mid-vitellogenic (MV) stage with yolk accumulation and the fish maintained their femaleness without sex reversal. Transcriptome analysis revealed up-regulation of pathways related to TGF-β signaling and endocytosis in the double mutant follicles. Interestingly, the expression of inhibin/activin βAa subunit (inhbaa) increased significantly in the double mutant ovary. Further knockout of inhbaa in the triple mutant (bmp15-/-;inha-/-;inhbaa-/-) resulted in the loss of yolk granules again. The serum levels of estradiol (E2) and vitellogenin (Vtg) both decreased significantly in bmp15 single mutant females (bmp15-/-), returned to normal in the double mutant (bmp15-/-;inha-/-), but reduced again significantly in the triple mutant (bmp15-/-;inha-/-;inhbaa-/-). E2 treatment could rescue the arrested follicles in bmp15-/-, and fadrozole (a nonsteroidal aromatase inhibitor) treatment blocked yolk accumulation in bmp15-/-;inha-/- fish. The loss of inhbaa also caused a reduction of Vtg receptor-like molecules (e.g., lrp1ab and lrp2a). In summary, the present study provided comprehensive genetic evidence that Bmp15 acts together with the activin-inhibin system in the follicle to control E2 production from the follicle, Vtg biosynthesis in the liver and its uptake by the developing oocytes.

Knock-out of vasotocin reduces reproductive success in female zebrafish, Danio rerio

The vertebrate nonapeptide vasotocin/vasopressin is evolutionarily highly conserved and acts as neuromodulator and endocrine/paracrine signaling molecule. Circumstantial and mechanistic evidence from pharmacological manipulations of the vasotocin system in several teleost fishes suggest sex- and species-specific reproductive roles of vasotocin. While effects of vasotocin on teleost reproductive physiology involve both courtship behaviors and the regulation of the hypothalamic-pituitary-gonadal (HPG) axes, comprehensive studies investigating behavioral and physiological reproductive consequences of genetic ablation of vasotocin in a genetically tractable fish model, such as the zebrafish, are currently lacking. Here, we report the generation of homozygous CRISPR/Cas9-based vasotocin gene knock-out zebrafish. Breeding pairs of vasotocin knock-out fish produce significantly fewer fertilized eggs per clutch compared to wildtype fish, an effect coincident with reduced female quivering courtship behavior. Crossbreeding experiments reveal that this reproductive phenotype is entirely female-dependent, as vasotocin-deficient males reproduce normally when paired with female wild-type fish. Histological analyses of vasotocin knock-out ovaries revealed an overall reduction in oocytes and differential distribution of oocyte maturation stages, demonstrating that the reproductive phenotype is linked to oocyte maturation and release. Ovarian hormone quantification and gene expression analysis in mutant fish indicated reduced synthesis of Prostaglandin F2α, a hormone involved in ovarian maturation, egg release and regulation of female courtship behavior in some cyprinids. However, acute injection of vasotocin did not rescue the female mutant reproductive phenotype, suggesting a contribution of organizational effects of vasotocin. Together, this study provides further support for emerging roles of vasotocin in female teleost reproduction in an important teleost model species.

Gadusol is a maternally provided sunscreen that protects fish embryos from DNA damage

Exposure to ultraviolet radiation (UVR) is harmful to living cells, leading organisms to evolve protective mechanisms against UVR-induced cellular damage and stress.1,2 UVR, particularly UVB (280-320 nm), can damage proteins and DNA, leading to errors during DNA repair and replication. Excessive UVR can induce cellular death. Aquatic organisms face risk of UV exposure as biologically harmful levels of UVB can penetrate >10 m in clear water.3 While melanin is the only known sunscreen in vertebrates, it often emerges late in embryonic development, rendering embryos of many species vulnerable during the earlier stages. Algae and microbes produce a class of sunscreening compounds known as mycosporine-like amino acids (MAAs).4 Fish eggs contain a similar compound called gadusol, whose role as a sunscreen has yet to be tested despite its discovery over 40 years ago.5 The recent finding that many vertebrate genomes contain a biosynthetic pathway for gadusol suggests that many fish may produce and use this molecule as a sunscreen.6 We generated a gadusol-deficient mutant zebrafish to investigate the role of gadusol in protecting fish embryos and larvae from UVR. Our results demonstrate that maternally provided gadusol is the primary sunscreen in embryonic and larval development, while melanin provides modest secondary protection. The gadusol biosynthetic pathway is retained in the vast majority of teleost genomes but is repeatedly lost in species whose young are no longer exposed to UVR. Our data demonstrate that gadusol is a maternally provided sunscreen that is critical for early-life survival in the most species-rich branch of the vertebrate phylogeny.

Defining the cellular complexity of the zebrafish bipotential gonad

Zebrafish are routinely used to model reproductive development, function, and disease, yet we still lack a clear understanding of the fundamental steps that occur during early bipotential gonad development, including when endothelial cells, pericytes, and macrophage cells arrive at the bipotential gonad to support gonad growth and differentiation. Here, we use a combination of transgenic reporters and single-cell sequencing analyses to define the arrival of different critical cell types to the larval zebrafish gonad. We determined that blood initially reaches the gonad via a vessel formed from the swim bladder artery, which we have termed the gonadal artery. We find that vascular and lymphatic development occurs concurrently in the bipotential zebrafish gonad and our data suggest that similar to what has been observed in developing zebrafish embryos, lymphatic endothelial cells in the gonad may be derived from vascular endothelial cells. We mined preexisting sequencing data sets to determine whether ovarian pericytes had unique gene expression signatures. We identified 215 genes that were uniquely expressed in ovarian pericytes that were not expressed in larval pericytes. Similar to what has been shown in the mouse ovary, our data suggest that pdgfrb+ pericytes may support the migration of endothelial tip cells during ovarian angiogenesis. Using a macrophage-driven photoconvertible protein, we found that macrophage established a nascent resident population as early as 12 dpf and can be observed removing cellular material during gonadal differentiation. This foundational information demonstrates that the early bipotential gonad contains complex cellular interactions, which likely shape the health and function of the mature, differentiated gonad.

Single cell transcriptome analyses of the developing zebrafish eye- perspectives and applications

Within a relatively short period of time, single cell transcriptome analyses (SCT) have become increasingly ubiquitous with transcriptomic research, uncovering plentiful details that boost our molecular understanding of various biological processes. Stemming from SCT analyses, the ever-growing number of newly assigned genetic markers increases our understanding of general function and development, while providing opportunities for identifying genes associated with disease. SCT analyses have been carried out using tissue from numerous organisms. However, despite the great potential of zebrafish as a model organism, other models are still preferably used. In this mini review, we focus on eye research as an example of the advantages in using zebrafish, particularly its usefulness for single cell transcriptome analyses of developmental processes. As studies have already shown, the unique opportunities offered by zebrafish, including similarities to the human eye, in combination with the possibility to analyze and extract specific cells at distinct developmental time points makes the model a uniquely powerful one. Particularly the practicality of collecting large numbers of embryos and therefore isolation of sufficient numbers of developing cells is a distinct advantage compared to other model organisms. Lastly, the advent of highly efficient genetic knockouts methods offers opportunities to characterize target gene function in a more cost-efficient way. In conclusion, we argue that the use of zebrafish for SCT approaches has great potential to further deepen our molecular understanding of not only eye development, but also many other organ systems.

Sex-biased gene expression at single-cell resolution: cause and consequence of sexual dimorphism

Gene expression differences between males and females are thought to be key for the evolution of sexual dimorphism, and sex-biased genes are often used to study the molecular footprint of sex-specific selection. However, gene expression is often measured from complex aggregations of diverse cell types, making it difficult to distinguish between sex differences in expression that are due to regulatory rewiring within similar cell types and those that are simply a consequence of developmental differences in cell-type abundance. To determine the role of regulatory versus developmental differences underlying sex-biased gene expression, we use single-cell transcriptomic data from multiple somatic and reproductive tissues of male and female guppies, a species that exhibits extensive phenotypic sexual dimorphism. Our analysis of gene expression at single-cell resolution demonstrates that nonisometric scaling between the cell populations within each tissue and heterogeneity in cell-type abundance between the sexes can influence inferred patterns of sex-biased gene expression by increasing both the false-positive and false-negative rates. Moreover, we show that, at the bulk level, the subset of sex-biased genes that are the product of sex differences in cell-type abundance can significantly confound patterns of coding-sequence evolution. Taken together, our results offer a unique insight into the effects of allometry and cellular heterogeneity on perceived patterns of sex-biased gene expression and highlight the power of single-cell RNA-sequencing in distinguishing between sex-biased genes that are the result of regulatory change and those that stem from sex differences in cell-type abundance, and hence are a consequence rather than a cause of sexual dimorphism.

RNA sequencing and expression analysis reveal a role for Lhx9 in the haploinsufficient adult mouse ovary

Understanding the molecular pathways that underpin ovarian development and function is vital for improving the research approaches to investigating fertility. Despite a significant improvement in our knowledge of molecular activity in the ovary, many questions remain unanswered in the quest to understand factors influencing fertility and ovarian pathologies such as cancer. Here, we present an investigation into the expression and function of the developmental transcription factor LIM Homeobox 9 (LHX9) in the adult mouse ovary. We have characterized Lhx9 expression in several cell types of the mature ovary across follicle stages. To evaluate possible LHX9 function in the adult ovary, we investigated ovarian anatomy and transcription in an Lhx9^+/- knockout mouse model displaying subfertility. Despite a lack of gross anatomical differences between genotypes, RNA-sequencing found that 90 differentially expressed genes between Lhx9^{+/ -} and Lhx9^+/+ mice. Gene ontology analyses revealed a reduced expression of genes with major roles in ovarian steroidogenesis and an increased expression of genes associated with ovarian cancer. Analysis of the ovarian epithelium revealed Lhx9^{+/ -} mice have a disorganized epithelial phenotype, corresponding to a significant increase in epithelial marker gene expression. These results provide an analysis of Lhx9 in the adult mouse ovary, suggesting a role in fertility and ovarian epithelial cancer.

The Evolution of Single-Cell RNA Sequencing Technology and Application: Progress and Perspectives

As an emerging sequencing technology, single-cell RNA sequencing (scRNA-Seq) has become a powerful tool for describing cell subpopulation classification and cell heterogeneity by achieving high-throughput and multidimensional analysis of individual cells and circumventing the shortcomings of traditional sequencing for detecting the average transcript level of cell populations. It has been applied to life science and medicine research fields such as tracking dynamic cell differentiation, revealing sensitive effector cells, and key molecular events of diseases. This review focuses on the recent technological innovations in scRNA-Seq, highlighting the latest research results with scRNA-Seq as the core technology in frontier research areas such as embryology, histology, oncology, and immunology. In addition, this review outlines the prospects for its innovative application in traditional Chinese medicine (TCM) research and discusses the key issues currently being addressed by scRNA-Seq and its great potential for exploring disease diagnostic targets and uncovering drug therapeutic targets in combination with multiomics technologies.

A chromosome-scale epigenetic map of the Hydra genome reveals conserved regulators of cell state

The epithelial and interstitial stem cells of the freshwater polyp Hydra are the best-characterized stem cell systems in any cnidarian, providing valuable insight into cell type evolution and the origin of stemness in animals. However, little is known about the transcriptional regulatory mechanisms that determine how these stem cells are maintained and how they give rise to their diverse differentiated progeny. To address such questions, a thorough understanding of transcriptional regulation in Hydra is needed. To this end, we generated extensive new resources for characterizing transcriptional regulation in Hydra, including new genome assemblies for Hydra oligactis and the AEP strain of Hydra vulgaris, an updated whole-animal single-cell RNA-seq atlas, and genome-wide maps of chromatin interactions, chromatin accessibility, sequence conservation, and histone modifications. These data revealed the existence of large kilobase-scale chromatin interaction domains in the Hydra genome that contain transcriptionally coregulated genes. We also uncovered the transcriptomic profiles of two previously molecularly uncharacterized cell types: isorhiza-type nematocytes and somatic gonad ectoderm. Finally, we identified novel candidate regulators of cell type-specific transcription, several of which have likely been conserved at least since the divergence of Hydra and the jellyfish Clytia hemisphaerica more than 400 million years ago.

wnt11f2 Zebrafish, an Animal Model for Development and New Insights in Bone Formation

Wnt signaling is a key regulator of osteoblast differentiation and mineralization in humans and animals, mediated by the canonical Wnt/β-catenin and noncanonical signaling pathways. Both pathways are crucial in regulating osteoblastogenesis and bone formation. The zebrafish silberblick (slb) carries a mutation in wnt11f2, a gene that contributes to embryonic morphogenesis; however, its role in bone morphology is unknown. wnt11f2 was originally known as wnt11; it was recently reclassified to avoid confusion in comparative genetics and disease modeling. The goal of this review is to summarize the characterization of the wnt11f2 zebrafish mutant and to deliver some new insights concerning its role in skeletal development. In addition to the previously described defects in early development in this mutant as well as craniofacial dysmorphia, we show an increase in tissue mineral density in the heterozygous mutant that points to a possible role of wnt11f2 in high bone mass phenotypes.

Single-Cell Transcriptomics Analysis Reveals a Cell Atlas and Cell Communication in Yak Ovary

Yaks (Bos grunniens) are the only bovine species that adapt well to the harsh high-altitude environment in the Qinghai-Tibetan plateau. However, the reproductive adaptation to the climate of the high elevation remains to be elucidated. Cell composition and molecular characteristics are the foundation of normal ovary function which determines reproductive performance. So, delineating ovarian characteristics at a cellular molecular level is conducive to elucidating the mechanism underlying the reproductive adaption of yaks. Here, the single-cell RNA-sequencing (scRNA-seq) was employed to depict an atlas containing different cell types with specific molecular signatures in the yak ovary. The cell types were identified on the basis of their specifically expressed genes and biological functions. As a result, a cellular atlas of yak ovary was established successfully containing theca cells, stromal cells, endothelial cells, smooth muscle cells, natural killer cells, macrophages, and proliferating cells. A cell-to-cell communication network between the distinct cell types was constructed. The theca cells were clustered into five subtypes based on their biological functions. Further, CYP11A1 was confirmed as a marker gene for the theca cells by immunofluorescence staining. Our work reveals an ovarian atlas at the cellular molecular level and contributes to providing insights into reproductive adaption in yaks.

Macrophage activation drives ovarian failure and masculinization

In humans, premature ovarian insufficiency (POI) is caused by autoimmunity and genetic factors, such as mutation of BMP15, a key ovarian determining gene. The cellular mechanisms associated with ovarian failure caused by BMP15 mutation and immune contributions to the disorder are not understood. BMP15's role in ovarian follicle development is conserved in vertebrates, including zebrafish. Using zebrafish, we established a causal link between macrophage activation and ovarian failure. We identified a germline-somatic gonadal cell-macrophage axis underlying ovarian atresia. Germline loss of Bmp15 triggers this axis that single-cell RNA sequencing and genetic analyses indicate involves activation of ovarian somatic cells that express conserved macrophage-activating ligands. Genetic ablation of macrophages blocks premature oocyte loss. Thus, the axis identified here represents potential therapeutic targets to preserve female fertility.

A single-cell transcriptomic atlas characterizes cell types and their molecular features in yak ovarian cortex

The ovary as one of the most dynamic organs produces steroids to orchestrate female secondary sexual characteristics, harbors ovarian reserve for oocytes, releases mature oocytes for fertilization, and maintains pregnancy. Yak (Bos grunniens) is the only bovid animal that can adapt to the harsh climatic conditions on the Qinghai-Tibetan Plateau (altitudes of over 3000 m above sea level). However, the cellular atlas is composed of oocytes and other somatic cells, and their individual molecular characteristics remain to be elucidated in the yak ovary. Here, single-cell RNA sequencing (scRNA-seq) was performed to delineate the molecular signature of various cell types in the yak ovarian cortex. A cellular atlas of yak ovarian cortex was constructed successfully on the basis of the differentially expressed genes (DEGs) from the distinct cell types and their functional enrichment analysis, comprising endothelial cells, nature kill cells, stromal cells, smooth muscle cells, oocytes, macrophages, epithelial cells, and granulosa cells. Meanwhile, the signature genes were determined based on their expression specificity in each cell type. A cell-to-cell communication network was built in light of the differentially overexpressed ligand and receptor genes from each cell type. Further, the oocytes were subdivided into four subtypes based on their individual DEGs and the functional enrichment of the DEGs. FST and TOP2A were identified as maker genes for oocytes by immunostaining in the yak ovarian cortex. The cellular atlas reveals the biological characteristics of the ovarian cortex at the cellular molecular level and provides insights into female reproductive biology via cellular communications in the yak.

Metzincin metalloproteases in PGC migration and gonadal sex conversion

Development of a functional gonad includes migration of primordial germ cells (PGCs), differentiations of somatic and germ cells, formation of primary follicles or spermatogenic cysts with somatic gonadal cells, development and maturation of gametes, and subsequent releasing of mature germ cells. These processes require extensive cellular and tissue remodeling, as well as broad alterations of the surrounding extracellular matrix (ECM). Metalloproteases, including MMPs (matrix metalloproteases), ADAMs (a disintegrin and metalloproteinases), and ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs), are suggested to have critical roles in the remodeling of the ECM during gonad development. However, few research articles and reviews are available on the functions and mechanisms of metalloproteases in remodeling gonadal ECM, gonadal development, or gonadal differentiation. Moreover, most studies focused on the roles of transcription and growth factors in early gonad development and primary sex determination, leaving a significant knowledge gap on how differentially expressed metalloproteases exert effects on the ECM, cell migration, development, and survival of germ cells during the development and differentiation of ovaries or testes. We will review gonad development with focus on the evidence of metalloprotease involvements, and with an emphasis on zebrafish as a model for studying gonadal sex differentiation and metalloprotease functions.

Dynamic regulation of gonadal transposon control across the lifespan of the naturally short-lived African turquoise killifish

Although germline cells are considered to be functionally "immortal," both the germline and supporting somatic cells in the gonad within an organism experience aging. With increased age at parenthood, the age-related decline in reproductive success has become an important biological issue for an aging population. However, molecular mechanisms underlying reproductive aging across sexes in vertebrates remain poorly understood. To decipher molecular drivers of vertebrate gonadal aging across sexes, we perform longitudinal characterization of the gonadal transcriptome throughout the lifespan in the naturally short-lived African turquoise killifish (Nothobranchius furzeri). By combining mRNA-seq and small RNA-seq from 26 individuals, we characterize the aging gonads of young-adult, middle-aged, and old female and male fish. We analyze changes in transcriptional patterns of genes, transposable elements (TEs), and piRNAs. We find that testes seem to undergo only marginal changes during aging. In contrast, in middle-aged ovaries, the time point associated with peak female fertility in this strain, PIWI pathway components are transiently down-regulated, TE transcription is elevated, and piRNA levels generally decrease, suggesting that egg quality may already be declining at middle-age. Furthermore, we show that piRNA ping-pong biogenesis declines steadily with age in ovaries, whereas it is maintained in aging testes. To our knowledge, this data set represents the most comprehensive transcriptomic data set for vertebrate gonadal aging. This resource also highlights important pathways that are regulated during reproductive aging in either ovaries or testes, which could ultimately be leveraged to help restore aspects of youthful reproductive function.

Genetic analysis of activin/inhibin β subunits in zebrafish development and reproduction

Activin and inhibin are both dimeric proteins sharing the same β subunits that belong to the TGF-β superfamily. They are well known for stimulating and inhibiting pituitary FSH secretion, respectively, in mammals. In addition, activin also acts as a mesoderm-inducing factor in frogs. However, their functions in development and reproduction of other species are poorly defined. In this study, we disrupted all three activin/inhibin β subunits (βAa, inhbaa; βAb, inhbab; and βB, inhbb) in zebrafish using CRISPR/Cas9. The loss of βAa/b but not βB led to a high mortality rate in the post-hatching stage. Surprisingly, the expression of fshb but not lhb in the pituitary increased in the female βA mutant together with aromatase (cyp19a1a) in the ovary. The single mutant of βAa/b showed normal folliculogenesis in young females; however, their double mutant (inhbaa-/-;inhbab-/-) showed delayed follicle activation, granulosa cell hypertrophy, stromal cell accumulation and tissue fibrosis. The ovary of inhbaa-/- deteriorated progressively after 180 dpf with reduced fecundity and the folliculogenesis ceased completely around 540 dpf. In addition, tumor- or cyst-like tissues started to appear in the inhbaa-/- ovary after about one year. In contrast to females, activin βAa/b mutant males showed normal spermatogenesis and fertility. As for activin βB subunit, the inhbb-/- mutant exhibited normal folliculogenesis, spermatogenesis and fertility in both sexes; however, the fecundity of mutant females decreased dramatically at 270 dpf with accumulation of early follicles. In summary, the activin-inhibin system plays an indispensable role in fish reproduction, in particular folliculogenesis and ovarian homeostasis.

A Germline-Specific Regulator of Mitochondrial Fusion is Required for Maintenance and Differentiation of Germline Stem and Progenitor Cells

Maintenance and differentiation of germline stem and progenitor cells (GSPCs) is important for sexual reproduction. Here, the authors identify zebrafish pld6 as a novel germline-specific gene by cross-analyzing different RNA sequencing results, and find that pld6 knockout mutants develop exclusively into infertile males. In pld6 mutants, GSPCs fail to differentiate and undergo apoptosis, leading to masculinization and infertility. Mitochondrial fusion in pld6-depleted GSPCs is severely impaired, and the mutants exhibit defects in piRNA biogenesis and transposon suppression. Overall, this work uncovers zebrafish Pld6 as a novel germline-specific regulator of mitochondrial fusion, and highlights its essential role in the maintenance and differentiation of GSPCs as well as gonadal development and gametogenesis.

Parental mutations influence wild-type offspring via transcriptional adaptation

Transgenerational epigenetic inheritance (TEI) is mostly discussed in the context of physiological or environmental factors. Here, we show intergenerational and transgenerational inheritance of transcriptional adaptation (TA), a process whereby mutant messenger RNA (mRNA) degradation affects gene expression, in nematodes and zebrafish. Wild-type offspring of animals heterozygous for mRNA-destabilizing alleles display increased expression of adapting genes. Notably, offspring of animals heterozygous for nontranscribing alleles do not display this response. Germline-specific mutations are sufficient to induce TA in wild-type offspring, indicating that, at least for some genes, mutations in somatic tissues are not necessary for this process. Microinjecting total RNA from germ cells of TA-displaying heterozygous zebrafish can trigger TA in wild-type embryos and in their progeny, suggesting a model whereby mutant mRNAs in the germline trigger a TA response that can be epigenetically inherited. In sum, this previously unidentified mode of TEI reveals a means by which parental mutations can modulate the offspring's transcriptome.