Multi-modal generative modeling for joint analysis of single-cell T cell receptor and gene expression data.

Recent advances in single-cell immune profiling have enabled the simultaneous measurement of transcriptome and T cell receptor (TCR) sequences, offering great potential for studying immune responses at the cellular level. However, integrating these diverse modalities across datasets is challenging due to their unique data characteristics and technical variations. Here, to address this, we develop the multimodal generative model mvTCR to fuse modality-specific information across transcriptome and TCR into a shared representation. Our analysis demonstrates the added value of multimodal over unimodal approaches to capture antigen specificity. Notably, we use mvTCR to distinguish T cell subpopulations binding to SARS-CoV-2 antigens from bystander cells. Furthermore, when combined with reference mapping approaches, mvTCR can map newly generated datasets to extensive T cell references, facilitating knowledge transfer. In summary, we envision mvTCR to enable a scalable analysis of multimodal immune profiling data and advance our understanding of immune responses.

Human SARS-CoV-2 challenge uncovers local and systemic response dynamics.

The COVID-19 pandemic is an ongoing global health threat, yet our understanding of the dynamics of early cellular responses to this disease remains limited1. Here in our SARS-CoV-2 human challenge study, we used single-cell multi-omics profiling of nasopharyngeal swabs and blood to temporally resolve abortive, transient and sustained infections in seronegative individuals challenged with pre-Alpha SARS-CoV-2. Our analyses revealed rapid changes in cell-type proportions and dozens of highly dynamic cellular response states in epithelial and immune cells associated with specific time points and infection status. We observed that the interferon response in blood preceded the nasopharyngeal response. Moreover, nasopharyngeal immune infiltration occurred early in samples from individuals with only transient infection and later in samples from individuals with sustained infection. High expression of HLA-DQA2 before inoculation was associated with preventing sustained infection. Ciliated cells showed multiple immune responses and were most permissive for viral replication, whereas nasopharyngeal T cells and macrophages were infected non-productively. We resolved 54 T cell states, including acutely activated T cells that clonally expanded while carrying convergent SARS-CoV-2 motifs. Our new computational pipeline Cell2TCR identifies activated antigen-responding T cells based on a gene expression signature and clusters these into clonotype groups and motifs. Overall, our detailed time series data can serve as a Rosetta stone for epithelial and immune cell responses and reveals early dynamic responses associated with protection against infection.

Enhanced CD95 and interleukin 18 signalling accompany T cell receptor Vß21.3+ activation in multi-inflammatory syndrome in children.

Multisystem inflammatory syndrome in children is a post-infectious presentation SARS-CoV-2 associated with expansion of the T cell receptor Vß21.3+ T-cell subgroup. Here we apply muti-single cell omics to compare the inflammatory process in children with acute respiratory COVID-19 and those presenting with non SARS-CoV-2 infections in children. Here we show that in Multi-Inflammatory Syndrome in Children (MIS-C), the natural killer cell and monocyte population demonstrate heightened CD95 (Fas) and Interleuking 18 receptor expression. Additionally, TCR Vß21.3+ CD4+ T-cells exhibit skewed differentiation towards T helper 1, 17 and regulatory T cells, with increased expression of the co-stimulation receptors ICOS, CD28 and interleukin 18 receptor. We observe no functional evidence for NLRP3 inflammasome pathway overactivation, though MIS-C monocytes show elevated active caspase 8. This, coupled with raised IL18 mRNA expression in CD16- NK cells on single cell RNA sequencing analysis, suggests interleukin 18 and CD95 signalling may trigger activation of TCR Vß21.3+ T-cells in MIS-C, driven by increased IL-18 production from activated monocytes and CD16- Natural Killer cells.

Dandelion uses the single-cell adaptive immune receptor repertoire to explore lymphocyte developmental origins.

Assessment of single-cell gene expression (single-cell RNA sequencing) and adaptive immune receptor (AIR) sequencing (scVDJ-seq) has been invaluable in studying lymphocyte biology. Here we introduce Dandelion, a computational pipeline for scVDJ-seq analysis. It enables the application of standard V(D)J analysis workflows to single-cell datasets, delivering improved V(D)J contig annotation and the identification of nonproductive and partially spliced contigs. We devised a strategy to create an AIR feature space that can be used for both differential V(D)J usage analysis and pseudotime trajectory inference. The application of Dandelion improved the alignment of human thymic development trajectories of double-positive T cells to mature single-positive CD4/CD8 T cells, generating predictions of factors regulating lineage commitment. Dandelion analysis of other cell compartments provided insights into the origins of human B1 cells and ILC/NK cell development, illustrating the power of our approach. Dandelion is available at https://www.github.com/zktuong/dandelion .

Atypical B cells and impaired SARS-CoV-2 neutralization following heterologous vaccination in the elderly.

Suboptimal responses to a primary vaccination course have been reported in the elderly, but there is little information regarding the impact of age on responses to booster third doses. Here, we show that individuals 70 years or older (median age 73, range 70-75) who received a primary two-dose schedule with AZD1222 and booster third dose with mRNA vaccine achieve significantly lower neutralizing antibody responses against SARS-CoV-2 spike pseudotyped virus compared with those younger than 70 (median age 66, range 54-69) at 1 month post booster. Impaired neutralization potency and breadth post third dose in the elderly is associated with circulating "atypical" spike-specific B cells expressing CD11c and FCRL5. However, when considering individuals who received three doses of mRNA vaccine, we did not observe differences in neutralization or enrichment in atypical B cells. This work highlights the finding that AdV and mRNA COVID-19 vaccine formats differentially instruct the memory B cell response.

Molecular architecture of the developing mouse brain.

The mammalian brain develops through a complex interplay of spatial cues generated by diffusible morphogens, cell-cell interactions and intrinsic genetic programs that result in probably more than a thousand distinct cell types. A complete understanding of this process requires a systematic characterization of cell states over the entire spatiotemporal range of brain development. The ability of single-cell RNA sequencing and spatial transcriptomics to reveal the molecular heterogeneity of complex tissues has therefore been particularly powerful in the nervous system. Previous studies have explored development in specific brain regions1-8, the whole adult brain9 and even entire embryos10. Here we report a comprehensive single-cell transcriptomic atlas of the embryonic mouse brain between gastrulation and birth. We identified almost eight hundred cellular states that describe a developmental program for the functional elements of the brain and its enclosing membranes, including the early neuroepithelium, region-specific secondary organizers, and both neurogenic and gliogenic progenitors. We also used in situ mRNA sequencing to map the spatial expression patterns of key developmental genes. Integrating the in situ data with our single-cell clusters revealed the precise spatial organization of neural progenitors during the patterning of the nervous system.