Single cell RNA-seq study of wild type and Hox9,10,11 mutant developing uterus

Secretion of WNT7A by UC-MSCs assist in promoting the endometrial epithelial regeneration

Stem cell therapy for intrauterine adhesions (IUAs) has been widely used in clinical treatment. However, intravenous injection lacks sufficient targeting capabilities, while in situ injection poses challenges in ensuring the effective survival of stem cells. Furthermore, the mechanism underlying the interaction between stem cells and endometrial cells in vivo remains poorly understood, and there is a lack of suitable in vitro models for studying these problems. Here, we designed an extracellular matrix (ECM)-adhesion mimic hydrogel for intrauterine administration, which was more effective than direct injection in treating IUAs. Additionally, we analyzed the epithelial-mesenchymal transition (EMT) and confirmed that the activation of endometrial epithelial stem cells is pivotal. Our findings demonstrated that umbilical cord mesenchymal stem cells (UC-MSCs) secrete WNT7A to activate endometrial epithelial stem cells, thereby accelerating regeneration of the endometrial epithelium. Concurrently, under transforming growth factor alpha (TGFA) stimulation secreted by the EMT epithelium, UC-MSCs upregulate E-cadherin while partially implanting into the endometrial epithelium.

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.

Hox genes and patterning the vertebrate body

The diversity of vertebrate body plans is dizzying, yet stunning for the many things they have in common. Vertebrates have inhabited virtually every part of the earth from its coldest to warmest climates. They locomote by swimming, flying, walking, slithering, or climbing, or combinations of these behaviors. And they exist in many different sizes, from the smallest of frogs, fish and lizards to giraffes, elephants, and blue whales. Despite these differences, vertebrates follow a remarkably similar blueprint for the establishment of their body plan. Within the relatively small amount of time required to complete gastrulation, the process through which the three germ layers, ectoderm, mesoderm, and endoderm are created, the embryo also generates its body axis and is simultaneously patterned. For the length of this axis, the genes that distinguish the neck from the rib cage or the trunk from the sacrum are the Hox genes. In vertebrates, there was evolutionary pressure to maintain this set of genes in the organism. Over the past decades, much has been learned regarding the regulatory mechanisms that ensure the appropriate expression of these genes along the main body axes. Genetic functions continue to be explored though much has been learned. Much less has been discerned on the identity of co-factors used by Hox proteins for the specificity of transcriptional regulation or what downstream targets and pathways are critical for patterning events, though there are notable exceptions. Current work in the field is demonstrating that Hox genes continue to function in many organs long after directing early patterning events. It is hopeful continued research will shed light on remaining questions regarding mechanisms used by this important and conserved set of transcriptional regulators.

Mapping uterine calcium dynamics during the ovulatory cycle in live mice

Uterine contraction patterns vary during the ovulatory cycle and throughout pregnancy but prior measurements have produced limited and conflicting information on these patterns. We combined a virally delivered genetically encoded calcium reporter (GCaMP8m) and ultra-widefield imaging in live nonpregnant mice to characterize uterine calcium dynamics at organ scale throughout the estrous cycle. Prior to ovulation (proestrus and estrus) uterine excitations primarily initiated in a region near the oviduct, but after ovulation (metestrus and diestrus), excitations initiated at loci homogeneously distributed throughout the organ. The frequency of excitation events was lowest in proestrus and estrus, higher in metestrus and highest in diestrus. These results establish a platform for mapping uterine activity, and show that the question of whether there is an anatomically localized trigger for uterine excitations depends on the estrous cycle phase.

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.

Single-cell insights into epithelial morphogenesis in the neonatal mouse uterus

The uterus is vital for successful reproduction in mammals, and two different types of epithelia (luminal and glandular) are essential for embryo implantation and pregnancy establishment. However, the essential cellular and molecular factors and pathways governing postnatal epithelium maturation, determination, and differentiation in developing uterus are yet to be elucidated. Here, the epithelium of the neonatal mouse uterus was isolated and subjected to single-cell transcriptome (scRNA-seq) analysis. Both the undifferentiated epithelium and determined luminal epithelium were heterogeneous and contained several different cell clusters based on single-cell transcription profiles. Substantial gene expression differences were evident as the epithelium matured and differentiated between postnatal days 1 to 15. Two new glandular epithelium-expressed genes (Gas6 and Cited4) were identified and validated by in situ hybridization. Trajectory analyses provided a framework for understanding epithelium maturation, lineage bifurcation, and differentiation. A candidate set of transcription factors and gene regulatory networks were identified that potentially direct epithelium lineage specification and morphogenesis. This atlas provides a foundation important to discover intrinsic cellular and molecular mechanisms directing uterine epithelium morphogenesis during a critical window of postnatal development.

Sctensor detects many-to-many cell-cell interactions from single cell RNA-sequencing data

BACKGROUND

Complex biological systems are described as a multitude of cell-cell interactions (CCIs). Recent single-cell RNA-sequencing studies focus on CCIs based on ligand-receptor (L-R) gene co-expression but the analytical methods are not appropriate to detect many-to-many CCIs.

RESULTS

In this work, we propose scTensor, a novel method for extracting representative triadic relationships (or hypergraphs), which include ligand-expression, receptor-expression, and related L-R pairs.

CONCLUSIONS

Through extensive studies with simulated and empirical datasets, we have shown that scTensor can detect some hypergraphs that cannot be detected using conventional CCI detection methods, especially when they include many-to-many relationships. scTensor is implemented as a freely available R/Bioconductor package.

Developmental estrogen exposure in mice disrupts uterine epithelial cell differentiation and causes adenocarcinoma via Wnt/β-catenin and PI3K/AKT signaling

Tissue development entails genetically programmed differentiation of immature cell types to mature, fully differentiated cells. Exposure during development to non-mutagenic environmental factors can contribute to cancer risk, but the underlying mechanisms are not understood. We used a mouse model of endometrial adenocarcinoma that results from brief developmental exposure to an estrogenic chemical, diethylstilbestrol (DES), to determine causative factors. Single-cell RNA sequencing (scRNAseq) and spatial transcriptomics of adult control uteri revealed novel markers of uterine epithelial stem cells (EpSCs), identified distinct luminal and glandular progenitor cell (PC) populations, and defined glandular and luminal epithelium (LE) cell differentiation trajectories. Neonatal DES exposure disrupted uterine epithelial cell differentiation, resulting in a failure to generate an EpSC population or distinguishable glandular and luminal progenitors or mature cells. Instead, the DES-exposed epithelial cells were characterized by a single proliferating PC population and widespread activation of Wnt/β-catenin signaling. The underlying endometrial stromal cells had dramatic increases in inflammatory signaling pathways and oxidative stress. Together, these changes activated phosphoinositide 3-kinase/AKT serine-threonine kinase signaling and malignant transformation of cells that were marked by phospho-AKT and the cancer-associated protein olfactomedin 4. Here, we defined a mechanistic pathway from developmental exposure to an endocrine disrupting chemical to the development of adult-onset cancer. These findings provide an explanation for how human cancers, which are often associated with abnormal activation of PI3K/AKT signaling, could result from exposure to environmental insults during development.

Hox genes in development and beyond

Hox genes encode evolutionarily conserved transcription factors that are essential for the proper development of bilaterian organisms. Hox genes are unique because they are spatially and temporally regulated during development in a manner that is dictated by their tightly linked genomic organization. Although their genetic function during embryonic development has been interrogated, less is known about how these transcription factors regulate downstream genes to direct morphogenetic events. Moreover, the continued expression and function of Hox genes at postnatal and adult stages highlights crucial roles for these genes throughout the life of an organism. Here, we provide an overview of Hox genes, highlighting their evolutionary history, their unique genomic organization and how this impacts the regulation of their expression, what is known about their protein structure, and their deployment in development and beyond.

Biology-inspired data-driven quality control for scientific discovery in single-cell transcriptomics

BACKGROUND

Quality control (QC) of cells, a critical first step in single-cell RNA sequencing data analysis, has largely relied on arbitrarily fixed data-agnostic thresholds applied to QC metrics such as gene complexity and fraction of reads mapping to mitochondrial genes. The few existing data-driven approaches perform QC at the level of samples or studies without accounting for biological variation.

RESULTS

We first demonstrate that QC metrics vary with both tissue and cell types across technologies, study conditions, and species. We then propose data-driven QC (ddqc), an unsupervised adaptive QC framework to perform flexible and data-driven QC at the level of cell types while retaining critical biological insights and improved power for downstream analysis. ddqc applies an adaptive threshold based on the median absolute deviation on four QC metrics (gene and UMI complexity, fraction of reads mapping to mitochondrial and ribosomal genes). ddqc retains over a third more cells when compared to conventional data-agnostic QC filters. Finally, we show that ddqc recovers biologically meaningful trends in gradation of gene complexity among cell types that can help answer questions of biological interest such as which cell types express the least and most number of transcripts overall, and ribosomal transcripts specifically.

CONCLUSIONS

ddqc retains cell types such as metabolically active parenchymal cells and specialized cells such as neutrophils which are often lost by conventional QC. Taken together, our work proposes a revised paradigm to quality filtering best practices-iterative QC, providing a data-driven QC framework compatible with observed biological diversity.

Single-Cell RNA-Sequencing Reveals Interactions between Endometrial Stromal Cells, Epithelial Cells, and Lymphocytes during Mouse Embryo Implantation

The decidualization of endometrial stromal cells (ESCs) is an essential process facilitating embryo implantation. However, the roles of non-decidualized and decidualized ESCs in regulating the microenvironment of a receptive endometrium remain unclear. We investigated single-cell transcriptomic changes in the uterus of a CD-1 mouse model at the post-implantation stage. The implantation and inter-implantation sites of the uteruses of pregnant mice at 4.5 and 5.5 days post-coitum were dissected for single-cell RNA sequencing. We identified eight cell types: epithelial cells, stromal cells, endothelial cells, mesothelial cells, lymphocytes, myocytes, myeloids, and pericytes. The ESC transcriptome suggests that the four ESC subtypes are involved in the extracellular remodeling during implantation. The trajectory plot of ESC subtypes indicates embryo implantation that involves a differentiation pathway from undifferentiated ESCs (ESC 1) to decidualized ESCs (DEC ESCs), with distinct signaling pathways between the ESC subtypes. Furthermore, the ligand-receptor analysis suggests that ESCs communicate with epithelial cells and immune cells through nectin and ICAM signaling. Collectively, both decidualized and non-decidualized ESCs may regulate the endometrial microenvironment for optimal endometrial receptivity and immune tolerance. This study provides insights on the molecular and cellular characteristics of mouse ESCs in modulating the epithelial and lymphocyte functions during early embryo implantation.

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.

The Estrogen Receptor α Cistrome in Human Endometrium and Epithelial Organoids

Endometrial health is affected by molecular processes that underlie estrogen responses. We assessed estrogen regulation of endometrial function by integrating the estrogen receptor α (ESR1) cistromes and transcriptomes of endometrial biopsies taken from the proliferative and mid-secretory phases of the menstrual cycle together with hormonally stimulated endometrial epithelial organoids. The cycle stage-specific ESR1 binding sites were determined by chromatin immunoprecipitation and next-generation sequencing and then integrated with changes in gene expression from RNA sequencing data to infer candidate ESR1 targets in normal endometrium. Genes with ESR1 binding in whole endometrium were enriched for chromatin modification and regulation of cell proliferation. The distribution of ESR1 binding sites in organoids was more distal from gene promoters when compared to primary endometrium and was more similar to the proliferative than the mid-secretory phase ESR1 cistrome. Inferred organoid estrogen/ESR1 candidate target genes affected formation of cellular protrusions and chromatin modification. Comparison of signaling effected by candidate ESR1 target genes in endometrium vs organoids reveals enrichment of both overlapping and distinct responses. Our analysis of the ESR1 cistromes and transcriptomes from endometrium and organoids provides important resources for understanding how estrogen affects endometrial health and function.

Contribution of the Wolffian duct mesenchyme to the formation of the female reproductive tract

The female reproductive tract develops from its embryonic precursor, the Müllerian duct. In close proximity to the Müllerian duct lies the precursor for the male reproductive tract, the Wolffian duct, which is eliminated in the female embryo during sexual differentiation. We discovered that a component of the Wolffian duct, its mesenchyme, is not eliminated after sexual differentiation. Instead, the Wolffian duct mesenchyme underwent changes in transcriptome and chromatin accessibility from male tract to female tract identity, and became a unique mesenchymal population in the female reproductive tract with localization and transcriptome distinct from the mesenchyme derived from the Müllerian duct. Partial ablation of the Wolffian duct mesenchyme stunted the growth of the fetal female reproductive tract in ex vivo organ culture. These findings reveal a new fetal origin of mesenchymal tissues for female reproductive tract formation and reshape our understanding of sexual differentiation of reproductive tracts.

Revisiting the Transcriptome Landscape of Pig Embryo Implantation Site at Single-Cell Resolution

Litter size is one of the most economically important traits in commercial pig farming. It has been estimated that approximately 30% of porcine embryos are lost during the peri-implantation period. Despite rapid advances over recent years, the molecular mechanism underlying embryo implantation in pigs remains poorly understood. In this study, the conceptus together with a small amount of its surrounding endometrial tissues at the implantation site was collected and subjected to single-cell RNA-seq using the 10x platform. Because embryo and maternal endometrium were genetically different, we successfully dissected embryonic cells from maternal endometrial cells in the data according to single nucleotide polymorphism information captured by single-cell RNA-seq. Undoubtedly, the interaction between trophoblast cells and uterine epithelial cells represents the key mechanism of embryo implantation. Using the CellChat tool, we revealed cell-cell communications between these 2 cell types in terms of secreted signaling, ECM-receptor interaction and cell-cell contact. Additionally, by analyzing the non-pregnant endometrium as control, we were able to identify global gene expression changes associated with embryo implantation in each cell type. Our data provide a valuable resource for deciphering the molecular mechanism of embryo implantation in pigs.

Regulation of posterior Hox genes by sex steroids explains vertebral variation in inbred mouse strains

A series of elegant embryo transfer experiments in the 1950s demonstrated that the uterine environment could alter vertebral patterning in inbred mouse strains. In the intervening decades, attention has tended to focus on the technical achievements involved and neglected the underlying biological question: how can genetically homogenous individuals have a heterogenous number of vertebrae? Here I revisit these experiments and, with the benefit of knowledge of the molecular-level processes of vertebral patterning gained over the intervening decades, suggest a novel hypothesis for homeotic transformation of the last lumbar vertebra to the adjacent sacral type through regulation of Hox genes by sex steroids. Hox genes are involved in both axial patterning and development of male and female reproductive systems and have been shown to be sensitive to sex steroids in vitro and in vivo. Regulation of these genes by sex steroids and resulting alterations to vertebral patterning may hint at a deep evolutionary link between the ribless lumbar region of mammals and the switch from egg-laying to embryo implantation. An appreciation of the impact of sex steroids on Hox genes may explain some puzzling aspects of human disease, and highlights the spine as a neglected target for in utero exposure to endocrine disruptors.

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.

Single-cell analysis of mouse uterus at the invasion phase of embryo implantation

Background

Embryo implantation into the uterus is a crucial step for human reproduction. A hypothesis has been proposed that the molecular circuit invented by trophoblasts for invasive embryo implantation during evolution might be misused by cancer cells to promote malignancy. Unfortunately, our current understanding of the molecular mechanism underlying embryo implantation is far from complete.

Results

Here we used the mouse as an animal model and generated a single-cell transcriptomic atlas of the embryo implantation site of mouse uterus at the invasion phase of embryo implantation on gestational day 6. We revealed 23 distinct cell clusters, including 5 stromal cell clusters, 2 epithelial cell clusters, 1 smooth muscle cell cluster, 2 pericyte clusters, 4 endothelial cell clusters, and 9 immune cell clusters. Through data analysis, we identified differentially expression changes in all uterine cell types upon embryo implantation. By integrated with single-cell RNA-seq data from E5.5 embryos, we predicted cell–cell crosstalk between trophoblasts and uterine cell types.

Conclusions

Our study provides a valuable resource for understanding of the molecular mechanism of embryo implantation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00749-y.

Deciphering mouse uterine receptivity for embryo implantation at single-cell resolution

Objectives

Mice are widely used as an animal model for studying human uterine receptivity for embryo implantation. Although transcriptional changes related to mouse uterine receptivity have been determined by using bulk RNA‐seq, the data are of limited value because the uterus is a complex organ consisting of many cell types. Here, we aimed to decipher mouse uterine receptivity for embryo implantation at single‐cell resolution.

Materials and methods

Single‐cell RNA sequencing was performed for the pre‐receptive and the receptive mouse uterus. Gene expression profiles in luminal epithelium and glandular epithelium were validated by comparing against a published laser capture microdissection (LCM)‐coupled microarray dataset.

Results

We revealed 19 distinct cell clusters, including 3 stromal cell clusters, 2 epithelial cell clusters, 1 smooth muscle cell cluster, 4 endothelial cell clusters and 8 immune cell clusters. We identified global gene expression changes associated with uterine receptivity in each cell type. Additionally, we predicted signalling interactions for distinct cell types to understand the crosstalk between the blastocyst and the receptive uterus.

Conclusion

Our data provide a valuable resource for deciphering the molecular mechanism underlying uterine receptivity in mice.

Very small embryonic-like stem cells (VSELs) regenerate whereas mesenchymal stromal cells (MSCs) rejuvenate diseased reproductive tissues

Compared to embryonic and induced pluripotent stem cells, mesenchymal stem/stromal cells (MSCs) have made their presence felt with good therapeutic promise and safety profile. Transplanting MSCs has successfully helped to reverse infertility and resulted in live births in animal models and also in humans. But the underlying mechanism for their therapeutic potential is not yet clear. MSCs are not pluripotent and hence lack plasticity to differentiate into multiple adult cell types. They rather act as 'paracrine providers' to the tissue-resident stem cells since similar beneficial effects are also observed when their secretome (microvesicles or exosomes) is transplanted. Cytokines, growth factors, signaling lipids, mRNAs, and miRNAs secreted by MSCs enables tissue-resident stem cells to undergo differentiation into specific cell types. Tissue-resident stem cells include pluripotent, very small embryonic-like stem cells (VSELs) and progenitors [spermatogonial (SSCs), ovarian (OSCs) and endometrial (EnSCs) stem cells in testes, ovary and uterus respectively] which function in a subtle manner to maintain life-long tissue homeostasis and regenerate damaged (non-functional) reproductive tissues by differentiating into sperm, oocytes and endometrial epithelial cells respectively. Similar to restoring spermatogenesis, primordial follicles numbers are increased upon transplanting MSCs. Published literature suggests that MSCs do not differentiate into epithelial cells in the endometrium. Nuclear OCT-4 positive VSELs and cytoplasmic OCT-4, AXIN2 and KERATIN-19 positive epithelial progenitors have a greater role during endometrial regeneration. We propose, transplantation of MSCs simply provides growth factors/cytokines essential for the tissue-resident stem/progenitor cells to undergo differentiation into sperm, eggs and endometrial epithelial cells in the reproductive tissues.

Spatial Transcriptomics analysis of uterine gene expression in enhancer of Zeste homolog 2 (Ezh2) conditional knockout mice

Histone proteins undergo various modifications that alter chromatin structure, including addition of methyl groups. Enhancer of homolog 2 (EZH2), is a histone methyltransferase that methylates lysine residue 27, and thereby, suppresses gene expression. EZH2 plays integral role in the uterus and other reproductive organs. We have previously shown that conditional deletion of uterine EZH2 results in increased proliferation of luminal and glandular epithelial cells, and RNAseq analyses reveal several uterine transcriptomic changes in Ezh2 conditional (c) knockout (KO) mice that can affect estrogen signaling pathways. To pinpoint the origin of such gene expression changes, we used the recently developed spatial transcriptomics (ST) method with the hypotheses that Ezh2cKO mice would predominantly demonstrate changes in epithelial cells and/or ablation of this gene would disrupt normal epithelial/stromal gene expression patterns. Uteri were collected from ovariectomized adult WT and Ezh2cKO mice and analyzed by ST. Asb4, Cxcl14, Dio2, and Igfbp5 were increased, Sult1d1, Mt3, and Lcn2 were reduced in Ezh2cKO uterine epithelium vs. WT epithelium. For Ezh2cKO uterine stroma, differentially expressed key hub genes included Cald1, Fbln1, Myh11, Acta2, and Tagln. Conditional loss of uterine Ezh2 also appears to shift the balance of gene expression profiles in epithelial vs. stromal tissue toward uterine epithelial cell and gland development and proliferation, consistent with uterine gland hyperplasia in these mice. Current findings provide further insight into how EZH2 may selectively affect uterine epithelial and stromal compartments. Additionally, these transcriptome data might provide the mechanistic understanding and valuable biomarkers for human endometrial disorders with epigenetic underpinnings.

Identification of Intercellular Crosstalk between Decidual Cells and Niche Cells in Mice

Decidualization is a crucial step for human reproduction, which is a prerequisite for embryo implantation, placentation and pregnancy maintenance. Despite rapid advances over recent years, the molecular mechanism underlying decidualization remains poorly understood. Here, we used the mouse as an animal model and generated a single-cell transcriptomic atlas of a mouse uterus during decidualization. By analyzing the undecidualized inter-implantation site of the uterus as a control, we were able to identify global gene expression changes associated with decidualization in each cell type. Additionally, we identified intercellular crosstalk between decidual cells and niche cells, including immune cells, endothelial cells and trophoblast cells. Our data provide a valuable resource for deciphering the molecular mechanism underlying decidualization.

Specific Transcriptomic Signatures and Dual Regulation of Steroidogenesis Between Fetal and Adult Mouse Leydig Cells

Leydig cells (LC) are the main testicular androgen-producing cells. In eutherian mammals, two types of LCs emerge successively during testicular development, fetal Leydig cells (FLCs) and adult Leydig cells (ALCs). Both display significant differences in androgen production and regulation. Using bulk RNA sequencing, we compared the transcriptomes of both LC populations to characterize their specific transcriptional and functional features. Despite similar transcriptomic profiles, a quarter of the genes show significant variations in expression between FLCs and ALCs. Non-transcriptional events, such as alternative splicing was also observed, including a high rate of intron retention in FLCs compared to ALCs. The use of single-cell RNA sequencing data also allowed the identification of nine FLC-specific genes and 50 ALC-specific genes. Expression of the corticotropin-releasing hormone 1 (Crhr1) receptor and the ACTH receptor melanocortin type 2 receptor (Mc2r) specifically in FLCs suggests a dual regulation of steroidogenesis. The androstenedione synthesis by FLCs is stimulated by luteinizing hormone (LH), corticotrophin-releasing hormone (CRH), and adrenocorticotropic hormone (ACTH) whereas the testosterone synthesis by ALCs is dependent exclusively on LH. Overall, our study provides a useful database to explore LC development and functions.

Studying Müllerian duct anomalies - from cataloguing phenotypes to discovering causation

Müllerian duct anomalies (MDAs) are developmental disorders of the Müllerian duct, the embryonic anlage of most of the female reproductive tract. The prevalence of MDAs is 6.7% in the general female population and 16.7% in women who exhibit recurrent miscarriages. Individuals affected by these anomalies suffer from high rates of infertility, first-trimester pregnancy losses, premature labour, placental retention, foetal growth retardation and foetal malpresentations. The aetiology of MDAs is complex and heterogeneous, displaying a range of clinical pictures that generally lack a direct genotype-phenotype correlation. De novo and familial cases sharing the same genomic lesions have been reported. The familial cases follow an autosomal-dominant inheritance, with reduced penetrance and variable expressivity. Furthermore, few genetic factors and molecular pathways underpinning Müllerian development and dysregulations causing MDAs have been identified. The current knowledge in this field predominantly derives from loss-of-function experiments in mouse and chicken models, as well as from human genetic association studies using traditional approaches, such as microarrays and Sanger sequencing, limiting the discovery of causal factors to few genetic entities from the coding genome. In this Review, we summarise the current state of the field, discuss limitations in the number of studies and patient samples that have stalled progress, and review how the development of new technologies provides a unique opportunity to overcome these limitations. Furthermore, we discuss how these new technologies can improve functional validation of potential causative alterations in MDAs.

Summary: Here, we review the current knowledge about Müllerian duct anomalies in the context of new high-throughput technologies and model systems and their implications in the prevention of these disorders.

Cell-Cell Communication at the Embryo Implantation Site of Mouse Uterus Revealed by Single-Cell Analysis

As a crucial step for human reproduction, embryo implantation is a low-efficiency process. Despite rapid advances in recent years, the molecular mechanism underlying embryo implantation remains poorly understood. Here, we used the mouse as an animal model and generated a single-cell transcriptomic atlas of embryo implantation sites. By analyzing inter-implantation sites of the uterus as control, we were able to identify global gene expression changes associated with embryo implantation in each cell type. Additionally, we predicted signaling interactions between uterine luminal epithelial cells and mural trophectoderm of blastocysts, which represent the key mechanism of embryo implantation. We also predicted signaling interactions between uterine epithelial-stromal crosstalk at implantation sites, which are crucial for post-implantation development. Our data provide a valuable resource for deciphering the molecular mechanism underlying embryo implantation.

Mechanistic Drivers of Müllerian Duct Development and Differentiation Into the Oviduct

The conduits of life; the animal oviducts and human fallopian tubes are of paramount importance for reproduction in amniotes. They connect the ovary with the uterus and are essential for fertility. They provide the appropriate environment for gamete maintenance, fertilization and preimplantation embryonic development. However, serious pathologies, such as ectopic pregnancy, malignancy and severe infections, occur in the oviducts. They can have drastic effects on fertility, and some are life-threatening. Despite the crucial importance of the oviducts in life, relatively little is known about the molecular drivers underpinning the embryonic development of their precursor structures, the Müllerian ducts, and their successive differentiation and maturation. The Müllerian ducts are simple rudimentary tubes comprised of an epithelial lumen surrounded by a mesenchymal layer. They differentiate into most of the adult female reproductive tract (FRT). The earliest sign of Müllerian duct formation is the thickening of the anterior mesonephric coelomic epithelium to form a placode of two distinct progenitor cells. It is proposed that one subset of progenitor cells undergoes partial epithelial-mesenchymal transition (pEMT), differentiating into immature Müllerian luminal cells, and another subset undergoes complete EMT to become Müllerian mesenchymal cells. These cells invaginate and proliferate forming the Müllerian ducts. Subsequently, pEMT would be reversed to generate differentiated epithelial cells lining the fully formed Müllerian lumen. The anterior Müllerian epithelial cells further specialize into the oviduct epithelial subtypes. This review highlights the key established molecular and genetic determinants of the processes involved in Müllerian duct development and the differentiation of its upper segment into oviducts. Furthermore, an extensive genome-wide survey of mouse knockout lines displaying Müllerian or oviduct phenotypes was undertaken. In addition to widely established genetic determinants of Müllerian duct development, our search has identified surprising associations between loss-of-function of several genes and high-penetrance abnormalities in the Müllerian duct and/or oviducts. Remarkably, these associations have not been investigated in any detail. Finally, we discuss future directions for research on Müllerian duct development and oviducts.

Downregulation of HOXA11 enhances endometrial cancer malignancy and cisplatin resistance via activating PTEN/AKT signaling pathway

PURPOSE

Endometrial cancer is the most common malignant tumor of female genital system worldwide. Homeobox A11 (HOXA11) is an evolutionarily conserved Homeobox gene closely implicated in carcinogenesis. However, the mechanisms of HOXA11 in the progression and cisplatin resistance of endometrial cancer remain unclear.

METHODS

The expression of HOXA11 was analyzed based on 548 endometrial cancer and 35 control tissues from The Cancer Genome Atlas (TCGA) database. Transwell assay was performed to investigate the effect of HOXA11 on endometrial cell migration and invasion. TUNEL staining was carried out to assay the role of HOXA11 in endometrial cell apoptosis. Western blot was employed to detect the protein levels of B cell lymphoma-2 (Bcl-2), Bcl-2 associated X (Bax), cleaved caspase-3, matrix metalloproteinase-2/9 (MMP/9), phosphatase and tensin homolog (PTEN), protein kinase B (AKT) and p-AKT.

RESULTS

TCGA data showed that HOXA11 expression was significantly down-regulated in endometrial cancer tissue samples. The overexpression of HOXA11 promoted the apoptosis, but inhibited the proliferation, migration and invasion of endometrial cancer cells. HOXA11 knockdown with small interfering RNA (siRNA) considerably repressed cell apoptosis, while promoted cell proliferation, migration, and invasion through PTEN/AKT signaling pathway. Interestingly, HOXA11 was lowly expressed in Ishikawa cells treated with cisplatin. In addition, HOXA11 knockdown increased the resistance of endometrial cancer to cisplatin through activating PTEN/AKT signaling pathway.

CONCLUSION

Low HOXA11 expression may promote the proliferation, migration, invasion of endometrial cancer cells, and increase their resistance to cisplatin through activating PTEN/AKT pathway.

Phenotyping of immune and endometrial epithelial cells in endometrial carcinomas revealed by single-cell RNA sequencing

Tumors are complex ecosystems harboring multiple cell types which might play a critical role in tumor progression and treatment response. The endometrial epithelial cell identities and immune microenvironment of endometrial carcinoma (ECC) are poorly characterized. In this study, a cellular map of endometrial carcinoma was generated by profiling 30,780 cells isolated from tumor and paratumor tissues from five patients using single-cell RNA sequencing. 7 cell types in lymphocytes, 7 types in myeloid cells and 3 types in endometrial epithelial cells were identified. Distinct CD8⁺ T cell states and different monocyte-macrophage populations were discovered, among which exhausted CD8⁺ T cells and macrophages were preferentially enriched in tumor. Both CD8⁺ T cells and macrophages comport with continuous activation model. Gene expression patterns examination and gene ontology enrichment analysis of endometrial epithelial cells revealed 3 subtypes: stem-like cells, secretory glandular cells and ciliated cells. Overall, our study presents a view of endometrial carcinoma at single-cell resolution that reveals the characteristics of endometrial epithelial cells in the endometrium, and provides a cellular landscape of the tumor immune microenvironment.

Single-Cell Transcriptomic Analysis Identifies a Unique Pulmonary Lymphangioleiomyomatosis Cell

Rationale: Lymphangioleiomyomatosis (LAM) is a metastatic neoplasm of reproductive-age women associated with mutations in tuberous sclerosis complex genes. LAM causes cystic remodeling of the lung and progressive respiratory failure. The sources and cellular characteristics of LAM cells underlying disease pathogenesis remain elusive.Objectives: Identification and characterization of LAM cells in human lung and uterus using a single-cell approach.Methods: Single-cell and single-nuclei RNA sequencing on LAM (n = 4) and control (n = 7) lungs, immunofluorescence confocal microscopy, ELISA, and aptamer proteomics were used to identify and validate LAMCORE cells and secreted biomarkers, predict cellular origins, and define molecular and cellular networks in LAM.Measurements and Main Results: A unique cell type termed LAMCORE was identified, which was distinct from, but closely related to, lung mesenchymal cells. LAMCORE cells expressing signature genes included known LAM markers such as PMEL, FIGF, CTSK, and MLANA and novel biomarkers validated by aptamer screening, ELISA, and immunofluorescence microscopy. LAM cells in lung and uterus are morphologically indistinguishable and share similar gene expression profiles and biallelic TSC2 mutations, supporting a potential uterine origin for the LAMCORE cell. Effects of LAM on resident pulmonary cell types indicated recruitment and activation of lymphatic endothelial cells.Conclusions: A unique population of LAMCORE cells was identified in lung and uterus of patients with LAM, sharing close transcriptomic identity. LAM cell selective markers, secreted biomarkers, and the predicted cellular molecular features provide new insights into the signaling and transcriptional programs that may serve as diagnostic markers and therapeutic targets to influence the pathogenesis of LAM.

Modeling Hypoxia-Induced Neuropathies Using a Fast and Scalable Human Motor Neuron Differentiation System

Human motor neuron (MN) diseases encompass a spectrum of disorders. A critical barrier to dissecting disease mechanisms is the lack of appropriate human MN models. Here, we describe a scalable, suspension-based differentiation system to generate functional human MN diseases in 3 weeks. Using this model, we translated recent findings that mRNA mis-localization plays a role in disease development to the human context by establishing a membrane-based system that allows efficient fractionation of MN cell soma and neurites. In response to hypoxia, used to mimic diabetic neuropathies, MNs upregulated mitochondrial transcripts in neurites; however, mitochondria were decreased. These data suggest that hypoxia may disrupt translation of mitochondrial mRNA, potentially leading to neurite damage and development of neuropathies. We report the development of a novel human MN model system to investigate mechanisms of disease affecting soma and/or neurites that facilitates the rapid generation and testing of patient-specific MN diseases.

Estrogen receptor α (ERα)-binding super-enhancers drive key mediators that control uterine estrogen responses in mice

Estrogen receptor α (ERα) modulates gene expression by interacting with chromatin regions that are frequently distal from the promoters of estrogen-regulated genes. Active chromatin-enriched "super-enhancer" (SE) regions, mainly observed in in vitro culture systems, often control production of key cell type-determining transcription factors. Here, we defined super-enhancers that bind to ERα in vivo within hormone-responsive uterine tissue in mice. We found that SEs are already formed prior to estrogen exposure at the onset of puberty. The genes at SEs encoded critical developmental factors, including retinoic acid receptor α (RARA) and homeobox D (HOXD). Using high-throughput chromosome conformation capture (Hi-C) along with DNA sequence analysis, we demonstrate that most SEs are located at a chromatin loop end and that most uterine genes in loop ends associated with these SEs are regulated by estrogen. Although the SEs were formed before puberty, SE-associated genes acquired optimal ERα-dependent expression after reproductive maturity, indicating that pubertal processes that occur after SE assembly and ERα binding are needed for gene responses. Genes associated with these SEs affected key estrogen-mediated uterine functions, including transforming growth factor β (TGFβ) and LIF interleukin-6 family cytokine (LIF) signaling pathways. To the best of our knowledge, this is the first identification of SE interactions that underlie hormonal regulation of genes in uterine tissue and optimal development of estrogen responses in this tissue.

Uterine Patterning, Endometrial Gland Development, and Implantation Failure in Mice Exposed Neonatally to Genistein

BACKGROUND

Embryo implantation relies on precise hormonal regulation, associated gene expression changes, and appropriate female reproductive tract tissue architecture. Female mice exposed neonatally to the phytoestrogen genistein (GEN) at doses similar to those in infants consuming soy-based infant formulas are infertile due in part to uterine implantation defects.

OBJECTIVES

Our goal was to determine the mechanisms by which neonatal GEN exposure causes implantation defects.

METHODS

Female mice were exposed to GEN on postnatal days (PND)1-5 and uterine tissues collected on PND5, PND22-26, and during pregnancy. Analysis of tissue weights, morphology, and gene expression was performed using standard histology, confocal imaging with three-dimensional analysis, real-time reverse transcription polymerase chain reaction (real-time RT-PCR), and microarrays. The response of ovariectomized adults to 17 β -estradiol (E2) and artificial decidualization were measured. Leukemia inhibitory factor (LIF) injections were given intraperitoneally and implantation sites visualized. Gene expression patterns were compared with curated data sets to identify upstream regulators.

RESULTS

GEN-exposed mice exhibited reduced uterine weight gain in response to E2 treatment or artificial decidualization compared with controls; however, expression of select hormone responsive genes remained similar between the two groups. Uteri from pregnant GEN-exposed mice were posteriorized and had reduced glandular epithelium. Implantation failure was not rescued by LIF administration. Microarray analysis of GEN-exposed uteri during early pregnancy revealed significant overlap with several conditional uterine knockout mouse models, including Foxa2, Wnt4, and Sox17. These models exhibit reduced endometrial glands, features of posteriorization and implantation failure. Expression of Foxa2, Wnt4, and Sox17, as well as genes important for neonatal uterine differentiation (Wnt7a, Hoxa10, and Msx2), were severely disrupted on PND5 in GEN-exposed mice.

DISCUSSION

Our findings suggest that neonatal GEN exposure in mice disrupts expression of genes important for uterine development, causing posteriorization and diminished gland function during pregnancy that contribute to implantation failure. These findings could have implications for women who consumed soy-based formulas as infants. https://doi.org/10.1289/EHP6336.

Mice lacking uterine enhancer of zeste homolog 2 have transcriptomic changes associated with uterine epithelial proliferation

Enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase that suppresses gene expression. Previously, we developed a conditional null model where EZH2 is knocked out in uterus. Deletion of uterine EZH2 increased proliferation of luminal and glandular epithelial cells. Herein, we used RNA-Seq in wild-type (WT) and EZH2 conditional knockout (Ezh2cKO) uteri to obtain mechanistic insights into the gene expression changes that underpin the pathogenesis observed in these mice. Ovariectomized adult Ezh2cKO mice were treated with vehicle (V) or 17β-estradiol (E2; 1 ng/g). Uteri were collected at postnatal day (PND) 75 for RNA-Seq or immunostaining for epithelial proliferation. Weighted gene coexpression network analysis was used to link uterine gene expression patterns and epithelial proliferation. In V-treated mice, 88 transcripts were differentially expressed (DEG) in Ezh2cKO mice, and Bmp5, Crabp2, Lgr5, and Sprr2f were upregulated. E2 treatment resulted in 40 DEG with Krt5, Krt15, Olig3, Crabp1, and Serpinb7 upregulated in Ezh2cKO compared with control mice. Transcript analysis relative to proliferation rates revealed two module eigengenes correlated with epithelial proliferation in WT V vs. Ezh2cKO V and WT E2 vs. Ezh2cKO E2 mice, with a positive relationship in the former and inverse in the latter. Notably, the ESR1, Wnt, and Hippo signaling pathways were among those functionally enriched in Ezh2cKO females. Current results reveal unique gene expression patterns in Ezh2cKO uterus and provide insight into how loss of this critical epigenetic regulator assumingly contributes to uterine abnormalities.

Molecular Studies on Pregnancy with Mouse Models

Pregnancy is a complex process that involves crosstalk among multiple cell types in both the endometrial and myometrial compartments at the maternal side to support the fetus. Genetic engineered mouse models have served as a major platform to dissect the convolute genetic interactions in a physiological context. Combining with various applications of next generation sequencing and genome editing, functional assays by mouse models have expanded the spectrum to include both coding and noncoding genome. The present review will highlight recent findings that are primarily based on studies of mouse models with emphasis on pathways for endometrial receptivity and myometrial contraction. Emerging novel technologies that may advance the research in these two aspects will also be discussed.

Self-renewing endometrial epithelial organoids of the human uterus

The human endometrium is essential in providing the site for implantation and maintaining the growth and survival of the conceptus. An unreceptive endometrium and disrupted maternal-conceptus interactions can cause infertility due to pregnancy loss or later pregnancy complications. Despite this, the role of uterine glands in first trimester human pregnancy is little understood. An established organoid protocol was used to generate and comprehensively analyze 3-dimensional endometrial epithelial organoid (EEO) cultures from human endometrial biopsies. The derived EEO expand long-term, are genetically stable, and can be cryopreserved. Using endometrium from 2 different donors, EEO were derived and then treated with estrogen (E2) for 2 d or E2 and medroxyprogesterone acetate (MPA) for 6 d. EEO cells were positive for the gland marker, FOXA2, and exhibited appropriate hormonal regulation of steroid hormone receptor expression. Real-time qPCR and bulk RNA-sequencing analysis revealed effects of hormone treatment on gene expression that recapitulated changes in proliferative and secretory phase endometrium. Single-cell RNA sequencing analysis revealed that several different epithelial cell types are present in the EEO whose proportion and gene expression changed with hormone treatment. The EEO model serves as an important platform for studying the physiology and pathology of the human endometrium.

A tale of two tracts: history, current advances, and future directions of research on sexual differentiation of reproductive tracts†

Alfred Jost's work in the 1940s laid the foundation of the current paradigm of sexual differentiation of reproductive tracts, which contends that testicular hormones drive the male patterning of reproductive tract system whereas the female phenotype arises by default. Once established, the sex-specific reproductive tracts undergo morphogenesis, giving rise to anatomically and functionally distinct tubular organs along the rostral-caudal axis. Impairment of sexual differentiation of reproductive tracts by genetic alteration and environmental exposure are the main causes of disorders of sex development, and infertility at adulthood. This review covers past and present work on sexual differentiation and morphogenesis of reproductive tracts, associated human disorders, and emerging technologies that have made impacts or could radically expand our knowledge in this field.

Single-cell sequencing of neonatal uterus reveals an Misr2+ endometrial progenitor indispensable for fertility

The Mullerian ducts are the anlagen of the female reproductive tract, which regress in the male fetus in response to MIS. This process is driven by subluminal mesenchymal cells expressing Misr2, which trigger the regression of the adjacent Mullerian ductal epithelium. In females, these Misr2+ cells are retained, yet their contribution to the development of the uterus remains unknown. Here, we report that subluminal Misr2+ cells persist postnatally in the uterus of rodents, but recede by week 37 of gestation in humans. Using single-cell RNA sequencing, we demonstrate that ectopic postnatal MIS administration inhibits these cells and prevents the formation of endometrial stroma in rodents, suggesting a progenitor function. Exposure to MIS during the first six days of life, by inhibiting specification of the stroma, dysregulates paracrine signals necessary for uterine development, eventually resulting in apoptosis of the Misr2+ cells, uterine hypoplasia, and complete infertility in the adult female.

In the womb, mammals possess all of the preliminary sexual structures necessary to become either male or female. This includes the Mullerian duct, which develops into the Fallopian tubes, uterus, cervix, and vagina in female fetuses. In male fetuses, the testis secretes a hormone called Mullerian inhibiting substance (MIS). This triggers the activity of a small group of cells, known as Misr2+ cells, that cause the Mullerian duct to degenerate, preventing males from developing female sexual organs.

It was not clear what happens to Misr2+ cells in female fetuses or if they affect how the uterus develops. Saatcioglu et al. now show that in newborn female mice and rats, a type of Misr2+ cell that sits within a thin inner layer of the developing uterus still responds to MIS. At this time, the uterus is in a critical early period of development. Treating the mice and rats with MIS protein during their first six days of life eventually caused the Misr2+ cells to die. The treatment also prevented a layer of connective tissue, known as the endometrial stroma, from forming in the uterus. As a result, the mice and rats were infertile and had severely underdeveloped uteri.

While the Misr2+ cells are present in newborn rats and mice, Saatcioglu et al. found that they disappeared before birth in humans. However, the overall results suggest that Misr2+ cells act as progenitor cells that develop into the cells of the endometrial stroma. Future work could investigate the roles these cells play in causing uterine developmental disorders and infertility disorders. Furthermore, the finding that MIS inhibits the Misr2+ cells could help researchers to develop treatments for uterine cancer and other conditions where the cells of the uterus grow and divide too much.