
==== Front
J Assist Reprod Genet
J Assist Reprod Genet
Journal of Assisted Reproduction and Genetics
1058-0468
1573-7330
Springer US New York

38995509
3173
10.1007/s10815-024-03173-y
Reproductive Physiology and Disease
Enrichment of cell cycle pathways in progesterone-treated endometrial organoids of infertile women compared to fertile women
http://orcid.org/0000-0003-3542-6215
Bui B. N. b.n.bui@umcutrecht.nl

1
http://orcid.org/0000-0001-8527-4585
Ardisasmita A. I. 2
van de Vliert F. H. 1
http://orcid.org/0009-0004-4350-3449
Abendroth M. S. 1
van Hoesel M. 1
http://orcid.org/0000-0001-7012-0436
Mackens S. 4
http://orcid.org/0000-0001-9147-2406
Fuchs S. A. 2
http://orcid.org/0000-0001-8094-1211
Nieuwenhuis E. E. S. 35
http://orcid.org/0000-0003-1366-2727
Broekmans F. J. M. 16
http://orcid.org/0000-0003-2652-0760
Steba G. S. 1
1 https://ror.org/0575yy874 grid.7692.a 0000 0000 9012 6352 Department of Gynaecology and Reproductive Medicine, University Medical Center Utrecht, Heidelberglaan 100, 3584 CX Utrecht, The Netherlands
2 https://ror.org/0575yy874 grid.7692.a 0000 0000 9012 6352 Department of Metabolic Diseases, University Medical Center Utrecht, Heidelberglaan 100, 3584 CX Utrecht, The Netherlands
3 grid.417100.3 0000 0004 0620 3132 Department of Pediatric Gastroenterology, Wilhelmina Children’s Hospital, University Medical Center Utrecht, Heidelberglaan 100, 3584 CX Utrecht, The Netherlands
4 grid.8767.e 0000 0001 2290 8069 Brussels IVF, Universitair Ziekenhuis Brussel, Vrije Universiteit Brussel, Laarbeeklaan 101, 1090 Brussels, Belgium
5 grid.5477.1 0000000120346234 Department of Science, University College Roosevelt, Lange Noordstraat 1, 4331 CB Middelburg, The Netherlands
6 Centre for Infertility Care, Dijklander Ziekenhuis, Purmerend, The Netherlands
12 7 2024
12 7 2024
9 2024
41 9 24052418
2 2 2024
13 6 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Purpose

To investigate whether the transcriptome profile differs between progesterone-treated infertile and fertile endometrial organoids.

Methods

Endometrial biopsies were obtained from 14 infertile and seven fertile women, after which organoids were generated from isolated epithelial cells. To mimic the secretory phase, organoids were sequentially treated with 17β-estradiol (E2) and progesterone (P4) and subjected to RNA sequencing. Differentially expressed genes (DEGs) were identified using DESeq2 (lfcThreshold = 0, log2 Fold Change ≥ 1.0 or ≤ −1.0), and a principal component analysis (PCA) plot was generated. Functional enrichment analysis was performed by overrepresentation analysis and Gene Set Enrichment Analysis (GSEA). To functionally assess proliferation, OrganoSeg surface measurements were performed before (T0) and after (T1) differentiation of organoids, and T1/T0 ratios were calculated to determine the proliferation rate.

Results

Although the PCA plot did not show clear clustering of the fertile and infertile samples, 363 significant DEGs (129 upregulated and 234 downregulated) were detected in infertile compared to fertile organoids. Mainly cell cycle processes were highly enriched in infertile organoids. Thus, we hypothesised that proliferative activity during differentiation may be higher in infertile organoids compared to fertile organoids. However, this could not be validated by cell surface measurements.

Conclusions

This study revealed that cell cycle processes were enriched in E2/P4-treated infertile endometrial organoids as compared to fertile organoids. This could reflect persistently higher proliferative activity of the endometrial epithelial cells in differentiated infertile organoids compared to fertile organoids. To confirm this hypothesis, further studies are warranted.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10815-024-03173-y.

Keywords

Endometrium
Organoids
Transcriptome
RNA-sequencing
Infertility
Fertility
issue-copyright-statement© Springer Science+Business Media, LLC, part of Springer Nature 2024
==== Body
pmcIntroduction

Infertility is defined as the failure to achieve a clinical pregnancy after 12 months or more of regular unprotected sexual intercourse [1]. Successful embryo implantation requires a viable blastocyst, a receptive endometrium, and a synchronised interaction between both [2]. In general, blastocysts are considered ‘viable’ when they demonstrate proper development in vitro and thus have the potential for successful pregnancy when transferred into the uterus. However, currently, there is no definitive method available to reliably predict whether an embryo will result in a successful pregnancy.

The endometrium is the inner lining of the uterus and plays a key role in reproduction [3]. The superficial functional layer sheds during menstruation and regenerates in the subsequent cycle, whereas the deeper basal layer harbours the endometrial progenitor stem cells from which regeneration of the endometrium occurs after each menstruation [4]. Histologically, the endometrium is composed of a single layer of luminal epithelial cells, covering a multicellular stromal layer with connective tissue, fibroblast-like stromal cells, tubular glands extending from the luminal epithelium, spiral arteries and immune cells [3]. The cyclic growth, differentiation and breakdown of the endometrium is regulated by the major ovarian sex steroids 17β-estradiol (E2) and progesterone (P4). The menstrual cycle starts with the proliferative phase, which is dominated by follicular E2 production and characterised by extensive proliferation of the endometrium. Following ovulation, P4 is produced by the corpus luteum, heralding the start of the secretory phase, in which the endometrium prepares for implantation [3, 5].

Previous endometrial transcriptome studies based on whole tissue biopsies have identified a vast number of differentially expressed genes (DEGs) in the endometrium of women with recurrent implantation failure (RIF) compared to fertile controls [6–10], suggesting differences in endometrial function. However, there is still an inadequate understanding of the potential role of the endometrium in implantation failure due to the inclusion of heterogeneous study populations, as there is no universally accepted definition of RIF [11]. Furthermore, difficulties to control for exogenous factors in vivo and heterogeneous cellular composition of the endometrial biopsy [12] may have contributed to poor concordance of findings across studies.

Functional endometrial in vivo studies are difficult due to ethical and practical reasons. Therefore, models mimicking the biology of the endometrium are indispensable. Although animal models have provided important insights [13], these do not completely resemble the reproductive tract of humans. Thus, in vitro models of the human endometrium are of great value to further study endometrial function and its role in implantation. The establishment of the three-dimensional (3D) endometrial organoid model has enabled researchers to study endometrial cells that more closely mimic in vivo physiology [14, 15]. Organoids are 3D in vitro structures that self-organise from tissue or stem cells under specific culture conditions and recapitulate the original organ’s key biological properties, such as its function and micro-anatomy [16, 17]. Endometrial organoids have proven to be robustly expandable while recapitulating the original donor’s characteristics and remaining phenotypically and genetically stable over multiple expansions [15], even after tissue cryopreservation [18].

The transcriptome of endometrial organoids has previously been analysed by both bulk and single-cell RNA-sequencing (RNA-seq), which revealed the effects of hormonal treatment on gene expression in endometrial organoids, the existence of several different epithelial cell types as well as important pathways determining the cell fate of these epithelial lineages [19–22]. There are no studies yet comparing the transcriptome of endometrial organoids derived from infertile women versus fertile women, which could provide insight into potential differences in reproductive endometrial function. The objective of this study was to compare the transcriptome profile, analysed by RNA-seq, of endometrial organoids of infertile and fertile women.

Materials and methods

Study population

Infertile women were defined as not having reached a clinical pregnancy (i.e. gestational sac visualised by ultrasound) after at least 12 months of regular unprotected intercourse. Endometrial tissue of infertile women was obtained within two randomised controlled trials (RCTs) on endometrial scratching (i.e. biopsy) (SCRaTCH-OFO trial and SCRaTCH trial) [23, 24], as well as an endometrium biobanking study (ENORM study). The full inclusion and exclusion criteria of the two RCTs have been described in detail elsewhere [23, 24]. To summarize, eligible participants for the SCRaTCH-OFO trial [24] were women aged between 18 and 38 years from couples diagnosed with unexplained infertility, who were trying to conceive spontaneously, with a predicted natural conception chance of ≥ 30% by the Hunault model [25]. Whereas for the SCRaTCH trial, eligible participants were women aged between 18 and 44 years who had failed implantation after one full in vitro fertilisation (IVF) or intracytoplasmic sperm injection (ICSI) cycle (i.e. after transfer of fresh and/or frozen embryos) with at least one embryo transfer, and who were planning a new IVF/ICSI cycle [23]. Within the ENORM study, endometrial tissue was obtained from women with RIF, who were defined as having had at least three failed embryo transfers.

Fertile women were defined as having conceived and reached a clinical pregnancy within 12 months of regular unprotected intercourse or intrauterine insemination (IUI) with donor sperm. Endometrial tissue of fertile women was also obtained within the ENORM study. Fertile women were eligible for this study if they fulfilled the following inclusion criteria: (1) aged between 18 and 38 years; (2) no previous infertility or miscarriage and (3) previous spontaneous conception within 12 months with delivery of a healthy baby at term; or undergoing IUI with donor sperm or preimplantation genetic testing (PGT) cycles without ever having conceived before, but not being classified as infertile yet (i.e. not having undergone at least 12 months of attempts) or undergoing IUI with donor sperm or PGT cycles with a previous spontaneous conception within 12 months with delivery of a healthy baby at term. Exclusion criteria were as follows: (1) a history of lower abdominal or pelvic infection; (2) an increased risk of intra-abdominal infection (e.g. due to intestinal surgery in the past or due to the presence of an immunological disease); (3) endometriosis stage III-IV according to the ASRM classification [26]; (4) previous caesarean section with niche formation; (5) the presence of untreated unilateral or bilateral hydrosalpinx; (6) increased risk of bleeding (e.g. the presence of severe coagulation disorders, such as haemophilia); (7) the presence of intrauterine adhesions; (8) previous intrauterine procedures (e.g. curettage, endometrial polyp resection, myomectomy, uterine septum resection); (9) having experienced one or more of the following complications during a previous pregnancy: preterm birth (i.e. before 37 weeks of gestational age), intrauterine growth restriction (i.e. birth weight below the 2.3th percentile), pregnancy induced hypertension, pre-eclampsia, HELLP syndrome, placenta accreta and manual removal of the placenta.

Endometrial biopsy, tissue processing and cryopreservation

Endometrial tissue was obtained by endometrial biopsy, using an endometrial biopsy catheter (Pipelle catheter). Our protocol for endometrial biopsy has been previously reported [27] and is briefly described below.

Endometrial biopsy was performed either in the mid-luteal phase of a natural cycle, 5 to 8 days after detection of the luteinising hormone (LH) surge by urinary tests (n = 16); or in the early luteal phase of a natural cycle (LH + 4; n = 1); or during contraceptive pill (n = 2) or vaginal ring use (n = 1); or during a mock artificial cycle on the 6th day of P4 supplementation (n = 1).

After insertion of a speculum, the cervix was cleaned with sterile water. The endometrial biopsy catheter was introduced through the cervix up to the uterine fundus. The piston of the catheter was drawn back to create a vacuum, and the catheter was slowly retracted in a time period of 1 to 2 min while constantly rotating 360°.

After the procedure, the tissue was cryopreserved for organoid development as previously described [18]. Briefly, the tissue was collected in collection medium (Advanced Dulbecco Modified Eagle’s Medium/Nutrient Mixture F-12 (DMEM/F-12) (Thermo Fisher Scientific, USA) with 5% foetal bovine serum (FBS) (Sigma–Aldrich, USA) and 2% penicillin–streptomycin (10,000 U/ml) (Thermo Fisher Scientific) and transported to the laboratory on ice for further processing in a laminar flow hood. After washing with Ca2+/Mg2+-free phosphate buffered saline (PBS) (Thermo Fisher Scientific) to remove as many blood clots as possible, the tissue was cut into small pieces of 1 mm3 and divided into two parts. Both tissue parts were collected in two cryovials containing 1 ml of cryopreservation medium each (Advanced DMEM/F-12 with 30% FBS and 10% dimethyl sulfoxide (DMSO) (Sigma–Aldrich) and slowly frozen to −80 °C overnight in a Mr. Frosty Freezing Container (Thermo Fisher Scientific). Subsequently, the two cryovials were stored in liquid nitrogen at − 196 °C until organoid development.

Endometrial organoid culture

Organoid culture from cryopreserved endometrial biopsies

The cryopreserved endometrial tissue was rapidly thawed at 37 °C and diluted in.

Advanced DMEM/F-12 with 10% FBS (1:10). After centrifugation (at 200 × g for 5 min) and removal of the supernatant, tissue dissociation, organoid culturing and passaging were performed as previously described [15].

In short, endometrial tissue was washed in PBS, further cut into finer pieces and incubated in 1 mg ml−1 collagenase IV (Thermo Fisher Scientific) with 10 μM Y-27632 (Merck Millipore, USA) at 37 °C for 1 to 2 h with gentle agitation. Tissue pieces were mechanically triturated every 20 min by pipetting up and down. The digestion was inactivated by 1:1 dilution with Advanced DMEM/F-12 with 10% FBS. The suspension with digested tissue pieces was run over a 40-µm cell strainer (Greiner Bio-One, Austria). The strainer was placed upside down in a petri dish and rinsed with Advanced DMEM/F-12 to remove the endometrial epithelial cells and glandular fragments. The medium was collected from the petri dish, and after centrifugation (200 × g, 5 min, 4 °C), the pellet of endometrial glandular fragments was resuspended in 70% Matrigel (Corning Incorporated, USA) and 30% DMEM/F-12 with 10 μM Y-27632. The cell-Matrigel suspension was plated in a 24-well cell culture plate (Corning Incorporated) with two drops of 15 µl per well. The plate was incubated at 37 °C and 5% CO2 for 20 min to allow the Matrigel drops to solidify. Pre-warmed endometrial organoid base culture medium (SupplementaryTable 1) was added with 10 μM Y-27632 once and refreshed every 48 h. After 96 h of culture, the organoid culture medium was supplemented with 1 nM 17β-estradiol (E2) (Sigma–Aldrich) (SupplementaryTable 1) and medium refreshments with E2 organoid culture medium continued every 48 h.

Passaging of endometrial organoids

Organoids were passaged every 10–14 days in a ratio of 1:4 to 1:5. Briefly, the organoid-containing Matrigel drops were dissolved in ice-cold Advanced DMEM/F-12, centrifugated (200 × g, 5 min, 4 °C), mechanically triturated, resuspended in 70% Matrigel and 30% Advanced DMEM/F-12, plated in a 24-wells cell culture plate and incubated at 37 °C and 5% CO2. Further maintenance of the organoid culture was performed as described above.

Hormonal treatment of endometrial organoids

Endometrial organoids between passages 2 and 7 were treated with hormones to mimic the secretory phase for further analysis by RNA-seq (Fig. 1a). During passaging, organoids of each organoid line were seeded in four wells of a 24-well cell culture plate, containing two drops of 15 µl per well. Organoids were first treated with organoid base culture medium without any hormones for 96 h, followed by base medium supplemented with 1 nM E2 for 96 h and finally base medium supplemented with 0.1 nM E2 and 200 ng ml−1 P4 (Sigma–Aldrich) for another 96 h (SupplementaryTable 1,Fig. 1a). Organoids were harvested and pooled in two conical tubes of 1.5 ml (Eppendorf, USA). After centrifugation at 200 × g for 5 min, the supernatant was removed and the cell pellets were stored at −80 °C.Fig. 1 Proliferation assay of infertile and fertile endometrial organoids. a Hormonal treatment of endometrial organoids. T−1 represented the start of estradiol (E2) supplementation, T0 the start of progesterone (P4) supplementation for differentiation of the organoids and T1 after differentiation. Photos were taken at the three-time points for evaluation in OrganoSeg. Two organoid lines received a second medium refreshment with E2 medium on days 8–9 and therefore P4 supplementation started 2–3 days later in these lines. BM, base medium; E2, estradiol medium; P4, progesterone medium. b Analysis of the two-dimensional surfaces of endometrial organoids in OrganoSeg, using the brightfield images taken at T0 and T1. The left image shows the original image, the middle image highlights the organoid borders and the right image delineates the measured organoid surfaces

RNA extraction

RNA extraction of the organoid cell pellets was performed by the Utrecht Sequencing Facility (USEQ) (Utrecht, The Netherlands). Cell pellets were thawed and homogenised in RLT Plus lysis buffer (Qiagen). Total RNA extraction was performed using the QIAsymphony SP and the QIAsymphony RNA Kit (Qiagen) following the manufacturer’s protocol. RNA was quantified using the Qubit™ RNA BR Assay Kit (Thermo Fisher Scientific) and the Qubit 4 Fluorometer (Thermo Fisher Scientific). RNA Integrity Number (RIN) was determined using the Agilent RNA 6000 Nano Kit (Agilent Technologies, USA) and the Agilent 2100 Bioanalyzer system (Agilent Technologies).

RNA-sequencing

RNA-seq was performed by USEQ in three runs. All endometrial organoid RNA samples (n = 21 lines in duplicate) had a RIN of ≥ 6.5. Of each sample, an input of 100 ng of total RNA was used to construct a library. Libraries were prepared using the TruSeq Stranded mRNA Kit (Illumina, USA) following the manufacturer’s protocol. xGen Dual Index UMI Adapters (Integrated DNA Technologies (IDT), USA) were used for indexing. Libraries had a final concentration of ≥ 1 ng/µl for sequencing, were pooled equimolarly and single-end (1 × 75 base pairs) sequenced on a NextSeq500 (Illumina).

RNA-sequencing analysis

Quality control of the raw sequence reads was performed with FastQC (v0.11.8) [28]. TrimGalore (v0.6.5) [29] was used for read trimming based on quality and adapter presence, after which quality control was performed again with FastQC. Ribosomal RNA reads were filtered out using SortMeRNA (v4.3.3) [30]. The resulting reads were aligned to the reference genome hg19 (GRCh37) using the STAR (v2.7.3a) aligner [31]. Follow-up quality control of the mapped (bam) files was performed using Sambamba (v0.7.0) [32], RSeQC (v3.0.1) [33] and PreSeq (v2.0.3) [34]. Read counts were generated using the Subread featureCounts module (v2.0.0) [35]. Batch effect adjustment was performed with ComBat-seq (v3.38.0) [36]. Differential gene expression analysis was performed with DESeq2 (v1.30.1) [37]. Normalised counts were generated by applying DESeq2 variance-stabilizing transformation (VST) to the read counts. To identify significantly differentially expressed genes (DEGs), a lfcThreshold of 0 was used in DESeq2 as a significance threshold, identifying all significantly expressed genes with a log2 Fold Change (Log2FC) of ≥ 1.0 or ≤ −1.0. Principal component analysis (PCA) was performed using normalised counts and plotted using ggplot2 (v3.3.3) [38].

Enrichment analysis

Significantly DEGs were subjected to overrepresentation analysis using Enrichr (https://maayanlab.cloud/Enrichr/) [39, 40], annotating Gene Ontology (GO) terms and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways to gene sets. Enrichr uses the Benjamini–Hochberg method to correct for multiple testing.

In addition, gene set enrichment analysis (GSEA) was performed using the GSEA software from the Broad Institute (v4.1.0, [41, 42]) in order to identify overrepresented GO gene sets. The pre-ranked enrichment method was applied, using a ranked list of DEGs as input for the GSEA-Pre-ranked tool. The DEGs were ranked based on the value obtained from the following formula: sign(log2foldchange) ×  − log10(pval). The analysis was performed using the following parameters: number of permutations = 1000; enrichment statistic = classic; max/min size = 500/15; normalisation method = meandiv. The Normalized Enrichment Score (NES) was used to indicate the degree to which the genes in a gene set were overrepresented.

Proliferation assay

The proliferation rate was functionally assessed using seven infertile and five fertile lines that had also been selected for RNA-seq. After passaging, the organoid fragments of each organoid line were seeded in 6 wells of a 24-well cell culture plate with a single drop of 20 µl in each well. Organoids were passaged in a ratio of 1:3 to 1:8, depending on the density of the organoids. During seeding of the drops, the wells of the upper and lower row of the cell culture plate were avoided to limit edge effects as much as possible. The organoids were cultured as described above under the “Organoid culture from cryopreserved endometrial biopsies” section.

Photos were taken of all wells using a brightfield microscope (Evos Cell Imaging, Thermo Fisher Scientific) with the 1.25 × objective at the following time points: prior to E2 supplementation (i.e. time point − 1, T−1), at the start of E2/P4 supplementation (i.e. T0) and at the end of E2/P4 differentiation (i.e. T1) (Fig. 1a). Photos were analysed using OrganoSeg, a software programme specifically developed to align spheroid figures in brightfield images of organoids [43] (Fig. 1b). OrganoSeg settings were uniformly configured for all individual wells within an organoid line. However, due to differences in organoid density across the different organoid lines, the settings differed between organoid lines. The surfaces of the organoids in pixels determined at T0 and T1 were used to calculate a ratio (T1/T0) to express the proliferation rate during differentiation.

Statistical analysis

Continuous data were statistically compared between two groups using the independent t-test for normally distributed data and the Wilcoxon rank sum test for non-normally distributed data. Categorical data were compared between groups using the chi-square test or Fisher’s exact test in case of an expected frequency of < 5. A p-value < 0.05 was considered to indicate a statistically significant difference. Statistical analysis was carried out using SPSS Statistics version 26 (IBM Corporation, USA).

Results

Demographics

Endometrial organoids were established from endometrial tissue obtained from seven fertile women and 14 infertile women. Clinical characteristics are summarized in Tables 1 and 2. The median age of both groups was respectively 35 and 34.5 years and did not differ significantly (p = 0.84). There were no significant differences in the median BMI and proportion of smokers between both groups (respectively p = 0.29 and p = 1.00). The infertile group consisted of women with RIF (n = 11) and infertile women who had not started ART yet, with thus a relatively good prognosis (n = 3) (Table 1). Within the RIF group, primary infertility was the most common type of infertility (81.9%), and the majority was diagnosed with unexplained infertility (45.5%), followed by male factor infertility (36.4%) (Table 1). Among the fertile women, all women had previously conceived within 12 months, and all but one woman had a previous term delivery. The exceptional case involved a woman who had a miscarriage at 10 weeks gestation. One woman (14.3%) had conceived spontaneously and six (85.7%) after donor insemination (Table 2). Table 1 Baseline characteristics of the infertile group

	Infertile (n = 14)	
Total infertile group	14	(100%)	
Median female age (IQR)—yrs	34.5	(31.5–39.0)	
Median female BMI (IQR)—kg/m2	22.3	(20.7–25.8)	
Median duration of infertility (IQR)—mo	31.0	(17.8–57.3)	
Smokers	1	(7.1%)	
Subgroup: infertile without ARTa	3	(21.4%)	
Median female age (IQR)—yrs	35.0	(32.0–35.5)	
Median female BMI (IQR)—kg/m2	25.7	(24.6–25.8)	
Median duration of infertility (IQR)—mo	18.0	(17.5–24.0)	
Smokers	0		
Type of infertility of the femaleb			
Primary	1	(33.3%)	
Secondary	2	(66.7%)	
Causes of infertility			
Unexplained	3	(100%)	
Parity			
Nulliparous	2	(66.7%)	
Primiparous	1	(33.3%)	
Subgroup: RIFc	11	(78.6%)	
Median female age (IQR)—yrs	34.0	(32.0–39.0)	
Median female BMI (IQR)—kg/m2	21.6	(20.6–24.2)	
Median duration of infertility (IQR)—mo	36.0	(18.0–58.0)	
Smokers	1	(9.1%)	
Type of infertility of the femaleb			
Primary	9	(81.9%)	
Secondary	2	(18.2%)	
Causes of infertility			
Unexplained	5	(45.5%)	
Male factor	4	(36.4%)	
Ovulatory disorder	1	(9.1%)	
Mixed	1	(9.1%)	
Parity			
Nulliparous	9	(81.8%)	
Primiparous	2	(18.2%)	
Data are presented as median (IQR) or number (%)

ART, assisted reproductive technology; BMI, body mass index; mo, months; RIF, recurrent implantation failure; yrs, years.

aDefined as women who did not reach a clinical pregnancy after at least 12 months of regular unprotected intercourse and who had not started ART yet.

bPrimary: female has never conceived before. Secondary: female has conceived before.

cDefined as women who have undergone IVF/ICSI with ≥3 failed embryo transfers without clinical pregnancy.

Table 2 Baseline characteristics of the fertile group

	Fertile (n = 7)	
Median female age (IQR)—yrs	35.0	(31.0–38.0)	
Median female BMI (IQR)—kg/m2	24.2	(22.0–30.1)	
Smokers	0	-	
Conceived by	
 Spontaneously	1	(14.3%)	
 Donor insemination	6	(85.7%)	
Parity	
 Nulliparous	1	(14.3%)	
 Primiparous	5	(71.4%)	
 Multiparous	1	(14.3%)	
Data are presented as median (IQR) or number (%)

BMI body mass index, yrs years

Transcriptome profile of infertile vs fertile endometrial organoids

Principal component analysis (PCA) did not show clustering of the samples based on clinical phenotype (i.e. fertile or infertile) (Fig. 2a). A total of 363 significant DEGs were identified, of which 129 were upregulated and 234 were downregulated in the infertile samples compared to fertile samples (Fig. 2b,Supplementary Table 2). Additionally, samples did not cluster based on timing of endometrial biopsy, infertility type or infertility diagnosis (Supplementary Figs. 1a-c).Fig. 2 RNA-seq gene expression analysis of infertile vs fertile endometrial organoids. a Principal component analysis plot of infertile and fertile endometrial organoid samples after batch effect correction. b Volcano plot of differentially expressed genes in infertile organoids compared to fertile organoids. In red are displayed the genes that are significantly differentially expressed with a log2 fold change (FC) ≥ 1.0 or ≤  − 1.0 and an adjusted p-value (padj) of < 0.05; in green the genes with a log2FC ≥ 1.0 or ≤  − 1.0, but with padj ≥ 0.05; in blue the genes with − 1.0 < log2FC < 1.0 and padj < 0.05; and in grey the genes that are not significantly differentially expressed (i.e. − 1.0 < log2FC < 1.0 and padj ≥ 0.05). NS not significant, FC fold change

Overrepresentation analysis of the significantly upregulated DEGs in infertile endometrial organoids compared to fertile endometrial organoids showed progesterone-mediated oocyte maturation, cell cycle, oocyte meiosis, cellular senescence and DNA replication as the significantly annotated pathways (Fig. 3a,Supplementary Table 3a). Genes that were significantly upregulated in multiple of these pathways were CCNA2, CCNB2, CCNB1, CDK1, AURKA, CDC20 and PCNA, which all had a function related to the cell cycle. The most significantly GO Biological Processes (BP) annotated to the significantly upregulated DEGs were mitotic spindle elongation, mitotic spindle midzone assembly and mitotic sister chromatid segregation, with PRC1, BIRC5, CDCA8 and KIF23 as overlapping DEGs, which were all linked to cell cycle processes (Fig. 3b,Supplementary Table 3b). Related significant GO Molecular Functions (MF) and Cellular Components (CC) are displayed in Supplementary Tables3c-d, which were also related to cell division. There were no significant pathways or GO BP and GO CC annotated to the significantly downregulated DEGs. As for GO MF, the significant molecular functions annotated to the significantly downregulated DEGs were androsterone dehydrogenase activity and oxidoreductase activity (Supplementary Table 3e).Fig. 3 Enrichment analysis annotating significantly differentially expressed genes (DEGs) to top 15 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) Biological Processes (BP). Pathways and BP are sorted by adjusted p-value. Terms with red bars have an adjusted p-value of < 0.1 a KEGG pathways annotated to significantly upregulated DEGs. b GO BP annotated to significantly upregulated DEGs. c KEGG pathways annotated to significantly downregulated DEGs. d GO BP annotated to significantly downregulated DEGs

GSEA of GO BP revealed that genes associated with chromosome organisation (NES = 3.81), protein-containing complex organisation (NES = 3.61), cell cycle (NES = 3.55) and mitotic cell cycle process (NES = 3.54) were the most significantly enriched, whereas genes associated with anatomical structure formation involved in morphogenesis (NES =  − 1.32), locomotion (NES =  − 0.73) and cellular response to endogenous stimulus (NES =  − 0.68) were the least enriched in samples of infertile women compared to those of fertile women (Fig. 4a, Supplementary Table 4a), though not significantly. With regards to GO MF, adenyl nucleotide binding (NES = 3.18) and ribonucleotide binding (NES = 2.65) were the most significantly enriched and RNA polymerase II transcription regulatory region sequence-specific DNA binding (NES =  − 0.90) and calcium ion binding (NES =  − 0.75) were the least (not significantly) enriched (Fig. 4b, Supplementary Table 4b). Finally, within the domain GO CC, chromosomal region (NES = 3.51) and chromosome (NES = 3.20) were most significantly enriched, and membrane protein complex (NES =  − 0.63) was not significantly the least enriched (Fig. 4c, Supplementary Table 4c).Fig. 4 Gene set enrichment analysis (GSEA) of Gene Ontology (GO) gene sets. Biological process (a), molecular function (b) and cellular component (c)

Proliferation assay

As cell cycle pathways were enriched in differentiated infertile endometrial organoids compared to fertile endometrial organoids, it was hypothesised that infertile differentiated organoids had higher proliferative activity than fertile differentiated organoids. To validate this hypothesis, the proliferation rate of endometrial organoids from both groups under E2/P4 exposure was assessed by OrganoSeg surface measurements at T1 and T0 (T1/T0 ratios) (Fig. 5a). The mean T1/T0 surface ratios were 1.55 and 1.60 in respectively the infertile and fertile group (Fig. 5b). No statistically significant differences in the T1/T0 ratios between infertile and fertile endometrial organoids (p = 0.69) were observed. In addition, no statistically significant differences were demonstrated in proliferation rate between infertile and fertile endometrial organoids prior to differentiation (T0/T−1 ratios), with a mean surface ratio of 2.37 in infertile organoids and 2.44 in fertile organoids (p = 0.85) (Fig. 5c).Fig. 5 Proliferation assay in infertile vs fertile endometrial organoids. a RNA-seq principal component analysis plot with the subset of organoid lines used for the proliferation assay indicated by the numbers (n = 7 infertile, n = 5 fertile). b, c Mean proliferation rates in infertile and fertile endometrial organoids expressed as T1/T0 (b) and T0/T−1 (c) OrganoSeg surface measurements

Discussion

Main findings

In this study, the transcriptome profiles of E2/P4-treated endometrial organoids of infertile and fertile women were compared. Although the PCA did not result in distinct separation of the samples from both groups, a total of 363 genes were significantly differentially expressed between infertile and fertile organoids. Analysis of significant DEGs by overrepresentation analysis as well as by GSEA showed mainly genes related to cell cycle processes to be highly enriched. Assessment of proliferation by surface measurements, right before and after differentiation with E2/P4, did not show any statistically significant differences between infertile and fertile organoids, and could not confirm our hypothesis of a higher proliferative activity during differentiation in infertile organoids compared to fertile organoids.

Interpretation with regard to previous studies

In accordance with the results of this study, significant upregulation of cell cycle processes in the mid-secretory phase endometrium of infertile RIF patients compared to (fertile) controls was reported in two previous studies in full endometrial biopsies [8, 44]. In concordance with this, increased proliferative activity, determined by Ki67 expression and reduced apoptosis were previously observed in the endometrial glandular epithelial cells of infertile women compared to fertile women [45]. While cell division is an important process during the proliferative phase of the menstrual cycle, proliferative activity is reduced in the in vivo stromal cells of the secretory phase endometrium and completely ceased in the endometrial epithelial cells [46]. In the endometrium of infertile women, the high proliferative activity may be sustained in the endometrial epithelial cells during the secretory phase, leading to inadequate differentiation in preparation to embryo implantation with defective endometrial receptivity, similarly to what has been previously suggested for the eutopic endometrium of women with endometriosis [47–49]. However, rather than upregulation, other studies reported downregulation of the cell cycle pathways in the mid-luteal phase endometrium of infertile RIF patients compared to controls [9, 50], which may suggest that cell cycle activity is altered in various ways in the infertile endometrium. These conflicting findings may also stem from substantial heterogeneity observed across previous endometrial transcriptome studies [51].

No differences in proliferation rate between infertile and fertile endometrial organoids were measured by OrganoSeg. There are several possible explanations for this result. Firstly, it could be possible that the method used was not sensitive enough to detect differences in proliferation. While OrganoSeg surface measurements allow morphometric analysis, this method is not validated yet for the assessment of cell proliferation. Other cell proliferation assays, based on DNA synthesis (e.g. BrdU, EdU assays); metabolic activity of cells (e.g. MTT, XTT assays) or staining of cell proliferation markers (e.g. Ki67) may have yielded different results. Another explanation could be the large heterogeneity within the samples as visualised in the PCA.

This is the first study comparing the transcriptome of infertile to fertile endometrial organoids. Recently, the apical protein secretions of endometrial organoids from infertile and fertile women were compared, showing 150 differentially secreted proteins between both groups, with enrichment in cell membrane movement pathways [52]. In addition, incubation of trophoblast progenitor spheroids with infertile organoid apical secretions significantly impaired the spheroids’ adhesion capacity to organoid epithelial cell monolayers compared to when using the secretions of fertile organoids [52]. These findings show that endometrial organoids are a valuable tool for identifying differences between the infertile and fertile endometrium. However, it is important to bear in mind that interactions with other cell types are also essential to prepare the endometrium for implantation and could modify epithelial cell function [53, 54]. Incorporation of multiple (patient-specific) cell types in the endometrial organoid model may help to mimic the in vivo endometrium but will concomitantly complicate this research system and the interpretation of results. These so-called co-culture models are currently being developed, but still face the challenge of long-term culture [55–58].

Strengths and limitations

The most important strength of this study is the novelty of comparing the transcriptomes of infertile to fertile endometrial organoids of a relatively large number of women. Endometrial organoids serve as a robust in vitro model, enabling the investigation of endometrial epithelial cells ex vivo under controlled conditions that are consistent across all cell lines. Additionally, gene expression profiling was performed by RNA-seq, a very powerful tool for this purpose due to its high reproducibility and low noise level [59].

Limitations of this study include a heterogeneous group of infertile women, including women with a male factor. This may have contributed to the variability in data within the infertile group and the lack of a clear separation of the infertile and fertile samples in our PCA. There was a similar level of variability within the fertile group, suggesting that the current terminologies to define both fertile and infertile groups may not capture underlying biology. Moreover, variability in endometrial biopsy timing may have contributed to the sample heterogeneity. In addition, endometrial organoids consist of different epithelial cell types [19, 22]. Differences may occur in subpopulations of epithelial cells both within and between organoid lines, which may not be discernible when assessing average gene expression through bulk RNA-seq.

Clinical implications and/or implications for future research

Previous studies have reported conflicting results concerning the activity of cell cycle pathways in infertile endometrium as compared to fertile endometrium [8, 9, 44, 50]. Endometrial organoids are a representative in vitro model to study the potential role of the endometrium in infertility in a controlled setting. Because epithelial endometrial organoids consist of multiple different epithelial cell types [19, 22], the use of single-cell RNA-seq will shed further light on potential differences in the function of specific endometrial cells between infertile and fertile organoids. A broad omics integrative analysis of both infertile and fertile (co-culture) endometrial organoids may further contribute to identifying biological mechanisms that underlie implantation failure to provide potential targets for future treatments.

Conclusion

Through transcriptome analysis, enriched cell cycle processes were observed in E2/P4-treated endometrial organoids of infertile women compared to those of fertile women, putatively reflecting higher proliferative activity of differentiated infertile organoids than fertile organoids. Since this could not functionally be confirmed by the OrganoSeg proliferation assay, this hypothesis requires further validation. Future integrative broad omics analysis on the single cell level, including endometrial organoid models, may contribute to improving our understanding of potential differences in the function of infertile and fertile endometrium.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary Material 1.

Supplementary Material 2.

Supplementary Material 3.

Acknowledgements

The authors are grateful to the women who donated tissue for research.

Furthermore, we thank Utrecht Sequencing Facility (USEQ) for providing sequencing service and data. USEQ is subsidised by the University Medical Center Utrecht, Hubrecht Institute, Utrecht University and The Netherlands X-omics Initiative (NWO project 184.034.019).

Author contribution

BB, SM, SF, EN, FB and GS contributed to the study conception and design. BB collected endometrial tissue. BB, FV, MA and MH cultured the endometrial organoids. BB and FV performed the proliferation assays. BB, AIA and FV analysed the data. BB and AIA drafted the manuscript. All authors critically appraised and revised the initial drafts of the manuscript, after which BB adjusted the manuscript for submission. All authors read and approved the submission of the final manuscript.

Funding

The two randomised controlled trials (SCRaTCH and SCRaTCH-OFO), in which endometrial tissue was obtained, were funded by the Netherlands Organisation for Health Research and Development, ‘ZonMw’ (ZonMW project numbers 843002601 and 843001808). The current endometrial organoid transcriptome project was funded by Merck.

Data availability

The RNA-sequencing data that support the findings of this study are openly available on Dataverse NL at 10.34894/IHHCKH.

Declarations

Ethics approval

Ethical approval was obtained from the Institutional Review Board of the University Medical Centre (UMC) Utrecht (15–495 (30 November 2015) and 17–592 (28 September 2017)) and the Biobank Research Ethics Committee of the UMC Utrecht (17–789 (15 March 2018)). Written informed consent was obtained from all study participants.

Competing interests

The SCRATCH and SCRaTCH-OFO trials and the current endometrial organoid transcriptome study received funding from public sources and Merck, but none of these sources had input in the study design and analysis. In addition, several authors report having received personal fees for consultations, conference attendance, and travel as follows: BB from Gedeon Richter Benelux and Guerbet; GS from Guerbet; and SM from Abbott Pharmaceuticals, Ferring, Merck, Oxolife and IBSA. FB is a paid member of the advisory boards of Merck and Ferring and also received speakers’ fees from Besins Healthcare, Cooper Surgical and Merck. All of these personal fees were unrelated to the study reported in this manuscript. The remaining authors declare no conflicts of interest relevant to this study.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Zegers-Hochschild F Adamson GD Dyer S Racowsky C De Mouzon J Sokol R The international glossary on infertility and fertility care, 2017 Hum Reprod 2017 32 1786 1801 10.1093/humrep/dex234 29117321
Zegers-Hochschild F, Adamson GD, Dyer S, Racowsky C, De Mouzon J, Sokol R, et al. The international glossary on infertility and fertility care, 2017. Hum Reprod. 2017;32:1786–801.29117321 10.1093/humrep/dex234
2. Diedrich K Fauser BCJM Devroey P Griesinger G The role of the endometrium and embryo in human implantation Hum Reprod Update 2007 13 365 377 10.1093/humupd/dmm011 17548368
Diedrich K, Fauser BCJM, Devroey P, Griesinger G. The role of the endometrium and embryo in human implantation. Hum Reprod Update. 2007;13:365–77.17548368 10.1093/humupd/dmm011
3. Critchley HOD, Maybin JA, Armstrong GM, Williams ARW. Physiology of the endometrium and regulation of menstruation. Physiol Rev. 2020;100:1149–79. Available from: https://journals.physiology.org/doi/10.1152/physrev.00031.2019
4. Gargett CE, Nguyen HPT, Ye L. Endometrial regeneration and endometrial stem/progenitor cells. Rev Endocr Metab Disord. 2012;13:235–51. Available from: http://link.springer.com/10.1007/s11154-012-9221-9
5. Noyes RW, Hertig AT, Rock J. Dating the endometrial biopsy. Fertil Steril. 1950;1:3–25. Available from: http://www.sciencedirect.com/science/article/pii/S0015028216300620
6. Shi C, Han HJ, Fan LJ, Guan J, Zheng XB, Chen X, et al. Diverse endometrial mRNA signatures during the window of implantation in patients with repeated implantation failure. Hum Fertil. 2018;21:183–94. Available from: 10.1080/14647273.2017.1324180
7. Zhou Q, Yan G, Ding L, Liu J, Yu X, Kong S, et al. EHD1 impairs decidualization by regulating the Wnt4/β-catenin signaling pathway in recurrent implantation failure. EBioMedicine . 2019;50:343–54. Available from: 10.1016/j.ebiom.2019.10.018
8. Bastu E Demiral I Gunel T Ulgen E Gumusoglu E Hosseini MK Potential marker pathways in the endometrium that may cause recurrent implantation failure Reprod Sci 2019 26 879 890 10.1177/1933719118792104 30081718
Bastu E, Demiral I, Gunel T, Ulgen E, Gumusoglu E, Hosseini MK, et al. Potential marker pathways in the endometrium that may cause recurrent implantation failure. Reprod Sci. 2019;26:879–90.30081718 10.1177/1933719118792104
9. Koler M Achache H Tsafrir A Smith Y Revel A Reich R Disrupted gene pattern in patients with repeated in vitro fertilization (IVF) failure Hum Reprod 2009 24 2541 2548 10.1093/humrep/dep193 19542175
Koler M, Achache H, Tsafrir A, Smith Y, Revel A, Reich R. Disrupted gene pattern in patients with repeated in vitro fertilization (IVF) failure. Hum Reprod. 2009;24:2541–8.19542175 10.1093/humrep/dep193
10. Tapia-Pizarro A Figueroa P Brito J Marín JC Munroe DJ Croxatto HB Endometrial gene expression reveals compromised progesterone signaling in women refractory to embryo implantation Reprod Biol Endocrinol 2014 12 1 15 10.1186/1477-7827-12-92 24383664
Tapia-Pizarro A, Figueroa P, Brito J, Marín JC, Munroe DJ, Croxatto HB. Endometrial gene expression reveals compromised progesterone signaling in women refractory to embryo implantation. Reprod Biol Endocrinol. 2014;12:1–15.24383664 10.1186/1477-7827-12-92
11. Cimadomo D, Craciunas L, Vermeulen N, Vomstein K, Toth B. Definition, diagnostic and therapeutic options in recurrent implantation failure: an international survey of clinicians and embryologists. Human Reproduction. 2021;36:305–17. Available from: https://academic.oup.com/humrep/article/36/2/305/6032160
12. Suhorutshenko M Kukushkina V Velthut-Meikas A Altmäe S Peters M Mägi R Endometrial receptivity revisited: endometrial transcriptome adjusted for tissue cellular heterogeneity Hum Reprod 2018 33 2074 2086 10.1093/humrep/dey301 30295736
Suhorutshenko M, Kukushkina V, Velthut-Meikas A, Altmäe S, Peters M, Mägi R, et al. Endometrial receptivity revisited: endometrial transcriptome adjusted for tissue cellular heterogeneity. Hum Reprod. 2018;33:2074–86.30295736 10.1093/humrep/dey301
13. Catalini L, Fedder J. Characteristics of the endometrium in menstruating species: lessons learned from the animal kingdom†. Biol Reprod. 2020;102:1160–9. Available from: https://academic.oup.com/biolreprod/article/102/6/1160/5775593
14. Turco MY, Gardner L, Hughes J, Cindrova-Davies T, Gomez MJ, Farrell L, et al. Long-term, hormone-responsive organoid cultures of human endometrium in a chemically defined medium. Nat Cell Biol. 2017;19:568–77. Available from: http://www.nature.com/articles/ncb3516
15. Boretto M Cox B Noben M Hendriks N Fassbender A Roose H Development of organoids from mouse and human endometrium showing endometrial epithelium physiology and long-term expandability Development (Cambridge) 2017 144 1775 1786 10.1242/dev.148478
Boretto M, Cox B, Noben M, Hendriks N, Fassbender A, Roose H, et al. Development of organoids from mouse and human endometrium showing endometrial epithelium physiology and long-term expandability. Development (Cambridge). 2017;144:1775–86.10.1242/dev.148478
16. Clevers H. Modeling development and disease with organoids. Cell. 2016;165:1586–97. Available from: 10.1016/j.cell.2016.05.082
17. Kretzschmar K Clevers H Organoids: modeling development and the stem cell niche in a dish Dev Cell 2016 38 590 600 10.1016/j.devcel.2016.08.014 27676432
Kretzschmar K, Clevers H. Organoids: modeling development and the stem cell niche in a dish. Dev Cell. 2016;38:590–600.27676432 10.1016/j.devcel.2016.08.014
18. Bui BN, Boretto M, Kobayashi H, van Hoesel M, Steba GS, van Hoogenhuijze N, et al. Organoids can be established reliably from cryopreserved biopsy catheter-derived endometrial tissue of infertile women. Reprod Biomed Online. 2020;00:1–9. Available from: 10.1016/j.rbmo.2020.03.019
19. Fitzgerald HC Dhakal P Behura SK Schust DJ Spencer TE Self-renewing endometrial epithelial organoids of the human uterus Proc Natl Acad Sci U S A 2019 116 23132 23142 10.1073/pnas.1915389116 31666317
Fitzgerald HC, Dhakal P, Behura SK, Schust DJ, Spencer TE. Self-renewing endometrial epithelial organoids of the human uterus. Proc Natl Acad Sci U S A. 2019;116:23132–42.31666317 10.1073/pnas.1915389116
20. Cochrane DR Campbell KR Greening K Ho GC Hopkins J Bui M Single cell transcriptomes of normal endometrial derived organoids uncover novel cell type markers and cryptic differentiation of primary tumours Journal of Pathology 2020 252 201 214 10.1002/path.5511 32686114
Cochrane DR, Campbell KR, Greening K, Ho GC, Hopkins J, Bui M, et al. Single cell transcriptomes of normal endometrial derived organoids uncover novel cell type markers and cryptic differentiation of primary tumours. Journal of Pathology. 2020;252:201–14.32686114 10.1002/path.5511
21. Hewitt SC, Wu S, Wang T, Ray M, Brolinson M, Young SL, et al. The estrogen receptor α cistrome in human endometrium and epithelial organoids. Endocrinology. 2022;163. Available from: https://academic.oup.com/endo/article/doi/10.1210/endocr/bqac116/6650761
22. Garcia-Alonso L Handfield L-F Roberts K Nikolakopoulou K Fernando RC Gardner L Mapping the temporal and spatial dynamics of the human endometrium in vivo and in vitro Nat Genet 2021 53 1698 1711 10.1038/s41588-021-00972-2 34857954
Garcia-Alonso L, Handfield L-F, Roberts K, Nikolakopoulou K, Fernando RC, Gardner L, et al. Mapping the temporal and spatial dynamics of the human endometrium in vivo and in vitro. Nat Genet. 2021;53:1698–711.34857954 10.1038/s41588-021-00972-2
23. van Hoogenhuijze NE, Mol F, Laven JSE, Groenewoud ER, Traas MAF, Janssen CAH, et al. Endometrial scratching in women with one failed IVF/ICSI cycle-outcomes of a randomised controlled trial (SCRaTCH). Hum Reprod. 2020;
24. Bui BN, Torrance HL, Janssen C, Cohlen B, de Bruin JP, den Hartog JE, et al. Does endometrial scratching increase the rate of spontaneous conception in couples with unexplained infertility and a good prognosis (Hunault > 30%)? Study protocol of the SCRaTCH-OFO trial: a randomized controlled trial. BMC Pregnancy Childbirth. 2018;18:511. Available from: https://bmcpregnancychildbirth.biomedcentral.com/articles/10.1186/s12884-018-2160-z
25. Hunault CC Habbema JDF Eijkemans MJC Collins JA Evers JLH te Velde ER Two new prediction rules for spontaneous pregnancy leading to live birth among subfertile couple, based on the synthesis of three previous models Hum Reprod 2004 19 2019 2026 10.1093/humrep/deh365 15192070
Hunault CC, Habbema JDF, Eijkemans MJC, Collins JA, Evers JLH, te Velde ER. Two new prediction rules for spontaneous pregnancy leading to live birth among subfertile couple, based on the synthesis of three previous models. Hum Reprod. 2004;19:2019–26.15192070 10.1093/humrep/deh365
26. American Society for Reproductive Medicine. Revised American Society for Reproductive Medicine classification of endometriosis: 1996. Fertil Steril. 1997;67:817–21. Available from: https://linkinghub.elsevier.com/retrieve/pii/S001502829781391X
27. van Hoogenhuijze NE Torrance HL Mol F Laven JSE Scheenjes E Traas MAF Endometrial scratching in women with implantation failure after a first IVF/ICSI cycle; does it lead to a higher live birth rate? The SCRaTCH study: A randomized controlled trial (NTR 5342) BMC Womens Health 2017 17 1 7 28049464
van Hoogenhuijze NE, Torrance HL, Mol F, Laven JSE, Scheenjes E, Traas MAF, et al. Endometrial scratching in women with implantation failure after a first IVF/ICSI cycle; does it lead to a higher live birth rate? The SCRaTCH study: A randomized controlled trial (NTR 5342). BMC Womens Health. 2017;17:1–7.28049464
28. Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. Available from: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/
29. Kreuger F. TrimGalore [Internet]. 2019. Available from: https://github.com/FelixKrueger/TrimGalore/releases
30. Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. 2012;28:3211–7. Available from: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/bts611
31. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21. Available from: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/bts635
32. Tarasov A, Vilella AJ, Cuppen E, Nijman IJ, Prins P. Sambamba: fast processing of NGS alignment formats. Bioinformatics. 2015;31:2032–4. Available from: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btv098
33. Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012;28:2184–5. Available from: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/bts356
34. Smith AD, Daley T, Deng C. PreSeq [Internet]. 2018. Available from: https://github.com/smithlabcode/preseq
35. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30. Available from: http://www.ncbi.nlm.nih.gov/pubmed/24227677
36. Zhang Y, Parmigiani G, Johnson WE. ComBat-seq: batch effect adjustment for RNA-seq count data. NAR Genom Bioinform [Internet]. 2020;2. Available from: https://academic.oup.com/nargab/article/doi/10.1093/nargab/lqaa078/5909519
37. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol [Internet]. 2014;15:550. Available from: http://genomebiology.biomedcentral.com/articles/10.1186/s13059-014-0550-8
38. Wickham H. ggplot2: elegant graphics for data analysis [Internet]. Springer-Verlag New York; 2016. Available from: https://ggplot2.tidyverse.org
39. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles G v., et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14.
40. Kuleshov M v., Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Enrichr: a comprehensive gene set enrichment analysis web server, update Nucleic Acids Res 2016 2016 44 W90 W97 10.1093/nar/gkw377
Kuleshov M, v., Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Enrichr: a comprehensive gene set enrichment analysis web server, et al. update. Nucleic Acids Res. 2016;2016(44):W90–7.10.1093/nar/gkw377
41. Mootha VK Lindgren CM Eriksson K-F Subramanian A Sihag S Lehar J PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes Nat Genet 2003 34 267 273 10.1038/ng1180 12808457
Mootha VK, Lindgren CM, Eriksson K-F, Subramanian A, Sihag S, Lehar J, et al. PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34:267–73.12808457 10.1038/ng1180
42. Subramanian A Tamayo P Mootha VK Mukherjee S Ebert BL Gillette MA Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles Proc Natl Acad Sci 2005 102 15545 15550 10.1073/pnas.0506580102 16199517
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102:15545–50.16199517 10.1073/pnas.0506580102
43. Borten MA, Bajikar SS, Sasaki N, Clevers H, Janes KA. Automated brightfield morphometry of 3D organoid populations by OrganoSeg. Sci Rep. 2018;8.
44. Devesa-Peiro A, Sebastian-Leon P, Garcia-Garcia F, Arnau V, Aleman A, Pellicer A, et al. Uterine disorders affecting female fertility: what are the molecular functions altered in endometrium? Fertil Steril. Elsevier Inc.; 2020. p. 1261–74.
45. Makker A, Goel MM, Nigam D, Mahdi AA, Das V, Agarwal A, et al. Aberrant Akt activation during implantation window in infertile women with intramural uterine fibroids. Reproductive Sciences [Internet]. 2018;25:1243–53. Available from: http://journals.sagepub.com/doi/10.1177/1933719117737844
46. Makieva S, Giacomini E, Ottolina J, Sanchez AM, Papaleo E, Viganò P. Inside the endometrial cell signaling subway: mind the gap(s). Int J Mol Sci. MDPI AG; 2018.
47. Lessey BA Kim JJ Endometrial receptivity in the eutopic endometrium of women with endometriosis: it is affected, and let me show you why Fertil Steril 2017 108 19 27 10.1016/j.fertnstert.2017.05.031 28602477
Lessey BA, Kim JJ. Endometrial receptivity in the eutopic endometrium of women with endometriosis: it is affected, and let me show you why. Fertil Steril. 2017;108:19–27.28602477 10.1016/j.fertnstert.2017.05.031
48. Park JS Lee JH Kim M Chang HJ Hwang KJ Chang KH Endometrium from women with endometriosis shows increased proliferation activity Fertil Steril 2009 92 1246 1249 10.1016/j.fertnstert.2009.04.025 19501353
Park JS, Lee JH, Kim M, Chang HJ, Hwang KJ, Chang KH. Endometrium from women with endometriosis shows increased proliferation activity. Fertil Steril. 2009;92:1246–9.19501353 10.1016/j.fertnstert.2009.04.025
49. Lessey B Lebovic D Taylor R Eutopic endometrium in women with endometriosis: ground zero for the study of implantation defects Semin Reprod Med 2013 31 109 124 10.1055/s-0032-1333476 23446858
Lessey B, Lebovic D, Taylor R. Eutopic endometrium in women with endometriosis: ground zero for the study of implantation defects. Semin Reprod Med. 2013;31:109–24.23446858 10.1055/s-0032-1333476
50. Koot YEM, van Hooff SR, Boomsma CM, van Leenen D, Koerkamp MJAG, Goddijn M, et al. An endometrial gene expression signature accurately predicts recurrent implantation failure after IVF. Sci Rep [Internet]. 2016;6:1–12. Available from: 10.1038/srep19411
51. Altmäe S Esteban FJ Stavreus-Evers A Simón C Giudice L Lessey BA Guidelines for the design, analysis and interpretation of “omics” data: focus on human endometrium Hum Reprod Update 2014 20 12 28 10.1093/humupd/dmt048 24082038
Altmäe S, Esteban FJ, Stavreus-Evers A, Simón C, Giudice L, Lessey BA, et al. Guidelines for the design, analysis and interpretation of “omics” data: focus on human endometrium. Hum Reprod Update. 2014;20:12–28.24082038 10.1093/humupd/dmt048
52. Zhou W, Barton S, Cui J, Santos LL, Yang G, Stern C, et al. Infertile human endometrial organoid apical protein secretions are dysregulated and impair trophoblast progenitor cell adhesion. Front Endocrinol (Lausanne). 2022;13.
53. Lee JY Lee M Lee SK Role of endometrial immune cells in implantation Clin Exp Reprod Med 2011 38 119 10.5653/cerm.2011.38.3.119 22384430
Lee JY, Lee M, Lee SK. Role of endometrial immune cells in implantation. Clin Exp Reprod Med. 2011;38:119.22384430 10.5653/cerm.2011.38.3.119
54. Evron A Goldman S Shalev E Effect of primary human endometrial stromal cells on epithelial cell receptivity and protein expression is dependent on menstrual cycle stage Hum Reprod 2011 26 176 190 10.1093/humrep/deq296 21098625
Evron A, Goldman S, Shalev E. Effect of primary human endometrial stromal cells on epithelial cell receptivity and protein expression is dependent on menstrual cycle stage. Hum Reprod. 2011;26:176–90.21098625 10.1093/humrep/deq296
55. Abbas Y Brunel LG Hollinshead MS Fernando RC Gardner L Duncan I Generation of a three-dimensional collagen scaffold-based model of the human endometrium Interface Focus 2020 10 20190079 10.1098/rsfs.2019.0079 32194932
Abbas Y, Brunel LG, Hollinshead MS, Fernando RC, Gardner L, Duncan I, et al. Generation of a three-dimensional collagen scaffold-based model of the human endometrium. Interface Focus. 2020;10:20190079.32194932 10.1098/rsfs.2019.0079
56. Heidari-Khoei H, Esfandiari F, Moini A, Yari S, Saber M, Ghaffari Novin M, et al. Derivation of hormone-responsive human endometrial organoids and stromal cells from cryopreserved biopsies. Exp Cell Res [Internet]. 2022;417:113205. Available from: 10.1016/j.yexcr.2022.113205
57. Rawlings TM, Makwana K, Taylor DM, Molè MA, Fishwick KJ, Tryfonos M, et al. Modelling the impact of decidual senescence on embryo implantation in human endometrial assembloids. Elife [Internet]. 2021;10. Available from: https://elifesciences.org/articles/69603
58. Wiwatpanit T, Murphy AR, Lu Z, Urbanek M, Burdette JE, Woodruff TK, et al. Scaffold-free endometrial organoids respond to excess androgens associated with polycystic ovarian syndrome. J Clin Endocrinol Metab [Internet]. 2020;105:769–80. Available from: https://academic.oup.com/jcem/article/105/3/769/5588091
59. Ching T Huang S Garmire LX Power analysis and sample size estimation for RNA-Seq differential expression RNA 2014 20 1684 1696 10.1261/rna.046011.114 25246651
Ching T, Huang S, Garmire LX. Power analysis and sample size estimation for RNA-Seq differential expression. RNA. 2014;20:1684–96.25246651 10.1261/rna.046011.114
