A Comprehensive Flow Cytometry Panel for Analysis of Idiopathic Subglottic Stenosis

OBJECTIVE

To present a comprehensive flow cytometry panel for idiopathic subglottic stenosis (iSGS).

STUDY DESIGN

Controlled ex vivo cohort study.

SETTING

Tertiary care academic hospital in a metropolitan area.

METHODS

Flow cytometry and single-cell RNA sequencing were performed on 9 paired normal and scar tissue samples from iSGS patients. Flow cytometry was used to assess the presence of myeloid (CD11b, CD14, CD15, Siglec8), lymphoid (CD3, CD4, CD8, gamma delta [γδ], FOXP3), endothelial (CD31), fibroblast (CD90, SMA), and epithelial (CD326, CK5) markers.

RESULTS

On flow cytometry, iSGS scar is characterized by an increased presence of myeloid, lymphoid, endothelial, and fibroblast cell types, but a decreased presence of epithelial cells. In the myeloid lineage, iSGS scar samples demonstrated increased CD11b+ monocytes (P < .001), Siglec8+ eosinophils (P = .03), and CD14+ monocytes (P = .02). In the lymphoid lineage, iSGS scar demonstrated increased CD3+ T-cells (P < .001), CD4+ helper T-cells (P < .001), γδ+ T-cells (P < .001), and FOXP3+ regulatory T-cells (P = .002). iSGS scar exhibited specific increases in CD90+ (P = .04) and SMA+ (P < .001) fibroblasts but decreased CD326+ (E-cadherin) epithelial cells (P = .01) relative to normal samples.

CONCLUSION

We present a comprehensive flow cytometry panel for iSGS. This flow panel may serve as a common platform among airway scientists to elucidate the cellular mechanisms underpinning iSGS and other upper airway pathologies. Scar iSGS samples demonstrate a distinct cellular profile relative to normal iSGS specimens, exhibiting increased fibroblast, endothelial, and inflammatory cell types but decreased epithelium.

Immune cell-related prognostic risk model and tumor immune environment modulation in esophageal carcinoma based on single-cell and bulk RNA sequencing

BACKGROUND

Immune cells play a pivotal role in the tumor microenvironment, exerting significant influence on tumor progression and patient outcomes, but the current biomarkers are insufficient to fully capture the complex and diverse tumor immune microenvironment and the impact of immunotherapy.

METHODS

The advent of single-cell sequencing allows us to explore the tumor microenvironment at an unprecedented resolution, enabling the identification and characterization of distinct subsets of immune cells, thereby paving the way for the development of prognostic models using immune cells. Leveraging single-cell data, our study deeply investigated the intricacies of immune microenvironment heterogeneity in esophageal carcinoma.

RESULTS

We elucidated the composition, functionality, evolution, and intercellular communication patterns of immune cells, culminating in the construction of an independent prognostic model at the single-cell level. Furthermore, we conducted a comprehensive analysis of disparities in immune infiltration and immune checkpoint expression between patients categorized into high- and low-risk groups, which may impact patient prognosis.

CONCLUSION

In summary, our study harnessed multiomics data to delineate the immune profile of esophageal carcinoma patients, provide a method for leveraging molecular signatures of immune cells to identify potential biomarkers, while concurrently providing evidence for the potential benefits of immunotherapy.

Pan-Cancer Single-Cell and Spatial-Resolved Profiling Reveals the Immunosuppressive Role of APOE+ Macrophages in Immune Checkpoint Inhibitor Therapy

The heterogeneity of macrophages influences the response to immune checkpoint inhibitor (ICI) therapy. However, few studies explore the impact of APOE+ macrophages on ICI therapy using single-cell RNA sequencing (scRNA-seq) and machine learning methods. The scRNA-seq and bulk RNA-seq data are Integrated to construct an M.Sig model for predicting ICI response based on the distinct molecular signatures of macrophage and machine learning algorithms. Comprehensive single-cell analysis as well as in vivo and in vitro experiments are applied to explore the potential mechanisms of the APOE+ macrophage in affecting ICI response. The M.Sig model shows clear advantages in predicting the efficacy and prognosis of ICI therapy in pan-cancer patients. The proportion of APOE+ macrophages is higher in ICI non-responders of triple-negative breast cancer compared with responders, and the interaction and longer distance between APOE+ macrophages and CD8+ exhausted T (Tex) cells affecting ICI response is confirmed by multiplex immunohistochemistry. In a 4T1 tumor-bearing mice model, the APOE inhibitor combined with ICI treatment shows the best efficacy. The M.Sig model using real-world immunotherapy data accurately predicts the ICI response of pan-cancer, which may be associated with the interaction between APOE+ macrophages and CD8+ Tex cells.

Immune features revealed by single-cell RNA and single-cell TCR/BCR sequencing in patients with rheumatoid arthritis receiving COVID-19 booster vaccination

Coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, have profoundly affected human health. Booster COVID-19 vaccines have demonstrated significant efficacy in reducing infection and severe cases. However, the effects of booster COVID-19 vaccines on key immune cell subsets and their responses in rheumatoid arthritis (RA) are not well understood. By using single-cell RNA sequencing (scRNA-seq) combined with scTCR/BCR-seq analysis, a total of 8 major and 27 minor cell clusters were identified from paired peripheral blood mononuclear cells (PBMCs) which were collected 1 week before and 4 weeks after booster vaccination in stable RA patients. Booster vaccination only had limited impact on the composition and proportions of PBMCs cell clusters. CD8+ cytotoxic T cells (CD8+T_CTL) showed a trend toward an increase after vaccination, while naive B cells and conventional dendritic cells (cDCs) showed a trend toward a decrease. Transcriptomic changes were observed after booster vaccination, primarily involving T/B cell receptor signaling pathways, phagosome, antigen processing and presenting, and viral myocarditis pathways. Interferon (IFN) and pro-inflammatory response gene sets were slightly upregulated across most major cell subpopulations in COVID-19 booster-vaccinated RA individuals. Plasma neutralizing antibody titers significantly increased after booster COVID-19 vaccination (p = 0.037). Single-cell TCR/BCR analysis revealed increased B cell clone expansion and repertoire diversity postvaccination, with no consistent alterations in T cells. Several clonotypes of BCRs and TCRs were identified to be significantly over-represented after vaccination, such as IGHV3-15 and TRBV28. Our study provided a comprehensive single-cell atlas of the peripheral immune response and TCR/BCR immune repertoire profiles to inactivated SARS-CoV-2 booster vaccination in RA patients, which helps us to understand vaccine-induced immune responses better.

Construction of a prognostic model with CAFs for predicting the prognosis and immunotherapeutic response of lung squamous cell carcinoma

Lung squamous cell carcinoma (LUSC) is one of the subtypes of lung cancer (LC) that contributes to approximately 25%-30% of its prevalence. Cancer-associated fibroblasts (CAFs) are key cellular components of the TME, and the large number of CAFs in tumour tissues creates a favourable environment for tumour development. However, the function of CAFs in the LUSC is complex and uncertain. First, we processed the scRNA-seq data and classified distinct types of CAFs. We also identified prognostic CAFRGs using univariate Cox analysis and conducted survival analysis. Additionally, we assessed immune cell infiltration in CAF clusters using ssGSEA. We developed a model with a significant prognostic correlation and verified the prognostic model. Furthermore, we explored the immune landscape of LUSC and further investigated the correlation between malignant features and LUSC. We identified CAFs and classified them into three categories: iCAFs, mCAFs and apCAFs. The survival analysis showed a significant correlation between apCAFs and iCAFs and LUSC patient prognosis. Kaplan-Meier analysis showed that patients in CAF cluster C showed a better survival probability compared to clusters A and B. In addition, we identified nine significant prognostic CAFRGs (CLDN1, TMX4, ALPL, PTX3, BHLHE40, TNFRSF12A, VKORC1, CST3 and ADD3) and subsequently employed multivariate Cox analysis to develop a signature and validate the model. Lastly, the correlation between CAFRG and malignant features indicates the potential role of CAFRG in promoting tumour angiogenesis, EMT and cell cycle alterations. We constructed a CAF prognostic signature for identifying potential prognostic CAFRGs and predicting the prognosis and immunotherapeutic response for LUSC. Our study may provide a more accurate prognostic assessment and immunotherapy targeting strategies for LUSC.

Single-cell transcriptomics reveal metastatic CLDN4+ cancer cells underlying the recurrence of malignant pleural effusion in patients with advanced non-small-cell lung cancer

BACKGROUND

Recurrent malignant pleural effusion (MPE) resulting from non-small-cell lung cancer (NSCLC) is easily refractory to conventional therapeutics and lacks predictive markers. The cellular or genetic signatures of recurrent MPE still remain largely uncertain.

METHODS

16 NSCLC patients with pleural effusions were recruited, followed by corresponding treatments based on primary tumours. Non-recurrent or recurrent MPE was determined after 3-6 weeks of treatments. The status of MPE was verified by computer tomography (CT) and cytopathology, and the baseline pleural fluids were collected for single-cell RNA sequencing (scRNA-seq). Samples were then integrated and profiled. Cellular communications and trajectories were inferred by bioinformatic algorithms. Comparative analysis was conducted and the results were further validated by quantitative polymerase chain reaction (qPCR) in a larger MPE cohort from the authors' centre (n = 64).

RESULTS

The scRNA-seq revealed that 33 590 cells were annotated as 7 major cell types and further characterized into 14 cell clusters precisely. The cell cluster C1, classified as Epithelial Cell Adhesion Molecule (EpCAM)+ metastatic cancer cell and correlated with activation of tight junction and adherence junction, was significantly enriched in the recurrent MPE group, in which Claudin-4 (CLDN4) was identified. The subset cell cluster C3 of C1, which was enriched in recurrent MPE and demonstrated a phenotype of ameboidal-type cell migration, also showed a markedly higher expression of CLDN4. Meanwhile, the expression of CLDN4 was positively correlated with E74 Like ETS Transcription Factor 3 (ELF3), EpCAM and Tumour Associated Calcium Signal Transducer 2 (TACSTD2), independent of driver-gene status. CLDN4 was also found to be associated with the expression of Hypoxia Inducible Factor 1 Subunit Alpha (HIF1A) and Vascular Endothelial Growth Factor A (VEGFA), and the cell cluster C1 was the major mediator in cellular communication of VEGFA signalling. In the extensive MPE cohort, a notably increased expression of CLDN4 in cells from pleural effusion among patients diagnosed with recurrent MPE was observed, compared with the non-recurrent group, which was also associated with a trend towards worse overall survival (OS).

CONCLUSIONS

CLDN4 could be considered as a predictive marker of recurrent MPE among patients with advanced NSCLC. Further validation for its clinical value in cohorts with larger sample size and in-depth mechanism studies on its biological function are warranted.

TRIAL REGISTRATION

Not applicable.

Integrating Single-Cell and Spatial Transcriptomics Reveals Heterogeneity of Early Pig Skin Development and a Subpopulation with Hair Placode Formation

The dermis and epidermis, crucial structural layers of the skin, encompass appendages, hair follicles (HFs), and intricate cellular heterogeneity. However, an integrated spatiotemporal transcriptomic atlas of embryonic skin has not yet been described and would be invaluable for studying skin-related diseases in humans. Here, single-cell and spatial transcriptomic analyses are performed on skin samples of normal and hairless fetal pigs across four developmental periods. The cross-species comparison of skin cells illustrated that the pig epidermis is more representative of the human epidermis than mice epidermis. Moreover, Phenome-wide association study analysis revealed that the conserved genes between pigs and humans are strongly associated with human skin-related diseases. In the epidermis, two lineage differentiation trajectories describe hair follicle (HF) morphogenesis and epidermal development. By comparing normal and hairless fetal pigs, it is found that the hair placode (Pc), the most characteristic initial structure in HFs, arises from progenitor-like OGN+/UCHL1+ cells. These progenitors appear earlier in development than the previously described early Pc cells and exhibit abnormal proliferation and migration during differentiation in hairless pigs. The study provides a valuable resource for in-depth insights into HF development, which may serve as a key reference atlas for studying human skin disease etiology using porcine models.

Integration of single-cell and RNA-seq data to explore the role of focal adhesion-related genes in osteoporosis

Integrin-based focal adhesion is one of the major mechanosensory in osteocytes. The aim of this study was to mine the hub genes associated with focal adhesion and investigate their roles in osteoporosis based on the data of single-cell RNA sequencing and RNA-sequencing. Two hub genes (FAM129A and RNF24) with the same expression trend and AUC values greater than 0.7 in both GSE56815 and GSE56116 cohorts were uncovered. The nomogram was created to predict the risk of OP based on two hub genes. Subsequently, the competing endogenous RNA network was established based on two hub genes, 14 microRNAs and five long noncoding RNAs. Meanwhile, transcription factors-hub gene network was established based on two hub genes and 14 TFs. Finally, 73 drugs were predicted, of which there were 13 drugs targeting FAM129A and 66 drugs targeting RNF24. In both mouse and human blood samples, FAM129A expression was decreased in granulocytes and RNF24 expression was increased in monocytes. In the mouse experiment, FAM129A and anti-RNF24 were found to partially alleviate the progression of osteoporosis. In conclusion, two hub genes related to focal adhesion were identified by combined scRNA-seq and RNA-seq analyses, which might supply a new insight for the treatment and evaluation of OP.

Distribution-Agnostic Deep Learning Enables Accurate Single-Cell Data Recovery and Transcriptional Regulation Interpretation

Single-cell RNA sequencing (scRNA-seq) is a robust method for studying gene expression at the single-cell level, but accurately quantifying genetic material is often hindered by limited mRNA capture, resulting in many missing expression values. Existing imputation methods rely on strict data assumptions, limiting their broader application, and lack reliable supervision, leading to biased signal recovery. To address these challenges, authors developed Bis, a distribution-agnostic deep learning model for accurately recovering missing sing-cell gene expression from multiple platforms. Bis is an optimal transport-based autoencoder model that can capture the intricate distribution of scRNA-seq data while addressing the characteristic sparsity by regularizing the cellular embedding space. Additionally, they propose a module using bulk RNA-seq data to guide reconstruction and ensure expression consistency. Experimental results show Bis outperforms other models across simulated and real datasets, showcasing superiority in various downstream analyses including batch effect removal, clustering, differential expression analysis, and trajectory inference. Moreover, Bis successfully restores gene expression levels in rare cell subsets in a tumor-matched peripheral blood dataset, revealing developmental characteristics of cytokine-induced natural killer cells within a head and neck squamous cell carcinoma microenvironment.

Comparative single-cell RNA sequencing analysis of immune response to inactivated vaccine and natural SARS-CoV-2 infection

Uncovering the immune response to an inactivated SARS-CoV-2 vaccine (In-Vac) and natural infection is crucial for comprehending COVID-19 immunology. Here we conducted an integrated analysis of single-cell RNA sequencing (scRNA-seq) data from serial peripheral blood mononuclear cell (PBMC) samples derived from 12 individuals receiving In-Vac compared with those from COVID-19 patients. Our study reveals that In-Vac induces subtle immunological changes in PBMC, including cell proportions and transcriptomes, compared with profound changes for natural infection. In-Vac modestly upregulates IFN-α but downregulates NF-κB pathways, while natural infection triggers hyperactive IFN-α and NF-κB pathways. Both In-Vac and natural infection alter T/B cell receptor repertoires, but COVID-19 has more significant change in preferential VJ gene, indicating a vigorous immune response. Our study reveals distinct patterns of cellular communications, including a selective activation of IL-15RA/IL-15 receptor pathway after In-Vac boost, suggesting its potential role in enhancing In-Vac-induced immunity. Collectively, our study illuminates multifaceted immune responses to In-Vac and natural infection, providing insights for optimizing SARS-CoV-2 vaccine efficacy.

An ultra high-throughput, massively multiplexable, single-cell RNA-seq platform in yeasts

Yeasts are naturally diverse, genetically tractable, and easy to grow such that researchers can investigate any number of genotypes, environments, or interactions thereof. However, studies of yeast transcriptomes have been limited by the processing capabilities of traditional RNA sequencing techniques. Here we optimize a powerful, high-throughput single-cell RNA sequencing (scRNAseq) platform, SPLiT-seq (Split Pool Ligation-based Transcriptome sequencing), for yeasts and apply it to 43,388 cells of multiple species and ploidies. This platform utilizes a combinatorial barcoding strategy to enable massively parallel RNA sequencing of hundreds of yeast genotypes or growth conditions at once. This method can be applied to most species or strains of yeast for a fraction of the cost of traditional scRNAseq approaches. Thus, our technology permits researchers to leverage "the awesome power of yeast" by allowing us to survey the transcriptome of hundreds of strains and environments in a short period of time and with no specialized equipment. The key to this method is that sequential barcodes are probabilistically appended to cDNA copies of RNA while the molecules remain trapped inside of each cell. Thus, the transcriptome of each cell is labeled with a unique combination of barcodes. Since SPLiT-seq uses the cell membrane as a container for this reaction, many cells can be processed together without the need to physically isolate them from one another in separate wells or droplets. Further, the first barcode in the sequence can be chosen intentionally to identify samples from different environments or genetic backgrounds, enabling multiplexing of hundreds of unique perturbations in a single experiment. In addition to greater multiplexing capabilities, our method also facilitates a deeper investigation of biological heterogeneity, given its single-cell nature. For example, in the data presented here, we detect transcriptionally distinct cell states related to cell cycle, ploidy, metabolic strategies, and so forth, all within clonal yeast populations grown in the same environment. Hence, our technology has two obvious and impactful applications for yeast research: the first is the general study of transcriptional phenotypes across many strains and environments, and the second is investigating cell-to-cell heterogeneity across the entire transcriptome.

Deciphering the cellular heterogeneity of the insect brain with single-cell RNA sequencing

Insects show highly complicated adaptive and sophisticated behaviors, including spatial orientation skills, learning ability, and social interaction. These behaviors are controlled by the insect brain, the central part of the nervous system. The tiny insect brain consists of millions of highly differentiated and interconnected cells forming a complex network. Decades of research has gone into an understanding of which parts of the insect brain possess particular behaviors, but exactly how they modulate these functional consequences needs to be clarified. Detailed description of the brain and behavior is required to decipher the complexity of cell types, as well as their connectivity and function. Single-cell RNA-sequencing (scRNA-seq) has emerged recently as a breakthrough technology to understand the transcriptome at cellular resolution. With scRNA-seq, it is possible to uncover the cellular heterogeneity of brain cells and elucidate their specific functions and state. In this review, we first review the basic structure of insect brains and the links to insect behaviors mainly focusing on learning and memory. Then the scRNA applications on insect brains are introduced by representative studies. Single-cell RNA-seq has allowed researchers to classify cell subpopulations within different insect brain regions, pinpoint single-cell developmental trajectories, and identify gene regulatory networks. These developments empower the advances in neuroscience and shed light on the intricate problems in understanding insect brain functions and behaviors.

Molecular mechanism analyses of post-traumatic epilepsy and hereditary epilepsy based on 10× single-cell transcriptome sequencing technology

BACKGROUND

Single-cell RNA sequencing analysis has been usually conducted on post-traumatic epilepsy (PET) and hereditary epilepsy (HE) patients; however, the transcriptome of patients with traumatic temporal lobe epilepsy has rarely been studied.

MATERIALS AND METHODS

Hippocampus tissues isolated from one patient with PTE and one patient with HE were used in the present study. Single cell isolates were prepared and captured using a 10× Genomics Chromium Single-Cell 3' kit (V3) according to the manufacturer's instructions. The libraries were sequenced on an Illumina NovaSeq 6000 sequencing system. Raw data were processed, and the cells were filtered and classified using the Seurat R package. Uniform Manifold Approximation and Projection was used for visualization. Differentially expressed genes (DEGs) were identified based on a p-value ≤0.01 and log fold change (FC) ≥0.25. Gene Ontology (GO, http://geneontology.org/) and KEGG (Kyoto Encyclopedia of Genes and Genomes, www.genome.jp/kegg) analyses were performed on the DEGs for enrichment analysis.

RESULTS

The reads obtained from the 10× genomic platform for PTE and HE were 39.56 M and 30.08 M, respectively. The Q30 score of the RNA reads was >91.6%. After filtering, 7479 PTE cells and 9357 HE cells remained for further study. More than 96.4% of the reads were mapped to GRCh38/GRCm38. The cells were differentially distributed in two groups, with higher numbers of oligodendrocytes (6522 vs. 2532) and astrocytes (133 vs. 52), and lower numbers of microglial cells (2242 vs. 3811), and neurons (3 vs. 203) present in the HE group than in the PTE group. The DEGs in four cell clusters were identified, with 25 being in oligodendrocytes (13 upregulated and 12 downregulated), 87 in microglia cells (42 upregulated and 45 downregulated), 222 in astrocytes (115 upregulated and 107 downregulated), and 393 in neurons (305 upregulated and 88 downregulated). The genes MTND1P23 (downregulated), XIST (downregulated), and RPS4Y1 (upregulated) were commonly expressed in all four cell clusters. The DEGs in microglial cells and astrocytes were enriched in the IL-17 signaling pathway.

CONCLUSION

Our study explored differences in cells found in a patient with PE compared to a patient with HE, and the transcriptome in the different cells was analyzed for the first time. Studying inflammatory and immune functions might be the best approach for investigating traumatic temporal lobe epilepsy in neurons.

Itraconazole inhibits tumor growth via CEBPB-mediated glycolysis in colorectal cancer

Advanced colorectal cancer (CRC) is characterized by a high recurrence and metastasis rate, which is the primary cause of patient mortality. Unfortunately, effective anti-cancer drugs for CRC are still lacking in clinical practice. We screened FDA-approved drugs by utilizing targeted organoid sequencing data and found that the antifungal drug itraconazole had a potential therapeutic effect on CRC tumors. However, the effect and mechanism of itraconazole on CRC tumors have not been investigated. A cell line-derived xenograft model in tumor-bearing mice was established and single-cell RNA sequencing was performed on tumor samples from four mice with or without itraconazole treatment. The proportion of cell populations and gene expression profiles was significantly different between the two groups. We found that itraconazole could inhibit tumor growth and glycolysis. We revealed that CEBPB was a new target for itraconazole, and that silencing CEBPB could repress CRC glycolysis and tumor growth by inhibiting ENO1 expression. Clinical analysis showed that CEBPB expression was obviously elevated in CRC patients, and was associated with poor survival. In summary, itraconazole treatment remodeled cell composition and gene expression profiles. Itraconazole inhibited cell glycolysis and tumor growth via the CEBPB-ENO1 axis. In this study, we illustrate a new energy metabolism mechanism for itraconazole on tumor growth in CRC that will provide a theoretical basis for CRC targeting/combination therapy.

Novel subsets of peripheral immune cells associated with promoting stroke recovery in mice

AIMS

Peripheral immune cells infiltrating into the brain trigger neuroinflammation after an ischemic stroke. Partial immune cells reprogram their function for neural repair. Which immune cells promote ischemic brain recovery needs further identification.

METHODS

We performed single-cell transcriptomic profiling of CD45high immune cells isolated from the ischemic hemisphere at subacute (5 days) and chronic (14 days) stages after ischemic stroke.

RESULTS

A subset of phagocytic macrophages was associated with neuron projection regeneration and tissue remodeling. We also identified a unique type of T cells with highly expressed macrophage markers, including C1q, Apoe, Hexb, and Fcer1g, which showed high abilities in tissue remodeling, myelination regulation, wound healing, and anti-neuroinflammation. Moreover, natural killer cells decreased cytotoxicity and increased energy and metabolic function in the chronic stage after ischemic stroke. Two subgroups of neutrophils upregulated CCL signals to recruit peripheral immune cells and released CXCL2 to keep self-recruiting at the chronic stage.

CONCLUSIONS

We identified subsets of peripheral immune cells that may provide potential therapeutic targets for promoting poststroke recovery.

Single-cell transcriptomics reveal the remodeling landscape of bladder in patients with obstruction-induced detrusor underactivity

Detrusor underactivity (DUA) is a common and thorny problem in urology, which severely impairs patients' bladder function and quality of life. However, its underlying pathophysiological mechanism remains unclear. Hence, we sequenced 69,973 cells from five controls and nine patients with bladder dysfunction using single-cell RNA sequencing. Twelve distinct cell types were identified and they showed high cellular and functional heterogeneity among each group. Among them, fibroblasts, macrophages, and epithelial cells had the most intercellular communications. Their aberrant gene expressions and altered intercellular interactions were mainly involved in extracellular matrix organization, inflammation/immune regulation, and cellular injury. Further re-cluster analysis revealed an accumulation of the RBFOX1+ fibroblasts and RIPOR2+ macrophages in dysfunctional bladder wall, which mediated bladder remodeling through dysfunctional extracellular matrix organization and inflammation/immune reaction. Besides, the subtype of the epithelial cells was significantly altered. They underwent an intricate process including inflammation, damage, and repair during bladder remodeling. Overall, this work constructed the first single-cell atlas for obstruction-induced DUA, which could provide a valuable resource for deciphering the cellular heterogeneity and function changes in DUA, as well as potential strategies for bladder function improvement.

Single-cell sequencing reveals the immune landscape of breast cancer patients with brain metastasis

BACKGROUND

Breast cancer has the highest incidence rate of cancer worldwide, and brain metastases (BrM) are among the most malignant cases. While some patients have benefited from immune checkpoint inhibitors (ICIs), the complex anatomical structure of the brain and the heterogeneity of metastatic tumors have made it difficult to characterize the tumor immune microenvironment (TME) of metastatic tumors.

METHODS

To address this, we used single-cell RNA sequencing (scRNA-seq) to analyze immune cells in the cerebrospinal fluid (CSF) of BrM patients with breast cancer, thereby providing a comprehensive view of the immune microenvironment landscape of BrM.

RESULTS

Based on canonical marker genes, we identified nine cell types, and further identified their subtypes through differential expression gene (DEG) analysis. We compared the changes in cells and functions in the immune microenvironment of patients with different prognoses. Our analysis revealed a series of genes that promote tumor immune function (CCR5, LYZ, IGKC, MS4A1, etc.) and inhibit tumor immune function (SCGB2A2, CD24, etc.).

CONCLUSIONS

The scRNA-seq in CSF provides a noninvasive method to describe the TME of breast cancer patients and guide immunotherapy.

Exploration of the prognostic prediction value of the PANoptosis-based risk score and its correlation with tumor immunity in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a common cancer with high mortality worldwide. PANoptosis is a novel inflammatory programmed cell death modality with the characteristics of pyroptosis, apoptosis and necroptosis. It is necessary to explore PANoptosis-related genes in LUAD patients and offer evidence for prognosis prediction and therapeutic strategies. Single-cell RNA sequencing data and RNA expression profiles of LUAD patients from The Cancer Genome Atlas and Gene Expression Omnibus databases are used to screen PANoptosis-related differential genes for the construction of a risk model. Fifteen PANoptosis-related markers with prognostic value were identified by Least Absolute Shrinkage and Selection Operator (LASSO)-Cox regression analysis. Kaplan-Meier analysis and receiver operating characteristic curve analysis further demonstrated the significant predictive capability. Immune infiltration, Single Nucleotide Variants (SNV) mutations, and clinical drug susceptibility were analyzed. In conclusion, a risk model of 15 PANoptosis-related genes has significant value in prognostic prediction for LUAD and has potential to direct clinical therapeutic strategies during the treatment.

Comparable CD8+ T-cell responses to SARS-CoV-2 vaccination in single-cell transcriptomics of recently allogeneic transplanted patients and healthy individuals

Despite extensive research on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccination responses in healthy individuals, there is comparatively little known beyond antibody titers and T-cell responses in the vulnerable cohort of patients after allogeneic hematopoietic stem cell transplantation (ASCT). In this study, we assessed the serological response and performed longitudinal multimodal analyses including T-cell functionality and single-cell RNA sequencing combined with T cell receptor (TCR)/B cell receptor (BCR) profiling in the context of BNT162b2 vaccination in ASCT patients. In addition, these data were compared to publicly available data sets of healthy vaccinees. Protective antibody titers were achieved in 40% of patients. We identified a distorted B- and T-cell distribution, a reduced TCR diversity, and increased levels of exhaustion marker expression as possible causes for the poorer vaccine response rates in ASCT patients. Immunoglobulin heavy chain gene rearrangement after vaccination proved to be highly variable in ASCT patients. Changes in TCRα and TCRβ gene rearrangement after vaccination differed from patterns observed in healthy vaccinees. Crucially, ASCT patients elicited comparable proportions of SARS-CoV-2 vaccine-induced (VI) CD8+ T-cells, characterized by a distinct gene expression pattern that is associated with SARS-CoV-2 specificity in healthy individuals. Our study underlines the impaired immune system and thus the lower vaccine response rates in ASCT patients. However, since protective vaccine responses and VI CD8+ T-cells can be induced in part of ASCT patients, our data advocate early posttransplant vaccination due to the high risk of infection in this vulnerable group.

Single-Cell Atlas of Neonatal Mouse Hearts Reveals an Unexpected Cardiomyocyte

BACKGROUND

Single-cell RNA sequencing is widely used in cancer research and organ development because of its powerful ability to analyze cellular heterogeneity. However, its application in cardiomyocytes is dissatisfactory mainly because the cardiomyocytes are too large and fragile to withstand traditional single-cell approaches.

METHODS AND RESULTS

Through designing the isolation procedure of neonatal mouse cardiac cells, we provide detailed cellular atlases of the heart at single-cell resolution across 4 different stages after birth. We have obtained 10 000 cardiomyocytes; to our knowledge, this is the most extensive reference framework to date. Moreover, we have discovered unexpected erythrocyte-like cardiomyocyte-terminal cardiomyocytes, comprising more than a third of all cardiomyocytes. Only a few genes are highly expressed in these cardiomyocytes. They are highly differentiated cardiomyocytes that function as contraction pumps. In addition, we have identified 2 cardiomyocyte-like conducting cells, lending support to the theory that the sinoatrial node pacemaker cells are specialized cardiomyocytes. Notably, we provide an initial blueprint for comprehensive interactions between cardiomyocytes and other cardiac cells.

CONCLUSIONS

This mouse cardiac cell atlas improves our understanding of cardiomyocyte heterogeneity and provides a valuable reference in response to varying physiological conditions and diseases.

Identification of RP11-770J1.4 as immune-related lncRNA regulating the CTXN1-cGAS-STING axis in histologically lower-grade glioma

Human gliomas are lethal brain cancers. Emerging evidence revealed the regulatory role of long noncoding RNAs (lncRNAs) in tumors. Here, we performed a comprehensive analysis of the expression profiles of RNAs in histologically lower-grade glioma (LGG). Enrichment analysis revealed that glioma is influenced by immune-related signatures. Survival analysis further established the close correlation between network features and glioma prognosis. Subsequent experiments showed lncRNA RP11-770J1.4 regulates CTXN1 expression through hsa-miR-124-3p. Correlation analysis identified lncRNA RP11-770J1.4 was immune related, specifically involved in the cytosolic DNA sensing pathway. Downregulated lncRNA RP11-770J1.4 resulted in increased spontaneous gene expression of the cGAS-STING pathway. Single-cell RNA sequencing analysis, along with investigations in a glioblastoma stem cell model and patient sample analysis, demonstrated the predominant localization of CTXN1 within tumor cores rather than peripheral regions. Immunohistochemistry staining established a negative correlation between CTXN1 expression and infiltration of CD8+ T cells. In vivo, Ctxn1 knockdown in GL261 cells led to decreased tumor burden and improved survival while increasing infiltration of CD8+ T cells. These findings unveil novel insights into the lncRNA RP11-770J1.4-CTXN1 as a potential immune regulatory axis, highlighting its therapeutic implications for histologically LGGs.

Single-cell RNA sequencing of retina revealed novel transcriptional landscape in high myopia and underlying cell-type-specific mechanisms

High myopia is a leading cause of blindness worldwide with increasing prevalence. Retina percepts visual information and triggers myopia development, but the underlying etiology is not fully understood because of cellular heterogeneity. In this study, single-cell RNA sequencing analysis was performed on retinas of mouse highly myopic and control eyes to dissect the involvement of each cell type during high myopia progression. For highly myopic photoreceptors, Hk2 inhibition underlying metabolic remodeling from aerobic glycolysis toward oxidative phosphorylation and excessive oxidative stress was identified. Importantly, a novel Apoe + rod subpopulation was specifically identified in highly myopic retina. In retinal neurons of highly myopic eyes, neurodegeneration was generally discovered, and the imbalanced ON/OFF signaling driven by cone-bipolar cells and the downregulated dopamine receptors in amacrine cells were among the most predominant findings, indicating the aberrant light processing in highly myopic eyes. Besides, microglia exhibited elevated expression of cytokines and TGF-β receptors, suggesting enhanced responses to inflammation and the growth-promoting states involved in high myopia progression. Furthermore, cell-cell communication network revealed attenuated neuronal interactions and increased glial/vascular interactions in highly myopic retinas. In conclusion, this study outlines the transcriptional landscape of highly myopic retina, providing novel insights into high myopia development and prevention.

Autoimmune uveitis in Behçet's disease and Vogt-Koyanagi-Harada disease differ in tissue immune infiltration and T cell clonality

Objectives

Non-infectious uveitis is often secondary to systemic autoimmune diseases, with Behçet's disease (BD) and Vogt-Koyanagi-Harada disease (VKHD) as the two most common causes. Uveitis in BD and VKHD can show similar clinical manifestations, but the underlying immunopathogenesis remains unclear.

Methods

To understand immune landscapes in inflammatory eye tissues, we performed single-cell RNA paired with T cell receptor (TCR) sequencing of immune cell infiltrates in aqueous humour from six patients with BD (N = 3) and VKHD (N = 3) uveitis patients.

Results

Although T cells strongly infiltrated in both types of autoimmune uveitis, myeloid cells only significantly presented in BD uveitis but not in VKHD uveitis. Conversely, VKHD uveitis but not BD uveitis showed an overwhelming dominance by CD4+ T cells (> 80%) within the T cell population due to expansion of CD4+ T cell clusters with effector memory (Tem) phenotypes. Correspondingly, VKHD uveitis demonstrated a selective expansion of CD4+ T cell clones which were enriched in pro-inflammatory Granzyme H+ CD4+ Tem cluster and showed TCR and Th1 pathway activation. In contrast, BD uveitis showed a preferential expansion of CD8+ T cell clones in pro-inflammatory Granzyme H+ CD8+ Tem cluster, and pathway activation for cytoskeleton remodelling, cellular adhesion and cytotoxicity.

Conclusion

Single-cell analyses of ocular tissues reveal distinct landscapes of immune cell infiltration and T-cell clonal expansions between VKHD and BD uveitis. Preferential involvements of pro-inflammatory CD4+ Th1 cells in VKHD and cytotoxic CD8+ T cells in BD suggest a difference in disease immunopathogenesis and can guide precision disease management.

Spatiotemporal single-cell transcriptomic profiling reveals inflammatory cell states in a mouse model of diffuse alveolar damage

Diffuse alveolar damage (DAD) triggers neutrophilic inflammation in damaged tissues of the lung, but little is known about the distinct roles of tissue structural cells in modulating the recruitment of neutrophils to damaged areas. Here, by combining single-cell and spatial transcriptomics, and using quantitative assays, we systematically analyze inflammatory cell states in a mouse model of DAD-induced neutrophilic inflammation after aerosolized intratracheal inoculation with ricin toxin. We show that homeostatic resident fibroblasts switch to a hyper-inflammatory state, and the subsequent occurrence of a CXCL1-CXCR2 chemokine axis between activated fibroblasts (AFib) as the signal sender and neutrophils as the signal receiver triggers further neutrophil recruitment. We also identify an anatomically localized inflamed niche (characterized by a close-knit spatial intercellular contact between recruited neutrophils and AFib) in peribronchial regions that facilitate the pulmonary inflammation outbreak. Our findings identify an intricate interplay between hyper-inflammatory fibroblasts and neutrophils and provide an overarching profile of dynamically changing inflammatory microenvironments during DAD progression.

Overview of single-cell RNA sequencing analysis and its application to spermatogenesis research

Background

Single-cell transcriptomics allows parallel analysis of multiple cell types in tissues. Because testes comprise somatic cells and germ cells at various stages of spermatogenesis, single-cell RNA sequencing is a powerful tool for investigating the complex process of spermatogenesis. However, single-cell RNA sequencing analysis needs extensive knowledge of experimental technologies and bioinformatics, making it difficult for many, particularly experimental biologists and clinicians, to use it.

Methods

Aiming to make single-cell RNA sequencing analysis familiar, this review article presents an overview of experimental and computational methods for single-cell RNA sequencing analysis with a history of transcriptomics. In addition, combining the PubMed search and manual curation, this review also provides a summary of recent novel insights into human and mouse spermatogenesis obtained using single-cell RNA sequencing analyses.

Main Findings

Single-cell RNA sequencing identified mesenchymal cells and type II innate lymphoid cells as novel testicular cell types in the adult mouse testes, as well as detailed subtypes of germ cells. This review outlines recent discoveries into germ cell development and subtypes, somatic cell development, and cell-cell interactions.

Conclusion

The findings on spermatogenesis obtained using single-cell RNA sequencing may contribute to a deeper understanding of spermatogenesis and provide new directions for male fertility therapy.

Global Characteristics and Dynamics of Single Immune Cells After Myocardial Infarction

Background Myocardial infarction (MI) is characterized by the emergence of dead or dying cardiomyocytes and excessive immune cell infiltration after coronary vessel occlusion. However, the complex transcriptional profile, pathways, cellular interactome, and transcriptional regulators of immune subpopulations after MI remain elusive. Methods and Results Here, male C57BL/6 mice were subjected to MI surgery and monitored for 1 day and 7 days, or sham surgery for 7 days, then cardiac CD45-positive immune cells were collected for single-cell RNA sequencing to determine immune heterogeneity. A total of 30 135 CD45⁺ immune cells were partitioned into macrophages, monocytes, neutrophils, dendritic cells, and T or B cells for further analysis. We showed that macrophages enriched for Olr1 and differentially expressed Gpnmb represented 2 crucial ischemia-associated macrophages with distinct proinflammatory and prophagocytic capabilities. In contrast to the proinflammatory subset of macrophages enriched for Olr1, Gpnmb-positive macrophages exhibited higher phagocytosis and fatty acid oxidation preference, which could be abolished by etomoxir treatment. In addition to macrophages, MI triggered prompt recruitment of neutrophils into murine hearts, which constituted the sequential cell-fate from naïve S100a4-positive, to activated Sell-high, to aging Icam1-high neutrophils. In silico tools predicted that the excessively expanded neutrophils at 1 day were attributed to chemokine C-C motif ligand/chemokine C-X-C motif ligand pathways, whereas CD80/inducible T-cell costimulator (ICOS) signaling was responsible for the immunosuppressive response at day 7 after MI. Finally, the Fos/AP-1 (activator protein 1) regulon was identified as the critical regulator of proinflammatory responses, which was significantly activated in patients with dilated cardiomyopathy and ischemic cardiomyopathy. We showed the enriched Fos/AP-1 target gene loci in genome-wide association study signals for coronary artery diseases and MI. Targeting Fos/AP-1 with the selective inhibitor T5224 blunted leukocyte infiltration and alleviated cardiac dysfunction in the preclinical murine MI model. Conclusions Taken together, this single-cell RNA sequencing data lay the groundwork for the understanding of immune cell heterogeneity and dynamics in murine ischemic hearts. Moreover, Fos/AP-1 inhibition mitigates inflammatory responses and cardiac dysfunction, which might provide potential therapeutic benefits for heart failure intervention after MI.

Single-cell transcriptomics reveals distinct cell response between acute and chronic pulmonary infection of Pseudomonas aeruginosa

Knowledge of the changes in the immune microenvironment during pulmonary bacterial acute and chronic infections is limited. The dissection of immune system may provide a basis for effective therapeutic strategies against bacterial infection. Here, we describe a single immune cell atlas of mouse lungs after acute and chronic Pseudomonas aeruginosa infection using single-cell transcriptomics, multiplex immunohistochemistry, and flow cytometry. Our single-cell transcriptomic analysis revealed large-scale comprehensive changes in immune cell composition and high variation in cell-cell interactions after acute and chronic P. aeruginosa infection. Bacterial infection reprograms the genetic architecture of immune cell populations. We identified specific immune cell types, including Cxcl2⁺ B cells and interstitial macrophages, in response to acute and chronic infection, such as their proportions, distribution, and functional status. Importantly, the patterns of immune cell response are drastically different between acute and chronic infection models. The distinct molecular signatures highlight the importance of the highly dynamic cell-cell interaction process in different pathological conditions, which has not been completely revealed previously. These findings provide a comprehensive and unbiased immune cellular landscape for respiratory pathogenesis in mice, enabling further understanding of immunologic mechanisms in infection and inflammatory diseases.

Single-Cell Transcriptome Identifies Drug-Resistance Signature and Immunosuppressive Microenvironment in Metastatic Small Cell Lung Cancer

Small cell lung cancer (SCLC) is a deadly neuroendocrine malignancy with high metastasis. However, the heterogeneity of metastatic SCLC at the single-cell level remains elusive. The single-cell transcriptome of a total of 24 081 cells in metastatic lymph node samples from seven SCLC patients via endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is examined. Genomic alterations are also examined by whole exome sequencing (WES) and the immune infiltration between SCLC and non-SCLC (NSCLC) is compared using public single-cell RNA sequencing (scRNA-seq) data. It is identified that malignant cells in lymph-node metastatic SCLC have inter-patient and intra-tumor heterogeneity characterized by distinct ASCL1 and NEUROD1 expression patterns. High expression of genes such as FZD8 in WNT pathway is associated with drug resistance in malignant cells. Compared to NSCLC, SCLC harbors a unique immunosuppressive tumor microenvironment. Malignant cells exhibit a pattern of attenuated MHC-I antigen presentation-related gene expression, which is associated with relatively low proportion of exhausted T cells. Natural killer (NK) cells display impaired antitumor function with high expression of TGFBR2. This work characterizes the inter-patient and intra-tumor heterogeneity of metastatic SCLC and uncovers the exhaustion signatures in NK cells, which may pave the way for novel treatments for SCLC including immune checkpoint blockade-based immunotherapy.

Tyrosine phosphatase SHP2 exacerbates psoriasis-like skin inflammation in mice via ERK5-dependent NETosis

Psoriasis is a chronic inflammatory skin disease, often accompanied by increased infiltration of immune cells, especially neutrophils. However, the detailed mechanism of the neutrophil function in psoriasis progression remains unclear. Here, we found that both Src homology-2 domain-containing protein tyrosine phosphatase-2 (SHP2) and neutrophils were highly correlated to developing psoriasis by single-cell ribonucleic acid (RNA) sequencing and experiment verification. The deficiency of SHP2 in neutrophils significantly alleviated psoriasis-like phenotype in an imiquimod-induced murine model. Interestingly, high levels of neutrophil extracellular traps (NETs) were produced in the inflamed lesions of psoriatic patients. In addition, imiquimod-induced psoriasis-like symptoms were remarkably ameliorated in peptidyl arginine deiminase 4 (PAD4) knockout mice, which cannot form NETs. Mechanistically, RNA-seq analysis revealed that SHP2 promoted the formation of NETs in neutrophils via the ERK5 pathway. Functionally, this mechanism resulted in the infiltration of pro-inflammatory cytokines such as TNF-α, IL-1β, IL-6, IL-17A, and CXCL-15, which enhances the inflammatory response in skin lesions and reinforces the cross-talk between neutrophils and keratinocytes, ultimately aggravating psoriasis. Our findings uncover a role for SHP2 in NET release and subsequent cell death known as NETosis in the progression of psoriasis and suggest that SHP2 may be a promising therapeutic target for psoriasis.

Single-cell RNA Seq reveals cellular landscape-specific characteristics and potential etiologies for adolescent idiopathic scoliosis

BACKGROUNDS

Abnormal vertebral growth and development have been found in adolescent idiopathic scoliosis (AIS) patients, and the proliferation and differentiation of bone development-related cells play important roles in its pathogenesis. However, a comprehensive single-cell-level differentiation roadmap in AIS has not been achieved.

METHODS

The present study compared the single-cell level cellular landscapes of spinal cancellous bone tissues between AIS patients and healthy subjects using high throughput single-cell RNA sequencing (scRNA-seq), which covers multiple cellular lineages including osteoblast, chondrocyte, osteoclast and related immunocytes. We constructed the differentiation trajectories of bone development-related cell lineages through pseudotime analysis, and the intercellular-communication networks between bone development-related cells and immunocytes were further developed.

RESULTS

A total of 11 distinct cell clusters were identified according to the genome-wide transcriptome profiles. t-Distributed stochastic neighbor embedding (t-SNE) analysis showed that mesenchymal stem cells (MSC) were classified into three subtypes: MSC-LOXL2, MSC-IGFBP5, and MSC-GJA1. Gene ontology (GO) analysis showed that MSC-GJA1 might possess greater osteoblast differentiation potential than the others. MSC-IGFBP5 was the specific MSC subtype observed only in AIS. There were two distinct gene expression clusters: OB-DPT and OB-OLFML2B, and the counts of osteoblasts derived from AIS was significantly less than that of non-AIS subjects. In AIS patients, MSC-IGFBP5 failed to differentiate into osteoblasts and exhibited negative regulation of cell proliferation and enhanced cell death. CPC-PCNA was found to be the specific chondrocyte progenitor cell (CPC) subtype observed only in AIS patients. The cell counts of OC-BIRC3 in AIS were less than those in controls. Pseudotime analysis suggested two possible distinct osteoclast differentiation patterns in AIS and control subjects. Monocytes in AIS mainly differentiated into OC-CRISP3.

CONCLUSIONS

Our single-cell analysis first revealed differences existed in the cellular states between AIS patients and healthy subjects and found the differentiation disruption of specific MSC and CPC clusters in AIS. Cell communication analysis provided the possible pathogenesis of osteoblast and chondrocyte differentiation dysfunction in AIS.

Intratumoral injection of caerin 1.1 and 1.9 peptides increases the efficacy of vaccinated TC-1 tumor-bearing mice with PD-1 blockade by modulating macrophage heterogeneity and the activation of CD8+ T cells in the tumor microenvironment

Objectives

Developing a vaccine formula that alters the tumor-infiltrating lymphocytes to be more immune active against a tumor is key to the improvement of clinical responses to immunotherapy. Here, we demonstrate that, in conjunction with E7 antigen-specific immunotherapy, and IL-10 and PD-1 blockade, intratumoral administration of caerin 1.1/1.9 peptides improves TC-1 tumor microenvironment (TME) to be more immune active than injection of a control peptide.

Methods

We compared the survival time of vaccinated TC-1 tumor-bearing mice with PD-1 and IL-10 blockade, in combination with a further injection of caerin 1.1/1.9 or control peptides. The tumor-infiltrating haematopoietic cells were examined by flow cytometry. Single-cell transcriptomics and proteomics were used to quantify changes in cellular activity across different cell types within the TME.

Results

The injection of caerin 1.1/1.9 increased the efficacy of vaccinated TC-1 tumor-bearing mice with anti-PD-1 treatment and largely expanded the populations of macrophages and NK cells with higher immune activation level, while reducing immunosuppressive macrophages. More activated CD8⁺ T cells were induced with higher populations of memory and effector-memory CD8⁺ T subsets. Computational integration of the proteome with the single-cell transcriptome supported activation of Stat1-modulated apoptosis and significant reduction in immune-suppressive B-cell function following caerin 1.1 and 1.9 treatment.

Conclusions

Caerin 1.1/1.9-containing treatment results in improved antitumor responses. Harnessing the novel candidate genes preferentially enriched in the immune active cell populations may allow further exploration of distinct macrophages, T cells and their functions in TC-1 tumors.

No evidence that plasmablasts transdifferentiate into developing neutrophils in severe COVID-19 disease

OBJECTIVES

A recent single-cell RNA sequencing study by Wilk et al. suggested that plasmablasts can transdifferentiate into 'developing neutrophils' in patients with severe COVID-19 disease. We explore the evidence for this.

METHODS

We downloaded the original data and code used by the authors in their study to replicate their findings and explore the possibility that regressing out variables may have led the authors to overfit their data.

RESULTS

The lineage relationship between plasmablasts and developing neutrophils breaks down when key features are not regressed out, and the data are not overfitted during the analysis.

CONCLUSION

Plasmablasts do not transdifferentiate into developing neutrophils. The single-cell RNA sequencing is a powerful technique for biological discovery and hypothesis generation. However, caution should be exercised in the bioinformatic analysis and interpretation of the data and findings cross-validated by orthogonal techniques.

SiglecF(HI) Marks Late-Stage Neutrophils of the Infarcted Heart: A Single-Cell Transcriptomic Analysis of Neutrophil Diversification

Background

Neutrophils are thought to be short‐lived first responders to tissue injuries such as myocardial infarction (MI), but little is known about their diversification or dynamics.

Methods and Results

We permanently ligated the left anterior descending coronary arteries of mice and performed single‐cell RNA sequencing and analysis of >28 000 neutrophil transcriptomes isolated from the heart, peripheral blood, and bone marrow of mice on days 1 to 4 after MI or at steady‐state. Unsupervised clustering of cardiac neutrophils revealed 5 major subsets, 3 of which originated in the bone marrow, including a late‐emerging granulocyte expressing SiglecF, a marker classically used to define eosinophils. SiglecF^HI neutrophils represented ≈25% of neutrophils on day 1 and grew to account for >50% of neutrophils by day 4 post‐MI. Validation studies using quantitative polymerase chain reaction of fluorescent‐activated cell sorter sorted Ly6G⁺SiglecF^HI and Ly6G⁺SiglecF^LO neutrophils confirmed the distinct nature of these populations. To confirm that the cells were neutrophils rather than eosinophils, we infarcted GATA‐deficient mice (∆dblGATA) and observed similar quantities of infiltrating Ly6G⁺SiglecF^HI cells despite marked reductions of conventional eosinophils. In contrast to other neutrophil subsets, Ly6G⁺SiglecF^HI neutrophils expressed high levels of Myc‐regulated genes, which are associated with longevity and are consistent with the persistence of this population on day 4 after MI.

Conclusions

Overall, our data provide a spatial and temporal atlas of neutrophil specialization in response to MI and reveal a dynamic proinflammatory cardiac Ly6G⁺SigF⁺(Myc⁺NFϰB⁺) neutrophil that has been overlooked because of negative selection.

Deciphering the regulatory landscape of fetal and adult γδ T-cell development at single-cell resolution

γδ T cells with distinct properties develop in the embryonic and adult thymus and have been identified as critical players in a broad range of infections, antitumor surveillance, autoimmune diseases, and tissue homeostasis. Despite their potential value for immunotherapy, differentiation of γδ T cells in the thymus is incompletely understood. Here, we establish a high-resolution map of γδ T-cell differentiation from the fetal and adult thymus using single-cell RNA sequencing. We reveal novel sub-types of immature and mature γδ T cells and identify an unpolarized thymic population which is expanded in the blood and lymph nodes. Our detailed comparative analysis reveals remarkable similarities between the gene networks active during fetal and adult γδ T-cell differentiation. By performing a combined single-cell analysis of Sox13, Maf, and Rorc knockout mice, we demonstrate sequential activation of these factors during IL-17-producing γδ T-cell (γδT17) differentiation. These findings substantially expand our understanding of γδ T-cell ontogeny in fetal and adult life. Our experimental and computational strategy provides a blueprint for comparing immune cell differentiation across developmental stages.

Multiplexed, Sequential Secretion Analysis of the Same Single Cells Reveals Distinct Effector Response Dynamics Dependent on the Initial Basal State

The effector response of immune cells dictated by an array of secreted proteins is a highly dynamic process, requiring sequential measurement of all relevant proteins from single cells. Herein, a microchip-based, 10-plexed, sequential secretion assay on the same single cells and at the scale of ≈5000 single cells measured simultaneously over 4 time points are shown. It is applied to investigating the time course of single human macrophage response to toll-like receptor 4 (TLR4) ligand lipopolysaccharide (LPS) and reveals four distinct activation modes for different proteins in single cells. Protein secretion dynamics classifies the cells into two major activation states dependent on the basal state of each cell. Single-cell RNA sequencing performed on the same samples at the matched time points further demonstrates the existence of two major activation states at the transcriptional level, which are enriched for translation versus inflammatory programs, respectively. These results show a cell-intrinsic heterogeneous response in a phenotypically homogeneous cell population. This work demonstrates the longitudinal tracking of protein secretion signature in thousands of single cells at multiple time points, providing dynamic information to better understand how individual immune cells react to pathogenic challenges over time and how they together constitute a population response.

Single-cell RNA sequencing reveals developmental heterogeneity among early lymphoid progenitors

Single-cell RNA sequencing is a powerful technology for assessing heterogeneity within defined cell populations. Here, we describe the heterogeneity of a B220⁺CD117^intCD19^-NK1.1^- uncommitted hematopoietic progenitor having combined lymphoid and myeloid potential. Phenotypic and functional assays revealed four subpopulations within the progenitor with distinct lineage developmental potentials. Among them, the Ly6D⁺SiglecH^-CD11c^- fraction was lymphoid-restricted exhibiting strong B-cell potential, whereas the Ly6D^-SiglecH^-CD11c^- fraction showed mixed lympho-myeloid potential. Single-cell RNA sequencing of these subsets revealed that the latter population comprised a mixture of cells with distinct lymphoid and myeloid transcriptional signatures and identified a subgroup as the potential precursor of Ly6D⁺SiglecH^-CD11c^- Subsequent functional assays confirmed that B220⁺CD117^intCD19^-NK1.1^- single cells are, with rare exceptions, not bipotent for lymphoid and myeloid lineages. A B-cell priming gradient was observed within the Ly6D⁺SiglecH^-CD11c^- subset and we propose a herein newly identified subgroup as the direct precursor of the first B-cell committed stage. Therefore, the apparent multipotency of B220⁺CD117^intCD19^-NK1.1^- progenitors results from underlying heterogeneity at the single-cell level and highlights the validity of single-cell transcriptomics for resolving cellular heterogeneity and developmental relationships among hematopoietic progenitors.

Alternative TSSs are co-regulated in single cells in the mouse brain

Alternative transcription start sites (TSSs) have been extensively studied genome‐wide for many cell types and have been shown to be important during development and to regulate transcript abundance between cell types. Likewise, single‐cell gene expression has been extensively studied for many cell types. However, how single cells use TSSs has not yet been examined. In particular, it is unknown whether alternative TSSs are independently expressed, or whether they are co‐activated or even mutually exclusive in single cells. Here, we use a previously published single‐cell RNA‐seq dataset, comprising thousands of cells, to study alternative TSS usage. We find that alternative TSS usage is a regulated process, and the correlation between two TSSs expressed in single cells of the same cell type is surprisingly high. Our findings indicate that TSSs generally are regulated by common factors rather than being independently regulated or stochastically expressed.