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

Abstract

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 dpf). Thus, forward genetic screens have identified very few genes required for gonad development. In addition, bulk RNA sequencing studies which identify genes expressed in the gonads do not have the resolution necessary to define minor cell populations that may play significant roles in 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.

Data availability

The raw and processed data reported in this paper are archived at NCBI GEO https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE191137) and in an interactively browsable forms at the Broad Institute Single-Cell Portal (https://singlecell.broadinstitute.org/single_cell/study/SCP928/40dpf-ovary-all-cells). Analysis code and objects are archived at github (https://github.com/yulongliu68/zeb_ov_ssRNAseq). Gene expression tables for the cell clusters identified are archived at Dryad: (https://datadryad.org/stash/share/CEd0Zs4oZKdinTWeJPKbWYjBq6hYq4QhVacQcFjf37E).

The following data sets were generated

Article and author information

Author details

  1. Yulong Liu

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
  2. Michelle E Kassack

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Matthew E McFaul

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Lana N Christensen

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Stefan Siebert

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
  6. Sydney R Wyatt

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
  7. Caramai N Kamei

    Mount Desert Island Biological Laboratory, Bar Harbor, United States
    Competing interests
    The authors declare that no competing interests exist.
  8. Samuel Horst

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
  9. Nayeli Arroyo

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
  10. Iain A Drummond

    Mount Desert Island Biological Laboratory, Bar Harbor, United States
    Competing interests
    The authors declare that no competing interests exist.
  11. Celina E Juliano

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4222-0987
  12. Bruce W Draper

    Department of Molecular and Cellular Biology, University of California, Davis, Davis, United States
    For correspondence
    bwdraper@ucdavis.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4397-7749

Funding

National Institutes of Health (R01 HD-081551)

  • Yulong Liu
  • Matthew E McFaul
  • Lana N Christensen
  • Bruce W Draper

National Science Foundation (IOS-1456737)

  • Michelle E Kassack
  • Lana N Christensen
  • Bruce W Draper

National Institutes of Health (T32 training grant ES-0070599)

  • Michelle E Kassack

National Institutes of Health (T32 training grant GM-007377)

  • Matthew E McFaul
  • Sydney R Wyatt

National Science Foundation (GRFP 2036201)

  • Sydney R Wyatt

National Institutes of Health (R35 GM133689)

  • Stefan Siebert
  • Celina E Juliano

National Institutes of Health (R01 DK126021)

  • Iain A Drummond

National Institutes of Health (R01)

  • Caramai N Kamei
  • Iain A Drummond

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Ethics

Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocols (#20200 and #20201) of the University of California, Davis.

Copyright

© 2022, Liu et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 6,525
    views
  • 1,007
    downloads
  • 63
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Yulong Liu
  2. Michelle E Kassack
  3. Matthew E McFaul
  4. Lana N Christensen
  5. Stefan Siebert
  6. Sydney R Wyatt
  7. Caramai N Kamei
  8. Samuel Horst
  9. Nayeli Arroyo
  10. Iain A Drummond
  11. Celina E Juliano
  12. Bruce W Draper
(2022)
Single-cell transcriptome reveals insights into the development and function of the zebrafish ovary
eLife 11:e76014.
https://doi.org/10.7554/eLife.76014

Share this article

https://doi.org/10.7554/eLife.76014

Further reading

    1. Developmental Biology
    Saira Amir, Olatunbosun Arowolo ... Alexander Suvorov
    Research Article

    Over the past several decades, a trend toward delayed childbirth has led to increases in parental age at the time of conception. Sperm epigenome undergoes age-dependent changes increasing risks of adverse conditions in offspring conceived by fathers of advanced age. The mechanism(s) linking paternal age with epigenetic changes in sperm remain unknown. The sperm epigenome is shaped in a compartment protected by the blood-testes barrier (BTB) known to deteriorate with age. Permeability of the BTB is regulated by the balance of two mTOR complexes in Sertoli cells where mTOR complex 1 (mTORC1) promotes the opening of the BTB and mTOR complex 2 (mTORC2) promotes its integrity. We hypothesized that this balance is also responsible for age-dependent changes in the sperm epigenome. To test this hypothesis, we analyzed reproductive outcomes, including sperm DNA methylation in transgenic mice with Sertoli cell-specific suppression of mTORC1 (Rptor KO) or mTORC2 (Rictor KO). mTORC2 suppression accelerated aging of the sperm DNA methylome and resulted in a reproductive phenotype concordant with older age, including decreased testes weight and sperm counts, and increased percent of morphologically abnormal spermatozoa and mitochondrial DNA copy number. Suppression of mTORC1 resulted in the shift of DNA methylome in sperm opposite to the shift associated with physiological aging – sperm DNA methylome rejuvenation and mild changes in sperm parameters. These results demonstrate for the first time that the balance of mTOR complexes in Sertoli cells regulates the rate of sperm epigenetic aging. Thus, mTOR pathway in Sertoli cells may be used as a novel target of therapeutic interventions to rejuvenate the sperm epigenome in advanced-age fathers.

    1. Cell Biology
    2. Developmental Biology
    Sarah Rubin, Ankit Agrawal ... Elazar Zelzer
    Research Article

    Chondrocyte columns, which are a hallmark of growth plate architecture, play a central role in bone elongation. Columns are formed by clonal expansion following rotation of the division plane, resulting in a stack of cells oriented parallel to the growth direction. In this work, we analyzed hundreds of Confetti multicolor clones in growth plates of mouse embryos using a pipeline comprising 3D imaging and algorithms for morphometric analysis. Surprisingly, analysis of the elevation angles between neighboring pairs of cells revealed that most cells did not display the typical stacking pattern associated with column formation, implying incomplete rotation of the division plane. Morphological analysis revealed that although embryonic clones were elongated, they formed clusters oriented perpendicular to the growth direction. Analysis of growth plates of postnatal mice revealed both complex columns, composed of ordered and disordered cell stacks, and small, disorganized clusters located in the outer edges. Finally, correlation between the temporal dynamics of the ratios between clusters and columns and between bone elongation and expansion suggests that clusters may promote expansion, whereas columns support elongation. Overall, our findings support the idea that modulations of division plane rotation of proliferating chondrocytes determines the formation of either clusters or columns, a multifunctional design that regulates morphogenesis throughout pre- and postnatal bone growth. Broadly, this work provides a new understanding of the cellular mechanisms underlying growth plate activity and bone elongation during development.