Expanding the Immunology Toolbox: Embracing Public-Data Reuse and Crowdsourcing

Crowdsourcing temporal transcriptomic coronavirus host infection data: Resources, guide, and novel insights

The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) reawakened the need to rapidly understand the molecular etiologies, pandemic potential, and prospective treatments of infectious agents. The lack of existing data on SARS-CoV-2 hampered early attempts to treat severe forms of coronavirus disease-2019 (COVID-19) during the pandemic. This study coupled existing transcriptomic data from severe acute respiratory syndrome-related coronavirus 1 (SARS-CoV-1) lung infection animal studies with crowdsourcing statistical approaches to derive temporal meta-signatures of host responses during early viral accumulation and subsequent clearance stages. Unsupervised and supervised machine learning approaches identified top dysregulated genes and potential biomarkers (e.g. CXCL10, BEX2, and ADM). Temporal meta-signatures revealed distinct gene expression programs with biological implications to a series of host responses underlying sustained Cxcl10 expression and Stat signaling. Cell cycle switched from G1/G0 phase genes, early in infection, to a G2/M gene signature during late infection that correlated with the enrichment of DNA damage response and repair genes. The SARS-CoV-1 meta-signatures were shown to closely emulate human SARS-CoV-2 host responses from emerging RNAseq, single cell, and proteomics data with early monocyte-macrophage activation followed by lymphocyte proliferation. The circulatory hormone adrenomedullin was observed as maximally elevated in elderly patients who died from COVID-19. Stage-specific correlations to compounds with potential to treat COVID-19 and future coronavirus infections were in part validated by a subset of twenty-four that are in clinical trials to treat COVID-19. This study represents a roadmap to leverage existing data in the public domain to derive novel molecular and biological insights and potential treatments to emerging human pathogens.

Results and lessons learned from the sbv IMPROVER metagenomics diagnostics for inflammatory bowel disease challenge

A growing body of evidence links gut microbiota changes with inflammatory bowel disease (IBD), raising the potential benefit of exploiting metagenomics data for non-invasive IBD diagnostics. The sbv IMPROVER metagenomics diagnosis for inflammatory bowel disease challenge investigated computational metagenomics methods for discriminating IBD and nonIBD subjects. Participants in this challenge were given independent training and test metagenomics data from IBD and nonIBD subjects, which could be wither either raw read data (sub-challenge 1, SC1) or processed Taxonomy- and Function-based profiles (sub-challenge 2, SC2). A total of 81 anonymized submissions were received between September 2019 and March 2020. Most participants' predictions performed better than random predictions in classifying IBD versus nonIBD, Ulcerative Colitis (UC) versus nonIBD, and Crohn's Disease (CD) versus nonIBD. However, discrimination between UC and CD remains challenging, with the classification quality similar to the set of random predictions. We analyzed the class prediction accuracy, the metagenomics features by the teams, and computational methods used. These results will be openly shared with the scientific community to help advance IBD research and illustrate the application of a range of computational methodologies for effective metagenomic classification.

Transcriptomic meta-analysis reveals unannotated long non-coding RNAs related to the immune response in sheep

Long non-coding RNAs (lncRNAs) are involved in several biological processes, including the immune system response to pathogens and vaccines. The annotation and functional characterization of lncRNAs is more advanced in humans than in livestock species. Here, we take advantage of the increasing number of high-throughput functional experiments deposited in public databases in order to uniformly analyse, profile unannotated lncRNAs and integrate 422 ovine RNA-seq samples from the ovine immune system. We identified 12302 unannotated lncRNA genes with support from independent CAGE-seq and histone modification ChIP-seq assays. Unannotated lncRNAs showed low expression levels and sequence conservation across other mammal species. There were differences in expression levels depending on the genomic location-based lncRNA classification. Differential expression analyses between unstimulated and samples stimulated with pathogen infection or vaccination resulted in hundreds of lncRNAs with changed expression. Gene co-expression analyses revealed immune gene-enriched clusters associated with immune system activation and related to interferon signalling, antiviral response or endoplasmic reticulum stress. Besides, differential co-expression networks were constructed in order to find condition-specific relationships between coding genes and lncRNAs. Overall, using a diverse set of immune system samples and bioinformatic approaches we identify several ovine lncRNAs associated with the response to an external stimulus. These findings help in the improvement of the ovine lncRNA catalogue and provide sheep-specific evidence for the implication in the general immune response for several lncRNAs.

OMiCC: An expanded and enhanced platform for meta-analysis of public gene expression data

OMiCC (OMics Compendia Commons) is a biologist-friendly web platform that facilitates data reuse and integration. Users can search over 40,000 publicly available gene expression studies, annotate and curate samples, and perform meta-analysis. Since the initial publication, we have incorporated RNA-seq datasets, compendia sharing, RESTful API support, and an additional meta-analysis method based on random effects. Here, we provide a step-by-step guide for using OMiCC.

For complete details on the use and execution of this protocol, please refer to Shah et al. (2016).

•

OMiCC (OMics Compendia Commons) is a free web-based tool for gene expression data reuse

•

Search publicly available studies to perform sample group comparisons to explore a disease

•

In meta-analysis, multiple studies are combined to identify coherent signals

•

OMiCC supports crowd-sharing and users can share their own analyses with the community

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Identification of key genes and pathways involved in vitiligo development based on integrated analysis

Vitiligo is a chronic skin condition lack of melanocytes. However, researches on the aetiology and pathogenesis of vitiligo are still under debate. This study aimed to explore the key genes and pathways associated with occurrence and development of vitiligo.Weighted gene coexpression network analysis (WGCNA) was applied to reanalyze the gene expression dataset GSE65127 systematically. Functional enrichments of these modules were carried out at gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set variation analysis (GSVA), and gene set enrichment analysis (GSEA). Then, a map of regulatory network was delineated according to pivot analysis and drug prediction. In addition, hub genes and crucial pathways were validated by an independent dataset GSE75819. The expressions of hub genes in modules were also tested by quantitative real-time polymerase chain reaction (qRT-PCR).Eight coexpressed modules were identified by WGCNA based on 5794 differentially expressed genes of vitiligo. Three modules were found to be significantly correlated with Lesional, Peri-Lesional, and Non-Lesional, respectively. The persistent maladjusted genes included 269 upregulated genes and 82 downregulated genes. The enrichments showed module genes were implicated in immune response, p53 signaling pathway, etc. According to GSEA and GSVA, dysregulated pathways were activated incessantly from Non-Lesional to Peri-Lesional and then to Lesional, 4 of which were verified by an independent dataset GSE75819. Finally, 42 transcription factors and 228 drugs were spotted. Focusing on the persistent maladjusted genes, a map of regulatory network was delineated. Hub genes (CACTIN, DCTN1, GPR143, HADH, MRPL47, NKTR, NUF2) and transcription factors (ITGAV, SYK, PDPK1) were validated by an independent dataset GSE75819. In addition, hub genes (CACTIN, DCTN1, GPR143, MRPL47, NKTR) were also confirmed by qRT-PCR.The present study, at least, might provide an integrated and in-depth insight for exploring the underlying mechanism of vitiligo and predicting potential diagnostic biomarkers and therapeutic targets.

Immune and Metabolic Signatures of COVID-19 Revealed by Transcriptomics Data Reuse

The current pandemic of coronavirus disease 19 (COVID-19) has affected millions of individuals and caused thousands of deaths worldwide. The pathophysiology of the disease is complex and mostly unknown. Therefore, identifying the molecular mechanisms that promote progression of the disease is critical to overcome this pandemic. To address such issues, recent studies have reported transcriptomic profiles of cells, tissues and fluids from COVID-19 patients that mainly demonstrated activation of humoral immunity, dysregulated type I and III interferon expression, intense innate immune responses and inflammatory signaling. Here, we provide novel perspectives on the pathophysiology of COVID-19 using robust functional approaches to analyze public transcriptome datasets. In addition, we compared the transcriptional signature of COVID-19 patients with individuals infected with SARS-CoV-1 and Influenza A (IAV) viruses. We identified a core transcriptional signature induced by the respiratory viruses in peripheral leukocytes, whereas the absence of significant type I interferon/antiviral responses characterized SARS-CoV-2 infection. We also identified the higher expression of genes involved in metabolic pathways including heme biosynthesis, oxidative phosphorylation and tryptophan metabolism. A BTM-driven meta-analysis of bronchoalveolar lavage fluid (BALF) from COVID-19 patients showed significant enrichment for neutrophils and chemokines, which were also significant in data from lung tissue of one deceased COVID-19 patient. Importantly, our results indicate higher expression of genes related to oxidative phosphorylation both in peripheral mononuclear leukocytes and BALF, suggesting a critical role for mitochondrial activity during SARS-CoV-2 infection. Collectively, these data point for immunopathological features and targets that can be therapeutically exploited to control COVID-19.

Systems biology and big data in asthma and allergy: recent discoveries and emerging challenges

Asthma is a common condition caused by immune and respiratory dysfunction, and it is often linked to allergy. A systems perspective may prove helpful in unravelling the complexity of asthma and allergy. Our aim is to give an overview of systems biology approaches used in allergy and asthma research. Specifically, we describe recent "omic"-level findings, and examine how these findings have been systematically integrated to generate further insight.Current research suggests that allergy is driven by genetic and epigenetic factors, in concert with environmental factors such as microbiome and diet, leading to early-life disturbance in immunological development and disruption of balance within key immuno-inflammatory pathways. Variation in inherited susceptibility and exposures causes heterogeneity in manifestations of asthma and other allergic diseases. Machine learning approaches are being used to explore this heterogeneity, and to probe the pathophysiological patterns or "endotypes" that correlate with subphenotypes of asthma and allergy. Mathematical models are being built based on genomic, transcriptomic and proteomic data to predict or discriminate disease phenotypes, and to describe the biomolecular networks behind asthma.The use of systems biology in allergy and asthma research is rapidly growing, and has so far yielded fruitful results. However, the scale and multidisciplinary nature of this research means that it is accompanied by new challenges. Ultimately, it is hoped that systems medicine, with its integration of omics data into clinical practice, can pave the way to more precise, personalised and effective management of asthma.

Genomic Circuitry Underlying Immunological Response to Pediatric Acute Respiratory Infection

Acute respiratory tract viral infections (ARTIs) cause significant morbidity and mortality. CD8 T cells are fundamental to host responses, but transcriptional alterations underlying anti-viral mechanisms and links to clinical characteristics remain unclear. CD8 T cell transcriptional circuitry in acutely ill pediatric patients with influenza-like illness was distinct for different viral pathogens. Although changes included expected upregulation of interferon-stimulated genes (ISGs), transcriptional downregulation was prominent upon exposure to innate immune signals in early IFV infection. Network analysis linked changes to severity of infection, asthma, sex, and age. An influenza pediatric signature (IPS) distinguished acute influenza from other ARTIs and outperformed other influenza prediction gene lists. The IPS allowed a deeper investigation of the connection between transcriptional alterations and clinical characteristics of acute illness, including age-based differences in circuits connecting the STAT1/2 pathway to ISGs. A CD8 T cell-focused systems immunology approach in pediatrics identified age-based alterations in ARTI host response pathways.

Low expression of KCNN3 may affect drug resistance in ovarian cancer

Drug resistance is a principal contributor to the poor prognosis of ovarian cancer (OC). Therefore, identifying factors that affect drug resistance in OC is critical. In the present study, 51 OC specimens from lab collections were immunohistochemically tested, public data for 489 samples from The Cancer Genome Atlas cohort and 1,656 samples from the Kaplan‑Meier Plotter were downloaded, and data were retrieved from Oncomine. It was identified that the mRNA and protein expression of the potassium calcium‑activated channel subfamily N member 3 (KCNN3) was markedly lower in OC tissues compared with normal tissues, and in drug‑resistant OC tissues compared with sensitive OC tissues. Low KCNN3 expression consistently predicted shorter disease‑free and overall survival (OS). Specifically, low KCNN3 expression predicted shorter OS in 395 patients with low expression levels of mucin‑16. There was additional evidence that KCNN3 expression is mediated by microRNA‑892b. Furthermore, text mining and analyses of protein and gene interactions indicated that KCNN3 affects drug resistance. To the best of the authors' knowledge, this is the first report to associate KCNN3 with poor prognosis and drug resistance in OC. The present findings indicated that KCNN3 is a potential prognostic marker and therapeutic target for OC.

Key genes and functional coexpression modules involved in the pathogenesis of systemic lupus erythematosus

We performed a systematic review of genome-wide gene expression datasets to identify key genes and functional modules involved in the pathogenesis of systemic lupus erythematosus (SLE) at a systems level. Genome-wide gene expression datasets involving SLE patients were searched in Gene Expression Omnibus and ArrayExpress databases. Robust rank aggregation (RRA) analysis was used to integrate those public datasets and identify key genes associated with SLE. The weighted gene coexpression network analysis (WGCNA) was adapted to identify functional modules involved in SLE pathogenesis, and the gene ontology enrichment analysis was utilized to explore their functions. The aberrant expressions of several randomly selected key genes were further validated in SLE patients through quantitative real-time polymerase chain reaction. Fifteen genome-wide gene expression datasets were finally included, which involved a total of 1,778 SLE patients and 408 healthy controls. A large number of significantly upregulated or downregulated genes were identified through RRA analysis, and some of those genes were novel SLE gene signatures and their molecular roles in etiology of SLE remained vague. WGCNA further successfully identified six main functional modules involved in the pathogenesis of SLE. The most important functional module involved in SLE included 182 genes and mainly enriched in biological processes, including defense response to virus, interferon signaling pathway, and cytokine-mediated signaling pathway. This study identifies a number of key genes and functional coexpression modules involved in SLE, which provides deepening insights into the molecular mechanism of SLE at a systems level and also provides some promising therapeutic targets.

Reuse of public, genome-wide, murine eosinophil expression data for hypotheses development

The eosinophil (Eos) surface phenotype and activation state is altered after recruitment into tissues and after exposure to pro-inflammatory cytokines. In addition, distinct Eos functional subsets have been described, suggesting that tissue-specific responses for Eos contribute to organ homeostasis. Understanding the mechanisms by which Eos subsets achieve their tissue-specific identity is currently an unmet goal for the eosinophil research community. Publicly archived expression data can be used to answer original questions, test and generate new hypotheses, and serve as a launching point for experimental design. With these goals in mind, we investigated the effect of genetic background, culture methods, and tissue residency on murine Eos gene expression using publicly available, genome-wide expression data. Eos differentiated from cultures have a gene expression profile that is distinct from that of native homeostatic Eos; thus, researchers can repurpose published expression data to aid in selecting the appropriate culture method to study their gene of interest. In addition, we identified Eos lung- and gastrointestinal-specific transcriptomes, highlighting the profound effect of local tissue environment on gene expression in a terminally differentiated granulocyte even at homeostasis. Expanding the "toolbox" of Eos researchers to include public-data reuse can reduce redundancy, increase research efficiency, and lead to new biological insights.

Rebooting Human Immunology

Recent progress in both conceptual and technological approaches to human immunology have rejuvenated a field that has long been in the shadow of the inbred mouse model. This is a healthy development both for the clinical relevance of immunology and for the fact that it is a way to gain access to the wealth of phenomenology in the many human diseases that involve the immune system. This is where we are likely to discover new immunological mechanisms and principals, especially those involving genetic heterogeneity or environmental influences that are difficult to model effectively in inbred mice. We also suggest that there are likely to be novel immunological mechanisms in long-lived, less fecund mammals such as human beings since they must remain healthy far longer than short-lived rodents in order for the species to survive.

Computational drug repositioning for rare diseases in the era of precision medicine

There are tremendous unmet needs in drug development for rare diseases. Computational drug repositioning is a promising approach and has been successfully applied to the development of treatments for diseases. However, how to utilize this knowledge and effectively conduct and implement computational drug repositioning approaches for rare disease therapies is still an open issue. Here, we focus on the means of utilizing accumulated genomic data for accelerating and facilitating drug repositioning for rare diseases. First, we summarize the current genome landscape of rare diseases. Second, we propose several promising bioinformatics approaches and pipelines for computational drug repositioning for rare diseases. Finally, we discuss recent regulatory incentives and other enablers in rare disease drug development and outline the remaining challenges.

A generic method for improving the spatial interoperability of medical and ecological databases

BACKGROUND

The availability of big data in healthcare and the intensive development of data reuse and georeferencing have opened up perspectives for health spatial analysis. However, fine-scale spatial studies of ecological and medical databases are limited by the change of support problem and thus a lack of spatial unit interoperability. The use of spatial disaggregation methods to solve this problem introduces errors into the spatial estimations. Here, we present a generic, two-step method for merging medical and ecological databases that avoids the use of spatial disaggregation methods, while maximizing the spatial resolution.

METHODS

Firstly, a mapping table is created after one or more transition matrices have been defined. The latter link the spatial units of the original databases to the spatial units of the final database. Secondly, the mapping table is validated by (1) comparing the covariates contained in the two original databases, and (2) checking the spatial validity with a spatial continuity criterion and a spatial resolution index.

RESULTS

We used our novel method to merge a medical database (the French national diagnosis-related group database, containing 5644 spatial units) with an ecological database (produced by the French National Institute of Statistics and Economic Studies, and containing with 36,594 spatial units). The mapping table yielded 5632 final spatial units. The mapping table's validity was evaluated by comparing the number of births in the medical database and the ecological databases in each final spatial unit. The median [interquartile range] relative difference was 2.3% [0; 5.7]. The spatial continuity criterion was low (2.4%), and the spatial resolution index was greater than for most French administrative areas.

CONCLUSIONS

Our innovative approach improves interoperability between medical and ecological databases and facilitates fine-scale spatial analyses. We have shown that disaggregation models and large aggregation techniques are not necessarily the best ways to tackle the change of support problem.

Environments Tune and Select Cellular Diversity

Technical advances in single-cell sequencing data and their application to greater samples is revealing substantial cell-to-cell variation in expression levels and propagation of this variation between molecules across cells. New quantitative approaches that apply mechanistic and statistical models in a systems-wide approach are illuminating the drivers of phenotypic diversity.

Meta-analysis of crowdsourced data compendia suggests pan-disease transcriptional signatures of autoimmunity

Background: The proliferation of publicly accessible large-scale biological data together with increasing availability of bioinformatics tools have the potential to transform biomedical research. Here we report a crowdsourcing Jamboree that explored whether a team of volunteer biologists without formal bioinformatics training could use OMiCC, a crowdsourcing web platform that facilitates the reuse and (meta-) analysis of public gene expression data, to compile and annotate gene expression data, and design comparisons between disease and control sample groups.

Methods: The Jamboree focused on several common human autoimmune diseases, including systemic lupus erythematosus (SLE), multiple sclerosis (MS), type I diabetes (DM1), and rheumatoid arthritis (RA), and the corresponding mouse models. Meta-analyses were performed in OMiCC using comparisons constructed by the participants to identify 1) gene expression signatures for each disease (disease versus healthy controls at the gene expression and biological pathway levels), 2) conserved signatures across all diseases within each species (pan-disease signatures), and 3) conserved signatures between species for each disease and across all diseases (cross-species signatures).

Results: A large number of differentially expressed genes were identified for each disease based on meta-analysis, with observed overlap among diseases both within and across species. Gene set/pathway enrichment of upregulated genes suggested conserved signatures (e.g., interferon) across all human and mouse conditions.

Conclusions: Our Jamboree exercise provides evidence that when enabled by appropriate tools, a "crowd" of biologists can work together to accelerate the pace by which the increasingly large amounts of public data can be reused and meta-analyzed for generating and testing hypotheses. Our encouraging experience suggests that a similar crowdsourcing approach can be used to explore other biological questions.