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.

Body Size, Diet Quality, and Epigenetic Aging: Cross-Sectional and Longitudinal Analyses

Epigenetic age is an emerging marker of health that is highly predictive of disease and mortality risk. There is a lack of evidence on whether lifestyle changes are associated with changes in epigenetic aging. We used data from 1 041 participants in the Melbourne Collaborative Cohort Study with blood DNA methylation measures at baseline (1990-1994, mean age: 57.4 years) and follow-up (2003-2007, mean age: 68.8 years). The Alternative Healthy Eating Index-2010 (AHEI-2010), the Mediterranean Dietary Score, and the Dietary Inflammatory Index were used as measures of diet quality, and weight, waist circumference, and waist-to-hip ratio as measures of body size. Five age-adjusted epigenetic aging measures were considered: GrimAge, PhenoAge, PCGrimAge, PCPhenoAge, and DunedinPACE. Multivariable linear regression models including restricted cubic splines were used to assess the cross-sectional and longitudinal associations of body size and diet quality with epigenetic aging. Associations between weight and epigenetic aging cross-sectionally at both time points were positive and appeared greater for DunedinPACE (per SD: β ~0.24) than for GrimAge and PhenoAge (β ~0.10). The longitudinal associations with weight change were markedly nonlinear (U-shaped) with stable weight being associated with the lowest epigenetic aging at follow-up, except for DunedinPACE, for which only weight gain showed a positive association. We found negative, linear associations for AHEI-2010 both cross-sectionally and longitudinally. Other adiposity measures and dietary scores showed similar results. In middle-aged to older adults, declining diet quality and weight gain may increase epigenetic age, while the association for weight loss may require further investigation. Our study sheds light on the potential of weight management and dietary improvement in slowing aging processes.

Breast and bowel cancers diagnosed in people 'too young to have cancer': A blueprint for research using family and twin studies

Young breast and bowel cancers (e.g., those diagnosed before age 40 or 50 years) have far greater morbidity and mortality in terms of years of life lost, and are increasing in incidence, but have been less studied. For breast and bowel cancers, the familial relative risks, and therefore the familial variances in age-specific log(incidence), are much greater at younger ages, but little of these familial variances has been explained. Studies of families and twins can address questions not easily answered by studies of unrelated individuals alone. We describe existing and emerging family and twin data that can provide special opportunities for discovery. We present designs and statistical analyses, including novel ideas such as the VALID (Variance in Age-specific Log Incidence Decomposition) model for causes of variation in risk, the DEPTH (DEPendency of association on the number of Top Hits) and other approaches to analyse genome-wide association study data, and the within-pair, ICE FALCON (Inference about Causation from Examining FAmiliaL CONfounding) and ICE CRISTAL (Inference about Causation from Examining Changes in Regression coefficients and Innovative STatistical AnaLysis) approaches to causation and familial confounding. Example applications to breast and colorectal cancer are presented. Motivated by the availability of the resources of the Breast and Colon Cancer Family Registries, we also present some ideas for future studies that could be applied to, and compared with, cancers diagnosed at older ages and address the challenges posed by young breast and bowel cancers.

Causation and familial confounding as explanations for the associations of polygenic risk scores with breast cancer: Evidence from innovative ICE FALCON and ICE CRISTAL analyses

A polygenic risk score (PRS) combines the associations of multiple genetic variants that could be due to direct causal effects, indirect genetic effects, or other sources of familial confounding. We have developed new approaches to assess evidence for and against causation by using family data for pairs of relatives (Inference about Causation from Examination of FAmiliaL CONfounding [ICE FALCON]) or measures of family history (Inference about Causation from Examining Changes in Regression coefficients and Innovative STatistical AnaLyses [ICE CRISTAL]). Inference is made from the changes in regression coefficients of relatives' PRSs or PRS and family history before and after adjusting for each other. We applied these approaches to two breast cancer PRSs and multiple studies and found that (a) for breast cancer diagnosed at a young age, for example, <50 years, there was no evidence that the PRSs were causal, while (b) for breast cancer diagnosed at later ages, there was consistent evidence for causation explaining increasing amounts of the PRS-disease association. The genetic variants in the PRS might be in linkage disequilibrium with truly causal variants and not causal themselves. These PRSs cause minimal heritability of breast cancer at younger ages. There is also evidence for nongenetic factors shared by first-degree relatives that explain breast cancer familial aggregation. Familial associations are not necessarily due to genes, and genetic associations are not necessarily causal.

Large-scale genome-wide association study of 398,238 women unveils seven novel loci associated with high-grade serous epithelial ovarian cancer risk

Background

Nineteen genomic regions have been associated with high-grade serous ovarian cancer (HGSOC). We used data from the Ovarian Cancer Association Consortium (OCAC), Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA), UK Biobank (UKBB), and FinnGen to identify novel HGSOC susceptibility loci and develop polygenic scores (PGS).

Methods

We analyzed >22 million variants for 398,238 women. Associations were assessed separately by consortium and meta-analysed. OCAC and CIMBA data were used to develop PGS which were trained on FinnGen data and validated in UKBB and BioBank Japan.

Results

Eight novel variants were associated with HGSOC risk. An interesting discovery biologically was finding that TP53 3'-UTR SNP rs78378222 was associated with HGSOC (per T allele relative risk (RR)=1.44, 95%CI:1.28-1.62, P=1.76×10-9). The optimal PGS included 64,518 variants and was associated with an odds ratio of 1.46 (95%CI:1.37-1.54) per standard deviation in the UKBB validation (AUROC curve=0.61, 95%CI:0.59-0.62).

Conclusions

This study represents the largest GWAS for HGSOC to date. The results highlight that improvements in imputation reference panels and increased sample sizes can identify HGSOC associated variants that previously went undetected, resulting in improved PGS. The use of updated PGS in cancer risk prediction algorithms will then improve personalized risk prediction for HGSOC.

Intratumoral presence of the genotoxic gut bacteria pks+ E. coli, Enterotoxigenic Bacteroides fragilis, and Fusobacterium nucleatum and their association with clinicopathological and molecular features of colorectal cancer

BACKGROUND

This study aimed to investigate clinicopathological and molecular tumour features associated with intratumoral pks+ Escherichia coli (pks+E.coli+), pks+E.coli- (non-E.coli bacteria harbouring the pks island), Enterotoxigenic Bacteroides fragilis (ETBF) and Fusobacterium nucleatum (F. nucleatum).

METHODS

We screened 1697 tumour-derived DNA samples from the Australasian Colorectal Cancer Family Registry, Melbourne Collaborative Cohort Study and the ANGELS study using targeted PCR.

RESULTS

Pks+E.coli+ was associated with male sex (P < 0.01) and APC:c.835-8 A > G somatic mutation (P = 0.03). The association between pks+E.coli+ and APC:c.835-8 A > G was specific to early-onset CRCs (diagnosed<45years, P = 0.02). The APC:c.835-A > G was not associated with pks+E.coli- (P = 0.36). F. nucleatum was associated with DNA mismatch repair deficiency (MMRd), BRAF:c.1799T>A p.V600E mutation, CpG island methylator phenotype, proximal tumour location, and high levels of tumour infiltrating lymphocytes (Ps < 0.01). In the stratified analysis by MMRd subgroups, F. nucleatum was associated with Lynch syndrome, MLH1 methylated and double MMR somatic mutated MMRd subgroups (Ps < 0.01).

CONCLUSION

Intratumoral pks+E.coli+ but not pks+E.coli- are associated with CRCs harbouring the APC:c.835-8 A > G somatic mutation, suggesting that this mutation is specifically related to DNA damage from colibactin-producing E.coli exposures. F. nucleatum was associated with both hereditary and sporadic MMRd subtypes, suggesting the MMRd tumour microenvironment is important for F. nucleatum colonisation irrespective of its cause.

Differences in polygenic score distributions in European ancestry populations: implications for breast cancer risk prediction

The 313-variant polygenic risk score (PRS313) provides a promising tool for breast cancer risk prediction. However, evaluation of the PRS313 across different European populations which could influence risk estimation has not been performed. Here, we explored the distribution of PRS313 across European populations using genotype data from 94,072 females without breast cancer, of European-ancestry from 21 countries participating in the Breast Cancer Association Consortium (BCAC) and 225,105 female participants from the UK Biobank. The mean PRS313 differed markedly across European countries, being highest in south-eastern Europe and lowest in north-western Europe. Using the overall European PRS313 distribution to categorise individuals leads to overestimation and underestimation of risk in some individuals from south-eastern and north-western countries, respectively. Adjustment for principal components explained most of the observed heterogeneity in mean PRS. Country-specific PRS distributions may be used to calibrate risk categories in individuals from different countries.

Genetic and Environmental Causes of Variation in an Automated Breast Cancer Risk Factor Based on Mammographic Textures

BACKGROUND

Cirrus is an automated risk predictor for breast cancer that comprises texture-based mammographic features and is mostly independent of mammographic density. We investigated genetic and environmental variance of variation in Cirrus.

METHODS

We measured Cirrus for 3,195 breast cancer-free participants, including 527 pairs of monozygotic (MZ) twins, 271 pairs of dizygotic (DZ) twins, and 1,599 siblings of twins. Multivariate normal models were used to estimate the variance and familial correlations of age-adjusted Cirrus as a function of age. The classic twin model was expanded to allow the shared environment effects to differ by zygosity. The SNP-based heritability was estimated for a subset of 2,356 participants.

RESULTS

There was no evidence that the variance or familial correlations depended on age. The familial correlations were 0.52 (SE, 0.03) for MZ pairs and 0.16(SE, 0.03) for DZ and non-twin sister pairs combined. Shared environmental factors specific to MZ pairs accounted for 20% of the variance. Additive genetic factors accounted for 32% (SE = 5%) of the variance, consistent with the SNP-based heritability of 36% (SE = 16%).

CONCLUSION

Cirrus is substantially familial due to genetic factors and an influence of shared environmental factors that was evident for MZ twin pairs only. The latter could be due to nongenetic factors operating in utero or in early life that are shared by MZ twins.

IMPACT

Early-life factors, shared more by MZ pairs than DZ/non-twin sister pairs, could play a role in the variation in Cirrus, consistent with early life being recognized as a critical window of vulnerability to breast carcinogens.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Germline Sequencing Analysis to Inform Clinical Gene Panel Testing for Aggressive Prostate Cancer

Importance

Germline gene panel testing is recommended for men with advanced prostate cancer (PCa) or a family history of cancer. While evidence is limited for some genes currently included in panel testing, gene panels are also likely to be incomplete and missing genes that influence PCa risk and aggressive disease.

Objective

To identify genes associated with aggressive PCa.

Design, Setting, and Participants

A 2-stage exome sequencing case-only genetic association study was conducted including men of European ancestry from 18 international studies. Data analysis was performed from January 2021 to March 2023. Participants were 9185 men with aggressive PCa (including 6033 who died of PCa and 2397 with confirmed metastasis) and 8361 men with nonaggressive PCa.

Exposure

Sequencing data were evaluated exome-wide and in a focused investigation of 29 DNA repair pathway and cancer susceptibility genes, many of which are included on gene panels.

Main Outcomes and Measures

The primary study outcomes were aggressive (category T4 or both T3 and Gleason score ≥8 tumors, metastatic PCa, or PCa death) vs nonaggressive PCa (category T1 or T2 and Gleason score ≤6 tumors without known recurrence), and metastatic vs nonaggressive PCa.

Results

A total of 17 546 men of European ancestry were included in the analyses; mean (SD) age at diagnosis was 65.1 (9.2) years in patients with aggressive PCa and 63.7 (8.0) years in those with nonaggressive disease. The strongest evidence of association with aggressive or metastatic PCa was noted for rare deleterious variants in known PCa risk genes BRCA2 and ATM (P ≤ 1.9 × 10-6), followed by NBN (P = 1.7 × 10-4). This study found nominal evidence (P < .05) of association with rare deleterious variants in MSH2, XRCC2, and MRE11A. Five other genes had evidence of greater risk (OR≥2) but carrier frequency differences between aggressive and nonaggressive PCa were not statistically significant: TP53, RAD51D, BARD1, GEN1, and SLX4. Deleterious variants in these 11 candidate genes were carried by 2.3% of patients with nonaggressive, 5.6% with aggressive, and 7.0% with metastatic PCa.

Conclusions and Relevance

The findings of this study provide further support for DNA repair and cancer susceptibility genes to better inform disease management in men with PCa and for extending testing to men with nonaggressive disease, as men carrying deleterious alleles in these genes are likely to develop more advanced disease.

Causal relationships between breast cancer risk factors based on mammographic features

BACKGROUND

Mammogram risk scores based on texture and density defined by different brightness thresholds are associated with breast cancer risk differently and could reveal distinct information about breast cancer risk. We aimed to investigate causal relationships between these intercorrelated mammogram risk scores to determine their relevance to breast cancer aetiology.

METHODS

We used digitised mammograms for 371 monozygotic twin pairs, aged 40-70 years without a prior diagnosis of breast cancer at the time of mammography, from the Australian Mammographic Density Twins and Sisters Study. We generated normalised, age-adjusted, and standardised risk scores based on textures using the Cirrus algorithm and on three spatially independent dense areas defined by increasing brightness threshold: light areas, bright areas, and brightest areas. Causal inference was made using the Inference about Causation from Examination of FAmilial CONfounding (ICE FALCON) method.

RESULTS

The mammogram risk scores were correlated within twin pairs and with each other (r = 0.22-0.81; all P < 0.005). We estimated that 28-92% of the associations between the risk scores could be attributed to causal relationships between the scores, with the rest attributed to familial confounders shared by the scores. There was consistent evidence for positive causal effects: of Cirrus, light areas, and bright areas on the brightest areas (accounting for 34%, 55%, and 85% of the associations, respectively); and of light areas and bright areas on Cirrus (accounting for 37% and 28%, respectively).

CONCLUSIONS

In a mammogram, the lighter (less dense) areas have a causal effect on the brightest (highly dense) areas, including through a causal pathway via textural features. These causal relationships help us gain insight into the relative aetiological importance of different mammographic features in breast cancer. For example our findings are consistent with the brightest areas being more aetiologically important than lighter areas for screen-detected breast cancer; conversely, light areas being more aetiologically important for interval breast cancer. Additionally, specific textural features capture aetiologically independent breast cancer risk information from dense areas. These findings highlight the utility of ICE FALCON and family data in decomposing the associations between intercorrelated disease biomarkers into distinct biological pathways.

"Out of the blue": A qualitative study exploring the experiences of women and next of kin receiving unexpected results from BRA-STRAP research gene panel testing

In the genomic era, the availability of gene panel and whole genome/exome sequencing is rapidly increasing. Opportunities for providing former patients with new genetic information are also increasing over time and recontacting former patients with new information is likely to become more common. Breast cancer Refined Analysis of Sequence Tests-Risk And Penetrance (BRA-STRAP) is an Australian study of individuals who had previously undertaken BRCA1 and BRCA2 genetic testing, with no pathogenic variants detected. Using a waiver of consent, stored DNA samples were retested using a breast/ovarian cancer gene panel and clinically significant results returned to the patient (or next of kin, if deceased). This qualitative study aimed to explore patient experiences, opinions, and expectations of recontacting in the Australian hereditary cancer setting. Participants were familial cancer clinic patients (or next of kin) who were notified of a new pathogenic variant identified via BRA-STRAP. In-depth, semi-structured interviews were conducted approximately 6 weeks post-result. Interviews were transcribed verbatim and analyzed using an inductive thematic approach. Thirty participants (all female; average age = 57; range 36-84) were interviewed. Twenty-five were probands, and five were next of kin. Most women reported initial shock upon being recontacted with unexpected news, after having obtained a sense of closure related to their initial genetic testing experiences and cancer diagnosis. For most, this initial distress was short-lived, followed by a process of readjustment, meaning-making and adaptation that was facilitated by perceived clinical and personal utility of the information. Women were overall satisfied with the waiver of consent approach and recontacting process. Results are in line with previous studies suggesting that patients have positive attitudes about recontacting. Women in this study valued new genetic information gained from retesting and were satisfied with the BRA-STRAP recontact model. Practice implications to facilitate readjustment and promote psychosocial adaptation were identified.

Variance of age-specific log incidence decomposition (VALID): a unifying model of measured and unmeasured genetic and non-genetic risks

BACKGROUND

The extent to which known and unknown factors explain how much people of the same age differ in disease risk is fundamental to epidemiology. Risk factors can be correlated in relatives, so familial aspects of risk (genetic and non-genetic) must be considered.

DEVELOPMENT

We present a unifying model (VALID) for variance in risk, with risk defined as log(incidence) or logit(cumulative incidence). Consider a normally distributed risk score with incidence increasing exponentially as the risk increases. VALID's building block is variance in risk, Δ2, where Δ = log(OPERA) is the difference in mean between cases and controls and OPERA is the odds ratio per standard deviation. A risk score correlated r between a pair of relatives generates a familial odds ratio of exp(rΔ2). Familial risk ratios, therefore, can be converted into variance components of risk, extending Fisher's classic decomposition of familial variation to binary traits. Under VALID, there is a natural upper limit to variance in risk caused by genetic factors, determined by the familial odds ratio for genetically identical twin pairs, but not to variation caused by non-genetic factors.

APPLICATION

For female breast cancer, VALID quantified how much variance in risk is explained-at different ages-by known and unknown major genes and polygenes, non-genomic risk factors correlated in relatives, and known individual-specific factors.

CONCLUSION

VALID has shown that, while substantial genetic risk factors have been discovered, much is unknown about genetic and familial aspects of breast cancer risk especially for young women, and little is known about individual-specific variance in risk.

Polymorphisms in genes of melatonin biosynthesis and signaling support the light-at-night hypothesis for breast cancer

Light-at-night triggers the decline of pineal gland melatonin biosynthesis and secretion and is an IARC-classified probable breast-cancer risk factor. We applied a large-scale molecular epidemiology approach to shed light on the putative role of melatonin in breast cancer. We investigated associations between breast-cancer risk and polymorphisms at genes of melatonin biosynthesis/signaling using a study population of 44,405 women from the Breast Cancer Association Consortium (22,992 cases, 21,413 population-based controls). Genotype data of 97 candidate single nucleotide polymorphisms (SNPs) at 18 defined gene regions were investigated for breast-cancer risk effects. We calculated adjusted odds ratios (ORs) and 95% confidence intervals (CI) by logistic regression for the main-effect analysis as well as stratified analyses by estrogen- and progesterone-receptor (ER, PR) status. SNP-SNP interactions were analyzed via a two-step procedure based on logic regression. The Bayesian false-discovery probability (BFDP) was used for all analyses to account for multiple testing. Noteworthy associations (BFDP < 0.8) included 10 linked SNPs in tryptophan hydroxylase 2 (TPH2) (e.g. rs1386492: OR = 1.07, 95% CI 1.02-1.12), and a SNP in the mitogen-activated protein kinase 8 (MAPK8) (rs10857561: OR = 1.11, 95% CI 1.04-1.18). The SNP-SNP interaction analysis revealed noteworthy interaction terms with TPH2- and MAPK-related SNPs (e.g. rs1386483R ∧ rs1473473D ∧ rs3729931D: OR = 1.20, 95% CI 1.09-1.32). In line with the light-at-night hypothesis that links shift work with elevated breast-cancer risks our results point to SNPs in TPH2 and MAPK-genes that may impact the intricate network of circadian regulation.

A likelihood ratio approach for utilizing case-control data in the clinical classification of rare sequence variants: application to BRCA1 and BRCA2

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 analyses 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 to classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and pre-formatted excel calculators for implementation of the method for rare variants in BRCA1, BRCA2 and other high-risk genes with known penetrance.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry

BACKGROUND

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

METHODS

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

RESULTS

Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94).

CONCLUSIONS

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

BACKGROUND

Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers.

AIM

To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS.

METHODS

Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death.

RESULTS

There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09-1.56)].

CONCLUSION

Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.

Methylation scores for smoking, alcohol consumption and body mass index and risk of seven types of cancer

Methylation marks of exposure to health risk factors may be useful markers of cancer risk as they might better capture current and past exposures than questionnaires, and reflect different individual responses to exposure. We used data from seven case-control studies nested within the Melbourne Collaborative Cohort Study of blood DNA methylation and risk of colorectal, gastric, kidney, lung, prostate and urothelial cancer, and B-cell lymphoma (N cases = 3123). Methylation scores (MS) for smoking, body mass index (BMI), and alcohol consumption were calculated based on published data as weighted averages of methylation values. Rate ratios (RR) and 95% confidence intervals for association with cancer risk were estimated using conditional logistic regression and expressed per SD increase of the MS, with and without adjustment for health-related confounders. The contribution of MS to discriminate cases from controls was evaluated using the area under the curve (AUC). After confounder adjustment, we observed: large associations (RR = 1.5-1.7) with lung cancer risk for smoking MS; moderate associations (RR = 1.2-1.3) with urothelial cancer risk for smoking MS and with mature B-cell neoplasm risk for BMI and alcohol MS; moderate to small associations (RR = 1.1-1.2) for BMI and alcohol MS with several cancer types and cancer overall. Generally small AUC increases were observed after inclusion of several MS in the same model (colorectal, gastric, kidney, urothelial cancers: +3%; lung cancer: +7%; B-cell neoplasms: +8%). Methylation scores for smoking, BMI and alcohol consumption show independent associations with cancer risk, and may provide some improvements in risk prediction.

Heritable methylation marks associated with prostate cancer risk

DNA methylation marks that are inherited from parents to offspring are known to play a role in cancer risk and could explain part of the familial risk for cancer. We therefore conducted a genome-wide search for heritable methylation marks associated with prostate cancer risk. Peripheral blood DNA methylation was measured for 133 of the 469 members of 25 multiple-case prostate cancer families, using the EPIC array. We used these families to systematically search the genome for methylation marks with Mendelian patterns of inheritance, then we tested the 1,000 most heritable marks for association with prostate cancer risk. After correcting for multiple testing, 41 heritable methylation marks were associated with prostate cancer risk. Separate analyses, based on 869 incident cases and 869 controls from a prospective cohort study, showed that 9 of these marks near the metastable epiallele VTRNA2-1 were also nominally associated with aggressive prostate cancer risk in the population.

Modifiable lifestyle risk factors and survival after diagnosis with multiple myeloma

BACKGROUND

While remaining incurable, median overall survival for MM now exceeds 5 years. Yet few studies have investigated how modifiable lifestyle factors influence survival. We investigate whether adiposity, diet, alcohol, or smoking are associated with MM-related fatality.

RESEARCH DESIGN AND METHODS

We recruited 760 incident cases of MM via cancer registries in two Australian states during 2010-2016. Participants returned questionnaires on health and lifestyle. Follow-up ended in 2020. Flexible parametric survival models were used to estimate hazard ratios (HR) and 95% confidence intervals (CI) for lifestyle exposures and risk of all-cause and MM-specific fatality.

RESULTS

Higher pre-diagnosis Alternative Healthy Eating Index (AHEI) scores were associated with reduced MM-specific fatality (per 10-unit score, HR = 0.84, 95%CI = 0.70-0.99). Pre-diagnosis alcohol consumption was inversely associated with MM-specific fatality, compared with nondrinkers (0.1-20 g per day, HR = 0.59, 95%CI = 0.39-0.90; >20 g per day, HR = 0.67, 95%CI = 0.40-1.13). Tobacco smoking was associated with increased all-cause fatality compared with never smoking (former smokers: HR = 1.44, 95%CI = 1.10-1.88; current smokers: HR = 1.30, 95%CI = 0.80-2.10). There was no association between pre-enrollment body mass index (BMI) and MM-specific or all-cause fatality.

CONCLUSIONS

Our findings support established recommendations for healthy diets and against smoking. Higher quality diet, as measured by the AHEI, may improve survival post diagnosis with MM.

Genome-wide Association Study of Bladder Cancer Reveals New Biological and Translational Insights

BACKGROUND

Genomic regions identified by genome-wide association studies (GWAS) for bladder cancer risk provide new insights into etiology.

OBJECTIVE

To identify new susceptibility variants for bladder cancer in a meta-analysis of new and existing genome-wide genotype data.

DESIGN, SETTING, AND PARTICIPANTS

Data from 32 studies that includes 13,790 bladder cancer cases and 343,502 controls of European ancestry were used for meta-analysis.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSES

Log-additive associations of genetic variants were assessed using logistic regression models. A fixed-effects model was used for meta-analysis of the results. Stratified analyses were conducted to evaluate effect modification by sex and smoking status. A polygenic risk score (PRS) was generated on the basis of known and novel susceptibility variants and tested for interaction with smoking.

RESULTS AND LIMITATIONS

Multiple novel bladder cancer susceptibility loci (6p.22.3, 7q36.3, 8q21.13, 9p21.3, 10q22.1, 19q13.33) as well as improved signals in three known regions (4p16.3, 5p15.33, 11p15.5) were identified, bringing the number of independent markers at genome-wide significance (p < 5 × 10-8) to 24. The 4p16.3 (FGFR3/TACC3) locus was associated with a stronger risk for women than for men (p-interaction = 0.002). Bladder cancer risk was increased by interactions between smoking status and genetic variants at 8p22 (NAT2; multiplicative p value for interaction [pM-I] = 0.004), 8q21.13 (PAG1; pM-I = 0.01), and 9p21.3 (LOC107987026/MTAP/CDKN2A; pM-I = 0.02). The PRS based on the 24 independent GWAS markers (odds ratio per standard deviation increase 1.49, 95% confidence interval 1.44-1.53), which also showed comparable results in two prospective cohorts (UK Biobank, PLCO trial), revealed an approximately fourfold difference in the lifetime risk of bladder cancer according to the PRS (e.g., 1st vs 10th decile) for both smokers and nonsmokers.

CONCLUSIONS

We report novel loci associated with risk of bladder cancer that provide clues to its biological underpinnings. Using 24 independent markers, we constructed a PRS to stratify lifetime risk. The PRS combined with smoking history, and other established risk factors, has the potential to inform future screening efforts for bladder cancer.

PATIENT SUMMARY

We identified new genetic markers that provide biological insights into the genetic causes of bladder cancer. These genetic risk factors combined with lifestyle risk factors, such as smoking, may inform future preventive and screening strategies for bladder cancer.

Spectrum and Frequency of Germline FANCM Protein-Truncating Variants in 44,803 European Female Breast Cancer Cases

FANCM germline protein truncating variants (PTVs) are moderate-risk factors for ER-negative breast cancer. We previously described the spectrum of FANCM PTVs in 114 European breast cancer cases. In the present, larger cohort, we report the spectrum and frequency of four common and 62 rare FANCM PTVs found in 274 carriers detected among 44,803 breast cancer cases. We confirmed that p.Gln1701* was the most common PTV in Northern Europe with lower frequencies in Southern Europe. In contrast, p.Gly1906Alafs*12 was the most common PTV in Southern Europe with decreasing frequencies in Central and Northern Europe. We verified that p.Arg658* was prevalent in Central Europe and had highest frequencies in Eastern Europe. We also confirmed that the fourth most common PTV, p.Gln498Thrfs*7, might be a founder variant from Lithuania. Based on the frequency distribution of the carriers of rare PTVs, we showed that the FANCM PTVs spectra in Southwestern and Central Europe were much more heterogeneous than those from Northeastern Europe. These findings will inform the development of more efficient FANCM genetic testing strategies for breast cancer cases from specific European populations.

Associations of height, body mass index, and weight gain with breast cancer risk in carriers of a pathogenic variant in BRCA1 or BRCA2: the BRCA1 and BRCA2 Cohort Consortium

INTRODUCTION

Height, body mass index (BMI), and weight gain are associated with breast cancer risk in the general population. It is unclear whether these associations also exist for carriers of pathogenic variants in the BRCA1 or BRCA2 genes.

PATIENTS AND METHODS

An international pooled cohort of 8091 BRCA1/2 variant carriers was used for retrospective and prospective analyses separately for premenopausal and postmenopausal women. Cox regression was used to estimate breast cancer risk associations with height, BMI, and weight change.

RESULTS

In the retrospective analysis, taller height was associated with risk of premenopausal breast cancer for BRCA2 variant carriers (HR 1.20 per 10 cm increase, 95% CI 1.04-1.38). Higher young-adult BMI was associated with lower premenopausal breast cancer risk for both BRCA1 (HR 0.75 per 5 kg/m2, 95% CI 0.66-0.84) and BRCA2 (HR 0.76, 95% CI 0.65-0.89) variant carriers in the retrospective analysis, with consistent, though not statistically significant, findings from the prospective analysis. In the prospective analysis, higher BMI and adult weight gain were associated with higher postmenopausal breast cancer risk for BRCA1 carriers (HR 1.20 per 5 kg/m2, 95% CI 1.02-1.42; and HR 1.10 per 5 kg weight gain, 95% CI 1.01-1.19, respectively).

CONCLUSION

Anthropometric measures are associated with breast cancer risk for BRCA1 and BRCA2 variant carriers, with relative risk estimates that are generally consistent with those for women from the general population.

Ovarian cancer pathology characteristics as predictors of variant pathogenicity in BRCA1 and BRCA2

BACKGROUND

The distribution of ovarian tumour characteristics differs between germline BRCA1 and BRCA2 pathogenic variant carriers and non-carriers. In this study, we assessed the utility of ovarian tumour characteristics as predictors of BRCA1 and BRCA2 variant pathogenicity, for application using the American College of Medical Genetics and the Association for Molecular Pathology (ACMG/AMP) variant classification system.

METHODS

Data for 10,373 ovarian cancer cases, including carriers and non-carriers of BRCA1 or BRCA2 pathogenic variants, were collected from unpublished international cohorts and consortia and published studies. Likelihood ratios (LR) were calculated for the association of ovarian cancer histology and other characteristics, with BRCA1 and BRCA2 variant pathogenicity. Estimates were aligned to ACMG/AMP code strengths (supporting, moderate, strong).

RESULTS

No histological subtype provided informative ACMG/AMP evidence in favour of BRCA1 and BRCA2 variant pathogenicity. Evidence against variant pathogenicity was estimated for the mucinous and clear cell histologies (supporting) and borderline cases (moderate). Refined associations are provided according to tumour grade, invasion and age at diagnosis.

CONCLUSIONS

We provide detailed estimates for predicting BRCA1 and BRCA2 variant pathogenicity based on ovarian tumour characteristics. This evidence can be combined with other variant information under the ACMG/AMP classification system, to improve classification and carrier clinical management.

Adiposity and plasma concentrations of kynurenine pathway metabolites and traditional markers of inflammation

AIM

The kynurenine pathway is increasingly recognised to play a role in inflammation and disease. We assessed the cross-sectional and longitudinal associations of adiposity measures (body mass index, waist-hip ratio, waist circumference and fat mass ratio) with plasma concentrations of kynurenine pathway metabolites and traditional markers of inflammation.

METHODS

We used data from 970 Melbourne Collaborative Cohort Study participants who had plasma markers measured at baseline (median age 59 years) and follow-up (median age 70 years). Linear regression was used to assess cross-sectional and longitudinal associations between four adiposity measures and concentrations of i) nine kynurenine pathway metabolites; ii) two derived markers; iii) eight traditional inflammatory markers.

RESULTS

Cross-sectionally, most kynurenine metabolites were strongly associated with adiposity measures at both time points; associations were generally stronger than for most inflammation markers except CRP (e.g. body mass index at baseline, quinolinic acid (per S.D. β = 0.30, 95%CI: 0.24-0.36, P = 10^-21), kynurenine (β = 0.25, 95%CI: 0.19-0.31, P = 10^-16) and CRP (β = 0.31, 95%CI: 0.25-0.37, P = 10^-24), and remained largely unchanged after adjustment for confounders. Longitudinally, changes in adiposity measures over approximately a decade were positively associated with changes in kynurenine metabolite concentrations (in particular for 3-hydroxyanthranilic acid, kynurenine and quinolinic acid), and more strongly so than for other markers of inflammation, including CRP.

CONCLUSIONS

In middle-aged and older adults, plasma concentrations of kynurenine metabolites are strongly associated with adiposity, both cross-sectionally and longitudinally. Our study demonstrates that kynurenine metabolites may be valuable markers to monitor the adverse consequences of obesity.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. We aimed to assess the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Analyses were based on 82,701 women diagnosed with invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations of treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR(95%CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR(95%CI) :1.30 (1.09-1.56)]. In conclusion, systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk. (Main MS: 3201 words).

Body mass index and molecular subtypes of colorectal cancer

BACKGROUND

Obesity is an established risk factor for colorectal cancer (CRC), but the evidence for the association is inconsistent across molecular subtypes of the disease.

METHODS

We pooled data on body mass index (BMI), tumor microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, and Jass classification types for 11 872 CRC cases and 11 013 controls from 11 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) adjusted for covariables.

RESULTS

Higher BMI was associated with increased CRC risk (OR per 5 kg/m2 = 1.18, 95% CI = 1.15 to 1.22). The positive association was stronger for men than women but similar across tumor subtypes defined by individual molecular markers. In analyses by Jass type, higher BMI was associated with elevated CRC risk for types 1-4 cases but not for type 5 CRC cases (considered familial-like/Lynch syndrome microsatellite instability-H, CpG island methylator phenotype-low or negative, BRAF-wild type, KRAS-wild type, OR = 1.04, 95% CI = 0.90 to 1.20). This pattern of associations for BMI and Jass types was consistent by sex and design of contributing studies (cohort or case-control).

CONCLUSIONS

In contrast to previous reports with fewer study participants, we found limited evidence of heterogeneity for the association between BMI and CRC risk according to molecular subtype, suggesting that obesity influences nearly all major pathways involved in colorectal carcinogenesis. The null association observed for the Jass type 5 suggests that BMI is not a risk factor for the development of CRC for individuals with Lynch syndrome.

FANCM missense variants and breast cancer risk: a case-control association study of 75,156 European women

Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07-2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08-1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.

Clonal hematopoiesis and risk of prostate cancer in large samples of European ancestry men

Little is known regarding the potential relationship between clonal hematopoiesis (CH) of indeterminate potential (CHIP), which is the expansion of hematopoietic stem cells with somatic mutations, and risk of prostate cancer, the fifth leading cause of cancer death of men worldwide. We evaluated the association of age-related CHIP with overall and aggressive prostate cancer risk in two large whole-exome sequencing studies of 75 047 European ancestry men, including 7663 prostate cancer cases, 2770 of which had aggressive disease, and 3266 men carrying CHIP variants. We found that CHIP, defined by over 50 CHIP genes individually and in aggregate, was not significantly associated with overall (aggregate HR = 0.93, 95% CI = 0.76-1.13, P = 0.46) or aggressive (aggregate OR = 1.14, 95% CI = 0.92-1.41, P = 0.22) prostate cancer risk. CHIP was weakly associated with genetic risk of overall prostate cancer, measured using a polygenic risk score (OR = 1.05 per unit increase, 95% CI = 1.01-1.10, P = 0.01). CHIP was not significantly associated with carrying pathogenic/likely pathogenic/deleterious variants in DNA repair genes, which have previously been found to be associated with aggressive prostate cancer. While findings from this study suggest that CHIP is likely not a risk factor for prostate cancer, it will be important to investigate other types of CH in association with prostate cancer risk.

Epigenome-wide association study of short-term temperature fluctuations based on within-sibship analyses in Australian females

BACKGROUND

Temperature fluctuations can affect human health independent of the effect of mean temperature. However, no study has evaluated whether short-term temperature fluctuations could affect DNA methylation.

METHODS

Peripheral blood DNA methylation for 479 female siblings of 130 families were analysed. Gridded daily temperatures data were obtained, linked to each participant's home address, and used to calculate nine different metrics of short-term temperature fluctuations: temperature variabilities (TVs) within the day of blood draw and preceding one to seven days (TV 0-1 to TV 0-7), diurnal temperature range (DTR), and temperature change between neighbouring days (TCN). Within-sibship design was used to perform epigenome-wide association analyses, adjusting for daily mean temperatures, and other important covariates (e.g., smoking, alcohol use, cell-type proportions). Differentially methylated regions (DMRs) were further identified. Multiple-testing comparisons with a significant threshold of 0.01 for cytosine-guanine dinucleotides (CpGs) and 0.05 for DMRs were applied.

RESULTS

Among 479 participants (mean age ± SD, 56.4 ± 7.9 years), we identified significant changes in methylation levels in 14 CpGs and 70 DMRs associated with temperature fluctuations. Almost all identified CpGs were associated with exposure to temperature fluctuations within three days. Differentially methylated signals were mapped to 68 genes that were linked to human diseases such as cancer (e.g., colorectal carcinoma, breast carcinoma, and metastatic neoplasms) and mental disorder (e.g., schizophrenia, mental depression, and bipolar disorder). The top three most significantly enriched gene ontology terms were Response to bacterium (TV 0-3), followed by Hydrolase activity, acting on ester bonds (TCN), and Oxidoreductase activity (TV 0-3).

CONCLUSIONS

Short-term temperature fluctuations were associated with differentially methylated signals across the human genome, which provides evidence on the potential biological mechanisms underlying the health impact of temperature fluctuations. Future studies are needed to further clarify the roles of DNA methylation in diseases associated with temperature fluctuations.

Wildfire-related PM2.5 and DNA methylation: An Australian twin and family study

BACKGROUND

Wildfire-related fine particulate matter (PM_2.5) has many adverse health impacts, but its impacts on human epigenome are unknown. We aimed to evaluate the associations between long-term exposure to wildfire-related PM_2.5 and blood DNA methylation, and whether the associations differ from those with non-wildfire-related PM_2.5.

METHODS

We studied 479 Australian women comprising 132 twin pairs and 215 of their sisters. Blood-derived DNA methylation was measured using the HumanMethylation450 BeadChip array. Data on 3-year (year of blood collection and previous two years) average wildfire-related and non-wildfire-related PM_2.5 at 0.01°×0.01° spatial resolution were created by combining information from satellite observations, chemical transport models, and ground-based observations. Exposure data were linked to each participant's home address, assuming the address did not change during the exposure window. For DNA methylation of each cytosine-guanine dinucleotide (CpG), and for global DNA methylation represented by the average of all measured CpGs or CpGs in repetitive elements, we evaluated their associations with wildfire- or non-wildfire-related PM_2.5 using a within-sibship analysis controlling for factors shared between siblings and other important covariates. Differentially methylated regions (DMRs) were defined by comb-p and DMRcate.

RESULTS

The 3-year average wildfire-related PM_2.5 (range: 0.3 to 7.6 µg/m³_, mean: 1.6 µg/m³) was negatively, but not significantly (p-values greater than 0.05) associated with all seven global DNA methylation measures. There were 26 CpGs and 33 DMRs associated with wildfire-related PM_2.5 (Bonferroni adjusted p-value < 0.05) mapped to 47 genes enriched for pathways related to inflammatory regulation and platelet activation. These genes have been related to many human diseases or phenotypes e.g., cancer, mental disorders, diabetes, obesity, asthma, blood pressure. These CpGs, DMRs and enriched pathways did not overlap with the 1 CpG and 7 DMRs associated with non-wildfire-related PM_2.5.

CONCLUSIONS

Long-term exposure to wildfire-related PM_2.5 was associated with various blood DNA methylation signatures in Australian women, and these were distinct from those associated with non-wildfire-related PM_2.5.

Australian genome-wide association study confirms higher female risk for adult glioma associated with variants in the region of CCDC26

BACKGROUND

Glioma accounts for approximately 80% of malignant adult brain cancer and its most common subtype, glioblastoma, has one of the lowest 5-year cancer survivals. Fifty risk-associated variants within 34 glioma genetic risk regions have been found by genome-wide association studies (GWAS) with a sex difference reported for 8q24.21 region. We conducted an Australian GWAS by glioma subtype and sex.

METHODS

We analysed genome-wide data from the Australian Genomics and Clinical Outcomes of Glioma (AGOG) consortium for 7,573,692 single nucleotide polymorphisms (SNPs) for 560 glioma cases and 2,237 controls of European ancestry. Cases were classified as glioblastoma, non-glioblastoma, astrocytoma or oligodendroglioma Logistic regression analysis was used to assess the associations of SNPs with glioma risk by subtype and by sex.

RESULTS

We replicated the previously reported glioma risk associations in the regions of 2q33.3 C2orf80, 2q37.3 D2HGDH, 5p15.33 TERT, 7p11.2 EGFR, 8q24.21 CCDC26, 9p21.3 CDKN2BAS, 11q21 MAML2, 11q23.3 PHLDB1, 15q24.2 ETFA, 16p13.3 RHBDF1, 16p13.3 LMF1, 17p13.1 TP53, 20q13.33 RTEL and 20q13.33 GMEB2 (P<0.05). We also replicated the previously reported sex difference at 8q24.21 CCDC26 (P=0.0024) with the association being nominally significant for both sexes (P<0.05).

CONCLUSIONS

Our study supports a stronger female risk association for the region 8q24.21 CCDC26 and highlights the importance of analysing glioma GWAS by sex. A better understanding of sex differences could provide biological insight into the cause of glioma with implications for prevention, risk prediction and treatment.

Distinct Reproductive Risk Profiles for Intrinsic-Like Breast Cancer Subtypes: Pooled Analysis of Population-Based Studies

BACKGROUND

Reproductive factors have been shown to be differentially associated with risk of estrogen receptor (ER)-positive and ER-negative breast cancer. However, their associations with intrinsic-like subtypes are less clear.

METHODS

Analyses included up to 23 353 cases and 71 072 controls pooled from 31 population-based case-control or cohort studies in the Breast Cancer Association Consortium across 16 countries on 4 continents. Polytomous logistic regression was used to estimate the association between reproductive factors and risk of breast cancer by intrinsic-like subtypes (luminal A-like, luminal B-like, luminal B-HER2-like, HER2-enriched-like, and triple-negative breast cancer) and by invasiveness. All statistical tests were 2-sided.

RESULTS

Compared with nulliparous women, parous women had a lower risk of luminal A-like, luminal B-like, luminal B-HER2-like, and HER2-enriched-like disease. This association was apparent only after approximately 10 years since last birth and became stronger with increasing time (odds ratio [OR] = 0.59, 95% confidence interval [CI] = 0.49 to 0.71; and OR = 0.36, 95% CI = 0.28 to 0.46 for multiparous women with luminal A-like tumors 20 to less than 25 years after last birth and 45 to less than 50 years after last birth, respectively). In contrast, parous women had a higher risk of triple-negative breast cancer right after their last birth (for multiparous women: OR = 3.12, 95% CI = 2.02 to 4.83) that was attenuated with time but persisted for decades (OR = 1.03, 95% CI = 0.79 to 1.34, for multiparous women 25 to less than 30 years after last birth). Older age at first birth (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) and breastfeeding (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) were associated with lower risk of triple-negative breast cancer but not with other disease subtypes. Younger age at menarche was associated with higher risk of all subtypes; older age at menopause was associated with higher risk of luminal A-like but not triple-negative breast cancer. Associations for in situ tumors were similar to luminal A-like.

CONCLUSIONS

This large and comprehensive study demonstrates a distinct reproductive risk factor profile for triple-negative breast cancer compared with other subtypes, with implications for the understanding of disease etiology and risk prediction.

Distinct germline genetic susceptibility profiles identified for common non-Hodgkin lymphoma subtypes

Lymphoma risk is elevated for relatives with common non-Hodgkin lymphoma (NHL) subtypes, suggesting shared genetic susceptibility across subtypes. To evaluate the extent of mutual heritability among NHL subtypes and discover novel loci shared among subtypes, we analyzed data from eight genome-wide association studies within the InterLymph Consortium, including 10,629 cases and 9505 controls. We utilized Association analysis based on SubSETs (ASSET) to discover loci for subsets of NHL subtypes and evaluated shared heritability across the genome using Genome-wide Complex Trait Analysis (GCTA) and polygenic risk scores. We discovered 17 genome-wide significant loci (P < 5 × 10-8) for subsets of NHL subtypes, including a novel locus at 10q23.33 (HHEX) (P = 3.27 × 10-9). Most subset associations were driven primarily by only one subtype. Genome-wide genetic correlations between pairs of subtypes varied broadly from 0.20 to 0.86, suggesting substantial heterogeneity in the extent of shared heritability among subtypes. Polygenic risk score analyses of established loci for different lymphoid malignancies identified strong associations with some NHL subtypes (P < 5 × 10-8), but weak or null associations with others. Although our analyses suggest partially shared heritability and biological pathways, they reveal substantial heterogeneity among NHL subtypes with each having its own distinct germline genetic architecture.

Does genetic predisposition modify the effect of lifestyle-related factors on DNA methylation?

Lifestyle-related phenotypes have been shown to be heritable and associated with DNA methylation. We aimed to investigate whether genetic predisposition to tobacco smoking, alcohol consumption, and higher body mass index (BMI) moderates the effect of these phenotypes on blood DNA methylation. We calculated polygenic scores (PGS) to quantify genetic predisposition to these phenotypes using training (N = 7,431) and validation (N = 4,307) samples. Using paired genetic-methylation data (N = 4,307), gene-environment interactions (i.e., PGS × lifestyle) were assessed using linear mixed-effects models with outcomes: 1) methylation at sites found to be strongly associated with smoking (1,061 CpGs), alcohol consumption (459 CpGs), and BMI (85 CpGs) and 2) two epigenetic ageing measures, PhenoAge and GrimAge. In the validation sample, PGS explained ~1.4% (P = 1 × 10-14), ~0.6% (P = 2 × 10-7), and ~8.7% (P = 7 × 10-87) of variance in smoking initiation, alcohol consumption, and BMI, respectively. Nominally significant interaction effects (P < 0.05) were found at 61, 14, and 7 CpGs for smoking, alcohol consumption, and BMI, respectively. There was strong evidence that all lifestyle-related phenotypes were positively associated with PhenoAge and GrimAge, except for alcohol consumption with PhenoAge. There was weak evidence that the association of smoking with GrimAge was attenuated in participants genetically predisposed to smoking (interaction term: -0.022, standard error [SE] = 0.012, P = 0.058) and that the association of alcohol consumption with PhenoAge was attenuated in those genetically predisposed to drink alcohol (interaction term: -0.030, SE = 0.015, P = 0.041). In conclusion, genetic susceptibility to unhealthy lifestyles did not strongly modify the association between observed lifestyle behaviour and blood DNA methylation. Potential associations were observed for epigenetic ageing measures, which should be replicated in additional studies.

PredictCBC-2.0: a contralateral breast cancer risk prediction model developed and validated in ~ 200,000 patients

BACKGROUND

Prediction of contralateral breast cancer (CBC) risk is challenging due to moderate performances of the known risk factors. We aimed to improve our previous risk prediction model (PredictCBC) by updated follow-up and including additional risk factors.

METHODS

We included data from 207,510 invasive breast cancer patients participating in 23 studies. In total, 8225 CBC events occurred over a median follow-up of 10.2 years. In addition to the previously included risk factors, PredictCBC-2.0 included CHEK2 c.1100delC, a 313 variant polygenic risk score (PRS-313), body mass index (BMI), and parity. Fine and Gray regression was used to fit the model. Calibration and a time-dependent area under the curve (AUC) at 5 and 10 years were assessed to determine the performance of the models. Decision curve analysis was performed to evaluate the net benefit of PredictCBC-2.0 and previous PredictCBC models.

RESULTS

The discrimination of PredictCBC-2.0 at 10 years was higher than PredictCBC with an AUC of 0.65 (95% prediction intervals (PI) 0.56-0.74) versus 0.63 (95%PI 0.54-0.71). PredictCBC-2.0 was well calibrated with an observed/expected ratio at 10 years of 0.92 (95%PI 0.34-2.54). Decision curve analysis for contralateral preventive mastectomy (CPM) showed the potential clinical utility of PredictCBC-2.0 between thresholds of 4 and 12% 10-year CBC risk for BRCA1/2 mutation carriers and non-carriers.

CONCLUSIONS

Additional genetic information beyond BRCA1/2 germline mutations improved CBC risk prediction and might help tailor clinical decision-making toward CPM or alternative preventive strategies. Identifying patients who benefit from CPM, especially in the general breast cancer population, remains challenging.

Segregation analysis of 17,425 population-based breast cancer families: Evidence for genetic susceptibility and risk prediction

Rare pathogenic variants in known breast cancer-susceptibility genes and known common susceptibility variants do not fully explain the familial aggregation of breast cancer. To investigate plausible genetic models for the residual familial aggregation, we studied 17,425 families ascertained through population-based probands, 86% of whom were screened for pathogenic variants in BRCA1, BRCA2, PALB2, CHEK2, ATM, and TP53 via gene-panel sequencing. We conducted complex segregation analyses and fitted genetic models in which breast cancer incidence depended on the effects of known susceptibility genes and other unidentified major genes and a normally distributed polygenic component. The proportion of familial variance explained by the six genes was 46% at age 20-29 years and decreased steadily with age thereafter. After allowing for these genes, the best fitting model for the residual familial variance included a recessive risk component with a combined genotype frequency of 1.7% (95% CI: 0.3%-5.4%) and a penetrance to age 80 years of 69% (95% CI: 38%-95%) for homozygotes, which may reflect the combined effects of multiple variants acting in a recessive manner, and a polygenic variance of 1.27 (95% CI: 0.94%-1.65), which did not vary with age. The proportion of the residual familial variance explained by the recessive risk component was 40% at age 20-29 years and decreased with age thereafter. The model predicted age-specific familial relative risks consistent with those observed by large epidemiological studies. The findings have implications for strategies to identify new breast cancer-susceptibility genes and improve disease-risk prediction, especially at a young age.

A validation of models for prediction of pathogenic variants in mismatch repair genes

PURPOSE

Models used to predict the probability of an individual having a pathogenic homozygous or heterozygous variant in a mismatch repair gene, such as MMRpro, are widely used. Recently, MMRpro was updated with new colorectal cancer penetrance estimates. The purpose of this study was to evaluate the predictive performance of MMRpro and other models for individuals with a family history of colorectal cancer.

METHODS

We performed a validation study of 4 models, Leiden, MMRpredict, PREMM5, and MMRpro, using 784 members of clinic-based families from the United States. Predicted probabilities were compared with germline testing results and evaluated for discrimination, calibration, and predictive accuracy. We analyzed several strategies to combine models and improve predictive performance.

RESULTS

MMRpro with additional tumor information (MMRpro+) and PREMM5 outperformed the other models in discrimination and predictive accuracy. MMRpro+ was the best calibrated with an observed to expected ratio of 0.98 (95% CI = 0.89-1.08). The combination models showed improvement over PREMM5 and performed similar to MMRpro+.

CONCLUSION

MMRpro+ and PREMM5 performed well in predicting the probability of having a pathogenic homozygous or heterozygous variant in a mismatch repair gene. They serve as useful clinical decision tools for identifying individuals who would benefit greatly from screening and prevention strategies.

Physical activity, sedentary time and breast cancer risk: a Mendelian randomisation study

OBJECTIVES

Physical inactivity and sedentary behaviour are associated with higher breast cancer risk in observational studies, but ascribing causality is difficult. Mendelian randomisation (MR) assesses causality by simulating randomised trial groups using genotype. We assessed whether lifelong physical activity or sedentary time, assessed using genotype, may be causally associated with breast cancer risk overall, pre/post-menopause, and by case-groups defined by tumour characteristics.

METHODS

We performed two-sample inverse-variance-weighted MR using individual-level Breast Cancer Association Consortium case-control data from 130 957 European-ancestry women (69 838 invasive cases), and published UK Biobank data (n=91 105-377 234). Genetic instruments were single nucleotide polymorphisms (SNPs) associated in UK Biobank with wrist-worn accelerometer-measured overall physical activity (nsnps=5) or sedentary time (nsnps=6), or accelerometer-measured (nsnps=1) or self-reported (nsnps=5) vigorous physical activity.

RESULTS

Greater genetically-predicted overall activity was associated with lower breast cancer overall risk (OR=0.59; 95% confidence interval (CI) 0.42 to 0.83 per-standard deviation (SD;~8 milligravities acceleration)) and for most case-groups. Genetically-predicted vigorous activity was associated with lower risk of pre/perimenopausal breast cancer (OR=0.62; 95% CI 0.45 to 0.87,≥3 vs. 0 self-reported days/week), with consistent estimates for most case-groups. Greater genetically-predicted sedentary time was associated with higher hormone-receptor-negative tumour risk (OR=1.77; 95% CI 1.07 to 2.92 per-SD (~7% time spent sedentary)), with elevated estimates for most case-groups. Results were robust to sensitivity analyses examining pleiotropy (including weighted-median-MR, MR-Egger).

CONCLUSION

Our study provides strong evidence that greater overall physical activity, greater vigorous activity, and lower sedentary time are likely to reduce breast cancer risk. More widespread adoption of active lifestyles may reduce the burden from the most common cancer in women.

Methylation-based markers of aging and lifestyle-related factors and risk of breast cancer: a pooled analysis of four prospective studies

BACKGROUND

DNA methylation in blood may reflect adverse exposures accumulated over the lifetime and could therefore provide potential improvements in the prediction of cancer risk. A substantial body of research has shown associations between epigenetic aging and risk of disease, including cancer. Here we aimed to study epigenetic measures of aging and lifestyle-related factors in association with risk of breast cancer.

METHODS

Using data from four prospective case-control studies nested in three cohorts of European ancestry participants, including a total of 1,655 breast cancer cases, we calculated three methylation-based measures of lifestyle factors (body mass index [BMI], tobacco smoking and alcohol consumption) and seven measures of epigenetic aging (Horvath-based, Hannum-based, PhenoAge and GrimAge). All measures were regression-adjusted for their respective risk factors and expressed per standard deviation (SD). Odds ratios (OR) and 95% confidence intervals (CI) were calculated using conditional or unconditional logistic regression and pooled using fixed-effects meta-analysis. Subgroup analyses were conducted by age at blood draw, time from blood sample to diagnosis, oestrogen receptor-positivity status and tumour stage.

RESULTS

None of the measures of epigenetic aging were associated with risk of breast cancer in the pooled analysis: Horvath 'age acceleration' (AA): OR per SD = 1.02, 95%CI: 0.95-1.10; AA-Hannum: OR = 1.03, 95%CI:0.95-1.12; PhenoAge: OR = 1.01, 95%CI: 0.94-1.09 and GrimAge: OR = 1.03, 95%CI: 0.94-1.12, in models adjusting for white blood cell proportions, body mass index, smoking and alcohol consumption. The BMI-adjusted predictor of BMI was associated with breast cancer risk, OR per SD = 1.09, 95%CI: 1.01-1.17. The results for the alcohol and smoking methylation-based predictors were consistent with a null association. Risk did not appear to substantially vary by age at blood draw, time to diagnosis or tumour characteristics.

CONCLUSION

We found no evidence that methylation-based measures of aging, smoking or alcohol consumption were associated with risk of breast cancer. A methylation-based marker of BMI was associated with risk and may provide insights into the underlying associations between BMI and breast cancer.

Inflammation and Epigenetic Aging Are Largely Independent Markers of Biological Aging and Mortality

Limited evidence exists on the link between inflammation and epigenetic aging. We aimed to (a) assess the cross-sectional and prospective associations of 22 inflammation-related plasma markers and a signature of inflammaging with epigenetic aging and (b) determine whether epigenetic aging and inflammaging are independently associated with mortality. Blood samples from 940 participants in the Melbourne Collaborative Cohort Study collected at baseline (1990-1994) and follow-up (2003-2007) were assayed for DNA methylation and 22 inflammation-related markers, including well-established markers (eg, interleukins and C-reactive protein) and metabolites of the tryptophan-kynurenine pathway. Four measures of epigenetic aging (PhenoAge, GrimAge, DunedinPoAm, and Zhang) and a signature of inflammaging were considered, adjusted for age, and transformed to Z scores. Associations were assessed using linear regression, and mortality hazard ratios (HR) and 95% confidence intervals (95% CI) were estimated using Cox regression. Cross-sectionally, most inflammation-related markers were associated with epigenetic aging measures, although with generally modest effect sizes (regression coefficients per SD ≤ 0.26) and explaining altogether between 1% and 11% of their variation. Prospectively, baseline inflammation-related markers were not, or only weakly, associated with epigenetic aging after 11 years of follow-up. Epigenetic aging and inflammaging were strongly and independently associated with mortality, for example, inflammaging: HR = 1.41, 95% CI = 1.27-1.56, p = 2 × 10-10, which was only slightly attenuated after adjustment for 4 epigenetic aging measures: HR = 1.35, 95% CI = 1.22-1.51, p = 7 × 10-9). Although cross-sectionally associated with epigenetic aging, inflammation-related markers accounted for a modest proportion of its variation. Inflammaging and epigenetic aging are essentially nonoverlapping markers of biological aging and may be used jointly to predict mortality.

Abstract 688: Multi-stage exome sequencing study of 17,546 aggressive and non-aggressive prostate cancer cases

Rare pathogenic variants in DNA repair genes have been found to influence risk of aggressive prostate cancer. We conducted a large case-only exome sequencing study to further understand the role of rare coding variation in aggressive prostate cancer in a study of 9,185 aggressive (prostate cancer death, metastatic disease, T4, or both T3 and Gleason≥8) and 8,361 non-aggressive cases (T1/T2 and Gleason≤6) of European ancestry from 19 international studies. Stage 1 samples (N=5,545) had whole exome-sequencing, and stage 2 samples (N=12,001) had targeted exome sequencing for 1,459 genes selected based on stage 1 results and previous evidence. Logistic regression models were used to evaluate gene-based tests and the aggregate effect of multiple genes to investigate whether carrying pathogenic/likely pathogenic/deleterious (P/LP/D) variants (18,759 identified) was associated with risk of aggressive prostate cancer, prostate cancer death (N=6,033), or metastatic disease (N=1,730) compared to non-aggressive disease. Gene-based tests were meta-analyzed across stages 1 and 2. BRCA2, ATM, and NBN had the most statistically significant gene-based results: BRCA2 P/LP/D variant carriers had 4.3-fold higher odds of aggressive disease (95% CI=3.15-5.86, P=4x10-20), 4.7-fold higher odds of prostate cancer death (95% CI=3.41-6.59, P=2x10-20), and 5.7-fold higher odds of metastatic disease (95% CI=3.71-8.76, P=2x10-15); ATM P/LP/D variant carriers had 2.2-fold higher odds of aggressive disease (95% CI=1.58-2.99, P=2x10-6), 2.2-fold higher odds of prostate cancer death (95% CI=1.52-3.05, P=2x10-5), and 3.0-fold higher odds of metastatic disease (95% CI=1.93-4.61, P=9x10-7); and NBN P/LP/D variant carriers had 5.9-fold higher odds of metastatic disease (95% CI=2.56-13.84, P=3x10-5). Among potentially novel genes with strong but not exome-wide significant statistical evidence were MMP19, involved in reproduction and metastasis, with carriers having 2.8-fold higher odds of prostate cancer death (95% CI=1.53-5.05, P=8x10-4); PKD2L2, involved in fertility, with carriers having 3.5-fold higher odds of prostate cancer death (95% CI=1.76-7.04, P=5x10-4); and SMPD1, involved in converting sphingomyelin to ceramide, with carriers having 5.3-fold higher odds of metastatic disease (95% CI=1.85-14.98, P=0.002). At least one P/LP/D variant within a subset of 24 previously curated candidate prostate cancer DNA repair genes was carried by 12.8% of aggressive cases (OR=1.48, 95% CI=1.34-1.64, P=3x10-14), 12.6% of cases who died due to prostate cancer (OR=1.47, 95% CI=1.31-1.65, P=3x10-11), and 15.1% of metastatic cases (OR=2.16, 95% CI=1.57-2.16, P=5x10-14) compared to 9.4% of non-aggressive cases. These findings support the importance of rare genetic variation in aggressive prostate cancer risk and may have important implications for prostate cancer risk stratification and screening.

Citation Format: Burcu F. Darst, Ed Saunders, Tokhir Dadaev, Xin Sheng, Peggy Wan, Loreall Pooler, Lucy Y. Xia, Stephen Chanock, Sonja I. Berndt, Susan M. Gapstur, Victoria Stevens, Demetrius Albanes, Stephanie J. Weinstein, Vincent Gnanapragasam, Graham G. Giles, Tu Nguyen-Dumont, Roger L. Milne, Mark M. Pomerantz, Julie A. Schmidt, Ruth C. Travis, Timothy J. Key, Konrad H. Stopsack, Lorelei A. Mucci, William J. Catalona, Beth Marosy, Kurt N. Hetrick, Kimberly F. Doheny, Robert J. MacInnis, Melissa C. Southey, Rosalind A. Eeles, Fredrik Wiklund, Zsofia Kote-Jarai, David V. Conti, Christopher A. Haiman. Multi-stage exome sequencing study of 17,546 aggressive and non-aggressive prostate cancer cases [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 688.

Genetic Aspects of Mammographic Density Measures Associated with Breast Cancer Risk

Cumulus, Altocumulus, and Cirrocumulus are measures of mammographic density defined at increasing pixel brightness thresholds, which, when converted to mammogram risk scores (MRSs), predict breast cancer risk. Twin and family studies suggest substantial variance in the MRSs could be explained by genetic factors. For 2559 women aged 30 to 80 years (mean 54 years), we measured the MRSs from digitized film mammograms and estimated the associations of the MRSs with a 313-SNP breast cancer polygenic risk score (PRS) and 202 individual SNPs associated with breast cancer risk. The PRS was weakly positively correlated (correlation coefficients ranged 0.05−0.08; all p < 0.04) with all the MRSs except the Cumulus-white MRS based on the “white but not bright area” (correlation coefficient = 0.04; p = 0.06). After adjusting for its association with the Altocumulus MRS, the PRS was not associated with the Cumulus MRS. There were MRS associations (Bonferroni-adjusted p < 0.04) with one SNP in the ATXN1 gene and nominally with some ESR1 SNPs. Less than 1% of the variance of the MRSs is explained by the genetic markers currently known to be associated with breast cancer risk. Discovering the genetic determinants of the bright, not white, regions of the mammogram could reveal substantial new genetic causes of breast cancer.

Breast cancer risks associated with missense variants in breast cancer susceptibility genes

BACKGROUND

Protein truncating variants in ATM, BRCA1, BRCA2, CHEK2, and PALB2 are associated with increased breast cancer risk, but risks associated with missense variants in these genes are uncertain.

METHODS

We analyzed data on 59,639 breast cancer cases and 53,165 controls from studies participating in the Breast Cancer Association Consortium BRIDGES project. We sampled training (80%) and validation (20%) sets to analyze rare missense variants in ATM (1146 training variants), BRCA1 (644), BRCA2 (1425), CHEK2 (325), and PALB2 (472). We evaluated breast cancer risks according to five in silico prediction-of-deleteriousness algorithms, functional protein domain, and frequency, using logistic regression models and also mixture models in which a subset of variants was assumed to be risk-associated.

RESULTS

The most predictive in silico algorithms were Helix (BRCA1, BRCA2 and CHEK2) and CADD (ATM). Increased risks appeared restricted to functional protein domains for ATM (FAT and PIK domains) and BRCA1 (RING and BRCT domains). For ATM, BRCA1, and BRCA2, data were compatible with small subsets (approximately 7%, 2%, and 0.6%, respectively) of rare missense variants giving similar risk to those of protein truncating variants in the same gene. For CHEK2, data were more consistent with a large fraction (approximately 60%) of rare missense variants giving a lower risk (OR 1.75, 95% CI (1.47-2.08)) than CHEK2 protein truncating variants. There was little evidence for an association with risk for missense variants in PALB2. The best fitting models were well calibrated in the validation set.

CONCLUSIONS

These results will inform risk prediction models and the selection of candidate variants for functional assays and could contribute to the clinical reporting of gene panel testing for breast cancer susceptibility.

Cancer Risks Associated With BRCA1 and BRCA2 Pathogenic Variants

PURPOSE

To provide precise age-specific risk estimates of cancers other than female breast and ovarian cancers associated with pathogenic variants (PVs) in BRCA1 and BRCA2 for effective cancer risk management.

METHODS

We used data from 3,184 BRCA1 and 2,157 BRCA2 families in the Consortium of Investigators of Modifiers of BRCA1/2 to estimate age-specific relative (RR) and absolute risks for 22 first primary cancer types adjusting for family ascertainment.

RESULTS

BRCA1 PVs were associated with risks of male breast (RR = 4.30; 95% CI, 1.09 to 16.96), pancreatic (RR = 2.36; 95% CI, 1.51 to 3.68), and stomach (RR = 2.17; 95% CI, 1.25 to 3.77) cancers. Associations with colorectal and gallbladder cancers were also suggested. BRCA2 PVs were associated with risks of male breast (RR = 44.0; 95% CI, 21.3 to 90.9), stomach (RR = 3.69; 95% CI, 2.40 to 5.67), pancreatic (RR = 3.34; 95% CI, 2.21 to 5.06), and prostate (RR = 2.22; 95% CI, 1.63 to 3.03) cancers. The stomach cancer RR was higher for females than males (6.89 v 2.76; P = .04). The absolute risks to age 80 years ranged from 0.4% for male breast cancer to approximately 2.5% for pancreatic cancer for BRCA1 carriers and from approximately 2.5% for pancreatic cancer to 27% for prostate cancer for BRCA2 carriers.

CONCLUSION

In addition to female breast and ovarian cancers, BRCA1 and BRCA2 PVs are associated with increased risks of male breast, pancreatic, stomach, and prostate (only BRCA2 PVs) cancers, but not with the risks of other previously suggested cancers. The estimated age-specific risks will refine cancer risk management in men and women with BRCA1/2 PVs.

Independent evaluation of melanoma polygenic risk scores in UK and Australian prospective cohorts

BACKGROUND

Previous studies suggest that polygenic risk scores (PRSs) may improve melanoma risk stratification. However, there has been limited independent validation of PRS-based risk prediction, particularly assessment of calibration (comparing predicted to observed risks).

OBJECTIVES

To evaluate PRS-based melanoma risk prediction in prospective UK and Australian cohorts with European ancestry.

METHODS

We analysed invasive melanoma incidence in the UK Biobank (UKB; n = 395 647, 1651 cases) and a case-cohort nested within the Melbourne Collaborative Cohort Study (MCCS, Australia; n = 4765, 303 cases). Three PRSs were evaluated: 68 single-nucleotide polymorphisms (SNPs) at 54 loci from a 2020 meta-analysis (PRS68), 50 SNPs significant in the 2020 meta-analysis excluding UKB (PRS50) and 45 SNPs at 21 loci known in 2018 (PRS45). Ten-year melanoma risks were calculated from population-level cancer registry data by age group and sex, with and without PRS adjustment.

RESULTS

Predicted absolute melanoma risks based on age and sex alone underestimated melanoma incidence in the UKB [ratio of expected/observed cases: E/O = 0·65, 95% confidence interval (CI) 0·62-0·68] and MCCS (E/O = 0·63, 95% CI 0·56-0·72). For UKB, calibration was improved by PRS adjustment, with PRS50-adjusted risks E/O = 0·91, 95% CI 0·87-0·95. The discriminative ability for PRS68- and PRS50-adjusted absolute risks was higher than for risks based on age and sex alone (Δ area under the curve 0·07-0·10, P < 0·0001), and higher than for PRS45-adjusted risks (Δ area under the curve 0·02-0·04, P < 0·001).

CONCLUSIONS

A PRS derived from a larger, more diverse meta-analysis improves risk prediction compared with an earlier PRS, and might help tailor melanoma prevention and early detection strategies to different risk levels. Recalibration of absolute risks may be necessary for application to specific populations.