Network-based prioritization and validation of regulators of vascular smooth muscle cell proliferation in disease.

Aberrant vascular smooth muscle cell (VSMC) homeostasis and proliferation characterize vascular diseases causing heart attack and stroke. Here we elucidate molecular determinants governing VSMC proliferation by reconstructing gene regulatory networks from single-cell transcriptomics and epigenetic profiling. We detect widespread activation of enhancers at disease-relevant loci in proliferation-predisposed VSMCs. We compared gene regulatory network rewiring between injury-responsive and nonresponsive VSMCs, which suggested shared transcription factors but differing target loci between VSMC states. Through in silico perturbation analysis, we identified and prioritized previously unrecognized regulators of proliferation, including RUNX1 and TIMP1. Moreover, we showed that the pioneer transcription factor RUNX1 increased VSMC responsiveness and that TIMP1 feeds back to promote VSMC proliferation through CD74-mediated STAT3 signaling. Both RUNX1 and the TIMP1-CD74 axis were expressed in human VSMCs, showing low levels in normal arteries and increased expression in disease, suggesting clinical relevance and potential as vascular disease targets.

Cellular mechanisms of oligoclonal vascular smooth muscle cell expansion in cardiovascular disease.

AIMS: Quiescent, differentiated adult vascular smooth muscle cells (VSMCs) can be induced to proliferate and switch phenotype. Such plasticity underlies blood vessel homeostasis and contributes to vascular disease development. Oligoclonal VSMC contribution is a hallmark of end-stage vascular disease. Here, we aim to understand cellular mechanisms underpinning generation of this VSMC oligoclonality. METHODS AND RESULTS: We investigate the dynamics of VSMC clone formation using confocal microscopy and single-cell transcriptomics in VSMC-lineage-traced animal models. We find that activation of medial VSMC proliferation occurs at low frequency after vascular injury and that only a subset of expanding clones migrate, which together drives formation of oligoclonal neointimal lesions. VSMC contribution in small atherosclerotic lesions is typically from one or two clones, similar to observations in mature lesions. Low frequency (<0.1%) of clonal VSMC proliferation is also observed in vitro. Single-cell RNA-sequencing revealed progressive cell state changes across a contiguous VSMC population at onset of injury-induced proliferation. Proliferating VSMCs mapped selectively to one of two distinct trajectories and were associated with cells showing extensive phenotypic switching. A proliferation-associated transitory state shared pronounced similarities with atypical SCA1+ VSMCs from uninjured mouse arteries and VSMCs in healthy human aorta. We show functionally that clonal expansion of SCA1+ VSMCs from healthy arteries occurs at higher rate and frequency compared with SCA1- cells. CONCLUSION: Our data suggest that activation of proliferation at low frequency is a general, cell-intrinsic feature of VSMCs. We show that rare VSMCs in healthy arteries display VSMC phenotypic switching akin to that observed in pathological vessel remodelling and that this is a conserved feature of mouse and human healthy arteries. The increased proliferation of modulated VSMCs from healthy arteries suggests that these cells respond more readily to disease-inducing cues and could drive oligoclonal VSMC expansion.

Brief encounters: The relationship between enhancer proximity and gene expression.

In a recent issue of Nature, Zuin et al. (2022) integrate a transcriptional enhancer into various genomic locations surrounding a reporter gene promoter in order to determine how enhancer-promoter proximity affects gene expression. Their findings suggest a non-linear relationship, which the authors explain via a quantitative biophysical model with intriguing biological implications.

Prioritisation of Candidate Genes Underpinning COVID-19 Host Genetic Traits Based on High-Resolution 3D Chromosomal Topology.

Genetic variants showing associations with specific biological traits and diseases detected by genome-wide association studies (GWAS) commonly map to non-coding DNA regulatory regions. Many of these regions are located considerable distances away from the genes they regulate and come into their proximity through 3D chromosomal interactions. We previously developed COGS, a statistical pipeline for linking GWAS variants with their putative target genes based on 3D chromosomal interaction data arising from high-resolution assays such as Promoter Capture Hi-C (PCHi-C). Here, we applied COGS to COVID-19 Host Genetic Consortium (HGI) GWAS meta-analysis data on COVID-19 susceptibility and severity using our previously generated PCHi-C results in 17 human primary cell types and SARS-CoV-2-infected lung carcinoma cells. We prioritise 251 genes putatively associated with these traits, including 16 out of 47 genes highlighted by the GWAS meta-analysis authors. The prioritised genes are expressed in a broad array of tissues, including, but not limited to, blood and brain cells, and are enriched for genes involved in the inflammatory response to viral infection. Our prioritised genes and pathways, in conjunction with results from other prioritisation approaches and targeted validation experiments, will aid in the understanding of COVID-19 pathology, paving the way for novel treatments.

Transcriptional enhancers and their communication with gene promoters.

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

Detecting chromosomal interactions in Capture Hi-C data with CHiCAGO and companion tools.

Capture Hi-C is widely used to obtain high-resolution profiles of chromosomal interactions involving, at least on one end, regions of interest such as gene promoters. Signal detection in Capture Hi-C data is challenging and cannot be adequately accomplished with tools developed for other chromosome conformation capture methods, including standard Hi-C. Capture Hi-C Analysis of Genomic Organization (CHiCAGO) is a computational pipeline developed specifically for Capture Hi-C analysis. It implements a statistical model accounting for biological and technical background components, as well as bespoke normalization and multiple testing procedures for this data type. Here we provide a step-by-step guide to the CHiCAGO workflow that is aimed at users with basic experience of the command line and R. We also describe more advanced strategies for tuning the key parameters for custom experiments and provide guidance on data preprocessing and downstream analysis using companion tools. In a typical experiment, CHiCAGO takes ~2-3 h to run, although pre- and postprocessing steps may take much longer.

Genetic perturbation of PU.1 binding and chromatin looping at neutrophil enhancers associates with autoimmune disease.

Neutrophils play fundamental roles in innate immune response, shape adaptive immunity, and are a potentially causal cell type underpinning genetic associations with immune system traits and diseases. Here, we profile the binding of myeloid master regulator PU.1 in primary neutrophils across nearly a hundred volunteers. We show that variants associated with differential PU.1 binding underlie genetically-driven differences in cell count and susceptibility to autoimmune and inflammatory diseases. We integrate these results with other multi-individual genomic readouts, revealing coordinated effects of PU.1 binding variants on the local chromatin state, enhancer-promoter contacts and downstream gene expression, and providing a functional interpretation for 27 genes underlying immune traits. Collectively, these results demonstrate the functional role of PU.1 and its target enhancers in neutrophil transcriptional control and immune disease susceptibility.

TOP1 inhibition therapy protects against SARS-CoV-2-induced lethal inflammation.

The ongoing pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is currently affecting millions of lives worldwide. Large retrospective studies indicate that an elevated level of inflammatory cytokines and pro-inflammatory factors are associated with both increased disease severity and mortality. Here, using multidimensional epigenetic, transcriptional, in vitro, and in vivo analyses, we report that topoisomerase 1 (TOP1) inhibition suppresses lethal inflammation induced by SARS-CoV-2. Therapeutic treatment with two doses of topotecan (TPT), an FDA-approved TOP1 inhibitor, suppresses infection-induced inflammation in hamsters. TPT treatment as late as 4 days post-infection reduces morbidity and rescues mortality in a transgenic mouse model. These results support the potential of TOP1 inhibition as an effective host-directed therapy against severe SARS-CoV-2 infection. TPT and its derivatives are inexpensive clinical-grade inhibitors available in most countries. Clinical trials are needed to evaluate the efficacy of repurposing TOP1 inhibitors for severe coronavirus disease 2019 (COVID-19) in humans.

Topoisomerase 1 inhibition therapy protects against SARS-CoV-2-induced inflammation and death in animal models.

The ongoing pandemic caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is currently affecting millions of lives worldwide. Large retrospective studies indicate that an elevated level of inflammatory cytokines and pro-inflammatory factors are associated with both increased disease severity and mortality. Here, using multidimensional epigenetic, transcriptional, in vitro and in vivo analyses, we report that Topoisomerase 1 (Top1) inhibition suppresses lethal inflammation induced by SARS-CoV-2. Therapeutic treatment with two doses of Topotecan (TPT), a FDA-approved Top1 inhibitor, suppresses infection-induced inflammation in hamsters. TPT treatment as late as four days post-infection reduces morbidity and rescues mortality in a transgenic mouse model. These results support the potential of Top1 inhibition as an effective host-directed therapy against severe SARS-CoV-2 infection. TPT and its derivatives are inexpensive clinical-grade inhibitors available in most countries. Clinical trials are needed to evaluate the efficacy of repurposing Top1 inhibitors for COVID-19 in humans.

Cohesin-Dependent and -Independent Mechanisms Mediate Chromosomal Contacts between Promoters and Enhancers.

It is currently assumed that 3D chromosomal organization plays a central role in transcriptional control. However, depletion of cohesin and CTCF affects the steady-state levels of only a minority of transcripts. Here, we use high-resolution Capture Hi-C to interrogate the dynamics of chromosomal contacts of all annotated human gene promoters upon degradation of cohesin and CTCF. We show that a majority of promoter-anchored contacts are lost in these conditions, but many contacts with distinct properties are maintained, and some new ones are gained. The rewiring of contacts between promoters and active enhancers upon cohesin degradation associates with rapid changes in target gene transcription as detected by SLAM sequencing (SLAM-seq). These results provide a mechanistic explanation for the limited, but consistent, effects of cohesin and CTCF depletion on steady-state transcription and suggest the existence of both cohesin-dependent and -independent mechanisms of enhancer-promoter pairing.

Functional effects of variation in transcription factor binding highlight long-range gene regulation by epromoters.

Identifying DNA cis-regulatory modules (CRMs) that control the expression of specific genes is crucial for deciphering the logic of transcriptional control. Natural genetic variation can point to the possible gene regulatory function of specific sequences through their allelic associations with gene expression. However, comprehensive identification of causal regulatory sequences in brute-force association testing without incorporating prior knowledge is challenging due to limited statistical power and effects of linkage disequilibrium. Sequence variants affecting transcription factor (TF) binding at CRMs have a strong potential to influence gene regulatory function, which provides a motivation for prioritizing such variants in association testing. Here, we generate an atlas of CRMs showing predicted allelic variation in TF binding affinity in human lymphoblastoid cell lines and test their association with the expression of their putative target genes inferred from Promoter Capture Hi-C and immediate linear proximity. We reveal >1300 CRM TF-binding variants associated with target gene expression, the majority of them undetected with standard association testing. A large proportion of CRMs showing associations with the expression of genes they contact in 3D localize to the promoter regions of other genes, supporting the notion of 'epromoters': dual-action CRMs with promoter and distal enhancer activity.

Chicdiff: a computational pipeline for detecting differential chromosomal interactions in Capture Hi-C data.

SUMMARY: Capture Hi-C is a powerful approach for detecting chromosomal interactions involving, at least on one end, DNA regions of interest, such as gene promoters. We present Chicdiff, an R package for robust detection of differential interactions in Capture Hi-C data. Chicdiff enhances a state-of-the-art differential testing approach for count data with bespoke normalization and multiple testing procedures that account for specific statistical properties of Capture Hi-C. We validate Chicdiff on published Promoter Capture Hi-C data in human Monocytes and CD4+ T cells, identifying multitudes of cell type-specific interactions, and confirming the overall positive association between promoter interactions and gene expression. AVAILABILITY AND IMPLEMENTATION: Chicdiff is implemented as an R package that is publicly available at https://github.com/RegulatoryGenomicsGroup/chicdiff. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Publisher Correction: Disease-relevant transcriptional signatures identified in individual smooth muscle cells from healthy mouse vessels.

The original version of this Article contained errors in the author affiliations.Martin R. Bennett was incorrectly associated with Nuclear Dynamics Programme, Babraham Institute, Babraham Research Campus, Cambridge, CB22 3AT, UK. This has now been corrected in both the PDF and HTML versions of the Article. Furthermore, Phoebe Oldach was incorrectly associated with Centre for Molecular Informatics, Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK.This has now been corrected in the HTML version of the Article. The PDF version of the Article was correct at the time of publication.

Disease-relevant transcriptional signatures identified in individual smooth muscle cells from healthy mouse vessels.

Vascular smooth muscle cells (VSMCs) show pronounced heterogeneity across and within vascular beds, with direct implications for their function in injury response and atherosclerosis. Here we combine single-cell transcriptomics with lineage tracing to examine VSMC heterogeneity in healthy mouse vessels. The transcriptional profiles of single VSMCs consistently reflect their region-specific developmental history and show heterogeneous expression of vascular disease-associated genes involved in inflammation, adhesion and migration. We detect a rare population of VSMC-lineage cells that express the multipotent progenitor marker Sca1, progressively downregulate contractile VSMC genes and upregulate genes associated with VSMC response to inflammation and growth factors. We find that Sca1 upregulation is a hallmark of VSMCs undergoing phenotypic switching in vitro and in vivo, and reveal an equivalent population of Sca1-positive VSMC-lineage cells in atherosclerotic plaques. Together, our analyses identify disease-relevant transcriptional signatures in VSMC-lineage cells in healthy blood vessels, with implications for disease susceptibility, diagnosis and prevention.

Author Correction: Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks.

In the original version of the Article, the gene symbol for tissue factor pathway inhibitor was inadvertently given as 'TFP1' instead of 'TFPI'. This has now been corrected in both the PDF and HTML versions of the Article.

Genome organization and chromatin analysis identify transcriptional downregulation of insulin-like growth factor signaling as a hallmark of aging in developing B cells.

BACKGROUND: Aging is characterized by loss of function of the adaptive immune system, but the underlying causes are poorly understood. To assess the molecular effects of aging on B cell development, we profiled gene expression and chromatin features genome-wide, including histone modifications and chromosome conformation, in bone marrow pro-B and pre-B cells from young and aged mice. RESULTS: Our analysis reveals that the expression levels of most genes are generally preserved in B cell precursors isolated from aged compared with young mice. Nonetheless, age-specific expression changes are observed at numerous genes, including microRNA encoding genes. Importantly, these changes are underpinned by multi-layered alterations in chromatin structure, including chromatin accessibility, histone modifications, long-range promoter interactions, and nuclear compartmentalization. Previous work has shown that differentiation is linked to changes in promoter-regulatory element interactions. We find that aging in B cell precursors is accompanied by rewiring of such interactions. We identify transcriptional downregulation of components of the insulin-like growth factor signaling pathway, in particular downregulation of Irs1 and upregulation of Let-7 microRNA expression, as a signature of the aged phenotype. These changes in expression are associated with specific alterations in H3K27me3 occupancy, suggesting that Polycomb-mediated repression plays a role in precursor B cell aging. CONCLUSIONS: Changes in chromatin and 3D genome organization play an important role in shaping the altered gene expression profile of aged precursor B cells. Components of the insulin-like growth factor signaling pathways are key targets of epigenetic regulation in aging in bone marrow B cell precursors.

Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks.

Long-range chromosomal interactions bring distal regulatory elements and promoters together to regulate gene expression in biological processes. By performing promoter capture Hi-C (PCHi-C) on human embryonic stem cell-derived cardiomyocytes (hESC-CMs), we show that such promoter interactions are a key mechanism by which enhancers contact their target genes after hESC-CM differentiation from hESCs. We also show that the promoter interactome of hESC-CMs is associated with expression quantitative trait loci (eQTLs) in cardiac left ventricular tissue; captures the dynamic process of genome reorganisation after hESC-CM differentiation; overlaps genome-wide association study (GWAS) regions associated with heart rate; and identifies new candidate genes in such regions. These findings indicate that regulatory elements in hESC-CMs identified by our approach control gene expression involved in ventricular conduction and rhythm of the heart. The study of promoter interactions in other hESC-derived cell types may be of utility in functional investigation of GWAS-associated regions.

Chromosome contacts in activated T cells identify autoimmune disease candidate genes.

BACKGROUND: Autoimmune disease-associated variants are preferentially found in regulatory regions in immune cells, particularly CD4+ T cells. Linking such regulatory regions to gene promoters in disease-relevant cell contexts facilitates identification of candidate disease genes. RESULTS: Within 4 h, activation of CD4+ T cells invokes changes in histone modifications and enhancer RNA transcription that correspond to altered expression of the interacting genes identified by promoter capture Hi-C. By integrating promoter capture Hi-C data with genetic associations for five autoimmune diseases, we prioritised 245 candidate genes with a median distance from peak signal to prioritised gene of 153 kb. Just under half (108/245) prioritised genes related to activation-sensitive interactions. This included IL2RA, where allele-specific expression analyses were consistent with its interaction-mediated regulation, illustrating the utility of the approach. CONCLUSIONS: Our systematic experimental framework offers an alternative approach to candidate causal gene identification for variants with cell state-specific functional effects, with achievable sample sizes.

Lineage-specific dynamic and pre-established enhancer-promoter contacts cooperate in terminal differentiation.

Chromosome conformation is an important feature of metazoan gene regulation; however, enhancer-promoter contact remodeling during cellular differentiation remains poorly understood. To address this, genome-wide promoter capture Hi-C (CHi-C) was performed during epidermal differentiation. Two classes of enhancer-promoter contacts associated with differentiation-induced genes were identified. The first class ('gained') increased in contact strength during differentiation in concert with enhancer acquisition of the H3K27ac activation mark. The second class ('stable') were pre-established in undifferentiated cells, with enhancers constitutively marked by H3K27ac. The stable class was associated with the canonical conformation regulator cohesin, whereas the gained class was not, implying distinct mechanisms of contact formation and regulation. Analysis of stable enhancers identified a new, essential role for a constitutively expressed, lineage-restricted ETS-family transcription factor, EHF, in epidermal differentiation. Furthermore, neither class of contacts was observed in pluripotent cells, suggesting that lineage-specific chromatin structure is established in tissue progenitor cells and is further remodeled in terminal differentiation.