Integrative analysis of COL6A3 in lupus nephritis: insights from single-cell transcriptomics and proteomics.

Introduction: Lupus nephritis (LN), a severe complication of systemic lupus erythematosus (SLE), presents significant challenges in patient management and treatment outcomes. The identification of novel LN-related biomarkers and therapeutic targets is critical to enhancing treatment outcomes and prognosis for patients. Methods: In this study, we analyzed single-cell expression data from LN (n=21) and healthy controls (n=3). A total of 143 differentially expressed genes were identified between the LN and control groups. Then, proteomics analysis of LN patients (n=9) and control (SLE patients without LN, n=11) revealed 55 differentially expressed genes among patients with LN and control group. We further utilizes protein-protein interaction network and functional enrichment analyses to elucidate the pivotal role of COL6A3 in key signaling pathways. Its diagnostic value is evaluate through its correlation with disease progression and renal function metrics, as well as Receiver Operating Characteristic Curve (ROC) analysis. Additionally, immunohistochemistry and qPCR experiments were performed to validate the expression of COL6A3 in LN. Results: By comparison of single-cell and proteomics data, we discovered that COL6A3 is significantly upregulated, highlighting it as a critical biomarker of LN. Our findings emphasize the substantial involvement of COL6A3 in the pathogenesis of LN, particularly noting its expression in mesangial cells. Through comprehensive protein-protein interaction network and functional enrichment analyses, we uncovered the pivotal role of COL6A3 in key signaling pathways including integrin-mediated signaling pathways, collagen-activated signaling pathways, and ECM-receptor interaction, suggesting potential therapeutic targets. The diagnostic utility is confirmed by its correlation with disease progression and renal function metrics of the glomerular filtration rate. ROC analysis further validates the diagnostic value of COL6A3, with the area under the ROC values of 0.879 in the in-house cohort, and 0.802 and 0.915 in tubular and glomerular external cohort samples, respectively. Furthermore, immunohistochemistry and qPCR experiments were consistent with those obtained from the single-cell RNA sequencing and proteomics studies. Discussion: These results proved that COL6A3 is a promising biomarker and therapeutic target, advancing personalized medicine strategies for LN.

Unveiling the molecular complexity of proliferative diabetic retinopathy through scRNA-seq, AlphaFold 2, and machine learning.

Background: Proliferative diabetic retinopathy (PDR), a major cause of blindness, is characterized by complex pathogenesis. This study integrates single-cell RNA sequencing (scRNA-seq), Non-negative Matrix Factorization (NMF), machine learning, and AlphaFold 2 methods to explore the molecular level of PDR. Methods: We analyzed scRNA-seq data from PDR patients and healthy controls to identify distinct cellular subtypes and gene expression patterns. NMF was used to define specific transcriptional programs in PDR. The oxidative stress-related genes (ORGs) identified within Meta-Program 1 were utilized to construct a predictive model using twelve machine learning algorithms. Furthermore, we employed AlphaFold 2 for the prediction of protein structures, complementing this with molecular docking to validate the structural foundation of potential therapeutic targets. We also analyzed protein-protein interaction (PPI) networks and the interplay among key ORGs. Results: Our scRNA-seq analysis revealed five major cell types and 14 subcell types in PDR patients, with significant differences in gene expression compared to those in controls. We identified three key meta-programs underscoring the role of microglia in the pathogenesis of PDR. Three critical ORGs (ALKBH1, PSIP1, and ATP13A2) were identified, with the best-performing predictive model demonstrating high accuracy (AUC of 0.989 in the training cohort and 0.833 in the validation cohort). Moreover, AlphaFold 2 predictions combined with molecular docking revealed that resveratrol has a strong affinity for ALKBH1, indicating its potential as a targeted therapeutic agent. PPI network analysis, revealed a complex network of interactions among the hub ORGs and other genes, suggesting a collective role in PDR pathogenesis. Conclusion: This study provides insights into the cellular and molecular aspects of PDR, identifying potential biomarkers and therapeutic targets using advanced technological approaches.

Advancing diabetes treatment: the role of mesenchymal stem cells in islet transplantation.

Diabetes mellitus, a prevalent global health challenge, significantly impacts societal and economic well-being. Islet transplantation is increasingly recognized as a viable treatment for type 1 diabetes that aims to restore endogenous insulin production and mitigate complications associated with exogenous insulin dependence. We review the role of mesenchymal stem cells (MSCs) in enhancing the efficacy of islet transplantation. MSCs, characterized by their immunomodulatory properties and differentiation potential, are increasingly seen as valuable in enhancing islet graft survival, reducing immune-mediated rejection, and supporting angiogenesis and tissue repair. The utilization of MSC-derived extracellular vesicles further exemplifies innovative approaches to improve transplantation outcomes. However, challenges such as MSC heterogeneity and the optimization of therapeutic applications persist. Advanced methodologies, including artificial intelligence (AI) and single-cell RNA sequencing (scRNA-seq), are highlighted as potential technologies for addressing these challenges, potentially steering MSC therapy toward more effective, personalized treatment modalities for diabetes. This review revealed that MSCs are important for advancing diabetes treatment strategies, particularly through islet transplantation. This highlights the importance of MSCs in the field of regenerative medicine, acknowledging both their potential and the challenges that must be navigated to fully realize their therapeutic promise.

Unraveling Key Factors for Hypoxia Tolerance in Contrasting Varieties of Cotton Rose by Comparative Morpho-physiological and Transcriptome Analysis.

The cotton rose (Hibiscus mutabilis) is a plant species commonly found in tropical and subtropical regions. It is remarkably resilient to waterlogging stress; however, the underlying mechanism behind this trait is yet unknown. This study used hypoxia-tolerant "Danbanhong" (DBH) and more hypoxia-sensitive "Yurui" (YR) genotypes and compared their morpho-physiological and transcriptional responses to hypoxic conditions. Notably, DBH had a higher number of adventitious roots (20.3) compared to YR (10.0), with longer adventitious roots in DBH (18.3 cm) than in YR (11.2 cm). Furthermore, the formation of aerenchyma was 3-fold greater in DBH compared to YR. Transcriptomic analysis revealed that DBH had more rapid transcriptional responses to hypoxia than YR. Identification of a greater number of differentially expressed genes (DEGs) for aerenchyma, adventitious root formation and development, and energy metabolism in DBH supported that DBH had better morphological and transcriptional adaptation than YR. DEG functional enrichment analysis indicated the involvement of variety-specific biological processes in adaption to hypoxia. Plant hormone signaling transduction, MAPK signaling pathway and carbon metabolism played more pronounced roles in DBH, whereas the ribosome genes were specifically induced in YR. These results show that effective multilevel coordination of adventitious root development and aerenchyma, in conjunction with plant hormone signaling and carbon metabolism, is required for increased hypoxia tolerance. This study provides new insights into the characterization of morpho-physiological and transcriptional responses to hypoxia in H. mutabilis, shedding light on the molecular mechanisms of its adaptation to hypoxic environments.

A brief review of the current status of pig islet xenotransplantation.

An estimated 1.5 million Americans suffer from Type I diabetes mellitus, and its incidence is increasing worldwide. Islet allotransplantation offers a treatment, but the availability of deceased human donor pancreases is limited. The transplantation of islets from gene-edited pigs, if successful, would resolve this problem. Pigs are now available in which the expression of the three known xenoantigens against which humans have natural (preformed) antibodies has been deleted, and in which several human 'protective' genes have been introduced. The transplantation of neonatal pig islets has some advantages over that of adult pig islets. Transplantation into the portal vein of the recipient results in loss of many islets from the instant blood-mediated inflammatory reaction (IBMIR) and so the search for an alternative site continues. The adaptive immune response can be largely suppressed by an immunosuppressive regimen based on blockade of the CD40/CD154 T cell co-stimulation pathway, whereas conventional therapy (e.g., based on tacrolimus) is less successful. We suggest that, despite the need for effective immunosuppressive therapy, the transplantation of 'free' islets will prove more successful than that of encapsulated islets. There are data to suggest that, in the absence of rejection, the function of pig islets, though less efficient than human islets, will be sufficient to maintain normoglycemia in diabetic recipients. Pig islets transplanted into immunosuppressed nonhuman primates have maintained normoglycemia for periods extending more than two years, illustrating the potential of this novel form of therapy.

Revolutionizing cancer immunotherapy in solid tumor: CAR engineering and single-cell sequencing insights.

The global increase in cancer incidence presents significant economic and societal challenges. While chimeric antigen receptor-modified T cell (CAR-T) therapy has demonstrated remarkable success in hematologic malignancies and has earned FDA approval, its translation to solid tumors encounters faces significant obstacles, primarily centered around identifying reliable tumor-associated antigens and navigating the complexities of the tumor microenvironment. Recent developments in single-cell RNA sequencing (scRNA-seq) have greatly enhanced our understanding of tumors by offering high-resolution, unbiased analysis of cellular heterogeneity and molecular patterns. These technologies have revolutionized our comprehension of tumor immunology and have led to notable progress in cancer immunotherapy. This mini-review explores the progress of chimeric antigen receptor (CAR) cell therapy in solid tumor treatment and the application of scRNA-seq at various stages following the administration of CAR cell products into the body. The advantages of scRNA-seq are poised to further advance the investigation of the biological characteristics of CAR cells in vivo, tumor immune evasion, the impact of different cellular components on clinical efficacy, the development of clinically relevant biomarkers, and the creation of new targeted drugs and combination therapy approaches. The integration of scRNA-seq with CAR therapy represents a promising avenue for future innovations in cancer immunotherapy. This synergy holds the potential to enhance the precision and efficacy of CAR cell therapies while expanding their applications to a broader range of malignancies.

Corrigendum: Single-cell RNA sequencing reveals the role of phosphorylation-related genes in hepatocellular carcinoma stem cells.

[This corrects the article DOI: 10.3389/fcell.2021.734287.].

FACS-based genome-wide CRISPR screens define key regulators of DNA damage signaling pathways.

DNA damage-activated signaling pathways are critical for coordinating multiple cellular processes, which must be tightly regulated to maintain genome stability. To provide a comprehensive and unbiased perspective of DNA damage response (DDR) signaling pathways, we performed 30 fluorescence-activated cell sorting (FACS)-based genome-wide CRISPR screens in human cell lines with antibodies recognizing distinct endogenous DNA damage signaling proteins to identify critical regulators involved in DDR. We discovered that proteasome-mediated processing is an early and prerequisite event for cells to trigger camptothecin- and etoposide-induced DDR signaling. Furthermore, we identified PRMT1 and PRMT5 as modulators that regulate ATM protein level. Moreover, we discovered that GNB1L is a key regulator of DDR signaling via its role as a co-chaperone specifically regulating PIKK proteins. Collectively, these screens offer a rich resource for further investigation of DDR, which may provide insight into strategies of targeting these DDR pathways to improve therapeutic outcomes.

Construction of a hepatocytes-related and protein kinase-related gene signature in HCC based on ScRNA-Seq analysis and machine learning algorithm.

With recent advancements in single-cell sequencing and machine learning methods, new insights into hepatocellular carcinoma (HCC) progression have been provided. Protein kinase-related genes (PKRGs) affect cell growth, differentiation, apoptosis, and signaling during HCC progression, making the predictive relevance of PKRGs in HCC highly necessary for personalized medicine. In this study, we analyzed single-cell data of HCC and used the machine learning method of LASSO regression to construct PKRG prediction models in six major cell types. CDK4 and AURKB were found to be the best PKRG prognostic signature for predicting the overall survival of HCC patients (including TCGA, ICGC, and GEO datasets) in hepatocytes. Independent clinical factors were further screened out using the COX regression method, and a nomogram combining PKRGs and cancer status was created. Treatment with Palbociclib (CDK4 Inhibitor) and Barasertib (AURKB Inhibitor) inhibited HCC cell migration. Patients classified as PKRG high- or low-risk groups showed different tumor mutation burdens, immune infiltrations, and gene enrichment. The PKRG high-risk group showed higher tumor mutation burdens and gene set enrichment analysis indicated that cell cycle, base excision repair, and RNA degradation pathways were more enriched in these patients. Additionally, the PKRG high-risk group demonstrated higher infiltration levels of Naïve CD8+ T cells, Endothelial cells, M2 macrophage, and Tregs than the low-risk group. In summary, this study established the hepatocytes-related PKRG signature for prognostic stratification at the single-cell level by using machine learning algorithms in HCC and identified potential HCC treatment targets based on the PKRG signature.

Single-cell RNA analysis to identify five cytokines signaling in immune-related genes for melanoma survival prognosis.

Melanoma is one of the deadliest skin cancers. Recently, developed single-cell sequencing has revealed fresh insights into melanoma. Cytokine signaling in the immune system is crucial for tumor development in melanoma. To evaluate melanoma patient diagnosis and treatment, the prediction value of cytokine signaling in immune-related genes (CSIRGs) is needed. In this study, the machine learning method of least absolute selection and shrinkage operator (LASSO) regression was used to establish a CSIRG prognostic signature of melanoma at the single-cell level. We discovered a 5-CSIRG signature that was substantially related to the overall survival of melanoma patients. We also constructed a nomogram that combined CSIRGs and clinical features. Overall survival of melanoma patients can be consistently predicted with good performance as well as accuracy by both the 5-CSIRG signature and nomograms. We compared the melanoma patients in the CSIRG high- and low-risk groups in terms of tumor mutation burden, infiltration of the immune system, and gene enrichment. High CSIRG-risk patients had a lower tumor mutational burden than low CSIRG-risk patients. The CSIRG high-risk patients had a higher infiltration of monocytes. Signaling pathways including oxidative phosphorylation, DNA replication, and aminoacyl tRNA biosynthesis were enriched in the high-risk group. For the first time, we constructed and validated a machine-learning model by single-cell RNA-sequencing datasets that have the potential to be a novel treatment target and might serve as a prognostic biomarker panel for melanoma. The 5-CSIRG signature may assist in predicting melanoma patient prognosis, biological characteristics, and appropriate therapy.

Construction of a lipid metabolism-related risk model for hepatocellular carcinoma by single cell and machine learning analysis.

One of the most common cancers is hepatocellular carcinoma (HCC). Numerous studies have shown the relationship between abnormal lipid metabolism-related genes (LMRGs) and malignancies. In most studies, the single LMRG was studied and has limited clinical application value. This study aims to develop a novel LMRG prognostic model for HCC patients and to study its utility for predictive, preventive, and personalized medicine. We used the single-cell RNA sequencing (scRNA-seq) dataset and TCGA dataset of HCC samples and discovered differentially expressed LMRGs between primary and metastatic HCC patients. By using the least absolute selection and shrinkage operator (LASSO) regression machine learning algorithm, we constructed a risk prognosis model with six LMRGs (AKR1C1, CYP27A1, CYP2C9, GLB1, HMGCS2, and PLPP1). The risk prognosis model was further validated in an external cohort of ICGC. We also constructed a nomogram that could accurately predict overall survival in HCC patients based on cancer status and LMRGs. Further investigation of the association between the LMRG model and somatic tumor mutational burden (TMB), tumor immune infiltration, and biological function was performed. We found that the most frequent somatic mutations in the LMRG high-risk group were CTNNB1, TTN, TP53, ALB, MUC16, and PCLO. Moreover, naïve CD8+ T cells, common myeloid progenitors, endothelial cells, granulocyte-monocyte progenitors, hematopoietic stem cells, M2 macrophages, and plasmacytoid dendritic cells were significantly correlated with the LMRG high-risk group. Finally, gene set enrichment analysis showed that RNA degradation, spliceosome, and lysosome pathways were associated with the LMRG high-risk group. For the first time, we used scRNA-seq and bulk RNA-seq to construct an LMRG-related risk score model, which may provide insights into more effective treatment strategies for predictive, preventive, and personalized medicine of HCC patients.

Construction of endothelial cell-related and autophagy-related prognostic models for hepatocellular carcinoma based on single-cell data.

BACKGROUND: We aimed to develop an autophagy-related prognostic model with single-cell RNA sequencing (ScRNA-Seq) data for hepatocellular carcinoma (HCC) patients. METHODS: ScRNA-Seq datasets of HCC patients were analyzed by Seurat. The expression of genes involved in canonical and noncanonical autophagy pathways in scRNA-seq data was also compared. Cox regression was applied to construct an AutRG risk prediction model. Subsequently, we examined the characteristics of AutRG high-risk and low-risk group patients. RESULTS: Six major cell types (hepatocytes, myeloid cells, T/NK cells, B cells, fibroblast cells, and endothelial cells) were identified in the scRNA-Seq dataset. The results showed that most of the canonical and noncanonical autophagy genes were highly expressed in hepatocytes, with the exception of MAP 1LC3B, SQSTM1, MAP 1LC3A, CYBB, and ATG3. Six AutRG risk prediction models originating from different cell types were constructed and compared. The AutRG prognostic signature (GAPDH, HSP90AA1, and TUBA1C) in endothelial cells had the best overall performance for predicting the overall survival of HCC patients, with 1-year, 3-year, and 5-year AUCs equal to 0.758, 0.68, and 0.651 in the training cohort and 0.760, 0.796, and 0.840 in the validation cohort, respectively. The different tumor mutation burden, immune infiltration, and gene set enrichment characteristics of the AutRG high-risk and low-risk group patients were identified. CONCLUSION: We constructed an endothelial cell-related and autophagy-related prognostic model of HCC patients using the ScRNA-Seq dataset for the first time. This model demonstrated the good calibration ability of HCC patients and provided a new understanding of the evaluation of prognosis.

Clinical and Prognostic Value of PPIA, SQSTM1, and CCL20 in Hepatocellular Carcinoma Patients by Single-Cell Transcriptome Analysis.

Hepatocellular carcinoma (HCC) is the most malignant and poor-prognosis subtype of primary liver cancer. The scRNA-seq approach provides unique insight into tumor cell behavior at the single-cell level. Cytokine signaling in the immune system plays an important role in tumorigenesis and has both pro-tumorigenic and anti-tumorigenic functions. A biomarker of cytokine signaling in immune-related genes (CSIRG) is urgently required to assess HCC patient diagnosis and treatment. By analyzing the expression profiles of HCC single cells, TCGA, and ICGC data, we discovered that three important CSIRG (PPIA, SQSTM1, and CCL20) were linked to the overall survival of HCC patients. Cancer status and three hub CSIRG were taken into account while creating a risk nomogram. The nomogram had a high level of predictability and accuracy. Based on the CSIRG risk score, a distinct pattern of somatic tumor mutational burden (TMB) was detected between the two groups. The enrichment of the pyrimidine metabolism pathway, purine metabolism pathway, and lysosome pathway in HCC was linked to the CSIRG high-risk scores. Overall, scRNA-seq and bulk RNA-seq were used to create a strong CSIRG signature for HCC diagnosis.

The transcriptome profile of RPE cells by the fullerenol against hydrogen peroxide stress.

Age-related macular degeneration (AMD) causes central vision impairment with increased incidence. In the pathogenesis of AMD, reactive oxygen species (ROS) are associated with RPE cell apoptosis. H2O2 is an oxidative toxicant and is used to establish the AMD in vitro model. However, the mechanisms of ROS in H2O2-induced AMD are still unclear. Fullerenol, a promising antioxidant of nanomaterials, protects RPE cells from ROS attack. In addition to working as a scavenger, little is known about the antioxidant mechanism of fullerenol in RPE cells. In this study, transcriptome sequencing was performed to examine the global changes in mRNA transcripts induced by H2O2 in human ARPE-19 cells. Moreover, we comprehensively investigated the protective effects of fullerenol against H2O2-induced oxidative injury by RNA sequencing. Gene Ontology enrichment analysis showed that those pathways related to the release of positive regulation of DNA-templated transcription and negative regulation of apoptotic process were affected. Finally, we found that 12 hub genes were related to the oxidative-protection function of fullerenol. In summary, H2O2 affected these hub genes and signaling pathways to regulate the senescence of RPE cells. Moreover, fullerenol is a potent nanomaterial that protects the RPE and would be a promising approach for AMD prevention.

Human IL-17 and TNF-α Additively or Synergistically Regulate the Expression of Proinflammatory Genes, Coagulation-Related Genes, and Tight Junction Genes in Porcine Aortic Endothelial Cells.

Immune rejection is the major limitation for porcine xenograft survival in primate recipients. Proinflammatory cytokines play important roles in immune rejection and have been found to mediate the pathological effects in various clinical and experimental transplantation trials. IL-17 and TNF-α play critical pathological roles in immune disorders, such as psoriasis and rheumatoid arthritis. However, the pathological roles of human IL-17 (hIL-17) and human TNF-α (hTNF-α) in xenotransplantation remain unclear. Here we found that hIL-17 and hTNF-α additively or synergistically regulate the expression of 697 genes in porcine aortic endothelial cells (PAECs). Overall, 415 genes were found to be synergistically regulated, while 282 genes were found to be additively regulated. Among these, 315 genes were upregulated and 382 genes were downregulated in PAECs. Furthermore, we found that hIL-17 and hTNF-α additively or synergistically induced the expression of various proinflammatory cytokines and chemokines (e.g., IL1α, IL6, and CXCL8) and decreased the expression of certain anti-inflammatory genes (e.g., IL10). Moreover, hIL-17 plus hTNF-α increased the expression of IL1R1 and IL6ST, receptors for IL1 and IL6, respectively, and decreased anti-inflammatory gene receptor expression (IL10R). hIL-17 and hTNF-α synergistically or additively induced CXCL8 and CCL2 expression and consequently promoted primary human neutrophil and human leukemia monocytic cell migration, respectively. In addition, hIL-17 and hTNF-α induced pro-coagulation gene (SERPINB2 and F3) expression and decreased anti-coagulation gene (TFPI, THBS1, and THBD) expression. Additionally, hIL-17 and hTNF-α synergistically decreased occludin expression and consequently promoted human antibody-mediated complement-dependent cytotoxicity. Interestingly, hTNF-α increased swine leukocyte antigen (SLA) class I expression; however, hIL-17 decreased TNF-α-mediated SLA-I upregulation. We concluded that hIL-17 and hTNF-α likely promote the inflammatory response, coagulation cascade, and xenoantibody-mediated cell injury. Thus, blockade of hIL-17 and hTNF-α together might be beneficial for xenograft survival in recipients.

Single-Cell Landscape of Mouse Islet Allograft and Syngeneic Graft.

Islet transplantation to treat the late stage of type 1 diabetic patient (T1DM) has recently made inspiring success in clinical trials. However, most patients experience a decline in islet graft function in one to three years due to immune rejection. Although the mechanisms of immune cells, including macrophages, dendritic cells (DCs), neutrophils, natural killer cells (NKs), B cells, and T cells, that mediate immune rejection have been investigated, the overall characteristics of immune infiltrates in islet allografts and syngeneic grafts remain unclear. Single-cell RNA sequencing (scRNA-seq) has provided us with new opportunities to study the complexity of the immune microenvironment in islet transplants. In the present study, we used scRNA-seq to comprehensively analyze the immune heterogeneity in the mouse model of islet transplantation. Our data revealed T lymphocytes and myeloid cells as the main immune components of grafts 7 days post-islet transplantation, especially in allografts. Moreover, our results indicated that allogeneic islet cells were transformed into antigen-presenting cell-like cells with highly expressed MHC class I molecules and genes involved in MHC class I-mediated antigen presentation. This transformation may dramatically facilitate the interaction with cytotoxic CD8+ T cells and promote the destruction of islet allografts. Our study provides insight into the transcriptomics and diverse microenvironment of islet grafts and their impacts on immune rejection.

Current Topics of Relevance to the Xenotransplantation of Free Pig Islets.

Pig islet xenotransplantation is a potential treatment for patients with type 1 diabetes. Current efforts are focused on identifying the optimal pig islet source and overcoming the immunological barrier. The optimal age of the pig donors remains controversial since both adult and neonatal pig islets have advantages. Isolation of adult islets using GMP grade collagenase has significantly improved the quantity and quality of adult islets, but neonatal islets can be isolated at a much lower cost. Certain culture media and coculture with mesenchymal stromal cells facilitate neonatal islet maturation and function. Genetic modification in pigs affords a promising strategy to prevent rejection. Deletion of expression of the three known carbohydrate xenoantigens (Gal, Neu5Gc, Sda) will certainly be beneficial in pig organ transplantation in humans, but this is not yet proven in islet transplantation, though the challenge of the '4th xenoantigen' may prove problematic in nonhuman primate models. Blockade of the CD40/CD154 costimulation pathway leads to long-term islet graft survival (of up to 965 days). Anti-CD40mAbs have already been applied in phase II clinical trials of islet allotransplantation. Fc region-modified anti-CD154mAbs successfully prevent the thrombotic complications reported previously. In this review, we discuss (I) the optimal age of the islet-source pig, (ii) progress in genetic modification of pigs, (iii) the immunosuppressive regimen for pig islet xenotransplantation, and (iv) the reduction in the instant blood-mediated inflammatory reaction.

The Lineage Differentiation and Dynamic Heterogeneity of Thymic Epithelial Cells During Thymus Organogenesis.

Although much progress has been made recently in revealing the heterogeneity of the thymic stromal components, the molecular programs of cell lineage divergency and temporal dynamics of thymic epithelial cell (TEC) development are largely elusive. Here, we constructed a single-cell transcriptional landscape of non-hematopoietic cells from mouse thymus spanning embryonic to adult stages, producing transcriptomes of 30,959 TECs. We resolved the transcriptional heterogeneity of developing TECs and highlighted the molecular nature of early TEC lineage determination and cortico-medullary thymic epithelial cell lineage divergency. We further characterized the differentiation dynamics of TECs by clarification of molecularly distinct cell states in the thymus developing trajectory. We also identified a population of Bpifa1+ Plet1+ mTECs that was preserved during thymus organogenesis and highly expressed tissue-resident adult stem cell markers. Finally, we highlighted the expression of Aire-dependent tissue-restricted antigens mainly in Aire+ Csn2+ mTECs and Spink5+ Dmkn+ mTECs in postnatal thymus. Overall, our data provided a comprehensive characterization of cell lineage differentiation, maturation, and temporal dynamics of thymic epithelial cells during thymus organogenesis.