Shared and organism-specific host responses to childhood diarrheal diseases revealed by whole blood transcript profiling

Systems immunology of transcriptional responses to viral infection identifies conserved antiviral pathways across macaques and humans

Viral pandemics and epidemics pose a significant global threat. While macaque models of viral disease are routinely used, it remains unclear how conserved antiviral responses are between macaques and humans. Therefore, we conducted a cross-species analysis of transcriptomic data from over 6,088 blood samples from macaques and humans infected with one of 31 viruses. Our findings demonstrate that irrespective of primate or viral species, there are conserved antiviral responses that are consistent across infection phase (acute, chronic, or latent) and viral genome type (DNA or RNA viruses). Leveraging longitudinal data from experimental challenges, we identify virus-specific response kinetics such as host responses to Coronaviridae and Orthomyxoviridae infections peaking 1-3 days earlier than responses to Filoviridae and Arenaviridae viral infections. Our results underscore macaque studies as a powerful tool for understanding viral pathogenesis and immune responses that translate to humans, with implications for viral therapeutic development and pandemic preparedness.

Transcriptional profiling of zebrafish identifies host factors controlling susceptibility to Shigella flexneri

Shigella flexneri is a human-adapted pathovar of Escherichia coli that can invade the intestinal epithelium, causing inflammation and bacillary dysentery. Although an important human pathogen, the host response to S. flexneri has not been fully described. Zebrafish larvae represent a valuable model for studying human infections in vivo. Here, we use a Shigella-zebrafish infection model to generate mRNA expression profiles of host response to Shigella infection at the whole-animal level. Immune response-related processes dominate the signature of early Shigella infection (6 h post-infection). Consistent with its clearance from the host, the signature of late Shigella infection (24 h post-infection) is significantly changed, and only a small set of immune-related genes remain differentially expressed, including acod1 and gpr84. Using mutant lines generated by ENU, CRISPR mutagenesis and F0 crispants, we show that acod1- and gpr84-deficient larvae are more susceptible to Shigella infection. Together, these results highlight the power of zebrafish to model infection by bacterial pathogens and reveal the mRNA expression of the early (acutely infected) and late (clearing) host response to Shigella infection.

Projecting genetic associations through gene expression patterns highlights disease etiology and drug mechanisms

Genes act in concert with each other in specific contexts to perform their functions. Determining how these genes influence complex traits requires a mechanistic understanding of expression regulation across different conditions. It has been shown that this insight is critical for developing new therapies. Transcriptome-wide association studies have helped uncover the role of individual genes in disease-relevant mechanisms. However, modern models of the architecture of complex traits predict that gene-gene interactions play a crucial role in disease origin and progression. Here we introduce PhenoPLIER, a computational approach that maps gene-trait associations and pharmacological perturbation data into a common latent representation for a joint analysis. This representation is based on modules of genes with similar expression patterns across the same conditions. We observe that diseases are significantly associated with gene modules expressed in relevant cell types, and our approach is accurate in predicting known drug-disease pairs and inferring mechanisms of action. Furthermore, using a CRISPR screen to analyze lipid regulation, we find that functionally important players lack associations but are prioritized in trait-associated modules by PhenoPLIER. By incorporating groups of co-expressed genes, PhenoPLIER can contextualize genetic associations and reveal potential targets missed by single-gene strategies.

Diagnosis of childhood febrile illness using a multi-class blood RNA molecular signature

BACKGROUND

Appropriate treatment and management of children presenting with fever depend on accurate and timely diagnosis, but current diagnostic tests lack sensitivity and specificity and are frequently too slow to inform initial treatment. As an alternative to pathogen detection, host gene expression signatures in blood have shown promise in discriminating several infectious and inflammatory diseases in a dichotomous manner. However, differential diagnosis requires simultaneous consideration of multiple diseases. Here, we show that diverse infectious and inflammatory diseases can be discriminated by the expression levels of a single panel of genes in blood.

METHODS

A multi-class supervised machine-learning approach, incorporating clinical consequence of misdiagnosis as a "cost" weighting, was applied to a whole-blood transcriptomic microarray dataset, incorporating 12 publicly available datasets, including 1,212 children with 18 infectious or inflammatory diseases. The transcriptional panel identified was further validated in a new RNA sequencing dataset comprising 411 febrile children.

FINDINGS

We identified 161 transcripts that classified patients into 18 disease categories, reflecting individual causative pathogen and specific disease, as well as reliable prediction of broad classes comprising bacterial infection, viral infection, malaria, tuberculosis, or inflammatory disease. The transcriptional panel was validated in an independent cohort and benchmarked against existing dichotomous RNA signatures.

CONCLUSIONS

Our data suggest that classification of febrile illness can be achieved with a single blood sample and opens the way for a new approach for clinical diagnosis.

FUNDING

European Union's Seventh Framework no. 279185; Horizon2020 no. 668303 PERFORM; Wellcome Trust (206508/Z/17/Z); Medical Research Foundation (MRF-160-0008-ELP-KAFO-C0801); NIHR Imperial BRC.

Identification of Pediatric Bacterial Gastroenteritis From Blood Counts and Interviews Based on Machine Learning

INTRODUCTION

 Differentiating between bacterial and viral gastroenteritis is crucial in pediatric enteritis practice. Our objective was to use machine learning (ML) to identify acute gastroenteritis (AG) caused by bacteria based on blood cell counts and interview findings.

METHODS

 ML was performed using a decision tree classifier based on data from previously published papers. We included 164 children between one and 108 months diagnosed with gastroenteritis, with 112 having bacterial AG and 52 having viral AG as subjects and controls. Feature selection was performed using least absolute shrinkage and selection operator (LASSO), and the classifier's performance was evaluated by five-fold cross-validation. Additionally, we presented a tree diagram of the decision tree classifier as a flowchart for practical applications.

RESULTS

 The area under curve (AUC) was 0.80, indicating a moderate model. Three important features in this model were platelet-lymphocyte ratio, eosinophil count, and leukocyte count.

CONCLUSIONS

 In conclusion, this study demonstrates that bacterial AG can be estimated from blood cell counts with moderate accuracy. These findings may be valuable in narrowing down bacterial AG in children with gastrointestinal symptoms.

Characterizing the landscape of gene expression variance in humans

Gene expression variance has been linked to organismal function and fitness but remains a commonly neglected aspect of molecular research. As a result, we lack a comprehensive understanding of the patterns of transcriptional variance across genes, and how this variance is linked to context-specific gene regulation and gene function. Here, we use 57 large publicly available RNA-seq data sets to investigate the landscape of gene expression variance. These studies cover a wide range of tissues and allowed us to assess if there are consistently more or less variable genes across tissues and data sets and what mechanisms drive these patterns. We show that gene expression variance is broadly similar across tissues and studies, indicating that the pattern of transcriptional variance is consistent. We use this similarity to create both global and within-tissue rankings of variation, which we use to show that function, sequence variation, and gene regulatory signatures contribute to gene expression variance. Low-variance genes are associated with fundamental cell processes and have lower levels of genetic polymorphisms, have higher gene-gene connectivity, and tend to be associated with chromatin states associated with transcription. In contrast, high-variance genes are enriched for genes involved in immune response, environmentally responsive genes, immediate early genes, and are associated with higher levels of polymorphisms. These results show that the pattern of transcriptional variance is not noise. Instead, it is a consistent gene trait that seems to be functionally constrained in human populations. Furthermore, this commonly neglected aspect of molecular phenotypic variation harbors important information to understand complex traits and disease.

Discrimination of bacterial and viral infection using host-RNA signatures integrated in a lab-on-chip platform

The unmet clinical need for accurate point-of-care (POC) diagnostic tests able to discriminate bacterial from viral infection demands a solution that can be used both within healthcare settings and in the field, and that can also stem the tide of antimicrobial resistance. Our approach to solve this problem combine the use of host gene signatures with our Lab-on-a-Chip (LoC) technology enabling low-cost POC expression analysis to detect Infectious Disease. Transcriptomics have been extensively investigated as a potential tool to be implemented in the diagnosis of infectious disease. On the other hand, LoC technologies using ion-sensitive field-effect transistor (ISFET), in conjunction with isothermal chemistries, are offering a promising alternative to conventional amplification instruments, owing to their portable and affordable nature. Currently, the data analysis of ISFET arrays are restricted to established methods by averaging the output of every sensor to give a single time-series. This simple approach makes unrealistic assumptions, leading to insufficient performance for applications that require accurate quantification such as Host-Transcriptomics. In order to reliably quantify transcripts on our LoC platform enabling the classification of infectious disease on-chip, we propose a novel data-driven algorithm for extracting time-to-positive values from ISFET arrays. The algorithm proposed correctly outputs a time-to-positive for all the reactions, with a high correlation to RT-qLAMP (0.85, R2 = 0.98, p < 0.01), resulting in a classification accuracy of 100% (CI, 95-100%). This work aims to bridge the gap between translating assays from microarray analysis to ISFET arrays providing benefits on tackling infectious disease and diagnostic testing in hard-to-reach areas of the world.

A multi-tissue study of immune gene expression profiling highlights the key role of the nasal epithelium in COVID-19 severity

Coronavirus Disease-19 (COVID-19) symptoms range from mild to severe illness; the cause for this differential response to infection remains unknown. Unravelling the immune mechanisms acting at different levels of the colonization process might be key to understand these differences. We carried out a multi-tissue (nasal, buccal and blood; n = 156) gene expression analysis of immune-related genes from patients affected by different COVID-19 severities, and healthy controls through the nCounter technology. Mild and asymptomatic cases showed a powerful innate antiviral response in nasal epithelium, characterized by activation of interferon (IFN) pathway and downstream cascades, successfully controlling the infection at local level. In contrast, weak macrophage/monocyte driven innate antiviral response and lack of IFN signalling activity were present in severe cases. Consequently, oral mucosa from severe patients showed signals of viral activity, cell arresting and viral dissemination to the lower respiratory tract, which ultimately could explain the exacerbated innate immune response and impaired adaptative immune responses observed at systemic level. Results from saliva transcriptome suggest that the buccal cavity might play a key role in SARS-CoV-2 infection and dissemination in patients with worse prognosis. Co-expression network analysis adds further support to these findings, by detecting modules specifically correlated with severity involved in the abovementioned biological routes; this analysis also provides new candidate genes that might be tested as biomarkers in future studies. We also found tissue specific severity-related signatures mainly represented by genes involved in the innate immune system and cytokine/chemokine signalling. Local immune response could be key to determine the course of the systemic response and thus COVID-19 severity. Our findings provide a framework to investigate severity host gene biomarkers and pathways that might be relevant to diagnosis, prognosis, and therapy.

Evaluating the Immunopathogenesis of Diabetes After Acute Pancreatitis in the Diabetes RElated to Acute Pancreatitis and Its Mechanisms Study: From the Type 1 Diabetes in Acute Pancreatitis Consortium

ABSTRACT

The association between acute pancreatitis (AP) and diabetes mellitus (DM) has long been established, with the initial descriptions of AP patients presenting with DM after a bout of AP published in the 1940s and 50s. However, the potential mechanisms involved, particularly those components related to the immune system, have not been well defined. The Diabetes RElated to Acute pancreatitis and its Mechanisms (DREAM) study is a multicenter clinical study designed to understand the frequency and phenotype of DM developing after AP. This article describes one objective of the DREAM study: to determine the immunologic mechanisms of DM after AP, including the contribution of β-cell autoimmunity. This component of the study will assess the presence of islet autoimmunity, as well as the magnitude and kinetics of the innate and adaptive immune response at enrollment and during longitudinal follow-up after 1 or more episodes of AP. Finally, DREAM will evaluate the relationship between immune features, DM development, and pancreatitis etiology and severity.

Transcript host-RNA signatures to discriminate bacterial and viral infections in febrile children

Traditional laboratory markers, such as white blood cell count, C-reactive protein, and procalcitonin, failed to discriminate viral and bacterial infections in children. The lack of an accurate diagnostic test has a negative impact on child's care, limiting the ability of early diagnosis and appropriate management of children. This, on the one hand, may lead to delayed recognition of sepsis and severe bacterial infections, which still represent the leading causes of child morbidity and mortality. On the other hand, this may lead to overuse of empiric antibiotic therapies, particularly for specific subgroups of patients, such as infants younger than 90 days of life or neutropenic patients. This approach has an adverse effect on costs, antibiotic resistance, and pediatric microbiota. Transcript host-RNA signatures are a new tool used to differentiate viral from bacterial infections by analyzing the transcriptional biosignatures of RNA in host leukocytes. In this systematic review, we evaluate the efficacy and the possible application of this new diagnostic method in febrile children, along with challenges in its implementation. Our review support the growing evidence that the application of these new tools can improve the characterization of the spectrum of bacterial and viral infections and optimize the use of antibiotics in children. IMPACT: Transcript host RNA signatures may allow to better characterize the spectrum of viral, bacterial, and inflammatory illnesses in febrile children and can be used with traditional diagnostic methods to determine if and when to start antibiotic therapy. This is the first review on the use of transcript RNA signatures in febrile children to distinguish viral from bacterial infections. Our review identified a wide variability of target populations and gold standards used to define sepsis and SBIs, limiting the generalization of our findings.

Combating Childhood Infections in LMICs: evaluating the contribution of Big Data Big data, biomarkers and proteomics: informing childhood diarrhoeal disease management in Low- and Middle-Income Countries

Despite efforts to reduce the global burden of childhood diarrhoea, 50% of all cases globally occur in children under five years in Low–Income and Middle- Income Countries (LMICs) and knowledge gaps remain regarding the aetiological diagnosis, introduction of diarrhoeal vaccines, and the role of environmental enteric dysfunction and severe acute malnutrition. Biomarkers may assist in understanding disease processes, from diagnostics, to management of childhood diarrhoea and the sequelae to vaccine development. Proteomics has the potential to assist in the identification of new biomarkers to understand the processes in the development of childhood diarrhoea and to aid in developing new vaccines. Centralised repositories that enable mining of large data sets to better characterise risk factors, the proteome of both the patient and the different diarrhoeal pathogens, and the environment, could inform patient management and vaccine development, providing a systems biological approach to address the burden of childhood diarrhoea in LMICs.

Inflated false discovery rate due to volcano plots: problem and solutions

MOTIVATION

Volcano plots are used to select the most interesting discoveries when too many discoveries remain after application of Benjamini-Hochberg's procedure (BH). The volcano plot suggests a double filtering procedure that selects features with both small adjusted $P$-value and large estimated effect size. Despite its popularity, this type of selection overlooks the fact that BH does not guarantee error control over filtered subsets of discoveries. Therefore the selected subset of features may include an inflated number of false discoveries.

RESULTS

In this paper, we illustrate the substantially inflated type I error rate of volcano plot selection with simulation experiments and RNA-seq data. In particular, we show that the feature with the largest estimated effect is a very likely false positive result. Next, we investigate two alternative approaches for multiple testing with double filtering that do not inflate the false discovery rate. Our procedure is implemented in an interactive web application and is publicly available.

Identification of a Minimal 3-Transcript Signature to Differentiate Viral from Bacterial Infection from Best Genome-Wide Host RNA Biomarkers: A Multi-Cohort Analysis

The fight against the spread of antibiotic resistance is one of the most important challenges facing health systems worldwide. Given the limitations of current diagnostic methods, the development of fast and accurate tests for the diagnosis of viral and bacterial infections would improve patient management and treatment, as well as contribute to reducing antibiotic misuse in clinical settings. In this scenario, analysis of host transcriptomics constitutes a promising target to develop new diagnostic tests based on the host-specific response to infections. We carried out a multi-cohort meta-analysis of blood transcriptomic data available in public databases, including 11 different studies and 1209 samples from virus- (n = 695) and bacteria- (n = 514) infected patients. We applied a Parallel Regularized Regression Model Search (PReMS) on a set of previously reported genes that distinguished viral from bacterial infection to find a minimum gene expression bio-signature. This strategy allowed us to detect three genes, namely BAFT, ISG15 and DNMT1, that clearly differentiate groups of infection with high accuracy (training set: area under the curve (AUC) 0.86 (sensitivity: 0.81; specificity: 0.87); testing set: AUC 0.87 (sensitivity: 0.82; specificity: 0.86)). BAFT and ISG15 are involved in processes related to immune response, while DNMT1 is related to the preservation of methylation patterns, and its expression is modulated by pathogen infections. We successfully tested this three-transcript signature in the 11 independent studies, demonstrating its high performance under different scenarios. The main advantage of this three-gene signature is the low number of genes needed to differentiate both groups of patient categories.

Host Transcriptomic Response Following Administration of Rotavirus Vaccine in Infants' Mimics Wild Type Infection

Background

Rotavirus (RV) is an enteric pathogen that has devastating impact on childhood morbidity and mortality worldwide. The immunologic mechanism underlying the protection achieved after RV vaccination is not yet fully understood.

Methods

We compared the transcriptome of children affected by community-acquired RV infection and children immunized with a live attenuated RV vaccine (RotaTeq®).

Results

RV vaccination mimics the wild type infection causing similar changes in children's transcriptome, including transcripts associated with cell cycle, diarrhea, nausea, vomiting, intussusception, and abnormal morphology of midgut. A machine learning approach allowed to detect a combination of nine-transcripts that differentiates vaccinated from convalescent-naturally infected children (AUC: 90%; 95%CI: 70-100) and distinguishes between acute-infected and healthy control children (in both cases, AUC: 100%; 95%CI: 100-100). We identified a miRNA hsa-mir-149 that seems to play a role in the host defense against viral pathogens and may have an antiviral role.

Discussion

Our findings might shed further light in the understanding of RV infection, its functional link to intussusception causes, as well as guide development of antiviral treatments and safer and more effective vaccines. The nine-transcript signature may constitute a marker of vaccine protection and helps to differentiate vaccinated from naturally infected or susceptible children.

Mining the Unmapped Reads in Bovine RNA-Seq Data Reveals the Prevalence of Bovine Herpes Virus-6 in European Dairy Cows and the Associated Changes in Their Phenotype and Leucocyte Transcriptome

Microbial RNA is detectable in host samples by aligning unmapped reads from RNA sequencing against taxon reference sequences, generating a score proportional to the microbial load. An RNA-Seq data analysis showed that 83.5% of leukocyte samples from six dairy herds in different EU countries contained bovine herpes virus-6 (BoHV-6). Phenotypic data on milk production, metabolic function, and disease collected during their first 50 days in milk (DIM) were compared between cows with low (1–200 and n = 114) or high (201–1175 and n = 24) BoHV-6 scores. There were no differences in milk production parameters, but high score cows had numerically fewer incidences of clinical mastitis (4.2% vs. 12.2%) and uterine disease (54.5% vs. 62.7%). Their metabolic status was worse, based on measurements of IGF-1 and various metabolites in blood and milk. A comparison of the global leukocyte transcriptome between high and low BoHV-6 score cows at around 14 DIM yielded 485 differentially expressed genes (DEGs). The top pathway from Gene Ontology (GO) enrichment analysis was the immune system process. Down-regulated genes in the high BoHV-6 cows included those encoding proteins involved in viral detection (DDX6 and DDX58), interferon response, and E3 ubiquitin ligase activity. This suggested that BoHV-6 may largely evade viral detection and that it does not cause clinical disease in dairy cows.

RNA-Seq Data-Mining Allows the Discovery of Two Long Non-Coding RNA Biomarkers of Viral Infection in Humans

There is a growing interest in unraveling gene expression mechanisms leading to viral host invasion and infection progression. Current findings reveal that long non-coding RNAs (lncRNAs) are implicated in the regulation of the immune system by influencing gene expression through a wide range of mechanisms. By mining whole-transcriptome shotgun sequencing (RNA-seq) data using machine learning approaches, we detected two lncRNAs (ENSG00000254680 and ENSG00000273149) that are downregulated in a wide range of viral infections and different cell types, including blood monocluclear cells, umbilical vein endothelial cells, and dermal fibroblasts. The efficiency of these two lncRNAs was positively validated in different viral phenotypic scenarios. These two lncRNAs showed a strong downregulation in virus-infected patients when compared to healthy control transcriptomes, indicating that these biomarkers are promising targets for infection diagnosis. To the best of our knowledge, this is the very first study using host lncRNAs biomarkers for the diagnosis of human viral infections.

Ancestry patterns inferred from massive RNA-seq data

There is a growing body of evidence suggesting that patterns of gene expression vary within and between human populations. However, the impact of this variation in human diseases has been poorly explored, in part owing to the lack of a standardized protocol to estimate biogeographical ancestry from gene expression studies. Here we examine several studies that provide new solid evidence indicating that the ancestral background of individuals impacts gene expression patterns. Next, we test a procedure to infer genetic ancestry from RNA-seq data in 25 data sets where information on ethnicity was reported. Genome data of reference continental populations retrieved from The 1000 Genomes Project were used for comparisons. Remarkably, only eight out of 25 data sets passed FastQC default filters. We demonstrate that, for these eight population sets, the ancestral background of donors could be inferred very efficiently, even in data sets including samples with complex patterns of admixture (e.g., American-admixed populations). For most of the gene expression data sets of suboptimal quality, ancestral inference yielded odd patterns. The present study thus brings a cautionary note for gene expression studies highlighting the importance to control for the potential confounding effect of ancestral genetic background.

A 2-transcript host cell signature distinguishes viral from bacterial diarrhea and it is influenced by the severity of symptoms

Recently, a biomarker signature consisting of 2-transcript host RNAs was proposed for discriminating bacterial from viral infections in febrile children. We evaluated the performance of this signature in a different disease scenario, namely a cohort of Mexican children (n = 174) suffering from acute diarrhea of different infectious etiologies. We first examined the admixed background of the patients, indicating that most of them have a predominantly Native American genetic ancestry with a variable amount of European background (ranging from 0% to 57%). The results confirm that the RNA test can discriminate between viral and bacterial causes of infection (t-test; P-value = 6.94×10⁻¹¹; AUC = 80%; sensitivity: 68% [95% CI: 55%–79%]; specificity: 84% [95% CI: 78%–90%]), but the strength of the signal differs substantially depending on the causal pathogen, with the stronger signal being that of Shigella (P-value = 3.14 × 10⁻¹²; AUC = 89; sensitivity: 70% [95% CI: 57%–83%]; specificity: 100% [95% CI: 100%–100%]). The accuracy of this test improves significantly when excluding mild cases (P-value = 2.13 × 10⁻⁶; AUC = 85%; sensitivity: 79% [95% CI: 58%–95%]; specificity: 78% [95% CI: 65%–88%]). The results broaden the scope of previous studies by incorporating different pathogens, variable levels of disease severity, and different ancestral background of patients, and add confirmatory support to the clinical utility of these 2-transcript biomarkers.