Analysis of the TGFβ-induced program in primary airway epithelial cells shows essential role of NF-κB/RelA signaling network in type II epithelial mesenchymal transition

Background

The airway epithelial cell plays a central role in coordinating the pulmonary response to injury and inflammation. Here, transforming growth factor-β (TGFβ) activates gene expression programs to induce stem cell-like properties, inhibit expression of differentiated epithelial adhesion proteins and express mesenchymal contractile proteins. This process is known as epithelial mesenchymal transition (EMT); although much is known about the role of EMT in cellular metastasis in an oncogene-transformed cell, less is known about Type II EMT, that occurring in normal epithelial cells. In this study, we applied next generation sequencing (RNA-Seq) in primary human airway epithelial cells to understand the gene program controlling Type II EMT and how cytokine-induced inflammation modifies it.

Results

Generalized linear modeling was performed on a two-factor RNA-Seq experiment of 6 treatments of telomerase immortalized human small airway epithelial cells (3 replicates). Using a stringent cut-off, we identified 3,478 differentially expressed genes (DEGs) in response to EMT. Unbiased transcription factor enrichment analysis identified three clusters of EMT regulators, one including SMADs/TP63 and another NF-κB/RelA. Surprisingly, we also observed 527 of the EMT DEGs were also regulated by the TNF-NF-κB/RelA pathway. This Type II EMT program was compared to Type III EMT in TGFβ stimulated A549 alveolar lung cancer cells, revealing significant functional differences. Moreover, we observe that Type II EMT modifies the outcome of the TNF program, reducing IFN signaling and enhancing integrin signaling. We confirmed experimentally that TGFβ-induced the NF-κB/RelA pathway by observing a 2-fold change in NF-κB/RelA nuclear translocation. A small molecule IKK inhibitor blocked TGFβ-induced core transcription factor (SNAIL1, ZEB1 and Twist1) and mesenchymal gene (FN1 and VIM) expression.

Conclusions

These data indicate that NF-κB/RelA controls a SMAD-independent gene network whose regulation is required for initiation of Type II EMT. Type II EMT dramatically affects the induction and kinetics of TNF-dependent gene networks.

Electronic supplementary material

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

The Multidisciplinary Translational Team (MTT) Model for Training and Development of Translational Research Investigators

Multiinstitutional research collaborations now form the most rapid and productive project execution structures in the health sciences. Effective adoption of a multidisciplinary team research approach is widely accepted as one mechanism enabling rapid translation of new discoveries into interventions in human health. Although the impact of successful team-based approaches facilitating innovation has been well-documented, its utility for training a new generation of scientists has not been thoroughly investigated. We describe the characteristics of how multidisciplinary translational teams (MTTs) promote career development of translational research scholars through competency building, interprofessional integration, and team-based mentoring approaches. Exploratory longitudinal and outcome assessments from our experience show that MTT membership had a positive effect on the development of translational research competencies, as determined by a self-report survey of 32 scholars. We also observed that all trainees produced a large number of collaborative publications that appeared to be associated with their CTSA association and participation with MTTs. We conclude that the MTT model provides a unique training environment for translational and team-based learning activities, for investigators at early stages of career development.

Abstract 458: Loss of Smooth Muscle α-actin in Mice Results in Thoracic Aortic Aneurysms via Increased Reactive Oxygen Species, Increased Nox4,and Increased Angiotensin II type 1 Receptor-Mediated Signaling

Objective: ACTA2 mutations cause 10-14% of familial thoracic aortic aneurysms and dissections. Mice deficient in smooth muscle α-actin (Acta2-/-) develop root and ascending thoracic aortic enlargement associated with thickening of the aortic media and fragmentation and disarray of elastic fibers. We hypothesized that blocking AT1 activation would block the aortic pathology and prevent aortic enlargement in Acta2-/- mice.

Methods and Results: Beginning at 4 weeks of age, Acta2-/- mice were treated with losartan or placebo (n≥10) for 6 months and echocardiograms were performed at baseline and every other month. The aortic root in Acta2-/- mice was found to undergo progressive dilatation. After 6 months of treatment, there was no difference in the diameter of the aortic root between wild-type (WT) mice and the losartan treated mice (p=0.44). Histologic analysis of Acta2-/- aortas demonstrated medial thickening and fragmentation of elastic fibers which was normalized by treatment with losartan. Gene expression of matrix metalloproteinase-2 and -9 (Mmp2, Mmp9), along with lumican and decorin, interleukin-6 (Il6) and phosphorylation of RelA (a subunit of nuclear factor κB, NF-κB) was increased in Acta2-/- aortas, and was corrected by treatment with losartan. NADPH oxidase 4 and AT1a mRNA was increased in Acta2-/- aortic smooth muscle cells (SMCs) and aortas. Increase of AT1a was blocked by lowering reactive oxygen species (ROS) with N-actetyl cysteine (NAC). Angiotensin II (AngII) dose response studies suggested Acta2-/- aortic SMCs had increased sensitivity to Ang II. Additionally, Acta2-/- SMCs had increased ROS compared to WT by flow cytometry (P<0.05). Increased ROS was also observed in aortic samples derived from Acta2-/- mice. Blocking ROS using NAC attenuated aneurysm formation by reducing the expression of Mmps and IL6 and reducing signaling through AT1 in Acta2-/- mice. Genetic deletion of AT1a attenuated ROS and expression of Mmps and IL6 levels in the Acta2-/-aortas but not to WT levels.

Conclusions: Our results demonstrate that complete loss of α-SMA leads to aortic dilation and pathologic changes by increasing SMC ROS levels, thus increasing sensitivity to AngII, which results in NF-κB activation and increased expression of Mmps and IL6.

Abstract 667: The Role of Aortic Wall-Nuclear Factor- kappa B (NF-kB) Signaling in Formation of Dissecting Aortic Aneurysms

Introduction: Abdominal aortic aneurysm is a devastating disease since it can lead to aortic rupture and instantaneous death. We previously demonstrated that IL-6 secreted from the aortic wall is necessary for the development of abdominal aortic aneurysm and dissection (AAD). Since IL-6 is a NF-kB/RelA dependant gene, we investigated the role of aortic wall- NF-kB/RelA signaling in the development of AAD.

Methods and Results: To test the role of aortic wall-RelA, we utilized Cre-Lox technology to delete RelA from aortic cells. Tamoxifen-inducible, Col1a2-promoter driven Cre mice (Col1a2-Cre) were crossed with mT/mG Cre-reporter mice to determine which aortic cells undergo genetic recombination after Cre activation. Flow cytometry analysis of the aortic wall indicated that 88% of the genetically recombined cells were SMCs and 8% were fibroblasts. Next, RelA floxed (RelA f/f) mice, generated in our lab, were crossed with Col1a2-Cre mice. RelA f/f Cre+ and RelA f/f Cre- were stimulated with tamoxifen for 10 days to generate aortic-RelA deficient (Ao-RelA-/-) or wild-type (Ao-RelA+/+) transgenics. Flow cytometry, qRT-PCR, and immunohistochemistry analysis suggested a depletion of aortic-RelA greater than 60%. To test the role of Ao-RelA in AAD, Ao-RelA -/- (n= 20) and Ao-RelA +/+ (n=14) mice were infused with angiotensin II for 7 days. Surprisingly, 20% of Ao-RelA-/- mice died from development of AAD and aortic rupture while no deaths were observed in Ao-RelA+/+ group. In addition, 40% of Ao-RelA-/- mice developed AAD compared to 14% of Ao-RelA+/+ mice. There was no significant difference in TUNEL staining or ERTR7+ fibroblast population between the two groups.

Conclusion: Our studies suggest that aortic wall-RelA may be necessary for protection from AAD.

Mixed-effects model of epithelial-mesenchymal transition reveals rewiring of signaling networks

The type II epithelial-mesenchymal transition (EMT) produces airway fibrosis and remodeling, contributing to the severity of asthma and chronic obstructive pulmonary disease. While numerous studies have been done on the mechanisms of the transition itself, few studies have investigated the system effects of EMT on signaling networks. Here, we use mixed effects modeling to develop a computational model of phospho-protein signaling data that compares human small airway epithelial cells (hSAECs) with their EMT-transformed counterparts across a series of perturbations with 8 ligands and 5 inhibitors, revealing previously uncharacterized changes in signaling in the EMT state. Strong couplings between menadione, TNFα and TGFβ and their known phospho-substrates were revealed after mixed effects modeling. Interestingly, the overall phospho-protein response was attenuated in EMT, with loss of Mena and TNFα coupling to heat shock protein (HSP)-27. These differences persisted after correction for EMT-induced changes in phospho-protein substrate abundance. Construction of network topology maps showed significant changes between the two cellular states, including a linkage between glycogen synthase kinase (GSK)-3α and small body size/mothers against decapentaplegic (SMAD)2. The model also predicted a loss of p38 mitogen activated protein kinase (p38MAPK)-independent HSP27 signaling, which we experimentally validated. We further characterized the relationship between HSP27 and signal transducers and activators of transcription (STAT)3 signaling, and determined that loss of HSP27 following EMT is only partially responsible for the downregulation of STAT3. These rewired connections represent therapeutic targets that could potentially reverse EMT and restore a normal phenotype to the respiratory mucosa.

Whole transcriptome analysis reveals an 8-oxoguanine DNA glycosylase-1-driven DNA repair-dependent gene expression linked to essential biological processes

Reactive oxygen species inflict oxidative modifications on various biological molecules, including DNA. One of the most abundant DNA base lesions, 8-oxo-7,8-dihydroguanine (8-oxoG) is repaired by 8-oxoguanine DNA glycosylase-1 (OGG1) during DNA base excision repair (OGG1-BER). 8-OxoG accumulation in DNA has been associated with various pathological and aging processes, although its role is unclear. The lack of OGG1-BER in Ogg1(-/-) mice resulted in decreased inflammatory responses and increased susceptibility to infections and metabolic disorders. Therefore, we proposed that OGG1 and/or 8-oxoG base may have a role in immune and homeostatic processes. To test our hypothesis, we challenged mouse lungs with OGG1-BER product 8-oxoG base and changes in gene expression were determined by RNA sequencing and data were analyzed by Gene Ontology and statistical tools. RNA-Seq analysis identified 1592 differentially expressed (≥ 3-fold change) transcripts. The upregulated mRNAs were related to biological processes, including homeostatic, immune-system, macrophage activation, regulation of liquid-surface tension, and response to stimulus. These processes were mediated by chemokines, cytokines, gonadotropin-releasing hormone receptor, integrin, and interleukin signaling pathways. Taken together, these findings point to a new paradigm showing that OGG1-BER plays a role in various biological processes that may benefit the host, but when in excess could be implicated in disease and/or aging processes.

Evolution of Multidisciplinary Translational Teams (MTTs): Insights for Accelerating Translational Innovations

There is growing consensus about the factors critical for development and productivity of multidisciplinary teams, but few studies have evaluated their longitudinal changes. We present a longitudinal study of 10 multidisciplinary translational teams (MTTs), based on team process and outcome measures, evaluated before and after 3 years of CTSA collaboration. Using a mixed methods approach, an expert panel of five judges (familiar with the progress of the teams) independently rated team performance based on four process and four outcome measures, and achieved a rating consensus. Although all teams made progress in translational domains, other process and outcome measures were highly variable. The trajectory profiles identified four categories of team performance. Objective bibliometric analysis of CTSA-supported MTTs with positive growth in process scores showed that these teams tended to have enhanced scientific outcomes and published in new scientific domains, indicating the conduct of innovative science. Case exemplars revealed that MTTs that experienced growth in both process and outcome evaluative criteria also experienced greater innovation, defined as publications in different areas of science. Of the eight evaluative criteria, leadership-related behaviors were the most resistant to the interventions introduced. Well-managed MTTs demonstrate objective productivity and facilitate innovation.

Ataxia telangiectasia mutated kinase mediates NF-κB serine 276 phosphorylation and interferon expression via the IRF7-RIG-I amplification loop in paramyxovirus infection

Respiratory syncytial virus (RSV) is a primary etiological agent of childhood lower respiratory tract disease. Molecular patterns induced by active infection trigger a coordinated retinoic acid-inducible gene I (RIG-I)-Toll-like receptor (TLR) signaling response to induce inflammatory cytokines and antiviral mucosal interferons. Recently, we discovered a nuclear oxidative stress-sensitive pathway mediated by the DNA damage response protein, ataxia telangiectasia mutated (ATM), in cytokine-induced NF-κB/RelA Ser 276 phosphorylation. Here we observe that ATM silencing results in enhanced single-strand RNA (ssRNA) replication of RSVand Sendai virus, due to decreased expression and secretion of type I and III interferons (IFNs), despite maintenance of IFN regulatory factor 3 (IRF3)-dependent IFN-stimulated genes (ISGs). In addition to enhanced oxidative stress, RSV replication enhances foci of phosphorylated histone 2AX variant (γH2AX), Ser 1981 phosphorylation of ATM, and IKKγ/NEMO-dependent ATM nuclear export, indicating activation of the DNA damage response. ATM-deficient cells show defective RSV-induced mitogen and stress-activated kinase 1 (MSK-1) Ser 376 phosphorylation and reduced RelA Ser 276 phosphorylation, whose formation is required for IRF7 expression. We observe that RelA inducibly binds the native IFN regulatory factor 7 (IRF7) promoter in an ATM-dependent manner, and IRF7 inducibly binds to the endogenous retinoic acid-inducible gene I (RIG-I) promoter. Ectopic IRF7 expression restores RIG-I expression and type I/III IFN expression in ATM-silenced cells. We conclude that paramyxoviruses trigger the DNA damage response, a pathway required for MSK1 activation of phospho Ser 276 RelA formation to trigger the IRF7-RIG-I amplification loop necessary for mucosal IFN production. These data provide the molecular pathogenesis for defects in the cellular innate immunity of patients with homozygous ATM mutations.

IMPORTANCE

RNA virus infections trigger cellular response pathways to limit spread to adjacent tissues. This "innate immune response" is mediated by germ line-encoded pattern recognition receptors that trigger activation of two, largely independent, intracellular NF-κB and IRF3 transcription factors. Downstream, expression of protective antiviral interferons is amplified by positive-feedback loops mediated by inducible interferon regulatory factors (IRFs) and retinoic acid inducible gene (RIG-I). Our results indicate that a nuclear oxidative stress- and DNA damage-sensing factor, ATM, is required to mediate a cross talk pathway between NF-κB and IRF7 through mediating phosphorylation of NF-κB. Our studies provide further information about the defects in cellular and innate immunity in patients with inherited ATM mutations.

Molecular classification of outcomes from dengue virus -3 infections

OBJECTIVES

Dengue virus (DENV) infection is a significant risk to over a third of the human population that causes a wide spectrum of illness, ranging from sub-clinical disease to intermediate syndrome of vascular complications called dengue fever complicated (DFC) and severe, dengue hemorrhagic fever (DHF). Methods for discriminating outcomes will impact clinical trials and understanding disease pathophysiology.

STUDY DESIGN

We integrated a proteomics discovery pipeline with a heuristics approach to develop a molecular classifier to identify an intermediate phenotype of DENV-3 infectious outcome.

RESULTS

121 differentially expressed proteins were identified in plasma from DHF vs dengue fever (DF), and informative candidates were selected using nonparametric statistics. These were combined with markers that measure complement activation, acute phase response, cellular leak, granulocyte differentiation and viral load. From this, we applied quantitative proteomics to select a 15 member panel of proteins that accurately predicted DF, DHF, and DFC using a random forest classifier. The classifier primarily relied on acute phase (A2M), complement (CFD), platelet counts and cellular leak (TPM4) to produce an 86% accuracy of prediction with an area under the receiver operating curve of >0.9 for DHF and DFC vs DF.

CONCLUSIONS

Integrating discovery and heuristic approaches to sample distinct pathophysiological processes is a powerful approach in infectious disease. Early detection of intermediate outcomes of DENV-3 will speed clinical trials evaluating vaccines or drug interventions.

Innate inflammation induced by the 8-oxoguanine DNA glycosylase-1-KRAS-NF-κB pathway

8-Oxoguanine-DNA glycosylase-1 (OGG1) is the primary enzyme for repairing 7,8-dihydro-8-oxoguanine (8-oxoG) via the DNA base excision repair pathway (OGG1-BER). Accumulation of 8-oxoG in the genomic DNA leads to genetic instability and carcinogenesis and is thought to contribute to the worsening of various inflammatory and disease processes. However, the disease mechanism is unknown. In this study, we proposed that the mechanistic link between OGG1-BER and proinflammatory gene expression is OGG1's guanine nucleotide exchange factor activity, acquired after interaction with the 8-oxoG base and consequent activation of the small GTPase RAS. To test this hypothesis, we used BALB/c mice expressing or deficient in OGG1 in their airway epithelium and various molecular biological approaches, including active RAS pulldown, reporter and Comet assays, small interfering RNA-mediated depletion of gene expression, quantitative RT-PCR, and immunoblotting. We report that the OGG1-initiated repair of oxidatively damaged DNA is a prerequisite for GDP → GTP exchange, KRAS-GTP-driven signaling via MAP kinases and PI3 kinases and mitogen-stress-related kinase-1 for NF-κB activation, proinflammatory chemokine/cytokine expression, and inflammatory cell recruitment to the airways. Mice deficient in OGG1-BER showed significantly decreased immune responses, whereas a lack of other Nei-like DNA glycosylases (i.e., NEIL1 and NEIL2) had no significant effect. These data unveil a previously unidentified role of OGG1-driven DNA BER in the generation of endogenous signals for inflammation in the innate signaling pathway.

8-Oxoguanine DNA glycosylase-1-mediated DNA repair is associated with Rho GTPase activation and α-smooth muscle actin polymerization

Reactive oxygen species (ROS) are activators of cell signaling and modify cellular molecules, including DNA. 8-Oxo-7,8-dihydroguanine (8-oxoG) is one of the prominent lesions in oxidatively damaged DNA, whose accumulation is causally linked to various diseases and aging processes, whereas its etiological relevance is unclear. 8-OxoG is repaired by the 8-oxoguanine DNA glycosylase-1 (OGG1)-initiated DNA base excision repair (BER) pathway. OGG1 binds free 8-oxoG and this complex functions as an activator of Ras family GTPases. Here we examined whether OGG1-initiated BER is associated with the activation of Rho GTPase and mediates changes in the cytoskeleton. To test this possibility, we induced OGG1-initiated BER in cultured cells and mouse lungs and used molecular approaches such as active Rho pull-down assays, siRNA ablation of gene expression, immune blotting, and microscopic imaging. We found that OGG1 physically interacts with Rho GTPase and, in the presence of 8-oxoG base, increases Rho-GTP levels in cultured cells and lungs, which mediates α-smooth muscle actin (α-SMA) polymerization into stress fibers and increases the level of α-SMA in insoluble cellular/tissue fractions. These changes were absent in cells lacking OGG1. These unexpected data and those showing that 8-oxoG repair is a lifetime process suggest that, via Rho GTPase, OGG1 could be involved in the cytoskeletal changes and organ remodeling observed in various chronic diseases.

TLR4 activation enhances the PD-L1-mediated tolerogenic capacity of colonic CD90+ stromal cells

Signaling via programmed death ligand-1 (PD-L1) and PD-L2 is crucial for maintaining peripheral tolerance. CD90(+) myofibroblasts/fibroblasts (CMFs) are major programmed cell death-1 (PD-1) ligand-expressing cells in normal human colonic mucosa. CMFs suppress activated CD4(+) T cell proliferation via PD-1 ligands. It is not known whether signaling through TLRs contribute to the regulation PD-1 ligands on CMFs upon colonic mucosal tolerance. In this study, we demonstrated that stimulation of TLR4 on human CMFs upregulates PD-L1, but not PD-L2, and reinforces CMF-mediated suppression of CD4(+) T cell proliferation and IFN-γ production. TLR4-mediated upregulation of PD-L1 on CMFs involved NF-κB pathways and was JAK2 and MyD88 dependent. MyD88-dependent stimulation of TLR1/2 and TLR5 also upregulated PD-L1 expression on CMFs in culture. PD-L1 expression was drastically decreased in vivo in the colonic mucosa of mice devoid of MyD88. Induction of MyD88 deficiency in CMFs in fibroblast-specific MyD88 conditional knockout mice resulted in a strong increase in a mucosal IFN-γ expression concomitantly with the abrogation of PD-L1 expression in CMFs under homeostasis and epithelial injury induced by dextran sodium sulfate. Together, these data suggest that MyD88-dependent TLR stimulation of CMFs in the normal colonic mucosa may reinforce these cells' anti-inflammatory capacity and thus contribute to the maintenance of mucosal tolerance.

ATM regulates NF-κB-dependent immediate-early genes via RelA Ser 276 phosphorylation coupled to CDK9 promoter recruitment

Ataxia-telangiectasia mutated (ATM), a member of the phosphatidylinositol 3 kinase-like kinase family, is a master regulator of the double strand DNA break-repair pathway after genotoxic stress. Here, we found ATM serves as an essential regulator of TNF-induced NF-kB pathway. We observed that TNF exposure of cells rapidly induced DNA double strand breaks and activates ATM. TNF-induced ROS promote nuclear IKKγ association with ubiquitin and its complex formation with ATM for nuclear export. Activated cytoplasmic ATM is involved in the selective recruitment of the E3-ubiquitin ligase β-TrCP to phospho-IκBα proteosomal degradation. Importantly, ATM binds and activates the catalytic subunit of protein kinase A (PKAc), ribosmal S6 kinase that controls RelA Ser 276 phosphorylation. In ATM knockdown cells, TNF-induced RelA Ser 276 phosphorylation is significantly decreased. We further observed decreased binding and recruitment of the transcriptional elongation complex containing cyclin dependent kinase-9 (CDK9; a kinase necessary for triggering transcriptional elongation) to promoters of NF-κB-dependent immediate-early cytokine genes, in ATM knockdown cells. We conclude that ATM is a nuclear damage-response signal modulator of TNF-induced NF-κB activation that plays a key scaffolding role in IκBα degradation and RelA Ser 276 phosphorylation. Our study provides a mechanistic explanation of decreased innate immune response associated with A-T mutation.

RIG-I enhanced interferon independent apoptosis upon Junin virus infection

Junin virus (JUNV) is the etiological agent of Argentine hemorrhagic fever (AHF), a human disease with a high case-fatality rate. It is widely accepted that arenaviral infections, including JUNV infections, are generally non-cytopathic. In contrast, here we demonstrated apoptosis induction in human lung epithelial carcinoma (A549), human hepatocarcinoma and Vero cells upon infection with the attenuated Candid#1 strain of, JUNV as determined by phosphatidylserine (PS) translocation, Caspase 3 (CASP3) activation, Poly (ADP-ribose) polymerase (PARP) cleavage and/or chromosomal DNA fragmentation. Moreover, as determined by DNA fragmentation, we found that the pathogenic Romero strain of JUNV was less cytopathic than Candid#1 in human hepatocarcinoma and Vero, but more apoptotic in A549 and Vero E6 cells. Additionally, we found that JUNV-induced apoptosis was enhanced by RIG-I signaling. Consistent with the previously reported role of RIG-I like helicase (RLH) signaling in initiating programmed cell death, we showed that cell death or DNA fragmentation of Candid#1-infected A549 cells was decreased upon siRNA or shRNA silencing of components of RIG-I pathway in spite of increased virus production. Similarly, we observed decreased DNA fragmentation in JUNV-infected human hepatocarcinoma cells deficient for RIG-I when compared with that of RIG-I-competent cells. In addition, DNA fragmentation detected upon Candid#1 infection of type I interferon (IFN)-deficient Vero cells suggested a type I IFN-independent mechanism of apoptosis induction in response to JUNV. Our work demonstrated for the first time apoptosis induction in various cells of mammalian origin in response to JUNV infection and partial mechanism of this cell death.

Systems biology approaches to understanding Epithelial Mesenchymal Transition (EMT) in mucosal remodeling and signaling in asthma

A pathological hallmark of asthma is chronic injury and repair, producing dysfunction of the epithelial barrier function. In this setting, increased oxidative stress, growth factor- and cytokine stimulation, together with extracellular matrix contact produces transcriptional reprogramming of the epithelial cell. This process results in epithelial-mesenchymal transition (EMT), a cellular state associated with loss of epithelial polarity, expression of mesenchymal markers, enhanced mobility and extracellular matrix remodeling. As a result, the cellular biology of the EMT state produces characteristic changes seen in severe, refractory asthma: myofibroblast expansion, epithelial trans-differentiation and subepithelial fibrosis. EMT also induces profound changes in epithelial responsiveness that affects innate immune signaling that may have impact on the adaptive immune response and effectiveness of glucocorticoid therapy in severe asthma. We discuss how this complex phenotype is beginning to be understood using systems biology-level approaches through perturbations coupled with high throughput profiling and computational modeling. Understanding the distinct changes induced by EMT at the systems level may provide translational strategies to reverse the altered signaling and physiology of refractory asthma.

Analysis and predictive modeling of asthma phenotypes

Molecular classification using robust biochemical measurements provides a level of diagnostic precision that is unattainable using indirect phenotypic measurements. Multidimensional measurements of proteins, genes, or metabolites (analytes) can identify subtle differences in the pathophysiology of patients with asthma in a way that is not otherwise possible using physiological or clinical assessments. We overview a method for relating biochemical analyte measurements to generate predictive models of discrete (categorical) clinical outcomes, a process referred to as "supervised classification." We consider problems inherent in wide (small n and large p) high-dimensional data, including the curse of dimensionality, collinearity and lack of information content. We suggest methods for reducing the data to the most informative features. We describe different approaches for phenotypic modeling, using logistic regression, classification and regression trees, random forest and nonparametric regression spline modeling. We provide guidance on post hoc model evaluation and methods to evaluate model performance using ROC curves and generalized additive models. The application of validated predictive models for outcome prediction will significantly impact the clinical management of asthma.

Multiplexed parallel reaction monitoring targeting histone modifications on the QExactive mass spectrometer

Histone acetylation and methylation play an important role in the regulation of gene expression. Irregular patterns of histone global acetylation and methylation have frequently been seen in various diseases. Quantitative analysis of these patterns is of high value for the evaluation of disease development and of outcomes from therapeutic treatment. Targeting histone acetylation and methylation by selected reaction monitoring (SRM) is one of the current quantitative methods. Here, we reported the use of the multiplexed parallel reaction monitoring (PRM) method on the QExactive mass spectrometer to target previously known lysine acetylation and methylation sites of histone H3 and H4 for the purpose of establishing precursor-product pairs for SRM. 55 modified peptides among which 29 were H3 K27/K36 modified peptides were detected from 24 targeted precursor ions included in the inclusion list. The identification was carried out directly from the trypsin digests of core histones that were separated without derivatization on a homemade capillary column packed with Waters YMC ODS-AQ reversed phase materials. Besides documenting the higher-energy c-trap dissociation (HCD) MS(2) spectra of previously known histone H3/H4 acetylated and methylated tryptic peptides, we identified novel H3 K18 methylation, H3 K27 monomethyl/acetyl duel modifications, H2B K23 acetylation, and H4 K20 acetylation in mammalian histones. The information gained from these experiments sets the foundation for quantification of histone modifications by targeted mass spectrometry methods directly from core histone samples.

Dynamic cross talk model of the epithelial innate immune response to double-stranded RNA stimulation: coordinated dynamics emerging from cell-level noise

We present an integrated dynamical cross-talk model of the epithelial innate immune response (IIR) incorporating RIG-I and TLR3 as the two major pattern recognition receptors (PRR) converging on the RelA and IRF3 transcriptional effectors. bioPN simulations reproduce biologically relevant gene-and protein abundance measurements in response to time course, gene silencing and dose-response perturbations both at the population and single cell level. Our computational predictions suggest that RelA and IRF3 are under auto- and cross-regulation. We predict, and confirm experimentally, that RIG-I mRNA expression is controlled by IRF7. We also predict the existence of a TLR3-dependent, IRF3-independent transcription factor (or factors) that control(s) expression of MAVS, IRF3 and members of the IKK family. Our model confirms the observed dsRNA dose-dependence of oscillatory patterns in single cells, with periods of 1-3 hr. Model fitting to time series, matched by knockdown data suggests that the NF-κB module operates in a different regime (with different coefficient values) than in the TNFα-stimulation experiments. In future studies, this model will serve as a foundation for identification of virus-encoded IIR antagonists and examination of stochastic effects of viral replication. Our model generates simulated time series, which reproduce the noisy oscillatory patterns of activity (with 1-3 hour period) observed in individual cells. Our work supports the hypothesis that the IIR is a phenomenon that emerged by evolution despite highly variable responses at an individual cell level.

Diagnostics for Statistical Variable Selection Methods for Prediction of Peptic Ulcer Disease in Helicobacter pylori Infection

BACKGROUND

The development of accurate classification models depends upon the methods used to identify the most relevant variables. The aim of this article is to evaluate variable selection methods to identify important variables in predicting a binary response using nonlinear statistical models. Our goals in model selection include producing non-overfitting stable models that are interpretable, that generate accurate predictions and have minimum bias. This work was motivated by data on clinical and laboratory features of Helicobacter pylori infections obtained from 60 individuals enrolled in a prospective observational study.

RESULTS

We carried out a comprehensive performance comparison of several nonlinear classification models over the H. pylori data set. We compared variable selection results by Multivariate Adaptive Regression Splines (MARS), Logistic Regression with regularization, Generalized Additive Models (GAMs) and Bayesian Variable Selection in GAMs. We found that the MARS model approach has the highest predictive power because the nonlinearity assumptions of candidate predictors are strongly satisfied, a finding demonstrated via deviance chi-square testing procedures in GAMs.

CONCLUSIONS

Our results suggest that the physiological free amino acids citrulline, histidine, lysine and arginine are the major features for predicting H. pylori peptic ulcer disease on the basis of amino acid profiling.

Assessing and evaluating multidisciplinary translational teams: a mixed methods approach

A case report illustrates how multidisciplinary translational teams can be assessed using outcome, process, and developmental types of evaluation using a mixed-methods approach. Types of evaluation appropriate for teams are considered in relation to relevant research questions and assessment methods. Logic models are applied to scientific projects and team development to inform choices between methods within a mixed-methods design. Use of an expert panel is reviewed, culminating in consensus ratings of 11 multidisciplinary teams and a final evaluation within a team-type taxonomy. Based on team maturation and scientific progress, teams were designated as (a) early in development, (b) traditional, (c) process focused, or (d) exemplary. Lessons learned from data reduction, use of mixed methods, and use of expert panels are explored.

Modulation of gene expression regulated by the transcription factor NF-κB/RelA

Modulators (Ms) are proteins that modify the activity of transcription factors (TFs) and influence expression of their target genes (TGs). To discover modulators of NF-κB/RelA, we first identified 365 NF-κB/RelA-binding proteins using liquid chromatography-tandem mass spectrometry (LC-MS/MS). We used a probabilistic model to infer 8349 (M, NF-κB/RelA, TG) triplets and their modes of modulatory action from our combined LC-MS/MS and ChIP-Seq (ChIP followed by next generation sequencing) data, published RelA modulators and TGs, and a compendium of gene expression profiles. Hierarchical clustering of the derived modulatory network revealed functional subnetworks and suggested new pathways modulating RelA transcriptional activity. The modulators with the highest number of TGs and most non-random distribution of action modes (measured by Shannon entropy) are consistent with published reports. Our results provide a repertoire of testable hypotheses for experimental validation. One of the NF-κB/RelA modulators we identified is STAT1. The inferred (STAT1, NF-κB/RelA, TG) triplets were validated by LC-selected reaction monitoring-MS and the results of STAT1 deletion in human fibrosarcoma cells. Overall, we have identified 562 NF-κB/RelA modulators, which are potential drug targets, and clarified mechanisms of achieving NF-κB/RelA multiple functions through modulators. Our approach can be readily applied to other TFs.

Identification of innate immune response endotypes in asthma: implications for personalized medicine

Asthma is an idiopathic disease characterized by episodic inflammation and reversible airway obstruction triggered by exposure to environmental agents. Because this disease is heterogeneous in onset, exacerbations, inflammatory states, and response to therapy, there is intense interest in developing personalized approaches to its management. Of focus in this review, the recognition that a component of the pathophysiology of asthma is mediated by inflammation has implications for understanding its etiology and individualizing its therapy. Despite understanding how Th2 polarization mediates asthma exacerbations by aeroallergen exposure, we do not yet fully understand how RNA virus infections produce asthmatic exacerbations. This review will summarize the explosion of information that has revealed how patterns produced by RNA virus infection trigger the innate immune response (IIR) in sentinel airway cells. When the IIR is triggered, these cells elaborate inflammatory cytokines and protective mucosal interferons whose actions activate long-lived adaptive immunity and limit organismal replication. Recent work has shown the multifaceted way that dysregulation of the IIR is linked to viral-induced exacerbation, steroid insensitivity, and T helper polarization of adaptive immunity. New developments in quantitative proteomics now enable accurate identification of subgroups of individuals that demonstrate activation of IIR ("innate endotype"). Potential applications to clinical research are proposed. Together, these developments open realistic prospects for how identification of the IIR endotype may inform asthma therapy in the future.

8-oxoguanine DNA glycosylase-1 augments proinflammatory gene expression by facilitating the recruitment of site-specific transcription factors

Among the insidious DNA base lesions, 8-oxo-7,8-dihydroguanine (8-oxoG) is one of the most abundant, a lesion that arises through the attack by reactive oxygen species on guanine, especially when located in cis-regulatory elements. 8-oxoG is repaired by the 8-oxoguanine glycosylase 1 (OGG1)-initiated DNA base excision repair pathway. In this study, we investigated whether 8-oxoG repair by OGG1 in promoter regions is compatible with a prompt gene expression and a host innate immune response. For this purpose, we used a mouse model of airway inflammation, supplemented with cell cultures, chromatin immunoprecipitation, small interfering RNA knockdown, real-time PCR, and comet and reporter transcription assays. Our data show that exposure of cells to TNF-α altered cellular redox, increased the 8-oxoG level in DNA, recruited OGG1 to promoter sequences, and transiently inhibited base excision repair of 8-oxoG. Promoter-associated OGG1 then enhanced NF-κB/RelA binding to cis-elements and facilitated recruitment of specificity protein 1, transcription initiation factor II-D, and p-RNA polymerase II, resulting in the rapid expression of chemokines/cytokines and inflammatory cell accumulation in mouse airways. Small interfering RNA depletion of OGG1 or prevention of guanine oxidation significantly decreased TNF-α-induced inflammatory responses. Taken together, these results show that nonproductive binding of OGG1 to 8-oxoG in promoter sequences could be an epigenetic mechanism to modulate gene expression for a prompt innate immune response.

IL-6 regulates extracellular matrix remodeling associated with aortic dilation in a fibrillin-1 hypomorphic mgR/mgR mouse model of severe Marfan syndrome

Background

Development of thoracic aortic aneurysms is the most significant clinical phenotype in patients with Marfan syndrome. An inflammatory response has been described in advanced stages of the disease. Because the hallmark of vascular inflammation is local interleukin‐6 (IL‐6) secretion, we explored the role of this proinflammatory cytokine in the formation of aortic aneurysms and rupture in hypomorphic fibrillin‐deficient mice (mgR/mgR).

Methods and Results

MgR/mgR mice developed ascending aortic aneurysms with significant dilation of the ascending aorta by 12 weeks (2.7±0.1 and 1.3±0.1 for mgR/mgR versus wild‐type mice, respectively; P<0.001). IL‐6 signaling was increased in mgR/mgR aortas measured by increases in IL‐6 and SOCS3 mRNA transcripts (P<0.05) and in cytokine secretion of IL‐6, MCP‐1, and GM‐CSF (P<0.05). To investigate the role of IL‐6 signaling, we generated mgR homozygous mice with IL‐6 deficiency (DKO). The extracellular matrix of mgR/mgR mice showed significant disruption of elastin and the presence of dysregulated collagen deposition in the medial‐adventitial border by second harmonic generation multiphoton autofluorescence microscopy. DKO mice exhibited less elastin and collagen degeneration than mgR/mgR mice, which was associated with decreased activity of matrix metalloproteinase‐9 and had significantly reduced aortic dilation (1.0±0.1 versus 1.6±0.2 mm change from baseline, DKO versus mgR/mgR, P<0.05) that did not affect rupture and survival.

Conclusion

Activation of IL‐6‐STAT3 signaling contributes to aneurysmal dilation in mgR/mgR mice through increased MMP‐9 activity, aggravating extracellular matrix degradation.

Measurement of the innate immune response in the airway

Asthma is an idiopathic disease associated with episodic inflammation and reversible airway obstruction that is triggered by environmental agents. Allergic and infectious agents trigger asthmatic exacerbations through the innate immune response (IIR). The IIR is activated by sentinel cells in the airways to elaborate inflammatory cytokines and protective mucosal interferons whose actions are designed to limit the spread of the organism, as well as to activate the adaptive immune response. We address the structure of the IIR pathway in sentinel cells of the airway and describe observations on its dysregulation. The IIR is triggered in a cell-type specific manner by germline-encoded pathogen recognition receptors (PPRs) including plasma membrane Toll-like receptors (TLRs) and the cytoplasmic Retinoic Acid-inducible Gene (RIG)-I-like RNA helicases, and protein kinase R (PKR). Although their mechanisms of intracellular signaling differ, both pathways converge on a small group of transcriptional effectors, nuclear factor-κB (NF-κB), IFN regulatory factor (IRF), and signal transducer and activator of transcription (STAT). We describe several distinct techniques to quantitate the IIR including assays based on quantitative real-time PCR (Q-RT-PCR) of NF-κB and IRF3-regulated genes, multiplex bead-based analysis of secreted proteins/cytokines and more recent developments in targeted, quantitative selected reaction monitoring (SRM)-mass spectrometry (MS). Application of these methods for quantitation of the IIR will further our understanding of the role of the IIR in asthma and its contribution to disease heterogeneity.

Innate immune response to arenaviral infection: a focus on the highly pathogenic New World hemorrhagic arenaviruses

Arenaviruses are enveloped, negative-stranded RNA viruses that belong to the family Arenaviridae. This diverse family can be further classified into OW (Old World) and NW (New World) arenaviruses based on their antigenicity, phylogeny, and geographical distribution. Many of the NW arenaviruses are highly pathogenic viruses that cause systemic human infections characterized by hemorrhagic fever and/or neurological manifestations, constituting public health problems in their endemic regions. NW arenavirus infection induces a variety of host innate immune responses, which could contribute to the viral pathogenesis and/or influence the final outcome of virus infection in vitro and in vivo. On the other hand, NW arenaviruses have also developed several strategies to counteract the host innate immune response. We will review current knowledge regarding the interplay between the host innate immune response and NW arenavirus infection in vitro and in vivo, with emphasis on viral-encoded proteins and their effect on the type I interferon response.

•

NW arenaviruses induce a variety of host innate immune responses.

•

The arenaviruses have several strategies to counteract host innate immune response.

•

We review the interplay between host innate immunity and the arenavirus infections.

Variable selection methods for developing a biomarker panel for prediction of dengue hemorrhagic fever

BACKGROUND

The choice of selection methods to identify important variables for binary classification modeling is critical to produce stable models that are interpretable, that generate accurate predictions and have minimum bias. This work is motivated by data on clinical and laboratory features of severe dengue infections (dengue hemorrhagic fever, DHF) obtained from 51 individuals enrolled in a prospective observational study of acute human dengue infections.

RESULTS

We carry out a comprehensive performance comparison using several classification models for DHF over the dengue data set. We compared variable selection results by Multivariate Adaptive Regression Splines, Learning Ensemble, Random Forest, Bayesian Moving Averaging, Stochastic Search Variable Selection, and Generalized Regularized Logistics Regression. Model averaging methods (bagging, boosting and ensemble learners) have higher accuracy, but the generalized regularized regression model has the highest predictive power because the linearity assumptions of candidate predictors are strongly satisfied via deviance chi-square testing procedures. Bootstrapping applications for evaluating predictive regression coefficients in regularized regression model are performed.

CONCLUSIONS

Feature reduction methods introduce inherent biases and therefore are data-type dependent. We propose that these limitations can be overcome using an exhaustive approach for searching feature space. Using this approach, our results suggest that IL-10, platelet and lymphocyte counts are the major features for predicting dengue DHF on the basis of blood chemistries and cytokine measurements.

Aortic remodeling after transverse aortic constriction in mice is attenuated with AT1 receptor blockade

OBJECTIVE

Although hypertension is the most common risk factor for thoracic aortic diseases, it is not understood how increased pressures on the ascending aorta lead to aortic aneurysms. We investigated the role of angiotensin II type 1 receptor activation in ascending aortic remodeling in response to increased biomechanical forces using a transverse aortic constriction (TAC) mouse model.

APPROACH AND RESULTS

Two weeks after TAC, the increased biomechanical pressures led to ascending aortic dilatation and thickening of the medial and adventitial layers of the aorta. There was significant adventitial hyperplasia and inflammatory responses in TAC ascending aortas were accompanied by increased adventitial collagen, elevated inflammatory and proliferative markers, and increased cell density attributable to accumulation of myofibroblasts and macrophages. Treatment with losartan significantly blocked TAC-induced vascular inflammation and macrophage accumulation. However, losartan only partially prevented TAC-induced adventitial hyperplasia, collagen accumulation, and ascending aortic dilatation. Increased Tgfb2 expression and phosphorylated-Smad2 staining in the medial layer of TAC ascending aortas were effectively blocked with losartan. In contrast, the increased Tgfb1 expression and adventitial phospho-Smad2 staining were only partially attenuated by losartan. In addition, losartan significantly blocked extracellular signal-regulated kinase activation and reactive oxygen species production in the TAC ascending aorta.

CONCLUSIONS

Inhibition of the angiotensin II type 1 receptor using losartan significantly attenuated the vascular remodeling associated with TAC but did not completely block the increased transforming growth factor-β1 expression, adventitial Smad2 signaling, and collagen accumulation. These results help to delineate the aortic transforming growth factor-β signaling that is dependent and independent of the angiotensin II type 1 receptor after TAC.

Use of theoretical peptide distributions in phosphoproteome analysis

The high mass accuracy and resolution of modern mass spectrometers provides new opportunities to employ theoretical peptide distributions in large-scale proteomic studies. We used theoretical distributions to study noise filtering and mass measurement errors and to examine mass-based differentiation of phosphorylated and nonphosphorylated peptides. Only the monoisotopic mass of the experimental precursor ion was necessary for this analysis. We found that peak deviations can be used to characterize the modification states of peptides in a sample. When applied to large-scale proteomic data sets, the peak deviation distribution can be used to filter chemical/electronic noise for singly charged species. Using peak deviation distributions, it is possible to separate the phosphorylated peptides from the nonphosphorylated peptides, enabling evaluation of the phosphoproteome content of a sample. Because this approach is simple, with light computational requirements, the analysis of theoretical peptide distributions has a significant potential for application to phosphoproteome analyses. For our studies we used publicly available data sets from three large-scale proteomic studies.

Applications of selected reaction monitoring (SRM)-mass spectrometry (MS) for quantitative measurement of signaling pathways

Quantitative measurement of the major regulatory proteins in signaling networks poses several technical challenges, including low abundance, the presence of post-translational modifications (PTMs), and the lack of suitable affinity detection reagents. Using the innate immune response (IIR) as a model signaling pathway, we illustrate the approach of stable isotope dilution (SID)-selected reaction monitoring (SRM)-mass spectrometry (MS) assays for quantification of low abundance signaling proteins. A work flow for SID-SRM-MS assay development is established for proteins with experimentally observed MS spectra and for those without. Using the interferon response factor (IRF)-3 transcription factor as an example, we illustrate the steps in high responding signature peptide identification, SID-SRM-MS assay optimization, and evaluation. SRM assays for normalization of IIR abundance to invariant housekeeping proteins are presented. We provide an example of SID-SRM assay development for post-translational modification (PTM) detection using an activating phospho-Ser modified NF-κB/RelA transcription factor, and describe challenges inherent in PTM-SID-SRM-MS assay development. Application of highly qualified quantitative, SID-SRM-MS assays will enable a systems-level approach to understanding the dynamics and kinetics of signaling in host cells, such as the IIR.

A probabilistic approach to learn chromatin architecture and accurate inference of the NF-κB/RelA regulatory network using ChIP-Seq

Using nuclear factor-κB (NF-κB) ChIP-Seq data, we present a framework for iterative learning of regulatory networks. For every possible transcription factor-binding site (TFBS)-putatively regulated gene pair, the relative distance and orientation are calculated to learn which TFBSs are most likely to regulate a given gene. Weighted TFBS contributions to putative gene regulation are integrated to derive an NF-κB gene network. A de novo motif enrichment analysis uncovers secondary TFBSs (AP1, SP1) at characteristic distances from NF-κB/RelA TFBSs. Comparison with experimental ENCODE ChIP-Seq data indicates that experimental TFBSs highly correlate with predicted sites. We observe that RelA-SP1-enriched promoters have distinct expression profiles from that of RelA-AP1 and are enriched in introns, CpG islands and DNase accessible sites. Sixteen novel NF-κB/RelA-regulated genes and TFBSs were experimentally validated, including TANK, a negative feedback gene whose expression is NF-κB/RelA dependent and requires a functional interaction with the AP1 TFBSs. Our probabilistic method yields more accurate NF-κB/RelA-regulated networks than a traditional, distance-based approach, confirmed by both analysis of gene expression and increased informativity of Genome Ontology annotations. Our analysis provides new insights into how co-occurring TFBSs and local chromatin context orchestrate activation of NF-κB/RelA sub-pathways differing in biological function and temporal expression patterns.

Short-term bed rest increases TLR4 and IL-6 expression in skeletal muscle of older adults

Bed rest induces significant loss of leg lean mass in older adults. Systemic and tissue inflammation also accelerates skeletal muscle loss, but it is unknown whether inflammation is associated to inactivity-induced muscle atrophy in healthy older adults. We determined if short-term bed rest increases toll-like receptor 4 (TLR4) signaling and pro-inflammatory markers in older adult skeletal muscle biopsy samples. Six healthy, older adults underwent seven consecutive days of bed rest. Muscle biopsies (vastus lateralis) were taken after an overnight fast before and at the end of bed rest. Serum cytokine expression was measured before and during bed rest. TLR4 signaling and cytokine mRNAs associated with pro- and anti-inflammation and anabolism were measured in muscle biopsy samples using Western blot analysis and qPCR. Participants lost ∼4% leg lean mass with bed rest. We found that after bed rest, muscle levels of TLR4 protein expression and interleukin-6 (IL-6), nuclear factor-κB1, interleukin-10, and 15 mRNA expression were increased after bed rest (P < 0.05). Additionally, the cytokines interferon-γ, and macrophage inflammatory protein-1β, were elevated in serum samples following bed rest (P < 0.05). We conclude that short-term bed rest in older adults modestly increased some pro- and anti-inflammatory cytokines in muscle samples while systemic changes in pro-inflammatory cytokines were mostly absent. Upregulation of TLR4 protein content suggests that bed rest in older adults increases the capacity to mount an exaggerated, and perhaps unnecessary, inflammatory response in the presence of specific TLR4 ligands, e.g., during acute illness.

Inducible tumor necrosis factor (TNF) receptor-associated factor-1 expression couples the canonical to the non-canonical NF-κB pathway in TNF stimulation

The NF-κB transcription factor mediates the inflammatory response through distinct (canonical and non-canonical) signaling pathways. The mechanisms controlling utilization of either of these pathways are largely unknown. Here we observe that TNF stimulation induces delayed NF-κB2/p100 processing and investigate the coupling mechanism. TNF stimulation induces TNF-associated factor-1 (TRAF-1) that directly binds NF-κB-inducing kinase (NIK) and stabilizes it from degradation by disrupting its interaction with TRAF2·cIAP2 ubiquitin ligase complex. We show that TRAF1 depletion prevents TNF-induced NIK stabilization and reduces p52 production. To further examine the interactions of TRAF1 and NIK with NF-κB2/p100 processing, we mathematically modeled TRAF1·NIK as a coupling signaling complex and validated computational inference by siRNA knockdown to show non-canonical pathway activation is dependent not only on TRAF1 induction but also NIK stabilization by forming TRAF1·NIK complex. Thus, these integrated computational-experimental studies of TNF-induced TRAF1 expression identified TRAF1·NIK as a central complex linking canonical and non-canonical pathways by disrupting the TRAF2-cIAP2 ubiquitin ligase complex. This feed-forward kinase pathway is essential for the activation of non-canonical pathway.

Interleukin-6-signal transducer and activator of transcription-3 signaling mediates aortic dissections induced by angiotensin II via the T-helper lymphocyte 17-interleukin 17 axis in C57BL/6 mice

OBJECTIVE

Dysregulated angiotensin II (Ang II) signaling induces local vascular interleukin-6 (IL-6) secretion, producing leukocyte infiltration and life-threatening aortic dissections. Precise mechanisms by which IL-6 signaling induces leukocyte recruitment remain unknown. T-helper 17 lymphocytes (Th17) have been implicated in vascular pathology, but their role in the development of aortic dissections is poorly understood. Here, we tested the relationship of IL-6-signal transducer and activator of transcription-3 signaling with Th17-induced inflammation in the formation of Ang II-induced dissections in C57BL/6 mice.

APPROACH AND RESULTS

Ang II infusion induced aortic dissections and CD4(+)-interleukin 17A (IL-17A)-expressing Th17 cell accumulation in C57BL/6 mice. A blunted local Th17 activation, macrophage recruitment, and reduced incidence of aortic dissections were seen in IL-6(-/-) mice. To determine the pathological roles of Th17 lymphocytes, we treated Ang II-infused mice with IL-17A-neutralizing antibody or infused Ang II in genetically deficient IL-17A mice and found decreased aortic chemokine monocytic chemotactic protein-1 production and macrophage recruitment, leading to a reduction in aortic dissections. This effect was independent of blood pressure in IL-17A-neutralizing antibody experiment. Application of a cell-permeable signal transducer and activator of transcription-3 inhibitor to downregulate the IL-6 pathway decreased aortic dilation and Th17 cell recruitment. We also observed increased aortic Th17 infiltration and IL-17 mRNA expression in patients with thoracic aortic dissections. Finally, we found that Ang II-mediated aortic dissections occurred independent of blood pressure changes.

CONCLUSIONS

Our results indicate that the IL-6-signal transducer and activator of transcription-3 signaling pathway converges on Th17 recruitment and IL-17A signaling upstream of macrophage recruitment, mediating aortic dissections.

The Culture of Translational Science Research: Participants' Stories

We apply a symbolic interactionist framework and a qualitative methodology to the examination of the everyday reality of translational science research (TSR). This is a growing scientific movement that aims to facilitate the efficient application of basic research to clinical service design and delivery. We describe the emerging culture of translational research at a mid-size medical center that received a Clinical and Translational Science Award from the National Institutes of Health. The stories related by scientists, clinicians, and students in interviews indicate that they make sense of the emerging inter- and cross-disciplinary, team-oriented culture of TSR through the refinement and redefinition of the significant symbols that inform their work while they attempt to master translational research by addressing the dilemmas it produces for them and their work. We see the strength, currency, adaptability, and energy of the core self-definition of "scientist" to be significant in shaping the emerging culture of translational research. We conclude by celebrating the value of interpretive ethnography for evaluation research.

Quantitation of the dynamic profiles of the innate immune response using multiplex selected reaction monitoring-mass spectrometry

The innate immune response (IIR) is a coordinated intracellular signaling network activated by the presence of pathogen-associated molecular patterns that limits pathogen spread and induces adaptive immunity. Although the precise temporal activation of the various arms of the IIR is a critical factor in the outcome of a disease, currently there are no quantitative multiplex methods for its measurement. In this study, we investigate the temporal activation pattern of the IIR in response to intracellular double-stranded RNA stimulation using a quantitative 10-plex stable isotope dilution-selected reaction monitoring-MS assay. We were able to observe rapid activation of both NF-κB and IRF3 signaling arms, with IRF3 demonstrating a transient response, whereas NF-κB underwent a delayed secondary amplification phase. Our measurements of the NF-κB-IκBα negative feedback loop indicate that about 20% of IκBα in the unstimulated cell is located within the nucleus and represents a population that is rapidly degraded in response to double-stranded RNA. Later in the time course of stimulation, the nuclear IκBα pool is repopulated first prior to its cytoplasmic accumulation. Examination of the IRF3 pathway components shows that double-stranded RNA induces initial consumption of the RIG-I PRR and the IRF3 kinase (TBK1). Stable isotope dilution-selected reaction monitoring-MS measurements after siRNA-mediated IRF3 or RelA knockdown suggests that a low nuclear threshold of NF-κB is required for inducing target gene expression, and that there is cross-inhibition of the NF-κB and IRF3 signaling arms. Finally, we were able to measure delayed noncanonical NF-κB activation by quantifying the abundance of the processed (52 kDa) NF-κB2 subunit in the nucleus. We conclude that quantitative proteomics measurement of the individual signaling arms of the IIR in response to system perturbations is significantly enabled by stable isotope dilution-selected reaction monitoring-MS-based quantification, and that this technique will reveal novel insights into the dynamics and connectivity of the IIR.

A structured approach to predictive modeling of a two-class problem using multidimensional data sets

Biological experiments in the post-genome era can generate a staggering amount of complex data that challenges experimentalists to extract meaningful information. Increasingly, the success of an appropriately controlled experiment relies on a robust data analysis pipeline. In this paper, we present a structured approach to the analysis of multidimensional data that relies on a close, two-way communication between the bioinformatician and experimentalist. A sequential approach employing data exploration (visualization, graphical and analytical study), pre-processing, feature reduction and supervised classification using machine learning is presented. This standardized approach is illustrated by an example from a proteomic data analysis that has been used to predict the risk of infectious disease outcome. Strategies for model selection and post hoc model diagnostics are presented and applied to the case illustration. We discuss some of the practical lessons we have learned applying supervised classification to multidimensional data sets, one of which is the importance of feature reduction in achieving optimal modeling performance.

The IL-6 trans-signaling-STAT3 pathway mediates ECM and cellular proliferation in fibroblasts from hypertrophic scar

The molecular mechanisms behind the pathogenesis of post-burn hypertrophic scar (HS) remain unclear. Here, we investigate the role of interleukin-6 (IL-6) trans-signaling-STAT3 pathway in HS fibroblasts (HSF) derived from burned-induced HS skin. HSF showed increased Tyr 705 STAT3 phosphorylation over normal fibroblast (NF) after IL-6•IL-6Rα stimulation by immunoassays. The endogenous STAT3 target gene, SOCS3, was upregulated in HSF and showed increased STAT3 binding on its promoter relative to NF in Chromatin Immunoprecipitation assay. We observed that the cell surface signaling transducer glycoprotein 130 is upregulated in HSF using Q-RT-PCR and flow cytometry. The production of excessive extracellular matrix (ECM), including the expression of alpha2 (1) procollagen (Col1A2) and fibronectin 1 (FN) were seen in HSFs. A STAT3 peptide inhibitor abrogated FN and Col1A2 gene expression in HSF indicating involvement of STAT3 in ECM production. The cellular proliferation markers Cyclin D1, Bcl-Xl and c-Myc were also upregulated in HSF and knockdown of STAT3 by siRNA attenuated c-Myc expression indicating the essential role of STAT3 in fibroblast proliferation. Taken together, our results suggest that the IL-6-trans-signaling-STAT3 pathway may play an integral role in HS pathogenesis and disruption of this pathway could be a potential therapeutic strategy for the treatment of burn-induced HS.

The CTSA as an exemplar framework for developing multidisciplinary translational teams

Translational science requires that scientists from multiple disciplines work together to improve the prevention, diagnosis, and treatment of human disease. Although a literature exists on the design and management of multidisciplinary teams, little has been written on multidisciplinary translational teams (MTTs). MTTs are distinct hybrid entities, with goals taken from both industry and academic models. We identified 30 design factors in 10 domains from a literature survey relevant to our MTT model: specific goals, structures, and processes. These dimensions were adapted to our own institutional environment in the selection and management of 11 MTTs that exploited resources of University of Texas Medical Branch (UTMB) Clinical and Translational Sciences Awards (CTSA). Case illustrations of two specific MTTs illustrate some of the challenges encountered and opportunities realized in terms of education and scientific advances. Network depiction of disciplinarity indicated that CTSA KRs and CTSA leadership contributed to discipline diversity especially in small (or nascent) MTTs. A separate depiction of MTT-KR utilization indicated that data analysis, translational technologies, and novel methods were heavily utilized by MTTs, whereas other KRs contributed significant effort to infrastructure development. We conclude that the CTSA can provide a rich infrastructural framework and scientific environment for the development of successful MTTs.

Junín virus infection activates the type I interferon pathway in a RIG-I-dependent manner

Junín virus (JUNV), an arenavirus, is the causative agent of Argentine hemorrhagic fever, an infectious human disease with 15-30% case fatality. The pathogenesis of AHF is still not well understood. Elevated levels of interferon and cytokines are reported in AHF patients, which might be correlated to the severity of the disease. However the innate immune response to JUNV infection has not been well evaluated. Previous studies have suggested that the virulent strain of JUNV does not induce IFN in human macrophages and monocytes, whereas the attenuated strain of JUNV was found to induce IFN response in murine macrophages via the TLR-2 signaling pathway. In this study, we investigated the interaction between JUNV and IFN pathway in human epithelial cells highly permissive to JUNV infection. We have determined the expression pattern of interferon-stimulated genes (ISGs) and IFN-β at both mRNA and protein levels during JUNV infection. Our results clearly indicate that JUNV infection activates the type I IFN response. STAT1 phosphorylation, a downstream marker of activation of IFN signaling pathway, was readily detected in JUNV infected IFN-competent cells. Our studies also demonstrated for the first time that RIG-I was required for IFN production during JUNV infection. IFN activation was detected during infection by either the virulent or attenuated vaccine strain of JUNV. Curiously, both virus strains were relatively insensitive to human IFN treatment. Our studies collectively indicated that JUNV infection could induce host type I IFN response and provided new insights into the interaction between JUNV and host innate immune system, which might be important in future studies on vaccine development and antiviral treatment.

Proteomics improves the prediction of burns mortality: results from regression spline modeling

Prediction of mortality in severely burned patients remains unreliable. Although clinical covariates and plasma protein abundance have been used with varying degrees of success, the triad of burn size, inhalation injury, and age remains the most reliable predictor. We investigated the effect of combining proteomics variables with these three clinical covariates on prediction of mortality in burned children. Serum samples were collected from 330 burned children (burns covering >25% of the total body surface area) between admission and the time of the first operation for clinical chemistry analyses and proteomic assays of cytokines. Principal component analysis revealed that serum protein abundance and the clinical covariates each provided independent information regarding patient survival. To determine whether combining proteomics with clinical variables improves prediction of patient mortality, we used multivariate adaptive regression splines, because the relationships between analytes and mortality were not linear. Combining these factors increased overall outcome prediction accuracy from 52% to 81% and area under the receiver operating characteristic curve from 0.82 to 0.95. Thus, the predictive accuracy of burns mortality is substantially improved by combining protein abundance information with clinical covariates in a multivariate adaptive regression splines classifier, a model currently being validated in a prospective study.

The dependence of expression of NF-κB-dependent genes: statistics and evolutionary conservation of control sequences in the promoter and in the 3' UTR

BACKGROUND

The NF-κB family plays a prominent role in the innate immune response, cell cycle activation or cell apoptosis. Upon stimulation by pathogen-associated patterns, such as viral RNA a kinase cascade is activated, which strips the NF-κB of its inhibitor IκBα molecule and allows it to translocate into the nucleus. Once in the nucleus, it activates transcription of approximately 90 genes whose kinetics of expression differ relative to when NF-κB translocates into the nucleus, referred to as Early, Middle and Late genes. It is not obvious what mechanism is responsible for segregation of the genes' timing of transcriptional response.

RESULTS

It is likely that the differences in timing are due, in part, to the number and type of transcription factor binding sites (TFBS), required for NF-κB itself as well as for the putative cofactors, in the Early versus Late genes. We therefore applied an evolutionary analysis of conserved TFBS. We also examined whether transcription dynamic was related to the presence of AU-rich elements (ARE) located in 3'UTR of the mRNA because recent studies have shown that the presence of AREs is associated with rapid gene induction. We found that Early genes were significantly enriched in NF-κB binding sites occurring in evolutionarily conserved domains compared to genes in the Late group. We also found that Early genes had significantly greater number of ARE sequences in the 3'UTR of the gene. The similarities observed among the Early genes were seen in comparison with distant species, while the Late genes promoter regions were much more diversified. Based on the promoter structure and ARE content, Middle genes can be divided into two subgroups which show similarities to Early and Late genes respectively.

CONCLUSIONS

Our data suggests that the rapid response of the NF-κB dependent Early genes may be due to both increased gene transcription due to NF-κB loading as well as the contribution of mRNA instability to the transcript profiles. Wider phylogenetic analysis of NF-κB dependent genes provides insight into the degree of cross-species similarity found in the Early genes, opposed to many differences in promoter structure that can be found among the Late genes. These data suggest that activation and expression of the Late genes is much more species-specific than of the Early genes.

A three-component biomarker panel for prediction of dengue hemorrhagic fever

Dengue virus infections are a major cause of morbidity in tropical countries. Early detection of dengue hemorrhagic fever (DHF) may help identify individuals that would benefit from intensive therapy. Predictive modeling was performed using 11 laboratory values of 51 individuals (38 DF and 13 DHF) obtained on initial presentation using logistic regression. We produced a robust model with an area under the curve of 0.9615 that retained IL-10 levels, platelets, and lymphocytes as the major predictive features. A classification and regression tree was developed on these features that were 86% accurate on cross-validation. The IL-10 levels and platelet counts were also identified as the most informative features associated with DHF using a Random Forest classifier. In the presence of polymerase chain reaction-proven acute dengue infections, we suggest a complete blood count and rapid measurement of IL-10 can assist in the triage of potential DHF cases for close follow-up or clinical intervention improving clinical outcome.

Improved mass defect model for theoretical tryptic peptides

Improvements in the mass accuracy and resolution of mass spectrometers have greatly aided mass spectrometry-based proteomics in profiling complex biological mixtures. With the use of innovative bioinformatics approaches, high mass accuracy and resolution information can be used for filtering chemical noise in mass spectral data. Using our recent algorithmic developments, we have generated the mass distributions of all theoretical tryptic peptides composed of 20 natural amino acids and with masses limited to 3.5 kDa. Peptide masses are distributed discretely, with well-defined peak clusters separated by empty or sparsely populated trough regions. Accurate models for peak centers and widths can be used to filter peptide signals from chemical noise. We modeled mass defects, the difference between monoisotopic and nominal masses, and peak centers and widths in the peptide mass distributions. We found that peak widths encompassing 95% of all peptide sequences are substantially smaller than previously thought. The result has implications for filtering out larger stretches of the mass axis. Mass defects of peptides exhibit an oscillatory behavior which is damped at high mass values. The periodicity of the oscillations is about 14 Da which is the most common difference between the masses of the 20 natural amino acids. To determine the effects of amino acid modifications on our findings, we examined the mass distributions of peptides composed of the 20 natural amino acids, oxidized Met, and phosphorylated Ser, Thr, and Tyr. We found that extension of the amino acid set by modifications increases the 95% peak width. Mass defects decrease, reflecting the fact that the average mass defect of natural amino acids is larger than that of oxidized Met. We propose that a new model for mass defects and peak widths of peptides may improve peptide identifications by filtering chemical noise in mass spectral data.

Discovery proteomics and nonparametric modeling pipeline in the development of a candidate biomarker panel for dengue hemorrhagic fever

Secondary dengue viral infection can produce capillary leakage associated with increased mortality known as dengue hemorrhagic fever (DHF). Because the mortality of DHF can be reduced by early detection and intensive support, improved methods for its detection are needed. We applied multidimensional protein profiling to predict outcomes in a prospective dengue surveillance study in South America. Plasma samples taken from initial clinical presentation of acute dengue infection were subjected to proteomics analyses using ELISA and a recently developed biofluid analysis platform. Demographics, clinical laboratory measurements, nine cytokines, and 419 plasma proteins collected at the time of initial presentation were compared between the DF and DHF outcomes. Here, the subject's gender, clinical parameters, two cytokines, and 42 proteins discriminated between the outcomes. These factors were reduced by multivariate adaptive regression splines (MARS) that a highly accurate classification model based on eight discriminant features with an area under the receiver operator curve (AUC) of 0.999. Model analysis indicated that the feature-outcome relationship were nonlinear. Although this DHF risk model will need validation in a larger cohort, we conclude that approaches to develop predictive biomarker models for disease outcome will need to incorporate nonparametric modeling approaches.

Two-step cross-linking for analysis of protein-chromatin interactions

Eukaryotic gene regulation is controlled, in part, by inducible transcription factor-binding regulatory sequences in a tissue-specific and hormone-responsive manner. The development of methods for the analysis of transcription factor interaction within native chromatin has been a significant advance for the systematic analyses of the timing of gene regulation and studies on the effects of chromatin modifying enzymes on promoter accessibility. Chromatin immunoprecipitation (ChIP) is a specific method involving formaldehyde mediated protein-chromatin fixation to preserve the interaction for subsequent target identification. However, the conventional single-step cross-linking technique does not preserve all protein-DNA interactions, especially for transcription factors in hyper-dynamic equilibrium with chromatin or for coactivator interactions. Here, we describe a versatile, efficient "two-step" XChIP method that involves sequential protein-protein fixation followed by protein-DNA fixation. This method has been used successfully for analysis of chromatin binding for transcription factors (NF-κB, STAT3), polymerases (RNA Pol II), coactivators (CBP/p300, CDK9), and chromatin structural proteins (modified histones). Modifications of DNA extraction and sonication suitable for downstream target identification by quantitative genomic PCR and next generation sequencing are described.

How cytokines co-occur across asthma patients: from bipartite network analysis to a molecular-based classification

Asthmatic patients are currently classified as either severe or non-severe based primarily on their response to glucocorticoids. However, because this classification is based on a post-hoc assessment of treatment response, it does not inform the rational staging of disease or therapy. Recent studies in other diseases suggest that a classification which includes molecular information could lead to more accurate diagnoses and prediction of treatment response. We therefore measured cytokine values in bronchoalveolar lavage (BAL) samples of the lower respiratory tract obtained from 83 asthma patients, and used bipartite network visualizations with associated quantitative measures to conduct an exploratory analysis of the co-occurrence of cytokines across patients. The analysis helped to identify three clusters of patients which had a complex but understandable interaction with three clusters of cytokines, leading to insights for a state-based classification of asthma patients. Furthermore, while the patient clusters were significantly different based on key pulmonary functions, they appeared to have no significant relationship to the current classification of asthma patients. These results suggest the need to define a molecular-based classification of asthma patients, which could improve the diagnosis and treatment of this disease.

High throughput short interfering RNA (siRNA) screening of the human kinome identifies novel kinases controlling the canonical nuclear factor-κB (NF-κB) activation pathway

Nuclear factor-κB (NF-κB) is an inducible cytoplasmic transcription factor that plays a role as a master regulator of airway mucosal inflammation. The prototypical ("canonical") NF-κB pathway controls cytoplasmic to nuclear translocation in response to stimulation by the mononuclear cytokine, TNF. Despite intensive investigation, the spectrum of kinases involved in the canonical NF-κB pathway has not yet been systematically determined. Here we have applied a high throughput siRNA-mediated loss-of-function screening assay to identify novel kinases important in TNF-induced NF-κB signaling. Type II A549 epithelial cells stably expressing an IL-8/luciferase reporter gene optimized for high throughput siRNA format (Z' score of 0.65) and siRNAs for 636 human kinases were reverse-transfected and screened in the assay. 36 candidate genes were identified that inhibited TNF signaling with a Z score deviation of <-1.3 in replicate plates. From this group, 11 kinases were selected for independent validation, of which eight were successfully silenced. Six kinases were validated, including ATM, CDK2, -5, and -7, CALM3, MAPAKP5, and MAP3K/MEKK3. The surprising function of ATM in TNF signaling was confirmed where reduced NF-κB/RelA translocation and Ser-276 phosphorylation were seen in ATM(-/-) mouse embryo fibroblasts. These data indicate that ATM is a key regulatory kinase that may control global NF-κB activation in the TNF-induced canonical pathway.

Sources of cell-to-cell variability in canonical nuclear factor-κB (NF-κB) signaling pathway inferred from single cell dynamic images

The canonical nuclear factor-κB (NF-κB) signaling pathway controls a gene network important in the cellular inflammatory response. Upon activation, NF-κB/RelA is released from cytoplasmic inhibitors, from where it translocates into the nucleus, subsequently activating negative feedback loops producing either monophasic or damped oscillatory nucleo-cytoplasmic dynamics. Although the population behavior of the NF-κB pathway has been extensively modeled, the sources of cell-to-cell variability are not well understood. We describe an integrated experimental-computational analysis of NF-κB/RelA translocation in a validated cell model exhibiting monophasic dynamics. Quantitative measures of cellular geometry and total cytoplasmic concentration and translocated RelA amounts were used as priors in Bayesian inference to estimate biophysically realistic parameter values based on dynamic live cell imaging studies of enhanced GFP-tagged RelA in stable transfectants. Bayesian inference was performed on multiple cells simultaneously, assuming identical reaction rate parameters, whereas cellular geometry and initial and total NF-κB concentration-related parameters were cell-specific. A subpopulation of cells exhibiting distinct kinetic profiles was identified that corresponded to differences in the IκBα translation rate. We conclude that cellular geometry, initial and total NF-κB concentration, IκBα translation, and IκBα degradation rates account for distinct cell-to-cell differences in canonical NF-κB translocation dynamics.