Different DNA methylome, transcriptome and histological features in uterine fibroids with and without MED12 mutations

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 tetratricopeptide 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.

Hypoxia in uterine fibroids: role in pathobiology and therapeutic opportunities

Uterine fibroids are the most common tumors in females affecting up to 70% of women world-wide, yet targeted therapeutic options are limited. Oxidative stress has recently surfaced as a key driver of fibroid pathogenesis and provides insights into hypoxia-induced cell transformation, extracellular matrix pathophysiology, hypoxic cell signaling cascades, and uterine biology. Hypoxia drives fibroid tumorigenesis through (1) promoting myometrial stem cell proliferation, (2) causing DNA damage propelling transformation of stem cells to tumor initiating cells, and (3) driving excess extracellular matrix (ECM) production. Common fibroid-associated DNA mutations include MED12 mutations, HMGA2 overexpression, and Fumarate hydratase loss of function. Evidence suggests an interaction between hypoxia signaling and these mutations. Fibroid development and growth are promoted by hypoxia-triggered cell signaling via various pathways including HIF-1, TGFβ, and Wnt/β-catenin. Fibroid-associated hypoxia persists due to antioxidant imbalance, ECM accumulation, and growth beyond adequate vascular supply. Current clinically available fibroid treatments do not take advantage of hypoxia-targeting therapies. Growing pre-clinical and clinical studies identify ROS inhibitors, anti-HIF-1 agents, Wnt/β-catenin inhibition, and TGFβ cascade inhibitors as agents that may reduce fibroid development and growth through targeting hypoxia.

Differential gene expression in decidualized human endometrial stromal cells induced by different stimuli

Decidualization can be induced by culturing human endometrial stromal cells (ESCs) with several decidualization stimuli, such as cAMP, medroxyprogesterone acetate (MPA) or Estradiol (E2). However, it has been unclear how decidualized cells induced by different stimuli are different. We compared transcriptomes and cellular functions of decidualized ESCs induced by different stimuli (MPA, E2 + MPA, cAMP, and cAMP + MPA). We also investigated which decidualization stimulus induces a closer in vivo decidualization. Differentially expressed genes (DEGs) and altered cellular functions by each decidualization stimuli were identified by RNA-sequence and gene-ontology analysis. DEGs was about two times higher for stimuli that use cAMP (cAMP and cAMP + MPA) than for stimuli that did not use cAMP (MPA and E2 + MPA). cAMP-using stimuli altered the cellular functions including angiogenesis, inflammation, immune system, and embryo implantation whereas MPA-using stimuli (MPA, E2 + MPA, and cAMP + MPA) altered the cellular functions associated with insulin signaling. A public single-cell RNA-sequence data of the human endometrium was utilized to analyze in vivo decidualization. The altered cellular functions by in vivo decidualization were close to those observed by cAMP + MPA-induced decidualization. In conclusion, decidualized cells induced by different stimuli have different transcriptome and cellular functions. cAMP + MPA may induce a decidualization most closely to in vivo decidualization.

Establishment of Noninvasive Prediction Models for the Diagnosis of Uterine Leiomyoma Subtypes

OBJECTIVE

To establish prediction models for the diagnosis of the subtypes of uterine leiomyomas by machine learning using magnetic resonance imaging (MRI) data.

METHODS

This is a prospective observational study. Ninety uterine leiomyoma samples were obtained from 51 patients who underwent surgery for uterine leiomyomas. Seventy-one samples (49 mediator complex subunit 12 [ MED12 ] mutation-positive and 22 MED12 mutation-negative leiomyomas) were assigned to the primary data set to establish prediction models. Nineteen samples (13 MED12 mutation-positive and 6 MED12 mutation-negative leiomyomas) were assigned to the unknown testing data set to validate the prediction model utility. The tumor signal intensity was quantified by seven MRI sequences (T2-weighted imaging, apparent diffusion coefficient, magnetic resonance elastography, T1 mapping, magnetization transfer contrast, T2* blood oxygenation level dependent, and arterial spin labeling) that can estimate the collagen and water contents of uterine leiomyomas. After surgery, the MED12 mutations were genotyped. These results were used to establish prediction models based on machine learning by applying support vector classification and logistic regression for the diagnosis of uterine leiomyoma subtypes. The performance of the prediction models was evaluated by cross-validation within the primary data set and then finally evaluated by external validation using the unknown testing data set.

RESULTS

The signal intensities of five MRI sequences (T2-weighted imaging, apparent diffusion coefficient, T1 mapping, magnetization transfer contrast, and T2* blood oxygenation level dependent) differed significantly between the subtypes. In cross-validation within the primary data set, both machine learning models (support vector classification and logistic regression) based on the five MRI sequences were highly predictive of the subtypes (area under the curve [AUC] 0.974 and 0.988, respectively). External validation with the unknown testing data set confirmed that both models were able to predict the subtypes for all samples (AUC 1.000, 100.0% accuracy). Our prediction models with T2-weighted imaging alone also showed high accuracy to discriminate the uterine leiomyoma subtypes.

CONCLUSION

We established noninvasive prediction models for the diagnosis of the subtypes of uterine leiomyomas by machine learning using MRI data.

Novel Approaches to Possible Targeted Therapies and Prophylaxis of Uterine Fibroids

Uterine leiomyomas are the most common benign tumors in women of childbearing age. They may lead to problems of conception or complications during the gestational period. The methods of treatment include surgical (myomectomy and hysterectomy, embolization of arteries) and therapeutic treatment (ulipristal acetate, leuprolide acetate, cetrorelix, goserelin, mifepristone). Both approaches are efficient but incompatible with pregnancy planning. Therefore, there is a call for medical practice to develop therapeutical means of preventing leiomyoma onset in patients planning on becoming pregnant. Based on the analysis of GWAS data on the search for mononucleotide polymorphisms associated with the risk of leiomyoma, in meta-transcriptomic and meta-methylomic studies, target proteins have been proposed. Prospective therapeutic treatments of leiomyoma may be based on chemical compounds, humanized recombinant antibodies, vaccines based on markers of the uterine leiomyoma cells that are absent in the adult organism, or DNA and RNA preparations. Three different nosological forms of the disease associated with driver mutations in the MED12, HMGA2, and FH genes should be considered when developing or prescribing drugs. For example, synthetic inhibitors and vaccines based on matrix metalloproteinases MMP11 and MMP16 are expected to be effective only for the prevention of the occurrence of MED12-dependent nodules.

MED12 mutation activates the tryptophan/kynurenine/AHR pathway to promote growth of uterine leiomyomas

Uterine leiomyomas cause heavy menstrual bleeding, anemia, and pregnancy loss in millions of women worldwide. Driver mutations in the transcriptional mediator complex subunit 12 (MED12) gene in uterine myometrial cells initiate 70% of leiomyomas that grow in a progesterone-dependent manner. We showed a distinct chromatin occupancy landscape of MED12 in mutant MED12 (mut-MED12) versus WT-MED12 leiomyomas. Integration of cistromic and transcriptomics data identified tryptophan 2,3-dioxygenase (TDO2) as the top mut-MED12 target gene that was significantly upregulated in mut-MED12 leiomyomas when compared with adjacent myometrium and WT-MED12 leiomyomas. TDO2 catalyzes the conversion of tryptophan to kynurenine, an aryl hydrocarbon receptor (AHR) ligand that we confirmed to be significantly elevated in mut-MED12 leiomyomas. Treatment of primary mut-MED12 leiomyoma cells with tryptophan or kynurenine stimulated AHR nuclear translocation, increased proliferation, inhibited apoptosis, and induced AHR-target gene expression, whereas blocking the TDO2/kynurenine/AHR pathway by siRNA or pharmacological treatment abolished these effects. Progesterone receptors regulated the expression of AHR and its target genes. In vivo, TDO2 expression positively correlated with the expression of genes crucial for leiomyoma growth. In summary, activation of the TDO2/kynurenine/AHR pathway selectively in mut-MED12 leiomyomas promoted tumor growth and may inform the future development of targeted treatments and precision medicine.

Developmental reprogramming of myometrial stem cells by endocrine disruptor linking to risk of uterine fibroids

BACKGROUND

The stage, when tissues and organs are growing, is very vulnerable to environmental influences, but it's not clear how exposure during this time causes changes to the epigenome and increases the risk of hormone-related illnesses like uterine fibroids (UFs).

METHODS

Developmental reprogramming of myometrial stem cells (MMSCs), the putative origin from which UFs originate, was investigated in vitro and in the Eker rat model by RNA-seq, ChIP-seq, RRBS, gain/loss of function analysis, and luciferase activity assays.

RESULTS

When exposed to the endocrine-disrupting chemical (EDC) diethylstilbestrol during Eker rat development, MMSCs undergo a reprogramming of their estrogen-responsive transcriptome. The reprogrammed genes in MMSCs are known as estrogen-responsive genes (ERGs) and are activated by mixed lineage leukemia protein-1 (MLL1) and DNA hypo-methylation mechanisms. Additionally, we observed a notable elevation in the expression of ERGs in MMSCs from Eker rats exposed to natural steroids after developmental exposure to EDC, thereby augmenting estrogen activity.

CONCLUSION

Our studies identify epigenetic mechanisms of MLL1/DNA hypo-methylation-mediated MMSC reprogramming. EDC exposure epigenetically targets MMSCs and leads to persistent changes in the expression of a subset of ERGs, imparting a hormonal imprint on the ERGs, resulting in a "hyper-estrogenic" phenotype, and increasing the hormone-dependent risk of UFs.

A View on Uterine Leiomyoma Genesis through the Prism of Genetic, Epigenetic and Cellular Heterogeneity

Uterine leiomyomas (ULs), frequent benign tumours of the female reproductive tract, are associated with a range of symptoms and significant morbidity. Despite extensive research, there is no consensus on essential points of UL initiation and development. The main reason for this is a pronounced inter- and intratumoral heterogeneity resulting from diverse and complicated mechanisms underlying UL pathobiology. In this review, we comprehensively analyse risk and protective factors for UL development, UL cellular composition, hormonal and paracrine signalling, epigenetic regulation and genetic abnormalities. We conclude the need to carefully update the concept of UL genesis in light of the current data. Staying within the framework of the existing hypotheses, we introduce a possible timeline for UL development and the associated key events—from potential prerequisites to the beginning of UL formation and the onset of driver and passenger changes.

Differential Expression of MED12-Associated Coding RNA Transcripts in Uterine Leiomyomas

Recent studies have demonstrated that somatic MED12 mutations in exon 2 occur at a frequency of up to 80% and have a functional role in leiomyoma pathogenesis. The objective of this study was to elucidate the expression profile of coding RNA transcripts in leiomyomas, with and without these mutations, and their paired myometrium. Next-generation RNA sequencing (NGS) was used to systematically profile the differentially expressed RNA transcripts from paired leiomyomas (n = 19). The differential analysis indicated there are 394 genes differentially and aberrantly expressed only in the mutated tumors. These genes were predominantly involved in the regulation of extracellular constituents. Of the differentially expressed genes that overlapped in the two comparison groups, the magnitude of change in gene expression was greater for many genes in tumors bearing MED12 mutations. Although the myometrium did not express MED12 mutations, there were marked differences in the transcriptome landscape of the myometrium from mutated and non-mutated specimens, with genes regulating the response to oxygen-containing compounds being most altered. In conclusion, MED12 mutations have profound effects on the expression of genes pivotal to leiomyoma pathogenesis in the tumor and the myometrium which could alter tumor characteristics and growth potential.

Combined histological and DNA methylome profiling approaches may provide insights into the pathophysiology of ovarian endometriomas

Purpose

To test the theory that invaginated ovarian surface epithelium and endometrial implants on the ovary form ovarian endometriomas.

Methods

Adhesion sites of ovarian endometrioma on the peritoneum and consecutive ovarian endometrioma cyst wall, called non-adhesion sites, were histologically examined. DNA methylomes of the adhesion sites, non-adhesion sites, and blueberry spots were compared with those of ovary, endometrium, and peritoneum.

Results

The non-adhesion sites showed an ovarian surface epithelium-like structure near the adhesion site, which continued to a columnar epithelium-like structure. Calretinin staining was strong in the ovarian surface epithelium-like structure but weak in the columnar epithelium-like structure. Estrogen receptors were absent in the ovarian surface epithelium-like structure, but present in the columnar epithelium-like structure. The adhesion sites had endometrial gland-like structures that expressed estrogen receptors. Analyses of DNA methylomes classified the non-adhesion sites and ovaries into the same group, suggesting that ovarian endometriomas originate from the ovarian surface epithelium. The adhesion sites, blueberry spots and peritoneum were classified in the same group, suggesting that the adhesion sites and blueberry spots originate from the peritoneum.

Conclusions

The present results support the invagination theory. Ovarian endometriomas consist of invaginated ovarian surface epithelium with celomic metaplasia and endometrium implants on the peritoneum.

PLP1 may serve as a potential diagnostic biomarker of uterine fibroids

Objective: We aim to identify the crucial genes or potential biomarkers associated with uterine fibroids (UFs), which may provide clinicians with evidence about the diagnostic biomarker of UFs and reveal the mechanism of its progression.

Methods: The gene expression and genome-wide DNA methylation profiles were obtained from Gene Expression Omnibus database (GEO). GSE45189, GSE31699, and GSE593 datasets were included. GEO2R and Venn diagrams were used to analyze the differentially expressed genes (DEGs) and extract the hub genes. Gene Ontology (GO) analysis was performed by the online tool Database for Annotation, Visualization, and Integrated Discovery (DAVID). The mRNA and protein expression of hub genes were validated by RT-qPCR, western blot, and immunohistochemistry. The receiver operating characteristic (ROC) curve was used to evaluate the diagnostic value.

Results: We detected 22 DEGs between UFs and normal myometrium, which were enriched in cell maturation, apoptotic process, hypoxia, protein binding, and cytoplasm for cell composition. By finding the intersection of the data between differentially expressed mRNA and DNA methylation profiles, 3 hub genes were identified, including transmembrane 4 L six family member 1 (TM4SF1), TNF superfamily member 10 (TNFSF10), and proteolipid protein 1 (PLP1). PLP1 was validated to be up-regulated significantly in UFs both at mRNA and protein levels. The area under the ROC curve (AUC) of PLP1 was 0.956, with a sensitivity of 79.2% and a specificity of 100%. Conclusion: Overall, our results indicate that PLP1 may be a potential diagnostic biomarker for uterine fibroids.