Accurate differential analysis of transcription factor activity from gene expression

Characterisation of the Circulating Transcriptomic Landscape in Inflammatory Bowel Disease Provides Evidence for Dysregulation of Multiple Transcription Factors Including NFE2, SPI1, CEBPB, and IRF2

AIM

To assess the pathobiological and translational importance of whole-blood transcriptomic analysis in inflammatory bowel disease [IBD].

METHODS

We analysed whole-blood expression profiles from paired-end sequencing in a discovery cohort of 590 Europeans recruited across six countries in the IBD Character initiative (newly diagnosed patients with Crohn's disease [CD; n = 156], ulcerative colitis [UC; n = 167], and controls [n = 267]), exploring differential expression [DESeq2], co-expression networks [WGCNA], and transcription factor involvement [EPEE, ChEA, DoRothEA]. Findings were validated by analysis of an independent replication cohort [99 CD, 100 UC, 95 controls]. In the discovery cohort, we also defined baseline expression correlates of future treatment escalation using cross-validated elastic-net and random forest modelling, along with a pragmatic ratio detection procedure.

RESULTS

Disease-specific transcriptomes were defined in IBD [8697 transcripts], CD [7152], and UC [8521], with the most highly significant changes in single genes, including CD177 (log2-fold change [LFC] = 4.63, p = 4.05 × 10-118), MCEMP1 [LFC = 2.45, p = 7.37 × 10-109], and S100A12 [LFC = 2.31, p = 2.15 × 10-93]. Significantly over-represented pathways included IL-1 [p = 1.58 × 10-11], IL-4, and IL-13 [p = 8.96 × 10-9]. Highly concordant results were obtained using multiple regulatory activity inference tools applied to the discovery and replication cohorts. These analyses demonstrated central roles in IBD for the transcription factors NFE2, SPI1 [PU.1], CEBPB, and IRF2, all regulators of cytokine signalling, based on a consistent signal across cohorts and transcription factor ranking methods. A number of simple transcriptome-based models were associated with the need for treatment escalation, including the binary CLEC5A/CDH2 expression ratio in UC (hazard ratio = 23.4, 95% confidence interval [CI] 5.3-102.0).

CONCLUSIONS

Transcriptomic analysis has allowed for a detailed characterisation of IBD pathobiology, with important potential translational implications.

Atg7 deficiency in microglia drives an altered transcriptomic profile associated with an impaired neuroinflammatory response

Microglia, resident immunocompetent cells of the central nervous system, can display a range of reaction states and thereby exhibit distinct biological functions across development, adulthood and under disease conditions. Distinct gene expression profiles are reported to define each of these microglial reaction states. Hence, the identification of modulators of selective microglial transcriptomic signature, which have the potential to regulate unique microglial function has gained interest. Here, we report the identification of ATG7 (Autophagy-related 7) as a selective modulator of an NF-κB-dependent transcriptional program controlling the pro-inflammatory response of microglia. We also uncover that microglial Atg7-deficiency was associated with reduced microglia-mediated neurotoxicity, and thus a loss of biological function associated with the pro-inflammatory microglial reactive state. Further, we show that Atg7-deficiency in microglia did not impact on their ability to respond to alternative stimulus, such as one driving them towards an anti-inflammatory/tumor supportive phenotype. The identification of distinct regulators, such as Atg7, controlling specific microglial transcriptional programs could lead to developing novel therapeutic strategies aiming to manipulate selected microglial phenotypes, instead of the whole microglial population with is associated with several pitfalls.

Determinants of correlated expression of transcription factors and their target genes

While transcription factors (TFs) are known to regulate the expression of their target genes (TGs), only a weak correlation of expression between TFs and their TGs has generally been observed. As lack of correlation could be caused by additional layers of regulation, the overall correlation distribution may hide the presence of a subset of regulatory TF-TG pairs with tight expression coupling. Using reported regulatory pairs in the plant Arabidopsis thaliana along with comprehensive gene expression information and testing a wide array of molecular features, we aimed to discern the molecular determinants of high expression correlation of TFs and their TGs. TF-family assignment, stress-response process involvement, short genomic distances of the TF-binding sites to the transcription start site of their TGs, few required protein-protein-interaction connections to establish physical interactions between the TF and polymerase-II, unambiguous TF-binding motifs, increased numbers of miRNA target-sites in TF-mRNAs, and a young evolutionary age of TGs were found particularly indicative of high TF-TG correlation. The modulating roles of post-transcriptional, post-translational processes, and epigenetic factors have been characterized as well. Our study reveals that regulatory pairs with high expression coupling are associated with specific molecular determinants.