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.

Age-specific nasal epithelial responses to SARS-CoV-2 infection.

Children infected with SARS-CoV-2 rarely progress to respiratory failure. However, the risk of mortality in infected people over 85 years of age remains high. Here we investigate differences in the cellular landscape and function of paediatric (<12 years), adult (30-50 years) and older adult (>70 years) ex vivo cultured nasal epithelial cells in response to infection with SARS-CoV-2. We show that cell tropism of SARS-CoV-2, and expression of ACE2 and TMPRSS2 in nasal epithelial cell subtypes, differ between age groups. While ciliated cells are viral replication centres across all age groups, a distinct goblet inflammatory subtype emerges in infected paediatric cultures and shows high expression of interferon-stimulated genes and incomplete viral replication. In contrast, older adult cultures infected with SARS-CoV-2 show a proportional increase in basaloid-like cells, which facilitate viral spread and are associated with altered epithelial repair pathways. We confirm age-specific induction of these cell types by integrating data from in vivo COVID-19 studies and validate that our in vitro model recapitulates early epithelial responses to SARS-CoV-2 infection.

Decoding chromatin states by proteomic profiling of nucleosome readers.

DNA and histone modifications combine into characteristic patterns that demarcate functional regions of the genome1,2. While many 'readers' of individual modifications have been described3-5, how chromatin states comprising composite modification signatures, histone variants and internucleosomal linker DNA are interpreted is a major open question. Here we use a multidimensional proteomics strategy to systematically examine the interaction of around 2,000 nuclear proteins with over 80 modified dinucleosomes representing promoter, enhancer and heterochromatin states. By deconvoluting complex nucleosome-binding profiles into networks of co-regulated proteins and distinct nucleosomal features driving protein recruitment or exclusion, we show comprehensively how chromatin states are decoded by chromatin readers. We find highly distinctive binding responses to different features, many factors that recognize multiple features, and that nucleosomal modifications and linker DNA operate largely independently in regulating protein binding to chromatin. Our online resource, the Modification Atlas of Regulation by Chromatin States (MARCS), provides in-depth analysis tools to engage with our results and advance the discovery of fundamental principles of genome regulation by chromatin states.

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 .

Early human lung immune cell development and its role in epithelial cell fate.

Studies of human lung development have focused on epithelial and mesenchymal cell types and function, but much less is known about the developing lung immune cells, even though the airways are a major site of mucosal immunity after birth. An unanswered question is whether tissue-resident immune cells play a role in shaping the tissue as it develops in utero. Here, we profiled human embryonic and fetal lung immune cells using scRNA-seq, smFISH, and immunohistochemistry. At the embryonic stage, we observed an early wave of innate immune cells, including innate lymphoid cells, natural killer cells, myeloid cells, and lineage progenitors. By the canalicular stage, we detected naive T lymphocytes expressing high levels of cytotoxicity genes and the presence of mature B lymphocytes, including B-1 cells. Our analysis suggests that fetal lungs provide a niche for full B cell maturation. Given the presence and diversity of immune cells during development, we also investigated their possible effect on epithelial maturation. We found that IL-1ß drives epithelial progenitor exit from self-renewal and differentiation to basal cells in vitro. In vivo, IL-1ß-producing myeloid cells were found throughout the lung and adjacent to epithelial tips, suggesting that immune cells may direct human lung epithelial development.

Acid ceramidase regulates innate immune memory.

Innate immune memory, also called "trained immunity," is a functional state of myeloid cells enabling enhanced immune responses. This phenomenon is important for host defense, but also plays a role in various immune-mediated conditions. We show that exogenously administered sphingolipids and inhibition of sphingolipid metabolizing enzymes modulate trained immunity. In particular, we reveal that acid ceramidase, an enzyme that converts ceramide to sphingosine, is a potent regulator of trained immunity. We show that acid ceramidase regulates the transcription of histone-modifying enzymes, resulting in profound changes in histone 3 lysine 27 acetylation and histone 3 lysine 4 trimethylation. We confirm our findings by identifying single-nucleotide polymorphisms in the region of ASAH1, the gene encoding acid ceramidase, that are associated with the trained immunity cytokine response. Our findings reveal an immunomodulatory effect of sphingolipids and identify acid ceramidase as a relevant therapeutic target to modulate trained immunity responses in innate immune-driven disorders.

Quantification of absolute transcription factor binding affinities in the native chromatin context using BANC-seq.

Transcription factor binding across the genome is regulated by DNA sequence and chromatin features. However, it is not yet possible to quantify the impact of chromatin context on transcription factor binding affinities. Here, we report a method called binding affinities to native chromatin by sequencing (BANC-seq) to determine absolute apparent binding affinities of transcription factors to native DNA across the genome. In BANC-seq, a concentration range of a tagged transcription factor is added to isolated nuclei. Concentration-dependent binding is then measured per sample to quantify apparent binding affinities across the genome. BANC-seq adds a quantitative dimension to transcription factor biology, which enables stratification of genomic targets based on transcription factor concentration and prediction of transcription factor binding sites under non-physiological conditions, such as disease-associated overexpression of (onco)genes. Notably, whereas consensus DNA binding motifs for transcription factors are important to establish high-affinity binding sites, these motifs are not always strictly required to generate nanomolar-affinity interactions in the genome.

A spatially resolved atlas of the human lung characterizes a gland-associated immune niche.

Single-cell transcriptomics has allowed unprecedented resolution of cell types/states in the human lung, but their spatial context is less well defined. To (re)define tissue architecture of lung and airways, we profiled five proximal-to-distal locations of healthy human lungs in depth using multi-omic single cell/nuclei and spatial transcriptomics (queryable at lungcellatlas.org ). Using computational data integration and analysis, we extend beyond the suspension cell paradigm and discover macro and micro-anatomical tissue compartments including previously unannotated cell types in the epithelial, vascular, stromal and nerve bundle micro-environments. We identify and implicate peribronchial fibroblasts in lung disease. Importantly, we discover and validate a survival niche for IgA plasma cells in the airway submucosal glands (SMG). We show that gland epithelial cells recruit B cells and IgA plasma cells, and promote longevity and antibody secretion locally through expression of CCL28, APRIL and IL-6. This new 'gland-associated immune niche' has implications for respiratory health.

Obesity Is Associated with Attenuated Tissue Immunity in COVID-19.

Rationale: Obesity affects 40% of U.S. adults, is associated with a proinflammatory state, and presents a significant risk factor for the development of severe coronavirus disease (COVID-19). To date, there is limited information on how obesity might affect immune cell responses in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Objectives: To determine the impact of obesity on respiratory tract immunity in COVID-19 across the human lifespan. Methods: We analyzed single-cell transcriptomes from BAL in three ventilated adult cohorts with (n = 24) or without (n = 9) COVID-19 from nasal immune cells in children with (n = 14) or without (n = 19) COVID-19, and from peripheral blood mononuclear cells in an independent adult COVID-19 cohort (n = 42), comparing obese and nonobese subjects. Measurements and Main Results: Surprisingly, we found that obese adult subjects had attenuated lung immune or inflammatory responses in SARS-CoV-2 infection, with decreased expression of IFN-α, IFN-γ, and TNF-α (tumor necrosis factor α) response gene signatures in almost all lung epithelial and immune cell subsets, and lower expression of IFNG and TNF in specific lung immune cells. Peripheral blood immune cells in an independent adult cohort showed a similar but less marked reduction in type-I IFN and IFNγ response genes, as well as decreased serum IFNα, in obese patients with SARS-CoV-2. Nasal immune cells from obese children with COVID-19 also showed reduced enrichment of IFN-α and IFN-γ response genes. Conclusions: These findings show blunted tissue immune responses in obese patients with COVID-19, with implications for treatment stratification, supporting the specific application of inhaled recombinant type-I IFNs in this vulnerable subset.

Liver Colonization by Colorectal Cancer Metastases Requires YAP-Controlled Plasticity at the Micrometastatic Stage.

Micrometastases of colorectal cancer can remain dormant for years prior to the formation of actively growing, clinically detectable lesions (i.e., colonization). A better understanding of this step in the metastatic cascade could help improve metastasis prevention and treatment. Here we analyzed liver specimens of patients with colorectal cancer and monitored real-time metastasis formation in mouse livers using intravital microscopy to reveal that micrometastatic lesions are devoid of cancer stem cells (CSC). However, lesions that grow into overt metastases demonstrated appearance of de novo CSCs through cellular plasticity at a multicellular stage. Clonal outgrowth of patient-derived colorectal cancer organoids phenocopied the cellular and transcriptomic changes observed during in vivo metastasis formation. First, formation of mature CSCs occurred at a multicellular stage and promoted growth. Conversely, failure of immature CSCs to generate more differentiated cells arrested growth, implying that cellular heterogeneity is required for continuous growth. Second, early-stage YAP activity was required for the survival of organoid-forming cells. However, subsequent attenuation of early-stage YAP activity was essential to allow for the formation of cell type heterogeneity, while persistent YAP signaling locked micro-organoids in a cellularly homogenous and growth-stalled state. Analysis of metastasis formation in mouse livers using single-cell RNA sequencing confirmed the transient presence of early-stage YAP activity, followed by emergence of CSC and non-CSC phenotypes, irrespective of the initial phenotype of the metastatic cell of origin. Thus, establishment of cellular heterogeneity after an initial YAP-controlled outgrowth phase marks the transition to continuously growing macrometastases. SIGNIFICANCE: Characterization of the cell type dynamics, composition, and transcriptome of early colorectal cancer liver metastases reveals that failure to establish cellular heterogeneity through YAP-controlled epithelial self-organization prohibits the outgrowth of micrometastases. See related commentary by LeBleu, p. 1870.

Local and systemic responses to SARS-CoV-2 infection in children and adults.

It is not fully understood why COVID-19 is typically milder in children1-3. Here, to examine the differences between children and adults in their response to SARS-CoV-2 infection, we analysed paediatric and adult patients with COVID-19 as well as healthy control individuals (total n = 93) using single-cell multi-omic profiling of matched nasal, tracheal, bronchial and blood samples. In the airways of healthy paediatric individuals, we observed cells that were already in an interferon-activated state, which after SARS-CoV-2 infection was further induced especially in airway immune cells. We postulate that higher paediatric innate interferon responses restrict viral replication and disease progression. The systemic response in children was characterized by increases in naive lymphocytes and a depletion of natural killer cells, whereas, in adults, cytotoxic T cells and interferon-stimulated subpopulations were significantly increased. We provide evidence that dendritic cells initiate interferon signalling in early infection, and identify epithelial cell states associated with COVID-19 and age. Our matching nasal and blood data show a strong interferon response in the airways with the induction of systemic interferon-stimulated populations, which were substantially reduced in paediatric patients. Together, we provide several mechanisms that explain the milder clinical syndrome observed in children.

Retinoic acid signaling drives differentiation toward the absorptive lineage in colorectal cancer.

Retinoic acid (RA) signaling is an important and conserved pathway that regulates cellular proliferation and differentiation. Furthermore, perturbed RA signaling is implicated in cancer initiation and progression. However, the mechanisms by which RA signaling contributes to homeostasis, malignant transformation, and disease progression in the intestine remain incompletely understood. Here, we report, in agreement with previous findings, that activation of the Retinoic Acid Receptor and the Retinoid X Receptor results in enhanced transcription of enterocyte-specific genes in mouse small intestinal organoids. Conversely, inhibition of this pathway results in reduced expression of genes associated with the absorptive lineage. Strikingly, this latter effect is conserved in a human organoid model for colorectal cancer (CRC) progression. We further show that RXR motif accessibility depends on progression state of CRC organoids. Finally, we show that reduced RXR target gene expression correlates with worse CRC prognosis, implying RA signaling as a putative therapeutic target in CRC.

Towards a Human Cell Atlas: Taking Notes from the Past.

Comprehensively characterizing the cellular composition and organization of tissues has been a long-term scientific challenge that has limited our ability to study fundamental and clinical aspects of human physiology. The Human Cell Atlas (HCA) is a global collaborative effort to create a reference map of all human cells as a basis for both understanding human health and diagnosing, monitoring, and treating disease. Many aspects of the HCA are analogous to the Human Genome Project (HGP), whose completion presents a major milestone in modern biology. To commemorate the HGP's 20-year anniversary of completion, we discuss the launch of the HCA in light of the HGP, and highlight recent progress by the HCA consortium.

Single-cell multi-omics analysis of the immune response in COVID-19.

Analysis of human blood immune cells provides insights into the coordinated response to viral infections such as severe acute respiratory syndrome coronavirus 2, which causes coronavirus disease 2019 (COVID-19). We performed single-cell transcriptome, surface proteome and T and B lymphocyte antigen receptor analyses of over 780,000 peripheral blood mononuclear cells from a cross-sectional cohort of 130 patients with varying severities of COVID-19. We identified expansion of nonclassical monocytes expressing complement transcripts (CD16+C1QA/B/C+) that sequester platelets and were predicted to replenish the alveolar macrophage pool in COVID-19. Early, uncommitted CD34+ hematopoietic stem/progenitor cells were primed toward megakaryopoiesis, accompanied by expanded megakaryocyte-committed progenitors and increased platelet activation. Clonally expanded CD8+ T cells and an increased ratio of CD8+ effector T cells to effector memory T cells characterized severe disease, while circulating follicular helper T cells accompanied mild disease. We observed a relative loss of IgA2 in symptomatic disease despite an overall expansion of plasmablasts and plasma cells. Our study highlights the coordinated immune response that contributes to COVID-19 pathogenesis and reveals discrete cellular components that can be targeted for therapy.

To NMD or Not To NMD: Nonsense-Mediated mRNA Decay in Cancer and Other Genetic Diseases.

The nonsense-mediated mRNA decay (NMD) pathway degrades some but not all mRNAs bearing premature termination codons (PTCs). Decades of work have elucidated the molecular mechanisms of NMD. More recently, statistical analyses of large genomic datasets have allowed the importance of known and novel 'rules of NMD' to be tested and combined into methods that accurately predict whether PTC-containing mRNAs are degraded or not. We discuss these genomic approaches and how they can be applied to identify diseases and individuals that may benefit from inhibition or activation of NMD. We also discuss the importance of NMD for gene editing and tumor evolution, and how inhibiting NMD may be an effective strategy to increase the efficacy of cancer immunotherapy.

Hydroxychloroquine Inhibits the Trained Innate Immune Response to Interferons.

Hydroxychloroquine is being investigated for a potential prophylactic effect in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, but its mechanism of action is poorly understood. Circulating leukocytes from the blood of coronavirus disease 2019 (COVID-19) patients show increased responses to Toll-like receptor ligands, suggestive of trained immunity. By analyzing interferon responses of peripheral blood mononuclear cells from healthy donors conditioned with heat-killed Candida, trained innate immunity can be modeled in vitro. In this model, hydroxychloroquine inhibits the responsiveness of these innate immune cells to virus-like stimuli and interferons. This is associated with a suppression of histone 3 lysine 27 acetylation and histone 3 lysine 4 trimethylation of inflammation-related genes, changes in the cellular lipidome, and decreased expression of interferon-stimulated genes. Our findings indicate that hydroxychloroquine inhibits trained immunity in vitro, which may not be beneficial for the antiviral innate immune response to SARS-CoV-2 infection in patients.