Rare missense mutations in RECQL and POLG associate with inherited predisposition to breast cancer.

Several known breast cancer susceptibility genes with moderate-to-high risk alleles encode proteins involved in DNA damage response (DDR). As these explain less than half of the hereditary breast cancer cases, additional predisposing alleles are likely to be discovered. Many of the previous studies utilizing massive parallel sequencing have focused on the protein-truncating variants, and the role of rare missense mutations has remained poorly addressed. To identify novel susceptibility factors, we have systematically analyzed the data from our parallel sequencing of 796 DDR genes in 189 Northern Finnish hereditary breast cancer patients for rare missense variants, predicted as deleterious. Thirty-five variants were studied here for the disease association using Finnish breast cancer case (n = 492-2,035) and control (n = 277-1,539) cohorts. As a result, two missense variants in genes involved in DNA replication, RECQL p.I156M and POLG p.L392V, the former involving genomic and the latter mitochondrial DNA replication, showed significant association with risk of breast cancer. Rare RECQL p.I156M allele was observed in breast cancer cases only (6/1,946, 0.3%, p = 0.043), whereas POLG p.L392V was two times more frequent in breast cancer cases (53/2,238, 2.4%) compared to controls (18/1,539, 1.2%, OR = 2.1, 95% CI 1.2-3.5, p = 0.010). Based on the current genetic data, both RECQL p.I156M and POLG p.L392V represent novel breast cancer predisposing alleles.

Genetic variation in the immunosuppression pathway genes and breast cancer susceptibility: a pooled analysis of 42,510 cases and 40,577 controls from the Breast Cancer Association Consortium.

Immunosuppression plays a pivotal role in assisting tumors to evade immune destruction and promoting tumor development. We hypothesized that genetic variation in the immunosuppression pathway genes may be implicated in breast cancer tumorigenesis. We included 42,510 female breast cancer cases and 40,577 controls of European ancestry from 37 studies in the Breast Cancer Association Consortium (2015) with available genotype data for 3595 single nucleotide polymorphisms (SNPs) in 133 candidate genes. Associations between genotyped SNPs and overall breast cancer risk, and secondarily according to estrogen receptor (ER) status, were assessed using multiple logistic regression models. Gene-level associations were assessed based on principal component analysis. Gene expression analyses were conducted using RNA sequencing level 3 data from The Cancer Genome Atlas for 989 breast tumor samples and 113 matched normal tissue samples. SNP rs1905339 (A>G) in the STAT3 region was associated with an increased breast cancer risk (per allele odds ratio 1.05, 95 % confidence interval 1.03-1.08; p value = 1.4 × 10(-6)). The association did not differ significantly by ER status. On the gene level, in addition to TGFBR2 and CCND1, IL5 and GM-CSF showed the strongest associations with overall breast cancer risk (p value = 1.0 × 10(-3) and 7.0 × 10(-3), respectively). Furthermore, STAT3 and IL5 but not GM-CSF were differentially expressed between breast tumor tissue and normal tissue (p value = 2.5 × 10(-3), 4.5 × 10(-4) and 0.63, respectively). Our data provide evidence that the immunosuppression pathway genes STAT3, IL5, and GM-CSF may be novel susceptibility loci for breast cancer in women of European ancestry.

Evidence of gene-environment interactions between common breast cancer susceptibility loci and established environmental risk factors.

Various common genetic susceptibility loci have been identified for breast cancer; however, it is unclear how they combine with lifestyle/environmental risk factors to influence risk. We undertook an international collaborative study to assess gene-environment interaction for risk of breast cancer. Data from 24 studies of the Breast Cancer Association Consortium were pooled. Using up to 34,793 invasive breast cancers and 41,099 controls, we examined whether the relative risks associated with 23 single nucleotide polymorphisms were modified by 10 established environmental risk factors (age at menarche, parity, breastfeeding, body mass index, height, oral contraceptive use, menopausal hormone therapy use, alcohol consumption, cigarette smoking, physical activity) in women of European ancestry. We used logistic regression models stratified by study and adjusted for age and performed likelihood ratio tests to assess gene-environment interactions. All statistical tests were two-sided. We replicated previously reported potential interactions between LSP1-rs3817198 and parity (Pinteraction = 2.4 × 10(-6)) and between CASP8-rs17468277 and alcohol consumption (Pinteraction = 3.1 × 10(-4)). Overall, the per-allele odds ratio (95% confidence interval) for LSP1-rs3817198 was 1.08 (1.01-1.16) in nulliparous women and ranged from 1.03 (0.96-1.10) in parous women with one birth to 1.26 (1.16-1.37) in women with at least four births. For CASP8-rs17468277, the per-allele OR was 0.91 (0.85-0.98) in those with an alcohol intake of <20 g/day and 1.45 (1.14-1.85) in those who drank ≥ 20 g/day. Additionally, interaction was found between 1p11.2-rs11249433 and ever being parous (Pinteraction = 5.3 × 10(-5)), with a per-allele OR of 1.14 (1.11-1.17) in parous women and 0.98 (0.92-1.05) in nulliparous women. These data provide first strong evidence that the risk of breast cancer associated with some common genetic variants may vary with environmental risk factors.

Identification of new genetic susceptibility loci for breast cancer through consideration of gene-environment interactions.

Genes that alter disease risk only in combination with certain environmental exposures may not be detected in genetic association analysis. By using methods accounting for gene-environment (G × E) interaction, we aimed to identify novel genetic loci associated with breast cancer risk. Up to 34,475 cases and 34,786 controls of European ancestry from up to 23 studies in the Breast Cancer Association Consortium were included. Overall, 71,527 single nucleotide polymorphisms (SNPs), enriched for association with breast cancer, were tested for interaction with 10 environmental risk factors using three recently proposed hybrid methods and a joint test of association and interaction. Analyses were adjusted for age, study, population stratification, and confounding factors as applicable. Three SNPs in two independent loci showed statistically significant association: SNPs rs10483028 and rs2242714 in perfect linkage disequilibrium on chromosome 21 and rs12197388 in ARID1B on chromosome 6. While rs12197388 was identified using the joint test with parity and with age at menarche (P-values = 3 × 10(-07)), the variants on chromosome 21 q22.12, which showed interaction with adult body mass index (BMI) in 8,891 postmenopausal women, were identified by all methods applied. SNP rs10483028 was associated with breast cancer in women with a BMI below 25 kg/m(2) (OR = 1.26, 95% CI 1.15-1.38) but not in women with a BMI of 30 kg/m(2) or higher (OR = 0.89, 95% CI 0.72-1.11, P for interaction = 3.2 × 10(-05)). Our findings confirm comparable power of the recent methods for detecting G × E interaction and the utility of using G × E interaction analyses to identify new susceptibility loci.

Annexin A1 expression in a pooled breast cancer series: association with tumor subtypes and prognosis.

BACKGROUND: Annexin A1 (ANXA1) is a protein related with the carcinogenesis process and metastasis formation in many tumors. However, little is known about the prognostic value of ANXA1 in breast cancer. The purpose of this study is to evaluate the association between ANXA1 expression, BRCA1/2 germline carriership, specific tumor subtypes and survival in breast cancer patients. METHODS: Clinical-pathological information and follow-up data were collected from nine breast cancer studies from the Breast Cancer Association Consortium (BCAC) (n = 5,752) and from one study of familial breast cancer patients with BRCA1/2 mutations (n = 107). ANXA1 expression was scored based on the percentage of immunohistochemical staining in tumor cells. Survival analyses were performed using a multivariable Cox model. RESULTS: The frequency of ANXA1 positive tumors was higher in familial breast cancer patients with BRCA1/2 mutations than in BCAC patients, with 48.6 % versus 12.4 %, respectively; P <0.0001. ANXA1 was also highly expressed in BCAC tumors that were poorly differentiated, triple negative, EGFR-CK5/6 positive or had developed in patients at a young age. In the first 5 years of follow-up, patients with ANXA1 positive tumors had a worse breast cancer-specific survival (BCSS) than ANXA1 negative (HRadj = 1.35; 95 % CI = 1.05-1.73), but the association weakened after 10 years (HRadj = 1.13; 95 % CI = 0.91-1.40). ANXA1 was a significant independent predictor of survival in HER2+ patients (10-years BCSS: HRadj = 1.70; 95 % CI = 1.17-2.45). CONCLUSIONS: ANXA1 is overexpressed in familial breast cancer patients with BRCA1/2 mutations and correlated with poor prognosis features: triple negative and poorly differentiated tumors. ANXA1 might be a biomarker candidate for breast cancer survival prediction in high risk groups such as HER2+ cases.

A polymorphism in the base excision repair gene PARP2 is associated with differential prognosis by chemotherapy among postmenopausal breast cancer patients.

BACKGROUND: Personalized therapy considering clinical and genetic patient characteristics will further improve breast cancer survival. Two widely used treatments, chemotherapy and radiotherapy, can induce oxidative DNA damage and, if not repaired, cell death. Since base excision repair (BER) activity is specific for oxidative DNA damage, we hypothesized that germline genetic variation in this pathway will affect breast cancer-specific survival depending on treatment. METHODS: We assessed in 1,408 postmenopausal breast cancer patients from the German MARIE study whether cancer specific survival after adjuvant chemotherapy, anthracycline chemotherapy, and radiotherapy is modulated by 127 Single Nucleotide Polymorphisms (SNPs) in 21 BER genes. For SNPs with interaction terms showing p<0.1 (likelihood ratio test) using multivariable Cox proportional hazard analyses, replication in 6,392 patients from nine studies of the Breast Cancer Association Consortium (BCAC) was performed. RESULTS: rs878156 in PARP2 showed a differential effect by chemotherapy (p=0.093) and was replicated in BCAC studies (p=0.009; combined analysis p=0.002). Compared to non-carriers, carriers of the variant G allele (minor allele frequency=0.07) showed better survival after chemotherapy (combined allelic hazard ratio (HR)=0.75, 95% 0.53-1.07) and poorer survival when not treated with chemotherapy (HR=1.42, 95% 1.08-1.85). A similar effect modification by rs878156 was observed for anthracycline-based chemotherapy in both MARIE and BCAC, with improved survival in carriers (combined allelic HR=0.73, 95% CI 0.40-1.32). None of the SNPs showed significant differential effects by radiotherapy. CONCLUSIONS: Our data suggest for the first time that a SNP in PARP2, rs878156, may together with other genetic variants modulate cancer specific survival in breast cancer patients depending on chemotherapy. These germline SNPs could contribute towards the design of predictive tests for breast cancer patients.

The role of genetic breast cancer susceptibility variants as prognostic factors.

Recent genome-wide association studies identified 11 single nucleotide polymorphisms (SNPs) associated with breast cancer (BC) risk. We investigated these and 62 other SNPs for their prognostic relevance. Confirmed BC risk SNPs rs17468277 (CASP8), rs1982073 (TGFB1), rs2981582 (FGFR2), rs13281615 (8q24), rs3817198 (LSP1), rs889312 (MAP3K1), rs3803662 (TOX3), rs13387042 (2q35), rs4973768 (SLC4A7), rs6504950 (COX11) and rs10941679 (5p12) were genotyped for 25 853 BC patients with the available follow-up; 62 other SNPs, which have been suggested as BC risk SNPs by a GWAS or as candidate SNPs from individual studies, were genotyped for replication purposes in subsets of these patients. Cox proportional hazard models were used to test the association of these SNPs with overall survival (OS) and BC-specific survival (BCS). For the confirmed loci, we performed an accessory analysis of publicly available gene expression data and the prognosis in a different patient group. One of the 11 SNPs, rs3803662 (TOX3) and none of the 62 candidate/GWAS SNPs were associated with OS and/or BCS at P<0.01. The genotypic-specific survival for rs3803662 suggested a recessive mode of action [hazard ratio (HR) of rare homozygous carriers=1.21; 95% CI: 1.09-1.35, P=0.0002 and HR=1.29; 95% CI: 1.12-1.47, P=0.0003 for OS and BCS, respectively]. This association was seen similarly in all analyzed tumor subgroups defined by nodal status, tumor size, grade and estrogen receptor. Breast tumor expression of these genes was not associated with prognosis. With the exception of rs3803662 (TOX3), there was no evidence that any of the SNPs associated with BC susceptibility were associated with the BC survival. Survival may be influenced by a distinct set of germline variants from those influencing susceptibility.

Comprehensive analysis of human tissues reveals unique expression and localization patterns of HSF1 and HSF2.

Heat shock factors (HSFs) are the main transcriptional regulators of the evolutionarily conserved heat shock response. Beyond cell stress, several studies have demonstrated that HSFs also contribute to a vast variety of human pathologies, ranging from metabolic diseases to cancer and neurodegeneration. Despite their evident role in mitigating cellular perturbations, the functions of HSF1 and HSF2 in physiological proteostasis have remained inconclusive. Here, we analyzed a comprehensive selection of paraffin-embedded human tissue samples with immunohistochemistry. We demonstrate that both HSF1 and HSF2 display distinct expression and subcellular localization patterns in benign tissues. HSF1 localizes to the nucleus in all epithelial cell types, whereas nuclear expression of HSF2 was limited to only a few cell types, especially the spermatogonia and the urothelial umbrella cells. We observed a consistent and robust cytoplasmic expression of HSF2 across all studied smooth muscle and endothelial cells, including the smooth muscle cells surrounding the vasculature and the high endothelial venules in lymph nodes. Outstandingly, HSF2 localized specifically at cell-cell adhesion sites in a broad selection of tissue types, such as the cardiac muscle, liver, and epididymis. To the best of our knowledge, this is the first study to systematically describe the expression and localization patterns of HSF1 and HSF2 in benign human tissues. Thus, our work expands the biological landscape of these factors and creates the foundation for the identification of specific roles of HSF1 and HSF2 in normal physiological processes.

DPYSL5 is highly expressed in treatment-induced neuroendocrine prostate cancer and promotes lineage plasticity via EZH2/PRC2.

Treatment-induced neuroendocrine prostate cancer (t-NEPC) is a lethal subtype of castration-resistant prostate cancer resistant to androgen receptor (AR) inhibitors. Our study unveils that AR suppresses the neuronal development protein dihydropyrimidinase-related protein 5 (DPYSL5), providing a mechanism for neuroendocrine transformation under androgen deprivation therapy. Our unique CRPC-NEPC cohort, comprising 135 patient tumor samples, including 55 t-NEPC patient samples, exhibits a high expression of DPYSL5 in t-NEPC patient tumors. DPYSL5 correlates with neuroendocrine-related markers and inversely with AR and PSA. DPYSL5 overexpression in prostate cancer cells induces a neuron-like phenotype, enhances invasion, proliferation, and upregulates stemness and neuroendocrine-related markers. Mechanistically, DPYSL5 promotes prostate cancer cell plasticity via EZH2-mediated PRC2 activation. Depletion of DPYSL5 decreases proliferation, induces G1 phase cell cycle arrest, reverses neuroendocrine phenotype, and upregulates luminal genes. In conclusion, DPYSL5 plays a critical role in regulating prostate cancer cell plasticity, and we propose the AR/DPYSL5/EZH2/PRC2 axis as a driver of t-NEPC progression.

Matriptase-2 gene (TMPRSS6) variants associate with breast cancer survival, and reduced expression is related to triple-negative breast cancer.

Matriptase-2 (TMPRSS6) has been identified as a breast cancer risk factor. Here, we examined relationships between TMPRSS6 genetic variations and breast cancer risk and survival, and determined the gene and protein expressions in breast tumors and assessed their clinical importance. Thirteen TMPRSS6 polymorphisms were genotyped in 462 invasive breast cancer cases and 458 controls. Gene expression was analyzed from 83 tumors and protein expression from 370 tumors. We then assessed the statistical significance of associations among genotypes, clinicopathological characteristics and survival. The TMPRSS6 variant rs2543519 was associated with breast cancer risk (p = 0.032). Multivariate analysis showed that four variants had effects on survival-rs2543519 (p = 0.017), rs2235324 (p = 0.038), rs14213212 (p = 0.044) and rs733655 (p = 0.021)-which were used to create a group variable that was associated with poorer prognosis correlating with more alleles related to reduced survival (p = 0.006; risk ratio, 2.375; 95% confidence interval, 1.287-4.382). Low gene expression was related to triple-negative breast cancer (p = 0.0001), and lower protein expression was detected in undifferentiated (p = 0.019), large (p = 0.014) and ductal or lobular tumors (p = 0.036). These results confirm the association of TMRRSS6 variants with breast cancer risk and survival. Matriptase-2 levels decrease with tumor progression, and lower gene expression is seen in poor-prognosis-related triple-negative breast cancers. Our study is the first to show that matriptase-2 gene variants are related to breast cancer prognosis, supporting matriptase-2 involvement in tumor development.

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

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.

Circulating Cell-Free DNA Reflects the Clonal Evolution of Breast Cancer Tumors.

Liquid biopsy of cell-free DNA (cfDNA) is proposed as a potential method for the early detection of breast cancer (BC) metastases and following the clonal evolution of BC. Though the use of liquid biopsy is a widely discussed topic in the field, only a few studies have demonstrated such usage so far. We sequenced the DNA of matched primary tumor and metastatic sites together with the matched cfDNA samples from 18 Eastern Finnish BC patients and investigated how well cfDNA reflected the clonal evolution of BC interpreted from tumor DNA. On average, liquid biopsy detected 56.2 ± 7.2% of the somatic variants that were present either in the matched primary tumor or metastatic sites. Despite the high discordance observed between matched samples, liquid biopsy was found to reflect the clonal evolution of BC and identify novel driver variants and therapeutic targets absent from the tumor DNA. Tumor-specific somatic variants were detected in cfDNA at the time of diagnosis and 8.4 ± 2.4 months prior to detection of locoregional recurrence or distant metastases. Our results demonstrate that the sequencing of cfDNA may be used for the early detection of locoregional and distant BC metastases. Observed discordance between tumor DNA sequencing and liquid biopsy supports the parallel sequencing of cfDNA and tumor DNA to yield the most comprehensive overview for the genetic landscape of BC.

Expression profiles of small non-coding RNAs in breast cancer tumors characterize clinicopathological features and show prognostic and predictive potential.

Precision medicine approaches are required for more effective therapies for cancer. As small non-coding RNAs (sncRNAs) have recently been suggested as intriguing candidates for cancer biomarkers and have shown potential also as novel therapeutic targets, we aimed at profiling the non-miRNA sncRNAs in a large sample set to evaluate their role in invasive breast cancer (BC). We used small RNA sequencing and 195 fresh-frozen invasive BC and 22 benign breast tissue samples to identify significant associations of small nucleolar RNAs, small nuclear RNAs, and miscellaneous RNAs with the clinicopathological features and patient outcome of BC. Ninety-six and five sncRNAs significantly distinguished (Padj < 0.01) invasive local BC from benign breast tissue and metastasized BC from invasive local BC, respectively. Furthermore, 69 sncRNAs significantly associated (Padj < 0.01) with the tumor grade, hormone receptor status, subtype, and/or tumor histology. Additionally, 42 sncRNAs were observed as candidates for prognostic markers and 29 for predictive markers for radiotherapy and/or tamoxifen response (P < 0.05). We discovered the clinical relevance of sncRNAs from each studied RNA type. By introducing new sncRNA biomarker candidates for invasive BC and validating the potential of previously described ones, we have guided the way for further research that is warranted for providing novel insights into BC biology.

MCF10CA Breast Cancer Cells Utilize Hyaluronan-Coated EV-Rich Trails for Coordinated Migration.

Invasion of tumor cells through the stroma is coordinated in response to migratory cues provided by the extracellular environment. One of the most abundant molecules in the tumor microenvironment is hyaluronan, a glycosaminoglycan known to promote many hallmarks of tumor progression, including the migratory potential of tumor cells. Strikingly, hyaluronan is also often found to coat extracellular vesicles (EVs) that originate from plasma membrane tentacles of tumor cells crucial for migration, such as filopodia, and are abundant in tumor niches. Thus, it is possible that hyaluronan and hyaluronan-coated EVs have a cooperative role in promoting migration. In this work, we compared the hyaluronan synthesis, EV secretion and migratory behavior of normal and aggressive breast cell lines from MCF10 series. Single live cell confocal imaging, electron microscopy and correlative light and electron microscopy experiments revealed that migrating tumor cells form EV-rich and hyaluronan -coated trails. These trails promote the pathfinding behavior of follower cells, which is dependent on hyaluronan. Specifically, we demonstrated that plasma membrane protrusions and EVs left behind by tumor cells during migration are strongly positive for CD9. Single cell tracking demonstrated a leader-follower behavior, which was significantly decreased upon removal of pericellular hyaluronan, indicating that hyaluronan promotes the pathfinding behavior of follower cells. Chick chorioallantoic membrane assays in ovo suggest that tumor cells behave similarly in 3D conditions. This study strengthens the important role of extracellular matrix production and architecture in coordinated tumor cell movements and validates the role of EVs as important components and regulators of tumor matrix. The results suggest that tumor cells can modify the extracellular niche by forming trails, which they subsequently follow coordinatively. Future studies will clarify in more detail the orchestrated role of hyaluronan, EVs and other extracellular cues in coordinated migration and pathfinding behavior of follower cells.

M1 Macrophages Induce Protumor Inflammation in Melanoma Cells through TNFR-NF-κB Signaling.

The tumor microenvironment, with distinctive cell types and a complex extracellular matrix has a tremendous impact on cancer progression. In this study, we investigated the effects of proinflammatory (M1) and immunosuppressive (M2) macrophages on hyaluronan (HA) matrix formation and inflammatory response in melanoma cells. Proinflammatory factors secreted from M1 macrophages stimulated the formation of a thick pericellular HA matrix in melanoma cells due to upregulation of HA synthase 2 (HAS2). HAS2 silencing reversed the effect of M1 conditioned medium on pericellular HA coat formation, and interestingly, it also partly downregulated the M1 conditioned medium‒induced upregulation of inflammation-related genes (IL1ß, IL6), as did the inhibitors for TNFR and IKKγ. Gene set enrichment analysis revealed that genes related to inflammatory responses and TNF-α signaling via NF-κB are enriched in the M1 conditioned medium‒treated melanoma cells. Moreover, the expression of matrix metalloproteinase 9 and three-dimensional cell invasion were induced in these cells, whereas M2 macrophages had no effect on HA synthesis, inflammatory response, or invasion. Our results indicate that the activation of TNFR-NF-κB signaling in M1 conditioned medium‒treated cells leads to HAS2 upregulation, which associates with a protumor inflammatory and invasive phenotype of melanoma cells.

Application of radiation omics in the development of adverse outcome pathway networks: an example of radiation-induced cardiovascular disease.

BACKGROUND: Epidemiological studies have indicated that exposure of the heart to doses of ionizing radiation as low as 0.5 Gy increases the risk of cardiac morbidity and mortality with a latency period of decades. The damaging effects of radiation to myocardial and endothelial structures and functions have been confirmed radiobiologically at high dose, but much less are known at low dose. Integration of radiation biology and epidemiology data is a recommended approach to improve the radiation risk assessment process. The adverse outcome pathway (AOP) framework offers a comprehensive tool to compile and translate mechanistic information into pathological endpoints which may be relevant for risk assessment at the different levels of a biological system. Omics technologies enable the generation of large volumes of biological data at various levels of complexity, from molecular pathways to functional organisms. Given the quality and quantity of available data across levels of biology, omics data can be attractive sources of information for use within the AOP framework. It is anticipated that radiation omics studies could improve our understanding of the molecular mechanisms behind the adverse effects of radiation on the cardiovascular system. In this review, we explored the available omics studies on radiation-induced cardiovascular disease (CVD) and their applicability to the proposed AOP for CVD. RESULTS: The results of 80 omics studies published on radiation-induced CVD over the past 20 years have been discussed in the context of the AOP of CVD proposed by Chauhan et al. Most of the available omics data on radiation-induced CVD are from proteomics, transcriptomics, and metabolomics, whereas few datasets were available from epigenomics and multi-omics. The omics data presented here show great promise in providing information for several key events (KEs) of the proposed AOP of CVD, particularly oxidative stress, alterations of energy metabolism, extracellular matrix (ECM), and vascular remodeling. CONCLUSIONS: The omics data presented here shows promise to inform the various levels of the proposed AOP of CVD. However, the data highlight the urgent need of designing omics studies to address the knowledge gap concerning different radiation scenarios, time after exposure, and experimental models. This review presents the evidence to build a qualitative omics-informed AOP and provides views on the potential benefits and challenges in using omics data to assess risk-related outcomes.

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

IMPORTANCE: Rare germline genetic variants in several genes are associated with increased breast cancer (BC) risk, but their precise contributions to different disease subtypes are unclear. This information is relevant to guidelines for gene panel testing and risk prediction. OBJECTIVE: To characterize tumors associated with BC susceptibility genes in large-scale population- or hospital-based studies. DESIGN, SETTING, AND PARTICIPANTS: The multicenter, international case-control analysis of the BRIDGES study included 42 680 patients and 46 387 control participants, comprising women aged 18 to 79 years who were sampled independently of family history from 38 studies. Studies were conducted between 1991 and 2016. Sequencing and analysis took place between 2016 and 2021. EXPOSURES: Protein-truncating variants and likely pathogenic missense variants in ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53. MAIN OUTCOMES AND MEASURES: The intrinsic-like BC subtypes as defined by estrogen receptor, progesterone receptor, and ERBB2 (formerly known as HER2) status, and tumor grade; morphology; size; stage; lymph node involvement; subtype-specific odds ratios (ORs) for carrying protein-truncating variants and pathogenic missense variants in the 9 BC susceptibility genes. RESULTS: The mean (SD) ages at interview (control participants) and diagnosis (cases) were 55.1 (11.9) and 55.8 (10.6) years, respectively; all participants were of European or East Asian ethnicity. There was substantial heterogeneity in the distribution of intrinsic subtypes by gene. RAD51C, RAD51D, and BARD1 variants were associated mainly with triple-negative disease (OR, 6.19 [95% CI, 3.17-12.12]; OR, 6.19 [95% CI, 2.99-12.79]; and OR, 10.05 [95% CI, 5.27-19.19], respectively). CHEK2 variants were associated with all subtypes (with ORs ranging from 2.21-3.17) except for triple-negative disease. For ATM variants, the association was strongest for the hormone receptor (HR)+ERBB2- high-grade subtype (OR, 4.99; 95% CI, 3.68-6.76). BRCA1 was associated with increased risk of all subtypes, but the ORs varied widely, being highest for triple-negative disease (OR, 55.32; 95% CI, 40.51-75.55). BRCA2 and PALB2 variants were also associated with triple-negative disease. TP53 variants were most strongly associated with HR+ERBB2+ and HR-ERBB2+ subtypes. Tumors occurring in pathogenic variant carriers were of higher grade. For most genes and subtypes, a decline in ORs was observed with increasing age. Together, the 9 genes were associated with 27.3% of all triple-negative tumors in women 40 years or younger. CONCLUSIONS AND RELEVANCE: The results of this case-control study suggest that variants in the 9 BC risk genes differ substantially in their associated pathology but are generally associated with triple-negative and/or high-grade disease. Knowing the age and tumor subtype distributions associated with individual BC genes can potentially aid guidelines for gene panel testing, risk prediction, and variant classification and guide targeted screening strategies.

Further evidence for the contribution of the RAD51C gene in hereditary breast and ovarian cancer susceptibility.

RAD51C, a RAD51 paralogue involved in homologous recombination, is a recently established Fanconi anemia and breast cancer predisposing factor. In the initial report, RAD51C mutations were shown to confer a high risk for both breast and ovarian tumors, but most of the replication studies published so far have failed to identify any additional susceptibility alleles. Here, we report a full mutation screening of the RAD51C gene in 147 Finnish familial breast cancer cases and in 232 unselected ovarian cancer cases originating from Finland and Sweden. In addition, in order to resolve whether common RAD51C SNPs are risk factors for breast cancer, we genotyped five tagging single nucleotide polymorphisms, rs12946522, rs304270, rs304283, rs17222691, and rs28363312, all located within the gene, from 993 Finnish breast cancer cases and 871 controls for cancer associated variants. Whereas, none of the studied common SNPs associated with breast cancer susceptibility, mutation analysis revealed two clearly pathogenic alterations. RAD51C c.-13_14del27 was observed in one familial breast cancer case and c.774delT in one unselected ovarian cancer case, thus confirming that RAD51C mutations are implicated in breast and ovarian cancer predisposition, although their overall frequency seems to be low. Independent identification of the very recently reported RAD51C c.774delT mutation in yet another patient originating from Sweden suggests that it might be a recurrent mutation in that population and should be studied further. The reliable estimation of the clinical implications of carrying a defective RAD51C allele still requires the identification of additional mutation positive families.

The debatable presence of PIWI-interacting RNAs in invasive breast cancer.

Numerous factors influence breast cancer (BC) prognosis, thus complicating the prediction of outcome. By identifying biomarkers that would distinguish the cases with poorer response to therapy already at the time of diagnosis, the rate of survival could be improved. Lately, Piwi-interacting RNAs (piRNAs) have been introduced as potential cancer biomarkers, however, due to the recently raised challenges in piRNA annotations, further evaluation of piRNAs' involvement in cancer is required. We performed small RNA sequencing in 227 fresh-frozen breast tissue samples from the Eastern Finnish Kuopio Breast Cancer Project material to study the presence of piRNAs in BC and their associations with the clinicopathological features and outcome of BC patients. We observed the presence of three small RNAs annotated as piRNA database entries (DQ596932, DQ570994, and DQ571955) in our samples. The actual species of these RNAs however remain uncertain. All three small RNAs were upregulated in grade III tumors and DQ596932 additionally in estrogen receptor negative tumors. Furthermore, patients with estrogen receptor positive BC and higher DQ571955 had shorter relapse-free survival and poorer BC-specific survival, thus indicating DQ571955 as a candidate predictive marker for radiotherapy response in estrogen receptor positive BC. DQ596932 showed possible prognostic value in BC, whereas DQ570994 was identified as a candidate predictive marker for tamoxifen and chemotherapy response. These three small RNAs appear as candidate biomarkers for BC, which could after further investigation provide novel approaches for the treatment of therapy resistant BC. Overall, our results indicate that the prevalence of piRNAs in cancer is most likely not as comprehensive as has been previously thought.