Trimodal single-cell profiling reveals a novel pediatric CD8αα+ T cell subset and broad age-related molecular reprogramming across the T cell compartment.

Age-associated changes in the T cell compartment are well described. However, limitations of current single-modal or bimodal single-cell assays, including flow cytometry, RNA-seq (RNA sequencing) and CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing), have restricted our ability to deconvolve more complex cellular and molecular changes. Here, we profile >300,000 single T cells from healthy children (aged 11-13 years) and older adults (aged 55-65 years) by using the trimodal assay TEA-seq (single-cell analysis of mRNA transcripts, surface protein epitopes and chromatin accessibility), which revealed that molecular programming of T cell subsets shifts toward a more activated basal state with age. Naive CD4+ T cells, considered relatively resistant to aging, exhibited pronounced transcriptional and epigenetic reprogramming. Moreover, we discovered a novel CD8αα+ T cell subset lost with age that is epigenetically poised for rapid effector responses and has distinct inhibitory, costimulatory and tissue-homing properties. Together, these data reveal new insights into age-associated changes in the T cell compartment that may contribute to differential immune responses.

Antibodies targeting surface membrane antigens in patients with chronic graft-versus-host disease.

Chronic graft-versus-host disease (cGVHD) after allogeneic hematopoietic stem cell transplant reflects a complex immune response resulting in chronic damage to multiple tissues. Previous studies indicated that donor B cells and the antibodies they produce play an important role in the development of cGVHD. To understand the pathogenic role of antibodies in cGVHD, we focused our studies on posttransplant production of immunoglobulin G antibodies targeting cell surface antigens expressed in multiple cGVHD affected tissues, due to their potential functional impact on living cells in vivo. Using plate-bound cell membrane proteins as targets, we detected a significantly higher level of antibodies reactive with these membrane antigens in patients who developed cGVHD, compared with those who did not and healthy donors. Plasma-reactive antibody levels increased significantly prior to the clinical diagnosis of cGVHD and were reduced following cGVHD therapies including prednisone, interleukin-2, or extracorporeal photophoresis. Using cell-based immunoprecipitation with plasma from cGVHD patients and mass spectrometry, we identified 43 membrane proteins targeted by these antibodies. The presence of antibodies in cGVHD patients' plasma that specifically target 6 of these proteins was validated. Antibodies reactive with these 6 antigens were more frequently detected in patients with cGVHD compared with patients without cGVHD and healthy donors. These results indicate that antibodies that target membrane antigens of living cells frequently develop in cGVHD patients and further support a role for B cells and antibodies in the development of cGVHD.

Abatacept increases T cell exhaustion in early RA individuals who carry HLA risk alleles.

Exhausted CD8 T cells (TEX) are associated with worse outcome in cancer yet better outcome in autoimmunity. Building on our past findings of increased TIGIT+KLRG1+ TEX with teplizumab therapy in type 1 diabetes (T1D), in the absence of treatment we found that the frequency of TIGIT+KLRG1+ TEX is stable within an individual but differs across individuals in both T1D and healthy control (HC) cohorts. This TIGIT+KLRG1+ CD8 TEX population shares an exhaustion-associated EOMES gene signature in HC, T1D, rheumatoid arthritis (RA), and cancer subjects, expresses multiple inhibitory receptors, and is hyporesponsive in vitro, together suggesting co-expression of TIGIT and KLRG1 may broadly define human peripheral exhausted cells. In HC and RA subjects, lower levels of EOMES transcriptional modules and frequency of TIGIT+KLRG1+ TEX were associated with RA HLA risk alleles (DR0401, 0404, 0405, 0408, 1001) even when considering disease status and cytomegalovirus (CMV) seropositivity. Moreover, the frequency of TIGIT+KLRG1+ TEX was significantly increased in RA HLA risk but not non-risk subjects treated with abatacept (CTLA4Ig). The DR4 association and selective modulation with abatacept suggests that therapeutic modulation of TEX may be more effective in DR4 subjects and TEX may be indirectly influenced by cellular interactions that are blocked by abatacept.

A comprehensive platform for analyzing longitudinal multi-omics data.

Longitudinal bulk and single-cell omics data is increasingly generated for biological and clinical research but is challenging to analyze due to its many intrinsic types of variations. We present PALMO ( https://github.com/aifimmunology/PALMO ), a platform that contains five analytical modules to examine longitudinal bulk and single-cell multi-omics data from multiple perspectives, including decomposition of sources of variations within the data, collection of stable or variable features across timepoints and participants, identification of up- or down-regulated markers across timepoints of individual participants, and investigation on samples of same participants for possible outlier events. We have tested PALMO performance on a complex longitudinal multi-omics dataset of five data modalities on the same samples and six external datasets of diverse background. Both PALMO and our longitudinal multi-omics dataset can be valuable resources to the scientific community.

Optimized workflow for human PBMC multiomic immunosurveillance studies.

Deep immune profiling is essential for understanding the human immune system in health and disease. Successful biological interpretation of this data requires consistent laboratory processing with minimal batch-to-batch variation. Here, we detail a robust pipeline for the profiling of human peripheral blood mononuclear cells by both high-dimensional flow cytometry and single-cell RNA-seq. These protocols reduce batch effects, generate reproducible data, and increase throughput. For complete details on the use and execution of this protocol, please refer to Savage et al. (2021).

Simultaneous trimodal single-cell measurement of transcripts, epitopes, and chromatin accessibility using TEA-seq.

Single-cell measurements of cellular characteristics have been instrumental in understanding the heterogeneous pathways that drive differentiation, cellular responses to signals, and human disease. Recent advances have allowed paired capture of protein abundance and transcriptomic state, but a lack of epigenetic information in these assays has left a missing link to gene regulation. Using the heterogeneous mixture of cells in human peripheral blood as a test case, we developed a novel scATAC-seq workflow that increases signal-to-noise and allows paired measurement of cell surface markers and chromatin accessibility: integrated cellular indexing of chromatin landscape and epitopes, called ICICLE-seq. We extended this approach using a droplet-based multiomics platform to develop a trimodal assay that simultaneously measures transcriptomics (scRNA-seq), epitopes, and chromatin accessibility (scATAC-seq) from thousands of single cells, which we term TEA-seq. Together, these multimodal single-cell assays provide a novel toolkit to identify type-specific gene regulation and expression grounded in phenotypically defined cell types.

Multimodal analysis for human ex vivo studies shows extensive molecular changes from delays in blood processing.

Multi-omic profiling of human peripheral blood is increasingly utilized to identify biomarkers and pathophysiologic mechanisms of disease. The importance of these platforms in clinical and translational studies led us to investigate the impact of delayed blood processing on the numbers and state of peripheral blood mononuclear cells (PBMC) and on the plasma proteome. Similar to previous studies, we show minimal effects of delayed processing on the numbers and general phenotype of PBMC up to 18 hours. In contrast, profound changes in the single-cell transcriptome and composition of the plasma proteome become evident as early as 6 hours after blood draw. These reflect patterns of cellular activation across diverse cell types that lead to progressive distancing of the gene expression state and plasma proteome from native in vivo biology. Differences accumulating during an overnight rest (18 hours) could confound relevant biologic variance related to many underlying disease states.

Longitudinal immune dynamics of mild COVID-19 define signatures of recovery and persistence.

SARS-CoV-2 has infected over 200 million and caused more than 4 million deaths to date. Most individuals (>80%) have mild symptoms and recover in the outpatient setting, but detailed studies of immune responses have focused primarily on moderate to severe COVID-19. We deeply profiled the longitudinal immune response in individuals with mild COVID-19 beginning with early time points post-infection (1-15 days) and proceeding through convalescence to >100 days after symptom onset. We correlated data from single cell analyses of peripheral blood cells, serum proteomics, virus-specific cellular and humoral immune responses, and clinical metadata. Acute infection was characterized by vigorous coordinated innate and adaptive immune activation that differed in character by age (young vs. old). We then characterized signals associated with recovery and convalescence to define and validate a new signature of inflammatory cytokines, gene expression, and chromatin accessibility that persists in individuals with post-acute sequelae of SARS-CoV-2 infection (PASC).

Genomic analyses of flow-sorted Hodgkin Reed-Sternberg cells reveal complementary mechanisms of immune evasion.

Classical Hodgkin lymphoma (cHL) is composed of rare malignant Hodgkin Reed-Sternberg (HRS) cells within an extensive, but ineffective, inflammatory/immune cell infiltrate. HRS cells exhibit near-universal somatic copy gains of chromosome 9p/9p24.1, which increase expression of the programmed cell death protein 1 (PD-1) ligands. To define genetic mechanisms of response and resistance to PD-1 blockade and identify complementary treatment targets, we performed whole-exome sequencing of flow cytometry-sorted HRS cells from 23 excisional biopsies of newly diagnosed cHLs, including 8 Epstein-Barr virus-positive (EBV+) tumors. We identified significantly mutated cancer candidate genes (CCGs) as well as somatic copy number alterations and structural variations and characterized their contribution to disease-defining immune evasion mechanisms and nuclear factor κB (NF-κB), JAK/STAT, and PI3K signaling pathways. EBV- cHLs had a higher prevalence of genetic alterations in the NF-κB and major histocompatibility complex class I antigen presentation pathways. In this young cHL cohort (median age, 26 years), we identified a predominant mutational signature of spontaneous deamination of cytosine- phosphate-guanines ("Aging"), in addition to apolipoprotein B mRNA editing catalytic polypeptide-like, activation-induced cytidine deaminase, and microsatellite instability (MSI)-associated hypermutation. In particular, the mutational burden in EBV- cHLs was among the highest reported, similar to that of carcinogen-induced tumors. Together, the overall high mutational burden, MSI-associated hypermutation, and newly identified genetic alterations represent additional potential bases for the efficacy of PD-1 blockade in cHL. Of note, recurrent cHL alterations, including B2M, TNFAIP3, STAT6, GNA13, and XPO1 mutations and 2p/2p15, 6p21.32, 6q23.3, and 9p/9p24.1 copy number alterations, were also identified in >20% of primary mediastinal B-cell lymphomas, highlighting shared pathogenetic mechanisms in these diseases.