DNA repair gene patterns as prognostic and predictive factors in molecular breast cancer subtypes

DNA Damage Response in Early Breast Cancer: A Phase III Cohort in the Phobos Study

We assessed the impact of DNA damage response and repair (DDR) biomarker expressions in 222 node-positive early breast cancer (BC) patients from a previous Phase III GOIM 9902 trial of adjuvant taxanes. At a median follow-up of 64 months, the original study showed no disease-free survival (DFS) or overall survival (OS) differences with the addition of docetaxel (D) to epirubicine-cyclophosphamide (EC). Immunohistochemistry was employed to assess the expression of DDR phosphoproteins (pATM, pATR, pCHK1, γH2AX, pRPA32, and pWEE1) in tumor tissue, and their association with clinical outcomes was evaluated through the Cox elastic net model. Over an extended follow-up of 234 months, we confirmed no significant differences in DFS or OS between patients treated with EC and those receiving D → EC. A DDR risk score, inversely driven by ATM and ATR expression, emerged as an independent prognostic factor for both DFS (HR = 0.41, p < 0.0001) and OS (HR = 0.61, p = 0.046). Further validation in a public adjuvant BC cohort was possible only for ATM, confirming its protective role. Overall, our findings confirm the potential role of the DDR pathway in BC prognostication and in shaping treatment strategies advocating for an integrated approach, combining molecular markers with clinical-pathological factors.

Specific and shared biological functions of PARP2 - is PARP2 really a lil' brother of PARP1?

PARP2, that belongs to the family of ADP-ribosyl transferase enzymes (ART), is a discovery of the millennium, as it was identified in 1999. Although PARP2 was described initially as a DNA repair factor, it is now evident that PARP2 partakes in the regulation or execution of multiple biological processes as inflammation, carcinogenesis and cancer progression, metabolism or oxidative stress-related diseases. Hereby, we review the involvement of PARP2 in these processes with the aim of understanding which processes are specific for PARP2, but not for other members of the ART family. A better understanding of the specific functions of PARP2 in all of these biological processes is crucial for the development of new PARP-centred selective therapies.

An image analysis approach to characterize micronuclei differences in different subtypes of breast cancer

BACKGROUND

Micronuclei (MN) have been used as screening, diagnostic and prognostic markers in multiple cancer types, including breast cancer (BC). However, the question that the MN present in all subtypes of BC are similar or different remains unanswered. We thus hypothesized that MN present in different subtypes of BC may differ in their contents which may be visible as differences in their morphologic and morphometric features. This study was thus carried out with the aim to identify the differences between MN morphometry, complexity, and texture in different subtypes of BC, such as estrogen and progesterone receptor-positive (ER+/PR+; MCF-7, T-47D), human epidermal growth factor receptor-positive (Her2 +;SKBR3) and triple-negative BC (TNBC; MDA-MB-231, MDA-MB-468) cell lines (CLs) by ImageJ software.

METHODS

For analysis of MN dimensions, MN irregularity, and texture, we used morphometry and two mathematical computer-assisted algorithms, i.e., fractal dimension (FD) and grey level co-occurrence matrix (GLCM) of ImageJ software.

RESULTS

MN area and perimeter values showed differences in the size of MN in different subtypes of BC, with the largest MN in TNBC CLs. GLCM parameters (entropy, angular second moment, inverse difference moment, contrast, and correlation) showed highly heterogenous texture in case of TNBC MN as compared to the others. FD analysis also revealed more complexity and irregularity in MN found in TNBC cells.

CONCLUSION

The study for the first time showed morphometric, architectural and texture related differences amongst MN present in different subtypes of BC. The above may reflect differences in their composition and contents. Further, these differences may point towards the distinct mechanisms involved in the formation of MN in different subtypes of BC that need to be explored further.

A comprehensive analysis of Fanconi anemia genes in Chinese patients with high-risk hereditary breast cancer

BACKGROUND

Four Fanconi anemia (FA) genes (BRCA1, BRCA2, PALB2 and RAD51C) are defined as breast cancer (BC) susceptibility genes. Other FA genes have been inconsistently associated with BC. Thus, the role of other FA genes in BC should be explored in specific populations.

METHODS

Mutations in 16 FA genes were screened with a 98-gene panel sequencing assay in a cohort of 1481 Chinese patients with high-risk hereditary BC. The association between mutations and clinicopathological characteristics as well as prognosis was analyzed. The risk of BC in carriers of FA gene mutations was assessed in the Genome Aggregation Database and the Westlake Biobank for Chinese cohort.

RESULTS

A total of 2.57% (38/1481) BC patients were identified who had 12 other FA gene germline mutations. Among them, the most frequently mutated gene was FANCA (8/1481, 0.54%). These 38 patients carried 35 distinct pathogenic/likely pathogenic variants, of which 21 were novel. We found one rare FANCB deleterious variant (c.1327-3dupT) in our cohort. There was a statistically significant difference in lymph node status between FA gene mutation carriers and non-carriers (p = 0.041). We observed a trend that mutation carriers had larger tumor sizes, lower estrogen receptor (ER) and progesterone receptor (PR) positivity rates, and lower 3.5-year invasive disease-free survival (iDFS) and distant recurrence-free survival (DRFS) rates than non-carriers (tumor size > 2 cm: 51.43% vs. 45.63%; ER positivity rates: 51.43% vs. 60.81%; PR positivity rates: 48.57% vs. 55.16%; 3.5-year iDFS rates: 58.8% vs. 66.7%; 3.5-year DRFS rates: 58.8% vs. 68.8%). The frequency of the mutations in FANCD2, FANCM and BRIP1 trended to be higher among BC cases than that in controls (p = 0.055, 0.08 and 0.08, respectively).

CONCLUSION

This study comprehensively estimated the prevalence, clinicopathological characteristics, prognosis and risk of BC associated with deleterious variants in FA genes in Chinese high-risk hereditary BC patients. It enriches our understanding of the role of FA genes with BC.

UBE2V2 promotes metastasis by regulating EMT and predicts a poor prognosis in lung adenocarcinoma

PURPOSE

As a member of the ubiquitin-conjugating enzyme (E2) family, UBE2V2 demonstrates significant tumorigenicity in many cancers. However, the relationship between UBE2V2 expression and the morbidity of lung adenocarcinoma (LUAD) is still unknown.

METHODS

We detected the mRNA and protein expression of UBE2V2 and analyzed its relationship with clinical parameters as well as survival prognosis based on bioinformatic and immunohistochemistry (IHC) in LUAD. The signaling pathway of UBE2V2 in the development of LUAD was obtained by GSEA. The TIMER database was used to investigate the association between UBE2V2 expression and the level of infiltration of different immune cells. Finally, we explored the effects of UBE2V2 knockdown on the proliferation, apoptosis, and migration of LUAD cells.

RESULTS

The results showed that UBE2V2 was a potential oncogene and might be considered an independent prognostic molecule for LUAD patients based on TCGA prediction (HR: 1.497 p = 0.012) and IHC (HR:1.864 p = 0.044). IHC showed that UBE2V2 was related to the following clinicopathological factors: gender (p = 0.043), stage (p = 0.042), and lymph node metastasis (p = 0.002). Finally, knockdown of UBE2V2 reduced the migration of LUAD cells by regulating EMT-related proteins. Knockdown of UBE2V2 induced LUAD cells to arrest in the G1 phase. Knockdown of UBE2V2 increased LUAD cell apoptosis and decreased proliferation, which might be related to the downregulation of PCNA and upregulation of P53 and ƳH2AX expression. Interestingly, UBE2V2 is negatively correlated with B cells, CD4+ T cells, macrophages, and dendritic cells.

CONCLUSION

UBE2V2 may be a valuable therapeutic target for lung cancer.

Identification of cell subpopulations associated with disease phenotypes from scRNA-seq data using PACSI

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) has revolutionized the transcriptomics field by advancing analyses from tissue-level to cell-level resolution. Despite the great advances in the development of computational methods for various steps of scRNA-seq analyses, one major bottleneck of the existing technologies remains in identifying the molecular relationship between disease phenotype and cell subpopulations, where "disease phenotype" refers to the clinical characteristics of each patient sample, and subpopulation refer to groups of single cells, which often do not correspond to clusters identified by standard single-cell clustering analysis. Here, we present PACSI, a method aimed at distinguishing cell subpopulations associated with disease phenotypes at the single-cell level.

RESULTS

PACSI takes advantage of the topological properties of biological networks to introduce a proximity-based measure that quantifies the correlation between each cell and the disease phenotype of interest. Applied to simulated data and four case studies, PACSI accurately identified cells associated with disease phenotypes such as diagnosis, prognosis, and response to immunotherapy. In addition, we demonstrated that PACSI can also be applied to spatial transcriptomics data and successfully label spots that are associated with poor survival of breast carcinoma.

CONCLUSIONS

PACSI is an efficient method to identify cell subpopulations associated with disease phenotypes. Our research shows that it has a broad range of applications in revealing mechanistic and clinical insights of diseases.

A novel role for Mms2 in the control of spontaneous mutagenesis and Pol3 abundance

Mms2 is a ubiquitin E2-variant protein with a very well-documented function in the tolerance pathway that protects both human and yeast cells from the lethal and mutagenic effects of DNA damage. Interestingly, a high expression level of human MMS2 is associated with poor survival prognosis in different cancer diseases. Here we have analyzed the physiological effects of Mms2 overproduction in yeast cells. We show that an increased level of this protein causes a spontaneous mutator effect independent of Ubc13, a cognate partner of Mms2 in the PCNA-polyubiquitinating complex responsible for the template switch. Instead, this new promutagenic role of Mms2 requires Ubc4 (E2) and two ubiquitin ligases of HECT and RING families, Rsp5 and Not4, respectively. We have established that the promutagenic activity of Mms2 is dependent on the activities of error-prone DNA polymerase ζ and Rev1. Additionally, it requires the ubiquitination of K164 in PCNA which facilitates recruitment of these translesion polymerases to the replication complex. Importantly, we have established also that the cellular abundance of Mms2 influences the cellular level of Pol3, the catalytic subunit of replicative DNA polymerase δ. Lack of Mms2 increases the Pol3 abundance, whereas in response to Mms2 overproduction the Pol3 level decreases. We hypothesize that increased levels of spontaneous mutagenesis may result from the Mms2-induced reduction in Pol3 accumulation leading to increased participation of error-prone polymerase ζ in the replication complex.

Molecular targets that sensitize cancer to radiation killing: From the bench to the bedside

Radiotherapy is a standard cytotoxic therapy against solid cancers. It uses ionizing radiation to kill tumor cells through damage to DNA, either directly or indirectly. Radioresistance is often associated with dysregulated DNA damage repair processes. Most radiosensitizers enhance radiation-mediated DNA damage and reduce the rate of DNA repair ultimately leading to accumulation of DNA damages, cell-cycle arrest, and cell death. Recently, agents targeting key signals in DNA damage response such as DNA repair pathways and cell-cycle have been developed. This new class of molecularly targeted radiosensitizing agents is being evaluated in preclinical and clinical studies to monitor their activity in potentiating radiation cytotoxicity of tumors and reducing normal tissue toxicity. The molecular pathways of DNA damage response are reviewed with a focus on the repair mechanisms, therapeutic targets under current clinical evaluation including ATM, ATR, CDK1, CDK4/6, CHK1, DNA-PKcs, PARP-1, Wee1, & MPS1/TTK and potential new targets (BUB1, and DNA LIG4) for radiation sensitization.

Exonucleases: Degrading DNA to Deal with Genome Damage, Cell Death, Inflammation and Cancer

Although DNA degradation might seem an unwanted event, it is essential in many cellular processes that are key to maintaining genomic stability and cell and organism homeostasis. The capacity to cut out nucleotides one at a time from the end of a DNA chain is present in enzymes called exonucleases. Exonuclease activity might come from enzymes with multiple other functions or specialized enzymes only dedicated to this function. Exonucleases are involved in central pathways of cell biology such as DNA replication, repair, and death, as well as tuning the immune response. Of note, malfunctioning of these enzymes is associated with immune disorders and cancer. In this review, we will dissect the impact of DNA degradation on the DNA damage response and its links with inflammation and cancer.

RecQ Helicase Somatic Alterations in Cancer

Named the “caretakers” of the genome, RecQ helicases function in several pathways to maintain genomic stability and repair DNA. This highly conserved family of enzymes consist of five different proteins in humans: RECQL1, BLM, WRN, RECQL4, and RECQL5. Biallelic germline mutations in BLM, WRN, and RECQL4 have been linked to rare cancer-predisposing syndromes. Emerging research has also implicated somatic alterations in RecQ helicases in a variety of cancers, including hematological malignancies, breast cancer, osteosarcoma, amongst others. These alterations in RecQ helicases, particularly overexpression, may lead to increased resistance of cancer cells to conventional chemotherapy. Downregulation of these proteins may allow for increased sensitivity to chemotherapy, and, therefore, may be important therapeutic targets. Here we provide a comprehensive review of our current understanding of the role of RecQ DNA helicases in cancer and discuss the potential therapeutic opportunities in targeting these helicases.

Observation of the cervical microbiome in the progression of cervical intraepithelial neoplasia

OBJECTIVE

Cervical microbial community in the cervical intraepithelial neoplasia and cervical cancer patients was analysed to study its composition, diversity and signalling pathways by high-throughput 16S rDNA sequencing,and the candidate genes associated with occurrence and progression of cervical intraepithelial neoplasia were screened out and the model was established to predict the evolution of cervical intraepithelial neoplasia malignant transformation from the cervical microbial genes aspect.

METHODS

Cervical tissues of normal, cervical intraepithelial neoplasia and cervical cancer patients without receiving any treatment were collected. The correlation between candidate genes and cervical intraepithelial neoplasia progression was initially determined by analyzing the microbial flora. Real-time fluorescence quantitative PCR was used to detect the expression of candidate genes in different cervical tissues, ROC curve and logistic regression was used to analyse and predict the risk factors related to the occurrence and progression of cervical intraepithelial neoplasia. Finally, the early warning model of cervical intraepithelial neoplasia occurrence and progression is established.

RESULTS

Cervical tissues from normal, cervical intraepithelial neoplasia and cervical cancer patients were collected for microbial community high-throughput 16S rDNA sequencing. The analysis revealed five different pathways related to cervical intraepithelial neoplasia. 10 candidate genes were selected by further bioinformatics analysis and preliminary screening. Real time PCR, ROC curve and Logistic regression analysis showed that human papillomavirus infection, TCT severity, ABCG2, TDG, PCNA were independent risk factors for cervical intraepithelial neoplasia. We used these indicators to establish a random forest model. Seven models were built through different combinations. The model 4 (ABCG2 + PCNA + TDG) was the best early warning model for the occurrence and progression of CIN.

CONCLUSIONS

A total of 5 differential pathways and 10 candidate genes related to occurrence and progression of cervical intraepithelial neoplasia were found in cervical microbial community. This study firstly identified the genes from cervical microbial community that play an important role in the occurrence and progression of cervical intraepithelial neoplasia. At the same time, the early warning model including ABCG2 + PCNA+TDG genes provided a new idea and target for clinical prediction and blocking the evolution of cervical intraepithelial neoplasia malignant transformation from the aspect of cervical microbiological related genes.

Analysis of potential genetic biomarkers and molecular mechanism of smoking-related postmenopausal osteoporosis using weighted gene co-expression network analysis and machine learning

OBJECTIVES

Smoking is a significant independent risk factor for postmenopausal osteoporosis, leading to genome variations in postmenopausal smokers. This study investigates potential biomarkers and molecular mechanisms of smoking-related postmenopausal osteoporosis (SRPO).

MATERIALS AND METHODS

The GSE13850 microarray dataset was downloaded from Gene Expression Omnibus (GEO). Gene modules associated with SRPO were identified using weighted gene co-expression network analysis (WGCNA), protein-protein interaction (PPI) analysis, and pathway and functional enrichment analyses. Feature genes were selected using two machine learning methods: support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF). The diagnostic efficiency of the selected genes was assessed by gene expression analysis and receiver operating characteristic curve.

RESULTS

Eight highly conserved modules were detected in the WGCNA network, and the genes in the module that was strongly correlated with SRPO were used for constructing the PPI network. A total of 113 hub genes were identified in the core network using topological network analysis. Enrichment analysis results showed that hub genes were closely associated with the regulation of RNA transcription and translation, ATPase activity, and immune-related signaling. Six genes (HNRNPC, PFDN2, PSMC5, RPS16, TCEB2, and UBE2V2) were selected as genetic biomarkers for SRPO by integrating the feature selection of SVM-RFE and RF.

CONCLUSION

The present study identified potential genetic biomarkers and provided a novel insight into the underlying molecular mechanism of SRPO.

UBE2V2 Positively Correlates With PD-L1 Expression and Confers Poor Patient Survival in Lung Adenocarcinoma

This research aims to explore the diagnostic and prognostic value of ubiquitin-conjugating enzyme E2 variant 2 (UBE2V2) in lung adenocarcinoma (LUAD). The expression of UBE2V2 in clinical specimens was evaluated by bioinformatics analyses and immunohistochemistry. Bioinformatics analyses relying on the The Cancer Genome Atlas (TCGA) database suggested the elevated UBE2V2 mRNA levels in LUAD in comparison to adjacent normal tissues. Gene set enrichment analyses and gene ontology term enrichment analyses further showed the involvement of UBE2V2 in the modulation of cell cycle and immune associated signaling. The correlation analyses in TCGA LUAD data set revealed the positive correlation between UBE2V2 and CCNE1, CCNE2, CCNA2, CCNB1, CCNB2, cyclin-dependent kinase (CDK)2, CDK4, and CDK1 at the mRNA level. Moreover, UBE2V2 mRNA levels were positively correlated with PD-L1 mRNA levels, the T classification, and poor survival of LUAD patients, and were negatively correlated with type II interferon response. Consistent with the results obtained from TCGA data mining, immunohistochemistry demonstrated that UBE2V2 protein levels were upregulated in LUAD in comparison to normal tissues and were positively associated with T classification. Intriguingly, a positive correlation between UBE2V2 protein levels and PD-L1 expression was also elucidated in clinical samples. Besides, UBE2V2 expression indicated a poor prognosis in LUAD patients. Our study found that UBE2V2 was identified as an independent prognostic indicator for LUAD and might serve as an alternative target for LUAD treatment.

The Fanconi anemia pathway and Breast Cancer: A comprehensive review of clinical data

The development of breast cancer depends on several risk factors, including environmental, lifestyle and genetic factors. Despite the evolution of DNA sequencing techniques and biomarker detection, the epidemiology and mechanisms of various breast cancer susceptibility genes have not been elucidated yet. Dysregulation of the DNA damage response causes genomic instability and increases the rate of mutagenesis and the risk of carcinogenesis. The Fanconi Anemia (FA) pathway is an important component of the DNA damage response and plays a critical role in the repair of DNA interstrand crosslinks and genomic stability. The FA pathway involves 22 recognized genes and specific mutations have been identified as the underlying defect in the majority of FA patients. A thorough understanding of the function and epidemiology of these genes in breast cancer is critical for the development and implementation of individualized therapies that target unique tumor profiles. Targeted therapies (PARP inhibitors) exploiting the FA pathway gene defects have been developed and have shown promising results. This narrative review summarizes the current literature on the involvement of FA genes in sporadic and familial breast cancer with a focus on clinical data derived from large cohorts.

An Aging-Related Gene Signature-Based Model for Risk Stratification and Prognosis Prediction in Lung Adenocarcinoma

Aging is an inevitable time-dependent process associated with a gradual decline in many physiological functions. Importantly, some studies have supported that aging may be involved in the development of lung adenocarcinoma (LUAD). However, no studies have described an aging-related gene (ARG)-based prognosis signature for LUAD. Accordingly, in this study, we analyzed ARG expression data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). After LASSO and Cox regression analyses, a six ARG-based signature (APOC3, EPOR, H2AFX, MXD1, PLCG2, and YWHAZ) was constructed using TCGA dataset that significantly stratified cases into high- and low-risk groups in terms of overall survival (OS). Cox regression analysis indicated that the ARG signature was an independent prognostic factor in LUAD. A nomogram based on the ARG signature and clinicopathological factors was developed in TCGA cohort and validated in the GEO dataset. Moreover, to visualize the prediction results, we established a web-based calculator yurong.shinyapps.io/ARGs_LUAD/. Calibration plots showed good consistency between the prediction of the nomogram and actual observations. Receiver operating characteristic curve and decision curve analyses indicated that the ARG nomogram had better OS prediction and clinical net benefit than the staging system. Taken together, these results established a genetic signature for LUAD based on ARGs, which may promote individualized treatment and provide promising novel molecular markers for immunotherapy.

Development and validation of a prognostic nomogram for patients with lung adenocarcinoma based on a novel 6-DNA repair-related gene signature

DNA repair-related genes (DRGs) have attracted much attention in the field of oncology. However, the prognostic role of DRGs and their biological function in lung adenocarcinoma (LUAD) remains rudimentary and inconclusive. In this study, 716 LUAD cases from two different cohorts were collected. Samples from The Cancer Genome Atlas (TCGA) were used as the training set, and data from Gene Expression Omnibus (GEO) datasets were used for validation. Using multivariate Cox analysis and LASSO regression, we constructed a DRG signature and used it, together with clinical indices, to develop a nomogram to predict 1-, 3-, and 5-year survival rates. We identified a six-DRG signature to estimate the survival of LUAD patients, which distinguished high-risk from low-risk patients with LUAD in both the training and validation cohorts. We also observed elevated levels of infiltrating CD4 memory activated T cells, resting NK cells, M0 and M1 macrophages, and activated mast cells in the high-risk group. Finally, a nomogram incorporating the signature and clinical parameters was superior to the American Joint Committee on Cancer (AJCC) staging system in predicting the survival of LUAD patients. The DRG prognostic signature and integrated nomogram could be a useful tool to predict prognosis in patients with LUAD.

Identification and external validation of a prognostic signature associated with DNA repair genes in gastric cancer

The aim of this study was to construct and validate a DNA repair-related gene signature for evaluating the overall survival (OS) of patients with gastric cancer (GC). Differentially expressed DNA repair genes between GC and normal gastric tissue samples obtained from the TCGA database were identified. Univariate Cox analysis was used to screen survival-related genes and multivariate Cox analysis was applied to construct a DNA repair-related gene signature. An integrated bioinformatics approach was performed to evaluate its diagnostic and prognostic value. The prognostic model and the expression levels of signature genes were validated using an independent external validation cohort. Two genes (CHAF1A, RMI1) were identified to establish the prognostic signature and patients ware stratified into high- and low-risk groups. Patients in high-risk group presented significant shorter survival time than patients in the low-risk group in both cohorts, which were verified by the ROC curves. Multivariate analysis showed that the prognostic signature was an independent predictor for patients with GC after adjustment for other known clinical parameters. A nomogram incorporating the signature and known clinical factors yielded better performance and net benefits in calibration plot and decision curve analyses. Further, the logistic regression classifier based on the two genes presented an excellent diagnostic power in differentiating early HCC and normal tissues with AUCs higher than 0.9. Moreover, Gene Set Enrichment Analysis revealed that diverse cancer-related pathways significantly clustered in the high-risk and low-risk groups. Immune cell infiltration analysis revealed that CHAF1A and RMI1 were correlated with several types of immune cell subtypes. A prognostic signature using CHAF1A and RMI1 was developed that effectively predicted different OS rates among patients with GC. This risk model provides new clinical evidence for the diagnostic accuracy and survival prediction of GC.

Non-Genomic Actions of Estrogens on the DNA Repair Pathways Are Associated With Chemotherapy Resistance in Breast Cancer

Estrogens have been implicated in the etiology of breast cancer for a long time. It has been stated that long-term exposure to estrogens is associated with a higher incidence of breast cancer, since estradiol (E₂) stimulates breast cell growth; however, its effect on DNA damage/repair is only starting to be investigated. Recent studies have documented that estrogens are able to modify the DNA damage response (DDR) and DNA repair mechanisms. On the other hand, it has been proposed that DDR machinery can be altered by estrogen signaling pathways, that can be related to cancer progression and chemoresistance. We have demonstrated that E₂ promotes c-Src activation and breast cancer cell motility, through a non-genomic pathway. This review discusses scientific evidence supporting this non-genomic mechanism where estrogen modifies the DNA repair pathways, and its relationship to potential causes of chemoresistance.

Multiple microarray analyses identify key genes associated with the development of Non-Small Cell Lung Cancer from Chronic Obstructive Pulmonary Disease

Introduction: Chronic obstructive pulmonary disease (COPD) is an independent risk factor of non-small cell lung cancer (NSCLC). This study aimed to analyze the key genes and potential molecular mechanisms that are involved in the development from COPD to NSCLC.

Methods: Expression profiles of COPD and NSCLC in GSE106899, GSE12472, and GSE12428 were downloaded from the Gene Expression Omnibus (GEO) database, followed by identification of the differentially expressed genes (DEGs) between COPD and NSCLC. Based on the identified DEGs, functional pathway enrichment and lung carcinogenesis-related networks analyses were performed and further visualized with Cytoscape software. Then, principal component analysis (PCA), cluster analysis, and support vector machines (SVM) verified the ability of the top modular genes to distinguish COPD from NSCLC. Additionally, the corrections between these key genes and clinical staging of NSCLC were studied using the UALCAN and HPA websites. Finally, a prognostic risk model was constructed based on multivariate Cox regression analysis. Kaplan-Meier survival curves of the top modular genes on the training and verification sets were generated.

Results: A total of 2350, 1914, and 1850 DEGs were obtained from GSE106899, GSE12472, and GSE12428 datasets, respectively. Following analysis of protein-protein interaction networks, the identified modular gene signatures containing H2AFX, MCM2, MCM3, MCM7, POLD1, and RPA1 were identified as markers for discrimination between COPD and NSCLC. The modular gene signatures were mainly enriched in the processes of DNA replication, cell cycle, mismatch repair, and others. Besides, the expression levels of these genes were significantly higher in NSCLC than in COPD, which was further verified by the immunohistochemistry. In addition, the high expression levels of H2AFX, MCM2, MCM7, and POLD1 correlate with poor prognosis of lung adenocarcinoma (LUAD). The Cox regression prognostic risk model showed the similar results and the predictive ability of this model is independent of other clinical variables.

Conclusions: This study revealed several key modules that closely relate to NSCLC with underlying disease COPD, which provide a deeper understanding of the potential mechanisms underlying the malignant development from COPD to NSCLC. This study provides valuable prognostic factors in high-risk lung cancer patients with COPD.

FAN1, a DNA Repair Nuclease, as a Modifier of Repeat Expansion Disorders

FAN1 encodes a DNA repair nuclease. Genetic deficiencies, copy number variants, and single nucleotide variants of FAN1 have been linked to karyomegalic interstitial nephritis, 15q13.3 microdeletion/microduplication syndrome (autism, schizophrenia, and epilepsy), cancer, and most recently repeat expansion diseases. For seven CAG repeat expansion diseases (Huntington's disease (HD) and certain spinocerebellar ataxias), modification of age of onset is linked to variants of specific DNA repair proteins. FAN1 variants are the strongest modifiers. Non-coding disease-delaying FAN1 variants and coding disease-hastening variants (p.R507H and p.R377W) are known, where the former may lead to increased FAN1 levels and the latter have unknown effects upon FAN1 functions. Current thoughts are that ongoing repeat expansions in disease-vulnerable tissues, as individuals age, promote disease onset. Fan1 is required to suppress against high levels of ongoing somatic CAG and CGG repeat expansions in tissues of HD and FMR1 transgenic mice respectively, in addition to participating in DNA interstrand crosslink repair. FAN1 is also a modifier of autism, schizophrenia, and epilepsy. Coupled with the association of these diseases with repeat expansions, this suggests a common mechanism, by which FAN1 modifies repeat diseases. Yet how any of the FAN1 variants modify disease is unknown. Here, we review FAN1 variants, associated clinical effects, protein structure, and the enzyme's attributed functional roles. We highlight how variants may alter its activities in DNA damage response and/or repeat instability. A thorough awareness of the FAN1 gene and FAN1 protein functions will reveal if and how it may be targeted for clinical benefit.

Increased Expression of UBE2T Predicting Poor Survival of Epithelial Ovarian Cancer: Based on Comprehensive Analysis of UBE2s, Clinical Samples, and the GEO Database

Ubiquitin-conjugating enzymes E2 (UBE2) have been reported in the microenvironment of various malignant tumors, but their correlation with ovarian cancer (OC) remains elusive. This study aimed to systematically analyze the expression patterns, prognostic value, genetic variation, and biological functions of 12 members of the UBE2 gene family in OC through the Oncomine, Gene Expression Profiling Interactive Analysis (GEPIA), Kaplan-Meier plotter, cBioPortal, and Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) databases, respectively. We found that the mRNA levels of UBE2C, UBE2N, UBE2S, and UBE2T were significantly upregulated in OC compared with those in normal ovarian tissue. In patients with serous ovarian cancer (SOC), UBE2A, UBE2B, UBE2C, UBE2G, UBE2R2, and UBE2T upregulation were associated with poor overall survival. Moreover, UBE2A, UBE2N, UBE2R2, and UBE2T upregulation and UBE2G downregulation were associated with poor progression-free survival. UBE2T exhibited a strong correlation with OC and was thus further examined. We found that UBE2T has a high diagnostic accuracy (area under the receiver operating characteristic curve = 0.969) in epithelial ovarian cancer (EOC). Immunohistochemical assays and the Gene Expression Omnibus (GEO) database revealed that UBE2T was significantly upregulated in EOC compared with that in borderline tumors, benign tumors, and normal ovarian tissues, and its high expression was associated with poor prognosis. The Cox model showed that UBE2T upregulation was an independent risk factor affecting the prognosis of EOC and SOC. Furthermore, UBE2T was associated with specific immune cells and was mainly involved in cell cycle-related events. Genomic analysis showed that TP53 and TTN mutations were associated with UBE2T expression. Gene copy number amplification and hypomethylation may be responsible for UBE2T upregulation in OC. In conclusion, UBE2 family members may play a role in the development of OC. Specifically, UBE2T could serve as a new prognostic marker and therapeutic target for this disease. (IRB Approval No. 2020PS533K).

Prediction of Overall Survival Among Female Patients With Breast Cancer Using a Prognostic Signature Based on 8 DNA Repair-Related Genes

IMPORTANCE

Breast cancer (BC), a common malignant tumor, ranks first among cancers in terms of morbidity and mortality among female patients. Currently, identifying effective prognostic models has a significant association with the prediction of the overall survival of patients with BC and guidance of clinicians in early diagnosis and treatment.

OBJECTIVES

To identify a potential DNA repair-related prognostic signature through a comprehensive evaluation and to further improve the accuracy of prediction of the overall survival of patients with BC.

DESIGN, SETTING, AND PARTICIPANTS

In this prognostic study, conducted from October 9, 2019, to February 3, 2020, the gene expression profiles and clinical data of patients with BC were collected from The Cancer Genome Atlas database. This study consisted of a training set from The Cancer Genome Atlas database and 2 validation cohorts from the Gene Expression Omnibus, which included 1096 patients with BC. A prognostic signature based on 8 DNA repair-related genes (DRGs) was developed to predict overall survival among female patients with BC.

MAIN OUTCOMES AND MEASURES

Primary screening prognostic biomarkers were analyzed using univariate Cox proportional hazards regression analysis and the least absolute shrinkage and selection operator Cox proportional hazards regression. A risk model was completely established through multivariate Cox proportional hazards regression analysis. Finally, a prognostic nomogram, combining the DRG signature and clinical characteristics of patients, was constructed. To examine the potential mechanisms of the DRGs, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed.

RESULTS

In this prognostic study based on samples from 1096 women with BC (mean [SD] age, 59.6 [13.1] years), 8 DRGs (MDC1, RPA3, MED17, DDB2, SFPQ, XRCC4, CYP19A1, and PARP3) were identified as prognostic biomarkers. The time-dependent receiver operating characteristic curve analysis suggested that the 8-gene signature had a good predictive accuracy. In the training cohort, the areas under the curve were 0.708 for 3-year survival and 0.704 for 5-year survival. In the validation cohort, the areas under the curve were 0.717 for 3-year survival and 0.772 for 5-year survival in the GSE9893 data set and 0.691 for 3-year survival and 0.718 for 5-year survival in the GSE42568 data set. This DRG signature mainly involved some regulation pathways of vascular endothelial cell proliferation.

CONCLUSIONS AND RELEVANCE

In this study, a prognostic signature using 8 DRGs was developed that successfully predicted overall survival among female patients with BC. This risk model provides new clinical evidence for the diagnostic accuracy and targeted treatment of BC.

MiRNAs Targeting Double Strand DNA Repair Pathways Lurk in Genomically Unstable Rare Fragile Sites and Determine Cancer Outcomes

Double strand break (DSB) repair mechanisms guard genome integrity and their deterioration causes genomic instability. Common and rare fragile sites (CFS and RFS, respectively) are particularly vulnerable to instability, and there is an inverse correlation between fragile site (FS) expression and DSB repair protein levels. Upon DSB repair dysfunction, genes residing at these sites are at greater risk of deregulation compared to genes located at non-FS. In this regard, it remains enigmatic why the incidence of miRNA genes at FS is higher compared to non-FS. Herein, using bioinformatics, we examined whether miRNA genes localized at FS inhibit components of DSB repair pathways and assessed their effects on cancer. We show that such miRNAs over-accumulate in RFS, and that FRAXA, which is expressed in Fragile X syndrome, is a conserved hotspot for miRNAs inhibiting DSB repair. Axes of FRAXA-residing miRNAs/DSB repair targets affect survival in a cancer type-specific manner. Moreover, copy number variations in the region encompassing these miRNA genes discriminate survival between male and female patients. Given that, thus far, only CFS have been considered relevant for carcinogenesis, our data are the first to associate RFS with cancer, through the impairment of DSB repair by the FRAXA-residing miRNAs.

A Hierarchical Machine Learning Model to Discover Gleason Grade-Specific Biomarkers in Prostate Cancer

(1) Background:One of the most common cancers that affect North American men and men worldwide is prostate cancer. The Gleason score is a pathological grading system to examine the potential aggressiveness of the disease in the prostate tissue. Advancements in computing and next-generation sequencing technology now allow us to study the genomic profiles of patients in association with their different Gleason scores more accurately and effectively. (2) Methods: In this study, we used a novel machine learning method to analyse gene expression of prostate tumours with different Gleason scores, and identify potential genetic biomarkers for each Gleason group. We obtained a publicly-available RNA-Seq dataset of a cohort of 104 prostate cancer patients from the National Center for Biotechnology Information's (NCBI) Gene Expression Omnibus (GEO) repository, and categorised patients based on their Gleason scores to create a hierarchy of disease progression. A hierarchical model with standard classifiers in different Gleason groups, also known as nodes, was developed to identify and predict nodes based on their mRNA or gene expression. In each node, patient samples were analysed via class imbalance and hybrid feature selection techniques to build the prediction model. The outcome from analysis of each node was a set of genes that could differentiate each Gleason group from the remaining groups. To validate the proposed method, the set of identified genes were used to classify a second dataset of 499 prostate cancer patients collected from cBioportal. (3) Results: The overall accuracy of applying this novel method to the first dataset was 93.3%; the method was further validated to have 87% accuracy using the second dataset. This method also identified genes that were not previously reported as potential biomarkers for specific Gleason groups. In particular, PIAS3 was identified as a potential biomarker for Gleason score 4 + 3 = 7, and UBE2V2 for Gleason score 6. (4) Insight: Previous reports show that the genes predicted by this newly proposed method strongly correlate with prostate cancer development and progression. Furthermore, pathway analysis shows that both PIAS3 and UBE2V2 share similar protein interaction pathways, the JAK/STAT signaling process.

A new risk score based on twelve hepatocellular carcinoma-specific gene expression can predict the patients' prognosis

A large panel of molecular biomarkers have been identified to predict the prognosis of hepatocellular carcinoma (HCC), yet with limited clinical application due to difficult extrapolation. We here generated a genetic risk score system comprised of 12 HCC-specific genes to better predict the prognosis of HCC patients. Four genomics profiling datasets (GSE5851, GSE28691, GSE15765 and GSE14323) were searched to seek HCC-specific genes by comparisons between cancer samples and normal liver tissues and between different subtypes of hepatic neoplasms. Univariate survival analysis screened HCC-specific genes associated with overall survival (OS) in the training dataset for next-step risk model construction. The prognostic value of the constructed HCC risk score system was then validated in the TCGA dataset. Stratified analysis indicated this scoring system showed better performance in elderly male patients with HBV infection and preoperative lower levels of creatinine, alpha-fetoprotein and platelet and higher level of albumin. Functional annotation of this risk model in high-risk patients revealed that pathways associated with cell cycle, cell migration and inflammation were significantly enriched. In summary, our constructed HCC-specific gene risk model demonstrated robustness and potentiality in predicting the prognosis of HCC patients, especially among elderly male patients with HBV infection and relatively better general conditions.

N-(2-hydroxyphenyl) acetamide (NA-2) elicits potent antitumor effect against human breast cancer cell line (MCF-7)

Breast cancer is the most dominating malignancy in females worldwide. Treatment with conventional chemotherapeutics is associated with severe adverse effects. Thus need of new compounds, with better therapeutic potential and lesser side effects still exist. In this context the present study is planned to investigate therapeutic potential of anti-inflammatory compound N-(2- hydroxyphenyl) acetamide (NA-2) against breast cancer cells (MCF-7). The compound was selected on the basis of its reported anti-inflammatory, anti-arthritic and anti-glioblastoma activities in our previous studies. MTT, Annexin-V-FITC and wound healing assays were used to analyze the effect of compound on growth inhibition, apoptosis and metastasis. While flow cytometry, RT-PCR and immunocytochemistry techniques were used to assess the effect of NA-2 on cell cycle arrest, and expression of apoptotic markers (Bax and Bcl-2) at both mRNA and protein level respectively. Data analysis revealed that NA-2 significantly inhibits growth of MCF-7 cells after 48 h treatment (IC₅₀ = 1.65 mM). NA-2 also delayed the wound healing process, arrested cell cycle at G0/G1 phase and induced apoptosis by enhancing Bax/Bcl-2 ratio. We concluded that NA-2 possesses strong anticancer activity against MCF-7 cells, which is mediated through different mechanisms, making it a useful molecule for the development of new antitumor drugs.

Mismatch repair-signature mutations activate gene enhancers across human colorectal cancer epigenomes

Commonly-mutated genes have been found for many cancers, but less is known about mutations in cis-regulatory elements. We leverage gains in tumor-specific enhancer activity, coupled with allele-biased mutation detection from H3K27ac ChIP-seq data, to pinpoint potential enhancer-activating mutations in colorectal cancer (CRC). Analysis of a genetically-diverse cohort of CRC specimens revealed that microsatellite instable (MSI) samples have a high indel rate within active enhancers. Enhancers with indels show evidence of positive selection, increased target gene expression, and a subset is highly recurrent. The indels affect short homopolymer tracts of A/T and increase affinity for FOX transcription factors. We further demonstrate that signature mismatch-repair (MMR) mutations activate enhancers using a xenograft tumor metastasis model, where mutations are induced naturally via CRISPR/Cas9 inactivation of MLH1 prior to tumor cell injection. Our results suggest that MMR signature mutations activate enhancers in CRC tumor epigenomes to provide a selective advantage.

Computational methods to identify bimodal gene expression and facilitate personalized treatment in cancer patients

Standard methods for detecting cancer-associated genes rely on comparison of sample means between cancer patients and healthy controls. While such methods have successfully identified several oncogenes and tumor suppressor genes, they neglect to account for heterogeneity within the cancer population. Genetic mutations, translocations, and amplifications are often inconsistent across tumors, and instead they often affect smaller subsets of patients. This concept gives rise to the idea of bimodally expressed genes, or genes that display two modes of expression within one population. Analysis of bimodal gene expression has been explored via a variety of techniques including test statistics and clustering. In this review, we summarize the methodologies used to quantify bimodal gene expression and address the utility of these genes in patient stratification and specialized therapeutics in breast and lung cancer. Finally we discuss the limitations and future directions for bimodal genes in the era of high-throughput sequencing and personalized medicine.

Hypomethylation of mismatch repair genes MLH1 and MSH2 is associated with chemotolerance of breast carcinoma: Clinical significance

BACKGROUND AND OBJECTIVES

The aim of the study was to understand the importance of mismatch repair genes MLH1 and MSH2 in chemotolerance and prognosis of breast carcinoma (BC).

METHODS

First, the alterations (deletion/methylation/expression) of MLH1 and MSH2 were analyzed in 45 neoadjuvant chemotherapy (NACT)-treated and 133 pretherapeutic BC samples. The chemotolerant BC cells were characterized by treating two BC cell lines MCF-7 and MDA MB 231 with two anthracycline antitumor antibiotics, doxorubicin and nogalamycin.

RESULTS

The deletion frequencies were 32% to 38% in MLH1/MSH2 genes and promoter methylation frequencies were 49% to 62% in MLH1 and 41% to 51% in MSH2 in both NACT-treated and pretherapeutic samples. The overall alteration of MLH1 and MSH2 was 58% to 71% in the samples. Reduced messenger RNA (mRNA) and protein expression were found in both the genes and it showed concordance with the molecular alterations. NACT-treated patients showed better prognosis. The chemotherapeutic drug induced increased mRNA/protein expression of the genes in BC cell lines was due to their promoter hypomethylation, as analyzed by quantitative methylation assay. This phenomenon was also evident in NACT-treated BC samples.

CONCLUSION

MLH1/MSH2 genes play a critical role in the development of BC. Hypomethylation of MLH1/MSH2 genes might be important in chemotolerance of the disease.

An integrative bioinformatics approach reveals coding and non-coding gene variants associated with gene expression profiles and outcome in breast cancer molecular subtypes

Background

Sequence variations in coding and non-coding regions of the genome can affect gene expression and signalling pathways, which in turn may influence disease outcome.

Methods

In this study, we integrated somatic mutations, gene expression and clinical data from 930 breast cancer patients included in the TCGA database. Genes associated with single mutations in molecular breast cancer subtypes were identified by the Mann-Whitney U-test and their prognostic value was evaluated by Kaplan-Meier and Cox regression analyses. Results were confirmed using gene expression profiles from the Metabric data set (n = 1988) and whole-genome sequencing data from the TCGA cohort (n = 117).

Results

The overall mutation rate in coding and non-coding regions were significantly higher in ER-negative/HER2-negative tumours (P = 2.8E–03 and P = 2.4E–07, respectively). Recurrent sequence variations were identified in non-coding regulatory regions of several cancer-associated genes, including NBPF1, PIK3CA and TP53. After multivariate regression analysis, gene signatures associated with three coding mutations (CDH1, MAP3K1 and TP53) and two non-coding variants (CRTC3 and STAG2) in cancer-related genes predicted prognosis in ER-positive/HER2-negative tumours.

Conclusions

These findings demonstrate that sequence alterations influence gene expression and oncogenic pathways, possibly affecting the outcome of breast cancer patients. Our data provide potential opportunities to identify non-coding variations with functional and clinical relevance in breast cancer.

Clonal evolution in relapsed and refractory diffuse large B-cell lymphoma is characterized by high dynamics of subclones

Little information is available about the role of certain mutations for clonal evolution and the clinical outcome during relapse in diffuse large B-cell lymphoma (DLBCL). Therefore, we analyzed formalin-fixed-paraffin-embedded tumor samples from first diagnosis, relapsed or refractory disease from 28 patients using next-generation sequencing of the exons of 104 coding genes. Non-synonymous mutations were present in 74 of the 104 genes tested. Primary tumor samples showed a median of 8 non-synonymous mutations (range: 0-24) with the used gene set. Lower numbers of non-synonymous mutations in the primary tumor were associated with a better median OS compared with higher numbers (28 versus 15 months, p=0.031). We observed three patterns of clonal evolution during relapse of disease: large global change, subclonal selection and no or minimal change possibly suggesting preprogrammed resistance. We conclude that targeted re-sequencing is a feasible and informative approach to characterize the molecular pattern of relapse and it creates novel insights into the role of dynamics of individual genes.

Bortezomib-induced miRNAs direct epigenetic silencing of locus genes and trigger apoptosis in leukemia

MicroRNAs (miRNAs) have been suggested to repress transcription via binding the 3′-untranslated regions of mRNAs. However, the involvement and details of miRNA-mediated epigenetic regulation, particularly in targeting genomic DNA and mediating epigenetic regulation, remain largely uninvestigated. In the present study, transcription factor CCAAT/enhancer binding protein delta (CEBPD) was responsive to the anticancer drug bortezomib, a clinical and highly selective drug for leukemia treatment, and contributed to bortezomib-induced cell death. Interestingly, following the identification of CEBPD-induced miRNAs, we found that miR-744, miR-3154 and miR-3162 could target CpG islands in the 5′-flanking region of the CEBPD gene. We previously demonstrated that the Yin Yang 1 (YY1)/polycomb group (PcG) protein/DNA methyltransferase (DNMT) complex is important for CCAAT/enhancer binding protein delta (CEBPD) gene inactivation; we further found that Argonaute 2 (Ago2) interacts with YY1 and binds to the CEBPD promoter. The miRNA/Ago2/YY1/PcG group protein/DNMT complex linked the inactivation of CEBPD and genes adjacent to its 5′-flanking region, including protein kinase DNA-activated catalytic polypeptide (PRKDC), minichromosome maintenance-deficient 4 (MCM4) and ubiquitin-conjugating enzyme E2 variant 2 (UBE2V2), upon bortezomib treatment. Moreover, we revealed that miRNA binding is necessary for YY1/PcG group protein/DNMT complex-mediated epigenetic gene silencing and is associated with bortezomib-induced methylation on genomic DNA. The present study successfully characterized the interactions of the miRNA/Ago2/YY1/PcG group protein/DNMT complex and provided new insights for miRNA-mediated epigenetic regulation in bortezomib-induced leukemic cell arrest and cell death.

Human RecQL4 helicase plays multifaceted roles in the genomic stability of normal and cancer cells

Human RecQ helicases that share homology with E. coli RecQ helicase play critical roles in diverse biological activities such as DNA replication, transcription, recombination and repair. Mutations in three of the five human RecQ helicases (RecQ1, WRN, BLM, RecQL4 and RecQ5) result in autosomal recessive syndromes characterized by accelerated aging symptoms and cancer incidence. Mutational inactivation of Werner (WRN) and Bloom (BLM) genes results in Werner syndrome (WS) and Bloom syndrome (BS) respectively. However, mutations in RecQL4 result in three human disorders: (I) Rothmund-Thomson syndrome (RTS), (II) RAPADILINO and (III) Baller-Gerold syndrome (BGS). Cells from WS, BS and RTS are characterized by a unique chromosomal anomaly indicating that each of the RecQ helicases performs specialized function(s) in a non-redundant manner. Elucidating the biological functions of RecQ helicases will enable us to understand not only the aging process but also to determine the cause for age-associated human diseases. Recent biochemical and molecular studies have given new insights into the multifaceted roles of RecQL4 that range from genomic stability to carcinogenesis and beyond. This review summarizes some of the existing and emerging knowledge on diverse biological functions of RecQL4 and its significance as a potential molecular target for cancer therapy.

Identification of evolutionarily conserved DNA damage response genes that alter sensitivity to cisplatin

Ovarian, head and neck, and other cancers are commonly treated with cisplatin and other DNA damaging cytotoxic agents. Altered DNA damage response (DDR) contributes to resistance of these tumors to chemotherapies, some targeted therapies, and radiation. DDR involves multiple protein complexes and signaling pathways, some of which are evolutionarily ancient and involve protein orthologs conserved from yeast to humans. To identify new regulators of cisplatin-resistance in human tumors, we integrated high throughput and curated datasets describing yeast genes that regulate sensitivity to cisplatin and/or ionizing radiation. Next, we clustered highly validated genes based on chemogenomic profiling, and then mapped orthologs of these genes in expanded genomic networks for multiple metazoans, including humans. This approach identified an enriched candidate set of genes involved in the regulation of resistance to radiation and/or cisplatin in humans. Direct functional assessment of selected candidate genes using RNA interference confirmed their activity in influencing cisplatin resistance, degree of γH2AX focus formation and ATR phosphorylation, in ovarian and head and neck cancer cell lines, suggesting impaired DDR signaling as the driving mechanism. This work enlarges the set of genes that may contribute to chemotherapy resistance and provides a new contextual resource for interpreting next generation sequencing (NGS) genomic profiling of tumors.

Base excision repair proteins couple activation-induced cytidine deaminase and endonuclease G during replication stress-induced MLL destabilization

The breakpoint cluster region of the MLL gene (MLLbcr) is frequently rearranged in therapy-related and infant acute leukaemia, but the destabilizing mechanism is poorly understood. We recently proposed that DNA replication stress results in MLLbcr cleavage via endonuclease G (EndoG) and represents the common denominator of genotoxic therapy-induced MLL destabilization. Here we performed a siRNA screen for new factors involved in replication stress-induced MLL rearrangements employing an enhanced green fluorescent protein-based reporter system. We identified 10 factors acting in line with EndoG in MLLbcr breakage or further downstream in the repair of the MLLbcr breaks, including activation-induced cytidine deaminase (AID), previously proposed to initiate MLLbcr rearrangements in an RNA transcription-dependent mechanism. Further analysis connected AID and EndoG in MLLbcr destabilization via base excision repair (BER) components. We show that replication stress-induced recruitment of EndoG to the MLLbcr and cleavage are AID/BER dependent. Notably, inhibition of the core BER factor Apurinic-apyrimidinic endonuclease 1 protects against MLLbcr cleavage in tumour and human cord blood-derived haematopoietic stem/progenitor cells, harbouring the cells of origin of leukaemia. We propose that off-target binding of AID to the MLLbcr initiates BER-mediated single-stranded DNA cleavage, which causes derailed EndoG activity ultimately resulting in leukaemogenic MLLbcr rearrangements.

Patterns of cell cycle checkpoint deregulation associated with intrinsic molecular subtypes of human breast cancer cells

Genomic instability is a hallmark of breast cancer, contributes to tumor heterogeneity, and influences chemotherapy resistance. Although Gap 2 and mitotic checkpoints are thought to prevent genomic instability, the role of these checkpoints in breast cancer is poorly understood. Here, we assess the Gap 2 and mitotic checkpoint functions of 24 breast cancer and immortalized mammary epithelial cell lines representing four of the six intrinsic molecular subtypes of breast cancer. We found that patterns of cell cycle checkpoint deregulation were associated with the intrinsic molecular subtype of breast cancer cell lines. Specifically, the luminal B and basal-like cell lines harbored two molecularly distinct Gap 2/mitosis checkpoint defects (impairment of the decatenation Gap 2 checkpoint and the spindle assembly checkpoint, respectively). All subtypes of breast cancer cell lines examined displayed aberrant DNA synthesis/Gap 2/mitosis progression and the basal-like and claudin-low cell lines exhibited increased percentages of chromatid cohesion defects. Furthermore, a decatenation Gap 2 checkpoint gene expression signature identified in the cell line panel correlated with clinical outcomes in breast cancer patients, suggesting that breast tumors may also harbor defects in decatenation Gap 2 checkpoint function. Taken together, these data imply that pharmacological targeting of signaling pathways driving these phenotypes may lead to the development of novel personalized treatment strategies for the latter two subtypes which currently lack targeted therapeutic options because of their triple negative breast cancer status.

Integrated MicroRNA-mRNA Profiling Identifies Oncostatin M as a Marker of Mesenchymal-Like ER-Negative/HER2-Negative Breast Cancer

MicroRNAs (miRNAs) simultaneously modulate different oncogenic networks, establishing a dynamic system of gene expression and pathway regulation. In this study, we analyzed global miRNA and messenger RNA (mRNA) expression profiles of 17 cell lines representing different molecular breast cancer subtypes. Spearman's rank correlation test was used to evaluate the correlation between miRNA and mRNA expression. Hierarchical clustering and pathway analysis were also performed. Publicly available gene expression profiles (n = 699) and tumor tissues (n = 80) were analyzed to assess the relevance of key miRNA-regulated pathways in human breast cancer. We identified 39 significantly deregulated miRNAs, and the integration between miRNA and mRNA data revealed the importance of immune-related pathways, particularly the Oncostatin M (OSM) signaling, associated with mesenchymal-like breast cancer cells. OSM levels correlated with genes involved in the inflammatory response, epithelial-to-mesenchymal transition (EMT), and epidermal growth factor (EGF) signaling in human estrogen receptor (ER)-negative/human epidermal growth factor receptor 2 (HER2)-negative breast cancer. Our results suggest that the deregulation of specific miRNAs may cooperatively impair immune and EMT pathways. The identification of the OSM inflammatory pathway as an important mediator of EMT in triple-negative breast cancer (TNBC) may provide a novel potential opportunity to improve therapeutic strategies.

Integrating Multi-omics Data to Dissect Mechanisms of DNA repair Dysregulation in Breast Cancer

DNA repair genes and pathways that are transcriptionally dysregulated in cancer provide the first line of evidence for the altered DNA repair status in tumours, and hence have been explored intensively as a source for biomarker discovery. The molecular mechanisms underlying DNA repair dysregulation, however, have not been systematically investigated in any cancer type. In this study, we performed a statistical analysis to dissect the roles of DNA copy number alteration (CNA), DNA methylation (DM) at gene promoter regions and the expression changes of transcription factors (TFs) in the differential expression of individual DNA repair genes in normal versus tumour breast samples. These gene-level results were summarised at pathway level to assess whether different DNA repair pathways are affected in distinct manners. Our results suggest that CNA and expression changes of TFs are major causes of DNA repair dysregulation in breast cancer, and that a subset of the identified TFs may exert global impacts on the dysregulation of multiple repair pathways. Our work hence provides novel insights into DNA repair dysregulation in breast cancer. These insights improve our understanding of the molecular basis of the DNA repair biomarkers identified thus far, and have potential to inform future biomarker discovery.

Individual and Combined Expression of DNA Damage Response Molecules PARP1, γH2AX, BRCA1, and BRCA2 Predict Shorter Survival of Soft Tissue Sarcoma Patients

DNA damage response (DDR) molecules are protective against genotoxic stresses. DDR molecules are also involved in the survival of cancer cells in patients undergoing anti-cancer therapies. Therefore, DDR molecules are potential markers of cancer progression in addition to being potential therapeutic targets. In this study, we evaluated the immunohistochemical expression of PARP1, γH2AX, BRCA1, and BRCA2 and their prognostic significance in 112 cases of soft tissue sarcoma (STS). The expression of PARP1, γH2AX, BRCA1, and BRCA2 were significantly associated with each other and were associated with higher tumor stage and presence of distant metastasis. The expression of PARP1, γH2AX, and BRCA2 were significantly associated with shorter disease-specific survival (DSS) and event-free survival (EFS) by univariate analysis. BRCA1 expression was associated with shorter DSS. Multivariate analysis revealed the expression of PARP1 and γH2AX to be independent indicators of poor prognosis of DSS and EFS. BRCA2 expression was an independent indicator of poor prognosis of DSS. In addition, the combined expressional patterns of PARP1, γH2AX, BRCA1, and BRCA2 (CSddrm) were independent prognostic predictors of DSS (P < 0.001) and EFS (P = 0.016). The ten-year DSS rate of the CSddrm-low, CSddrm-intermediate, and CSddrm-high subgroups were 81%, 26%, and 0%, respectively. In conclusion, this study demonstrates that the individual and combined expression patterns of the DDR molecules PARP1, γH2AX, BRCA1, and BRCA2 could be predictive of the prognosis of STS patients and suggests that controlling the activity of these DDR molecules could be employed in new therapeutic stratagems for the treatment of STS.

Deciphering and Targeting Oncogenic Mutations and Pathways in Breast Cancer

: Advances in DNA and RNA sequencing revealed substantially greater genomic complexity in breast cancer than simple models of a few driver mutations would suggest. Only very few, recurrent mutations or copy-number variations in cancer-causing genes have been identified. The two most common alterations in breast cancer are TP53 (affecting the majority of triple-negative breast cancers) and PIK3CA (affecting almost half of estrogen receptor-positive cancers) mutations, followed by a long tail of individually rare mutations affecting <1%-20% of cases. Each cancer harbors from a few dozen to a few hundred potentially high-functional impact somatic variants, along with a much larger number of potentially high-functional impact germline variants. It is likely that it is the combined effect of all genomic variations that drives the clinical behavior of a given cancer. Furthermore, entirely new classes of oncogenic events are being discovered in the noncoding areas of the genome and in noncoding RNA species driven by errors in RNA editing. In light of this complexity, it is not unexpected that, with the exception of HER2 amplification, no robust molecular predictors of benefit from targeted therapies have been identified. In this review, we summarize the current genomic portrait of breast cancer, focusing on genetic aberrations that are actively being targeted with investigational drugs.

IMPLICATIONS FOR PRACTICE

Next-generation sequencing is now widely available in the clinic, but interpretation of the results is challenging, and its impact on treatment selection is often limited. This work provides an overview of frequently encountered molecular abnormalities in breast cancer and discusses their potential therapeutic implications. This review emphasizes the importance of administering investigational targeted therapies, or off-label use of approved targeted drugs, in the context of a formal clinical trial or registry programs to facilitate learning about the clinical utility of tumor target profiling.

Differential Potential of Pharmacological PARP Inhibitors for Inhibiting Cell Proliferation and Inducing Apoptosis in Human Breast Cancer Cells

BRCA1/2-mutant cells are hypersensitive to inactivation of poly(ADP-ribose) polymerase 1 (PARP-1). We recently showed that inhibition of PARP-1 by NU1025 is strongly cytotoxic for BRCA1-positive BT-20 cells, but not BRCA1-deficient SKBr-3 cells. These results raised the possibility that other PARP-1 inhibitors, particularly those tested in clinical trials, may be more efficacious against BRCA1-deficient SKBr-3 breast cancer cells than NU1025. Thus, in the presented study the cytotoxicity of four PARP inhibitors under clinical evaluation (olaparib, rucaparib, iniparib and AZD2461) was examined and compared to that of NU1025. The sensitivity of breast cancer cells to the PARP-1 inhibition strongly varied. Remarkably, BRCA-1-deficient SKBr-3 cells were almost completely insensitive to NU1025, olaparib and rucaparib, whereas BRCA1-expressing BT-20 cells were strongly affected by NU1025 even at low doses. In contrast, iniparib and AZD2461 were cytotoxic for both BT-20 and SKBr-3 cells. Of the four tested PARP-1 inhibitors only AZD2461 strongly affected cell cycle progression. Interestingly, the anti-proliferative and pro-apoptotic potential of the tested PARP-1 inhibitors clearly correlated with their capacity to damage DNA. Further analyses revealed that proteomic signatures of the two studied breast cancer cell lines strongly differ, and a set of 197 proteins was differentially expressed in NU1025-treated BT-20 cancer cells. These results indicate that BT-20 cells may harbor an unknown defect in DNA repair pathway(s) rendering them sensitive to PARP-1 inhibition. They also imply that therapeutic applicability of PARP-1 inhibitors is not limited to BRCA mutation carriers but can be extended to patients harboring deficiencies in other components of the pathway(s).

Genomic Integration of High-Risk HPV Alters Gene Expression in Oropharyngeal Squamous Cell Carcinoma

High-risk HPV (hrHPV) is the leading etiologic factor in oropharyngeal cancer. HPV-positive oropharyngeal tumors generally respond well to therapy, with complete recovery in approximately 80% of patients. However, it remains unclear why some patients are nonresponsive to treatment, with 20% of patients recurring within 5 years. In this study, viral factors were examined for possible clues to differences in tumor behavior. Oropharynx tumors that responded well to therapy were compared with those that persisted and recurred. Viral oncogene alternate transcripts were assessed, and cellular sites of viral integration were mapped and sequenced. Effects of integration on gene expression were assessed by transcript analysis at the integration sites. All of the tumors demonstrated active viral oncogenesis, indicated by expression of HPV E6 and E7 oncogenes and alternate E6 splicing. In the responsive tumors, HPV integration occurred exclusively in intergenic chromosome regions, except for one tumor with viral integration into TP63. Each recurrent tumor exhibited complex HPV integration patterns into cancer-associated genes, including TNFRSF13B, SCN2A, SH2B1, UBE2V2, SMOC1, NFIA, and SEMA6D Disrupted cellular transcripts were identified in the region of integration in four of the seven affected genes.

IMPLICATIONS

Integration of transcriptionally active hrHPV into cellular intergenic regions associates with tumor behavior by altering gene expression. Mol Cancer Res; 14(10); 941-52. ©2016 AACR.

Gene expression in "young adult type" breast cancer: a retrospective analysis

Background: Young women with breast cancer experience inferior outcome and commonly manifest aggressive biological subtypes. Data is controversial regarding biological differences between breast tumors in young (diagnosed at <40 years of age) versus older women. We hypothesize there may be age-related expression differences in key genes for proliferation, invasion and metastasis within and across breast cancer subtypes, and that these differences correlate with outcome.

Methods: Using clinically-annotated gene expression data from 778 breast tumors from three public databases, we compared clinico-pathologic characteristics, mRNA expression of 17 selected genes, and outcome, as a function of age (< 40 years vs. ≥ 40 years).

Results: 14 of 17 genes were differentially expressed in tumors of young vs. older women, 4 of which persisted after correction for subtype and grade (p ≤0.05). BUB1, KRT5, and MYCN were overexpressed and CXCL2 underexpressed in young women. In multivariate analysis, overexpression of cytokeratin genes predicted inferior DFS only for young women. Overexpression of ANGPTL4 strongly predicted inferior DFS in basal but not HER2-enriched tumors in young women. Overexpression of cytokeratin genes and MYBL2 and low SNAI1 expression correlated with inferior DFS in HER2-enriched tumors in younger women. Kaplan-Meier analysis within the basal and HER2-enriched subgroups showed that overexpression of cytokeratin genes was associated with inferior DFS for young, but not older women.

Conclusions: This preliminary study reveals age- and subtype-related differences in expression of key breast cancer genes for proliferation, invasion and metastasis, which correlate with prognostic differences in young women and suggest targeted therapies.

Estrogen Receptor Expression Is Associated with DNA Repair Capacity in Breast Cancer

Estrogen-receptor-positive (ER+) tumors employ complex signaling that engages in crosstalk with multiple pathways through genomic and non-genomic regulation. A greater understanding of these pathways is important for developing improved biomarkers that can better determine treatment choices, risk of recurrence and cancer progression. Deficiencies in DNA repair capacity (DRC) is a hallmark of breast cancer (BC); therefore, in this work we tested whether ER signaling influences DRC. We analyzed the association between ER positivity (% receptor activation) and DRC in 270 BC patients, then further stratified our analysis by HER2 receptor status. Our results show that among HER2 negative, the likelihood of having low DRC values among ER- women is 1.92 (95% CI: 1.03, 3.57) times the likelihood of having low DRC values among ER+ women, even adjusting for different potential confounders (p<0.05); however, a contrary pattern was observed among HER2 positives women. In conclusion, there is an association between DRC levels and ER status, and this association is modified by HER2 receptor status. Adding a DNA repair capacity test to hormone receptor testing may provide new information on defective DNA repair phenotypes, which could better stratify BC patients who have ER+ tumors. ER+/HER2- tumors are heterogeneous, incompletely defined, and clinically challenging to treat; the addition of a DRC test could better characterize and classify these patients as well as help clinicians select optimal therapies, which could improve outcomes and reduce recurrences.

Human Helicase RECQL4 Drives Cisplatin Resistance in Gastric Cancer by Activating an AKT-YB1-MDR1 Signaling Pathway

Elevation of the DNA-unwinding helicase RECQL4, which participates in various DNA repair pathways, has been suggested to contribute to the pathogenicity of various human cancers, including gastric cancer. In this study, we addressed the prognostic and chemotherapeutic significance of RECQL4 in human gastric cancer, which has yet to be determined. We observed significant increases in RECQL4 mRNA or protein in >70% of three independent sets of human gastric cancer specimens examined, relative to normal gastric tissues. Strikingly, high RECQL4 expression in primary tumors correlated well with poor survival and gastric cancer lines with high RECQL4 expression displayed increased resistance to cisplatin treatment. Mechanistic investigations revealed a novel role for RECQL4 in transcriptional regulation of the multidrug resistance gene MDR1, through a physical interaction with the transcription factor YB1. Notably, ectopic expression of RECQL4 in cisplatin-sensitive gastric cancer cells with low endogenous RECQL4 was sufficient to render them resistant to cisplatin, in a manner associated with YB1 elevation and MDR1 activation. Conversely, RECQL4 silencing in cisplatin-resistant gastric cancer cells with high endogenous RECQL4 suppressed YB1 phosphorylation, reduced MDR1 expression, and resensitized cells to cisplatin. In establishing RECQL4 as a critical mediator of cisplatin resistance in gastric cancer cells, our findings provide a therapeutic rationale to target RECQL4 or the downstream AKT-YB1-MDR1 axis to improve gastric cancer treatment. Cancer Res; 76(10); 3057-66. ©2016 AACR.

Meta-Analyses of Microarray Datasets Identifies ANO1 and FADD as Prognostic Markers of Head and Neck Cancer

The head and neck squamous cell carcinoma (HNSCC) transcriptome has been profiled extensively, nevertheless, identifying biomarkers that are clinically relevant and thereby with translational benefit, has been a major challenge. The objective of this study was to use a meta-analysis based approach to catalog candidate biomarkers with high potential for clinical application in HNSCC. Data from publically available microarray series (N = 20) profiled using Agilent (4X44K G4112F) and Affymetrix (HGU133A, U133A_2, U133Plus 2) platforms was downloaded and analyzed in a platform/chip-specific manner (GeneSpring software v12.5, Agilent, USA). Principal Component Analysis (PCA) and clustering analysis was carried out iteratively for segregating outliers; 140 normal and 277 tumor samples from 15 series were included in the final analysis. The analyses identified 181 differentially expressed, concordant and statistically significant genes; STRING analysis revealed interactions between 122 of them, with two major gene clusters connected by multiple nodes (MYC, FOS and HSPA4). Validation in the HNSCC-specific database (N = 528) in The Cancer Genome Atlas (TCGA) identified a panel (ECT2, ANO1, TP63, FADD, EXT1, NCBP2) that was altered in 30% of the samples. Validation in treatment naïve (Group I; N = 12) and post treatment (Group II; N = 12) patients identified 8 genes significantly associated with the disease (Area under curve>0.6). Correlation with recurrence/re-recurrence showed ANO1 had highest efficacy (sensitivity: 0.8, specificity: 0.6) to predict failure in Group I. UBE2V2, PLAC8, FADD and TTK showed high sensitivity (1.00) in Group I while UBE2V2 and CRYM were highly sensitive (>0.8) in predicting re-recurrence in Group II. Further, TCGA analysis showed that ANO1 and FADD, located at 11q13, were co-expressed at transcript level and significantly associated with overall and disease-free survival (p<0.05). The meta-analysis approach adopted in this study has identified candidate markers correlated with disease outcome in HNSCC; further validation in a larger cohort of patients will establish their clinical relevance.

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.