Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma.
Animals
Cell Line, Tumor
Disease Models, Animal
Gene Expression Regulation
Genes, Reporter
/ genetics
Genetic Predisposition to Disease
/ genetics
Genome-Wide Association Study
/ methods
HEK293 Cells
Humans
Melanocytes
/ metabolism
Melanoma
/ genetics
Mutation
Myxovirus Resistance Proteins
/ genetics
Polymorphism, Single Nucleotide
Proto-Oncogene Proteins B-raf
/ genetics
Quantitative Trait Loci
/ genetics
Zebrafish
/ genetics
Journal
Nature communications
ISSN: 2041-1723
Titre abrégé: Nat Commun
Pays: England
ID NLM: 101528555
Informations de publication
Date de publication:
01 06 2020
01 06 2020
Historique:
received:
06
06
2019
accepted:
12
05
2020
entrez:
3
6
2020
pubmed:
3
6
2020
medline:
18
8
2020
Statut:
epublish
Résumé
Genome-wide association studies (GWAS) have identified ~20 melanoma susceptibility loci, most of which are not functionally characterized. Here we report an approach integrating massively-parallel reporter assays (MPRA) with cell-type-specific epigenome and expression quantitative trait loci (eQTL) to identify susceptibility genes/variants from multiple GWAS loci. From 832 high-LD variants, we identify 39 candidate functional variants from 14 loci displaying allelic transcriptional activity, a subset of which corroborates four colocalizing melanocyte cis-eQTL genes. Among these, we further characterize the locus encompassing the HIV-1 restriction gene, MX2 (Chr21q22.3), and validate a functional intronic variant, rs398206. rs398206 mediates the binding of the transcription factor, YY1, to increase MX2 levels, consistent with the cis-eQTL of MX2 in primary human melanocytes. Melanocyte-specific expression of human MX2 in a zebrafish model demonstrates accelerated melanoma formation in a BRAF
Identifiants
pubmed: 32483191
doi: 10.1038/s41467-020-16590-1
pii: 10.1038/s41467-020-16590-1
pmc: PMC7264232
doi:
Substances chimiques
MX2 protein, human
0
Myxovirus Resistance Proteins
0
Proto-Oncogene Proteins B-raf
EC 2.7.11.1
Types de publication
Journal Article
Research Support, N.I.H., Extramural
Research Support, N.I.H., Intramural
Research Support, Non-U.S. Gov't
Langues
eng
Sous-ensembles de citation
IM
Pagination
2718Subventions
Organisme : NCI NIH HHS
ID : R01 CA083115
Pays : United States
Organisme : NCI NIH HHS
ID : R01 CA103846
Pays : United States
Organisme : NIAAA NIH HHS
ID : R01 AA010384
Pays : United States
Organisme : Cancer Research UK
ID : C588/A19167
Pays : United Kingdom
Organisme : NHGRI NIH HHS
ID : R01 HG009120
Pays : United States
Organisme : Howard Hughes Medical Institute
Pays : United States
Organisme : NCI NIH HHS
ID : P01 CA163222
Pays : United States
Références
Law, M. H. et al. Genome-wide meta-analysis identifies five new susceptibility loci for cutaneous malignant melanoma. Nat. Genet. 47, 987–995 (2015).
pubmed: 26237428
pmcid: 4557485
doi: 10.1038/ng.3373
Gudbjartsson, D. F. et al. ASIP and TYR pigmentation variants associate with cutaneous melanoma and basal cell carcinoma. Nat. Genet. 40, 886–891 (2008).
pubmed: 18488027
doi: 10.1038/ng.161
Brown, K. M. et al. Common sequence variants on 20q11.22 confer melanoma susceptibility. Nat. Genet. 40, 838–840 (2008).
pubmed: 18488026
pmcid: 2755512
doi: 10.1038/ng.163
Bishop, D. T. et al. Genome-wide association study identifies three loci associated with melanoma risk. Nat. Genet. 41, 920–925 (2009).
pubmed: 19578364
pmcid: 2741419
doi: 10.1038/ng.411
Falchi, M. et al. Genome-wide association study identifies variants at 9p21 and 22q13 associated with development of cutaneous nevi. Nat. Genet. 41, 915–919 (2009).
pubmed: 19578365
doi: 10.1038/ng.410
Iles, M. M. et al. A variant in FTO shows association with melanoma risk not due to BMI. Nat. Genet. 45, 428–32, 432e1 (2013).
pubmed: 23455637
pmcid: 3640814
doi: 10.1038/ng.2571
Rafnar, T. et al. Sequence variants at the TERT-CLPTM1L locus associate with many cancer types. Nat. Genet. 41, 221–227 (2009).
pubmed: 19151717
pmcid: 4525478
doi: 10.1038/ng.296
Macgregor, S. et al. Genome-wide association study identifies a new melanoma susceptibility locus at 1q21.3. Nat. Genet. 43, 1114–1118 (2011).
pubmed: 21983785
pmcid: 3227560
doi: 10.1038/ng.958
Barrett, J. H. et al. Genome-wide association study identifies three new melanoma susceptibility loci. Nat. Genet. 43, 1108–1113 (2011).
pubmed: 21983787
pmcid: 3251256
doi: 10.1038/ng.959
Ransohoff, K. J. et al. Two-stage genome-wide association study identifies a novel susceptibility locus associated with melanoma. Oncotarget 8, 17586–17592 (2017).
pubmed: 28212542
pmcid: 5392271
doi: 10.18632/oncotarget.15230
Duffy, D. L. et al. Multiple pigmentation gene polymorphisms account for a substantial proportion of risk of cutaneous malignant melanoma. J. Invest. Dermatol. 130, 520–528 (2010).
pubmed: 19710684
doi: 10.1038/jid.2009.258
Beaumont, K. A. et al. Receptor function, dominant negative activity and phenotype correlations for MC1R variant alleles. Hum. Mol. Genet. 16, 2249–2260 (2007).
pubmed: 17616515
doi: 10.1093/hmg/ddm177
Tsetskhladze, Z. R. et al. Functional assessment of human coding mutations affecting skin pigmentation using zebrafish. PLoS ONE 7, e47398 (2012).
pubmed: 23071798
pmcid: 3468441
doi: 10.1371/journal.pone.0047398
Halaban, R. et al. Endoplasmic reticulum retention is a common defect associated with tyrosinase-negative albinism. Proc. Natl Acad. Sci. USA 97, 5889–5894 (2000).
pubmed: 10823941
doi: 10.1073/pnas.97.11.5889
Liu, F. et al. Genetics of skin color variation in Europeans: genome-wide association studies with functional follow-up. Hum. Genet. 134, 823–835 (2015).
pubmed: 25963972
pmcid: 4495261
doi: 10.1007/s00439-015-1559-0
Duffy, D. L. et al. Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways. Nat. Commun. 9, 4774 (2018).
pubmed: 30429480
pmcid: 6235897
doi: 10.1038/s41467-018-06649-5
Duffy, D. L. et al. IRF4 variants have age-specific effects on nevus count and predispose to melanoma. Am. J. Hum. Genet. 87, 6–16 (2010).
pubmed: 20602913
pmcid: 2896771
doi: 10.1016/j.ajhg.2010.05.017
Choi, J. et al. A common intronic variant of PARP1 confers melanoma risk and mediates melanocyte growth via regulation of MITF. Nat. Genet. 49, 1326–1335 (2017).
pubmed: 28759004
doi: 10.1038/ng.3927
Fang, J. et al. Functional characterization of a multi-cancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148. Nat. Commun. 8, 15034 (2017).
pubmed: 28447668
pmcid: 5414179
doi: 10.1038/ncomms15034
Ulirsch, J. C. et al. Systematic functional dissection of common genetic variation affecting red blood cell traits. Cell 165, 1530–1545 (2016).
pubmed: 27259154
pmcid: 4893171
doi: 10.1016/j.cell.2016.04.048
Kheradpour, P. et al. Systematic dissection of regulatory motifs in 2000 predicted human enhancers using a massively parallel reporter assay. Genome Res. 23, 800–811 (2013).
pubmed: 23512712
pmcid: 3638136
doi: 10.1101/gr.144899.112
Melnikov, A. et al. Systematic dissection and optimization of inducible enhancers in human cells using a massively parallel reporter assay. Nat. Biotechnol. 30, 271–277 (2012).
pubmed: 22371084
pmcid: 3297981
doi: 10.1038/nbt.2137
GTEx Consortium et al. Genetic effects on gene expression across human tissues. Nature 550, 204–213 (2017).
pubmed: 5776756
pmcid: 5776756
doi: 10.1038/nature24277
Zhang, T. et al. Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes. Genome Res. 28, 1621–1635 (2018).
pubmed: 30333196
pmcid: 6211648
doi: 10.1101/gr.233304.117
Hoggart, C. J., Whittaker, J. C., De Iorio, M. & Balding, D. J. Simultaneous analysis of all SNPs in genome-wide and re-sequencing association studies. PLoS Genet. 4, e1000130 (2008).
pubmed: 18654633
pmcid: 2464715
doi: 10.1371/journal.pgen.1000130
Verfaillie, A. et al. Decoding the regulatory landscape of melanoma reveals TEADS as regulators of the invasive cell state. Nat. Commun. 6, 6683 (2015).
pubmed: 25865119
pmcid: 4403341
doi: 10.1038/ncomms7683
Melnikov, A., Zhang, X., Rogov, P., Wang, L. & Mikkelsen, T. S. Massively parallel reporter assays in cultured mammalian cells. J. Vis. Exp. 90, e51719 (2014).
Coetzee, S. G., Coetzee, G. A. & Hazelett, D. J. motifbreakR: an R/Bioconductor package for predicting variant effects at transcription factor binding sites. Bioinformatics 31, 3847–3849 (2015).
pubmed: 26272984
pmcid: 4653394
Kichaev, G. et al. Improved methods for multi-trait fine mapping of pleiotropic risk loci. Bioinformatics 33, 248–255 (2017).
pubmed: 27663501
doi: 10.1093/bioinformatics/btw615
Li, Y. I. et al. Annotation-free quantification of RNA splicing using LeafCutter. Nat. Genet. 50, 151–158 (2018).
pubmed: 29229983
doi: 10.1038/s41588-017-0004-9
Baylin, S. B. & Jones, P. A. A decade of exploring the cancer epigenome-biological and translational implications. Nat. Rev. Cancer 11, 726–734 (2011).
pubmed: 21941284
pmcid: 3307543
doi: 10.1038/nrc3130
Kulis, M. et al. Epigenomic analysis detects widespread gene-body DNA hypomethylation in chronic lymphocytic leukemia. Nat. Genet. 44, 1236–1242 (2012).
pubmed: 23064414
doi: 10.1038/ng.2443
Maunakea, A. K. et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature 466, 253–257 (2010).
pubmed: 20613842
pmcid: 3998662
doi: 10.1038/nature09165
Varley, K. E. et al. Dynamic DNA methylation across diverse human cell lines and tissues. Genome Res. 23, 555–567 (2013).
pubmed: 23325432
pmcid: 3589544
doi: 10.1101/gr.147942.112
Wolff, E. M. et al. Hypomethylation of a LINE-1 promoter activates an alternate transcript of the MET oncogene in bladders with cancer. PLoS Genet. 6, e1000917 (2010).
pubmed: 20421991
pmcid: 2858672
doi: 10.1371/journal.pgen.1000917
Gordon, S., Akopyan, G., Garban, H. & Bonavida, B. Transcription factor YY1: structure, function, and therapeutic implications in cancer biology. Oncogene 25, 1125–1142 (2006).
pubmed: 16314846
doi: 10.1038/sj.onc.1209080
Li, J. et al. YY1 regulates melanocyte development and function by cooperating with MITF. PLoS Genet. 8, e1002688 (2012).
pubmed: 22570637
pmcid: 3342948
doi: 10.1371/journal.pgen.1002688
Yeo, N. C. et al. An enhanced CRISPR repressor for targeted mammalian gene regulation. Nat. Methods 15, 611–616 (2018).
pubmed: 30013045
pmcid: 6129399
doi: 10.1038/s41592-018-0048-5
Weintraub, A. S. et al. YY1 is a structural regulator of enhancer-promoter loops. Cell 171, 1573–1588 e28 (2017).
pubmed: 29224777
pmcid: 5785279
doi: 10.1016/j.cell.2017.11.008
Goujon, C. et al. Human MX2 is an interferon-induced post-entry inhibitor of HIV-1 infection. Nature 502, 559–562 (2013).
pubmed: 24048477
doi: 10.1038/nature12542
Kane, M. et al. MX2 is an interferon-induced inhibitor of HIV-1 infection. Nature 502, 563–566 (2013).
pubmed: 24121441
pmcid: 3912734
doi: 10.1038/nature12653
King, M. C., Raposo, G. & Lemmon, M. A. Inhibition of nuclear import and cell-cycle progression by mutated forms of the dynamin-like GTPase MxB. Proc. Natl Acad. Sci. USA 101, 8957–8962 (2004).
pubmed: 15184662
doi: 10.1073/pnas.0403167101
Ceol, C. J. et al. The histone methyltransferase SETDB1 is recurrently amplified in melanoma and accelerates its onset. Nature 471, 513–517 (2011).
pubmed: 21430779
pmcid: 3348545
doi: 10.1038/nature09806
Chundru, V. K. et al. Examining the impact of imputation errors on fine-mapping using DNA methylation QTL as a model trait. Genetics 212, 577–586 (2019).
pubmed: 31040117
pmcid: 6614908
doi: 10.1534/genetics.118.301861
Schaid, D. J., Chen, W. & Larson, N. B. From genome-wide associations to candidate causal variants by statistical fine-mapping. Nat. Rev. Genet. 19, 491–504 (2018).
pubmed: 29844615
pmcid: 6050137
doi: 10.1038/s41576-018-0016-z
Amit, I. et al. A module of negative feedback regulators defines growth factor signaling. Nat. Genet. 39, 503–512 (2007).
pubmed: 17322878
doi: 10.1038/ng1987
Katsuoka, F. & Yamamoto, M. Small Maf proteins (MafF, MafG, MafK): history, structure and function. Gene 586, 197–205 (2016).
pubmed: 27058431
pmcid: 4911266
doi: 10.1016/j.gene.2016.03.058
Denat, L., Kadekaro, A. L., Marrot, L., Leachman, S. A. & Abdel-Malek, Z. A. Melanocytes as instigators and victims of oxidative stress. J. Invest. Dermatol. 134, 1512–1518 (2014).
pubmed: 24573173
pmcid: 4418514
doi: 10.1038/jid.2014.65
Olson, O. C. & Joyce, J. A. Cysteine cathepsin proteases: regulators of cancer progression and therapeutic response. Nat. Rev. Cancer 15, 712–729 (2015).
pubmed: 26597527
doi: 10.1038/nrc4027
McIlwain, D. R., Berger, T. & Mak, T. W. Caspase functions in cell death and disease. Cold Spring Harb. Perspect. Biol. 5, a008656 (2013).
pubmed: 23545416
pmcid: 3683896
doi: 10.1101/cshperspect.a008656
Michailidou, K. et al. Association analysis identifies 65 new breast cancer risk loci. Nature 551, 92–94 (2017).
pubmed: 29059683
pmcid: 5798588
doi: 10.1038/nature24284
Stacey, S. N. et al. New basal cell carcinoma susceptibility loci. Nat. Commun. 6, 6825 (2015).
pubmed: 25855136
pmcid: 4403348
doi: 10.1038/ncomms7825
Juraleviciute, M. et al. MX 2 is a novel regulator of cell cycle in melanoma cells. Pigment Cell Melanoma Res. 33, 446–457 (2019).
doi: 10.1111/pcmr.12837
Barrett, J. H. et al. Fine mapping of genetic susceptibility loci for melanoma reveals a mixture of single variant and multiple variant regions. Int. J. Cancer 136, 1351–1360 (2015).
pubmed: 25077817
doi: 10.1002/ijc.29099
Stephens, M. & Balding, D. J. Bayesian statistical methods for genetic association studies. Nat. Rev. Genet. 10, 681–690 (2009).
pubmed: 19763151
doi: 10.1038/nrg2615
Vignal, C. M., Bansal, A. T. & Balding, D. J. Using penalised logistic regression to fine map HLA variants for rheumatoid arthritis. Ann. Hum. Genet. 75, 655–664 (2011).
pubmed: 21883133
doi: 10.1111/j.1469-1809.2011.00670.x
Long, J. S. & Ervin, L. H. Using heteroscedasticity consistent standard errors in the linear regression model. Am. Statistician 54, 217–224 (2000).
Kulakovskiy, I. V. et al. HOCOMOCO: towards a complete collection of transcription factor binding models for human and mouse via large-scale ChIP-Seq analysis. Nucleic Acids Res. 46, D252–D259 (2018).
pubmed: 29140464
doi: 10.1093/nar/gkx1106
Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinforma. 12, 323 (2011).
doi: 10.1186/1471-2105-12-323
Ongen, H., Buil, A., Brown, A. A., Dermitzakis, E. T. & Delaneau, O. Fast and efficient QTL mapper for thousands of molecular phenotypes. Bioinformatics 32, 1479–1485 (2016).
pubmed: 26708335
doi: 10.1093/bioinformatics/btv722
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
pubmed: 23104886
pmcid: 23104886
doi: 10.1093/bioinformatics/bts635
Assenov, Y. et al. Comprehensive analysis of DNA methylation data with RnBeads. Nat. Methods 11, 1138–1140 (2014).
pubmed: 25262207
pmcid: 4216143
doi: 10.1038/nmeth.3115
Aryee, M. J. et al. Minfi: a flexible and comprehensive bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 30, 1363–1369 (2014).
pubmed: 24478339
pmcid: 4016708
doi: 10.1093/bioinformatics/btu049
Spruijt, C. G., Baymaz, H. I. & Vermeulen, M. Identifying specific protein-DNA interactions using SILAC-based quantitative proteomics. Methods Mol. Biol. 977, 137–157 (2013).
pubmed: 23436359
doi: 10.1007/978-1-62703-284-1_11
Makowski, M. M. et al. An interaction proteomics survey of transcription factor binding at recurrent TERT promoter mutations. Proteomics 16, 417–426 (2016).
pubmed: 26553150
doi: 10.1002/pmic.201500327
Cox, J. & Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 26, 1367–1372 (2008).
pubmed: 19029910
pmcid: 19029910
doi: 10.1038/nbt.1511
Chari, R., Yeo, N. C., Chavez, A. & Church, G. M. sgRNA Scorer 2.0: a species-independent model to predict CRISPR/Cas9 activity. ACS Synth. Biol. 6, 902–904 (2017).
pubmed: 28146356
pmcid: 5793212
doi: 10.1021/acssynbio.6b00343
Lagonigro, M. S. et al. CTAB-urea method purifies RNA from melanin for cDNA microarray analysis. Pigment Cell Res. 17, 312–315 (2004).
pubmed: 15140079
doi: 10.1111/j.1600-0749.2004.00155.x
Liao, Y., Smyth, G. K. & Shi, W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 47, e47 (2019).
pubmed: 30783653
pmcid: 6486549
doi: 10.1093/nar/gkz114
Varet, H., Brillet-Gueguen, L., Coppee, J. Y. & Dillies, M. A. SARTools: a DESeq2- and EdgeR-Based R pipeline for comprehensive differential analysis of RNA-Seq Data. PLoS ONE 11, e0157022 (2016).
pubmed: 27280887
pmcid: 4900645
doi: 10.1371/journal.pone.0157022
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140 (2010).
pubmed: 19910308
doi: 10.1093/bioinformatics/btp616
Ewels, P., Magnusson, M., Lundin, S. & Kaller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–3048 (2016).
pubmed: 27312411
pmcid: 5039924
doi: 10.1093/bioinformatics/btw354
Patro, R., Duggal, G., Love, M. I., Irizarry, R. A. & Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419 (2017).
pubmed: 28263959
pmcid: 28263959
doi: 10.1038/nmeth.4197
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
pubmed: 25516281
pmcid: 25516281
doi: 10.1186/s13059-014-0550-8