Cicero Predicts cis-Regulatory DNA Interactions from Single-Cell Chromatin Accessibility Data

Molecular Cell

Volumn 71, Issue 5, 2018, Pages 858-871.e8

  Pliner, Hannah A ^a   Packer, Jonathan S ^a   McFaline Figueroa, José L ^a   Cusanovich, Darren A ^a   Daza, Riza M ^a   Aghamirzaie, Delasa ^a   Srivatsan, Sanjay ^a   Qiu, Xiaojie ^a   Jackson, Dana ^a   Minkina, Anna ^a   Adey, Andrew C ^b   Steemers, Frank J ^c   Shendure, Jay ^a,d,e   Trapnell, Cole ^a,e  

^a UNIVERSITY OF WASHINGTON   (United States)

^b OREGON HEALTH AND SCIENCE UNIVERSITY   (United States)

^c ILLUMINA INC   (United States)

^d HOWARD HUGHES MEDICAL INSTITUTE   (United States)

^e BROTMAN BATY INSTITUTE FOR PRECISION MEDICINE   (United States)

Author keywords

ATAC seq; chromatin accessibility; co accessibility; gene regulation; machine learning; myoblast differentiation; single cell

Indexed keywords

E1A ASSOCIATED P300 PROTEIN; HISTONE; HOMEOBOX PROTEIN MEIS1; MYOD1 PROTEIN; PRE B CELL LEUKEMIA TRANSCRIPTION FACTOR 1; TRANSCRIPTION FACTOR; DNA;

ADOLESCENT; ALGORITHM; ARTICLE; CELL DIFFERENTIATION; CHROMATIN; CONTROLLED STUDY; DNA ELEMENT; DNA STRUCTURE; FEMALE; GENE CONTROL; GENE EXPRESSION; GENE TARGETING; HISTONE MODIFICATION; HUMAN; HUMAN CELL; INTERACTIONS WITH DNA; NORMAL HUMAN; PREDICTION; PROMOTER REGION; PROTEIN DNA INTERACTION; REGULATORY DNA SEQUENCE; SINGLE CELL ANALYSIS; SKELETAL MYOBLAST; TEMPORAL ANALYSIS; CHROMATIN ASSEMBLY AND DISASSEMBLY; ENHANCER REGION; GENE EXPRESSION REGULATION; GENETICS; HOMEOBOX; MYOBLAST; PHYSIOLOGY;

ADOLESCENT; CELL DIFFERENTIATION; CHROMATIN; CHROMATIN ASSEMBLY AND DISASSEMBLY; DNA; ENHANCER ELEMENTS, GENETIC; FEMALE; GENE EXPRESSION REGULATION; GENES, HOMEOBOX; HISTONES; HUMANS; MYOBLASTS; TRANSCRIPTION FACTORS;

EID: 85053831649     PISSN: 10972765     EISSN: 10974164     Source Type: Journal    
DOI: 10.1016/j.molcel.2018.06.044     Document Type: Article

References (62)

1. Amini, S., Pushkarev, D., Christiansen, L., Kostem, E., Royce, T., Turk, C., Pignatelli, N., Adey, A., Kitzman, J.O., Vijayan, K., et al. Haplotype-resolved whole-genome sequencing by contiguity-preserving transposition and combinatorial indexing. Nat Genet 46 (2014), 1343–1349.
2. Benezra, R., Davis, R.L., Lockshon, D., Turner, D.L., Weintraub, H., The protein Id: a negative regulator of helix-loop-helix DNA binding proteins. Cell 61 (1990), 49–59.
3. Berkes, C.A., Bergstrom, D.A., Penn, B.H., Seaver, K.J., Knoepfler, P.S., Tapscott, S.J., Pbx marks genes for activation by MyoD indicating a role for a homeodomain protein in establishing myogenic potential. Mol. Cell 14 (2004), 465–477.
4. Beygelzimer, A., Kakadet, S., Langford, J., Arya, S., Mount, D., and Li, S. (2013). FNN: fast nearest neighbor search algorithms and applications. R Package Version 1.
5. Blondel, V.D., Guillaume, J.-L., Lambiotte, R., Lefebvre, E., Fast unfolding of communities in large networks. J. Stat. Mech., 2008, 10.1088/1742-5468/2008/10/P10008.
6. Bolger, A.M., Lohse, M., Usadel, B., Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30 (2014), 2114–2120.
7. Breiman, L., Bagging predictors. Mach. Learn. 24 (1996), 123–140.
8. Budden, D.M., Hurley, D.G., Crampin, E.J., Predictive modelling of gene expression from transcriptional regulatory elements. Brief. Bioinform. 16 (2015), 616–628.
9. Buenrostro, J.D., Giresi, P.G., Zaba, L.C., Chang, H.Y., Greenleaf, W.J., Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat. Methods 10 (2013), 1213–1218.
10. Cairns, J., Freire-Pritchett, P., Wingett, S.W., Várnai, C., Dimond, A., Plagnol, V., Zerbino, D., Schoenfelder, S., Javierre, B.-M., Osborne, C., et al. CHiCAGO: robust detection of DNA looping interactions in Capture Hi-C data. Genome Biol., 17, 2016, 127.
11. Cao, Y., Yao, Z., Sarkar, D., Lawrence, M., Sanchez, G.J., Parker, M.H., MacQuarrie, K.L., Davison, J., Morgan, M.T., Ruzzo, W.L., et al. Genome-wide MyoD binding in skeletal muscle cells: a potential for broad cellular reprogramming. Dev. Cell 18 (2010), 662–674.
12. Cusanovich, D.A., Daza, R., Adey, A., Pliner, H.A., Christiansen, L., Gunderson, K.L., Steemers, F.J., Trapnell, C., Shendure, J., Multiplex single cell profiling of chromatin accessibility by combinatorial cellular indexing. Science 348 (2015), 910–914.
13. Cusanovich, D.A., Reddington, J.P., Garfield, D.A., Daza, R.M., Aghamirzaie, D., Marco-Ferreres, R., Pliner, H.A., Christiansen, L., Qiu, X., Steemers, F.J., et al. The cis-regulatory dynamics of embryonic development at single-cell resolution. Nature 555 (2018), 538–542.
14. de la Serna, I.L., Ohkawa, Y., Berkes, C.A., Bergstrom, D.A., Dacwag, C.S., Tapscott, S.J., Imbalzano, A.N., MyoD targets chromatin remodeling complexes to the myogenin locus prior to forming a stable DNA-bound complex. Mol. Cell. Biol. 25 (2005), 3997–4009.
15. de Laat, W., Grosveld, F., Spatial organization of gene expression: the active chromatin hub. Chromosome Res. 11 (2003), 447–459.
16. Dekker, J., Marti-Renom, M.A., Mirny, L.A., Exploring the three-dimensional organization of genomes: interpreting chromatin interaction data. Nat. Rev. Genet. 14 (2013), 390–403.
17. Dilworth, F.J., Seaver, K.J., Fishburn, A.L., Htet, S.L., Tapscott, S.J., In vitro transcription system delineates the distinct roles of the coactivators pCAF and p300 during MyoD/E47-dependent transactivation. Proc. Natl. Acad. Sci. USA 101 (2004), 11593–11598.
18. ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature 489 (2012), 57–74.
19. Felsenfeld, G., Boyes, J., Chung, J., Clark, D., Studitsky, V., Chromatin structure and gene expression. Proc. Natl. Acad. Sci. USA 93 (1996), 9384–9388.
20. Fong, A.P., Yao, Z., Zhong, J.W., Johnson, N.M., Farr, G.H. 3rd, Maves, L., Tapscott, S.J., Conversion of MyoD to a neurogenic factor: binding site specificity determines lineage. Cell Rep. 10 (2015), 1937–1946.
21. Fortin, J.-P., Hansen, K.D., Reconstructing A/B compartments as revealed by Hi-C using long-range correlations in epigenetic data. Genome Biol., 16, 2015, 180.
22. Friedman, J., Hastie, T., Tibshirani, R., Sparse inverse covariance estimation with the graphical lasso. Biostatistics 9 (2008), 432–441.
23. Goldberg, A.V., and Tarjan, R.E. (1986). A new approach to the maximum flow problem. In Proceedings of the Eighteenth Annual ACM Symposium on Theory of Computing, pp. 136–146.
24. Grant, C.E., Bailey, T.L., Noble, W.S., FIMO: scanning for occurrences of a given motif. Bioinformatics 27 (2011), 1017–1018.
25. Gusev, A., Lee, S.H., Trynka, G., Finucane, H., Vilhjálmsson, B.J., Xu, H., Zang, C., Ripke, S., Bulik-Sullivan, B., Stahl, E., et al. Schizophrenia Working Group of the Psychiatric Genomics Consortium SWE-SCZ Consortium Schizophrenia Working Group of the Psychiatric Genomics Consortium, SWE-SCZ Consortium. Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases. Am. J. Hum. Genet. 95 (2014), 535–552.
26. Harrow, J., Frankish, A., Gonzalez, J.M., Tapanari, E., Diekhans, M., Kokocinski, F., Aken, B.L., Barrell, D., Zadissa, A., Searle, S., et al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res 22 (2012), 1760–1774.
27. Killick, R., Fearnhead, P., Eckley, I.A., Optimal detection of changepoints with a linear computational cost. J. Am. Stat. Assoc. 107 (2012), 1590–1598.
28. Knoepfler, P.S., Bergstrom, D.A., Uetsuki, T., Dac-Korytko, I., Sun, Y.H., Wright, W.E., Tapscott, S.J., Kamps, M.P., A conserved motif N-terminal to the DNA-binding domains of myogenic bHLH transcription factors mediates cooperative DNA binding with pbx-Meis1/Prep1. Nucleic Acids Res. 27 (1999), 3752–3761.
29. Krom, Y.D., Dumonceaux, J., Mamchaoui, K., den Hamer, B., Mariot, V., Negroni, E., Geng, L.N., Martin, N., Tawil, R., Tapscott, S.J., et al. Generation of isogenic D4Z4 contracted and noncontracted immortal muscle cell clones from a mosaic patient: a cellular model for FSHD. Am. J. Pathol. 181 (2012), 1387–1401.
30. Langmead, B., Salzberg, S.L., Fast gapped-read alignment with Bowtie 2. Nat. Methods 9 (2012), 357–359.
31. Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G., Abecasis, G., Durbin, R., 1000 Genome Project Data Processing Subgroup. The sequence alignment/map format and SAMtools. Bioinformatics 25 (2009), 2078–2079.
32. Libbrecht, M.W., Rodriguez, O., Weng, Z., Hoffman, M., Bilmes, J.A., Noble, W.S., A unified encyclopedia of human functional DNA elements through fully automated annotation of 164 human cell types. bioRxiv, 2016, 10.1101/086025.
33. Love, M.I., Huber, W., Anders, S., Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol., 15, 2014, 550.
34. MacQuarrie, K.L., Yao, Z., Fong, A.P., Diede, S.J., Rudzinski, E.R., Hawkins, D.S., Tapscott, S.J., Comparison of genome-wide binding of MyoD in normal human myogenic cells and rhabdomyosarcomas identifies regional and local suppression of promyogenic transcription factors. Mol. Cell. Biol. 33 (2013), 773–784.
35. Mathelier, A., Fornes, O., Arenillas, D.J., Chen, C.-Y., Denay, G., Lee, J., Shi, W., Shyr, C., Tan, G., Worsley-Hunt, R., et al. JASPAR 2016: a major expansion and update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. 44:D1 (2016), D110–D115.
36. Maurano, M.T., Humbert, R., Rynes, E., Thurman, R.E., Haugen, E., Wang, H., Reynolds, A.P., Sandstrom, R., Qu, H., Brody, J., et al. Systematic localization of common disease-associated variation in regulatory DNA. Science 337 (2012), 1190–1195.
37. Maves, L., Waskiewicz, A.J., Paul, B., Cao, Y., Tyler, A., Moens, C.B., Tapscott, S.J., Pbx homeodomain proteins direct Myod activity to promote fast-muscle differentiation. Development 134 (2007), 3371–3382.
38. Merico, D., Isserlin, R., Stueker, O., Emili, A., Bader, G.D., Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS ONE, 5, 2010, e13984.
39. Mifsud, B., Tavares-Cadete, F., Young, A.N., Sugar, R., Schoenfelder, S., Ferreira, L., Wingett, S.W., Andrews, S., Grey, W., Ewels, P.A., et al. Mapping long-range promoter contacts in human cells with high-resolution capture Hi-C. Nat. Genet. 47 (2015), 598–606.
40. Molkentin, J.D., Black, B.L., Martin, J.F., Olson, E.N., Cooperative activation of muscle gene expression by MEF2 and myogenic bHLH proteins. Cell 83 (1995), 1125–1136.
41. Philipot, O., Joliot, V., Ait-Mohamed, O., Pellentz, C., Robin, P., Fritsch, L., Ait-Si-Ali, S., The core binding factor CBF negatively regulates skeletal muscle terminal differentiation. PLoS ONE, 5, 2010, e9425.
42. Puri, P.L., Sartorelli, V., Yang, X.J., Hamamori, Y., Ogryzko, V.V., Howard, B.H., Kedes, L., Wang, J.Y., Graessmann, A., Nakatani, Y., Levrero, M., Differential roles of p300 and PCAF acetyltransferases in muscle differentiation. Mol. Cell 1 (1997), 35–45.
43. Qiu, X., Mao, Q., Tang, Y., Wang, L., Chawla, R., Pliner, H.A., Trapnell, C., Reversed graph embedding resolves complex single-cell trajectories. Nat. Methods 14 (2017), 979–982.
44. Qiu, X., Hill, A., Packer, J., Lin, D., Ma, Y.-A., Trapnell, C., Single-cell mRNA quantification and differential analysis with Census. Nat. Methods 14 (2017), 309–315.
45. Quinlan, A.R., Hall, I.M., BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26 (2010), 841–842.
46. Rao, S.S.P., Huntley, M.H., Durand, N.C., Stamenova, E.K., Bochkov, I.D., Robinson, J.T., Sanborn, A.L., Machol, I., Omer, A.D., Lander, E.S., Aiden, E.L., A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell 159 (2014), 1665–1680.
47. Sanborn, A.L., Rao, S.S.P., Huang, S.-C., Durand, N.C., Huntley, M.H., Jewett, A.I., Bochkov, I.D., Chinnappan, D., Cutkosky, A., Li, J., et al. Chromatin extrusion explains key features of loop and domain formation in wild-type and engineered genomes. Proc. Natl. Acad. Sci. USA 112 (2015), E6456–E6465.
48. Sartorelli, V., Huang, J., Hamamori, Y., Kedes, L., Molecular mechanisms of myogenic coactivation by p300: direct interaction with the activation domain of MyoD and with the MADS box of MEF2C. Mol. Cell. Biol. 17 (1997), 1010–1026.
49. Scrucca, L., Fop, M., Murphy, T.B., Raftery, A.E., mclust 5: clustering, classification and density estimation using Gaussian finite mixture models. R J. 8 (2016), 289–317.
50. Serra, C., Palacios, D., Mozzetta, C., Forcales, S.V., Morantte, I., Ripani, M., Jones, D.R., Du, K., Jhala, U.S., Simone, C., Puri, P.L., Functional interdependence at the chromatin level between the MKK6/p38 and IGF1/PI3K/AKT pathways during muscle differentiation. Mol. Cell 28 (2007), 200–213.
51. Simone, C., Forcales, S.V., Hill, D.A., Imbalzano, A.N., Latella, L., Puri, P.L., p38 pathway targets SWI-SNF chromatin-remodeling complex to muscle-specific loci. Nat. Genet. 36 (2004), 738–743.
52. Simpson, E.H., The interpretation of interaction in contingency tables. J. R. Stat. Soc. Series B Stat. Methodol. 13 (1951), 238–241.
53. Snider, L., Geng, L.N., Lemmers, R.J.L.F., Kyba, M., Ware, C.B., Nelson, A.M., Tawil, R., Filippova, G.N., van der Maarel, S.M., Tapscott, S.J., Miller, D.G., Facioscapulohumeral dystrophy: incomplete suppression of a retrotransposed gene. PLoS Genet., 6, 2010, e1001181.
54. Tang, Z., Luo, O.J., Li, X., Zheng, M., Zhu, J.J., Szalaj, P., Trzaskoma, P., Magalska, A., Wlodarczyk, J., Ruszczycki, B., et al. CTCF-mediated human 3D genome architecture reveals chromatin topology for transcription. Cell 163 (2015), 1611–1627.
55. Thurman, R.E., Rynes, E., Humbert, R., Vierstra, J., Maurano, M.T., Haugen, E., Sheffield, N.C., Stergachis, A.B., Wang, H., Vernot, B., et al. The accessible chromatin landscape of the human genome. Nature 489 (2012), 75–82.
56. Tolhuis, B., Palstra, R.-J., Splinter, E., Grosveld, F., de Laat, W., Looping and interaction between hypersensitive sites in the active β-globin locus. Mol. Cell 10 (2002), 1453–1465.
57. Trapnell, C., Defining cell types and states with single-cell genomics. Genome Res. 25 (2015), 1491–1498.
58. Trapnell, C., Cacchiarelli, D., Grimsby, J., Pokharel, P., Li, S., Morse, M., Lennon, N.J., Livak, K.J., Mikkelsen, T.S., Rinn, J.L., The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat. Biotechnol. 32 (2014), 381–386.
59. Väremo, L., Nielsen, J., Nookaew, I., Enriching the gene set analysis of genome-wide data by incorporating directionality of gene expression and combining statistical hypotheses and methods. Nucleic Acids Res. 41 (2013), 4378–4391.
60. Yoshida, N., Yoshida, S., Koishi, K., Masuda, K., Nabeshima, Y., Cell heterogeneity upon myogenic differentiation: down-regulation of MyoD and Myf-5 generates ‘reserve cells’. J. Cell Sci. 111 (1998), 769–779.
61. Zhang, Y., Liu, T., Meyer, C.A., Eeckhoute, J., Johnson, D.S., Bernstein, B.E., Nusbaum, C., Myers, R.M., Brown, M., Li, W., Liu, X.S., Model-based analysis of ChIP-seq (MACS). Genome Biol., 9, 2008, R137.
62. Zou, H., Hastie, T., Regularization and variable selection via the elastic net. J. R. Stat. Soc. Series B Stat. Methodol. 67 (2005), 301–320.