Multiomic analysis of uterine leiomyomas in self-described Black and White women: molecular insights into health disparities.

BACKGROUND: Black women are at an increased risk of developing uterine leiomyomas and experiencing worse disease prognosis than White women. Epidemiologic and molecular factors have been identified as underlying these disparities, but there remains a paucity of deep, multiomic analysis investigating molecular differences in uterine leiomyomas from Black and White patients. OBJECTIVE: To identify molecular alterations within uterine leiomyoma tissues correlating with patient race by multiomic analyses of uterine leiomyomas collected from cohorts of Black and White women. STUDY DESIGN: We performed multiomic analysis of uterine leiomyomas from Black (42) and White (47) women undergoing hysterectomy for symptomatic uterine leiomyomata. In addition, our analysis included the application of orthogonal methods to evaluate fibroid biomechanical properties, such as second harmonic generation microscopy, uniaxial compression testing, and shear-wave ultrasonography analyses. RESULTS: We found a greater proportion of MED12 mutant uterine leiomyomas from Black women (>35% increase; Mann-Whitney U, P<.001). MED12 mutant tumors exhibited an elevated abundance of extracellular matrix proteins, including several collagen isoforms, involved in the regulation of the core matrisome. Histologic analysis of tissue fibrosis using trichrome staining and secondary harmonic generation microscopy confirmed that MED12 mutant tumors are more fibrotic than MED12 wild-type tumors. Using shear-wave ultrasonography in a prospectively collected cohort, Black patients had fibroids that were firmer than White patients, even when similar in size. In addition, these analyses uncovered ancestry-linked expression quantitative trait loci with altered allele frequencies in African and European populations correlating with differential abundance of several proteins in uterine leiomyomas independently of MED12 mutation status, including tetracoidpeptide repeat protein 38. CONCLUSION: Our study shows that Black women have a higher prevalence of uterine leiomyomas harboring mutations in MED12 and that this mutational status correlates with increased tissue fibrosis compared with wild-type uterine leiomyomas. Our study provides insights into molecular alterations correlating with racial disparities in uterine leiomyomas and improves our understanding of the molecular etiology underlying uterine leiomyoma development within these populations.

ProteoMixture: A cell type deconvolution tool for bulk tissue proteomic data.

Numerous multi-omic investigations of cancer tissue have documented varying and poor pairwise transcript:protein quantitative correlations, and most deconvolution tools aiming to predict cell type proportions (cell admixture) have been developed and credentialed using transcript-level data alone. To estimate cell admixture using protein abundance data, we analyzed proteome and transcriptome data generated from contrived admixtures of tumor, stroma, and immune cell models or those selectively harvested from the tissue microenvironment by laser microdissection from high grade serous ovarian cancer (HGSOC) tumors. Co-quantified transcripts and proteins performed similarly to estimate stroma and immune cell admixture (r ≥ 0.63) in two commonly used deconvolution algorithms, ESTIMATE or ConsensusTME. We further developed and optimized protein-based signatures estimating cell admixture proportions and benchmarked these using bulk tumor proteomic data from over 150 patients with HGSOC. The optimized protein signatures supporting cell type proportion estimates from bulk tissue proteomic data are available at https://lmdomics.org/ProteoMixture/.

Proteogenomic analysis of enriched HGSOC tumor epithelium identifies prognostic signatures and therapeutic vulnerabilities.

We performed a deep proteogenomic analysis of bulk tumor and laser microdissection enriched tumor cell populations from high-grade serous ovarian cancer (HGSOC) tissue specimens spanning a broad spectrum of purity. We identified patients with longer progression-free survival had increased immune-related signatures and validated proteins correlating with tumor-infiltrating lymphocytes in 65 tumors from an independent cohort of HGSOC patients, as well as with overall survival in an additional 126 HGSOC patient cohort. We identified that homologous recombination deficient (HRD) tumors are enriched in pathways associated with metabolism and oxidative phosphorylation that we validated in independent patient cohorts. We further identified that polycomb complex protein BMI-1 is elevated in HR proficient (HRP) tumors, that elevated BMI-1 correlates with poor overall survival in HRP but not HRD HGSOC patients, and that HRP HGSOC cells are uniquely sensitive to BMI-1 inhibition.

Mapping three-dimensional intratumor proteomic heterogeneity in uterine serous carcinoma by multiregion microsampling.

BACKGROUND: Although uterine serous carcinoma (USC) represents a small proportion of all uterine cancer cases, patients with this aggressive subtype typically have high rates of chemotherapy resistance and disease recurrence that collectively result in a disproportionately high death rate. The goal of this study was to provide a deeper view of the tumor microenvironment of this poorly characterized uterine cancer variant through multi-region microsampling and quantitative proteomics. METHODS: Tumor epithelium, tumor-involved stroma, and whole "bulk" tissue were harvested by laser microdissection (LMD) from spatially resolved levels from nine USC patient tumor specimens and underwent proteomic analysis by mass spectrometry and reverse phase protein arrays, as well as transcriptomic analysis by RNA-sequencing for one patient's tumor. RESULTS: LMD enriched cell subpopulations demonstrated varying degrees of relatedness, indicating substantial intratumor heterogeneity emphasizing the necessity for enrichment of cellular subpopulations prior to molecular analysis. Known prognostic biomarkers were quantified with stable levels in both LMD enriched tumor and stroma, which were shown to be highly variable in bulk tissue. These USC data were further used in a comparative analysis with a data generated from another serous gynecologic malignancy, high grade serous ovarian carcinoma, and have been added to our publicly available data analysis tool, the Heterogeneity Analysis Portal ( https://lmdomics.org/ ). CONCLUSIONS: Here we identified extensive three-dimensional heterogeneity within the USC tumor microenvironment, with disease-relevant biomarkers present in both the tumor and the stroma. These data underscore the critical need for upfront enrichment of cellular subpopulations from tissue specimens for spatial proteogenomic analysis.

The murine metastatic microenvironment of experimental brain metastases of breast cancer differs by host age in vivo: a proteomic study.

Breast cancer in young patients is known to exhibit more aggressive biological behavior and is associated with a less favorable prognosis than the same disease in older patients, owing in part to an increased incidence of brain metastases. The mechanistic explanations behind these findings remain poorly understood. We recently reported that young mice, in comparison to older mice, developed significantly greater brain metastases in four mouse models of triple-negative and luminal B breast cancer. Here we have performed a quantitative mass spectrometry-based proteomic analysis to identify proteins potentially contributing to age-related disparities in the development of breast cancer brain metastases. Using a mouse hematogenous model of brain-tropic triple-negative breast cancer (MDA-MB-231BR), we harvested subpopulations of tumor metastases, the tumor-adjacent metastatic microenvironment, and uninvolved brain tissues via laser microdissection followed by quantitative proteomic analysis using high resolution mass spectrometry to characterize differentially abundant proteins potentially contributing to age-dependent rates of brain metastasis. Pathway analysis revealed significant alterations in signaling pathways, particularly in the metastatic microenvironment, modulating tumorigenesis, metabolic processes, inflammation, and neuronal signaling. Tenascin C (TNC) was significantly elevated in all laser microdissection (LMD) enriched compartments harvested from young mice relative to older hosts, which was validated and confirmed by immunoblot analysis of whole brain lysates. Additional in vitro studies including migration and wound-healing assays demonstrated TNC as a positive regulator of tumor cell migration. These results provide important new insights regarding microenvironmental factors, including TNC, as mechanisms contributing to the increased brain cancer metastatic phenotype observed in young breast cancer patients.

Metronomic dosing of ovarian cancer cells with the ATR inhibitor AZD6738 leads to loss of CDC25A expression and resistance to ATRi treatment.

OBJECTIVE: ATR kinase inhibitors promote cell killing by inducing replication stress and through potentiation of genotoxic agents in gynecologic cancer cells. To explore mechanisms of acquired resistance to ATRi in ovarian cancer, we characterized ATRi-resistant ovarian cancer cells generated by metronomic dosing with the clinical ATR inhibitor AZD6738. METHODS: ATRi-resistant ovarian cancer cells (OVCAR3 and OV90) were generated by dosing with AZD6738 and assessed for sensitivity to Chk1i (LY2603618), PARPi (Olaparib) and combination with cisplatin or a CDK4/6 inhibitor (Palbociclib). Models were characterized by diverse methods including silencing CDC25A in OV90 cells and assessing impact on ATRi response. Serum proteomic analysis of ATRi-resistant OV90 xenografts was performed to identify circulating biomarker candidates of ATRi-resistance. RESULTS: AZD6738-resistant cell lines are refractory to LY2603618, but not to Olaparib or combinations with cisplatin. Cell cycle analyses showed ATRi-resistant cells exhibit G1/S arrest following AZD6738 treatment. Accordingly, combination with Palbociclib confers resistance to AZD6738. AZD6738-resistant cells exhibit altered abundances of G1/S phase regulatory proteins, including loss of CDC25A in AZD6738-resistant OV90 cells. Silencing of CDC25A in OV90 cells confers resistance to AZD6738. Serum proteomics from AZD6738-resistant OV90 xenografts identified Vitamin D-Binding Protein (GC), Apolipoprotein E (APOE) and A1 (APOA1) as significantly elevated in AZD6738-resistant backgrounds. CONCLUSIONS: We show that metronomic dosing of ovarian cancer cells with AZD6738 results in resistance to ATR/ Chk1 inhibitors, that loss of CDC25A expression represents a mechanism of resistance to ATRi treatment in ovarian cancer cells and identify several circulating biomarker candidates of CDC25A low, AZD6738-resistant ovarian cancer cells.

Multiomic analysis of homologous recombination-deficient end-stage high-grade serous ovarian cancer.

High-grade serous ovarian cancer (HGSC) is frequently characterized by homologous recombination (HR) DNA repair deficiency and, while most such tumors are sensitive to initial treatment, acquired resistance is common. We undertook a multiomics approach to interrogate molecular diversity in end-stage disease, using multiple autopsy samples collected from 15 women with HR-deficient HGSC. Patients had polyclonal disease, and several resistance mechanisms were identified within most patients, including reversion mutations and HR restoration by other means. We also observed frequent whole-genome duplication and global changes in immune composition with evidence of immune escape. This analysis highlights diverse evolutionary changes within HGSC that evade therapy and ultimately overwhelm individual patients.

Molecular Correlates of Venous Thromboembolism (VTE) in Ovarian Cancer.

BACKGROUND: The incidence of venous thromboembolism (VTE) in patients with ovarian cancer is higher than most solid tumors, ranging between 10-30%, and a diagnosis of VTE in this patient population is associated with worse oncologic outcomes. The tumor-specific molecular factors that may lead to the development of VTE are not well understood. OBJECTIVES: The aim of this study was to identify molecular features present in ovarian tumors of patients with VTE compared to those without. METHODS: We performed a multiplatform omics analysis incorporating RNA and DNA sequencing, quantitative proteomics, as well as immune cell profiling of high-grade serous ovarian carcinoma (HGSC) samples from a cohort of 32 patients with or without VTE. RESULTS: Pathway analyses revealed upregulation of both inflammatory and coagulation pathways in the VTE group. While DNA whole-exome sequencing failed to identify significant coding alterations between the groups, the results of an integrated proteomic and RNA sequencing analysis indicated that there is a relationship between VTE and the expression of platelet-derived growth factor subunit B (PDGFB) and extracellular proteins in tumor cells, namely collagens, that are correlated with the formation of thrombosis. CONCLUSIONS: In this comprehensive analysis of HGSC tumor tissues from patients with and without VTE, we identified markers unique to the VTE group that could contribute to development of thrombosis. Our findings provide additional insights into the molecular alterations underlying the development of VTE in ovarian cancer patients and invite further investigation into potential predictive biomarkers of VTE in ovarian cancer.

Peptide ancestry informative markers in uterine neoplasms from women of European, African, and Asian ancestry.

Characterization of ancestry-linked peptide variants in disease-relevant patient tissues represents a foundational step to connect patient ancestry with disease pathogenesis. Nonsynonymous single-nucleotide polymorphisms encoding missense substitutions within tryptic peptides exhibiting high allele frequencies in European, African, and East Asian populations, termed peptide ancestry informative markers (pAIMs), were prioritized from 1000 genomes. In silico analysis identified that as few as 20 pAIMs can determine ancestry proportions similarly to >260K SNPs (R2 = 0.99). Multiplexed proteomic analysis of >100 human endometrial cancer cell lines and uterine leiomyoma tissues combined resulted in the quantitation of 62 pAIMs that correlate with patient race and genotype-confirmed ancestry. Candidates include a D451E substitution in GC vitamin D-binding protein previously associated with altered vitamin D levels in African and European populations. pAIMs will support generalized proteoancestry assessment as well as efforts investigating the impact of ancestry on the human proteome and how this relates to the pathogenesis of uterine neoplasms.

Proteogenomic landscape of uterine leiomyomas from hereditary leiomyomatosis and renal cell cancer patients.

Pathogenic mutations in fumarate hydratase (FH) drive hereditary leiomyomatosis and renal cell cancer (HLRCC) and increase the risk of developing uterine leiomyomas (ULMs). An integrated proteogenomic analysis of ULMs from HLRCC (n = 16; FH-mutation confirmed) and non-syndromic (NS) patients (n = 12) identified a significantly higher protein:transcript correlation in HLRCC (R = 0.35) vs. NS ULMs (R = 0.242, MWU p = 0.0015). Co-altered proteins and transcripts (228) included antioxidant response element (ARE) target genes, such as thioredoxin reductase 1 (TXNRD1), and correlated with activation of NRF2-mediated oxidative stress response signaling in HLRCC ULMs. We confirm 185 transcripts previously described as altered between HLRCC and NS ULMs, 51 co-altered at the protein level and several elevated in HLRCC ULMs are involved in regulating cellular metabolism and glycolysis signaling. Furthermore, 367 S-(2-succino)cysteine peptides were identified in HLRCC ULMs, of which sixty were significantly elevated in HLRCC vs. NS ULMs (LogFC = 1.86, MWU p < 0.0001). These results confirm and define novel proteogenomic alterations in uterine leiomyoma tissues collected from HLRCC patients and underscore conserved molecular alterations correlating with inactivation of the FH tumor suppressor gene.