Gene Expression Variability in Mammalian Embryonic Stem Cells Using Single Cell RNA-seq Data

Overview

Journal Comput Biol Chem

Publisher Elsevier

Date 2016 Mar 9

PMID 26951854

Citations 17

Authors

Anna Mantsoki

Guillaume Devailly

Anagha Joshi

Affiliations

Soon will be listed here.

Abstract

Background: Gene expression heterogeneity contributes to development as well as disease progression. Due to technological limitations, most studies to date have focused on differences in mean expression across experimental conditions, rather than differences in gene expression variance. The advent of single cell RNA sequencing has now made it feasible to study gene expression heterogeneity and to characterise genes based on their coefficient of variation.

Methods: We collected single cell gene expression profiles for 32 human and 39 mouse embryonic stem cells and studied correlation between diverse characteristics such as network connectivity and coefficient of variation (CV) across single cells. We further systematically characterised properties unique to High CV genes.

Results: Highly expressed genes tended to have a low CV and were enriched for cell cycle genes. In contrast, High CV genes were co-expressed with other High CV genes, were enriched for bivalent (H3K4me3 and H3K27me3) marked promoters and showed enrichment for response to DNA damage and DNA repair.

Conclusions: Taken together, this analysis demonstrates the divergent characteristics of genes based on their CV. High CV genes tend to form co-expression clusters and they explain bivalency at least in part.

Citing Articles

Self-organization underlies developmental robustness in plants.

Kong S, Zhu M, Roeder A Cells Dev. 2024; :203936.

PMID: 38960068 PMC: 11688513. DOI: 10.1016/j.cdev.2024.203936.

Somatosensory neurons express specific sets of lincRNAs, and lincRNA CLAP promotes itch sensation in mice.

Wang B, Jiang B, Li G, Dong F, Luo Z, Cai B EMBO Rep. 2022; 24(2):e54313.

PMID: 36524339 PMC: 9900349. DOI: 10.15252/embr.202154313.

Stochastic expression of invasion genes in Plasmodium falciparum schizonts.

Tripathi J, Zhu L, Nayak S, Stoklasa M, Bozdech Z Nat Commun. 2022; 13(1):3004.

PMID: 35637187 PMC: 9151791. DOI: 10.1038/s41467-022-30605-z.

Identification of a ZC3H12D-regulated competing endogenous RNA network for prognosis of lung adenocarcinoma at single-cell level.

Chen W, Guo Z, Wu J, Lin G, Chen S, Lin Q BMC Cancer. 2022; 22(1):115.

PMID: 35090416 PMC: 8796579. DOI: 10.1186/s12885-021-08992-1.

DUBStepR is a scalable correlation-based feature selection method for accurately clustering single-cell data.

Ranjan B, Sun W, Park J, Mishra K, Schmidt F, Xie R Nat Commun. 2021; 12(1):5849.

PMID: 34615861 PMC: 8494900. DOI: 10.1038/s41467-021-26085-2.

References

Singh A, Chappell J, Trost R, Lin L, Wang T, Tang J . Cell-cycle control of developmentally regulated transcription factors accounts for heterogeneity in human pluripotent cells. Stem Cell Reports. 2013; 1(6):532-44. PMC: 3871385. DOI: 10.1016/j.stemcr.2013.10.009. View

Wu Z, Yang M, Liu H, Guo H, Wang Y, Cheng H . Role of nuclear receptor coactivator 3 (Ncoa3) in pluripotency maintenance. J Biol Chem. 2012; 287(45):38295-304. PMC: 3488098. DOI: 10.1074/jbc.M112.373092. View

Karolchik D, Barber G, Casper J, Clawson H, Cline M, Diekhans M . The UCSC Genome Browser database: 2014 update. Nucleic Acids Res. 2013; 42(Database issue):D764-70. PMC: 3964947. DOI: 10.1093/nar/gkt1168. View

Tang F, Barbacioru C, Bao S, Lee C, Nordman E, Wang X . Tracing the derivation of embryonic stem cells from the inner cell mass by single-cell RNA-Seq analysis. Cell Stem Cell. 2010; 6(5):468-78. PMC: 2954317. DOI: 10.1016/j.stem.2010.03.015. View

Streets A, Zhang X, Cao C, Pang Y, Wu X, Xiong L . Microfluidic single-cell whole-transcriptome sequencing. Proc Natl Acad Sci U S A. 2014; 111(19):7048-53. PMC: 4024894. DOI: 10.1073/pnas.1402030111. View

Xu H, Baroukh C, Dannenfelser R, Chen E, Tan C, Kou Y . ESCAPE: database for integrating high-content published data collected from human and mouse embryonic stem cells. Database (Oxford). 2013; 2013:bat045. PMC: 3689438. DOI: 10.1093/database/bat045. View

Mar J, Matigian N, Mackay-Sim A, Mellick G, Sue C, Silburn P . Variance of gene expression identifies altered network constraints in neurological disease. PLoS Genet. 2011; 7(8):e1002207. PMC: 3154954. DOI: 10.1371/journal.pgen.1002207. View

Pearson M, Carbone R, Sebastiani C, Cioce M, Fagioli M, Saito S . PML regulates p53 acetylation and premature senescence induced by oncogenic Ras. Nature. 2000; 406(6792):207-10. DOI: 10.1038/35018127. View

Mason E, Mar J, Laslett A, Pera M, Quackenbush J, Wolvetang E . Gene expression variability as a unifying element of the pluripotency network. Stem Cell Reports. 2014; 3(2):365-77. PMC: 4175554. DOI: 10.1016/j.stemcr.2014.06.008. View

10.

Trapnell C, Pachter L, Salzberg S . TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009; 25(9):1105-11. PMC: 2672628. DOI: 10.1093/bioinformatics/btp120. View

11.

Raj A, Rifkin S, Andersen E, van Oudenaarden A . Variability in gene expression underlies incomplete penetrance. Nature. 2010; 463(7283):913-8. PMC: 2836165. DOI: 10.1038/nature08781. View

12.

Xue Z, Huang K, Cai C, Cai L, Jiang C, Feng Y . Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing. Nature. 2013; 500(7464):593-7. PMC: 4950944. DOI: 10.1038/nature12364. View

13.

Lercher M, Urrutia A, Hurst L . Clustering of housekeeping genes provides a unified model of gene order in the human genome. Nat Genet. 2002; 31(2):180-3. DOI: 10.1038/ng887. View

14.

Vousden K, Prives C . Blinded by the Light: The Growing Complexity of p53. Cell. 2009; 137(3):413-31. DOI: 10.1016/j.cell.2009.04.037. View

15.

Bernstein B, Mikkelsen T, Xie X, Kamal M, Huebert D, Cuff J . A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006; 125(2):315-26. DOI: 10.1016/j.cell.2006.02.041. View

16.

Trapnell C, Williams B, Pertea G, Mortazavi A, Kwan G, van Baren M . Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010; 28(5):511-5. PMC: 3146043. DOI: 10.1038/nbt.1621. View

17.

Karwacki-Neisius V, Goke J, Osorno R, Halbritter F, Ng J, Weisse A . Reduced Oct4 expression directs a robust pluripotent state with distinct signaling activity and increased enhancer occupancy by Oct4 and Nanog. Cell Stem Cell. 2013; 12(5):531-45. PMC: 3650585. DOI: 10.1016/j.stem.2013.04.023. View

18.

Guberman J, Ai J, Arnaiz O, Baran J, Blake A, Baldock R . BioMart Central Portal: an open database network for the biological community. Database (Oxford). 2011; 2011:bar041. PMC: 3263598. DOI: 10.1093/database/bar041. View

19.

Schulze A, Downward J . Navigating gene expression using microarrays--a technology review. Nat Cell Biol. 2001; 3(8):E190-5. DOI: 10.1038/35087138. View

20.

Dixon J, Selvaraj S, Yue F, Kim A, Li Y, Shen Y . Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012; 485(7398):376-80. PMC: 3356448. DOI: 10.1038/nature11082. View