BCL9/9L proteins enhance the transcriptional output of the β-catenin/TCF transcriptional complex and contribute critically to upholding the high WNT signaling level required for stemness maintenance in the intestinal epithelium. Here we show that a BCL9/9L-dependent gene signature derived from independent mouse colorectal cancer (CRC) models unprecedentedly separates patient subgroups with regard to progression free and overall survival. We found that this effect was by and large attributable to stemness related gene sets. Remarkably, this signature proved associated with recently described poor prognosis CRC subtypes exhibiting high stemness and/or epithelial-to-mesenchymal transition (EMT) traits. Consistent with the notion that high WNT signaling is required for stemness maintenance, ablating Bcl9/9l-β-catenin in murine oncogenic intestinal organoids provoked their differentiation and completely abrogated their tumorigenicity, while not affecting their proliferation. Therapeutic strategies aimed at targeting WNT responses may be limited by intestinal toxicity. Our findings suggest that attenuating WNT signaling to an extent that affects stemness maintenance without disturbing intestinal renewal might be well tolerated and prove sufficient to reduce CRC recurrence and dramatically improve disease outcome.

While T cell-based immunotherapies are steadily improving, there are still many patients who progress, despite T cell-infiltrated tumors. Emerging evidence suggests that T cells themselves may provoke immune escape of cancer cells. Here, we describe a well-controlled co-culture system for studying the dynamic T cell - cancer cell interplay, using human melanoma as a model. We explain starting material, controls, and culture parameters to establish reproducible and comparable cultures with highly heterogeneous tumor cells. Low passage melanoma cell lines and melanoma-specific CD8+ T cell clones generated from patient blood were cultured together for up to 3 days. Living melanoma cells were isolated from the co-culture system by fluorescence-activated cell sorting. We demonstrate that the characterization of isolated melanoma cells is feasible using flow cytometry for protein expression analysis as well as an Agilent whole human genome microarray and the NanoString technology for differential gene expression analysis. In addition, we identify five genes (ALG12, GUSB, RPLP0, KRBA2, and ADAT2) that are stably expressed in melanoma cells independent of the presence of T cells or the T cell-derived cytokines IFNγ and TNFα. These genes are essential for correct normalization of gene expression data by NanoString. Further to the characterization of melanoma cells after exposure to CTLs, this experimental system might be suitable to answer a series of questions, including how the affinity of CTLs for their target antigen influences the melanoma cell response and whether CTL-induced gene expression changes in melanoma cells are reversible. Taken together, our human T cell - melanoma cell culture system is well suited to characterize immune-related mechanisms in cancer cells.

Microsatellite instability (MSI) occurs in 10-20% of colorectal tumours and is associated with good prognosis. Here we describe the development and validation of a genomic signature that identifies colorectal cancer patients with MSI caused by DNA mismatch repair deficiency with high accuracy. Microsatellite status for 276 stage II and III colorectal tumours has been determined. Full-genome expression data was used to identify genes that correlate with MSI status. A subset of these samples (n = 73) had sequencing data for 615 genes available. An MSI gene signature of 64 genes was developed and validated in two independent validation sets: the first consisting of frozen samples from 132 stage II patients; and the second consisting of FFPE samples from the PETACC-3 trial (n = 625). The 64-gene MSI signature identified MSI patients in the first validation set with a sensitivity of 90.3% and an overall accuracy of 84.8%, with an AUC of 0.942 (95% CI, 0.888-0.975). In the second validation, the signature also showed excellent performance, with a sensitivity 94.3% and an overall accuracy of 90.6%, with an AUC of 0.965 (95% CI, 0.943-0.988). Besides correct identification of MSI patients, the gene signature identified a group of MSI-like patients that were MSS by standard assessment but MSI by signature assessment. The MSI-signature could be linked to a deficient MMR phenotype, as both MSI and MSI-like patients showed a high mutation frequency (8.2% and 6.4% of 615 genes assayed, respectively) as compared to patients classified as MSS (1.6% mutation frequency). The MSI signature showed prognostic power in stage II patients (n = 215) with a hazard ratio of 0.252 (p = 0.0145). Patients with an MSI-like phenotype had also an improved survival when compared to MSS patients. The MSI signature was translated to a diagnostic microarray and technically and clinically validated in FFPE and frozen samples.

Cancer omics data are exponentially created and associated with clinical variables, and important findings can be extracted based on bioinformatics approaches which can then be experimentally validated. Many of these findings are related to a specific class of non-coding RNA molecules called microRNAs (miRNAs) (post-transcriptional regulators of mRNA expression). The related research field is quite heterogeneous and bioinformaticians, clinicians, statisticians and biologists, as well as data miners and engineers collaborate to cure stored data and on new impulses coming from the output of the latest Next Generation Sequencing technologies. Here we review the main research findings on miRNA of the first 10 years in colon cancer research with an emphasis on possible uses in clinical practice. This review intends to provide a road map in the jungle of publications of miRNA in colorectal cancer, focusing on data availability and new ways to generate biologically relevant information out of these huge amounts of data.

Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) experiments are widely used to determine, within entire genomes, the occupancy sites of any protein of interest, including, for example, transcription factors, RNA polymerases, or histones with or without various modifications. In addition to allowing the determination of occupancy sites within one cell type and under one condition, this method allows, in principle, the establishment and comparison of occupancy maps in various cell types, tissues, and conditions. Such comparisons require, however, that samples be normalized. Widely used normalization methods that include a quantile normalization step perform well when factor occupancy varies at a subset of sites, but may miss uniform genome-wide increases or decreases in site occupancy. We describe a spike adjustment procedure (SAP) that, unlike commonly used normalization methods intervening at the analysis stage, entails an experimental step prior to immunoprecipitation. A constant, low amount from a single batch of chromatin of a foreign genome is added to the experimental chromatin. This "spike" chromatin then serves as an internal control to which the experimental signals can be adjusted. We show that the method improves similarity between replicates and reveals biological differences including global and largely uniform changes.

Hepatoblastoma is the most common malignant pediatric liver cancer. Histological evaluation of tumor biopsies is used to distinguish among the different subtypes of hepatoblastoma, with fetal and embryonal representing the two main epithelial components. With frequent CTNNB1 mutations, hepatoblastoma is a Wnt/β-catenin-driven malignancy. Considering that Wnt activation has been associated with tumor metabolic reprogramming, we characterized the metabolic profile of cells from hepatoblastoma and compared it to cells from hepatocellular carcinoma. First, we demonstrated that glucose transporter GLUT3 is a direct TCF4/β-catenin target gene. RNA sequencing enabled to identify molecular and metabolic features specific to hepatoblastoma and revealed that several glycolytic enzymes are overexpressed in embryonal-like compared to fetal-like tumor cells. This led us to implement successfully three biomarkers to distinguish embryonal from fetal components by immunohistochemistry from a large panel of human hepatoblastoma samples. Functional analyses demonstrated that embryonal-like hepatoblastoma cells are highly glycolytic and sensitive to hexokinase-1 silencing. Altogether, our findings reveal a new, metabolic classification of human hepatoblastoma, with potential future implications for patients' diagnosis and treatment.

The criteria for choosing relevant cell lines among a vast panel of available intestinal-derived lines exhibiting a wide range of functional properties are still ill-defined. The objective of this study was, therefore, to establish objective criteria for choosing relevant cell lines to assess their appropriateness as tumor models as well as for drug absorption studies.

Down syndrome (DS) is characterized by extensive phenotypic variability, with most traits occurring in only a fraction of affected individuals. Substantial gene-expression variation is present among unaffected individuals, and this variation has a strong genetic component. Since DS is caused by genomic-dosage imbalance, we hypothesize that gene-expression variation of human chromosome 21 (HSA21) genes in individuals with DS has an impact on the phenotypic variability among affected individuals. We studied gene-expression variation in 14 lymphoblastoid and 17 fibroblast cell lines from individuals with DS and an equal number of controls. Gene expression was assayed using quantitative real-time polymerase chain reaction on 100 and 106 HSA21 genes and 23 and 26 non-HSA21 genes in lymphoblastoid and fibroblast cell lines, respectively. Surprisingly, only 39% and 62% of HSA21 genes in lymphoblastoid and fibroblast cells, respectively, showed a statistically significant difference between DS and normal samples, although the average up-regulation of HSA21 genes was close to the expected 1.5-fold in both cell types. Gene-expression variation in DS and normal samples was evaluated using the Kolmogorov-Smirnov test. According to the degree of overlap in expression levels, we classified all genes into 3 groups: (A) nonoverlapping, (B) partially overlapping, and (C) extensively overlapping expression distributions between normal and DS samples. We hypothesize that, in each cell type, group A genes are the most dosage sensitive and are most likely involved in the constant DS traits, group B genes might be involved in variable DS traits, and group C genes are not dosage sensitive and are least likely to participate in DS pathological phenotypes. This study provides the first extensive data set on HSA21 gene-expression variation in DS and underscores its role in modulating the outcome of gene-dosage imbalance.

The elucidation of mechanisms underlying telencephalic neural development has been limited by the lack of knowledge regarding the molecular and cellular aspects of the ganglionic eminence (GE), an embryonic structure that supplies the brain with diverse sets of GABAergic neurons. Here, we report a comprehensive transcriptomic analysis of this structure including its medial (MGE), lateral (LGE) and caudal (CGE) subdivisions and its temporal dynamics in 12.5 to 16 day-old rat embryos. Surprisingly, comparison across subdivisions showed that CGE gene expression was the most unique providing unbiased genetic evidence for its differentiation from MGE and LGE. The molecular signature of the CGE comprised a large set of genes, including Rwdd3, Cyp26b1, Nr2f2, Egr3, Cpta1, Slit3, and Hod, of which several encode cell signaling and migration molecules such as WNT5A, DOCK9, VSNL1 and PRG1. Temporal analysis of the MGE revealed differential expression of unique sets of cell specification and migration genes, with early expression of Hes1, Lhx2, Ctgf and Mdk, and late enrichment of Olfm3, SerpinE2 and Wdr44. These GE profiles reveal new candidate regulators of spatiotemporally governed GABAergic neuronogenesis.

Molecular characterization of breast and other cancers by gene expression profiling has corroborated existing classifications and revealed novel subtypes. Most profiling studies are based on fresh frozen (FF) tumor material which is available only for a limited number of samples while thousands of tumor samples exist as formalin-fixed, paraffin-embedded (FFPE) blocks. Unfortunately, RNA derived of FFPE material is fragmented and chemically modified impairing expression measurements by standard procedures. Robust protocols for isolation of RNA from FFPE material suitable for stable and reproducible measurement of gene expression (e.g. by quantitative reverse transcriptase PCR, QPCR) remain a major challenge.

Gene expression data obtained in large studies hold great promises for discovering disease signatures or subtypes through data analysis. It is also prone to technical variation, whose removal is essential to avoid spurious discoveries. Because this variation is not always known and can be confounded with biological signals, its removal is a challenging task. Here we provide a step-wise procedure and comprehensive analysis of the MINDACT microarray dataset. The MINDACT trial enrolled 6693 breast cancer patients and prospectively validated the gene expression signature MammaPrint for outcome prediction. The study also yielded a full-transcriptome microarray for each tumor. We show for the first time in such a large dataset how technical variation can be removed while retaining expected biological signals. Because of its unprecedented size, we hope the resulting adjusted dataset will be an invaluable tool to discover or test gene expression signatures and to advance our understanding of breast cancer.

High T-cell infiltration in colorectal cancer (CRC) correlates with a favorable disease outcome and immunotherapy response. This, however, is only observed in a small subset of CRC patients. A better understanding of the factors influencing tumor T-cell responses in CRC could inspire novel therapeutic approaches to achieve broader immunotherapy responsiveness. Here, we investigated T cell-suppressive properties of different myeloid cell types in an inducible colon tumor mouse model. The most potent inhibitors of T-cell activity were tumor-infiltrating neutrophils. Gene expression analysis and combined in vitro and in vivo tests indicated that T-cell suppression is mediated by neutrophil-secreted metalloproteinase activation of latent TGFβ. CRC patient neutrophils similarly suppressed T cells via TGFβ in vitro, and public gene expression datasets suggested that T-cell activity is lowest in CRCs with combined neutrophil infiltration and TGFβ activation. Thus, the interaction of neutrophils with a TGFβ-rich tumor microenvironment may represent a conserved immunosuppressive mechanism in CRC.

Virus-specific PD1+ Tcf1+ memory-like CD8+ T cells (TMLs) maintain the CD8+ T cell response during chronic viral infection. However, the fate of these cells following cessation of persistent antigen exposure has been unclear. Here, we find that TMLs persist upon transfer into antigen-free hosts and form memory following recall stimulation. Phenotypic, functional, and transcriptome analyses show that TML-derived memory cells resemble those arising in response to acute, resolved infection, but they retain features of chronically stimulated cells, including elevated PD-1 and Tox and reduced cytokine expression. This chronic infection imprint is largely accounted for by constitutive Tox expression. Virus-specific Tcf1+ CD8+ T cells that persist after clearance of systemic infection also display a chronic infection imprint. Notwithstanding, renewed virus exposure induces a recall response, which controls virus infection in part. Thus, cessation of chronic antigen exposure yields a memory CD8+ T cell compartment that reflects prior stimulation.

The mechanisms regulating exhaustion of tumor-infiltrating lymphocytes (TIL) and responsiveness to PD-1 blockade remain partly unknown. In human ovarian cancer, we show that tumor-specific CD8+ TIL accumulate in tumor islets, where they engage antigen and upregulate PD-1, which restrains their functions. Intraepithelial PD-1+CD8+ TIL can be, however, polyfunctional. PD-1+ TIL indeed exhibit a continuum of exhaustion states, with variable levels of CD28 costimulation, which is provided by antigen-presenting cells (APC) in intraepithelial tumor myeloid niches. CD28 costimulation is associated with improved effector fitness of exhausted CD8+ TIL and is required for their activation upon PD-1 blockade, which also requires tumor myeloid APC. Exhausted TIL lacking proper CD28 costimulation in situ fail to respond to PD-1 blockade, and their response may be rescued by local CTLA-4 blockade and tumor APC stimulation via CD40L.

The understanding of the role of B cells in patients with solid tumors remains insufficient. We found that circulating B cells produced TNFα and/or IL-6, associated with unresponsiveness and poor overall survival of melanoma patients treated with anti-CTLA4 antibody. Transcriptome analysis of B cells from melanoma metastases showed enriched expression of inflammatory response genes. Publicly available single B cell data from the tumor microenvironment revealed a negative correlation between TNFα expression and response to immune checkpoint blockade. These findings suggest that B cells contribute to tumor growth via the production of inflammatory cytokines. Possibly, these B cells are different from tertiary lymphoid structure-associated B cells, which have been described to correlate with favorable clinical outcome of cancer patients. Further studies are required to identify and characterize B cell subsets and their functions promoting or counteracting tumor growth, with the aim to identify biomarkers and novel treatment targets.

The mechanisms maintaining adult lymphatic vascular specialization throughout life and their role in coordinating inter-organ communication to sustain homeostasis remain elusive. We report that inactivation of the mechanosensitive transcription factor Foxc2 in adult lymphatic endothelium leads to a stepwise intestine-to-lung systemic failure. Foxc2 loss compromised the gut epithelial barrier, promoted dysbiosis and bacterial translocation to peripheral lymph nodes, and increased circulating levels of purine metabolites and angiopoietin-2. Commensal microbiota depletion dampened systemic pro-inflammatory cytokine levels, corrected intestinal lymphatic dysfunction, and improved survival. Foxc2 loss skewed the specialization of lymphatic endothelial subsets, leading to populations with mixed, pro-fibrotic identities and to emergence of lymph node-like endothelial cells. Our study uncovers a cross-talk between lymphatic vascular function and commensal microbiota, provides single-cell atlas of lymphatic endothelial subtypes, and reveals organ-specific and systemic effects of dysfunctional lymphatics. These effects potentially contribute to the pathogenesis of diseases, such as inflammatory bowel disease, cancer, or lymphedema.

The reason why most individuals with COVID-19 have relatively limited symptoms while other develop respiratory distress with life-threatening complications remains unknown. Increasing evidence suggests that COVID-19 associated adverse outcomes mainly rely on dysregulated immunity. Here, we compared transcriptomic profiles of blood cells from 103 patients with different severity levels of COVID-19 with that of 27 healthy and 22 influenza-infected individuals. Data provided a complete overview of SARS-CoV-2-induced immune signature, including a dramatic defect in IFN responses, a reduction of toxicity-related molecules in NK cells, an increased degranulation of neutrophils, a dysregulation of T cells, a dramatic increase in B cell function and immunoglobulin production, as well as an important over-expression of genes involved in metabolism and cell cycle in patients infected with SARS-CoV-2 compared to those infected with influenza viruses. These features also differed according to COVID-19 severity. Overall and specific gene expression patterns across groups can be visualized on an interactive website (https://bix.unil.ch/covid/). Collectively, these transcriptomic host responses to SARS-CoV-2 infection are discussed in the context of current studies, thereby improving our understanding of COVID-19 pathogenesis and shaping the severity level of COVID-19.

MicroRNAs (miRNAs) have important regulatory functions in cellular processes and have shown promising potential as prognostic markers for disease outcome in patients with cancer. The aim of the present study was to find miRNA expression profiles in whole blood that were prognostic for overall survival (OS) in patients with metastatic colorectal cancer (mCRC) treated with cetuximab and irinotecan.

The reciprocal interaction between cancer cells and the tissue-specific stroma is critical for primary and metastatic tumor growth progression. Prostate cancer cells colonize preferentially bone (osteotropism), where they alter the physiological balance between osteoblast-mediated bone formation and osteoclast-mediated bone resorption, and elicit prevalently an osteoblastic response (osteoinduction). The molecular cues provided by osteoblasts for the survival and growth of bone metastatic prostate cancer cells are largely unknown. We exploited the sufficient divergence between human and mouse RNA sequences together with redefinition of highly species-specific gene arrays by computer-aided and experimental exclusion of cross-hybridizing oligonucleotide probes. This strategy allowed the dissection of the stroma (mouse) from the cancer cell (human) transcriptome in bone metastasis xenograft models of human osteoinductive prostate cancer cells (VCaP and C4-2B). As a result, we generated the osteoblastic bone metastasis-associated stroma transcriptome (OB-BMST). Subtraction of genes shared by inflammation, wound healing and desmoplastic responses, and by the tissue type-independent stroma responses to a variety of non-osteotropic and osteotropic primary cancers generated a curated gene signature ("Core" OB-BMST) putatively representing the bone marrow/bone-specific stroma response to prostate cancer-induced, osteoblastic bone metastasis. The expression pattern of three representative Core OB-BMST genes (PTN, EPHA3 and FSCN1) seems to confirm the bone specificity of this response. A robust induction of genes involved in osteogenesis and angiogenesis dominates both the OB-BMST and Core OB-BMST. This translates in an amplification of hematopoietic and, remarkably, prostate epithelial stem cell niche components that may function as a self-reinforcing bone metastatic niche providing a growth support specific for osteoinductive prostate cancer cells. The induction of this combinatorial stem cell niche is a novel mechanism that may also explain cancer cell osteotropism and local interference with hematopoiesis (myelophthisis). Accordingly, these stem cell niche components may represent innovative therapeutic targets and/or serum biomarkers in osteoblastic bone metastasis.

Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms.