skip to main content

Title: Transcriptomic and proteomic signatures of stemness and differentiation in the colon crypt

Intestinal stem cells are non-quiescent, dividing epithelial cells that rapidly differentiate into progenitor cells of the absorptive and secretory cell lineages. The kinetics of this process is rapid such that the epithelium is replaced weekly. To determine how the transcriptome and proteome keep pace with rapid differentiation, we developed a new cell sorting method to purify mouse colon epithelial cells. Here we show that alternative mRNA splicing and polyadenylation dominate changes in the transcriptome as stem cells differentiate into progenitors. In contrast, as progenitors differentiate into mature cell types, changes in mRNA levels dominate the transcriptome. RNA processing targets regulators of cell cycle, RNA, cell adhesion, SUMOylation, and Wnt and Notch signaling. Additionally, global proteome profiling detected >2,800 proteins and revealed RNA:protein patterns of abundance and correlation. Paired together, these data highlight new potentials for autocrine and feedback regulation and provide new insights into cell state transitions in the crypt.

more » « less
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Communications Biology
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Background

    Alternative RNA splicing is widely dysregulated in cancers including lung adenocarcinoma, where aberrant splicing events are frequently caused by somatic splice site mutations or somatic mutations of splicing factor genes. However, the majority of mis-splicing in cancers is unexplained by these known mechanisms. We hypothesize that the aberrant Ras signaling characteristic of lung cancers plays a role in promoting the alternative splicing observed in tumors.


    We recently performed transcriptome and proteome profiling of human lung epithelial cells ectopically expressing oncogenic KRAS and another cancer-associated Ras GTPase, RIT1. Unbiased analysis of phosphoproteome data identified altered splicing factor phosphorylation in KRAS-mutant cells, so we performed differential alternative splicing analysis using rMATS to identify significantly altered isoforms in lung epithelial cells. To determine whether these isoforms were uniquely regulated by KRAS, we performed a large-scale splicing screen in which we generated over 300 unique RNA sequencing profiles of isogenic A549 lung adenocarcinoma cells ectopically expressing 75 different wild-type or variant alleles across 28 genes implicated in lung cancer.


    Mass spectrometry data showed widespread downregulation of splicing factor phosphorylation in lung epithelial cells expressing mutant KRAS compared to cells expressing wild-type KRAS. We observed alternative splicing in the same cells, with 2196 and 2416 skipped exon events in KRASG12Vand KRASQ61Hcells, respectively, 997 of which were shared (p < 0.001 by hypergeometric test). In the high-throughput splicing screen, mutant KRAS induced the greatest number of differential alternative splicing events, second only to the RNA binding protein RBM45 and its variant RBM45M126I. We identified ten high confidence cassette exon events across multiple KRAS variants and cell lines. These included differential splicing of the Myc Associated Zinc Finger (MAZ). As MAZ regulates expression of KRAS, this splice variant may be a mechanism for the cell to modulate wild-type KRAS levels in the presence of oncogenic KRAS.


    Proteomic and transcriptomic profiling of lung epithelial cells uncovered splicing factor phosphorylation and mRNA splicing events regulated by oncogenic KRAS. These data suggest that in addition to widespread transcriptional changes, the Ras signaling pathway can promote post-transcriptional splicing changes that may contribute to oncogenic processes.

    more » « less
  2. In NSCLC, there is a pressing need for immunotherapy predictive biomarkers. The processes underlying B-cell dysfunction, as well as their prognostic importance in NSCLC, are unknown. Tumor-specific B-cell gene co-expression networks were constructed by comparing the Boolean implication modeling of single-cell RNA sequencing of NSCLC tumor B cells and normal B cells. Proliferation genes were selected from the networks using in vitro CRISPR-Cas9/RNA interfering (RNAi) screening data in more than 92 human NSCLC epithelial cell lines. The prognostic and predictive evaluation was performed using public NSCLC transcriptome and proteome profiles. A B cell proliferation and prognostic gene co-expression network was present only in normal lung B cells and missing in NSCLC tumor B cells. A nine-gene signature was identified from this B cell network that provided accurate prognostic stratification using bulk NSCLC tumor transcriptome (n = 1313) and proteome profiles (n = 103). Multiple genes (HLA-DRA, HLA-DRB1, OAS1, and CD74) differentially expressed in NSCLC B cells, peripheral blood lymphocytes, and tumor T cells had concordant prognostic indications at the mRNA and protein expression levels. The selected genes were associated with drug sensitivity/resistance to 10 commonly used NSCLC therapeutic regimens. Lestaurtinib was discovered as a potential repositioning drug for treating NSCLC. 
    more » « less

    Cytosolic Ca2+ is a highly dynamic, tightly regulated and broadly conserved cellular signal. Ca2+ dynamics have been studied widely in cellular monocultures, yet organs in vivo comprise heterogeneous populations of stem and differentiated cells. Here, we examine Ca2+ dynamics in the adult Drosophila intestine, a self-renewing epithelial organ in which stem cells continuously produce daughters that differentiate into either enteroendocrine cells or enterocytes. Live imaging of whole organs ex vivo reveals that stem-cell daughters adopt strikingly distinct patterns of Ca2+ oscillations after differentiation: enteroendocrine cells exhibit single-cell Ca2+ oscillations, whereas enterocytes exhibit rhythmic, long-range Ca2+ waves. These multicellular waves do not propagate through immature progenitors (stem cells and enteroblasts), of which the oscillation frequency is approximately half that of enteroendocrine cells. Organ-scale inhibition of gap junctions eliminates Ca2+ oscillations in all cell types – even, intriguingly, in progenitor and enteroendocrine cells that are surrounded only by enterocytes. Our findings establish that cells adopt fate-specific modes of Ca2+ dynamics as they terminally differentiate and reveal that the oscillatory dynamics of different cell types in a single, coherent epithelium are paced independently.

    more » « less
  4. Abstract

    Inorganic phosphate is an essential nutrient acquired by cells from their environment. Here, we characterize the adaptative responses of fission yeast to chronic phosphate starvation, during which cells enter a state of quiescence, initially fully reversible upon replenishing phosphate after 2 days but resulting in gradual loss of viability during 4 weeks of starvation. Time-resolved analyses of changes in mRNA levels revealed a coherent transcriptional program in which phosphate dynamics and autophagy were upregulated, while the machineries for rRNA synthesis and ribosome assembly, and for tRNA synthesis and maturation, were downregulated in tandem with global repression of genes encoding ribosomal proteins and translation factors. Consistent with the transcriptome changes, proteome analysis highlighted global depletion of 102 ribosomal proteins. Concomitant with this ribosomal protein deficit, 28S and 18S rRNAs became vulnerable to site-specific cleavages that generated temporally stable rRNA fragments. The finding that Maf1, a repressor of RNA polymerase III transcription, was upregulated during phosphate starvation prompted a hypothesis that its activity might prolong lifespan of the quiescent cells by limiting production of tRNAs. Indeed, we found that deletion of maf1 results in precocious death of phosphate-starved cells via a distinctive starvation-induced pathway associated with tRNA overproduction and dysfunctional tRNA biogenesis.

    more » « less
  5. Abstract

    The anterior pituitary controls key biological processes, including growth, metabolism, reproduction, and stress responses through distinct cell types that each secrete specific hormones. The anterior pituitary cells show a remarkable level of cell type plasticity that mediates the shifts in hormone-producing cell populations that are required to meet organismal needs. The molecular mechanisms underlying pituitary cell plasticity are not well understood. Recent work has implicated the pituitary stem cell populations and specifically, the mRNA binding proteins of the Musashi family in control of pituitary cell type identity. In this study we have identified the target mRNAs that mediate Musashi function in the adult mouse pituitary and demonstrate the requirement for Musashi function in vivo. Using Musashi RNA immunoprecipitation, we identify a cohort of 1184 mRNAs that show specific Musashi binding. Identified Musashi targets include the Gnrhr mRNA, which encodes the gonadotropin-releasing hormone receptor (GnRHR), and the Fshb mRNA, encoding follicle-stimulating hormone (FSH). Reporter assays reveal that Musashi functions to exert repression of translation of the Fshb mRNA, in addition to the previously observed repression of the Gnrhr mRNA. Importantly, mice engineered to lack Musashi in gonadotropes demonstrate a failure to repress translation of the endogenous Gnrhr and Fshb mRNAs during the estrous cycle and display a significant heterogeneity in litter sizes. The range of identified target mRNAs suggests that, in addition to these key gonadotrope proteins, Musashi may exert broad regulatory control over the pituitary proteome in a cell type–specific manner.

    more » « less