Co-observation of germline pathogenic variants in breast cancer predisposition genes: Results from analysis of the BRIDGES sequencing dataset.

Co-observation of a gene variant with a pathogenic variant in another gene that explains the disease presentation has been designated as evidence against pathogenicity for commonly used variant classification guidelines. Multiple variant curation expert panels have specified, from consensus opinion, that this evidence type is not applicable for the classification of breast cancer predisposition gene variants. Statistical analysis of sequence data for 55,815 individuals diagnosed with breast cancer from the BRIDGES sequencing project was undertaken to formally assess the utility of co-observation data for germline variant classification. Our analysis included expected loss-of-function variants in 11 breast cancer predisposition genes and pathogenic missense variants in BRCA1, BRCA2, and TP53. We assessed whether co-observation of pathogenic variants in two different genes occurred more or less often than expected under the assumption of independence. Co-observation of pathogenic variants in each of BRCA1, BRCA2, and PALB2 with the remaining genes was less frequent than expected. This evidence for depletion remained after adjustment for age at diagnosis, study design (familial versus population-based), and country. Co-observation of a variant of uncertain significance in BRCA1, BRCA2, or PALB2 with a pathogenic variant in another breast cancer gene equated to supporting evidence against pathogenicity following criterion strength assignment based on the likelihood ratio and showed utility in reclassification of missense BRCA1 and BRCA2 variants identified in BRIDGES. Our approach has applicability for assessing the value of co-observation as a predictor of variant pathogenicity in other clinical contexts, including for gene-specific guidelines developed by ClinGen Variant Curation Expert Panels.

Cancer risks among first-degree relatives of women with a genetic predisposition to breast cancer.

BACKGROUND: Associations between germline alterations in women and cancer risks among their relatives are largely unknown. METHODS: We identified women from 2 Swedish cohorts Karolinska Mammography Project for Risk Prediction of Breast Cancer (KARMA) and prevalent KARMA (pKARMA), including 28 362 women with genotyping data and 13 226 with sequencing data. Using Swedish Multi-Generation Register, we linked these women to 133 389 first-degree relatives. Associations between protein-truncating variants in 8 risk genes and breast cancer polygenic risk score in index women and cancer risks among their relatives were modeled via Cox regression. RESULTS: Female relatives of index women who were protein-truncating variant carriers in any of the 8 risk genes had an increased breast cancer risk compared with those of noncarriers (hazard ratio [HR] = 1.85, 95% confidence interval [CI] = 1.52 to 2.27), with the strongest association found for protein-truncating variants in BRCA1 and 2. These relatives had a statistically higher risk of early onset than late-onset breast cancer (P = .001). Elevated breast cancer risk was also observed in female relatives of index women with higher polygenic risk score (HR per SD = 1.28, 95% CI = 1.23 to 1.32). The estimated lifetime risk was 22.3% for female relatives of protein-truncating variant carriers and 14.4% for those related to women in the top polygenic risk score quartile. Moreover, relatives of index women with protein-truncating variant presence (HR = 1.30, 95% CI = 1.06 to 1.59) or higher polygenic risk score (HR per SD = 1.04, 95% CI = 1.01 to 1.07) were also at higher risk of nonbreast hereditary breast and ovary cancer syndrome-related cancers. CONCLUSIONS: Protein-truncating variants of risk genes and higher polygenic risk score in index women are associated with an increased risk of breast and other hereditary breast and ovary syndrome-related cancers among relatives.

Investigation of Genetic Alterations Associated With Interval Breast Cancer.

Importance: Breast cancers (BCs) diagnosed between 2 screening examinations are called interval cancers (ICs), and they have worse clinicopathological characteristics and poorer prognosis than screen-detected cancers (SDCs). However, the association of rare germline genetic variants with IC have not been studied. Objective: To evaluate whether rare germline deleterious protein-truncating variants (PTVs) can be applied to discriminate between IC and SDC while considering mammographic density. Design, Setting, and Participants: This population-based genetic association study was based on women aged 40 to 76 years who were attending mammographic screening in Sweden. All women with a diagnosis of BC between January 2001 and January 2016 were included, together with age-matched controls. Patients with BC were followed up for survival until 2021. Statistical analysis was performed from September 2021 to December 2022. Exposure: Germline PTVs in 34 BC susceptibility genes as analyzed by targeted sequencing. Main Outcomes and Measures: Odds ratios (ORs) were used to compare IC with SDC using logistic regression. Hazard ratios were used to investigate BC-specific survival using Cox regression. Results: All 4121 patients with BC (IC, n = 1229; SDC, n = 2892) were female, with a mean (SD) age of 55.5 (7.1) years. There were 5631 age-matched controls. The PTVs of the ATM, BRCA1, BRCA2, CHEK2, and PALB2 genes were more common in patients with IC compared with SDC (OR, 1.48; 95% CI, 1.06-2.05). This association was primarily influenced by BRCA1/2 and PALB2 variants. A family history of BC together with PTVs of any of these genes synergistically increased the probability of receiving a diagnosis of IC rather than SDC (OR, 3.95; 95% CI, 1.97-7.92). Furthermore, 10-year BC-specific survival revealed that if a patient received a diagnosis of an IC, carriers of PTVs in any of these 5 genes had significantly worse survival compared with patients not carrying any of them (hazard ratio, 2.04; 95% CI, 1.06-3.92). All of these associations were further pronounced in a subset of patients with IC who had a low mammographic density at prior screening examination. Conclusions and Relevance: The results of this study may be helpful in future optimizations of screening programs that aim to lower mortality as well as the clinical treatment of patients with BC.

The impact of circulating protein levels identified by affinity proteomics on short-term, overall breast cancer risk.

OBJECTIVE: Current breast cancer risk prediction scores and algorithms can potentially be further improved by including molecular markers. To this end, we studied the association of circulating plasma proteins using Proximity Extension Assay (PEA) with incident breast cancer risk. SUBJECTS: In this study, we included 1577 women participating in the prospective KARMA mammographic screening cohort. RESULTS: In a targeted panel of 164 proteins, we found 8 candidates nominally significantly associated with short-term breast cancer risk (P < 0.05). Similarly, in an exploratory panel consisting of 2204 proteins, 115 were found nominally significantly associated (P < 0.05). However, none of the identified protein levels remained significant after adjustment for multiple testing. This lack of statistically significant findings was not due to limited power, but attributable to the small effect sizes observed even for nominally significant proteins. Similarly, adding plasma protein levels to established risk factors did not improve breast cancer risk prediction accuracy. CONCLUSIONS: Our results indicate that the levels of the studied plasma proteins captured by the PEA method are unlikely to offer additional benefits for risk prediction of short-term overall breast cancer risk but could provide interesting insights into the biological basis of breast cancer in the future.

Evaluation of circulating plasma proteins in breast cancer using Mendelian randomisation.

Biomarkers for early detection of breast cancer may complement population screening approaches to enable earlier and more precise treatment. The blood proteome is an important source for biomarker discovery but so far, few proteins have been identified with breast cancer risk. Here, we measure 2929 unique proteins in plasma from 598 women selected from the Karolinska Mammography Project to explore the association between protein levels, clinical characteristics, and gene variants, and to identify proteins with a causal role in breast cancer. We present 812 cis-acting protein quantitative trait loci for 737 proteins which are used as instruments in Mendelian randomisation analyses of breast cancer risk. Of those, we present five proteins (CD160, DNPH1, LAYN, LRRC37A2 and TLR1) that show a potential causal role in breast cancer risk with confirmatory results in independent cohorts. Our study suggests that these proteins should be further explored as biomarkers and potential drug targets in breast cancer.

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.

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.

Prediction of breast cancer risk for sisters of women attending screening.

BACKGROUND: Risk assessment is important for breast cancer prevention and early detection. We aimed to examine whether common risk factors, mammographic features, and breast cancer risk prediction scores of a woman were associated with breast cancer risk for her sisters. METHODS: We included 53 051 women from the Karolinska Mammography Project for Risk Prediction of Breast Cancer (KARMA) study. Established risk factors were derived using self-reported questionnaires, mammograms, and single nucleotide polymorphism genotyping. Using the Swedish Multi-Generation Register, we identified 32 198 sisters of the KARMA women (including 5352 KARMA participants and 26 846 nonparticipants). Cox models were used to estimate the hazard ratios of breast cancer for both women and their sisters, respectively. RESULTS: A higher breast cancer polygenic risk score, a history of benign breast disease, and higher breast density in women were associated with an increased risk of breast cancer for both women and their sisters. No statistically significant association was observed between breast microcalcifications and masses in women and breast cancer risk for their sisters. Furthermore, higher breast cancer risk scores in women were associated with an increased risk of breast cancer for their sisters. Specifically, the hazard ratios for breast cancer per 1 standard deviation increase in age-adjusted KARMA, Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA), and Tyrer-Cuzick risk scores were 1.16 (95% confidence interval [CI] = 1.07 to 1.27), 1.23 (95% CI = 1.12 to 1.35), and 1.21 (95% CI = 1.11 to 1.32), respectively. CONCLUSION: A woman's breast cancer risk factors are associated with her sister's breast cancer risk. However, the clinical utility of these findings requires further investigation.

Risk of estrogen receptor-specific breast cancer by family history of estrogen receptor subtypes and other cancers.

BACKGROUND: The extent to which the risk of estrogen receptor (ER)-specific breast cancer is associated with ER status of breast cancer and other cancers among first-degree relatives is unclear. METHODS: This population-based cohort included 464 707 cancer-free women in Stockholm, Sweden, during 1978-2019. For ER-negative and ER-positive breast cancers, we estimated hazard ratios (HRs) associated with ER status of female first-degree relatives with breast cancer and of other cancers in all first-degree relatives. Associations between ER-negative and ER-positive status by family cancer history were estimated using logistic regression in a case-only design. RESULTS: Women with familial ER-positive breast cancer had 1.87 times (95% confidence interval [CI] = 1.77 to 1.97) higher risk of ER-positive subtype, whereas the corresponding hazard ratio for ER-negative was 2.54 (95% CI = 2.08 to 3.10) when having familial ER-negative breast cancer. The risk increased with an increasing number of female first-degree relatives having concordant subtypes and younger age at diagnosis (Ptrend <.001 for both). Nonbreast cancers among first-degree relatives were associated with both ER-positive (HR = 1.14, 95% CI = 1.10 to 1.17) and ER-negative (HR = 1.08, 95% CI = 1.01 to 1.16) breast cancers. Compared with women with ER-positive breast cancer, women with ER-negative breast cancer were more likely to have family history of liver (odds ratio [OR] = 1.33, 95% CI = 1.05 to 1.67), ovary (OR = 1.28, 95% CI = 1.01 to 1.61), and testicle cancer (OR = 1.79, 95% CI = 1.01 to 3.16) but less likely to have family history of endometrial cancer (OR = 0.77, 95% CI = 0.60 to 1.00) and leukemia (OR = 0.72, 95% CI = 0.56 to 0.91). CONCLUSIONS: Risk of ER-specific breast cancer differs according to ER status of female first-degree relatives with breast cancer and some other cancers of first-degree relatives. This family history information should be considered in the individual risk prediction for ER subtypes.

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).

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.

Percutaneous transhepatic or endoscopic ultrasound-guided biliary drainage in malignant distal bile duct obstruction using a self-expanding metal stent: Study protocol for a prospective European multicenter trial (PUMa trial).

BACKGROUND: Endoscopic ultrasound-guided biliary drainage (EUS-BD) was associated with better clinical success and a lower rate of adverse events (AEs) than fluoroscopy-guided percutaneous transhepatic biliary drainage (PTBD) in recent single center studies with mainly retrospective design and small case numbers (< 50). The aim of this prospective European multicenter study is to compare both drainage procedures using ultrasound-guidance and primary metal stent implantation in patients with malignant distal bile duct obstruction (PUMa Trial). METHODS: The study is designed as a non-randomized, controlled, parallel group, non-inferiority trial. Each of the 16 study centers performs the procedure with the best local expertise (PTBD or EUS-BD). In PTBD, bile duct access is performed by ultrasound guidance. EUS-BD is performed as an endoscopic ultrasound (EUS)-guided hepaticogastrostomy (EUS-HGS), EUS-guided choledochoduodenostomy (EUS-CDS) or EUS-guided antegrade stenting (EUS-AGS). Insertion of a metal stent is intended in both procedures in the first session. Primary end point is technical success. Secondary end points are clinical success, duration pf procedure, AEs graded by severity, length of hospital stay, re-intervention rate and survival within 6 months. The target case number is 212 patients (12 calculated dropouts included). DISCUSSION: This study might help to clarify whether PTBD is non-inferior to EUS-BD concerning technical success, and whether one of both interventions is superior in terms of efficacy and safety in one or more secondary endpoints. Randomization is not provided as both procedures are rarely used after failed endoscopic biliary drainage and study centers usually prefer one of both procedures that they can perform best. TRIAL REGISTRATION: ClinicalTrials.gov ID: NCT03546049 (22.05.2018).

Risk of heart disease following treatment for breast cancer - results from a population-based cohort study.

Background: There is a rising concern about treatment-associated cardiotoxicities in breast cancer patients. This study aimed to determine the time- and treatment-specific incidence of arrhythmia, heart failure, and ischemic heart disease in women diagnosed with breast cancer. Methods: A register-based matched cohort study was conducted including 8015 breast cancer patients diagnosed from 2001 to 2008 in the Stockholm-Gotland region and followed up until 2017. Time-dependent risks of arrhythmia, heart failure, and ischemic heart disease in breast cancer patients were assessed using flexible parametric models as compared to matched controls from general population. Treatment-specific effects were estimated in breast cancer patients using Cox model. Results: Time-dependent analyses revealed long-term increased risks of arrhythmia and heart failure following breast cancer diagnosis. Hazard ratios (HRs) within the first year of diagnosis were 2.14 (95% CI = 1.63-2.81) for arrhythmia and 2.71 (95% CI = 1.70-4.33) for heart failure. HR more than 10 years following diagnosis was 1.42 (95% CI = 1.21-1.67) for arrhythmia and 1.28 (95% CI = 1.03-1.59) for heart failure. The risk for ischemic heart disease was significantly increased only during the first year after diagnosis (HR = 1.45, 95% CI = 1.03-2.04). Trastuzumab and anthracyclines were associated with increased risk of heart failure. Aromatase inhibitors, but not tamoxifen, were associated with risk of ischemic heart disease. No increased risk of heart disease was identified following locoregional radiotherapy. Conclusions: Administration of systemic adjuvant therapies appears to be associated with increased risks of heart disease. The risk estimates observed in this study may aid adjuvant therapy decision-making and patient counseling in oncology practices. Funding: This work was supported by the Swedish Research Council (grant no: 2018-02547); Swedish Cancer Society (grant no: CAN-19-0266); and FORTE (grant no: 2016-00081).

Interval breast cancer is associated with interferon immune response.

BACKGROUND: The aggressive nature of breast cancers detected between planned mammographic screens, so-called interval cancers, remains elusive. Here, we aim to characterise underlying molecular features of interval cancer. METHODS: From 672 patients with invasive breast cancer, we analysed gene expression differences between 90 'true' interval cancer cases (i.e. women with low-dense breasts defined as per cent mammographic density <25%) and 310 screen-detected tumours while accounting for PAM50 subtypes and thus overall tumour aggressiveness. We computed an interval cancer gene expression profile (IC-Gx) in a total of 2270 breast tumours (regardless of interval cancer status) and tested for association with expression-based immune subtypes in breast cancer. In addition, we investigated the contribution of inherited and somatic genetic variants in distinct features of interval cancer. RESULTS: We identified 8331 genes nominally associated with interval cancer (P-value < 0.05, fold-change > 1.5). Gene set enrichment analysis showed immune-related pathways as key processes altered in interval cancer. Our IC-Gx, based on 47 genes with the strongest associations (false discovery rate < 0.05), was found to be associated mainly with immune subtypes involving interferon response. We isolated an interaction network of interval cancer and interferon genes for which a significant load of somatic and germline variants in class I interferon genes was observed. CONCLUSION: We identified novel molecular features of interval breast cancer highlighting interferon pathways as a potential target for prevention or treatment.

Mammographic features are associated with cardiometabolic disease risk and mortality.

AIMS: In recent years, microcalcifications identified in routine mammograms were found to be associated with cardiometabolic disease in women. Here, we aimed to systematically evaluate the association of microcalcifications and other mammographic features with cardiometabolic disease risk and mortality in a large screening cohort and to understand a potential genetic contribution. METHODS AND RESULTS: This study included 57 867 women from a prospective mammographic screening cohort in Sweden (KARMA) and 49 583 sisters. Cardiometabolic disease diagnoses and mortality and medication were extracted by linkage to Swedish population registries with virtually no missing data. In the cardiometabolic phenome-wide association study, we found that a higher number of microcalcifications were associated with increased risk for multiple cardiometabolic diseases, particularly in women with pre-existing cardiometabolic diseases. In contrast, dense breasts were associated with a lower incidence of cardiometabolic diseases. Importantly, we observed similar associations in sisters of KARMA women, indicating a potential genetic overlap between mammographic features and cardiometabolic traits. Finally, we observed that the presence of microcalcifications was associated with increased cardiometabolic mortality in women with pre-existing cardiometabolic diseases (hazard ratio and 95% confidence interval: 1.79 [1.24-2.58], P = 0.002) while we did not find such effects in women without cardiometabolic diseases. CONCLUSIONS: We found that mammographic features are associated with cardiometabolic risk and mortality. Our results strengthen the notion that a combination of mammographic features and other breast cancer risk factors could be a novel and affordable tool to assess cardiometabolic health in women attending mammographic screening.

Lymph node metastases in breast cancer: Investigating associations with tumor characteristics, molecular subtypes and polygenic risk score using a continuous growth model.

We investigate the association between rate of breast cancer lymph node spread and grade, estrogen receptor (ER) status, progesteron receptor status, decision tree derived PAM50 molecular subtype and a polygenic risk score (PRS), using data on 10 950 women included from two different data sources. Lymph node spread was analyzed using a novel continuous tumor progression model that adjusts for tumor volume in a biologically motivated way and that incorporates covariates of interest. Grades 2 and 3 tumors, respectively, were associated with 1.63 and 2.17 times faster rates of lymph node spread than Grade 1 tumors (P < 10-16 ). ER/PR negative breast cancer was associated with a 1.25/1.19 times faster spread than ER/PR positive breast cancer, respectively (P = .0011 and .0012). Among the molecular subtypes luminal A, luminal B, Her2-enriched and basal-like, Her2-enriched breast cancer was associated with 1.53 times faster spread than luminal A cancer (P = .00072). PRS was not associated with the rate of lymph node spread. Continuous growth models are useful for quantifying associations between lymph node spread and tumor characteristics. These may be useful for building realistic progression models for microsimulation studies used to design individualized screening programs.

Deciphering the genetic and epidemiological landscape of mitochondrial DNA abundance.

Mitochondrial (MT) dysfunction is a hallmark of aging and has been associated with most aging-related diseases as well as immunological processes. However, little is known about aging, lifestyle and genetic factors influencing mitochondrial DNA (mtDNA) abundance. In this study, mtDNA abundance was estimated from the weighted intensities of probes mapping to the MT genome in 295,150 participants from the UK Biobank. We found that the abundance of mtDNA was significantly elevated in women compared to men, was negatively correlated with advanced age, higher smoking exposure, greater body-mass index, higher frailty index as well as elevated red and white blood cell count and lower mortality. In addition, several biochemistry markers in blood-related to cholesterol metabolism, ion homeostasis and kidney function were found to be significantly associated with mtDNA abundance. By performing a genome-wide association study, we identified 50 independent regions genome-wide significantly associated with mtDNA abundance which harbour multiple genes involved in the immune system, cancer as well as mitochondrial function. Using mixed effects models, we estimated the SNP-heritability of mtDNA abundance to be around 8%. To investigate the consequence of altered mtDNA abundance, we performed a phenome-wide association study and found that mtDNA abundance is involved in risk for leukaemia, hematologic diseases as well as hypertension. Thus, estimating mtDNA abundance from genotyping arrays has the potential to provide novel insights into age- and disease-relevant processes, particularly those related to immunity and established mitochondrial functions.

Association between breast cancer risk and disease aggressiveness: Characterizing underlying gene expression patterns.

The association between breast cancer risk defined by the Tyrer-Cuzick score (TC) and disease prognosis is not well established. Here, we investigated the relationship between 5-year TC and disease aggressiveness and then characterized underlying molecular processes. In a case-only study (n = 2474), we studied the association of TC with molecular subtypes and tumor characteristics. In a subset of patients (n = 672), we correlated gene expression to TC and computed a low-risk TC gene expression (TC-Gx) profile, that is, a profile expected to be negatively associated with risk, which we used to test for association with disease aggressiveness. We performed enrichment analysis to pinpoint molecular processes likely to be altered in low-risk tumors. A higher TC was found to be inversely associated with more aggressive surrogate molecular subtypes and tumor characteristics (P < .05) including Ki-67 proliferation status (P < 5 × 10-07 ). Our low-risk TC-Gx, based on the weighted sum of 37 expression values of genes strongly correlated with TC, was associated with basal-like (P < 5 × 10-13 ), HER2-enriched subtype (P < 5 × 10-07 ) and worse 10-year breast cancer-specific survival (log-rank P < 5 × 10-04 ). Associations between low-risk TC-Gx and more aggressive molecular subtypes were replicated in an independent cohort from The Cancer Genome Atlas database (n = 975). Gene expression that correlated with low TC was enriched in proliferation and oncogenic signaling pathways (FDR < 0.05). Moreover, higher proliferation was a key factor explaining the association with worse survival. Women who developed breast cancer despite having a low risk were diagnosed with more aggressive tumors and had a worse prognosis, most likely driven by increased proliferation. Our findings imply the need to establish risk factors associated with more aggressive breast cancer subtypes.

A systems genomics approach to uncover the molecular properties of cancer genes.

Genes involved in cancer are under constant evolutionary pressure, potentially resulting in diverse molecular properties. In this study, we explore 23 omic features from publicly available databases to define the molecular profile of different classes of cancer genes. Cancer genes were grouped according to mutational landscape (germline and somatically mutated genes), role in cancer initiation (cancer driver genes) or cancer survival (survival genes), as well as being implicated by genome-wide association studies (GWAS genes). For each gene, we also computed feature scores based on all omic features, effectively summarizing how closely a gene resembles cancer genes of the respective class. In general, cancer genes are longer, have a lower GC content, have more isoforms with shorter exons, are expressed in more tissues and have more transcription factor binding sites than non-cancer genes. We found that germline genes more closely resemble single tissue GWAS genes while somatic genes are more similar to pleiotropic cancer GWAS genes. As a proof-of-principle, we utilized aggregated feature scores to prioritize genes in breast cancer GWAS loci and found that top ranking genes were enriched in cancer related pathways. In conclusion, we have identified multiple omic features associated with different classes of cancer genes, which can assist prioritization of genes in cancer gene discovery.