Material and methods
Mice
Seven-to eight -week-old CBA/J female mice were purchased from Janvier Laboratory and housed
under standard conditions (ad libitum food and water, 12h light -dark cycle) in the animal care facility
of Institut Cochin.
Endometriosis model
The endometriosis model consisted of a syngeneic graft of uterine horns to generate endometriosis-
like lesions. In details, donor mice were sacrificed by cervical dislocation , and the uterine horns were
surgically extracted and transferred to a Petri dish containing PBS. The uterine horns were opened
longitudinally with micro scissors and 3 to 5 -mm-length samples were prepared for grafting on the
internal face of the peritoneum of the recipient mice (each horn fragment was weighed). A
preoperative gavage of all donor mice with 100µg/kg/day of 17β-estradiol was performed for two days
before sacrifice (to synchronize estral cycles). Recipient mice were anesthetized using isoflurane. An
incision was made on the ventral midline, and one or two donor horn fragments were sutured onto
the parietal peritoneum with two 7/0 polypropylene stitches (Prolen®, Ethicon, Somerville, NJ , USA).
In all mice, tissue samples were sutured at identical positions in the abdominal wall to ensure that the
host tissue sites exhibited comparable vascularization. The incision was then sutured with a 6/0 nylon
thread.
Immune training and endometriosis
Three independent experiments were performed using 5 –10 mice per group using the previously
described endometriosis model. The control group underwent a sham surgery with a midline incision
and peritoneal stitches. For the endometriosis group, where we evaluated in vivo LPSlow innate immune
training (IMM-EDT), mice were administered daily peritoneal injections of low doses (0.1 mg/kg) of LPS
for 5 days the week before endometriosis induction. The PBS-EDT group received PBS injection instead.
Lesion size was measured in the PBS -EDT and IMM -EDT groups using a high-frequency ultrasound
imaging system (Vevo® 2100 VisualSonics; Toronto, CA) at 1 and 3 weeks after surgery.
Ethical statements
The animal experiment protocol was approved with the ethical approval number DAP20 -104 and
authorization APAFIS#30327 by the Ministry of Higher Education and Scientific Research.
Mice pregnancy follow-up
Three weeks after endometriosis induction, CBA/J female mice from the three groups (Sham, PBS-EDT,
and IMM-EDT) were mated with DBA/2 male mice. The presence of a vaginal plug was considered
embryonic day 0.5 (E0.5) . A high -frequency ultrasound imaging system (Vevo® 2100 VisualSonics;
Toronto, CA) was used to assess the number of implantation sites (between E7.5 and E10.5) and the
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
number of alive/dead/resorbed fetuses (between E11.5 and E14.5). For ultrasound imaging, pregnant
mice were anesthetized with isoflurane and restrained on a heated stage. The fur was removed from
the abdomen using hair removal agents, and pre-warmed ultrasound contact gel was applied to the
shaved abdomen. The beating heart was detected in the living fetuses . Dead fetuses displayed no
beating heart but had visible organ structures. Resorbed embryos/fetuses displayed an echogenic dot
with no discernible organ structure. Gestational age was confirmed by observing key developmental
features (E8.5 heart, head , and whole embryo; E9.5 amniotic membrane, yolk sac , and cerebral
ventricles; E10.5 umbilical cord, placenta, and eyes; E12.5 spine; and E13.5 face, skull bones, and ribs).
At the end of gestation (E18.5), the mice were sacrificed to harvest the fetuses, placentas, uterus, and
lesion tissues . The number of fetuses, their weights, and the weights of the placentas and
endometriosis lesions were recorded. The r esorption rate was calculated as (number of implantation
sites at the first ultrasound - number of live pups at sacrifice) / number of implantation sites at the first
ultrasound × 100. The lesion volume was calculated using measurements obtained by ultrasound
(length × width × width).
Statistical analyses
All data from the mouse experiments were analyzed using GraphPad Prism 8.0 software ( GraphPad,
San Diego, CA, USA). For the comparison of two groups of data, we used the Student t -test with
variables following a normal distribution or the Mann-Whitney U test otherwise. Statistical significance
was set at p < 0.05.
RNA sequencing
Mouse tissues sampling
The three groups of CBA/J mice were obtained again as described previously and mated with DBA mice
three weeks after surgery. Pregnancy was assessed by the presence of a vaginal plug (E0.5 stage), and
ultrasound was performed at E7.5 -E8.5 to confirm pregnancy. Mice were sacrificed at E9.5 , and the
feto-maternal interfaces were retrieved from the uterus and either snap-frozen in nitrogen (for snRNA-
seq) or digested to isolate immune cells (for CITE -seq). Fetal -maternal interfaces showing signs of
abnormal development during ultrasound and/or dissection were not harvested to avoid non-specific
signals (considered as dying structures).
Nuclei isolation and single nucleus RNA sequencing (snRNAseq)
Nuclei were isolated following a published method 22. Briefly, materno-fetal interface tissues were
lysed for 10 min at 37°C with the lysis buffer described in the method. Samples were dounced with 10
strokes and filtered with a 70µm strainer. The filtrate was centrifuged and filtered again using a 40µm
strainer. The filtrate was centrifuged again , and the pellet was resuspended in the staining buffer
described in the method with specific antibodies. Nuclei from each placenta were tagged using
BioLegend ® TotalSeq ™ antibodies. Nuclei were processed using the Single Cell 3′ Gene Expression kit
v3 (10× Chromium, 1000076) according to the manufacturer’s instructions. In brief, 10000 -16,000
nuclei per sample were loaded onto Single Cell Chips B to recover as many nuclei as possible (targeting
3,000–10,000 nuclei per sample) while limiting the potential for doublets. Using a Chromium
Controller, Gel Bead -In Emulsions were generated, and samples were subsequently processed to
isolate and amplify cDNA , and ultimately construct libraries. The quality and concentration of cDNA
were evaluated using an Agilent 2100 Bioanalyzer. The quality and concentration of libraries were
evaluated by qPCR and on an Agilent 2200 TapeStation, and libraries were sequenced on an Illumina
NovaSeq 6000 through the Genom’IC platform of Institut Cochin.
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
Immune cell isolation, surface protein barcoding, single cell RNA and epitope sequencing
(CITE-seq)
Previously snap -frozen E9.5 feto -maternal interfaces were digested using a Multi Tissue Tissue
Dissociation Kit (Miltenyi Biotec 130-110-201) and a gentleMACS Dissociator (Miltenyi Biotec 130-134-
029) at 37°C for 15 min. Samples were filtered using a 100µm strainer and centrifuged to obtain the
pellet. The pellet was resuspended in red blood cell lysis buffer (BioLegend # 420301) for 2 min at room
temperature. Samples were centrifuged, and the pellet was resuspended in RPMI 2% FCS. Samples
were incubated w ith TruStain FcX PLUS Blocking Reagent (BioLegend 156603) for 10 min at 4°C.
Samples were centrifuged, and the pellet was incubated for 30 min at 4°C with a mix of antibodies
containing CD45.2 and Zombie NIR (dead/live marker) for sorting, and several epitope markers listed
below. Each sample was incubated with a specific antibody tagged with a hashtag (TotalSeq B
antimouse), which is listed below. Two wash cycles (centrifugation and addition of RPMI 2% PCF) were
performed for each sample. Samples were sorted using a BD FACSAria ™ III Cell Sorter keeping all the
CD45+ and Zo mbieNIR- cells. The sorted samples were then processed using the Single Cell 3′ Gene
Expression kit v3 (10× Chromium, #1000076) and following the same sequencing protocol as previously
described.
Antibody Barcode Sequence Reference (Cat #)
B0001 (anti-mouse CD4) AACAAGACCCTTGAG 100573
B0002 (anti-mouse CD8a) TACCCGTAATAGCGT 100783
B0182 (anti-mouse CD3) GTATGTCCGCTCGAT 100257
B0184 (anti-mouse CD335 (NKp46)) CCCTTTCACCTCGAA 137641
B0114 (anti-mouse F4/80) TTAACTTCAGCCCGT 123155
B0013 (anti-mouse Ly-6C) AAGTCGTGAGGCATG 128053
B0839 (anti-mouse Ly49H) CCAGTAGGCTTATTA 144721
B0301 Hashtag 1 ACCCACCAGTAAGAC 155831
B0302 Hashtag 2 GGTCGAGAGCATTCA 155833
B0303 Hashtag 3 CTTGCCGCATGTCAT 155835
B0304 Hashtag 4 AAAGCATTCTTCACG 155837
Sequences quality control and alignment
Base calling was completed using Illumina® NovaSeq 6000 RTA v3.4.4 software, and BCL base call files
were converted to FASTQ files using bcl2fastq conversion software (v2.20). Using the Cell Ranger
software suite (v4.0.0) (10X Genomics®), FASTQ files were aligned to the GRCm38 mouse reference
genome (modification steps provided by 10X Genomics), and gene -barcode count matrices were
generated for all samples.
Computational analysis of snRNA-seq and CITE seq data
All basic analyses were performed in R using Seurat 4.1.1 package included in a Snakemake pipeline
with the following rules (with key parameters described below): quality control, demultiplexing,
normalization, integration , and clustering. Muscat package was used to perform differential gene
expression analysis. Scripts are available at the bkheira/doridot-sncell GitHub repository.
Quality control & demultiplexing
All raw gene-barcode count matrices were converted into Seurat objects in R for initial quality control
filtering. Nuclei/cells were filtered using their RNA molecule count (>200 and 500 and < 10000) , and percentage of mitochondrial sequenced genes (<10%). Since samples
were tagged with antibodies presenting a barcode (TotalSeq anti-mouse hashtags), we could separate
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
them in our analysis by reading these barcodes. An assay “HTO” was created in the Seurat Object and
contained barcode data. HTO assay data were normalized and demultiplexed using the “HTOdemux”
function of Seurat. Doublet and negative cells for the barcodes were removed from each Seurat object.
Normalization & integration
Seurat objects were normalized using SCTransform. To identify the cell populations present in all
groups, the Seurat object of each group was integrated with the others. Integration features were
selected using the “SelectIntegrationFeatures” function. The SCTransformed data were prepared using
the “PrepSCTIntegration ” function. Integration anchors were found using the
“FindIntegrationAnchors” function, and Seurat objects were integrated using “IntegrateData”.
Clustering
Clustering was performed using the “FindNeighbors“ and “FindClusters” functions of Seurat with
different clustering resolutions . The optimum resolution was chosen using Clustree package in R.
Uniform Manifold Approximation and Projection (UMAP) was performed using the RunUMAP function.
DEG identification & Gene Set Enrichment Analysis
Differentially expressed genes were identified using Muscat 23 package in R. Muscat is a package in R
that allows pseudobulk analysis by pooling single-cell data in groups of unique clusters per sample. We
followed the authors instructions to perform a differential analysis of the mouse data. Data were then
filtered : we considered differentially expressed a gene with a p -value adjusted to multiple
comparisons smaller than 0.05, an absolute log2FC greater than 0.6 and a minimal expression (logCPM)
of 5.
We used the FGSEA package in R using the tutorial available on the Biostatsquid blog24. GSEA was
performed using the Hallmark database. We added some visualization plots using the package ggplot2.
Illustrations
Figures 1A and 2A were created using BioRender. BOUZID, K. (2025) https://BioRender.com/trkr2lh
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
Results
of three similar experiments. The number of mice is indicated in the figure by n, with a precision of
fetuses (F) and gestant mice (M), if relevant. (B) Number of implantations per mouse (C) Fetal resorption rate
per mouse. For (B) and (C), data are presented as median ± interquartile range. (D) Fetal weight per pup at
E18.5. *: p-value < 0.05 (Mann-Whitney U test).
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
Histological evaluation of E18.5 placentas did not reveal any obvious anomalies associated with
endometriosis (Supplementary Figure 2).
Pre-conceptional introduction of low-dose LPS reduces gestational complications
associated with endometriosis in a mouse model
As we previously showed that immunomodulatory (IMM by repeated exposure to low doses of LPS)
treatment can reduce the lesion size in endometriosis mice 15, we wanted to see if this
immunomodulatory treatment could also reduce the deleterious effect of endometriosis on
pregnancy. We conducted an experiment (Figure 2A) in which endometriosis induction was preceded
by peritoneal injections of repeated low doses of LPS (IMM) or PBS. Three weeks after surgery, the
endometriosis mice were mated with DBA/2 males, and gestation was followed up as previously
described. We previously showed that IMM has no effect on the resorption rate in mice without
endometriosis21.
Figure 2: Immunomodulation (IMM) treatment partially corrects gestational complications in endometriosis
mice. (A) Experimental protocol: repeated injections of PBS or low-dose LPS were administered before inducing
endometriosis by surgery, followed by the same protocol as shown in Figure 1A. (B) Evolution of the lesion
volume as observed on ultrasound. For (C), (D ) and (E), data are presented as mean ± SEM. PBS-EDT: group
receiving PBS and two lesions of endometriosis. IMM-EDT: group receiving LPS at low doses and with two lesions
of endometriosis. (C) Total weight of the two lesions at the time of sacrifice. Student ’s t-test (D) Number of
implantations per mouse, Student ’s t -test (E) Fetal resorption per mouse, Student ’s t -test (F) Correlation
between total lesion weight and resorption rate per mouse, PBS -EDT Pearson r = 0.72 (p = 0.0013), LPS -EDT
Pearson r = -0.21 (p = 0.38). (G) Fetal weights. *: P-value < 0.05, **: P-value < 0.01.
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
The volume of the lesions was measured during ultrasounds following gestation, and IMM treatment
showed a trend toward reducing the volume of the lesions (Figure 2B, p=0.15 by mixed effect analysis).
At the end of gestation, the lesion weight of the IMM -EDT group was significantly lower than that of
the PBS-EDT control group (Figure 2C, p=0.01). The number of implantations was significantly higher
in the IMM-treated group (20% increase, Figure 2D, p = 0.048). The fetal resorption rate was 38% in
the PBS-treated group and 24% in the IMM-treated group (Figure 2E, p = 0.06), which is closer to that
observed in mice without endometriosis (sham group in Figure 1C, 19%). The significant correlation
between lesion size and resorption rate in the endometriosis group (r=0,73, p=0.0013) was lost with
the IMM (Figure 2F). Fetal weights at E18.5 were equivalent in all groups (Figure 2G).
SnRNAseq reveals an increase of Gata4 and a decrease of P rap1 expression in decidual
cells in endometriotic mice
Since there are a decreased mean number of implantation sites and an increased resorption rate in
the endometriosis group, we hypothesized that there is a possible dysfunction of the decidua and/or
the placenta that leads to an altered or interrupted embryo development. To validate this hypothesis,
feto-maternal interfaces at E9.5 without visible signs of ongoing resorption were collected from mice
with sham surgery or with two lesions of endometriosis, and their transcriptomes were sequenced at
the single-nuclei level.
We identified 13 cell clusters in our dataset that were distributed homogeneously among the samples
(Figure 3A, 3B). Wt1, a known marker of decidualized cells25, was found in seven clusters (Figure 3D).
Among them, there are two decidual fibroblast/stromal clusters, identified by the expression of
collagen genes ( Col3a1) and several genes involved in the remodelling of the extracellular matrix of
the decidua during implantation and placentation (Sulf1, Adamtsl1)26,27 ,28. Interestingly, these decidual
fibroblasts expressed Il15 and its receptor Il15ra, which are involved in the activation of uterine Natural
Killer (uNK) 29,30. The other Wt1+ clusters were identified as decidual stromal cells (DSC) as they
expressed several known markers of decidual stromal cells at different levels ( Hand2, Pgr, and Esr1),
suggesting that the clusters may represent different stages of differentiation or different subtypes of
decidual stromal cells. All clusters’ markers are available in Supplementary Table 1.
The expression of the maternal gene Xist in placental samples associated with a male embryo
segregated the cells of embryonic origin from those of maternal origin (Figure 3C). We confirmed their
trophoblast origin by the expression of specific genes that mark them as trophoblast progenitors, such
as Epcam and Ror2 expression in embryonic cells 31, 32 and trophoblast cells with Mct1/Sl16a1,
CD147/Bsg, Glut1/Slc2a133, and Nr6a134 (Figure 3F). The expression of Pecam1 was used to identify
endothelial cells within the placenta 35,36. Cells of hematopoietic origin were identified by Ptprc (CD45
gene, Figure 3E). Pericytes have been identified using the Pdgfrb gene.
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
Through differential expression analysis, we identified several differentially expressed genes (DEG) in
all clusters (Supplementary Table 2). In the different clusters of decidual stromal cells, we observed
two genes that were consistently differentially expressed: Gata4 was upregulated and Prap1 was
downregulated in the group with endometriosis -like lesions. Interestingly, these 2 genes differential
expression were normalized by the immunomodulatory treatment, with a decrease of Gata4 and an
increase in Prap1 in the IMM-EDT group compared to the PBS-EDT group.
Figure 3: Single -nucleus RNA sequencing of E9.5 fetomaternal interface reveals its composition and
characteristics (A) UMAP of cell transcriptomes (B) Graph representing the proportion of a cell type (cluster) in
each sample (top) and the number of cells retained for analyses after quality control filtering (bottom) (C)
Expression of Xist gene in materno-fetal interfaces with male embryos (D) Expression of Wt1 (E) Expression of
Ptprc (CD45 gene) (F) Dotplot of cluster identifying genes. DSC: Decidual stromal cells
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
For the next analysis, we grouped the clusters by big cell types (decidual stromal cells, trophoblast
cells, endothelial cells, immune cells, and pericytes) and compared the groups of mice (sham, PBS-EDT,
IMM-EDT) by performing differential expression analysis and gene set enrichment analysis (GSEA)
(Figure 4, Supplementary Table 3 ). When plotting a heatmap of the top (up to 10) DEG in each cell
Figure 4: Differential expression analysis reveals alterations driven by endometriosis and partially corrected
by the IMM treatment
(A) Number of differentially expressed genes per group of cells (B) Heatmap of the top 10 differentially
expressed genes between the sham and endometriosis groups. Cell clusters were grouped by cell type as
shown in Figure 3B. (C) Volcano plot showing the differentially expressed genes within the decidual stromal
cells (D) Dotplot of Gene Set Enrichment Analysis: top 10 significantly enriched pathways from Sham vs PBS -
EDT analysis and how they are in the PBS -EDT vs IMM -EDT in the pooled DSC group. Pathways involved in
proliferation are surrounded by red rectangles. * Significant enrichment after adjustment.
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
group, we observed that the transcriptomic expression of the IMM-EDT group was often between the
expression of the control and endometriosis groups, indicating a partial correction of the gene
expression differences induced by endometriosis by IMM treatment. We observed transcriptomic
changes induced by endometriosis in all the cell types ( Figure 4A, B). We then chose to focus on
decidual stromal cells, which are the group of cells with the most differentially expressed genes
induced by the presence of endo metriotic lesions ( Figure 4A ). The differential expression analysis
showed an overexpression of Gata4 (log2FC of 2.27 and adjusted p -value of 1.7 × 10−5) and a
downregulation of Prap1 (log2FC of -2.37 and adjusted p-value of 1.63 × 10−5 ) in the PBS-EDT group
compared to the Sham group (Figure 4C). IMM led to the correction of these expressions (Figure 4B).
GSEA showed that pathways related to proliferation (Myc targets, E2F targets, G2M checkpoints, and
KRAS signalling) were negatively enriched in the endometriosis group. This downregulation was
significantly corrected by IMM treatment, except for the G2M checkpoint (Figure 4D).
CITEseq analysis reveals transcriptomic differences in all immune cell types and an
inflammation in endometriotic placentas
After immunomodulation, it is relevant to examine the consequences o n immune cells directly. As
these cells account for a small proportion of the maternal fetal interface and thus are not well
represented in the whole interface snRNAseq dataset, we performed Cellular Indexing of
Transcriptomes and Epitopes by Sequencing (CITEse q or proteotranscriptomics ) on isolated immune
cells (CD45+) in the E9.5 materno-fetal interface.
We identified 11 clusters of cells in this dataset (Figure 5A ). Cluster markers are available in
Supplementary Table 5. All the cell types were equally represented in the samples (Figure 5B). A major
group of myeloid cells was recognized by Cd68 and Itgam (CD11b) expression ( Figure 5G ) and the
expression of the monocyte marker Ly6C ( Figure 5D). Within this group, we identified dendritic cells
expressing Cd209a and a small cluster of plasmacytoid dendritic cells expressing Siglech. Neutrophils
were identified by the expression of calprotein genes ( S100a9 and S100a8). Macrophages, identified
here with the antibody F4/80 (Figure 5C), had several marker genes, such as Wwp1, Ccl8, Stab1, Gas6,
and Fcrls. Three other clusters were identified as monocytes (all expressing Ly6C and/or Itgam), and
specific genes differentiating these subcategories of monocytes were identified.
The other major cell type was natural killer (NK) cells and innate lymphoid cells (ILC) identified by the
different granzyme genes (Gzmb, Gzmc, Gzme) and perforin (Prf1) and with the surface protein data
of the CITEseq, using the NKp46 antibody ( Figure 5E). T cells were identified by the surface proteins
CD3 (figure 5F), CD4, and CD8a, and at the transcriptomic level, they had specific expressions of Itk and
Il7r. B cells were identified with immunoglobulin genes such as Igkc or Ighm, as well as genes like Ebf1,
Bank1, Cd79a, and Cd79b.
Finally, in the CD45+ cells, we also had a small cluster of endothelial cells recognized by the expression
of known genes such as Flt1 and Egfl7.
We could see in the PBS-EDT group compared to the Sham group few transcriptomic differences in all
the cell types ( Figure 6A). The IMM treatment induced a higher number of significant transcriptomic
differences in all cell types compared to the Sham group ( Supplementary Figure 3A) and to the PBS -
EDT group ( Figure 6A ). Indeed, immunomodulation of the immune system leads to transcriptional
changes37. A gene that was consistently downregulated among the cell groups (except macrophages)
was Slc15a2 (Figure 6B), which encodes the peptide transporter PEPT2 (see Supplementary Table 7).
To determine the pathways in which these genes are involved, we performed GSEA. Notably, we
observed an upregulation of TNF α signalling in macrophages in endometriosis, which tended to be
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
corrected by IMM treatment ( Figure 6C). This upregulation is also observed in endometriosis in the
pooled group of monocytes (Supplementary figure 3), in T and B cells, but only corrected by the IMM
in B cells (see Supplementary Table 8). In NK and ILC, endometriosis led to a negative enrichment of
IFNγ and IFN α response pathways, and IMM tended to normalize this pathway ( Figure 6D ). This
negative enrichment of IFN γ response in endometriosis is also found in decidual stromal cells,
endothelial cells, and trophoblasts.
Figure 5: Composition of E9.5 fetomaternal interface by single cell proteotranscriptomics sequencing of
immune cells
(A) UMAP of cell transcriptomes (B) Graph representing the proportion of a cell type (cluster) in each sample
(top) and the number of cells retained for analyses after quality control filtering (bottom) (C) Featureplot of
the surface antibody F4/80 staining macrophages (D) Ly-6C (E) Nkp46 (F) CD3 (G) Dotplot of cluster identifying
genes and antibodies
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
Discussion
In this study, we first examined the impact of endometriosis-like lesions on gestation in CBA/JxDBA/2N
crossing, with fewer implantations and a higher proportion of fetal resorptions than in the sham group.
Other studies have shown the consequences of endometriosis on gestation in mouse models, including
a reduction in the pregnancy rate 9,10,12 (gestant mice per group), increase in the resorption rate 11,12,
Figure 6: Differential expression analysis in immune cells reveals alterations driven by endometriosis and
partially corrected by the IMM treatment
(A) Number of differentially expressed genes in each cell group. (B) Heatmap of the top 3 differentially
expressed genes between the sham and endometriosis groups. Cell clusters were grouped by cell types as
shown in Figure 5A. (C) Dotplot of Gene Set Enrichment Analysis : top (up to) 10 significantly enriched
pathways from Sham vs PBS-EDT analysis and how they are in the PBS-EDT vs IMM-EDT in Macrophages and
(D) in NK and ILC. * Significant enrichment after adjustment.
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
decrease in litter size12, and reduced birth weight for the pups11. In our study, we found no significant
difference in fetal weight, which may be due to the fact that mice were sacrificed before parturition;
however, other studies also showed no difference in pup weight at E18.510. Currently, most of these
studies on mouse models have been used as preclinical models to test the effects of potentially
beneficial molecules, but the underlying mechanisms of these complications are yet to be understood.
Our study suggests that the presence of lesions and/or inflammation induced by endometriosis leads
to a deleterious environment for pregnancy development.
This is further substantiated by the improvement in pregnancy outcomes with IMM treatment. In
accordance with previous findings of our team, inducing immunotolerance reduced the evolution and
final size of endometriosis lesions and had a beneficial effect on pregnancy complications, with an
increase in the number of implantations and a decrease in the fetal resorption rate. In the IMM
condition, the correlation between lesion size and fetal resorption was lost, suggesting that the
beneficial effect of IMM o n pregnancy may be due to a reduction in inflammation rather than a
reduction in lesion size.
Due to the increased risk of miscarriage in endometriosis, there may be defects in early pregnancy. To
study this, we used our mouse model again and retrieved the feto -maternal interface at E9.5, a stage
at which the structure of the placenta is developing. Transcriptomic sequencing has shown that
endometriosis induces changes in all cell types of the developing placenta, even in trophoblasts, which
are exclusively fetal. These fetal cells are derived from oocytes exposed to chronic inflammation
induced by endometriosis, which may explain the consequences on trophoblastic cells. These changes
in fetal cells of the placenta raise questions about the transmission of endometriosis consequences to
the fetus via epigenetic signals. Indeed, endometriosis has a genetic heritability of 50%, and until now,
several GWAS have only been able to explain 5% 38 of the genetics of endometriosis. In addition to
potential rare variants not detected by GWAS, missing heritability could also be due to epigenetics and
should be studied in further research. Our model may allow for the exploration of its potential impact
on the next generation.
The major cell type at E9.5 in the developing placenta is decidual stromal cells (DSC), which are derived
from the decidualization of endometrial stromal cells. We observed most transcriptomic changes in
the presence of endometriosis in this group of cells. Gata4 gene expression was upregulated by
endometriosis and normalized by IMM treatment. This upregulation of Gata4 has already been found
in endometriosis within ectopic and eutopic endometrium in patients 39, making our mouse model
relevant for studying this upregulation. The involvement of Gata4 in gestation has not been studied,
but its family members GATA2 and GATA3 are important during preimplantation and early post -
implantation development, with knockouts of both Gata2 and Gata3 leading to placental defects40 in
trophoblasts.
DSC also showed downregulation of Prap1 in endometriosis, which was corrected by IMM treatment.
In mice, Prap1 appears to be involved in establishing uterine receptivity to the embryo 41 and is an
indicator of successful implantation 42. Therefore, its downregulation may participate in the reduced
number of implantations and increased resorption observed in endometriotic mice. Prap1 is a gene
positively regulated by the progesterone receptor PGR in the mouse endometrium 43. Pgr expression
in our dataset was non-significantly downregulated by endometriosis (logFC of -0.40, p -value of
0.0044, and adjusted p-value of 0.077) and may participate in the downregulation of Prap1 expression.
This downregulation of Pgr, although not significant in our dataset when we adjusted the p -value for
multiple testing, should be further studied, and our model may be interesting to study the
progesterone resistance found in endometriosis in human44.
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
In our GSEA of mouse DSCs at E9.5, most of the downregulated pathways in the endometriosis group
were related to proliferation. A defect in proliferation during early gestation could explain miscarriages
or further placental pathologies found to be increased in endometriosis.
We identified two clusters of decidual stromal cells that expressed extracellular matrix remodelling
genes and named them decidual fibroblasts. These decidual fibroblasts expressed Il15, which encodes
IL-15, an important factor for the differentiation and proliferation of uNK cells, allowing uNK cells to
support embryo implantation, spiral artery remodelling, and immune tolerance during pregnancy 45.
Until now, it was known that IL -15 signalling was coming from decidual stromal cells 30, but here we
show that only a specific subtype of DSC has this role. Interestingly, these cells also express Il15ra, a
gene encoding a receptor for IL-15, suggesting the possibility of an autocrine signal within these cells.
When sorting the immune cells from the E9.5 feto -maternal interfaces, we identified a majority of
myeloid cells and up to 20 -25% of lymphoid cells , concordant with previous littereature 46. This is
different from human early placentas, where NK cells represent 70% of immune cells 47 and
macrophages represent 20 and 30%48. Previous studies have shown that less than 1% of immune cells
are of fetal origin at E10.5 46, a later stage than that in our study; therefore, it is likely that the vast
majority of cells in our dataset are of maternal origin.
Endometriosis induced transcriptomic changes in all immune cells. When focusing on macrophages in
endometriosis, GSEA showed an upregulation of an inflammatory pathway, TNF α via NF κB. In
endometriosis, it has been shown that macrophages are more activated and participate in
inflammation within the peritoneal cavity 49. Our results suggest that macrophages are inflammatory,
even at the maternal-fetal interface. We can imagine that inflammation at an early stage of gestation
could lead to a poor establishment of the maternal -fetal interface and thus induce
miscarriages/resorptions. The upregulation of inflammatory pathways was no longer observed when
endometriosis mice were previously trained with low doses of LPS (IMM). We have already shown that
peritoneal macrophages in a mouse model of endometriosis become more immuno tolerant with LPS
training at low doses15, and we show here that this immunotolerance is also induced in the maternal -
fetal interface.
The other major immune cells during pregnancy are NK cells. In these cells in endometriosis mice, the
IFNγ and IFNα response pathways were downregulated, whereas the Myc target V1 pathway, which is
involved in proliferation, was upregulated. uNK cells are the main source of IFN γ in the murine
materno-fetal interface50. IFNγ is essential for the development of gestation in mice and participates
in NK maturation51. Thus, NK cells, which may be less able to respond to this signal, may not mature
well and may disrupt the development of pregnancy. Mice that are IFNγ-/- present a high rate of fetal
resorptions, no spiral artery remodelling, and decidua necrosis. They also have more uNK at the MFI,
as if they are trying to compensate for the lack of IFN γ by having more producing cells 50. The
upregulation of a proliferation pathway in NK cells in our dataset may be a consequence of the
downregulated IFNγ response. IFN-α and IFN-β are regulators of NK cell activation and their production
of IFNγ52, and in mice that are Ifna-/-, the same phenotype as Ifng-/- is observed at the materno-fetal
interface, suggesting that NK cell IFN-α response is also essential for the development of the MFI. This
downregulated response to IFN γ was also found in the decidual cells, endothelial cells, and
trophoblasts of our dataset, and only the downregulation of IFN α was found in decidual cells and
trophoblasts. Given the importance of IFN γ signaling in vascular remodelling and maintenance of the
decidua, this general d ownregulation of the response may be a cause of disruption of gestational
development.
Of note, B cells are among the immune cells with most differentially expressed genes in endometriosis
compared to sham group. GSEA showed enrichment of several pathways involved in inflammation and
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
immune response. B cells involvement in endometriosis is still incompletely defined and controversial
with evidence supporting increased B cell activation and auto -antibody production 53. Our mouse
model could be used to further study B cells in endometriosis.
Our results show that a murine model is pertinent for studying gestational complications induced by
endometriosis and reveals mechanisms affecting decidual stromal and immune cells at the beginning
of gestation. Our findings suggest that targeting the immune system may be a therapeutic strategy to
improve pregnancy outcomes and reduce the risk of miscarriage.
Acknowledgments and funding:
We would like to thank the different platforms at the Institut Cochin and Université Paris Cité that
helped us obtain our results: GENOM’IC, CYBIO, HIST’IM, BIOINFORMAT’IC, and IPOP-UP.
We would like to thank Université Paris Cité for the IDEX funding of the PlacentAtlas project (ANR-18-
IDEX-0001).
K. B. received a PhD fellowship from BioSPC doctoral School at Université Paris Cité. C. M. is supported
by ANR-20-CE14-0004. L.D. is funded by the European Union (European Research Council (ERC) starting
grant No 101078556). Views and opinions expressed are however those of the author(s) only and do
not necessarily reflect those of the European Union. Neither the European Union nor the granting
authority can be held responsible for them.
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
1. Zhang, T., De Carolis, C., Man, G. C. W. & Wang, C. C. The link between immunity, autoimmunity
and endometriosis: a literature update. Autoimmun. Rev. 17, 945–955 (2018).
2. Endometriosis-associated infertility: aspects of pathophysiological mechanisms and treatment
options - Tanbo - 2017 - Acta Obstetricia et Gynecologica Scandinavica - Wiley Online Library.
https://obgyn-onlinelibrary-wiley-com.ezproxy.u-paris.fr/doi/10.1111/aogs.13082.
3. Lalani, S. et al. Endometriosis and adverse maternal, fetal and neonatal outcomes, a systematic
review and meta-analysis. Hum. Reprod. 33, 1854–1865 (2018).
4. Brosens, I. Endometriosis and the outcome of in vitro fertilization. Fertil. Steril. 81, 1198–1200
(2004).
5. Huang, Y. et al. Miscarriage on Endometriosis and Adenomyosis in Women by Assisted
Reproductive Technology or with Spontaneous Conception: A Systematic Review and Meta -
Analysis. BioMed Res. Int. 2020, 4381346 (2020).
6. Borghese, B. et al. Low Birth Weight Is Strongly Associated with the Risk of Deep Infiltrating
Endometriosis: Results of a 743 Case-Control Study. PLOS ONE 10, e0117387 (2015).
7. Jenabi, E. & Fereidooni, B. The association between endometriosis and the risk of LBW: A meta -
analysis. J. Endometr. Pelvic Pain Disord. 10, 72–78 (2018).
8. Gao, M., Allebeck, P., Mishra, G. D. & Koupil, I. Developmental origins of endometriosis: a Swedish
cohort study. J Epidemiol Community Health 73, 353–359 (2019).
9. Park, Y. et al. Oleuropein suppresses endometriosis progression and improves the fertility of mice
with endometriosis. J. Biomed. Sci. 29, 1–20 (2022).
10. Bilotas, M. A. et al. Interplay between Endometriosis and Pregnancy in a Mouse Model. PLOS ONE
10, e0124900 (2015).
11. França, P. R. de C. et al. R-954, a bradykinin B1 receptor antagonist, as a potential therapy in a
preclinical endometriosis model. Peptides 181, 171294 (2024).
12. Lai, Z.-Z. et al. Protopanaxadiol improves endometriosis associated infertility and miscarriage in
sex hormones receptors -dependent and independent manners. Int. J. Biol. Sci. 17, 1878 –1894
(2021).
13. Elsherbini, M. et al. Impact of Chronic Exposure to Endometriosis on Perinatal Outcomes:
Establishment of a Mouse Model. Biomedicines 10, 2627 (2022).
14. Soni, U. K. et al. MCP-1 promotes ILK phosphorylation at Ser246 during endometriosis
development and affects the pregnancy outcome. Mol. Hum. Reprod. 31, gaaf004 (2025).
15. Jeljeli, M. et al. Macrophage Immune Memory Controls Endometriosis in Mice and Humans. Cell
Rep. 33, 108325 (2020).
16. Ander, S. E., Diamond, M. S. & Coyne, C. B. Immune responses at the maternal-fetal interface. Sci.
Immunol. 4, eaat6114 (2019).
17. Izumi, G. et al. Involvement of immune cells in the pathogenesis of endometriosis. J. Obstet.
Gynaecol. Res. 44, 191–198 (2018).
18. Chuang, P. -C., Wu, M. -H., Shoji, Y. & Tsai, S. -J. Downregulation of CD36 results in reduced
phagocytic ability of peritoneal macrophages of women with endometriosis. J. Pathol. 219, 232–
241 (2009).
19. Jeljeli, M. et al. Trained immunity modulates inflammation -induced fibrosis. Nat. Commun. 10,
5670 (2019).
20. Mujal, A. M., Delconte, R. B. & Sun, J. C. Natural Killer Cells: From Innate to Adaptive Features.
Annu. Rev. Immunol. 39, 417–447 (2021).
21. Dang, Y. et al. BCG-trained innate immunity leads to fetal growth restriction by altering immune
cell profile in the mouse developing placenta. J. Leukoc. Biol. 111, 1009–1020 (2022).
22. Gaublomme, J. T. et al. Nuclei multiplexing with barcoded antibodies for single-nucleus genomics.
Nat. Commun. 10, 2907 (2019).
23. Crowell, H. L. et al. muscat detects subpopulation -specific state transitions from multi -sample
multi-condition single-cell transcriptomics data. Nat. Commun. 11, 6077 (2020).
24. Twomey, L. Easy Gene Set Enrichment Analysis in R with fgsea() - biostatsquid.com.
https://biostatsquid.com/fgsea-tutorial-gsea/ (2023).
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
25. Zhou, J., Rauscher, F. J. & Bondy, C. Wilms’ tumor (WT1) gene expression in rat decidual
differentiation. Differentiation 54, 109–114 (1993).
26. Lum, D. H., Tan, J., Rosen, S. D. & Werb, Z. Gene Trap Disruption of the Mouse Heparan Sulfate 6-
O-Endosulfatase Gene, Sulf2. Mol. Cell. Biol. 27, 678–688 (2007).
27. Kirn-Safran, C. B., D’Souza, S. S. & Carson, D. D. Heparan sulfate proteoglycans and their binding
proteins in embryo implantation and placentation. Semin. Cell Dev. Biol. 19, 187–193 (2008).
28. Hirohata, S. et al. Punctin, a Novel ADAMTS-like Molecule, ADAMTSL-1, in Extracellular Matrix*. J.
Biol. Chem. 277, 12182–12189 (2002).
29. Brighton, P. J. et al. Clearance of senescent decidual cells by uterine natural killer cells in cycling
human endometrium. eLife 6, e31274 (2017).
30. Sharma, S., Godbole, G. & Modi, D. Decidual Control of Trophoblast Invasion. Am. J. Reprod.
Immunol. 75, 341–350 (2016).
31. Ueno, M. et al. c-Met-Dependent Multipotent Labyrinth Trophoblast Progenitors Establish
Placental Exchange Interface. Dev. Cell 27, 373–386 (2013).
32. Marsh, B. & Blelloch, R. Single nuclei RNA-seq of mouse placental labyrinth development. eLife 9,
e60266 (2020).
33. Nagai, A., Takebe, K., Nio -Kobayashi, J., Takahashi -Iwanaga, H. & Iwanaga, T. Cellular Expression
of the Monocarboxylate Transporter (MCT) Family in the Placenta of Mice. Placenta 31, 126–133
(2010).
34. Mehta, D. V., Kim, Y.-S., Dixon, D. & Jetten, A. M. Characterization of the Expression of the Retinoid-
related, Testis-associated Receptor (RTR) in Trophoblasts. Placenta 23, 281–287 (2002).
35. van Mourik, J. A., Leeksma, O. C., Reinders, J. H., de Groot, P. G. & Zandbergen -Spaargaren, J.
Vascular endothelial cells synthesize a plasma membrane protein indistinguishable from the
platelet membrane glycoprotein IIa. J. Biol. Chem. 260, 11300–11306 (1985).
36. Jakeman, L. B., Winer, J., Bennett, G. L., Altar, C. A. & Ferrara, N. Binding sites for vascular
endothelial growth factor are localized on endothelial cells in adult rat tissues. J. Clin. Invest. 89,
244–253 (1992).
37. Netea, M. G. et al. Trained immunity: A program of innate immune memory in health and disease.
Science 352, aaf1098 (2016).
38. Rahmioglu, N. et al. The genetic basis of endometriosis and comorbidity with other pain and
inflammatory conditions. Nat. Genet. 55, 423–436 (2023).
39. Fouquet, B., Santulli, P., Noel, J.-C. & Misrahi, M. Ovarian-like differentiation in eutopic and ectopic
endometrioses with aberrant FSH receptor, INSL3 and GATA4/6 expression. BBA Clin. 6, 143–152
(2016).
40. Home, P. et al. Genetic redundancy of GATA factors in the extraembryonic trophoblast lineage
ensures the progression of preimplantation and postimplantation mammalian development.
Development 144, 876–888 (2017).
41. Xiong, G. -F. et al. Estradiol-regulated proline -rich acid protein 1 is repressed by class I histone
deacetylase and functions in peri -implantation mouse uterus. Mol. Cell. Endocrinol. 331, 23–33
(2011).
42. Diao, H., Xiao, S., Zhao, F. & Ye, X. Uterine luminal epithelium –specific proline-rich acidic protein
1 (PRAP1) as a marker for successful embryo implantation. Fertil. Steril. 94, 2808-2811.e1 (2010).
43. Kim, T. H. & Jeong, J.-W. Proline-Rich Acidic Protein 1 (PRAP1) is a Target of ARID1A and PGR in the
Murine Uterus. Dev. Reprod. 23, 277–284 (2019).
44. Zhang, P. & Wang, G. Progesterone Resistance in Endometriosis: Current Evidence and Putative
Mechanisms. Int. J. Mol. Sci. 24, 6992 (2023).
45. Mori, M., Bogdan, A., Balassa, T., Csabai, T. & Szekeres-Bartho, J. The decidua—the maternal bed
embracing the embryo—maintains the pregnancy. Semin. Immunopathol. 38, 635–649 (2016).
46. Moore, A. R. et al. Gestationally dependent immune organization at the maternal-fetal interface.
Cell Rep. 41, 111651 (2022).
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint
47. King, A., Balendran, N., Wooding, P., Carter, N. P. & Loke, Y. W. CD3 - leukocytes present in the
human uterus during early placentation: phenotypic and morphologic characterization of the
CD56++ population. Dev. Immunol. 1, 169–190 (1991).
48. Lessin, D. L., Hunt, J. S., King, C. R. & Wood, G. W. Antigen Expression by Cells Near the Maternal -
Fetal Interface. Am. J. Reprod. Immunol. Microbiol. 16, 1–7 (1988).
49. Riccio, L. da G. C. et al. Immunology of endometriosis. Best Pract. Res. Clin. Obstet. Gynaecol. 50,
39–49 (2018).
50. Ashkar, Ali. A. & Croy, B. A. Interferon-γ Contributes to the Normalcy of Murine Pregnancy1. Biol.
Reprod. 61, 493–502 (1999).
51. Ashkar, A. A., Di Santo, J. P. & Croy, B. A. Interferon γ Contributes to Initiation of Uterine Vascular
Modification, Decidual Integrity, and Uterine Natural Killer Cell Maturation during Normal Murine
Pregnancy. J. Exp. Med. 192, 259–270 (2000).
52. Matikainen, S. et al. IFN-α and IL-18 synergistically enhance IFN -γ production in human NK cells:
differential regulation of Stat4 activation and IFN -γ gene expression by IFN -α and IL -12. Eur. J.
Immunol. 31, 2236–2245 (2001).
53. Riccio, L. G. C., Baracat, E. C., Chapron, C., Batteux, F. & Abrão, M. S. The role of the B lymphocytes
in endometriosis: A systematic review. J. Reprod. Immunol. 123, 29–34 (2017).
.CC-BY-NC-ND 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted January 2, 2026. ; https://doi.org/10.64898/2026.01.02.695166doi: bioRxiv preprint