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.

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.

Environmentally relevant uptake, elimination, and metabolic changes following early embryonic exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin in zebrafish.

Dioxin and dioxin-like compounds are ubiquitous environmental contaminants that induce toxicity by binding to the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor. The zebrafish model has been used to define the developmental toxicity observed following exposure to exogenous AHR ligands such as the potent agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (dioxin, TCDD). While the model has successfully identified cellular targets of TCDD and molecular mechanisms mediating TCDD-induced phenotypes, fundamental information such as the body burden produced by standard exposure models is still unknown. We performed targeted gas chromatography (GC) high-resolution mass spectrometry (HRMS) in tandem with non-targeted liquid chromatography (LC) HRMS to quantify TCDD uptake, model the elimination dynamics of TCDD, and determine how TCDD exposure affects the zebrafish metabolome. We found that 50 ppt, 10 ppb, and 1 ppb waterborne exposures to TCDD during early embryogenesis produced environmentally relevant body burdens: 38 ± 4.34, 26.6 ± 1.2, and 8.53 ± 0.341 pg/embryo, respectively, at 24 hours post fertilization. TCDD exposure was associated with the dysregulation of metabolic pathways that are associated with the AHR signaling pathway as well as pathways shown to be affected in mammals following TCDD exposure. In addition, we discovered that TCDD exposure affected several metabolic pathways that are critical for brain development and function including glutamate metabolism, chondroitin sulfate biosynthesis, and tyrosine metabolism. Together, these data demonstrate that existing exposure methods produce environmentally relevant body burdens of TCDD in zebrafish and provide insight into the biochemical pathways impacted by toxicant-induced AHR activation.

Proper modulation of AHR signaling is necessary for establishing neural connectivity and oligodendrocyte precursor cell development in the embryonic zebrafish brain.

2,3,7,8-tetrachlorodibenzo-[p]-dioxin (TCDD) is a persistent global pollutant that exhibits a high affinity for the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor. Epidemiological studies have associated AHR agonist exposure with multiple human neuropathologies. Consistent with the human data, research studies using laboratory models have linked pollutant-induced AHR activation to disruptions in learning and memory as well as motor impairments. Our understanding of endogenous AHR functions in brain development is limited and, correspondingly, scientists are still determining which cell types and brain regions are sensitive to AHR modulation. To identify novel phenotypes resulting from pollutant-induced AHR activation and ahr2 loss of function, we utilized the optically transparent zebrafish model. Early embryonic TCDD exposure impaired embryonic brain morphogenesis, resulted in ventriculomegaly, and disrupted neural connectivity in the optic tectum, habenula, cerebellum, and olfactory bulb. Altered neural network formation was accompanied by reduced expression of synaptic vesicle 2. Loss of ahr2 function also impaired nascent network development, but did not affect gross brain or ventricular morphology. To determine whether neural AHR activation was sufficient to disrupt connectivity, we used the Gal4/UAS system to express a constitutively active AHR specifically in differentiated neurons and observed disruptions only in the cerebellum; thus, suggesting that the phenotypes resulting from global AHR activation likely involve multiple cell types. Consistent with this hypothesis, we found that TCDD exposure reduced the number of oligodendrocyte precursor cells and their derivatives. Together, our findings indicate that proper modulation of AHR signaling is necessary for the growth and maturation of the embryonic zebrafish brain.

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

Zebrafish are an established research organism that has made many contributions to our understanding of vertebrate tissue and organ development, yet there are still significant gaps in our understanding of the genes that regulate gonad development, sex, and reproduction. Unlike the development of many organs, such as the brain and heart that form during the first few days of development, zebrafish gonads do not begin to form until the larval stage (≥5 days post-fertilization). Thus, forward genetic screens have identified very few genes required for gonad development. In addition, bulk RNA-sequencing studies that identify genes expressed in the gonads do not have the resolution necessary to define minor cell populations that may play significant roles in the development and function of these organs. To overcome these limitations, we have used single-cell RNA sequencing to determine the transcriptomes of cells isolated from juvenile zebrafish ovaries. This resulted in the profiles of 10,658 germ cells and 14,431 somatic cells. Our germ cell data represents all developmental stages from germline stem cells to early meiotic oocytes. Our somatic cell data represents all known somatic cell types, including follicle cells, theca cells, and ovarian stromal cells. Further analysis revealed an unexpected number of cell subpopulations within these broadly defined cell types. To further define their functional significance, we determined the location of these cell subpopulations within the ovary. Finally, we used gene knockout experiments to determine the roles of foxl2l and wnt9b for oocyte development and sex determination and/or differentiation, respectively. Our results reveal novel insights into zebrafish ovarian development and function, and the transcriptome profiles will provide a valuable resource for future studies.

Genetic regulation of sex determination and maintenance in zebrafish (Danio rerio).

Over the last several decades zebrafish (Danio rerio) has become a major model organism for the study of vertebrate development and physiology. Given this, it may be surprising how little is known about the mechanism that zebrafish use to determine sex. While zebrafish are a gonochoristic species (having two sexes) that do not switch sex as adults, it was appreciated early on that sex ratios obtained from breeding lab domesticated lines were not typically a 1:1 ratio of male and female, suggesting that sex was not determined by a strict chromosomal mechanism. Here we will review the recent progress toward defining the genetic mechanism for sex determination in both wild and domesticated zebrafish.

Female Sex Development and Reproductive Duct Formation Depend on Wnt4a in Zebrafish.

In laboratory strains of zebrafish, sex determination occurs in the absence of a typical sex chromosome and it is not known what regulates the proportion of animals that develop as males or females. Many sex determination and gonad differentiation genes that act downstream of a sex chromosome are well conserved among vertebrates, but studies that test their contribution to this process have mostly been limited to mammalian models. In mammals, WNT4 is a signaling ligand that is essential for ovary and Müllerian duct development, where it antagonizes the male-promoting FGF9 signal. Wnt4 is well conserved across all vertebrates, but it is not known if Wnt4 plays a role in sex determination and/or the differentiation of sex organs in nonmammalian vertebrates. This question is especially interesting in teleosts, such as zebrafish, because they lack an Fgf9 ortholog. Here we show that wnt4a is the ortholog of mammalian Wnt4, and that wnt4b was present in the last common ancestor of humans and zebrafish, but was lost in mammals. We show that wnt4a loss-of-function mutants develop predominantly as males and conclude that wnt4a activity promotes female sex determination and/or differentiation in zebrafish. Additionally, both male and female wnt4a mutants are sterile due to defects in reproductive duct development. Together these results strongly argue that Wnt4a is a conserved regulator of female sex determination and reproductive duct development in mammalian and nonmammalian vertebrates.