Search

Abstract Background Aging is a complex, heterogeneous process that has multiple causes. Knowledge on genomic, epigenomic and transcriptomic changes during the aging process shed light on understanding the aging mechanism. A recent breakthrough in biotechnology, single cell RNAseq, is revolutionizing aging study by providing gene expression profile of the entire transcriptome of individual cells. Many interesting information could be inferred from this new type of data with the help of novel computational methods. Results In this manuscript a novel statistical method, penalized Latent Dirichlet Allocation (pLDA), is applied to an aging mouse blood scRNA-seq data set. A pipeline is built for cell type and aging prediction. The sequence of models in the pipeline take scRNA-seq expression counts as input, preprocess the data using pLDA and predict the cell type and aging status. Conclusions pLDA learns a dimension reduced representation of the expression profile. This representation allows identification of cell types and has predictability of the age of cells.

The ability of the CD8+HLA-DR+ regulatory T lymphocytes population to regulate the immune response was first described several years ago. It is known that the suppressive effects of these cells depend on intercellular interactions and are mediated by the expression of checkpoint inhibitor molecules such as CTLA-4, TIM-3, PD-1 and LAG-3. The CD8+HLA-DR+ regulatory T cells also share some properties with conventional CD4+ regulatory T lymphocytes. Nevertheless, the characteristics and function of this subpopulation remain poorly understood. Furthermore, studying the properties of CD8+HLA-DR+ regulatory T cells becomes relevant in the light of general age-associated changes in the human immune system and the increased sensitivity of CD8+T lymphocytes to these changes. Therefore, the aim of this study was to investigate the age dynamics and search for transcriptional signatures of the CD8+HLA-DR+ regulatory T cells. For this purpose, flow cytometric analysis of peripheral blood mononuclear cells from 18 donors aged 21 to 85 years was performed. Bioinformatic analysis of single-cell RNA sequencing data was carried out to search for signatures. It was found that CD8+HLA-DR+ regulatory T cells accumulate with age. The transcriptional signatures of this population consist of genes involved in antigen presentation and cytotoxicity, along with a decrease in the expression of genes encoding proteins of activating protein 1 complex. These data suggest mechanisms of suppressor function of CD8+HLA-DR+ regulatory T lymphocytes associated with the ability of these cells to present antigens and perform cytotoxic activity against effector T lymphocytes. The accumulation of the studied cells may imply a potential influence of CD8+HLA-DR+ regulatory T lymphocytes on the efficiency of adaptive immune response in the aging. Further studies of this population may provide insights into its role in age-related changes in the immune system and develop strategies to improve the immune response in the elderly.

Pheochromocytomas (PC) and paragangliomas (PG) are rare neuroendocrine tumors with varied genetic makeup and are associated with high cardiovascular morbidity and a variable risk of malignancy. The source of the transcriptional heterogeneity of the disease and the underlying biological processes that determine the outcome of PCPG remain largely unclear. We focused on PCPG tumors with germline SDHB and RET mutations, which represent distinct prognostic groups with worse or better prognoses, respectively. We applied single-nuclei RNA sequencing (snRNA-seq) to tissue samples from 11 patients and found high patient-to-patient transcriptome heterogeneity in neuroendocrine tumor cells. The tumor microenvironment also showed heterogeneous profiles, mainly contributed bymacrophages of the immune cell clusters and Schwann cells of the stroma. By performing non-negative matrix factorization, we identified common transcriptional programs active in RET and SDHB, as well as distinct modules, including neuronal development, hormone synthesis and secretion, and DNA replication. Similarities between the transcriptomes of the tumor cells and those of the chromaffin- and precursor cell types suggests different developmental stages at which PC and PG tumors appear to be arrested. [ABSTRACT FROM AUTHOR]

Due to lack of targetable receptors and intertumoral heterogeneity, triple negative breast cancer (TNBC) remains particularly difficult to treat. Doxorubicin (DOX) is typically used as nonselective neoadjuvant chemotherapy, but the diversity of treatment efficacy remains unclear. Comparable to variability in clinical response, an experimental model of TNBC using a 4T1 syngeneic mouse model was found to elicit a differential response to a seven-day treatment regimen of DOX. Single-cell RNA sequencing identified an increase in T cells in tumors that responded to DOX treatment compared to tumors that continued to grow uninhibited. Additionally, compared to resistant tumors, DOX sensitive tumors contained significantly more CD4 T helper cells (339%), γδ T cells (727%), Naïve T cells (278%), and activated CD8 T cells (130%). Furthermore, transcriptional profiles of tumor infiltrated T cells in DOX responsive tumors revealed decreased exhaustion, increased chemokine/cytokine expression, and increased activation and cytotoxic activity. γδ T cell derived IL-17A was identified to be highly abundant in the sensitive tumor microenvironment. IL-17A was also found to directly increase sensitivity of TNBC cells in combination with DOX treatment. In TNBC tumors sensitive to DOX, increased IL-17A levels lead to a direct effect on cancer cell responsiveness and chronic stimulation of tumor infiltrated T cells leading to improved chemotherapeutic efficacy. IL-17A's role as a chemosensitive cytokine in TNBC may offer new opportunities for treating chemoresistant breast tumors and other cancer types.

Neurofibromatosis Type 1 (NF1) is one of the most common genetically inherited disorders that affects 1 in 3000 children annually. Clinical manifestations vary widely but nearly always include the development of cutaneous, plexiform and diffuse neurofibromas that are managed over many years. Recent single-cell transcriptomics profiling efforts of neurofibromas have begun to reveal cell signaling processes. However, the cell signaling networks in mature, non- cutaneous neurofibromas remain unexplored. Here, we present insights into the cellular composition and signaling within mature neurofibromas, contrasting with normal adjacent tissue, in a porcine model of NF1 using single-cell RNA sequencing (scRNA-seq) analysis and histopathological characterization. These neurofibromas exhibited classic diffuse-type histologic morphology and expected patterns of S100, SOX10, GFAP, and CD34 immunohistochemistry. The porcine mature neurofibromas closely resemble human neurofibromas histologically and contain all known cellular components of their human counterparts. The scRNA-seq confirmed the presence of all expected cell types within these neurofibromas and identified novel populations of fibroblasts and immune cells, which may contribute to the tumor microenvironment by suppressing inflammation, promoting M2 macrophage polarization, increasing fibrosis, and driving the proliferation of Schwann cells. Notably, we identified tumor-associated IDO1+/CD274+ (PD-L1)+ dendritic cells, which represent the first such observation in any NF1 animal model and suggest the role of the upregulation of immune checkpoints in mature neurofibromas. Finally, we observed that cell types in the tumor microenvironment are poised to promote immune evasion, extracellular matrix reconstruction, and nerve regeneration. [ABSTRACT FROM AUTHOR]

Introduction: Inflammatory epidermolysis bullosa acquisita (EBA) is characterized by a neutrophilic response to anti-type VII collagen (COL7) antibodies resulting in the development of skin inflammation and blistering. The antibody transfer model of EBA closely mirrors this EBA phenotype. Methods: To better understand the changes induced in neutrophils upon recruitment from peripheral blood into lesional skin in EBA, we performed single-cell RNA-sequencing of whole blood and skin dissociate to capture minimally perturbed neutrophils and characterize their transcriptome. Results: Through this approach, we identified clear distinctions between circulating activated neutrophils and intradermal neutrophils. Most strikingly, the gene expression of multiple C-type lectin receptors, which have previously been reported to orchestrate host defense against fungi and select bacteria, were markedly dysregulated. After confirming the upregulation of Clec4n, Clec4d, and Clec4e in experimental EBA as well as in lesional skin from patients with inflammatory EBA, we performed functional studies in globally deficient Clec4e−/ − and Clec4d−/ − mice as well as in neutrophil-specific Clec4n−/ − mice. Deficiency in these genes did not reduce disease in the EBA model. Discussion: Collectively, our results suggest that while the upregulation of Clec4n, Clec4d, and Clec4e is a hallmark of activated dermal neutrophil populations, their individual contribution to the pathogenesis of EBA is dispensable. [ABSTRACT FROM AUTHOR]

COVID-19 has posed a significant threat to global health. Early data has revealed that IL-6, a key regulatory cytokine, plays an important role in the cytokine storm of COVID-19. Multiple trials are therefore looking at the effects of Tocilizumab, an IL-6 receptor antibody that inhibits IL-6 activity, on treatment of COVID-19, with promising findings. As part of a clinical trial looking at the effects of Tocilizumab treatment on kidney transplant recipients with subclinical rejection, we performed single-cell RNA sequencing of comparing stimulated PBMCs before and after Tocilizumab treatment. We leveraged this data to create an in vitro cytokine storm model, to better understand the effects of Tocilizumab in the presence of inflammation. Tocilizumab-treated cells had reduced expression of inflammatory-mediated genes and biologic pathways, particularly amongst monocytes. These results support the hypothesis that Tocilizumab may hinder the cytokine storm of COVID-19, through a demonstration of biologic impact at the single-cell level.

Background: Human embryonic stem cell-derived cardiomyocytes (hESC-CMs) can be used as a source for cell delivery to remuscularize the heart after myocardial infarction. Despite their therapeutic potential, the emergence of ventricular arrhythmias has limited their application. We previously developed a double reporter hESC line to isolate first heart field (FHF: TBX5 + NKX2-5 +) and second heart field (SHF: TBX5 - NKX2-5 + ) CMs. Herein, we explore the role of TBX5 and its effects on underlying gene regulatory networks driving phenotypical and functional differences between these two populations. Methods: We used a combination of tools and techniques for rapid and unsupervised profiling of FHF and SHF populations at the transcriptional, translational, and functional level including single cell RNA (scRNA) and bulk RNA sequencing, atomic force and quantitative phase microscopy, respirometry, and electrophysiology. Results: Gene ontology analysis revealed three biological processes attributed to TBX5 expression: sarcomeric structure, oxidative phosphorylation, and calcium ion handling. Interestingly, migratory pathways were enriched in SHF population. SHF-like CMs display less sarcomeric organization compared to FHF-like CMs, despite prolonged in vitro culture. Atomic force and quantitative phase microscopy showed increased cellular stiffness and decreased mass distribution over time in FHF compared to SHF populations, respectively. Electrophysiological studies showed longer plateau in action potentials recorded from FHF-like CMs, consistent with their increased expression of calcium handling genes. Interestingly, both populations showed nearly identical respiratory profiles with the only significant functional difference being higher ATP generation-linked oxygen consumption rate in FHF-like CMs. Our findings suggest that FHF-like CMs display more mature features given their enhanced sarcomeric alignment, calcium handling, and decreased migratory characteristics. Finally, pseudotime analyses revealed a closer association of the FHF population to human fetal CMs along the developmental trajectory. Conclusion: Our studies reveal that distinguishing FHF and SHF populations based on TBX5 expression leads to a significant impact on their downstream functional properties. FHF CMs display more mature characteristics such as enhanced sarcomeric organization and improved calcium handling, with closer positioning along the differentiation trajectory to human fetal hearts. These data suggest that the FHF CMs may be a more suitable candidate for cardiac regeneration.

Peripheral neuropathic pain (PNP), neuropathic pain that arises from a damage or disease affecting the peripheral nervous system, is associated with an extremely large disease burden, and there is an increasing and urgent need for new therapies for treating this disorder. In this review we have highlighted therapeutic targets that may be translated into disease modifying therapies for PNP associated with peripheral neuropathy. We have also discussed how genetic studies and novel technologies, such as optogenetics, chemogenetics and single-cell RNAsequencing, have been increasingly successful in revealing novel mechanisms underlying PNP. Additionally, consideration of the role of non-neuronal cells and communication between the skin and sensory afferents is presented to highlight the potential use of drug treatment that could be applied topically, bypassing drug side effects. We conclude by discussing the current difficulties to the development of effective new therapies and, most importantly, how we might improve the translation of targets for peripheral neuropathic pain identified from studies in animal models to the clinic. [ABSTRACT FROM AUTHOR]

COVID-19 has posed a significant threat to global health. Early data has revealed that IL-6, a key regulatory cytokine, plays an important role in the cytokine storm of COVID-19. Multiple trials are therefore looking at the effects of Tocilizumab, an IL-6 receptor antibody that inhibits IL-6 activity, on treatment of COVID-19, with promising findings. As part of a clinical trial looking at the effects of Tocilizumab treatment on kidney transplant recipients with subclinical rejection, we performed single-cell RNA sequencing of comparing stimulated PBMCs before and after Tocilizumab treatment. We leveraged this data to create an in vitro cytokine storm model, to better understand the effects of Tocilizumab in the presence of inflammation. Tocilizumab-treated cells had reduced expression of inflammatory-mediated genes and biologic pathways, particularly amongst monocytes. These results support the hypothesis that Tocilizumab may hinder the cytokine storm of COVID-19, through a demonstration of biologic impact at the single-cell level.

Background: T cell lymphomas (TCL) are a group of heterogeneous diseases with over 40 subtypes. In this study, we identified a novel TCL subtype which was featured by a unique T cell receptor (TCR) presentation, α, β and γ chains were co-existing in a single malignant T cell. Case presentation: A 45-year-old male patient was diagnosed T cell lymphoma after 2-month of abdominal distension and liver enlargement. Combining histology review, PET-CT scanning and immunophenotyes, the patient was not classified to any existing TCL subtypes. To better understand this unclassified TCL case, we performed single cell RNA sequencing paired with TCR sequencing on the patient's PBMC and bone marrow samples. To our surprise, we identified that the malignant T cells had a very rare TCR combination, by expressing two α chains, one β chain and one γ chain simultaneously. We further studied the molecular pathogenesis and tumor cell heterogeneity of this rare TCL subtype. A set of potential therapeutic targets were identified from the transcriptome data, such as CCL5, KLRG1 and CD38. Conclusions: We identified the first TCL case co-expressing α, β and γ chains and dissected its molecular pathogenesis, providing valuable information for precision medicine options for this novel TCL subtype. [ABSTRACT FROM AUTHOR]

Translationally relevant animal models are essential for the successful translation of basic science findings into clinical medicine. While rodent models are widely accessible, there are numerous limitations that prevent the extrapolation of findings to human medicine. One approach to overcome these limitations is to use animal models that are genetically diverse and naturally develop disease. For example, pet dogs spontaneously develop diseases that recapitulate the natural progression seen in humans and live in similar environments alongside humans. Thus, dogs represent a useful animal model for many areas of research. Despite the value of the canine model, species specific reagent limitations have hampered in depth characterization of canine immune cells, which constrains the conclusions that can be drawn from canine immunotherapy studies. To address this need, we used single-cell RNA sequencing to characterize the heterogeneity of circulating leukocytes in healthy dogs (n = 7) and osteosarcoma (OS) affected dogs (n = 10). We present a cellular atlas of leukocytes in healthy dogs, then employ the dataset to investigate the impact of primary OS tumors on the transcriptome of circulating leukocytes. We identified 36 unique cell populations amongst dog circulating leukocytes, with a remarkable amount of heterogeneity in CD4 T cell subtypes. In our comparison of healthy dogs and dogs with OS, we identified relative increases in the abundances of polymorphonuclear (PMN-) and monocytic (M-) myeloid-derived suppressor cells (MDSCs), as well as aberrations in gene expression within myeloid cells. Overall, this study provides a detailed atlas of canine leukocytes and investigates how the presence of osteosarcoma alters the transcriptional profiles of circulating immune cells. [ABSTRACT FROM AUTHOR]

Cortical inhibitory interneurons form a broad spectrum of subtypes. This diversity suggests a division of labor, in which each cell type supports a distinct function. In the present era of optimisation-based algorithms, it is tempting to speculate that these functions were the evolutionary or developmental driving force for the spectrum of interneurons we see in the mature mammalian brain. In this study, we evaluated this hypothesis using the two most common interneuron types, parvalbumin (PV) and somatostatin (SST) expressing cells, as examples. PV and SST interneurons control the activity in the cell bodies and the apical dendrites of excitatory pyramidal cells, respectively, due to a combination of anatomical and synaptic properties. But was this compartment-specific inhibition indeed the function for which PV and SST cells originally evolved? Does the compartmental structure of pyramidal cells shape the diversification of PV and SST interneurons over development? To address these questions, we reviewed and reanalyzed publicly available data on the development and evolution of PV and SST interneurons on one hand, and pyramidal cell morphology on the other. These data speak against the idea that the compartment structure of pyramidal cells drove the diversification into PV and SST interneurons. In particular, pyramidal cells mature late, while interneurons are likely committed to a particular fate (PV vs. SST) during early development. Moreover, comparative anatomy and single cell RNA-sequencing data indicate that PV and SST cells, but not the compartment structure of pyramidal cells, existed in the last common ancestor of mammals and reptiles. Specifically, turtle and songbird SST cells also express the Elfn1 and Cbln4 genes that are thought to play a role in compartment-specific inhibition in mammals. PV and SST cells therefore evolved and developed the properties that allow them to provide compartment-specific inhibition before there was selective pressure for this function. This suggest that interneuron diversity originally resulted from a different evolutionary driving force and was only later co-opted for the compartment-specific inhibition it seems to serve in mammals today. Future experiments could further test this idea using our computational reconstruction of ancestral Elfn1 protein sequences. [ABSTRACT FROM AUTHOR]

Neurofibromatosis Type 1 (NF1) is one of the most common genetically inherited disorders that affects 1 in 3000 children annually. Clinical manifestations vary widely but nearly always include the development of cutaneous, plexiform and diffuse neurofibromas that are managed over many years. Recent single-cell transcriptomics profiling efforts of neurofibromas have begun to reveal cell signaling processes. However, the cell signaling networks in mature, non-cutaneous neurofibromas remain unexplored. Here, we present insights into the cellular composition and signaling within mature neurofibromas, contrasting with normal adjacent tissue, in a porcine model of NF1 using single-cell RNA sequencing (scRNA-seq) analysis and histopathological characterization. These neurofibromas exhibited classic diffuse-type histologic morphology and expected patterns of S100, SOX10, GFAP, and CD34 immunohistochemistry. The porcine mature neurofibromas closely resemble human neurofibromas histologically and contain all known cellular components of their human counterparts. The scRNA-seq confirmed the presence of all expected cell types within these neurofibromas and identified novel populations of fibroblasts and immune cells, which may contribute to the tumor microenvironment by suppressing inflammation, promoting M2 macrophage polarization, increasing fibrosis, and driving the proliferation of Schwann cells. Notably, we identified tumor-associated IDO1+/CD274+ (PD-L1)+ dendritic cells, which represent the first such observation in any NF1 animal model and suggest the role of the upregulation of immune checkpoints in mature neurofibromas. Finally, we observed that cell types in the tumor microenvironment are poised to promote immune evasion, extracellular matrix reconstruction, and nerve regeneration.

IntroductionInflammatory epidermolysis bullosa acquisita (EBA) is characterized by a neutrophilic response to anti-type VII collagen (COL7) antibodies resulting in the development of skin inflammation and blistering. The antibody transfer model of EBA closely mirrors this EBA phenotype.MethodsTo better understand the changes induced in neutrophils upon recruitment from peripheral blood into lesional skin in EBA, we performed single-cell RNA-sequencing of whole blood and skin dissociate to capture minimally perturbed neutrophils and characterize their transcriptome.ResultsThrough this approach, we identified clear distinctions between circulating activated neutrophils and intradermal neutrophils. Most strikingly, the gene expression of multiple C-type lectin receptors, which have previously been reported to orchestrate host defense against fungi and select bacteria, were markedly dysregulated. After confirming the upregulation of Clec4n, Clec4d, and Clec4e in experimental EBA as well as in lesional skin from patients with inflammatory EBA, we performed functional studies in globally deficient Clec4e−/− and Clec4d−/− mice as well as in neutrophil-specific Clec4n−/− mice. Deficiency in these genes did not reduce disease in the EBA model.DiscussionCollectively, our results suggest that while the upregulation of Clec4n, Clec4d, and Clec4e is a hallmark of activated dermal neutrophil populations, their individual contribution to the pathogenesis of EBA is dispensable.

Peripheral neuropathic pain (PNP), neuropathic pain that arises from a damage or disease affecting the peripheral nervous system, is associated with an extremely large disease burden, and there is an increasing and urgent need for new therapies for treating this disorder. In this review we have highlighted therapeutic targets that may be translated into disease modifying therapies for PNP associated with peripheral neuropathy. We have also discussed how genetic studies and novel technologies, such as optogenetics, chemogenetics and single-cell RNA-sequencing, have been increasingly successful in revealing novel mechanisms underlying PNP. Additionally, consideration of the role of non-neuronal cells and communication between the skin and sensory afferents is presented to highlight the potential use of drug treatment that could be applied topically, bypassing drug side effects. We conclude by discussing the current difficulties to the development of effective new therapies and, most importantly, how we might improve the translation of targets for peripheral neuropathic pain identified from studies in animal models to the clinic.

Translationally relevant animal models are essential for the successful translation of basic science findings into clinical medicine. While rodent models are widely accessible, there are numerous limitations that prevent the extrapolation of findings to human medicine. One approach to overcome these limitations is to use animal models that are genetically diverse and naturally develop disease. For example, pet dogs spontaneously develop diseases that recapitulate the natural progression seen in humans and live in similar environments alongside humans. Thus, dogs represent a useful animal model for many areas of research. Despite the value of the canine model, species specific reagent limitations have hampered in depth characterization of canine immune cells, which constrains the conclusions that can be drawn from canine immunotherapy studies. To address this need, we used single-cell RNA sequencing to characterize the heterogeneity of circulating leukocytes in healthy dogs (n = 7) and osteosarcoma (OS) affected dogs (n = 10). We present a cellular atlas of leukocytes in healthy dogs, then employ the dataset to investigate the impact of primary OS tumors on the transcriptome of circulating leukocytes. We identified 36 unique cell populations amongst dog circulating leukocytes, with a remarkable amount of heterogeneity in CD4 T cell subtypes. In our comparison of healthy dogs and dogs with OS, we identified relative increases in the abundances of polymorphonuclear (PMN-) and monocytic (M-) myeloid-derived suppressor cells (MDSCs), as well as aberrations in gene expression within myeloid cells. Overall, this study provides a detailed atlas of canine leukocytes and investigates how the presence of osteosarcoma alters the transcriptional profiles of circulating immune cells.

Cortical inhibitory interneurons form a broad spectrum of subtypes. This diversity suggests a division of labor, in which each cell type supports a distinct function. In the present era of optimisation-based algorithms, it is tempting to speculate that these functions were the evolutionary or developmental driving force for the spectrum of interneurons we see in the mature mammalian brain. In this study, we evaluated this hypothesis using the two most common interneuron types, parvalbumin (PV) and somatostatin (SST) expressing cells, as examples. PV and SST interneurons control the activity in the cell bodies and the apical dendrites of excitatory pyramidal cells, respectively, due to a combination of anatomical and synaptic properties. But was this compartment-specific inhibition indeed the function for which PV and SST cells originally evolved? Does the compartmental structure of pyramidal cells shape the diversification of PV and SST interneurons over development? To address these questions, we reviewed and reanalyzed publicly available data on the development and evolution of PV and SST interneurons on one hand, and pyramidal cell morphology on the other. These data speak against the idea that the compartment structure of pyramidal cells drove the diversification into PV and SST interneurons. In particular, pyramidal cells mature late, while interneurons are likely committed to a particular fate (PV vs. SST) during early development. Moreover, comparative anatomy and single cell RNA-sequencing data indicate that PV and SST cells, but not the compartment structure of pyramidal cells, existed in the last common ancestor of mammals and reptiles. Specifically, turtle and songbird SST cells also express the Elfn1 and Cbln4 genes that are thought to play a role in compartment-specific inhibition in mammals. PV and SST cells therefore evolved and developed the properties that allow them to provide compartment-specific inhibition before there was selective pressure for this function. This suggest that interneuron diversity originally resulted from a different evolutionary driving force and was only later co-opted for the compartment-specific inhibition it seems to serve in mammals today. Future experiments could further test this idea using our computational reconstruction of ancestral Elfn1 protein sequences.

Defining the immunological landscape of human tissue is an important area of research, but challenges include the impact of tissue disaggregation on cell phenotypes and the low abundance of immune cells in many tissues. Here, we describe methods to troubleshoot and standardize Cellular Indexing of Transcriptomes and Epitopes by sequencing (CITE-seq) for studies involving enzymatic digestion of human tissue. We tested epitope susceptibility of 92 antibodies commonly used to differentiate immune lineages and cell states on human peripheral blood mononuclear cells following treatment with an enzymatic digestion cocktail used to isolate islets. We observed CD4, CD8a, CD25, CD27, CD120b, CCR4, CCR6, and PD1 display significant sensitivity to enzymatic treatment, effects that often could not be overcome with alternate antibodies. Comparison of flow cytometry-based CITE-seq antibody titrations and sequencing data supports that for the majority of antibodies, flow cytometry accurately predicts optimal antibody concentrations for CITE-seq. Comparison by CITE-seq of immune cells in enzymatically digested islet tissue and donor-matched spleen not treated with enzymes revealed little digestion-induced epitope cleavage, suggesting increased sensitivity of CITE-seq and/or that the islet structure may protect resident immune cells from enzymes. Within islets, CITE-seq identified immune cells difficult to identify by transcriptional signatures alone, such as distinct tissue-resident T cell subsets, mast cells, and innate lymphoid cells (ILCs). Collectively this study identifies strategies for the rational design and testing of CITE-seq antibodies for single-cell studies of immune cells within islets and other tissues. [ABSTRACT FROM AUTHOR]

Introduction: High stress in the compartment surrounded by the deep fascia can cause acute compartment syndrome (ACS) that may result in necrosis of the limbs. The study aims to investigate the cellular heterogeneity of the deep fascia in ACS patients by single-cell RNA sequencing (scRNA-seq). Methods: We collected deep fascia samples from patients with ACS (high-stress group, HG, n=3) and patients receiving thigh amputation due to osteosarcoma (normal-stress group, NG, n=3). We utilized ultrasound and scanning electron microscopy to observe the morphologic change of the deep fascia, used multiplex staining and multispectral imaging to explore immune cell infiltration, and applied scRNA-seq to investigate the cellular heterogeneity of the deep fascia and to identify differentially expressed genes. Results: Notably, we identified GZMK+ interferon-act CD4 central memory T cells as a specific high-stress compartment subcluster expressing interferon-related genes. Additionally, the changes in the proportions of inflammation-related subclusters, such as the increased proportion of M2 macrophages and decreased proportion of M1 macrophages, may play crucial roles in the balance of pro-inflammatory and anti-inflammatory in the development of ACS. Furthermore, we found that heat shock protein genes were highly expressed but metal ion-related genes (S100 family and metallothionein family) were down-regulated in various subpopulations under high stress. Conclusions: We identified a high stress-specific subcluster and variations in immune cells and fibroblast subclusters, as well as their differentially expressed genes, in ACS patients. Our findings reveal the functions of the deep fascia in the pathophysiology of ACS, providing new approaches for its treatment and prevention. [ABSTRACT FROM AUTHOR]

Defining the immunological landscape of human tissue is an important area of research, but challenges include the impact of tissue disaggregation on cell phenotypes and the low abundance of immune cells in many tissues. Here, we describe methods to troubleshoot and standardize Cellular Indexing of Transcriptomes and Epitopes by sequencing (CITE-seq) for studies involving enzymatic digestion of human tissue. We tested epitope susceptibility of 92 antibodies commonly used to differentiate immune lineages and cell states on human peripheral blood mononuclear cells following treatment with an enzymatic digestion cocktail used to isolate islets. We observed CD4, CD8a, CD25, CD27, CD120b, CCR4, CCR6, and PD1 display significant sensitivity to enzymatic treatment, effects that often could not be overcome with alternate antibodies. Comparison of flow cytometry-based CITE-seq antibody titrations and sequencing data supports that for the majority of antibodies, flow cytometry accurately predicts optimal antibody concentrations for CITE-seq. Comparison by CITE-seq of immune cells in enzymatically digested islet tissue and donor-matched spleen not treated with enzymes revealed little digestion-induced epitope cleavage, suggesting increased sensitivity of CITE-seq and/or that the islet structure may protect resident immune cells from enzymes. Within islets, CITE-seq identified immune cells difficult to identify by transcriptional signatures alone, such as distinct tissue-resident T cell subsets, mast cells, and innate lymphoid cells (ILCs). Collectively this study identifies strategies for the rational design and testing of CITE-seq antibodies for single-cell studies of immune cells within islets and other tissues.

Necroptosis is a programmed cell death playing a significant role in cancer. Although necroptosis has been related to tumor immune environment (TIME) remodeling and cancer prognosis, however, the role of necroptosis-related genes (NRGs) in glioma is still elusive. In this study, a total of 159 NRGs were obtained, and parameters such as mutation rate, copy number variation (CNV), and relative expression level were assessed. Then, we constructed an 18-NRGs-based necroptosis-related signature (NRS) in the TCGA dataset, which could predict the patient's prognosis and was validated in two external CGGA datasets. We also explored the correlation between NRS and glioma TIME, chemotherapy sensitivity, and certain immunotherapy-related factors. The two necroptosis-related subtypes were discovered and could also distinguish the patients' prognosis. Through the glioblastoma (GBM) scRNAseq data analysis, NRGs' expression levels in different GBM patient tissue cell subsets were investigated and the relative necroptosis status of different cell subsets was assessed, with the microglia score culminating among all. Moreover, we found a high infiltration level of immunosuppressive cells in glioma TIME, which was associated with poor prognosis in the high-NRS glioma patient group. Finally, the necroptosis suppressor CASP8 exhibited a high expression in glioma and was associated with poor prognosis. Subsequent experiments were performed in human glioma cell lines and patients' tissue specimens to verify the bioinformatic analytic findings about CASP8. Altogether, this study provides comprehensive evidence revealing a prognostic value of NRGs in glioma, which is associated with TIME regulation. [ABSTRACT FROM AUTHOR]

CD11c+ T-bet+ atypical B cells (ABCs) have been identified in the context of vaccination, acute and chronic infections and autoimmune disease. However, the origins and functions of ABCs remain elusive. A major obstacle in the study of ABCs, and human MBCs more generally, has been the use of different phenotypic markers in different contexts to identify what appear to be phenotypically similar cells. Advances in single-cell RNA sequencing (scRNAseq) technology have allowed researchers to accurately identify ABCs in different immune contexts such as diseases and tissues. Notably, recent studies utilizing single cell techniques have demonstrated ABCs are a highly conserved memory B cell lineage. This analysis has also revealed that ABCs are more abundant in ostensibly healthy donors than previously thought. Nonetheless, the normal function of these cells remains elusive. In this review, we will focus on scRNA-seq studies to discuss recent advances in our understanding about the development and functions of ABCs. [ABSTRACT FROM AUTHOR]

The freshwater snail Biomphalaria glabrata is an intermediate host of Schistosoma mansoni, the agent of human intestinal schistosomiasis. However, much is to be discovered about its innate immune system that appears as a complex black box, in which the immune cells (called hemocytes) play a major role in both cellular and humoral response towards pathogens. Until now, hemocyte classification has been based exclusively on cell morphology and ultrastructural description and depending on the authors considered from 2 to 5 hemocyte populations have been described. In this study, we proposed to evaluate the hemocyte heterogeneity at the transcriptomic level. To accomplish this objective, we used single cell RNA sequencing (scRNAseq) technology coupled to a droplet-based system to separate hemocytes and analyze their transcriptome at a unique cell level in naive Biomphalaria glabrata snails. We were able to demonstrate the presence of 7 hemocyte transcriptomic populations defined by the expression of specific marker genes. As a result, scRNAseq approach showed a high heterogeneity within hemocytes, but provides a detailed description of the different hemocyte transcriptomic populations in B. glabrata supported by distinct cellular functions and lineage trajectory. As a main result, scRNAseq revealed the 3 main population as a super-group of hemocyte diversity but, on the contrary, a great hemocytes plasticity with a probable capacity of hemocytes to engage to different activation pathways. This work opens a new field of research to understand the role of hemocytes particularly in response to pathogens, and towards S. mansoni parasites. [ABSTRACT FROM AUTHOR]

Pheochromocytoma, neuroendocrine tumor, single cell RNA-sequencing, transcriptome, heterogeneity, SDHB, RET, paraganglinoma; Pheochromocytomas (PC) and paragangliomas (PG) are rare neuroendocrine tumors with varied genetic makeup and are associated with high cardiovascular morbidity and a variable risk of malignancy. The source of the transcriptional heterogeneity of the disease and the underlying biological processes that determine the outcome of PCPG remain largely unclear. We focused on PCPG tumors with germline SDHB and RET mutations, which represent distinct prognostic groups with worse or better prognoses, respectively. We applied single-nuclei RNA sequencing (snRNA-seq) to tissue samples from 11 patients and found high patient-to-patient transcriptome heterogeneity in neuroendocrine tumor cells. The tumor microenvironment also showed heterogeneous profiles, mainly contributed by macrophages of the immune cell clusters and Schwann cells of the stroma. By performing non-negative matrix factorization, we identified common transcriptional programs active in RET and SDHB, as well as distinct modules, including neuronal development, hormone synthesis and secretion, and DNA replication. Similarities between the transcriptomes of the tumor cells and those of the chromaffin- and precursor cell types suggests different developmental stages at which PC and PG tumors appear to be arrested. [ABSTRACT FROM AUTHOR]

Due to lack of targetable receptors and intertumoral heterogeneity, triple negative breast cancer (TNBC) remains particularly difficult to treat. Doxorubicin (DOX) is typically used as nonselective neoadjuvant chemotherapy, but the diversity of treatment efficacy remains unclear. Comparable to variability in clinical response, an experimental model of TNBC using a 4T1 syngeneic mouse model was found to elicit a differential response to a seven-day treatment regimen of DOX. Single-cell RNA sequencing identified an increase in T cells in tumors that responded to DOX treatment compared to tumors that continued to grow uninhibited. Additionally, compared to resistant tumors, DOX sensitive tumors contained significantly more CD4 T helper cells (339%), γδ T cells (727%), Naïve T cells (278%), and activated CD8 T cells (130%). Furthermore, transcriptional profiles of tumor infiltrated T cells in DOX responsive tumors revealed decreased exhaustion, increased chemokine/cytokine expression, and increased activation and cytotoxic activity. γδ T cell derived IL-17A was identified to be highly abundant in the sensitive tumor microenvironment. IL-17A was also found to directly increase sensitivity of TNBC cells in combination with DOX treatment. In TNBC tumors sensitive to DOX, increased IL-17A levels lead to a direct effect on cancer cell responsiveness and chronic stimulation of tumor infiltrated T cells leading to improved chemotherapeutic efficacy. IL-17A’s role as a chemosensitive cytokine in TNBC may offer new opportunities for treating chemoresistant breast tumors and other cancer types [ABSTRACT FROM AUTHOR]

The freshwater snail Biomphalaria glabrata is an intermediate host of Schistosoma mansoni, the agent of human intestinal schistosomiasis. However, much is to be discovered about its innate immune system that appears as a complex black box, in which the immune cells (called hemocytes) play a major role in both cellular and humoral response towards pathogens. Until now, hemocyte classification has been based exclusively on cell morphology and ultrastructural description and depending on the authors considered from 2 to 5 hemocyte populations have been described. In this study, we proposed to evaluate the hemocyte heterogeneity at the transcriptomic level. To accomplish this objective, we used single cell RNA sequencing (scRNAseq) technology coupled to a droplet-based system to separate hemocytes and analyze their transcriptome at a unique cell level in naive Biomphalaria glabrata snails. We were able to demonstrate the presence of 7 hemocyte transcriptomic populations defined by the expression of specific marker genes. As a result, scRNAseq approach showed a high heterogeneity within hemocytes, but provides a detailed description of the different hemocyte transcriptomic populations in B. glabrata supported by distinct cellular functions and lineage trajectory. As a main result, scRNAseq revealed the 3 main population as a super-group of hemocyte diversity but, on the contrary, a great hemocytes plasticity with a probable capacity of hemocytes to engage to different activation pathways. This work opens a new field of research to understand the role of hemocytes particularly in response to pathogens, and towards S. mansoni parasites.

CD11c+ T-bet+ atypical B cells (ABCs) have been identified in the context of vaccination, acute and chronic infections and autoimmune disease. However, the origins and functions of ABCs remain elusive. A major obstacle in the study of ABCs, and human MBCs more generally, has been the use of different phenotypic markers in different contexts to identify what appear to be phenotypically similar cells. Advances in single-cell RNA sequencing (scRNA-seq) technology have allowed researchers to accurately identify ABCs in different immune contexts such as diseases and tissues. Notably, recent studies utilizing single cell techniques have demonstrated ABCs are a highly conserved memory B cell lineage. This analysis has also revealed that ABCs are more abundant in ostensibly healthy donors than previously thought. Nonetheless, the normal function of these cells remains elusive. In this review, we will focus on scRNA-seq studies to discuss recent advances in our understanding about the development and functions of ABCs.

Pheochromocytomas (PC) and paragangliomas (PG) are rare neuroendocrine tumors with varied genetic makeup and are associated with high cardiovascular morbidity and a variable risk of malignancy. The source of the transcriptional heterogeneity of the disease and the underlying biological processes that determine the outcome of PCPG remain largely unclear. We focused on PCPG tumors with germline SDHB and RET mutations, which represent distinct prognostic groups with worse or better prognoses, respectively. We applied single-nuclei RNA sequencing (snRNA-seq) to tissue samples from 11 patients and found high patient-to-patient transcriptome heterogeneity in neuroendocrine tumor cells. The tumor microenvironment also showed heterogeneous profiles, mainly contributed by macrophages of the immune cell clusters and Schwann cells of the stroma. By performing non-negative matrix factorization, we identified common transcriptional programs active in RET and SDHB, as well as distinct modules, including neuronal development, hormone synthesis and secretion, and DNA replication. Similarities between the transcriptomes of the tumor cells and those of the chromaffin- and precursor cell types suggests different developmental stages at which PC and PG tumors appear to be arrested.

This document is an erratum published in the journal Frontiers in Oncology. The erratum corrects an error in the abstract of a previous article titled "Extensive patient-to-patient single nucleus transcriptome heterogeneity in pheochromocytomas and paragangliomas." The correction provides a revised abstract that describes the genetic makeup and transcriptional heterogeneity of pheochromocytomas and paragangliomas, rare neuroendocrine tumors. The correction also mentions the identification of common and distinct transcriptional programs in these tumors. The publisher apologizes for the mistake and states that the original article has been updated. [Extracted from the article]

Plasmacytoid dendritic cells [pDCs] represent a rare innate immune subset uniquely endowed with the capacity to produce substantial amounts of type-I interferons. This function of pDCs is critical for effective antiviral defenses and has been implicated in autoimmunity. While IFN-I and select cytokines have been recognized as pDC secreted products, a comprehensive agnostic profiling of the pDC secretome in response to a physiologic stimulus has not been reported. We applied LC-MS/MS to catalogue the repertoire of proteins secreted by pDCs in the unperturbed condition and in response to challenge with influenza H1N1. We report the identification of a baseline pDC secretome, and the repertoire of virus-induced proteins including most type-I interferons, various cytokines, chemokines and granzyme B. Additionally, using single-cell RNA-seq [scRNA-seq], we perform multidimensional analyses of pDC transcriptional diversity immediately ex vivo and following stimulation. Our data evidence preexisting pDC heterogeneity, with subsequent highly specialized roles within the pDC population upon stimulation ranging from dedicated cytokine super-producers to cells with APC-like traits. Dynamic expression of transcription factors and surface markers characterize subclusters within activated pDCs. Integrating the proteomic and transcriptomic datasets confirms the pDC-subcluster origin of the proteins identified in the secretome. Our findings represent the most comprehensive molecular characterization of primary human pDCs at baseline, and in response to influenza virus, reported to date. [ABSTRACT FROM AUTHOR]

Background: Treatment of many types of hearing instability in humans, including sudden sensorineural hearing loss, Meniere's disease, and autoimmune inner ear disease, rely heavily on the utilization of corticosteroids delivered both by oral and transtympanic routes. Despite this use, there is heterogeneity in the response to treatment with corticosteroids in humans with these diseases. The mechanisms by which corticosteroids exert their effect and the cell types in which they exert their effects in the inner ear remain poorly characterized. In this study, we localize steroid-responsive genes to cochlear cell types using previously published transcriptome datasets from the mammalian cochlea. Methods: Steroid-responsive genes were localized to specific cochlear cell types using existing transcriptome datasets from wild-type mammalian cochlea exposed to systemic and transtympanic steroids, as well as previously published single-cell and single-nucleus RNA-sequencing datasets from the mammalian cochlea. Gene ontology (GO) analysis of differentially expressed genes (DEGs) was performed using PANTHER to investigate cellular processes implicated in transtympanic vs. systemic steroid action in the cochlea. Results: Steroid-responsive genes were localized to specific cell types and regions in the cochlea including the stria vascularis, organ of Corti, and spiral ganglion neurons (SGN). Analyses demonstrate differential prevalence of steroid-responsive genes. GO analysis demonstrated steroid-responsive DEGs in the SGN to be associated with angiogenesis, apoptosis, and cytokine-mediated anti-inflammatory pathways. Conclusions: Single-cell and single-nucleus transcriptome datasets localize steroid-responsive genes to specific regions in the cochlea. Further study of these regionally-specific steroid-responsive genes may provide insight into the mechanisms of and clinical response to corticosteroids in diseases of hearing instability. [ABSTRACT FROM AUTHOR]

Plasmacytoid dendritic cells [pDCs] represent a rare innate immune subset uniquely endowed with the capacity to produce substantial amounts of type-I interferons. This function of pDCs is critical for effective antiviral defenses and has been implicated in autoimmunity. While IFN-I and select cytokines have been recognized as pDC secreted products, a comprehensive agnostic profiling of the pDC secretome in response to a physiologic stimulus has not been reported. We applied LC-MS/MS to catalogue the repertoire of proteins secreted by pDCs in the unperturbed condition and in response to challenge with influenza H1N1. We report the identification of a baseline pDC secretome, and the repertoire of virus-induced proteins including most type-I interferons, various cytokines, chemokines and granzyme B. Additionally, using single-cell RNA-seq [scRNA-seq], we perform multidimensional analyses of pDC transcriptional diversity immediately ex vivo and following stimulation. Our data evidence preexisting pDC heterogeneity, with subsequent highly specialized roles within the pDC population upon stimulation ranging from dedicated cytokine super-producers to cells with APC-like traits. Dynamic expression of transcription factors and surface markers characterize subclusters within activated pDCs. Integrating the proteomic and transcriptomic datasets confirms the pDC-subcluster origin of the proteins identified in the secretome. Our findings represent the most comprehensive molecular characterization of primary human pDCs at baseline, and in response to influenza virus, reported to date.

BackgroundTreatment of many types of hearing instability in humans, including sudden sensorineural hearing loss, Meniere's disease, and autoimmune inner ear disease, rely heavily on the utilization of corticosteroids delivered both by oral and transtympanic routes. Despite this use, there is heterogeneity in the response to treatment with corticosteroids in humans with these diseases. The mechanisms by which corticosteroids exert their effect and the cell types in which they exert their effects in the inner ear remain poorly characterized. In this study, we localize steroid-responsive genes to cochlear cell types using previously published transcriptome datasets from the mammalian cochlea.MethodsSteroid-responsive genes were localized to specific cochlear cell types using existing transcriptome datasets from wild-type mammalian cochlea exposed to systemic and transtympanic steroids, as well as previously published single-cell and single-nucleus RNA-sequencing datasets from the mammalian cochlea. Gene ontology (GO) analysis of differentially expressed genes (DEGs) was performed using PANTHER to investigate cellular processes implicated in transtympanic vs. systemic steroid action in the cochlea.ResultsSteroid-responsive genes were localized to specific cell types and regions in the cochlea including the stria vascularis, organ of Corti, and spiral ganglion neurons (SGN). Analyses demonstrate differential prevalence of steroid-responsive genes. GO analysis demonstrated steroid-responsive DEGs in the SGN to be associated with angiogenesis, apoptosis, and cytokine-mediated anti-inflammatory pathways.ConclusionsSingle-cell and single-nucleus transcriptome datasets localize steroid-responsive genes to specific regions in the cochlea. Further study of these regionally-specific steroid-responsive genes may provide insight into the mechanisms of and clinical response to corticosteroids in diseases of hearing instability.

Background: Aging is a complex, heterogeneous process that has multiple causes. Knowledge on genomic, epigenomic and transcriptomic changes during the aging process shed light on understanding the aging mechanism. A recent breakthrough in biotechnology, single cell RNAseq, is revolutionizing aging study by providing gene expression profile of the entire transcriptome of individual cells. Many interesting information could be inferred from this new type of data with the help of novel computational methods., Results: In this manuscript a novel statistical method, penalized Latent Dirichlet Allocation (pLDA), is applied to an aging mouse blood scRNA-seq data set. A pipeline is built for cell type and aging prediction. The sequence of models in the pipeline take scRNA-seq expression counts as input, preprocess the data using pLDA and predict the cell type and aging status., Conclusions: pLDA learns a dimension reduced representation of the expression profile. This representation allows identification of cell types and has predictability of the age of cells., (© 2024. The Author(s).)

Solid organ transplant recipients require long-term immunosuppression for prevention of rejection. Calcineurin inhibitor (CNI)-based immunosuppressive regimens have remained the primary means for immunosuppression for four decades now, yet little is known about their effects on graft resident and infiltrating immune cell populations. Similarly, the understanding of rejection biology under specific types of immunosuppression remains to be defined. Furthermore, development of innovative, rationally designed targeted therapeutics for mitigating or preventing rejection requires a fundamental understanding of the immunobiology that underlies the rejection process. The established use of microarray technologies in transplantation has provided great insight into gene transcripts associated with allograft rejection but does not characterize rejection on a single cell level. Therefore, the development of novel genomics tools, such as single cell sequencing techniques, combined with powerful bioinformatics approaches, has enabled characterization of immune processes at the single cell level. This can provide profound insights into the rejection process, including identification of resident and infiltrating cell transcriptomes, cell-cell interactions, and T cell receptor α/β repertoires. In this review, we discuss genomic analysis techniques, including microarray, bulk RNAseq (bulkSeq), single-cell RNAseq (scRNAseq), and spatial transcriptomic (ST) techniques, including considerations of their benefits and limitations. Further, other techniques, such as chromatin analysis via assay for transposase-accessible chromatin sequencing (ATACseq), bioinformatic regulatory network analyses, and protein-based approaches are also examined. Application of these tools will play a crucial role in redefining transplant rejection with single cell resolution and likely aid in the development of future immunomodulatory therapies in solid organ transplantation. [ABSTRACT FROM AUTHOR]

Background: Human embryonic stem cell-derived cardiomyocytes (hESC-CMs) can be used as a source for cell delivery to remuscularize the heart after myocardial infarction. Despite their therapeutic potential, the emergence of ventricular arrhythmias has limited their application. We previously developed a double reporter hESC line to isolate first heart field (FHF: TBX5+NKX2-5+) and second heart field (SHF: TBX5-NKX2-5+) CMs. Herein, we explore the role of TBX5 and its effects on underlying gene regulatory networks driving phenotypical and functional differences between these two populations.Methods: We used a combination of tools and techniques for rapid and unsupervised profiling of FHF and SHF populations at the transcriptional, translational, and functional level including single cell RNA (scRNA) and bulk RNA sequencing, atomic force and quantitative phase microscopy, respirometry, and electrophysiology.Results: Gene ontology analysis revealed three biological processes attributed to TBX5 expression: sarcomeric structure, oxidative phosphorylation, and calcium ion handling. Interestingly, migratory pathways were enriched in SHF population. SHF-like CMs display less sarcomeric organization compared to FHF-like CMs, despite prolonged in vitro culture. Atomic force and quantitative phase microscopy showed increased cellular stiffness and decreased mass distribution over time in FHF compared to SHF populations, respectively. Electrophysiological studies showed longer plateau in action potentials recorded from FHF-like CMs, consistent with their increased expression of calcium handling genes. Interestingly, both populations showed nearly identical respiratory profiles with the only significant functional difference being higher ATP generation-linked oxygen consumption rate in FHF-like CMs. Our findings suggest that FHF-like CMs display more mature features given their enhanced sarcomeric alignment, calcium handling, and decreased migratory characteristics. Finally, pseudotime analyses revealed a closer association of the FHF population to human fetal CMs along the developmental trajectory.Conclusion: Our studies reveal that distinguishing FHF and SHF populations based on TBX5 expression leads to a significant impact on their downstream functional properties. FHF CMs display more mature characteristics such as enhanced sarcomeric organization and improved calcium handling, with closer positioning along the differentiation trajectory to human fetal hearts. These data suggest that the FHF CMs may be a more suitable candidate for cardiac regeneration.

Large collections of annotated single-cell RNA sequencing (scRNA-seq) experiments are being generated across different organs, conditions and organisms on different platforms. The diversity, volume and complexity of this aggregated data requires new analysis techniques to extract actionable knowledge. Fundamental to most analysis are key abilities such as: identification of similar cells across different experiments and transferring annotations from an annotated dataset to an unannotated one. There have been many strategies explored in achieving these goals, and they focuses primarily on aligning and re-clustering datasets of interest. In this work, we are interested in exploring the applicability of deep metric learning methods as a form of distance function to capture similarity between cells and facilitate the transfer of cell type annotation for similar cells across different experiments. Toward this aim, we developed MapCell, a few-shot training approach using Siamese Neural Networks (SNNs) to learn a generalizable distance metric that can differentiate between single cell types. Requiring only a small training set, we demonstrated that SNN derived distance metric can perform accurate transfer of annotation across different scRNA-seq platforms, batches, species and also aid in flagging novel cell types.

Solid organ transplant recipients require long-term immunosuppression for prevention of rejection. Calcineurin inhibitor (CNI)-based immunosuppressive regimens have remained the primary means for immunosuppression for four decades now, yet little is known about their effects on graft resident and infiltrating immune cell populations. Similarly, the understanding of rejection biology under specific types of immunosuppression remains to be defined. Furthermore, development of innovative, rationally designed targeted therapeutics for mitigating or preventing rejection requires a fundamental understanding of the immunobiology that underlies the rejection process. The established use of microarray technologies in transplantation has provided great insight into gene transcripts associated with allograft rejection but does not characterize rejection on a single cell level. Therefore, the development of novel genomics tools, such as single cell sequencing techniques, combined with powerful bioinformatics approaches, has enabled characterization of immune processes at the single cell level. This can provide profound insights into the rejection process, including identification of resident and infiltrating cell transcriptomes, cell-cell interactions, and T cell receptor α/β repertoires. In this review, we discuss genomic analysis techniques, including microarray, bulk RNAseq (bulkSeq), single-cell RNAseq (scRNAseq), and spatial transcriptomic (ST) techniques, including considerations of their benefits and limitations. Further, other techniques, such as chromatin analysis via assay for transposase-accessible chromatin sequencing (ATACseq), bioinformatic regulatory network analyses, and protein-based approaches are also examined. Application of these tools will play a crucial role in redefining transplant rejection with single cell resolution and likely aid in the development of future immunomodulatory therapies in solid organ transplantation.

The catalog of the Drosophila immune cells was until recently limited to three major cell types, based on morphology, function and few molecular markers. Three recent single cell studies highlight the presence of several subgroups, revealing a large diversity in the molecular signature of the larval immune cells. Since these studies rely on somewhat different experimental and analytical approaches, we here compare the datasets and identify eight common, robust subgroups associated to distinct functions such as proliferation, immune response, phagocytosis or secretion. Similar comparative analyses with datasets from different stages and tissues disclose the presence of larval immune cells resembling embryonic hemocyte progenitors and the expression of specific properties in larval immune cells associated with peripheral tissues.

COVID-19 has posed a significant threat to global health. Early data has revealed that IL-6, a key regulatory cytokine, plays an important role in the cytokine storm of COVID-19. Multiple trials are therefore looking at the effects of Tocilizumab, an IL-6 receptor antibody that inhibits IL-6 activity, on treatment of COVID-19, with promising findings. As part of a clinical trial looking at the effects of Tocilizumab treatment on kidney transplant recipients with subclinical rejection, we performed single-cell RNA sequencing of comparing stimulated PBMCs before and after Tocilizumab treatment. We leveraged this data to create an in vitro cytokine storm model, to better understand the effects of Tocilizumab in the presence of inflammation. Tocilizumab-treated cells had reduced expression of inflammatory-mediated genes and biologic pathways, particularly amongst monocytes. These results support the hypothesis that Tocilizumab may hinder the cytokine storm of COVID-19, through a demonstration of biologic impact at the single-cell level.

Background: Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that mainly transfers from human to human via respiratory and gastrointestinal routes. The S-glycoprotein in the virus is the key factor for the entry of SARS-CoV-2 into the cell, which contains two functional domains: S1 is an angiotensin-converting enzyme 2 (ACE2) receptor binding domain, and S2 is necessary for fusion of the coronavirus and cell membranes. Moreover, it has been reported that ACE2 is likely to be the receptor for SARS-CoV-2. In addition, mRNA level expression of Furin enzyme and ACE2 receptor had been reported in airway epithelia, cardiac tissue, and enteric canals. However, the expression patterns of ACE2 and Furin in different cell types of oral tissues are still unclear.Methods: In order to investigate the potential infective channel of the new coronavirus via the oropharyngeal cavity, we analyze the expression of ACE2 and Furin in human oral mucosa using the public single-cell sequence datasets. Furthermore, immunohistochemistry was performed in mucosal tissue from different oral anatomical sites to confirm the expression of ACE2 and Furin at the protein level.Results: The bioinformatics results indicated the differential expression of ACE2 and Furin on epithelial cells from different oral anatomical sites. Immunohistochemistry results revealed that both the ACE2-positive and Furin-positive cells in the target tissues were mainly positioned in the epithelial layers, partly expressed in fibroblasts, further confirming the bioinformatics results.Conclusions: Based on these findings, we speculated that SARS-CoV-2 could invade oral mucosal cells through two possible routes: binding to the ACE2 receptor and fusion with cell membrane activated by Furin protease. Our results indicated that oral mucosa tissues are susceptible to SARS-CoV-2 that could facilitate COVID-19 infection via respiratory and fecal–oral routes.

The gonads of vertebrate embryos are unique among organs because they have a developmental choice; ovary or testis formation. Given the importance of proper gonad formation for sexual development and reproduction, considerable research has been conducted over the years to elucidate the genetic and cellular mechanisms of gonad formation and sexual differentiation. While the molecular trigger for gonadal sex differentiation into ovary of testis can vary among vertebrates, from egg temperature to sex-chromosome linked master genes, the downstream molecular pathways are largely conserved. The cell biology of gonadal formation and differentiation has long thought to also be conserved. However, recent discoveries point to divergent mechanisms of gonad formation, at least among birds and mammals. In this mini-review, we provide an overview of cell lineage allocation during gonadal sex differentiation in the mouse model, focusing on the key supporting and steroidogenic cells and drawing on recent insights provided by single cell RNA-sequencing. We compare this data with emerging information in the chicken model. We highlight surprising differences in cell lineage specification between species and identify gaps in our current understanding of the cell biology underlying gonadogenesis.

Objectives: To identify sub-populations of intestinal preproglucagon-expressing (PPG) cells producing Glucagon-like Peptide-1, and their associated expression profiles of sensory receptors, thereby enabling the discovery of therapeutic strategies that target these cell populations for the treatment of diabetes and obesity. Methods: We performed single cell RNA sequencing of PPG-cells purified by flow cytometry from the upper small intestine of 3 GLU-Venus mice. Cells from 2 mice were sequenced at low depth, and from the third mouse at high depth. High quality sequencing data from 234 PPG-cells were used to identify clusters by tSNE analysis. qPCR was performed to compare the longitudinal and crypt/villus locations of cluster-specific genes. Immunofluorescence and mass spectrometry were used to confirm protein expression. Results: PPG-cells formed 3 major clusters: a group with typical characteristics of classical L-cells, including high expression of Gcg and Pyy (comprising 51% of all PPG-cells); a cell type overlapping with Gip-expressing K-cells (14%); and a unique cluster expressing Tph1 and Pzp that was predominantly located in proximal small intestine villi and co-produced 5-HT (35%). Expression of G-protein coupled receptors differed between clusters, suggesting the cell types are differentially regulated and would be differentially targetable. Conclusions: Our findings support the emerging concept that many enteroendocrine cell populations are highly overlapping, with individual cells producing a range of peptides previously assigned to distinct cell types. Different receptor expression profiles across the clusters highlight potential drug targets to increase gut hormone secretion for the treatment of diabetes and obesity.