Recent technical advancements have facilitated the mapping of epigenomes at single-cell resolution; however, the throughput and quality of these methods have limited their widespread adoption. Here ...we describe a high-quality (10
nuclear fragments per cell) droplet-microfluidics-based method for single-cell profiling of chromatin accessibility. We use this approach, named 'droplet single-cell assay for transposase-accessible chromatin using sequencing' (dscATAC-seq), to assay 46,653 cells for the unbiased discovery of cell types and regulatory elements in adult mouse brain. We further increase the throughput of this platform by combining it with combinatorial indexing (dsciATAC-seq), enabling single-cell studies at a massive scale. We demonstrate the utility of this approach by measuring chromatin accessibility across 136,463 resting and stimulated human bone marrow-derived cells to reveal changes in the cis- and trans-regulatory landscape across cell types and under stimulatory conditions at single-cell resolution. Altogether, we describe a total of 510,123 single-cell profiles, demonstrating the scalability and flexibility of this droplet-based platform.
Widespread linkage disequilibrium and incomplete annotation of cell-to-cell state variation represent substantial challenges to elucidating mechanisms of trait-associated genetic variation. Here we ...perform genetic fine-mapping for blood cell traits in the UK Biobank to identify putative causal variants. These variants are enriched in genes encoding proteins in trait-relevant biological pathways and in accessible chromatin of hematopoietic progenitors. For regulatory variants, we explore patterns of developmental enhancer activity, predict molecular mechanisms, and identify likely target genes. In several instances, we localize multiple independent variants to the same regulatory element or gene. We further observe that variants with pleiotropic effects preferentially act in common progenitor populations to direct the production of distinct lineages. Finally, we leverage fine-mapped variants in conjunction with continuous epigenomic annotations to identify trait-cell type enrichments within closely related populations and in single cells. Our study provides a comprehensive framework for single-variant and single-cell analyses of genetic associations.
Regulatory networks that maintain functional, differentiated cell states are often dysregulated in tumor development. Here, we use single-cell epigenomics to profile chromatin state transitions in a ...mouse model of lung adenocarcinoma (LUAD). We identify an epigenomic continuum representing loss of cellular identity and progression toward a metastatic state. We define co-accessible regulatory programs and infer key activating and repressive chromatin regulators of these cell states. Among these co-accessibility programs, we identify a pre-metastatic transition, characterized by activation of RUNX transcription factors, which mediates extracellular matrix remodeling to promote metastasis and is predictive of survival across human LUAD patients. Together, these results demonstrate the power of single-cell epigenomics to identify regulatory programs to uncover mechanisms and key biomarkers of tumor progression.
Display omitted
•Cancer progression is marked by a continuum of heterogeneous epigenomic states•Lung cancer cells adopt features of other cell identities across tumor evolution•RUNX2 transcription factor activity is associated with a pre-metastatic cell state•Regulatory programs defined from mouse models are useful predictive biomarkers
Using a generalizable framework to leverage single-cell chromatin accessibility data to investigate cell state transitions across tumor evolution in a mouse model of lung adenocarcinoma, LaFave et al. elucidate a pre-metastatic cell state that arises in primary tumors prior to metastasis.
Transcriptional dysregulation has long been recognized as central to the pathogenesis of Huntington's disease (HD). MicroRNAs (miRNAs) represent a major system of post-transcriptional regulation, by ...either preventing translational initiation or by targeting transcripts for storage or for degradation. Using next-generation miRNA sequencing in prefrontal cortex (Brodmann Area 9) of twelve HD and nine controls, we identified five miRNAs (miR-10b-5p, miR-196a-5p, miR-196b-5p, miR-615-3p and miR-1247-5p) up-regulated in HD at genome-wide significance (FDR q-value<0.05). Three of these, miR-196a-5p, miR-196b-5p and miR-615-3p, were expressed at near zero levels in control brains. Expression was verified for all five miRNAs using reverse transcription quantitative PCR and all but miR-1247-5p were replicated in an independent sample (8HD/8C). Ectopic miR-10b-5p expression in PC12 HTT-Q73 cells increased survival by MTT assay and cell viability staining suggesting increased expression may be a protective response. All of the miRNAs but miR-1247-5p are located in intergenic regions of Hox clusters. Total mRNA sequencing in the same samples identified fifteen of 55 genes within the Hox cluster gene regions as differentially expressed in HD, and the Hox genes immediately adjacent to the four Hox cluster miRNAs as up-regulated. Pathway analysis of mRNA targets of these miRNAs implicated functions for neuronal differentiation, neurite outgrowth, cell death and survival. In regression models among the HD brains, huntingtin CAG repeat size, onset age and age at death were independently found to be inversely related to miR-10b-5p levels. CAG repeat size and onset age were independently inversely related to miR-196a-5p, onset age was inversely related to miR-196b-5p and age at death was inversely related to miR-615-3p expression. These results suggest these Hox-related miRNAs may be involved in neuroprotective response in HD. Recently, miRNAs have shown promise as biomarkers for human diseases and given their relationship to disease expression, these miRNAs are biomarker candidates in HD.
Carcinoma associated fibroblasts (CAFs) form the main constituents of tumor stroma and play an important role in tumor growth and invasion. The presence of CAFs is a strong predictor of poor ...prognosis of head and neck squamous cell carcinoma. Despite significant progress in determining the role of CAFs in tumor progression, the mechanisms contributing to their activation remain poorly characterized, in part due to fibroblast heterogeneity and the scarcity of reliable fibroblast surface markers. To search for such markers in oral squamous cell carcinoma (OSCC), we applied a novel approach that uses RNA-sequencing data derived from the cancer genome atlas (TCGA). Specifically, our strategy allowed for an unbiased identification of genes whose expression was closely associated with a set of bona fide stroma-specific transcripts, namely the interstitial collagens COL1A1, COL1A2, and COL3A1. Among the top hits were genes involved in cellular matrix remodeling and tumor invasion and migration, including platelet-derived growth factor receptor beta (PDGFRβ), which was found to be the highest-ranking receptor protein genome-wide. Similar analyses performed on ten additional TCGA cancer datasets revealed that other tumor types shared CAF markers with OSCC, including PDGFRβ, which was found to significantly correlate with the reference collagen expression in ten of the 11 cancer types tested. Subsequent immunostaining of OSCC specimens demonstrated that PDGFRβ was abundantly expressed in stromal fibroblasts of all tested cases (12/12), while it was absent in tumor cells, with greater specificity than other known markers such as alpha smooth muscle actin or podoplanin (3/11). Overall, this study identified PDGFRβ as a novel marker of stromal activation in OSCC, and further characterized a list of promising candidate CAF markers that may be relevant to other carcinomas. Our novel approach provides for a fast and accurate method to identify CAF markers without the need for large-scale immunostaining experiments.
Oral squamous cell carcinoma (OSCC) is a prevalent form of cancer that develops from the epithelium of the oral cavity. OSCC is on the rise worldwide, and death rates associated with the disease are ...particularly high. Despite progress in understanding the mutational and expression landscape associated with OSCC, advances in deciphering these alterations for the development of therapeutic strategies have been limited. Further insight into the molecular cues that contribute to OSCC is therefore required. Here, we show that the transcriptional regulators YAP (YAP1) and TAZ (WWTR1), which are key effectors of the Hippo pathway, drive protumorigenic signals in OSCC. Regions of premalignant oral tissues exhibit aberrant nuclear YAP accumulation, suggesting that dysregulated YAP activity contributes to the onset of OSCC. Supporting this premise, we determined that nuclear YAP and TAZ activity drives OSCC cell proliferation, survival, and migration in vitro, and is required for OSCC tumor growth and metastasis in vivo. Global gene expression profiles associated with YAP and TAZ knockdown revealed changes in the control of gene expression implicated in protumorigenic signaling, including those required for cell cycle progression and survival. Notably, the transcriptional signature regulated by YAP and TAZ significantly correlates with gene expression changes occurring in human OSCCs identified by The Cancer Genome Atlas (TCGA), emphasizing a central role for YAP and TAZ in OSCC biology.
This study defines a YAP/TAZ-regulated transcriptional program in OSCC and reveals novel roles for nuclear YAP/TAZ activity in the onset and progression of this cancer.
Inflammation can trigger lasting phenotypes in immune and non-immune cells. Whether and how human infections and associated inflammation can form innate immune memory in hematopoietic stem and ...progenitor cells (HSPC) has remained unclear. We found that circulating HSPC, enriched from peripheral blood, captured the diversity of bone marrow HSPC, enabling investigation of their epigenomic reprogramming following coronavirus disease 2019 (COVID-19). Alterations in innate immune phenotypes and epigenetic programs of HSPC persisted for months to 1 year following severe COVID-19 and were associated with distinct transcription factor (TF) activities, altered regulation of inflammatory programs, and durable increases in myelopoiesis. HSPC epigenomic alterations were conveyed, through differentiation, to progeny innate immune cells. Early activity of IL-6 contributed to these persistent phenotypes in human COVID-19 and a mouse coronavirus infection model. Epigenetic reprogramming of HSPC may underlie altered immune function following infection and be broadly relevant, especially for millions of COVID-19 survivors.
Display omitted
•Severe COVID-19 programs durable epigenetic changes and hyper-activation in monocytes•Circulating HSPC capture post-COVID-19 changes in hematopoiesis and stem cell programs•Post-COVID-19 HSPC convey epigenetic and transcriptional memory to mature progeny cells•IL-6 contributes to epigenetic reprogramming of mouse and human HSPC and myeloid cells
Severe COVID-19 can reprogram hematopoiesis and establish epigenetic memory in hematopoietic stem and progenitor cells (HSPC) and progeny myeloid cells for up to 1 year. These durable alterations, which could affect post-infection immune responses and equilibrium, are controlled in part by the activity of IL-6 during acute disease.
Cells require coordinated control over gene expression when responding to environmental stimuli. Here we apply scATAC-seq and single-cell RNA sequencing (scRNA-seq) in resting and stimulated human ...blood cells. Collectively, we generate ∼91,000 single-cell profiles, allowing us to probe the cis-regulatory landscape of the immunological response across cell types, stimuli, and time. Advancing tools to integrate multi-omics data, we develop functional inference of gene regulation (FigR), a framework to computationally pair scATAC-seq with scRNA-seq cells, connect distal cis-regulatory elements to genes, and infer gene-regulatory networks (GRNs) to identify candidate transcription factor (TF) regulators. Utilizing these paired multi-omics data, we define domains of regulatory chromatin (DORCs) of immune stimulation and find that cells alter chromatin accessibility and gene expression at timescales of minutes. Construction of the stimulation GRN elucidates TF activity at disease-associated DORCs. Overall, FigR enables elucidation of regulatory interactions across single-cell data, providing new opportunities to understand the function of cells within tissues.
Display omitted
•A multi-omics resource of resting and immunologically stimulated human blood•Data integration identifies expression-linked regulatory elements across cells•FigR reveals TF activator and repressor drivers of chromatin accessibility•Gene-regulatory modeling identifies disease-associated regulatory interactions
Single-cell methods for measuring chromatin accessibility (ATAC-seq) and gene expression (RNA-seq) are rapidly evolving, but tools to integrate data and infer gene-regulatory relationships remain limited. Here we generate multi-omics data of resting and stimulated human blood cells and present a new computational framework for constructing gene-regulatory networks (GRNs). Specifically, we describe functional inference of gene regulation (FigR), a workflow to (1) pair scATAC-seq with scRNA-seq, (2) connect cis-regulatory elements to target genes, and (3) identify TF-gene relationships.
Head and neck squamous cell carcinoma (HNSCC) is an aggressive malignancy characterized by tumor heterogeneity, locoregional metastases, and resistance to existing treatments. Although a number of ...genomic and molecular alterations associated with HNSCC have been identified, they have had limited impact on the clinical management of this disease. To date, few targeted therapies are available for HNSCC, and only a small fraction of patients have benefited from these treatments. A frequent feature of HNSCC is the inappropriate activation of β-catenin that has been implicated in cell survival and in the maintenance and expansion of stem cell-like populations, thought to be the underlying cause of tumor recurrence and resistance to treatment. However, the therapeutic value of targeting β-catenin activity in HNSCC has not been explored.
We utilized a combination of computational and experimental profiling approaches to examine the effects of blocking the interaction between β-catenin and cAMP-responsive element binding (CREB)-binding protein (CBP) using the small molecule inhibitor ICG-001. We generated and annotated in vitro treatment gene expression signatures of HNSCC cells, derived from human oral squamous cell carcinomas (OSCCs), using microarrays. We validated the anti-tumorigenic activity of ICG-001 in vivo using SCC-derived tumor xenografts in murine models, as well as embryonic zebrafish-based screens of sorted stem cell-like subpopulations. Additionally, ICG-001-inhibition signatures were overlaid with RNA-sequencing data from The Cancer Genome Atlas (TCGA) for human OSCCs to evaluate its association with tumor progression and prognosis.
ICG-001 inhibited HNSCC cell proliferation and tumor growth in cellular and murine models, respectively, while promoting intercellular adhesion and loss of invasive phenotypes. Furthermore, ICG-001 preferentially targeted the ability of subpopulations of stem-like cells to establish metastatic tumors in zebrafish. Significantly, interrogation of the ICG-001 inhibition-associated gene expression signature in the TCGA OSCC human cohort indicated that the targeted β-catenin/CBP transcriptional activity tracked with tumor status, advanced tumor grade, and poor overall patient survival.
Collectively, our results identify β-catenin/CBP interaction as a novel target for anti-HNSCC therapy and provide evidence that derivatives of ICG-001 with enhanced inhibitory activity may serve as an effective strategy to interfere with aggressive features of HNSCC.
Lysine-specific demethylase 1 (LSD1) is a nuclear histone demethylase and a member of the amine oxidase (AO) family. LSD1 is a flavin-containing AO that specifically catalyzes the demethylation of ...mono- and di-methylated histone H3 lysine 4 through an FAD-dependent oxidative reaction. LSD1 is inappropriately upregulated in lung, liver, brain and esophageal cancers, where it promotes cancer initiation, progression, and metastasis. However, unlike other lysine-specific demethylases, the role and specific targets of LSD1 in oral squamous cell carcinoma (OSCC) pathogenesis remain unknown. We show that LSD1 protein expression was increased in malignant OSCC tissues in a clinical tissue microarray, and its expression correlated with progressive tumor stages. In an orthotopic oral cancer mouse model, LSD1 overexpression in aggressive HSC-3 cells promoted metastasis whereas knockdown of LSD1 inhibited tumor spread, suggesting that LSD1 is a key regulator of OSCC metastasis. Pharmacological inhibition of LSD1 using a specific small molecule inhibitor, GSK-LSD1, down-regulated EGF signaling pathway. Further, GSK-LSD1 attenuates CTGF/CCN2, MMP13, LOXL4 and vimentin expression but increased E-cadherin expression in pre-existing, patient-derived tonsillar OSCC xenografts. Similarly, GSK-LSD1 inhibited proliferation and CTGF expression in mesenchymal cells, including myoepithelial cells and osteosarcoma cells. In addition, gene set enrichment analysis revealed that GSK-LSD1 increased p53 expression and apoptosis while inhibiting c-myc, β-catenin and YAP-induced oncogenic transcriptional networks. These data reveal that aberrant LSD1 activation regulates key OSCC microenvironment and EMT promoting factors, including CTGF, LOXL4 and MMP13.