Clonally related uterine leiomyomas are common and display branched tumor evolution

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.

3'RNA and whole-genome sequencing of archival uterine leiomyomas reveal a tumor subtype with chromosomal rearrangements affecting either HMGA2, HMGA1, or PLAG1

Uterine leiomyomas, or fibroids, are very common smooth muscle tumors that arise from the myometrium. They can be divided into distinct molecular subtypes. We have previously shown that 3'RNA-sequencing is highly effective in classifying archival formalin-fixed paraffin-embedded (FFPE) leiomyomas according to the underlying mutation. In this study, we performed 3'RNA-sequencing with 111 FFPE leiomyomas previously classified as negative for driver alterations in mediator complex subunit 12 (MED12), high mobility group AT-hook 2 (HMGA2), and fumarate hydratase (FH) by Sanger sequencing and immunohistochemistry. This revealed 43 tumors that displayed expression features typically seen in HMGA2-positive tumors, including overexpression of PLAG1. We explored 12 such leiomyomas by whole-genome sequencing to identify their underlying genomic drivers and to evaluate the feasibility of detecting chromosomal driver alterations from FFPE material. Four tumors with significant HMGA2 overexpression at the protein-level served as controls. We identified chromosomal rearrangements targeting either HMGA2, HMGA1, or PLAG1 in all 16 tumors, demonstrating that it is possible to detect chromosomal driver alterations in archival leiomyoma specimens as old as 18 years. Furthermore, two tumors displayed biallelic loss of DEPDC5 and one tumor harbored a COL4A5-COL4A6 deletion. These observations suggest that instead of only HMGA2-positive leiomyomas, a distinct leiomyoma subtype is characterized by rearrangements targeting either HMGA2, HMGA1, or PLAG1. The results indicate that the frequency of HMGA2-positive leiomyomas may be higher than estimated in previous studies where immunohistochemistry has been used. This study also demonstrates the feasibility of detecting chromosomal driver alterations from archival FFPE material.

Comparison of 2SC, AKR1B10, and FH Antibodies as Potential Biomarkers for FH-deficient Uterine Leiomyomas

Hereditary leiomyomatosis and renal cell cancer (HLRCC) is a tumor predisposition syndrome caused by germline fumarate hydratase (FH) mutations and characterized by uterine and cutaneous leiomyomas and renal cell cancer. Currently, there is no generally approved method to differentiate FH-deficient uterine leiomyomas from other leiomyomas. Here, we analyzed 3 antibodies (S-(2-succino)-cysteine [2SC], aldo-keto reductase family 1, member B10 [AKR1B10], and FH) as potential biomarkers. The study consisted of 2 sample series. The first series included 155 formalin-fixed paraffin-embedded uterine leiomyomas, of which 90 were from HLRCC patients and 65 were sporadic. The second series included 1590 unselected fresh frozen leiomyomas. Twenty-seven tumors were from known HLRCC patients, while the FH status for the remaining 1563 tumors has been determined by copy number analysis and Sanger sequencing revealing 45 tumors with monoallelic (n=33) or biallelic (n=12) FH loss. Altogether 197 samples were included in immunohistochemical analyses: all 155 samples from series 1 and 42 available corresponding formalin-fixed paraffin-embedded samples from series 2 (15 tumors with monoallelic and 7 with biallelic FH loss, 20 with no FH deletion). Results show that 2SC performed best with 100% sensitivity and specificity. Scoring was straightforward with unambiguously positive or negative results. AKR1B10 identified most tumors accurately with 100% sensitivity and 99% specificity. FH was 100% specific but showed slightly reduced 91% sensitivity. Both FH and AKR1B10 displayed also intermediate staining intensities. We suggest that when patient's medical history and/or histopathologic tumor characteristics indicate potential FH-deficiency, the tumor's FH status is determined by 2SC staining. When aberrant staining is observed, the patient can be directed to genetic counseling and mutation screening.

Cytogenomic Profile of Uterine Leiomyoma: In Vivo vs. In Vitro Comparison

We performed a comparative cytogenomic analysis of cultured and uncultured uterine leiomyoma (UL) samples. The experimental approach included karyotyping, aCGH, verification of the detected chromosomal abnormalities by metaphase and interphase FISH, MED12 mutation analysis and telomere measurement by Q-FISH. An abnormal karyotype was detected in 12 out of 32 cultured UL samples. In five karyotypically abnormal ULs, MED12 mutations were found. The chromosomal abnormalities in ULs were present mostly by complex rearrangements, including chromothripsis. In both karyotypically normal and abnormal ULs, telomeres were ~40% shorter than in the corresponding myometrium, being possibly prerequisite to chromosomal rearrangements. The uncultured samples of six karyotypically abnormal ULs were checked for the detected chromosomal abnormalities through interphase FISH with individually designed DNA probe sets. All chromosomal abnormalities detected in cultured ULs were found in corresponding uncultured samples. In all tumors, clonal spectra were present by the karyotypically abnormal cell clone/clones which coexisted with karyotypically normal ones, suggesting that chromosomal abnormalities acted as drivers, rather than triggers, of the neoplastic process. In vitro propagation did not cause any changes in the spectrum of the cell clones, but altered their ratio compared to uncultured sample. The alterations were unique for every UL. Compared to its uncultured counterpart, the frequency of chromosomally abnormal cells in the cultured sample was higher in some ULs and lower in others. To summarize, ULs are characterized by both inter- and intratumor genetic heterogeneity. Regardless of its MED12 status, a tumor may be comprised of clones with and without chromosomal abnormalities. In contrast to the clonal spectrum, which is unique and constant for each UL, the clonal frequency demonstrates up or down shifts under in vitro conditions, most probably determined by the unequal ability of cells with different genetic aberrations to exist outside the body.

Endometriosis and Uterine Fibroids (Leiomyomata): Comorbidity, Risks and Implications

Uterine Fibroids (leiomyomata) and endometriosis affect millions of women world-wide. Although aetiology and natural history of the conditions are markedly different, symptoms can overlap and make differential diagnoses necessary, often using invasive methods such as laparoscopy. Considerable comorbidity exists between the two conditions and needs to be taken into account when treating fibroids and/or endometriosis. The genetic foundations of both uterine fibroids and endometriosis remain to be fully understood but recent evidence suggest common underpinnings. Here, we discuss the comorbidity of uterine fibroids and endometriosis and the implications for diagnosis, treatment and risks.

MIPUP: minimum perfect unmixed phylogenies for multi-sampled tumors via branchings and ILP

Motivation

Discovering the evolution of a tumor may help identify driver mutations and provide a more comprehensive view on the history of the tumor. Recent studies have tackled this problem using multiple samples sequenced from a tumor, and due to clinical implications, this has attracted great interest. However, such samples usually mix several distinct tumor subclones, which confounds the discovery of the tumor phylogeny.

Results

We study a natural problem formulation requiring to decompose the tumor samples into several subclones with the objective of forming a minimum perfect phylogeny. We propose an Integer Linear Programming formulation for it, and implement it into a method called MIPUP. We tested the ability of MIPUP and of four popular tools LICHeE, AncesTree, CITUP, Treeomics to reconstruct the tumor phylogeny. On simulated data, MIPUP shows up to a 34% improvement under the ancestor-descendant relations metric. On four real datasets, MIPUP’s reconstructions proved to be generally more faithful than those of LICHeE.

Availability and implementation

MIPUP is available at https://github.com/zhero9/MIPUP as open source.

Supplementary information

Supplementary data are available at Bioinformatics online.

Evidence of biomechanical and collagen heterogeneity in uterine fibroids

Objective

Uterine fibroids (leiomyomas) are common benign tumors of the myometrium but their molecular pathobiology remains elusive. These stiff and often large tumors contain abundant extracellular matrix (ECM), including large amounts of collagen, and can lead to significant morbidities. After observing structural multiformities of uterine fibroids, we aimed to explore this heterogeneity by focusing on collagen and tissue stiffness.

Methods

For 19 fibroids, ranging in size from 3 to 11 centimeters, from eight women we documented gross appearance and evaluated collagen content by Masson trichrome staining. Collagen types were determined in additional samples by serial extraction and gel electrophoresis. Biomechanical stiffness was evaluated by rheometry.

Results

Fibroid slices displayed different gross morphology and some fibroids had characteristics of two or more patterns: classical whorled (n = 8); nodular (n = 9); interweaving trabecular (n = 9); other (n = 1). All examined fibroids contained at least 37% collagen. Tested samples included type I, III, and V collagen of different proportions. Fibroid stiffness was not correlated with the overall collagen content (correlation coefficient 0.22). Neither stiffness nor collagen content was correlated with fibroid size. Stiffness among fibroids ranged from 3028 to 14180 Pa (CV 36.7%; p<0.001, one-way ANOVA). Stiffness within individual fibroids was also not uniform and variability ranged from CV 1.6 to 42.9%.

Conclusions

The observed heterogeneity in structure, collagen content, and stiffness highlights that fibroid regions differ in architectural status. These differences might be associated with variations in local pressure, biomechanical signaling, and altered growth. We conclude the design of all fibroid studies should account for such heterogeneity because samples from different regions have different characteristics. Our understanding of fibroid pathophysiology will greatly increase through the investigation of the complexity of the chemical and biochemical signaling in fibroid development, the correlation of collagen content and mechanical properties in uterine fibroids, and the mechanical forces involved in fibroid development as affected by the various components of the ECM.

Genetic predisposition to uterine leiomyoma is determined by loci for genitourinary development and genome stability

Uterine leiomyomas (ULs) are benign tumors that are a major burden to women’s health. A genome-wide association study on 15,453 UL cases and 392,628 controls was performed, followed by replication of the genomic risk in six cohorts. Effects of the risk alleles were evaluated in view of molecular and clinical characteristics. 22 loci displayed a genome-wide significant association. The likely predisposition genes could be grouped to two biological processes. Genes involved in genome stability were represented by TERT, TERC, OBFC1 - highlighting the role of telomere maintenance - TP53 and ATM. Genes involved in genitourinary development, WNT4, WT1, SALL1, MED12, ESR1, GREB1, FOXO1, DMRT1 and uterine stem cell marker antigen CD44, formed another strong subgroup. The combined risk contributed by the 22 loci was associated with MED12 mutation-positive tumors. The findings link genes for uterine development and genetic stability to leiomyomagenesis, and in part explain the more frequent occurrence of UL in women of African origin.

Fibroids – also known as uterine leiomyomas, or myomas – are a very common form of benign tumor that grows in the muscle wall of the uterus. As many as 70% of women develop fibroids in their lifetime. About a fifth of women report symptoms including severe pain, heavy bleeding during periods and complications in pregnancy. In the United States, the cost of treating fibroids is estimated to be $34 billion each year.

Despite the prevalence of fibroids in women, there are few treatments available. Drugs to target them have limited effect and often an invasive procedure such as surgery is needed to remove the tumors. However, a better understanding of the genetics of fibroids could lead to a way to develop better treatment options.

Välimäki, Kuisma et al. used a genome-wide association study to seek out DNA variations that are more common in people with fibroids. Using data from the UK Biobank, the genomes of over 15,000 women with fibroids were analyzed against a control population of over 392,000 individuals. The analysis revealed 22 regions of the genome that were associated with fibroids. These regions included genes that may well contribute to fibroid development, such as the gene TP53, which influences the stability of the genome, and ESR1, which codes for a receptor for estrogen – a hormone known to play a role in the growth of fibroids. Variation in a set of genes known to control development of the female reproductive organs was also identified in women with fibroids.

The findings are the result of the largest genome-wide association study on fibroids, revealing a set of genes that could influence the development of fibroids. Studying these genes could lead to more effective drug development to treat fibroids. Revealing this group of genes could also help to identify women at high risk of developing fibroids and help to prevent or manage the condition.

Factors targeting MED12 to drive tumorigenesis?

Mediator Complex Subunit 12 (MED12) is part of the transcriptional preinitiation machinery. Mutations of its gene predominantly occur in two types of highly frequent benign tumors, uterine leiomyomas and fibroadenomas of the breast, where they apparently act as driver mutations. Nevertheless, their presence is not restricted to benign tumors having been found at considerable frequencies in uterine leiomyosarcomas, malignant phyllodes tumors, and chronic lymphocytic leukemia also. Most of the mutations are located within exon 2 of the gene but in rare cases the intron 1/exon 2 boundary or exon 1 are affected. As to their type, predominantly single nucleotide exchanges with a hotspot in one codon are found, but small deletions clustering around that hotspot also are not uncommon. These latter deletions are leaving the open reading frame intact. As to the types of mutations, so far no apparent differences between the tumor entities affected have emerged. Interestingly, this pattern with small deletions clustered around the hotspot of single nucleotide exchanges resembles that seen as a result of targeted gene editing. In contrast to other driver mutations the percentage of MED12-mutation positive tumors of independent clonal origin increases with the number of tumors per patient suggesting unknown etiological factors supporting site specific mutagenesis.  These factors may act by inducing simultaneous site-specific double strand breaks the erroneous repair of which may lead to corresponding mutations. As inducers of DNA damage and its repair such as foreign nucleic acids of the microbiome displaying sequence homology to the putative target site might play a role. Interestingly, a 16 base pair homology of the hotspot to a putative terminator base-paired hairpin sequence of a Staphylococcus aureus tRNA gene cluster has been noted which might form R-loop like structures with its target sequence thus inducing said changes.

Clinical, pathologic, cytogenetic, and molecular profiling in self-identified black women with uterine leiomyomata

Black women are disproportionately affected by uterine leiomyomata (UL), or fibroids, compared to other racial groups, having a greater lifetime risk of developing UL and an earlier age of diagnosis. In order to elucidate molecular and genetic mechanisms responsible for the increased prevalence and morbidity associated with UL in black women, clinical, pathologic, cytogenetic, and select molecular profiling (MED12 mutation analysis) of 75 self-reported black women undergoing surgical treatment for UL was performed. Our observations are broadly representative of previous cytogenetic studies of UL: karyotypically abnormal tumors were detected in 30.7% of women and 17.4% of analyzed tumors. No notable association was observed between race and increased occurrence of cytogenetic abnormalities that might contribute to any population-specific morbidity or prevalence rate. Our data on MED12 mutation analyses (73.2% of tumors harbored a MED12 mutation) provide additional support for a significant role of MED12 in tumorigenesis. Although the effect of MED12-mediated tumorigenesis appears significant irrespective of race, other genetic events such as the distribution of karyotypic abnormalities appear differently in black women. This case series indicates that presently recognized genetic and molecular characteristics of UL do not appear to explain the increased prevalence and morbidity of UL in black women.

Somatic evolutionary timings of driver mutations

BACKGROUND

A unified analysis of DNA sequences from hundreds of tumors concluded that the driver mutations primarily occur in the earliest stages of cancer formation, with relatively few driver mutation events detected in the late-arising subclones. However, emerging evidence from the sequencing of multiple tumors and tumor regions per individual suggests that late-arising subclones with additional driver mutations are underestimated in single-sample analyses.

METHODS

To test whether driver mutations generally map to early tumor development, we examined multi-regional tumor sequencing data from 101 individuals reported in 11 published studies. Following previous studies, we annotated mutations as early-arising when all tumors/regions had those mutations (ubiquitous). We then inferred the fraction of mutations occurring early and compared it with late-arising mutations that were found in only single tumors/regions.

RESULTS

While a large fraction of driver mutations in tumors occurred relatively early in cancers, later driver mutations occurred at least as frequently as the early drivers in a substantial number of patients. This result was robust to many different approaches to annotate driver mutations. The relative frequency of early and late driver mutations varied among patients of the same cancer type and in different cancer types. We found that previous reports of the preponderance of early driver mutations were primarily informed by analysis of single tumor variant allele profiles, with which it is challenging to clearly distinguish between early and late drivers.

CONCLUSIONS

The origin and preponderance of new driver mutations are not limited to early stages of tumor evolution, with different tumors and regions showing distinct driver mutations and, consequently, distinct characteristics. Therefore, tumors with extensive intratumor heterogeneity appear to have many newly acquired drivers.

Subtype-Specific Tumor-Associated Fibroblasts Contribute to the Pathogenesis of Uterine Leiomyoma

Recent genomic studies have identified subtypes of uterine leiomyoma (LM) with distinctive genetic alterations. Here, we report the elucidation of the biological characteristics of the two most prevalent uterine leiomyoma subtypes, MED12-mutant (MED12-LM) and HMGA2-overexpressing (HMGA2-LM) uterine leiomyomas. Because each tumor carries only one genetic alteration, both subtypes are considered to be monoclonal. Approximately 90% of cells in HMGA2-uterine leiomyoma were smooth muscle cells (SMC) with HMGA2 overexpression. In contrast, MED12-LM consisted of similar numbers of SMC and non-SMC, which were mostly tumor-associated fibroblasts (TAF). Paradoxically, TAF carried no mutations in MED12, suggesting an interaction between SMC and TAF to coordinate their growth. The higher amount of extracellular matrix in MED12-LM than HMGA2-LM was partially due to the high concentration of collagen-producing TAF. SMC growth in a xenograft assay was driven by progesterone in both uterine leiomyoma subtypes. In contrast, TAF in MED12-LM proliferated in response to estradiol, whereas progesterone had no effect. The high concentration of estrogen-responsive TAF in MED12-LM explains the inconsistent discoveries between in vivo and in vitro studies on the mitogenic effect of estrogen and raises questions regarding the accuracy of previous studies utilizing MED12-LM cell culture. In addition, the differential effects of estradiol and progesterone on these uterine leiomyoma subtypes emphasize the importance of subtypes and genotypes in designing nonsurgical therapeutic strategies for uterine leiomyoma. Cancer Res; 77(24); 6891-901. ©2017 AACR.

Prevalence and clinical significance of mediator complex subunit 12 mutations in 362 Han Chinese samples with uterine leiomyoma

Uterine leiomyomas (ULs) are the most common gynecological benign tumors originating from the myometrium. Prevalent mutations in the mediator complex subunit 12 (MED12) gene have been identified in ULs, and functional evidence has revealed that these mutations may promote the development of ULs. However, whether MED12 mutations are associated with certain clinical characteristics in ULs remains largely unknown. In the present study, the potential mutations of MED12 and its paralogous gene, mediator complex subunit 12-like (MED12L), were screened in 362 UL tumors from Han Chinese patients. A total of 158 out of 362 UL tumors (43.6%) were identified as harboring MED12 somatic mutations, and the majority of these mutations were restricted to the 44th residue. MED12 mutations were also observed in 2 out of 145 (1.4%) adjacent control myometrium. Furthermore, the mutation spectrum of MED12 in the concurrent leiomyomas was noticeably different. Correlation analysis of MED12 mutations with the available clinical features indicated that patients with mutated MED12 tended to have smaller cervical diameters. By contrast, no MED12L mutation was identified in the present samples. In summary, the present study demonstrated the presence of prevalent MED12 somatic mutations in UL samples, and the MED12 mutation was associated with smaller cervical diameters. The low mutation frequency of MED12 in adjacent control myometrium indicated that MED12 mutation may be an early event in the pathogenesis of ULs. Furthermore, MED12 mutation status in concurrent tumors from multiple leiomyomas supported several prior observations that the majority of these tumors arose independently.

Integrated data analysis reveals uterine leiomyoma subtypes with distinct driver pathways and biomarkers

Uterine leiomyomas are common benign smooth muscle tumors that impose a major burden on women's health. Recent sequencing studies have revealed recurrent and mutually exclusive mutations in leiomyomas, suggesting the involvement of molecularly distinct pathways. In this study, we explored transcriptional differences among leiomyomas harboring different genetic drivers, including high mobility group AT-hook 2 (HMGA2) rearrangements, mediator complex subunit 12 (MED12) mutations, biallelic inactivation of fumarate hydratase (FH), and collagen, type IV, alpha 5 and collagen, type IV, alpha 6 (COL4A5-COL4A6) deletions. We also explored the transcriptional consequences of 7q22, 22q, and 1p deletions, aiming to identify possible target genes. We investigated 94 leiomyomas and 60 corresponding myometrial tissues using exon arrays, whole genome sequencing, and SNP arrays. This integrative approach revealed subtype-specific expression changes in key driver pathways, including Wnt/β-catenin, Prolactin, and insulin-like growth factor (IGF)1 signaling. Leiomyomas with HMGA2 aberrations displayed highly significant up-regulation of the proto-oncogene pleomorphic adenoma gene 1 (PLAG1), suggesting that HMGA2 promotes tumorigenesis through PLAG1 activation. This was supported by the identification of genetic PLAG1 alterations resulting in expression signatures as seen in leiomyomas with HMGA2 aberrations. RAD51 paralog B (RAD51B), the preferential translocation partner of HMGA2, was up-regulated in MED12 mutant lesions, suggesting a role for this gene in the genesis of leiomyomas. FH-deficient leiomyomas were uniquely characterized by activation of nuclear factor erythroid 2-related factor 2 (NRF2) target genes, supporting the hypothesis that accumulation of fumarate leads to activation of the oncogenic transcription factor NRF2. This study emphasizes the need for molecular stratification in leiomyoma research and possibly in clinical practice as well. Further research is needed to determine whether the candidate biomarkers presented herein can provide guidance for managing the millions of patients affected by these lesions.