Targeting metastasis-initiating cells through the fatty acid receptor CD36.

The fact that the identity of the cells that initiate metastasis in most human cancers is unknown hampers the development of antimetastatic therapies. Here we describe a subpopulation of CD44bright cells in human oral carcinomas that do not overexpress mesenchymal genes, are slow-cycling, express high levels of the fatty acid receptor CD36 and lipid metabolism genes, and are unique in their ability to initiate metastasis. Palmitic acid or a high-fat diet specifically boosts the metastatic potential of CD36+ metastasis-initiating cells in a CD36-dependent manner. The use of neutralizing antibodies to block CD36 causes almost complete inhibition of metastasis in immunodeficient or immunocompetent orthotopic mouse models of human oral cancer, with no side effects. Clinically, the presence of CD36+ metastasis-initiating cells correlates with a poor prognosis for numerous types of carcinomas, and inhibition of CD36 also impairs metastasis, at least in human melanoma- and breast cancer-derived tumours. Together, our results indicate that metastasis-initiating cells particularly rely on dietary lipids to promote metastasis.

Integrating gene expression and epidemiological data for the discovery of genetic interactions associated with cancer risk.

Dozens of common genetic variants associated with cancer risk have been identified through genome-wide association studies (GWASs). However, these variants only explain a modest fraction of the heritability of disease. The missing heritability has been attributed to several factors, among them the existence of genetic interactions (G × G). Systematic screens for G × G in model organisms have revealed their fundamental influence in complex phenotypes. In this scenario, G × G overlap significantly with other types of gene and/or protein relationships. Here, by integrating predicted G × G from GWAS data and complex- and context-defined gene coexpression profiles, we provide evidence for G × G associated with cancer risk. G × G predicted from a breast cancer GWAS dataset identified significant overlaps [relative enrichments (REs) of 8-36%, empirical P values < 0.05 to 10(-4)] with complex (non-linear) gene coexpression in breast tumors. The use of gene or protein data not specific for breast cancer did not reveal overlaps. According to the predicted G × G, experimental assays demonstrated functional interplay between lipoma-preferred partner and transforming growth factor-ß signaling in the MCF10A non-tumorigenic mammary epithelial cell model. Next, integration of pancreatic tumor gene expression profiles with pancreatic cancer G × G predicted from a GWAS corroborated the observations made for breast cancer risk (REs of 25-59%). The method presented here can potentially support the identification of genetic interactions associated with cancer risk, providing novel mechanistic hypotheses for carcinogenesis.

Identification of candidate susceptibility genes for colorectal cancer through eQTL analysis.

In this study, we aim to identify the genes responsible for colorectal cancer risk behind the loci identified in genome-wide association studies (GWAS). These genes may be candidate targets for developing new strategies for prevention or therapy. We analyzed the association of genotypes for 26 GWAS single nucleotide polymorphisms (SNPs) with the expression of genes within a 2 Mb region (cis-eQTLs). Affymetrix Human Genome U219 expression arrays were used to assess gene expression in two series of samples, one of healthy colonic mucosa (n = 47) and other of normal mucosa adjacent to colon cancer (n = 97, total 144). Paired tumor tissues (n = 97) were also analyzed but did not provide additional findings. Partial Pearson correlation (r), adjusted for sample type, was used for the analysis. We have found Bonferroni-significant cis-eQTLs in three loci: rs3802842 in 11q23.1 associated to C11orf53, COLCA1 (C11orf92) and COLCA2 (C11orf93; r = 0.60); rs7136702 in 12q13.12 associated to DIP2B (r = 0.63) and rs5934683 in Xp22.3 associated to SHROOM2 and GPR143 (r = 0.47). For loci in chromosomes 11 and 12, we have found other SNPs in linkage disequilibrium that are more strongly associated with the expression of the identified genes and are better functional candidates: rs7130173 for 11q23.1 (r = 0.66) and rs61927768 for 12q13.12 (r = 0.86). These SNPs are located in DNA regions that may harbor enhancers or transcription factor binding sites. The analysis of trans-eQTLs has identified additional genes in these loci that may have common regulatory mechanisms as shown by the analysis of protein-protein interaction networks.

Aberrant gene expression in mucosa adjacent to tumor reveals a molecular crosstalk in colon cancer.

BACKGROUND: A colorectal tumor is not an isolated entity growing in a restricted location of the body. The patient's gut environment constitutes the framework where the tumor evolves and this relationship promotes and includes a complex and tight correlation of the tumor with inflammation, blood vessels formation, nutrition, and gut microbiome composition. The tumor influence in the environment could both promote an anti-tumor or a pro-tumor response. METHODS: A set of 98 paired adjacent mucosa and tumor tissues from colorectal cancer (CRC) patients and 50 colon mucosa from healthy donors (246 samples in total) were included in this work. RNA extracted from each sample was hybridized in Affymetrix chips Human Genome U219. Functional relationships between genes were inferred by means of systems biology using both transcriptional regulation networks (ARACNe algorithm) and protein-protein interaction networks (BIANA software). RESULTS: Here we report a transcriptomic analysis revealing a number of genes activated in adjacent mucosa from CRC patients, not activated in mucosa from healthy donors. A functional analysis of these genes suggested that this active reaction of the adjacent mucosa was related to the presence of the tumor. Transcriptional and protein-interaction networks were used to further elucidate this response of normal gut in front of the tumor, revealing a crosstalk between proteins secreted by the tumor and receptors activated in the adjacent colon tissue; and vice versa. Remarkably, Slit family of proteins activated ROBO receptors in tumor whereas tumor-secreted proteins transduced a cellular signal finally activating AP-1 in adjacent tissue. CONCLUSIONS: The systems-level approach provides new insights into the micro-ecology of colorectal tumorogenesis. Disrupting this intricate molecular network of cell-cell communication and pro-inflammatory microenvironment could be a therapeutic target in CRC patients.

COLONOMICS - integrative omics data of one hundred paired normal-tumoral samples from colon cancer patients.

Colonomics is a multi-omics dataset that includes 250 samples: 50 samples from healthy colon mucosa donors and 100 paired samples from colon cancer patients (tumor/adjacent). From these samples, Colonomics project includes data from genotyping, DNA methylation, gene expression, whole exome sequencing and micro-RNAs (miRNAs) expression. It also includes data from copy number variation (CNV) from tumoral samples. In addition, clinical data from all these samples is available. The aims of the project were to explore and integrate these datasets to describe colon cancer at molecular level and to compare normal and tumoral tissues. Also, to improve screening by finding biomarkers for the diagnosis and prognosis of colon cancer. This project has its own website including four browsers allowing users to explore Colonomics datasets. Since generated data could be reuse for the scientific community for exploratory or validation purposes, here we describe omics datasets included in the Colonomics project as well as results from multi-omics layers integration.

Adherence to the Mediterranean diet and risk of coronary heart disease in the Spanish EPIC Cohort Study.

No known cohort study has investigated whether the Mediterranean diet can reduce incident coronary heart disease (CHD) events in a Mediterranean population. This study examined the relation between Mediterranean diet adherence and risk of incident CHD events in the 5 Spanish centers of the European Prospective Investigation into Cancer and Nutrition. Analysis included 41,078 participants aged 29-69 years, recruited in 1992-1996 and followed up until December 2004 (mean follow-up:10.4 years). Confirmed incident fatal and nonfatal CHD events were analyzed according to Mediterranean diet adherence, measured by using an 18-unit relative Mediterranean diet score. A total of 609 participants (79% male) had a fatal or nonfatal confirmed acute myocardial infarction (n = 468) or unstable angina requiring revascularization (n = 141). After stratification by center and age and adjustment for recognized CHD risk factors, high compared with low relative Mediterranean diet score was associated with a significant reduction in CHD risk (hazard ratio = 0.60, 95% confidence interval: 0.47, 0.77). A 1-unit increase in relative Mediterranean diet score was associated with a 6% reduced risk of CHD (95% confidence interval: 0.91, 0.97), with similar risk reductions by sex. Mediterranean diet adherence was associated with a significantly reduced CHD risk in this Mediterranean country, supporting its role in primary prevention of CHD in healthy populations.

Assessing associations between the AURKA-HMMR-TPX2-TUBG1 functional module and breast cancer risk in BRCA1/2 mutation carriers.

While interplay between BRCA1 and AURKA-RHAMM-TPX2-TUBG1 regulates mammary epithelial polarization, common genetic variation in HMMR (gene product RHAMM) may be associated with risk of breast cancer in BRCA1 mutation carriers. Following on these observations, we further assessed the link between the AURKA-HMMR-TPX2-TUBG1 functional module and risk of breast cancer in BRCA1 or BRCA2 mutation carriers. Forty-one single nucleotide polymorphisms (SNPs) were genotyped in 15,252 BRCA1 and 8,211 BRCA2 mutation carriers and subsequently analyzed using a retrospective likelihood approach. The association of HMMR rs299290 with breast cancer risk in BRCA1 mutation carriers was confirmed: per-allele hazard ratio (HR) = 1.10, 95% confidence interval (CI) 1.04-1.15, p = 1.9 x 10(-4) (false discovery rate (FDR)-adjusted p = 0.043). Variation in CSTF1, located next to AURKA, was also found to be associated with breast cancer risk in BRCA2 mutation carriers: rs2426618 per-allele HR = 1.10, 95% CI 1.03-1.16, p = 0.005 (FDR-adjusted p = 0.045). Assessment of pairwise interactions provided suggestions (FDR-adjusted pinteraction values > 0.05) for deviations from the multiplicative model for rs299290 and CSTF1 rs6064391, and rs299290 and TUBG1 rs11649877 in both BRCA1 and BRCA2 mutation carriers. Following these suggestions, the expression of HMMR and AURKA or TUBG1 in sporadic breast tumors was found to potentially interact, influencing patients' survival. Together, the results of this study support the hypothesis of a causative link between altered function of AURKA-HMMR-TPX2-TUBG1 and breast carcinogenesis in BRCA1/2 mutation carriers.

A monotonic and prognostic genomic signature from fibroblasts for colorectal cancer initiation, progression, and metastasis.

UNLABELLED: The differential gene expression patterns between normal colonic fibroblasts (NCF), carcinoma-associated fibroblasts from primary tumors (CAF-PT), and CAFs from hepatic metastasis (CAF-LM) are hypothesized to be useful for predicting relapse in primary tumors. A transcriptomic profile of NCF (n = 9), CAF-PT (n = 14), and CAF-LM (n = 11) was derived. Prediction Analysis of Microarrays (PAM) was used to obtain molecular details for each fibroblast class, and differentially expressed transcripts were used to classify patients according to recurrence status. A number of transcripts (n = 277) were common to all three types of fibroblasts and whose expression level was sequentially deregulated according to the transition: NCF→CAF-PT→CAF-LM. Importantly, the gene signature was able to accurately classify patients with primary tumors according to their prognosis. This capacity was exploited to obtain a refined 19-gene classifier that predicted recurrence with high accuracy in two independent datasets of patients with colorectal cancer and correlates with fibroblast migratory potential. The prognostic power of this genomic signature is strong evidence of the link between the tumor-stroma microenvironment and cancer progression. Furthermore, the 19-gene classifier was able to identify low-risk patients very accurately, which is of particular importance for stage II patients, who would benefit from the omission of chemotherapy, especially T4N0 patients, who are clinically classified as being at high risk. IMPLICATIONS: A defined stromal gene expression signature predicts relapse in patients with colorectal cancer.

Differences between CAFs and their paired NCF from adjacent colonic mucosa reveal functional heterogeneity of CAFs, providing prognostic information.

Little is known about the difference in gene expression between carcinoma-associated fibroblasts (CAFs) and paired normal colonic fibroblasts (NCFs) in colorectal cancer. Paired CAFs and NCFs were isolated from eight primary human colorectal carcinoma specimens. In culture conditions, soluble factors secreted by CAFs in the conditioned media increased clonogenicity and migration of epithelial cancer cells lines to a greater extent than did NCF. In vivo, CAFs were more competent as tumour growth enhancers than paired NCFs when co-inoculated with colorectal cell lines. Gene expression analysis of microarrays of CAF and paired NCF populations enabled us to identify 108 deregulated genes (38 upregulated and 70 downregulated genes). Most of those genes are fibroblast-specific. This has been validated in silico in dataset GSE39396 and by qPCR in selected genes. GSEA analysis revealed a differential transcriptomic profile of CAFs, mainly involving the Wnt signallingsignalling pathway, focal adhesion and cell cycle. Both deregulated genes and biological processes involved depicted a considerable degree of overlap with deregulated genes reported in breast, lung, oesophagus and prostate CAFs. These observations suggest that similar transcriptomic programs may be active in the transition from normal fibroblast in adjacent tissues to CAFs, independently of their anatomic demarcation. Additionally NCF already depicted an activated pattern associated with inflammation. The deregulated genes signature score seemed to correlate with CAF tumour promoter abilities in vitro, suggesting a high degree of heterogeneity between CAFs, and it has also prognostic value in two independent datasets. Further characterization of the roles these biomarkers play in cancer will reveal how CAFs provide cancer cells with a suitable microenvironment and may help in the development of new therapeutic targets for cancer treatment.

A 5-gene classifier from the carcinoma-associated fibroblast transcriptomic profile and clinical outcome in colorectal cancer.

Based on 108 differentially expressed genes between carcinoma-associated fibroblasts (CAFs) and paired normal colonic fibroblasts we recently reported, a 5-gene classifier for relapse prediction in Stage II/III colorectal cancer (CRC ) was developed. Its predictive value was validated in datasets GSE17538, GSE33113 and GSE14095. An additional validation was performed in a metacohort (n=317) and 142 CRC patients by means of RT-PCR. The 5-gene classifier was significantly associated with increased relapse risk and death from CRC across all validation series of Stage II/III patients used. Multivariate Cox regression analyses confirmed the independent prognostic value of the stromal classifier (HR=2.67; P=0.002). Post-test probabilities provided evidence of the suitability of the 5-gene classifier in clinical practice, identifying a subgroup of Stage-II patients who were at high risk of relapse. Moreover, the a priory worst prognosis mesenchymal subtype of tumours can be stratified according to the physiological status of their carcinoma-associated fibroblasts. In conclusion the CAFs-derived 5-gene classifier provides more accurate information about outcome than conventional clinicopathological criteria and it could be useful to take clinical decisions, especially in Stage II. Additionally, the classifier put into relevance the CAF's intratumoral heterogeneity and might contribute to find relevant targets for depleting adequate CAFS subtypes.

Tools for protein-protein interaction network analysis in cancer research.

As cancer is a complex disease, the representation of a malignant cell as a protein-protein interaction network (PPIN) and its subsequent analysis can provide insight into the behaviour of cancer cells and lead to the discovery of new biomarkers. The aim of this review is to help life-science researchers without previous computer programming skills to extract meaningful biological information from such networks, taking advantage of easy-to-use, public bioinformatics tools. It is structured in four parts: the first section describes the pipeline of consecutive steps from network construction to biological hypothesis generation. The second part provides a repository of public, user-friendly tools for network construction, visualisation and analysis. Two different and complementary approaches of network analysis are presented: the topological approach studies the network as a whole by means of structural graph theory, whereas the global approach divides the PPIN into sub-graphs, or modules. In section three, some concepts and tools regarding heterogeneous molecular data integration through a PPIN are described. Finally, the fourth part is an example of how to extract meaningful biological information from a colorectal cancer PPIN using some of the described tools.

Cancer develops, progresses and responds to therapies through restricted perturbation of the protein-protein interaction network.

The products of genes mutated or differentially expressed in cancer tend to occupy central positions within the network of protein-protein interactions, or the interactome network. Integration of different types of gene and protein relationships has considerably increased the understanding of the mechanisms of carcinogenesis, while also enhancing the applicability of expression signatures. In this scenario, however, it remains unknown how cancer develops, progresses and responds to therapies in a potentially controlled manner at the systems level. Here, by applying the concepts of load transfer and cascading failures in power grids, we examine the impact and transmission of cancer-related gene expression changes in the interactome network. Relative to random perturbations, this study reveals topological robustness associated with all cancer conditions. In addition, experimental perturbation of a central cancer node, which consists of over-expression of the α-synuclein (SNCA) protein in MCF7 breast cancer cells, also reveals robustness. Conversely, a search for proteins with an opposite topological impact identifies the autophagy pathway. Mechanistically, the existence of smaller shortest paths among cancer-related proteins appears to be a topological feature that partially contributes to the restricted perturbation of the network. Together, the results of this study suggest that cancer develops, progresses and responds to therapies following controlled, restricted perturbation of the interactome network.

Clinical value of prognosis gene expression signatures in colorectal cancer: a systematic review.

INTRODUCTION: The traditional staging system is inadequate to identify those patients with stage II colorectal cancer (CRC) at high risk of recurrence or with stage III CRC at low risk. A number of gene expression signatures to predict CRC prognosis have been proposed, but none is routinely used in the clinic. The aim of this work was to assess the prediction ability and potential clinical usefulness of these signatures in a series of independent datasets. METHODS: A literature review identified 31 gene expression signatures that used gene expression data to predict prognosis in CRC tissue. The search was based on the PubMed database and was restricted to papers published from January 2004 to December 2011. Eleven CRC gene expression datasets with outcome information were identified and downloaded from public repositories. Random Forest classifier was used to build predictors from the gene lists. Matthews correlation coefficient was chosen as a measure of classification accuracy and its associated p-value was used to assess association with prognosis. For clinical usefulness evaluation, positive and negative post-tests probabilities were computed in stage II and III samples. RESULTS: Five gene signatures showed significant association with prognosis and provided reasonable prediction accuracy in their own training datasets. Nevertheless, all signatures showed low reproducibility in independent data. Stratified analyses by stage or microsatellite instability status showed significant association but limited discrimination ability, especially in stage II tumors. From a clinical perspective, the most predictive signatures showed a minor but significant improvement over the classical staging system. CONCLUSIONS: The published signatures show low prediction accuracy but moderate clinical usefulness. Although gene expression data may inform prognosis, better strategies for signature validation are needed to encourage their widespread use in the clinic.

Exploring the link between germline and somatic genetic alterations in breast carcinogenesis.

Recent genome-wide association studies (GWASs) have identified candidate genes contributing to cancer risk through low-penetrance mutations. Many of these genes were unexpected and, intriguingly, included well-known players in carcinogenesis at the somatic level. To assess the hypothesis of a germline-somatic link in carcinogenesis, we evaluated the distribution of somatic gene labels within the ordered results of a breast cancer risk GWAS. This analysis suggested frequent influence on risk of genetic variation in loci encoding for "driver kinases" (i.e., kinases encoded by genes that showed higher somatic mutation rates than expected by chance and, therefore, whose deregulation may contribute to cancer development and/or progression). Assessment of these predictions using a population-based case-control study in Poland replicated the association for rs3732568 in EPHB1 (odds ratio (OR) = 0.79; 95% confidence interval (CI): 0.63-0.98; P(trend) = 0.031). Analyses by early age at diagnosis and by estrogen receptor α (ERα) tumor status indicated potential associations for rs6852678 in CDKL2 (OR = 0.32, 95% CI: 0.10-1.00; P(recessive) = 0.044) and rs10878640 in DYRK2 (OR = 2.39, 95% CI: 1.32-4.30; P(dominant) = 0.003), and for rs12765929, rs9836340, rs4707795 in BMPR1A, EPHA3 and EPHA7, respectively (ERα tumor status P(interaction)<0.05). The identification of three novel candidates as EPH receptor genes might indicate a link between perturbed compartmentalization of early neoplastic lesions and breast cancer risk and progression. Together, these data may lay the foundations for replication in additional populations and could potentially increase our knowledge of the underlying molecular mechanisms of breast carcinogenesis.

Biological convergence of cancer signatures.

Gene expression profiling has identified cancer prognostic and predictive signatures with superior performance to conventional histopathological or clinical parameters. Consequently, signatures are being incorporated into clinical practice and will soon influence everyday decisions in oncology. However, the slight overlap in the gene identity between signatures for the same cancer type or condition raises questions about their biological and clinical implications. To clarify these issues, better understanding of the molecular properties and possible interactions underlying apparently dissimilar signatures is needed. Here, we evaluated whether the signatures of 24 independent studies are related at the genome, transcriptome or proteome levels. Significant associations were consistently observed across these molecular layers, which suggest the existence of a common cancer cell phenotype. Convergence on cell proliferation and death supports the pivotal involvement of these processes in prognosis, metastasis and treatment response. In addition, functional and molecular associations were identified with the immune response in different cancer types and conditions that complement the contribution of cell proliferation and death. Examination of additional, independent, cancer datasets corroborated our observations. This study proposes a comprehensive strategy for interpreting cancer signatures that reveals common design principles and systems-level properties.

Biological processes, properties and molecular wiring diagrams of candidate low-penetrance breast cancer susceptibility genes.

BACKGROUND: Recent advances in whole-genome association studies (WGASs) for human cancer risk are beginning to provide the part lists of low-penetrance susceptibility genes. However, statistical analysis in these studies is complicated by the vast number of genetic variants examined and the weak effects observed, as a result of which constraints must be incorporated into the study design and analytical approach. In this scenario, biological attributes beyond the adjusted statistics generally receive little attention and, more importantly, the fundamental biological characteristics of low-penetrance susceptibility genes have yet to be determined. METHODS: We applied an integrative approach for identifying candidate low-penetrance breast cancer susceptibility genes, their characteristics and molecular networks through the analysis of diverse sources of biological evidence. RESULTS: First, examination of the distribution of Gene Ontology terms in ordered WGAS results identified asymmetrical distribution of Cell Communication and Cell Death processes linked to risk. Second, analysis of 11 different types of molecular or functional relationships in genomic and proteomic data sets defined the "omic" properties of candidate genes: i/ differential expression in tumors relative to normal tissue; ii/ somatic genomic copy number changes correlating with gene expression levels; iii/ differentially expressed across age at diagnosis; and iv/ expression changes after BRCA1 perturbation. Finally, network modeling of the effects of variants on germline gene expression showed higher connectivity than expected by chance between novel candidates and with known susceptibility genes, which supports functional relationships and provides mechanistic hypotheses of risk. CONCLUSION: This study proposes that cell communication and cell death are major biological processes perturbed in risk of breast cancer conferred by low-penetrance variants, and defines the common omic properties, molecular interactions and possible functional effects of candidate genes and proteins.

Changes in food supply in Mediterranean countries from 1961 to 2001.

OBJECTIVE: To describe geographical differences and time trends in the supply of the most important food components of the traditional Mediterranean diet. DESIGN: Food supply data collected from national food balance sheets for the period 1961-2001. SETTING: Selected Mediterranean countries: Spain, Italy, France, Greece, Algeria, Morocco, Tunisia and Turkey. RESULTS: Differences of almost 30-fold and five-fold were found in the supply of olive oil and fruits and vegetables, respectively, among the Mediterranean countries studied during the 1960s. A favourable increasing trend for the supply of fruit and vegetables was observed in most Mediterranean countries. However, an increase in the supply of meats and dairy products and a decrease in the supply of cereals and wine were observed in European Mediterranean countries from 1961 until 2001. Only in African and Asiatic Mediterranean countries were cereals the base of food supply. During the 1990s, Greece's food supply pattern was closest to the traditional Mediterranean diet, while Italy and Spain maintained a high availability of fruits, vegetables and olive oil, but were losing the other typical components. Among African and Asiatic Mediterranean countries, only Turkey presented a traditional Mediterranean dietary pattern except with respect to olive oil, the supply of which was very low. France showed a Western dietary pattern, with a high supply of animal products and a low supply of olive oil. CONCLUSIONS: Dietary supplies in the Mediterranean area were quite heterogeneous in the 1960s and have experienced a process of Westernization, especially in European Mediterranean countries.