Identification of a gene signature in cell cycle pathway for breast cancer prognosis using gene expression profiling data

NNBGWO-BRCA marker: Neural Network and binary grey wolf optimization based Breast cancer biomarker discovery framework using multi-omics dataset

BACKGROUND AND OBJECTIVE

Breast cancer is a multifaceted condition characterized by diverse features and a substantial mortality rate, underscoring the imperative for timely detection and intervention. The utilization of multi-omics data has gained significant traction in recent years to identify biomarkers and classify subtypes in breast cancer. This kind of research idea from part to whole will also be an inevitable trend in future life science research. Deep learning can integrate and analyze multi-omics data to predict cancer subtypes, which can further drive targeted therapies. However, there are few articles leveraging the nature of deep learning for feature selection. Therefore, this paper proposes a Neural Network and Binary grey Wolf Optimization based BReast CAncer bioMarker (NNBGWO-BRCAMarker) discovery framework using multi-omics data to obtain a series of biomarkers for precise classification of breast cancer subtypes.

METHODS

NNBGWO-BRCAMarker consists of two phases: in the first phase, relevant genes are selected using the weights obtained from a trained feedforward neural network; in the second phase, the binary grey wolf optimization algorithm is leveraged to further screen the selected genes, resulting in a set of potential breast cancer biomarkers.

RESULTS

The SVM classifier with RBF kernel achieved a classification accuracy of 0.9242 ± 0.03 when trained using the 80 biomarkers identified by NNBGWO-BRCAMarker, as evidenced by the experimental results. We conducted a comprehensive gene set analysis, prognostic analysis, and druggability analysis, unveiling 25 druggable genes, 16 enriched pathways strongly linked to specific subtypes of breast cancer, and 8 genes linked to prognostic outcomes.

CONCLUSIONS

The proposed framework successfully identified 80 biomarkers from the multi-omics data, enabling accurate classification of breast cancer subtypes. This discovery may offer novel insights for clinicians to pursue in further studies.

Differential role of glucocorticoid receptor based on its cell type specific expression on tumor cells and infiltrating lymphocytes

BACKGROUND

The glucocorticoid receptor (GR) is frequently expressed in breast cancer (BC), and its prognostic implications are contingent on estrogen receptor (ER) status. To address conflicting reports and explore therapeutic potential, a GR signature (GRsig) independent of ER status was developed. We also investigated cell type-specific GR protein expression in BC tumor epithelial cells and infiltrating lymphocytes.

METHODS

GRsig was derived from Dexamethasone treated cell lines through a bioinformatic pipeline. Immunohistochemistry assessed GR protein expression. Associations between GRsig and tumor phenotypes (proliferation, cytolytic activity (CYT), immune cell distribution, and epithelial-to-mesenchymal transition (EMT) were explored in public datasets. Single-cell RNA sequencing data evaluated context-dependent GR roles, and a cell type-specific prognostic role was assessed in an independent BC cohort.

RESULTS

High GRsig levels were associated with a favorable prognosis across BC subtypes. Tumor-specific high GRsig correlated with lower proliferation, increased CYT, and anti-tumorigenic immune cells. Single-cell data analysis revealed higher GRsig expression in immune cells, negatively correlating with EMT while a positive correlation was observed with EMT primarily in tumor and stromal cells. Univariate and multivariate analyses demonstrated the robust and independent predictive capability of GRsig for favorable prognosis. GR protein expression on immune cells in triple-negative tumors indicated a favorable prognosis.

CONCLUSION

This study underscores the cell type-specific role of GR, where its expression on tumor cells is associated with aggressive features like EMT, while in infiltrating lymphocytes, it predicts a better prognosis, particularly within TNBC tumors. The GRsig emerges as a promising independent prognostic indicator across diverse BC subtypes.

Identification of immune-related gene signature for non-small cell lung cancer patients with immune checkpoint inhibitors

Background

The utilization of immune checkpoint inhibitors (ICIs) has become the established protocol for treating advanced non-small cell lung cancer (NSCLC). This work aimed to identify the immune-related gene signature that can predict the prognosis of NSCLC patients receiving ICI treatment.

Methods

The ImmPort database was queried to obtain a list of immune-related genes (IRGs). Differentially expressed IRGs in NSCLC patients were identified using the TCGA database. RNA-seq data and clinical information from NSCLC patients receiving immunotherapy were obtained from the GEO database (GSE93157 and ////). A gene signature was generated through multivariate Cox and LASSO regression analyses. The prognostic value and function of this gene signature were thoroughly investigated using comprehensive bioinformatics analyses.

Results

A total of 6 prognostic-related genes were identified from 617 differentially expressed genes, and two prognostic-related differentially expressed genes (CAMP and IL17A) were determined to construct gene signature. Our gene signature demonstrated superior performance compared to other clinicopathological parameters in predicting the prognosis of NSCLC patients receiving immunotherapy, with an area under the ROC curve (AUC) of 0.812. Furthermore, immune infiltration analysis indicated that the high-risk group was enriched with resting CD4 T cell memory, while the low-risk group showed a "hot" tumor microenvironment that promotes anti-tumor immunity in NSCLC patients.

Conclusion

Gene signatures based on immune-related genes exhibited excellent indicator performance of prognosis and immune infiltration, which has the potential to be an effective biomarker for NSCLC with ICI treatment.

AMEERA-4: a randomized, preoperative window-of-opportunity study of amcenestrant versus letrozole in early breast cancer

BACKGROUND

Window-of-opportunity (WOO) studies provide insights into the clinical activity of new drugs in breast cancer.

METHODS

AMEERA-4 (NCT04191382) was a WOO study undertaken to compare the pharmacodynamic effects of amcenestrant, a selective estrogen receptor degrader, with those of letrozole in postmenopausal women with newly diagnosed, operable estrogen receptor-positive, human epidermal growth factor receptor 2-negative (ER+/HER2-) breast cancer. Women were randomized (1:1:1) to receive amcenestrant 400 mg, amcenestrant 200 mg, or letrozole 2.5 mg once daily for 14 days before breast surgery. The primary endpoint was change in Ki67 between baseline and Day 15 (i.e., day of surgery).

RESULTS

Enrollment was stopped early because of slow recruitment, in the context of the COVID-19 pandemic. The modified intent-to-treat population consisted of 95 study participants with baseline and post-treatment Ki67 values, whereas the safety population included 104 participants who had received at least one dose of study medication. Relative change from baseline in Ki67 was - 75.9% (95% confidence interval [CI] - 81.9 to - 67.9) for amcenestrant 400 mg, - 68.2% (- 75.7 to - 58.4) for amcenestrant 200 mg, and - 77.7% (- 83.4 to - 70.0) for letrozole (geometric least-squares mean [LSM] estimates). Absolute change in ER H-score from baseline (LSM estimate) was - 176.7 in the amcenestrant 400 mg arm, - 202.9 in the amcenestrant 200 mg arm, and - 32.5 in the letrozole arm. There were no Grade ≥ 3 treatment-related adverse events.

CONCLUSIONS

Both amcenestrant and letrozole demonstrated antiproliferative activity in postmenopausal women with previously untreated, operable ER+/HER2- breast cancer and had good overall tolerability.

TRIAL REGISTRATION

ClinicalTrials.gov, NCT04191382 https://clinicaltrials.gov/ct2/show/NCT04191382 . Registered 9 December 2019.

EARN: an ensemble machine learning algorithm to predict driver genes in metastatic breast cancer

BACKGROUND

Today, there are a lot of markers on the prognosis and diagnosis of complex diseases such as primary breast cancer. However, our understanding of the drivers that influence cancer aggression is limited.

METHODS

In this work, we study somatic mutation data consists of 450 metastatic breast tumor samples from cBio Cancer Genomics Portal. We use four software tools to extract features from this data. Then, an ensemble classifier (EC) learning algorithm called EARN (Ensemble of Artificial Neural Network, Random Forest, and non-linear Support Vector Machine) is proposed to evaluate plausible driver genes for metastatic breast cancer (MBCA). The decision-making strategy for the proposed ensemble machine is based on the aggregation of the predicted scores obtained from individual learning classifiers to be prioritized homo sapiens genes annotated as protein-coding from NCBI.

RESULTS

This study is an attempt to focus on the findings in several aspects of MBCA prognosis and diagnosis. First, drivers and passengers predicted by SVM, ANN, RF, and EARN are introduced. Second, biological inferences of predictions are discussed based on gene set enrichment analysis. Third, statistical validation and comparison of all learning methods are performed by some evaluation metrics. Finally, the pathway enrichment analysis (PEA) using ReactomeFIVIz tool (FDR < 0.03) for the top 100 genes predicted by EARN leads us to propose a new gene set panel for MBCA. It includes HDAC3, ABAT, GRIN1, PLCB1, and KPNA2 as well as NCOR1, TBL1XR1, SIRT4, KRAS, CACNA1E, PRKCG, GPS2, SIN3A, ACTB, KDM6B, and PRMT1. Furthermore, we compare results for MBCA to other outputs regarding 983 primary tumor samples of breast invasive carcinoma (BRCA) obtained from the Cancer Genome Atlas (TCGA). The comparison between outputs shows that ROC-AUC reaches 99.24% using EARN for MBCA and 99.79% for BRCA. This statistical result is better than three individual classifiers in each case.

CONCLUSIONS

This research using an integrative approach assists precision oncologists to design compact targeted panels that eliminate the need for whole-genome/exome sequencing. The schematic representation of the proposed model is presented as the Graphic abstract.

A Clinicogenetic Prognostic Classifier for Prediction of Recurrence and Survival in Asian Breast Cancer Patients

Background

Several prognostic factors affect the recurrence of breast cancer in patients who undergo mastectomy. Assays of the expression profiles of multiple genes increase the probability of overexpression of certain genes and thus can potentially characterize the risk of metastasis.

Methods

We propose a 20-gene classifier for predicting patients with high/low risk of recurrence within 5 years. Gene expression levels from a quantitative PCR assay were used to screen 473 luminal breast cancer patients treated at Taiwan Hospital (positive for estrogen and progesterone receptors, negative for human epidermal growth factor receptor 2). Gene expression scores, along with clinical information (age, tumor stage, and nodal stage), were evaluated for risk prediction. The classifier could correctly predict patients with and without relapse (logistic regression, P<0.05).

Results

A Cox proportional hazards regression analysis showed that the 20-gene panel was prognostic with hazard ratios of 5.63 (95% confidence interval 2.77-11.5, univariate) and 5.56 (2.62-11.8, multivariate) for the “genetic” model, and of 8.02 (3.52-18.3, univariate) and 19.8 (5.96-65.87, multivariate) for the “clinicogenetic” model during a 5-year follow-up.

Conclusions

The proposed 20-gene classifier can successfully separate the patients into two risk groups, and the two risk group had significantly different relapse rate and prognosis. This 20-gene classifier can provide better estimation of prognosis, which can help physicians to make better personalized treatment plans.

SARS-CoV-2 early infection signature identified potential key infection mechanisms and drug targets

Background

The ongoing COVID-19 outbreak has caused devastating mortality and posed a significant threat to public health worldwide. Despite the severity of this illness and 2.3 million worldwide deaths, the disease mechanism is mostly unknown. Previous studies that characterized differential gene expression due to SARS-CoV-2 infection lacked robust validation. Although vaccines are  now available, effective treatment options are still out of reach.

Results

To characterize the transcriptional activity of SARS-CoV-2 infection, a gene signature consisting of 25 genes was generated using a publicly available RNA-Sequencing (RNA-Seq) dataset of cultured cells infected with SARS-CoV-2. The signature estimated infection level accurately in bronchoalveolar lavage fluid (BALF) cells and peripheral blood mononuclear cells (PBMCs) from healthy and infected patients (mean 0.001 vs. 0.958; P < 0.0001). These signature genes were investigated in their ability to distinguish the severity of SARS-CoV-2 infection in a single-cell RNA-Sequencing dataset. TNFAIP3, PPP1R15A, NFKBIA, and IFIT2 had shown bimodal gene expression in various immune cells from severely infected patients compared to healthy or moderate infection cases. Finally, this signature was assessed using the publicly available ConnectivityMap database to identify potential disease mechanisms and drug repurposing candidates. Pharmacological classes of tricyclic antidepressants, SRC-inhibitors, HDAC inhibitors, MEK inhibitors, and drugs such as atorvastatin, ibuprofen, and ketoconazole showed strong negative associations (connectivity score < − 90), highlighting the need for further evaluation of these candidates for their efficacy in treating SARS-CoV-2 infection.

Conclusions

Thus, using the 25-gene SARS-CoV-2 infection signature, the SARS-CoV-2 infection status was captured in BALF cells, PBMCs and postmortem lung biopsies. In addition, candidate SARS-CoV-2 therapies with known safety profiles were identified. The signature genes could potentially also be used to characterize the COVID-19 disease severity in patients’ expression profiles of BALF cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07433-4.

Elevated EXO1 expression is associated with breast carcinogenesis and poor prognosis

Background

Breast cancer is the most common cancer and leading cause of cancer mortality in women worldwide. Exonuclease 1 (EXO1), a protein with 5' to 3' exonuclease and RNase H activity, could be involved in mismatch repair and recombination. This study aims to investigate the prognostic value of EXO1 in breast cancer and explore the association between EXO1 expression and breast carcinogenesis.

Methods

The data of 1,215 breast cancer susceptibility gene (BRCA) samples were obtained from The Cancer Genome Atlas (TCGA). Real-time quantitative polymerase chain reaction (RT-qPCR) further verified the elevated mRNA expression level of EXO1 in human BRCA cells MDA-MB231 compared with that in human breast epithelial cells MCF-10A. EXO1 copy number was proved to be correlated with its expression level. Besides, Kaplan-Meier analysis, differentially expressed genes and function enrichment analysis were performed.

Results

Analysis of data from The Cancer Genome Atlas (TCGA) revealed that the EXO1 expression level in breast cancer tissues was significantly increased. Real-time quantitative polymerase chain reaction (RT-qPCR) supported the elevated mRNA expression level of EXO1 in human breast cancer cells MDA-MB231 compared with that in human breast epithelial cells MCF-10A. EXO1 copy number was shown to be correlated with its expression level. Kaplan-Meier analysis showed that elevated EXO1 was an indicator of poor breast cancer prognosis. Furthermore, differentially expressed genes and function enrichment analysis indicated that the cell cycle pathway and cardiac muscle contraction pathway were activated and inhibited respectively in breast cancer samples with high EXO1 expression.

Conclusions

Therefore, this study shows that elevated EXO1 expression is associated with carcinogenesis and poor prognosis in breast cancer, and might be a biomarker for breast cancer treatment.

Tuberculosis comorbidity with rheumatoid arthritis: Gene signatures, associated biomarkers, and screening

Rheumatoid arthritis (RA) is known to be related to an elevated risk of infections because of its pathobiology and the use of immunosuppressive therapies. Reactivation of latent tuberculosis (TB) infection is a serious issue in patients with RA, especially after receiving anti-TNFs therapy. TNF blocking reinforces the TB granuloma formation and maintenance and the growth of Mycobacterium tuberculosis (Mtb). After intercurrent of TB infection, the standard recommendation is that the treatment with TNF inhibitors to be withheld despite its impressive effect on suppression of inflammation until the infection has resolved. Knowing pathways and mechanisms that are common between two diseases might help to find the mechanistic basis of this comorbidity, as well as provide us a new approach to apply them as therapeutic targets or diagnostic biomarkers. Also, screening for latent TB before initiation of an anti-TNF therapy can minimize complications. This review summarizes the shared gene signature between TB and RA and discusses the biomarkers for early detection of this infection, and screening procedures as well.

An integrative multi-omics network-based approach identifies key regulators for breast cancer

Although genome-wide association studies (GWASs) have successfully identified thousands of risk variants for human complex diseases, understanding the biological function and molecular mechanisms of the associated SNPs involved in complex diseases is challenging. Here we developed a framework named integrative multi-omics network-based approach (IMNA), aiming to identify potential key genes in regulatory networks by integrating molecular interactions across multiple biological scales, including GWAS signals, gene expression-based signatures, chromatin interactions and protein interactions from the network topology. We applied this approach to breast cancer, and prioritized key genes involved in regulatory networks. We also developed an abnormal gene expression score (AGES) signature based on the gene expression deviation of the top 20 rank-ordered genes in breast cancer. The AGES values are associated with genetic variants, tumor properties and patient survival outcomes. Among the top 20 genes, RNASEH2A was identified as a new candidate gene for breast cancer. Thus, our integrative network-based approach provides a genetic-driven framework to unveil tissue-specific interactions from multiple biological scales and reveal potential key regulatory genes for breast cancer. This approach can also be applied in other complex diseases such as ovarian cancer to unravel underlying mechanisms and help for developing therapeutic targets.

Pathway Mutations in Breast Cancer Using Whole-Exome Sequencing

The genomic landscape of breast cancer (BC) is complex. The purpose of this study was to decipher the mutational profiles of Taiwanese patients with BC using next-generation sequencing. We performed whole-exome sequencing on DNA from 24 tumor tissue specimens from BC patients. Sanger sequencing was used to validate the identified variants. Sanger sequencing was also performed on paired adjacent nontumor tissues. After genotype calling and algorithmic annotations, we identified 49 deleterious variants in canonical cancer-related genes in our BC cohort. The most frequently mutated genes were PIK3CA (16.67%), FKBP9 (12.5%), TP53 (12.5%), ATM (8.33%), CHEK2 (8.33%), FOXO3 (8.33%), NTRK1 (8.33%), and NUTM2B (8.33%). Seven mutated variants (ATR p.V1581fs, CSF1R p.R579Q, GATA3 p.T356delinsTMKS, LRP5 p.W389*, MAP3K1 p.T918fs, MET p.K1161fs, and MTR p.P1178S) were novel variants that are not present in any gene mutation database. After grouping the samples according to molecular subtype, we found that the cell cycle, MAPK, and chemokine signaling pathways in the luminal A subtype of BC; the focal adhesion, axon guidance, and endocytosis pathways in the luminal B subtype; and amyotrophic lateral sclerosis in the basal-like subtype were exclusively altered. Survival curve analysis showed that the presence of the MAPK signaling pathway and endocytosis mutations were correlated with a poor prognosis. These survival data were consistent with cBioPortal analyses of 2,051 BC cases. We discovered novel mutations in patients with BC. These results have implications for developing strategic, adjuvant, and gene-targeted therapies.

Weighted gene expression profiles identify diagnostic and prognostic genes for lung adenocarcinoma and squamous cell carcinoma

OBJECTIVE

To construct a diagnostic signature to distinguish lung adenocarcinoma from lung squamous cell carcinoma and a prognostic signature to predict the risk of death for patients with nonsmall-cell lung cancer, with satisfactory predictive performances, good stabilities, small sizes and meaningful biological implications.

METHODS

Pathway-based feature selection methods utilize pathway information as a priori to provide insightful clues on potential biomarkers from the biological perspective, and such incorporation may be realized by adding weights to test statistics or gene expression values. In this study, weighted gene expression profiles were generated using the GeneRank method and then the LASSO method was used to identify discriminative and prognostic genes.

RESULTS

The five-gene diagnostic signature including keratin 5 (KRT5), mucin 1 (MUC1), triggering receptor expressed on myeloid cells 1 (TREM1), complement C3 (C3) and transmembrane serine protease 2 (TMPRSS2) achieved a predictive error of 12.8% and a Generalized Brier Score of 0.108, while the five-gene prognostic signature including alcohol dehydrogenase 1C (class I), gamma polypeptide (ADH1C), alpha-2-glycoprotein 1, zinc-binding (AZGP1), clusterin (CLU), cyclin dependent kinase 1 (CDK1) and paternally expressed 10 (PEG10) obtained a log-rank P-value of 0.03 and a C-index of 0.622 on the test set.

CONCLUSIONS

Besides good predictive capacity, model parsimony and stability, the identified diagnostic and prognostic genes were highly relevant to lung cancer. A large-sized prospective study to explore the utilization of these genes in a clinical setting is warranted.

Gene relevance based on multiple evidences in complex networks

MOTIVATION

Multi-omics approaches offer the opportunity to reconstruct a more complete picture of the molecular events associated with human diseases, but pose challenges in data analysis. Network-based methods for the analysis of multi-omics leverage the complex web of macromolecular interactions occurring within cells to extract significant patterns of molecular alterations. Existing network-based approaches typically address specific combinations of omics and are limited in terms of the number of layers that can be jointly analysed. In this study, we investigate the application of network diffusion to quantify gene relevance on the basis of multiple evidences (layers).

RESULTS

We introduce a gene score (mND) that quantifies the relevance of a gene in a biological process taking into account the network proximity of the gene and its first neighbours to other altered genes. We show that mND has a better performance over existing methods in finding altered genes in network proximity in one or more layers. We also report good performances in recovering known cancer genes. The pipeline described in this article is broadly applicable, because it can handle different types of inputs: in addition to multi-omics datasets, datasets that are stratified in many classes (e.g., cell clusters emerging from single cell analyses) or a combination of the two scenarios.

AVAILABILITY AND IMPLEMENTATION

The R package 'mND' is available at URL: https://www.itb.cnr.it/mnd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genes and functions from breast cancer signatures

BACKGROUND

Breast cancer is a heterogeneous disease and personalized medicine is the hope for the improvement of the clinical outcome. Multi-gene signatures for breast cancer stratification have been extensively studied in the past decades and more than 30 different signatures have been reported. A major concern is the minimal overlap of genes among the reported signatures. We investigated the breast cancer signature genes to address our hypothesis that the genes of different signature may share common functions, as well as to use these previously reported signature genes to build better prognostic models.

METHODS

A total of 33 signatures and the corresponding gene lists were investigated. We first examined the gene frequency and the gene overlap in these signatures. Then the gene functions of each signature gene list were analysed and compared by the KEGG pathways and gene ontology (GO) terms. A classifier built using the common genes was tested using the METABRIC (Molecular Taxonomy of Breast Cancer International Consortium) data. The common genes were also tested for building the Yin Yang gene mean expression ratio (YMR) signature using public datasets (GSE1456 and GSE2034).

RESULTS

Among a total of 2239 genes collected from the 33 breast cancer signatures, only 238 genes overlapped in at least two signatures; while from a total of 1979 function terms enriched in the 33 signature gene lists, 429 terms were common in at least two signatures. Most of the common function terms were involved in cell cycle processes. While there is almost no common overlapping genes between signatures developed for ER-positive (e.g. 21-gene signature) and those developed for ER-negative (e.g. basal signatures) tumours, they have common function terms such as cell death, regulation of cell proliferation. We used the 62 genes that were common in at least three signatures as a classifier and subtyped 1141 METABRIC cases including 144 normal samples into nine subgroups. These subgroups showed different clinical outcome. Among the 238 common genes, we selected those genes that are more highly expressed in normal breast tissue than in tumours as Yang genes and those more highly expressed in tumours than in normal as Yin genes and built a YMR model signature. This YMR showed significance in risk stratification in two datasets (GSE1456 and GSE2034).

CONCLUSIONS

The lack of significant numbers of overlapping genes among most breast cancer signatures can be partially explained by our discovery that these signature genes represent groups with similar functions. The genes collected from these previously reported signatures are valuable resources for new model development. The subtype classifier and YMR signature built from the common genes showed promising results.

Functional characterization of breast cancer using pathway profiles

BACKGROUND

The molecular characteristics of human diseases are often represented by a list of genes termed "signature genes". A significant challenge facing this approach is that of reproducibility: signatures developed on a set of patients may fail to perform well on different sets of patients. As diseases are resulted from perturbed cellular functions, irrespective of the particular genes that contribute to the function, it may be more appropriate to characterize diseases based on these perturbed cellular functions.

METHODS

We proposed a profile-based approach to characterize a disease using a binary vector whose elements indicate whether a given function is perturbed based on the enrichment analysis of expression data between normal and tumor tissues. Using breast cancer and its four primary clinically relevant subtypes as examples, this approach is evaluated based on the reproducibility, accuracy and resolution of the resulting pathway profiles.

RESULTS

Pathway profiles for breast cancer and its subtypes are constructed based on data obtained from microarray and RNA-Seq data sets provided by The Cancer Genome Atlas (TCGA), and an additional microarray data set provided by The European Genome-phenome Archive (EGA). An average reproducibility of 68% is achieved between different data sets (TCGA microarray vs. EGA microarray data) and 67% average reproducibility is achieved between different technologies (TCGA microarray vs. TCGA RNA-Seq data). Among the enriched pathways, 74% of them are known to be associated with breast cancer or other cancers. About 40% of the identified pathways are enriched in all four subtypes, with 4, 2, 4, and 7 pathways enriched only in luminal A, luminal B, triple-negative, and HER2+ subtypes, respectively. Comparison of profiles between subtypes, as well as other diseases, shows that luminal A and luminal B subtypes are more similar to the HER2+ subtype than to the triple-negative subtype, and subtypes of breast cancer are more likely to be closer to each other than to other diseases.

CONCLUSIONS

Our results demonstrate that pathway profiles can successfully characterize both common and distinct functional characteristics of four subtypes of breast cancer and other related diseases, with acceptable reproducibility, high accuracy and reasonable resolution.

Chemotherapy-induced dynamic gene expression changes in vivo are prognostic in ovarian cancer

BACKGROUND

The response of ovarian cancer patients to carboplatin and paclitaxel is variable, necessitating identification of biomarkers that can reliably predict drug sensitivity and resistance. In this study, we sought to identify dynamically controlled genes and pathways associated with drug response and its time dependence.

METHODS

Gene expression was assessed for 14 days post-treatment with carboplatin or carboplatin-paclitaxel in xenografts from two ovarian cancer models: platinum-sensitive serous adenocarcinoma-derived OV1002 and a mixed clear cell/endometrioid carcinoma-derived HOX424 with reduced sensitivity to platinum.

RESULTS

Tumour volume reduction was observed in both xenografts, but more dominantly in OV1002. Upregulated genes in OV1002 were involved in DNA repair, cell cycle and apoptosis, whereas downregulated genes were involved in oxygen-consuming metabolic processes and apoptosis control. Carboplatin-paclitaxel triggered a more comprehensive response than carboplatin only in both xenografts. In HOX424, apoptosis and cell cycle were upregulated, whereas Wnt signalling was inhibited. Genes downregulated after day 7 from both xenografts were predictive of overall survival. Overrepresented pathways were also predictive of outcome.

CONCLUSIONS

Late expressed genes are prognostic in ovarian tumours in a dynamic manner. This longitudinal gene expression study further elucidates chemotherapy response in two models, stressing the importance of delayed biomarker detection and guiding optimal timing of biopsies.

HIF-1 is involved in the negative regulation of AURKA expression in breast cancer cell lines under hypoxic conditions

Numerous microarray-based gene expression studies performed on several types of solid tumors revealed significant changes in key genes involved in progression and regulation of the cell cycle, including AURKA that is known to be overexpressed in many types of human malignancies. Tumor hypoxia is associated with poor prognosis in several cancer types, including breast cancer (BC). Since hypoxia is a condition that influences the expression of many genes involved in tumorigenesis, proliferation, and cell cycle regulation, we performed a microarray-based gene expression analysis in order to identify differentially expressed genes in BC cell lines exposed to hypoxia. This analysis showed that hypoxia induces a down-regulation of AURKA expression. Although hypoxia is a tumor feature, the molecular mechanisms that regulate AURKA expression in response to hypoxia in BC are still unknown. For the first time, we demonstrated that HIF-1 activation downstream of hypoxia could drive AURKA down-regulation in BC cells. In fact, we found that siRNA-mediated knockdown of HIF-1α significantly reduces the AURKA down-regulation in BC cells under hypoxia. The aim of our study was to obtain new insights into AURKA transcriptional regulation in hypoxic conditions. Luciferase reporter assays showed a reduction of AURKA promoter activity in hypoxia. Unlike the previous findings, we hypothesize a new possible mechanism where HIF-1, rather than inducing transcriptional activation, could promote the AURKA down-regulation via its binding to hypoxia-responsive elements into the proximal region of the AURKA promoter. The present study shows that hypoxia directly links HIF-1 with AURKA expression, suggesting a possible pathophysiological role of this new pathway in BC and confirming HIF-1 as an important player linking an environmental signal to the AURKA promoter. Since AURKA down-regulation overrides the estrogen-mediated growth and chemoresistance in BC cells, these findings could be important for the development of new possible therapies against BC.

The inference of breast cancer metastasis through gene regulatory networks

Understanding the mechanisms of gene regulation during breast cancer is one of the most difficult problems among oncologists because this regulation is likely comprised of complex genetic interactions. Given this complexity, a computational study using the Bayesian network technique has been employed to construct a gene regulatory network from microarray data. Although the Bayesian network has been notified as a prominent method to infer gene regulatory processes, learning the Bayesian network structure is NP hard and computationally intricate. Therefore, we propose a novel inference method based on low-order conditional independence that extends to the case of the Bayesian network to deal with a large number of genes and an insufficient sample size. This method has been evaluated and compared with full-order conditional independence and different prognostic indices on a publicly available breast cancer data set. Our results suggest that the low-order conditional independence method will be able to handle a large number of genes in a small sample size with the least mean square error. In addition, this proposed method performs significantly better than other methods, including the full-order conditional independence and the St. Gallen consensus criteria. The proposed method achieved an area under the ROC curve of 0.79203, whereas the full-order conditional independence and the St. Gallen consensus criteria obtained 0.76438 and 0.73810, respectively. Furthermore, our empirical evaluation using the low-order conditional independence method has demonstrated a promising relationship between six gene regulators and two regulated genes and will be further investigated as potential breast cancer metastasis prognostic markers.

Predictive Power Estimation Algorithm (PPEA)--a new algorithm to reduce overfitting for genomic biomarker discovery

Toxicogenomics promises to aid in predicting adverse effects, understanding the mechanisms of drug action or toxicity, and uncovering unexpected or secondary pharmacology. However, modeling adverse effects using high dimensional and high noise genomic data is prone to over-fitting. Models constructed from such data sets often consist of a large number of genes with no obvious functional relevance to the biological effect the model intends to predict that can make it challenging to interpret the modeling results. To address these issues, we developed a novel algorithm, Predictive Power Estimation Algorithm (PPEA), which estimates the predictive power of each individual transcript through an iterative two-way bootstrapping procedure. By repeatedly enforcing that the sample number is larger than the transcript number, in each iteration of modeling and testing, PPEA reduces the potential risk of overfitting. We show with three different cases studies that: (1) PPEA can quickly derive a reliable rank order of predictive power of individual transcripts in a relatively small number of iterations, (2) the top ranked transcripts tend to be functionally related to the phenotype they are intended to predict, (3) using only the most predictive top ranked transcripts greatly facilitates development of multiplex assay such as qRT-PCR as a biomarker, and (4) more importantly, we were able to demonstrate that a small number of genes identified from the top-ranked transcripts are highly predictive of phenotype as their expression changes distinguished adverse from nonadverse effects of compounds in completely independent tests. Thus, we believe that the PPEA model effectively addresses the over-fitting problem and can be used to facilitate genomic biomarker discovery for predictive toxicology and drug responses.

Integration and analysis of genome-scale data from gliomas

Primary brain tumors are a leading cause of cancer-related mortality among young adults and children. The most common primary malignant brain tumor, glioma, carries a median survival of only 14 months. Two major multi-institutional programs, the Glioma Molecular Diagnostic Initiative and The Cancer Genome Atlas, have pursued a comprehensive genomic characterization of a large number of clinical glioma samples using a variety of technologies to measure gene expression, chromosomal copy number alterations, somatic and germline mutations, DNA methylation, microRNA, and proteomic changes. Classification of gliomas on the basis of gene expression has revealed six major subtypes and provided insights into the underlying biology of each subtype. Integration of genome-wide data from different technologies has been used to identify many potential protein targets in this disease, while increasing the reliability and biological interpretability of results. Mapping genomic changes onto both known and inferred cellular networks represents the next level of analysis, and has yielded proteins with key roles in tumorigenesis. Ultimately, the information gained from these approaches will be used to create customized therapeutic regimens for each patient based on the unique genomic signature of the individual tumor. In this Review, we describe efforts to characterize gliomas using genomic data, and consider how insights gained from these analyses promise to increase understanding of the biological underpinnings of the disease and lead the way to new therapeutic strategies.

Genetic tumor profiling and genetically targeted cancer therapy

OBJECTIVES

To discuss how understanding and manipulation of tumor genetics information and technology shapes cancer care today and what changes might be expected in the near future.

DATA SOURCES

Published articles, web resources, clinical practice.

CONCLUSIONS

Advances in our understanding of genes and their regulation provide a promise of more personalized cancer care, allowing selection of the most safe and effective therapy in an individual situation.

IMPLICATIONS FOR NURSING PRACTICE

Rapid progress in the technology of tumor profiling and targeted cancer therapies challenges nurses to keep up-to-date to provide quality patient education and care.

Dynamic changes in gene expression in vivo predict prognosis of tamoxifen-treated patients with breast cancer

INTRODUCTION

Tamoxifen is the most widely prescribed anti-estrogen treatment for patients with estrogen receptor (ER)-positive breast cancer. However, there is still a need for biomarkers that reliably predict endocrine sensitivity in breast cancers and these may well be expressed in a dynamic manner.

METHODS

In this study we assessed gene expression changes at multiple time points (days 1, 2, 4, 7, 14) after tamoxifen treatment in the ER-positive ZR-75-1 xenograft model that displays significant changes in apoptosis, proliferation and angiogenesis within 2 days of therapy.

RESULTS

Hierarchical clustering identified six time-related gene expression patterns, which separated into three groups: two with early/transient responses, two with continuous/late responses and two with variable response patterns. The early/transient response represented reductions in many genes that are involved in cell cycle and proliferation (e.g. BUB1B, CCNA2, CDKN3, MKI67, UBE2C), whereas the continuous/late changed genes represented the more classical estrogen response genes (e.g. TFF1, TFF3, IGFBP5). Genes and the proteins they encode were confirmed to have similar temporal patterns of expression in vitro and in vivo and correlated with reduction in tumour volume in primary breast cancer. The profiles of genes that were most differentially expressed on days 2, 4 and 7 following treatment were able to predict prognosis, whereas those most changed on days 1 and 14 were not, in four tamoxifen treated datasets representing a total of 404 patients.

CONCLUSIONS

Both early/transient/proliferation response genes and continuous/late/estrogen-response genes are able to predict prognosis of primary breast tumours in a dynamic manner. Temporal expression of therapy-response genes is clearly an important factor in characterising the response to endocrine therapy in breast tumours which has significant implications for the timing of biopsies in neoadjuvant biomarker studies.

Spectral Clustering of Microarray Data Elucidates the Roles of Microenvironment Remodeling and Immune Responses in Survival of Head and Neck Squamous Cell Carcinoma

Purpose

To identify functionally related prognostic gene sets for head and neck squamous cell carcinoma (HNSCC) by unsupervised statistical analysis of microarray data.

Patients and Methods

Microarray analysis was performed on 14 normal oral epithelium and 71 HNSCCs from patients with outcome data. Spectral clustering (SC) analysis of the data set identified multiple vectors representing distinct aspects of gene expression heterogeneity between samples. Gene ontology (GO) analysis of vector gene lists identified gene sets significantly enriched within defined biologic pathways. The prognostic significance of these was established by Cox survival analysis.

Results

The most influential SC vectors were V2 and V3. V2 separated normal from tumor samples. GO analysis of V2 gene lists identified pathways with heterogeneous expression between HNSCCs, notably focal adhesion (FA)/extracellular matrix remodeling and cytokine-cytokine receptor (CR) interactions. Similar analysis of V3 gene lists identified further heterogeneity in CR pathways. V2CR genes represent an innate immune response, whereas high expression of V3CR genes represented an adaptive immune response that was not dependent on human papillomavirus status. Survival analysis demonstrated that the FA gene set was prognostic of poor outcome, whereas classification for adaptive immune response by the CR gene set was prognostic of good outcome. A combined FA&CR model dramatically exceeded the performance of current clinical classifiers (P < .001 in our cohort and, importantly, P = .007 in an independent cohort of 60 HNSCCs).

Conclusion

The application of SC and GO algorithms to HNSCC microarray data identified gene sets highly significant for predicting patient outcome. Further large-scale studies will establish the usefulness of these gene sets in the clinical management of HNSCC.

[Genome-wide expression profiling as a clinical tool: are we there yet?]

Breast cancer is a heterogeneous disease, encompassing a plethora of histological types and clinical courses. Current histopathological classification systems for breast cancer are based on descriptive entities that are of prognostic significance. Few prognostic markers beyond those offered by histopathological analysis are available. Furthermore, a very limited armamentarium of predictive biomarkers has been introduced in clinical practice. High throughput molecular technologies are reshaping our understanding of breast cancer, of which microarray-based gene expression has received the most attention. This method has been successfully used to derive a molecular taxonomy for breast cancer, which has provided interesting insights into the biology of the disease. Microarray-based class prediction studies have generated a multitude of prognostic/predictive signatures. Although these signatures have not been fully translated to clinical practice as yet, they herald the promise of an improvement in breast cancer treatment decision-making. It should be noted, however, that most of the signatures developed to date seem to have discriminatory power almost restricted to oestrogen receptor-positive disease. This review addresses the contribution of gene expression profiling to our understanding of breast cancer and its clinical management and what has yet to be done for these classifiers to be incorporated in clinical practice.

Microarray-based Gene Expression Profiling as a Clinical Tool for Breast Cancer Management: Are We There Yet?

Breast cancer is a heterogeneous disease, encompassing several histological types and clinical behaviors. Current histopathological classification systems are based on descriptive entities with prognostic significance. Few prognostic and predictive markers beyond those offered by histopathological analysis are available. High-throughput molecular technologies are reshaping our understanding of breast cancer, of which microarray-based gene expression has received most attention. This method has been used to derive a molecular taxonomy for breast cancer, which has provided interesting insights into the biology of the disease. Class prediction studies have generated a multitude of prognostic/predictive signatures, which herald the promise for an improvement in treatment decision making. However, most of the signatures developed to date seem to have discriminatory power almost restricted to estrogen receptor—positive disease. This review addresses the contribution of gene expression profiling to our understanding of breast cancer and its clinical management.