The evolution of ovarian somatic cells characterized by transcriptome and chromatin accessibility across rodents, monkeys, and humans

Abstract

The ovary plays a crucial role in the reproductive system of female mammals by producing mature oocytes through folliculogenesis. Non-human model organisms are extensively utilized in research on human ovarian biology, thus necessitating the investigation of conservation and divergence in molecular mechanisms across species. In this study, we employed integrative single-cell analysis of transcriptome and chromatin accessibility to identify the evolutionary conservation and divergence patterns of ovaries among humans, monkeys, mice, rats, and rabbits. Our analyses revealed that theca cells exhibited the most significant changes during evolution based on scRNA-seq and scATAC-seq datasets. Furthermore, we discovered common cis-regulatory architectures in theca cells across species by conducting joint analyses of scRNA-seq and scATAC-seq datasets. These findings have potential applications in non-human biomedical and genetic research to validate molecular mechanisms found in human organisms. Additionally, our investigation into non-coding genomic regions identified intergenic highly transcribed regions (igHTRs) that may contribute to the evolution of species-specific phenotypic traits. Overall, our study provides valuable insights into understanding the molecular characteristics of adult ovaries while offering new perspectives for studying human ovarian physiology and diseases.

Keywords

single-cell multi-omics / ovarian somatic cells / evolutionary conservation and divergence / theca cells

References

[1] Richards JS , Pangas SA . The ovary: basic biology and clinical implications. J Clin Invest 2010; 120: 963- 72.

[2] Baerwald AR , Adams GP , Pierson RA . Ovarian antral folliculogenesis during the human menstrual cycle: a review. Hum Reprod Update 2012; 18: 73- 91.

[3] Hsueh AJ , Kawamura K , Cheng Y , et al. Intraovarian control of early folliculogenesis. Endocr Rev 2015; 36: 1- 24.

[4] Cavalcanti GS , Carvalho KC , Ferreira CS , et al. Granulosa cells and follicular development: a brief review. Rev Assoc Med Bras (1992) 2023; 69: e20230175.

[5] Havelock JC , Rainey WE , Carr BR . Ovarian granulosa cell lines. Mol Cell Endocrinol 2004; 228: 67- 78.

[6] Skinner MK . Regulation of primordial follicle assembly and development. Hum Reprod Update 2005; 11: 461- 71.

[7] Richards JS , Ren YA , Candelaria N , et al. Ovarian follicular theca cell recruitment, differentiation, and impact on fertility: 2017 update. Endocr Rev 2018; 39: 1- 20.

[8] Young JM , McNeilly AS . Theca: the forgotten cell of the ovarian follicle. Reproduction 2010; 140: 489- 504.

[9] Knight PG , Glister C . TGF-beta superfamily members and ovarian follicle development. Reproduction 2006; 132: 191- 206.

[10] Ferrer MJ , Silva AF , Abruzzese GA , et al. Lipid metabolism and relevant disorders to female reproductive health. Curr Med Chem 2021; 28: 5625- 47.

[11] Seli E , Babayev E , Collins SC , et al. Minireview: metabolism of female reproduction: regulatory mechanisms and clinical implications. Mol Endocrinol 2014; 28: 790- 804.

[12] Chen M , Long X , Chen M , et al. Integration of single-cell transcriptome and chromatin accessibility of early gonads development among goats, pigs, macaques, and humans. Cell Rep 2022; 41: 111587.

[13] Fan X , Bialecka M , Moustakas I , et al. Single-cell reconstruction of follicular remodeling in the human adult ovary. Nat Commun 2019; 10: 3164.

[14] Guahmich NL , Man L , Wang J , et al. Human theca arises from ovarian stroma and is comprised of three discrete subtypes. Commun Biol 2023; 6: 7.

[15] Lengyel E , Li Y , Weigert M , et al. A molecular atlas of the human postmenopausal fallopian tube and ovary from single-cell RNA and ATAC sequencing. Cell Rep 2022; 41: 111838.

[16] Long X , Yang Q , Qian J , et al. Obesity modulates cell-cell interactions during ovarian folliculogenesis. iScience 2022; 25: 103627.

[17] Morris ME , Meinsohn M-C , Chauvin M , et al. A single-cell atlas of the cycling murine ovary. Elife 2022; 11: e77239.

[18] Niu W , Spradling AC . Two distinct pathways of pregranulosa cell differentiation support follicle formation in the mouse ovary. Proc Natl Acad Sci U S A 2020; 117: 20015- 26.

[19] Wamaitha SE , Nie X , Pandolfi EC , et al. Single-cell analysis of the developing human ovary defines distinct insights into ovarian somatic and germline progenitors. Dev Cell 2023; 58: 2097- 111.e3.

[20] Wang JJ , Ge W , Zhai Q-Y , et al. Single-cell transcriptome landscape of ovarian cells during primordial follicle assembly in mice. PLoS Biol 2020; 18: e3001025.

[21] Wang S , Zheng Y , Li J , et al. Single-cell transcriptomic atlas of primate ovarian aging. Cell 2020; 180: 585- 600.e19.

[22] Jakimiuk AJ , Weitsman SR , Navab A , et al. Luteinizing hormone receptor, steroidogenesis acute regulatory protein, and steroidogenic enzyme messenger ribonucleic acids are overexpressed in thecal and granulosa cells from polycystic ovaries. J Clin Endocrinol Metab 2001; 86: 1318- 23.

[23] Ortega I , Cress AB , Wong DH , et al. Simvastatin reduces steroidogenesis by inhibiting Cyp17a1 gene expression in rat ovarian theca-interstitial cells. Biol Reprod 2012; 86: 1- 9.

[24] Wickenheisser JK , Nelson-Degrave VL , McAllister JM . Dysregulation of cytochrome P450 17alpha-hydroxylase messenger ribonucleic acid stability in theca cells isolated from women with polycystic ovary syndrome. J Clin Endocrinol Metab 2005; 90: 1720- 7.

[25] Adam M , Saller S , Ströbl S , et al. Decorin is a part of the ovarian extracellular matrix in primates and may act as a signaling molecule. Hum Reprod 2012; 27: 3249- 58.

[26] Rotgers E , Jorgensen A , Yao HH . At the crossroads of fatesomatic cell lineage specification in the fetal gonad. Endocr Rev 2018; 39: 739- 59.

[27] Harlow CR , Davidson L , Burns KH , et al. FSH and TGF-beta superfamily members regulate granulosa cell connective tissue growth factor gene expression in vitro and in vivo. Endocrinology 2002; 143: 3316- 25.

[28] Jorgez CJ , Klysik M , Jamin SP , et al. Granulosa cell-specific inactivation of follistatin causes female fertility defects. Mol Endocrinol 2004; 18: 953- 67.

[29] Nagaraja AK , Middlebrook BS , Rajanahally S , et al. Defective gonadotropin-dependent ovarian folliculogenesis and granulosa cell gene expression in inhibin-deficient mice. Endocrinology 2010; 151: 4994- 5006.

[30] Stein-O'Brien GL , Clark BS , Sherman T , et al. Decomposing cell identity for transfer learning across cellular measurements, platforms, tissues, and species. Cell Syst 2019; 8: 395- 411.e8.

[31] Wang J , Wu DC , Guo H , et al. Hyperandrogenemia and insulin resistance: the chief culprit of polycystic ovary syndrome. Life Sci 2019; 236: 116940.

[32] Li YE , Preissl S , Miller M , et al. A comparative atlas of single-cell chromatin accessibility in the human brain. Science 2023; 382: eadf7044.

[33] Morabito S , Miyoshi E , Michael N , et al. Single-nucleus chromatin accessibility and transcriptomic characterization of Alzheimer's disease. Nat Genet 2021; 53: 1143- 55.

[34] Klein JD , Qu C , Yang X , et al. c-Fos repression by Piwi regulates drosophila ovarian germline formation and tissue morphogenesis. PLoS Genet 2016; 12: e1006281.

[35] Ge T , Wen Y , Li B , et al. Single-cell sequencing reveals the reproductive variations between primiparous and multiparous Hu ewes. J Anim Sci Biotechnol 2023; 14: 144.

[36] Russell DL , Doyle KM , Gonzales-Robayna I , et al. Egr-1 induction in rat granulosa cells by follicle-stimulating hormone and luteinizing hormone: combinatorial regulation by transcription factors cyclic adenosine 3',5'-monophosphate regulatory element binding protein, serum response factor, sp1, and early growth response factor-1. Mol Endocrinol 2003; 17: 520- 33.

[37] Han P , Guerrero-Netro H , Estienne A , et al. Effects of fibroblast growth factors and the transcription factor, early growth response 1, on bovine theca cells. Mol Cell Endocrinol 2018; 476: 96- 102.

[38] Minnoye L , Taskiran II , Mauduit D , et al. Cross-species analysis of enhancer logic using deep learning. Genome Res 2020; 30: 1815- 34.

[39] LaVoie HA . Transcriptional control of genes mediating ovarian follicular growth, differentiation, and steroidogenesis in pigs. Mol Reprod Dev 2017; 84: 788- 801.

[40] Pusalkar M , Meherji P , Gokral J , et al. CYP11A1 and CYP17 promoter polymorphisms associate with hyperandrogenemia in polycystic ovary syndrome. Fertil Steril 2009; 92: 653- 9.

[41] Du H , Sarno J , Taylor HS . HOXA10 inhibits Kruppel-like factor 9 expression in the human endometrial epithelium. Biol Reprod 2010; 83: 205- 11.

[42] Baysoy A , Bai Z , Satija R , et al. The technological landscape and applications of single-cell multi-omics. Nat Rev Mol Cell Biol 2023; 24: 695- 713.

[43] Yan R , Gu C , You D , et al. Decoding dynamic epigenetic landscapes in human oocytes using single-cell multi-omics sequencing. Cell Stem Cell 2021; 28: 1641- 56.e7.

[44] Stuart T , Butler A , Hoffman P , et al. Comprehensive integration of single-cell data. Cell 2019; 177: 1888- 902.e21.

[45] Stuart T , Srivastava A , Madad S , et al. Single-cell chromatin state analysis with Signac. Nat Methods 2021; 18: 1333- 41.

[46] Wu T , Hu E , Xu S , et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation (Camb) 2021; 2: 100141.

[47] Granja JM , Corces MR , Pierce SE , et al. ArchR is a scalable software package for integrative single-cell chromatin accessibility analysis. Nat Genet 2021; 53: 403- 11.

[48] Zhang Y , Liu T , Meyer CA , et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol 2008; 9: R137.

RIGHTS & PERMISSIONS

The Author(s). Published by Oxford University Press on behalf of Higher Education Press.