Search

Hypoxia-inducible transcription factor 1 (HIF-1) has been shown to enhance microbial killing and ameliorate the course of bacterial infections. While the impact of HIF-1 on inflammatory diseases of the gut has been studied intensively, its function in bacterial infections of the gastrointestinal tract remains largely elusive. With the help of a publicly available gene expression data set, we inferred significant activation of HIF-1 after oral infection of mice with Salmonella enterica serovar Typhimurium. Immunohistochemistry and Western blot analyses confirmed marked HIF-1α protein stabilization, especially in the intestinal epithelium. This prompted us to analyze conditional Hif1a-deficient mice to examine cell type-specific functions of HIF-1 in this model. Our results demonstrate enhanced noncanonical induction of HIF-1 activity upon Salmonella infection in the intestinal epithelium as well as in macrophages. Surprisingly, Hif1a deletion in intestinal epithelial cells did not impact inflammatory gene expression, bacterial spread, or disease outcomes. In contrast, Hif1a deletion in myeloid cells enhanced intestinal Cxcl2 expression and reduced the cecal Salmonella load. In vitro, HIF-1α-deficient macrophages showed overall impaired transcription of mRNA encoding proinflammatory factors; however, the intracellular survival of Salmonella was not impacted by HIF-1α deficiency.

Nucleic acids symposium series 47(19), 10010-10026 (2019). doi:10.1093/nar/gkz805, Published by Oxford Univ. Press22313, Oxford

Background Transcriptomic studies have contributed to fundamental knowledge of myocardial remodeling in human heart failure (HF). However, the key HF genes reported are often inconsistent between studies, and systematic efforts to integrate evidence from multiple patient cohorts are lacking. Here, we aimed to provide a framework for comprehensive comparison and analysis of publicly available data sets resulting in an unbiased consensus transcriptional signature of human end‐stage HF. Methods and Results We curated and uniformly processed 16 public transcriptomic studies of left ventricular samples from 263 healthy and 653 failing human hearts. First, we evaluated the degree of consistency between studies by using linear classifiers and overrepresentation analysis. Then, we meta‐analyzed the deregulation of 14 041 genes to extract a consensus signature of HF. Finally, to functionally characterize this signature, we estimated the activities of 343 transcription factors, 14 signaling pathways, and 182 micro RNAs, as well as the enrichment of 5998 biological processes. Machine learning approaches revealed conserved disease patterns across all studies independent of technical differences. These consistent molecular changes were prioritized with a meta‐analysis, functionally characterized and validated on external data. We provide all results in a free public resource ( https://saezlab.shinyapps.io/reheat/ ) and exemplified usage by deciphering fetal gene reprogramming and tracing the potential myocardial origin of the plasma proteome markers in patients with HF. Conclusions Even though technical and sampling variability confound the identification of differentially expressed genes in individual studies, we demonstrated that coordinated molecular responses during end‐stage HF are conserved. The presented resource is crucial to complement findings in independent studies and decipher fundamental changes in failing myocardium.

The hypoxia-inducible transcription factor 1 (HIF-1) has been shown to ameliorate different bacterial infections through enhancement of microbial killing. While the impact of HIF-1 on inflammatory diseases of the gut has been studied intensively, its function in bacterial infections of the intestine remains largely elusive. With the help of a publicly available gene expression data set, we could infer significant activation of the HIF-1 transcription factor after oral infection of mice with Salmonella Typhimurium. This prompted us to apply lineage-restricted deletion of the Hif1a locus in mice to examine cell type-specific functions of HIF-1 in this model. We show hypoxia-independent induction of HIF-1 activity upon Salmonella infection in the intestinal epithelium as well as in macrophages. Surprisingly, Hif1a deletion in intestinal epithelial cells impacted neither disease outcome nor inflammatory activity. The conditional knockout of Hif1a in myeloid cells enhanced the mRNA expression of the largely pro-inflammatory chemokine Cxcl2, revealing a potentially inflammatory effect of HIF-1 deficiency in myeloid cells in the gut in vivo. Again, the disease outcome was not affected. In vitro HIF-1-deficient macrophages showed an overall impaired transcription of pro-inflammatory factors, however, Salmonella bypassed direct intracellular, bactericidal HIF-1-dependent mechanisms in a Salmonella pathogenicity island (SPI)-2 independent manner. Taken together, our data suggest that HIF-1 in intestinal epithelial and myeloid cells is either dispensable or compensable in the immune defense against Salmonella Typhimurium.

Mouse models are frequently used to study chronic liver diseases (CLDs). To assess their translational relevance, we quantified the similarity of commonly used mouse models to human CLDs based on transcriptome data. Gene‐expression data from 372 patients were compared with data from acute and chronic mouse models consisting of 227 mice, and additionally to nine published gene sets of chronic mouse models. Genes consistently altered in humans and mice were mapped to liver cell types based on single‐cell RNA‐sequencing data and validated by immunostaining. Considering the top differentially expressed genes, the similarity between humans and mice varied among the mouse models and depended on the period of damage induction. The highest recall (0.4) and precision (0.33) were observed for the model with 12‐months damage induction by CCl4 and by a Western diet, respectively. Genes consistently up‐regulated between the chronic CCl4 model and human CLDs were enriched in inflammatory and developmental processes, and mostly mapped to cholangiocytes, macrophages, and endothelial and mesenchymal cells. Down‐regulated genes were enriched in metabolic processes and mapped to hepatocytes. Immunostaining confirmed the regulation of selected genes and their cell type specificity. Genes that were up‐regulated in both acute and chronic models showed higher recall and precision with respect to human CLDs than exclusively acute or chronic genes. Conclusion: Similarly regulated genes in human and mouse CLDs were identified. Despite major interspecies differences, mouse models detected 40% of the genes significantly altered in human CLD. The translational relevance of individual genes can be assessed at https://saezlab.shinyapps.io/liverdiseaseatlas/., In the present study, we observed that – although major interspecies differences remain – improved mouse models show up to 40% of the gene expression changes seen in liver tissue of human NAFLD or NASH, which is much higher than previously reported. Moreover, we identified gene sets consistently regulated in human and mouse chronic liver disease. Based on single‐cell RNA‐sequencing data set we mapped liver cell types to those genes and validated them by immunostaining.

Background & Aims In chronic liver diseases, inflammation induces oxidative stress and thus may contribute to the progression of liver injury, fibrosis, and carcinogenesis. The KEAP1/NRF2 axis is a major regulator of cellular redox balance. In the present study, we investigated whether the KEAP1/NRF2 system is involved in liver disease progression in humans and mice. Methods The clinical relevance of oxidative stress was investigated by liver RNA sequencing in a well-characterized cohort of patients with non-alcoholic fatty liver disease (n = 63) and correlated with histological and clinical parameters. For functional analysis, hepatocyte-specific Nemo knockout (NEMOΔhepa) mice were crossed with hepatocyte-specific Keap1 knockout (KEAP1Δhepa) mice. Results Immunohistochemical analysis of human liver sections showed increased oxidative stress and high NRF2 expression in patients with chronic liver disease. RNA sequencing of liver samples in a human pediatric NAFLD cohort revealed a significant increase of NRF2 activation correlating with the grade of inflammation, but not with the grade of steatosis, which could be confirmed in a second adult NASH cohort. In mice, microarray analysis revealed that Keap1 deletion induces NRF2 target genes involved in glutathione metabolism and xenobiotic stress (e.g., Nqo1). Furthermore, deficiency of one of the most important antioxidants, glutathione (GSH), in NEMOΔhepa livers was rescued after deleting Keap1. As a consequence, NEMOΔhepa/KEAP1Δhepa livers showed reduced apoptosis compared to NEMOΔhepa livers as well as a dramatic downregulation of genes involved in cell cycle regulation and DNA replication. Consequently, NEMOΔhepa/KEAP1Δhepa compared to NEMOΔhepa livers displayed decreased fibrogenesis, lower tumor incidence, reduced tumor number, and decreased tumor size. Conclusions NRF2 activation in patients with non-alcoholic steatohepatitis correlates with the grade of inflammation, but not steatosis. Functional analysis in mice demonstrated that NRF2 activation in chronic liver disease is protective by ameliorating fibrogenesis, initiation and progression of hepatocellular carcinogenesis. Lay summary The KEAP1 (Kelch-like ECH-associated protein-1)/NRF2 (erythroid 2-related factor 2) axis has a major role in regulating cellular redox balance. Herein, we show that NRF2 activity correlates with the grade of inflammation in patients with non-alcoholic steatohepatitis. Functional studies in mice actually show that NRF2 activation, resulting from KEAP1 deletion, protects against fibrosis and cancer.

2.AbstractAimsTranscriptomic studies have contributed to fundamental knowledge of myocardial remodeling in human heart failure (HF). However, the agreement on key regulated genes in HF is limited and systematic efforts to integrate evidence from multiple patient cohorts are lacking. Here we aimed to provide an unbiased consensus transcriptional signature of human end-stage HF by comprehensive comparison and analysis of publicly available datasets.Methods and ResultsWe curated and uniformly processed 16 public transcriptomic studies of left ventricular samples from 263 healthy and 653 failing human hearts. Transfer learning approaches revealed conserved disease patterns across all studies independent of technical differences. We meta-analyzed the dysregulation of 14041 genes to extract a consensus signature of HF. Estimation of the activities of 343 transcription factors, 14 signalling pathways, and 182 micro RNAs, as well as the enrichment of 5998 biological processes confirmed the established aspects of the functional landscape of the disease and revealed novel ones. We provide all results in a free public resource https://saezlab.shinyapps.io/reheat/ to facilitate further use and interpretation of the results. We exemplify usage by deciphering fetal gene reprogramming and tracing myocardial origin of the plasma proteome biomarkers in HF patients.ConclusionWe demonstrated the feasibility of combining transcriptional studies from different HF patient cohorts. This compendium provides a robust and consistent collection of molecular markers of end-stage HF that may guide the identification of novel targets with diagnostic or therapeutic relevance.

BACKGROUND: Many functional analysis tools have been developed to extract functional and mechanistic insight from bulk transcriptome data. With the advent of single-cell RNA sequencing (scRNA-seq), it is in principle possible to do such an analysis for single cells. However, scRNA-seq data has characteristics such as drop-out events and low library sizes. It is thus not clear if functional TF and pathway analysis tools established for bulk sequencing can be applied to scRNA-seq in a meaningful way. RESULTS: To address this question, we perform benchmark studies on simulated and real scRNA-seq data. We include the bulk-RNA tools PROGENy, GO enrichment, and DoRothEA that estimate pathway and transcription factor (TF) activities, respectively, and compare them against the tools SCENIC/AUCell and metaVIPER, designed for scRNA-seq. For the in silico study, we simulate single cells from TF/pathway perturbation bulk RNA-seq experiments. We complement the simulated data with real scRNA-seq data upon CRISPR-mediated knock-out. Our benchmarks on simulated and real data reveal comparable performance to the original bulk data. Additionally, we show that the TF and pathway activities preserve cell type-specific variability by analyzing a mixture sample sequenced with 13 scRNA-seq protocols. We also provide the benchmark data for further use by the community. CONCLUSIONS: Our analyses suggest that bulk-based functional analysis tools that use manually curated footprint gene sets can be applied to scRNA-seq data, partially outperforming dedicated single-cell tools. Furthermore, we find that the performance of functional analysis tools is more sensitive to the gene sets than to the statistic used. CHH is supported by the German Federal Ministry of Education and Research (BMBF)-funded project Systems Medicine of the Liver (LiSyM, FKZ: 031 L0049). MPK, BAJ, and DAL are supported by NIH Grant U54-CA217377. BS is supported by the Premium Postdoctoral Fellowship Program of the Hungarian Academy of Sciences. HH is a Miguel Servet (CP14/00229) researcher funded by the Spanish Institute of Health Carlos III (ISCIII). This work has received funding from the Ministerio de Ciencia, Innovación y Universidades (SAF2017-89109-P; AEI/FEDER, UE)

Introduction To develop effective therapies and identify novel early biomarkers for chronic kidney disease, an understanding of the molecular mechanisms orchestrating it is essential. We here set out to understand how differences in chronic kidney disease (CKD) origin are reflected in gene expression. To this end, we integrated publicly available human glomerular microarray gene expression data for 9 kidney disease entities that account for most of CKD worldwide. Our primary goal was to demonstrate the possibilities and potential on data analysis and integration to the nephrology community. Methods We integrated data from 5 publicly available studies and compared glomerular gene expression profiles of disease with that of controls from nontumor parts of kidney cancer nephrectomy tissues. A major challenge was the integration of the data from different sources, platforms, and conditions that we mitigated with a bespoke stringent procedure. Results We performed a global transcriptome-based delineation of different kidney disease entities, obtaining a transcriptomic diffusion map of their similarities and differences based on the genes that acquire a consistent differential expression between each kidney disease entity and nephrectomy tissue. We derived functional insights by inferring the activity of signaling pathways and transcription factors from the collected gene expression data and identified potential drug candidates based on expression signature matching. We validated representative findings by immunostaining in human kidney biopsies indicating, for example, that the transcription factor FOXM1 is significantly and specifically expressed in parietal epithelial cells in rapidly progressive glomerulonephritis (RPGN) whereas not expressed in control kidney tissue. Furthermore, we found drug candidates by matching the signature on expression of drugs to that of the CKD entities, in particular, the Food and Drug Administration–approved drug nilotinib. Conclusion These results provide a foundation to comprehend the specific molecular mechanisms underlying different kidney disease entities that can pave the way to identify biomarkers and potential therapeutic targets. To facilitate further use, we provide our results as a free interactive Web application: https://saezlab.shinyapps.io/ckd_landscape/. However, because of the limitations of the data and the difficulties in its integration, any specific result should be considered with caution. Indeed, we consider this study rather an illustration of the value of functional genomics and integration of existing data., Graphical abstract

Mouse models of non-alcoholic fatty liver disease (NAFLD) are required to define therapeutic targets, but detailed time-resolved studies to establish a sequence of events are lacking. Here, we fed male C57Bl/6N mice a Western or standard diet over 48 weeks. Multiscale time-resolved characterization was performed using RNA-seq, histopathology, immunohistochemistry, intravital imaging, and blood chemistry, the results were compared to human disease. Acetaminophen toxicity and ammonia metabolism were additionally analyzed as functional readouts. We identified a sequence of eight key events: formation of lipid droplets, inflammatory foci, lipogranulomas, zonal reorganization, cell death and replacement proliferation, ductular reaction, fibrogenesis, and hepatocellular cancer. Functional changes included resistance to acetaminophen and altered nitrogen metabolism. The transcriptomic landscape was characterized by two large clusters of monotonously increasing or decreasing genes, and a smaller number of ‘rest-and-jump genes’ that initially remained unaltered but became differentially expressed only at week 12 or later. Approximately 30% of the genes altered in human NAFLD are also altered in the present mouse model and an increasing overlap with genes altered in human HCC occurred at weeks 30–48. In conclusion, the observed sequence of events recapitulates many features of human disease and offers a basis for the identification of therapeutic targets.

Genome research 29(8), 1363-1375 (2019). doi:10.1101/gr.240663.118, Published by HighWire Press, Stanford, Calif.

Many tools have been developed to extract functional and mechanistic insight from bulk transcriptome profiling data. With the advent of single-cell RNA sequencing (scRNA-seq), it is in principle possible to do such an analysis for single cells. However, scRNA-seq data has characteristics such as drop-out events, low library sizes and a comparatively large number of samples/cells. It is thus not clear if functional genomics tools established for bulk sequencing can be applied to scRNA-seq in a meaningful way. To address this question, we performed benchmark studies on in silico and in vitro single-cell RNA-seq data. We included the bulk-RNA tools PROGENy, GO enrichment and DoRothEA that estimate pathway and transcription factor (TF) activities, respectively, and compared them against the tools AUCell and metaVIPER, designed for scRNA-seq. For the in silico study we simulated single cells from TF/pathway perturbation bulk RNA-seq experiments. Our simulation strategy guarantees that the information of the original perturbation is preserved while resembling the characteristics of scRNA-seq data. We complemented the in silico data with in vitro scRNA-seq data upon CRISPR-mediated knock-out. Our benchmarks on both the simulated and real data revealed comparable performance to the original bulk data. Additionally, we showed that the TF and pathway activities preserve cell-type specific variability by analysing a mixture sample sequenced with 13 scRNA-seq different protocols. Our analyses suggest that bulk functional genomics tools can be applied to scRNA-seq data, outperforming dedicated single cell tools. Furthermore we provide a benchmark for further methods development by the community.

The prediction of transcription factor (TF) activities from the gene expression of their targets (i.e., TF regulon) is becoming a widely used approach to characterize the functional status of transcriptional regulatory circuits. Several strategies and data sets have been proposed to link the target genes likely regulated by a TF, each one providing a different level of evidence. The most established ones are (1) manually curated repositories, (2) interactions derived from ChIP-seq binding data, (3) in silico prediction of TF binding on gene promoters, and (4) reverse-engineered regulons from large gene expression data sets. However, it is not known how these different sources of regulons affect the TF activity estimations and, thereby, downstream analysis and interpretation. Here we compared the accuracy and biases of these strategies to define human TF regulons by means of their ability to predict changes in TF activities in three reference benchmark data sets. We assembled a collection of TF-target interactions for 1541 human TFs and evaluated how different molecular and regulatory properties of the TFs, such as the DNA-binding domain, specificities, or mode of interaction with the chromatin, affect the predictions of TF activity. We assessed their coverage and found little overlap on the regulons derived from each strategy and better performance by literature-curated information followed by ChIP-seq data. We provide an integrated resource of all TF-target interactions derived through these strategies, with confidence scores, as a resource for enhanced prediction of TF activities.

Transcriptome profiling followed by differential gene expression analysis often leads to lists of genes that are hard to analyze and interpret. Functional genomics tools are powerful approaches for downstream analysis, as they summarize the large and noisy gene expression space into a smaller number of biological meaningful features. In particular, methods that estimate the activity of processes by mapping transcripts level to process members are popular. However, footprints of either a pathway or transcription factor (TF) on gene expression show superior performance over mapping-based gene sets. These footprints are largely developed for humans and their usability in the broadly-used model organism Mus musculus is uncertain. Evolutionary conservation of the gene regulatory system suggests that footprints of human pathways and TFs can functionally characterize mice data. In this paper we analyze this hypothesis. We perform a comprehensive benchmark study exploiting two state-of-the-art footprint methods, DoRothEA and an extended version of PROGENy. These methods infer TF and pathway activity, respectively. Our results show that both can recover mouse perturbations, confirming our hypothesis that footprints are conserved between mice and humans. Subsequently, we illustrate the usability of PROGENy and DoRothEA by recovering pathway/TF-disease associations from newly generated disease sets. Additionally, we provide pathway and TF activity scores for a large collection of human and mouse perturbation and disease experiments (2374). We believe that this resource, available for interactive exploration and download ( https://saezlab.shinyapps.io/footprint_scores/ ), can have broad applications including the study of diseases and therapeutics.

Transcriptome profiling followed by differential gene expression analysis often leads to unclear lists of genes which are hard to analyse and interpret. Functional genomic tools are powerful approaches for downstream analysis, as they summarize the large and noisy gene expression space in a smaller number of biological meaningful features. In particular, methods that estimate the activity of processes by mapping transcripts level to process members are popular. However, footprints of either a pathway or transcription factor (TF) on gene expression show superior performance over mapping-based gene sets. These footprints are largely developed for human and their usability in the broadly-used model organism Mus musculus is uncertain. Evolutionary conservation of the gene regulatory system suggests that footprints of human pathways and TFs can functionally characterize mice data. In this paper we analyze this hypothesis. We perform a comprehensive benchmark study exploiting two state-of-the-art footprint methods, DoRothEA and an extended version of PROGENy. These methods infer TF and pathway activity, respectively. Our results show that both can recover mouse perturbations, confirming our hypothesis that footprints are conserved between mice and humans. Subsequently, we illustrate the usability of PROGENy and DoRothEA by recovering pathway/TF-disease associations from newly generated disease sets. Additionally, we provide pathway and TF activity scores for a large collection of human and mouse perturbation and disease experiments (2,374). We believe that this resource, available for interactive exploration and download (https://saezlab.shinyapps.io/footprint_scores/), can have broad applications including the study of diseases and therapeutics.

To develop efficient therapies and identify novel early biomarkers for chronic kidney disease an understanding of the molecular mechanisms orchestrating it is essential. We here set out to understand how differences in CKD origin are reflected in gene expression. To this end, we integrated publicly available human glomerular microarray gene expression data for nine kidney disease entities that account for a majority of CKD worldwide. We included data from five distinct studies and compared glomerular gene expression profiles to that of non-tumor parts of kidney cancer nephrectomy tissues. A major challenge was the integration of the data from different sources, platforms and conditions, that we mitigated with a bespoke stringent procedure. This allowed us to perform a global transcriptome-based delineation of different kidney disease entities, obtaining a landscape of their similarities and differences based on the genes that acquire a consistent differential expression between each kidney disease entity and nephrectomy tissue. Furthermore, we derived functional insights by inferring activity of signaling pathways and transcription factors from the collected gene expression data, and identified potential drug candidates based on expression signature matching. We validated representative findings by immunostaining in human kidney biopsies indicating e.g. that the transcription factor FOXM1 is significantly and specifically expressed in parietal epithelial cells in RPGN whereas not expressed in control kidney tissue. These results provide a foundation to comprehend the specific molecular mechanisms underlying different kidney disease entities, that can pave the way to identify biomarkers and potential therapeutic targets. To facilitate this, we provide our results as a free interactive web application: https://saezlab.shinyapps.io/ckd_landscape/.Translational StatementChronic kidney disease is a combination of entities with different etiologies. We integrate and analyse transcriptomics analysis of glomerular from different entities to dissect their different pathophysiology, what might help to identify novel entity-specific therapeutic targets.