Single-Cell Gene Expression Analysis Revealed Immune Cell Signatures of Delta COVID-19.

The coronavirus disease 2019 (COVID-19) is accompanied by a cytokine storm with the release of many proinflammatory factors and development of respiratory syndrome. Several SARS-CoV-2 lineages have been identified, and the Delta variant (B.1.617), linked with high mortality risk, has become dominant in many countries. Understanding the immune responses associated with COVID-19 lineages may therefore aid the development of therapeutic and diagnostic strategies. Multiple single-cell gene expression studies revealed innate and adaptive immunological factors and pathways correlated with COVID-19 severity. Additional investigations covering host-pathogen response characteristics for infection caused by different lineages are required. Here, we performed single-cell transcriptome profiling of blood mononuclear cells from the individuals with different severity of the COVID-19 and virus lineages to uncover variant specific molecular factors associated with immunity. We identified significant changes in lymphoid and myeloid cells. Our study highlights that an abundant population of monocytes with specific gene expression signatures accompanies Delta lineage of SARS-CoV-2 and contributes to COVID-19 pathogenesis inferring immune components for targeted therapy.

Coupled scRNA-Seq and Intracellular Protein Activity Reveal an Immunosuppressive Role of TREM2 in Cancer.

Cell function and activity are regulated through integration of signaling, epigenetic, transcriptional, and metabolic pathways. Here, we introduce INs-seq, an integrated technology for massively parallel recording of single-cell RNA sequencing (scRNA-seq) and intracellular protein activity. We demonstrate the broad utility of INs-seq for discovering new immune subsets by profiling different intracellular signatures of immune signaling, transcription factor combinations, and metabolic activity. Comprehensive mapping of Arginase 1-expressing cells within tumor models, a metabolic immune signature of suppressive activity, discovers novel Arg1+ Trem2+ regulatory myeloid (Mreg) cells and identifies markers, metabolic activity, and pathways associated with these cells. Genetic ablation of Trem2 in mice inhibits accumulation of intra-tumoral Mreg cells, leading to a marked decrease in dysfunctional CD8+ T cells and reduced tumor growth. This study establishes INs-seq as a broadly applicable technology for elucidating integrated transcriptional and intra-cellular maps and identifies the molecular signature of myeloid suppressive cells in tumors.

ACE-inhibition induces a cardioprotective transcriptional response in the metabolic syndrome heart.

Cardiovascular disease associated with metabolic syndrome has a high prevalence, but the mechanistic basis of metabolic cardiomyopathy remains poorly understood. We characterised the cardiac transcriptome in a murine metabolic syndrome (MetS) model (LDLR-/-; ob/ob, DKO) relative to the healthy, control heart (C57BL/6, WT) and the transcriptional changes induced by ACE-inhibition in those hearts. RNA-Seq, differential gene expression and transcription factor analysis identified 288 genes differentially expressed between DKO and WT hearts implicating 72 pathways. Hallmarks of metabolic cardiomyopathy were increased activity in integrin-linked kinase signalling, Rho signalling, dendritic cell maturation, production of nitric oxide and reactive oxygen species in macrophages, atherosclerosis, LXR-RXR signalling, cardiac hypertrophy, and acute phase response pathways. ACE-inhibition had a limited effect on gene expression in WT (55 genes, 23 pathways), and a prominent effect in DKO hearts (1143 genes, 104 pathways). In DKO hearts, ACE-I appears to counteract some of the MetS-specific pathways, while also activating cardioprotective mechanisms. We conclude that MetS and control murine hearts have unique transcriptional profiles and exhibit a partially specific transcriptional response to ACE-inhibition.

Identification of cis-regulatory mutations generating de novo edges in personalized cancer gene regulatory networks.

The identification of functional non-coding mutations is a key challenge in the field of genomics. Here we introduce µ-cisTarget to filter, annotate and prioritize cis-regulatory mutations based on their putative effect on the underlying "personal" gene regulatory network. We validated µ-cisTarget by re-analyzing the TAL1 and LMO1 enhancer mutations in T-ALL, and the TERT promoter mutation in melanoma. Next, we re-sequenced the full genomes of ten cancer cell lines and used matched transcriptome data and motif discovery to identify master regulators with de novo binding sites that result in the up-regulation of nearby oncogenic drivers. µ-cisTarget is available from http://mucistarget.aertslab.org .

Multiplex enhancer-reporter assays uncover unsophisticated TP53 enhancer logic.

Transcription factors regulate their target genes by binding to regulatory regions in the genome. Although the binding preferences of TP53 are known, it remains unclear what distinguishes functional enhancers from nonfunctional binding. In addition, the genome is scattered with recognition sequences that remain unoccupied. Using two complementary techniques of multiplex enhancer-reporter assays, we discovered that functional enhancers could be discriminated from nonfunctional binding events by the occurrence of a single TP53 canonical motif. By combining machine learning with a meta-analysis of TP53 ChIP-seq data sets, we identified a core set of more than 1000 responsive enhancers in the human genome. This TP53 cistrome is invariably used between cell types and experimental conditions, whereas differences among experiments can be attributed to indirect nonfunctional binding events. Our data suggest that TP53 enhancers represent a class of unsophisticated cell-autonomous enhancers containing a single TP53 binding site, distinct from complex developmental enhancers that integrate signals from multiple transcription factors.

Identification of High-Impact cis-Regulatory Mutations Using Transcription Factor Specific Random Forest Models.

Cancer genomes contain vast amounts of somatic mutations, many of which are passenger mutations not involved in oncogenesis. Whereas driver mutations in protein-coding genes can be distinguished from passenger mutations based on their recurrence, non-coding mutations are usually not recurrent at the same position. Therefore, it is still unclear how to identify cis-regulatory driver mutations, particularly when chromatin data from the same patient is not available, thus relying only on sequence and expression information. Here we use machine-learning methods to predict functional regulatory regions using sequence information alone, and compare the predicted activity of the mutated region with the reference sequence. This way we define the Predicted Regulatory Impact of a Mutation in an Enhancer (PRIME). We find that the recently identified driver mutation in the TAL1 enhancer has a high PRIME score, representing a "gain-of-target" for MYB, whereas the highly recurrent TERT promoter mutation has a surprisingly low PRIME score. We trained Random Forest models for 45 cancer-related transcription factors, and used these to score variations in the HeLa genome and somatic mutations across more than five hundred cancer genomes. Each model predicts only a small fraction of non-coding mutations with a potential impact on the function of the encompassing regulatory region. Nevertheless, as these few candidate driver mutations are often linked to gains in chromatin activity and gene expression, they may contribute to the oncogenic program by altering the expression levels of specific oncogenes and tumor suppressor genes.

Decoding the regulatory landscape of melanoma reveals TEADS as regulators of the invasive cell state.

Transcriptional reprogramming of proliferative melanoma cells into a phenotypically distinct invasive cell subpopulation is a critical event at the origin of metastatic spreading. Here we generate transcriptome, open chromatin and histone modification maps of melanoma cultures; and integrate this data with existing transcriptome and DNA methylation profiles from tumour biopsies to gain insight into the mechanisms underlying this key reprogramming event. This shows thousands of genomic regulatory regions underlying the proliferative and invasive states, identifying SOX10/MITF and AP-1/TEAD as regulators, respectively. Knockdown of TEADs shows a previously unrecognized role in the invasive gene network and establishes a causative link between these transcription factors, cell invasion and sensitivity to MAPK inhibitors. Using regulatory landscapes and in silico analysis, we show that transcriptional reprogramming underlies the distinct cellular states present in melanoma. Furthermore, it reveals an essential role for the TEADs, linking it to clinically relevant mechanisms such as invasion and resistance.

iRegulon: from a gene list to a gene regulatory network using large motif and track collections.

Identifying master regulators of biological processes and mapping their downstream gene networks are key challenges in systems biology. We developed a computational method, called iRegulon, to reverse-engineer the transcriptional regulatory network underlying a co-expressed gene set using cis-regulatory sequence analysis. iRegulon implements a genome-wide ranking-and-recovery approach to detect enriched transcription factor motifs and their optimal sets of direct targets. We increase the accuracy of network inference by using very large motif collections of up to ten thousand position weight matrices collected from various species, and linking these to candidate human TFs via a motif2TF procedure. We validate iRegulon on gene sets derived from ENCODE ChIP-seq data with increasing levels of noise, and we compare iRegulon with existing motif discovery methods. Next, we use iRegulon on more challenging types of gene lists, including microRNA target sets, protein-protein interaction networks, and genetic perturbation data. In particular, we over-activate p53 in breast cancer cells, followed by RNA-seq and ChIP-seq, and could identify an extensive up-regulated network controlled directly by p53. Similarly we map a repressive network with no indication of direct p53 regulation but rather an indirect effect via E2F and NFY. Finally, we generalize our computational framework to include regulatory tracks such as ChIP-seq data and show how motif and track discovery can be combined to map functional regulatory interactions among co-expressed genes. iRegulon is available as a Cytoscape plugin from http://iregulon.aertslab.org.