Genome-wide association study of germline variants and breast cancer-specific mortality

Additional prognostic value of polymorphisms within the 3'-untranslated region of programmed cell death pathway genes in early-stage breast cancer

Introduction

The programmed cell death (PCD) pathway plays an important role in restricting cancer cell survival and proliferation. However, limited studies have investigated the association between genetic variants in the 3'-untranslated region of the PCD pathway genes and breast cancer outcomes.

Methods

In this study, we genotyped 28 potentially functional single nucleotide polymorphisms (SNPs) in 23 PCD pathway genes in 1,177 patients with early-stage breast cancer (EBC) from a Han Chinese population. The median follow-up period was 174 months.

Results

Among all the candidate SNPs, four independent SNPs (rs4900321 and rs7150025 in ATG2B, rs6753785 in BCL2L11, and rs2213181 in c-Kit) were associated with invasive disease-free survival (iDFS), distant disease-free survival (DDFS), breast cancer-specific survival (BCSS) and overall survival (OS), respectively. Further combined genotypes of these four SNPs revealed that the survival decreased as the number of unfavorable genotypes increased (Ptrend = 1.0 × 10-6, 8.5 × 10-8, 3.6 × 10-4, and 1.3 × 10-4 for iDFS, DDFS, BCSS, and OS, respectively). Receiver operating characteristic curve analysis demonstrated that incorporating unfavorable genotypes and clinicopathological variables improved the ability to predict EBC survival (P = 0.006, 0.004, 0.029, and 0.019 for iDFS, DDFS, BCSS, and OS, respectively). Additionally, rs6753785 and rs2213181 were associated with BCL2L11 and c-Kit mRNA expression, respectively.

Conclusions

Our results suggest that these four SNPs may act as novel biomarkers for EBC survival, possibly by modulating the expression of the corresponding genes.

Machine learning reveals genetic modifiers of the immune microenvironment of cancer

Heritability in the immune tumor microenvironment (iTME) has been widely observed yet remains largely uncharacterized. Here, we developed a machine learning approach to map iTME modifiers within loci from genome-wide association studies (GWASs) for breast cancer (BrCa) incidence. A random forest model was trained on a positive set of immune-oncology (I-O) targets, and then used to assign I-O target probability scores to 1,362 candidate genes in linkage disequilibrium with 155 BrCa GWAS loci. Cluster analysis of the most probable candidates revealed two subfamilies of genes related to effector functions and adaptive immune responses, suggesting that iTME modifiers impact multiple aspects of anticancer immunity. Two of the top ranking BrCa candidates, LSP1 and TLR1, were orthogonally validated as iTME modifiers using BrCa patient biopsies and comparative mapping studies, respectively. Collectively, these data demonstrate a robust and flexible framework for functionally fine-mapping GWAS risk loci to identify translatable therapeutic targets.

A Mendelian Randomization Analysis of 55 Genetically Predicted Metabolic Traits with Breast Cancer Survival Outcomes in the Pathways Study

Previous studies suggest associations of metabolic syndromes with breast cancer prognosis, yet the evidence is mixed. In recent years, the maturation of genome-wide association study findings has led to the development of polygenic scores (PGS) for many common traits, making it feasible to use Mendelian randomization to examine associations between metabolic traits and breast cancer outcomes. In the Pathways Study of 3,902 patients and a median follow-up time of 10.5 years, we adapted a Mendelian randomization approach to calculate PGS for 55 metabolic traits and tested their associations with seven survival outcomes. Multivariable Cox proportional hazards models were used to derive HRs and 95% confidence intervals (CI) with adjustment for covariates. The highest tertile (T3) of PGS for cardiovascular disease was associated with shorter overall survival (HR = 1.34, 95% CI = 1.11-1.61) and second primary cancer-free survival (HR = 1.31, 95% CI = 1.12-1.53). PGS for hypertension (T3) was associated with shorter overall survival (HR = 1.20, 95% CI = 1.00-1.43), second primary cancer-free survival (HR = 1.24, 95% CI = 1.06-1.45), invasive disease-free survival (HR = 1.18, 95% CI = 1.01-1.38), and disease-free survival (HR = 1.21, 95% CI = 1.04-1.39). PGS for serum cystatin C levels (T3) was associated with longer disease-free survival (HR = 0.82, 95% CI = 0.71-0.95), breast event-free survival (HR = 0.74, 95% CI = 0.61-0.91), and breast cancer-specific survival (HR = 0.72, 95% CI = 0.54-0.95). The above associations were significant at a nominal P < 0.05 level but not after correcting for multiple testing (Bonferroni P < 0.0009). Our analyses revealed notable associations of PGS for cardiovascular disease, hypertension, and cystatin C levels with breast cancer survival outcomes. These findings implicate metabolic traits in breast cancer prognosis.

Significance

To our knowledge, this is the largest study of PGS for metabolic traits with breast cancer prognosis. The findings revealed significant associations of PGS for cardiovascular disease, hypertension, and cystatin C levels with several breast cancer survival outcomes. These findings implicate an underappreciated role of metabolic traits in breast cancer prognosis that would warrant further exploration.

Association between breast cancer and thyroid cancer risk: a two-sample Mendelian randomization study

Background

Breast and thyroid cancer are increasingly prevalent, but it remains unclear whether the observed associations are due to heightened medical surveillance or intrinsic etiological factors. Observational studies are vulnerable to residual confounding, reverse causality, and bias, which can compromise causal inference. In this study, we employed a two-sample Mendelian randomization (MR) analysis to establish a causal link between breast cancer and heightened thyroid cancer risk.

Methods

We obtained the single nucleotide polymorphisms (SNPs) associated with breast cancer from a genome-wide association study (GWAS) conducted by the Breast Cancer Association Consortium (BCAC). The FinnGen consortium's latest and largest accessible GWAS thyroid cancer data at the summary level. We performed four MR analyses, including the inverse-variance-weighted (IVW), weighted median, MR-Egger regression, and weighted mode, to evaluate the potential causal connection between genetically predicted breast cancer and higher risk for thyroid cancer. Sensitivity analysis, heterogeneity and pleiotropy tests were used to ensure the reliability of our findings.

Results

Our study revealed causal relationship between genetically predicted breast cancer and thyroid cancer (IVW method, odds ratio (OR) = 1.135, 95% confidence interval (CI): 1.006 to 1.279, P = 0.038). However, there was no causal association between genetically predicted triple-negative breast cancer and thyroid cancer (OR = 0.817, 95% CI: 0.610 to 1.095, P = 0.177). There was no directional pleiotropy or horizontal pleiotropy in the present study.

Conclusion

This two-sample MR study supports a causal link between ER-positive breast cancer and heightened the risk of thyroid cancer. Our analysis did not reveal a direct correlation between triple-negative breast cancer and thyroid cancer.

Development and internal-external validation of statistical and machine learning models for breast cancer prognostication: cohort study

OBJECTIVE

To develop a clinically useful model that estimates the 10 year risk of breast cancer related mortality in women (self-reported female sex) with breast cancer of any stage, comparing results from regression and machine learning approaches.

DESIGN

Population based cohort study.

SETTING

QResearch primary care database in England, with individual level linkage to the national cancer registry, Hospital Episodes Statistics, and national mortality registers.

PARTICIPANTS

141 765 women aged 20 years and older with a diagnosis of invasive breast cancer between 1 January 2000 and 31 December 2020.

MAIN OUTCOME MEASURES

Four model building strategies comprising two regression (Cox proportional hazards and competing risks regression) and two machine learning (XGBoost and an artificial neural network) approaches. Internal-external cross validation was used for model evaluation. Random effects meta-analysis that pooled estimates of discrimination and calibration metrics, calibration plots, and decision curve analysis were used to assess model performance, transportability, and clinical utility.

RESULTS

During a median 4.16 years (interquartile range 1.76-8.26) of follow-up, 21 688 breast cancer related deaths and 11 454 deaths from other causes occurred. Restricting to 10 years maximum follow-up from breast cancer diagnosis, 20 367 breast cancer related deaths occurred during a total of 688 564.81 person years. The crude breast cancer mortality rate was 295.79 per 10 000 person years (95% confidence interval 291.75 to 299.88). Predictors varied for each regression model, but both Cox and competing risks models included age at diagnosis, body mass index, smoking status, route to diagnosis, hormone receptor status, cancer stage, and grade of breast cancer. The Cox model's random effects meta-analysis pooled estimate for Harrell's C index was the highest of any model at 0.858 (95% confidence interval 0.853 to 0.864, and 95% prediction interval 0.843 to 0.873). It appeared acceptably calibrated on calibration plots. The competing risks regression model had good discrimination: pooled Harrell's C index 0.849 (0.839 to 0.859, and 0.821 to 0.876, and evidence of systematic miscalibration on summary metrics was lacking. The machine learning models had acceptable discrimination overall (Harrell's C index: XGBoost 0.821 (0.813 to 0.828, and 0.805 to 0.837); neural network 0.847 (0.835 to 0.858, and 0.816 to 0.878)), but had more complex patterns of miscalibration and more variable regional and stage specific performance. Decision curve analysis suggested that the Cox and competing risks regression models tested may have higher clinical utility than the two machine learning approaches.

CONCLUSION

In women with breast cancer of any stage, using the predictors available in this dataset, regression based methods had better and more consistent performance compared with machine learning approaches and may be worthy of further evaluation for potential clinical use, such as for stratified follow-up.

Associations of a Breast Cancer Polygenic Risk Score With Tumor Characteristics and Survival

PURPOSE

A polygenic risk score (PRS) consisting of 313 common genetic variants (PRS313) is associated with risk of breast cancer and contralateral breast cancer. This study aimed to evaluate the association of the PRS313 with clinicopathologic characteristics of, and survival following, breast cancer.

METHODS

Women with invasive breast cancer were included, 98,397 of European ancestry and 12,920 of Asian ancestry, from the Breast Cancer Association Consortium (BCAC), and 683 women from the European MINDACT trial. Associations between PRS313 and clinicopathologic characteristics, including the 70-gene signature for MINDACT, were evaluated using logistic regression analyses. Associations of PRS313 (continuous, per standard deviation) with overall survival (OS) and breast cancer-specific survival (BCSS) were evaluated with Cox regression, adjusted for clinicopathologic characteristics and treatment.

RESULTS

The PRS313 was associated with more favorable tumor characteristics. In BCAC, increasing PRS313 was associated with lower grade, hormone receptor-positive status, and smaller tumor size. In MINDACT, PRS313 was associated with a low risk 70-gene signature. In European women from BCAC, higher PRS313 was associated with better OS and BCSS: hazard ratio (HR) 0.96 (95% CI, 0.94 to 0.97) and 0.96 (95% CI, 0.94 to 0.98), but the association disappeared after adjustment for clinicopathologic characteristics (and treatment): OS HR, 1.01 (95% CI, 0.98 to 1.05) and BCSS HR, 1.02 (95% CI, 0.98 to 1.07). The results in MINDACT and Asian women from BCAC were consistent.

CONCLUSION

An increased PRS313 is associated with favorable tumor characteristics, but is not independently associated with prognosis. Thus, PRS313 has no role in the clinical management of primary breast cancer at the time of diagnosis. Nevertheless, breast cancer mortality rates will be higher for women with higher PRS313 as increasing PRS313 is associated with an increased risk of disease. This information is crucial for modeling effective stratified screening programs.

The impact of coding germline variants on contralateral breast cancer risk and survival

Evidence linking coding germline variants in breast cancer (BC)-susceptibility genes other than BRCA1, BRCA2, and CHEK2 with contralateral breast cancer (CBC) risk and breast cancer-specific survival (BCSS) is scarce. The aim of this study was to assess the association of protein-truncating variants (PTVs) and rare missense variants (MSVs) in nine known (ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53) and 25 suspected BC-susceptibility genes with CBC risk and BCSS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated with Cox regression models. Analyses included 34,401 women of European ancestry diagnosed with BC, including 676 CBCs and 3,449 BC deaths; the median follow-up was 10.9 years. Subtype analyses were based on estrogen receptor (ER) status of the first BC. Combined PTVs and pathogenic/likely pathogenic MSVs in BRCA1, BRCA2, and TP53 and PTVs in CHEK2 and PALB2 were associated with increased CBC risk [HRs (95% CIs): 2.88 (1.70-4.87), 2.31 (1.39-3.85), 8.29 (2.53-27.21), 2.25 (1.55-3.27), and 2.67 (1.33-5.35), respectively]. The strongest evidence of association with BCSS was for PTVs and pathogenic/likely pathogenic MSVs in BRCA2 (ER-positive BC) and TP53 and PTVs in CHEK2 [HRs (95% CIs): 1.53 (1.13-2.07), 2.08 (0.95-4.57), and 1.39 (1.13-1.72), respectively, after adjusting for tumor characteristics and treatment]. HRs were essentially unchanged when censoring for CBC, suggesting that these associations are not completely explained by increased CBC risk, tumor characteristics, or treatment. There was limited evidence of associations of PTVs and/or rare MSVs with CBC risk or BCSS for the 25 suspected BC genes. The CBC findings are relevant to treatment decisions, follow-up, and screening after BC diagnosis.

A Swedish Familial Genome-Wide Haplotype Analysis Identified Five Novel Breast Cancer Susceptibility Loci on 9p24.3, 11q22.3, 15q11.2, 16q24.1 and Xq21.31

Most breast cancer heritability is unexplained. We hypothesized that analysis of unrelated familial cases in a GWAS context could enable the identification of novel susceptibility loci. In order to examine the association of a haplotype with breast cancer risk, we performed a genome-wide haplotype association study using a sliding window analysis of window sizes 1-25 SNPs in 650 familial invasive breast cancer cases and 5021 controls. We identified five novel risk loci on 9p24.3 (OR 3.4; p 4.9 × 10^-11), 11q22.3 (OR 2.4; p 5.2 × 10^-9), 15q11.2 (OR 3.6; p 2.3 × 10^-8), 16q24.1 (OR 3; p 3 × 10^-8) and Xq21.31 (OR 3.3; p 1.7 × 10^-8) and confirmed three well-known loci on 10q25.13, 11q13.3, and 16q12.1. In total, 1593 significant risk haplotypes and 39 risk SNPs were distributed on the eight loci. In comparison with unselected breast cancer cases from a previous study, the OR was increased in the familial analysis in all eight loci. Analyzing familial cancer cases and controls enabled the identification of novel breast cancer susceptibility loci.

Genome-Wide Analysis of Rare Haplotypes Associated with Breast Cancer Risk

Numerous common genetic variants have been linked to breast cancer risk, but they only partially explain the total breast cancer heritability. Inference from Nordic population-based twin data indicates rare high-risk loci as the chief determinant of breast cancer risk. Here, we use haplotypes, rather than single variants, to identify rare high-risk loci for breast cancer. With computationally phased genotypes from 181,034 white British women in the UK Biobank, a genome-wide haplotype-breast cancer association analysis was conducted using sliding windows of 5 to 500 consecutive array-genotyped variants. In the discovery stage, haplotype-breast cancer associations were evaluated retrospectively in the prestudy-enrollment data including 5,487 breast cancer cases. Breast cancer hazard ratios (HR) for additive haplotypic effects were estimated using Cox regression. The replication analysis included a prospective cohort of women free of breast cancer at enrollment, of whom 3,524 later developed breast cancer. This two-stage analysis detected 13 rare loci (frequency <1%), each associated with an appreciable breast cancer-risk increase (discovery: HRs = 2.84-6.10, P < 5 × 10-8; replication: HRs = 2.08-5.61, P < 0.01). In contrast, the variants that formed these rare haplotypes individually exhibited much smaller effects. Functional annotation revealed extensive cis-regulatory DNA elements in breast cancer-related cells underlying the replicated rare haplotypes. Using phased, imputed genotypes from 30,064 cases and 25,282 controls in the DRIVE OncoArray case-control study, 6 of the 13 rare-loci associations were found generalizable (odds ratio estimates: 1.48-7.67, P < 0.05). This study demonstrates the complementary advantage of utilizing rare haplotypes to capture novel risk loci and suggests the potential for the discovery of more genetic elements contributing to cancer heritability as large data sets of germline whole-genome sequencing become available.

SIGNIFICANCE

A genome-wide two-stage haplotype analysis identifies rare haplotypes associated with breast cancer risk and suggests that the rare risk haplotypes represent long-range interactions with regulatory consequences influencing cancer risk.

Prognostic and diagnostic values of non-coding RNAs as biomarkers for breast cancer: An umbrella review and pan-cancer analysis

Background: Breast cancer (BC) is the most common cancer in women. The incidence and morbidity of BC are expected to rise rapidly. The stage at which BC is diagnosed has a significant impact on clinical outcomes. When detected early, an overall 5-year survival rate of up to 90% is possible. Although numerous studies have been conducted to assess the prognostic and diagnostic values of non-coding RNAs (ncRNAs) in breast cancer, their overall potential remains unclear. In this field of study, there are various systematic reviews and meta-analysis studies that report volumes of data. In this study, we tried to collect all these systematic reviews and meta-analysis studies in order to re-analyze their data without any restriction to breast cancer or non-coding RNA type, to make it as comprehensive as possible. Methods: Three databases, namely, PubMed, Scopus, and Web of Science (WoS), were searched to find any relevant meta-analysis studies. After thoroughly searching, the screening of titles, abstracts, and full-text and the quality of all included studies were assessed using the AMSTAR tool. All the required data including hazard ratios (HRs), sensitivity (SENS), and specificity (SPEC) were extracted for further analysis, and all analyses were carried out using Stata. Results: In the prognostic part, our initial search of three databases produced 10,548 articles, of which 58 studies were included in the current study. We assessed the correlation of non-coding RNA (ncRNA) expression with different survival outcomes in breast cancer patients: overall survival (OS) (HR = 1.521), disease-free survival (DFS) (HR = 1.33), recurrence-free survival (RFS) (HR = 1.66), progression-free survival (PFS) (HR = 1.71), metastasis-free survival (MFS) (HR = 0.90), and disease-specific survival (DSS) (HR = 0.37). After eliminating low-quality studies, the results did not change significantly. In the diagnostic part, 22 articles and 30 datasets were retrieved from 8,453 articles. The quality of all studies was determined. The bivariate and random-effects models were used to assess the diagnostic value of ncRNAs. The overall area under the curve (AUC) of ncRNAs in differentiated patients is 0.88 (SENS: 80% and SPEC: 82%). There was no difference in the potential of single and combined ncRNAs in differentiated BC patients. However, the overall potential of microRNAs (miRNAs) is higher than that of long non-coding RNAs (lncRNAs). No evidence of publication bias was found in the current study. Nine miRNAs, four lncRNAs, and five gene targets showed significant OS and RFS between normal and cancer patients based on pan-cancer data analysis, demonstrating their potential prognostic value. Conclusion: The present umbrella review showed that ncRNAs, including lncRNAs and miRNAs, can be used as prognostic and diagnostic biomarkers for breast cancer patients, regardless of the sample sources, ethnicity of patients, and subtype of breast cancer.

The association between plasma chemokines and breast cancer risk and prognosis: A mendelian randomization study

Background: Despite the potential role of several chemokines in the migration of cytotoxic immune cells to prohibit breast cancer cell proliferation, a comprehensive view of chemokines and the risk and prognosis of breast cancer is scarce, and little is known about their causal associations. Methods: With a two-sample Mendelian randomization (MR) approach, genetic instruments associated with 30 plasma chemokines were created. Their genetic associations with breast cancer and its survival by molecular subtypes were extracted from the recent genome-wide association study of 133,384 breast cancer cases and 113,789 controls, with available survival information for 96,661 patients. We further tested the associations between the polygenic risk score (PRS) for chemokines and breast cancer in the UK Biobank cohort using logistic regression models, while the association with breast cancer survival was tested using Cox regression models. In addition, the association between chemokine expression in tumors and breast cancer survival was also analyzed in the TCGA cohort using Cox regression models. Results: Plasma CCL5 was causally associated with breast cancer in the MR analysis, which was significant in the luminal and HER-2 enriched subtypes and further confirmed using PRS analysis (OR = 0.94, 95% CI = 0.89-1.00). A potential causal association with breast cancer survival was only found for plasma CCL19, especially for ER-positive patients. Although not replicated in the UK Biobank, we still found an inverse association between CCL19 expression in tumors and breast cancer overall and relapse-free survival in the TCGA cohort (HR = 0.58, 95% CI = 0.35-0.95). Conclusion: We observed an inverse association between genetic predisposition to CCL5 and breast cancer, while CCL19 was associated with breast cancer survival. These associations suggested the potential of these chemokines as tools for breast cancer prevention and treatment.

The impact of personality on the risk and survival of breast cancer: a Mendelian randomization analysis

BACKGROUND

It has long been hypothesized that personality plays a causative role in incidence and outcome of breast cancer (BC), but epidemiological evidence of association between personality and BC is inconsistent.

METHOD

We used two-sample Mendelian randomization analysis to estimate the impact of personality on the risk and survival of BC. In total, 109 single nucleotide polymorphisms (SNPs) were utilized as instruments of neuroticism from a large-scale Genome-Wide Association Studies (GWAS), and five SNPs were utilized as instruments of extraversion from Genetic of Personality Consortium and 23andMe. Genetic association with the risk and survival of overall and individual subtype BC were obtained from the Breast Cancer Association Consortium.

RESULT

Neuroticism is significantly associated with the risk of overall BC [odds ratio (OR) 1.06; 95% confidence interval (CI) 1.01-1.11; p = 0.015] and the risk of luminal A BC (OR 1.09; 95% CI 1.03-1.16; p = 0.004). Extraversion is not associated with the risk of BC. None of neuroticism or extraversion is associated with the survival of BC.

CONCLUSION

Neuroticism was associated with a modest increased risk of BC and particularly luminal A BC.

Exosomal miR-1304-3p promotes breast cancer progression in African Americans by activating cancer-associated adipocytes

Breast cancer displays disparities in mortality between African Americans and Caucasian Americans. However, the exact molecular mechanisms remain elusive. Here, we identify miR-1304-3p as the most upregulated microRNA in African American patients. Importantly, its expression significantly correlates with poor progression-free survival in African American patients. Ectopic expression of miR-1304 promotes tumor progression in vivo. Exosomal miR-1304-3p activates cancer-associated adipocytes that release lipids and enhance cancer cell growth. Moreover, we identify the anti-adipogenic gene GATA2 as the target of miR-1304-3p. Notably, a single nucleotide polymorphism (SNP) located in the miR-1304 stem-loop region shows a significant difference in frequencies of the G allele between African and Caucasian American groups, which promotes the maturation of miR-1304-3p. Therefore, our results reveal a mechanism of the disparity in breast cancer progression and suggest a potential utility of miR-1304-3p and the associated SNP as biomarkers for predicting the outcome of African American patients.

Rat Mammary carcinoma susceptibility 3 (Mcs3) pleiotropy, socioenvironmental interaction, and comparative genomics with orthologous human 15q25.1-25.2

Genome-wide association studies of breast cancer susceptibility have revealed risk-associated genetic variants and nominated candidate genes; however, the identification of causal variants and genes is often undetermined by genome-wide association studies. Comparative genomics, utilizing Rattus norvegicus strains differing in susceptibility to mammary tumor development, is a complimentary approach to identify breast cancer susceptibility genes. Mammary carcinoma susceptibility 3 (Mcs3) is a Copenhagen (COP/NHsd) allele that confers resistance to mammary carcinomas when introgressed into a mammary carcinoma susceptible Wistar Furth (WF/NHsd) genome. Here, Mcs3 was positionally mapped to a 7.2-Mb region of RNO1 spanning rs8149408 to rs107402736 (chr1:143700228-150929594, build 6.0/rn6) using WF.COP congenic strains and 7,12-dimethylbenz(a)anthracene-induced mammary carcinogenesis. Male and female WF.COP-Mcs3 rats had significantly lower body mass compared to the Wistar Furth strain. The effect on female body mass was observed only when females were raised in the absence of males indicating a socioenvironmental interaction. Furthermore, female WF.COP-Mcs3 rats, raised in the absence of males, did not develop enhanced lobuloalveolar morphologies compared to those observed in the Wistar Furth strain. Human 15q25.1-25.2 was determined to be orthologous to rat Mcs3 (chr15:80005820-82285404 and chr15:83134545-84130720, build GRCh38/hg38). A public database search of 15q25.1-25.2 revealed genome-wide significant and nominally significant associations for body mass traits and breast cancer risk. These results support the existence of a breast cancer risk-associated allele at human 15q25.1-25.2 and warrant ultrafine mapping of rat Mcs3 and human 15q25.1-25.2 to discover novel causal genes and variants.

Update Breast Cancer 2022 Part 3 - Early-Stage Breast Cancer

This review summarizes recent developments in the prevention and treatment of patients with early-stage breast cancer. The individual disease risk for different molecular subtypes was investigated in a large epidemiological study. With regard to treatment, new data are available from long-term follow-up of the Aphinity study, as well as new data on neoadjuvant therapy with atezolizumab in HER2-positive patients. Biomarkers, such as residual cancer burden, were investigated in the context of pembrolizumab therapy. A Genomic Grade Index study in elderly patients is one of a group of studies investigating the use of modern multigene tests to identify patients with an excellent prognosis in whom chemotherapy may be avoided. These and other aspects of the latest developments in the diagnosis and treatment of breast cancer are described in this review.

Integrative multi-omic analysis identifies genetically influenced DNA methylation biomarkers for breast and prostate cancers

Aberrant DNA methylation has emerged as a hallmark in several cancers and contributes to risk, oncogenesis, progression, and prognosis. In this study, we performed imputation-based and conventional methylome-wide association analyses for breast cancer (BrCa) and prostate cancer (PrCa). The imputation-based approach identified DNA methylation at cytosine-phosphate-guanine sites (CpGs) associated with BrCa and PrCa risk utilising genome-wide association summary statistics (N_BrCa = 228,951, N_PrCa = 140,254) and prebuilt methylation prediction models, while the conventional approach identified CpG associations utilising TCGA and GEO experimental methylation data (N_BrCa = 621, N_PrCa = 241). Enrichment analysis of the association results implicated 77 and 81 genetically influenced CpGs for BrCa and PrCa, respectively. Furthermore, analysis of differential gene expression around these CpGs suggests a genome-epigenome-transcriptome mechanistic relationship. Conditional analyses identified multiple independent secondary SNP associations (P_cond < 0.05) around 28 BrCa and 22 PrCa CpGs. Cross-cancer analysis identified eight common CpGs, including a strong therapeutic target in SREBF1 (17p11.2)—a key player in lipid metabolism. These findings highlight the utility of integrative analysis of multi-omic cancer data to identify robust biomarkers and understand their regulatory effects on cancer risk.

Methylome-wide association studies identify genetically-influenced CpGs associated with breast and prostate cancer risk and (epi)genome-transcriptome mechanistic relationships, with lipid metabolism genes implicated as potential therapeutic targets.

UACA locus is associated with breast cancer chemoresistance and survival

Few germline genetic variants have been robustly linked with breast cancer outcomes. We conducted trans-ethnic meta genome-wide association study (GWAS) of overall survival (OS) in 3973 breast cancer patients from the Pathways Study, one of the largest prospective breast cancer survivor cohorts. A locus spanning the UACA gene, a key regulator of tumor suppressor Par-4, was associated with OS in patients taking Par-4 dependent chemotherapies, including anthracyclines and anti-HER2 therapy, at a genome-wide significance level (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 1.27 \times 10^{ - 9}$$\end{document}P=1.27×10−9). This association was confirmed in meta-analysis across four independent prospective breast cancer cohorts (combined hazard ratio = 1.84, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 1.28 \times 10^{ - 11}$$\end{document}P=1.28×10−11). Transcriptome-wide association study revealed higher UACA gene expression was significantly associated with worse OS (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 4.68 \times 10^{ - 7}$$\end{document}P=4.68×10−7). Our study identified the UACA locus as a genetic predictor of patient outcome following treatment with anthracyclines and/or anti-HER2 therapy, which may have clinical utility in formulating appropriate treatment strategies for breast cancer patients based on their genetic makeup.

The current status of risk-stratified breast screening

Apart from high-risk scenarios such as the presence of highly penetrant genetic mutations, breast screening typically comprises mammography or tomosynthesis strategies defined by age. However, age-based screening ignores the range of breast cancer risks that individual women may possess and is antithetical to the ambitions of personalised early detection. Whilst screening mammography reduces breast cancer mortality, this is at the risk of potentially significant harms including overdiagnosis with overtreatment, and psychological morbidity associated with false positives. In risk-stratified screening, individualised risk assessment may inform screening intensity/interval, starting age, imaging modality used, or even decisions not to screen. However, clear evidence for its benefits and harms needs to be established. In this scoping review, the authors summarise the established and emerging evidence regarding several critical dependencies for successful risk-stratified breast screening: risk prediction model performance, epidemiological studies, retrospective clinical evaluations, health economic evaluations and qualitative research on feasibility and acceptability. Family history, breast density or reproductive factors are not on their own suitable for precisely estimating risk and risk prediction models increasingly incorporate combinations of demographic, clinical, genetic and imaging-related parameters. Clinical evaluations of risk-stratified screening are currently limited. Epidemiological evidence is sparse, and randomised trials only began in recent years.

A Swedish Genome-Wide Haplotype Association Analysis Identifies a Novel Breast Cancer Susceptibility Locus in 8p21.2 and Characterizes Three Loci on Chromosomes 10, 11 and 16

(1) Background: The heritability of breast cancer is partly explained but much of the genetic contribution remains to be identified. Haplotypes are often used as markers of ethnicity as they are preserved through generations. We have previously demonstrated that haplotype analysis, in addition to standard SNP association studies, could give novel and more detailed information on genetic cancer susceptibility. (2) Methods: In order to examine the association of a SNP or a haplotype to breast cancer risk, we performed a genome wide haplotype association study, using sliding window analysis of window sizes 1−25 and 50 SNPs, in 3200 Swedish breast cancer cases and 5021 controls. (3) Results: We identified a novel breast cancer susceptibility locus in 8p21.1 (OR 2.08; p 3.92 × 10−8), confirmed three known loci in 10q26.13, 11q13.3, 16q12.1-2 and further identified novel subloci within these three loci. Altogether 76 risk SNPs, 3302 risk haplotypes of window size 2−25 and 113 risk haplotypes of window size 50 at p < 5 × 10−8 on chromosomes 8, 10, 11 and 16 were identified. In the known loci haplotype analysis reached an OR of 1.48 in overall breast cancer and in familial cases OR 1.68. (4) Conclusions: Analyzing haplotypes, rather than single variants, could detect novel susceptibility loci even in small study populations but the method requires a fairly homogenous study population.

A robust method for collider bias correction in conditional genome-wide association studies

Estimated genetic associations with prognosis, or conditional on a phenotype (e.g. disease incidence), may be affected by collider bias, whereby conditioning on the phenotype induces associations between causes of the phenotype and prognosis. We propose a method, ‘Slope-Hunter’, that uses model-based clustering to identify and utilise the class of variants only affecting the phenotype to estimate the adjustment factor, assuming this class explains more variation in the phenotype than any other variant classes. Simulation studies show that our approach eliminates the bias and outperforms alternatives even in the presence of genetic correlation. In a study of fasting blood insulin levels (FI) conditional on body mass index, we eliminate paradoxical associations of the underweight loci: COBLLI; PPARG with increased FI, and reveal an association for the locus rs1421085 (FTO). In an analysis of a case-only study for breast cancer mortality, a single region remains associated with more pronounced results.

Genetic associations can be biased by conditioning on a phenotype. This study presents ‘Slope-Hunter’, a method which uses model-based clustering to correct this bias, even in the presence of genetic correlation, assuming the class of SNPs affecting only the collider explains more variation in the collider than any other class of SNPs.

Rare germline copy number variants (CNVs) and breast cancer risk

Germline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Mendelian Randomization: Concepts and Scope

Mendelian randomization (MR) is a method of studying the causal effects of modifiable exposures (i.e., potential risk factors) on health, social, and economic outcomes using genetic variants associated with the specific exposures of interest. MR provides a more robust understanding of the influence of these exposures on outcomes because germline genetic variants are randomly inherited from parents to offspring and, as a result, should not be related to potential confounding factors that influence exposure-outcome associations. The genetic variant can therefore be used as a tool to link the proposed risk factor and outcome, and to estimate this effect with less confounding and bias than conventional epidemiological approaches. We describe the scope of MR, highlighting the range of applications being made possible as genetic data sets and resources become larger and more freely available. We outline the MR approach in detail, covering concepts, assumptions, and estimation methods. We cover some common misconceptions, provide strategies for overcoming violation of assumptions, and discuss future prospects for extending the clinical applicability, methodological innovations, robustness, and generalizability of MR findings.

Identification and verification of immune-related gene prognostic signature based on ssGSEA for breast cancer

Introduction

Breast cancer (BC) is the most common cancer in women worldwide and has a high mortality rate. The fact that the tumor microenvironment affects clinical outcomes of all types of cancers underlines the involvement of various immune-related genes (IRGs). Therefore, this study aimed to establish an IRGs-based signature for the prognosis of BC patients.

Material and methods

In this study, 12 immune cell infiltrating degrees in 1,102 BC cases from The Cancer Genome Atlas (TCGA) database were assessed, and RNA-sequencing (RNA-seq) data of these samples were analyzed by single-sample gene set enrichment analysis (ssGSEA). Based on the results, high, low, and middle immune infiltrating clusters were constructed. A total of 138 overlapped differentially expressed genes (DEGs) were identified in the high and low infiltrating clusters, as well as in normal and BC samples. Univariate Cox regression and LASSO analyses were also performed. Furthermore, GSEA suggested some highly enriched pathways in the different immune infiltrating clusters, leading to a better understanding of potential mechanisms of immune infiltration in BC.

Results

Finally, 19 immune-related genes were identified that could be utilized as a potential prognostic biomarker for BC. Kaplan-Meier plot and ROC curve, univariate as well as multivariate Cox analyses were carried out, which suggested that the 19-IRG-based signature is a significant prognosis factor independent of clinical features. Based on the analysis of protein-protein interactions (PPI), the three hub genes were identified.

Conclusions

These results provide a new method to predict the prognosis and survival of BC based on the three genes' features.

Gene-Level Germline Contributions to Clinical Risk of Recurrence Scores in Black and White Patients with Breast Cancer

Continuous risk of recurrence scores (CRS) based on tumor gene expression are vital prognostic tools for breast cancer. Studies have shown that Black women (BW) have higher CRS than White women (WW). Although systemic injustices contribute substantially to breast cancer disparities, evidence of biological and germline contributions is emerging. In this study, we investigated germline genetic associations with CRS and CRS disparity using approaches modeled after transcriptome-wide association studies (TWAS). In the Carolina Breast Cancer Study, using race-specific predictive models of tumor expression from germline genetics, we performed race-stratified (N = 1,043 WW, 1,083 BW) linear regressions of three CRS (ROR-S: PAM50 subtype score; proliferation score; ROR-P: ROR-S plus proliferation score) on imputed tumor genetically regulated tumor expression (GReX). Bayesian multivariate regression and adaptive shrinkage tested GReX-prioritized genes for associations with tumor PAM50 expression and subtype to elucidate patterns of germline regulation underlying GReX-CRS associations. At FDR-adjusted P < 0.10, 7 and 1 GReX prioritized genes among WW and BW, respectively. Among WW, CRS were positively associated with MCM10, FAM64A, CCNB2, and MMP1 GReX and negatively associated with VAV3, PCSK6, and GNG11 GReX. Among BW, higher MMP1 GReX predicted lower proliferation score and ROR-P. GReX-prioritized gene and PAM50 tumor expression associations highlighted potential mechanisms for GReX-prioritized gene to CRS associations. Among patients with breast cancer, differential germline associations with CRS were found by race, underscoring the need for larger, diverse datasets in molecular studies of breast cancer. These findings also suggest possible germline trans-regulation of PAM50 tumor expression, with potential implications for CRS interpretation in clinical settings. SIGNIFICANCE: This study identifies race-specific genetic associations with breast cancer risk of recurrence scores and suggests mediation of these associations by PAM50 subtype and expression, with implications for clinical interpretation of these scores.

Whole-Genome Genotyping Using DNA Microarrays for Population Genetics

The field of population genetics has exploded in the last two decades following the sequencing of the human genome in 2001 (Green et al. Nature 526:29-31, 2015). Tools to measure genetic variation have matured significantly throughout this advancement in knowledge (Lenoir and Giannella. J Biomed Discov Collab 1:11, 2006; Marzancola et al. Methods Mol Biol 1368:161-178, 2016). In this chapter, the focus is on the laboratory methods developed to perform genome-wide genotyping utilizing DNA microarrays, which is one of the most commonly used molecular techniques to assess global genetic variation (Heller MJ, Annu Rev Biomed Eng 4:129-153, 2002). DNA microarrays allow for the interrogation of hundreds of thousands of SNPs (single nucleotide polymorphisms) at once utilizing array-based technology in conjunction with fluorescent molecular labels in a process referred to as genotyping (Marzancola et al. Methods Mol Biol 1368:161-178, 2016). Genotype data can be utilized to associate certain phenotypes in relation with specific genetic variants within a population in a process known as genome-wide association studies or GWAS (Charlesworth and Charlesworth. Heredity (Edinb) 118(1):2-9, 2017; Casillas and Barbadilla. Genetics 205(3):1003-1035, 2017). This experimental technique is a multiple-day process involving the combination of DNA extraction, amplification, fragmentation, binding, and staining (Illumina Infinium HTS Assay Protocol Guide, 2013). Many vendors supply platforms and products to assess global genetic variation using DNA microarrays (Illumina Infinium HTS Assay Protocol Guide, 2013). In this chapter, the focus is on the methods utilized to generate high-quality genotype data with the Illumina^® Infinium Global Screening Array. Although data analysis and quality control are not the focus for this chapter, they are also briefly addressed.

Germline variants and breast cancer survival in patients with distant metastases at primary breast cancer diagnosis

Breast cancer metastasis accounts for most of the deaths from breast cancer. Identification of germline variants associated with survival in aggressive types of breast cancer may inform understanding of breast cancer progression and assist treatment. In this analysis, we studied the associations between germline variants and breast cancer survival for patients with distant metastases at primary breast cancer diagnosis. We used data from the Breast Cancer Association Consortium (BCAC) including 1062 women of European ancestry with metastatic breast cancer, 606 of whom died of breast cancer. We identified two germline variants on chromosome 1, rs138569520 and rs146023652, significantly associated with breast cancer-specific survival (P = 3.19 × 10-8 and 4.42 × 10-8). In silico analysis suggested a potential regulatory effect of the variants on the nearby target genes SDE2 and H3F3A. However, the variants showed no evidence of association in a smaller replication dataset. The validation dataset was obtained from the SNPs to Risk of Metastasis (StoRM) study and included 293 patients with metastatic primary breast cancer at diagnosis. Ultimately, larger replication studies are needed to confirm the identified associations.

Association of germline genetic variants with breast cancer-specific survival in patient subgroups defined by clinic-pathological variables related to tumor biology and type of systemic treatment

BACKGROUND

Given the high heterogeneity among breast tumors, associations between common germline genetic variants and survival that may exist within specific subgroups could go undetected in an unstratified set of breast cancer patients.

METHODS

We performed genome-wide association analyses within 15 subgroups of breast cancer patients based on prognostic factors, including hormone receptors, tumor grade, age, and type of systemic treatment. Analyses were based on 91,686 female patients of European ancestry from the Breast Cancer Association Consortium, including 7531 breast cancer-specific deaths over a median follow-up of 8.1 years. Cox regression was used to assess associations of common germline variants with 15-year and 5-year breast cancer-specific survival. We assessed the probability of these associations being true positives via the Bayesian false discovery probability (BFDP < 0.15).

RESULTS

Evidence of associations with breast cancer-specific survival was observed in three patient subgroups, with variant rs5934618 in patients with grade 3 tumors (15-year-hazard ratio (HR) [95% confidence interval (CI)] 1.32 [1.20, 1.45], P = 1.4E-08, BFDP = 0.01, per G allele); variant rs4679741 in patients with ER-positive tumors treated with endocrine therapy (15-year-HR [95% CI] 1.18 [1.11, 1.26], P = 1.6E-07, BFDP = 0.09, per G allele); variants rs1106333 (15-year-HR [95% CI] 1.68 [1.39,2.03], P = 5.6E-08, BFDP = 0.12, per A allele) and rs78754389 (5-year-HR [95% CI] 1.79 [1.46,2.20], P = 1.7E-08, BFDP = 0.07, per A allele), in patients with ER-negative tumors treated with chemotherapy.

CONCLUSIONS

We found evidence of four loci associated with breast cancer-specific survival within three patient subgroups. There was limited evidence for the existence of associations in other patient subgroups. However, the power for many subgroups is limited due to the low number of events. Even so, our results suggest that the impact of common germline genetic variants on breast cancer-specific survival might be limited.

Survival Differences in Chinese Versus White Women With Breast Cancer in the United States: A SEER-Based Analysis

PURPOSE

The affect of race on breast cancer prognosis is not well understood. We compared crude and adjusted breast cancer survival rates of Chinese women versus White women in the United States.

METHODS

We conducted a cohort study of Chinese and White women with breast cancer diagnosed between 2004 to 2015 in the SEER 18 registries database. We abstracted information on age at diagnosis, tumor size, grade, lymph node status, receptor status, surgical treatment, receipt of radiotherapy and chemotherapy, and death. We compared crude breast cancer–specific mortality between the two ethnic groups. We calculated adjusted hazard ratios (HRs) in a propensity-matched design using the Cox proportional hazards model. P < .05 was considered statistically significant.

RESULTS

There were 7,553 Chinese women (1.8%) and 414,618 White women (98.2%) with stage I-IV breast cancer in the SEER database. There were small differences in demographics, nodal burden, and clinical stage between Chinese and White women. Ten-year breast cancer–specific survival was 88.8% for Chinese women and 85.6% for White women (HR, 0.73; 95% CI, 0.67 to 0.80; P < .0001). In a propensity-matched analysis among women with stage I–IIIC breast cancer, the HR was 0.71 (95% CI, 0.62 to 0.81; P < .0001). Annual mortality rates in White women exceeded those in Chinese women for the first 9 years after diagnosis.

CONCLUSION

Chinese women in the United States have superior breast cancer–specific survival compared with White women. The reason for the observed difference is not clear. Differences in demographic and tumor features between Chinese and White women with breast cancer may contribute to the disparity, as may the possibility of intrinsic biologic differences.

The role of functional polymorphisms in oxidative stress-related genes on early-stage breast cancer survival

BACKGROUND

Genetic variations in oxidative stress-related genes may alter the coded protein level and impact the pathogenesis of breast cancer.

METHODS

The current study investigated the associations of functional single nucleotide polymorphisms in the NFE2L2, HMOX1, P21, TXNRD2, and ATF3 genes with the early-stage breast cancer clinicopathological characteristics and disease-free survival, metastasis-free survival, and overall survival. A total of 202 Eastern European (Lithuanian) women with primary I-II stage breast cancer were involved. Genotyping of the single nucleotide polymorphisms was performed using TaqMan single nucleotide polymorphisms genotyping assays.

RESULTS

The CA+AA genotypes of P21 rs1801270 were significantly less frequent in patients with lymph node metastasis and larger tumor size (P=0.041 and P=0.022, respectively). The TT genotype in ATF3 rs3125289 had significantly lower risk of estrogen receptor (ER), progesterone receptor (PR) negative, and human epidermal growth factor receptor 2 (HER2) positive status (P=0.023, P=0.046, and P=0.040, respectively). In both, univariate and multivariate Cox analysis, TXNRD2 rs1139793 GG genotype vs. GA+AA was a negative prognostic factor for disease-free survival (multivariate hazard ratio (HR) 2.248; P=0.025) and overall survival (multivariate HR 2.248; P=0.029). The ATF3 rs11119982 CC genotype in the genotype model was a negative prognostic factor for disease-free survival (multivariate HR 5.878; P=0.006), metastasis-free survival (multivariate HR 4.759; P=0.018), and overall survival (multivariate HR 3.280; P=0.048).

CONCLUSION

Our findings suggest that P21 rs1801270 is associated with lymph node metastasis and larger tumor size, and ATF3 rs3125289 is associated with ER, PR, and HER2 status. Two potential, novel, early-stage breast cancer survival biomarkers, TXNRD2 rs1139793 and ATF3 rs11119982, were detected. Further investigations are needed to confirm the results of the current study.

Lack of association of CD44-rs353630 and CHI3L2-rs684559 with pancreatic ductal adenocarcinoma survival

Although pancreatic ductal adenocarcinoma (PDAC) survival is poor, there are differences in patients' response to the treatments. Detection of predictive biomarkers explaining these differences is of the utmost importance. In a recent study two genetic markers (CD44-rs353630 and CHI3L2-rs684559) were reported to be associated with survival after PDAC resection. We attempted to replicate the associations in 1856 PDAC patients (685 resected with stage I/II) from the PANcreatic Disease ReseArch (PANDoRA) consortium. We also analysed the combined effect of the two genotypes in order to compare our results with what was previously reported. Additional stratified analyses considering TNM stage of the disease and whether the patients received surgery were also performed. We observed no statistically significant associations, except for the heterozygous carriers of CD44-rs353630, who were associated with worse OS (HR = 5.01; 95% CI 1.58-15.88; p = 0.006) among patients with stage I disease. This association is in the opposite direction of those reported previously, suggesting that data obtained in such small subgroups are hardly replicable and should be considered cautiously. The two polymorphisms combined did not show any statistically significant association. Our results suggest that the effect of CD44-rs353630 and CHI3L2-rs684559 cannot be generalized to all PDAC patients.

Landscape of somatic mutations in breast cancer: new opportunities for targeted therapies in Saudi Arabian patients

Breast cancer (BCa) ranks first in incidence rate among cancers in Arab females. The association between genetic polymorphisms in tumor suppressor genes and the risk of BCa has been studied in many ethnic populations with conflicting conclusions while Arab females and Saudi Arabian studies are still lacking. We screened a cohort of Saudi BCa patients by NGS using a bespoke gene panel to clarify the genetic landscape of this population, correlating and assessing genetic findings with clinical outcomes. We identified a total of 263 mutations spanning 51 genes, including several frequently mutated. Among the genes analyzed, the highest mutation rates were found in PIK3CA (12.9%), BRCA2 (11.7%), BRCA1 (10.2%), TP53 (6.0%), MSH2 (3.8%), PMS2 (3.8%), BARD1 (3.8%), MLH1 (3.4%), CDH1 (3.0%), RAD50 (3.0%), MSH6 (3.0%), NF1 (2.6%), in addition to others. We identified multiple common recurrent variants and previously reported mutations. We also identified 46 novel variants in 22 genes that were predicted to have a pathogenic effect. Survival analysis according to the four most common mutations (BRCA1, BRCA2, TP53, and PIK3CA) showed reduced survival in BRCA1 and BRCA2-mutant patients compared to total patients. Moreover, BRCA2 was demonstrated as an independent predictor of reduced survival using independent Cox proportional hazard models.

We reveal the landscape of the mutations associated with BCa in Saudi women, highlighting the importance of routine genetic sequencing in implementation of precision therapies in KSA.

A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

Breast cancer (BC) risk for BRCA1 and BRCA2 mutation carriers varies by genetic and familial factors. About 50 common variants have been shown to modify BC risk for mutation carriers. All but three, were identified in general population studies. Other mutation carrier-specific susceptibility variants may exist but studies of mutation carriers have so far been underpowered. We conduct a novel case-only genome-wide association study comparing genotype frequencies between 60,212 general population BC cases and 13,007 cases with BRCA1 or BRCA2 mutations. We identify robust novel associations for 2 variants with BC for BRCA1 and 3 for BRCA2 mutation carriers, P < 10-8, at 5 loci, which are not associated with risk in the general population. They include rs60882887 at 11p11.2 where MADD, SP11 and EIF1, genes previously implicated in BC biology, are predicted as potential targets. These findings will contribute towards customising BC polygenic risk scores for BRCA1 and BRCA2 mutation carriers.

DNA damage response as a prognostic indicator in metastatic breast cancer via mutational analysis

Background

High tumor heterogeneity contributes to breast cancer recurrence and metastasis. However, the lack of indicators to serve as precise and reliable means of predicting breast cancer prognosis has yet to be addressed. This study aims to reveal the prognostic relevance of mutations in metastatic breast cancer (MBC) by large-scale circulating tumor DNA (ctDNA) analysis in China.

Methods

We performed ctDNA panel-captured sequencing of 958 blood samples from MBC patients including 494 hormone receptor (HR)-positive cases, 130 human epidermal growth factor receptor 2-positive cases, and 177 triple-negative breast cancer (TNBC) cases. The somatic mutations and potential targets were assessed. Progression-free survival (PFS) was analyzed using the Kaplan-Meier method.

Results

In 801 of the 958 MBC blood samples, 663 mutated genes and 5,829 nonsynonymous alterations were identified. Mutated genes of the highest frequency were tumor protein p53 (TP53, 54%), phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha (PIK3CA, 41%), estrogen receptor 1 (ESR1, 12%), myeloid/lymphoid or mixed-lineage leukemia protein 3 (MLL3, 11%), DNA (cytosine-5)-methyltransferase 3A (DNMT3A, 10%), erb-b2 receptor tyrosine kinase 2 (ERBB2, 10%), GATA binding protein 3 (GATA3, 8%), FAT atypical cadherin 1 (FAT1, 7%), phosphatase and tensin homolog (PTEN, 6%), and mitogen-activated protein kinase kinase kinase 1 (MAP3K1, 6%). Enriched mutations and driver genes in MBC varied across stages and in multiple subtypes. Moreover, TP53, ERBB2, or coexisting TP53/PIK3CA mutations in MBC were remarkably related with shorter PFS. Mutated DNA damage response (DDR) genes were significantly associated with tumor mutation burden and mutant-allele tumor heterogeneity score, as well as with worse clinical outcome.

Conclusions

Our findings indicate that the mutations of TP53, PIK3CA, ERBB2, and in particular, DDR genes, in MBC might be relevant indicators of unfavorable prognosis in MBC.

Comprehensive assessments of germline deletion structural variants reveal the association between prognostic MUC4 and CEP72 deletions and immune response gene expression in colorectal cancer patients

Background

Functional disruptions by large germline genomic structural variants in susceptible genes are known risks for cancer. We used deletion structural variants (DSVs) generated from germline whole-genome sequencing (WGS) and DSV immune-related association tumor microenvironment (TME) to predict cancer risk and prognosis.

Methods

We investigated the contribution of germline DSVs to cancer susceptibility and prognosis by silicon and causal inference models. DSVs in germline WGS data were generated from the blood samples of 192 cancer and 499 non-cancer subjects. Clinical information, including family cancer history (FCH), was obtained from the National Cheng Kung University Hospital and Taiwan Biobank. Ninety-nine colorectal cancer (CRC) patients had immune response gene expression data. We used joint calling tools and an attention-weighted model to build the cancer risk predictive model and identify DSVs in familial cancer. The survival support vector machine (survival-SVM) was used to select prognostic DSVs.

Results

We identified 671 DSVs that could predict cancer risk. The area under the curve (AUC) of the receiver operating characteristic curve (ROC) of the attention-weighted model was 0.71. The 3 most frequent DSV genes observed in cancer patients were identified as ADCY9, AURKAPS1, and RAB3GAP2 (p < 0.05). The DSVs in SGSM2 and LHFPL3 were relevant to colorectal cancer. We found a higher incidence of FCH in cancer patients than in non-cancer subjects (p < 0.05). SMYD3 and NKD2DSV genes were associated with cancer patients with FCH (p < 0.05). We identified 65 immune-associated DSV markers for assessing cancer prognosis (p < 0.05). The functional protein of MUC4 DSV gene interacted with MAGE1 expression, according to the STRING database. The causal inference model showed that deleting the CEP72 DSV gene affect the recurrence-free survival (RFS) of IFIT1 expression.

Conclusions

We established an explainable attention-weighted model for cancer risk prediction and used the survival-SVM for prognostic stratification by using germline DSVs and immune gene expression datasets. Comprehensive assessments of germline DSVs can predict the cancer risk and clinical outcome of colon cancer patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-020-00302-3.

Breast cancer risk factors and their effects on survival: a Mendelian randomisation study

BACKGROUND

Observational studies have investigated the association of risk factors with breast cancer prognosis. However, the results have been conflicting and it has been challenging to establish causality due to potential residual confounding. Using a Mendelian randomisation (MR) approach, we aimed to examine the potential causal association between breast cancer-specific survival and nine established risk factors for breast cancer: alcohol consumption, body mass index, height, physical activity, mammographic density, age at menarche or menopause, smoking, and type 2 diabetes mellitus (T2DM).

METHODS

We conducted a two-sample MR analysis on data from the Breast Cancer Association Consortium (BCAC) and risk factor summary estimates from the GWAS Catalog. The BCAC data included 86,627 female patients of European ancestry with 7054 breast cancer-specific deaths during 15 years of follow-up. Of these, 59,378 were estrogen receptor (ER)-positive and 13,692 were ER-negative breast cancer patients. For the significant association, we used sensitivity analyses and a multivariable MR model. All risk factor associations were also examined in a model adjusted by other prognostic factors.

RESULTS

Increased genetic liability to T2DM was significantly associated with worse breast cancer-specific survival (hazard ratio [HR] = 1.10, 95% confidence interval [CI] = 1.03-1.17, P value [P] = 0.003). There were no significant associations after multiple testing correction for any of the risk factors in the ER-status subtypes. For the reported significant association with T2DM, the sensitivity analyses did not show evidence for violation of the MR assumptions nor that the association was due to increased BMI. The association remained significant when adjusting by other prognostic factors.

CONCLUSIONS

This extensive MR analysis suggests that T2DM may be causally associated with worse breast cancer-specific survival and therefore that treating T2DM may improve prognosis.

Recognition of Immune Microenvironment Landscape and Immune-Related Prognostic Genes in Breast Cancer

Background

Breast cancer (BC) is the most common malignant tumor in women. The immunophenotype of tumor microenvironment (TME) has shown great therapeutic potential in tumor.

Method

The transcriptome was obtained from TCGA and GEO data. Immune infiltration was analyzed by single-sample gene set enrichment (ssGSEA). The immune feature was constructed by Cox regression analysis. In addition, the coexpression of differential expression genes (DEGs) was identified. Through enrichment analysis, the function and pathway of module genes were identified. The somatic mutations related to immune characteristics were analyzed by Maftools. By using the consistency clustering algorithm, the molecular subtypes were constructed, and the overall survival time (OS) was predicted.

Results

Immune landscape can be divided into low immune infiltration and high immune infiltration. Cox regression analysis identified seven immune cells as protective factors of BC. In the coexpression modules for DEGs of high and low immune infiltration, module 1 was related to T cells and high immune infiltration. In particular, the area under the curve (AUC) value of hub gene SASH3 was the highest, and the correlation with T cells was stronger in the high immune infiltration. Enrichment analysis found that oxidative stress, T cell aggregation, and apoptosis were observed in high immune infiltration. In addition, TP53 was identified as the most important somatic gene mutation related to immune characteristics. Importantly, we also constructed seven immune cell-based breast cancer subtypes to predict OS.

Conclusion

We evaluated the immune landscape of BC and constructed the gene characteristics related to the immune landscape. The potential of T cells and SASH3 in immunotherapy of BC was revealed, which may guide the development of new clinical treatment strategies.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma.

Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness.

To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness.

The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry.

We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10⁻⁸) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

Association of germline variation with the survival of women with BRCA1/2 pathogenic variants and breast cancer

Germline genetic variation has been suggested to influence the survival of breast cancer patients independently of tumor pathology. We have studied survival associations of genetic variants in two etiologically unique groups of breast cancer patients, the carriers of germline pathogenic variants in BRCA1 or BRCA2 genes. We found that rs57025206 was significantly associated with the overall survival, predicting higher mortality of BRCA1 carrier patients with estrogen receptor-negative breast cancer, with a hazard ratio 4.37 (95% confidence interval 3.03-6.30, P = 3.1 × 10-9). Multivariable analysis adjusted for tumor characteristics suggested that rs57025206 was an independent survival marker. In addition, our exploratory analyses suggest that the associations between genetic variants and breast cancer patient survival may depend on tumor biological subgroup and clinical patient characteristics.

What Is the Heritability of Periodontitis? A Systematic Review

The aim of this study was to systematically appraise the existing literature on the yet-unclear heritability of gingivitis and periodontitis. This review was conducted following the PRISMA guidelines. A search was conducted through the electronic databases Medline, Embase, LILACS, Cochrane Library, Open Grey, Google Scholar, and Research Gate, as complemented by a hand search, for human studies reporting measures of heritability of gingivitis and periodontitis. A total of 9,037 papers were initially identified from combined databases and 10,810 on Google Scholar. After full-text reading, 28 articles met the inclusion criteria and were carried forward to data abstraction. The reviewed data included information from >50,000 human subjects. Meta-analyses were performed by grouping studies based on design and outcome. Heritability ( H²) of periodontitis was estimated at 0.38 (95% CI, 0.34 to 0.43; I² = 12.9%) in twin studies, 0.15 (95% CI, 0.06 to 0.24; I² = 0%) in other family studies, and 0.29 (95% CI, 0.21 to 0.38; I² = 61.2%) when twin and other family studies were combined. Genome-wide association studies detected a lower heritability estimate of 0.07 (95% CI, -0.02 to 0.15) for combined definitions of periodontitis, increasing with disease severity and when the interaction with smoking was included. Furthermore, heritability tended to be lower among older age groups. Heritability for the self-reported gingivitis trait was estimated at 0.29 (95% CI, 0.22 to 0.36; I² = 37.6%), while it was not statistically significant for clinically measured gingivitis. This systematic review brings forward summary evidence to confirm that up to a third of the periodontitis variance in the population is due to genetic factors. This seems consistent across the different studied populations and increases with disease severity. In summary, up to a third of the variance of periodontitis in the population is due to genetic factors, with higher heritability for more severe disease.

Association of angiogenesis and inflammation-related gene functional polymorphisms with early-stage breast cancer prognosis

Genetic variations in inflammation- and angiogenesis-related genes may alter the coded protein level and impact the pathogenesis of breast cancer (BC). The present study investigated the association of functional single nucleotide polymorphisms (SNPs) in the VEGFA, IL-1β, IL-1α and IL-6 genes with the early-stage BC phenotype and survival. Genomic DNA and clinical data were collected for 202 adult Eastern European (Lithuanian) women with primary I-II stage BC. Genotyping of the SNPs was performed using TaqMan SNP genotyping assays. Nine VEGFA, IL-1β, IL-1α and IL-6 polymorphisms were analysed. The VEGFA and IL-6 haplotypes were inferred using Phase software. Patients were prospectively followed-up for recurrence, occurrence of metastasis and mortality until April 30, 2019. All studied genotypes were in Hardy-Weinberg equilibrium and had the same distribution as the 1,000 Genomes project Phase 3 dataset for European population. Significant associations of the studied SNPs with clinicopathologic variables were observed between IL-1α rs1800587 C allele and larger primary tumour size; IL-6 rs1800797 A allele, rs1800797 GA genotype, rs1800795 C allele, IL-6 (rs1800797-re1800795) AC diplotype and hormonal receptor-positive disease; IL-6 rs1800797 A allele and HER2 negative status. In univariate Cox survival analysis, IL-1α rs1800587 CC and IL-6 rs1800797 GG genotype carriers exhibited worse disease-free survival (DFS), metastasis-free survival (MFS) and overall survival (OS). The IL-6 rs1800795 GG genotype was associated with worse OS. IL-6 (rs1800797, rs1800795) GG/GG diplotype carriers had shorter MFS and OS. Multivariate Cox survival analysis revealed that the IL-1α rs1800587 CC genotype was an independent negative prognostic factor for DFS, MFS and OS, and the IL6 GG/GG diplotype was an independent negative prognostic factor for MFS and OS. According to the present study, functional SNPs in the IL-1α and IL-6 genes may contribute to the identification of patients at higher risk of BC recurrence, development of metastases and worse OS among early-stage patients with BC.

A framework for transcriptome-wide association studies in breast cancer in diverse study populations

BACKGROUND

The relationship between germline genetic variation and breast cancer survival is largely unknown, especially in understudied minority populations who often have poorer survival. Genome-wide association studies (GWAS) have interrogated breast cancer survival but often are underpowered due to subtype heterogeneity and clinical covariates and detect loci in non-coding regions that are difficult to interpret. Transcriptome-wide association studies (TWAS) show increased power in detecting functionally relevant loci by leveraging expression quantitative trait loci (eQTLs) from external reference panels in relevant tissues. However, ancestry- or race-specific reference panels may be needed to draw correct inference in ancestrally diverse cohorts. Such panels for breast cancer are lacking.

RESULTS

We provide a framework for TWAS for breast cancer in diverse populations, using data from the Carolina Breast Cancer Study (CBCS), a population-based cohort that oversampled black women. We perform eQTL analysis for 406 breast cancer-related genes to train race-stratified predictive models of tumor expression from germline genotypes. Using these models, we impute expression in independent data from CBCS and TCGA, accounting for sampling variability in assessing performance. These models are not applicable across race, and their predictive performance varies across tumor subtype. Within CBCS (N = 3,828), at a false discovery-adjusted significance of 0.10 and stratifying for race, we identify associations in black women near AURKA, CAPN13, PIK3CA, and SERPINB5 via TWAS that are underpowered in GWAS.

CONCLUSIONS

We show that carefully implemented and thoroughly validated TWAS is an efficient approach for understanding the genetics underpinning breast cancer outcomes in diverse populations.

A network analysis to identify mediators of germline-driven differences in breast cancer prognosis

Identifying the underlying genetic drivers of the heritability of breast cancer prognosis remains elusive. We adapt a network-based approach to handle underpowered complex datasets to provide new insights into the potential function of germline variants in breast cancer prognosis. This network-based analysis studies ~7.3 million variants in 84,457 breast cancer patients in relation to breast cancer survival and confirms the results on 12,381 independent patients. Aggregating the prognostic effects of genetic variants across multiple genes, we identify four gene modules associated with survival in estrogen receptor (ER)-negative and one in ER-positive disease. The modules show biological enrichment for cancer-related processes such as G-alpha signaling, circadian clock, angiogenesis, and Rho-GTPases in apoptosis.

Factors influencing harmonized health data collection, sharing and linkage in Denmark and Switzerland: A systematic review

Introduction

The digitalization of medicine has led to a considerable growth of heterogeneous health datasets, which could improve healthcare research if integrated into the clinical life cycle. This process requires, amongst other things, the harmonization of these datasets, which is a prerequisite to improve their quality, re-usability and interoperability. However, there is a wide range of factors that either hinder or favor the harmonized collection, sharing and linkage of health data.

Objective

This systematic review aims to identify barriers and facilitators to health data harmonization—including data sharing and linkage—by a comparative analysis of studies from Denmark and Switzerland.

Methods

Publications from PubMed, Web of Science, EMBASE and CINAHL involving cross-institutional or cross-border collection, sharing or linkage of health data from Denmark or Switzerland were searched to identify the reported barriers and facilitators to data harmonization.

Results

Of the 345 projects included, 240 were single-country and 105 were multinational studies. Regarding national projects, a Swiss study reported on average more barriers and facilitators than a Danish study. Barriers and facilitators of a technical nature were most frequently reported.

Conclusion

This systematic review gathered evidence from Denmark and Switzerland on barriers and facilitators concerning data harmonization, sharing and linkage. Barriers and facilitators were strictly interrelated with the national context where projects were carried out. Structural changes, such as legislation implemented at the national level, were mirrored in the projects. This underlines the impact of national strategies in the field of health data. Our findings also suggest that more openness and clarity in the reporting of both barriers and facilitators to data harmonization constitute a key element to promote the successful management of new projects using health data and the implementation of proper policies in this field. Our study findings are thus meaningful beyond these two countries.