Spurdle, Amanda B.

Abstract Background Breast cancer is etiologically heterogeneous, but which risk factors differ in their associations across tumor subtypes remains unclear. We conducted a large, pooled analysis to evaluate independent, dose-response associations between breast cancer risk factors and quantitative tumor features. Methods Analyses of 15,731 invasive breast cancers from 24 studies evaluated associations (p-trend) between reproductive and hormonal factors, body mass index (BMI), alcohol, smoking, and family history in relation to quantitative immunohistochemistry measures on tissue microarrays (ER, PR, HER2, KI67, TP53) and tumor grade. Analyses in a subset of 10 population-based studies estimated subtype-specific odds ratios (ORs) comparing cases to controls. A Bayesian False Discovery Probability (BFDP) <0.2 was used to identify associations with strong statistical evidence. Results Nulliparity and later age at menopause were associated with higher ER-positivity (p-trend=0.021 and 0.001, respectively), with corresponding OR[ER+] (95% CI) = 1.49 (1.16-1.90) for nulliparous vs. parous and 1.07 (1.03-1.11) per 5 years. Current combined menopausal hormone therapy (MHT) use was associated with lower grade (p-trend<0.001), with OR [grade1] = 3.37 (2.69-4.21) for current vs. never users. Higher BMI was associated with lower ER-positivity and higher grade in premenopausal women (p-trend<0.001 and <0.001), with OR[ER+] = 0.80 (0.74-0.87) and OR[grade1] = 0.75 (0.63-0.88) per 5 units, and with higher PR-positivity and higher grade in postmenopausal women (p-trend<0.001 and <0.001), with OR[PR+] = 1.08 (1.03-1.14) and OR[grade 3] = 1.10 (1.04-1.17) per 5 units. Conclusion This pooled analysis of 15,731 cases showed that nulliparity, age at menopause, MHT, and BMI have independent, dose-response associations with ER, PR, and grade, clarifying patterns of etiologic heterogeneity. Associations with HER2, KI67 and TP53, or other risk factors did not meet our threshold for strong evidence.

Abstract The success of precision cancer prevention and treatment hinges on accurate discrimination between benign and pathogenic germline alleles in cancer susceptibility genes. Variants of uncertain clinical significance (VUS) present a challenge for cancer risk assessment and functional data will be essential to resolve many VUS. To address this, we developed a cloud-based environment to collate, curate, integrate and analyze all published functional data related to BRCA1 and BRCA2 missense VUS. Then, functional evidence criteria for pathogenicity were determined for each assay using the American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) classification guidelines and ClinGen criteria. Our approach involved retrieving published articles reporting one or more functional assays of BRCA1/2 missense variants reporting variant impact on various biochemical and cell biological functions. Functional results were recorded using original authors' thresholds and categories, then converted to ordinal variables to harmonize: [0 = no functional impact], [1 = intermediate impact], and [2 = functional impact]. For BRCA1, we updated our previous published integration (Genet Med. 2021;23:306-315) with 14 additional publications, collating results from 53 instances of functional assays reporting functional data on 3,219 unique missense variants. For BRCA2, 28 publications were identified reporting data from 141 instances of functional assays on 4,759 unique missense variants. Utilizing a panel of 542 and 396 known pathogenic and benign reference variants for BRCA1 and BRCA2, respectively, we determined the sensitivity, specificity, and ACMG/AMP odds of pathogenicity for each instance of the assay. The remaining variants were assigned a relevant functional assay code weight for each assay instance. Our study successfully derived ACMG/AMP evidence strength from functional data for 3,040 BRCA1 (94% of 3,219) and 2,704 BRCA2 (57% of 4,759) missense variants. This comprehensive work underscores the utility of functional data in informing classification of BRCA1 and BRCA2 missense VUS. Citation Format: Paulo Lyra, Lucas Dalcolmo, Michael Parsons, Thales Nepomuceno, Samuel Brito, Nam Phuong N. Nguyen, Geise de Oliveira, Joao Paulo da Siva, Laura Caleca, Tanisha Taneja, Chunling Hu, Marcy Richardson, Maria Rossing, Aditi Hazra, Alexandra Martins, Sandrine Caputo, Gael A. Millot, Joanne NgeowYuen Yie, Marcelo A. Carvalho, Melissa Cline, Paolo Radice, Rachael Carlsen, Romy Mesman, Valentina Zampiga, Rehan Villani, Vijay Josef, Shyam Sharan, Kyriaki Michailidou, Amanda B. Spurdle, Fergus Couch, Maaike P.G. Vreeswijk, Alvaro N. Monteiro. Integration of functional data to classify BRCA1/2 missense variants: An ENIGMA project [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7325.

A large number of variants identified through clinical genetic testing in disease susceptibility genes are of uncertain significance (VUS). Following the recommendations of the American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP), the frequency in case-control datasets (PS4 criterion) can inform their interpretation. We present a novel case-control likelihood ratio-based method that incorporates gene-specific age-related penetrance. We demonstrate the utility of this method in the analysis of simulated and real datasets. In the analysis of simulated data, the likelihood ratio method was more powerful compared to other methods. Likelihood ratios were calculated for a case-control dataset of BRCA1 and BRCA2 variants from the Breast Cancer Association Consortium (BCAC) and compared with logistic regression results. A larger number of variants reached evidence in favor of pathogenicity, and a substantial number of variants had evidence against pathogenicity—findings that would not have been reached using other case-control analysis methods. Our novel method provides greater power to classify rare variants compared with classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and preformatted Excel calculators for implementation of the method for rare variants in BRCA1, BRCA2, and other high-risk genes with known penetrance.

Abstract Purpose: Germline pathogenic variants in CHEK2 confer moderately elevated breast cancer risk (odds ratio, OR ∼ 2.5), qualifying carriers for enhanced breast cancer screening. Besides pathogenic variants, dozens of missense CHEK2 variants of uncertain significance (VUS) have been identified, hampering the clinical utility of germline genetic testing (GGT). Experimental Design: We collected 460 CHEK2 missense VUS identified by the ENIGMA consortium in 15 countries. Their functional characterization was performed using CHEK2-complementation assays quantifying KAP1 phosphorylation and CHK2 autophosphorylation in human RPE1–CHEK2-knockout cells. Concordant results in both functional assays were used to categorize CHEK2 VUS from 12 ENIGMA case–control datasets, including 73,048 female patients with breast cancer and 88,658 ethnicity-matched controls. Results: A total of 430/460 VUS were successfully analyzed, of which 340 (79.1%) were concordant in both functional assays and categorized as functionally impaired (N = 102), functionally intermediate (N = 12), or functionally wild-type (WT)–like (N = 226). We then examined their association with breast cancer risk in the case–control analysis. The OR and 95% CI (confidence intervals) for carriers of functionally impaired, intermediate, and WT-like variants were 2.83 (95% CI, 2.35–3.41), 1.57 (95% CI, 1.41–1.75), and 1.19 (95% CI, 1.08–1.31), respectively. The meta-analysis of population-specific datasets showed similar results. Conclusions: We determined the functional consequences for the majority of CHEK2 missense VUS found in patients with breast cancer (3,660/4,436; 82.5%). Carriers of functionally impaired missense variants accounted for 0.5% of patients with breast cancer and were associated with a moderate risk similar to that of truncating CHEK2 variants. In contrast, 2.2% of all patients with breast cancer carried functionally wild-type/intermediate missense variants with no clinically relevant breast cancer risk in heterozygous carriers.

Abstract Background The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci. Methods We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci. Results We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci. Conclusions Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Abstract Background Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors. Methods We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set. Results We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls. Conclusions In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Abstract Background: Accumulating evidence suggests a relationship between endometrial cancer and ovarian cancer. Independent genome-wide association studies (GWAS) for endometrial cancer and ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. We aimed to identify joint endometrial and ovarian cancer risk loci by performing a meta-analysis of GWAS summary statistics from these two cancers. Methods: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data. Results: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). We found seven loci associated with risk for both cancers (PBonferroni < 2.4 × 10−9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P < 5 × 10−7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation. Conclusions: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis. Impact: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.

Abstract Blood lipids have been associated with the development of a range of cancers, including breast, lung and colorectal cancer. For endometrial cancer, observational studies have reported inconsistent associations between blood lipids and cancer risk. To reduce biases from unmeasured confounding, we performed a bidirectional, two‐sample Mendelian randomization analysis to investigate the relationship between levels of three blood lipids (low‐density lipoprotein [LDL] and high‐density lipoprotein [HDL] cholesterol, and triglycerides) and endometrial cancer risk. Genetic variants associated with each of these blood lipid levels ( P < 5 × 10 −8 ) were identified as instrumental variables, and assessed using genome‐wide association study data from the Endometrial Cancer Association Consortium (12 906 cases and 108 979 controls) and the Global Lipids Genetic Consortium (n = 188 578). Mendelian randomization analyses found genetically raised LDL cholesterol levels to be associated with lower risks of endometrial cancer of all histologies combined, and of endometrioid and non‐endometrioid subtypes. Conversely, higher genetically predicted HDL cholesterol levels were associated with increased risk of non‐endometrioid endometrial cancer. After accounting for the potential confounding role of obesity (as measured by genetic variants associated with body mass index), the association between genetically predicted increased LDL cholesterol levels and lower endometrial cancer risk remained significant, especially for non‐endometrioid endometrial cancer. There was no evidence to support a role for triglycerides in endometrial cancer development. Our study supports a role for LDL and HDL cholesterol in the development of non‐endometrioid endometrial cancer. Further studies are required to understand the mechanisms underlying these findings.

Abstract Background: A polygenic hazard score (PHS), the weighted sum of 54 SNP genotypes, was previously validated for association with clinically significant prostate cancer and for improved prostate cancer screening accuracy. Here, we assess the potential impact of PHS-informed screening. Methods: United Kingdom population incidence data (Cancer Research United Kingdom) and data from the Cluster Randomized Trial of PSA Testing for Prostate Cancer were combined to estimate age-specific clinically significant prostate cancer incidence (Gleason score ≥7, stage T3–T4, PSA ≥10, or nodal/distant metastases). Using HRs estimated from the ProtecT prostate cancer trial, age-specific incidence rates were calculated for various PHS risk percentiles. Risk-equivalent age, when someone with a given PHS percentile has prostate cancer risk equivalent to an average 50-year-old man (50-year-standard risk), was derived from PHS and incidence data. Positive predictive value (PPV) of PSA testing for clinically significant prostate cancer was calculated using PHS-adjusted age groups. Results: The expected age at diagnosis of clinically significant prostate cancer differs by 19 years between the 1st and 99th PHS percentiles: men with PHS in the 1st and 99th percentiles reach the 50-year-standard risk level at ages 60 and 41, respectively. PPV of PSA was higher for men with higher PHS-adjusted age. Conclusions: PHS provides individualized estimates of risk-equivalent age for clinically significant prostate cancer. Screening initiation could be adjusted by a man's PHS. Impact: Personalized genetic risk assessments could inform prostate cancer screening decisions.

AbstractPrevious transcriptome‐wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome‐wide association studies (GWAS), but analyses of breast cancer subtype‐specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta‐analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER−). We further compared associations with ER+ and ER− subtypes, using a case‐only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER– breast cancer. We further identified 30 TWAS‐significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast‐cancer gene in three of six regions harboring multiple TWAS‐significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER− breast cancer.