Cytokine systems approach demonstrates differences in innate and pro-inflammatory host responses between genetically distinct MERS-CoV isolates.

BACKGROUND: The recent emergence of a novel coronavirus in the Middle East (designated MERS-CoV) is a reminder of the zoonotic and pathogenic potential of emerging coronaviruses in humans. Clinical features of Middle East respiratory syndrome (MERS) include atypical pneumonia and progressive respiratory failure that is highly reminiscent of severe acute respiratory syndrome (SARS) caused by SARS-CoV. The host response is a key component of highly pathogenic respiratory virus infection. Here, we computationally analyzed gene expression changes in a human airway epithelial cell line infected with two genetically distinct MERS-CoV strains obtained from human patients, MERS-CoV SA 1 and MERS-CoV Eng 1. RESULTS: Using topological techniques, including persistence homology and filtered clustering, we performed a comparative transcriptional analysis of human Calu-3 cell host responses to the different MERS-CoV strains, with MERS-CoV Eng 1 inducing early kinetic changes, between 3 and 12 hours post infection, compared to MERS-CoV SA 1. Robust transcriptional changes distinguished the two MERS-CoV strains predominantly at the late time points. Combining statistical analysis of infection and cytokine-stimulated Calu-3 transcriptomics, we identified differential innate responses, including up-regulation of extracellular remodeling genes following MERS-CoV Eng 1 infection and differential pro-inflammatory responses. CONCLUSIONS: Through our genomics-based approach, we found topological differences in the kinetics and magnitude of the host response to MERS-CoV SA 1 and MERS-CoV Eng 1, with differential expression of innate immune and pro-inflammatory responsive genes as a result of IFN, TNF and IL-1α signaling. Predicted activation for STAT3 mediating gene expression relevant for epithelial cell-to-cell adherens and junction signaling in MERS-CoV Eng 1 infection suggest that these transcriptional differences may be the result of amino acid differences in viral proteins known to modulate innate immunity during MERS-CoV infection.

Activator-independent transcription of Snf1-dependent genes in mutants lacking histone tails.

Transcriptional regulation of Snf1-dependent genes occurs in part by histone-acetylation-dependent binding of the transcription factor Adr1. Analysis of previously published microarray data indicated unscheduled transcription of a large number of Snf1- and Adr1-dependent genes when either the histone H3 or H4 tail was deleted. Quantitative real-time PCR confirmed that the tails were important to preserve stringent transcriptional repression of Snf1-dependent genes when glucose was present. The absence of the tails allowed Adr1 and RNA Polymerase II to bind promoters in normally inhibitory conditions. The promoters escaped glucose repression to a limited extent and the weak constitutive ADH2 transcription induced by deletion of the histone tails was transcription factor- and Snf1-independent. These effects were apparently due to a permissive chromatin structure that allowed transcription in the absence of repression mediated by the histone tails. Deleting REG1, and thus activating Snf1 in the H3 tail mutant enhanced transcription in repressing conditions, indicating that Snf1 and the H3 tail influence transcription independently. Deleting REG1 in the histone H4 tail mutant appeared to be lethal, even in the absence of Snf1, suggesting that Reg1 and the H4 tail have redundant functions that are important for cell viability.

Host regulatory network response to infection with highly pathogenic H5N1 avian influenza virus.

During the last decade, more than half of humans infected with highly pathogenic avian influenza (HPAI) H5N1 viruses have died, yet virus-induced host signaling has yet to be clearly elucidated. Airway epithelia are known to produce inflammatory mediators that contribute to HPAI H5N1-mediated pathogenicity, but a comprehensive analysis of the host response in this cell type is lacking. Here, we leveraged a system approach to identify and statistically validate signaling subnetworks that define the dynamic transcriptional response of human bronchial epithelial cells after infection with influenza A/Vietnam/1203/2004 (H5N1, VN1203). Importantly, we validated a subset of transcripts from one subnetwork in both Calu-3 cells and mice. A more detailed examination of two subnetworks involved in the immune response and keratinization processes revealed potential novel mediators of HPAI H5N1 pathogenesis and host response signaling. Finally, we show how these results compare to those for a less virulent strain of influenza virus. Using emergent network properties, we provide fresh insight into the host response to HPAI H5N1 virus infection and identify novel avenues for perturbation studies and potential therapeutic interventions for fatal HPAI H5N1 disease.

Snf1-independent, glucose-resistant transcription of Adr1-dependent genes in a mediator mutant of Saccharomyces cerevisiae.

Glucose represses transcription of a network of co-regulated genes in Saccharomyces cerevisiae, ensuring that it is utilized before poorer carbon sources are metabolized. Adr1 is a glucose-regulated transcription factor whose promoter binding and activity require Snf1, the yeast homologue of the AMP-activated protein kinase in higher eukaryotes. In this study we found that a temperature-sensitive allele of MED14, a Mediator middle subunit that tethers the tail to the body, allowed a low level of Adr1-independent ADH2 expression that can be enhanced by Adr1 in a dose-dependent manner. A low level of TATA-independent ADH2 expression was observed in the med14-truncated strain and transcription of ADH2 and other Adr1-dependent genes occurred in the absence of Snf1 and chromatin remodeling coactivators. Loss of ADH2 promoter nucleosomes had occurred in the med14 strain in repressing conditions and did not require ADR1. A global analysis of transcription revealed that loss of Med14 function was associated with both up- and down- regulation of several groups of co-regulated genes, with ADR1-dependent genes being the most highly represented in the upregulated class. Expression of most genes was not significantly affected by the loss of Med14 function.

Strong TH1-biased CD4 T cell responses are associated with diminished SIV vaccine efficacy.

Activated CD4 T cells are a major target of HIV infection. Results from the STEP HIV vaccine trial highlighted a potential role for total activated CD4 T cells in promoting HIV acquisition. However, the influence of vaccine insert-specific CD4 T cell responses on HIV acquisition is not known. Here, using the data obtained from four macaque studies, we show that the DNA prime/modified vaccinia Ankara boost vaccine induced interferon γ (IFNγ+) CD4 T cells [T helper 1 (TH1) cells] rapidly migrate to multiple tissues including colon, cervix, and vaginal mucosa. These mucosal TH1 cells persisted at higher frequencies and expressed higher density of CCR5, a viral coreceptor, compared to cells in blood. After intravaginal or intrarectal simian immunodeficiency virus (SIV)/simian-human immunodeficiency virus (SHIV) challenges, strong vaccine protection was evident only in animals that had lower frequencies of vaccine-specific TH1 cells but not in animals that had higher frequencies of TH1 cells, despite comparable vaccine-induced humoral and CD8 T cell immunity in both groups. An RNA transcriptome signature in blood at 7 days after priming immunization from one study was associated with induction of fewer TH1-type CD4 cells and enhanced protection. These results demonstrate that high and persisting frequencies of HIV vaccine-induced TH1-biased CD4 T cells in the intestinal and genital mucosa can mitigate beneficial effects of protective antibodies and CD8 T cells, highlighting a critical role of priming immunization and vaccine adjuvants in modulating HIV vaccine efficacy.

Next-generation sequencing reveals a controlled immune response to Zaire Ebola virus challenge in cynomolgus macaques immunized with vesicular stomatitis virus expressing Zaire Ebola virus glycoprotein (VSVΔG/EBOVgp).

Vesicular stomatitis virus expressing Zaire Ebola virus (EBOV) glycoprotein (VSVΔG/EBOVgp) could be used as a vaccine to meet the 2014 Ebola virus outbreak. To characterize the host response to this vaccine, we used mRNA sequencing to analyze peripheral blood mononuclear cells (PBMCs) from cynomolgus macaques after VSVΔG/EBOVgp immunization and subsequent EBOV challenge. We found a controlled transcriptional response that transitioned to immune regulation as the EBOV was cleared. This observation supports the safety of the vaccine.

Multiple low-dose challenges in a rhesus macaque AIDS vaccine trial result in an evolving host response that affects protective outcome.

Using whole-blood transcriptional profiling, we investigated differences in the host response to vaccination and challenge in a rhesus macaque AIDS vaccine trial. Samples were collected from animals prior to and after vaccination with live, irradiated vaccine cells secreting the modified endoplasmic reticulum chaperone gp96-Ig loaded with simian immunodeficiency virus (SIV) peptides, either alone or in combination with a SIV-gp120 protein boost. Additional samples were collected following multiple low-dose rectal challenges with SIVmac251. Animals in the boosted group had a 73% reduced risk of infection. Surprisingly, few changes in gene expression were observed during the vaccination phase. Focusing on postchallenge comparisons, in particular for protected animals, we identified a host response signature of protection comprised of strong interferon signaling after the first challenge, which then largely abated after further challenges. We also identified a host response signature, comprised of early macrophage-mediated inflammatory responses, in animals with undetectable viral loads 5 days after the first challenge but with unusually high viral titers after subsequent challenges. Statistical analysis showed that prime-boost vaccination significantly lowered the probability of infection in a time-consistent manner throughout several challenges. Given that humoral responses in the prime-boost group were highly significant prechallenge correlates of protection, the strong innate signaling after the first challenge suggests that interferon signaling may enhance vaccine-induced antibody responses and is an important contributor to protection from infection during repeated low-dose exposure to SIV.

The draft genome sequence of the ferret (Mustela putorius furo) facilitates study of human respiratory disease.

The domestic ferret (Mustela putorius furo) is an important animal model for multiple human respiratory diseases. It is considered the 'gold standard' for modeling human influenza virus infection and transmission. Here we describe the 2.41 Gb draft genome assembly of the domestic ferret, constituting 2.28 Gb of sequence plus gaps. We annotated 19,910 protein-coding genes on this assembly using RNA-seq data from 21 ferret tissues. We characterized the ferret host response to two influenza virus infections by RNA-seq analysis of 42 ferret samples from influenza time-course data and showed distinct signatures in ferret trachea and lung tissues specific to 1918 or 2009 human pandemic influenza virus infections. Using microarray data from 16 ferret samples reflecting cystic fibrosis disease progression, we showed that transcriptional changes in the CFTR-knockout ferret lung reflect pathways of early disease that cannot be readily studied in human infants with cystic fibrosis disease.

A comprehensive collection of systems biology data characterizing the host response to viral infection.

The Systems Biology for Infectious Diseases Research program was established by the U.S. National Institute of Allergy and Infectious Diseases to investigate host-pathogen interactions at a systems level. This program generated 47 transcriptomic and proteomic datasets from 30 studies that investigate in vivo and in vitro host responses to viral infections. Human pathogens in the Orthomyxoviridae and Coronaviridae families, especially pandemic H1N1 and avian H5N1 influenza A viruses and severe acute respiratory syndrome coronavirus (SARS-CoV), were investigated. Study validation was demonstrated via experimental quality control measures and meta-analysis of independent experiments performed under similar conditions. Primary assay results are archived at the GEO and PeptideAtlas public repositories, while processed statistical results together with standardized metadata are publically available at the Influenza Research Database (www.fludb.org) and the Virus Pathogen Resource (www.viprbrc.org). By comparing data from mutant versus wild-type virus and host strains, RNA versus protein differential expression, and infection with genetically similar strains, these data can be used to further investigate genetic and physiological determinants of host responses to viral infection.

Drug repurposing: a better approach for infectious disease drug discovery?

The advent of publicly available databases containing system-wide phenotypic data of the host response to both drugs and pathogens, in conjunction with bioinformatics and computational methods now allows for in silico predictions of FDA-approved drugs as treatments against infection diseases. This systems biology approach captures the complexity of both the pathogen and drug host response in the form of expression patterns or molecular interaction networks without having to understand the underlying mechanisms of action. These drug repurposing techniques have been successful in identifying new drug candidates for several types of cancers and were recently used to identify potential therapeutics against influenza, the newly discovered Middle Eastern Respiratory Syndrome coronavirus and several parasitic diseases. These new approaches have the potential to significantly reduce both the time and cost for infectious diseases drug discovery.

Systems virology: host-directed approaches to viral pathogenesis and drug targeting.

High-throughput molecular profiling and computational biology are changing the face of virology, providing a new appreciation of the importance of the host in viral pathogenesis and offering unprecedented opportunities for better diagnostics, therapeutics and vaccines. Here, we provide a snapshot of the evolution of systems virology, from global gene expression profiling and signatures of disease outcome, to geometry-based computational methods that promise to yield novel therapeutic targets, personalized medicine and a deeper understanding of how viruses cause disease. To realize these goals, pipettes and Petri dishes need to join forces with the powers of mathematics and computational biology.

Mechanisms of severe acute respiratory syndrome coronavirus-induced acute lung injury.

UNLABELLED: Systems biology offers considerable promise in uncovering novel pathways by which viruses and other microbial pathogens interact with host signaling and expression networks to mediate disease severity. In this study, we have developed an unbiased modeling approach to identify new pathways and network connections mediating acute lung injury, using severe acute respiratory syndrome coronavirus (SARS-CoV) as a model pathogen. We utilized a time course of matched virologic, pathological, and transcriptomic data within a novel methodological framework that can detect pathway enrichment among key highly connected network genes. This unbiased approach produced a high-priority list of 4 genes in one pathway out of over 3,500 genes that were differentially expressed following SARS-CoV infection. With these data, we predicted that the urokinase and other wound repair pathways would regulate lethal versus sublethal disease following SARS-CoV infection in mice. We validated the importance of the urokinase pathway for SARS-CoV disease severity using genetically defined knockout mice, proteomic correlates of pathway activation, and pathological disease severity. The results of these studies demonstrate that a fine balance exists between host coagulation and fibrinolysin pathways regulating pathological disease outcomes, including diffuse alveolar damage and acute lung injury, following infection with highly pathogenic respiratory viruses, such as SARS-CoV. IMPORTANCE: Severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in 2002 and 2003, and infected patients developed an atypical pneumonia, acute lung injury (ALI), and acute respiratory distress syndrome (ARDS) leading to pulmonary fibrosis and death. We identified sets of differentially expressed genes that contribute to ALI and ARDS using lethal and sublethal SARS-CoV infection models. Mathematical prioritization of our gene sets identified the urokinase and extracellular matrix remodeling pathways as the most enriched pathways. By infecting Serpine1-knockout mice, we showed that the urokinase pathway had a significant effect on both lung pathology and overall SARS-CoV pathogenesis. These results demonstrate the effective use of unbiased modeling techniques for identification of high-priority host targets that regulate disease outcomes. Similar transcriptional signatures were noted in 1918 and 2009 H1N1 influenza virus-infected mice, suggesting a common, potentially treatable mechanism in development of virus-induced ALI.

Analysis of nucleosome positioning using a nucleosome-scanning assay.

The nucleosome-scanning assay (NuSA) couples isolation of mononucleosomal DNA after micrococcal nuclease (MNase) digestion with quantitative real-time PCR (qPCR) to map nucleosome positions in chromatin. It is a relatively simple, rapid procedure that can produce a high-resolution map of nucleosome location and occupancy and thus is suitable for analyzing individual promoters in great detail. The analysis can also quantify the protection of DNA sequences due to interaction with proteins other than nucleosomes and show how this protection varies when conditions change. When coupled with chromatin immunoprecipitation (ChIP), NuSA can identify histone variants and modifications associated with specific nucleosomes.

Next-generation sequencing reveals HIV-1-mediated suppression of T cell activation and RNA processing and regulation of noncoding RNA expression in a CD4+ T cell line.

UNLABELLED: Next-generation sequencing (NGS) enables the highly sensitive measurement of whole transcriptomes. We report the first application to our knowledge of this technology to the analysis of RNA from a CD4+ T cell line infected with intact HIV. We sequenced the total mRNA from infected cells and detected differences in the expression of both host and viral mRNA. Viral reads represented a large portion of the total mapped sequencing reads: approximately 20% at 12 h postinfection (hpi) and 40% at 24 hpi. We also detected a small but significant suppression of T cell activation-related genes at 12 hpi. This suppression persisted and expanded by 24 hpi, providing new possible markers of virus-induced T cell cytopathology. By 24 hpi, the expression of over 50% of detectable host loci was also altered, indicating widespread alteration of host processes, including RNA processing, splicing, and transport to an extent not previously reported. In addition, next-generation sequencing provided insights into alternative viral RNA splice events and the expression of noncoding RNAs, including microRNA host genes. IMPORTANCE: Recent advances in sequencing technology now allow the measurement of effectively all the RNA in a cell. This approach is especially useful for studying models of virus infection, as it allows the simultaneous measurement of both host and viral RNA. Using next-generation sequencing (NGS), we measured changes in total mRNA from a HIV-infected T cell line. To our knowledge, this is the first application of this technology to the investigation of HIV-host interactions involving intact HIV. We directly measured the amount of viral mRNA in infected cells and detected novel viral RNA splice variants and changes in the host expression of noncoding RNA species. We also detected small changes in T cell activation and other host processes during the early stages of viral replication that increased near the peak of viral replication, providing new candidate biomarkers of T cell death.

A systems biology approach to infectious disease research: innovating the pathogen-host research paradigm.

The twentieth century was marked by extraordinary advances in our understanding of microbes and infectious disease, but pandemics remain, food and waterborne illnesses are frequent, multidrug-resistant microbes are on the rise, and the needed drugs and vaccines have not been developed. The scientific approaches of the past-including the intense focus on individual genes and proteins typical of molecular biology-have not been sufficient to address these challenges. The first decade of the twenty-first century has seen remarkable innovations in technology and computational methods. These new tools provide nearly comprehensive views of complex biological systems and can provide a correspondingly deeper understanding of pathogen-host interactions. To take full advantage of these innovations, the National Institute of Allergy and Infectious Diseases recently initiated the Systems Biology Program for Infectious Disease Research. As participants of the Systems Biology Program, we think that the time is at hand to redefine the pathogen-host research paradigm.

Adr1 and Cat8 mediate coactivator recruitment and chromatin remodeling at glucose-regulated genes.

BACKGROUND: Adr1 and Cat8 co-regulate numerous glucose-repressed genes in S. cerevisiae, presenting a unique opportunity to explore their individual roles in coactivator recruitment, chromatin remodeling, and transcription. METHODOLOGY/PRINCIPAL FINDINGS: We determined the individual contributions of Cat8 and Adr1 on the expression of a cohort of glucose-repressed genes and found three broad categories: genes that need both activators for full derepression, genes that rely mostly on Cat8 and genes that require only Adr1. Through combined expression and recruitment data, along with analysis of chromatin remodeling at two of these genes, ADH2 and FBP1, we clarified how these activators achieve this wide range of co-regulation. We find that Adr1 and Cat8 are not intrinsically different in their abilities to recruit coactivators but rather, promoter context appears to dictate which activator is responsible for recruitment to specific genes. These promoter-specific contributions are also apparent in the chromatin remodeling that accompanies derepression: ADH2 requires both Adr1 and Cat8, whereas, at FBP1, significant remodeling occurs with Cat8 alone. Although over-expression of Adr1 can compensate for loss of Cat8 at many genes in terms of both activation and chromatin remodeling, this over-expression cannot complement all of the cat8Delta phenotypes. CONCLUSIONS/SIGNIFICANCE: Thus, at many of the glucose-repressed genes, Cat8 and Adr1 appear to have interchangeable roles and promoter architecture may dictate the roles of these activators.

The transcriptional coactivators SAGA, SWI/SNF, and mediator make distinct contributions to activation of glucose-repressed genes.

The paradigm of activation via ordered recruitment has evolved into a complicated picture as the influence of coactivators and chromatin structures on gene regulation becomes understood. We present here a comprehensive study of many elements of activation of ADH2 and FBP1, two glucose-regulated genes. We identify SWI/SNF as the major chromatin-remodeling complex at these genes, whereas SAGA (Spt-Ada-Gcn5-acetyltransferase complex) is required for stable recruitment of other coactivators. Mediator plays a crucial role in expression of both genes but does not affect chromatin remodeling. We found that Adr1 bound unaided by coactivators to ADH2, but Cat8 binding depended on coactivators at FBP1. Taken together, our results suggest that commonly regulated genes share many aspects of activation, but that gene-specific regulators or elements of promoter architecture may account for small differences in the mechanism of activation. Finally, we found that activator overexpression can compensate for the loss of SWI/SNF but not for the loss of SAGA.

A poised initiation complex is activated by SNF1.

Snf1, the yeast AMP kinase homolog, is essential for derepression of glucose-repressed genes that are activated by Adr1. Although required for Adr1 DNA binding, the precise role of Snf1 is unknown. Deletion of histone deacetylase genes allowed constitutive promoter binding of Adr1 and Cat8, another activator of glucose-repressed genes. In repressed conditions, at the Adr1-and Cat8-dependent ADH2 promoter, partial chromatin remodeling had occurred, and the activators recruited a partial preinitiation complex that included RNA polymerase II. Transcription did not occur, however, unless Snf1 was activated, suggesting a Snf1-dependent event that occurs after RNA polymerase II recruitment. Glucose regulation persisted because shifting to low glucose increased expression. Glucose repression could be completely relieved by combining the three elements of 1) chromatin perturbation by mutation of histone deacetylases, 2) activation of Snf1, and 3) the addition of an Adr1 mutant that by itself confers only weak constitutive activity.

PSGL-1 and mTOR regulate translation of ROCK-1 and physiological functions of macrophages.

Rho-associated kinases (ROCKs) are critical molecules involved in the physiological functions of macrophages, such as chemotaxis and phagocytosis. We demonstrate that macrophage adherence promotes rapid changes in physiological functions that depend on translational upregulation of preformed ROCK-1 mRNA, but not ROCK-2 mRNA. Before adherence, both ROCK mRNAs were present in the cytoplasm of macrophages, whereas ROCK proteins were undetectable. Macrophage adherence promoted signaling through P-selectin glycoprotein ligand-1 (PSGL-1)/Akt/mTOR that resulted in synthesis of ROCK-1, but not ROCK-2. Following synthesis, ROCK-1 was catalytically active. In addition, there was a rapamycin/sirolimus-sensitive enhanced loading of ribosomes on preformed ROCK-1 mRNAs. Inhibition of mTOR by rapamycin abolished ROCK-1 synthesis in macrophages resulting in an inhibition of chemotaxis and phagocytosis. Macrophages from PSGL-1-deficient mice recapitulated pharmacological inhibitor studies. These results indicate that receptor-mediated regulation at the level of translation is a component of a rapid set of mechanisms required to direct the macrophage phenotype upon adherence and suggest a mechanism for the immunosuppressive and anti-inflammatory effects of rapamycin/sirolimus.

The undertranslated transcriptome reveals widespread translational silencing by alternative 5' transcript leaders.

BACKGROUND: Translational efficiencies in Saccharomyces cerevisiae vary from transcript to transcript by approximately two orders of magnitude. Many of the poorly translated transcripts were found to respond to the appropriate external stimulus by recruiting ribosomes. Unexpectedly, a high frequency of these transcripts showed the appearance of altered 5' leaders that coincide with increased ribosome loading. RESULTS: Of the detectable transcripts in S. cerevisiae, 8% were found to be underloaded with ribosomes. Gene ontology categories of responses to stress or external stimuli were overrepresented in this population of transcripts. Seventeen poorly loaded transcripts involved in responses to pheromone, nitrogen starvation, and osmotic stress were selected for detailed study and were found to respond to the appropriate environmental signal with increased ribosome loading. Twelve of these regulated transcripts exhibited structural changes in their 5' transcript leaders in response to the environmental signal. In many of these the coding region remained intact, whereas regulated shortening of the 5' end truncated the open reading frame in others. Colinearity between the gene and transcript sequences eliminated regulated splicing as a mechanism for these alterations in structure. CONCLUSION: Frequent occurrence of coordinated changes in transcript structure and translation efficiency, in at least three different gene regulatory networks, suggests a widespread phenomenon. It is likely that many of these altered 5' leaders arose from changes in promoter usage. We speculate that production of translationally silenced transcripts may be one mechanism for allowing low-level transcription activity necessary for maintaining an open chromatin structure while not allowing inappropriate protein production.