Background: Rupture and erosion of advanced atherosclerotic lesions with a resultant myocardial infarction or stroke are the leading worldwide cause of death. However, we have a limited understanding of the identity, origin, and function of many cells that make up late-stage atherosclerotic lesions, as well as the mechanisms by which they control plaque stability. Methods: We conducted a comprehensive single-cell RNA sequencing of advanced human carotid endarterectomy samples and compared these with single-cell RNA sequencing from murine microdissected advanced atherosclerotic lesions with smooth muscle cell (SMC) and endothelial lineage tracing to survey all plaque cell types and rigorously determine their origin. We further used chromatin immunoprecipitation sequencing (ChIP-seq), bulk RNA sequencing, and an innovative dual lineage tracing mouse to understand the mechanism by which SMC phenotypic transitions affect lesion pathogenesis. Results: We provide evidence that SMC-specific Klf4- versus Oct4-knockout showed virtually opposite genomic signatures, and their putative target genes play an important role regulating SMC phenotypic changes. Single-cell RNA sequencing revealed remarkable similarity of transcriptomic clusters between mouse and human lesions and extensive plasticity of SMC- and endothelial cell-derived cells including 7 distinct clusters, most negative for traditional markers. In particular, SMC contributed to a Myh11 ^- , Lgals3 ⁺ population with a chondrocyte-like gene signature that was markedly reduced with SMC- Klf4 knockout. We observed that SMCs that activate Lgals3 compose up to two thirds of all SMC in lesions. However, initial activation of Lgals3 in these cells does not represent conversion to a terminally differentiated state, but rather represents transition of these cells to a unique stem cell marker gene–positive, extracellular matrix-remodeling, “pioneer” cell phenotype that is the first to invest within lesions and subsequently gives rise to at least 3 other SMC phenotypes within advanced lesions, including Klf4-dependent osteogenic phenotypes likely to contribute to plaque calcification and plaque destabilization. Conclusions: Taken together, these results provide evidence that SMC-derived cells within advanced mouse and human atherosclerotic lesions exhibit far greater phenotypic plasticity than generally believed, with Klf4 regulating transition to multiple phenotypes including Lgals3 ⁺ osteogenic cells likely to be detrimental for late-stage atherosclerosis plaque pathogenesis.

We have devised a method for isolating virtually pure and comprehensive libraries of restriction fragments that contained replication initiation sites (bubbles) in vivo. We have now sequenced and mapped the bubble-containing fragments from GM06990, a near-normal EBV-transformed lymphoblastoid cell line, and have compared origin distributions with a comprehensive replication timing study recently published for this cell line. We find that early-firing origins, which represent ∼32% of all origins, overwhelmingly represent zones, associate only marginally with active transcription units, are localized within large domains of open chromatin, and are significantly associated with DNase I hypersensitivity. Origin “density” falls from early- to mid-S-phase, but rises again in late S-phase to levels only 17% lower than in early S-phase. Unexpectedly, late origin density calculated on the 1-Mb scale increases as a function of increasing chromatin compaction. Furthermore, the median efficiency of origins in late-replicating, heterochromatic domains is only 25% lower than in early-replicating euchromatic loci. Thus, the activation of early- and late-firing origins must be regulated by quintessentially different mechanisms. The aggregate data can be unified into a model in which initiation site selection is driven almost entirely by epigenetic factors that fashion both the long-range and local chromatin environments, with underlying DNA sequence and local transcriptional activity playing only minor roles. Importantly, the comprehensive origin map we have prepared for GM06990 overlaps moderately well with origin maps recently reported for the genomes of four different human cell lines based on the distributions of small nascent strands.

Capturing Binding Location and Speed: Transcription factor–binding sites in chromatin can be mapped by the chromatin immunoprecipitation (ChIP) assay, which analyzes formaldehyde-fixed chromatin fragments obtained from cells. However, the standard ChIP assay does not provide information about how stable the inter-actions are. Other approaches, including live-cell imaging, can reveal aspects of the dynamic behavior of transcription factors but are limited either in location precision or time resolution. Poorey et al. (p. 369 , published online 3 October) developed a model to explain how the ChIP signal relates to formaldehyde cross-linking time, and they developed a method to measure chromatin site–specific binding dynamics with high temporal resolution.

We have used a novel bubble-trapping procedure to construct nearly pure and comprehensive human origin libraries from early S- and log-phase HeLa cells, and from log-phase GM06990, a karyotypically normal lymphoblastoid cell line. When hybridized to ENCODE tiling arrays, these libraries illuminated 15.3%, 16.4%, and 21.8% of the genome in the ENCODE regions, respectively. Approximately half of the origin fragments cluster into zones, and their signals are generally higher than those of isolated fragments. Interestingly, initiation events are distributed about equally between genic and intergenic template sequences. While only 13.2% and 14.0% of genes within the ENCODE regions are actually transcribed in HeLa and GM06990 cells, 54.5% and 25.6% of zonal origin fragments overlap transcribed genes, most with activating chromatin marks in their promoters. Our data suggest that cell synchronization activates a significant number of inchoate origins. In addition, HeLa and GM06990 cells activate remarkably different origin populations. Finally, there is only moderate concordance between the log-phase HeLa bubble map and published maps of small nascent strands for this cell line.

Background: In the last decade, biochemical studies have revealed that epigenetic modifications including histone modifications, histone variants and DNA methylation form a complex network that regulate the state of chromatin and processes that depend on it including transcription and DNA replication. Currently, a large number of these epigenetic modifications are being mapped in a variety of cell lines at different stages of development using high throughput sequencing by members of the ENCODE consortium, the NIH Roadmap Epigenomics Program and the Human Epigenome Project. An extremely promising and underexplored area of research is the application of machine learning methods, which are designed to construct predictive network models, to these large-scale epigenomic data sets.Results: Using a ChIP-Seq data set of 20 histone lysine and arginine methylations and histone variant H2A.Z in human CD4⁺T-cells, we built predictive models of gene expression as a function of histone modification/variant levels using Multilinear (ML) Regression and Multivariate Adaptive Regression Splines (MARS). Along with extensive crosstalk among the 20 histone methylations, we found H4R3me2 was the most and second most globally repressive histone methylation among the 20 studied in the ML and MARS models, respectively. In support of our finding, a number of experimental studies show that PRMT5-catalyzed symmetric dimethylation of H4R3 is associated with repression of gene expression. This includes a recent study, which demonstrated that H4R3me2 is required for DNMT3A-mediated DNA methylation--a known global repressor of gene expression.Conclusion: In stark contrast to univariate analysis of the relationship between H4R3me2 and gene expression levels, our study showed that the regulatory role of some modifications like H4R3me2 is masked by confounding variables, but can be elucidated by multivariate/systems-level approaches.

Soluble growth factors and cytokines within the tumor microenvironment aid in the induction of the epithelial-to-mesenchymal transition (EMT). Although EMT promotes the development of cancer-initiating cells (CIC), cellular mechanisms by which cancer cells maintain mesenchymal phenotypes remain poorly understood. Work presented here indicates that induction of EMT stimulates non–small cell lung cancer (NSCLC) to secrete soluble factors that function in an autocrine fashion. Using gene expression profiling of all annotated and predicted secreted gene products, we find that NF-κB activity is required to upregulate INHBA/Activin, a morphogen in the TGFβ superfamily. INHBA is capable of inducing and maintaining mesenchymal phenotypes, including the expression of EMT master-switch regulators and self-renewal factors that sustain CIC phenotypes and promote lung metastasis. Our work demonstrates that INHBA mRNA and protein expression are commonly elevated in primary human NSCLC and provide evidence that INHBA is a critical autocrine factor that maintains mesenchymal properties of CICs to promote metastasis in NSCLC. Cancer Res; 75(2); 426–35. ©2014 AACR.