Multi-ancestry genome-wide meta-analysis identifies novel basal cell carcinoma loci and shared genetic effects with squamous cell carcinoma.

Basal cell carcinoma (BCC) is one of the most common malignancies worldwide, yet its genetic determinants are incompletely defined. We perform a European ancestry genome-wide association (GWA) meta-analysis and a Hispanic/Latino ancestry GWA meta-analysis and meta-analyze both in a multi-ancestry GWAS meta-analysis of BCC, totaling 50,531 BCC cases and 762,234 controls from four cohorts (GERA, Mass-General Brigham Biobank, UK Biobank, and 23andMe research cohort). Here we identify 122 BCC-associated loci, of which 36 were novel, and subsequently fine-mapped these associations. We also identify an association of the well-known pigment gene SLC45A2 as well as associations at RCC2 and CLPTM1L with BCC in Hispanic/Latinos. We examine these BCC loci for association with cutaneous squamous cell carcinoma (cSCC) in 16,407 SCC cases and 762,486 controls of European ancestry, and 33 SNPs show evidence of association. Our study findings provide important insights into the genetic basis of BCC and cSCC susceptibility.

Association of Behavioral and Clinical Risk Factors With Cataract: A Two-Sample Mendelian Randomization Study.

Purpose: To investigate the association of genetically determined primary open-angle glaucoma (POAG), myopic refractive error (RE), type 2 diabetes (T2D), blood pressure (BP), body mass index (BMI), cigarette smoking, and alcohol consumption with the risk of age-related cataract. Methods: To assess potential causal effects of clinical or behavioral factors on cataract risk, we conducted two-sample Mendelian randomization analyses. Genetic instruments, based on common genetic variants associated with risk factors at genome-wide significance (P < 5 × 10-8), were derived from published genome-wide association studies (GWAS). For age-related cataract, we used GWAS summary statistics from our previous GWAS conducted in the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (28,092 cataract cases and 50,487 controls; all non-Hispanic whites) or in the UK Biobank (31,852 cataract cases and 428,084 controls; all European-descent individuals). We used the inverse-variance weighted (IVW) method as our primary source of Mendelian randomization estimates and conducted common sensitivity analyses. Results: We found that genetically determined POAG and mean spherical equivalent RE were significantly associated with cataract risk (IVW model: odds ratio [OR] = 1.04; 95% confidence interval [CI], 1.01-1.08; P = 0.018; per diopter more hyperopic: OR = 0.92; 95% CI, 0.89-0.93; P = 6.51 × 10-13, respectively). In contrast, genetically determined T2D, BP, BMI, cigarette smoking, or alcohol consumption were not associated with cataract risk (P > 0.05). Conclusions: Our results provide evidence that genetic risks for POAG and myopia may be causal risk factors for age-related cataract. These results are consistent with previous observational studies reporting associations of myopia with cataract risk. This information may support population cataract risk stratification and screening strategies.

Rare coding variants in CHRNB2 reduce the likelihood of smoking.

Human genetic studies of smoking behavior have been thus far largely limited to common variants. Studying rare coding variants has the potential to identify drug targets. We performed an exome-wide association study of smoking phenotypes in up to 749,459 individuals and discovered a protective association in CHRNB2, encoding the ß2 subunit of the α4ß2 nicotine acetylcholine receptor. Rare predicted loss-of-function and likely deleterious missense variants in CHRNB2 in aggregate were associated with a 35% decreased odds for smoking heavily (odds ratio (OR) = 0.65, confidence interval (CI) = 0.56-0.76, P = 1.9 × 10-8). An independent common variant association in the protective direction ( rs2072659 ; OR = 0.96; CI = 0.94-0.98; P = 5.3 × 10-6) was also evident, suggesting an allelic series. Our findings in humans align with decades-old experimental observations in mice that ß2 loss abolishes nicotine-mediated neuronal responses and attenuates nicotine self-administration. Our genetic discovery will inspire future drug designs targeting CHRNB2 in the brain for the treatment of nicotine addiction.

MiXcan: a framework for cell-type-aware transcriptome-wide association studies with an application to breast cancer.

Human bulk tissue samples comprise multiple cell types with diverse roles in disease etiology. Conventional transcriptome-wide association study approaches predict genetically regulated gene expression at the tissue level, without considering cell-type heterogeneity, and test associations of predicted tissue-level expression with disease. Here we develop MiXcan, a cell-type-aware transcriptome-wide association study approach that predicts cell-type-level expression, identifies disease-associated genes via combination of cell-type-level association signals for multiple cell types, and provides insight into the disease-critical cell type. As a proof of concept, we conducted cell-type-aware analyses of breast cancer in 58,648 women and identified 12 transcriptome-wide significant genes using MiXcan compared with only eight genes using conventional approaches. Importantly, MiXcan identified genes with distinct associations in mammary epithelial versus stromal cells, including three new breast cancer susceptibility genes. These findings demonstrate that cell-type-aware transcriptome-wide analyses can reveal new insights into the genetic and cellular etiology of breast cancer and other diseases.

Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Assessment of genetic susceptibility to multiple primary cancers through whole-exome sequencing in two large multi-ancestry studies.

BACKGROUND: Up to one of every six individuals diagnosed with one cancer will be diagnosed with a second primary cancer in their lifetime. Genetic factors contributing to the development of multiple primary cancers, beyond known cancer syndromes, have been underexplored. METHODS: To characterize genetic susceptibility to multiple cancers, we conducted a pan-cancer, whole-exome sequencing study of individuals drawn from two large multi-ancestry populations (6429 cases, 165,853 controls). We created two groupings of individuals diagnosed with multiple primary cancers: (1) an overall combined set with at least two cancers across any of 36 organ sites and (2) cancer-specific sets defined by an index cancer at one of 16 organ sites with at least 50 cases from each study population. We then investigated whether variants identified from exome sequencing were associated with these sets of multiple cancer cases in comparison to individuals with one and, separately, no cancers. RESULTS: We identified 22 variant-phenotype associations, 10 of which have not been previously discovered and were significantly overrepresented among individuals with multiple cancers, compared to those with a single cancer. CONCLUSIONS: Overall, we describe variants and genes that may play a fundamental role in the development of multiple primary cancers and improve our understanding of shared mechanisms underlying carcinogenesis.

Genome-wide association study of actinic keratosis identifies new susceptibility loci implicated in pigmentation and immune regulation pathways.

Actinic keratosis (AK) is a common precancerous cutaneous neoplasm that arises on chronically sun-exposed skin. AK susceptibility has a moderate genetic component, and although a few susceptibility loci have been identified, including IRF4, TYR, and MC1R, additional loci have yet to be discovered. We conducted a genome-wide association study of AK in non-Hispanic white participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (n = 63,110, discovery cohort), with validation in the Mass-General Brigham (MGB) Biobank cohort (n = 29,130). We identified eleven loci (P < 5 × 10-8), including seven novel loci, of which four novel loci were validated. In a meta-analysis (GERA + MGB), one additional novel locus, TRPS1, was identified. Genes within the identified loci are implicated in pigmentation (SLC45A2, IRF4, BNC2, TYR, DEF8, RALY, HERC2, and TRPS1), immune regulation (FOXP1 and HLA-DQA1), and cell signaling and tissue remodeling (MMP24) pathways. Our findings provide novel insight into the genetics and pathogenesis of AK susceptibility.

Exome sequencing and analysis of 454,787 UK Biobank participants.

A major goal in human genetics is to use natural variation to understand the phenotypic consequences of altering each protein-coding gene in the genome. Here we used exome sequencing1 to explore protein-altering variants and their consequences in 454,787 participants in the UK Biobank study2. We identified 12 million coding variants, including around 1 million loss-of-function and around 1.8 million deleterious missense variants. When these were tested for association with 3,994 health-related traits, we found 564 genes with trait associations at P ≤ 2.18 × 10-11. Rare variant associations were enriched in loci from genome-wide association studies (GWAS), but most (91%) were independent of common variant signals. We discovered several risk-increasing associations with traits related to liver disease, eye disease and cancer, among others, as well as risk-lowering associations for hypertension (SLC9A3R2), diabetes (MAP3K15, FAM234A) and asthma (SLC27A3). Six genes were associated with brain imaging phenotypes, including two involved in neural development (GBE1, PLD1). Of the signals available and powered for replication in an independent cohort, 81% were confirmed; furthermore, association signals were generally consistent across individuals of European, Asian and African ancestry. We illustrate the ability of exome sequencing to identify gene-trait associations, elucidate gene function and pinpoint effector genes that underlie GWAS signals at scale.

A large multiethnic GWAS meta-analysis of cataract identifies new risk loci and sex-specific effects.

Cataract is the leading cause of blindness among the elderly worldwide and cataract surgery is one of the most common operations performed in the United States. As the genetic etiology of cataract formation remains unclear, we conducted a multiethnic genome-wide association meta-analysis, combining results from the GERA and UK Biobank cohorts, and tested for replication in the 23andMe research cohort. We report 54 genome-wide significant loci, 37 of which were novel. Sex-stratified analyses identified CASP7 as an additional novel locus specific to women. We show that genes within or near 80% of the cataract-associated loci are significantly expressed and/or enriched-expressed in the mouse lens across various spatiotemporal stages as per iSyTE analysis. Furthermore, iSyTE shows 32 candidate genes in the associated loci have altered gene expression in 9 different gene perturbation mouse models of lens defects/cataract, suggesting their relevance to lens biology. Our work provides further insight into the complex genetic architecture of cataract susceptibility.

Cigarette smoking behaviors and the importance of ethnicity and genetic ancestry.

Cigarette smoking contributes to numerous diseases and is one of the leading causes of death in the United States. Smoking behaviors vary widely across race/ethnicity, but it is not clear why. Here, we examine the contribution of genetic ancestry to variation in two smoking-related traits in 43,485 individuals from four race/ethnicity groups (non-Hispanic white, Hispanic/Latino, East Asian, and African American) from a single U.S. healthcare plan. Smoking prevalence was the lowest among East Asians (22.7%) and the highest among non-Hispanic whites (38.5%). We observed significant associations between genetic ancestry and smoking-related traits. Within East Asians, we observed higher smoking prevalence with greater European (versus Asian) ancestry (P = 9.95 × 10-12). Within Hispanic/Latinos, higher cigarettes per day (CPD) was associated with greater European ancestry (P = 3.34 × 10-25). Within non-Hispanic whites, the lowest number of CPD was observed for individuals of southeastern European ancestry (P = 9.06 × 10-5). These associations remained after considering known smoking-associated loci, education, socioeconomic factors, and marital status. Our findings support the role of genetic ancestry and socioeconomic factors in cigarette smoking behaviors in non-Hispanic whites, Hispanic/Latinos, and East Asians.

Cross-cancer evaluation of polygenic risk scores for 16 cancer types in two large cohorts.

Even distinct cancer types share biological hallmarks. Here, we investigate polygenic risk score (PRS)-specific pleiotropy across 16 cancers in European ancestry individuals from the Genetic Epidemiology Research on Adult Health and Aging cohort (16,012 cases, 50,552 controls) and UK Biobank (48,969 cases, 359,802 controls). Within cohorts, each PRS is evaluated in multivariable logistic regression models against all other cancer types. Results are then meta-analyzed across cohorts. Ten positive and one inverse cross-cancer associations are found after multiple testing correction. Two pairs show bidirectional associations; the melanoma PRS is positively associated with oral cavity/pharyngeal cancer and vice versa, whereas the lung cancer PRS is positively associated with oral cavity/pharyngeal cancer, and the oral cavity/pharyngeal cancer PRS is inversely associated with lung cancer. Overall, we validate known, and uncover previously unreported, patterns of pleiotropy that have the potential to inform investigations of risk prediction, shared etiology, and precision cancer prevention strategies.

A Large-Scale Association Study Detects Novel Rare Variants, Risk Genes, Functional Elements, and Polygenic Architecture of Prostate Cancer Susceptibility.

To identify rare variants associated with prostate cancer susceptibility and better characterize the mechanisms and cumulative disease risk associated with common risk variants, we conducted an integrated study of prostate cancer genetic etiology in two cohorts using custom genotyping microarrays, large imputation reference panels, and functional annotation approaches. Specifically, 11,984 men (6,196 prostate cancer cases and 5,788 controls) of European ancestry from Northern California Kaiser Permanente were genotyped and meta-analyzed with 196,269 men of European ancestry (7,917 prostate cancer cases and 188,352 controls) from the UK Biobank. Three novel loci, including two rare variants (European ancestry minor allele frequency < 0.01, at 3p21.31 and 8p12), were significant genome wide in a meta-analysis. Gene-based rare variant tests implicated a known prostate cancer gene (HOXB13), as well as a novel candidate gene (ILDR1), which encodes a receptor highly expressed in prostate tissue and is related to the B7/CD28 family of T-cell immune checkpoint markers. Haplotypic patterns of long-range linkage disequilibrium were observed for rare genetic variants at HOXB13 and other loci, reflecting their evolutionary history. In addition, a polygenic risk score (PRS) of 188 prostate cancer variants was strongly associated with risk (90th vs. 40th-60th percentile OR = 2.62, P = 2.55 × 10-191). Many of the 188 variants exhibited functional signatures of gene expression regulation or transcription factor binding, including a 6-fold difference in log-probability of androgen receptor binding at the variant rs2680708 (17q22). Rare variant and PRS associations, with concomitant functional interpretation of risk mechanisms, can help clarify the full genetic architecture of prostate cancer and other complex traits. SIGNIFICANCE: This study maps the biological relationships between diverse risk factors for prostate cancer, integrating different functional datasets to interpret and model genome-wide data from over 200,000 men with and without prostate cancer.See related commentary by Lachance, p. 1637.

Genetic ancestry, skin pigmentation, and the risk of cutaneous squamous cell carcinoma in Hispanic/Latino and non-Hispanic white populations.

Although cutaneous squamous cell carcinoma (cSCC) is one of the most common malignancies in individuals of European ancestry, the incidence of cSCC in Hispanic/Latinos is also increasing. cSCC has both a genetic and environmental etiology. Here, we examine the role of genetic ancestry, skin pigmentation, and sun exposure in Hispanic/Latinos and non-Hispanic whites on cSCC risk. We observe an increased cSCC risk with greater European ancestry (P = 1.27 × 10-42) within Hispanic/Latinos and with greater northern (P = 2.38 × 10-65) and western (P = 2.28 × 10-49) European ancestry within non-Hispanic whites. These associations are significantly, but not completely, attenuated after considering skin pigmentation-associated loci, history of actinic keratosis, and sun-protected versus sun-exposed anatomical sites. We also report an association of the well-known pigment variant Ala111Thr (rs1426654) at SLC24A5 with cSCC in Hispanic/Latinos. These findings demonstrate a strong correlation of northwestern European genetic ancestry with cSCC risk in both Hispanic/Latinos and non-Hispanic whites, largely but not entirely mediated through its impact on skin pigmentation.

Identification of 31 loci for mammographic density phenotypes and their associations with breast cancer risk.

Mammographic density (MD) phenotypes are strongly associated with breast cancer risk and highly heritable. In this GWAS meta-analysis of 24,192 women, we identify 31 MD loci at P < 5 × 10-8, tripling the number known to 46. Seventeen identified MD loci also are associated with breast cancer risk in an independent meta-analysis (P < 0.05). Mendelian randomization analyses show that genetic estimates of dense area (DA), nondense area (NDA), and percent density (PD) are all significantly associated with breast cancer risk (P < 0.05). Pathway analyses reveal distinct biological processes involving DA, NDA and PD loci. These findings provide additional insights into the genetic basis of MD phenotypes and their associations with breast cancer risk.

Pan-cancer study detects genetic risk variants and shared genetic basis in two large cohorts.

Deciphering the shared genetic basis of distinct cancers has the potential to elucidate carcinogenic mechanisms and inform broadly applicable risk assessment efforts. Here, we undertake genome-wide association studies (GWAS) and comprehensive evaluations of heritability and pleiotropy across 18 cancer types in two large, population-based cohorts: the UK Biobank (408,786 European ancestry individuals; 48,961 cancer cases) and the Kaiser Permanente Genetic Epidemiology Research on Adult Health and Aging cohorts (66,526 European ancestry individuals; 16,001 cancer cases). The GWAS detect 21 genome-wide significant associations independent of previously reported results. Investigations of pleiotropy identify 12 cancer pairs exhibiting either positive or negative genetic correlations; 25 pleiotropic loci; and 100 independent pleiotropic variants, many of which are regulatory elements and/or influence cross-tissue gene expression. Our findings demonstrate widespread pleiotropy and offer further insight into the complex genetic architecture of cross-cancer susceptibility.

Genetic and environmental factors underlying keratinocyte carcinoma risk.

Recent large-scale GWAS and large epidemiologic studies have accelerated the discovery of genes and environmental factors that contribute to the risk of keratinocyte carcinoma (KC), which includes basal cell carcinoma (BCC) and squamous cell carcinoma (SCC). This Review summarizes the genomic regions associated with SCC and BCC risk, examines the genetic overlap between SCC and BCC, and discusses biological pathways involved in SCC and BCC development. Next, we review environmental factors that are associated with KC risk, including those that are shared between SCC and BCC as well as others that associated with only one type of KC. We conclude with a critical appraisal of current research and potential directions for future research.

Functional validity, role, and implications of heavy alcohol consumption genetic loci.

High alcohol consumption is a risk factor for morbidity and mortality, yet few genetic loci have been robustly associated with alcohol intake. Here, we use U.K. Biobank (n = 125,249) and GERA (n = 47,967) datasets to determine genetic factors associated with extreme population-level alcohol consumption and examine the functional validity of outcomes using model organisms and in silico techniques. We identified six loci attaining genome-wide significant association with alcohol consumption after meta-analysis and meeting our criteria for replication: ADH1B (lead SNP: rs1229984), KLB (rs13130794), BTF3P13 (rs144198753), GCKR (rs1260326), SLC39A8 (rs13107325), and DRD2 (rs11214609). A conserved role in phenotypic responses to alcohol was observed for all genetic targets available for investigation (ADH1B, GCKR, SLC39A8, and KLB) in Caenorhabditis elegans. Evidence of causal links to lung cancer, and shared genetic architecture with gout and hypertension was also found. These findings offer insight into genes, pathways, and relationships for disease risk associated with high alcohol consumption.

Genome-Wide Association and Functional Studies Reveal Novel Pharmacological Mechanisms for Allopurinol.

Allopurinol, which lowers uric acid (UA) concentration, is increasingly being recognized for its benefits in cardiovascular and renal disease. However, response to allopurinol is variable. We gathered samples from 4,446 multiethnic subjects for a genome-wide association study of allopurinol response. Consistent with previous studies, we observed that the Q141K variant in ABCG2 (rs2231142), which encodes the efflux pump breast cancer resistance protein (BCRP), associated with worse response to allopurinol. However, for the first time this association reached genome-wide level significance (P = 8.06 × 10-11 ). Additionally, we identified a novel association with a variant in GREM2 (rs1934341, P = 3.22 × 10-6 ). In vitro studies identified oxypurinol, the active metabolite of allopurinol, as an inhibitor of the UA transporter GLUT9, suggesting that oxypurinol may modulate UA reabsorption. These results provide strong evidence for a role of BCRP Q141K in allopurinol response, and suggest that allopurinol may have additional hypouricemic effects beyond xanthine oxidase inhibition.

Targeted Genomic Profiling of Acral Melanoma.

BACKGROUND: Acral melanoma is a rare type of melanoma that affects world populations irrespective of skin color and has worse survival than other cutaneous melanomas. It has relatively few single nucleotide mutations without the UV signature of cutaneous melanomas, but instead has a genetic landscape characterized by structural rearrangements and amplifications. BRAF mutations are less common than in other cutaneous melanomas, and knowledge about alternative therapeutic targets is incomplete. METHODS: To identify alternative therapeutic targets, we performed targeted deep-sequencing on 122 acral melanomas. We confirmed the loss of the tumor suppressors p16 and NF1 by immunohistochemistry in select cases. RESULTS: In addition to BRAF (21.3%), NRAS (27.9%), and KIT (11.5%) mutations, we identified a broad array of MAPK pathway activating alterations, including fusions of BRAF (2.5%), NTRK3 (2.5%), ALK (0.8%), and PRKCA (0.8%), which can be targeted by available inhibitors. Inactivation of NF1 occurred in 18 cases (14.8%). Inactivation of the NF1 cooperating factor SPRED1 occurred in eight cases (6.6%) as an alternative mechanism of disrupting the negative regulation of RAS. Amplifications recurrently affected narrow loci containing PAK1 and GAB2 (n = 27, 22.1%), CDK4 (n = 27, 22.1%), CCND1 (n = 24, 19.7%), EP300 (n = 20, 16.4%), YAP1 (n = 15, 12.3%), MDM2 (n = 13, 10.7%), and TERT (n = 13, 10.7%) providing additional and possibly complementary therapeutic targets. Acral melanomas with BRAFV600E mutations harbored fewer genomic amplifications and were more common in patients with European ancestry. CONCLUSION: Our findings support a new, molecularly based subclassification of acral melanoma with potential therapeutic implications: BRAFV600E mutant acral melanomas with characteristics similar to nonacral melanomas that could benefit from BRAF inhibitor therapy, and non-BRAFV600E mutant acral melanomas. Acral melanomas without BRAFV600E mutations harbor a broad array of therapeutically relevant alterations. Expanded molecular profiling would increase the detection of potentially targetable alterations for this subtype of acral melanoma.