Computational identification of transcriptional regulators in human endotoxemia

Rapid induction of expression by LPS is accompanied by favorable chromatin and rapid binding of c-Jun

The response to infection is managed in mammals by a coordinated immune response. Innate responses are rapid and hard wired and have been demonstrated to be regulated at the level of chromatin accessibility. This study examined primary human monocyte responses to LPS as a model of innate responses to bacteria. We utilized inhibitors of chromatin modifying enzymes to understand the inter-relationships of the chromatin complexes regulating transcription. Multiplex digital gene detection was utilized to quantitate changes in mRNA levels for genes induced by LPS. In the first 30 min, genes that were highly induced by LPS as a group exhibited minimal effect of the chemical inhibitors of chromatin modifications. At 60 min, the more highly expressed genes were markedly more inhibitable. The effects of the inhibitors were almost entirely concordant in spite of different mechanisms of action. Two focus groups of genes with either high LPS inducibility at 30 min or high LPS inducibility at 60 min (but not at 30 min) were further examined by ChIP assay. NFκB p65 binding was increased at the promoters of 30- and 60-min highly inducible genes equivalently. Binding of c-Jun was increased after LPS in the 30-min inducible gene set but not the 60-min inducible gene set. H3K4me3 and H4ac were not detectably altered by LPS stimulation. Baseline H3K4me3 and H4ac were higher in the 30-min highly inducible gene set compared to the 60-min highly inducible gene set. NFκB and JNK inhibitors led to diminished H4ac after LPS. The effects of DRB and C646 were greater for LPS-induced IL6 transcription at 30 min and LPS-stimulated H4ac compared to TNF where transcription was largely unaffected by the inhibitors. In conclusion, genes with very rapidly induced expression after LPS exhibited more favorable chromatin characteristics at baseline and were less inhibitable than genes induced at the later time points.

On the analysis of complex biological supply chains: From Process Systems Engineering to Quantitative Systems Pharmacology

The use of models in biology has become particularly relevant as it enables investigators to develop a mechanistic framework for understanding the operating principles of living systems as well as in quantitatively predicting their response to both pathological perturbations and pharmacological interventions. This application has resulted in a synergistic convergence of systems biology and pharmacokinetic-pharmacodynamic modeling techniques that has led to the emergence of quantitative systems pharmacology (QSP). In this review, we discuss how the foundational principles of chemical process systems engineering inform the progressive development of more physiologically-based systems biology models.

Understanding Physiology in the Continuum: Integration of Information from Multiple -Omics Levels

In this paper, we discuss approaches for integrating biological information reflecting diverse physiologic levels. In particular, we explore statistical and model-based methods for integrating transcriptomic, proteomic and metabolomics data. Our case studies reflect responses to a systemic inflammatory stimulus and in response to an anti-inflammatory treatment. Our paper serves partly as a review of existing methods and partly as a means to demonstrate, using case studies related to human endotoxemia and response to methylprednisolone (MPL) treatment, how specific questions may require specific methods, thus emphasizing the non-uniqueness of the approaches. Finally, we explore novel ways for integrating -omics information with PKPD models, toward the development of more integrated pharmacology models.

Revealing complex function, process and pathway interactions with high-throughput expression and biological annotation data

The biological relationships both between and within the functions, processes and pathways that operate within complex biological systems are only poorly characterized, making the interpretation of large scale gene expression datasets extremely challenging. Here, we present an approach that integrates gene expression and biological annotation data to identify and describe the interactions between biological functions, processes and pathways that govern a phenotype of interest. The product is a global, interconnected network, not of genes but of functions, processes and pathways, that represents the biological relationships within the system. We validated our approach on two high-throughput expression datasets describing organismal and organ development. Our findings are well supported by the available literature, confirming that developmental processes and apoptosis play key roles in cell differentiation. Furthermore, our results suggest that processes related to pluripotency and lineage commitment, which are known to be critical for development, interact mainly indirectly, through genes implicated in more general biological processes. Moreover, we provide evidence that supports the relevance of cell spatial organization in the developing liver for proper liver function. Our strategy can be viewed as an abstraction that is useful to interpret high-throughput data and devise further experiments.

Quantitative Systems Pharmacology: A Framework for Context

Quantitative Systems Pharmacology (QSP) is receiving increased attention. As the momentum builds and the expectations grow it is important to (re)assess and formalize the basic concepts and approaches. In this short review, I argue that QSP, in addition to enabling the rational integration of data and development of complex models, maybe more importantly, provides the foundations for developing an integrated framework for the assessment of drugs and their impact on disease within a broader context expanding the envelope to account in great detail for physiology, environment and prior history. I articulate some of the critical enablers, major obstacles and exciting opportunities manifesting themselves along the way. Charting such overarching themes will enable practitioners to identify major and defining factors as the field progressively moves towards personalized and precision health care delivery.

Targeting Androgen/Estrogen Receptors Crosstalk in Cancer

The actions of estrogens are mediated by estrogen receptors, ERα and ERβ. Recent genomic landscaping of ERα- and ERβ-binding sites has revealed important distinctions regarding their transcriptional activity. ERβ and its isoforms have been correlated with endocrine treatment responsiveness in breast tumors, while post-translational modifications, receptor dimerization patterns, and subcellular localization are increasingly recognized as crucial modulators in prostate carcinogenesis. Androgen receptor (AR) is essential for the development and progression of prostate cancer as well as of certain breast cancer types. The balance between the activity of these two hormone receptors and their molecular interactions in different clinical settings is influenced by several coregulators. This comprises a dynamic regulatory network enhancing or limiting the activity of AR-directed treatments in breast and prostate tumorigenesis. In this review, we discuss the molecular background regarding the therapeutic targeting of androgen/estrogen receptor crosstalk in breast and prostate cancer.

Integrated transcriptional and metabolic profiling in human endotoxemia

In this meta-study, we aimed to integrate biological insights gained from two levels of -omics analyses on the response to systemic inflammation induced by lipopolysaccharide in humans. We characterized the interplay between plasma metabolite compositions and transcriptional response of leukocytes through integration of transcriptomics with plasma metabonomics. We hypothesized that the drastic changes in the immediate environment of the leukocytes might have an adaptive effect on shaping their transcriptional response in conjunction with the initial inflammatory stimuli. Indeed, we observed that leukocytes, most notably, tune the activity of lipid- and protein-associated processes at the transcriptional level in accordance with the fluctuations in metabolite compositions of surrounding plasma. A closer look into the transcriptional control of only metabolic pathways uncovered alterations in bioenergetics and defenses against oxidative stress closely associated with mitochondrial dysfunction and shifts in energy production observed during inflammatory processes.

Human metabolic response to systemic inflammation: assessment of the concordance between experimental endotoxemia and clinical cases of sepsis/SIRS

Introduction

Two recent, independent, studies conducted novel metabolomics analyses relevant to human sepsis progression; one was a human model of endotoxin (lipopolysaccharide (LPS)) challenge (experimental endotoxemia) and the other was community acquired pneumonia and sepsis outcome diagnostic study (CAPSOD). The purpose of the present study was to assess the concordance of metabolic responses to LPS and community-acquired sepsis.

Methods

We tested the hypothesis that the patterns of metabolic response elicited by endotoxin would agree with those in clinical sepsis. Alterations in the plasma metabolome of the subjects challenged with LPS were compared with those of sepsis patients who had been stratified into two groups: sepsis patients with confirmed infection and non-infected patients who exhibited systemic inflammatory response syndrome (SIRS) criteria. Common metabolites between endotoxemia and both these groups were individually identified, together with their direction of change and functional classifications.

Results

Response to endotoxemia at the metabolome level elicited characteristics that agree well with those observed in sepsis patients despite the high degree of variability in the response of these patients. Moreover, some distinct features of SIRS have been identified. Upon stratification of sepsis patients based on 28-day survival, the direction of change in 21 of 23 metabolites was the same in endotoxemia and sepsis survival groups.

Conclusions

The observed concordance in plasma metabolomes of LPS-treated subjects and sepsis survivors strengthens the relevance of endotoxemia to clinical research as a physiological model of community-acquired sepsis, and gives valuable insights into the metabolic changes that constitute a homeostatic response. Furthermore, recapitulation of metabolic differences between sepsis non-survivors and survivors in LPS-treated subjects can enable further research on the development and assessment of rational clinical therapies to prevent sepsis mortality. Compared with earlier studies which focused exclusively on comparing transcriptional dynamics, the distinct metabolomic responses to systemic inflammation with or without confirmed infection, suggest that the metabolome is much better at differentiating these pathophysiologies. Finally, the metabolic changes in the recovering patients shift towards the LPS-induced response pattern strengthening the notion that the metabolic, as well as transcriptional responses, characteristic to the endotoxemia model represent necessary and “healthy” responses to infectious stimuli.

Electronic supplementary material

The online version of this article (doi:10.1186/s13054-015-0783-2) contains supplementary material, which is available to authorized users.

Temporal metabolic profiling of plasma during endotoxemia in humans

Endotoxemia induced by the administration of low-dose lipopolysaccharide (LPS) to healthy human volunteers is a well-established experimental protocol and has served as a reproducible platform for investigating the responses to systemic inflammation. Because metabolic composition of a tissue or body fluid is uniquely altered by stimuli and provides information about the dominant regulatory mechanisms at various cellular processes, understanding the global metabolic response to systemic inflammation constitutes a major part in this investigation complementing the studies undertaken so far in both clinical and systems biology fields. This article communicates the first proof-of-principle metabonomic analysis, which comprised global biochemical profiles in human plasma samples from healthy subjects given intravenous endotoxin at 2 ng/kg. Concentrations of a total of 366 plasma biochemicals were determined in archived blood samples collected from 15 endotoxin-treated subjects at five time points within 24 h after treatment and compared with control samples collected from four saline-treated subjects. Principal component analysis within this data set determined the sixth hour as a critical time point separating development and recovery phases of the LPS-induced metabolic changes. Consensus clustering of the differential metabolites identified two distinct subsets of metabolites that displayed common coherent profiles with opposing directionality. The first group of metabolites, which were mostly associated with pathways related to lipid metabolism, was upregulated within the first 6 h and downregulated by the 24th hour following LPS administration. The second group of metabolites, in contrast, was first downregulated until the sixth hour, then upregulated. Metabolites in this group were predominantly amino acids or their derivatives. In summary, nontargeted biochemical profiling and unsupervised multivariate analyses highlighted the prominent roles of lipid and protein metabolism in regulating the response to systemic inflammation while also revealing their dynamics in opposite directions.

Identification of the key differential transcriptional responses of human whole blood following TLR2 or TLR4 ligation in-vitro

The use of human whole blood for transcriptomic analysis has potential advantages over the use of isolated immune cells for studying the transcriptional response to pathogens and their products. Whole blood stimulation can be carried out in a laboratory without the expertise or equipment to isolate immune cells from blood, with the added advantage of being able to undertake experiments using very small volumes of blood. Toll like receptors (TLRs) are a family of pattern recognition receptors which recognise highly conserved microbial products. Using the TLR2 ligand (Pam3CSK4) and the TLR4 ligand (LPS), human whole blood was stimulated for 0, 1, 3, 6, 12 or 24 hours at which times mRNA was isolated and a comparative microarray was undertaken. A common NFκB transcriptional programme was identified following both TLR2 and TLR4 ligation which peaked at between 3 to 6 hours including upregulation of many of the NFκB family members. In contrast an interferon transcriptional response was observed following TLR4 but not TLR2 ligation as early as 1 hour post stimulation and peaking at 6 hours. These results recapitulate the findings observed in previously published studies using isolated murine and human myeloid cells indicating that in vitro stimulated human whole blood can be used to interrogate the early transcriptional kinetic response of innate cells to TLR ligands. Our study demonstrates that a transcriptomic analysis of mRNA isolated from human whole blood can delineate both the temporal response and the key transcriptional differences following TLR2 and TLR4 ligation.

Bioinformatics analysis of transcriptional regulation of circadian genes in rat liver

Background

The circadian clock is a critical regulator of biological functions controlling behavioral, physiological and biochemical processes. Because the liver is the primary regulator of metabolites within the mammalian body and the disruption of circadian rhythms in liver is associated with severe illness, circadian regulators would play a strong role in maintaining liver function. However, the regulatory structure that governs circadian dynamics within the liver at a transcriptional level remains unknown. To explore this aspect, we analyzed hepatic transcriptional dynamics in Sprague-Dawley rats over a period of 24 hours to assess the genome-wide responses.

Results

Using an unsupervised consensus clustering method, we identified four major gene expression clusters, corresponding to central carbon and nitrogen metabolism, membrane integrity, immune function, and DNA repair, all of which have dynamics which suggest regulation in a circadian manner. With the assumption that transcription factors (TFs) that are differentially expressed and contain CLOCK:BMAL1 binding sites on their proximal promoters are likely to be clock-controlled TFs, we were able to use promoter analysis to putatively identify additional clock-controlled TFs besides PARF and RORA families. These TFs are both functionally and temporally related to the clusters they regulate. Furthermore, we also identified significant sets of clock TFs that are potentially transcriptional regulators of gene clusters.

Conclusions

All together, we were able to propose a regulatory structure for circadian regulation which represents alternative paths for circadian control of different functions within the liver. Our prediction has been affirmed by functional and temporal analyses which are able to extend for similar studies.

Tissue-specific gene expression and regulation in liver and muscle following chronic corticosteroid administration

Although corticosteroids (CSs) affect gene expression in multiple tissues, the array of genes that are regulated by these catabolic steroids is diverse, highly tissue specific, and depends on their functions in the tissue. Liver has many important functions in performing and regulating diverse metabolic processes. Muscle, in addition to its mechanical role, is critical in maintaining systemic energy homeostasis and accounts for about 80% of insulin-directed glucose disposal. Consequently, a better understanding of CS pharmacogenomic effects in these tissues would provide valuable information regarding the tissue-specificity of transcriptional dynamics, and would provide insights into the underlying molecular mechanisms of action for both beneficial and detrimental effects. We performed an integrated analysis of transcriptional data from liver and muscle in response to methylprednisolone (MPL) infusion, which included clustering and functional annotation of clustered gene groups, promoter extraction and putative transcription factor (TF) identification, and finally, regulatory closeness (RC) identification. This analysis allowed the identification of critical transcriptional responses and CS-responsive functions in liver and muscle during chronic MPL administration, the prediction of putative transcriptional regulators relevant to transcriptional responses of CS-affected genes which are also potential secondary bio-signals altering expression levels of target-genes, and the exploration of the tissue-specificity and biological significance of gene expression patterns, CS-responsive functions, and transcriptional regulation. The analysis provided an integrated description of the genomic and functional effects of chronic MPL infusion in liver and muscle.

Inhibitor of nuclear factor-κB, SN50, attenuates lipopolysaccharide-induced lung injury in an isolated and perfused rat lung model

NF-κB cell permeable inhibitory peptide (SN50) inhibits translocation of nuclear factor-κB (NF-κB) and production of inflammatory cytokines that are implicated in lipopolysaccharide (LPS)-induced lung injury (LPSLI). However, the protective effect of SN50 in LPSLI is unclear. We explored the cellular and molecular mechanisms of SN50 treatment in LPSLI. LPSLI was induced by intratracheal instillation of 10 mg/kg LPS using an isolated and perfused rat lung model. SN50 was administered in the perfusate 15 minutes before LPS was administered. Hemodynamics, lung histologic change, inflammatory responses, and activation of apoptotic pathways were evaluated. After LPSLI, increased pulmonary vascular permeability and lung weight gain was observed. The levels of interleukin (IL)-1β, tumor necrosis factor (TNF)-α, myeloperoxidase, and macrophage inflammatory protein-2 increased in bronchoalveolar lavage fluids. Lung-tissue expression of TNF-α, IL-1β, mitogen-activated protein kinases (MAPKs), caspase-3, p-AKT (serine-threonine kinase, also known as protein kinase B), and plasminogen activator inhibitor-1 (PAI-1) was greater in the LPS group compared with controls. Upregulation and activation of NF-κB was associated with increased lung injury in LPSLI. SN50 attenuated the inflammatory responses, including expression of IL-1β, TNF-α, myeloperoxidase, MAPKs, PAI-1, and NF-κB; downregulation of apoptosis indicated by caspase-3 and p-AKT expression was also observed. In addition, SN50 mitigated the increase in the lung weight, pulmonary vascular permeability, and lung injury. In conclusion, LPSLI is associated with inflammatory responses, apoptosis, and coagulation. NF-κB is an important therapeutic target in the treatment of LPSLI. SN50 inhibits translocation of NF-κB and attenuates LPSLI.

Dynamics of hepatic gene expression profile in a rat cecal ligation and puncture model

BACKGROUND

Sepsis remains a major clinical challenge in intensive care units. The difficulty in developing new and more effective treatments for sepsis exemplifies our incomplete understanding of the underlying pathophysiology of it. One of the more widely used rodent models for studying polymicrobial sepsis is cecal ligation and puncture (CLP). While a number of CLP studies investigated the ensuing systemic inflammatory response, they usually focus on a single time point post-CLP and therefore fail to describe the dynamics of the response. Furthermore, previous studies mostly use surgery without infection (herein referred to as sham CLP, SCLP) as a control for the CLP model, however, SCLP represents an aseptic injurious event that also stimulates a systemic inflammatory response. Thus, there is a need to better understand the dynamics and expression patterns of both injury- and sepsis-induced gene expression alterations to identify potential regulatory targets. In this direction, we characterized the response of the liver within the first 24 h in a rat model of SCLP and CLP using a time series of microarray gene expression data.

METHODS

Rats were randomly divided into three groups: sham, SCLP, and CLP. Rats in SCLP group are subjected to laparotomy, cecal ligation, and puncture while those in CLP group are subjected to the similar procedures without cecal ligation and puncture. Animals were saline resuscitated and sacrificed at defined time points (0, 2, 4, 8, 16, and 24 h). Liver tissues were explanted and analyzed for their gene expression profiles using microarray technology. Unoperated animals (Sham) serve as negative controls. After identifying differentially expressed probesets between sham and SCLP or CLP conditions over time, the concatenated data sets corresponding to these differentially expressed probesets in sham and SCLP or CLP groups were combined and analyzed using a "consensus clustering" approach. Promoters of genes that share common characteristics were extracted and compared with gene batteries comprised of co-expressed genes to identify putatative transcription factors, which could be responsible for the co-regulation of those genes.

RESULTS

The SCLP/CLP genes whose expression patterns significantly changed compared with sham over time were identified, clustered, and finally analyzed for pathway enrichment. Our results indicate that both CLP and SCLP triggered the activation of a proinflammatory response, enhanced synthesis of acute-phase proteins, increased metabolism, and tissue damage markers. Genes triggered by CLP, which can be directly linked to bacteria removal functions, were absent in SCLP injury. In addition, genes relevant to oxidative stress induced damage were unique to CLP injury, which may be due to the increased severity of CLP injury versus SCLP injury. Pathway enrichment identified pathways with similar functionality but different dynamics in the two injury models, indicating that the functions controlled by those pathways are under the influence of different transcription factors and regulatory mechanisms. Putatively identified transcription factors, notably including cAMP response element-binding (CREB), nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), and signal transducer and activator of transcription (STAT), were obtained through analysis of the promoter regions in the SCLP/CLP genes. Our results show that while transcription factors such as NF-κB, homeodomain transcription factor (HOMF), and GATA transcription factor (GATA) were common in both injuries for the IL-6 signaling pathway, there were many other transcription factors associated with that pathway which were unique to CLP, including forkhead (FKHD), hairy/enhancer of split family (HESF), and interferon regulatory factor family (IRFF). There were 17 transcription factors that were identified as important in at least two pathways in the CLP injury, but only seven transcription factors with that property in the SCLP injury. This also supports the hypothesis of unique regulatory modules that govern the pathways present in both the CLP and SCLP response.

CONCLUSIONS

By using microarrays to assess multiple genes in a high throughput manner, we demonstrate that an inflammatory response involving different dynamics and different genes is triggered by SCLP and CLP. From our analysis of the CLP data, the key characteristics of sepsis are a proinflammatory response, which drives hypermetabolism, immune cell activation, and damage from oxidative stress. This contrasts with SCLP, which triggers a modified inflammatory response leading to no immune cell activation, decreased detoxification potential, and hyper metabolism. Many of the identified transcription factors that drive the CLP-induced response are not found in the SCLP group, suggesting that SCLP and CLP induce different types of inflammatory responses via different regulatory pathways.