Multiomics integration of 22 immune-mediated monogenic diseases reveals an emergent axis of human immune health.

Monogenic diseases are often studied in isolation due to their rarity. Here we utilize multiomics to assess 22 monogenic immune-mediated conditions with age- and sex-matched healthy controls. Despite clearly detectable disease-specific and "pan-disease" signatures, individuals possess stable personal immune states over time. Temporally stable differences among subjects tend to dominate over differences attributable to disease conditions or medication use. Unsupervised principal variation analysis of personal immune states and machine learning classification distinguishing between healthy controls and patients converge to a metric of immune health (IHM). The IHM discriminates healthy from multiple polygenic autoimmune and inflammatory disease states in independent cohorts, marks healthy aging, and is a pre-vaccination predictor of antibody responses to influenza vaccination in the elderly. We identified easy-to-measure circulating protein biomarker surrogates of the IHM that capture immune health variations beyond age. Our work provides a conceptual framework and biomarkers for defining and measuring human immune health.

Influenza vaccination reveals sex dimorphic imprints of prior mild COVID-19.

Acute viral infections can have durable functional impacts on the immune system long after recovery, but how they affect homeostatic immune states and responses to future perturbations remain poorly understood1-4. Here we use systems immunology approaches, including longitudinal multimodal single-cell analysis (surface proteins, transcriptome and V(D)J sequences) to comparatively assess baseline immune statuses and responses to influenza vaccination in 33 healthy individuals after recovery from mild, non-hospitalized COVID-19 (mean, 151 days after diagnosis) and 40 age- and sex-matched control individuals who had never had COVID-19. At the baseline and independent of time after COVID-19, recoverees had elevated T cell activation signatures and lower expression of innate immune genes including Toll-like receptors in monocytes. Male individuals who had recovered from COVID-19 had coordinately higher innate, influenza-specific plasmablast, and antibody responses after vaccination compared with healthy male individuals and female individuals who had recovered from COVID-19, in part because male recoverees had monocytes with higher IL-15 responses early after vaccination coupled with elevated prevaccination frequencies of 'virtual memory'-like CD8+ T cells poised to produce more IFNγ after IL-15 stimulation. Moreover, the expression of the repressed innate immune genes in monocytes increased by day 1 to day 28 after vaccination in recoverees, therefore moving towards the prevaccination baseline of the healthy control individuals. By contrast, these genes decreased on day 1 and returned to the baseline by day 28 in the control individuals. Our study reveals sex-dimorphic effects of previous mild COVID-19 and suggests that viral infections in humans can establish new immunological set-points that affect future immune responses in an antigen-agnostic manner.

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).

Early human B cell signatures of the primary antibody response to mRNA vaccination.

Messenger RNA (mRNA) vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are highly effective at inducing protective immunity. However, weak antibody responses are seen in some individuals, and cellular correlates of immunity remain poorly defined, especially for B cells. Here we used unbiased approaches to longitudinally dissect primary antibody, plasmablast, and memory B cell (MBC) responses to the two-dose mRNA-1273 vaccine in SARS-CoV-2-naive adults. Coordinated immunoglobulin A (IgA) and IgG antibody responses were preceded by bursts of spike-specific plasmablasts after both doses but earlier and more intensely after dose 2. While antibody and B cell cellular responses were generally robust, they also varied within the cohort and decreased over time after a dose-2 peak. Both antigen-nonspecific postvaccination plasmablast frequency after dose 1 and their spike-specific counterparts early after dose 2 correlated with subsequent antibody levels. This correlation between early plasmablasts and antibodies remained for titers measured at 6 months after vaccination. Several distinct antigen-specific MBC populations emerged postvaccination with varying kinetics, including two MBC populations that correlated with 2- and 6-month antibody titers. Both were IgG-expressing MBCs: one less mature, appearing as a correlate after the first dose, while the other MBC correlate showed a more mature and resting phenotype, emerging as a correlate later after dose 2. This latter MBC was also a major contributor to the sustained spike-specific MBC response observed at month 6. Thus, these plasmablasts and MBCs that emerged after both the first and second doses with distinct kinetics are potential determinants of the magnitude and durability of antibodies in response to mRNA-based vaccination.

Influenza vaccination and single cell multiomics reveal sex dimorphic immune imprints of prior mild COVID-19.

Viral infections can have profound and durable functional impacts on the immune system. There is an urgent need to characterize the long-term immune effects of SARS-CoV-2 infection given the persistence of symptoms in some individuals and the continued threat of novel variants. Here we use systems immunology, including longitudinal multimodal single cell analysis (surface proteins, transcriptome, and V(D)J sequences) from 33 previously healthy individuals after recovery from mild, non-hospitalized COVID-19 and 40 age- and sex-matched healthy controls with no history of COVID-19 to comparatively assess the post-infection immune status (mean: 151 days after diagnosis) and subsequent innate and adaptive responses to seasonal influenza vaccination. Identification of both sex-specific and -independent temporally stable changes, including signatures of T-cell activation and repression of innate defense/immune receptor genes (e.g., Toll-like receptors) in monocytes, suggest that mild COVID-19 can establish new post-recovery immunological set-points. COVID-19-recovered males had higher innate, influenza-specific plasmablast, and antibody responses after vaccination compared to healthy males and COVID-19-recovered females, partly attributable to elevated pre-vaccination frequencies of a GPR56 expressing CD8+ T-cell subset in male recoverees that are "poised" to produce higher levels of IFNγ upon inflammatory stimulation. Intriguingly, by day 1 post-vaccination in COVID-19-recovered subjects, the expression of the repressed genes in monocytes increased and moved towards the pre-vaccination baseline of healthy controls, suggesting that the acute inflammation induced by vaccination could partly reset the immune states established by mild COVID-19. Our study reveals sex-dimorphic immune imprints and in vivo functional impacts of mild COVID-19 in humans, suggesting that prior COVID-19, and possibly respiratory viral infections in general, could change future responses to vaccination and in turn, vaccines could help reset the immune system after COVID-19, both in an antigen-agnostic manner.

Pre-vaccination and early B cell signatures predict antibody response to SARS-CoV-2 mRNA vaccine.

SARS-CoV-2 mRNA vaccines are highly effective, although weak antibody responses are seen in some individuals with correlates of immunity that remain poorly understood. Here we longitudinally dissected antibody, plasmablast, and memory B cell (MBC) responses to the two-dose Moderna mRNA vaccine in SARS-CoV-2-uninfected adults. Robust, coordinated IgA and IgG antibody responses were preceded by bursts of spike-specific plasmablasts after both doses, but earlier and more intensely after dose two. Distinct antigen-specific MBC populations also emerged post-vaccination with varying kinetics. We identified antigen non-specific pre-vaccination MBC and post-vaccination plasmablasts after dose one and their spike-specific counterparts early after dose two that correlated with subsequent antibody levels. These baseline and response signatures can thus provide early indicators of serological efficacy and explain response variability in the population.

Time-resolved systems immunology reveals a late juncture linked to fatal COVID-19.

COVID-19 exhibits extensive patient-to-patient heterogeneity. To link immune response variation to disease severity and outcome over time, we longitudinally assessed circulating proteins as well as 188 surface protein markers, transcriptome, and T cell receptor sequence simultaneously in single peripheral immune cells from COVID-19 patients. Conditional-independence network analysis revealed primary correlates of disease severity, including gene expression signatures of apoptosis in plasmacytoid dendritic cells and attenuated inflammation but increased fatty acid metabolism in CD56dimCD16hi NK cells linked positively to circulating interleukin (IL)-15. CD8+ T cell activation was apparent without signs of exhaustion. Although cellular inflammation was depressed in severe patients early after hospitalization, it became elevated by days 17-23 post symptom onset, suggestive of a late wave of inflammatory responses. Furthermore, circulating protein trajectories at this time were divergent between and predictive of recovery versus fatal outcomes. Our findings stress the importance of timing in the analysis, clinical monitoring, and therapeutic intervention of COVID-19.

Environment Tunes Propagation of Cell-to-Cell Variation in the Human Macrophage Gene Network.

Cell-to-cell variation in gene expression and the propagation of such variation (PoV or "noise propagation") from one gene to another in the gene network, as reflected by gene-gene correlation across single cells, are commonly observed in single-cell transcriptomic studies and can shape the phenotypic diversity of cell populations. While gene network "rewiring" is known to accompany cellular adaptation to different environments, how PoV changes between environments and its underlying regulatory mechanisms are less understood. Here, we systematically explored context-dependent PoV among genes in human macrophages, utilizing different cytokines as natural perturbations of multiple molecular parameters that may influence PoV. Our single-cell, epigenomic, computational, and stochastic simulation analyses reveal that environmental adaptation can tune PoV to potentially shape cellular heterogeneity by changing parameters such as the degree of phosphorylation and transcription factor-chromatin interactions. This quantitative tuning of PoV may be a widespread, yet underexplored, property of cellular adaptation to distinct environments.

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

New technologies have been propelling dramatic increases in the volume and diversity of large-scale public data, which can potentially be reused to answer questions beyond those originally envisioned. However, this often requires computational and statistical skills beyond the reach of most bench scientists. The development of educational and accessible computational tools is thus critical, as are crowdsourcing efforts that utilize the community's expertise to curate public data for hypothesis generation and testing. Here we review the history of public-data reuse and argue for greater incorporation of computational and statistical sciences into the biomedical education curriculum and the development of biologist-friendly crowdsourcing tools. Finally, we provide a resource list for the reuse of public data and highlight an illustrative crowdsourcing exercise to explore public gene-expression data of human autoimmune diseases and corresponding mouse models. Through education, tool development, and community engagement, immunologists will be poised to transform public data into biological insights.

Humoral Fingerprinting of Immune Responses: 'Super-Resolution', High-Dimensional Serology.

In a recent study, Chung et al. report the development of a high-dimensional approach to assess humoral responses to immune perturbation that goes beyond antibody neutralization and titers. This approach enables the identification of potentially novel correlates and mechanisms of protective immunity to HIV vaccination, thus offering a glimpse of how dense phenotyping of serological responses coupled with bioinformatics analysis could lead to much-sought-after markers of protective vaccination responses.

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.

Three-Dimensional Spot Detection in Ratiometric Fluorescence Imaging For Measurement of Subcellular Organelles.

Lysosomes are subcellular organelles playing a vital role in the endocytosis process of the cell. Lysosomal acidity is an important factor in assuring proper functioning of the enzymes within the organelle, and can be assessed by labeling the lysosomes with pH-sensitive fluorescence probes. To enhance our understanding of the acidification mechanisms, the goal of this work is to develop a method that can accurately detect and characterize the acidity of each lysosome captured in ratiometric fluorescence images. We present an algorithm that utilizes the h-dome transformation and reconciles spots detected independently from two wavelength channels. We evaluated our algorithm using simulated images for which the exact locations were known. The h-dome algorithm achieved an f-score as high as 0.890. We also computed the fluorescence ratios from lysosomes in live HeLa cell images with known lysosomal pHs. Using leave-one-out cross-validation, we demonstrated that the new algorithm was able to achieve much better pH prediction accuracy than the conventional method.

Overview of BioCreative II gene normalization.

BACKGROUND: The goal of the gene normalization task is to link genes or gene products mentioned in the literature to biological databases. This is a key step in an accurate search of the biological literature. It is a challenging task, even for the human expert; genes are often described rather than referred to by gene symbol and, confusingly, one gene name may refer to different genes (often from different organisms). For BioCreative II, the task was to list the Entrez Gene identifiers for human genes or gene products mentioned in PubMed/MEDLINE abstracts. We selected abstracts associated with articles previously curated for human genes. We provided 281 expert-annotated abstracts containing 684 gene identifiers for training, and a blind test set of 262 documents containing 785 identifiers, with a gold standard created by expert annotators. Inter-annotator agreement was measured at over 90%. RESULTS: Twenty groups submitted one to three runs each, for a total of 54 runs. Three systems achieved F-measures (balanced precision and recall) between 0.80 and 0.81. Combining the system outputs using simple voting schemes and classifiers obtained improved results; the best composite system achieved an F-measure of 0.92 with 10-fold cross-validation. A 'maximum recall' system based on the pooled responses of all participants gave a recall of 0.97 (with precision 0.23), identifying 763 out of 785 identifiers. CONCLUSION: Major advances for the BioCreative II gene normalization task include broader participation (20 versus 8 teams) and a pooled system performance comparable to human experts, at over 90% agreement. These results show promise as tools to link the literature with biological databases.

Introducing meta-services for biomedical information extraction.

We introduce the first meta-service for information extraction in molecular biology, the BioCreative MetaServer (BCMS; http://bcms.bioinfo.cnio.es/). This prototype platform is a joint effort of 13 research groups and provides automatically generated annotations for PubMed/Medline abstracts. Annotation types cover gene names, gene IDs, species, and protein-protein interactions. The annotations are distributed by the meta-server in both human and machine readable formats (HTML/XML). This service is intended to be used by biomedical researchers and database annotators, and in biomedical language processing. The platform allows direct comparison, unified access, and result aggregation of the annotations.

Rule-based human gene normalization in biomedical text with confidence estimation.

The ability to identify gene mentions in text and normalize them to the proper unique identifiers is crucial for "down-stream" text mining applications in bioinformatics. We have developed a rule-based algorithm that divides the normalization task into two steps. The first step includes pattern matching for gene symbols and an approximate term searching technique for gene names. Next, the algorithm measures several features based on morphological, statistical, and contextual information to estimate the level of confidence that the correct identifier is selected for a potential mention. Uniqueness, inverse distance, and coverage are three novel features we quantified. The algorithm was evaluated against the BioCreAtIvE datasets. The feature weights were tuned by the Nealder-Mead simplex method. An F-score of .7622 and an AUC (area under the recall-precision curve) of .7461 were achieved on the test data using the set of weights optimized to the training data.