Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers.

Squamous cell carcinomas (SqCC) of the aerodigestive tract have similar etiological risk factors. Although genetic risk variants for individual cancers have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. To identify novel and pleotropic SqCC risk variants, we performed a meta-analysis of GWAS data on lung SqCC (LuSqCC), oro/pharyngeal SqCC (OSqCC), laryngeal SqCC (LaSqCC) and esophageal SqCC (ESqCC) cancers, totaling 13,887 cases and 61,961 controls of European ancestry. We identified one novel genome-wide significant (Pmeta<5x10-8) aerodigestive SqCC susceptibility loci in the 2q33.1 region (rs56321285, TMEM273). Additionally, three previously unknown loci reached suggestive significance (Pmeta<5x10-7): 1q32.1 (rs12133735, near MDM4), 5q31.2 (rs13181561, TMEM173) and 19p13.11 (rs61494113, ABHD8). Multiple previously identified loci for aerodigestive SqCC also showed evidence of pleiotropy in at least another SqCC site, these include: 4q23 (ADH1B), 6p21.33 (STK19), 6p21.32 (HLA-DQB1), 9p21.33 (CDKN2B-AS1) and 13q13.1(BRCA2). Gene-based association and gene set enrichment identified a set of 48 SqCC-related genes rel to DNA damage and epigenetic regulation pathways. Our study highlights the importance of cross-cancer analyses to identify pleiotropic risk loci of histology-related cancers arising at distinct anatomical sites.

Integration of multiomic annotation data to prioritize and characterize inflammation and immune-related risk variants in squamous cell lung cancer.

Clinical trial results have recently demonstrated that inhibiting inflammation by targeting the interleukin-1ß pathway can offer a significant reduction in lung cancer incidence and mortality, highlighting a pressing and unmet need to understand the benefits of inflammation-focused lung cancer therapies at the genetic level. While numerous genome-wide association studies (GWAS) have explored the genetic etiology of lung cancer, there remains a large gap between the type of information that may be gleaned from an association study and the depth of understanding necessary to explain and drive translational findings. Thus, in this study we jointly model and integrate extensive multiomics data sources, utilizing a total of 40 genome-wide functional annotations that augment previously published results from the International Lung Cancer Consortium (ILCCO) GWAS, to prioritize and characterize single nucleotide polymorphisms (SNPs) that increase risk of squamous cell lung cancer through the inflammatory and immune responses. Our work bridges the gap between correlative analysis and translational follow-up research, refining GWAS association measures in an interpretable and systematic manner. In particular, reanalysis of the ILCCO data highlights the impact of highly associated SNPs from nuclear factor-κB signaling pathway genes as well as major histocompatibility complex mediated variation in immune responses. One consequence of prioritizing likely functional SNPs is the pruning of variants that might be selected for follow-up work by over an order of magnitude, from potentially tens of thousands to hundreds. The strategies we introduce provide informative and interpretable approaches for incorporating extensive genome-wide annotation data in analysis of genetic association studies.

Cancer-driving H3G34V/R/D mutations block H3K36 methylation and H3K36me3-MutSα interaction.

Somatic mutations on glycine 34 of histone H3 (H3G34) cause pediatric cancers, but the underlying oncogenic mechanism remains unknown. We demonstrate that substituting H3G34 with arginine, valine, or aspartate (H3G34R/V/D), which converts the non-side chain glycine to a large side chain-containing residue, blocks H3 lysine 36 (H3K36) dimethylation and trimethylation by histone methyltransferases, including SETD2, an H3K36-specific trimethyltransferase. Our structural analysis reveals that the H3 "G33-G34" motif is recognized by a narrow substrate channel, and that H3G34/R/V/D mutations impair the catalytic activity of SETD2 due to steric clashes that impede optimal SETD2-H3K36 interaction. H3G34R/V/D mutations also block H3K36me3 from interacting with mismatch repair (MMR) protein MutSα, preventing the recruitment of the MMR machinery to chromatin. Cells harboring H3G34R/V/D mutations display a mutator phenotype similar to that observed in MMR-defective cells. Therefore, H3G34R/V/D mutations promote genome instability and tumorigenesis by inhibiting MMR activity.

Evaluation of associations between genetically predicted circulating protein biomarkers and breast cancer risk.

A small number of circulating proteins have been reported to be associated with breast cancer risk, with inconsistent results. Herein, we attempted to identify novel protein biomarkers for breast cancer via the integration of genomics and proteomics data. In the Breast Cancer Association Consortium (BCAC), with 122,977 cases and 105,974 controls of European descendants, we evaluated the associations of the genetically predicted concentrations of >1,400 circulating proteins with breast cancer risk. We used data from a large-scale protein quantitative trait loci (pQTL) analysis as our study instrument. Summary statistics for these pQTL variants related to breast cancer risk were obtained from the BCAC and used to estimate odds ratios (OR) for each protein using the inverse-variance weighted method. We identified 56 proteins significantly associated with breast cancer risk by instrumental analysis (false discovery rate <0.05). Of these, the concentrations of 32 were influenced by variants close to a breast cancer susceptibility locus (ABO, 9q34.2). Many of these proteins, such as insulin receptor, insulin-like growth factor receptor 1 and other membrane receptors (OR: 0.82-1.18, p values: 6.96 × 10-4 -3.28 × 10-8 ), are linked to insulin resistance and estrogen receptor signaling pathways. Proteins identified at other loci include those involved in biological processes such as alcohol and lipid metabolism, proteolysis, apoptosis, immune regulation and cell motility and proliferation. Consistent associations were observed for 22 proteins in the UK Biobank data (p < 0.05). The study identifies potential novel biomarkers for breast cancer, but further investigation is needed to replicate our findings.

Genome-wide association study of INDELs identified four novel susceptibility loci associated with lung cancer risk.

Genome-wide association studies (GWAS) have identified 45 susceptibility loci associated with lung cancer. Only less than SNPs, small insertions and deletions (INDELs) are the second most abundant genetic polymorphisms in the human genome. INDELs are highly associated with multiple human diseases, including lung cancer. However, limited studies with large-scale samples have been available to systematically evaluate the effects of INDELs on lung cancer risk. Here, we performed a large-scale meta-analysis to evaluate INDELs and their risk for lung cancer in 23,202 cases and 19,048 controls. Functional annotations were performed to further explore the potential function of lung cancer risk INDELs. Conditional analysis was used to clarify the relationship between INDELs and SNPs. Four new risk loci were identified in genome-wide INDEL analysis (1p13.2: rs5777156, Insertion, OR = 0.92, p = 9.10 × 10-8 ; 4q28.2: rs58404727, Deletion, OR = 1.19, p = 5.25 × 10-7 ; 12p13.31: rs71450133, Deletion, OR = 1.09, p = 8.83 × 10-7 ; and 14q22.3: rs34057993, Deletion, OR = 0.90, p = 7.64 × 10-8 ). The eQTL analysis and functional annotation suggested that INDELs might affect lung cancer susceptibility by regulating the expression of target genes. After conducting conditional analysis on potential causal SNPs, the INDELs in the new loci were still nominally significant. Our findings indicate that INDELs could be potentially functional genetic variants for lung cancer risk. Further functional experiments are needed to better understand INDEL mechanisms in carcinogenesis.

Transcriptome-wide association study reveals candidate causal genes for lung cancer.

We have recently completed the largest GWAS on lung cancer including 29,266 cases and 56,450 controls of European descent. The goal of our study has been to integrate the complete GWAS results with a large-scale expression quantitative trait loci (eQTL) mapping study in human lung tissues (n = 1,038) to identify candidate causal genes for lung cancer. We performed transcriptome-wide association study (TWAS) for lung cancer overall, by histology (adenocarcinoma, squamous cell carcinoma and small cell lung cancer) and smoking subgroups (never- and ever-smokers). We performed replication analysis using lung data from the Genotype-Tissue Expression (GTEx) project. DNA damage assays were performed in human lung fibroblasts for selected TWAS genes. As expected, the main TWAS signal for all histological subtypes and ever-smokers was on chromosome 15q25. The gene most strongly associated with lung cancer at this locus using the TWAS approach was IREB2 (pTWAS = 1.09E-99), where lower predicted expression increased lung cancer risk. A new lung adenocarcinoma susceptibility locus was revealed on 9p13.3 and associated with higher predicted expression of AQP3 (pTWAS = 3.72E-6). Among the 45 previously described lung cancer GWAS loci, we mapped candidate target gene for 17 of them. The association AQP3-adenocarcinoma on 9p13.3 was replicated using GTEx (pTWAS = 6.55E-5). Consistent with the effect of risk alleles on gene expression levels, IREB2 knockdown and AQP3 overproduction promote endogenous DNA damage. These findings indicate genes whose expression in lung tissue directly influences lung cancer risk.

Impact of individual level uncertainty of lung cancer polygenic risk score (PRS) on risk stratification.

BACKGROUND: Although polygenic risk score (PRS) has emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS), the individual-level accuracy of lung cancer PRS and the extent to which its impact on subsequent clinical applications remains largely unexplored. METHODS: Lung cancer PRSs and confidence/credible interval (CI) were constructed using two statistical approaches for each individual: (1) the weighted sum of 16 GWAS-derived significant SNP loci and the CI through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes), among 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium. Individuals were classified into different genetic risk subgroups based on the relationship between their own PRS mean/PRS CI and the population level threshold. RESULTS: Considerable variances in PRS point estimates at the individual level were observed for both methods, with an average standard deviation (s.d.) of 0.12 for PRS-16-CV and a much larger s.d. of 0.88 for PRS-Bayes. Using PRS-16-CV, only 25.0% of individuals with PRS point estimates in the lowest decile of PRS and 16.8% in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individuals. Only 19% of the individuals were concordantly identified as having high genetic risk (> 90th percentile) using the two PRS estimators. An increased relative risk of lung cancer comparing the highest PRS percentile to the lowest was observed when taking the CI into account (OR = 2.73, 95% CI: 2.12-3.50, P-value = 4.13 × 10-15) compared to using PRS-16-CV mean (OR = 2.23, 95% CI: 1.99-2.49, P-value = 5.70 × 10-46). Improved risk prediction performance with higher AUC was consistently observed in individuals identified by PRS-16-CV CI, and the best performance was achieved by incorporating age, gender, and detailed smoking pack-years (AUC: 0.73, 95% CI = 0.72-0.74). CONCLUSIONS: Lung cancer PRS estimates using different methods have modest correlations at the individual level, highlighting the importance of considering individual-level uncertainty when evaluating the practical utility of PRS.

A Large-Scale Genome-Wide Gene-Gene Interaction Study of Lung Cancer Susceptibility in Europeans With a Trans-Ethnic Validation in Asians.

INTRODUCTION: Although genome-wide association studies have been conducted to investigate genetic variation of lung tumorigenesis, little is known about gene-gene (G × G) interactions that may influence the risk of non-small cell lung cancer (NSCLC). METHODS: Leveraging a total of 445,221 European-descent participants from the International Lung Cancer Consortium OncoArray project, Transdisciplinary Research in Cancer of the Lung and UK Biobank, we performed a large-scale genome-wide G × G interaction study on European NSCLC risk by a series of analyses. First, we used BiForce to evaluate and rank more than 58 billion G × G interactions from 340,958 single-nucleotide polymorphisms (SNPs). Then, the top interactions were further tested by demographically adjusted logistic regression models. Finally, we used the selected interactions to build lung cancer screening models of NSCLC, separately, for never and ever smokers. RESULTS: With the Bonferroni correction, we identified eight statistically significant pairs of SNPs, which predominantly appeared in the 6p21.32 and 5p15.33 regions (e.g., rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.17, p = 6.57 × 10-13; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.17, p = 2.43 × 10-13; rs2858859HLA-DQA1 and rs9275572HLA-DQA2, ORinteraction = 1.15, p = 2.84 × 10-13; rs2853668TERT and rs62329694CLPTM1L, ORinteraction = 0.73, p = 2.70 × 10-13). Notably, even with much genetic heterogeneity across ethnicities, three pairs of SNPs in the 6p21.32 region identified from the European-ancestry population remained significant among an Asian population from the Nanjing Medical University Global Screening Array project (rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.13, p = 0.008; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.11, p = 5.23 × 10-4; rs3135369BTNL2 and rs9271300HLA-DQA1, ORinteraction = 0.89, p = 0.006). The interaction-empowered polygenetic risk score that integrated classical polygenetic risk score and G × G information score was remarkable in lung cancer risk stratification. CONCLUSIONS: Important G × G interactions were identified and enriched in the 5p15.33 and 6p21.32 regions, which may enhance lung cancer screening models.

Comprehensive functional annotation of susceptibility variants identifies genetic heterogeneity between lung adenocarcinoma and squamous cell carcinoma.

Although genome-wide association studies have identified more than eighty genetic variants associated with non-small cell lung cancer (NSCLC) risk, biological mechanisms of these variants remain largely unknown. By integrating a large-scale genotype data of 15 581 lung adenocarcinoma (AD) cases, 8350 squamous cell carcinoma (SqCC) cases, and 27 355 controls, as well as multiple transcriptome and epigenomic databases, we conducted histology-specific meta-analyses and functional annotations of both reported and novel susceptibility variants. We identified 3064 credible risk variants for NSCLC, which were overrepresented in enhancer-like and promoter-like histone modification peaks as well as DNase I hypersensitive sites. Transcription factor enrichment analysis revealed that USF1 was AD-specific while CREB1 was SqCC-specific. Functional annotation and gene-based analysis implicated 894 target genes, including 274 specifics for AD and 123 for SqCC, which were overrepresented in somatic driver genes (ER = 1.95, P = 0.005). Pathway enrichment analysis and Gene-Set Enrichment Analysis revealed that AD genes were primarily involved in immune-related pathways, while SqCC genes were homologous recombination deficiency related. Our results illustrate the molecular basis of both well-studied and new susceptibility loci of NSCLC, providing not only novel insights into the genetic heterogeneity between AD and SqCC but also a set of plausible gene targets for post-GWAS functional experiments.

The underrecognized progressive nature of N370S Gaucher disease and assessment of cancer risk in 403 patients.

Mutations in GBA1 gene that encodes lysosomal glucocerebrosidase result in Type 1 Gaucher Disease (GD), the commonest lysosomal storage disorder; the most prevalent disease mutation is N370S. We investigated the heterogeneity and natural course of N370S GD in 403 patients. Demographic, clinical, and genetic characteristics of GD at presentation were examined in a cross-sectional study. In addition, the relative risk (RR) of cancer in patients compared with age-, sex-, and ethnic-group adjusted national rates of cancer was determined. Of the 403 patients, 54% of patients were homozygous (N370S/N370S) and 46% were compound heterozygous for the N370S mutation (N370S/other). The majority of N370S/N370S patients displayed a phenotype characterized by late onset, predominantly skeletal disease, whereas the majority of N370S/other patients displayed early onset, predominantly visceral/hematologic disease, P < 0.0001. There was a striking increase in lifetime risk of multiple myeloma in the entire cohort (RR 25, 95% CI 9.17-54.40), mostly confined to N370S homozygous patients. The risk of other hematologic malignancies (RR 3.45, 95% CI 1.49-6.79), and overall cancer risk (RR 1.80, 95% CI 1.32-2.40) was increased. Homozygous N370S GD leads to adult-onset progressive skeletal disease with relative sparing of the viscera, a strikingly high risk of multiple myeloma, and an increased risk of other cancers. High incidence of gammopathy suggests an important role of the adaptive immune system in the development of GD. Adult patients with GD should be monitored for skeletal disease and cancers including multiple myeloma.

Fine mapping of MHC region in lung cancer highlights independent susceptibility loci by ethnicity.

Lung cancer has several genetic associations identified within the major histocompatibility complex (MHC); although the basis for these associations remains elusive. Here, we analyze MHC genetic variation among 26,044 lung cancer patients and 20,836 controls densely genotyped across the MHC, using the Illumina Illumina OncoArray or Illumina 660W SNP microarray. We impute sequence variation in classical HLA genes, fine-map MHC associations for lung cancer risk with major histologies and compare results between ethnicities. Independent and novel associations within HLA genes are identified in Europeans including amino acids in the HLA-B*0801 peptide binding groove and an independent HLA-DQB1*06 loci group. In Asians, associations are driven by two independent HLA allele sets that both increase risk in HLA-DQB1*0401 and HLA-DRB1*0701; the latter better represented by the amino acid Ala-104. These results implicate several HLA-tumor peptide interactions as the major MHC factor modulating lung cancer susceptibility.

Pleiotropy of genetic variants on obesity and smoking phenotypes: Results from the Oncoarray Project of The International Lung Cancer Consortium.

Obesity and cigarette smoking are correlated through complex relationships. Common genetic causes may contribute to these correlations. In this study, we selected 241 loci potentially associated with body mass index (BMI) based on the Genetic Investigation of ANthropometric Traits (GIANT) consortium data and calculated a BMI genetic risk score (BMI-GRS) for 17,037 individuals of European descent from the Oncoarray Project of the International Lung Cancer Consortium (ILCCO). Smokers had a significantly higher BMI-GRS than never-smokers (p = 0.016 and 0.010 before and after adjustment for BMI, respectively). The BMI-GRS was also positively correlated with pack-years of smoking (p<0.001) in smokers. Based on causal network inference analyses, seven and five of 241 SNPs were classified to pleiotropic models for BMI/smoking status and BMI/pack-years, respectively. Among them, three and four SNPs associated with smoking status and pack-years (p<0.05), respectively, were followed up in the ever-smoking data of the Tobacco, Alcohol and Genetics (TAG) consortium. Among these seven candidate SNPs, one SNP (rs11030104, BDNF) achieved statistical significance after Bonferroni correction for multiple testing, and three suggestive SNPs (rs13021737, TMEM18; rs11583200, ELAVL4; and rs6990042, SGCZ) achieved a nominal statistical significance. Our results suggest that there is a common genetic component between BMI and smoking, and pleiotropy analysis can be useful to identify novel genetic loci of complex phenotypes.

Breast cancer risk following bilateral oophorectomy in BRCA1 and BRCA2 mutation carriers: an international case-control study.

PURPOSE: The purpose of this study was to estimate the extent of protection offered against breast cancer by prophylactic oophorectomy in carriers of BRCA1 or BRCA2 mutations and to determine to what extent risk reduction varies with age at oophorectomy, age at diagnosis, and time elapsed since surgery. PATIENTS AND METHODS: We analyzed 1,439 patients with breast cancer and 1,866 matched controls derived from a registry of BRCA1 and BRCA2 carriers. We estimated odds ratios (ORs) of breast cancer for having had a bilateral oophorectomy, using conditional logistic regression, matched for parity and for oral contraceptive use. RESULTS: A previous history of oophorectomy was associated with a significant reduction in breast cancer risk of 56% for BRCA1 carriers (OR = 0.44; 95% CI, 0.29 to 0.66) and of 46% for BRCA2 carriers (OR = 0.57; 95% CI, 0.28 to 1.15). The risk reduction was greater if the oophorectomy was performed before age 40 (OR = 0.36; 95% CI, 0.20 to 0.64 for BRCA1 carriers) than after age 40 (OR = 0.53; 95% CI, 0.30 to 0.91). The protective effect was evident for 15 years post-oophorectomy (OR = 0.39; 95% CI, 0.26 to 0.57). CONCLUSION: Oophorectomy is an effective means of reducing the risk of breast cancer in carriers of BRCA1 mutations. The data suggest oophorectomy is protective in BRCA2 carriers as well, but needs to be confirmed in other studies.

A Germline Variant on Chromosome 4q31.1 Associates with Susceptibility to Developing Colon Cancer Metastasis.

We tested for germline variants showing association to colon cancer metastasis using a genome-wide association study that compared Ashkenazi Jewish individuals with stage IV metastatic colon cancers versus those with stage I or II non-metastatic colon cancers. In a two-stage study design, we demonstrated significant association to developing metastatic disease for rs60745952, that in Ashkenazi discovery and validation cohorts, respectively, showed an odds ratio (OR) = 2.3 (P = 2.73E-06) and OR = 1.89 (P = 8.05E-04) (exceeding validation threshold of 0.0044). Significant association to metastatic colon cancer was further confirmed by a meta-analysis of rs60745952 in these datasets plus an additional Ashkenazi validation cohort (OR = 1.92; 95% CI: 1.28-2.87), and by a permutation test that demonstrated a significantly longer haplotype surrounding rs60745952 in the stage IV samples. rs60745952, located in an intergenic region on chromosome 4q31.1, and not previously associated with cancer, is, thus, a germline genetic marker for susceptibility to developing colon cancer metastases among Ashkenazi Jews.