Integrative single-cell expression and functional studies unravels a sensitization to cytarabine-based chemotherapy through HIF pathway inhibition in AML leukemia stem cells.

Relapse remains a major challenge in the clinical management of acute myeloid leukemia (AML) and is driven by rare therapy-resistant leukemia stem cells (LSCs) that reside in specific bone marrow niches. Hypoxia signaling maintains cells in a quiescent and metabolically relaxed state, desensitizing them to chemotherapy. This suggests the hypothesis that hypoxia contributes to the chemoresistance of AML-LSCs and may represent a therapeutic target to sensitize AML-LSCs to chemotherapy. Here, we identify HIFhigh and HIFlow specific AML subgroups (inv(16)/t(8;21) and MLLr, respectively) and provide a comprehensive single-cell expression atlas of 119,000 AML cells and AML-LSCs in paired diagnostic-relapse samples from these molecular subgroups. The HIF/hypoxia pathway signature is attenuated in AML-LSCs compared with more differentiated AML cells but is more expressed than in healthy hematopoietic cells. Importantly, chemical inhibition of HIF cooperates with standard-of-care chemotherapy to impair AML growth and to substantially eliminate AML-LSCs in vitro and in vivo. These findings support the HIF pathway in the stem cell-driven drug resistance of AML and unravel avenues for combinatorial targeted and chemotherapy-based approaches to specifically eliminate AML-LSCs.

Identification of a Twelve-microRNA Signature with Prognostic Value in Stage II Microsatellite Stable Colon Cancer.

We aimed to identify and validate a set of miRNAs that could serve as a prognostic signature useful to determine the recurrence risk for patients with COAD. Small RNAs from tumors of 100 stage II, untreated, MSS colon cancer patients were sequenced for the discovery step. For this purpose, we built an miRNA score using an elastic net Cox regression model based on the disease-free survival status. Patients were grouped into high or low recurrence risk categories based on the median value of the score. We then validated these results in an independent sample of stage II microsatellite stable tumor tissues, with a hazard ratio of 3.24, (CI95% = 1.05-10.0) and a 10-year area under the receiver operating characteristic curve of 0.67. Functional analysis of the miRNAs present in the signature identified key pathways in cancer progression. In conclusion, the proposed signature of 12 miRNAs can contribute to improving the prediction of disease relapse in patients with stage II MSS colorectal cancer, and might be useful in deciding which patients may benefit from adjuvant chemotherapy.

Tumor suppressor role of RBM22 in prostate cancer acting as a dual-factor regulating alternative splicing and transcription of key oncogenic genes.

Prostate cancer (PCa) is one of the leading causes of cancer-related deaths among men. Consequently, the identification of novel molecular targets for treatment is urgently needed to improve patients' outcomes. Our group recently reported that some elements of the cellular machinery controlling alternative-splicing might be useful as potential novel therapeutic tools against advanced PCa. However, the presence and functional role of RBM22, a key spliceosome component, in PCa remains unknown. Therefore, RBM22 levels were firstly interrogated in 3 human cohorts and 2 preclinical mouse models (TRAMP/Pbsn-Myc). Results were validated in in silico using 2 additional cohorts. Then, functional effects in response to RBM22 overexpression (proliferation, migration, tumorspheres/colonies formation) were tested in PCa models in vitro (LNCaP, 22Rv1, and PC-3 cell-lines) and in vivo (xenograft). High throughput methods (ie, RNA-seq, nCounter PanCancer Pathways Panel) were performed in RBM22 overexpressing cells and xenograft tumors. We found that RBM22 levels were down-regulated (mRNA and protein) in PCa samples, and were inversely associated with key clinical aggressiveness features. Consistently, a gradual reduction of RBM22 from non-tumor to poorly differentiated PCa samples was observed in transgenic models (TRAMP/Pbsn-Myc). Notably, RBM22 overexpression decreased aggressiveness features in vitro, and in vivo. These actions were associated with the splicing dysregulation of numerous genes and to the downregulation of critical upstream regulators of cell-cycle (i.e., CDK1/CCND1/EPAS1). Altogether, our data demonstrate that RBM22 plays a critical pathophysiological role in PCa and invites to suggest that targeting negative regulators of RBM22 expression/activity could represent a novel therapeutic strategy to tackle this disease.

A convergent malignant phenotype in B-cell acute lymphoblastic leukemia involving the splicing factor SRRM1.

A significant proportion of infant B-cell acute lymphoblastic leukemia (B-ALL) patients remains with a dismal prognosis due to yet undetermined mechanisms. We performed a comprehensive multicohort analysis of gene expression, gene fusions, and RNA splicing alterations to uncover molecular signatures potentially linked to the observed poor outcome. We identified 87 fusions with significant allele frequency across patients and shared functional impacts, suggesting common mechanisms across fusions. We further identified a gene expression signature that predicts high risk independently of the gene fusion background and includes the upregulation of the splicing factor SRRM1. Experiments in B-ALL cell lines provided further evidence for the role of SRRM1 on cell survival, proliferation, and invasion. Supplementary analysis revealed that SRRM1 potentially modulates splicing events associated with poor outcomes through protein-protein interactions with other splicing factors. Our findings reveal a potential convergent mechanism of aberrant RNA processing that sustains a malignant phenotype independently of the underlying gene fusion and that could potentially complement current clinical strategies in infant B-ALL.

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.

Colon-specific eQTL analysis to inform on functional SNPs.

BACKGROUND: Genome-wide association studies on colorectal cancer have identified more than 60 susceptibility loci, but for most of them there is no clear knowledge of functionality or the underlying gene responsible for the risk modification. Expression quantitative trail loci (eQTL) may provide functional information for such single nucleotide polymorphisms (SNPs). METHODS: We have performed detailed eQTL analysis specific for colon tissue on a series of 97 colon tumours, their paired adjacent normal mucosa and 47 colon mucosa samples donated by healthy individuals. R package MatrixEQTL was used to search for genome-wide cis-eQTL and trans-eQTL fitting linear models adjusted for age, gender and tissue type to rank transformed expression data. RESULTS: The cis-eQTL analyses has revealed 29,073 SNP-gene associations with permutation-adjusted P-values < 0.01. These correspond to 363 unique genes. The trans-eQTL analysis identified 10,665 significant SNP-gene associations, most of them in the same chromosome, further than 1 Mb of the gene. We provide a web tool to search for specific SNPs or genes. The tool calculates Pearson or Spearman correlation, and allows to select tissue type for analysis. Data and plots can be exported. CONCLUSIONS: This resource should be useful to prioritise SNPs for further functional studies and to identify relevant genes behind identified loci.

Telomere length alterations in microsatellite stable colorectal cancer and association with the immune response.

Telomeres are repetitive sequences (TTAGGG) located at the end of chromosomes. Telomeres progressively shorten with each cell replication cycle, ultimately leading to chromosomal instability and loss of cell viability. Telomere length anomaly appears to be one of the earliest and most prevalent genetic alterations in malignant transformation. Here we aim to estimate telomere length from whole-exome sequencing data in colon tumors and normal colonic mucosa, and to analyze the potential association of telomere length with clinical factors and gene expression in colon cancer. Reads containing at least five repetitions of the telomere sequence (TTAGGG) were extracted from the raw sequences of 42 adjacent normal-tumor paired samples. The number of reads from the tumor sample was normalized to build the Tumor Telomere Length Ratio (TTLR), considered an estimation of telomere length change in the tumor compared to the paired normal tissue. We evaluated the associations between TTLR and clinical factors, gene expression and copy number (CN) aberrations measured in the same tumor samples. Colon tumors showed significantly shorter telomeres than their paired normal samples. No significant association was observed between TTLR and gender, age, tumor location, prognosis, stromal infiltration or molecular subtypes. The functional gene set enrichment analysis showed pathways related to immune response significantly associated with TLLR. By extracting a relative measure of telomere length from whole-exome sequencing data, we have assessed that colon tumor cells predominantly shorten telomeres, and this alteration is associated with expression changes in genes related to immune response and inflammation in tumor cells.

Exome Sequencing Reveals AMER1 as a Frequently Mutated Gene in Colorectal Cancer.

PURPOSE: Somatic mutations occur at early stages of adenoma and accumulate throughout colorectal cancer progression. The aim of this study was to characterize the mutational landscape of stage II tumors and to search for novel recurrent mutations likely implicated in colorectal cancer tumorigenesis. EXPERIMENTAL DESIGN: The exomic DNA of 42 stage II, microsatellite-stable colon tumors and their paired mucosae were sequenced. Other molecular data available in the discovery dataset [gene expression, methylation, and copy number variations (CNV)] were used to further characterize these tumors. Additional datasets comprising 553 colorectal cancer samples were used to validate the discovered mutations. RESULTS: As a result, 4,886 somatic single-nucleotide variants (SNV) were found. Almost all SNVs were private changes, with few mutations shared by more than one tumor, thus revealing tumor-specific mutational landscapes. Nevertheless, these diverse mutations converged into common cellular pathways, such as cell cycle or apoptosis. Among this mutational heterogeneity, variants resulting in early stop codons in the AMER1 (also known as FAM123B or WTX) gene emerged as recurrent mutations in colorectal cancer. Losses of AMER1 by other mechanisms apart from mutations such as methylation and copy number aberrations were also found. Tumors lacking this tumor suppressor gene exhibited a mesenchymal phenotype characterized by inhibition of the canonical Wnt pathway. CONCLUSIONS: In silico and experimental validation in independent datasets confirmed the existence of functional mutations in AMER1 in approximately 10% of analyzed colorectal cancer tumors. Moreover, these tumors exhibited a characteristic phenotype.

miR-1269 promotes metastasis and forms a positive feedback loop with TGF-ß.

As patient survival drops precipitously from early-stage cancers to late-stage and metastatic cancers, microRNAs that promote relapse and metastasis can serve as prognostic and predictive markers as well as therapeutic targets for chemoprevention. Here we show that miR-1269a promotes colorectal cancer (CRC) metastasis and forms a positive feedback loop with TGF-ß signalling. miR-1269a is upregulated in late-stage CRCs, and long-term monitoring of 100 stage II CRC patients revealed that miR-1269a expression in their surgically removed primary tumours is strongly associated with risk of CRC relapse and metastasis. Consistent with clinical observations, miR-1269a significantly increases the ability of CRC cells to invade and metastasize in vivo. TGF-ß activates miR-1269 via Sox4, while miR-1269a enhances TGF-ß signalling by targeting Smad7 and HOXD10, hence forming a positive feedback loop. Our findings suggest that miR-1269a is a potential marker to inform adjuvant chemotherapy decisions for CRC patients and a potential therapeutic target to deter metastasis.

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.