Ablation of hematopoietic stem cell derived adipocytes reduces tumor burden in syngeneic mouse models of high-grade serous carcinoma.

In this study we examined the influence of hematopoietic stem cell-derived adipocytes (HSCDAs) on the proliferation and metastasis of high-grade serous carcinoma (HGSC) - the most common type of ovarian cancer. HSCDAs are a subtype of adipocytes that differentiate from myeloid precursors that traffic from bone marrow to adipose tissue and accumulate therein. These are distinct from conventional mesenchymal adipocytes (CMAs), which are derived from mesenchymal precursors. We hypothesized that HSCDAs promote HGSC progression and establish a pro-tumoral niche within peritoneal adipose tissues such as the omentum. Primary human white adipose tissue samples were obtained via biopsy and then sorted into myeloid and mesenchymal populations through flow cytometry. These adipose precursors were then differentiated in vitro into mature HSCDAs and CMAs, respectively. Transcriptomic analysis showed that HSCDAs have a distinct transcriptional profile from CMAs, including downregulation of cell cycle and upregulation of multiple metabolic and adipogenic pathways. Using ELISA, we found that HSCDAs secreted greater amounts of inflammatory cytokines IL-6 and IL-8 than CMAs. Next, we incubated HGSC cells in conditioned media from HSCDAs and CMAs and performed proliferation and protein expression profiling. HGSC cells in HSCDA media, compared to those in CMA media, had elevated expression of protein markers related to epithelial to mesenchymal plasticity, including fibronectin, as well as increased serine phosphorylation of pro-survival AKT1/2. Conversely, HGSC cells in HSCDA media exhibited comparably downregulated expression of tumor suppressors including the Wnt regulator GSK3ß. Depending on the cell line and adipose donor, HGSC cells also showed altered growth rates in conditioned media. We next investigated the role of HSCDAs in HGSC progression and metastasis in vivo . We generated immunocompetent mice that were either HSCDA Proficient (can make both adipocyte subtypes) or Deficient (can only make CMAs). Using these models, we conducted two independent tumor studies using the ID8 ( Tp53-/- , Brca2-/- ) and SO ( Tp53-/- , Brca1/2 wild-type, Hras and Myc amplified) syngeneic models. Overall tumor burden was lower in HSCDA Deficient mice in both models. In the ID8 model, omental tumors from HSCDA Deficient mice showed reduced proliferation (Ki67) and apoptosis (cleaved caspase 3) relative to those from Proficient mice. Transcriptionally, omental ID8 tumors from HSCDA Deficient downregulated oxidative phosphorylation, adipogenesis, and fatty acid metabolism relative to tumors from HSCDA Proficient mice. These pathways were enriched in HSCDA cells in vitro , suggesting that ablation of HSCDAs had a significant influence on the tumor metabolic environment. Reduced inflammatory pathways in ID8 tumors from HSCDA Deficient mice were also observed leading us to interrogate immune cell infiltration into omental tumors. Compared to HSCDA Proficient mice, tumors from HSCDA Deficient mice showed reduced densities of dendritic cells (DC) and natural killer (NK) cells, as well as fewer DCs, NKs, and B-cells in proximity to tumor cells, as determined by spatial analysis. Overall, our data suggest that HSCDAs promote HGSC survival and plasticity while downregulating expression of tumor suppressors and altering the peritoneal immune and metabolic environment to promote HGSC progression.

Chemotherapy drives tertiary lymphoid structures that correlate with ICI-responsive TCF1+CD8+ T cells in metastatic ovarian cancer.

PURPOSE: Patients with high-grade serous ovarian carcinoma (HGSOC) are virtually insensitive to immune checkpoint inhibitors (ICIs) employed as standalone therapeutics, at least in part reflecting microenvironmental immunosuppression. Thus, conventional chemotherapeutics and targeted anticancer agents that not only mediate cytotoxic effects but also promote the recruitment of immune effector cells to the HGSOC microenvironment stand out as promising combinatorial partners for ICIs in this oncological indication. EXPERIMENTAL DESIGN: We harnessed a variety of transcriptomic, spatial and functional assays to characterize the differential impact of neo-adjuvant paclitaxel-carboplatin on the immunological configuration of paired primary and metastatic HGSOC biopsies as compared to NACT-naïve HGSOC samples from 5 independent patient cohorts. RESULTS: We found neo-adjuvant chemotherapy (NACT)-driven endoplasmic reticulum stress and calreticulin exposure in metastatic HGSOC lesions culminates with the establishment of a dense immune infiltrate including follicular T cells (TFH cells), a prerequisite for mature tertiary lymphoid structure (TLS) formation. In this context, TLS maturation was associated with an increased intratumoral density of ICI-sensitive TCF1+PD-1+ CD8+ T cells over their ICI-insensitive TIM-3+PD-1+ counterparts. Consistent with this notion, chemotherapy coupled with a PD-1-targeting ICI provided a significant survival benefit over either therapeutic approach in syngeneic models of HGSOC bearing high (but not low) tumor mutational burden. CONCLUSION: Altogether, our findings suggest that NACT promotes TLS formation and maturation in HGSOC lesions, de facto preserving an intratumoral ICI-sensitive T-cell phenotype. These observations emphasize the role of rational design, especially relative to the administration schedule, for clinical trials testing chemotherapy plus ICIs in patients with HGSOC.

Spatiotemporal architecture of immune cells and cancer-associated fibroblasts in high-grade serous ovarian carcinoma.

High-grade serous ovarian carcinoma (HGSOC), the deadliest form of ovarian cancer, is typically diagnosed after it has metastasized and often relapses after standard-of-care platinum-based chemotherapy, likely due to advanced tumor stage, heterogeneity, and immune evasion and tumor-promoting signaling from the tumor microenvironment. To understand how spatial heterogeneity contributes to HGSOC progression and early relapse, we profiled an HGSOC tissue microarray of patient-matched longitudinal samples from 42 patients. We found spatial patterns associated with early relapse, including changes in T cell localization, malformed tertiary lymphoid structure (TLS)-like aggregates, and increased podoplanin-positive cancer-associated fibroblasts (CAFs). Using spatial features to compartmentalize the tissue, we found that plasma cells distribute in two different compartments associated with TLS-like aggregates and CAFs, and these distinct microenvironments may account for the conflicting reports about the role of plasma cells in HGSOC prognosis.

Metabolic targeting of cancer associated fibroblasts overcomes T-cell exclusion and chemoresistance in soft-tissue sarcomas.

T cell-based immunotherapies have exhibited promising outcomes in tumor control; however, their efficacy is limited in immune-excluded tumors. Cancer-associated fibroblasts (CAFs) play a pivotal role in shaping the tumor microenvironment and modulating immune infiltration. Despite the identification of distinct CAF subtypes using single-cell RNA-sequencing (scRNA-seq), their functional impact on hindering T-cell infiltration remains unclear, particularly in soft-tissue sarcomas (STS) characterized by low response rates to T cell-based therapies. In this study, we characterize the STS microenvironment using murine models (in female mice) with distinct immune composition by scRNA-seq, and identify a subset of CAFs we termed glycolytic cancer-associated fibroblasts (glyCAF). GlyCAF rely on GLUT1-dependent expression of CXCL16 to impede cytotoxic T-cell infiltration into the tumor parenchyma. Targeting glycolysis decreases T-cell restrictive glyCAF accumulation at the tumor margin, thereby enhancing T-cell infiltration and augmenting the efficacy of chemotherapy. These findings highlight avenues for combinatorial therapeutic interventions in sarcomas and possibly other solid tumors. Further investigations and clinical trials are needed to validate these potential strategies and translate them into clinical practice.

Tertiary lymphoid structures and B cells determine clinically relevant T cell phenotypes in ovarian cancer.

Intratumoral tertiary lymphoid structures (TLSs) have been associated with improved outcome in various cohorts of patients with cancer, reflecting their contribution to the development of tumor-targeting immunity. Here, we demonstrate that high-grade serous ovarian carcinoma (HGSOC) contains distinct immune aggregates with varying degrees of organization and maturation. Specifically, mature TLSs (mTLS) as forming only in 16% of HGSOCs with relatively elevated tumor mutational burden (TMB) are associated with an increased intratumoral density of CD8+ effector T (TEFF) cells and TIM3+PD1+, hence poorly immune checkpoint inhibitor (ICI)-sensitive, CD8+ T cells. Conversely, CD8+ T cells from immunologically hot tumors like non-small cell lung carcinoma (NSCLC) are enriched in ICI-responsive TCF1+ PD1+ T cells. Spatial B-cell profiling identifies patterns of in situ maturation and differentiation associated with mTLSs. Moreover, B-cell depletion promotes signs of a dysfunctional CD8+ T cell compartment among tumor-infiltrating lymphocytes from freshly isolated HGSOC and NSCLC biopsies. Taken together, our data demonstrate that - at odds with NSCLC - HGSOC is associated with a low density of follicular helper T cells and thus develops a limited number of mTLS that might be insufficient to preserve a ICI-sensitive TCF1+PD1+ CD8+ T cell phenotype. These findings point to key quantitative and qualitative differences between mTLSs in ICI-responsive vs ICI-irresponsive neoplasms that may guide the development of alternative immunotherapies for patients with HGSOC.

INHBA(+) cancer-associated fibroblasts generate an immunosuppressive tumor microenvironment in ovarian cancer.

Effective targeting of cancer-associated fibroblasts (CAFs) is hindered by the lack of specific biomarkers and a poor understanding of the mechanisms by which different populations of CAFs contribute to cancer progression. While the role of TGFß in CAFs is well-studied, less attention has been focused on a structurally and functionally similar protein, Activin A (encoded by INHBA). Here, we identified INHBA(+) CAFs as key players in tumor promotion and immunosuppression. Spatiotemporal analyses of patient-matched primary, metastatic, and recurrent ovarian carcinomas revealed that aggressive metastatic tumors enriched in INHBA(+) CAFs were also enriched in regulatory T cells (Tregs). In ovarian cancer mouse models, intraperitoneal injection of the Activin A neutralizing antibody attenuated tumor progression and infiltration with pro-tumorigenic subsets of myofibroblasts and macrophages. Downregulation of INHBA in human ovarian CAFs inhibited pro-tumorigenic CAF functions. Co-culture of human ovarian CAFs and T cells revealed the dependence of Treg differentiation on direct contact with INHBA(+) CAFs. Mechanistically, INHBA/recombinant Activin A in CAFs induced the autocrine expression of PD-L1 through SMAD2-dependent signaling, which promoted Treg differentiation. Collectively, our study identified an INHBA(+) subset of immunomodulatory pro-tumoral CAFs as a potential therapeutic target in advanced ovarian cancers which typically show a poor response to immunotherapy.

Tumor-Associated Macrophages Expand Chemoresistant, Ovarian Cancer Stem-Like Cells.

The persistence of ovarian cancer stem-like cells (OvCSCs) after chemotherapy resistance has been implicated in relapse. However, the ability of these relatively quiescent cells to produce the robust tumor regrowth necessary for relapse remains an enigma. Since normal stem cells exist in a niche, and tumor-associated macrophages (TAMs) are the highest abundance immune cell within ovarian tumors, we hypothesized that TAMs may influence OvCSC proliferation. To test this, we optimized OvCSC enrichment by sphere culture and in vitro polarization of monocytes to a TAM-like M2 phenotype. Using cocultures that permitted the exchange of only soluble factors, we found that M2 macrophages increased the proliferation of sphere cells. Longer-term exposure (5-7 days) to soluble TAM factors led to retention of some stem cell features by OvCSCs but loss of others, suggesting that TAMs may support an intermediate stemness phenotype in OvCSCs. Although TAM coculture decreased the percentage of OvCSCs surviving chemotherapy, it increased the overall number. We therefore sought to determine the influence of this interaction on chemotherapy efficacy in vivo and found that inhibiting macrophages improved chemotherapy response. Comparing the gene expression changes in OvCSCs cocultured with TAMs to publicly available patient data identified 34 genes upregulated in OvCSCs by exposure to soluble TAM factors whose expression correlates with outcome. Overall, these data suggest that TAMs may influence OvCSC proliferation and impact therapeutic response.

Metabolic dependencies and targets in ovarian cancer.

Reprogramming of cellular metabolism is a hallmark of cancer. Cancer cells undergo metabolic adaptations to maintain tumorigenicity and survive under the attack of immune cells and chemotherapy in the tumor microenvironment. Metabolic alterations in ovarian cancer in part overlap with findings from other solid tumors and in part reflect unique traits. Altered metabolic pathways not only facilitate ovarian cancer cells' survival and proliferation but also endow them to metastasize, acquire resistance to chemotherapy, maintain cancer stem cell phenotype and escape the effects of anti-tumor immune defense. In this review, we comprehensively review the metabolic signatures of ovarian cancer and their impact on cancer initiation, progression, and resistance to treatment. We highlight novel therapeutic strategies targeting metabolic pathways under development.

Ciliated Cells in Ovarian Cancer Decrease with Increasing Tumor Grade and Disease Progression.

Ciliated cell markers expressed in epithelial ovarian cancers (EOC) are associated with improved survival. We examined the distribution of cells expressing ciliated cell markers in various EOC histologies and stages. Immunohistochemistry and/or multiplex immunofluorescence were used to determine the expression of FOXJ1 and/or CAPS (ciliated cell markers) in tissue microarrays including 4 normal fallopian tubes, 6 normal endometria, 16 cystadenomas, 25 borderline tumors, 21 low-grade carcinomas, and 118 high-grade carcinomas (HGSOC) (46 serous, 29 endometrioid, 30 clear cell, 13 mucinous). CAPS+ cells were observed in normal fallopian tubes and endometria and in ~85% of serous benign and borderline tumors and low-grade carcinomas but only in <40% of HGSOC. mRNA data from an independent cohort showed higher FOXJ1 and CAPS expression in serous borderline tumors and low-grade carcinomas compared to HGSOC. In HGSOC, ciliated cell-positive markers were observed in 52% primary tumors compared to 26% of patient-matched synchronous metastases, and 24% metachronous metastases (p = 0.009). mRNA data from an independent HGSOC cohort showed lower levels of CAPS in metastases than in primary tumors (p = 0.03). Overall, the study revealed that ciliated cells were less common in mucinous EOC, the percentage of ciliated cell marker-positive cases decreased with increasing grade, and the percentage of ciliated cells decreased in HGSOC metastases compared to patient-matched primary tumors.

Dynamic Changes in the Extracellular Matrix in Primary, Metastatic, and Recurrent Ovarian Cancers.

Cancer-associated fibroblasts (CAFs) and their extracellular matrix are active participants in cancer progression. While it is known that functionally different subpopulations of CAFs co-exist in ovarian cancer, it is unclear whether certain CAF subsets are enriched during metastatic progression and/or chemotherapy. Using computational image analyses of patient-matched primary high-grade serous ovarian carcinomas, synchronous pre-chemotherapy metastases, and metachronous post-chemotherapy metastases from 42 patients, we documented the dynamic spatiotemporal changes in the extracellular matrix, fibroblasts, epithelial cells, immune cells, and CAF subsets expressing different extracellular matrix components. Among the different CAF subsets, COL11A1+ CAFs were associated with linearized collagen fibers and exhibited the greatest enrichment in pre- and post-chemotherapy metastases compared to matched primary tumors. Although pre- and post-chemotherapy metastases were associated with increased CD8+ T cell infiltration, the infiltrate was not always evenly distributed between the stroma and cancer cells, leading to an increased frequency of the immune-excluded phenotype where the majority of CD8+ T cells are present in the tumor stroma but absent from the tumor parenchyma. Overall, most of the differences in the tumor microenvironment were observed between primary tumors and metastases, while fewer differences were observed between pre- and post-treatment metastases. These data suggest that the tumor microenvironment is largely determined by the primary vs. metastatic location of the tumor while chemotherapy does not have a significant impact on the host microenvironment.

FOXC2 Promotes Vasculogenic Mimicry in Ovarian Cancer.

FOXC2 is a forkhead family transcription factor that plays a critical role in specifying mesenchymal cell fate during embryogenesis. FOXC2 expression is associated with increased metastasis and poor survival in various solid malignancies. Using in vitro and in vivo assays in mouse ovarian cancer cell lines, we confirmed the previously reported mechanisms by which FOXC2 could promote cancer growth, metastasis, and drug resistance, including epithelial-mesenchymal transition, stem cell-like differentiation, and resistance to anoikis. In addition, we showed that FOXC2 expression is associated with vasculogenic mimicry in mouse and human ovarian cancers. FOXC2 overexpression increased the ability of human ovarian cancer cells to form vascular-like structures in vitro, while inhibition of FOXC2 had the opposite effect. Thus, we present a novel mechanism by which FOXC2 might contribute to cancer aggressiveness and poor patient survival.

Ovarian cancer cell fate regulation by the dynamics between saturated and unsaturated fatty acids.

Fatty acids are an important source of energy and a key component of phospholipids in membranes and organelles. Saturated fatty acids (SFAs) are converted into unsaturated fatty acids (UFAs) by stearoyl Co-A desaturase (SCD), an enzyme active in cancer. Here, we studied how the dynamics between SFAs and UFAs regulated by SCD impacts ovarian cancer cell survival and tumor progression. SCD depletion or inhibition caused lower levels of UFAs vs. SFAs and altered fatty acyl chain plasticity, as demonstrated by lipidomics and stimulated Raman scattering (SRS) microscopy. Further, increased levels of SFAs resulting from SCD knockdown triggered endoplasmic reticulum (ER) stress response with brisk activation of IRE1α/XBP1 and PERK/eIF2α/ATF4 axes. Disorganized ER membrane was visualized by electron microscopy and SRS imaging in ovarian cancer cells in which SCD was knocked down. The induction of long-term mild ER stress or short-time severe ER stress by the increased levels of SFAs and loss of UFAs led to cell death. However, ER stress and apoptosis could be readily rescued by supplementation with UFAs and reequilibration of SFA/UFA levels. The effects of SCD knockdown or inhibition observed in vitro translated into suppression of intraperitoneal tumor growth in ovarian cancer xenograft models. Furthermore, a combined intervention using an SCD inhibitor and an SFA-enriched diet initiated ER stress in tumors growing in vivo and potently blocked their dissemination. In all, our data support SCD as a key regulator of the cancer cell fate under metabolic stress and point to treatment strategies targeting the lipid balance.

Computational pathology in ovarian cancer.

Histopathologic evaluations of tissue sections are key to diagnosing and managing ovarian cancer. Pathologists empirically assess and integrate visual information, such as cellular density, nuclear atypia, mitotic figures, architectural growth patterns, and higher-order patterns, to determine the tumor type and grade, which guides oncologists in selecting appropriate treatment options. Latent data embedded in pathology slides can be extracted using computational imaging. Computers can analyze digital slide images to simultaneously quantify thousands of features, some of which are visible with a manual microscope, such as nuclear size and shape, while others, such as entropy, eccentricity, and fractal dimensions, are quantitatively beyond the grasp of the human mind. Applications of artificial intelligence and machine learning tools to interpret digital image data provide new opportunities to explore and quantify the spatial organization of tissues, cells, and subcellular structures. In comparison to genomic, epigenomic, transcriptomic, and proteomic patterns, morphologic and spatial patterns are expected to be more informative as quantitative biomarkers of complex and dynamic tumor biology. As computational pathology is not limited to visual data, nuanced subvisual alterations that occur in the seemingly "normal" pre-cancer microenvironment could facilitate research in early cancer detection and prevention. Currently, efforts to maximize the utility of computational pathology are focused on integrating image data with other -omics platforms that lack spatial information, thereby providing a new way to relate the molecular, spatial, and microenvironmental characteristics of cancer. Despite a dire need for improvements in ovarian cancer prevention, early detection, and treatment, the ovarian cancer field has lagged behind other cancers in the application of computational pathology. The intent of this review is to encourage ovarian cancer research teams to apply existing and/or develop additional tools in computational pathology for ovarian cancer and actively contribute to advancing this important field.

Focal Serous Tubal Intra-Epithelial Carcinoma Lesions Are Associated With Global Changes in the Fallopian Tube Epithelia and Stroma.

Objective: Serous tubal intra-epithelial carcinoma (STIC) lesions are thought to be precursors to high-grade serous ovarian cancer (HGSOC), but HGSOC is not always accompanied by STIC. Our study was designed to determine if there are global visual and subvisual microenvironmental differences between fallopian tubes with and without STIC lesions. Methods: Computational image analyses were used to identify potential morphometric and topologic differences in stromal and epithelial cells in samples from three age-matched groups of fallopian tubes. The Benign group comprised normal fallopian tubes from women with benign conditions while the STIC and NoSTIC groups consisted of fallopian tubes from women with HGSOC, with and without STIC lesions, respectively. For the morphometric feature extraction and analysis of the stromal architecture, the image tiles in the STIC group were further divided into the stroma away from the STIC (AwaySTIC) and the stroma near the STIC (NearSTIC). QuPath software was used to identify and quantitate secretory and ciliated epithelial cells. A secretory cell expansion (SCE) or a ciliated cell expansion (CCE) was defined as a monolayered contiguous run of >10 secretory or ciliated cells uninterrupted by the other cell type. Results: Image analyses of the tubal stroma revealed gradual architectural differences from the Benign to NoSTIC to AwaySTIC to NearSTIC groups. In the epithelial topology analysis, the relative number of SCE and the average number of cells within SCE were higher in the STIC group than in the Benign and NoSTIC groups. In addition, aging was associated with an increased relative number of SCE and a decreased relative number of CCE. ROC analysis determined that an average of 15 cells within SCE was the optimal cutoff value indicating the presence of a STIC lesion in the tubal epithelium. Conclusions: Our findings suggest that global stromal alterations and age-associated reorganization of tubal secretory and ciliated cells are associated with STIC lesions. Further studies will need to determine if these alterations precede STIC lesions and provide permissible conditions for the formation of STIC.

Integrated single-cell transcriptome analysis reveals heterogeneity of esophageal squamous cell carcinoma microenvironment.

The tumor microenvironment is a highly complex ecosystem of diverse cell types, which shape cancer biology and impact the responsiveness to therapy. Here, we analyze the microenvironment of esophageal squamous cell carcinoma (ESCC) using single-cell transcriptome sequencing in 62,161 cells from blood, adjacent nonmalignant and matched tumor samples from 11 ESCC patients. We uncover heterogeneity in most cell types of the ESCC stroma, particularly in the fibroblast and immune cell compartments. We identify a tumor-specific subset of CST1+ myofibroblasts with prognostic values and potential biological significance. CST1+ myofibroblasts are also highly tumor-specific in other cancer types. Additionally, a subset of antigen-presenting fibroblasts is revealed and validated. Analyses of myeloid and T lymphoid lineages highlight the immunosuppressive nature of the ESCC microenvironment, and identify cancer-specific expression of immune checkpoint inhibitors. This work establishes a rich resource of stromal cell types of the ESCC microenvironment for further understanding of ESCC biology.

Human iPSC-derived fallopian tube organoids with BRCA1 mutation recapitulate early-stage carcinogenesis.

Germline pathogenic mutations in BReast CAncer (BRCA1) genes are thought to drive normal fallopian tube epithelial (FTE) cell transformation to high-grade serous ovarian cancer. No human models capture the sequence of events for disease initiation and progression. Here, we generate induced pluripotent stem cells (iPSCs) from healthy individuals and young ovarian cancer patients with germline pathogenic BRCA1 mutations (BRCA1mut). Following differentiation into FTE organoids, BRCA1mut lines exhibit cellular abnormalities consistent with neoplastic transformation compared to controls. BRCA1mut organoids show an increased production of cancer-specific proteins and survival following transplantation into mice. Organoids from women with the most aggressive ovarian cancer show the greatest pathology, indicating the potential value to predict clinical severity prior to disease onset. These human FTE organoids from BRCA1mut carriers provide a faithful physiological in vitro model of FTE lesion generation and early carcinogenesis. This platform can be used for personalized mechanistic and drug screening studies.

Immunological configuration of ovarian carcinoma: features and impact on disease outcome.

Epithelial ovarian carcinoma (EOC) is a relatively rare malignancy but is the fifth-leading cause of cancer-related death in women, largely reflecting early, prediagnosis dissemination of malignant disease to the peritoneum. At odds with other neoplasms, EOC is virtually insensitive to immune checkpoint inhibitors, correlating with a tumor microenvironment that exhibits poor infiltration by immune cells and active immunosuppression. Here, we comparatively summarize the humoral and cellular features of primary and metastatic EOC, comparatively analyze their impact on disease outcome, and propose measures to alter them in support of treatment sensitivity and superior patient survival.

Identification of Infertility-Associated Topologically Important Genes Using Weighted Co-expression Network Analysis.

Endometriosis has been associated with a high risk of infertility. However, the underlying molecular mechanism of infertility in endometriosis remains poorly understood. In our study, we aimed to discover topologically important genes related to infertility in endometriosis, based on the structure network mining. We used microarray data from the Gene Expression Omnibus (GEO) database to construct a weighted gene co-expression network for fertile and infertile women with endometriosis and to identify gene modules highly correlated with clinical features of infertility in endometriosis. Additionally, the protein-protein interaction network analysis was used to identify the potential 20 hub messenger RNAs (mRNAs) while the network topological analysis was used to identify nine candidate long non-coding RNAs (lncRNAs). Functional annotations of clinically significant modules and lncRNAs revealed that hub genes might be involved in infertility in endometriosis by regulating G protein-coupled receptor signaling (GPCR) activity. Gene Set Enrichment Analysis showed that the phospholipase C-activating GPCR signaling pathway is correlated with infertility in patients with endometriosis. Taken together, our analysis has identified 29 hub genes which might lead to infertility in endometriosis through the regulation of the GPCR network.