Decidual remodelling of midluteal endometrium leads to a short implantation window after which the uterine mucosa either breaks down or is transformed into a robust matrix that accommodates the placenta throughout pregnancy. To gain insights into the underlying mechanisms, we established and characterised endometrial assembloids, consisting of gland-like organoids and primary stromal cells. Single-cell transcriptomics revealed that decidualized assembloids closely resemble midluteal endometrium, harbouring differentiated and senescent subpopulations in both glands and stroma. We show that acute senescence in glandular epithelium drives secretion of multiple canonical implantation factors, whereas in the stroma it calibrates the emergence of anti-inflammatory decidual cells and pro-inflammatory senescent decidual cells. Pharmacological inhibition of stress responses in pre-decidual cells accelerated decidualization by eliminating the emergence of senescent decidual cells. In co-culture experiments, accelerated decidualization resulted in entrapment of collapsed human blastocysts in a robust, static decidual matrix. By contrast, the presence of senescent decidual cells created a dynamic implantation environment, enabling embryo expansion and attachment, although their persistence led to gradual disintegration of assembloids. Our findings suggest that decidual senescence controls endometrial fate decisions at implantation and highlight how endometrial assembloids may accelerate the discovery of new treatments to prevent reproductive failure.

STUDY QUESTION Can the accuracy of timing of luteal phase endometrial biopsies based on urinary ovulation testing be improved by measuring the expression of a small number of genes and a continuous, non-categorical modelling approach? SUMMARY ANSWER Measuring the expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3 and SLC15A2) is sufficient to obtain substantially more accurate timing estimates and to assess the reliability of timing estimates for each sample. WHAT IS KNOWN ALREADY Commercially available endometrial timing approaches based on gene expression require large gene sets and use a categorical approach that classifies samples as pre-receptive, receptive or post-receptive. STUDY DESIGN, SIZE, DURATION Gene expression was measured by RTq-PCR in different sample sets, comprising a total of 664 endometrial biopsies obtained 4–12 days after a self-reported positive home ovulation test. A further 36 endometrial samples were profiled by RTq-PCR as well as RNA-sequencing. PARTICIPANTS/MATERIALS, SETTING, METHODS A computational procedure, named ‘EndoTime’, was established that models the temporal profile of each gene and estimates the timing of each sample. Iterating these steps, temporal profiles are gradually refined as sample timings are being updated, and confidence in timing estimates is increased. After convergence, the method reports updated timing estimates for each sample while preserving the overall distribution of time points. MAIN RESULTS AND THE ROLE OF CHANCE The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3: P < 0.005; CXCL14: P < 2.7e−6; DPP4: P < 3.7e−13). Pearson correlation between EndoTime estimates for the same sample set but based on RTq-PCR or RNA-sequencing data showed a high degree of congruency between the two (P = 8.6e−10, R2 = 0.687). Estimated timings did not differ significantly between control subjects and patients with recurrent pregnancy loss or recurrent implantation failure (P > 0.05). LARGE SCALE DATA The RTq-PCR data files are available via the GitHub repository for the EndoTime software at https://github.com/AE-Mitchell/EndoTime, as is the code used for pre-processing of RTq-PCR data. The RNA-sequencing data are available on GEO (accession GSE180485). LIMITATIONS, REASONS FOR CAUTION Timing estimates are informed by glandular gene expression and will only represent the temporal state of other endometrial cell types if in synchrony with the epithelium. Methods that estimate the day of ovulation are still required as these data are essential inputs in our method. Our approach, in its current iteration, performs batch correction such that larger sample batches impart greater accuracy to timing estimations. In theory, our method requires endometrial samples obtained at different days in the luteal phase. In practice, however, this is not a concern as timings based on urinary ovulation testing are associated with a sufficient level of noise to ensure that a variety of time points will be sampled. WIDER IMPLICATIONS OF THE FINDINGS Our method is the first to assay the temporal state of luteal-phase endometrial samples on a continuous domain. It is freely available with fully shared data and open-source software. EndoTime enables accurate temporal profiling of any gene in luteal endometrial samples for a wide range of research applications and, potentially, clinical use. STUDY FUNDING/COMPETING INTEREST(S) This study was supported by a Wellcome Trust Investigator Award (Grant/Award Number: 212233/Z/18/Z) and the Tommy's National Miscarriage Research Centre. None of the authors have any competing interests. J.L. was funded by the Biotechnology and Biological Sciences Research Council (UK) through the Midlands Integrative Biology Training Partnership (MIBTP, BB/M01116X/1).

Despite advances in assisted reproductive techniques in the four decades since the first human birth after in vitro fertilisation, 1-2 % of couples experience recurrent implantation failure, and some will never achieve a successful pregnancy even in the absence of a confirmed dysfunction. Furthermore, 1-2 % of couples who do conceive, either naturally or with assistance, will experience recurrent early loss of karyotypically normal pregnancies. In both cases, embryo-endometrial interaction is a clear candidate for exploration. The impossibility of studying implantation processes within the human body has necessitated the use of animal models and cell culture approaches. Recent advances in 3-dimensional modelling techniques, namely the advent of organoids, present an exciting opportunity to elucidate the unanswerable within human reproduction. In this review, we will explore the ontogeny of implantation modelling and propose a roadmap to application and discovery.

Decidual remodelling of midluteal endometrium leads to a short implantation window after which the uterine mucosa either breaks down or is transformed into a robust matrix that accommodates the placenta throughout pregnancy. To gain insights into the underlying mechanisms, we established and characterised endometrial assembloids, consisting of gland organoids and primary stromal cells. Single-cell transcriptomics revealed that decidualized assembloids closely resemble midluteal endometrium, harbouring differentiated and senescent subpopulations in both glands and stroma. We show that acute senescence in glandular epithelium drives secretion of multiple canonical implantation factors, whereas in the stroma it calibrates the emergence of anti-inflammatory decidual cells and pro-inflammatory senescent decidual cells. Pharmacological inhibition of stress responses in pre-decidual cells accelerated decidualization by inhibiting senescence and mesenchymal-epithelial transition, processes involved in endometrial breakdown and regeneration, respectively. Accelerated decidualization resulted in entrapment of co-cultured human blastocysts in a largely static decidual matrix. By contrast, the presence of senescent decidual cells created a dynamic implantation environment, enabling embryo expansion and attachment, although their persistence led to gradual disintegration of assembloids. Our findings demonstrate that senescence controls endometrial fate decisions at implantation and highlight how endometrial assembloids may accelerate the discovery of new treatments to prevent reproductive failure.

Study question Can expression measurements of a small panel of genes be used to develop a continuous, non-categorical model for the improvement of endometrial biopsy timing accuracy? Summary answer Measuring expression levels of six genes (IL2RB, IGFBP1, CXCL14, DPP4, GPX3, and SLC15A2) is sufficient to obtain and assess substantially more accurate timing estimates. What is known already Commercially available endometrial timing approaches based on gene expression require much larger gene sets and use a categorical approach that classifies samples as pre-receptive, receptive, or post-receptive. Study design, size, duration Gene expression was measured by RT-qPCR in different sample sets, comprising a total of 664 endometrial biopsies obtained 4 to 12 days after a self-reported positive home ovulation test. A further 36 endometrial samples were profiled by RT-qPCR as well as RNA-sequencing. Participants/materials, setting, methods A computational procedure, named ‘EndoTime’, was established that models the temporal profile of each gene and estimates the timing of each sample. Iterating these steps, temporal profiles are gradually refined as sample timings are being updated, and confidence in timing estimates is increased. After convergence, the method reports updated timing estimates for each sample while preserving the overall distribution of time points. Main results and the role of chance The Wilcoxon rank-sum test was used to confirm that ordering samples by EndoTime estimates yields sharper temporal expression profiles for held-out genes (not used when determining sample timings) than ordering the same expression values by patient-reported times (GPX3 : p < 0.005; CXCL14 : p < 2.7e-6; DPP4 : p < 3.7e-13). Pearson correlation between EndoTime estimates for the same sample set but based on RT-qPCR or RNA-sequencing data showed high degree of congruency between the two (p = 8.6e-10, R-squared = 0.687). Estimated timings did not differ significantly between control subjects and patients with recurrent pregnancy loss or recurrent implantation failure (p > 0.05). Limitations, reasons for caution Timing estimates are informed by glandular gene expression and will only represent the temporal state of other endometrial cell types if in synchrony with the epithelium. Additionally, methods that estimate the day of ovulation are still required as these data are essential inputs in our method. Wider implications of the findings Our method is the first to assay the temporal state of endometrial samples on a continuous domain, enabling accurate temporal profiling of any gene in luteal phase biopsies for a wide range of research applications and, potentially, clinical use. It is freely available, open-source software including supporting data sets. Trial registration number Not applicable