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.

; ; ; ; ; ; ; ; ; ; ; ;
Publication Date:
Journal Name:
Communications Biology
Nature Publishing Group
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 2416more »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.

    « 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.

    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.

  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 tRNAmore »biogenesis.

    « less
  5. Abstract Background

    RNA sequencing is a powerful approach to quantify the genome-wide distribution of mRNA molecules in a population to gain deeper understanding of cellular functions and phenotypes. However, unlike eukaryotic cells, mRNA sequencing of bacterial samples is more challenging due to the absence of a poly-A tail that typically enables efficient capture and enrichment of mRNA from the abundant rRNA molecules in a cell. Moreover, bacterial cells frequently contain 100-fold lower quantities of RNA compared to mammalian cells, which further complicates mRNA sequencing from non-cultivable and non-model bacterial species. To overcome these limitations, we report EMBR-seq (Enrichment of mRNA by Blocked rRNA), a method that efficiently depletes 5S, 16S and 23S rRNA using blocking primers to prevent their amplification.


    EMBR-seq results in 90% of the sequenced RNA molecules from anE. coliculture deriving from mRNA. We demonstrate that this increased efficiency provides a deeper view of the transcriptome without introducing technical amplification-induced biases. Moreover, compared to recent methods that employ a large array of oligonucleotides to deplete rRNA, EMBR-seq uses a single or a few oligonucleotides per rRNA, thereby making this new technology significantly more cost-effective, especially when applied to varied bacterial species. Finally, compared to existing commercial kits for bacterialmore »rRNA depletion, we show that EMBR-seq can be used to successfully quantify the transcriptome from more than 500-fold lower starting total RNA.


    EMBR-seq provides an efficient and cost-effective approach to quantify global gene expression profiles from low input bacterial samples.

    « less