Single-cell immunophenotyping of the skin lesion erythema migrans identifies IgM memory B cells

Single-cell sequencing reveals distinct immune cell features in cutaneous lesions of pemphigus vulgaris and bullous pemphigoid

Bullous pemphigoid (BP) and pemphigus vulgaris (PV) are two common subtypes of autoimmune bullous disease (AIBD). The key role of circulating autoreactive immune cells contributing to skin damage of AIBD has been widely recognized. Nevertheless, the immune characteristics in cutaneous lesions remain unclear. Here, we performed single-cell RNA sequencing (scRNA-seq) and single-cell VDJ sequencing (scRNA-seq) to generate transcriptional profiles for cells and T/B cell clonetype in skin lesions of BP and PV. We found that the proportions of NK&T, macrophages/ dendritic cells, B cells, and mast cells increased in BP and PV lesions. Then, BP and PV cells constituted over 75% of all myeloid cell subtypes, CD4+ T cell subtypes and CD8+ T cell subtypes. Strikingly, CD8+ Trm was identified to be expanded in PV, and located in the intermediate state of the pseudotime trajectory from CD8+ Tm to CD8+ Tem. Interestingly, CD8+ Tem and CD4+ Treg highly expressed exhaustion-related genes, especially in BP lesions. Moreover, the enhanced cell communication between stromal cells and immune cells like B cells and macrophages/ dendritic cells was also identified in BP and PV lesions. Finally, clone expansion was observed in T cells of BP and PV compared with HC, while CD8+ Trm represented the highest ratio of hyperexpanded TCR clones among all T cell subtypes. Our study generally depicts a large and comprehensive single-cell landscape of cutaneous lesions and highlights immune cell features in BP and PV. This offers potential research targets for further investigation.

Inferring B Cell Phylogenies from Paired H and L Chain BCR Sequences with Dowser

Abs are vital to human immune responses and are composed of genetically variable H and L chains. These structures are initially expressed as BCRs. BCR diversity is shaped through somatic hypermutation and selection during immune responses. This evolutionary process produces B cell clones, cells that descend from a common ancestor but differ by mutations. Phylogenetic trees inferred from BCR sequences can reconstruct the history of mutations within a clone. Until recently, BCR sequencing technologies separated H and L chains, but advancements in single-cell sequencing now pair H and L chains from individual cells. However, it is unclear how these separate genes should be combined to infer B cell phylogenies. In this study, we investigated strategies for using paired H and L chain sequences to build phylogenetic trees. We found that incorporating L chains significantly improved tree accuracy and reproducibility across all methods tested. This improvement was greater than the difference between tree-building methods and persisted even when mixing bulk and single-cell sequencing data. However, we also found that many phylogenetic methods estimated significantly biased branch lengths when some L chains were missing, such as when mixing single-cell and bulk BCR data. This bias was eliminated using maximum likelihood methods with separate branch lengths for H and L chain gene partitions. Thus, we recommend using maximum likelihood methods with separate H and L chain partitions, especially when mixing data types. We implemented these methods in the R package Dowser: https://dowser.readthedocs.io.

Dissociation of murine oral mucosal tissues for single cell applications

Single-cell RNA sequencing and flow cytometry approaches have been instrumental in understanding cellular states within various tissues and organs. However, tissue dissociation methods can potentially alter results and create bias due to preferential recovery of particular cell types. Here we present efforts to optimize methods for dissociation of murine oral mucosal tissues and provide three different protocols that can be utilized to isolate major cell populations in the oral mucosa. These methods can be used both in health and in states of inflammation, such as periodontitis. The optimized protocols use different enzymatic approaches (collagenase II, collagenase IV and the Miltenyi whole skin dissociation kit) and yield preferential recovery of immune, stromal and epithelial cells, respectively. We suggest choosing the dissociation method based on the cell population of interest to study, while understanding the limitations of each approach.

Cellular and transcriptome signatures unveiled by single-cell RNA-Seq following ex vivo infection of murine splenocytes with Borrelia burgdorferi

Introduction

Lyme disease, the most common tick-borne infectious disease in the US, is caused by a spirochetal pathogen Borrelia burgdorferi (Bb). Distinct host responses are observed in susceptible and resistant strains of inbred of mice following infection with Bb reflecting a subset of inflammatory responses observed in human Lyme disease. The advent of post-genomic methodologies and genomic data sets enables dissecting the host responses to advance therapeutic options for limiting the pathogen transmission and/or treatment of Lyme disease.

Methods

In this study, we used single-cell RNA-Seq analysis in conjunction with mouse genomics exploiting GFP-expressing Bb to sort GFP+ splenocytes and GFP- bystander cells to uncover novel molecular and cellular signatures that contribute to early stages of immune responses against Bb.

Results

These data decoded the heterogeneity of splenic neutrophils, macrophages, NK cells, B cells, and T cells in C3H/HeN mice in response to Bb infection. Increased mRNA abundance of apoptosis-related genes was observed in neutrophils and macrophages clustered from GFP+ splenocytes. Moreover, complement-mediated phagocytosis-related genes such as C1q and Ficolin were elevated in an inflammatory macrophage subset, suggesting upregulation of these genes during the interaction of macrophages with Bb-infected neutrophils. In addition, the role of DUSP1 in regulating the expression of Casp3 and pro-inflammatory cytokines Cxcl1, Cxcl2, Il1b, and Ccl5 in Bb-infected neutrophils were identified.

Discussion

These findings serve as a growing catalog of cell phenotypes/biomarkers among murine splenocytes that can be exploited for limiting spirochetal burden to limit the transmission of the agent of Lyme disease to humans via reservoir hosts.

Bulk and single-nucleus RNA sequencing highlight immune pathways induced in individuals during an Ixodes scapularis tick bite

Ticks are hematophagous arthropods that use a complex mixture of salivary proteins to evade host defenses while taking a blood meal. Little is known about the immunological and physiological consequences of tick feeding on humans. Here, we performed the first bulk and single-nucleus RNA sequencing (snRNA-seq) of skin and blood of four persons presenting with naturally acquired, attached Ixodes scapularis ticks. Pathways and individual genes associated with innate and adaptive immunity were identified based on bulk RNA sequencing, including interleukin-17 signaling and platelet activation pathways at the site of tick attachment or in peripheral blood. snRNA-seq further revealed that the Hippo signaling, cell adhesion, and axon guidance pathways were involved in the response to an I. scapularis bite in humans. Features of the host response in these individuals also overlapped with that of laboratory guinea pigs exposed to I. scapularis and which acquired resistance to ticks. These findings offer novel insights for the development of new biomarkers for I. scapularis exposure and anti-tick vaccines for human use.

Inferring B cell phylogenies from paired heavy and light chain BCR sequences with Dowser

Antibodies are vital to human immune responses and are composed of genetically variable heavy and light chains. These structures are initially expressed as B cell receptors (BCRs). BCR diversity is shaped through somatic hypermutation and selection during immune responses. This evolutionary process produces B cell clones, cells that descend from a common ancestor but differ by mutations. Phylogenetic trees inferred from BCR sequences can reconstruct the history of mutations within a clone. Until recently, BCR sequencing technologies separated heavy and light chains, but advancements in single cell sequencing now pair heavy and light chains from individual cells. However, it is unclear how these separate genes should be combined to infer B cell phylogenies. In this study, we investigated strategies for using paired heavy and light chain sequences to build phylogenetic trees. We found incorporating light chains significantly improved tree accuracy and reproducibility across all methods tested. This improvement was greater than the difference between tree building methods and persisted even when mixing bulk and single cell sequencing data. However, we also found that many phylogenetic methods estimated significantly biased branch lengths when some light chains were missing, such as when mixing single cell and bulk BCR data. This bias was eliminated using maximum likelihood methods with separate branch lengths for heavy and light chain gene partitions. Thus, we recommend using maximum likelihood methods with separate heavy and light chain partitions, especially when mixing data types. We implemented these methods in the R package Dowser: https://dowser.readthedocs.io.

Gene set predictor for post-treatment Lyme disease

Lyme disease (LD) is tick-borne disease whose post-treatment sequelae are not well understood. For this study, we enrolled 152 individuals with symptoms of post-treatment LD (PTLD) to profile their peripheral blood mononuclear cells (PBMCs) with RNA sequencing (RNA-seq). Combined with RNA-seq data from 72 individuals with acute LD and 44 uninfected controls, we investigated differences in differential gene expression. We observe that most individuals with PTLD have an inflammatory signature that is distinguished from the acute LD group. By distilling gene sets from this study with gene sets from other sources, we identify a subset of genes that are highly expressed in the cohorts but are not already established as biomarkers for inflammatory response or other viral or bacterial infections. We further reduce this gene set by feature importance to establish an mRNA biomarker set capable of distinguishing healthy individuals from those with acute LD or PTLD as a candidate for translation into an LD diagnostic.

Tick feeding modulates the human skin immune landscape to facilitate tick-borne pathogen transmission

During cutaneous tick attachment, the feeding cavity becomes a site of transmission for tick salivary compounds and tick-borne pathogens. However, the immunological consequences of tick feeding for human skin remain unclear. Here, we assessed human skin and blood samples upon tick bite and developed a human skin explant model mimicking Ixodes ricinus bites and tick-borne pathogen infection. Following tick attachment, we observed rapidly occurring patterns of immunomodulation, including increases in neutrophils and cutaneous B and T cells. T cells upregulated tissue residency markers, while lymphocytic cytokine production was impaired. In early stages of Borrelia burgdorferi model infections, we detected strain-specific immune responses and close spatial relationships between macrophages and spirochetes. Preincubation of spirochetes with tick salivary gland extracts hampered accumulation of immune cells and increased spirochete loads. Collectively, we showed that tick feeding exerts profound changes on the skin immune network that interfere with the primary response against tick-borne pathogens.

Research progress on application of single-cell TCR/BCR sequencing technology to the tumor immune microenvironment, autoimmune diseases, and infectious diseases

Single-cell omics is the profiling of individual cells through sequencing and other technologies including high-throughput analysis for single-cell resolution, cell classification, and identification as well as time series analyses. Unlike multicellular studies, single-cell omics overcomes the problem of cellular heterogeneity. It provides new methods and perspectives for in-depth analyses of the behavior and mechanism of individual cells in the cell population and their relationship with the body, and plays an important role in basic research and precision medicine. Single-cell sequencing technologies mainly include single-cell transcriptome sequencing, single-cell assay for transposase accessible chromatin with high-throughput sequencing, single-cell immune profiling (single-cell T-cell receptor [TCR]/B-cell receptor [BCR] sequencing), and single-cell transcriptomics. Single-cell TCR/BCR sequencing can be used to obtain a large amount of single-cell gene expression and immunomics data at one time, and combined with transcriptome sequencing and TCR/BCR diversity data, can resolve immune cell heterogeneity. This paper summarizes the progress in applying single-cell TCR/BCR sequencing technology to the tumor immune microenvironment, autoimmune diseases, infectious diseases, immunotherapy, and chronic inflammatory diseases, and discusses its shortcomings and prospects for future application.

Single-cell transcriptomics in human skin research: available technologies, technical considerations, and disease applications

Single‐cell technologies have revolutionized research in the last decade, including for skin biology. Single‐cell RNA sequencing has emerged as a powerful tool allowing the dissection of human disease pathophysiology at unprecedented resolution by assessing cell‐to‐cell variation, facilitating identification of rare cell populations and elucidating cellular heterogeneity. In dermatology, this technology has been widely applied to inflammatory skin disorders, fibrotic skin diseases, wound healing complications and cutaneous neoplasms. Here, we discuss the available technologies and technical considerations of single‐cell RNA sequencing and describe its applications to a broad spectrum of dermatological diseases.