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.

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

Breast cancer includes several subtypes with distinct characteristic biological, pathological, and clinical features. Elucidating subtype-specific genetic etiology could provide insights into the heterogeneity of breast cancer to facilitate development of improved prevention and treatment approaches. Here, we conducted pairwise case-case comparisons among five breast cancer subtypes by applying a case-case GWAS (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 to luminal cancer. Functional investigation further identified rs16867605 at 8q22.3 as a SNP that modulates enhancer activity of GRHL2. Subtype-informative polygenic risk scores (PRS) were derived, and patients with a high subtype-informative PRS had an up to 2-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 breast cancer patients. Overall, these findings have advanced our understanding of the genetic etiology of breast cancer subtypes, particularly for TNBC.

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.

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.

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. Herein, we conducted large alternative polyadenylation-wide association studies (APA-WAS) to investigate associations of APA levels with 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, including breast, ovary, prostate, colorectum, lung, and pancreas. 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 co-localization 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 post-transcriptional activities of their target genes compared to 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.

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.

Large-scale alternative polyadenylation (APA)-wide association studies to identify putative susceptibility genes in human common cancers.

Alternative polyadenylation (APA) modulates mRNA processing in the 3' untranslated regions (3'UTR), which affect mRNA stability and translation efficiency. Here, we build genetic models to predict APA levels in multiple tissues using sequencing data of 1,337 samples from the Genotype-Tissue Expression, and apply these models to assess associations between genetically predicted APA levels and cancer risk with data from large genome-wide association studies of six common cancers, including breast, ovary, prostate, colorectum, lung, and pancreas among European-ancestry populations. At a Bonferroni-corrected P □<□0.05, we identify 58 risk genes, including seven in newly identified loci. Using luciferase reporter assays, we demonstrate that risk alleles of 3'UTR variants, rs324015 ( STAT6 ), rs2280503 ( DIP2B ), rs1128450 ( FBXO38 ) and rs145220637 ( LDAH ), could significantly increase post-transcriptional activities of their target genes compared to reference alleles. Further gene knockdown experiments confirm their oncogenic roles. Our study provides additional insight into the genetic susceptibility of these common cancers.

Integrating genomics and proteomics data to identify candidate plasma biomarkers for lung cancer risk among European descendants.

BACKGROUND: Plasma proteins are potential biomarkers for complex diseases. We aimed to identify plasma protein biomarkers for lung cancer. METHODS: We investigated genetically predicted plasma levels of 1130 proteins in association with lung cancer risk among 29,266 cases and 56,450 controls of European descent. For proteins significantly associated with lung cancer risk, we evaluated associations of genetically predicted expression of their coding genes with the risk of lung cancer. RESULTS: Nine proteins were identified with genetically predicted plasma levels significantly associated with overall lung cancer risk at a false discovery rate (FDR) of <0.05. Proteins C2, MICA, AIF1, and CTSH were associated with increased lung cancer risk, while proteins SFTPB, HLA-DQA2, MICB, NRP1, and GMFG were associated with decreased lung cancer risk. Stratified analyses by histological types revealed the cross-subtype consistency of these nine associations and identified an additional protein, ICAM5, significantly associated with lung adenocarcinoma risk (FDR < 0.05). Coding genes of NRP1 and ICAM5 proteins are located at two loci that have never been reported by previous GWAS. Genetically predicted blood levels of genes C2, AIF1, and CTSH were associated with lung cancer risk, in directions consistent with those shown in protein-level analyses. CONCLUSION: Identification of novel plasma protein biomarkers provided new insights into the biology of lung cancer.

Identification of target proteins for breast cancer genetic risk loci and blood risk biomarkers in a large study by integrating genomic and proteomic data.

Genome-wide association studies (GWAS) have identified around 200 loci associated with breast cancer risk. However, protein targets for these loci remain largely unknown. Identifying protein targets and biomarkers can improve the understanding of cancer biology and etiology and identify high-risk individuals for cancer prevention. In this study, we investigated genetically predicted levels of 1142 circulating proteins with breast cancer risk in 133 384 cases and 113 789 controls of European ancestry included in the Breast Cancer Association Consortium (BCAC). We identified 22 blood protein biomarkers associated with the risk of overall breast cancer at a false discovery rate (FDR) <0.05, including nine proteins encoded by genes located at least 500 kb away from previously reported risk variants for breast cancer. Analyses focusing on 124 encoding genes located at GWAS-identified breast cancer risk loci found 20 proteins associated with overall breast cancer risk and one protein associated with triple-negative breast cancer risk at FDR <0.05. Adjustment for the GWAS-identified risk variants significantly attenuated the association for 13 of these proteins, suggesting that these proteins may be the targets of these GWAS-identified risk loci. The identified proteins are involved in various biological processes, including glutathione conjugation, STAT5 signaling and NF-κB signaling pathways. Our study identified novel protein targets and risk biomarkers for breast cancer risk.

Serum Lipid Profiles and Cholesterol-Lowering Medication Use in Relation to Subsequent Risk of Colorectal Cancer in the UK Biobank Cohort.

BACKGROUND: Dyslipidemia is closely associated with metabolic syndrome, a known risk factor for colorectal cancer. However, the association of dyslipidemia with colorectal cancer risk is controversial. Most previous studies did not consider cholesterol-lowering medication use at the time of lipid measurements, which could bias findings. METHODS: We analyzed data from 384,862 UK Biobank participants to disentangle the associations between blood lipids and colorectal cancer risk. Serum levels of total cholesterol, high- and low-density lipoprotein cholesterol (HDL-C, LDL-C), and triglyceride were measured at study baseline. Multivariable-adjusted Cox models were used to estimate HRs and 95% confidence intervals (CI). RESULTS: During a median follow-up time of 8.2 years, 3,150 incident primary colorectal cancer cases were identified. Triglyceride levels were positively, while HDL-C levels were inversely associated with colorectal cancer risk (both Ptrend < 0.005). No significant associations were found for total cholesterol and LDL-C. However, among nonusers of cholesterol-lowering medications, a high total cholesterol level (> 6.7 mmol/L, HR = 1.11; 95% CI, 1.00-1.24) and LDL-C level (>4.1 mmol/L, HR = 1.11; 95% CI, 0.99-1.23) was associated with an increased colorectal cancer risk compared with the referent group (5.2-6.2 mmol/L and 2.6-3.4 mmol/L for total and LDL cholesterol, respectively). Compared with nonusers, cholesterol-lowering medication users had 15% increased colorectal cancer risk (HR = 1.15; 95% CI, 1.04-1.26). CONCLUSIONS: Circulating total cholesterol, LDL-C, HDL-C and triglyceride were modestly associated with colorectal cancer risk. IMPACT: Our findings call for careful consideration of cholesterol-lowering medication use in future studies of blood lipid-colorectal cancer associations.

Genome- and transcriptome-wide association studies of 386,000 Asian and European-ancestry women provide new insights into breast cancer genetics.

By combining data from 160,500 individuals with breast cancer and 226,196 controls of Asian and European ancestry, we conducted genome- and transcriptome-wide association studies of breast cancer. We identified 222 genetic risk loci and 137 genes that were associated with breast cancer risk at a p < 5.0 × 10-8 and a Bonferroni-corrected p < 4.6 × 10-6, respectively. Of them, 32 loci and 15 genes showed a significantly different association between ER-positive and ER-negative breast cancer after Bonferroni correction. Significant ancestral differences in risk variant allele frequencies and their association strengths with breast cancer risk were identified. Of the significant associations identified in this study, 17 loci and 14 genes are located 1Mb away from any of the previously reported breast cancer risk variants. Pathways analyses including 221 putative risk genes identified multiple signaling pathways that may play a significant role in the development of breast cancer. Our study provides a comprehensive understanding of and new biological insights into the genetics of this common malignancy.

Incorporating Polygenic Risk Scores and Nongenetic Risk Factors for Breast Cancer Risk Prediction Among Asian Women.

Importance: Polygenic risk scores (PRSs) have shown promise in breast cancer risk prediction; however, limited studies have been conducted among Asian women. Objective: To develop breast cancer risk prediction models for Asian women incorporating PRSs and nongenetic risk factors. Design, Setting, and Participants: This diagnostic study included women of Asian ancestry from the Asia Breast Cancer Consortium. PRSs were developed using data from genomewide association studies (GWASs) of breast cancer conducted among 123 041 women with Asian ancestry (including 18 650 women with breast cancer) using 3 approaches: (1) reported PRS for women with European ancestry; (2) breast cancer-associated single-nucleotide variations (SNVs) identified by fine-mapping of GWAS-identified risk loci; and (3) genomewide risk prediction algorithms. A nongenetic risk score (NGRS) was built, including 7 well-established nongenetic risk factors, using data of 416 case participants and 1558 control participants from a prospective cohort study. PRSs were initially validated in an independent data set including 1426 case participants and 1323 control participants and further evaluated, along with the NGRS, in the second data set including 368 case participants and 736 control participants nested within a prospective cohort study. Main Outcomes and Measures: Logistic regression was used to examine associations of risk scores with breast cancer risk to estimate odds ratios (ORs) with 95% CIs and area under the receiver operating characteristic curve (AUC). Results: A total of 126 894 women of Asian ancestry were included; 20 444 (16.1%) had breast cancer. The mean (SD) age ranged from 49.1 (10.8) to 54.4 (10.4) years for case participants and 50.6 (9.5) to 54.0 (7.4) years for control participants among studies that provided demographic characteristics. In the prospective cohort, a PRS with 111 SNVs developed using the fine-mapping approach (PRS111) showed a prediction performance comparable with a genomewide PRS that included more than 855 000 SNVs. The OR per SD increase of PRS111 score was 1.67 (95% CI, 1.46-1.92), with an AUC of 0.639 (95% CI, 0.604-0.674). The NGRS had a limited predictive ability (AUC, 0.565; 95% CI, 0.529-0.601). Compared with the average risk group (40th-60th percentile), women in the top 5% of PRS111 and NGRS were at a 3.84-fold (95% CI, 2.30-6.46) and 2.10-fold (95% CI, 1.22-3.62) higher risk of breast cancer, respectively. The prediction model including both PRS111 and NGRS achieved the highest prediction accuracy (AUC, 0.648; 95% CI, 0.613-0.682). Conclusions and Relevance: In this study, PRSs derived using breast cancer risk-associated SNVs had similar predictive performance in Asian and European women. Including nongenetic risk factors in models further improved prediction accuracy. These findings support the utility of these models in developing personalized screening and prevention strategies.

Associations of genetic susceptibility to 16 cancers with risk of breast cancer overall and by intrinsic subtypes.

Certain genetic variants are associated with risks of multiple cancers. We investigated breast cancer risk with overall genetic susceptibility to each of 16 other cancers. We constructed polygenic risk scores (PRS) for 16 cancers using risk variants identified by genome-wide association studies. We evaluated the associations of these PRSs with breast cancer risk (overall and by subtypes) using Breast Cancer Association Consortium data, including 106,278 cases and 91,477 controls of European ancestry. Odds ratios (OR) and 95% confidence intervals (CIs) were estimated to measure the association of each PRS with breast cancer risk. Data from the UK Biobank, including 4,337 cases and 209,983 non-cases, were used to replicate the findings. A 5%-8% significantly elevated risk of overall breast cancer was associated with per unit increase of the PRS for glioma and cancers of the corpus uteri, stomach, or colorectum. Analyses by subtype revealed that the PRS for corpus uteri cancer (OR = 1.09; 95% CI, 1.03-1.15) and stomach cancer (OR = 1.07; 95% CI, 1.03-1.12) were associated with estrogen receptor-positive breast cancer, while ovarian cancer PRS was associated with triple-negative breast cancer (OR = 1.25; 95% CI, 1.01-1.55). UK Biobank data supported the positive associations of overall breast cancer risk with PRS for melanoma and cancers of the stomach, colorectum, and ovary. Our study provides strong evidence for shared genetic susceptibility of breast cancer with several other cancers. Results from our study help uncover the genetic basis for breast and other cancers and identify individuals at high risk for multiple cancers.