miRNA-mRNA integrated analysis reveals roles for miRNAs in primary breast tumors

Introduction

Few studies have performed expression profiling of both miRNA and mRNA from the same primary breast carcinomas. In this study we present and analyze data derived from expression profiling of 799 miRNAs in 101 primary human breast tumors, along with genome-wide mRNA profiles and extensive clinical information.

Methods

We investigate the relationship between these molecular components, in terms of their correlation with each other and with clinical characteristics. We use a systems biology approach to examine the correlative relationship between miRNA and mRNAs using statistical enrichment methods.

Results

We identify statistical significant differential expression of miRNAs between molecular intrinsic subtypes, and between samples with different levels of proliferation. Specifically, we point to miRNAs significantly associated with TP53 and ER status. We also show that several cellular processes, such as proliferation, cell adhesion and immune response, are strongly associated with certain miRNAs. We validate the role of miRNAs in regulating proliferation using high-throughput lysate-microarrays on cell lines and point to potential drivers of this process.

Conclusion

This study provides a comprehensive dataset as well as methods and system-level results that jointly form a basis for further work on understanding the role of miRNA in primary breast cancer.     

Gene Expression-Based Classification of Non-Small Cell Lung Carcinomas and Survival Prediction

Background

Current clinical therapy of non-small cell lung cancer depends on histo-pathological classification. This approach poorly predicts clinical outcome for individual patients. Gene expression profiling holds promise to improve clinical stratification, thus paving the way for individualized therapy.

Methodology and Principal Findings

A genome-wide gene expression analysis was performed on a cohort of 91 patients. We used 91 tumor- and 65 adjacent normal lung tissue samples. We defined sets of predictor genes (probe sets) with the expression profiles. The power of predictor genes was evaluated using an independent cohort of 96 non-small cell lung cancer- and 6 normal lung samples. We identified a tumor signature of 5 genes that aggregates the 156 tumor and normal samples into the expected groups. We also identified a histology signature of 75 genes, which classifies the samples in the major histological subtypes of non-small cell lung cancer. Correlation analysis identified 17 genes which showed the best association with post-surgery survival time. This signature was used for stratification of all patients in two risk groups. Kaplan-Meier survival curves show that the two groups display a significant difference in post-surgery survival time (p = 5.6E-6). The performance of the signatures was validated using a patient cohort of similar size (Duke University, n = 96). Compared to previously published prognostic signatures for NSCLC, the 17 gene signature performed well on these two cohorts.

Conclusions

The gene signatures identified are promising tools for histo-pathological classification of non-small cell lung cancer, and may improve the prediction of clinical outcome.     

A Predictive Phosphorylation Signature of Lung Cancer

Background

Aberrant activation of signaling pathways drives many of the fundamental biological processes that accompany tumor initiation and progression. Inappropriate phosphorylation of intermediates in these signaling pathways are a frequently observed molecular lesion that accompanies the undesirable activation or repression of pro- and anti-oncogenic pathways. Therefore, methods which directly query signaling pathway activation via phosphorylation assays in individual cancer biopsies are expected to provide important insights into the molecular “logic” that distinguishes cancer and normal tissue on one hand, and enables personalized intervention strategies on the other.

Results

We first document the largest available set of tyrosine phosphorylation sites that are, individually, differentially phosphorylated in lung cancer, thus providing an immediate set of drug targets. Next, we develop a novel computational methodology to identify pathways whose phosphorylation activity is strongly correlated with the lung cancer phenotype. Finally, we demonstrate the feasibility of classifying lung cancers based on multi-variate phosphorylation signatures.

Conclusions

Highly predictive and biologically transparent phosphorylation signatures of lung cancer provide evidence for the existence of a robust set of phosphorylation mechanisms (captured by the signatures) present in the majority of lung cancers, and that reliably distinguish each lung cancer from normal. This approach should improve our understanding of cancer and help guide its treatment, since the phosphorylation signatures highlight proteins and pathways whose phosphorylation should be inhibited in order to prevent unregulated proliferation.     

Characterizing tyrosine phosphorylation signaling in lung cancer using SH2 profiling

Background

Tyrosine kinases drive the proliferation and survival of many human cancers. Thus profiling the global state of tyrosine phosphorylation of a tumor is likely to provide a wealth of information that can be used to classify tumors for prognosis and prediction. However, the comprehensive analysis of tyrosine phosphorylation of large numbers of human cancer specimens is technically challenging using current methods.

Methodology/Principal Findings

We used a phosphoproteomic method termed SH2 profiling to characterize the global state of phosphotyrosine (pTyr) signaling in human lung cancer cell lines. This method quantifies the phosphorylated binding sites for SH2 domains, which are used by cells to respond to changes in pTyr during signaling. Cells could be grouped based on SH2 binding patterns, with some clusters correlated with EGF receptor (EGFR) or K-RAS mutation status. Binding of specific SH2 domains, most prominently RAS pathway activators Grb2 and ShcA, correlated with EGFR mutation and sensitivity to the EGFR inhibitor erlotinib. SH2 binding patterns also reflected MET activation and could identify cells driven by multiple kinases. The pTyr responses of cells treated with kinase inhibitors provided evidence of distinct mechanisms of inhibition.

Conclusions/Significance

This study illustrates the potential of modular protein domains and their proteomic binding profiles as powerful molecular diagnostic tools for tumor classification and biomarker identification.     

DNA Methylation Profiles and Their Relationship with Cytogenetic Status in Adult Acute Myeloid Leukemia

Background

Aberrant promoter DNA methylation has been shown to play a role in acute myeloid leukemia (AML) pathophysiology. However, further studies to discuss the prognostic value and the relationship of the epigenetic signatures with defined genomic rearrangements in acute myeloid leukemia are required.

Methodology/Principal Findings

We carried out high-throughput methylation profiling on 116 de novo AML cases and we validated the significant biomarkers in an independent cohort of 244 AML cases. Methylation signatures were associated with the presence of a specific cytogenetic status. In normal karyotype cases, aberrant methylation of the promoter of DBC1 was validated as a predictor of the disease-free and overall survival. Furthermore, DBC1 expression was significantly silenced in the aberrantly methylated samples. Patients with chromosome rearrangements showed distinct methylation signatures. To establish the role of fusion proteins in the epigenetic profiles, 20 additional samples of human hematopoietic stem/progenitor cells (HSPC) transduced with common fusion genes were studied and compared with patient samples carrying the same rearrangements. The presence of MLL rearrangements in HSPC induced the methylation profile observed in the MLL-positive primary samples. In contrast, fusion genes such as AML1/ETO or CBFB/MYH11 failed to reproduce the epigenetic signature observed in the patients.

Conclusions/Significance

Our study provides a comprehensive epigenetic profiling of AML, identifies new clinical markers for cases with a normal karyotype, and reveals relevant biological information related to the role of fusion proteins on the methylation signature.     

Similar Source of Differential Blood mRNAs in Lung Cancer and Pulmonary Inflammatory Diseases: Calls for Improved Strategy for Identifying Cancer-Specific Biomarkers

Background

Many studies try to identify cancer diagnostic biomarkers by comparing peripheral whole blood (PWB) of cancer samples and healthy controls, explicitly or implicitly assuming that such biomarkers are potential candidate biomarkers for distinguishing cancer from nonmalignant inflammation-associated diseases.

Methods

Multiple PWB gene expression profiles for lung cancer/inflammation-associated pulmonary diseases were used for differential mRNAs identification and comparison and for proportion estimation of PWB cell subtypes.

Results

The differentially expressed genes (DE genes) between lung cancer/inflammation-associated pulmonary patients and healthy controls were reproducibly identified in different datasets. For these DE genes observed in lung cancer/inflammation-associated pulmonary diseases, more than 90.2% were differentially expressed between myeloid cells and lymphoid cells, with at least 96.8% having consistent directions of regulation (up- or down-regulations) in myeloid cells compared to lymphoid cells, explainable by the shifted populations of PWB cell subtypes under the disease conditions. The comparison of DE genes for lung cancer and inflammation-associated pulmonary diseases showed that the overlapping genes were 100% consistent in the sense of direction of regulation.

Conclusions

The differential blood mRNAs observed in lung cancer and in inflammation-associated pulmonary diseases were similar, both mainly reflecting the difference between myeloid cells and lymphoid cells predominantly determined by PWB cell population shifts. Thus, the strategy of comparing cancer with healthy controls may provide little information of the ability of the identified candidate biomarkers in discriminating cancer from inflammation-associated pulmonary diseases.     

Experimental validation of methods for differential gene expression analysis and sample pooling in RNA-seq

Background

Massively parallel cDNA sequencing (RNA-seq) experiments are gradually superseding microarrays in quantitative gene expression profiling. However, many biologists are uncertain about the choice of differentially expressed gene (DEG) analysis methods and the validity of cost-saving sample pooling strategies for their RNA-seq experiments. Hence, we performed experimental validation of DEGs identified by Cuffdiff2, edgeR, DESeq2 and Two-stage Poisson Model (TSPM) in a RNA-seq experiment involving mice amygdalae micro-punches, using high-throughput qPCR on independent biological replicate samples. Moreover, we sequenced RNA-pools and compared their results with sequencing corresponding individual RNA samples.

Results

False-positivity rate of Cuffdiff2 and false-negativity rates of DESeq2 and TSPM were high. Among the four investigated DEG analysis methods, sensitivity and specificity of edgeR was relatively high. We documented the pooling bias and that the DEGs identified in pooled samples suffered low positive predictive values.

Conclusions

Our results highlighted the need for combined use of more sensitive DEG analysis methods and high-throughput validation of identified DEGs in future RNA-seq experiments. They indicated limited utility of sample pooling strategies for RNA-seq in similar setups and supported increasing the number of biological replicate samples.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1767-y) contains supplementary material, which is available to authorized users.     

Biological effects of a de novo designed myxoma virus peptide analogue: evaluation of cytotoxicity on tumor cells

Background

The Resonant Recognition Model (RRM) is a physico-mathematical model that interprets protein sequence linear information using digital signal processing methods. In this study the RRM concept was employed for structure-function analysis of myxoma virus (MV) proteins and the design of a short bioactive therapeutic peptide with MV-like antitumor/cytotoxic activity.

Methodology/Principal Findings

The analogue RRM-MV was designed by RRM as a linear 18 aa 2.3 kDa peptide. The biological activity of this computationally designed peptide analogue against cancer and normal cell lines was investigated. The cellular cytotoxicity effects were confirmed by confocal immunofluorescence microscopy, by measuring the levels of cytoplasmic lactate dehydrogenase (LDH) and by Prestoblue cell viability assay for up to 72 hours in peptide treated and non-treated cell cultures. Our results revealed that RRM-MV induced a significant dose and time-dependent cytotoxic effect on murine and human cancer cell lines. Yet, when normal murine cell lines were similarly treated with RRM-MV, no cytotoxic effects were observed. Furthermore, the non-bioactive RRM designed peptide RRM-C produced negligible cytotoxic effects on these cancer and normal cell lines when used at similar concentrations. The presence/absence of phosphorylated Akt activity in B16F0 mouse melanoma cells was assessed to indicate the possible apoptosis signalling pathway that could be affected by the peptide treatment. So far, Akt activity did not seem to be significantly affected by RRM-MV as is the case for the original viral protein.

Conclusions/Significance

Our findings indicate the successful application of the RRM concept to design a bioactive peptide analogue (RRM-MV) with cytotoxic effects on tumor cells only. This 2.345 kDa peptide analogue to a 49 kDa viral protein may be suitable to be developed as a potential cancer therapeutic. These results also open a new direction to the rational design of therapeutic agents for future cancer treatment.     

Wnt Pathway Activity in Breast Cancer Sub-Types and Stem-Like Cells

Introduction

Wnt signalling has been implicated in stem cell regulation however its role in breast cancer stem cell regulation remains unclear.

Methods

We used a panel of normal and breast cancer cell lines to assess Wnt pathway gene and protein expression, and for the investigation of Wnt signalling within stem cell-enriched populations, mRNA and protein expression was analysed after the selection of anoikis-resistant cells. Finally, cell lines and patient-derived samples were used to investigate Wnt pathway effects on stem cell activity in vitro.

Results

Wnt pathway signalling increased in cancer compared to normal breast and in both cell lines and patient samples, expression of Wnt pathway genes correlated with estrogen receptor (ER) expression. Furthermore, specific Wnt pathway genes were predictive for recurrence within subtypes of breast cancer. Canonical Wnt pathway genes were increased in breast cancer stem cell-enriched populations in comparison to normal breast stem cell-enriched populations. Furthermore in cell lines, the ligand Wnt3a increased whilst the inhibitor DKK1 reduced mammosphere formation with the greatest inhibitory effects observed in ER+ve breast cancer cell lines. In patient-derived metastatic breast cancer samples, only ER-ve mammospheres were responsive to the ligand Wnt3a. However, the inhibitor DKK1 efficiently inhibited both ER+ve and ER-ve breast cancer but not normal mammosphere formation, suggesting that the Wnt pathway is aberrantly activated in breast cancer mammospheres.

Conclusions

Collectively, these data highlight differential Wnt signalling in breast cancer subtypes and activity in patient-derived metastatic cancer stem-like cells indicating a potential for Wnt-targeted treatment in breast cancers.     

Altered expression of insulin receptor isoforms in breast cancer

Purpose

Insulin-like growth factor (IGF) signaling through human insulin receptor isoform A (IR-A) contributes to tumorigenesis and intrinsic resistance to anti-IGF1R therapy. In the present study, we (a) developed quantitative TaqMan real time-PCR-based assays (qRT-PCR) to measure human insulin receptor isoforms with high specificity, (b) evaluated isoform expression levels in molecularly-defined breast cancer subtypes, and (c) identified the IR-A:IR-B mRNA ratio as a potential biomarker guiding patient stratification for anti-IGF therapies.

Experimental Design

mRNA expression levels of IR-A and IR-B were measured in 42 primary breast cancers and 19 matched adjacent normal tissues with TaqMan qRT-PCR assays. The results were further confirmed in 165 breast cancers. The tumor samples were profiled using whole genome microarrays and subsequently subtyped using the PAM50 breast cancer gene signature. The relationship between the IR-A:IR-B ratio and cancer subtype, as well as markers of proliferation were characterized.

Results

The mRNA expression levels of IR-A in the breast tumors were similar to those observed in the adjacent normal tissues, while the mRNA levels of IR-B were significantly decreased in tumors. The IR-A:IR-B ratio was significantly higher in luminal B breast cancer than in luminal A. Strong concordance between the IR-A:IR-B ratio and the composite Oncotype DX proliferation score was observed for stratifying the latter two breast cancer subtypes.

Conclusions

The reduction in IR-B expression is the key to the altered IR-A:IR-B ratio observed in breast cancer. The IR-A:IR-B ratio may have biomarker utility in guiding a patient stratification strategy for an anti-IGF therapeutic.     

Plasma free amino acid profiling of five types of cancer patients and its application for early detection

Background

Recently, rapid advances have been made in metabolomics-based, easy-to-use early cancer detection methods using blood samples. Among metabolites, profiling of plasma free amino acids (PFAAs) is a promising approach because PFAAs link all organ systems and have important roles in metabolism. Furthermore, PFAA profiles are known to be influenced by specific diseases, including cancers. Therefore, the purpose of the present study was to determine the characteristics of the PFAA profiles in cancer patients and the possibility of using this information for early detection.

Methods and Findings

Plasma samples were collected from approximately 200 patients from multiple institutes, each diagnosed with one of the following five types of cancer: lung, gastric, colorectal, breast, or prostate cancer. Patients were compared to gender- and age- matched controls also used in this study. The PFAA levels were measured using high-performance liquid chromatography (HPLC)–electrospray ionization (ESI)–mass spectrometry (MS). Univariate analysis revealed significant differences in the PFAA profiles between the controls and the patients with any of the five types of cancer listed above, even those with asymptomatic early-stage disease. Furthermore, multivariate analysis clearly discriminated the cancer patients from the controls in terms of the area under the receiver-operator characteristics curve (AUC of ROC >0.75 for each cancer), regardless of cancer stage. Because this study was designed as case-control study, further investigations, including model construction and validation using cohorts with larger sample sizes, are necessary to determine the usefulness of PFAA profiling.

Conclusions

These findings suggest that PFAA profiling has great potential for improving cancer screening and diagnosis and understanding disease pathogenesis. PFAA profiles can also be used to determine various disease diagnoses from a single blood sample, which involves a relatively simple plasma assay and imposes a lower physical burden on subjects when compared to existing screening methods.     

Overexpression and small molecule-triggered downregulation of CIP2A in lung cancer

Background

Lung cancer is the leading cause of cancer deaths worldwide, with a five-year overall survival rate of only 15%. Cancerous inhibitor of PP2A (CIP2A) is a human oncoprotein inhibiting PP2A in many human malignancies. However, whether CIP2A can be a new drug target for lung cancer is largely unclear.

Methodology/Principal Findings

Normal and malignant lung tissues were derived from 60 lung cancer patients from southern China. RT-PCR, Western blotting and immunohistochemistry were used to evaluate the expression of CIP2A. We found that among the 60 patients, CIP2A was undetectable or very low in paratumor normal tissues, but was dramatically elevated in tumor samples in 38 (63.3%) patients. CIP2A overexpression was associated with cigarette smoking. Silencing CIP2A by siRNA inhibited the proliferation and clonogenic activity of lung cancer cells. Intriguingly, we found a natural compound, rabdocoetsin B which is extracted from a Traditional Chinese Medicinal herb Rabdosia coetsa, could induce down-regulation of CIP2A and inactivation of Akt pathway, and inhibit proliferation and induce apoptosis in a variety of lung cancer cells.

Conclusions/Significance

Our findings strongly indicate that CIP2A could be an effective target for lung cancer drug development, and the therapeutic potentials of CIP2A-targeting agents warrant further investigation.     

An integrated quantification method to increase the precision,robustness, and resolution of protein measurement in human plasma samples

Background

Current quantification methods for mass spectrometry (MS)-based proteomics either do not provide sufficient control of variability or are difficult to implement for routine clinical testing.

Results

We present here an integrated quantification (InteQuan) method that better controls pre-analytical and analytical variability than the popular quantification method using stable isotope-labeled standard peptides (SISQuan). We quantified 16 lung cancer biomarker candidates in human plasma samples in three assessment studies, using immunoaffinity depletion coupled with multiple reaction monitoring (MRM) MS. InteQuan outperformed SISQuan in precision in all three studies and tolerated a two-fold difference in sample loading. The three studies lasted over six months and encountered major changes in experimental settings. Nevertheless, plasma proteins in low ng/ml to low μg/ml concentrations were measured with a median technical coefficient of variation (CV) of 11.9% using InteQuan. The corresponding median CV using SISQuan was 15.3% after linear fitting. Furthermore, InteQuan surpassed SISQuan in measuring biological difference among clinical samples and in distinguishing benign versus cancer plasma samples.

Conclusions

We demonstrated that InteQuan is a simple yet robust quantification method for MS-based quantitative proteomics, especially for applications in biomarker research and in routine clinical testing.

Electronic supplementary material

The online version of this article (doi:10.1186/1559-0275-12-3) contains supplementary material, which is available to authorized users.     

Blood peptidome-degradome profile of breast cancer

Background

Cancer invasion and metastasis are closely associated with activities within the degradome; however, little is known about whether these activities can be detected in the blood of cancer patients.

Methodology and Principal Findings

The peptidome-degradome profiles of pooled blood plasma sampled from 15 breast cancer patients (BCP) and age, race, and menopausal status matched control healthy persons (HP) were globally characterized using advanced comprehensive separations combined with tandem Fourier transform mass spectrometry and new data analysis approaches that facilitated top-down peptidomic analysis. The BCP pool displayed 71 degradome protein substrates that encompassed 839 distinct peptidome peptides. In contrast, the HP 50 degradome substrates found encompassed 425 peptides. We find that the ratios of the peptidome peptide relative abundances can vary as much as >4000 fold between BCP and HP. The experimental results also show differential degradation of substrates in the BCP sample in their functional domains, including the proteolytic and inhibitory sites of the plasmin-antiplasmin and thrombin-antithrombin systems, the main chains of the extracellular matrix protection proteins, the excessive degradation of innate immune system key convertases and membrane attack complex components, as well as several other cancer suppressor proteins.

Conclusions

Degradomics-peptidomics profiling of blood plasma is highly sensitive to changes not evidenced by conventional bottom-up proteomics and potentially provides unique signatures of possible diagnostic utility.     

Ultra-deep sequencing reveals the microRNA expression pattern of the human stomach

Background

While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia.

Methodology/Principal Findings

A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue.

Conclusions/Significance

This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.     

Aberrant Mucin5B expression in lung adenocarcinomas detected by iTRAQ labeling quantitative proteomics and immunohistochemistry

Background

Lung cancer is the number one cause of cancer-related deaths in the United States and worldwide. The complex protein changes and/or signature of protein expression in lung cancer, particularly in non-small cell lung cancer (NSCLC) has not been well defined. Although several studies have investigated the protein profile in lung cancers, the knowledge is far from complete. Among early studies, mucin5B (MUC5B) has been suggested to play an important role in the tumor progression. MUC5B is the major gel-forming mucin in the airway. In this study, we investigated the overall protein profile and MUC5B expression in lung adenocarcinomas, the most common type of NSCLCs.

Methods

Lung adenocarcinoma tissue in formalin-fixed paraffin-embedded (FFPE) blocks was collected and microdissected. Peptides from 8 tumors and 8 tumor-matched normal lung tissue were extracted and labeled with 8-channel iTRAQ reagents. The labeled peptides were identified and quantified by LC-MS/MS using an LTQ Orbitrap Velos mass spectrometer. MUC5B expression identified by iTRAQ labeling was further validated using immunohistochemistry (IHC) on tumor tissue microarray (TMA).

Results

A total of 1288 peptides from 210 proteins were identified and quantified in tumor tissues. Twenty-two proteins showed a greater than 1.5-fold differences between tumor and tumor-matched normal lung tissues. Fifteen proteins, including MUC5B, showed significant changes in tumor tissues. The aberrant expression of MUC5B was further identified in 71.1% of lung adenocarcinomas in the TMA.

Discussions

A subset of tumor-associated proteins was differentially expressed in lung adenocarcinomas. The differential expression of MUC5B in lung adenocarcinomas suggests its role as a potential biomarker in the detection of adenocarcinomas.