Longitudinal peripheral blood transcriptional analysis of a patient with severe Ebola virus disease

Supervised learning approaches for predicting Ebola-Human Protein-Protein interactions

The goal of this research work is to predict protein-protein interactions (PPIs) between the Ebola virus and the host who is at risk of infection. Since there are very limited databases available on the Ebola virus; we have prepared a comprehensive database of all the PPIs between the Ebola virus and human proteins (EbolaInt). Our work focuses on the finding of some new protein-protein interactions between humans and the Ebola virus using some state- of-the-arts machine learning techniques. However, it is basically a two-class problem with a positive interacting dataset and a negative non-interacting dataset. These datasets contain various sequence-based human protein features such as structure of amino acid and conjoint triad and domain-related features. In this research, we have briefly discussed and used some well-known supervised learning approaches to predict PPIs between human proteins and Ebola virus proteins, including K-nearest neighbours (KNN), random forest (RF), support vector machine (SVM), and deep feed-forward multi-layer perceptron (DMLP) etc. We have validated our prediction results using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Our goal with this prediction is to compare all other models' accuracy, precision, recall, and f1-score for predicting these PPIs. In the result section, DMLP is giving the highest accuracy along with the prediction of 2655 potential human target proteins.

Machine learning-based analysis of Ebola virus' impact on gene expression in nonhuman primates

Introduction

This study introduces the Supervised Magnitude-Altitude Scoring (SMAS) methodology, a novel machine learning-based approach for analyzing gene expression data from non-human primates (NHPs) infected with Ebola virus (EBOV). By focusing on host-pathogen interactions, this research aims to enhance the understanding and identification of critical biomarkers for Ebola infection.

Methods

We utilized a comprehensive dataset of NanoString gene expression profiles from Ebola-infected NHPs. The SMAS system combines gene selection based on both statistical significance and expression changes. Employing linear classifiers such as logistic regression, the method facilitates precise differentiation between RT-qPCR positive and negative NHP samples.

Results

The application of SMAS led to the identification of IFI6 and IFI27 as key biomarkers, which demonstrated perfect predictive performance with 100% accuracy and optimal Area Under the Curve (AUC) metrics in classifying various stages of Ebola infection. Additionally, genes including MX1, OAS1, and ISG15 were significantly upregulated, underscoring their vital roles in the immune response to EBOV.

Discussion

Gene Ontology (GO) analysis further elucidated the involvement of these genes in critical biological processes and immune response pathways, reinforcing their significance in Ebola pathogenesis. Our findings highlight the efficacy of the SMAS methodology in revealing complex genetic interactions and response mechanisms, which are essential for advancing the development of diagnostic tools and therapeutic strategies.

Conclusion

This study provides valuable insights into EBOV pathogenesis, demonstrating the potential of SMAS to enhance the precision of diagnostics and interventions for Ebola and other viral infections.

Long-term Sudan Virus Ebola Survivors Maintain Multiple Antiviral Defense Mechanisms

BACKGROUND

The critical issues of sustained memory immunity following ebolavirus disease among long-term survivors are still unclear.

METHODS

Here, we examine virus-specific immune and inflammatory responses following in vitro challengd in 12 Sudan virus (SUDV) long-term survivors from Uganda's 2000-2001 Gulu outbreak, 15 years after recovery. Total RNA from isolated SUDV-stimulated and unstimulated peripheral blood mononuclear cells was extracted and analyzed. Matched serum samples were also collected to determine SUDV IgG levels and functionality.

RESULTS

We detected persistent humoral (58%, 7 of 12) and cellular (33%, 4 of 12) immune responses in SUDV long-term survivors and identified critical molecular mechanisms of innate and adaptive immunity. Gene expression in immune pathways, the interferon signaling system, antiviral defense response, and activation and regulation of T- and B-cell responses were observed. SUDV long-term survivors also maintained robust virus-specific IgG antibodies capable of polyfunctional responses, including neutralizing and innate Fc effector functions.

CONCLUSIONS

Data integration identified significant correlations among humoral and cellular immune responses and pinpointed a specific innate and adaptive gene expression signature associated with long-lasting immunity. This could help identify natural and vaccine correlates of protection against ebolavirus disease.

Peripheral immune responses to filoviruses in a reservoir versus spillover hosts reveal transcriptional correlates of disease

Several filoviruses, including Marburg virus (MARV), cause severe disease in humans and nonhuman primates (NHPs). However, the Egyptian rousette bat (ERB, Rousettus aegyptiacus), the only known MARV reservoir, shows no overt illness upon natural or experimental infection, which, like other bat hosts of zoonoses, is due to well-adapted, likely species-specific immune features. Despite advances in understanding reservoir immune responses to filoviruses, ERB peripheral blood responses to MARV and how they compare to those of diseased filovirus-infected spillover hosts remain ill-defined. We thus conducted a longitudinal analysis of ERB blood gene responses during acute MARV infection. These data were then contrasted with a compilation of published primate blood response studies to elucidate gene correlates of filovirus protection versus disease. Our work expands on previous findings in MARV-infected ERBs by supporting both host resistance and disease tolerance mechanisms, offers insight into the peripheral immunocellular repertoire during infection, and provides the most direct known cross-examination between reservoir and spillover hosts of the most prevalently-regulated response genes, pathways and activities associated with differences in filovirus pathogenesis and pathogenicity.

A metric for evaluating biological information in gene sets and its application to identify co-expressed gene clusters in PBMC

Recent technological advances have made the gathering of comprehensive gene expression datasets a commodity. This has shifted the limiting step of transcriptomic studies from the accumulation of data to their analyses and interpretation. The main problem in analyzing transcriptomics data is that the number of independent samples is typically much lower (<100) than the number of genes whose expression is quantified (typically >14,000). To address this, it would be desirable to reduce the gathered data’s dimensionality without losing information. Clustering genes into discrete modules is one of the most commonly used tools to accomplish this task. While there are multiple clustering approaches, there is a lack of informative metrics available to evaluate the resultant clusters’ biological quality. Here we present a metric that incorporates known ground truth gene sets to quantify gene clusters’ biological quality derived from standard clustering techniques. The GECO (Ground truth Evaluation of Clustering Outcomes) metric demonstrates that quantitative and repeatable scoring of gene clusters is not only possible but computationally lightweight and robust. Unlike current methods, it allows direct comparison between gene clusters generated by different clustering techniques. It also reveals that current cluster analysis techniques often underestimate the number of clusters that should be formed from a dataset, which leads to fewer clusters of lower quality. As a test case, we applied GECO combined with k-means clustering to derive an optimal set of co-expressed gene modules derived from PBMC, which we show to be superior to previously generated modules generated on whole-blood. Overall, GECO provides a rational metric to test and compare different clustering approaches to analyze high-dimensional transcriptomic data.

Next-generation sequencing has spurred the creation of many techniques that attempt to distill large datasets down to a more manageable size without losing valuable information, more simply referred to as dimensionality reduction. We have sought to contribute to this effort by focusing not directly on dimensionality reduction but on interpreting the results of the most common technique used for dimensionality reduction of sequencing data: gene clustering. While methods to generate gene clusters have been well explored, the evaluation of cluster quality has not, i.e., answering the question "Have we made biologically significant clusters?" We have developed a metric that can be used to answer this question. Our metric incorporates prior biological knowledge about the data to determine if the clustering process was optimal by looking at how genes are grouped in gene clusters and determine if they make sense biologically. Our metric can also be used to provide a discrete range of values that indicate how to generate clusters with the highest potential biological information content. This metric can be utilized by any -omics level study to generate study-specific gene clusters while reducing the time spent validating gene clusters and improving confidence in the resultant clusters.

Interferon-Induced HERC5 Inhibits Ebola Virus Particle Production and Is Antagonized by Ebola Glycoprotein

Survival following Ebola virus (EBOV) infection correlates with the ability to mount an early and robust interferon (IFN) response. The host IFN-induced proteins that contribute to controlling EBOV replication are not fully known. Among the top genes with the strongest early increases in expression after infection in vivo is IFN-induced HERC5. Using a transcription- and replication-competent VLP system, we showed that HERC5 inhibits EBOV virus-like particle (VLP) replication by depleting EBOV mRNAs. The HERC5 RCC1-like domain was necessary and sufficient for this inhibition and did not require zinc finger antiviral protein (ZAP). Moreover, we showed that EBOV (Zaire) glycoprotein (GP) but not Marburg virus GP antagonized HERC5 early during infection. Our data identify a novel ‘protagonist–antagonistic’ relationship between HERC5 and GP in the early stages of EBOV infection that could be exploited for the development of novel antiviral therapeutics.

Distinct transcriptional responses to fatal Ebola virus infection in cynomolgus and rhesus macaques suggest species-specific immune responses

Ebola virus (EBOV) is a negative single-stranded RNA virus within the Filoviridae family and the causative agent of Ebola virus disease (EVD). Nonhuman primates (NHPs), including cynomolgus and rhesus macaques, are considered the gold standard animal model to interrogate mechanisms of EBOV pathogenesis. However, despite significant genetic similarity (>90%), NHP species display different clinical presentation following EBOV infection, notably a ∼1-2 days delay in disease progression. Consequently, evaluation of therapeutics is generally conducted in rhesus macaques, whereas cynomolgus macaques are utilized to determine efficacy of preventative treatments, notably vaccines. This observation is in line with reported differences in disease severity and host responses between these two NHP following infection with simian varicella virus, influenza A and SARS-CoV-2. However, the molecular underpinnings of these differential outcomes following viral infections remain poorly defined. In this study, we compared published transcriptional profiles obtained from cynomolgus and rhesus macaques infected with the EBOV-Makona Guinea C07 using bivariate and regression analyses to elucidate differences in host responses. We report the presence of a shared core of differentially expressed genes (DEGs) reflecting EVD pathology, including aberrant inflammation, lymphopenia, and coagulopathy. However, the magnitudes of change differed between the two macaque species. These findings suggest that the differential clinical presentation of EVD in these two species is mediated by altered transcriptional responses.

Frontline Science: CD40 signaling restricts RNA virus replication in Mϕs, leading to rapid innate immune control of acute virus infection

Many acute viral infections target tissue Mϕs, yet the mechanisms of Mϕ-mediated control of viruses are poorly understood. Here, we report that CD40 expressed by peritoneal Mϕs restricts early infection of a broad range of RNA viruses. Loss of CD40 expression enhanced virus replication as early as 12-24 h of infection and, conversely, stimulation of CD40 signaling with an agonistic Ab blocked infection. With peritoneal cell populations infected with the filovirus, wild-type (WT) Ebola virus (EBOV), or a BSL2 model virus, recombinant vesicular stomatitis virus encoding Ebola virus glycoprotein (rVSV/EBOV GP), we examined the mechanism conferring protection. Here, we demonstrate that restricted virus replication in Mϕs required CD154/CD40 interactions that stimulated IL-12 production through TRAF6-dependent signaling. In turn, IL-12 production resulted in IFN-γ production, which induced proinflammatory polarization of Mϕs, protecting the cells from infection. These CD40-dependent events protected mice against virus challenge. CD40-/- mice were exquisitely sensitive to intraperitoneal challenge with a dose of rVSV/EBOV GP that was sublethal to CD40+/+ mice, exhibiting viremia within 12 h of infection and rapidly succumbing to infection. This study identifies a previously unappreciated role for Mϕ-intrinsic CD40 signaling in controlling acute virus infection.

Long-lasting severe immune dysfunction in Ebola virus disease survivors

Long-term follow up studies from Ebola virus disease (EVD) survivors (EBOV_S) are lacking. Here, we evaluate immune and gene expression profiles in 35 Guinean EBOV_S from the last West African outbreak, a median of 23 months (IQR [18–25]) after discharge from treatment center. Compared with healthy donors, EBOV_S exhibit increases of blood markers of inflammation, intestinal tissue damage, T cell and B cell activation and a depletion of circulating dendritic cells. All survivors have EBOV-specific IgG antibodies and robust and polyfunctional EBOV-specific memory T-cell responses. Deep sequencing of the genes expressed in blood reveals an enrichment in ‘inflammation’ and ‘antiviral’ pathways. Integrated analyses identify specific immune markers associated with the persistence of clinical symptoms. This study identifies a set of biological and genetic markers that could be used to define a signature of “chronic Ebola virus disease (CEVD)”.

Patients who have recovered from Ebola virus can have ongoing health problems. Here, the authors show that 35 Guinean survivors of the last West African Ebola epidemic have a chronic disease with high inflammatory cytokine expression and other markers of immune activation as well as evidence of intestinal tissue damage nearly two years after their release from hospital.

Longitudinal Human Antibody Repertoire against Complete Viral Proteome from Ebola Virus Survivor Reveals Protective Sites for Vaccine Design

Evolution of antibody repertoire against the Ebola virus (EBOV) proteome was characterized in an acutely infected patient receiving supportive care alone to elucidate virus-host interactions over time. Differential kinetics are observed for IgM-IgG-IgA epitope diversity, antibody binding, and affinity maturation to EBOV proteins. During acute illness, antibodies predominate to VP40 and glycoprotein (GP). At day 13 of clinical illness, a marked increase in antibody titers to most EBOV proteins and affinity maturation to GP is associated with rapid decline in viral replication and illness severity. At one year, despite undetectable virus, a diverse IgM repertoire against VP40 and GP epitopes is observed suggesting occult viral persistence. Rabbit immunization experiments identify key immunodominant sites of GP, while challenge studies in mice found these epitopes induce EBOV-neutralizing antibodies and protect against lethal EBOV challenge. This study reveals markers of viral persistence and provides promising approaches for development and evaluation of vaccines and therapeutics.

Modeling Challenges of Ebola Virus-Host Dynamics during Infection and Treatment

Mathematical modeling of Ebola virus (EBOV)-host dynamics during infection and treatment in vivo is in its infancy due to few studies with frequent viral kinetic data, lack of approved antiviral therapies, and limited insight into the timing of EBOV infection of cells and tissues throughout the body. Current in-host mathematical models simplify EBOV infection by assuming a single homogeneous compartment of infection. In particular, a recent modeling study assumed the liver as the largest solid organ targeted by EBOV infection and predicted that nearly all cells become refractory to infection within seven days of initial infection without antiviral treatment. We compared our observations of EBOV kinetics in multiple anatomic compartments and hepatocellular injury in a critically ill patient with Ebola virus disease (EVD) with this model's predictions. We also explored the model's predictions, with and without antiviral therapy, by recapitulating the model using published inputs and assumptions. Our findings highlight the challenges of modeling EBOV-host dynamics and therapeutic efficacy and emphasize the need for iterative interdisciplinary efforts to refine mathematical models that might advance understanding of EVD pathogenesis and treatment.

Cytokine Effects on the Entry of Filovirus Envelope Pseudotyped Virus-Like Particles into Primary Human Macrophages

Macrophages are one of the first and also a major site of filovirus replication and, in addition, are a source of multiple cytokines, presumed to play a critical role in the pathogenesis of the viral infection. Some of these cytokines are known to induce macrophage phenotypic changes in vitro, but how macrophage polarization may affect the cell susceptibility to filovirus entry remains largely unstudied. We generated different macrophage subsets using cytokine pre-treatment and subsequently tested their ability to fuse with beta-lactamase containing virus-like particles (VLP), pseudotyped with the surface glycoprotein of Ebola virus (EBOV) or the glycoproteins of other clinically relevant filovirus species. We found that pre-incubation of primary human monocyte-derived macrophages (MDM) with interleukin-10 (IL-10) significantly enhanced filovirus entry into cells obtained from multiple healthy donors, and the IL-10 effect was preserved in the presence of pro-inflammatory cytokines found to be elevated during EBOV disease. In contrast, fusion of IL-10-treated macrophages with influenza hemagglutinin/neuraminidase pseudotyped VLPs was unchanged or slightly reduced. Importantly, our in vitro data showing enhanced virus entry are consistent with the correlation established between elevated serum IL-10 and increased mortality in filovirus infected patients and also reveal a novel mechanism that may account for the IL-10-mediated increase in filovirus pathogenicity.

Ebola Immunity: Gaining a Winning Position in Lightning Chess

Zaire ebolavirus (EBOV), one of five species in the genus Ebolavirus, is the causative agent of the hemorrhagic fever disease epidemic that claimed more than 11,000 lives from 2014 to 2016 in West Africa. The combination of EBOV's ability to disseminate broadly and rapidly within the host and its high pathogenicity pose unique challenges to the human immune system postinfection. Potential transmission from apparently healthy EBOV survivors reported in the recent epidemic raises questions about EBOV persistence and immune surveillance mechanisms. Clinical, virological, and immunological data collected since the West Africa epidemic have greatly enhanced our knowledge of host-virus interactions. However, critical knowledge gaps remain in our understanding of what is necessary for an effective host immune response for protection against, or for clearance of, EBOV infection. This review provides an overview of immune responses against EBOV and discusses those associated with the success or failure to control EBOV infection.

Inflammatory and Humoral Immune Response during Ebola Virus Infection in Survivor and Fatal Cases Occurred in Sierra Leone during the 2014⁻2016 Outbreak in West Africa

Ebola virus (EBOV) infection is characterized by an excessive inflammatory response, a loss of lymphocytes and a general paralysis of the immune system, however pathophysiological mechanisms are not fully understood. In a cohort of 23 fatal and 21 survivors of ebola virus disease (EVD) cases admitted to the Emergency Ebola-Treatment-Center in Goderich (Freetown, Sierra Leone) during the 2014 to 2016 EBOV epidemic in Western Africa, we analyzed the pathway-focused gene expression profile of secreted proteins involved in the immune response and the levels of specific anti-EBOV IgM and IgG from the time of admission till discharge or death. We observed a dysregulated inflammatory response in fatal patients as compared to survivors, mainly consisting of the upregulation of inflammatory mediators, whose extent directly correlated with viremia levels. The upregulation persisted and intensified during the late phase of infection. Relevant differences were also found in humoral immunity, as an earlier and more robust EBOV antibody response was observed in survivor patients.

A conserved transcriptional response to intranasal Ebola virus exposure in nonhuman primates prior to onset of fever

Ebola virus disease (EVD), caused by Ebola virus (EBOV), is a severe illness characterized by case fatality rates of up to 90%. The sporadic nature of outbreaks in resource-limited areas has hindered the ability to characterize the pathogenesis of EVD at all stages of infection but particularly early host responses. Pathogenesis is often studied in nonhuman primate (NHP) models of disease that replicate major aspects of human EVD. Typically, NHP models use a large infectious dose, are carried out through intramuscular or aerosol exposure, and have a fairly uniform disease course. By contrast, we report our analysis of the host response to EBOV after intranasal exposure. Twelve cynomolgus macaques were infected with 100 plaque-forming units of EBOV/Makona through intranasal exposure and presented with varying times to onset of EVD. We used RNA sequencing and a newly developed NanoString CodeSet to monitor the host response via changes in RNA transcripts over time. When individual animal gene expression data were phased based on the onset of sustained fever, the first clinical sign of severe disease, mathematical models indicated that interferon-stimulated genes appeared as early as 4 days before fever onset. This demonstrates that lethal EVD has a uniform and predictable response to infection regardless of time to onset. Furthermore, expression of a subset of genes could predict disease development before other host-based indications of infection such as fever.

Multi-platform 'Omics Analysis of Human Ebola Virus Disease Pathogenesis

The pathogenesis of human Ebola virus disease (EVD) is complex. EVD is characterized by high levels of virus replication and dissemination, dysregulated immune responses, extensive virus- and host-mediated tissue damage, and disordered coagulation. To clarify how host responses contribute to EVD pathophysiology, we performed multi-platform 'omics analysis of peripheral blood mononuclear cells and plasma from EVD patients. Our results indicate that EVD molecular signatures overlap with those of sepsis, imply that pancreatic enzymes contribute to tissue damage in fatal EVD, and suggest that Ebola virus infection may induce aberrant neutrophils whose activity could explain hallmarks of fatal EVD. Moreover, integrated biomarker prediction identified putative biomarkers from different data platforms that differentiated survivors and fatalities early after infection. This work reveals insight into EVD pathogenesis, suggests an effective approach for biomarker identification, and provides an important community resource for further analysis of human EVD severity.

Comparison of Transcriptomic Platforms for Analysis of Whole Blood from Ebola-Infected Cynomolgus Macaques

Ebola virus disease (EVD) is a serious illness with mortality rates of 20-90% in various outbreaks. EVD is characterized by robust virus replication and strong host inflammatory response. Analyzing host immune responses has increasingly involved multimodal approaches including transcriptomics to profile gene expression. We studied cynomolgus macaques exposed to Ebola virus Makona via different routes with the intent of comparing RNA-Seq to a NanoString nCounter codeset targeting 769 non-human primate (NHP) genes. RNA-Seq analysis of serial blood samples showed different routes led to the same overall transcriptional response seen in previously reported EBOV-exposed NHP studies. Both platforms displayed a strong correlation in gene expression patterns, including a strong induction of innate immune response genes at early times post-exposure, and neutrophil-associated genes at later time points. A 41-gene classifier was tested in both platforms for ability to cluster samples by infection status. Both NanoString and RNA-Seq could be used to predict relative abundances of circulating immune cell populations that matched traditional hematology. This demonstrates the complementarity of RNA-Seq and NanoString. Moreover, the development of an NHP-specific NanoString codeset should augment studies of filoviruses and other high containment infectious diseases without the infrastructure requirements of RNA-Seq technology.

Transcriptome Analysis of Circulating Immune Cell Subsets Highlight the Role of Monocytes in Zaire Ebola Virus Makona Pathogenesis

Existing models of Ebola virus disease (EVD) suggest antigen-presenting cells are initial targets of Zaire ebolavirus (ZEBOV). In vitro studies have shown that ZEBOV infection of monocytes and macrophages results in the production of inflammatory mediators, which may cause lymphocyte apoptosis. However, these findings have not been corroborated by in vivo studies. In this study, we report the first longitudinal analysis of transcriptional changes in purified monocytes, T-cells, and B-cells isolated from cynomolgus macaques following infection with ZEBOV-Makona. Our data reveal monocytes as one of the major immune cell subsets that supports ZEBOV replication in vivo. In addition, we report a marked increase in the transcription of genes involved in inflammation, coagulation, and vascular disease within monocytes, suggesting that monocytes contribute to EVD manifestations. Further, genes important for antigen presentation and regulation of immunity were downregulated, potentially subverting development of adaptive immunity. In contrast, lymphocytes, which do not support ZEBOV replication, showed transcriptional changes limited to a small number of interferon-stimulated genes (ISGs) and a failure to upregulate genes associated with an antiviral effector immune response. Collectively, these data suggest that ZEBOV-infected monocytes play a significant role in ZEBOV-Makona pathogenesis and strategies to suppress virus replication or modify innate responses to infection in these cells should be a priority for therapeutic intervention.

Evidence of Ebola Virus Replication and High Concentration in Semen of a Patient During Recovery

In one patient over time, we found that concentration of Ebola virus RNA in semen during recovery is remarkably higher than blood at peak illness. Virus in semen is replication-competent with no change in viral genome over time. Presence of sense RNA suggests replication in cells present in semen.

Host Transcriptional Response to Ebola Virus Infection

Ebola virus disease (EVD) is a serious illness that causes severe disease in humans and non-human primates (NHPs) and has mortality rates up to 90%. EVD is caused by the Ebolavirus and currently there are no licensed therapeutics or vaccines to treat EVD. Due to its high mortality rates and potential as a bioterrorist weapon, a better understanding of the disease is of high priority. Multiparametric analysis techniques allow for a more complete understanding of a disease and the host response. Analysis of RNA species present in a sample can lead to a greater understanding of activation or suppression of different states of the immune response. Transcriptomic analyses such as microarrays and RNA-Sequencing (RNA-Seq) have been important tools to better understand the global gene expression response to EVD. In this review, we outline the current knowledge gained by transcriptomic analysis of EVD.