Search

ObjectiveThere is an established yet unexplained link between interferon (IFN) and systemic lupus erythematosus (SLE). The expression of sequences derived from transposable elements (TEs) may contribute to SLE phenotypes, specifically production of type I IFNs and generation of autoantibodies.MethodsWe profiled cell-sorted RNA-sequencing data (CD4+ T cells, CD14+ monocytes, CD19+ B cells, and natural killer cells) from peripheral blood mononuclear cells of 120 patients with SLE and quantified TE expression identifying 27,135 TEs. We tested for differential TE expression across 10 SLE phenotypes, including autoantibody production and disease activity.ResultsWe found 731 differentially expressed (DE) TEs across all SLE phenotypes that were mostly cell specific and phenotype specific. DE TEs were enriched for specific families and open reading frames of viral genes encoded in TE sequences. Increased expression of DE TEs was associated with genes involved in antiviral activity, such as LY6E, ISG15, and TRIM22, and pathways such as IFN signaling.ConclusionThese findings suggest that expression of TEs contributes to activation of SLE-related mechanisms in a cell-specific manner, which can impact disease diagnostics and therapeutics.

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.

BACKGROUNDPatients hospitalized for COVID-19 exhibit diverse clinical outcomes, with outcomes for some individuals diverging over time even though their initial disease severity appears similar to that of other patients. A systematic evaluation of molecular and cellular profiles over the full disease course can link immune programs and their coordination with progression heterogeneity.METHODSWe performed deep immunophenotyping and conducted longitudinal multiomics modeling, integrating 10 assays for 1,152 Immunophenotyping Assessment in a COVID-19 Cohort (IMPACC) study participants and identifying several immune cascades that were significant drivers of differential clinical outcomes.RESULTSIncreasing disease severity was driven by a temporal pattern that began with the early upregulation of immunosuppressive metabolites and then elevated levels of inflammatory cytokines, signatures of coagulation, formation of neutrophil extracellular traps, and T cell functional dysregulation. A second immune cascade, predictive of 28-day mortality among critically ill patients, was characterized by reduced total plasma Igs and B cells and dysregulated IFN responsiveness. We demonstrated that the balance disruption between IFN-stimulated genes and IFN inhibitors is a crucial biomarker of COVID-19 mortality, potentially contributing to failure of viral clearance in patients with fatal illness.CONCLUSIONOur longitudinal multiomics profiling study revealed temporal coordination across diverse omics that potentially explain the disease progression, providing insights that can inform the targeted development of therapies for patients hospitalized with COVID-19, especially those who are critically ill.TRIAL REGISTRATIONClinicalTrials.gov NCT04378777.FUNDINGNIH (5R01AI135803-03, 5U19AI118608-04, 5U19AI128910-04, 4U19AI090023-11, 4U19AI118610-06, R01AI145835-01A1S1, 5U19AI062629-17, 5U19AI057229-17, 5U19AI125357-05, 5U19AI128913-03, 3U19AI077439-13, 5U54AI142766-03, 5R01AI104870-07, 3U19AI089992-09, 3U19AI128913-03, and 5T32DA018926-18); NIAID, NIH (3U19AI1289130, U19AI128913-04S1, and R01AI122220); and National Science Foundation (DMS2310836).

Rheumatoid arthritis (RA) management leans toward achieving remission or low disease activity. In this study, we conducted single-cell RNA sequencing (scRNA-Seq) of peripheral blood mononuclear cells (PBMCs) from 36 individuals (18 patients with RA and 18 matched controls, accounting for age, sex, race, and ethnicity), to identify disease-relevant cell subsets and cell type-specific signatures associated with disease activity. Our analysis revealed 18 distinct PBMC subsets, including an IFN-induced transmembrane 3-overexpressing (IFITM3-overexpressing) IFN-activated monocyte subset. We observed an increase in CD4+ T effector memory cells in patients with moderate-high disease activity (DAS28-CRP ≥ 3.2) and a decrease in nonclassical monocytes in patients with low disease activity or remission (DAS28-CRP < 3.2). Pseudobulk analysis by cell type identified 168 differentially expressed genes between RA and matched controls, with a downregulation of proinflammatory genes in the γδ T cell subset, alteration of genes associated with RA predisposition in the IFN-activated subset, and nonclassical monocytes. Additionally, we identified a gene signature associated with moderate-high disease activity, characterized by upregulation of proinflammatory genes such as TNF, JUN, EGR1, IFIT2, MAFB, and G0S2 and downregulation of genes including HLA-DQB1, HLA-DRB5, and TNFSF13B. Notably, cell-cell communication analysis revealed an upregulation of signaling pathways, including VISTA, in both moderate-high and remission-low disease activity contexts. Our findings provide valuable insights into the systemic cellular and molecular mechanisms underlying RA disease activity.

The glycosylation of IgG plays a critical role during human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, activating immune cells and inducing cytokine production. However, the role of IgM N-glycosylation has not been studied during human acute viral infection. The analysis of IgM N-glycosylation from healthy controls and hospitalized coronavirus disease 2019 (COVID-19) patients reveals increased high-mannose and sialylation that correlates with COVID-19 severity. These trends are confirmed within SARS-CoV-2-specific immunoglobulin N-glycan profiles. Moreover, the degree of total IgM mannosylation and sialylation correlate significantly with markers of disease severity. We link the changes of IgM N-glycosylation with the expression of Golgi glycosyltransferases. Lastly, we observe antigen-specific IgM antibody-dependent complement deposition is elevated in severe COVID-19 patients and modulated by exoglycosidase digestion. Taken together, this work links the IgM N-glycosylation with COVID-19 severity and highlights the need to understand IgM glycosylation and downstream immune function during human disease.

Post-acute sequelae of SARS-CoV-2 (PASC) is a significant public health concern. We describe Patient Reported Outcomes (PROs) on 590 participants prospectively assessed from hospital admission for COVID-19 through one year after discharge. Modeling identified 4 PRO clusters based on reported deficits (minimal, physical, mental/cognitive, and multidomain), supporting heterogenous clinical presentations in PASC, with sub-phenotypes associated with female sex and distinctive comorbidities. During the acute phase of disease, a higher respiratory SARS-CoV-2 viral burden and lower Receptor Binding Domain and Spike antibody titers were associated with both the physical predominant and the multidomain deficit clusters. A lower frequency of circulating B lymphocytes by mass cytometry (CyTOF) was observed in the multidomain deficit cluster. Circulating fibroblast growth factor 21 (FGF21) was significantly elevated in the mental/cognitive predominant and the multidomain clusters. Future efforts to link PASC to acute anti-viral host responses may help to better target treatment and prevention of PASC.

Altered myeloid inflammation and lymphopenia are hallmarks of severe infections. We identified the upregulated EN-RAGE gene program in airway and blood myeloid cells from patients with acute lung injury from SARS-CoV-2 or other causes across 7 cohorts. This program was associated with greater clinical severity and predicted future mechanical ventilation and death. EN-RAGEhi myeloid cells express features consistent with suppressor cell functionality, including low HLA-DR and high PD-L1. Sustained EN-RAGE program expression in airway and blood myeloid cells correlated with clinical severity and increasing expression of T cell dysfunction markers. IL-6 upregulated many EN-RAGE program genes in monocytes in vitro. IL-6 signaling blockade by tocilizumab in a placebo-controlled clinical trial led to rapid normalization of EN-RAGE and T cell gene expression. This identifies IL-6 as a key driver of myeloid dysregulation associated with worse clinical outcomes in COVID-19 patients and provides insights into shared pathophysiological mechanisms in non-COVID-19 ARDS.

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.

The IMPACC cohort, composed of >1,000 hospitalized COVID-19 participants, contains five illness trajectory groups (TGs) during acute infection (first 28 days), ranging from milder (TG1-3) to more severe disease course (TG4) and death (TG5). Here, we report deep immunophenotyping, profiling of >15,000 longitudinal blood and nasal samples from 540 participants of the IMPACC cohort, using 14 distinct assays. These unbiased analyses identify cellular and molecular signatures present within 72 h of hospital admission that distinguish moderate from severe and fatal COVID-19 disease. Importantly, cellular and molecular states also distinguish participants with more severe disease that recover or stabilize within 28 days from those that progress to fatal outcomes (TG4 vs. TG5). Furthermore, our longitudinal design reveals that these biologic states display distinct temporal patterns associated with clinical outcomes. Characterizing host immune responses in relation to heterogeneity in disease course may inform clinical prognosis and opportunities for intervention.

Ancestral SARS coronavirus-2 (SARS-CoV-2) and variants of concern (VOC) caused a global pandemic with a spectrum of disease severity. The mechanistic explaining variations related to airway epithelium are relatively understudied. Here, we biobanked airway organoids (AO) by preserving stem cell function. We optimized viral infection with H1N1/PR8 and comprehensively characterized epithelial responses to SARS-CoV-2 infection in phenotypically stable AO from 20 different subjects. We discovered Tetraspanin-8 (TSPAN8) as a facilitator of SARS-CoV-2 infection. TSPAN8 facilitates SARS-CoV-2 infection rates independently of ACE2-Spike interaction. In head-to-head comparisons with Ancestral SARS-CoV-2, Delta and Omicron VOC displayed lower overall infection rates of AO but triggered changes in epithelial response. All variants shared highest tropism for ciliated and goblet cells. TSPAN8-blocking antibodies diminish SARS-CoV-2 infection and may spur novel avenues for COVID-19 therapy.

In the past decade, high-dimensional single-cell technologies have revolutionized basic and translational immunology research and are now a key element of the toolbox used by scientists to study the immune system. However, analysis of the data generated by these approaches often requires clustering algorithms and dimensionality reduction representation, which are computationally intense and difficult to evaluate and optimize. Here, we present Cytometry Clustering Optimization and Evaluation (Cyclone), an analysis pipeline integrating dimensionality reduction, clustering, evaluation, and optimization of clustering resolution, and downstream visualization tools facilitating the analysis of a wide range of cytometry data. We benchmarked and validated Cyclone on mass cytometry (CyTOF), full-spectrum fluorescence-based cytometry, and multiplexed immunofluorescence (IF) in a variety of biological contexts, including infectious diseases and cancer. In each instance, Cyclone not only recapitulates gold standard immune cell identification but also enables the unsupervised identification of lymphocytes and mononuclear phagocyte subsets that are associated with distinct biological features. Altogether, the Cyclone pipeline is a versatile and accessible pipeline for performing, optimizing, and evaluating clustering on a variety of cytometry datasets, which will further power immunology research and provide a scaffold for biological discovery.

Rationale: Autopsy and biomarker studies suggest that endotheliopathy contributes to coronavirus disease (COVID-19)-associated acute respiratory distress syndrome. However, the effects of COVID-19 on the lung endothelium are not well defined. We hypothesized that the lung endotheliopathy of COVID-19 is caused by circulating host factors and direct endothelial infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Objectives: We aimed to determine the effects of SARS-CoV-2 or sera from patients with COVID-19 on the permeability and inflammatory activation of lung microvascular endothelial cells. Methods: Human lung microvascular endothelial cells were treated with live SARS-CoV-2; inactivated viral particles; or sera from patients with COVID-19, patients without COVID-19, and healthy volunteers. Permeability was determined by measuring transendothelial resistance to electrical current flow, where decreased resistance signifies increased permeability. Inflammatory mediators were quantified in culture supernatants. Endothelial biomarkers were quantified in patient sera. Measurements and Main Results: Viral PCR confirmed that SARS-CoV-2 enters and replicates in endothelial cells. Live SARS-CoV-2, but not dead virus or spike protein, induces endothelial permeability and secretion of plasminogen activator inhibitor 1 and vascular endothelial growth factor. There was substantial variability in the effects of SARS-CoV-2 on endothelial cells from different donors. Sera from patients with COVID-19 induced endothelial permeability, which correlated with disease severity. Serum levels of endothelial activation and injury biomarkers were increased in patients with COVID-19 and correlated with severity of illness. Conclusions: SARS-CoV-2 infects and dysregulates endothelial cell functions. Circulating factors in patients with COVID-19 also induce endothelial cell dysfunction. Our data point to roles for both systemic factors acting on lung endothelial cells and viral infection of endothelial cells in COVID-19-associated endotheliopathy.

BackgroundBetter understanding of the association between characteristics of patients hospitalized with coronavirus disease 2019 (COVID-19) and outcome is needed to further improve upon patient management.MethodsImmunophenotyping Assessment in a COVID-19 Cohort (IMPACC) is a prospective, observational study of 1164 patients from 20 hospitals across the United States. Disease severity was assessed using a 7-point ordinal scale based on degree of respiratory illness. Patients were prospectively surveyed for 1 year after discharge for post-acute sequalae of COVID-19 (PASC) through quarterly surveys. Demographics, comorbidities, radiographic findings, clinical laboratory values, SARS-CoV-2 PCR and serology were captured over a 28-day period. Multivariable logistic regression was performed.FindingsThe median age was 59 years (interquartile range [IQR] 20); 711 (61%) were men; overall mortality was 14%, and 228 (20%) required invasive mechanical ventilation. Unsupervised clustering of ordinal score over time revealed distinct disease course trajectories. Risk factors associated with prolonged hospitalization or death by day 28 included age ≥ 65 years (odds ratio [OR], 2.01; 95% CI 1.28-3.17), Hispanic ethnicity (OR, 1.71; 95% CI 1.13-2.57), elevated baseline creatinine (OR 2.80; 95% CI 1.63- 4.80) or troponin (OR 1.89; 95% 1.03-3.47), baseline lymphopenia (OR 2.19; 95% CI 1.61-2.97), presence of infiltrate by chest imaging (OR 3.16; 95% CI 1.96-5.10), and high SARS-CoV2 viral load (OR 1.53; 95% CI 1.17-2.00). Fatal cases had the lowest ratio of SARS-CoV-2 antibody to viral load levels compared to other trajectories over time (p=0.001). 589 survivors (51%) completed at least one survey at follow-up with 305 (52%) having at least one symptom consistent with PASC, most commonly dyspnea (56% among symptomatic patients). Female sex was the only associated risk factor for PASC.InterpretationIntegration of PCR cycle threshold, and antibody values with demographics, comorbidities, and laboratory/radiographic findings identified risk factors for 28-day outcome severity, though only female sex was associated with PASC. Longitudinal clinical phenotyping offers important insights, and provides a framework for immunophenotyping for acute and long COVID-19.FundingNIH.

While studies have elucidated many pathophysiological elements of COVID-19, little is known about immunological changes during COVID-19 resolution. We analyzed immune cells and phosphorylated signaling states at single-cell resolution from longitudinal blood samples of patients hospitalized with COVID-19, pneumonia and/or sepsis, and healthy individuals by mass cytometry. COVID-19 patients showed distinct immune compositions and an early, coordinated, and elevated immune cell signaling profile associated with early hospital discharge. Intra-patient longitudinal analysis revealed changes in myeloid and T cell frequencies and a reduction in immune cell signaling across cell types that accompanied disease resolution and discharge. These changes, together with increases in regulatory T cells and reduced signaling in basophils, also accompanied recovery from respiratory failure and were associated with better outcomes at time of admission. Therefore, although patients have heterogeneous immunological baselines and highly variable disease courses, a core immunological trajectory exists that defines recovery from severe SARS-CoV-2 infection.

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.

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.

Cancers display significant heterogeneity with respect to tissue of origin, driver mutations, and other features of the surrounding tissue. It is likely that individual tumors engage common patterns of the immune system-here "archetypes"-creating prototypical non-destructive tumor immune microenvironments (TMEs) and modulating tumor-targeting. To discover the dominant immune system archetypes, the University of California, San Francisco (UCSF) Immunoprofiler Initiative (IPI) processed 364 individual tumors across 12 cancer types using standardized protocols. Computational clustering of flow cytometry and transcriptomic data obtained from cell sub-compartments uncovered dominant patterns of immune composition across cancers. These archetypes were profound insofar as they also differentiated tumors based upon unique immune and tumor gene-expression patterns. They also partitioned well-established classifications of tumor biology. The IPI resource provides a template for understanding cancer immunity as a collection of dominant patterns of immune organization and provides a rational path forward to learn how to modulate these to improve therapy.

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.

Alzheimer's Disease (AD) is a complex neurodegenerative disease that gravely affects patients and imposes an immense burden on caregivers. Apolipoprotein E4 (APOE4) has been identified as the most common genetic risk factor for AD, yet the molecular mechanisms connecting APOE4 to AD are not well understood. Past transcriptomic analyses in AD have revealed APOE genotype-specific transcriptomic differences; however, these differences have not been explored at a single-cell level. To elucidate more complex APOE genotype-specific disease-relevant changes masked by the bulk analysis, we leverage the first two single-nucleus RNA sequencing AD datasets from human brain samples, including nearly 55,000 cells from the prefrontal and entorhinal cortices. In each brain region, we performed a case versus control APOE genotype-stratified differential gene expression analysis and pathway network enrichment in astrocytes, microglia, neurons, oligodendrocytes, and oligodendrocyte progenitor cells. We observed more global transcriptomic changes in APOE4 positive AD cells and identified differences across APOE genotypes primarily in glial cell types. Our findings highlight the differential transcriptomic perturbations of APOE isoforms at a single-cell level in AD pathogenesis and have implications for precision medicine development in the diagnosis and treatment of AD.

BackgroundSince the beginning of the coronavirus disease 2019 (COVID-19) pandemic, there has been increasing urgency to identify pathophysiological characteristics leading to severe clinical course in patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Human leukocyte antigen alleles (HLA) have been suggested as potential genetic host factors that affect individual immune response to SARS-CoV-2. We sought to evaluate this hypothesis by conducting a multicenter study using HLA sequencing.MethodsWe analyzed the association between COVID-19 severity and HLAs in 435 individuals from Germany (n = 135), Spain (n = 133), Switzerland (n = 20) and the United States (n = 147), who had been enrolled from March 2020 to August 2020. This study included patients older than 18 years, diagnosed with COVID-19 and representing the full spectrum of the disease. Finally, we tested our results by meta-analysing data from prior genome-wide association studies (GWAS).FindingsWe describe a potential association of HLA-C*04:01 with severe clinical course of COVID-19. Carriers of HLA-C*04:01 had twice the risk of intubation when infected with SARS-CoV-2 (risk ratio 1.5 [95% CI 1.1-2.1], odds ratio 3.5 [95% CI 1.9-6.6], adjusted p-value = 0.0074). These findings are based on data from four countries and corroborated by independent results from GWAS. Our findings are biologically plausible, as HLA-C*04:01 has fewer predicted bindings sites for relevant SARS-CoV-2 peptides compared to other HLA alleles.InterpretationHLA-C*04:01 carrier state is associated with severe clinical course in SARS-CoV-2. Our findings suggest that HLA class I alleles have a relevant role in immune defense against SARS-CoV-2.FundingFunded by Roche Sequencing Solutions, Inc.

Neutralizing autoantibodies against type I interferons (IFNs) have been found in some patients with critical coronavirus disease 2019 (COVID-19), the disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, the prevalence of these antibodies, their longitudinal dynamics across the disease severity scale, and their functional effects on circulating leukocytes remain unknown. Here, in 284 patients with COVID-19, we found type I IFN–specific autoantibodies in peripheral blood samples from 19% of patients with critical disease and 6% of patients with severe disease. We found no type I IFN autoantibodies in individuals with moderate disease. Longitudinal profiling of over 600,000 peripheral blood mononuclear cells using multiplexed single-cell epitope and transcriptome sequencing from 54 patients with COVID-19 and 26 non–COVID-19 controls revealed a lack of type I IFN–stimulated gene (ISG-I) responses in myeloid cells from patients with critical disease. This was especially evident in dendritic cell populations isolated from patients with critical disease producing type I IFN–specific autoantibodies. Moreover, we found elevated expression of the inhibitory receptor leukocyte-associated immunoglobulin-like receptor 1 (LAIR1) on the surface of monocytes isolated from patients with critical disease early in the disease course. LAIR1 expression is inversely correlated with ISG-I expression response in patients with COVID-19 but is not expressed in healthy controls. The deficient ISG-I response observed in patients with critical COVID-19 with and without type I IFN–specific autoantibodies supports a unifying model for disease pathogenesis involving ISG-I suppression through convergent mechanisms.

Acute myeloid leukemia patients refractory to induction therapy or relapsed within one year have poor outcomes. Autocrine production of hepatocyte growth factor by myeloid blasts drives leukemogenesis in pre-clinical models. A phase Ib trial evaluated ficlatuzumab, a first-in-class anti-HGF antibody, in combination with cytarabine in this high-risk population. Dose-limiting toxicities were not observed, and 20 mg/kg was established as the recommended phase II dose. The most frequent treatment-related adverse event was febrile neutropenia. Among 17 evaluable patients, the overall response rate was 53%, all complete remissions. Phospho-proteomic mass cytometry showed potent on-target suppression of p-MET after ficlatuzumab treatment and that attenuation of p-S6 was associated with clinical response. Multiplexed single cell RNA sequencing using prospectively acquired patient specimens identified interferon response genes as adverse predictive factors. The ficlatuzumab and cytarabine combination is well-tolerated with favorable efficacy. High-dimensional analyses at single-cell resolution represent promising approaches for identifying biomarkers of response and mechanisms of resistance in prospective clinical studies.

SummaryA 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 implementationdittoSeq is an R package available through Bioconductor via an open source MIT license.Supplementary informationSupplementary data are available at Bioinformatics online.

Although infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has pleiotropic and systemic effects in some individuals1-3, many others experience milder symptoms. Here, to gain a more comprehensive understanding of the distinction between severe and mild phenotypes in the pathology of coronavirus disease 2019 (COVID-19) and its origins, we performed a whole-blood-preserving single-cell analysis protocol to integrate contributions from all major immune cell types of the blood-including neutrophils, monocytes, platelets, lymphocytes and the contents of the serum. Patients with mild COVID-19 exhibit a coordinated pattern of expression of interferon-stimulated genes (ISGs)3 across every cell population, whereas these ISG-expressing cells are systemically absent in patients with severe disease. Paradoxically, individuals with severe COVID-19 produce very high titres of anti-SARS-CoV-2 antibodies and have a lower viral load compared to individuals with mild disease. Examination of the serum from patients with severe COVID-19 shows that these patients uniquely produce antibodies that functionally block the production of the ISG-expressing cells associated with mild disease, by activating conserved signalling circuits that dampen cellular responses to interferons. Overzealous antibody responses pit the immune system against itself in many patients with COVID-19, and perhaps also in individuals with other viral infections. Our findings reveal potential targets for immunotherapies in patients with severe COVID-19 to re-engage viral defence.

To study disease development, an inventory of an organ's cell types and understanding of physiologic function is paramount. Here, we performed single-cell RNA-sequencing to examine heterogeneity of murine pancreatic duct cells, pancreatobiliary cells, and intrapancreatic bile duct cells. We describe an epithelial-mesenchymal transitory axis in our three pancreatic duct subpopulations and identify osteopontin as a regulator of this fate decision as well as human duct cell dedifferentiation. Our results further identify functional heterogeneity within pancreatic duct subpopulations by elucidating a role for geminin in accumulation of DNA damage in the setting of chronic pancreatitis. Our findings implicate diverse functional roles for subpopulations of pancreatic duct cells in maintenance of duct cell identity and disease progression and establish a comprehensive road map of murine pancreatic duct cell, pancreatobiliary cell, and intrapancreatic bile duct cell homeostasis.

Influenza is a significant cause of morbidity and mortality worldwide. Here we show changes in the abundance and activation states of more than 50 immune cell subsets in 35 individuals over 11 time points during human A/California/2009 (H1N1) virus challenge monitored using mass cytometry along with other clinical assessments. Peak change in monocyte, B cell, and T cell subset frequencies coincided with peak virus shedding, followed by marked activation of T and NK cells. Results led to the identification of CD38 as a critical regulator of plasmacytoid dendritic cell function in response to influenza virus. Machine learning using study-derived clinical parameters and single-cell data effectively classified and predicted susceptibility to infection. The coordinated immune cell dynamics defined in this study provide a framework for identifying novel correlates of protection in the evaluation of future influenza 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.

MotivationHigh-parameter single-cell technologies can reveal novel cell populations of interest, but studying or validating these populations using lower-parameter methods remains challenging.ResultsHere, we present GateFinder, an algorithm that enriches high-dimensional cell types with simple, stepwise polygon gates requiring only two markers at a time. A series of case studies of complex cell types illustrates how simplified enrichment strategies can enable more efficient assays, reveal novel biomarkers and clarify underlying biology.Availability and implementationThe GateFinder algorithm is implemented as a free and open-source package for BioConductor: https://nalab.stanford.edu/gatefinder.Supplementary informationSupplementary 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.

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.

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.