Search

Background: Uveal melanoma (UM), the most common adult intraocular tumor, is characterized by high malignancy and poor prognosis in advanced stages. Angiogenesis is critical for UM development, however, not only the role of vascular endothelial dysfunction in UM remains unknown, but also their analysis at the single-cell level has been lacking. A comprehensive analysis is essential to clarify the role of the endothelium in the development of UM. Methods: By using single-cell RNA transcriptomics data of 11 cases of primary and liver metastasis UM, we analyzed the endothelial cell status. In addition, we analyzed and validated ECs in the in vitro model and collected clinical specimens. Subsequently, we explored the impact of endothelial dysfunction on UM cell migration and explored the mechanisms responsible for the endothelial cell abnormalities and the reasons for their peripheral effects. Results: UM metastasis has a significantly higher percentage of vascular endothelial cells compared to in situ tumors, and endothelial cells in metastasis show significant senescence. Senescent endothelial cells in metastatic tumors showed significant Krüppel-like factor 4 (KLF4) upregulation, overexpression of KLF4 in normal endothelial cells induced senescence, and knockdown of KLF4 in senescent endothelium inhibited senescence, suggesting that KLF4 is a driver gene for endothelial senescence. KLF4-induced endothelial senescence drove tumor cell migration through a senescence-associated secretory phenotype (SASP), of which the most important component of the effector was CXCL12 (C-X-C motif chemokine ligand 12), and participated in the composition of the immunosuppressive microenvironment. Conclusion: This study provides an undesirable insight of senescent endothelial cells in promoting UM metastasis. [ABSTRACT FROM AUTHOR]

Hepatitis B surface antigen (HBsAg) loss or functional cure (FC) is considered the optimal therapeutic outcome for patients with chronic hepatitis B (CHB). However, the immune-pathological biomarkers and underlying mechanisms of FC remain unclear. In this study we comprehensively interrogate disease-associated cell states identified within intrahepatic tissue and matched PBMCs (peripheral blood mononuclear cells) from patients with CHB or after FC, at the resolution of single cells, to provide novel insights into putative mechanisms underlying FC. We combined single-cell transcriptomics (single-cell RNA sequencing) with multiparametric flow cytometry-based immune phenotyping, and multiplexed immunofluorescence to elucidate the immunopathological cell states associated with CHB vs. FC. We found that the intrahepatic environment in CHB and FC displays specific cell identities and molecular signatures that are distinct from those found in matched PBMCs. FC is associated with the emergence of an altered adaptive immune response marked by CD4 cytotoxic T lymphocytes, and an activated innate response represented by liver-resident natural killer cells, specific Kupffer cell subtypes and marginated neutrophils. Surprisingly, we found MHC class II-expressing hepatocytes in patients achieving FC, as well as low but persistent levels of covalently closed circular DNA and pregenomic RNA, which may play an important role in FC. Our study provides conceptually novel insights into the immuno-pathological control of HBV cure, and opens exciting new avenues for clinical management, biomarker discovery and therapeutic development. We believe that the discoveries from this study, as it relates to the activation of an innate and altered immune response that may facilitate sustained, low-grade inflammation, may have broader implications in the resolution of chronic viral hepatitis. This study dissects the immuno-pathological cell states associated with functionally cured chronic hepatitis B (defined by the loss of HBV surface antigen or HBsAg). We identified the sustained presence of very low viral load, accessory antigen-presenting hepatocytes, adaptive-memory-like natural killer cells, and the emergence of helper CD4 T cells with cytotoxic or effector-like signatures associated with functional cure, suggesting previously unsuspected alterations in the adaptive immune response, as well as a key role for the innate immune response in achieving or maintaining functional cure. Overall, the insights generated from this study may provide new avenues for the development of alternative therapies as well as patient surveillance for better clinical management of chronic hepatitis B. [Display omitted] • An altered adaptive immune response – associated with CD4 cytotoxic T lymphocytes – emerges in functionally cured patients. • Increased recruitment and margination of neutrophils in the liver suggests a low grade immune active state linked to functional cure. • MHC class II-expressing hepatocytes were identified as accessory antigen-presenting cells. • Alterations in the immuno-pathological cell states identified in the liver environment are not always reflected in peripheral circulation. [ABSTRACT FROM AUTHOR]

Xiaoyan Huang,1,2 Xiangrong Zhao,1,2 Yaping Li,1,2 Yangmeng Feng,1,2 Guoan Zhang,3 Qiyu Wang,4 Cuixiang Xu1,2 1Shaanxi Provincial Key Laboratory of Infection and Immune Diseases, Shaanxi Provincial People’s Hospital, Xi’an, People’s Republic of China; 2Shaanxi Engineering Research Center of Cell Immunology, Shaanxi Provincial People’s Hospital, Xi’an, People’s Republic of China; 3Department of Cardiovascular Surgery, Shaanxi Provincial People’s Hospital, Xi’an, People’s Republic of China; 4Department of Graduate School, Yan’an University, Yan’an, People’s Republic of ChinaCorrespondence: Cuixiang Xu, Email xucuixiang1129@163.comBackground: Dilated cardiomyopathy (DCM) is the second leading cause of heart failure, with intricate pathophysiological underpinnings. In order to shed fresh light on the mechanistic research of DCM, we combined bulk RNA-seq and single-cell RNA-seq (scRNA-seq) data to examine significant cells and genes implicated in the disease.Methods: This analysis employed publicly accessible bulk RNA-seq and scRNA-seq DCM datasets. The scRNA-seq data underwent normalization, principal component, and t-distribution stochastic neighbor embedding analysis. Cell-to-cell communication networks and activity analysis were conducted using CellChat. Utilizing enrichment analysis, the marker genes’ role in the active cells was evaluated. After screening by limma software and weighted gene co-expression network analysis, the differentially expressed genes (DEGs) served as hub genes. Furthermore, these hub genes were subjected to immunological studies, transcription factor expression, and gene set enrichment. Lastly, the expression of the four hub genes and their connection to DCM were verified using the rat models.Results: Fibroblasts and monocytes were chosen as hub cells from among the eight identified cell clusters; their marker genes intersected with DEGs to yield six hub genes. In addition, the six hub genes and the essential module genes intersected to yield four essential genes (ASPN, SFRP4, LUM, and FRZB) that were connected to the Wnt signaling pathway and highly expressed in fibroblast. The four hub DEGs had an expression pattern in the DCM rat model experiment results that was in line with the findings of the bioinformatics study. Additionally, there was a strong correlation between decreased cardiac function and the up-regulation of ASPN, SFRP4, LUM, and FRZB.Conclusion: Ultimately, bulk RNA-seq and scRNA-seq data identified fibroblasts and monocytes as the main cell types implicated in DCM. The highly expressed genes ASPN, FRZB, LUM, and SFRP4 in fibroblasts may aid in the mechanistic investigation of DCM.Keywords: dilated cardiomyopathy, single cell RNA-seq, bulk RNA-seq, weighted gene co-expression network analysis, wnt signaling pathway, fibroblasts

Abstract Single-cell RNA-seq (scRNA-seq) is widely used for transcriptome profiling, but most analyses focus on gene-level events, with less attention devoted to alternative splicing. Here, we present scASfind, a novel computational method to allow for quantitative analysis of cell type-specific splicing events using full-length scRNA-seq data. ScASfind utilizes an efficient data structure to store the percent spliced-in value for each splicing event. This makes it possible to exhaustively search for patterns among all differential splicing events, allowing us to identify marker events, mutually exclusive events, and events involving large blocks of exons that are specific to one or more cell types.

Abstract Background Uveal melanoma (UM), the most common adult intraocular tumor, is characterized by high malignancy and poor prognosis in advanced stages. Angiogenesis is critical for UM development, however, not only the role of vascular endothelial dysfunction in UM remains unknown, but also their analysis at the single-cell level has been lacking. A comprehensive analysis is essential to clarify the role of the endothelium in the development of UM. Methods By using single-cell RNA transcriptomics data of 11 cases of primary and liver metastasis UM, we analyzed the endothelial cell status. In addition, we analyzed and validated ECs in the in vitro model and collected clinical specimens. Subsequently, we explored the impact of endothelial dysfunction on UM cell migration and explored the mechanisms responsible for the endothelial cell abnormalities and the reasons for their peripheral effects. Results UM metastasis has a significantly higher percentage of vascular endothelial cells compared to in situ tumors, and endothelial cells in metastasis show significant senescence. Senescent endothelial cells in metastatic tumors showed significant Krüppel-like factor 4 (KLF4) upregulation, overexpression of KLF4 in normal endothelial cells induced senescence, and knockdown of KLF4 in senescent endothelium inhibited senescence, suggesting that KLF4 is a driver gene for endothelial senescence. KLF4-induced endothelial senescence drove tumor cell migration through a senescence-associated secretory phenotype (SASP), of which the most important component of the effector was CXCL12 (C-X-C motif chemokine ligand 12), and participated in the composition of the immunosuppressive microenvironment. Conclusion This study provides an undesirable insight of senescent endothelial cells in promoting UM metastasis.

Abstract Objective Proliferative nodular formation represents a characteristic pathological feature of benign prostatic hyperplasia (BPH) and serves as the primary cause for prostate volume enlargement and consequent lower urinary tract symptoms (LUTS). Its specific mechanism is largely unknown, although several cellular processes have been reported to be involved in BPH initiation and development and highlighted the crucial role of epithelial cells in proliferative nodular formation. However, the technological limitations hinder the in vivo investigation of BPH patients. Methods The robust cell type decomposition (RCTD) method was employed to integrate spatial transcriptomics and single cell RNA sequencing profiles, enabling the elucidation of epithelial cell alterations during nodular formation. Immunofluorescent and immunohistochemical staining was performed for verification. Results The alterations of epithelial cells during the formation of nodules in BPH was observed, and a distinct subgroup of basal epithelial (BE) cells, referred to as BE5, was identified to play a crucial role in driving this progression through the hypoxia-induced epithelial-mesenchymal transition (EMT) signaling pathway. BE5 served as both the initiating cell during nodular formation and the transitional cell during the transformation from luminal epithelial (LE) to BE cells. A distinguishing characteristic of the BE5 cell subgroup in patients with BPH was its heightened hypoxia and upregulated expression of FOS. Histological verification results confirmed a significant association between c-Fos expression and key biological processes such as hypoxia and cell proliferation, as well as the close relationship between hypoxia and EMT in BPH tissues. Furthermore, a strong link between c-Fos expression and the progression of BPH was also been validated. Additionally, notable functional differences were observed in glandular and stromal nodules regarding BE5 cells, with BE5 in glandular nodules exhibiting enhanced capacities for EMT and cell proliferation characterized by club-like cell markers. Conclusions This study elucidated the comprehensive landscape of epithelial cells during in vivo nodular formation in patients, thereby offering novel insights into the initiation and progression of BPH.

Interactions between host and bacterial cells are integral to human physiology. The complexity of host–microbe interactions extends to different cell types, spatial aspects, and phenotypic heterogeneity, requiring high‐resolution approaches to capture their full complexity. The latest breakthroughs in single‐cell RNA sequencing (scRNA‐seq) have opened up a new era of studies in host–pathogen interactions. Here, we first report a high‐throughput cross‐species dual scRNA‐seq technology by using random primers to simultaneously capture both eukaryotic and bacterial RNAs (scRandom‐seq). Using reference cells, scRandom‐seq can detect individual eukaryotic and bacterial cells with high throughput and high specificity. Acinetobacter baumannii (A.b) is a highly opportunistic and nosocomial pathogen that displays resistance to many antibiotics, posing a significant threat to human health, calling for discoveries and treatment. In the A.b infection model, scRandom‐seq witnessed polarization of THP‐1 derived‐macrophages and the intracellular A.b‐induced ferroptosis‐stress in host cells. The inhibition of ferroptosis by Ferrostatin‐1 (Fer‐1) resulted in the improvement of cell vitality and resistance to A.b infection, indicating the potential to resist related infections. scRandom‐seq provides a high‐throughput cross‐species dual single‐cell RNA profiling tool that will facilitate future discoveries in unraveling the complex interactions of host–microbe interactions in infection systems and tumor micro‐environments. [ABSTRACT FROM AUTHOR]

Objective: Proliferative nodular formation represents a characteristic pathological feature of benign prostatic hyperplasia (BPH) and serves as the primary cause for prostate volume enlargement and consequent lower urinary tract symptoms (LUTS). Its specific mechanism is largely unknown, although several cellular processes have been reported to be involved in BPH initiation and development and highlighted the crucial role of epithelial cells in proliferative nodular formation. However, the technological limitations hinder the in vivo investigation of BPH patients. Methods: The robust cell type decomposition (RCTD) method was employed to integrate spatial transcriptomics and single cell RNA sequencing profiles, enabling the elucidation of epithelial cell alterations during nodular formation. Immunofluorescent and immunohistochemical staining was performed for verification. Results: The alterations of epithelial cells during the formation of nodules in BPH was observed, and a distinct subgroup of basal epithelial (BE) cells, referred to as BE5, was identified to play a crucial role in driving this progression through the hypoxia-induced epithelial-mesenchymal transition (EMT) signaling pathway. BE5 served as both the initiating cell during nodular formation and the transitional cell during the transformation from luminal epithelial (LE) to BE cells. A distinguishing characteristic of the BE5 cell subgroup in patients with BPH was its heightened hypoxia and upregulated expression of FOS. Histological verification results confirmed a significant association between c-Fos expression and key biological processes such as hypoxia and cell proliferation, as well as the close relationship between hypoxia and EMT in BPH tissues. Furthermore, a strong link between c-Fos expression and the progression of BPH was also been validated. Additionally, notable functional differences were observed in glandular and stromal nodules regarding BE5 cells, with BE5 in glandular nodules exhibiting enhanced capacities for EMT and cell proliferation characterized by club-like cell markers. Conclusions: This study elucidated the comprehensive landscape of epithelial cells during in vivo nodular formation in patients, thereby offering novel insights into the initiation and progression of BPH. [ABSTRACT FROM AUTHOR]

Circular RNAs (circRNAs) are a crucial yet relatively unexplored class of transcripts known for their tissue‐ and cell‐type‐specific expression patterns. Despite the advances in single‐cell and spatial transcriptomics, these technologies face difficulties in effectively profiling circRNAs due to inherent limitations in circRNA sequencing efficiency. To address this gap, a deep learning model, CIRI‐deep, is presented for comprehensive prediction of circRNA regulation on diverse types of RNA‐seq data. CIRI‐deep is trained on an extensive dataset of 25 million high‐confidence circRNA regulation events and achieved high performances on both test and leave‐out data, ensuring its accuracy in inferring differential events from RNA‐seq data. It is demonstrated that CIRI‐deep and its adapted version enable various circRNA analyses, including cluster‐ or region‐specific circRNA detection, BSJ ratio map visualization, and trans and cis feature importance evaluation. Collectively, CIRI‐deep's adaptability extends to all major types of RNA‐seq datasets including single‐cell and spatial transcriptomic data, which will undoubtedly broaden the horizons of circRNA research. [ABSTRACT FROM AUTHOR]

Glioblastoma (GBM) cells have the potential to switch from being “proliferative cells” to peritumoral “invasive cells”. Peritumoral GBM cells have highly invasive properties that allow them to survive surgery, leading to recurrence. The mechanisms underlying the manner in which the tumor microenvironment (TME) regulates the invasiveness of GBM remain unclear. Single-cell RNA sequencing analysis revealed heterogeneity in GBM cells, microglia and macrophages. In this study, the Oncostatin M receptor (OSMR) and leukemia inhibitory factor receptor (LIFR) expression indicated higher invasiveness in core GBM cells. Under environmental stress, the expression of OSMR and LIFR were up-regulated with the effect of hypoxic, acidic, and low-glucose conditions in vitro. Functional experiments revealed that TME stress significantly influences the proliferation, migration and invasion of GBM cells. The differences in core/peripheral TMEs in GBM affected the invasive properties, indicating the significant role of OSMR expression within the TME in tumor progression and postoperative therapy.

Abstract Circular RNAs (circRNAs) are a crucial yet relatively unexplored class of transcripts known for their tissue‐ and cell‐type‐specific expression patterns. Despite the advances in single‐cell and spatial transcriptomics, these technologies face difficulties in effectively profiling circRNAs due to inherent limitations in circRNA sequencing efficiency. To address this gap, a deep learning model, CIRI‐deep, is presented for comprehensive prediction of circRNA regulation on diverse types of RNA‐seq data. CIRI‐deep is trained on an extensive dataset of 25 million high‐confidence circRNA regulation events and achieved high performances on both test and leave‐out data, ensuring its accuracy in inferring differential events from RNA‐seq data. It is demonstrated that CIRI‐deep and its adapted version enable various circRNA analyses, including cluster‐ or region‐specific circRNA detection, BSJ ratio map visualization, and trans and cis feature importance evaluation. Collectively, CIRI‐deep's adaptability extends to all major types of RNA‐seq datasets including single‐cell and spatial transcriptomic data, which will undoubtedly broaden the horizons of circRNA research.

Background: Although PD-1 blockade is effective for treating several types of cancer, the efficacy of this agent in glioblastoma is largely limited. To overcome non-responders and the immunosuppressive tumor microenvironment, combinational immunotherapeutic strategies with anti-PD-1 need to be considered. Here, we developed IL-12-secreting mesenchymal stem cells (MSC_IL-12) with glioblastoma tropism and evaluated the therapeutic effects of anti-PD-1, MSC_IL-12, and their combination against glioblastoma. Methods: Therapeutic responses were evaluated using an immunocompetent mouse orthotopic model. Tumor-infiltrating lymphocytes (TILs) were analyzed using immunofluorescent imaging. Single-cell transcriptome was obtained from mouse brains after treatments. Results: Anti-PD-1 and MSC_IL-12 showed complete tumor remission in 25.0% (4/16) and 23.1% (3/13) of glioblastoma-implanted mice, respectively, and their combination yielded synergistic antitumor efficacy indicated by 50.0% (6/12) of complete tumor remission. Analyses of TILs revealed that anti-PD-1 increased CD8+ T cells, while MSC_IL-12 led to infiltration of CD4+ T cells and NK cells. Both therapies reduced frequencies of Tregs. All these aspects observed in each monotherapy group were superimposed in the combination group. Notably, no tumor growth was observed upon rechallenge in cured mice, indicating long-term immunity against glioblastoma provoked by the therapies. Single-cell RNA-seq data confirmed these results and revealed that the combined treatment led to immune-favorable tumor microenvironment–CD4+, CD8+ T cells, effector memory T cells, and activated microglia were increased, whereas exhausted T cells, Tregs, and M2 polarized microglia were reduced. Conclusion: Anti-PD-1 and MSC_IL-12 monotherapies show long-term therapeutic responses, and their combination further enhances antitumor efficacy against glioblastoma via inducing immune-favorable tumor microenvironment.

The emergence of definitive human haematopoietic stem cells (HSCs) from Carnegie Stage (CS) 14 to CS17 in the aorta-gonad-mesonephros (AGM) region is a tightly regulated process. Previously, we conducted spatial transcriptomic analysis of the human AGM region at the end of this period (CS16/CS17) and identified secreted factors involved in HSC development. Here, we extend ouranalysis to investigate the progression of dorso-ventral polarised signalling around the dorsal aorta over the entire period of HSC emergence. Our results reveal a dramatic increase in ventral signalling complexity from the CS13-CS14 transition, coinciding with the first appearance of definitive HSCs. We further observe stage-specific changes in signalling up to CS17, which may underpin the step-wise maturation of HSCs described in the mouse model. The data-rich resource is also presented in an online interface enabling in silico analysis of molecular interactions between spatially defined domains of the AGM region. This resource will be of particular interest for researchers studying mechanisms underlying human HSC development as well as those developing in vitro methods for the generation of clinically relevant HSCs from pluripotent stem cells. [ABSTRACT FROM AUTHOR]

Conventional dendritic cells (cDCs) are antigen-presenting cells (APCs) that play a central role in linking innate and adaptive immunity. cDCs have been well described in a number of different mammalian species, but remain poorly characterised in the chicken. In this study, we use previously described chicken cDC specific reagents, a novel gene-edited chicken line and single-cell RNA sequencing (scRNAseq) to characterise chicken splenic cDCs. In contrast to mammals, scRNAseq analysis indicates that the chicken spleen contains a single, chemokine receptor XCR1 expressing, cDC subset. By sexual maturity the XCR1+ cDC population is the most abundant mononuclear phagocyte cell subset in the chicken spleen. scRNAseq analysis revealed substantial heterogeneity within the chicken splenic XCR1+ cDC population. Immature MHC class II (MHCII)LOW XCR1+ cDCs expressed a range of viral resistance genes. Maturation to MHCIIHIGH XCR1+ cDCs was associated with reduced expression of anti-viral gene expression and increased expression of genes related to antigen presentation via the MHCII and cross-presentation pathways. To visualise and transiently ablate chicken XCR1+ cDCs in situ, we generated XCR1-iCaspase9-RFP chickens using a CRISPR-Cas9 knockin transgenesis approach to precisely edit the XCR1 locus, replacing the XCR1 coding region with genes for a fluorescent protein (TagRFP), and inducible Caspase 9. After inducible ablation, the chicken spleen is initially repopulated by immature CD1.1+ XCR1+ cDCs. XCR1+ cDCs are abundant in the splenic red pulp, in close association with CD8+ T-cells. Knockout of XCR1 prevented this clustering of cDCs withCD8+ T-cells. Taken together these data indicate a conserved role for chicken and mammalian XCR1+ cDCs in driving CD8+ T-cells responses. [ABSTRACT FROM AUTHOR]

The liver is an important organ in the human body which is in charge of many key physiological processes, including metabolism, detoxification, and immunity. Normally, homeostasis of macrophages in the liver maintains liver function. Disruption of macrophage homeostasis, evolution among different cell subsets, and altered function contribute to the progression of non ⁃ alcoholic fatty liver diseases (NAFLD). Single cell RNA⁃seq technology is used to identify specific macrophage subsets for the treatment of non⁃ alcoholic fatty liver diseases. This article will review the classification, heterogeneity, and role of hepatic macrophages in non⁃alcoholic fatty liver diseases [ABSTRACT FROM AUTHOR]

Colorectal cancer (CRC), also known as colon cancer, is the third most common cancer and the fourth most cause of cancer-related death in the world. CRC can be classified into two major subtypes, including microsatellite instability (MSI) and microsatellite stability (MSS), which showed different characteristics in immunotherapy. Low sensitivity of diagnostic biomarkers and metastasis are still the principal cause of mortality, especially in MSI. Here, applying computational programs, we identified recurring expression programs based on single cell RNA sequencing (scRNA-Seq) data of CRC cell lines. Notably, three MSI specific recurring modules were identified by non-negative matrix factorization (NMF). High NMF score genes enriched in the function of metabolism and inflammatory response. Focusing on top specific active transcription factor (TF), RUNX3 (Runt-related transcription factor 3), our results suggest that T cell infiltration was increased in RUNX3 high MSI CRC samples. Unbiased Gene Set Enrichment Analysis (GSEA) showed that RUNX3 was strongly associated with immune and metastasis related functions, such as Interferon Gamma (IFN-γ) and EPITHELIAL MESENCHYMAL TRANSITION (EMT). In addition, RUNX3 shows specific highly activated at epigenetic level in MSI compared with other gastrointestinal carcinomas. Positive correlation between RUNX3 and most immune checkpoints further confirmed RUNX3 might have crucial roles in MSI cancer progression and immunotherapy. Taken together, these results indicate significant tumor heterogeneity of two CRC subtypes at single-cell level and epigenetic modification level. These results also linked transcriptional dysregulation with immune infiltration at single-cell level in MSI, which may advance the application of scRNA-Seq technology in immunotherapy and contribute to developing novel biomarkers of this malignancy.

Background: Malignant tumours, particularly non-small cell lung cancer (NSCLC), pose a significant threat to human health due to their prevalence and lethality. Treatment methods for NSCLC vary greatly among individuals, making it crucial to identify predictive markers. Moreover, during tumour initiation and progression, tumour cells can release signaling molecules to induce polarization of macrophages towards a more tumour friendly M2 phenotype, which can promote tumour growth, metastasis, and drug resistance. Methods: We employed a comprehensive approach, combining bulk RNA-seq and single-cell sequencing analysis. Results: In our study, we used bulk RNA-seq and single-cell sequencing methods to analyze differential cells in NSCLC and adjacent tissues, searching for relevant marker genes that can predict prognosis and drug efficacy. We scrutinized biological phenomena such as macrophage-related gene methylation, copy number variation, and alternative splicing. Additionally, we utilized a co-culture technique of immune and tumour cells to explore the role of these genes in macrophage polarization. Our findings revealed distinct differences in macrophages between cancerous and adjacent tissues. We identified ANP32A, CCL20, ERAP2, MYD88, TMEM126B, TUBB6, and ZNF655 as macrophage-related genes that correlate with NSCLC patient prognosis and immunotherapy efficacy. Notably, ERAP2, TUBB6, CCL20, and TMEM126B can induce macrophage M0 to M2 polarization, promoting tumour proliferation. Conclusion: These findings significantly contribute to our understanding of the NSCLC tumour immune microenvironment. They pave the way for further research into the potential of these genes as targets for regulating tumour occurrence and development.

Abstract Differential gene expression in bulk transcriptomics data can reflect change of transcript abundance within a cell type and/or change in the proportions of cell types. Expression deconvolution methods can help differentiate these scenarios. BEDwARS is a Bayesian deconvolution method designed to address differences between reference signatures of cell types and corresponding true signatures underlying bulk transcriptomic profiles. BEDwARS is more robust to noisy reference signatures and outperforms leading in-class methods for estimating cell type proportions and signatures. Application of BEDwARS to dihydropyridine dehydrogenase deficiency identified the possible involvement of ciliopathy and impaired translational control in the etiology of the disorder.

Abstract As an important part of the stomach, gastric antrum secretes gastrin which can regulate acid secretion and gastric emptying. Although most cell types in the gastric antrum are identified, the comparison of cell composition and gene expression in the gastric antrum among different species are not explored. In this study, we collected antrum epithelial tissues from human, pig, rat and mouse for scRNA-seq and compared cell types and gene expression among species. In pig antral epithelium, we identified a novel cell cluster, which is marked by high expression of AQP5, F3, CLCA1 and RRAD. We also discovered that the porcine antral epithelium has stronger immune function than the other species. Further analysis revealed that this may be due to the insufficient function of porcine immune cells. Together, our results replenish the information of multiple species of gastric antral epithelium at the single cell level and provide resources for understanding the homeostasis maintenance and regeneration of gastric antrum epithelium.

Hemocytes of the colonial tunicate Botrylloides leachii play important roles throughout the animal’s life span, including transport and storage of nutrients, respiration, regeneration, budding, sexual reproduction, defense responses, and tunic generation. Nevertheless, very little is known about the involvement of hemocytes or their functions in the remarkable torpor phenomenon characteristic of this species. Changes in water temperature result in rapid and dramatic morphological changes in which the entire colony degenerates, leaving remnants devoid of feeding and reproductive organs that form compacted and opaque masses of lacunae filled with hemocyte cells. Here, we study hemocyte populations in active and hibernating colonies at the morphological, cellular, and molecular levels using histological and transmission electron microscopy observations, primary cell culture observations, and single-cell transcriptomics. This study defines and captures the different cell types in torpor and control stages and further highlights torpor-associated cell types. Multinucleated cells (MNCs) appear only in torpor stages and under in vitro conditions and are most likely formed from cell fusions. Bacteria-carrying phagocytes are cells specific to a torpor state with yet unknown function. Single-cell sequencing analysis revealed 14 transcriptionally distinct cell clusters. The hibernating colonies had a low throughput of cells, yet all but two transcriptional clusters were present in hibernating colonies. A comparison of gene expressions in the same cell clusters revealed torpor-specific transcriptional modalities in seven of the cell types. Single-cell sequencing generates an enormous amount of valuable data that can serve researchers in future studies of ascidians and torpor phenomena and provide opportunities for future meta-analysis studies. However, the technical challenges create bottlenecks for the full exploitation of single-cell RNA sequencing (scRNA-seq) data. We developed an interactive, searchable, and intuitive cloud-based database where researchers can easily explore the single-cell transcriptomics data generated by us from active and torpid B. leachii colonies. The data can be interrogated and downloaded for further analysis and comparative studies. The Dashboard is available at: http://bleachii.tauberbioinformatics.org/seurat?id=Bleachii-v1. By elucidating hemocyte populations during hibernation, these results provide the basis for future studies of hibernation at the cellular and molecular levels in B. leachii and comparative studies of hibernation phenomena in other organisms.

The liver is critical for the digestive and immune systems. Although the physiology and pathology of liver have been well studied and many scRNA‐seq data are generated, a database and landscape for characterizing cell types and gene expression in different liver diseases or developmental stages at single‐cell resolution are lacking. Hence, scLiverDB is developed, a specialized database for human and mouse liver transcriptomes to unravel the landscape of liver cell types, cell heterogeneity and gene expression at single‐cell resolution across various liver diseases/cell types/developmental stages. To date, 62 datasets including 9,050 samples and 1,741,734 cells is curated. A uniform workflow is used, which included quality control, dimensional reduction, clustering, and cell‐type annotation to analyze datasets on the same platform; integrated manual and automatic methods for accurate cell‐type identification and provided a user‐friendly web interface with multiscale functions. There are two case studies to show the usefulness of scLiverDB, which identified the LTB (lymphotoxin Beta) gene as a potential biomarker of lymphoid cells differentiation and showed the expression changes of Foxa3 (forkhead box A3) in liver chronic progressive diseases. This work provides a crucial resource to resolve molecular and cellular information in normal, diseased, and developing human and mouse livers. [ABSTRACT FROM AUTHOR]

As an important part of the stomach, gastric antrum secretes gastrin which can regulate acid secretion and gastric emptying. Although most cell types in the gastric antrum are identified, the comparison of cell composition and gene expression in the gastric antrum among different species are not explored. In this study, we collected antrum epithelial tissues from human, pig, rat and mouse for scRNA-seq and compared cell types and gene expression among species. In pig antral epithelium, we identified a novel cell cluster, which is marked by high expression of AQP5, F3, CLCA1 and RRAD. We also discovered that the porcine antral epithelium has stronger immune function than the other species. Further analysis revealed that this may be due to the insufficient function of porcine immune cells. Together, our results replenish the information of multiple species of gastric antral epithelium at the single cell level and provide resources for understanding the homeostasis maintenance and regeneration of gastric antrum epithelium. [ABSTRACT FROM AUTHOR]

Conventional dendritic cells (cDCs) are antigen-presenting cells (APCs) that play a central role in linking innate and adaptive immunity. cDCs have been well described in a number of different mammalian species, but remain poorly characterised in the chicken. In this study, we use previously described chicken cDC specific reagents, a novel gene-edited chicken line and single-cell RNA sequencing (scRNAseq) to characterise chicken splenic cDCs. In contrast to mammals, scRNAseq analysis indicates that the chicken spleen contains a single, chemokine receptor XCR1 expressing, cDC subset. By sexual maturity the XCR1+ cDC population is the most abundant mononuclear phagocyte cell subset in the chicken spleen. scRNAseq analysis revealed substantial heterogeneity within the chicken splenic XCR1+ cDC population. Immature MHC class II (MHCII)LOW XCR1+ cDCs expressed a range of viral resistance genes. Maturation to MHCIIHIGH XCR1+ cDCs was associated with reduced expression of anti-viral gene expression and increased expression of genes related to antigen presentation via the MHCII and cross-presentation pathways. To visualise and transiently ablate chicken XCR1+ cDCs in situ, we generated XCR1-iCaspase9-RFP chickens using a CRISPR-Cas9 knockin transgenesis approach to precisely edit the XCR1 locus, replacing the XCR1 coding region with genes for a fluorescent protein (TagRFP), and inducible Caspase 9. After inducible ablation, the chicken spleen is initially repopulated by immature CD1.1+ XCR1+ cDCs. XCR1+ cDCs are abundant in the splenic red pulp, in close association with CD8+ T-cells. Knockout of XCR1 prevented this clustering of cDCs with CD8+ T-cells. Taken together these data indicate a conserved role for chicken and mammalian XCR1+ cDCs in driving CD8+ T-cells responses.

IntroductionSeveral efforts have been made to describe the complexity of T cell heterogeneity during the COVID-19 disease; however, there remain gaps in our understanding in terms of the granularity within.MethodsFor this attempt, we performed a single-cell transcriptomic analysis of 33 individuals (4 healthy, 16 COVID-19 positive patients, and 13 COVID-19 recovered individuals).ResultsWe found CD8+ T cell-biased lymphopenia in COVID-19 patients compared to healthy and recovered individuals. We also found an optimal Th1/Th2 ratio, indicating an effective immune response during COVID-19. Expansion of activated CD4+ T and NK T was detected in the COVID-19-positive individuals. Surprisingly, we found cellular and metal ion homeostasis pathways enriched in the COVID-19-positive individuals compared to the healthy and recovered in the CD8+ T cell populations (CD8+ TCM and CD8+ TEM) as well as activated CD4+ T cells.DiscussionIn summary, the COVID-19-positive individuals exhibit a dynamic T cell mediated response. This response may have a possible association with the dysregulation of non-canonical pathways, including housekeeping functions in addition to the conventional antiviral immune response mediated by the T cell subpopulation. These findings considerably extend our insights into the heterogeneity of T cell response during and post-SARS-CoV-2 infection.

In response to the increasing demand for lung transplantation, ex vivo lung perfusion (EVLP) has extended the number of suitable donor lungs by rehabilitating marginal organs. However despite an expanding use in clinical practice, the responses of the different lung cell types to EVLP are not known. In order to advance our mechanistic understanding and establish a refine tool for improvement of EVLP, we conducted a pioneer study involving single cell RNA-seq on human lungs declined for transplantation. Functional enrichment analyses were performed upon integration of data sets generated at 4 h (clinical duration) and 10 h (prolonged duration) from two human lungs processed to EVLP. Pathways related to inflammation were predicted activated in epithelial and blood endothelial cells, in monocyte-derived macrophages and temporally at 4 h in alveolar macrophages. Pathways related to cytoskeleton signaling/organization were predicted reduced in most cell types mainly at 10 h. We identified a division of labor between cell types for the selected expression of cytokine and chemokine genes that varied according to time. Immune cells including CD4+ and CD8+ T cells, NK cells, mast cells and conventional dendritic cells displayed gene expression patterns indicating blunted activation, already at 4 h in several instances and further more at 10 h. Therefore despite inducing inflammatory responses, EVLP appears to dampen the activation of major lung immune cell types, what may be beneficial to the outcome of transplantation. Our results also support that therapeutics approaches aiming at reducing inflammation upon EVLP should target both the alveolar and vascular compartments. [ABSTRACT FROM AUTHOR]

Abstract Background Infectious Salmon Anaemia Virus (ISAV) is an Orthomixovirus that represents a large problem for salmonid aquaculture worldwide. Current prevention and treatment methods are only partially effective. Genetic selection and genome engineering have the potential to develop ISAV resistant salmon stocks. Both strategies can benefit from an improved understanding of the genomic regulation of ISAV pathogenesis. Here, we used single-cell RNA sequencing of an Atlantic salmon cell line to provide the first high dimensional insight into the transcriptional landscape that underpins host-virus interaction during early ISAV infection. Results Salmon head kidney (SHK-1) cells were single-cell RNA sequenced at 24, 48 and 96 h post-ISAV challenge. At 24 h post infection, cells showed expression signatures consistent with viral entry, with genes such as PI3K, FAK or JNK being upregulated relative to uninfected cells. At 48 and 96 h, infected cells showed a clear anti-viral response, characterised by the expression of IFNA2 or IRF2. Uninfected bystander cells at 48 and 96 h also showed clear transcriptional differences, potentially suggesting paracrine signalling from infected cells. These bystander cells expressed pathways such as mRNA sensing, RNA degradation, ubiquitination or proteasome; and up-regulation of mitochondrial ribosome genes also seemed to play a role in the host response to the infection. Correlation between viral and host genes revealed novel genes potentially key for this fish-virus interaction. Conclusions This study has increased our understanding of the cellular response of Atlantic salmon during ISAV infection and revealed host-virus interactions at the cellular level. Our results highlight various potential key genes in this host-virus interaction, which can be manipulated in future functional studies to increase the resistance of Atlantic salmon to ISAV.

Abstract Ovarian cancer is a highly heterogeneous gynecological malignancy that seriously affects the survival and prognosis of female patients. Single-cell sequencing and transcriptome analysis can effectively characterize tumor heterogeneity to better study the mechanism of occurrence and development. In this study, we identified differentially expressed genes with different differentiation outcomes of tumor cells by analyzing a single-cell dataset. Based on the differentially expressed genes, we explored the differences in function and tumor microenvironment among clusters via consensus clustering. Meanwhile, WGCNA was employed to obtain key genes related to ovarian cancer. On the basis of the TCGA and GEO datasets, we constructed a risk model consisting of 7 genes using the LASSO regression model, and successfully verified that the model was characterized as an independent prognostic factor, efficiently predicting the survival prognosis of patients. In addition, immune signature analysis showed that patients in the high-risk group exhibited lower anti-tumor immune cell infiltration and immunosuppressive status, and had poorer responsiveness to chemotherapeutic drugs and immunotherapy. In conclusion, our study provided a 7-gene prognostic model based on the heterogeneity of OC cells for ovarian cancer patients, which could effectively predict the prognosis of patients and identify the immune microenvironment status of patients.

Abstract Background Glioma is the most prevalent primary tumor of the central nervous system. Glioblastoma multiforme (GBM) is the most malignant form of glioma with an extremely poor prognosis. A novel, regulated cell death form of copper-induced cell death called “cuproptosis” provides a new prospect for cancer treatment by regulating cuproptosis. Methods Data from bulk RNA sequencing (RNA-seq) analysis (The Cancer Genome Atlas cohort and Chinese Glioma Genome Atlas cohort) and single cell RNA-seq (scRNA-seq) analysis were integrated to reveal their relationships. A scoring system was constructed according to the cuproptosis-related gene expression, and core genes were experimentally verified using real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), Western blot (WB), immunohistochemistry (IHC), and immunofluorescence (IF). Moreover, cell counting kit-8 (CCK8), colony formation, 5-ethynyl-2’-deoxyuridine (EdU) incorporation, transwell, and flow cytometry cell cycle were performed to evaluate cell proliferation, invasion, and migration. Results The Cuproptosis Activation Scoring (CuAS) Model has stable and independent prognostic efficacy, as verified by two CGGA datasets. Epiregulin (EREG), the gene of the model has the most contributions in the principal component analysis (PCA), is an onco-immunological gene that can affect immunity by influencing the expression of programmed death-ligand 1 (PD-L1) and mediate the process of cuproptosis by influencing the expression of ferredoxin 1 (FDX1). Single cell transcriptome analysis revealed that high CuAS GBM cells are found in vascular endothelial growth factor A (VEGFA) + malignant cells. Oligodendrocyte transcription factor 1 (OLIG1) + malignant is the original clone, and VEGF and CD99 are the differential pathways of specific cell communication between the high and low CuAS groups. This was also demonstrated by immunofluorescence in the tissue sections. Furthermore, CuAS has therapeutic potential for immunotherapy, and we predict that many drugs (methotrexate, NU7441, KU -0063794, GDC-0941, cabozantinib, and NVP-BEZ235) may be used in patients with high CuAS. Conclusion EREG is the core onco-immunological biomarker of CuAS and modulates the cross-talk between VEGF and CD99 signaling in glioblastoma, and CuAS may provide support for immunotherapy and chemotherapy.

Nrf2, encoded by the gene Nfe2l2, is a broadly expressed transcription factor that regulates gene expression in response to reactive oxygen species (ROS) and oxidative stress. It is commonly referred to as a ubiquitous pathway, but this generalization overlooks work indicating that Nrf2 is essentially unexpressed in some neuronal populations. To explore whether this pattern extends throughout the central nervous system (CNS), we quantified Nfe2l2 expression and chromatin accessibility at the Nfe2l2 locus across multiple single cell datasets. In both the mouse and human CNS, Nfe2l2 was repressed in almost all mature neurons, but highly expressed in non-neuronal support cells, and this pattern was robust across multiple human CNS diseases. A subset of key Nrf2 target genes, like Slc7a11, also remained low in neurons. Thus, these data suggest that while most cells express Nfe2l2, with activity determined by ROS levels, neurons actively avoid Nrf2 activity by keeping Nfe2l2 expression low.

In response to the increasing demand for lung transplantation, ex vivo lung perfusion (EVLP) has extended the number of suitable donor lungs by rehabilitating marginal organs. However despite an expanding use in clinical practice, the responses of the different lung cell types to EVLP are not known. In order to advance our mechanistic understanding and establish a refine tool for improvement of EVLP, we conducted a pioneer study involving single cell RNA-seq on human lungs declined for transplantation. Functional enrichment analyses were performed upon integration of data sets generated at 4 h (clinical duration) and 10 h (prolonged duration) from two human lungs processed to EVLP. Pathways related to inflammation were predicted activated in epithelial and blood endothelial cells, in monocyte-derived macrophages and temporally at 4 h in alveolar macrophages. Pathways related to cytoskeleton signaling/organization were predicted reduced in most cell types mainly at 10 h. We identified a division of labor between cell types for the selected expression of cytokine and chemokine genes that varied according to time. Immune cells including CD4+ and CD8+ T cells, NK cells, mast cells and conventional dendritic cells displayed gene expression patterns indicating blunted activation, already at 4 h in several instances and further more at 10 h. Therefore despite inducing inflammatory responses, EVLP appears to dampen the activation of major lung immune cell types, what may be beneficial to the outcome of transplantation. Our results also support that therapeutics approaches aiming at reducing inflammation upon EVLP should target both the alveolar and vascular compartments.

Introduction: Emodin (EMO), a natural derivative of the anthraquinone family mainly extracted from rhubarb (Rheum palmatum), has previously been demonstrated to possess superior anti-inflammatory properties from a single target or pathway. In order to explore the underlying mechanism of action of EMO against rheumatoid arthritis (RA), a network pharmacology approach was employed. Methods: A gene expression profile from GSE55457 available from the Gene Expression Omnibus (GEO) database was used to identify the targets of EMO action. Further, single cell RNA sequencing data from GEO database of RA patients (GSE159117) were downloaded and analysed. To further investigate the anti-RA effect of EMO on MH7A cells, the expression of IL-6 and IL-1ß were monitored. Finally, RNA-seq analyses were conducted on synovial fibroblasts from EMOtreated. Result: We screened the key targets of EMO against RA using network pharmacology methods, including HMGB1, STAT1, EGR1, NR3C1, EGFR, MAPK14, CASP3, CASP1, IL4, IL13, IKBKB and FN1, and their reliability was verified using ROC curve. Single-cell RNA sequencing data analysis showed that these core target proteins mainly played a role by modulating monocytes. The anti-RA effect of EMO was further verified with MH7A cells, which showed that EMO could block cell differentiation and reduce the expression of IL-6 and IL-1ß. WB experiments confirmed that EMO could affect the expression of COX2, HMBG1 and the phosphorylation of p38. Finally, sequencing of synovial fibroblasts from rats treated with EMO showed consistent results with those predicted and verified, further proving the anti-inflammatory effect of EMO. Conclusion: Our research shows that EMO inhibits inflammatory response of rheumatoid arthritis (RA) by targeting HMGB1, STAT1, EGR1, NR3C1, EGFR, MAPK14, CASP3, CASP1, IL4, IL13, IKBKB, FN1 and Monocytes/macrophages. [ABSTRACT FROM AUTHOR]

Background: Myocardial infarction (MI) is a common cardiac condition with a high incidence of morbidity and mortality. Despite extensive medical treatment for MI, the development and outcomes of post-MI heart failure (HF) continue to be major factors contributing to poor post-MI prognosis. Currently, there are few predictors of post-MI heart failure. Methods: In this study, we re-examined single-cell RNA sequencing and bulk RNA sequencing datasets derived from the peripheral blood samples of patients with myocardial infarction, including patients who developed heart failure and those who did not develop heart failure after myocardial infarction. Using marker genes of the relevant cell subtypes, a signature was generated and validated using relevant bulk datasets and human blood samples. Results: We identified a subtype of immune-activated B cells that distinguished post-MI HF patients from non-HF patients. Polymerase chain reaction was used to confirm these findings in independent cohorts. By combining the specific marker genes of B cell subtypes, we developed a prediction model of 13 markers that can predict the risk of HF in patients after myocardial infarction, providing new ideas and tools for clinical diagnosis and treatment. Conclusion: Sub-cluster B cells may play a significant role in post-MI HF. We found that the STING1, HSPB1, CCL5, ACTN1, and ITGB2 genes in patients with post-MI HF showed the same trend of increase as those without post-MI HF. [ABSTRACT FROM AUTHOR]

Pericytes (PCs) are essential components of the blood brain barrier. Brain PCs are critical for dynamically regulating blood flow, for maintaining vascular integrity and their dysregulation is associated with a myriad of disorders such as Alzheimer's disease. To understand their physiological and molecular functions, studies have increasingly focused on primary brain PC isolation and culture. Multiple methods for PC culture have been developed over the years, however, it is still unclear how primary PCs compare to their in vivo counterparts. To address this question, we compared cultured brain PCs at passage 5 and 20 to adult and embryonic brain PCs directly isolated from mouse brains via single cell RNA-seq. Cultured PCs were highly homogeneous, and were most similar to embryonic PCs, while displaying a significantly different transcriptional profile to adult brain PCs. Cultured PCs downregulated canonical PC markers and extracellular matrix (ECM) genes. Importantly, expression of PC markers and ECM genes could be improved by co-culture with brain endothelial cells, showing the importance of the endothelium in maintaining PC identity and function. Taken together, these results highlight key transcriptional differences between cultured and in vivo PCs which should be considered when performing in vitro experiments with brain PCs. [ABSTRACT FROM AUTHOR]

Background: Tuberculosis (TB) is the deadliest communicable disease in the world with the exception of the ongoing COVID-19 pandemic. Programmed cell death (PCD) patterns play key roles in the development and progression of many disease states such that they may offer value as effective biomarkers or therapeutic targets that can aid in identifying and treating TB patients. Materials and methods: The Gene Expression Omnibus (GEO) was used to gather TB-related datasets after which immune cell profiles in these data were analyzed to examine the potential TB-related loss of immune homeostasis. Profiling of differentially expressed PCD-related genes was performed, after which candidate hub PCD-associated genes were selected via a machine learning approach. TB patients were then stratified into two subsets based on the expression of PCD-related genes via consensus clustering. The potential roles of these PCD-associated genes in other TB-related diseases were further examined. Results: In total, 14 PCD-related differentially expressed genes (DEGs) were identified and highly expressed in TB patient samples and significantly correlated with the abundance of many immune cell types. Machine learning algorithms enabled the selection of seven hub PCD-related genes that were used to establish PCD-associated patient subgroups, followed by the validation of these subgroups in independent datasets. These findings, together with GSVA results, indicated that immune-related pathways were significantly enriched in TB patients exhibiting high levels of PCD-related gene expression, whereas metabolic pathways were significantly enriched in the other patient group. Single cell RNA-seq (scRNA-seq) further highlighted significant differences in the immune status of these different TB patient samples. Furthermore, we used CMap to predict five potential drugs for TB-related diseases. Conclusion: These results highlight clear enrichment of PCD-related gene expression in TB patients and suggest that this PCD activity is closely associated with immune cell abundance. This thus indicates that PCD may play a role in TB progression through the induction or dysregulation of an immune response. These findings provide a foundation for further research aimed at clarifying the molecular drivers of TB, the selection of appropriate diagnostic biomarkers, and the design of novel therapeutic interventions aimed at treating this deadly infectious disease. [ABSTRACT FROM AUTHOR]

Radiotherapy of head and neck cancers frequently leads to irreversible hypofunction of salivary glands, which severely compromises the quality of life and is extremely difficult to treat. We found recently that salivary gland resident macrophages are sensitive to radiation and interact with epithelial progenitors and endothelial cells through homeostatic paracrine factors. Heterogeneous subpopulations of resident macrophages are present in other organs with distinct functions, whereas subpopulations of salivary gland resident macrophages with distinct functions or transcriptional profiles have not been reported yet. Using single-cell RNA sequencing, we found that mouse submandibular glands (SMGs) contain 2 distinct self-renewing resident macrophage subsets, an MHC-IIhi subset present in many other organs and an uncommon Csf2r+ subset. The main source of Csf2 in SMGs are innate lymphoid cells (ILCs) that rely on IL15 for maintenance, while the main source of IL15 protein is Csf2r+ resident macrophages, indicating a homeostatic paracrine interaction between these cells. Csf2r+ resident macrophages are the major source of hepatocyte growth factor (Hgf) that regulates homeostasis of SMG epithelial progenitors. Meanwhile, Csf2r+ resident macrophages are responsive to Hedgehog signaling that can rescue salivary function impaired by radiation. Consistently, irradiation persistently decreased numbers of ILCs and levels of IL15 and Csf2 in SMGs, which were all recovered by transient activation of Hedgehog signaling after radiation. Csf2r+ resident macrophages and MHC-IIhi resident macrophages share transcriptome profiles of perivascular macrophages and macrophages associated with nerves and/or epithelial cells in other organs, respectively, and such niche preferences were supported by lineage tracing and immunofluorescent staining. These findings reveal an uncommon resident macrophage subset that regulates the homeostasis of the salivary gland and is promising as the target to restore salivary gland function impaired by radiation. [ABSTRACT FROM AUTHOR]

Single-cell RNA-seq data contains a lot of dropouts hampering downstream analyses due to the low number and inefficient capture of mRNAs in individual cells. Here, we present Epi-Impute, a computational method for dropout imputation by reconciling expression and epigenomic data. Epi-Impute leverages single-cell ATAC-seq data as an additional source of information about gene activity to reduce the number of dropouts. We demonstrate that Epi-Impute outperforms existing methods, especially for very sparse single-cell RNA-seq data sets, significantly reducing imputation error. At the same time, Epi-Impute accurately captures the primary distribution of gene expression across cells while preserving the gene-gene and cell-cell relationship in the data. Moreover, Epi-Impute allows for the discovery of functionally relevant cell clusters as a result of the increased resolution of scRNA-seq data due to imputation. [ABSTRACT FROM AUTHOR]

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease in which motor neurons throughout the brain and spinal cord progressively degenerate resulting in muscle atrophy, paralysis and death. Recent studies using animal models of ALS implicate multiple cell-types (e.g., astrocytes and microglia) in ALS pathogenesis in the spinal motor systems. To ascertain cellular vulnerability and cell-type specific mechanisms of ALS in the brainstem that orchestrates oral-motor functions, we conducted parallel single cell RNA sequencing (scRNA-seq) analysis using the high-throughput Drop-seq method. We isolated 1894 and 3199 cells from the brainstem of wildtype and mutant SOD1 symptomatic mice respectively, at postnatal day 100. We recovered major known cell types and neuronal subpopulations, such as interneurons and motor neurons, and trigeminal ganglion (TG) peripheral sensory neurons, as well as, previously uncharacterized interneuron subtypes. We found that the majority of the cell types displayed transcriptomic alterations in ALS mice. Differentially expressed genes (DEGs) of individual cell populations revealed cell-type specific alterations in numerous pathways, including previously known ALS pathways such as inflammation (in microglia), stress response (ependymal and an uncharacterized cell population), neurogenesis (astrocytes, oligodendrocytes, neurons), synapse organization and transmission (microglia, oligodendrocyte precursor cells, and neuronal subtypes), and mitochondrial function (uncharacterized cell populations). Other cell-type specific processes altered in SOD1 mutant brainstem include those from motor neurons (axon regeneration, voltage-gated sodium and potassium channels underlying excitability, potassium ion transport), trigeminal sensory neurons (detection of temperature stimulus involved in sensory perception), and cellular response to toxic substances (uncharacterized cell populations). DEGs consistently altered across cell types (e.g., Malat1), as well as cell-type specific DEGs, were identified. Importantly, DEGs from various cell types overlapped with known ALS genes from the literature and with top hits from an existing human ALS genome-wide association study (GWAS), implicating the potential cell types in which the ALS genes function with ALS pathogenesis. Our molecular investigation at single cell resolution provides comprehensive insights into the cell types, genes and pathways altered in the brainstem in a widely used ALS mouse model.

The development of single-cell RNA sequencing (scRNA-seq) has allowed high-resolution analysis of cell-type diversity and transcriptional networks controlling cell-fate specification. To identify the transcriptional networks governing human retinal development, we performed scRNA-seq analysis on 16 time points from developing retina as well as four early stages of retinal organoid differentiation. We identified evolutionarily conserved patterns of gene expression during retinal progenitor maturation and specification of all seven major retinal cell types. Furthermore, we identified gene-expression differences between developing macula and periphery and between distinct populations of horizontal cells. We also identified species-specific patterns of gene expression during human and mouse retinal development. Finally, we identified an unexpected role for ATOH7 expression in regulation of photoreceptor specification during late retinogenesis. These results provide a roadmap to future studies of human retinal development and may help guide the design of cell-based therapies for treating retinal dystrophies.

Low back pain is one of the most serious public health problems worldwide and the major clinical manifestation of intervertebral disc degeneration (IVDD). The key pathological change during IVDD is dysfunction of the annulus fibrosus (AF). However, due to the lack of an in-depth understanding of AF biology, the methods to reconstruct the AF are very limited. In this study, the mice AF cell atlas were decoded by single-cell RNA sequencing to provide a guide for AF reconstruction. The results first identify a new population of AF cells, fibrochondrocyte-like AF cells, which synthesize both collagen I and collagen II and are potential functional cells for AF reconstruction. According to the dual features of the AF extracellular matrix, a composite hydrogel based on the acylation of methacrylated silk fibroin with methacrylated hyaluronic acid was produced. To obtain the ability to stimulate differentiation, the composite hydrogels were combined with a fibrochondrocyte-inducing supplement. Finally, reconstruction of the AF defects, by the novel AF stem cell-loaded composite hydrogel, could be observed, its amount of chondroid matrices recovered to 31.7% of AF aera which is significantly higher than that in other control groups. In summary, this study decodes the AF cell atlas, based on which a novel strategy for AF reconstruction is proposed.

Abstract Background Single-cell RNA sequencing (scRNA-seq) allows the detection of rare cell types in complex tissues. The detection of markers for rare cell types is useful for further biological analysis of, for example, flow cytometry and imaging data sets for either physical isolation or spatial characterization of these cells. However, only a few computational approaches consider the problem of selecting specific marker genes from scRNA-seq data. Results Here, we propose sc2marker, which is based on the maximum margin index and a database of proteins with antibodies, to select markers for flow cytometry or imaging. We evaluated the performances of sc2marker and competing methods in ranking known markers in scRNA-seq data of immune and stromal cells. The results showed that sc2marker performed better than the competing methods in accuracy, while having a competitive running time.

Introduction: The underlying pathophysiological mechanisms of cerebral ischemia reperfusion injury (CIRI) is intricate, and current studies suggest that neuron, astrocyte, microglia, endothelial cell, and pericyte all have different phenotypic changes of specific cell types after ischemic stroke. And microglia account for the largest proportion after CIRI. Previous transcriptomic studies of ischemic stroke have typically focused on the 24 hours after CIRI, obscuring the dynamics of cellular subclusters throughout the disease process. Therefore, traditional methods for identifying cell types and their subclusters may not be sufficient to fully unveil the complexity of single-cell transcriptional profile dynamics caused by an ischemic stroke. Methods: In this study, to explore the dynamic transcriptional profile of single cells after CIRI, we used single-cell State Transition Across-samples of RNA-seq data (scSTAR), a new bioinformatics method, to analyze the single-cell transcriptional profile of day 1, 3, and 7 of transient middle cerebral artery occlusion (tMCAO) mice. Combining our bulk RNA sequences and proteomics data, we found the importance of the integrin beta 2 (Itgb2) gene in postmodeling. And microglia of Itgb2+ and Itgb2- were clustered by the scSTAR method. Finally, the functions of the subpopulations were defined by Matescape, and three different time points after tMCAO were found to exhibit specific functions. Results: Our analysis revealed a dynamic transcriptional profile of single cells in microglia after tMCAO and explored the important role of Itgb2 contributed to microglia by combined transcriptomics and proteomics analysis after modeling. Frontiers in Our further analysis revealed that the Itgb2+ microglia subcluster was mainly involved in energy metabolism, cell cycle, angiogenesis, neuronal myelin formation, and repair at 1, 3, and 7 days after tMCAO, respectively. Discussion: Our results suggested that Itgb2+ microglia act as a time-specific multifunctional immunomodulatory subcluster during CIRI, and the underlying mechanisms remain to be further investigated. [ABSTRACT FROM AUTHOR]

Background: Infectious Salmon Anaemia Virus (ISAV) is an Orthomixovirus that represents a large problem for salmonid aquaculture worldwide. Current prevention and treatment methods are only partially effective. Genetic selection and genome engineering have the potential to develop ISAV resistant salmon stocks. Both strategies can benefit from an improved understanding of the genomic regulation of ISAV pathogenesis. Here, we used single-cell RNA sequencing of an Atlantic salmon cell line to provide the first high dimensional insight into the transcriptional landscape that underpins host-virus interaction during early ISAV infection. Results: Salmon head kidney (SHK-1) cells were single-cell RNA sequenced at 24, 48 and 96 h post-ISAV challenge. At 24 h post infection, cells showed expression signatures consistent with viral entry, with genes such as PI3K, FAK or JNK being upregulated relative to uninfected cells. At 48 and 96 h, infected cells showed a clear anti-viral response, characterised by the expression of IFNA2 or IRF2. Uninfected bystander cells at 48 and 96 h also showed clear transcriptional differences, potentially suggesting paracrine signalling from infected cells. These bystander cells expressed pathways such as mRNA sensing, RNA degradation, ubiquitination or proteasome; and up-regulation of mitochondrial ribosome genes also seemed to play a role in the host response to the infection. Correlation between viral and host genes revealed novel genes potentially key for this fish-virus interaction. Conclusions: This study has increased our understanding of the cellular response of Atlantic salmon during ISAV infection and revealed host-virus interactions at the cellular level. Our results highlight various potential key genes in this host-virus interaction, which can be manipulated in future functional studies to increase the resistance of Atlantic salmon to ISAV. [ABSTRACT FROM AUTHOR]