Cell-type specific analysis of physiological action of estrogen in mouse oviducts

Progesterone signaling in oviductal epithelial cells modulates the immune response to support preimplantation embryonic development

More than 60% of pregnancy losses occur during the first trimester, highlighting the need to understand the role of the oviduct in early pregnancy. In this study, we conditionally ablated the classical progesterone receptor (Pgr) in oviductal epithelial cells, called the Pgrd/d mouse model. We found that 40% of embryos collected from Pgrd/d females were nonviable or developmentally delayed, indicating that epithelial PGR expression is crucial for embryonic development. Single-cell RNA sequencing revealed up-regulation of proinflammatory genes, including interleukin-22 (IL-22), in the epithelial cells of Pgrd/d females. Pharmacological inhibition of inflammation using nonsteroidal anti-inflammatory drugs significantly reduced IL-22 levels in the oviducts and rescued embryonic developmental rates in Pgrd/d females. Coculture of wild-type zygotes with IL-22 significantly decreased the number of expanded blastocysts. Our findings suggest that progesterone signaling is vital for immunoregulation and normal preimplantation development, potentially providing insights for developing diagnostic tools and therapeutic strategies to address pregnancy failures.

Multi-omics analyses and machine learning prediction of oviductal responses in the presence of gametes and embryos

The oviduct is the site of fertilization and preimplantation embryo development in mammals. Evidence suggests that gametes alter oviductal gene expression. To delineate the adaptive interactions between the oviduct and gamete/embryo, we performed a multi-omics characterization of oviductal tissues utilizing bulk RNA-sequencing (RNA-seq), single-cell RNA-sequencing (scRNA-seq), and proteomics collected from distal and proximal at various stages after mating in mice. We observed robust region-specific transcriptional signatures. Specifically, the presence of sperm induces genes involved in pro-inflammatory responses in the proximal region at 0.5 days post-coitus (dpc). Genes involved in inflammatory responses were produced specifically by secretory epithelial cells in the oviduct. At 1.5 and 2.5 dpc, genes involved in pyruvate and glycolysis were enriched in the proximal region, potentially providing metabolic support for developing embryos. Abundant proteins in the oviductal fluid were differentially observed between naturally fertilized and superovulated samples. RNA-seq data were used to identify transcription factors predicted to influence protein abundance in the proteomic data via a novel machine learning model based on transformers of integrating transcriptomics and proteomics data. The transformers identified influential transcription factors and correlated predictive protein expressions in alignment with the in vivo-derived data. Lastly, we found some differences between inflammatory responses in sperm-exposed mouse oviducts compared to hydrosalpinx Fallopian tubes from patients. In conclusion, our multi-omics characterization and subsequent in vivo confirmation of proteins/RNAs indicate that the oviduct is adaptive and responsive to the presence of sperm and embryos in a spatiotemporal manner.

Decoding Müllerian Duct Epithelial Regionalization

Müllerian ducts (MD), also known as paramesonephric ducts, are the primordial anlage of the female reproductive tract organs including the oviduct, uterus, cervix and upper vagina along the craniocaudal axis. Although the general architecture of MD-derived organs is conserved, each organ possesses their unique epithelial structures and cell types to confer their region-specific functions, which collectively coordinate successful fertilization and pregnancy. MD epithelial fate decisions and differentiation along the craniocaudal axis is dependent on spatiotemporal regulation of intrinsic transcription factors and extrinsic signals derived from the mesenchyme. Findings from genetic mouse models, single-cell sequencing studies, and organoid cultures have significantly advanced our understanding of the cellular and molecular mechanisms of MD regionalization. In this review, we first discuss the diversity of epithelial morphologies and cell types in the female reproductive tract organs. Then, we discuss the roles of key transcription factors (Hox, transcriptional cascade driving multiciliogenesis, Foxa2, and P63), signaling pathways (estrogen/ESR1, Wnt/β-catenin, hedgehog, and retinoic acid), and epigenetic factors (microRNAs, chromatin remodeling factors, and histone modification enzymes) in region-specific MD differentiation. Further deciphering molecular mechanisms of MD craniocaudal patterning will open new avenues to improve our strategies for prevention, diagnosis, and treatment of Müllerian anomalies and female reproductive tract disorders.

FOXA1 expression and its association with mucin expression and KRAS mutation in ovarian mucinous tumors: implications for tumor progression and differentiation

This study aimed to investigate forkhead box A1 (FOXA1) expression in ovarian mucinous tumors and its association with mucin expression and KRAS mutation status to clarify its role in tumor progression and differentiation. We analyzed 57 normal tissues or benign ovarian lesions, 110 mucinous ovarian tumors, including mucinous carcinomas, and 214 other ovarian epithelial carcinomas using immunohistochemistry for FOXA1, MUC2, MUC5AC, and MUC6. We also performed KRAS mutation analysis. Strong nuclear staining of FOXA1 was observed in Walthard nests, Brenner tumors, and fallopian tube ciliated epithelium. FOXA1 expression was significantly associated with mucinous histology in ovarian epithelial carcinomas (P < 0.001). In mucinous tumors, FOXA1 was expressed in 73.6% of cystadenomas, 91.4% of borderline tumors, 100% of borderline tumors with intraepithelial carcinomas, and 87.5% of carcinomas. MUC2 expression progressively increased from mucinous cystadenomas to borderline tumors (P < 0.050) and significantly correlated with FOXA1 expression (P = 0.024). The prevalence of KRAS mutations also tended to increase with the malignancy of mucinous tumors (P < 0.050); however, KRAS mutations were significantly enriched in FOXA1-negative cystadenomas compared with FOXA1-positive cystadenomas (P < 0.050). A stepwise increase was noted in the percentage of both KRAS mutations and FOXA1 expression from cystadenoma to carcinoma. Mucinous ovarian tumors commonly express FOXA1. The co-occurrence of KRAS mutations and FOXA1 expression may be important for driving the progression and intestinal differentiation of mucinous ovarian tumors.

Multi-omics analyses and machine learning prediction of oviductal responses in the presence of gametes and embryos

The oviduct is the site of fertilization and preimplantation embryo development in mammals. Evidence suggests that gametes alter oviductal gene expression. To delineate the adaptive interactions between the oviduct and gamete/embryo, we performed a multi-omics characterization of oviductal tissues utilizing bulk RNA-sequencing (RNA-seq), single-cell RNA-sequencing (scRNA-seq), and proteomics collected from distal and proximal at various stages after mating in mice. We observed robust region-specific transcriptional signatures. Specifically, the presence of sperm induces genes involved in pro-inflammatory responses in the proximal region at 0.5 days post-coitus (dpc). Genes involved in inflammatory responses were produced specifically by secretory epithelial cells in the oviduct. At 1.5 and 2.5 dpc, genes involved in pyruvate and glycolysis were enriched in the proximal region, potentially providing metabolic support for developing embryos. Abundant proteins in the oviductal fluid were differentially observed between naturally fertilized and superovulated samples. RNA-seq data were used to identify transcription factors predicted to influence protein abundance in the proteomic data via a novel machine learning model based on transformers of integrating transcriptomics and proteomics data. The transformers identified influential transcription factors and correlated predictive protein expressions in alignment with the in vivo-derived data. Lastly, we found some differences between inflammatory responses in sperm-exposed mouse oviducts compared to hydrosalpinx fallopian tubes from patients. In conclusion, our multi-omics characterization and subsequent in vivo confirmation of proteins/RNAs indicate that the oviduct is adaptive and responsive to the presence of sperm and embryos in a spatiotemporal manner.

Morphometric Evaluation of the Ampullae of the Oviducts in Prepubertal Gilts Treated With Chorionic Gonadotropins

This study assessed morphometric traits of the ampulla of the oviducts in prepubertal gilts treated with chorionic gonadotropins. With the day of slaughter as D0, gilts were assigned to four treatments (n = 8 each): control (untreated), eCG (200 IU eCG on D3), eCG+hCG (1200 IU eCG on D6 plus 500 IU hCG on D3), and eCG+hCG+AI (the previous treatment plus artificial insemination on D1). Blood and ampullae samples were collected at slaughter. Serum progesterone concentrations were higher for gilts treated with hCG than for those in the eCG and control treatments (p < 0.001), but estradiol concentrations did not differ (p > 0.05). The epithelium, muscle and lumen areas and the inner and larger ampullae diameters did not differ across treatments (p > 0.05). Therefore, treatment with chorionic gonadotropins did not alter the ampullae morphometry of prepubertal gilts.

Single-cell transcriptomic profiling of the neonatal oviduct and uterus reveals new insights into upper Müllerian duct regionalization

The upper Müllerian duct (MD) is patterned and specified into two morphologically and functionally distinct organs, the oviduct and uterus. It is known that this regionalization process is instructed by inductive signals from the adjacent mesenchyme. However, the interaction landscape between epithelium and mesenchyme during upper MD development remains largely unknown. Here, we performed single-cell transcriptomic profiling of mouse neonatal oviducts and uteri at the initiation of MD epithelial differentiation (postnatal day 3). We identified major cell types including epithelium, mesenchyme, pericytes, mesothelium, endothelium, and immune cells in both organs with established markers. Moreover, we uncovered region-specific epithelial and mesenchymal subpopulations and then deduced region-specific ligand-receptor pairs mediating mesenchymal-epithelial interactions along the craniocaudal axis. Unexpectedly, we discovered a mesenchymal subpopulation marked by neurofilaments with specific localizations at the mesometrial pole of both the neonatal oviduct and uterus. Lastly, we analyzed and revealed organ-specific signature genes of pericytes and mesothelial cells. Taken together, our study enriches our knowledge of upper MD development, and provides a manageable list of potential genes, pathways, and region-specific cell subtypes for future functional studies.

Metabolic Profiling in Plasma and Brain Induced by 17β-Estradiol Supplementation in Ovariectomized Mice

17β-Estradiol is an ovarian hormone that regulates energy circulation and storage by acting on the central nervous system. However, the metabolic differences between the blood and brain when stimulated by 17β-estradiol are poorly understood. Moreover, research using menopause-induced models to investigate primary metabolites in the blood and brain is limited. Thus, this study aimed to identify metabolic changes in the plasma and brain resulting from 17β-estradiol supplementation in an estrogen-deficient mouse model. Three groups of mice were utilized: sham-operated mice (Sham), ovariectomized mice (OVX), and ovariectomized mice that received a weekly supplementation of 17β-estradiol (E2). Plasma and brain samples from these mice were subjected to metabolic analysis using gas chromatography-time-of-flight-mass spectrometry. Compared with the plasma samples from the Sham and OVX groups, the plasma samples from the E2 group contained higher contents of branched-chain amino acids (BCAAs), such as valine, isoleucine, and leucine. Meanwhile, the brain samples from the E2 group contained higher contents of most metabolites, including BCAAs, neurotransmitters, tricarboxylic acid cycle intermediates, and fatty acids, than those from the two other groups. This study is the first to reveal differences in energy metabolism induced by 17β-estradiol supplementation through brain metabolic profiling of ovariectomized mice, emphasizing the importance of brain metabolic profiling in menopausal hormone research.

Sperm in the Mammalian Female Reproductive Tract: Surfing Through the Tract to Try to Beat the Odds

Mammalian sperm are deposited in the vagina or the cervix/uterus at coitus or at artificial insemination, and the fertilizing sperm move through the female reproductive tract to the ampulla of the oviduct, the site of fertilization. But the destination of most sperm is not the oviduct. Most sperm are carried by retrograde fluid flow to the vagina, are phagocytosed, and/or do not pass barriers on the pathway to the oviduct. The sperm that reach the site of fertilization are the exceptions and winners of one of the most stringent selection processes in nature. This review discusses the challenges sperm encounter and how the few sperm that reach the site of fertilization overcome them. The sperm that reach the goal must navigate viscoelastic fluid, swim vigorously and cooperatively along the walls of the female tract, avoid the innate immune system, and respond to potential cues to direct their movement.

Convergent differentiation of multiciliated cells

Multiciliated cells (MCCs) are epithelial cells that control body fluid flow and contribute to the clearance of pathogenic microbes and other particles from the airways, egg transport in oviducts, and circulation of cerebrospinal fluid in the central nervous system. Although MCCs have shared functions to control fluid flow via coordinated motility of multiple ciliary structures, they are found in multiple mammalian tissues originating from distinct germ layers and differentiate via distinct developmental pathways. To understand the similarities and differences of MCCs in multiple tissues, we investigated single-cell transcriptome data of nasal epithelial cells, bronchial tubes, fallopian tubes, and ependymal cells in the subventricular zone from humans and mice by cross-species data integration. Expression of cilia-associated genes was indistinguishable between these MCCs, although cell populations had unique properties by the species and tissue, demonstrating that they share the same final differentiation status for ciliary functions. We further analyzed the final differentiation step of MCCs from their distinctive progenitors and confirmed their convergent gene set expression for ciliogenesis at the final step. These results may provide new insight into understanding ciliogenesis during the developmental process.

Single-cell transcriptomic profiling of the neonatal oviduct and uterus reveals new insights into upper Müllerian duct regionalization

The upper Müllerian duct (MD) is patterned and specified into two morphologically and functionally distinct organs, the oviduct and uterus. It is known that this regionalization process is instructed by inductive signals from the adjacent mesenchyme. However, the interaction landscape between epithelium and mesenchyme during upper MD development remains largely unknown. Here, we performed single-cell transcriptomic profiling of mouse neonatal oviducts and uteri at the initiation of MD epithelial differentiation (postnatal day 3). We identified major cell types including epithelium, mesenchyme, pericytes, mesothelium, endothelium, and immune cells in both organs with established markers. Moreover, we uncovered region-specific epithelial and mesenchymal subpopulations and then deduced region-specific ligand-receptor pairs mediating mesenchymal-epithelial interactions along the craniocaudal axis. Unexpectedly, we discovered a mesenchymal subpopulation marked by neurofilaments with specific localizations at the mesometrial pole of both the neonatal oviduct and uterus. Lastly, we analyzed and revealed organ-specific signature genes of pericytes and mesothelial cells. Taken together, our study enriches our knowledge of upper Müllerian duct development, and provides a manageable list of potential genes, pathways, and region-specific cell subtypes for future functional studies.

Classical Estrogen Signaling in Ciliated Epithelial Cells of the Oviduct Is Nonessential for Fertility in Female Mice

Ciliary action performs a critical role in the oviduct (Fallopian tube) during pregnancy establishment through sperm and egg transport. The disruption of normal ciliary function in the oviduct affects oocyte pick-up and is a contributing factor to female infertility. Estrogen is an important regulator of ciliary action in the oviduct and promotes ciliogenesis in several species. Global loss of estrogen receptor α (ESR1) leads to infertility. We have previously shown that ESR1 in the oviductal epithelial cell layer is required for female fertility. Here, we assessed the role of estrogen on transcriptional regulation of ciliated epithelial cells of the oviduct using single-cell RNA-sequencing analysis. We observed minor variations in ciliated cell genes in the proximal region (isthmus and uterotubal junction) of the oviduct. However, 17β-estradiol treatment had little impact on the gene expression profile of ciliated epithelial cells. We also conditionally ablated Esr1 from ciliated epithelial cells of the oviduct (called ciliated Esr1d/d mice). Our studies showed that ciliated Esr1d/d females had fertility rates comparable to control females, did not display any disruptions in preimplantation embryo development or embryo transport to the uterus, and had comparable cilia formation to control females. However, we observed some incomplete deletion of Esr1 in the ciliated epithelial cells, especially in the ampulla region. Nevertheless, our data suggest that ESR1 expression in ciliated cells of the oviduct is dispensable for ciliogenesis and nonessential for female fertility in mice.

Chromatin remodeling of prostaglandin signaling in smooth muscle enables mouse embryo passage through the female reproductive tract

Early mammalian development occurs during embryo transit of the female reproductive tract. Transport is orchestrated by secreted oviduct fluid, unidirectional beating of epithelial cilia, and smooth muscle contractions. Using gene-edited mice, we document that conditional disruption of a component of the SWI/SNF chromatin remodeling complex in smooth muscle cells prevents transport through the oviduct without perturbing embryogenesis. Analysis with RNA sequencing (RNA-seq), transposase-accessible chromatin with sequencing (ATAC-seq), chromatin immunocleavage sequencing (ChIC-seq), and pharmacologic rescue experiments implicated prostaglandin signaling pathways. In comparison with controls, gene-edited mice had compromised chromatin accessibility at enhancer/promoters of Ptgs2, Pla2g16, Pla2r1, and Ptger3 (EP3) as well as decreased enhancer-promoter interactive looping critical for Ptgs2 (aka Cox-2) expression in a SWI/SNF complex-dependent manner. Treatment of wild-type mice with prostaglandin inhibitors phenocopied the genetically induced defect.

Morphogenesis of the female reproductive tract along antero-posterior and dorso-ventral axes is dependent on Amhr2+ mesenchyme in mice†

Morphogenesis of the female reproductive tract is regulated by the mesenchyme. However, the identity of the mesenchymal lineage that directs the morphogenesis of the female reproductive tract has not been determined. Using in vivo genetic cell ablation, we identified Amhr2+ mesenchyme as an essential mesenchymal population in patterning the female reproductive tract. After partial ablation of Amhr2+ mesenchymal cells, the oviduct failed to develop its characteristic coiling due to decreased epithelial proliferation and tubule elongation during development. The uterus displayed a reduction in size and showed decreased cellular proliferation in both epithelial and mesenchymal compartments. More importantly, in the uterus, partial ablation of Amhr2+ mesenchyme caused abnormal lumen shape and altered the direction of its long axis from the dorsal-ventral axis to the left-right axis (i.e., perpendicular to the dorsal-ventral axis). Despite these morphological defects, epithelia underwent normal differentiation into secretory and ciliated cells in the oviduct and glandular epithelial cells in the uterus. These results demonstrated that Amhr2+ mesenchyme can direct female reproductive tract morphogenesis by regulating epithelial proliferation and lumen shape without affecting the differentiation of epithelial cell types.

Spatial transcriptomic profiles of mouse uterine microenvironments at pregnancy day 7.5†

Uterine dysfunctions lead to fertility disorders and pregnancy complications. Normal uterine functions at pregnancy depend on crosstalk among multiple cell types in uterine microenvironments. Here, we performed the spatial transcriptomics and single-cell RNA-seq assays to determine local gene expression profiles at the embryo implantation site of the mouse uterus on pregnancy day 7.5 (D7.5). The spatial transcriptomic annotation identified 11 domains of distinct gene signatures, including a mesometrial myometrium, an anti-mesometrial myometrium, a mesometrial decidua enriched with natural killer cells, a vascular sinus zone for maternal vessel remodeling, a fetal-maternal interface, a primary decidual zone, a transition decidual zone, a secondary decidual zone, undifferentiated stroma, uterine glands, and the embryo. The scRNA-Seq identified 12 types of cells in the D7.5 uterus including three types of stromal fibroblasts with differentiated and undifferentiated markers, one cluster of epithelium including luminal and glandular epithelium, mesothelium, endothelia, pericytes, myelomonocytic cell, natural killer cells, and lymphocyte B. These single-cell RNA signatures were then utilized to deconvolute the cell-type compositions of each individual uterine microenvironment. Functional annotation assays on spatial transcriptomic data revealed uterine microenvironments with distinguished metabolic preferences, immune responses, and various cellular behaviors that are regulated by region-specific endocrine and paracrine signals. Global interactome among regions is also projected based on the spatial transcriptomic data. This study provides high-resolution transcriptome profiles with locality information at the embryo implantation site to facilitate further investigations on molecular mechanisms for normal pregnancy progression.

Progesterone and Inflammatory Response in the Oviduct during Physiological and Pathological Conditions

Progesterone has been shown to be a potent suppressor of several inflammatory pathways. During pregnancy, progesterone levels increase, allowing for normal pregnancy establishment and maintenance. The dysregulation of progesterone, as well as inflammation, leads to poor pregnancy outcomes. However, it is unclear how progesterone imbalance could impact inflammatory responses in the oviduct and subsequently result in early pregnancy loss. Therefore, in this review, we describe the role of progesterone signaling in regulating the inflammatory response, with a focus on the oviduct and pathological conditions in the Fallopian tubes.

Spatiotemporal transcriptional dynamics of the cycling mouse oviduct

Female fertility in mammals requires iterative remodeling of the entire adult female reproductive tract across the menstrual/estrous cycle. However, while transcriptome dynamics across the estrous cycle have been reported in human and bovine models, no global analysis of gene expression across the estrous cycle has yet been reported for the mouse. Here, we examined the cellular composition and global transcriptional dynamics of the mouse oviduct along the anteroposterior axis and across the estrous cycle. We observed robust patterns of differential gene expression along the anteroposterior axis, but we found surprisingly few changes in gene expression across the estrous cycle. Notable gene expression differences along the anteroposterior axis included a surprising enrichment for genes related to embryonic development, such as Hox and Wnt genes. The relatively stable transcriptional dynamics across the estrous cycle differ markedly from other mammals, leading us to speculate that this is an evolutionarily derived state that may reflect the extremely rapid five-day mouse estrous cycle. This dataset fills a critical gap by providing an important genomic resource for a highly tractable genetic model of mammalian female reproduction.