Genome-wide Association Study: Not so benign

Susceptibility to uterine fibroids, benign tumors that affect the health of many women, is linked to genes that are responsible for preserving genome integrity and promoting genitourinary development.
  1. Barbara Rivera  Is a corresponding author
  1. McGill University, Canada
  2. Jewish General Hospital, Canada

The established consensus is that benign tumors do not become malignant and do not metastasize, but this does not mean they are not clinically important. Uterine leiomyomas are a perfect example. These tumors, which are commonly referred to as fibroids or myomas, affect up to 70% of women before menopause (Stewart et al., 2017), and cause multiple clinical problems, including excessive uterine bleeding, pelvic pain, infertility and miscarriage. Fibroids are also the leading cause of hysterectomy, and therefore incur substantial medical costs around the world (Soliman et al., 2015).

Even in this cancer genomics era, our knowledge of the processes underlying the growth of these tumors is limited. Although they derive from a single smooth muscle cell in the uterine wall, fibroids are a heterogeneous tumor group. Broadly, they can be divided into three main groups, the most common being fibroids with alterations in MED12 (a gene found on the X chromosome), followed by lesions that are deficient in an enzyme called fumarate hydratase, and lesions in which a protein called HMGA2 is overexpressed. Numerous risk factors contribute to fibroid development, including age, obesity, hormones, hypertension and ancestry, with African-American women at the top of the prevalence chart (Dvorská et al., 2017).

Multiple fibroids can occur as part of a hereditary syndrome called familial leiomyomatosis that is caused by germline mutations in the gene coding for fumarate hydratase. Moreover, and independent of familial leiomyomatosis, the first-degree relatives of affected women have a 2.5 times greater chance of developing fibroids than the general population. So, are there loci in the genome that are linked to a general susceptibility to leiomyomas? A handful of studies have explored this question and identified a few loci that are associated with a small increase in the risk of leiomyomas. However, these results are restricted to specific populations or have limited biological significance (Cha et al., 2011; Edwards et al., 2013; Hellwege et al., 2017). Now, in eLife, Lauri Aaltonen of the University of Helsinki and co-workers – including Niko Välimäki and Heli Kuisma as joint first authors – report how they have performed a genome-wide association study that sheds light on the development of leiomyomas (Välimäki et al., 2018).

The study was performed in five stages (Figure 1). In the first phase the researchers analyzed cases and controls from the UK Biobank, which resulted in a list of 50 candidate variants from 22 loci in the genome that might be associated with fibroids. These variants were single nucleotide polymorphisms (SNPs): that is, variants in which one nucleotide had been replaced by a different nucleotide. The researchers then computed the polygenic genomic risk score (GRS), which is a sum of the risk associated to each candidate variant, and validated the association of this score to fibroid development in an independent cohort from Helsinki. In the second phase the researchers combined the two populations in a case-control meta-analysis, which led to another seven candidate variants being identified (bringing the total to 57). The genomic risk score was recalculated in this second phase.

A genome-wide association study on leiomyomas.

The study performed by Välimäki et al. had five phases. (1) Discovery: regions along the genome were compared in 15,453 women with leiomyomas (the cases) and a control sample of 392,628 women from the UK Biobank to identify locations where the two sets of genomes differed. 50 variants (in the form of a SNPs; single nucleotide polymorphisms) at 22 loci were identified and used to compute a genomic risk score (GRS). (2) Meta-analysis: Välimäki et al. then validated the risk score obtained in the first phase by studying cases and controls from Helsinki. Combining the UK Biobank and Helsinki samples lead to the identification of another seven variants. The risk was also recalculated at the end of this phase. (3) Validation: The risk score was then validated in six different ethnic groups (Northern Finland, Caribbean, African, Irish, Indian and 'other white'). When analyzed individually, only one SNP out of the total of 57 gave a significant signal in three different cohorts (UK Biobank, Helsinki and Northern Finland). (4) Biological significance: The genes located within the loci associated with leiomyomas can be divided into two groups: genes related to genome integrity and genes commonly altered in cancers (left), and genes related to the development of the genitourinary system and genes involved in hormonal processes (right). (5) Phenotypic associations: The risk score was associated with the number of leiomyomas per patient, with a higher prevalence of leiomyomas in African women, and with the presence of MED12 mutations in the tumors. A SNP located within the TERT gene was associated with shorter telomeres.

In the third phase Välimäki et al. validated their candidates in six different ethnic groups. The genomic risk score was linked to the development of uterine leiomyomas in all six groups, confirming the involvement of these loci in the disease, independent of ethnicity. The researchers also identified an elevated risk in Black Africans compared to Irish (Caucasian), which could explain the differences in the leiomyoma prevalence.

In the fourth phase Välimäki et al. explored the biological significance of the loci they had identified by searching for genes located within these loci. The genes they identified could be classified into two groups: i) genes implicated in preserving the integrity of the genome and genes commonly implicated in cancer (such as TP53, which is sometimes called the 'guardian of the genome', and TERT, which helps to ensure that the ends of chromosomes remain intact during cell divisions); ii) genes related to the development of the genitourinary system and genes involved in hormonal processes, some of which are also deregulated in cancers (such as the gene for the estrogen receptor gene ESR1).

In the fifth phase the researchers showed that the risk score was positively associated with the number of leiomyomas per patient and a higher prevalence of the disease in African women compared to women with Caucasian ancestry. Overall the results are consistent with a recent study, by a different group, that compared cases and controls from the Icelandic population and Caucasians from the UK Biobank (Rafnar et al., 2018). Surprisingly, the results of Välimäki et al. showed that a variant called rs5937008, which is located in chromosome X near the MED12 gene locus, predisposed women to the subtype of leiomyomas in which the MED12 gene is altered. This result suggests that carriers of the rs5937008 variant have an increased risk of developing MED12-mutated leiomyomas in particular.

The study by Aaltonen, Välimäki, Kuisma and co-workers – who are based the University of Helsinki, Helsinki University Hospital, the University of Oulu and the Karolinska Institute – raises questions that are specific to leiomyomas. For example, what mechanisms orchestrate the process between susceptibility and tumorigenesis? There are also more general questions that remain unanswered: why, for example, do benign tumors – which can carry mutations similar to those seen in malignant tumors – remain benign (Crago et al., 2015; Agnihotri et al., 2016)? Answering such questions will benefit both research into leiomyomas and cancer research in general.

References

Article and author information

Author details

  1. Barbara Rivera

    Barbara Rivera is in the Gerald Bronfman Department of Oncology, McGill University, Montreal, Canada, and the Lady Davis Institute for Medical Research, Jewish General Hospital, Montreal, Canada

    For correspondence
    barbara.riverapolo@mcgill.ca
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9434-6288

Publication history

  1. Version of Record published: October 26, 2018 (version 1)

Copyright

© 2018, Rivera

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,007
    Page views
  • 77
    Downloads
  • 0
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Barbara Rivera
(2018)
Genome-wide Association Study: Not so benign
eLife 7:e42181.
https://doi.org/10.7554/eLife.42181

Further reading

    1. Genetics and Genomics
    2. Immunology and Inflammation
    Jean-David Larouche, Céline M Laumont ... Claude Perreault
    Research Article

    Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

    1. Genetics and Genomics
    Pianpian Zhao, Zhifeng Sheng ... Hou-Feng Zheng
    Research Article

    The ‘diabetic bone paradox’ suggested that type 2 diabetes (T2D) patients would have higher areal bone mineral density (BMD) but higher fracture risk than individuals without T2D. In this study, we found that the genetically predicted T2D was associated with higher BMD and lower risk of fracture in both weighted genetic risk score (wGRS) and two-sample Mendelian randomization (MR) analyses. We also identified ten genomic loci shared between T2D and fracture, with the top signal at SNP rs4580892 in the intron of gene RSPO3. And the higher expression in adipose subcutaneous and higher protein level in plasma of RSPO3 were associated with increased risk of T2D, but decreased risk of fracture. In the prospective study, T2D was observed to be associated with higher risk of fracture, but BMI mediated 30.2% of the protective effect. However, when stratified by the T2D-related risk factors for fracture, we observed that the effect of T2D on the risk of fracture decreased when the number of T2D-related risk factors decreased, and the association became non-significant if the T2D patients carried none of the risk factors. In conclusion, the genetically determined T2D might not be associated with higher risk of fracture. And the shared genetic architecture between T2D and fracture suggested a top signal around RSPO3 gene. The observed effect size of T2D on fracture risk decreased if the T2D-related risk factors could be eliminated. Therefore, it is important to manage the complications of T2D to prevent the risk of fracture.