Single-cell RNA sequencing methods have led to improved understanding of the heterogeneity and transcriptomic states present in complex biological systems. Recently, the development of novel single-cell technologies for assaying additional modalities, specifically genomic, epigenomic, proteomic, and spatial data, allows for unprecedented insight into cellular biology. While certain technologies collect multiple measurements from the same cells simultaneously, even when modalities are separately assayed in different cells, we can apply novel computational methods to integrate these data. The application of computational integration methods to multimodal paired and unpaired data results in rich information about the identities of the cells present and the interactions between different levels of biology, such as between genetic variation and transcription. In this review, we both discuss the single-cell technologies for measuring these modalities and describe and characterize a variety of computational integration methods for combining the resulting data to leverage multimodal information toward greater biological insight.

Summary

A visualization suite for major forms of bulk and single-cell RNAseq data in R. dittoSeq is color blindness-friendly by default, robustly documented to power ease-of-use and allows highly customizable generation of both daily-use and publication-quality figures.

Availability and implementation

dittoSeq is an R package available through Bioconductor via an open source MIT license.

Supplementary information

Supplementary data are available at Bioinformatics online.

Sulfomenaquinone (SMK) is a recently identified metabolite that is unique to the Mycobacterium tuberculosis (M. tuberculosis) complex and is shown to modulate its virulence. Here, we report the identification of the SMK biosynthetic operon that, in addition to a previously identified sulfotransferase stf3, includes a putative cytochrome P450 gene (cyp128) and a gene of unknown function, rv2269c. We demonstrate that cyp128 and stf3 are sufficient for the biosynthesis of SMK from menaquinone and rv2269c exhibits promoter activity in M. tuberculosis. Loss of Stf3 expression, but not that of Cyp128, is correlated with elevated levels of menaquinone-9, an essential component in the electron-transport chain in M. tuberculosis. Finally, we showed in a mouse model of infection that the loss of cyp128 exhibits a hypervirulent phenotype similar to that in previous studies of the stf3 mutant. These findings provide a platform for defining the molecular basis of SMK's role in M. tuberculosis pathogenesis.

Assessing the health and competence of the immune system is central to evaluating vaccination responses, autoimmune conditions, cancer prognosis, and treatment. With an increasing number of studies examining immune dysregulation, there is a growing need for a curated reference of variation in immune parameters in healthy individuals. We used mass cytometry (CyTOF) to profile blood from 86 humans in response to 15 ex vivo immune stimuli. We present reference ranges for cell-specific immune markers and highlight differences that appear across sex and age. We identified modules of immune features that suggest there exists an underlying structure to the immune system based on signaling pathway responses across cell types. We observed increased MAPK signaling in inflammatory pathways in innate immune cells and greater overall coordination of immune cell responses in females. In contrast, males exhibited stronger pSTAT1 and pTBK1 responses. These reference data are publicly available as a resource for immune profiling studies.

Animal models are an integral part of the drug development and evaluation process. However, they are unsurprisingly imperfect reflections of humans, and the extent and nature of many immunological differences are unknown. With the rise of targeted and biological therapeutics, it is increasingly important that we understand the molecular differences in the immunological behavior of humans and model organisms. However, very few antibodies are raised against non-human primate antigens, and databases of cross-reactivity between species are incomplete. Thus, we screened 332 antibodies in five immune cell populations in blood from humans and four non-human primate species generating a comprehensive cross-reactivity catalog that includes cell type-specificity. We used this catalog to create large mass cytometry universal cross-species phenotyping and signaling panels for humans, along with three of the model organisms most similar to humans: rhesus and cynomolgus macaques and African green monkeys; and one of the mammalian models most widely used in drug development: C57BL/6 mice. As a proof-of-principle, we measured immune cell signaling responses across all five species to an array of 15 stimuli using mass cytometry. We found numerous instances of different cellular phenotypes and immune signaling events occurring within and between species, and detailed three examples (double-positive T cell frequency and signaling; granulocyte response to Bacillus anthracis antigen; and B cell subsets). We also explore the correlation of herpes simian B virus serostatus on the immune profile. Antibody panels and the full dataset generated are available online as a resource to enable future studies comparing immune responses across species during the evaluation of therapeutics.

High-dimensional single-cell technologies present new opportunities for biological discovery, but the complex nature of the resulting datasets makes it challenging to perform comprehensive analysis. One particular challenge is the analysis of single-cell time course datasets: how to identify unique cell populations and track how they change across time points. To facilitate this analysis, we developed FLOW-MAP, a graphical user interface (GUI)-based software tool that uses graph layout analysis with sequential time ordering to visualize cellular trajectories in high-dimensional single-cell datasets obtained from flow cytometry, mass cytometry or single-cell RNA sequencing (scRNAseq) experiments. Here we provide a detailed description of the FLOW-MAP algorithm and how to use the open-source R package FLOWMAPR via its GUI or with text-based commands. This approach can be applied to many dynamic processes, including in vitro stem cell differentiation, in vivo development, oncogenesis, the emergence of drug resistance and cell signaling dynamics. To demonstrate our approach, we perform a step-by-step analysis of a single-cell mass cytometry time course dataset from mouse embryonic stem cells differentiating into the three germ layers: endoderm, mesoderm and ectoderm. In addition, we demonstrate FLOW-MAP analysis of a previously published scRNAseq dataset. Using both synthetic and experimental datasets for comparison, we perform FLOW-MAP analysis side by side with other single-cell analysis methods, to illustrate when it is advantageous to use the FLOW-MAP approach. The protocol takes between 30 min and 1.5 h to complete.

The study of traditional populations provides a view of human-associated microbes unperturbed by industrialization, as well as a window into the microbiota that co-evolved with humans. Here we discuss our recent work characterizing the microbiota from the Hadza hunter-gatherers of Tanzania. We found seasonal shifts in bacterial taxa, diversity, and carbohydrate utilization by the microbiota. When compared to the microbiota composition from other populations around the world, the Hadza microbiota shares bacterial families with other traditional societies that are rare or absent from microbiotas of industrialized nations. We present additional observations from the Hadza microbiota and their lifestyle and environment, including microbes detected on hands, water, and animal sources, how the microbiota varies with sex and age, and the short-term effects of introducing agricultural products into the diet. In the context of our previously published findings and of these additional observations, we discuss a path forward for future work.

Alzheimer's disease (AD) is a pervasive neurodegenerative disorder that disproportionately affects women. Since neural anatomy and disease pathophysiology differ by sex, investigating sex-specific mechanisms in AD pathophysiology can inform new therapeutic approaches for both sexes. Previous bulk human brain RNA sequencing studies have revealed sex differences in dysregulated molecular pathways related to energy production, neuronal function, and immune response; however, the sex differences in disease mechanisms are yet to be examined comprehensively on a single-cell level. We leveraged nearly 74,000 cells from human prefrontal and entorhinal cortex samples from the first two publicly available single-cell RNA sequencing AD datasets to perform a case versus control sex-stratified differential gene expression analysis and pathway network enrichment in a cell type-specific manner for each brain region. Our examination at the single-cell level revealed sex differences in AD prominently in glial cells of the prefrontal cortex. In the entorhinal cortex, we observed the same genes and networks to be perturbed in opposing directions between sexes in AD relative to healthy state. Our findings contribute to growing evidence of sex differences in AD-related transcriptomic changes, which can fuel the development of therapies that may prove more effective at reversing AD pathophysiology.

Juvenile dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I IFN response and autoantibodies. Treatment options are limited due to an incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of patients with JDM at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment toward an immature naive state as a hallmark of JDM at diagnosis. Furthermore, we find that these changes in B cells are paralleled by T cell signatures suggestive of Th2-mediated inflammation that persist despite disease quiescence. We applied network analysis to reveal that hyperactivation of the type I IFN response in all immune populations is coordinated with previously masked cell states including dysfunctional protein processing in CD4+ T cells and regulation of cell death programming in NK cells, CD8+ T cells, and γδ T cells. Together, these findings unveil the coordinated immune dysregulation underpinning JDM and provide insight into strategies for restoring balance in immune function.

Immune cells function in an interacting hierarchy that coordinates the activities of various cell types according to genetic and environmental contexts. We developed graphical approaches to construct an extensible immune reference map from mass cytometry data of cells from different organs, incorporating landmark cell populations as flags on the map to compare cells from distinct samples. The maps recapitulated canonical cellular phenotypes and revealed reproducible, tissue-specific deviations. The approach revealed influences of genetic variation and circadian rhythms on immune system structure, enabled direct comparisons of murine and human blood cell phenotypes, and even enabled archival fluorescence-based flow cytometry data to be mapped onto the reference framework. This foundational reference map provides a working definition of systemic immune organization to which new data can be integrated to reveal deviations driven by genetics, environment, or pathology.