Single-cell RNA-seq identified novel genes involved in primordial follicle formation

Single-cell RNA-seq identified novel genes involved in primordial follicle formation Hang-Jing Tan Zi-Heng Deng Hui Shen Hong-Wen Deng Hong-Mei Xiao Introduction

The number of primordial follicles (PFs) in mammals determines the ovarian reserve, and impairment of primordial follicle formation (PFF) will cause premature ovarian insufficiency (POI).
Methods

By analyzing public single-cell RNA sequencing performed during PFF on mice and human ovaries, we identified novel functional genes and novel ligand-receptor interaction during PFF. Based on immunofluorescence and in vitro ovarian culture, we confirmed mechanisms of genes and ligand-receptor interaction in PFF. We also applied whole exome sequencing (WES) in 93 cases with POI and whole genome sequencing (WGS) in 465 controls. Variants in POI patients were further investigated by in silico analysis and functional verification.
Results

We revealed ANXA7 (annexin A7) and GTF2F1 (general transcription factor IIF subunit 1) in germ cells to be novel potentially genes in promoting PFF. Ligand Mdk (midkine) in germ cells and its receptor Sdc1 (syndecan 1) in granulosa cells are novel interaction crucial for PFF. Based on immunofluorescence, we confirmed significant up-regulation of ANXA7 in PFs compared with germline cysts, and uniform expression of GTF2F1, MDK and SDC1 during PFF, in 25 weeks human fetal ovary. In vitro investigation indicated that Anxa7 and Gtf2f1 are vital for mice PFF by regulating Jak/Stat3 and Jnk signaling pathways, respectively. Ligand-receptor (Mdk-Sdc1) are crucial for PFF by regulating Pi3k-akt signaling pathway. Two heterozygous variants in GTF2F1, and one heterozygous variants in SDC1 were identified in cases, but no variant were identified in controls. The protein level of GTF2F1 or SDC1 in POI cases are significantly lower than that of controls, indicating the pathogenic effects of the two genes on ovarian function were dosage dependent.
Discussion

Our study identified novel genes and novel ligand-receptor interaction during PFF, and further expanding the genetic architecture of POI.
12 11 2023 1285667 10.3389/fendo.2023.1285667 1 10.3389/crossmark-policy frontiersin.org true https://creativecommons.org/licenses/by/4.0/ 10.3389/fendo.2023.1285667 https://www.frontiersin.org/articles/10.3389/fendo.2023.1285667/full https://www.frontiersin.org/articles/10.3389/fendo.2023.1285667/full PloS Biol Wang 18 2020 Single-cell transcriptome landscape of ovarian cells during primordial follicle assembly in mice 10.1371/journal.pbio.3001025 Dev Biol Cui 383 2013 Cytoplasmic polyadenylation is a major mRNA regulator during oogenesis and egg activation in Drosophila 10.1016/j.ydbio.2013.08.013 Aging Cell He 20 2021 Single-cell RNA-Seq reveals a highly coordinated transcriptional program in mouse germ cells during primordial follicle formation 10.1111/acel.13424 Hum Reprod Update Skinner 11 2005 Regulation of primordial follicle assembly and development 10.1093/humupd/dmi020 J Assist Reprod Genet Kirshenbaum 36 2019 Premature ovarian insufficiency (POI) and autoimmunity-an update appraisal 10.1007/s10815-019-01572-0 FERTIL STERIL Qin 101 2014 Novel variants in the SOHLH2 gene are implicated in human premature ovarian failure 10.1016/j.fertnstert.2014.01.001 Am J Hum Genet Zhao 82 2008 Transcription factor FIGLA is mutated in patients with premature ovarian failure 10.1016/j.ajhg.2008.04.018 Proc Natl Acad Sci U.S.A. Niu 117 2020 Two distinct pathways of pregranulosa cell differentiation support follicle formation in the mouse ovary 10.1073/pnas.2005570117 Cell Stem Cell Li 20 2017 Single-cell RNA-seq analysis maps development of human germline cells and gonadal niche interactions 10.1016/j.stem.2017.03.007 Nat Biotechnol Butler 36 2018 Integrating single-cell transcriptomic data across different conditions, technologies, and species 10.1038/nbt.4096 Nat Methods Qiu 14 2017 Reversed graph embedding resolves complex single-cell trajectories 10.1038/nmeth.4402 BMC Bioinf Langfelder 9 559 2008 WGCNA: an R package for weighted correlation network analysis 10.1186/1471-2105-9-559 Nat Methods Aibar 14 2017 SCENIC: single-cell regulatory network inference and clustering 10.1038/nmeth.4463 Nat Commun Jin 12 1088 2021 Inference and analysis of cell-cell communication using CellChat 10.1038/s41467-021-21246-9 J Cachexia Sarcopenia Muscle Lv 12 2021 Human gut microbiome impacts skeletal muscle mass via gut microbial synthesis of the short-chain fatty acid butyrate among healthy menopausal women 10.1002/jcsm.12788 Nat Methods Schwarz 7 2010 MutationTaster evaluates disease-causing potential of sequence alterations 10.1038/nmeth0810-575 Nucleic Acids Res Rentzsch 47 2019 CADD: predicting the deleteriousness of variants throughout the human genome 10.1093/nar/gky1016 BIOINFORMATICS Quang 31 2015 DANN: a deep learning approach for annotating the pathogenicity of genetic variants 10.1093/bioinformatics/btu703 J Clin Endocrinol Metab Albader 107 2022 Insights of noncanonical splice-site variants on RNA splicing in patients with congenital hypothyroidism 10.1210/clinem/dgab737 BIOINFORMATICS Ramani 35 2019 PhastWeb: a web interface for evolutionary conservation scoring of multiple sequence alignments using phastCons and phyloP 10.1093/bioinformatics/bty966 PloS Genet Huber 16 2020 Population genetic models of GERP scores suggest pervasive turnover of constrained sites across mammalian evolution 10.1371/journal.pgen.1008827 Nucleic Acids Res Ng 31 2003 SIFT: Predicting amino acid changes that affect protein function 10.1093/nar/gkg509 Curr Protoc Hum Genet Adzhubei t7 2013 Predicting functional effect of human missense mutations using PolyPhen-2 10.1002/0471142905.hg0720s76 Sci Rep Lu 5 10576 2015 A statistical framework to predict functional non-coding regions in the human genome through integrated analysis of annotation data 10.1038/srep10576 Cell Death Dis Zhou 12 540 2021 Long non-coding RNA Xist regulates oocyte loss via suppressing miR-23b-3p/miR-29a-3p maturation and upregulating STX17 in perinatal mouse ovaries 10.1038/s41419-021-03831-4 Nucleic Acids Res Wang 48 2020 LHX8 and SOHLH1 transcription factor networks regulate mouse oocyte growth and differentiation 10.1093/nar/gkaa101 J Mol Cell Biol Cai 12 2020 SP1 governs primordial folliculogenesis by regulating pregranulosa cell development in mice 10.1093/jmcb/mjz059 Cell Death Dis Zhang 12 559 2021 HDAC6 regulates primordial follicle activation through mTOR signaling pathway 10.1038/s41419-021-03842-1 Development Niu 143 2016 JNK signaling regulates E-cadherin junctions in germline cysts and determines primordial follicle formation in mice 10.1242/dev.132175 Stem Cells Childs 28 2010 BMP signaling in the human fetal ovary is developmentally regulated and promotes primordial germ cell apoptosis 10.1002/stem.440 Mol Reprod Dev Workman 90 295 2023 RNA sequencing and expression analysis reveal a role for Lhx9 in the haploinsufficient adult mouse ovary 10.1002/mrd.23686 Biol Reprod Rusovici 69 64 2003 Expression and distribution of AP-1 transcription factors in the porcine ovary 10.1095/biolreprod.102.013995 Appl Immunohistochem Mol Morphol Oktay 12 2004 Immunohistochemical analysis of tyrosine phosphorylation and AP-1 transcription factors c-Jun, Jun D, and Fos family during early ovarian follicle development in the mouse 10.1097/00129039-200412000-00014 Biol Reprod Padua 91 24 2014 Simultaneous gene deletion of gata4 and gata6 leads to early disruption of follicular development and germ cell loss in the murine ovary 10.1095/biolreprod.113.117002 Dev Biol Efimenko 381 2013 The transcription factor GATA4 is required for follicular development and normal ovarian function 10.1016/j.ydbio.2013.06.004 Am J Reprod Immunol Mo 82 2019 ANXA7 regulates trophoblast proliferation and apoptosis in preeclampsia 10.1111/aji.13183 Biol Open Huang 7 2018 JAK signaling regulates germline cyst breakdown and primordial follicle formation in mice 10.1242/bio.029470 Biochem Biophys Res Commun Chen 453 2014 The JNK1/JNK3 interactome–contributions by the JNK3 unique N-terminus and JNK common docking site residues 10.1016/j.bbrc.2014.09.122 Exp Cell Res Stepp 316 2010 Syndecan-1 regulates cell migration and fibronectin fibril assembly 10.1016/j.yexcr.2010.05.020 Hum Cell Yu 33 2020 SDC1 promotes cisplatin resistance in hepatic carcinoma cells via PI3K-AKT pathway 10.1007/s13577-020-00362-6 Reprod Biol Endocrinol Gong 18 121 2020 Growth hormone activates PI3K/Akt signaling and inhibits ROS accumulation and apoptosis in granulosa cells of patients with polycystic ovary syndrome 10.1186/s12958-020-00677-x BMC Bioinf Venselaar 11 548 2010 Protein structure analysis of mutations causing inheritable diseases. An e-Science approach with life scientist friendly interfaces 10.1186/1471-2105-11-548 Int J Clin Exp Pathol Huang 8 2015 Role of midkine-progranulin interaction during angiogenesis of hepatocellular carcinoma Hum Reprod Update Qin 21 787 2015 Genetics of primary ovarian insufficiency: new developments and opportunities 10.1093/humupd/dmv036 DataSheet_1.zip 10.3389/fendo.2023.1285667.s001 https://www.frontiersin.org/articles/10.3389/fendo.2023.1285667/supplementary-material/10.3389/fendo.2023.1285667.s001