Integrating muti-omics data to identify tissue-specific DNA methylation biomarkers for cancer risk.

The relationship between tissue-specific DNA methylation and cancer risk remains inadequately elucidated. Leveraging resources from the Genotype-Tissue Expression consortium, here we develop genetic models to predict DNA methylation at CpG sites across the genome for seven tissues and apply these models to genome-wide association study data of corresponding cancers, namely breast, colorectal, renal cell, lung, ovarian, prostate, and testicular germ cell cancers. At Bonferroni-corrected P < 0.05, we identify 4248 CpGs that are significantly associated with cancer risk, of which 95.4% (4052) are specific to a particular cancer type. Notably, 92 CpGs within 55 putative novel loci retain significant associations with cancer risk after conditioning on proximal signals identified by genome-wide association studies. Integrative multi-omics analyses reveal 854 CpG-gene-cancer trios, suggesting that DNA methylation at 309 distinct CpGs might influence cancer risk through regulating the expression of 205 unique cis-genes. These findings substantially advance our understanding of the interplay between genetics, epigenetics, and gene expression in cancer etiology.

Association of ESR1 Germline Variants with TP53 Somatic Variants in Breast Tumors in a Genome-wide Study.

In breast tumors, somatic mutation frequencies in TP53 and PIK3CA vary by tumor subtype and ancestry. Emerging data suggest tumor mutation status is associated with germline variants and genetic ancestry. We aimed to identify germline variants that are associated with somatic TP53 or PIK3CA mutation status in breast tumors. A genome-wide association study was conducted in 2,850 women of European ancestry with breast cancer using TP53 and PIK3CA mutation status (positive or negative) as well as specific functional categories [e.g., TP53 gain-of-function (GOF) and loss-of-function, PIK3CA activating] as phenotypes. Germline variants showing evidence of association were selected for validation analyses and tested in multiple independent datasets. Discovery association analyses found five variants associated with TP53 mutation status with P values <1 × 10-6 and 33 variants with P values <1 × 10-5. Forty-four variants were associated with PIK3CA mutation status with P values <1 × 10-5. In validation analyses, only variants at the ESR1 locus were associated with TP53 mutation status after multiple comparisons corrections. Combined analyses in European and Malaysian populations found ESR1 locus variants rs9383938 and rs9479090 associated with the presence of TP53 mutations overall (P values 2 × 10-11 and 4.6 × 10-10, respectively). rs9383938 also showed association with TP53 GOF mutations (P value 6.1 × 10-7). rs9479090 showed suggestive evidence (P value 0.02) for association with TP53 mutation status in African ancestry populations. No other variants were significantly associated with TP53 or PIK3CA mutation status. Larger studies are needed to confirm these findings and determine if additional variants contribute to ancestry-specific differences in mutation frequency. SIGNIFICANCE: Emerging data show ancestry-specific differences in TP53 and PIK3CA mutation frequency in breast tumors suggesting that germline variants may influence somatic mutational processes. This study identified variants near ESR1 associated with TP53 mutation status and identified additional loci with suggestive association which may provide biological insight into observed differences.

Dietary intake, obesity, and physical activity in association with biliary tract cancer risk: Results from meta-analyses of individual-level data from prospective cohort studies of 723,326 adults.

Biliary tract cancer (BTC) is a rare and aggressive malignancy with increasing incidence. Most BTC cases are diagnosed with metastatic disease which carries a 5-year survival rate of <5%. Physical activity, diet, and obesity might be associated with BTC risk, but studies have been limited particularly in African descendants. We addressed this knowledge gap by evaluating associations of BTC risk with obesity, physical activity, and dietary intakes in 723,326 adult participants in four cohort studies conducted in China, the United Kingdom, and the United States. Cox regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) in each cohort; results were combined using meta-analysis. All cohorts had ≥11 median follow-up years with 839 incident BTC cases combined. BTC risk was positively associated with body mass index (BMI) and waist-to-hip ratio (WHR) whereas physical activity, fruit intake, and fish intake were inversely associated. HR and (95% CI) comparing BMI >35.0 to 18.5-24.9: 1.71 (1.26, 2.31), p-trend <.0001; comparing BMI-adjusted WHR top to bottom quartile: 1.20 (0.94, 1.53), p-trend = .05; comparing ≥15-0 metabolic equivalent task-hours/week 0.76 (0.61, 0.94), p-trend = .009; comparing highest to lowest intake tertile for fruit and fish 0.79 (0.66, 0.95), p-trend = .01; 0.82 (0.68, 0.98), p-trend = .04, respectively. Associations were, in general, similar across ancestry groups. Our study provides strong evidence for important roles of obesity, diet, and physical activity in BTC etiology and stresses the need for lifestyle modification to combat the rising incidence of this fatal malignancy.

Case-Case Genome-Wide Analyses Identify Subtype-Informative Variants That Confer Risk for Breast Cancer.

Breast cancer includes several subtypes with distinct characteristic biological, pathologic, and clinical features. Elucidating subtype-specific genetic etiology could provide insights into the heterogeneity of breast cancer to facilitate the development of improved prevention and treatment approaches. In this study, we conducted pairwise case-case comparisons among five breast cancer subtypes by applying a case-case genome-wide association study (CC-GWAS) approach to summary statistics data of the Breast Cancer Association Consortium. The approach identified 13 statistically significant loci and eight suggestive loci, the majority of which were identified from comparisons between triple-negative breast cancer (TNBC) and luminal A breast cancer. Associations of lead variants in 12 loci remained statistically significant after accounting for previously reported breast cancer susceptibility variants, among which, two were genome-wide significant. Fine mapping implicated putative functional/causal variants and risk genes at several loci, e.g., 3q26.31/TNFSF10, 8q22.3/NACAP1/GRHL2, and 8q23.3/LINC00536/TRPS1, for TNBC as compared with luminal cancer. Functional investigation further identified rs16867605 at 8q22.3 as a SNP that modulates the enhancer activity of GRHL2. Subtype-informative polygenic risk scores (PRS) were derived, and patients with a high subtype-informative PRS had an up to two-fold increased risk of being diagnosed with TNBC instead of luminal cancers. The CC-GWAS PRS remained statistically significant after adjusting for TNBC PRS derived from traditional case-control GWAS in The Cancer Genome Atlas and the African Ancestry Breast Cancer Genetic Consortium. The CC-GWAS PRS was also associated with overall survival and disease-specific survival among patients with breast cancer. Overall, these findings have advanced our understanding of the genetic etiology of breast cancer subtypes, particularly for TNBC. Significance: The discovery of subtype-informative genetic risk variants for breast cancer advances our understanding of the etiologic heterogeneity of breast cancer, which could accelerate the identification of targets and personalized strategies for prevention and treatment.

Using genome and transcriptome data from African-ancestry female participants to identify putative breast cancer susceptibility genes.

African-ancestry (AA) participants are underrepresented in genetics research. Here, we conducted a transcriptome-wide association study (TWAS) in AA female participants to identify putative breast cancer susceptibility genes. We built genetic models to predict levels of gene expression, exon junction, and 3' UTR alternative polyadenylation using genomic and transcriptomic data generated in normal breast tissues from 150 AA participants and then used these models to perform association analyses using genomic data from 18,034 cases and 22,104 controls. At Bonferroni-corrected P < 0.05, we identified six genes associated with breast cancer risk, including four genes not previously reported (CTD-3080P12.3, EN1, LINC01956 and NUP210L). Most of these genes showed a stronger association with risk of estrogen-receptor (ER) negative or triple-negative than ER-positive breast cancer. We also replicated the associations with 29 genes reported in previous TWAS at P < 0.05 (one-sided), providing further support for an association of these genes with breast cancer risk. Our study sheds new light on the genetic basis of breast cancer and highlights the value of conducting research in AA populations.

Large-Scale Alternative Polyadenylation-Wide Association Studies to Identify Putative Cancer Susceptibility Genes.

Alternative polyadenylation (APA) modulates mRNA processing in the 3'-untranslated regions (3' UTR), affecting mRNA stability and translation efficiency. Research into genetically regulated APA has the potential to provide insights into cancer risk. In this study, we conducted large APA-wide association studies to investigate associations between APA levels and cancer risk. Genetic models were built to predict APA levels in multiple tissues using genotype and RNA sequencing data from 1,337 samples from the Genotype-Tissue Expression project. Associations of genetically predicted APA levels with cancer risk were assessed by applying the prediction models to data from large genome-wide association studies of six common cancers among European ancestry populations: breast, ovarian, prostate, colorectal, lung, and pancreatic cancers. A total of 58 risk genes (corresponding to 76 APA sites) were associated with at least one type of cancer, including 25 genes previously not linked to cancer susceptibility. Of the identified risk APAs, 97.4% and 26.3% were supported by 3'-UTR APA quantitative trait loci and colocalization analyses, respectively. Luciferase reporter assays for four selected putative regulatory 3'-UTR variants demonstrated that the risk alleles of 3'-UTR variants, rs324015 (STAT6), rs2280503 (DIP2B), rs1128450 (FBXO38), and rs145220637 (LDHA), significantly increased the posttranscriptional activities of their target genes compared with reference alleles. Furthermore, knockdown of the target genes confirmed their ability to promote proliferation and migration. Overall, this study provides insights into the role of APA in the genetic susceptibility to common cancers. Significance: Systematic evaluation of associations of alternative polyadenylation with cancer risk reveals 58 putative susceptibility genes, highlighting the contribution of genetically regulated alternative polyadenylation of 3'UTRs to genetic susceptibility to cancer.

Genome-wide association analyses of breast cancer in women of African ancestry identify new susceptibility loci and improve risk prediction.

We performed genome-wide association studies of breast cancer including 18,034 cases and 22,104 controls of African ancestry. Genetic variants at 12 loci were associated with breast cancer risk (P < 5 × 10-8), including associations of a low-frequency missense variant rs61751053 in ARHGEF38 with overall breast cancer (odds ratio (OR) = 1.48) and a common variant rs76664032 at chromosome 2q14.2 with triple-negative breast cancer (TNBC) (OR = 1.30). Approximately 15.4% of cases with TNBC carried six risk alleles in three genome-wide association study-identified TNBC risk variants, with an OR of 4.21 (95% confidence interval = 2.66-7.03) compared with those carrying fewer than two risk alleles. A polygenic risk score (PRS) showed an area under the receiver operating characteristic curve of 0.60 for the prediction of breast cancer risk, which outperformed PRS derived using data from females of European ancestry. Our study markedly increases the population diversity in genetic studies for breast cancer and demonstrates the utility of PRS for risk prediction in females of African ancestry.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Lifestyle factors, genetic susceptibility to obesity and their interactions on coronary artery disease risk: A cohort study in the UK Biobank.

OBJECTIVE: We aimed to evaluate potential modifying effects of genetic susceptibility to obesity on the association of lifestyle factors with coronary artery disease (CAD) risk. METHODS: A total of 328,606 participants (54% women) were included using data from the UK Biobank. We evaluated the risk of developing CAD associated with obesity-related polygenic scores (PGSs) and healthy lifestyle scores (HLSs). HLSs were constructed using six lifestyle factors. Obesity PGSs were created using genetic variants identified by genome-wide association studies, including 941 variants for body mass index (BMI) and 457 for waist-to-hip ratio (WHR). Both HLSs and PGSs were categorized into three groups. RESULTS: During a 9-year median follow-up, 14,541 participants developed CAD. An unhealthy lifestyle was significantly associated with an increased CAD risk (hazard ratio [HR] = 2.24, 95% confidence interval [CI] = 2.09-2.40). High BMI and WHR PGSs were each significantly associated with an increased CAD risk (HRBMI = 1.23, 1.17-1.29; HRWHR = 1.15, 1.09-1.21). Lifestyle factors explained 41% (95% CI = 38%-45%) of CAD, while genetic variants for BMI explained only 10% (7%-14%). Risks of CAD were increased with poorer HLS independent of obesity-related PGSs. Individuals with the most unhealthy lifestyle and highest BMI PGS had the highest risk of CAD risk (HR = 2.59, 95% CI = 2.26-2.97), compared with participants with the healthiest lifestyle and lowest BMI PGS. CONCLUSIONS: While the observational nature of the study precludes the establishment of causality, our study provides supports for a causal association between obesity and CAD risk and the importance of lifestyle modification in the prevention of CAD.

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia.

BACKGROUND: Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs. METHODS: We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG). RESULTS: In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16). CONCLUSION: The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

Evaluating interactions of polygenic risk scores and NAT2 genotypes with tobacco smoking in bladder cancer risk.

Tobacco smoking is the most important risk factor for bladder cancer. Previous studies have identified the N-acetyltransferase (NAT2) gene in association with bladder cancer risk. The NAT2 gene encodes an enzyme that metabolizes aromatic amines, carcinogens commonly found in tobacco smoke. In our study, we evaluated potential interactions of tobacco smoking with NAT2 genotypes and polygenic risk score (PRS) for bladder cancer, using data from the UK Biobank, a large prospective cohort study. We used Cox proportional hazards models to measure the strength of the association. The PRS was derived using genetic risk variants identified by genome-wide association studies for bladder cancer. With an average of 10.1 years of follow-up of 390 678 eligible participants of European descent, 769 incident bladder cancer cases were identified. Current smokers with a PRS in the highest tertile had a higher risk of developing bladder cancer (HR: 6.45, 95% CI: 4.51-9.24) than current smokers with a PRS in the lowest tertile (HR: 2.41, 95% CI: 1.52-3.84; P for additive interaction = <.001). A similar interaction was found for genetically predicted metabolizing NAT2 phenotype and tobacco smoking where current smokers with the slow NAT2 phenotype had an increased risk of developing bladder cancer (HR: 5.70, 95% CI: 2.64-12.30) than current smokers with the fast NAT2 phenotype (HR: 3.61, 95% CI: 1.14-11.37; P for additive interaction = .100). Our study provides support for considering both genetic and lifestyle risk factors in developing prevention measures for bladder cancer.

Association of ESR1 germline variants with TP53 somatic variants in breast tumors in a genome-wide study.

Background: In breast tumors, somatic mutation frequencies in TP53 and PIK3CA vary by tumor subtype and ancestry. HER2 positive and triple negative breast cancers (TNBC) have a higher frequency of TP53 somatic mutations than other subtypes. PIK3CA mutations are more frequently observed in hormone receptor positive tumors. Emerging data suggest tumor mutation status is associated with germline variants and genetic ancestry. We aimed to identify germline variants that are associated with somatic TP53 or PIK3CA mutation status in breast tumors. Methods: A genome-wide association study was conducted using breast cancer mutation status of TP53 and PIK3CA and functional mutation categories including TP53 gain of function (GOF) and loss of function mutations and PIK3CA activating/hotspot mutations. The discovery analysis consisted of 2850 European ancestry women from three datasets. Germline variants showing evidence of association with somatic mutations were selected for validation analyses based on predicted function, allele frequency, and proximity to known cancer genes or risk loci. Candidate variants were assessed for association with mutation status in a multi-ancestry validation study, a Malaysian study, and a study of African American/Black women with TNBC. Results: The discovery Germline x Mutation (GxM) association study found five variants associated with one or more TP53 phenotypes with P values <1×10-6, 33 variants associated with one or more TP53 phenotypes with P values <1×10-5, and 44 variants associated with one or more PIK3CA phenotypes with P values <1×10-5. In the multi-ancestry and Malaysian validation studies, germline ESR1 locus variant, rs9383938, was associated with the presence of TP53 mutations overall (P values 6.8×10-5 and 9.8×10-8, respectively) and TP53 GOF mutations (P value 8.4×10-6). Multiple variants showed suggestive evidence of association with PIK3CA mutation status in the validation studies, but none were significant after correction for multiple comparisons. Conclusions: We found evidence that germline variants were associated with TP53 and PIK3CA mutation status in breast cancers. Variants near the estrogen receptor alpha gene, ESR1, were significantly associated with overall TP53 mutations and GOF mutations. Larger multi-ancestry studies are needed to confirm these findings and determine if these variants contribute to ancestry-specific differences in mutation frequency.