Epigenetic memory of coronavirus infection in innate immune cells and their progenitors.

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.

Chromatin Potential Identified by Shared Single-Cell Profiling of RNA and Chromatin.

Cell differentiation and function are regulated across multiple layers of gene regulation, including modulation of gene expression by changes in chromatin accessibility. However, differentiation is an asynchronous process precluding a temporal understanding of regulatory events leading to cell fate commitment. Here we developed simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq), a highly scalable approach for measurement of chromatin accessibility and gene expression in the same single cell, applicable to different tissues. Using 34,774 joint profiles from mouse skin, we develop a computational strategy to identify cis-regulatory interactions and define domains of regulatory chromatin (DORCs) that significantly overlap with super-enhancers. During lineage commitment, chromatin accessibility at DORCs precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for lineage commitment. We computationally infer chromatin potential as a quantitative measure of chromatin lineage-priming and use it to predict cell fate outcomes. SHARE-seq is an extensible platform to study regulatory circuitry across diverse cells in tissues.

MicroRNAs located in the Hox gene clusters are implicated in huntington's disease pathogenesis.

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.

Single-cell multi-scale footprinting reveals the modular organization of DNA regulatory elements.

Cis-regulatory elements control gene expression and are dynamic in their structure, reflecting changes to the composition of diverse effector proteins over time1-3. Here we sought to connect the structural changes at cis-regulatory elements to alterations in cellular fate and function. To do this we developed PRINT, a computational method that uses deep learning to correct sequence bias in chromatin accessibility data and identifies multi-scale footprints of DNA-protein interactions. We find that multi-scale footprints enable more accurate inference of TF and nucleosome binding. Using PRINT with single-cell multi-omics, we discover wide-spread changes to the structure and function of candidate cis-regulatory elements (cCREs) across hematopoiesis, wherein nucleosomes slide, expose DNA for TF binding, and promote gene expression. Activity segmentation using the co-variance across cell states identifies "sub-cCREs" as modular cCRE subunits of regulatory DNA. We apply this single-cell and PRINT approach to characterize the age-associated alterations to cCREs within hematopoietic stem cells (HSCs). Remarkably, we find a spectrum of aging alterations among HSCs corresponding to a global gain of sub-cCRE activity while preserving cCRE accessibility. Collectively, we reveal the functional importance of cCRE structure across cell states, highlighting changes to gene regulation at single-cell and single-base-pair resolution.

Functional inference of gene regulation using single-cell multi-omics.

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 scA-TAC-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.

Smarca4 Inactivation Promotes Lineage-Specific Transformation and Early Metastatic Features in the Lung.

SMARCA4/BRG1 encodes for one of two mutually exclusive ATPases present in mammalian SWI/SNF chromatin remodeling complexes and is frequently mutated in human lung adenocarcinoma. However, the functional consequences of SMARCA4 mutation on tumor initiation, progression, and chromatin regulation in lung cancer remain poorly understood. Here, we demonstrate that loss of Smarca4 sensitizes club cell secretory protein-positive cells within the lung in a cell type-dependent fashion to malignant transformation and tumor progression, resulting in highly advanced dedifferentiated tumors and increased metastatic incidence. Consistent with these phenotypes, Smarca4-deficient primary tumors lack lung lineage transcription factor activities and resemble a metastatic cell state. Mechanistically, we show that Smarca4 loss impairs the function of all three classes of SWI/SNF complexes, resulting in decreased chromatin accessibility at lung lineage motifs and ultimately accelerating tumor progression. Thus, we propose that the SWI/SNF complex via Smarca4 acts as a gatekeeper for lineage-specific cellular transformation and metastasis during lung cancer evolution. SIGNIFICANCE: We demonstrate cell-type specificity in the tumor-suppressive functions of SMARCA4 in the lung, pointing toward a critical role of the cell-of-origin in driving SWI/SNF-mutant lung adenocarcinoma. We further show the direct effects of SMARCA4 loss on SWI/SNF function and chromatin regulation that cause aggressive malignancy during lung cancer evolution.This article is highlighted in the In This Issue feature, p. 275.

ß-Catenin/CBP inhibition alters epidermal growth factor receptor fucosylation status in oral squamous cell carcinoma.

Epidermal growth factor receptor (EGFR) is a major driver of head and neck cancer, a devastating malignancy with a major sub-site in the oral cavity manifesting as oral squamous cell carcinoma (OSCC). EGFR is a glycoprotein receptor tyrosine kinase (RTK) whose activity is upregulated in >80% OSCC. Current anti-EGFR therapy relies on the use of cetuximab, a monoclonal antibody against EGFR, although it has had only a limited response in patients. Here, we uncover a novel mechanism regulating EGFR activity, identifying a role of the nuclear branch of the Wnt/ß-catenin signaling pathway, the ß-catenin/CBP axis, in control of post-translational modification of N-glycans on the EGFR. Genomic and structural analyses reveal that ß-catenin/CBP signaling represses fucosylation on the antennae of N-linked glycans on EGFR. By employing nUPLC-MS/MS, we determined that malignant human OSCC cells harbor EGFR with a paucity of N-glycan antennary fucosylation, while indolent cells display higher levels of fucosylation at sites N420 and N579. Additionally, treatment with either ICG-001 or E7386, which are both small molecule inhibitors of ß-catenin/CBP signaling, leads to increased transcriptional expression of fucosyltransferases FUT2 and FUT3, with a concomitant increase in EGFR N-glycan antennary fucosylation. In order to discover which fucosylated glycan epitopes are involved in the observed effect, we performed in-depth characterization of multiply-fucosylated N-glycans via tandem mass spectrometry analysis of the EGFR tryptic glycopeptides. Data are available via ProteomeXchange with identifier PXD017060. We propose that ß-catenin/CBP signaling promotes EGFR oncogenic activity in OSCC by inhibiting its N-glycan antennary fucosylation through transcriptional repression of FUT2 and FUT3.

Epigenomic State Transitions Characterize Tumor Progression in Mouse Lung Adenocarcinoma.

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.

CaDrA: A Computational Framework for Performing Candidate Driver Analyses Using Genomic Features.

The identification of genetic alteration combinations as drivers of a given phenotypic outcome, such as drug sensitivity, gene or protein expression, and pathway activity, is a challenging task that is essential to gaining new biological insights and to discovering therapeutic targets. Existing methods designed to predict complementary drivers of such outcomes lack analytical flexibility, including the support for joint analyses of multiple genomic alteration types, such as somatic mutations and copy number alterations, multiple scoring functions, and rigorous significance and reproducibility testing procedures. To address these limitations, we developed Candidate Driver Analysis or CaDrA, an integrative framework that implements a step-wise heuristic search approach to identify functionally relevant subsets of genomic features that, together, are maximally associated with a specific outcome of interest. We show CaDrA's overall high sensitivity and specificity for typically sized multi-omic datasets using simulated data, and demonstrate CaDrA's ability to identify known mutations linked with sensitivity of cancer cells to drug treatment using data from the Cancer Cell Line Encyclopedia (CCLE). We further apply CaDrA to identify novel regulators of oncogenic activity mediated by Hippo signaling pathway effectors YAP and TAZ in primary breast cancer tumors using data from The Cancer Genome Atlas (TCGA), which we functionally validate in vitro. Finally, we use pan-cancer TCGA protein expression data to show the high reproducibility of CaDrA's search procedure. Collectively, this work demonstrates the utility of our framework for supporting the fast querying of large, publicly available multi-omics datasets, including but not limited to TCGA and CCLE, for potential drivers of a given target profile of interest.

Droplet-based combinatorial indexing for massive-scale single-cell chromatin accessibility.

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 (105 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.

Interrogation of human hematopoiesis at single-cell and single-variant resolution.

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.

Functional and genomic analyses reveal therapeutic potential of targeting ß-catenin/CBP activity in head and neck cancer.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Inhibition of LSD1 epigenetically attenuates oral cancer growth and metastasis.

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.

PDGFRß Is a Novel Marker of Stromal Activation in Oral Squamous Cell Carcinomas.

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.

A YAP/TAZ-Regulated Molecular Signature Is Associated with Oral Squamous Cell Carcinoma.

UNLABELLED: 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. IMPLICATIONS: 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.

miR-10b-5p expression in Huntington's disease brain relates to age of onset and the extent of striatal involvement.

BACKGROUND: MicroRNAs (miRNAs) are small non-coding RNAs that recognize sites of complementarity of target messenger RNAs, resulting in transcriptional regulation and translational repression of target genes. In Huntington's disease (HD), a neurodegenerative disease caused by a trinucleotide repeat expansion, miRNA dyregulation has been reported, which may impact gene expression and modify the progression and severity of HD. METHODS: We performed next-generation miRNA sequence analysis in prefrontal cortex (Brodmann Area 9) from 26 HD, 2 HD gene positive, and 36 control brains. Neuropathological information was available for all HD brains, including age at disease onset, CAG-repeat size, Vonsattel grade, and Hadzi-Vonsattel striatal and cortical scores, a continuous measure of the extent of neurodegeneration. Linear models were performed to examine the relationship of miRNA expression to these clinical features, and messenger RNA targets of associated miRNAs were tested for gene ontology term enrichment. RESULTS: We identified 75 miRNAs differentially expressed in HD brain (FDR q-value <0.05). Among the HD brains, nine miRNAs were significantly associated with Vonsattel grade of neuropathological involvement and three of these, miR-10b-5p, miR-10b-3p, and miR-302a-3p, significantly related to the Hadzi-Vonsattel striatal score (a continuous measure of striatal involvement) after adjustment for CAG length. Five miRNAs (miR-10b-5p, miR-196a-5p, miR-196b-5p, miR-10b-3p, and miR-106a-5p) were identified as having a significant relationship to CAG length-adjusted age of onset including miR-10b-5p, the mostly strongly over-expressed miRNA in HD cases. Although prefrontal cortex was the source of tissue profiled in these studies, the relationship of miR-10b-5p expression to striatal involvement in the disease was independent of cortical involvement. Correlation of miRNAs to the clinical features clustered by direction of effect and the gene targets of the observed miRNAs showed association to processes relating to nervous system development and transcriptional regulation. CONCLUSIONS: These results demonstrate that miRNA expression in cortical BA9 provides insight into striatal involvement and support a role for these miRNAs, particularly miR-10b-5p, in HD pathogenicity. The miRNAs identified in our studies of postmortem brain tissue may be detectable in peripheral fluids and thus warrant consideration as accessible biomarkers for disease stage, rate of progression, and other important clinical characteristics of HD.