Two-stage genome-wide association study identifies a novel susceptibility locus associated with melanoma

A blood- and brain-based EWAS of smoking

DNA methylation offers an objective method to assess the impact of smoking. In this work, we conduct a Bayesian EWAS of smoking pack years (n = 17,865, ~850k sites, Illumina EPIC array) and extend it by analysing whole genome data of smokers and non-smokers from Generation Scotland (n = 46, ~4-21 million sites via TWIST and Oxford Nanopore sequencing). We develop mCigarette, an epigenetic biomarker of smoking, and test it in two British cohorts. Results of brain- and blood-based EWAS (nbrain=14, nblood = 882, >450k sites, Illumina arrays) reveal several loci with near-perfect discrimination of smoking status, but which do not overlap across tissues. Furthermore, we perform a GWAS of epigenetic smoking, identifying several smoking-related loci. Overall, we improve smoking-related biomarker accuracy and enhance the understanding of the effects of smoking by integrating DNA methylation data from multiple tissues and cohorts.

Knowledge graph and frontier trends in melanoma-associated ncRNAs: a bibliometric analysis from 2006 to 2023

Objectives

Malignant melanoma (MM) is a highly malignant skin tumor. Although research on non-coding RNAs (ncRNAs) of MM has advanced swiftly in recent years, no specific bibliometric analyses have been conducted on this topic. The present study aims to summarize the knowledge graphs and frontier trends and to provide new perspectives and direction of collaboration for researchers.

Method

Research data on melanoma and ncRNA published from January 1, 2006 to October 9, 2023 were retrieved and extracted from Web of Science. R Studio (Version 4.3.1), Scimago Graphica (Version 1.0.36), VOSviewer version (1.6.19), and Citespace (6.2.4R) were used to analyze the publications, countries, journals, institutions, authors, keywords, references, and other relevant data and to build collaboration network graphs and co-occurrence network graphs accordingly.

Results

A total of 1,222 articles were retrieved, involving 4,894 authors, 385 journals, 43,220 references, 2413 keywords, and 1,651 institutions in 47 countries. The average annual growth rate in the number of articles was 25.02% from 2006 to 2023; among all the journals, Plos One had the highest number of publications and citations, which are 42 publications and 2,228 citations, respectively. Chinese researchers were the most prolific publishers in this field, having published a total of 657 articles, among which 42 were published by Shanghai Jiao Tong University, which was the most productive institution. In recent years, the most explored keywords included long non-coding RNAs, immunotherapy, and exosm. According to the timeline chart of reference co-citation, "functional role" has been the most explored hotspot since 2015, and human cancer is a newly emerged hotspot after 2021.

Conclusion

Through a bibliometric analysis, this study included all publications on ncRNAs and melanoma that were published in English from 2006 to 2023 in Web of Science to analyze the trends in the number of publications, international research focuses, and the direction of collaboration. The results of this study may provide information on knowledge graph, frontier trends and identify research topics in melanoma. More current research proved that ncRNA plays a crucial role in the biological behavior of melanoma including proliferation, invasion, metastasis, drug resistance, etc. With the development of research on ncRNA and melanoma, ncRNA may great potential in development of early diagnosis, targeted therapy and efficacy evaluation in the future. The results of this study also provide new perspectives and research partners for researchers in this field.

Mapping chromatin interactions at melanoma susceptibility loci and cell-type specific dataset integration uncovers distant gene targets of cis-regulation

Genome-wide association studies (GWAS) of melanoma risk have identified 68 independent signals at 54 loci. For most loci, specific functional variants and their respective target genes remain to be established. Capture-HiC is an assay that links fine-mapped risk variants to candidate target genes by comprehensively mapping cell-type specific chromatin interactions. We performed a melanoma GWAS region-focused capture-HiC assay in human primary melanocytes to identify physical interactions between fine-mapped risk variants and potential causal melanoma susceptibility genes. Overall, chromatin interaction data alone nominated potential causal genes for 61 of the 68 melanoma risk signals, identifying many candidates beyond those reported by previous studies. We further integrated these data with cell-type specific epigenomic (chromatin state, accessibility), gene expression (eQTL/TWAS), DNA methylation (meQTL/MWAS), and massively parallel reporter assay (MPRA) data to prioritize potentially cis-regulatory variants and their respective candidate gene targets. From the set of fine-mapped variants across these loci, we identified 140 prioritized candidate causal variants linked to 195 candidate genes at 42 risk signals. In addition, we developed an integrative scoring system to facilitate candidate gene prioritization, integrating melanocyte and melanoma datasets. Notably, at several GWAS risk signals we observed long-range chromatin connections (500 kb to >1 Mb) with distant candidate target genes. We validated several such cis-regulatory interactions using CRISPR inhibition, providing evidence for known cancer driver genes MDM4 and CBL, as well as the SRY-box transcription factor SOX4, as likely melanoma risk genes.

Identification of genomic-wide genetic links between cutaneous melanoma and obesity-related physical traits via cFDR

BACKGROUND

Both epidemiological and clinical studies have suggested the comorbidity between cutaneous melanoma (CM) and obesity-related physical traits. However, it remains unclear about their shared genetic architecture.

OBJECTIVE

To determine the shared genetic architecture between CM and obesity-related physical traits through conditional false discovery rate (cFDR) analysis.

METHOD

Quantile-quantile plots were firstly built to assess the pleiotropic enrichment of shared single nucleotide polymorphisms between CM and each trait. Then, cFDR and conjunctional cFDR (ccFDR) were used to identify the shared risk loci between CM and each trait. Moreover, the functional evaluation of shared risk genes was carried out through analyses of expression quantitative trait loci (eQTL), Kyoto Encyclopedia of Genes and Genomes and gene ontology, respectively. Finally, single-cell sequence analysis was performed to locate the expression of eQTL-mapped genes in tissues.

RESULTS

Successive pleiotropic enrichment was found between CM and 5 obesity-related traits or height. 24 shared risk loci were identified between CM and 13 traits except appendicular lean mass using ccFDR analysis, with 17 novel and 4 validated loci. The functions of ccFDR-identified and eQTL-mapped genes were revealed to be mainly involved in cellular senescence, proliferation, meiotic nuclear division, cell cycle, and the metabolism of lipid, cholesterol and glucose. Single-cell sequence analysis showed that keratinocytes contribute to the occurrence and aggressiveness of CM through secreting paracrine cytokines.

CONCLUSION

Our findings demonstrate the significant genetic correlation between CM and obesity-related physical traits, which may provide a novel genetical basis for the pathogenesis and treatment of CM.

Regulatory miRNAs and lncRNAs in Skin Cancer: A Narrative Review

Non-coding RNAs (ncRNAs) have a significant regulatory role in the pathogenesis of skin cancer, despite the fact that protein-coding genes have generally been the focus of research efforts in the field. We comment on the actions of long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) in the current review with an eye toward potential therapeutic treatments. LncRNAs are remarkably adaptable, acting as scaffolding, guides, or decoys to modify key signaling pathways (i.e., the Wnt/β-catenin pathway) and gene expression. As post-transcriptional gatekeepers, miRNAs control gene expression by attaching to messenger RNAs and causing their degradation or suppression during translation. Cell cycle regulation, cellular differentiation, and immunological responses are all affected by the dysregulation of miRNAs observed in skin cancer. NcRNAs also show promise as diagnostic biomarkers and prognostic indicators. Unraveling the complexity of the regulatory networks governed by ncRNAs in skin cancer offers unprecedented opportunities for groundbreaking targeted therapies, revolutionizing the landscape of dermatologic care.

Patterns of item nonresponse behaviour to survey questionnaires are systematic and associated with genetic loci

Response to survey questionnaires is vital for social and behavioural research, and most analyses assume full and accurate response by participants. However, nonresponse is common and impedes proper interpretation and generalizability of results. We examined item nonresponse behaviour across 109 questionnaire items in the UK Biobank (N = 360,628). Phenotypic factor scores for two participant-selected nonresponse answers, 'Prefer not to answer' (PNA) and 'I don't know' (IDK), each predicted participant nonresponse in follow-up surveys (incremental pseudo-R2 = 0.056), even when controlling for education and self-reported health (incremental pseudo-R2 = 0.046). After performing genome-wide association studies of our factors, PNA and IDK were highly genetically correlated with one another (rg = 0.73 (s.e. = 0.03)) and with education (rg,PNA = -0.51 (s.e. = 0.03); rg,IDK = -0.38 (s.e. = 0.02)), health (rg,PNA = 0.51 (s.e. = 0.03); rg,IDK = 0.49 (s.e. = 0.02)) and income (rg,PNA = -0.57 (s.e. = 0.04); rg,IDK = -0.46 (s.e. = 0.02)), with additional unique genetic associations observed for both PNA and IDK (P < 5 × 10-8). We discuss how these associations may bias studies of traits correlated with item nonresponse and demonstrate how this bias may substantially affect genome-wide association studies. While the UK Biobank data are deidentified, we further protected participant privacy by avoiding exploring non-response behaviour to single questions, assuring that no information can be used to associate results with any particular respondents.

Explainable multi-task learning improves the parallel estimation of polygenic risk scores for many diseases through shared genetic basis

Many complex diseases share common genetic determinants and are comorbid in a population. We hypothesized that the co-occurrences of diseases and their overlapping genetic etiology can be exploited to simultaneously improve multiple diseases' polygenic risk scores (PRS). This hypothesis was tested using a multi-task learning (MTL) approach based on an explainable neural network architecture. We found that parallel estimations of the PRS for 17 prevalent cancers in a pan-cancer MTL model were generally more accurate than independent estimations for individual cancers in comparable single-task learning (STL) models. Such performance improvement conferred by positive transfer learning was also observed consistently for 60 prevalent non-cancer diseases in a pan-disease MTL model. Interpretation of the MTL models revealed significant genetic correlations between the important sets of single nucleotide polymorphisms used by the neural network for PRS estimation. This suggested a well-connected network of diseases with shared genetic basis.

BASP1 promotes high glucose-induced endothelial apoptosis in diabetes via activation of EGFR signaling

AIMS

Diabetes mellitus is a common chronic disease of glucose metabolism. Endothelial dysfunction is an early event in diabetes complicated by cardiovascular disease. This study aimed to reveal the expression of BASP1 and its biological roles in endothelial cell dysfunction in diabetes complicated by cardiovascular disease.

MATERIALS AND METHODS

By analyzing the databases related to diabetes complicated with coronary heart disease, BASP1 was screened out as an upregulated gene. Human umbilical vein endothelial cells (HUVECs) and primary mouse aortic endothelial cells were treated with high glucose to establish cell models of diabetes-related endothelial dysfunction, and the expression changes of BASP1 were verified by RT-qPCR, western blot, and immunofluorescence. BASP1 was silenced or overexpressed by siRNA or overexpression plasmid, and its effects on cell migration, apoptosis, tube formation, inflammatory response, and ROS were detected. The possible signaling pathway of BASP1 was found and the mechanism of BASP1 on promoting the progression of endothelial dysfunction was explored using the EGFR inhibitor, gefitinib.

RESULTS

Bioinformatics analysis indicated that the expression of BASP1 in patients with diabetes mellitus and concomitant coronary heart disease was increased. High glucose induced the upregulation of BASP1 expression in endothelial cells, and showed a time-dependent relationship. Silencing of BASP1 alleviated the damage of high glucose to endothelial cells. BASP1 regulated EGFR positively. The promoting effect of BASP1 on endothelial cell apoptosis may be achieved by regulating the EGFR pathway.

CONCLUSION

BASP1 promotes endothelial cell injury induced by high glucose in patients with diabetes, which may be activated by activating the EGFR pathway.

A genome-wide association study of germline variation and melanoma prognosis

Background

The high mortality of cutaneous melanoma (CM) is partly due to unpredictable patterns of disease progression in patients with early-stage lesions. The reliable prediction of advanced disease risk from early-stage CM, is an urgent clinical need, especially given the recent expansion of immune checkpoint inhibitor therapy to the adjuvant setting. In our study, we comprehensively investigated the role of germline variants as CM prognostic markers.

Methods

We performed a genome-wide association analysis in two independent cohorts of N=551 (discovery), and N=550 (validation) early-stage immunotherapy-naïve melanoma patients. A multivariable Cox proportional hazard regression model was used to identify associations with overall survival in the discovery group, followed by a validation analysis. Transcriptomic profiling and survival analysis were used to elucidate the biological relevance of candidate genes associated with CM progression.

Results

We found two independent associations of germline variants with melanoma prognosis. The alternate alleles of these two SNPs were both associated with an increased risk of death [rs60970102 in MELK: HR=3.14 (2.05-4.81), p=1.48×10-7; and rs77480547 in SH3BP4: HR=3.02 (2.02-4.52), p=7.58×10-8, both in the pooled cohort]. The addition of the combined risk alleles (CRA) of the identified variants into the prognostic model improved the predictive power, as opposed to a model of clinical covariates alone.

Conclusions

Our study provides suggestive evidence of novel melanoma germline prognostic markers, implicating two candidate genes: an oncogene MELK and a tumor suppressor SH3BP4, both previously suggested to affect CM progression. Pending further validation, these findings suggest that the genetic factors may improve the prognostic stratification of high-risk early-stage CM patients, and propose putative biological insights for potential therapeutic investigation of these targets to prevent aggressive outcome from early-stage melanoma.

The performance of multi-gene panels for breast/ovarian cancer predisposition

BRCA1 and BRCA2 are the most mutated genes in breast cancer. We analyzed 48 breast cancer subjects using two methods that differ in terms of number of genes investigated and strategy used (primers: Panel A - 12 genes - vs probes: Panel B - 48 genes). Both the panels and procedures identified "pathogenic" or "likely pathogenic" variants in TP53, ATM, CHEK2 and BARD1 besides BRCA1 and BRCA2. Panel B identified two other putatively pathogenic variants in RNASEL and in RAD50. Identification of variants other than the BRCA genes can be useful in patient management. A total of 121 variants were distributed within the 12 genes and were correctly detected by both panels. However, the number of calls without divergence, namely ± 0.10 difference of allelic frequency, was 78.3%, while calls with a divergence below 0.10 was 16.7%, thus indicating that only 5% (n = 275) of 5,412 calls had a divergence above 0.10. Although these panels differ from each other, both are useful in different situations, particularly when patients should be tested for genes other than BRCA1/2 (as occurs in patients affected by a so called hereditary syndrome) or for therapeutic purposes.

A multi-phenotype analysis reveals 19 susceptibility loci for basal cell carcinoma and 15 for squamous cell carcinoma

Basal cell carcinoma and squamous cell carcinoma are the most common skin cancers, and have genetic overlap with melanoma, pigmentation traits, autoimmune diseases, and blood biochemistry biomarkers. In this multi-trait genetic analysis of over 300,000 participants from Europe, Australia and the United States, we reveal 78 risk loci for basal cell carcinoma (19 previously unknown and replicated) and 69 for squamous cell carcinoma (15 previously unknown and replicated). The previously unknown risk loci are implicated in cancer development and progression (e.g. CDKL1), pigmentation (e.g. TPCN2), cardiometabolic (e.g. FADS2), and immune-regulatory pathways for innate immunity (e.g. IFIH1), and HIV-1 viral load modulation (e.g. CCR5). We also report an optimised polygenic risk score for effective risk stratification for keratinocyte cancer in the Canadian Longitudinal Study of Aging (794 cases and 18139 controls), which could facilitate skin cancer surveillance e.g. in high risk subpopulations such as transplantees.

Implementation of individualised polygenic risk score analysis: a test case of a family of four

BACKGROUND

Polygenic risk scores (PRS) have been widely applied in research studies, showing how population groups can be stratified into risk categories for many common conditions. As healthcare systems consider applying PRS to keep their populations healthy, little work has been carried out demonstrating their implementation at an individual level.

CASE PRESENTATION

We performed a systematic curation of PRS sources from established data repositories, selecting 15 phenotypes, comprising an excess of 37 million SNPs related to cancer, cardiovascular, metabolic and autoimmune diseases. We tested selected phenotypes using whole genome sequencing data for a family of four related individuals. Individual risk scores were given percentile values based upon reference distributions among 1000 Genomes Iberians, Europeans, or all samples. Over 96 billion allele effects were calculated in order to obtain the PRS for each of the individuals analysed here.

CONCLUSIONS

Our results highlight the need for further standardisation in the way PRS are developed and shared, the importance of individual risk assessment rather than the assumption of inherited averages, and the challenges currently posed when translating PRS into risk metrics.

Genome-Wide Association Study Suggests the Variant rs7551288*A within the DHCR24 Gene Is Associated with Poor Overall Survival in Melanoma Patients

Simple Summary

The aim of this work was to investigate prognostic genetic factors in melanoma patients. Phenotypic and disease data as well as biomaterial were collected after informed consent from patients followed up in a Skin Cancer Center of a University clinic. Genome-wide analysis (GWAS) was performed with survival data of 556 melanoma patients and genetic data including more than 300,000 common polymorphisms. The SNP rs7551288 reached suggestive genome-wide significance (p = 2 × 10⁻⁶). This intronic variant of the DHCR24 gene is involved in the cholesterol synthesis pathway. Further analyses and a literature review suggest an important role of this locus for the clinical course of disease in melanoma patients.

Abstract

Melanoma incidence rates are high among individuals with fair skin and multiple naevi. Established prognostic factors are tumour specific, and less is known about prognostic host factors. A total of 556 stage I to stage IV melanoma patients from Germany with phenotypic and disease-specific data were analysed; 64 of these patients died of melanoma after a median follow-up time of 8 years. Germline DNA was assessed by the HumanCoreExome BeadChip and data of 356,384 common polymorphisms distributed over all 23 chromosomes were used for a genome-wide analysis. A suggestive genome-wide significant association of the intronic allele rs7551288*A with diminished melanoma-specific survival was detected (p = 2 × 10⁻⁶). The frequency of rs7551288*A was 0.43 and was not associated with melanoma risk, hair and eye colour, tanning and total naevus count. Cox regression multivariate analyses revealed a 5.31-fold increased risk of melanoma-specific death for patients with the rs7551288 A/A genotype, independent of tumour thickness, ulceration and stage of disease at diagnoses. The variant rs7551288 belongs to the DHCR24 gene, which encodes Seladin-1, an enzyme involved in the biosynthesis of cholesterol. Further investigations are needed to confirm this genetic variant as a novel prognostic biomarker and to explore whether specific treatment strategies for melanoma patients might be derived from it.

Coding and regulatory variants are associated with serum protein levels and disease

Circulating proteins can be used to diagnose and predict disease-related outcomes. A deep serum proteome survey recently revealed close associations between serum protein networks and common disease. In the current study, 54,469 low-frequency and common exome-array variants were compared to 4782 protein measurements in the serum of 5343 individuals from the AGES Reykjavik cohort. This analysis identifies a large number of serum proteins with genetic signatures overlapping those of many diseases. More specifically, using a study-wide significance threshold, we find that 2021 independent exome array variants are associated with serum levels of 1942 proteins. These variants reside in genetic loci shared by hundreds of complex disease traits, highlighting serum proteins' emerging role as biomarkers and potential causative agents of a wide range of diseases.

Melanoma susceptibility: an update on genetic and epigenetic findings

Malignant melanoma is one of the most highly ranked cancers in terms of years of life lost. Hereditary melanoma with its increased familial susceptibility is thought to affect up to 12% of all melanoma patients. In the past, only a few high-penetrance genes associated with familial melanoma, such as CDKN2A and CDK4, have been clinically tested. However, findings now indicate that melanoma is a cancer most likely to develop not only due to high-penetrance variants but also due to polygenic inheritance patterns, leaving no clear division between the hereditary and sporadic development of malignant melanoma. Various pathogenic low-penetrance variants were recently discovered through genome-wide association studies, and are now translated into polygenic risk scores. These can show superior sensitivity rates for the prediction of melanoma susceptibility and related mixed cancer syndromes than risk scores based on phenotypic traits of the patients, with odds ratios of up to 5.7 for patients in risk groups. In addition to describing genetic findings, we also review the first results of epigenetic research showing constitutional methylation changes that alter the susceptibility to cutaneous melanoma and its risk factors.

Genomic Risk Score for Melanoma in a Prospective Study of Older Individuals

BACKGROUND

Recent genome-wide association meta-analysis for melanoma doubled the number of previously identified variants. We assessed the performance of an updated polygenic risk score (PRS) in a population of older individuals, where melanoma incidence and cumulative ultraviolet radiation exposure is greatest.

METHODS

We assessed a PRS for cutaneous melanoma comprising 55 variants in a prospective study of 12 712 individuals in the ASPirin in Reducing Events in the Elderly Trial. We evaluated incident melanomas diagnosed during the trial and prevalent melanomas diagnosed preenrolment (self-reported). Multivariable models examined associations between PRS as a continuous variable (per SD) and categorical (low-risk [0%-20%], medium-risk [21%-80%], high-risk [81%-100%] groups) with incident melanoma. Logistic regression examined the association between PRS and prevalent melanoma.

RESULTS

At baseline, mean participant age was 75 years; 55.0% were female, and 528 (4.2%) had prevalent melanomas. During follow-up (median = 4.7 years), 120 (1.0%) incident cutaneous melanomas occurred, 98 of which were in participants with no history. PRS was associated with incident melanoma (hazard ratio = 1.46 per SD, 95% confidence interval [CI] = 1.20 to 1.77) and prevalent melanoma (odds ratio [OR] = 1.55 per SD, 95% CI = 1.42 to 1.69). Participants in the highest-risk PRS group had increased risk compared with the low-risk group for incident melanoma (OR = 2.51, 95% CI = 1.28 to 4.92) and prevalent melanoma (OR = 3.66, 95% CI = 2.69 to 5.05). When stratifying by sex, only males had an association between the PRS and incident melanoma, whereas both sexes had an association between the PRS and prevalent melanoma.

CONCLUSIONS

A genomic risk score is associated with melanoma risk in older individuals and may contribute to targeted surveillance.

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

The association between citrus consumption and melanoma risk in the UK Biobank

BACKGROUND

Melanoma incidence has been dramatically increasing worldwide. Psoralen, a known photocarcinogen, is naturally abundant in citrus products, leading to the hypothesis that high citrus consumption may increase melanoma risk.

OBJECTIVES

To investigate the association between total citrus consumption and melanoma risk, and the association between individual citrus products and melanoma risk, and to test for interactions between total citrus intake and established melanoma risk factors.

METHODS

Logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the association between citrus consumption and melanoma risk among 1592 cases and 197 372 controls from the UK Biobank cohort. Citrus consumption data were collected via five rounds of 24-h recall questionnaires. International Classification of Diseases codes were used to determine melanoma outcome.

RESULTS

After adjusting for potential confounders, participants in the highest category of total citrus intake (> 2 servings per day) had a significantly increased risk of melanoma (OR 1·63, 95% CI 1·24-2·12) relative to those with no consumption. For individual citrus products, participants with the most orange and orange juice consumption (> 1 serving per day) had a significantly increased melanoma risk relative to those with no consumption (OR 1·79, 95% CI 1·07-2·78 and OR 1·54, 95% CI 1·10-2·10, respectively). Fair- or very fair-skinned participants with high citrus consumption had an even greater melanoma risk (OR 1·75, 95% CI 1·31-2·29).

CONCLUSIONS

High citrus consumption was associated with an increased risk of melanoma in a large, prospective, population-based cohort. Further validation of these findings could lead to improved melanoma prevention strategies.

A custom capture sequence approach for oculocutaneous albinism identifies structural variant alleles at the OCA2 locus

Oculocutaneous albinism (OCA) is a heritable disorder of pigment production that manifests as hypopigmentation and altered eye development. Exon sequencing of known OCA genes is unsuccessful in producing a complete molecular diagnosis for a significant number of affected individuals. We sequenced the DNA of individuals with OCA using short-read custom capture sequencing that targeted coding, intronic, and noncoding regulatory regions of known OCA genes, and genome-wide association study-associated pigmentation loci. We identified an OCA2 complex structural variant (CxSV), defined by a 143 kb inverted segment reintroduced in intron 1, upstream of the native location. The corresponding CxSV junctions were observed in 11/390 probands screened. The 143 kb CxSV presents in one family as a copy number variant duplication for the 143 kb region. In the remaining 10/11 families, the 143 kb CxSV acquired an additional 184 kb deletion across the same region, restoring exons 3-19 of OCA2 to a copy-number neutral state. Allele-associated haplotype analysis found rare SNVs rs374519281 and rs139696407 are linked with the 143 kb CxSV in both OCA2 alleles. For individuals in which customary molecular evaluation does not reveal a biallelic OCA diagnosis, we recommend preliminary screening for these haplotype-associated rare variants, followed by junction-specific validation for the OCA2 143 kb CxSV.

Advancing drug discovery using the power of the human genome

Human genetics plays an increasingly important role in drug development and population health. Here we review the history of human genetics in the context of accelerating the discovery of therapies, present examples of how human genetics evidence supports successful drug targets, and discuss how polygenic risk scores could be beneficial in various clinical settings. We highlight the value of direct-to-consumer platforms in the era of fast-paced big data biotechnology, and how diverse genetic and health data can benefit society. © 2021 23andMe, Inc. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

A genome-wide association study identifies novel gene associations with facial skin wrinkling and mole count in Latin Americans

BACKGROUND

Genome-wide association studies (GWASs) have identified genes influencing skin ageing and mole count in Europeans, but little is known about the relevance of these (or other genes) in non-Europeans.

OBJECTIVES

To conduct a GWAS for facial skin ageing and mole count in adults < 40 years old, of mixed European, Native American and African ancestry, recruited in Latin America.

METHODS

Skin ageing and mole count scores were obtained from facial photographs of over 6000 individuals. After quality control checks, three wrinkling traits and mole count were retained for genetic analyses. DNA samples were genotyped with Illumina's HumanOmniExpress chip. Association testing was performed on around 8 703 729 single-nucleotide polymorphisms (SNPs) across the autosomal genome.

RESULTS

Genome-wide significant association was observed at four genome regions: two were associated with wrinkling (in 1p13·3 and 21q21·2), one with mole count (in 1q32·3) and one with both wrinkling and mole count (in 5p13·2). Associated SNPs in 5p13·2 and in 1p13·3 are intronic within SLC45A2 and VAV3, respectively, while SNPs in 1q32·3 are near the SLC30A1 gene, and those in 21q21·2 occur in a gene desert. Analyses of SNPs in IRF4 and MC1R are consistent with a role of these genes in skin ageing.

CONCLUSIONS

We replicate the association of wrinkling with variants in SLC45A2, IRF4 and MC1R reported in Europeans. We identify VAV3 and SLC30A1 as two novel candidate genes impacting on wrinkling and mole count, respectively. We provide the first evidence that SLC45A2 influences mole count, in addition to variants in this gene affecting melanoma risk in Europeans.

Integrative Analysis of Omics Data Reveals Regulatory Network of CDK10 in Vitiligo Risk

Vitiligo is a multifactorial polygenic disorder, characterized by acquired depigmented skin and overlying hair resulting from the destruction of melanocytes. Genome-wide association studies (GWASs) of vitiligo have identified approximately 100 genetic variants. However, the identification of functional genes and their regulatory elements remains a challenge. To prioritize putative functional genes and DNAm sites, we performed a Summary data-based Mendelian Randomization (SMR) and heterogeneity in dependent instruments (HEIDI) test to integrate omics summary statistics from GWAS, expression quantitative trait locus (eQTL), and methylation quantitative trait loci (meQTL) analysis of large sample size. By integrating omics data, we identified two newly putative functional genes (SPATA2L and CDK10) associated with vitiligo and further validated CDK10 by qRT-PCR in independent samples. We also identified 17 vitiligo-associated DNA methylation (DNAm) sites in Chr16, of which cg05175606 was significantly associated with the expression of CDK10 and vitiligo. Colocalization analyses detected transcript of CDK10 in the blood and skin colocalizing with cg05175606 at single nucleotide polymorphism (SNP) rs77651727. Our findings revealed that a shared genetic variant rs77651727 alters the cg05175606 as well as up-regulates gene expression of CDK10 and further decreases the risk of vitiligo.

Cancer PRSweb: An Online Repository with Polygenic Risk Scores for Major Cancer Traits and Their Evaluation in Two Independent Biobanks

To facilitate scientific collaboration on polygenic risk scores (PRSs) research, we created an extensive PRS online repository for 35 common cancer traits integrating freely available genome-wide association studies (GWASs) summary statistics from three sources: published GWASs, the NHGRI-EBI GWAS Catalog, and UK Biobank-based GWASs. Our framework condenses these summary statistics into PRSs using various approaches such as linkage disequilibrium pruning/p value thresholding (fixed or data-adaptively optimized thresholds) and penalized, genome-wide effect size weighting. We evaluated the PRSs in two biobanks: the Michigan Genomics Initiative (MGI), a longitudinal biorepository effort at Michigan Medicine, and the population-based UK Biobank (UKB). For each PRS construct, we provide measures on predictive performance and discrimination. Besides PRS evaluation, the Cancer-PRSweb platform features construct downloads and phenome-wide PRS association study results (PRS-PheWAS) for predictive PRSs. We expect this integrated platform to accelerate PRS-related cancer research.

Overexpression of BASP1 Indicates a Poor Prognosis in Head and Neck Squamous Cell Carcinoma

Objective:

Brain abundant membrane attached signal protein 1 (BASP1) was originally identified as a membrane and cytoplasmic protein. Recent studies have shown that BASP1 highly expressed in cancer and promoted the proliferation of cancer. However, the role of BASP1 in head and neck squamous cell carcinoma (HNSCC) is largely unknown. Here, we performed a systematic data analysis to examine whether BASP1 can function as prognostic marker in HNSCC.

Methods:

In this study, we used Oncomine, and UALCAN, databases to analyze the expression of BASP1 in HNSCC. We used Kaplan-Meier plotter to evaluate the effect of BASP1 on clinical prognosis. In addition, we also analyzed genetic alterations, interaction network, and functional enrichment of BASP1.

Results:

BASP1 mRNA expression level was remarkably increased in HNSCC than in normal tissues (P=1.624e-12). Moreover, high BASP1 expression was significantly related to poor survival (p=0.00056) in HNSCC patients. In addition, BASP1 gene amplified in 5% of HNSCC patients which contributes to the overexpression of BASP1.

Conclusions:

These findings suggest that BASP1 was frequently amplified which contributes to the overexpression of BASP1, thereby promoting HNSCC progression. Thus, these results indicate that BASP1 might serve as a biomarker to predict the progression and prognosis of HNSCC patients.

Multiplex melanoma families are enriched for polygenic risk

Cancers, including cutaneous melanoma, can cluster in families. In addition to environmental etiological factors such as ultraviolet radiation, cutaneous melanoma has a strong genetic component. Genetic risks for cutaneous melanoma range from rare, high-penetrance mutations to common, low-penetrance variants. Known high-penetrance mutations account for only about half of all densely affected cutaneous melanoma families, and the causes of familial clustering in the remainder are unknown. We hypothesize that some clustering is due to the cumulative effect of a large number of variants of individually small effect. Common, low-penetrance genetic risk variants can be combined into polygenic risk scores. We used a polygenic risk score for cutaneous melanoma to compare families without known high-penetrance mutations with unrelated melanoma cases and melanoma-free controls. Family members had significantly higher mean polygenic load for cutaneous melanoma than unrelated cases or melanoma-free healthy controls (Bonferroni-corrected t-test P = 1.5 × 10-5 and 6.3 × 10-45, respectively). Whole genome sequencing of germline DNA from 51 members of 21 families with low polygenic risk for melanoma identified a CDKN2A p.G101W mutation in a single family but no other candidate high-penetrance melanoma susceptibility genes. This work provides further evidence that melanoma, like many other common complex disorders, can arise from the joint action of multiple predisposing factors, including rare high-penetrance mutations, as well as via a combination of large numbers of alleles of small effect.

Inherited variations in human pigmentation-related genes modulate cutaneous melanoma risk and clinicopathological features in Brazilian population

Ultraviolet light exposure and cutaneous pigmentation are important host risk factors for cutaneous melanoma (CM), and it is well known that inherited ability to produce melanin varies in humans. The study aimed to identify single-nucleotide variants (SNVs) on pigmentation-related genes with importance in risk and clinicopathological aspects of CM. The study was conducted in two stages. In stage 1, 103 CM patients and 103 controls were analyzed using Genome-Wide Human SNV Arrays in order to identify SNVs in pigmentation-related genes, and the most important SNVs were selected for data validation in stage 2 by real-time polymerase-chain reaction in 247 CM patients and 280 controls. ADCY3 c.675+9196T>G, CREB1 c.303+373G>A, and MITF c.938-325G>A were selected for data validation among 74 SNVs. Individuals with CREB1 GA or AA genotype and allele "A" were under 1.79 and 1.47-fold increased risks of CM than others, respectively. Excesses of CREB1 AA and MITF AA genotype were seen in patients with tumors at Clark levels III to V (27.8% versus 13.7%) and at III or IV stages (46.1% versus 24.9%) compared to others, respectively. When compared to others, patients with ADCY3 TT had 1.89 more chances of presenting CM progression, and those with MITF GA or AA had 2.20 more chances of evolving to death by CM. Our data provide, for the first time, preliminary evidence that inherited abnormalities in ADCY3, CREB1, and MITF pigmentation-related genes, not only can increase the risk to CM, but also influence CM patients' clinicopathological features.

Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma

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^V600E background. Our integrative approach streamlines GWAS follow-up studies and highlights a pleiotropic function of MX2 in melanoma susceptibility.

There are more than 20 known melanoma susceptibility genes. Here, using a massively parallel reporter assay, the authors identify risk-associated variants that alter gene transcription, and demonstrate that expression of one such gene, MX2, leads to the promotion of melanoma in a zebrafish model.

Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility

Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 × 10-8) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.

Evolution of approaches to identify melanoma missing heritability

Introduction: Around 10% of melanoma patients have a positive family history of melanoma and/or related cancers. Although a germline pathogenic variant in a high-risk gene can be identified in up to 40% of these patients, the remaining part of melanoma heritability remains largely unexplained.Areas covered: The aim of this review is to provide an overview of the impact that new technologies and new research approaches had and are having on finding more efficient ways to unravel the missing heritability in melanoma.Expert opinion: High-throughput sequencing technologies have been crucial in increasing the number of genes/loci that might be implicated in melanoma predisposition. However, results from these approaches may have been inferior to the expectations, due to an increase in quantitative information which hasn't been followed at the same speed by an improvement of the methods to correctly interpret these data. Optimal approaches for improving our knowledge on melanoma heritability are currently based on segregation analysis coupled with functional assessment of candidate genes. An improvement of computational methods to infer genotype-phenotype correlations could help address the issue of missing heritability.

Exome Sequencing Analysis Identifies Rare Variants in ATM and RPL8 That Are Associated With Shorter Telomere Length

Telomeres are important for maintaining genomic stability. Telomere length has been associated with aging, disease, and mortality and is highly heritable (∼82%). In this study, we aimed to identify rare genetic variants associated with telomere length using whole-exome sequence data. We studied 1,303 participants of the Erasmus Rucphen Family (ERF) study, 1,259 of the Rotterdam Study (RS), and 674 of the British Heart Foundation Family Heart Study (BHF-FHS). We conducted two analyses, first we analyzed the family-based ERF study and used the RS and BHF-FHS for replication. Second, we combined the summary data of the three studies in a meta-analysis. Telomere length was measured by quantitative polymerase chain reaction in blood. We identified nine rare variants significantly associated with telomere length (p-value < 1.42 × 10^–7, minor allele frequency of 0.2–0.5%) in the ERF study. Eight of these variants (in C11orf65, ACAT1, NPAT, ATM, KDELC2, and EXPH5) were located on chromosome 11q22.3 that contains ATM, a gene involved in telomere maintenance. Although we were unable to replicate the variants in the RS and BHF-FHS (p-value ≥ 0.21), segregation analysis showed that all variants segregate with shorter telomere length in a family. In the meta-analysis of all studies, a nominally significant association with LTL was observed with a rare variant in RPL8 (p-value = 1.48 × 10⁻⁶), which has previously been associated with age. Additionally, a novel rare variant in the known RTEL1 locus showed suggestive evidence for association (p-value = 1.18 × 10^–4) with LTL. To conclude, we identified novel rare variants associated with telomere length. Larger samples size are needed to confirm these findings and to identify additional variants.

Identification of a Gene-Related Risk Signature in Melanoma Patients Using Bioinformatic Profiling

Introduction

Gene signature has been used to predict prognosis in melanoma patients. Meanwhile, the efficacy of immunotherapy was correlated with particular genes expression or mutation. In this study, we systematically explored the gene expression pattern in the melanoma-immune microenvironment and its relationship with prognosis.

Methods

A cohort of 122 melanoma cases with whole-genome microarray expression data were enrolled from the Gene Expression Omnibus (GEO) database. The findings were validated using The Cancer Genome Atlas (TCGA) database. A principal component analysis (PCA), gene set enrichment analysis (GSEA), and gene oncology (GO) analysis were performed to explore the bioinformatic implications.

Results

Different gene expression patterns were identified according to the clinical stage. All eligible gene sets were analyzed, and the 8 genes (GPR87, KIT, SH3GL3, PVRL1, ATP1B1, CDAN1, FAU, and TNFSF14) with the greatest prognostic impact on melanoma. A gene-related risk signature was developed to distinguish patients with a high or low risk of an unfavorable outcome, and this signature was validated using the TCGA database. Furthermore, the prognostic significance of the signature between the classified subgroups was verified as an independent prognostic predictor of melanoma. Additionally, the low-risk melanoma patients presented an enhanced immune phenotype compared to that of the high-risk gene signature patients.

Conclusions

The gene pattern differences in melanoma were profiled, and a gene signature that could independently predict melanoma patients with a high risk of poor survival was established, highlighting the relationship between prognosis and the local immune response.

The brain acid-soluble protein 1 (BASP1) interferes with the oncogenic capacity of MYC and its binding to calmodulin

The MYC protein is a transcription factor with oncogenic potential controlling fundamental cellular processes such as cell proliferation, metabolism, differentiation, and apoptosis. The MYC gene is a major cancer driver, and elevated MYC protein levels are a hallmark of most human cancers. We have previously shown that the brain acid-soluble protein 1 gene (BASP1) is specifically downregulated by the v-myc oncogene and that ectopic BASP1 expression inhibits v-myc-induced cell transformation. The 11-amino acid effector domain of the BASP1 protein interacts with the calcium sensor calmodulin (CaM) and is mainly responsible for this inhibitory function. We also reported recently that CaM interacts with all MYC variant proteins and that ectopic CaM increases the transactivation and transformation potential of the v-Myc protein. Here, we show that the presence of excess BASP1 or of a synthetic BASP1 effector domain peptide leads to displacement of v-Myc from CaM. The protein stability of v-Myc is decreased in cells co-expressing v-Myc and BASP1, which may account for the inhibition of v-Myc. Furthermore, suppression of v-Myc-triggered transcriptional activation and cell transformation is compensated by ectopic CaM, suggesting that BASP1-mediated withdrawal of CaM from v-Myc is a crucial event in the inhibition. In view of the tumor-suppressive role of BASP1 which was recently also reported for human cancer, small compounds or peptides based on the BASP1 effector domain could be used in drug development strategies aimed at tumors with high MYC expression.

Combined analysis of keratinocyte cancers identifies novel genome-wide loci

The keratinocyte cancers (KC), basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are the most common cancers in fair-skinned people. KC treatment represents the second highest cancer healthcare expenditure in Australia. Increasing our understanding of the genetic architecture of KC may provide new avenues for prevention and treatment. We first conducted a series of genome-wide association studies (GWAS) of KC across three European ancestry datasets from Australia, Europe and USA, and used linkage disequilibrium (LD) Score regression (LDSC) to estimate their pairwise genetic correlations. We employed a multiple-trait approach to map genes across the combined set of KC GWAS (total N = 47 742 cases, 634 413 controls). We also performed meta-analyses of BCC and SCC separately to identify trait specific loci. We found substantial genetic correlations (generally 0.5-1) between BCC and SCC suggesting overlapping genetic risk variants. The multiple trait combined KC GWAS identified 63 independent genome-wide significant loci, 29 of which were novel. Individual separate meta-analyses of BCC and SCC identified an additional 13 novel loci not found in the combined KC analysis. Three new loci were implicated using gene-based tests. New loci included common variants in BRCA2 (distinct to known rare high penetrance cancer risk variants), and in CTLA4, a target of immunotherapy in melanoma. We found shared and trait specific genetic contributions to BCC and SCC. Considering both, we identified a total of 79 independent risk loci, 45 of which are novel.

Genome-wide association studies and polygenic risk scores for skin cancer: clinically useful yet?

BACKGROUND

Genome-wide association studies (GWAS) have identified thousands of susceptibility variants, although most have been associated with small individual risk estimates that offer little predictive value. However, combining multiple variants into polygenic risk scores (PRS) may be more informative. Multiple studies have developed PRS composed of GWAS-identified variants for cutaneous cancers. This review highlights data from these studies.

OBJECTIVES

To review published GWAS and PRS studies for melanoma, cutaneous squamous cell carcinoma (cSCC) and basal cell carcinoma (BCC), and discuss their potential clinical utility.

METHODS

We searched PubMed and the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue to identify relevant studies.

RESULTS

Results from 21 GWAS (11 melanoma, 3 cSCC, 7 BCC) and 11 PRS studies are summarized. Six loci in pigmentation genes overlap between these three cancers (ASIP/RALY, IRF4, MC1R, OCA2, SLC45A2 and TYR). Additional loci overlap for cSCC/BCC and BCC/melanoma, but no other loci are shared between cSCC and melanoma. PRS for melanoma show roughly two-to-threefold increases in risk and modest improvements in risk prediction (2-7% increases). PRS are associated with twofold and threefold increases in risk of cSCC and BCC, respectively, with small improvements (2% increase) in predictive ability.

CONCLUSIONS

Existing data indicate that PRS may offer small, but potentially meaningful, improvements to risk prediction. Additional research is needed to clarify the potential utility of PRS in cutaneous carcinomas. Clinical translation will require well-powered validation studies incorporating known risk factors to evaluate PRS as tools for screening. What's already known about this topic? Over 50 susceptibility loci for melanoma, basal cell carcinoma (BCC) and cutaneous squamous cell carcinoma (cSCC) have been identified in genome-wide association studies (GWAS). Polygenic risk scores (PRS) using variants identified from GWAS have also been developed for melanoma, BCC and cSCC, and investigated with respect to clinical risk prediction. What does this study add? This review provides an overview of GWAS findings and the potential clinical utility of PRS for melanoma, BCC and cSCC.

A Trans-Ethnic Genome-Wide Association Study of Uterine Fibroids

Uterine fibroids affect up to 77% of women by menopause and account for up to $34 billion in healthcare costs each year. Although fibroid risk is heritable, genetic risk for fibroids is not well understood. We conducted a two-stage case-control meta-analysis of genetic variants in European and African ancestry women with and without fibroids classified by a previously published algorithm requiring pelvic imaging or confirmed diagnosis. Women from seven electronic Medical Records and Genomics (eMERGE) network sites (3,704 imaging-confirmed cases and 5,591 imaging-confirmed controls) and women of African and European ancestry from UK Biobank (UKB, 5,772 cases and 61,457 controls) were included in the discovery genome-wide association study (GWAS) meta-analysis. Variants showing evidence of association in Stage I GWAS (P < 1 × 10-5) were targeted in an independent replication sample of African and European ancestry individuals from the UKB (Stage II) (12,358 cases and 138,477 controls). Logistic regression models were fit with genetic markers imputed to a 1000 Genomes reference and adjusted for principal components for each race- and site-specific dataset, followed by fixed-effects meta-analysis. Final analysis with 21,804 cases and 205,525 controls identified 326 genome-wide significant variants in 11 loci, with three novel loci at chromosome 1q24 (sentinel-SNP rs14361789; P = 4.7 × 10-8), chromosome 16q12.1 (sentinel-SNP rs4785384; P = 1.5 × 10-9) and chromosome 20q13.1 (sentinel-SNP rs6094982; P = 2.6 × 10-8). Our statistically significant findings further support previously reported loci including SNPs near WT1, TNRC6B, SYNE1, BET1L, and CDC42/WNT4. We report evidence of ancestry-specific findings for sentinel-SNP rs10917151 in the CDC42/WNT4 locus (P = 1.76 × 10-24). Ancestry-specific effect-estimates for rs10917151 were in opposite directions (P-Het-between-groups = 0.04) for predominantly African (OR = 0.84) and predominantly European women (OR = 1.16). Genetically-predicted gene expression of several genes including LUZP1 in vagina (P = 4.6 × 10-8), OBFC1 in esophageal mucosa (P = 8.7 × 10-8), NUDT13 in multiple tissues including subcutaneous adipose tissue (P = 3.3 × 10-6), and HEATR3 in skeletal muscle tissue (P = 5.8 × 10-6) were associated with fibroids. The finding for HEATR3 was supported by SNP-based summary Mendelian randomization analysis. Our study suggests that fibroid risk variants act through regulatory mechanisms affecting gene expression and are comprised of alleles that are both ancestry-specific and shared across continental ancestries.

Exploring various polygenic risk scores for skin cancer in the phenomes of the Michigan genomics initiative and the UK Biobank with a visual catalog: PRSWeb

Polygenic risk scores (PRS) are designed to serve as single summary measures that are easy to construct, condensing information from a large number of genetic variants associated with a disease. They have been used for stratification and prediction of disease risk. The primary focus of this paper is to demonstrate how we can combine PRS and electronic health records data to better understand the shared and unique genetic architecture and etiology of disease subtypes that may be both related and heterogeneous. PRS construction strategies often depend on the purpose of the study, the available data/summary estimates, and the underlying genetic architecture of a disease. We consider several choices for constructing a PRS using data obtained from various publicly-available sources including the UK Biobank and evaluate their abilities to predict not just the primary phenotype but also secondary phenotypes derived from electronic health records (EHR). This study was conducted using data from 30,702 unrelated, genotyped patients of recent European descent from the Michigan Genomics Initiative (MGI), a longitudinal biorepository effort within Michigan Medicine. We examine the three most common skin cancer subtypes in the USA: basal cell carcinoma, cutaneous squamous cell carcinoma, and melanoma. Using these PRS for various skin cancer subtypes, we conduct a phenome-wide association study (PheWAS) within the MGI data to evaluate PRS associations with secondary traits. PheWAS results are then replicated using population-based UK Biobank data and compared across various PRS construction methods. We develop an accompanying visual catalog called PRSweb that provides detailed PheWAS results and allows users to directly compare different PRS construction methods.

A Unique Family of Neuronal Signaling Proteins Implicated in Oncogenesis and Tumor Suppression

The neuronal proteins GAP43 (neuromodulin), MARCKS, and BASP1 are highly expressed in the growth cones of nerve cells where they are involved in signal transmission and cytoskeleton organization. Although their primary structures are unrelated, these signaling proteins share several structural properties like fatty acid modification, and the presence of cationic effector domains. GAP43, MARCKS, and BASP1 bind to cell membrane phospholipids, a process reversibly regulated by protein kinase C-phosphorylation or by binding to the calcium sensor calmodulin (CaM). GAP43, MARCKS, and BASP1 are also expressed in non-neuronal cells, where they may have important functions to manage cytoskeleton architecture, and in case of MARCKS and BASP1 to act as cofactors in transcriptional regulation. During neoplastic cell transformation, the proteins reveal differential expression in normal vs. tumor cells, and display intrinsic tumor promoting or tumor suppressive activities. Whereas GAP43 and MARCKS are oncogenic, tumor suppressive functions have been ascribed to BASP1 and in part to MARCKS depending on the cell type. Like MARCKS, the myristoylated BASP1 protein is localized both in the cytoplasm and in the cell nucleus. Nuclear BASP1 participates in gene regulation converting the Wilms tumor transcription factor WT1 from an oncoprotein into a tumor suppressor. The BASP1 gene is downregulated in many human tumor cell lines particularly in those derived from leukemias, which display elevated levels of WT1 and of the major cancer driver MYC. BASP1 specifically inhibits MYC-induced cell transformation in cultured cells. The tumor suppressive functions of BASP1 and MARCKS could be exploited to expand the spectrum of future innovative therapeutic approaches to inhibit growth and viability of susceptible human tumors.

Systematic evaluation of cancer-specific genetic risk score for 11 types of cancer in The Cancer Genome Atlas and Electronic Medical Records and Genomics cohorts

BACKGROUND

Genetic risk score (GRS) is an odds ratio (OR)-weighted and population-standardized method for measuring cumulative effect of multiple risk-associated single nucleotide polymorphisms (SNPs). We hypothesize that GRS is a valid tool for risk assessment of most common cancers.

METHODS

Utilizing genotype and phenotype data from The Cancer Genome Atlas (TCGA) and Electronic Medical Records and Genomics (eMERGE), we tested 11 cancer-specific GRSs (bladder, breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, prostate, renal, and thyroid cancer) for association with the respective cancer type. Cancer-specific GRSs were calculated, for the first time in these cohorts, based on previously published risk-associated SNPs using the Caucasian subjects in these two cohorts.

RESULTS

Mean cancer-specific GRS in the population controls of eMERGE approximated the expected value of 1.00 (between 0.98 and 1.02) for all 11 types of cancer. Mean cancer-specific GRS was consistently higher in respective cancer patients than controls for all 11 types of cancer (P < 0.05). When subjects were categorized into low-, average-, and high-risk groups based on cancer-specific GRS (<0.5, 0.5-1.5, and >1.5, respectively), significant dose-response associations of higher cancer-specific GRS with higher OR of respective type of cancer were found for nine types of cancer (P_-trend  < 0.05). More than 64% subjects in the population controls of eMERGE can be classified as high risk for at least one type of these cancers.

CONCLUSION

Validity of GRS for predicting cancer risk is demonstrated for most types of cancer. If confirmed in larger studies, cancer-specific GRS may have the potential for developing personalized cancer screening strategy.

Familial Risks and Mortality in Second Primary Cancers in Melanoma

Background

Malignant melanoma (MM) patients are at increasing risk of developing second primary cancers (SPCs). We assessed mortality and risk of SPCs in MM patients with siblings or parents affected with same cancer compared with that of the general population.

Methods

We used the Swedish Family-Cancer Database to assess relative risks (RRs) and causes of death in SPCs until 2015 in patients with a MM diagnosis between 1958 and 2015. We identified 35 451patients with MM among whom 3212 received a subsequent diagnosis of SPC. RRs of SPCs after MM diagnosis were calculated stratifying over concordant family history of cancer in first-degree relatives.

Results

Familial RRs were increased for second melanoma (RR = 19.28, 95% CI = 16.71 to 22.25), squamous cell skin cancer (RR = 7.58, 95% CI = 5.57 to 10.29), leukemia (RR = 5.69, 95% CI = 2.96 to 10.94), bladder (RR = 4.15, 95% CI = 2.50 to 6.89), ovarian (RR = 3.89, 95% CI = 1.46 to 10.37), kidney cancer (RR = 3.77, 95% CI = 1.57 to 9.06), cancer of unknown primary (RR = 3.67, 95% CI = 1.65 to 8.16), nervous system (RR = 2.88, 95% CI = 1.20 to 6.93), breast (RR = 2.34, 95% CI = 1.92 to 2.84), lung (RR = 2.24, 95% CI = 1.50 to 3.35), and prostate cancer (RR = 2.22, 95% CI = 1.89 to 2.61) with statistical significance. For all cancers, familial RR was in excess (2.09, 95% CI = 2.02 to 2.16 vs 1.78, 95% CI = 1.69 to 1.87; P_trend < .0001). Cause of death in MM patients with SPC is shown to be dependent on the cancer site though SPCs contributed to majority of deaths.

Conclusions

SPCs appear higher with prior family history of cancer and contribute to mortality. SPC was the most common cause of death in patients with SPC and is almost uniformly the major contributing cause of death for all cancer sites. For improved survival in MM patients, prevention and early detection of SPCs would be important.

Sex-specific gene and pathway modeling of inherited glioma risk

Background

To date, genome-wide association studies (GWAS) have identified 25 risk variants for glioma, explaining 30% of heritable risk. Most histologies occur with significantly higher incidence in males, and this difference is not explained by currently known risk factors. A previous GWAS identified sex-specific glioma risk variants, and this analysis aims to further elucidate risk variation by sex using gene- and pathway-based approaches.

Methods

Results from the Glioma International Case-Control Study were used as a testing set, and results from 3 GWAS were combined via meta-analysis and used as a validation set. Using summary statistics for nominally significant autosomal SNPs (P < 0.01 in a previous meta-analysis) and nominally significant X-chromosome SNPs (P < 0.01), 3 algorithms (Pascal, BimBam, and GATES) were used to generate gene scores, and Pascal was used to generate pathway scores. Results were considered statistically significant in the discovery set when P < 3.3 × 10-6 and in the validation set when P < 0.001 in 2 of 3 algorithms.

Results

Twenty-five genes within 5 regions and 19 genes within 6 regions reached statistical significance in at least 2 of 3 algorithms in males and females, respectively. EGFR was significantly associated with all glioma and glioblastoma in males only and a female-specific association in TERT, all of which remained nominally significant after conditioning on known risk loci. There were nominal associations with the BioCarta telomeres pathway in both males and females.

Conclusions

These results provide additional evidence that there may be differences by sex in genetic risk for glioma. Additional analyses may further elucidate the biological processes through which this risk is conferred.

An atlas of genetic associations in UK Biobank

Genome-wide association studies (GWAS) have identified many loci contributing to variation in complex traits, yet the majority of loci that contribute to the heritability of complex traits remain elusive. Large study populations with sufficient statistical power are required to detect the small effect sizes of the yet unidentified genetic variants. However, the analysis of huge cohorts, like UK Biobank, is challenging. Here, we present an atlas of genetic associations for 118 non-binary and 660 binary traits of 452,264 UK Biobank participants of European ancestry. Results are compiled in a publicly accessible database that allows querying genome-wide association results for 9,113,133 genetic variants, as well as downloading GWAS summary statistics for over 30 million imputed genetic variants (>23 billion phenotype-genotype pairs). Our atlas of associations (GeneATLAS, http://geneatlas.roslin.ed.ac.uk ) will help researchers to query UK Biobank results in an easy and uniform way without the need to incur high computational costs.

Off-target phenotypes in forensic DNA phenotyping and biogeographic ancestry inference: A resource

With recent advances in DNA sequencing technologies it has become feasible and cost effective to genotype larger marker sets for forensic purposes. Two technologies that make use of the larger marker sets have come into focus in forensic research and applications; inference of biogeographic ancestry (BGA) and forensic DNA phenotyping (FDP). These methods hold the promise to reveal information about a yet unknown perpetrator from a DNA sample. In contrast, DNA-profiling, that is a standard practice in case work, relies on matching DNA-profiles between crime scene material and suspects on a database of DNA-profiles. Markers for DNA-profiling were developed under the premise to reveal as little additional information about the human source of the profile as possible, the rationale being that personal privacy rights have to be balanced against the public interest in solving a crime. The same argument holds for markers used in BGA and FDP; these markers might also reveal information on off-target phenotypes (OTPs), that go beyond BGA and the phenotypes targeted in FDP. In particular, health related OTPs might shift the balance between privacy protection and public interest. However, to our knowledge, there is currently no convenient resource available to incorporate knowledge on OTPs in BGA and FDP assay design and application. In order to provide such a resource, we performed a systematic search for OTPs associated with a comprehensive set of markers (1766 SNPs) used or suggested to be used for BGA inference and FDP. In this set, we identified a relatively small number of 27 SNPs (1.53%) that convey information on diverse health related OTPs such as cancer risk, induced asthma, or risk of alcoholism. Some of these SNPs are commonly used for FDP and BGA across different marker sets. We conclude that the effects of SNP markers used in FDP and BGA on OTPs are currently limited, with few exceptions that should be considered in a balanced decision on assay design and application.

Genetic polymorphisms may influence the vertical growth rate of melanoma

Background: Identification of new predictive markers in melanoma is of great clinical importance. This study was aimed to analyze association between selected common variants in the cancer susceptibility genes and melanoma progression at the time of diagnosis. Material and Method: The study included 243 consecutive patients with melanoma. Genotyping was performed using real-time PCR. Results: Our data revealed modest association between xeroderma pigmentosum complementation group D (XPD) codon 312 polymorphism and tumor thickness (as defined by Breslow score; XPD D312N CC: 3.00 ± 3.78mm, CT: 1.71 ± 2.48mm, TT: 2,53 ± 3,24mm, P=0.023). The CT genotype in XPD D312N polymorphism was more frequently represented in non-invasive melanomas compared to deeply penetrating tumors. None of the common SNPs in cyclin dependent kinase inhibitor 2A (CDKN2A), vitamin D receptor (VDR), melanocortin 1 receptor (MC1R) were associated with Breslow depth. Conclusion: These findings suggest that genetic alteration in XPD contributes to melanoma progression and may be a potential diagnostic and molecular prognostic marker.

IDPpi: Protein-Protein Interaction Analyses of Human Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are characterized by the lack of a fixed tertiary structure and are involved in the regulation of key biological processes via binding to multiple protein partners. IDPs are malleable, adapting to structurally different partners, and this flexibility stems from features encoded in the primary structure. The assumption that universal sequence information will facilitate coverage of the sparse zones of the human interactome motivated us to explore the possibility of predicting protein-protein interactions (PPIs) that involve IDPs based on sequence characteristics. We developed a method that relies on features of the interacting and non-interacting protein pairs and utilizes machine learning to classify and predict IDP PPIs. Consideration of both sequence determinants specific for conformational organizations and the multiplicity of IDP interactions in the training phase ensured a reliable approach that is superior to current state-of-the-art methods. By applying a strict evaluation procedure, we confirm that our method predicts interactions of the IDP of interest even on the proteome-scale. This service is provided as a web tool to expedite the discovery of new interactions and IDP functions with enhanced efficiency.

Genome-wide association study in 176,678 Europeans reveals genetic loci for tanning response to sun exposure

The skin’s tendency to sunburn rather than tan is a major risk factor for skin cancer. Here we report a large genome-wide association study of ease of skin tanning in 176,678 subjects of European ancestry. We identify significant association with tanning ability at 20 loci. We confirm previously identified associations at six of these loci, and report 14 novel loci, of which ten have never been associated with pigmentation-related phenotypes. Our results also suggest that variants at the AHR/AGR3 locus, previously associated with cutaneous malignant melanoma the underlying mechanism of which is poorly understood, might act on disease risk through modulation of tanning ability.

The skin’s tanning response to sun exposure shows great interindividual variability. Here, Visconti et al. perform a genome-wide association study for ease of skin tanning and identify 20 genetic loci, ten of which had not previously been associated with pigmentation-related traits.

Methylation-associated silencing of BASP1 contributes to leukemogenesis in t(8;21) acute myeloid leukemia

The AML1-ETO fusion protein (A/E), which results from the t(8;21) translocation, is considered to be a leukemia-initiating event. Identifying the mechanisms underlying the oncogenic activity of A/E remains a major challenge. In this study, we identified a specific down-regulation of brain acid-soluble protein 1 (BASP1) in t(8;21) acute myeloid leukemia (AML). A/E recognized AML1-binding sites and recruited DNA methyltransferase 3a (DNMT3a) to the BASP1 promoter sequence, which triggered DNA methylation-mediated silencing of BASP1. Ectopic expression of BASP1 inhibited proliferation and the colony-forming ability of A/E-positive AML cell lines and led to apoptosis and cell cycle arrest. The DNMT inhibitor decitabine up-regulated the expression of BASP1 in A/E-positive AML cell lines. In conclusion, our data suggest that BASP1 silencing via promoter methylation may be involved in A/E-mediated leukemogenesis and that BASP1 targeting may be an actionable therapeutic strategy in t(8;21) AML.

A chromosomal rearrangement commonly observed in certain leukemias selectively inactivates a gene that otherwise thwarts cancerous growth. Between 7 and 12% of acute myeloid leukemia cases exhibit a dramatic alteration in chromosomal structure that results in the production of an abnormal fusion protein. Researchers led by Li Yu at the General Hospital of Shenzen University in China have learned that this protein promotes disease progression by switching off an important tumor suppressor. Yu and colleagues showed that it binds a genomic sequence that regulates the gene encoding a second protein called BASP1, dramatically reducing its production. This gene silencing facilitates tumor growth. Chemicals that reactivated BASP1 production slowed proliferation and initiated ‘self-destruct’ mechanisms in leukemia cells. These findings suggest that BASP1-oriented therapies could offer a fruitful avenue of treatment for some patients.