Computationally restoring the potency of a clinical antibody against Omicron.

The COVID-19 pandemic underscored the promise of monoclonal antibody-based prophylactic and therapeutic drugs1-3 and revealed how quickly viral escape can curtail effective options4,5. When the SARS-CoV-2 Omicron variant emerged in 2021, many antibody drug products lost potency, including Evusheld and its constituent, cilgavimab4-6. Cilgavimab, like its progenitor COV2-2130, is a class 3 antibody that is compatible with other antibodies in combination4 and is challenging to replace with existing approaches. Rapidly modifying such high-value antibodies to restore efficacy against emerging variants is a compelling mitigation strategy. We sought to redesign and renew the efficacy of COV2-2130 against Omicron BA.1 and BA.1.1 strains while maintaining efficacy against the dominant Delta variant. Here we show that our computationally redesigned antibody, 2130-1-0114-112, achieves this objective, simultaneously increases neutralization potency against Delta and subsequent variants of concern, and provides protection in vivo against the strains tested: WA1/2020, BA.1.1 and BA.5. Deep mutational scanning of tens of thousands of pseudovirus variants reveals that 2130-1-0114-112 improves broad potency without increasing escape liabilities. Our results suggest that computational approaches can optimize an antibody to target multiple escape variants, while simultaneously enriching potency. Our computational approach does not require experimental iterations or pre-existing binding data, thus enabling rapid response strategies to address escape variants or lessen escape vulnerabilities.

A compendium of multi-omics data illuminating host responses to lethal human virus infections.

Human infections caused by viral pathogens trigger a complex gamut of host responses that limit disease, resolve infection, generate immunity, and contribute to severe disease or death. Here, we present experimental methods and multi-omics data capture approaches representing the global host response to infection generated from 45 individual experiments involving human viruses from the Orthomyxoviridae, Filoviridae, Flaviviridae, and Coronaviridae families. Analogous experimental designs were implemented across human or mouse host model systems, longitudinal samples were collected over defined time courses, and global multi-omics data (transcriptomics, proteomics, metabolomics, and lipidomics) were acquired by microarray, RNA sequencing, or mass spectrometry analyses. For comparison, we have included transcriptomics datasets from cells treated with type I and type II human interferon. Raw multi-omics data and metadata were deposited in public repositories, and we provide a central location linking the raw data with experimental metadata and ready-to-use, quality-controlled, statistically processed multi-omics datasets not previously available in any public repository. This compendium of infection-induced host response data for reuse will be useful for those endeavouring to understand viral disease pathophysiology and network biology.

Genetic loci regulate Sarbecovirus pathogenesis: A comparison across mice and humans.

Coronavirus (CoV) cause considerable morbidity and mortality in humans and other mammals, as evidenced by the emergence of Severe Acute Respiratory CoV (SARS-CoV) in 2003, Middle East Respiratory CoV (MERS-CoV) in 2012, and SARS-CoV-2 in 2019. Although poorly characterized, natural genetic variation in human and other mammals modulate virus pathogenesis, as reflected by the spectrum of clinical outcomes ranging from asymptomatic infections to lethal disease. Using multiple human epidemic and zoonotic Sarbecoviruses, coupled with murine Collaborative Cross genetic reference populations, we identify several dozen quantitative trait loci that regulate SARS-like group-2B CoV pathogenesis and replication. Under a Chr4 QTL, we deleted a candidate interferon stimulated gene, Trim14 which resulted in enhanced SARS-CoV titers and clinical disease, suggesting an antiviral role during infection. Importantly, about 60 % of the murine QTL encode susceptibility genes identified as priority candidates from human genome-wide association studies (GWAS) studies after SARS-CoV-2 infection, suggesting that similar selective forces have targeted analogous genes and pathways to regulate Sarbecovirus disease across diverse mammalian hosts. These studies provide an experimental platform in rodents to investigate the molecular-genetic mechanisms by which potential cross mammalian susceptibility loci and genes regulate type-specific and cross-SARS-like group 2B CoV replication, immunity, and pathogenesis in rodent models. Our study also provides a paradigm for identifying susceptibility loci for other highly heterogeneous and virulent viruses that sporadically emerge from zoonotic reservoirs to plague human and animal populations.

Intestinal helminth infection impairs vaccine-induced T cell responses and protection against SARS-CoV-2.

Although vaccines have reduced COVID-19 disease burden, their efficacy in helminth infection endemic areas is not well characterized. We evaluated the impact of infection by Heligmosomoides polygyrus bakeri (Hpb), a murine intestinal hookworm, on the efficacy of an mRNA vaccine targeting the Wuhan-1 spike protein of SARS-CoV-2. Although immunization generated similar B cell responses in Hpb-infected and uninfected mice, polyfunctional CD4+ and CD8+ T cell responses were markedly reduced in Hpb-infected mice. Hpb-infected and mRNA vaccinated mice were protected against the ancestral SARS-CoV-2 strain WA1/2020, but control of lung infection was diminished against an Omicron variant compared to animals immunized without Hpb infection. Helminth mediated suppression of spike-specific CD8+ T cell responses occurred independently of STAT6 signaling, whereas blockade of IL-10 rescued vaccine-induced CD8+ T cell responses. In mice, intestinal helminth infection impairs vaccine induced T cell responses via an IL-10 pathway and compromises protection against antigenically shifted SARS-CoV-2 variants.

Domain-based mRNA vaccines encoding spike protein N-terminal and receptor binding domains confer protection against SARS-CoV-2.

With the success of messenger RNA (mRNA) vaccines against coronavirus disease 2019, strategies can now focus on improving vaccine potency, breadth, and stability. We designed and evaluated domain-based mRNA vaccines encoding the wild-type spike protein receptor binding domain (RBD) or N-terminal domain (NTD) alone or in combination. An NTD-RBD-linked candidate vaccine, mRNA-1283, showed improved antigen expression, antibody responses, and stability at refrigerated temperatures (2° to 8°C) compared with the clinically available mRNA-1273, which encodes the full-length spike protein. In BALB/c mice administered mRNA-1283 as a primary series, booster, or variant-specific booster, similar or greater immune responses from viral challenge were observed against wild-type, beta, delta, or omicron (BA.1) viruses compared with mRNA-1273-immunized mice, especially at lower vaccine dosages. K18-hACE2 mice immunized with mRNA-1283 or mRNA-1273 as a primary series demonstrated similar degrees of protection from challenge with SARS-CoV-2 Delta and Omicron variants at all vaccine dosages. These results support clinical assessment of mRNA-1283, which has now entered clinical trials (NCT05137236).

An Antiviral Role for TRIM14 in Ebola Virus Infection.

Ebola virus (EBOV) is a highly pathogenic virus that encodes 7 multifunctional structural proteins. Multiple host factors have been reported to interact with the EBOV proteins. Here, we found that tripartite motif-containing 14 (TRIM14), an interferon-stimulated gene that mediates cellular signaling pathways associated with type I interferon and inflammatory cytokine production, interacts with EBOV nucleoprotein to enhance interferon-ß (IFN-ß) and nuclear factor-κB (NF-κB) promotor activation. Moreover, TRIM14 overexpression reduced viral replication in an infectious but biologically contained EBOVΔVP30 system by approximately 10-fold without affecting viral protein expression. Furthermore, TRM14-deficient mice were more susceptible to mouse-adapted EBOV infection than wild-type mice. Our data suggest that TRIM14 is a host factor with anti-EBOV activity that limits EBOV pathogenesis.

Computationally restoring the potency of a clinical antibody against SARS-CoV-2 Omicron subvariants.

The COVID-19 pandemic underscored the promise of monoclonal antibody-based prophylactic and therapeutic drugs1-3, but also revealed how quickly viral escape can curtail effective options4,5. With the emergence of the SARS-CoV-2 Omicron variant in late 2021, many clinically used antibody drug products lost potency, including Evusheld™ and its constituent, cilgavimab4,6. Cilgavimab, like its progenitor COV2-2130, is a class 3 antibody that is compatible with other antibodies in combination4 and is challenging to replace with existing approaches. Rapidly modifying such high-value antibodies with a known clinical profile to restore efficacy against emerging variants is a compelling mitigation strategy. We sought to redesign COV2-2130 to rescue in vivo efficacy against Omicron BA.1 and BA.1.1 strains while maintaining efficacy against the contemporaneously dominant Delta variant. Here we show that our computationally redesigned antibody, 2130-1-0114-112, achieves this objective, simultaneously increases neutralization potency against Delta and many variants of concern that subsequently emerged, and provides protection in vivo against the strains tested, WA1/2020, BA.1.1, and BA.5. Deep mutational scanning of tens of thousands pseudovirus variants reveals 2130-1-0114-112 improves broad potency without incurring additional escape liabilities. Our results suggest that computational approaches can optimize an antibody to target multiple escape variants, while simultaneously enriching potency. Because our approach is computationally driven, not requiring experimental iterations or pre-existing binding data, it could enable rapid response strategies to address escape variants or pre-emptively mitigate escape vulnerabilities.

Bivalent SARS-CoV-2 mRNA vaccines increase breadth of neutralization and protect against the BA.5 Omicron variant in mice.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants in the Omicron lineage has resulted in diminished Coronavirus Disease 2019 (COVID-19) vaccine efficacy and persistent transmission. In this study, we evaluated the immunogenicity and protective efficacy of two, recently authorized, bivalent COVID-19 vaccines that contain two mRNAs encoding Wuhan-1 and either BA.1 (mRNA-1273.214) or BA.4/5 (mRNA-1273.222) spike proteins. As a primary two-dose immunization series in mice, both bivalent vaccines induced greater neutralizing antibody responses against Omicron variants than the parental, monovalent mRNA-1273 vaccine. When administered to mice as a booster at 7 months after the primary vaccination series with mRNA-1273, the bivalent vaccines induced broadly neutralizing antibody responses. Whereas most anti-Omicron receptor binding domain antibodies in serum induced by mRNA-1273, mRNA-1273.214 and mRNA-1273.222 boosters cross-reacted with the antecedent Wuhan-1 spike antigen, the mRNA-1273.214 and mRNA-1273.222 bivalent vaccine boosters also induced unique BA.1-specific and BA.4/5-specific responses, respectively. Although boosting with parental or bivalent mRNA vaccines substantially improved protection against BA.5 compared to mice receiving two vaccine doses, the levels of infection, inflammation and pathology in the lung were lowest in animals administered the bivalent mRNA vaccines. Thus, boosting with bivalent Omicron-based mRNA-1273.214 or mRNA-1273.222 vaccines enhances immunogenicity and confers protection in mice against a currently circulating SARS-CoV-2 strain.

Development of SARS-CoV-2 mRNA vaccines encoding spike N-terminal and receptor binding domains.

With the success of mRNA vaccines against coronavirus disease 2019 (COVID-19), strategies can now focus on improving vaccine potency, breadth, and stability. We present the design and preclinical evaluation of domain-based mRNA vaccines encoding the wild-type spike-protein receptor-binding (RBD) and/or N-terminal domains (NTD). An NTD-RBD linked candidate vaccine, mRNA-1283, showed improved antigen expression, antibody responses, and stability at refrigerated temperatures (2-8°C) compared with the clinically available mRNA-1273, which encodes the full-length spike protein. In mice administered mRNA-1283 as a primary series, booster, or variant-specific booster, similar or greater immune responses and protection from viral challenge were observed against wild-type, beta, delta, or omicron (BA. 1) compared with mRNA-1273 immunized mice, especially at lower vaccine dosages. These results support clinical assessment of mRNA-1283 ( NCT05137236 ). One Sentence Summary: A domain-based mRNA vaccine, mRNA-1283, is immunogenic and protective against SARS-CoV-2 and emerging variants in mice.

Bivalent SARS-CoV-2 mRNA vaccines increase breadth of neutralization and protect against the BA.5 Omicron variant.

The emergence of SARS-CoV-2 variants in the Omicron lineage with large numbers of substitutions in the spike protein that can evade antibody neutralization has resulted in diminished vaccine efficacy and persistent transmission. One strategy to broaden vaccine-induced immunity is to administer bivalent vaccines that encode for spike proteins from both historical and newly-emerged variant strains. Here, we evaluated the immunogenicity and protective efficacy of two bivalent vaccines that recently were authorized for use in Europe and the United States and contain two mRNAs encoding Wuhan-1 and either BA.1 (mRNA-1273.214) or BA.4/5 (mRNA-1273.222) spike proteins. As a primary immunization series in BALB/c mice, both bivalent vaccines induced broader neutralizing antibody responses than the constituent monovalent vaccines (mRNA-1273 [Wuhan-1], mRNA-1273.529 [BA.1], and mRNA-1273-045 [BA.4/5]). When administered to K18-hACE2 transgenic mice as a booster at 7 months after the primary vaccination series with mRNA-1273, the bivalent vaccines induced greater breadth and magnitude of neutralizing antibodies compared to an mRNA-1273 booster. Moreover, the response in bivalent vaccine-boosted mice was associated with increased protection against BA.5 infection and inflammation in the lung. Thus, boosting with bivalent Omicron-based mRNA-1273.214 or mRNA-1273.222 vaccines enhances immunogenicity and protection against currently circulating SARS-CoV-2 strains.

Kathryn V. Holmes: A Career of Contributions to the Coronavirus Field.

Over the past two years, scientific research has moved at an unprecedented rate in response to the COVID-19 pandemic. The rapid development of effective vaccines and therapeutics would not have been possible without extensive background knowledge on coronaviruses developed over decades by researchers, including Kathryn (Kay) Holmes. Kay's research team discovered the first coronavirus receptors for mouse hepatitis virus and human coronavirus 229E and contributed a wealth of information on coronaviral spike glycoproteins and receptor interactions that are critical determinants of host and tissue specificity. She collaborated with several research laboratories to contribute knowledge in additional areas, including coronaviral pathogenesis, epidemiology, and evolution. Throughout her career, Kay was an extremely dedicated and thoughtful mentor to numerous graduate students and post-doctoral fellows. This article provides a review of her contributions to the coronavirus field and her exemplary mentoring.

mRNA-1273 and Ad26.COV2.S vaccines protect against the B.1.621 variant of SARS-CoV-2.

BACKGROUND: Since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 2019, viral variants with greater transmissibility or immune-evasion properties have arisen, which could jeopardize recently deployed vaccine- and antibody-based countermeasures. METHODS: Here, we evaluated in mice and hamsters the efficacy of a pre-clinical version of the Moderna mRNA vaccine (mRNA-1273) and the Johnson & Johnson recombinant adenoviral-vectored vaccine (Ad26.COV2.S) against the B.1.621 (Mu) variant of SARS-CoV-2, which contains spike mutations T95I, Y144S, Y145N, R346K, E484K, N501Y, D614G, P681H, and D950N. FINDINGS: Immunization of 129S2 and K18-human ACE2 transgenic mice with the mRNA-1273 vaccine protected against weight loss, lung infection, and lung pathology after challenge with the B.1.621 or WA1/2020 N501Y/D614G SARS-CoV-2 strain. Similarly, immunization of 129S2 mice and Syrian hamsters with a high dose of Ad26.COV2.S reduced lung infection after B.1.621 virus challenge. CONCLUSIONS: Thus, immunity induced by the mRNA-1273 or Ad26.COV2.S vaccine can protect against the B.1.621 variant of SARS-CoV-2 in multiple animal models. FUNDING: This study was supported by the NIH (R01 AI157155 and U01 AI151810), NIAID Centers of Excellence for Influenza Research and Response [CEIRR] contracts 75N93021C00014 and 75N93021C00016, and the Collaborative Influenza Vaccine Innovation Centers [CIVIC] contract 75N93019C00051. It was also supported, in part, by the National Institutes of Allergy and Infectious Diseases Center for Research on Influenza Pathogenesis (HHSN272201400008C) and the Japan Program for Infectious Diseases Research and Infrastructure (JP21wm0125002) from the Japan Agency for Medical Research and Development (AMED).

Boosting with variant-matched or historical mRNA vaccines protects against Omicron infection in mice.

The large number of spike substitutions in Omicron lineage variants (BA.1, BA.1.1., and BA.2) could jeopardize the efficacy of SARS-CoV-2 vaccines. We evaluated in mice the protective efficacy of the Moderna mRNA-1273 vaccine against BA.1 before or after boosting. Whereas two doses of mRNA-1273 vaccine induced high levels of neutralizing antibodies against historical WA1/2020 strains, lower levels against BA.1 were associated with breakthrough infection and inflammation in the lungs. A primary vaccination series with mRNA-1273.529, an Omicron-matched vaccine, potently neutralized BA.1 but inhibited historical or other SARS-CoV-2 variants less effectively. However, boosting with either mRNA-1273 or mRNA-1273.529 vaccines increased neutralizing titers and protection against BA.1 and BA.2 infection. Nonetheless, the neutralizing antibody titers were higher, and lung viral burden and cytokines were slightly lower in mice boosted with mRNA-1273.529 and challenged with BA.1. Thus, boosting with mRNA-1273 or mRNA-1273.529 enhances protection against Omicron infection with limited differences in efficacy measured.

Boosting with Omicron-matched or historical mRNA vaccines increases neutralizing antibody responses and protection against B.1.1.529 infection in mice.

The B.1.1.529 Omicron variant jeopardizes vaccines designed with early pandemic spike antigens. Here, we evaluated in mice the protective activity of the Moderna mRNA-1273 vaccine against B.1.1.529 before or after boosting with preclinical mRNA-1273 or mRNA-1273.529, an Omicron-matched vaccine. Whereas two doses of mRNA-1273 vaccine induced high levels of serum neutralizing antibodies against historical WA1/2020 strains, levels were lower against B.1.1.529 and associated with infection and inflammation in the lung. A primary vaccination series with mRNA-1273.529 potently neutralized B.1.1.529 but showed limited inhibition of historical or other SARS-CoV-2 variants. However, boosting with mRNA-1273 or mRNA-1273.529 vaccines increased serum neutralizing titers and protection against B.1.1.529 infection. Nonetheless, the levels of inhibitory antibodies were higher, and viral burden and cytokines in the lung were slightly lower in mice given the Omicron-matched mRNA booster. Thus, in mice, boosting with mRNA-1273 or mRNA-1273.529 enhances protection against B.1.1.529 infection with limited differences in efficacy measured.

Protective activity of mRNA vaccines against ancestral and variant SARS-CoV-2 strains.

Although mRNA vaccines encoding the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) prevent COVID-19, the emergence of new viral variants jeopardizes their efficacy. Here, we assessed the immunogenicity and protective activity of historical (mRNA-1273, designed for Wuhan-1 spike protein) or modified (mRNA-1273.351, designed for B.1.351 spike protein) Moderna mRNA vaccines in 129S2 and K18-hACE2 mice. Mice were immunized with either high-dose or low-dose formulations of the mRNA vaccines, where low-dose vaccination modeled suboptimal immune responses. Immunization with formulations at either dose induced neutralizing antibodies in serum against ancestral SARS-CoV-2 WA1/2020 and several virus variants, although serum titers were lower against the B.1.617.2 (Delta) virus. Protection against weight loss and lung pathology was observed with all high-dose vaccines against all viruses. However, low-dose formulations of the vaccines, which produced lower magnitude antibody and T cell responses, showed breakthrough lung infections with B.1.617.2 and development of pneumonia in K18-hACE2 mice. Thus, in individuals with reduced immunity after mRNA vaccination, breakthrough infection and disease may occur with some SARS-CoV-2 variants.

Protective activity of mRNA vaccines against ancestral and variant SARS-CoV-2 strains.

Although mRNA vaccines prevent COVID-19, variants jeopardize their efficacy as immunity wanes. Here, we assessed the immunogenicity and protective activity of historical (mRNA-1273, designed for Wuhan-1 spike) or modified (mRNA-1273.351, designed for B.1.351 spike) preclinical Moderna mRNA vaccines in 129S2 and K18-hACE2 mice. Immunization with high or low dose formulations of mRNA vaccines induced neutralizing antibodies in serum against ancestral SARS-CoV-2 and several variants, although levels were lower particularly against the B.1.617.2 (Delta) virus. Protection against weight loss and lung pathology was observed with all high-dose vaccines against all viruses. Nonetheless, low-dose formulations of the vaccines, which produced lower magnitude antibody and T cell responses, and serve as a possible model for waning immunity, showed breakthrough lung infection and pneumonia with B.1.617.2. Thus, as levels of immunity induced by mRNA vaccines decline, breakthrough infection and disease likely will occur with some SARS-CoV-2 variants, suggesting a need for additional booster regimens.

Mesalamine Reduces Intestinal ACE2 Expression Without Modifying SARS-CoV-2 Infection or Disease Severity in Mice.

INTRODUCTION: Coronavirus Disease 2019 (COVID-19) is an ongoing public health crisis that has sickened or precipitated death in millions. The etiologic agent of COVID-19, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), infects the intestinal epithelium, and can induce GI symptoms similar to the human inflammatory bowel diseases (IBD). An international surveillance epidemiology study (SECURE-IBD) reported that the standardized mortality ratio trends higher in IBD patients (1.5-1.8) and that mesalamine/sulfasalazine therapy correlates with poor outcome. The goal of our study was to experimentally address the relationship between mesalamine and SARS-CoV-2 entry, replication, and/or pathogenesis. METHODS: Viral infection was performed with a chimeric vesicular stomatitis virus expressing SARS-CoV-2 spike protein and EGFP (VSV-SARS-CoV-2) and SARS-CoV-2 virus derived from an infectious cDNA clone of 2019n-CoV/USA_WA1/2020. Primary human ileal spheroids derived from healthy donors were grown as 3D spheroids or on 2D transwells. We assessed the effect of 10 mM mesalamine (Millipore Sigma) on viral RNA levels, as well as the expression of the SARS-CoV-2 receptor angiotensin II-converting enzyme 2 (ACE2), Transmembrane Serine Protease 2 (TMPRSS2), TMPRSS4, Cathepsin B (CTSB) and CTSL by qRT-PCR. 8-12 week old K18-ACE2 were treated orally with PBS or mesalamine at 200 mg/kg daily. Mice were inoculated intranasally with 1Ã-10 3 FFU of SARS-CoV-2. Mice were weighed daily and viral titers were determined 7 days post infection (dpi) by qRT-PCR. For the intestinal viral entry model, VSV-SARS-CoV-2 was injected into a ligated intestinal loop of anesthetized K18-ACE2 mice and tissues were harvested 6 hours post-infection. RESULTS: We found no change in viral RNA levels in human intestinal epithelial cells in response to mesalamine. Expression of ACE2 was reduced following mesalamine treatment in enteroids, while CTSL expression was increased. Mice receiving mesalamine lost weight at similar rates compared to mice receiving vehicle control. Mesalamine treatment did not change viral load in the lung, heart, or intestinal tissues harvested at 7 dpi. Pretreatment with mesalamine did not modulate intestinal entry of the chimeric VSV-SARS-CoV-2 in K18-ACE2 mice. CONCLUSIONS: Mesalamine did not alter viral entry, replication, or pathogenesis in vitro or in mouse models. Mesalamine treatment reduced expression of the viral receptor ACE2 while concurrently increasing CTSL expression in human ileum organoids.

Helminth-virus interactions: determinants of coinfection outcomes.

Viral infections are often studied in model mammalian organisms under specific pathogen-free conditions. However, in nature, coinfections are common, and infection with one organism can alter host susceptibility to infection with another. Helminth parasites share a long coevolutionary history with mammalian hosts and have shaped host physiology, metabolism, immunity, and the composition of the microbiome. Published studies suggest that helminth infection can either be beneficial or detrimental during viral infection. Here, we discuss coinfection studies in mouse models and use them to define key determinants that impact outcomes, including the type of antiviral immunity, the tissue tropism of both the helminth and the virus, and the timing of viral infection in relation to the helminth lifecycle. We also explore the current mechanistic understanding of how helminth-virus coinfection impacts host immunity and viral pathogenesis. While much attention has been placed on the impact of the gut bacterial microbiome on immunity to infection, we suggest that enteric helminths, as a part of the eukaryotic macrobiome, also represent an important modulator of disease pathogenesis and severity following virus infection.

Hypergraph models of biological networks to identify genes critical to pathogenic viral response.

BACKGROUND: Representing biological networks as graphs is a powerful approach to reveal underlying patterns, signatures, and critical components from high-throughput biomolecular data. However, graphs do not natively capture the multi-way relationships present among genes and proteins in biological systems. Hypergraphs are generalizations of graphs that naturally model multi-way relationships and have shown promise in modeling systems such as protein complexes and metabolic reactions. In this paper we seek to understand how hypergraphs can more faithfully identify, and potentially predict, important genes based on complex relationships inferred from genomic expression data sets. RESULTS: We compiled a novel data set of transcriptional host response to pathogenic viral infections and formulated relationships between genes as a hypergraph where hyperedges represent significantly perturbed genes, and vertices represent individual biological samples with specific experimental conditions. We find that hypergraph betweenness centrality is a superior method for identification of genes important to viral response when compared with graph centrality. CONCLUSIONS: Our results demonstrate the utility of using hypergraphs to represent complex biological systems and highlight central important responses in common to a variety of highly pathogenic viruses.