Identification of biomarkers predictive of metastasis development in early-stage colorectal cancer using network-based regularization.

Colorectal cancer (CRC) is the third most common cancer and the second most deathly worldwide. It is a very heterogeneous disease that can develop via distinct pathways where metastasis is the primary cause of death. Therefore, it is crucial to understand the molecular mechanisms underlying metastasis. RNA-sequencing is an essential tool used for studying the transcriptional landscape. However, the high-dimensionality of gene expression data makes selecting novel metastatic biomarkers problematic. To distinguish early-stage CRC patients at risk of developing metastasis from those that are not, three types of binary classification approaches were used: (1) classification methods (decision trees, linear and radial kernel support vector machines, logistic regression, and random forest) using differentially expressed genes (DEGs) as input features; (2) regularized logistic regression based on the Elastic Net penalty and the proposed iTwiner-a network-based regularizer accounting for gene correlation information; and (3) classification methods based on the genes pre-selected using regularized logistic regression. Classifiers using the DEGs as features showed similar results, with random forest showing the highest accuracy. Using regularized logistic regression on the full dataset yielded no improvement in the methods' accuracy. Further classification using the pre-selected genes found by different penalty factors, instead of the DEGs, significantly improved the accuracy of the binary classifiers. Moreover, the use of network-based correlation information (iTwiner) for gene selection produced the best classification results and the identification of more stable and robust gene sets. Some are known to be tumor suppressor genes (OPCML-IT2), to be related to resistance to cancer therapies (RAC1P3), or to be involved in several cancer processes such as genome stability (XRCC6P2), tumor growth and metastasis (MIR602) and regulation of gene transcription (NME2P2). We show that the classification of CRC patients based on pre-selected features by regularized logistic regression is a valuable alternative to using DEGs, significantly increasing the models' predictive performance. Moreover, the use of correlation-based penalization for biomarker selection stands as a promising strategy for predicting patients' groups based on RNA-seq data.

Genetic and microenvironmental intra-tumor heterogeneity impacts colorectal cancer evolution and metastatic development.

Colorectal cancer (CRC) is a highly diverse disease, where different genomic instability pathways shape genetic clonal diversity and tumor microenvironment. Although intra-tumor heterogeneity has been characterized in primary tumors, its origin and consequences in CRC outcome is not fully understood. Therefore, we assessed intra- and inter-tumor heterogeneity of a prospective cohort of 136 CRC samples. We demonstrate that CRC diversity is forged by asynchronous forms of molecular alterations, where mutational and chromosomal instability collectively boost CRC genetic and microenvironment intra-tumor heterogeneity. We were able to depict predictor signatures of cancer-related genes that can foresee heterogeneity levels across the different tumor consensus molecular subtypes (CMS) and primary tumor location. Finally, we show that high genetic and microenvironment heterogeneity are associated with lower metastatic potential, whereas late-emerging copy number variations favor metastasis development and polyclonal seeding. This study provides an exhaustive portrait of the interplay between genetic and microenvironment intra-tumor heterogeneity across CMS subtypes, depicting molecular events with predictive value of CRC progression and metastasis development.

HERVs establish a distinct molecular subtype in stage II/III colorectal cancer with poor outcome.

Colorectal cancer (CRC) is one of the most lethal malignancies. The extreme heterogeneity in survival rate is driving the need for new prognostic biomarkers. Human endogenous retroviruses (hERVs) have been suggested to influence tumor progression, oncogenesis and elicit an immune response. We examined multiple next-generation sequencing (NGS)-derived biomarkers in 114 CRC patients with paired whole-exome and whole-transcriptome sequencing (WES and WTS, respectively). First, we demonstrate that the median expression of hERVs can serve as a potential biomarker for prognosis, relapse, and resistance to chemotherapy in stage II and III CRC. We show that hERV expression and CD8+ tumor-infiltrating T-lymphocytes (TILs) synergistically stratify overall and relapse-free survival (OS and RFS): the median OS of the CD8-/hERV+ subgroup was 29.8 months compared with 37.5 months for other subgroups (HR = 4.4, log-rank P < 0.001). Combing NGS-based biomarkers (hERV/CD8 status) with clinicopathological factors provided a better prediction of patient survival compared to clinicopathological factors alone. Moreover, we explored the association between genomic and transcriptomic features of tumors with high hERV expression and establish this subtype as distinct from previously described consensus molecular subtypes of CRC. Overall, our results underscore a previously unknown role for hERVs in leading to a more aggressive subtype of CRC.

Epigenetic silencing of miR-125b is required for normal B-cell development.

Deregulation of several microRNAs (miRs) can influence critical developmental checkpoints during hematopoiesis as well as cell functions, eventually leading to the development of autoimmune disease or cancer. We found that miR-125b is expressed in bone marrow multipotent progenitors and myeloid cells but shut down in the B-cell lineage, and the gene encoding miR-125b lacked transcriptional activation markers in B cells. To understand the biological importance of the physiological silencing of miR-125b expression in B cells, we drove its expression in the B-cell lineage and found that dysregulated miR-125b expression impaired egress of immature B cells from the bone marrow to peripheral blood. Such impairment appeared to be mediated primarily by inhibited expression of the sphingosine-1-phosphate receptor 1 (S1PR1). Enforced expression of S1PR1 or clustered regularly interspaced short palindromic repeats/Cas9-mediated genome editing of the miR-125b targeting site in the S1PR1 3' untranslated region rescued the miR-125b-mediated defect in B-cell egress. In addition to impaired B-cell egress, miR-125b dysregulation initially reduced pre-B-cell output but later induced pre-B-cell lymphoma/leukemia in mice. Genetic deletion of IRF4 was found in miR-125b-induced B-cell cancer, but its role in oncogenic miR-125b-induced B-cell transformation is still unknown. Here, we further demonstrated an interaction of the effects of miR-125b and IRF4 in cancer induction by showing that miR125b-induced B-cell leukemia was greatly accelerated in IRF4 homozygous mutant mice. Thus, we conclude that physiological silencing of miR-125b is required for normal B-cell development and also acts as a mechanism of cancer suppression.

Dual mechanisms by which miR-125b represses IRF4 to induce myeloid and B-cell leukemias.

The oncomir microRNA-125b (miR-125b) is upregulated in a variety of human neoplastic blood disorders and constitutive upregulation of miR-125b in mice can promote myeloid and B-cell leukemia. We found that miR-125b promotes myeloid and B-cell neoplasm by inducing tumorigenesis in hematopoietic progenitor cells. Our study demonstrates that miR-125b induces myeloid leukemia by enhancing myeloid progenitor output from stem cells as well as inducing immortality, self-renewal, and tumorigenesis in myeloid progenitors. Through functional and genetic analyses, we demonstrated that miR-125b induces myeloid and B-cell leukemia by inhibiting interferon regulatory factor 4 (IRF4) but through distinct mechanisms; it induces myeloid leukemia through repressing IRF4 at the messenger RNA (mRNA) level without altering the genomic DNA and induces B-cell leukemia via genetic deletion of the gene encoding IRF4.

Heme-mediated SPI-C induction promotes monocyte differentiation into iron-recycling macrophages.

Splenic red pulp macrophages (RPM) degrade senescent erythrocytes and recycle heme-associated iron. The transcription factor SPI-C is selectively expressed by RPM and is required for their development, but the physiologic stimulus inducing Spic is unknown. Here, we report that Spic also regulated the development of F4/80(+)VCAM1(+) bone marrow macrophages (BMM) and that Spic expression in BMM and RPM development was induced by heme, a metabolite of erythrocyte degradation. Pathologic hemolysis induced loss of RPM and BMM due to excess heme but induced Spic in monocytes to generate new RPM and BMM. Spic expression in monocytes was constitutively inhibited by the transcriptional repressor BACH1. Heme induced proteasome-dependent BACH1 degradation and rapid Spic derepression. Furthermore, cysteine-proline dipeptide motifs in BACH1 that mediate heme-dependent degradation were necessary for Spic induction by heme. These findings are the first example of metabolite-driven differentiation of a tissue-resident macrophage subset and provide new insights into iron homeostasis.

The Yin and Yang of microRNAs: leukemia and immunity.

Yin and Yang are two complementary forces that together describe the nature of real-world elements. Yin is the dark side; Yang is the light side. We describe microRNAs having both Yin and Yang characteristics because they can contribute to normal function (Yang) but also to autoimmunity, myeloproliferation, and cancer (Yin). We have been working on a number of microRNAs that have these dual characteristics and here we focus on two, miR-125b and miR-146a. We have concentrated on these two RNAs because we have very extensive knowledge of them, much of it from our laboratory, and also because they provide a strong contrast: the effects of overexpression of miR-125b are rapid, suggesting that it acts directly, whereas the effects of miR-146a are slow to develop, suggesting that they arise from chronic alterations in cellular behavior.

Regulation of APC development, immune response, and autoimmunity by Bach1/HO-1 pathway in mice.

APCs are essential for innate and adaptive immunity as well as self-immune tolerance. Here, we show that the Cap'n'collar member Bach1 regulates the generation of APCs, specifically macrophages and dendritic cells, in mice. The impaired APC development in Bach1(-/-) mice was accompanied by defects in downstream T-cell responses and partial protection from experimental autoimmune encephalomyelitis. Genomewide analyses identified a panel of Bach1 target genes and ablation of the direct Bach1 target gene HO-1 exacerbated the impaired APC development observed in Bach1(-/-) mice. This was attributed to the impaired ability of HO-1(-/-)Bach1(-/-) double mutants to produce upstream APC progenitor cells, including common myeloid progenitor (CMP)-Flk2(+). By contrast, we observed an increase in hematopoietic stem-progenitor cells (HSPCs) in these mice, suggesting a developmental block in the progression of HSPCs to CMP-Flk2(+) and subsequently APCs.

Oncomir miR-125b regulates hematopoiesis by targeting the gene Lin28A.

MicroRNA-125b (miR-125b) is up-regulated in patients with leukemia. Overexpression of miR-125b alone in mice causes a very aggressive, transplantable myeloid leukemia. Before leukemia, these mice do not display elevation of white blood cells in the spleen or bone marrow; rather, the hematopoietic compartment shows lineage-skewing, with myeloid cell numbers dramatically increased and B-cell numbers severely diminished. miR-125b exerts this effect by up-regulating the number of common myeloid progenitors while inhibiting development of pre-B cells. We applied a miR-125b sponge loss of function system in vivo to show that miR-125b physiologically regulates hematopoietic development. Investigating the mechanism by which miR-125b regulates hematopoiesis, we found that, among a panel of candidate targets, the mRNA for Lin28A, an induced pluripotent stem cell gene, was most repressed by miR-125b in mouse hematopoietic stem and progenitor cells. Overexpressing Lin28A in the mouse hematopoietic system mimicked the phenotype observed on inhibiting miR-125b function, leading to a decrease in hematopoietic output. Relevant to the miR-125b overexpression phenotype, we also found that knockdown of Lin28A led to hematopoietic lineage-skewing, with increased myeloid and decreased B-cell numbers. Thus, the miR-125b target Lin28A is an important regulator of hematopoiesis and a primary target of miR-125b in the hematopoietic system.

MicroRNA-125b potentiates macrophage activation.

MicroRNA (miR)-125b expression is modulated in macrophages in response to stimulatory cues. In this study, we report a functional role of miR-125b in macrophages. We found that miR-125b is enriched in macrophages compared with lymphoid cells and whole immune tissues. Enforced expression of miR-125b drives macrophages to adapt an activated morphology that is accompanied by increased costimulatory factor expression and elevated responsiveness to IFN-γ, whereas anti-miR-125b treatment decreases CD80 surface expression. To determine whether these alterations in cell signaling, gene expression, and morphology have functional consequences, we examined the ability of macrophages with enhanced miR-125b expression to present Ags and found that they better stimulate T cell activation than control macrophages. Further indicating increased function, these macrophages were more effective at killing EL4 tumor cells in vitro and in vivo. Moreover, miR-125b repressed IFN regulatory factor 4 (IRF4), and IRF4 knockdown in macrophages mimicked the miR-125b overexpression phenotype. In summary, our evidence suggests that miR-125b is at least partly responsible for generating the activated nature of macrophages, at least partially by reducing IRF4 levels, and potentiates the functional role of macrophages in inducing immune responses.

A coactivator role of CARM1 in the dysregulation of ß-catenin activity in colorectal cancer cell growth and gene expression.

Aberrant activation of Wnt/ß-catenin signaling, resulting in the expression of Wnt-regulated oncogenes, is recognized as a critical factor in the etiology of colorectal cancer. Occupancy of ß-catenin at promoters of Wnt target genes drives transcription, but the mechanism of ß-catenin action remains poorly understood. Here, we show that CARM1 (coactivator-associated arginine methyltransferase 1) interacts with ß-catenin and positively modulates ß-catenin-mediated gene expression. In colorectal cancer cells with constitutively high Wnt/ß-catenin activity, depletion of CARM1 inhibits expression of endogenous Wnt/ß-catenin target genes and suppresses clonal survival and anchorage-independent growth. We also identified a colorectal cancer cell line (RKO) with a low basal level of ß-catenin, which is dramatically elevated by treatment with Wnt3a. Wnt3a also increased the expression of a subset of endogenous Wnt target genes, and CARM1 was required for the Wnt-induced expression of these target genes and the accompanying dimethylation of arginine 17 of histone H3. Depletion of ß-catenin from RKO cells diminished the Wnt-induced occupancy of CARM1 on a Wnt target gene, indicating that CARM1 is recruited to Wnt target genes through its interaction with ß-catenin and contributes to transcriptional activation by mediating events (including histone H3 methylation) that are downstream from the actions of ß-catenin. Therefore, CARM1 is an important positive modulator of Wnt/ß-catenin transcription and neoplastic transformation, and may thereby represent a novel target for therapeutic intervention in cancers involving aberrantly activated Wnt/ß-catenin signaling.

Glucocorticoid regulation of the circadian clock modulates glucose homeostasis.

Circadian clock genes are regulated by glucocorticoids; however, whether this regulation is a direct or secondary effect and the physiological consequences of this regulation were unknown. Here, we identified glucocorticoid response elements (GREs) at multiple clock genes and showed that 3 were directly regulated by the glucocorticoid receptor. We determined that a GRE within the core clock gene Per2 was continuously occupied during rhythmic expression and essential for glucocorticoid regulation of that gene in vivo. We further demonstrated that mice with a genomic deletion spanning this GRE expressed elevated leptin levels and were protected from glucose intolerance and insulin resistance on glucocorticoid treatment but not from muscle wasting. We conclude that Per2 is an integral component of a particular glucocorticoid regulatory pathway and that glucocorticoid regulation of the peripheral clock is selectively required for some actions of glucocorticoids.

Angiopoietin-like 4 (ANGPTL4, fasting-induced adipose factor) is a direct glucocorticoid receptor target and participates in glucocorticoid-regulated triglyceride metabolism.

Glucocorticoids are important regulators of lipid homeostasis, and chronically elevated glucocorticoid levels induce hypertriglyceridemia, hepatic steatosis, and visceral obesity. The occupied glucocorticoid receptor (GR) is a transcription factor. However, those genes regulating lipid metabolism under GR control are not fully known. Angiopoietin-like 4 (ANGPTL4, fasting-induced adipose factor), a protein inhibitor of lipoprotein lipase, is synthesized and secreted during fasting, when circulating glucocorticoid levels are physiologically increased. We therefore tested whether the ANGPTL4 gene (Angptl4) is transcriptionally controlled by GR. We show that treatment with the synthetic glucocorticoid dexamethasone increased Angptl4 mRNA levels in primary hepatocytes and adipocytes (2-3-fold) and in the livers and white adipose tissue of mice (approximately 4-fold). We tested the mechanism of this increase in H4IIE hepatoma cells and found that dexamethasone treatment increased the transcriptional rate of Angptl4. Using bioinformatics and chromatin immunoprecipitation, we identified a GR binding site within the rat Angptl4 sequence. A reporter plasmid containing this site was markedly activated by dexamethasone, indicative of a functional glucocorticoid response element. Dexamethasone treatment also increased histone H4 acetylation and DNase I accessibility in genomic regions near this site, further supporting that it is a glucocorticoid response element. Glucocorticoids promote the flux of triglycerides from white adipose tissue to liver. We found that mice lacking ANGPTL4 (Angptl4(-/-)) had reductions in dexamethasone-induced hypertriglyceridemia and hepatic steatosis, suggesting that ANGPTL4 is required for this flux. Overall, we establish that ANGPTL4 is a direct GR target that participates in glucocorticoid-regulated triglyceride metabolism.

DNA binding site sequence directs glucocorticoid receptor structure and activity.

Genes are not simply turned on or off, but instead their expression is fine-tuned to meet the needs of a cell. How genes are modulated so precisely is not well understood. The glucocorticoid receptor (GR) regulates target genes by associating with specific DNA binding sites, the sequences of which differ between genes. Traditionally, these binding sites have been viewed only as docking sites. Using structural, biochemical, and cell-based assays, we show that GR binding sequences, differing by as little as a single base pair, differentially affect GR conformation and regulatory activity. We therefore propose that DNA is a sequence-specific allosteric ligand of GR that tailors the activity of the receptor toward specific target genes.

The unfolded protein response during prostate cancer development.

Accumulation of misfolded proteins in the endoplasmic reticulum (ER) induces the unfolded protein response (UPR). The UPR promotes cell survival by adjusting ER protein folding capacity but if homeostasis cannot be re-established, apoptosis is induced. The execution of life/death decisions is regulated by the three UPR branches (IRE1, PERK, ATF6) and their downstream effectors. Events that offset the balance of the UPR branches can have devastating consequences, and UPR misregulation has been correlated with various diseases, including metabolic and neurodegenerative diseases and cancer. In cancer, upregulation of the UPR is thought to provide a growth advantage to tumor cells. In contrast to this prevailing view, we report here an analysis of data obtained by others indicating that all three UPR branches appear selectively down-regulated in mouse models of prostate tumorigenesis.

Conservation analysis predicts in vivo occupancy of glucocorticoid receptor-binding sequences at glucocorticoid-induced genes.

The glucocorticoid receptor (GR) interacts with specific GR-binding sequences (GBSs) at glucocorticoid response elements (GREs) to orchestrate transcriptional networks. Although the sequences of the GBSs are highly variable among different GREs, the precise sequence within an individual GRE is highly conserved. In this study, we examined whether sequence conservation of sites resembling GBSs is sufficient to predict GR occupancy of GREs at genes responsive to glucocorticoids. Indeed, we found that the level of conservation of these sites at genes up-regulated by glucocorticoids in mouse C3H10T1/2 mesenchymal stem-like cells correlated directly with the extent of occupancy by GR. In striking contrast, we failed to observe GR occupancy of GBSs at genes repressed by glucocorticoids, despite the occurrence of these sites at a frequency similar to that of the induced genes. Thus, GR occupancy of the GBS motif correlates with induction but not repression, and GBS conservation alone is sufficient to predict GR occupancy and GRE function at induced genes.

Cell- and gene-specific regulation of primary target genes by the androgen receptor.

The androgen receptor (AR) mediates the physiologic and pathophysiologic effects of androgens including sexual differentiation, prostate development, and cancer progression by binding to genomic androgen response elements (AREs), which influence transcription of AR target genes. The composition and context of AREs differ between genes, thus enabling AR to confer multiple regulatory functions within a single nucleus. We used expression profiling of an immortalized human prostate epithelial cell line to identify 205 androgen-responsive genes (ARGs), most of them novel. In addition, we performed chromatin immunoprecipitation to identify 524 AR binding regions and validated in reporter assays the ARE activities of several such regions. Interestingly, 67% of our AREs resided within approximately 50 kb of the transcription start sites of 84% of our ARGs. Indeed, most ARGs were associated with two or more AREs, and ARGs were sometimes themselves linked in gene clusters containing up to 13 AREs and 12 ARGs. AREs appeared typically to be composite elements, containing AR binding sequences adjacent to binding motifs for other transcriptional regulators. Functionally, ARGs were commonly involved in prostate cell proliferation, communication, differentiation, and possibly cancer progression. Our results provide new insights into cell- and gene-specific mechanisms of transcriptional regulation of androgen-responsive gene networks.

Determinants of cell- and gene-specific transcriptional regulation by the glucocorticoid receptor.

The glucocorticoid receptor (GR) associates with glucocorticoid response elements (GREs) and regulates selective gene transcription in a cell-specific manner. Native GREs are typically thought to be composite elements that recruit GR as well as other regulatory factors into functional complexes. We assessed whether GR occupancy is commonly a limiting determinant of GRE function as well as the extent to which core GR binding sequences and GRE architecture are conserved at functional loci. We surveyed 100-kb regions surrounding each of 548 known or potentially glucocorticoid-responsive genes in A549 human lung cells for GR-occupied GREs. We found that GR was bound in A549 cells predominately near genes responsive to glucocorticoids in those cells and not at genes regulated by GR in other cells. The GREs were positionally conserved at each responsive gene but across the set of responsive genes were distributed equally upstream and downstream of the transcription start sites, with 63% of them >10 kb from those sites. Strikingly, although the core GR binding sequences across the set of GREs varied extensively around a consensus, the precise sequence at an individual GRE was conserved across four mammalian species. Similarly, sequences flanking the core GR binding sites also varied among GREs but were conserved at individual GREs. We conclude that GR occupancy is a primary determinant of glucocorticoid responsiveness in A549 cells and that core GR binding sequences as well as GRE architecture likely harbor gene-specific regulatory information.