A mixed-model approach for powerful testing of genetic associations with cancer risk incorporating tumor characteristics

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

Cis- and trans-eQTL TWAS of breast and ovarian cancer identify more than 100 risk associated genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWAS) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have only considered regulatory effects of risk associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWAS of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents a first look into the role of trans-eQTLs in the complex molecular mechanisms underlying these diseases.

The sequence kernel association test for multicategorical outcomes

Disease heterogeneity is ubiquitous in biomedical and clinical studies. In genetic studies, researchers are increasingly interested in understanding the distinct genetic underpinning of subtypes of diseases. However, existing set-based analysis methods for genome-wide association studies are either inadequate or inefficient to handle such multicategorical outcomes. In this paper, we proposed a novel set-based association analysis method, sequence kernel association test (SKAT)-MC, the sequence kernel association test for multicategorical outcomes (nominal or ordinal), which jointly evaluates the relationship between a set of variants (common and rare) and disease subtypes. Through comprehensive simulation studies, we showed that SKAT-MC effectively preserves the nominal type I error rate while substantially increases the statistical power compared to existing methods under various scenarios. We applied SKAT-MC to the Polish breast cancer study (PBCS), and identified gene FGFR2 was significantly associated with estrogen receptor (ER)+ and ER- breast cancer subtypes. We also investigated educational attainment using UK Biobank data (N = 127 , 127 $N=127,127$ ) with SKAT-MC, and identified 21 significant genes in the genome. Consequently, SKAT-MC is a powerful and efficient analysis tool for genetic association studies with multicategorical outcomes. A freely distributed R package SKAT-MC can be accessed at https://github.com/Zhiwen-Owen-Jiang/SKATMC.

Nitrated Polycyclic Aromatic Hydrocarbon (Nitro-PAH) Signatures and Somatic Mutations in Diesel Exhaust-Exposed Bladder Tumors

BACKGROUND

Diesel exhaust is a complex mixture, including polycyclic aromatic hydrocarbons (PAH) and nitrated PAHs (nitro-PAH), many of which are potent mutagens and possible bladder carcinogens. To explore the association between diesel exposure and bladder carcinogenesis, we examined the relationship between exposure and somatic mutations and mutational signatures in bladder tumors.

METHODS

Targeted sequencing was conducted in bladder tumors from the New England Bladder Cancer Study. Using data on 797 cases and 1,418 controls, two-stage polytomous logistic regression was used to evaluate etiologic heterogeneity between bladder cancer subtypes and quantitative, lifetime estimates of respirable elemental carbon (REC), a surrogate for diesel exposure. Poisson regression was used to evaluate associations between REC and mutational signatures.

RESULTS

We observed significant heterogeneity in the diesel-bladder cancer risk relationship, with a strong positive association among cases with high-grade, nonmuscle invasive TP53-mutated tumors compared with controls [ORTop Tertile vs.Unexposed, 4.8; 95% confidence interval (CI), 2.2-10.5; Ptrend < 0.001; Pheterogeneity = 0.002]. In muscle-invasive tumors, we observed a positive association between diesel exposure and the nitro-PAH signatures of 1,6-dintropyrene (RR, 1.93; 95% CI, 1.28-2.92) and 3-nitrobenzoic acid (RR, 1.97; 95% CI, 1.33-2.92).

CONCLUSIONS

The relationship between diesel exhaust and bladder cancer was heterogeneous based on the presence of TP53 mutations in tumors, further supporting the link between PAH exposure and TP53 mutations in carcinogenesis. Future studies that can identify nitro-PAH signatures in exposed tumors are warranted to add human data supporting the link between diesel and bladder cancer.

IMPACT

This study provides additional insight into the etiology and possible mechanisms related to diesel exhaust-induced bladder cancer.

Pathology of Tumors Associated With Pathogenic Germline Variants in 9 Breast Cancer Susceptibility Genes

IMPORTANCE

Rare germline genetic variants in several genes are associated with increased breast cancer (BC) risk, but their precise contributions to different disease subtypes are unclear. This information is relevant to guidelines for gene panel testing and risk prediction.

OBJECTIVE

To characterize tumors associated with BC susceptibility genes in large-scale population- or hospital-based studies.

DESIGN, SETTING, AND PARTICIPANTS

The multicenter, international case-control analysis of the BRIDGES study included 42 680 patients and 46 387 control participants, comprising women aged 18 to 79 years who were sampled independently of family history from 38 studies. Studies were conducted between 1991 and 2016. Sequencing and analysis took place between 2016 and 2021.

EXPOSURES

Protein-truncating variants and likely pathogenic missense variants in ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53.

MAIN OUTCOMES AND MEASURES

The intrinsic-like BC subtypes as defined by estrogen receptor, progesterone receptor, and ERBB2 (formerly known as HER2) status, and tumor grade; morphology; size; stage; lymph node involvement; subtype-specific odds ratios (ORs) for carrying protein-truncating variants and pathogenic missense variants in the 9 BC susceptibility genes.

RESULTS

The mean (SD) ages at interview (control participants) and diagnosis (cases) were 55.1 (11.9) and 55.8 (10.6) years, respectively; all participants were of European or East Asian ethnicity. There was substantial heterogeneity in the distribution of intrinsic subtypes by gene. RAD51C, RAD51D, and BARD1 variants were associated mainly with triple-negative disease (OR, 6.19 [95% CI, 3.17-12.12]; OR, 6.19 [95% CI, 2.99-12.79]; and OR, 10.05 [95% CI, 5.27-19.19], respectively). CHEK2 variants were associated with all subtypes (with ORs ranging from 2.21-3.17) except for triple-negative disease. For ATM variants, the association was strongest for the hormone receptor (HR)+ERBB2- high-grade subtype (OR, 4.99; 95% CI, 3.68-6.76). BRCA1 was associated with increased risk of all subtypes, but the ORs varied widely, being highest for triple-negative disease (OR, 55.32; 95% CI, 40.51-75.55). BRCA2 and PALB2 variants were also associated with triple-negative disease. TP53 variants were most strongly associated with HR+ERBB2+ and HR-ERBB2+ subtypes. Tumors occurring in pathogenic variant carriers were of higher grade. For most genes and subtypes, a decline in ORs was observed with increasing age. Together, the 9 genes were associated with 27.3% of all triple-negative tumors in women 40 years or younger.

CONCLUSIONS AND RELEVANCE

The results of this case-control study suggest that variants in the 9 BC risk genes differ substantially in their associated pathology but are generally associated with triple-negative and/or high-grade disease. Knowing the age and tumor subtype distributions associated with individual BC genes can potentially aid guidelines for gene panel testing, risk prediction, and variant classification and guide targeted screening strategies.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

METHODS

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

RESULTS

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

CONCLUSION

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

TP53 Pathway Function, Estrogen Receptor Status, and Breast Cancer Risk Factors in the Carolina Breast Cancer Study

BACKGROUND

TP53 and estrogen receptor (ER) both play essential roles in breast cancer development and progression, with recent research revealing cross-talk between TP53 and ER signaling pathways. Although many studies have demonstrated heterogeneity of risk factor associations across ER subtypes, associations by TP53 status have been inconsistent.

METHODS

This case-case analysis included incident breast cancer cases (47% Black) from the Carolina Breast Cancer Study (1993-2013). Formalin-fixed paraffin-embedded tumor samples were classified for TP53 functional status (mutant-like/wild-type-like) using a validated RNA signature. For IHC-based TP53 status, mutant-like was classified as at least 10% positivity. We used two-stage polytomous logistic regression to evaluate risk factor heterogeneity due to RNA-based TP53 and/or ER, adjusting for each other and for PR, HER2, and grade. We then compared this with the results when using IHC-based TP53 classification.

RESULTS

The RNA-based classifier identified 55% of tumors as TP53 wild-type-like and 45% as mutant-like. Several hormone-related factors (oral contraceptive use, menopausal status, age at menopause, and pre- and postmenopausal body mass index) were associated with TP53 mutant-like status, whereas reproductive factors (age at first birth and parity) and smoking were associated with ER status. Multiparity was associated with both TP53 and ER. When classifying TP53 status using IHC methods, no associations were observed with TP53. Associations observed with RNA-based TP53 remained after accounting for basal-like subtype.

CONCLUSIONS

This case-case study found breast cancer risk factors associated with RNA-based TP53 and ER.

IMPACT

RNA-based TP53 and ER represent an emerging etiologic schema of interest in breast cancer prevention research.

Genome-wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analyses

Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype1-3. To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 × 10-8), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.

Diesel exhaust and bladder cancer risk by pathologic stage and grade subtypes

BACKGROUND

The International Agency for Research on Cancer (IARC) classifies diesel engine exhaust as carcinogenic to humans based on sufficient evidence for lung cancer. IARC noted, however, an increased risk of bladder cancer (based on limited evidence).

OBJECTIVE

To evaluate the association between quantitative, lifetime occupational diesel exhaust exposure and risk of urothelial cell carcinoma of the bladder (UBC) overall and according to pathological subtypes.

METHODS

Data from personal interviews with 1944 UBC cases, as well as formalin-fixed paraffin-embedded tumor tissue blocks, and 2135 controls were pooled from two case-control studies conducted in the U.S. and Spain. Lifetime occupational histories combined with exposure-oriented questions were used to estimate cumulative exposure to respirable elemental carbon (REC), a primary surrogate for diesel exhaust. Unconditional logistic regression and two-stage polytomous logistic regression were used to calculate odds ratios (ORs) and 95% confidence intervals (CIs), adjusting for smoking and other risk factors.

RESULTS

Exposure to cumulative REC was associated with an increased risk of UBC; workers with cumulative REC >396 μg/m³-years had an OR of 1.61 (95% CI, 1.08-2.40). At this level of cumulative exposure, similar results were observed in the U.S. and Spain, OR = 1.75 (95% CI, 0.97-3.15) and OR = 1.54 (95% CI, 0.89-2.68), respectively. In lagged analysis, we also observed a consistent increased risk among workers with cumulative REC >396 μg/m³-years (range of ORs = 1.52-1.93) for all lag intervals evaluated (5-40 years). When we accounted for tumor subtypes defined by stage and grade, a significant association between diesel exhaust exposure and UBC was apparent (global test for association p = 0.0019).

CONCLUSIONS

Combining data from two large epidemiologic studies, our results provide further evidence that diesel exhaust exposure increases the risk of UBC.