Recent advances in spatial transcriptomics (ST) have brought unprecedented opportunities to understand tissue organization and function in spatial context. However, it is still challenging to precisely dissect spatial domains with similar gene expression and histology in situ. Here, we present DeepST, an accurate and universal deep learning framework to identify spatial domains, which performs better than the existing state-of-the-art methods on benchmarking datasets of the human dorsolateral prefrontal cortex. Further testing on a breast cancer ST dataset, we showed that DeepST can dissect spatial domains in cancer tissue at a finer scale. Moreover, DeepST can achieve not only effective batch integration of ST data generated from multiple batches or different technologies, but also expandable capabilities for processing other spatial omics data. Together, our results demonstrate that DeepST has the exceptional capacity for identifying spatial domains, making it a desirable tool to gain novel insights from ST studies.

Given the chronic inflammatory nature of Parkinson’s disease (PD), T cell immunity may be important for disease onset. Here, we performed single-cell transcriptome and TCR sequencing, and conducted integrative analyses to decode composition, function and lineage relationship of T cells in the blood and cerebrospinal fluid of PD. Combined expression and TCR-based lineage tracking, we discovered a large population of CD8+ T cells showing continuous progression from central memory to terminal effector T cells in PD patients. Additionally, we identified a group of cytotoxic CD4+ T cells (CD4 CTLs) remarkably expanded in PD patients, which derived from Th1 cells by TCR-based fate decision. Finally, we screened putative TCR–antigen pairs that existed in both blood and cerebrospinal fluid of PD patients to provide potential evidence for peripheral T cells to participate in neuronal degeneration. Our study provides valuable insights and rich resources for understanding the adaptive immune response in PD.

The world is facing a pandemic of Corona Virus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Adaptive immune responses are essential for SARS-CoV-2 virus clearance. Although a large body of studies have been conducted to investigate the immune mechanism in COVID-19 patients, we still lack a comprehensive understanding of the BCR repertoire in patients. In this study, we used the single-cell V(D)J sequencing to characterize the BCR repertoire across convalescent COVID-19 patients. We observed that the BCR diversity was significantly reduced in disease compared with healthy controls. And BCRs tend to skew toward different V gene segments in COVID-19 and healthy controls. The CDR3 sequences of heavy chain in clonal BCRs in patients were more convergent than that in healthy controls. In addition, we discovered increased IgG and IgA isotypes in the disease, including IgG1, IgG3 and IgA1. In all clonal BCRs, IgG isotypes had the most frequent class switch recombination events and the highest somatic hypermutation rate, especially IgG3. Moreover, we found that an IgG3 cluster from different clonal groups had the same IGHV, IGHJ and CDR3 sequences (IGHV4-4-CARLANTNQFYDSSSYLNAMDVW-IGHJ6). Overall, our study provides a comprehensive characterization of the BCR repertoire in COVID-19 patients, which contributes to the understanding of the mechanism for the immune response to SARS-CoV-2 infection.

Immune system plays important roles in the pathogenesis of Parkinson’s disease (PD). However, the role of B cells in this complex disease are still not fully understood. B cells produce antibodies but can also regulate immune responses. In order to decode the relative contribution of peripheral B cell subtypes to the etiology of PD, we performed single cell RNA and BCR sequencing for 10,466 B cells from 8 PD patients and 6 age-matched healthy controls. We observed significant increased memory B cells and significant decreased naïve B cells in PD patients compared to healthy controls. Notably, we also discovered increased IgG and IgA isotypes and more frequent class switch recombination events in PD patients. Moreover, we identified preferential V and J gene segments of B cell receptors in PD patients as the evidence of convergent selection in PD. Finally, we found a marked clonal expanded memory B cell population in PD patients, up-regulating both MHC II genes (HLA-DRB5, HLA-DQA2 and HLA-DPB1) and transcription factor activator protein 1 (AP-1), suggesting that the antigen presentation capacity of B cells was enhanced and B cells were activated in PD patients. Overall, this study conducted a comprehensive analysis of peripheral B cell characteristics of PD patients, which provided novel insights into the humoral immune response in the pathogenesis of PD.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing an outbreak of coronavirus disease 2019 (COVID-19), has been undergoing various mutations. The analysis of the structural and energetic effects of mutations on protein-protein interactions between the receptor binding domain (RBD) of SARS-CoV-2 and angiotensin converting enzyme 2 (ACE2) or neutralizing monoclonal antibodies will be beneficial for epidemic surveillance, diagnosis, and optimization of neutralizing agents. According to the molecular dynamics simulation, a key mutation N439K in the SARS-CoV-2 RBD region created a new salt bridge which resulted in greater electrostatic complementarity. Furthermore, the N439K-mutated RBD bound hACE2 with a higher affinity than wild-type, which may lead to more infectious. In addition, the N439K-mutated RBD was markedly resistant to the SARS-CoV-2 neutralizing antibody REGN10987, which may lead to the failure of neutralization. These findings would offer guidance on the development of neutralizing antibodies and the prevention of COVID-19.