The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Tianjin Central Hospital of Obstetrics and Gynecology (Approval No. CDPK991890). The requirement for written informed consent was formally waived due to the retrospective observational design. All retrieved clinical data were strictly anonymized and de-identified prior to analysis to ensure patient confidentiality and data security. The research tools used in the protocol are listed in the Table of Materials.
1. Patient selection and cohort definition
Electronic medical records were reviewed to identify patients who underwent fertility-sparing surgical treatment for pathologically confirmed endometriosis at Tianjin Central Hospital of Obstetrics and Gynecology between January 2020 and December 2024, thereby establishing a single-center retrospective cohort. Eligible patients were women of reproductive age with intraoperative and pathological confirmation of endometriosis who underwent fertility-preserving surgery with retention of the uterus and at least one ovary.
Patients were required to have complete baseline clinical characteristics, disease severity metrics including revised American Society for Reproductive Medicine (rASRM) staging, lesion phenotype classification, operative details, and at least 12 months of continuous postoperative follow-up. Cases were classified according to the dominant disease phenotype, including isolated ovarian endometrioma, deep infiltrating endometriosis (DIE) with or without endometrioma, adenomyosis-associated disease, and superficial peritoneal or adhesion-dominant disease. These classifications were used to support phenotype-aware interpretation of pooled outcome estimates.
Patients without definitive pathological confirmation of endometriosis, those with concurrent systemic malignancies, those who underwent hysterectomy or bilateral oophorectomy, and those lost to follow-up or lacking primary fertility outcome data were excluded. All reporting was conducted in accordance with the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines. The patient screening workflow and final cohort assembly process are shown in Figure 1.

Figure 1: Patient selection flowchart and study design. A total of 312 patients with surgically suspected endometriosis were screened between January 2020 and December 2024. After application of exclusion criteria, 219 patients with pathologically confirmed endometriosis were included in the final analysis and classified into laparoscopic surgery (LS; n = 150) and laparotomic surgery (LT; n = 69) cohorts. The figure also summarizes phenotype classification, collected clinical variables, postoperative management variables, and fertility and recurrence outcomes evaluated during follow-up. Please click here to view a larger version of this figure.
2. Surgical intervention
- Laparoscopic approach
Patients undergoing laparoscopic surgery received general anesthesia and were positioned in the modified lithotomy position. Pneumoperitoneum was established using either a Veress needle or an open Hasson technique at the umbilicus, with intra-abdominal pressure maintained at 12–14 mmHg. A primary 10 mm trocar was inserted for the laparoscope, and three accessory 5 mm trocars were placed in the lower abdomen under direct visualization to facilitate pelvic access and operative triangulation.
The pelvic cavity was systematically inspected to determine the extent of endometriotic disease and to assign the rASRM stage. Adhesiolysis was performed using cold scissors or ultrasonic energy devices to mobilize adherent structures and restore pelvic anatomy. The extent of lesion excision was determined by intraoperative assessment of tissue planes, with priority given to preserving critical neurovascular structures in anatomically distorted regions.
Ovarian endometriomas were managed using a cystectomy stripping technique. The cleavage plane between the cyst capsule and the normal ovarian cortex was identified whenever feasible, and hemostasis was achieved with selective bipolar coagulation, minimizing thermal injury to preserve ovarian tissue and ovarian reserve.
Deep-infiltrating endometriosis (DIE) lesions were resected according to intraoperative anatomical findings and surgical feasibility. When deep nodular involvement was identified, ureterolysis and pararectal space dissection were performed as appropriate. Multidisciplinary surgical collaboration was undertaken when bowel or urinary tract resection was required.
- Laparotomic approach
Patients undergoing laparotomic surgery received general anesthesia and underwent either a Pfannenstiel or midline vertical abdominal incision to access the pelvic cavity. A self-retaining retractor was used to optimize visualization of pelvic structures.
Complex pelvic adhesions, particularly in cases with complete cul-de-sac obliteration, were identified and carefully dissected using a combination of sharp and blunt techniques to expose critical neurovascular and ureteral landmarks. Extensive endometriotic lesions and large ovarian cysts were resected using conventional surgical instrumentation. Hemostasis was achieved primarily through meticulous suture ligation, with minimal thermal coagulation. Following completion of the procedure, the pelvic cavity was irrigated with sterile saline, hemostasis was confirmed, and the abdominal wall was closed in anatomical layers.
3. Postoperative management and fertility guidance
Postoperative care was individualized according to clinical recovery and disease characteristics. Analgesia was administered as indicated, and recovery parameters, including time to first flatus and length of hospital stay, were recorded throughout the postoperative period.
Postoperative hormonal suppression therapy was initiated within four weeks after surgery according to lesion severity, residual disease status, reproductive goals, and clinician judgment. Gonadotropin-releasing hormone agonists or oral progestins were administered for three to six months in selected patients with advanced disease phenotypes to suppress residual disease activity and reduce postoperative inflammatory burden.
Comprehensive fertility counseling was provided on the basis of postoperative pelvic anatomy, ovarian reserve, patient age, fertility timeline, and individual reproductive preferences. Natural conception was encouraged when clinically appropriate. Patients with diminished ovarian reserve, advanced reproductive age, severe disease, or prolonged infertility were counseled regarding timely referral for assisted reproductive technology (ART) rather than assuming surgery alone would represent the preferred fertility strategy.
4. Outcome evaluation and follow-up
Patients underwent structured postoperative follow-up assessments at 3, 6, and 12 months after surgery. Reproductive outcomes were documented systematically, including clinical pregnancy (CP) and time to pregnancy (TTP).
Clinical pregnancy was defined by ultrasonographic visualization of an intrauterine gestational sac. Time to pregnancy was calculated from the initiation of conception attempts to the first recognized clinical pregnancy.
Disease recurrence was monitored using structured symptom assessment and pelvic ultrasonography. Suspected recurrence events were independently reviewed by two experienced gynecologists. Imaging recurrence was defined as the detection of newly developed cystic lesions that met the established morphological criteria for ovarian endometrioma. Operational definitions for fertility outcomes and recurrence endpoints are summarized in Table 1.
Table 1: Definitions of fertility outcomes, recurrence endpoints, and follow-up schedule. Clinical definitions, diagnostic criteria, treatment variables, and scheduled follow-up time points applied throughout the analysis are summarized. Abbreviations: ART = assisted reproductive technology; TTP = time to pregnancy; VAS = visual analog scale; MRI = magnetic resonance imaging; GnRH-a = gonadotropin-releasing hormone agonist; COC = combined oral contraceptive; RFS = recurrence-free survival; rASRM = revised American Society for Reproductive Medicine. Please click here to download this Table.
5. Data extraction and statistical analysis
Demographic, surgical, and therapeutic variables were extracted from the institutional hospital database. Key analytical covariates included age, body mass index (BMI), duration of infertility, and preoperative cancer antigen 125 (CA125) levels.
Continuous variables were summarized as means with standard deviations (SDs) or medians with interquartile ranges (IQRs), depending on data distribution. Categorical variables were reported as frequencies and percentages. Missing data were addressed using multiple imputation when appropriate, and sensitivity analyses were conducted to assess the potential impact of incomplete records on the study findings.
Baseline imbalance between surgical groups was quantified using standardized mean differences (SMDs), with values < 0.10 indicating acceptable balance. Multivariable regression served as the primary adjusted analytical approach for the aggregate cohort. Propensity-score matching or inverse probability of treatment weighting (IPTW) was specified only as a sensitivity analysis and was to be reported only if patient-level covariate overlap was adequate, including matched or weighted sample size and post-adjustment SMDs. Because the present submission reports the aggregate multivariable model rather than a matched or weighted cohort, propensity-based findings should not be interpreted unless those diagnostics are supplied.
Adjusted reproductive outcomes were evaluated using multivariable logistic regression models. Time-to-event outcomes, including recurrence-free survival (RFS), were assessed using Kaplan–Meier survival analysis and Cox proportional hazards regression.
All statistical tests were two-sided, and P < 0.05 was considered statistically significant. Subgroup analyses were interpreted as exploratory and hypothesis-generating due to limited subgroup sample sizes and the potential for multiple comparisons. The primary recurrence-free survival analysis was restricted to the prespecified 12-month postoperative period. Any analyses extending beyond 12 months were considered exploratory.
The reported aggregate regression models and survival analyses were performed using R software (version 4.3.0; R Foundation for Statistical Computing). Analytical workflows were verified through standardized code execution to support methodological reproducibility.