Loss of Hepatic Leucine-Rich Repeat-Containing G-Protein Coupled Receptors 4 and 5 Promotes Nonalcoholic Fatty Liver Disease.

The roof plate-specific spondin-leucine-rich repeat-containing G-protein coupled receptor 4/5 (LGR4/5)-zinc and ring finger 3 (ZNRF3)/ring finger protein 43 (RNF43) module is a master regulator of hepatic Wnt/ß-catenin signaling and metabolic zonation. However, its impact on nonalcoholic fatty liver disease (NAFLD) remains unclear. The current study investigated whether hepatic epithelial cell-specific loss of the Wnt/ß-catenin modulator Lgr4/5 promoted NAFLD. The 3- and 6-month-old mice with hepatic epithelial cell-specific deletion of both receptors Lgr4/5 (Lgr4/5dLKO) were compared with control mice fed with normal diet (ND) or high-fat diet (HFD). Six-month-old HFD-fed Lgr4/5dLKO mice developed hepatic steatosis and fibrosis but the control mice did not. Serum cholesterol-high-density lipoprotein and total cholesterol levels in 3- and 6-month-old HFD-fed Lgr4/5dLKO mice were decreased compared with those in control mice. An ex vivo primary hepatocyte culture assay and a comprehensive bile acid (BA) characterization in liver, plasma, bile, and feces demonstrated that ND-fed Lgr4/5dLKO mice had impaired BA secretion, predisposing them to develop cholestatic characteristics. Lipidome and RNA-sequencing analyses demonstrated severe alterations in several lipid species and pathways controlling lipid metabolism in the livers of Lgr4/5dLKO mice. In conclusion, loss of hepatic Wnt/ß-catenin activity by Lgr4/5 deletion led to loss of BA secretion, cholestatic features, altered lipid homeostasis, and deregulation of lipoprotein pathways. Both BA and intrinsic lipid alterations contributed to the onset of NAFLD.

Unraveling the role of low-frequency mutated genes in breast cancer.

Motivation: Breast cancer is the most commonly diagnosed malignancy in women and the second cause of cancer death in developed countries. While advancements in early detection and therapeutic options have led to a significant decrease in mortality, response to treatment is affected by the genetic heterogeneity of the disease. Recent genome-wide DNA mutation analyses revealed the existence of hundreds of low-frequency mutated genes, in addition to known cancer drivers: a finding that is prompting research into the impact of these genes on the pathogenesis of the disease. Results: Herein, we describe a strategy towards the characterization of the role of low-frequency mutated genes in breast cancer. Through the combined analyses of publicly available gene expression and mutational datasets, we identified several Cancer Gene Modules (CMs) that we re-organized in Gene Regulatory Networks (GRN) enriched in low-frequency mutated genes. Importantly, these low-frequency mutated genes were mutually exclusive with known cancer drivers. Finally, we provide evidence that gene expression analysis of these mutated GRNs can predict resistance/sensitivity to chemotherapeutic drugs for breast cancer treatment. Availability and implementation: Datasets are available at https://www.ncbi.nlm.nih.gov/geo/ and at https://www.ebi.ac.uk/ega/datasets/. Molecular signatures and GSEA software are available at http://www.gsea-msigdb.org/gsea/index.jsp. Source codes are available at https://github.com/EleonoraLusito/Reverse_Engineering_BC_GRNs. Supplementary information: Supplementary data are available at Bioinformatics online.

Screening of Intestinal Crypt Organoids: A Simple Readout for Complex Biology.

Oral and intestinal mucositis is a debilitating side effect of radiation treatment. A mouse model of radiation-induced mucositis leads to weight loss and tissue damage, reflecting the human ailment as it responds to keratinocyte growth factor (KGF), the standard-of-care treatment. Cultured intestinal crypt organoids allowed the development of an assay monitoring the effect of treatments of intestinal epithelium to radiation-induced damage. This in vitro assay resembles the mouse model as KGF and roof plate-specific spondin-1 (RSPO1) enhanced crypt organoid recovery following radiation. Screening identified compounds that increased the survival of organoids postradiation. Testing of these compounds revealed that the organoids changed their responses over time. Unbiased transcriptome analysis was performed on crypt organoid cultures at various time points in culture to investigate this adaptive behavior. A number of genes and pathways were found to be modulated over time, providing a rationale for the altered sensitivity of the organoid cultures. This report describes an in vitro assay that reflects aspects of human disease. The assay was used to identify bioactive compounds, which served as probes to interrogate the biology of crypt organoids over prolonged culture. The pathways that are changing over time may offer potential targets for treatment of mucositis.

A Vulnerability of a Subset of Colon Cancers with Potential Clinical Utility.

BRAF(V600E) mutant colon cancers (CCs) have a characteristic gene expression signature that is also found in some tumors lacking this mutation. Collectively, they are referred to as "BRAF-like" tumors and represent some 20% of CCs. We used a shRNA-based genetic screen focused on genes upregulated in BRAF(V600E) CCs to identify vulnerabilities of this tumor subtype that might be exploited therapeutically. Here, we identify RANBP2 (also known as NUP358) as essential for survival of BRAF-like, but not for non-BRAF-like, CC cells. Suppression of RANBP2 results in mitotic defects only in BRAF-like CC cells, leading to cell death. Mechanistically, RANBP2 silencing reduces microtubule outgrowth from the kinetochores, thereby inducing spindle perturbations, providing an explanation for the observed mitotic defects. We find that BRAF-like CCs display far greater sensitivity to the microtubule poison vinorelbine both in vitro and in vivo, suggesting that vinorelbine is a potential tailored treatment for BRAF-like CCs.

Mining cancer gene expression databases for latent information on intronic microRNAs.

Around 50% of all human microRNAs reside within introns of coding genes and are usually co-transcribed. Gene expression datasets, therefore, should contain a wealth of miRNA-relevant latent information, exploitable for many basic and translational research aims. The present study was undertaken to investigate this possibility. We developed an in silico approach to identify intronic-miRNAs relevant to breast cancer, using public gene expression datasets. This led to the identification of a miRNA signature for aggressive breast cancer, and to the characterization of novel roles of selected miRNAs in cancer-related biological phenotypes. Unexpectedly, in a number of cases, expression regulation of the intronic-miRNA was more relevant than the expression of their host gene. These results provide a proof of principle for the validity of our intronic miRNA mining strategy, which we envision can be applied not only to cancer research, but also to other biological and biomedical fields.

Gene expression patterns unveil a new level of molecular heterogeneity in colorectal cancer.

The recognition that colorectal cancer (CRC) is a heterogeneous disease in terms of clinical behaviour and response to therapy translates into an urgent need for robust molecular disease subclassifiers that can explain this heterogeneity beyond current parameters (MSI, KRAS, BRAF). Attempts to fill this gap are emerging. The Cancer Genome Atlas (TGCA) reported two main CRC groups, based on the incidence and spectrum of mutated genes, and another paper reported an EMT expression signature defined subgroup. We performed a prior free analysis of CRC heterogeneity on 1113 CRC gene expression profiles and confronted our findings to established molecular determinants and clinical, histopathological and survival data. Unsupervised clustering based on gene modules allowed us to distinguish at least five different gene expression CRC subtypes, which we call surface crypt-like, lower crypt-like, CIMP-H-like, mesenchymal and mixed. A gene set enrichment analysis combined with literature search of gene module members identified distinct biological motifs in different subtypes. The subtypes, which were not derived based on outcome, nonetheless showed differences in prognosis. Known gene copy number variations and mutations in key cancer-associated genes differed between subtypes, but the subtypes provided molecular information beyond that contained in these variables. Morphological features significantly differed between subtypes. The objective existence of the subtypes and their clinical and molecular characteristics were validated in an independent set of 720 CRC expression profiles. Our subtypes provide a novel perspective on the heterogeneity of CRC. The proposed subtypes should be further explored retrospectively on existing clinical trial datasets and, when sufficiently robust, be prospectively assessed for clinical relevance in terms of prognosis and treatment response predictive capacity. Original microarray data were uploaded to the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress/) under Accession Nos E-MTAB-990 and E-MTAB-1026.

DUSP 4 expression identifies a subset of colorectal cancer tumors that differ in MAPK activation, regardless of the genotype.

As dual-specificity phosphatase (DUSP) expression has been correlated to sensitivity to MEK inhibitors, DUSP expression levels may indicate activation of the mitogen-activated protein kinase (MAPK) pathway in many tumor types. In this study, we investigate if DUSP levels can indicate different levels of MAPK activation within colorectal cancer (CRC) patients. In three different CRC patient microarray datasets, we analyzed the expression of DUSP1. DUSP4 and DUSP6 according to mutational status, their correlation with survival and their association with different clinical characteristics. DUSP4 was significantly differentially expressed between all mutational subgroups with the highest expression in BRAF mutated tumors. Moreover, high DUSP4 expression was associated with a worse overall survival and with clinical characteristics typical for BRAF mutant patients. The use of DUSP expression as a predictive biomarker towards MAPK targeted therapy in CRC patients needs further investigation.

A snapshot of the physical and functional wiring of the Eps15 homology domain network in the nematode.

Protein interaction modules coordinate the connections within and the activity of intracellular signaling networks. The Eps15 Homology (EH) module, a protein-protein interaction domain that is a key feature of the EH-network, was originally identified in a few proteins involved in endocytosis and vesicle trafficking, and has subsequently also been implicated in actin reorganization, nuclear shuttling, and DNA repair. Here we report an extensive characterization of the physical connections and of the functional wirings of the EH-network in the nematode. Our data show that one of the major physiological roles of the EH-network is in neurotransmission. In addition, we found that the proteins of the network intersect, and possibly coordinate, a number of "territories" of cellular activity including endocytosis/recycling/vesicle transport, actin dynamics, general metabolism and signal transduction, ubiquitination/degradation of proteins, DNA replication/repair, and miRNA biogenesis and processing.

Integrated analysis of molecular and clinical prognostic factors in stage II/III colon cancer.

BACKGROUND: The prognostic potential of individual clinical and molecular parameters in stage II/III colon cancer has been investigated, but a thorough multivariable assessment of their relative impact is missing. METHODS: Tumors from patients (N = 1404) in the PETACC3 adjuvant chemotherapy trial were examined for BRAF and KRAS mutations, microsatellite instability (MSI), chromosome 18q loss of heterozygosity (18qLOH), and SMAD4 expression. Their importance in predicting relapse-free survival (RFS) and overall survival (OS) was assessed by Kaplan-Meier analyses, Cox regression models, and recursive partitioning trees. All statistical tests were two-sided. RESULTS: MSI-high status and SMAD4 focal loss of expression were identified as independent prognostic factors with better RFS (hazard ratio [HR] of recurrence = 0.54, 95% CI = 0.37 to 0.81, P = .003) and OS (HR of death = 0.43, 95% CI = 0.27 to 0.70, P = .001) for MSI-high status and worse RFS (HR = 1.47, 95% CI = 1.19 to 1.81, P < .001) and OS (HR = 1.58, 95% CI = 1.23 to 2.01, P < .001) for SMAD4 loss. 18qLOH did not have any prognostic value in RFS or OS. Recursive partitioning identified refinements of TNM into new clinically interesting prognostic subgroups. Notably, T3N1 tumors with MSI-high status and retained SMAD4 expression had outcomes similar to stage II disease. CONCLUSIONS: Concomitant assessment of molecular and clinical markers in multivariable analysis is essential to confirm or refute their independent prognostic value. Including molecular markers with independent prognostic value might allow more accurate prediction of prognosis than TNM staging alone.

A comprehensive characterization of genome-wide copy number aberrations in colorectal cancer reveals novel oncogenes and patterns of alterations.

To develop a comprehensive overview of copy number aberrations (CNAs) in stage-II/III colorectal cancer (CRC), we characterized 302 tumors from the PETACC-3 clinical trial. Microsatellite-stable (MSS) samples (n = 269) had 66 minimal common CNA regions, with frequent gains on 20 q (72.5%), 7 (41.8%), 8 q (33.1%) and 13 q (51.0%) and losses on 18 (58.6%), 4 q (26%) and 21 q (21.6%). MSS tumors have significantly more CNAs than microsatellite-instable (MSI) tumors: within the MSI tumors a novel deletion of the tumor suppressor WWOX at 16 q23.1 was identified (p<0.01). Focal aberrations identified by the GISTIC method confirmed amplifications of oncogenes including EGFR, ERBB2, CCND1, MET, and MYC, and deletions of tumor suppressors including TP53, APC, and SMAD4, and gene expression was highly concordant with copy number aberration for these genes. Novel amplicons included putative oncogenes such as WNK1 and HNF4A, which also showed high concordance between copy number and expression. Survival analysis associated a specific patient segment featured by chromosome 20 q gains to an improved overall survival, which might be due to higher expression of genes such as EEF1B2 and PTK6. The CNA clustering also grouped tumors characterized by a poor prognosis BRAF-mutant-like signature derived from mRNA data from this cohort. We further revealed non-random correlation between CNAs among unlinked loci, including positive correlation between 20 q gain and 8 q gain, and 20 q gain and chromosome 18 loss, consistent with co-selection of these CNAs. These results reinforce the non-random nature of somatic CNAs in stage-II/III CRC and highlight loci and genes that may play an important role in driving the development and outcome of this disease.

Abrogation of junctional adhesion molecule-A expression induces cell apoptosis and reduces breast cancer progression.

Intercellular junctions promote homotypic cell to cell adhesion and transfer intracellular signals which control cell growth and apoptosis. Junctional adhesion molecule-A (JAM-A) is a transmembrane immunoglobulin located at tight junctions of normal epithelial cells of mammary ducts and glands. In the present paper we show that JAM-A acts as a survival factor for mammary carcinoma cells. JAM-A null mice expressing Polyoma Middle T under MMTV promoter develop significantly smaller mammary tumors than JAM-A positive mice. Angiogenesis and inflammatory or immune infiltrate were not statistically modified in absence of JAM-A but tumor cell apoptosis was significantly increased. Tumor cells isolated from JAM-A null mice or 4T1 cells incubated with JAM-A blocking antibodies showed reduced growth and increased apoptosis which paralleled altered junctional architecture and adhesive function. In a breast cancer clinical data set, tissue microarray data show that JAM-A expression correlates with poor prognosis. Gene expression analysis of mouse tumor samples showed a correlation between genes enriched in human G3 tumors and genes over expressed in JAM-A +/+ mammary tumors. Conversely, genes enriched in G1 human tumors correlate with genes overexpressed in JAM-A-/- tumors. We conclude that down regulation of JAM-A reduces tumor aggressive behavior by increasing cell susceptibility to apoptosis. JAM-A may be considered a negative prognostic factor and a potential therapeutic target.

Interplay between oncogene-induced DNA damage response and heterochromatin in senescence and cancer.

Two major mechanisms have been causally implicated in the establishment of cellular senescence: the activation of the DNA damage response (DDR) pathway and the formation of senescence-associated heterochromatic foci (SAHF). Here we show that in human fibroblasts resistant to premature p16(INK4a) induction, SAHF are preferentially formed following oncogene activation but are not detected during replicative cellular senescence or on exposure to a variety of senescence-inducing stimuli. Oncogene-induced SAHF formation depends on DNA replication and ATR (ataxia telangiectasia and Rad3-related). Inactivation of ATM (ataxia telangiectasia mutated) or p53 allows the proliferation of oncogene-expressing cells that retain increased heterochromatin induction. In human cancers, levels of heterochromatin markers are higher than in normal tissues, and are independent of the proliferative index or stage of the tumours. Pharmacological and genetic perturbation of heterochromatin in oncogene-expressing cells increase DDR signalling and lead to apoptosis. In vivo, a histone deacetylase inhibitor (HDACi) causes heterochromatin relaxation, increased DDR, apoptosis and tumour regression. These results indicate that heterochromatin induced by oncogenic stress restrains DDR and suggest that the use of chromatin-modifying drugs in cancer therapies may benefit from the study of chromatin and DDR status of tumours.

Alterations of the Notch pathway in lung cancer.

Notch signaling regulates cell specification and homeostasis of stem cell compartments, and it is counteracted by the cell fate determinant Numb. Both Numb and Notch have been implicated in human tumors. Here, we show that Notch signaling is altered in approximately one third of non-small-cell lung carcinomas (NSCLCs), which are the leading cause of cancer-related deaths: in approximately 30% of NSCLCs, loss of Numb expression leads to increased Notch activity, while in a smaller fraction of cases (around 10%), gain-of-function mutations of the NOTCH-1 gene are present. Activation of Notch correlates with poor clinical outcomes in NSCLC patients without TP53 mutations. Finally, primary epithelial cell cultures, derived from NSCLC harboring constitutive activation of the Notch pathway, are selectively killed by inhibitors of Notch (gamma-secretase inhibitors), showing that the proliferative advantage of these tumors is dependent upon Notch signaling. Our results show that the deregulation of the Notch pathway is a relatively frequent event in NSCLCs and suggest that it might represent a possible target for molecular therapies in these tumors.

Gene expression time-series analysis of camptothecin effects in U87-MG and DBTRG-05 glioblastoma cell lines.

BACKGROUND: The clinical efficacy of camptothecin (CPT), a drug specifically targeting topoisomerase I (TopoI), is under evaluation for the treatment of malignant gliomas. Due to the high unresponsiveness of these tumours to chemotherapy, it would be very important to study the signalling network that drives camptothecin outcome in this type of cancer cells. To address this issue, we had previously compared the expression profile of human U87-MG glioblastoma cells with that of a CPT-resistant counterpart, giving evidence that the development of a robust inflammatory response was the main transcriptional effect associated with CPT resistance. Here we report time-related changes and cell line specific patterns of gene expression after CPT treatment by using two p53 wild-type glioblastoma cell lines, U87-MG and DBTRG-05, with different sensitivities to TopoI inhibition. RESULTS: First, we demonstrated that CPT treatment brings the two cell lines to completely different outcomes: accelerated senescence in U87-MG and apoptosis in DBTRG-05 cells. Then, to understand the different susceptibility to CPT, we used oligo-microarray to identify the genes whose expression was regulated during a time-course treatment, ranging from 2 h to 72 h. The statistical analysis of microarray data by MAANOVA (MicroArray ANalysis Of VAriance) showed much less modulated genes in apoptotic DBTRG-05 cells (155) with respect to the senescent U87-MG cells (3168), where the number of down-regulated genes largely exceeded that of the up-regulated ones (80% vs. 20%). Despite this great difference, the two data-sets showed a large overlapping (60% circa) mainly due to the expression of early stress responsive genes. The use of High-Throughput GoMINER and EASE tools, for functional analysis of significantly enriched GO terms, highlighted common cellular processes and showed that U87-MG and DBTRG-05 cells shared many GO terms, which are related to the down-regulation of cell cycle and mitosis and to the up-regulation of cell growth inhibition and DNA damage.Furthermore, the down-regulation of MYC and DP1 genes, which act as key transcription factors in cell growth control, together with the inhibition of BUB1, BUB3 and MAD2 mRNAs, which are known to be involved in the spindle checkpoint pathway, were specifically associated with the execution of senescence in U87-MG cells and addressed as critical factors that could drive the choice between different CPT-inducible effectors programs. In U87-MG cells we also found inflammation response and IL1-beta induction, as late transcriptional effects of Topo I treatment but these changes were only partially involved in the senescence development, as shown by IL1-beta gene silencing. CONCLUSION: By comparing the transcription profile of two glioblastoma cell lines treated with camptothecin, we were able to identify the common cellular pathways activated upon Topo I inhibition. Moreover, our results helped in identifying some key genes whose expression seemed to be associated with the execution of senescence or apoptosis in U87-MG and DBTRG-05 cells, respectively.