Differential expression of miRNAs in pancreatobiliary type of periampullary adenocarcinoma and its associated stroma

Transcriptomic-Based Microenvironment Classification Reveals Precision Medicine Strategies for Pancreatic Ductal Adenocarcinoma

BACKGROUND & AIMS

The complex tumor microenvironment (TME) of pancreatic ductal adenocarcinoma (PDAC) has hindered the development of reliable predictive biomarkers for targeted therapy and immunomodulatory strategies. A comprehensive characterization of the TME is necessary to advance precision therapeutics in PDAC.

METHODS

A transcriptomic profiling platform for TME classification based on functional gene signatures was applied to 14 publicly available PDAC datasets (n = 1657) and validated in a clinically annotated independent cohort of patients with PDAC (n = 79). Four distinct subtypes were identified using unsupervised clustering and assessed to evaluate predictive and prognostic utility.

RESULTS

TME classification using transcriptomic profiling identified 4 biologically distinct subtypes based on their TME immune composition: immune enriched (IE); immune enriched, fibrotic (IE/F); fibrotic (F); and immune depleted (D). The IE and IE/F subtypes demonstrated a more favorable prognosis and potential for response to immunotherapy compared with the F and D subtypes. Most lung metastases and liver metastases were subtypes IE and D, respectively, indicating the role of clonal phenotype and immune milieu in developing personalized therapeutic strategies. In addition, distinct TMEs with potential therapeutic implications were identified in treatment-naive primary tumors compared with tumors that underwent neoadjuvant therapy.

CONCLUSIONS

This novel approach defines a distinct subgroup of PADC patients that may benefit from immunotherapeutic strategies based on their TME subtype and provides a framework to select patients for prospective clinical trials investigating precision immunotherapy in PDAC. Further, the predictive utility and real-world clinical applicability espoused by this transcriptomic-based TME classification approach will accelerate the advancement of precision medicine in PDAC.

HtrA2 Independently Predicts Poor Prognosis and Correlates with Immune Cell Infiltration in Hepatocellular Carcinoma

High-temperature requirement protein A2 (HtrA2), a mitochondrial protein, is related to apoptosis regulation. However, the role of HtrA2 in hepatocellular carcinoma (HCC) remains unclear. In the present study, we explored the prognostic value and expression pattern of HtrA2 in HCC and confirmed its independent value for predicting outcomes via Cox analyses. LinkedOmics and GEPIA2 were used to construct the coexpression and functional networks of HtrA2. Additionally, the data obtained from TCGA was analyzed to investigate the relationship between the infiltration of immune cells and HtrA2 mRNA expression. Finally, the expression pattern of HtrA2 in HCC was confirmed by wet-lab experiments. The results showed high HtrA2 expression (P < 0.001) presented in tumor tissues in TCGA-HCC. Moreover, high HtrA2 expression was confirmed to be associated with poor HCC patient survival (P < 0.05). HtrA2 has also been recognized as an essential risk factor for overall survival (P=0.01, HR = 1.654, 95% CI 1.128-2.425), disease-specific survival (P=0.004, HR = 2.204, 95% CI 1.294-3.753), and progression-free interval (P=0.007, HR = 1.637, 95% CI 1.145-2.341) of HCC. HCC patients with low HtrA2 methylation had worse overall survival than patients with high methylation (P=0.0019). Functional network analysis suggests that HtrA2 regulates mitochondrial homeostasis through pathways involving multiple microRNAs and transcription factors in HCC. In addition, HtrA2 expression correlated with infiltrating levels of multiple immune cell populations. At last, increased expression of HtrA2 in HCC was confirmed using wet-lab experiments. Our study provides evidence that the upregulation of HtrA2 in HCC is an independent predictor of prognosis. Our results provide the foundation for further study on the roles of HtrA2 in HCC tumorigenesis.

PLK2 Single Nucleotide Variant in Gastric Cancer Patients Affects miR-23b-5p Binding

PURPOSE

Chromosomal instability is a hallmark of gastric cancer (GC). It can be driven by single nucleotide variants (SNVs) in cell cycle genes. We investigated the associations between SNVs in candidate genes, PLK2, PLK3, and ATM, and GC risk and clinicopathological features.

MATERIALS AND METHODS

The genotyping study included 542 patients with GC and healthy controls. Generalized linear models were used for the risk and clinicopathological association analyses. Survival analysis was performed using the Kaplan-Meier method. The binding of candidate miRs was analyzed using a luciferase reporter assay.

RESULTS

The PLK2 Crs15009-Crs963615 haplotype was under-represented in the GC group compared to that in the control group (Pcorr=0.050). Male patients with the PLK2 rs963615 CT genotype had a lower risk of GC, whereas female patients had a higher risk (P=0.023; P=0.026). The PLK2 rs963615 CT genotype was associated with the absence of vascular invasion (P=0.012). The PLK3 rs12404160 AA genotype was associated with a higher risk of GC in the male population (P=0.015). The ATM Trs228589-Ars189037-Grs4585 haplotype was associated with a higher risk of GC (P<0.001). The ATM rs228589, rs189037, and rs4585 genotypes TA+AA, AG+GG, and TG+GG were associated with the absence of perineural invasion (P=0.034). In vitro analysis showed that the cancer-associated miR-23b-5p mimic specifically bound to the PLK2 rs15009 G allele (P=0.0097). Moreover, low miR-23b expression predicted longer 10-year survival (P=0.0066) in patients with GC.

CONCLUSIONS

PLK2, PLK3, and ATM SNVs could potentially be helpful for the prediction of GC risk and clinicopathological features. PLK2 rs15009 affects the binding of miR-23b-5p. MiR-23b-5p expression status could serve as a prognostic marker for survival in patients with GC.

Dynamic contract-enhanced CT-based radiomics for differentiation of pancreatobiliary-type and intestinal-type periampullary carcinomas

AIM

To investigate whether computed tomography (CT) radiomics can differentiate pancreatobiliary-type from intestinal-type periampullary carcinomas.

MATERIALS AND METHODS

CT radiomics of 96 patients (54 pancreatobiliary type and 42 intestinal type) with surgically confirmed periampullary carcinoma were assessed retrospectively. Volumes of interest (VOIs) were delineated manually. Radiomic features were extracted from preoperative CT images. A single-phase model and combined-phase model were constructed. Five-fold cross-validation and five machine-learning algorithms were utilised for model construction. The diagnostic performance of the models was evaluated by receiver operating characteristic (ROC) curves, and indicators included area under the curve (AUC), accuracy, sensitivity, specificity, and precision. ROC curves were compared using DeLong's test.

RESULTS

A total of 788 features were extracted on each phase. After feature selection using least absolute shrinkage and selection operator (LASSO) algorithm, the number of selected optimal feature was 18 (plain scan), nine (arterial phase), two (venous phase), 23 (delayed phase), 15 (three enhanced phases), and 29 (all phases), respectively. For the single-phase model, the delayed-phase model using the logistic regression (LR) algorithm showed the best prediction performance with AUC, accuracy, sensitivity, specificity, and precision of 0.89, 0.83, 0.80, 0.88, and 0.93, respectively. Two combined-phase models showed better results than the single-phase models. The model of all phases using the LR algorithm showed the best prediction performance with AUC, accuracy, sensitivity, specificity, and precision of 0.96, 0.88, 0.90, 0.93, and 0.92, respectively.

CONCLUSION

Radiomic models based on preoperative CT images can differentiate pancreatobiliary-type from intestinal-type periampullary carcinomas, in particular, the model of all phases using the LR algorithm.

miR-194-5p down-regulates tumor cell PD-L1 expression and promotes anti-tumor immunity in pancreatic cancer

Pancreatic cancer is a highly malignant cancer of the digestive tract. Studies have shown that in some types of cancer, a high level of microRNA-194-5p (miR-194-5p) is beneficial for controlling tumor progression, while in other cancers it plays a completely opposite role. However, how miR-194-5p affects anti-tumor immunity of pancreatic cancer remains unclear. In this study, we found that high expression of miR-194-5p in human pancreatic cancer patients is associated with a better survival rate, while increased expression of programmed cell death ligand 1 (PD-L1) in human pancreatic cancer patients is associated with a worse survival rate. In pancreatic cancer, the expression level of PD-L1 is negatively correlated with the expression level of miR-194-5p, and we identified that PD-L1 was target gene of miR-194-5p. In addition, we found that overexpression of miR-194-5p inhibited the migration, invasion and proliferation of pancreatic cancer cells in vitro. The orthotopic mouse model of pancreatic cancer shown that miR-194-5p suppressed the progression of pancreatic cancer, promoted the infiltration of CD8⁺ T cells in tumor immune microenvironments, and enhanced the IFN-γ production of CD8⁺ T cells. Consistently, the co-culture experiments showed that overexpression of miR-194-5p in tumor cell enhanced IFN-γ production by CD8⁺ T cells. In conclusion, miR-194-5p may serve as a novel immunotherapeutic target for pancreatic ductal adenocarcinoma (PDAC) by inhibiting the expression of PD-L1, and play important roles in inhibiting the progression of pancreatic cancer and boosting the anti-tumor effect of CD8⁺ T cells.

An Integrated Data Analysis of mRNA, miRNA and Signaling Pathways in Pancreatic Cancer

Although many genes and miRNAs have been reported for various cancers, pancreatic cancer's specific genes or miRNAs have not been studied precisely yet. Therefore, we have analyzed the gene and miRNA expression profile of pancreatic cancer data in the gene expression omnibus (GEO) database. The microarray-derived miRNAs and mRNAs were annotated by gene ontology (GO) and signaling pathway analysis. We also recognized mRNAs that were targeted by miRNA through the mirDIP database. An integrated analysis of the microarray revealed that only 6 out of 43 common miRNAs had significant differences in their expression profiles between the tumor and normal groups (P value < 0.05 and |log Fold Changes (logFC)|> 1). The hsa-miR-210 had upregulation, whereas hsa-miR-375, hsa-miR-216a, hsa-miR-217, hsa-miR-216b and hsa-miR-634 had downregulation in pancreatic cancer (PC). The analysis results also revealed 109 common mRNAs by microarray and mirDIP 4.1 databases. Pathway analysis showed that amoebiasis, axon guidance, PI3K-Akt signaling pathway, absorption and focal adhesion, adherens junction, platelet activation, protein digestion, human papillomavirus infection, extracellular matrix (ECM) receptor interaction, and riboflavin metabolism played important roles in pancreatic cancer. GO analysis revealed the significant enrichment in the three terms of biological process, cellular component, and molecular function, which were identified as the most important processes associated strongly with pancreatic cancer. In conclusion, DTL, CDH11, COL5A1, ITGA2, KIF14, SMC4, VCAN, hsa-mir-210, hsa-mir-217, hsa-mir-216a, hsa-mir-216b, hsa-mir-375 and hsa-mir-634 can be reported as the novel diagnostic or even therapeutic markers for the future studies. Also, the hsa-mir-107 and hsa-mir-125a-5p with COL5A1, CDH11 and TGFBR1 genes can be introduced as major miRNA and genes on the miRNA-drug-mRNA network. The new regulatory network created in our study could give a deeper knowledge of the pancreatic cancer.

Integrated analysis of circular RNA-associated ceRNA network in pancreatic ductal adenocarcinoma

Circular RNAs (circRNAs) have displayed dysregulated expression in several types of cancer. However, the functions of the majority of circRNAs in pancreatic ductal adenocarcinoma (PDAC) remain unknown. The present study aimed to investigate the expression, functions and molecular mechanisms of circRNAs in PDAC. The circRNAs, mRNAs and the microRNA (miRNAs) expression profiles were obtained from three Gene Expression Omnibus microarray datasets, and a circRNA-miRNA-mRNA and circRNA-miRNA-hubgene network was established. The interactions between proteins were analyzed using the Search Tool for the Retrieval of Interacting Genes/Proteins database, and hubgenes were identified using the MCODE plugin. A total of eight differentially expressed circRNAs (DEcircRNAs), 44 differentially expressed miRNAs (DEmiRNAs), and 2,052 differentially expressed mRNAs (DEmRNAs) were identified. The present study successfully constructed a circRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network based on four circRNAs, six miRNAs and 111 mRNAs in PDAC. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathways analyses indicated that DEmRNAs may participate in the pathogenesis and progression of PDAC. The protein-protein interaction network and module analysis identified six hubgenes (THBS1, FN1, TIMP3, TGFB2, ITGA1 and ITGA3). Furthermore, the circRNA-miRNA-hubgene regulatory modules were constructed based on the three DEcircRNAs, one DEmiRNAs and five DEmRNAs. In conclusion, the results of the present study improve the current understanding of the pathogenesis of PDAC.

Differential expression of the miR-17-92 cluster and miR-17 family in breast cancer according to tumor type; results from the Norwegian Women and Cancer (NOWAC) study

Background

MicroRNAs (miRNAs) are promising biomarkers due to their structural stability and distinct expression profile in various cancers. We wanted to explore the miRNA expression in benign breast tissue and breast cancer subgroups in the Norwegian Women and Cancer study.

Methods

Specimens and histopathological data from study participants in Northern Norway diagnosed with breast cancer, and benign tissue from breast reduction surgery were collected. Main molecular subtypes were based on surrogate markers; luminal A (ER+ and/or PR+, HER2− and Ki67 ≤ 30%), luminal B (ER+ and/or PR+, HER2− and Ki67 > 30% or ER+ and/or PR+ and HER2+), HER2 positive (ER− and PR− and HER2+) and triple-negative (ER−, PR− and HER2−). RNA was extracted from formalin-fixed paraffin-embedded (FFPE) tissue, and miRNAs were successfully analyzed in 102 cancers and 36 benign controls using the 7th generation miRCURY LNA microarray containing probes targeting all human miRNAs as annotated in miRBASE version 19.0. Validation with RT-qPCR was performed.

Results

On average, 450 miRNAs were detected in each sample, and 304 miRNAs were significantly different between malignant and benign tissue. Subgroup analyses of cancer cases revealed 23 miRNAs significantly different between ER+ and ER− tumors, and 47 miRNAs different between tumors stratified according to grade. Significantly higher levels were found in high grade tumors for miR-17-5p (p = 0.006), miR-20a-5p (p = 0.007), miR-106b-5p (p = 0.007), miR-93-5p (p = 0.007) and miR-25-3p (p = 0.015) from the paralogous clusters miR-17-92 and miR-106b-25. Expression of miR-17-5p (p = 0.0029), miR-20a-5p (p = 0.0021), miR-92a-3p (p = 0.011) and miR-106b-5p (p = 0.021) was significantly higher in triple-negative tumors compared to the rest, and miR-17-5p and miR-20a-5p were significantly lower in luminal A tumors.

Conclusions

miRNA expression profiles were significantly different between malignant and benign tissue and between cancer subgroups according to ER− status, grade and molecular subtype. miRNAs in the miR-17-92 cluster and miR-17 family were overexpressed in high grade and triple-negative tumors associated with aggressive behavior. The expression and functional role of these miRNAs should be further studied in breast cancer to explore their potential as biomarkers in diagnostic pathology and clinical oncology.

Differences in microRNA expression between melanoma and healthy adjacent skin

Background

The tumor microenvironment is composed of cancer-associated fibroblasts, tumor-associated macrophages, endothelial cells, immune cells, signaling molecules and extracellular matrix structures, which closelycommunicate with the tumor via multiple mechanisms. MicroRNAs are paracrine regulators that provide a direct interaction between the microenvironment and cancer cells. In the presentstudy, we aimed to identify the microRNA expression profile in melanoma compared with thatin healthy adjacent skin, with a further assessment of altered microRNA signaling pathways and target genes.

Methods

Formalin-fixed paraffin-embedded (FFPE) melanoma tissue samples were separated by dissection into tumor and surrounding health tissue fragments. MicroRNA expression profiles were obtained by microarray using Gene Atlas Microarray System (Affymetrix, California, USA). To confirm microarray results real-time PCR was carried out. Bioinformatic analysis was performed using the DIANA-miRPath v.3.0 database. Target genes for miR-146a-5p were determined using three algorithms: TargetScan 7.0, miRWalk 2.0 and miRTarBase v.4.5.

Results

A microarray profiling revealed 143 microRNAs asdifferent in tumor versus adjacent tissues. Expression level of hsa-miR-146a-5p showedto be higher in melanoma cells as compared to thehealthy adjacent skin. The bioinformatic study has determined several signaling cascades associated with miR-146a-5p:Toll-like receptor pathway, NF-κB pathway, ErB pathway, and measles signaling pathway. The 38 target genes have been shown for miR-146a-5p of which NRAS gene is known asone of the most frequent mutated in melanoma.

Conclusions

Elucidation of the role of miR-146-a-5p in complex interactions between the tumor and the cells of healthy adjacent skin is necessary for our understanding of the mechanisms oftumor progression. Significant differences found between cancer cells and adjacent tissues in microRNA expression profile corresponding to divergent mRNA/protein levels in these structures should be taken into account when tumor samples characterization estimatedby high-throughput methods.

Electronic supplementary material

The online version of this article (10.1186/s12895-018-0081-1) contains supplementary material, which is available to authorized users.

A systems medicine approach for finding target proteins affecting treatment outcomes in patients with non-Hodgkin lymphoma

Autoantibody profiling with a systems medicine approach can help identify critical dysregulated signaling pathways (SPs) in cancers. In this way, immunoglobulins G (IgG) purified from the serum samples of 92 healthy controls, 10 pre-treated (PR) non-Hodgkin lymphoma (NHL) patients, and 20 NHL patients who underwent chemotherapy (PS) were screened with a phage-displayed random peptide library. Protein-protein interaction networks of the PR and PS groups were analyzed and visualized by Gephi. The results indicated AXIN2, SENP2, TOP2A, FZD6, NLK, HDAC2, HDAC1, and EHMT2, in addition to CAMK2A, PLCG1, PLCG2, GRM5, GRIN2B, GRIN2D, CACNA2D3, and SPTAN1 as hubs in 11 and 7 modules of PR and PS networks, respectively. PR- and PS-specific hubs were evaluated in the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome databases. The PR-specific hubs were involved in Wnt SP, signaling by Notch1 in cancer, telomere maintenance, and transcriptional misregulation. In contrast, glutamate receptor SP, Fc receptor-related pathways, growth factors-related SPs, and Wnt SP were statistically significant enriched pathways, based on the pathway analysis of PS hubs. The results revealed that the most PR-specific proteins were associated with events involved in tumor development, while chemotherapy in the PS group was associated with side effects of drugs and/or cancer recurrence. As the findings demonstrated, PR- and PS-specific proteins in this study can be promising therapeutic targets in future studies.

No Functional Role for microRNA-342 in a Mouse Model of Pancreatic Acinar Carcinoma

The intronic microRNA (miR)-342 has been proposed as a potent tumor-suppressor gene. miR-342 is found to be downregulated or epigenetically silenced in multiple different tumor sites, and this loss of expression permits the upregulation of several key oncogenic pathways. In several different cell lines, lower miR-342 expression results in enhanced proliferation and metastasis potential, both in vitro and in xenogenic transplant conditions. Here, we sought to determine the function of miR-342 in an in vivo spontaneous cancer model, using the Ela1-TAg transgenic model of pancreatic acinar carcinoma. Through longitudinal magnetic resonance imaging monitoring of Ela1-TAg transgenic mice, either wild-type or knockout for miR-342, we found no role for miR-342 in the development, growth rate, or pathogenicity of pancreatic acinar carcinoma. These results indicate the importance of assessing miR function in the complex physiology of in vivo model systems and indicate that further functional testing of miR-342 is required before concluding it is a bona fide tumor-suppressor-miR.

Tissue MicroRNA profiles as diagnostic and prognostic biomarkers in patients with resectable pancreatic ductal adenocarcinoma and periampullary cancers

Background

The aim of this study was to validate previously described diagnostic and prognostic microRNA expression profiles in tissue samples from patients with pancreatic cancer and other periampullary cancers.

Methods

Expression of 46 selected microRNAs was studied in formalin-fixed paraffin-embedded tissue from patients with resected pancreatic ductal adenocarcinoma (n = 165), ampullary cancer (n=59), duodenal cancer (n = 6), distal common bile duct cancer (n = 21), and gastric cancer (n = 20); chronic pancreatitis (n = 39); and normal pancreas (n = 35). The microRNAs were analyzed by PCR using the Fluidigm platform.

Results

Twenty-two microRNAs were significantly differently expressed in patients with pancreatic cancer when compared to healthy controls and chronic pancreatitis patients; 17 miRNAs were upregulated (miR-21-5p, −23a-3p, −31-5p, −34c-5p, −93-3p, −135b-3p, −155-5p, −186-5p, −196b-5p, −203, −205-5p, −210, −222-3p, −451, −492, −614, and miR-622) and 5 were downregulated (miR-122-5p, −130b-3p, −216b, −217, and miR-375). MicroRNAs were grouped into diagnostic indices of varying complexity. Ten microRNAs associated with prognosis were identified (let-7 g, miR-29a-5p, −34a-5p, −125a-3p, −146a-5p, −187, −205-5p, −212-3p, −222-5p, and miR-450b-5p). Prognostic indices based on differences in expression of 2 different microRNAs were constructed for pancreatic and ampullary cancer combined and separately (30, 5, and 21 indices).

Conclusion

The study confirms that pancreatic cancer tissue has a microRNA expression profile that is different from that of other periampullary cancers, chronic pancreatitis, and normal pancreas. We identified prognostic microRNAs and microRNA indices that were associated with shorter overall survival in patients with radically resected pancreatic cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s40364-017-0087-6) contains supplementary material, which is available to authorized users.

Updated therapeutic outcome for patients with periampullary and pancreatic cancer related to recent translational research

Chemotherapy with improved effect in patients with metastatic pancreatic cancer has recently been established, launching a new era for patients with this very aggressive disease. FOLFIRINOX and gemcitabine plus nab-paclitaxel are different regimens, both capable of stabilizing the disease, thus increasing the number of patients who can reach second line and even third line of treatment. Concurrently, new windows of opportunity open for nutritional support and other therapeutic interventions, improving quality of life. Also pancreatic surgery has changed significantly during the latest years. Extended operations, including vascular/multivisceral resections are frequently performed in specialized centers, pushing borders of resectability. Potentially curative treatment including neoadjuvant and adjuvant chemotherapy is offered new patient groups. Translational research is the basis for the essential understanding of the ongoing development. Even thou biomarkers for clinical management of patients with periampullary tumors have almost been lacking, biomarker driven trials are now in progress. New insight is constantly made available for clinicians; one recent example is selection of patients for gemcitabine treatment based on the expression level of the human equilibrium nucleoside transporter 1. An example of new diagnostic tools is identification of early pancreatic cancer patients by a three-biomarker panel in urine: The proteins lymphatic vessel endothelial hyaluronan receptor 1, regenerating gene 1 alpha and translation elongation factor 1 alpha. Requirement of treatment guideline revisions is intensifying, as combined chemotherapy regimens result in unexpected advantages. The European Study Group for Pancreatic Cancer 4 trial outcome is an illustration: Addition of capecitabine in the adjuvant setting improved overall survival more than expected from the effect in advanced disease. Rapid implementation of new treatment options is mandatory when progress finally extends to patients with this serious disease.

Differential expression of miRNAs in pancreatobiliary type of periampullary adenocarcinoma and its associated stroma

Periampullary adenocarcinomas can be of two histological subtypes, intestinal or pancreatobiliary. The latter is more frequent and aggressive, and characterized by a prominent desmoplastic stroma, which is tightly related to the biology of the cancer, including its poor response to chemotherapy. Whereas miRNAs are known to regulate various cellular processes and interactions between cells, their exact role in periampullary carcinoma remains to be characterized, especially with respect to the prominent stromal component of pancreatobiliary type cancers. The present study aimed at elucidating this role by miRNA expression profiling of the carcinomatous and stromal component in twenty periampullary adenocarcinomas of pancreatobiliary type. miRNA expression profiles were compared between carcinoma cells, stromal cells and normal tissue samples. A total of 43 miRNAs were found to be differentially expressed between carcinoma and stroma of which 11 belong to three miRNA families (miR-17, miR-15 and miR-515). The levels of expression of miRNAs miR-17, miR-20a, miR-20b, miR-223, miR-10b, miR-2964a and miR-342 were observed to be higher and miR-519e to be lower in the stromal component compared to the carcinomatous and normal components. They follow a trend where expression in stroma is highest followed by carcinoma and then normal tissue. Pathway analysis revealed that pathways regulating tumor-stroma interactions such as ECM interaction remodeling, epithelial-mesenchymal transition, focal adhesion pathway, TGF-beta, MAPK signaling, axon guidance and endocytosis were differently regulated. The miRNA-mRNA mediated interactions between carcinoma and stromal cells add new knowledge regarding tumor-stroma interactions.