Interactions between brain-resident and peripheral infiltrated immune cells are thought to contribute to neuroplasticity after cerebral ischemia. However, conventional bulk sequencing makes it challenging to depict this complex immune network. Using single-cell RNA sequencing, we mapped compositional and transcriptional features of peri-infarct immune cells. Microglia were the predominant cell type in the peri-infarct region, displaying a more diverse activation pattern than the typical pro- and anti-inflammatory state, with axon tract-associated microglia (ATMs) being associated with neuronal regeneration. Trajectory inference suggested that infiltrated monocyte-derived macrophages (MDMs) exhibited a gradual fate trajectory transition to activated MDMs. Inter-cellular crosstalk between MDMs and microglia orchestrated anti-inflammatory and repair-promoting microglia phenotypes and promoted post-stroke neurogenesis, with SOX2 and related Akt/CREB signaling as the underlying mechanisms. This description of the brain's immune landscape and its relationship with neurogenesis provides new insight into promoting neural repair by regulating neuroinflammatory responses.
Humans ; Ischemic Stroke ; Brain/metabolism* ; Macrophages ; Brain Ischemia/metabolism* ; Microglia/metabolism* ; Gene Expression Profiling ; Anti-Inflammatory Agents ; Neuronal Plasticity/physiology* ; Infarction/metabolism*

Tumor progression is closely related to tumor tissue metabolism and reshaping of the microenvironment. Oral squamous cell carcinoma (OSCC), a representative hypoxic tumor, has a heterogeneous internal metabolic environment. To clarify the relationship between different metabolic regions and the tumor immune microenvironment (TME) in OSCC, Single cell (SC) and spatial transcriptomics (ST) sequencing of OSCC tissues were performed. The proportion of TME in the ST data was obtained through SPOTlight deconvolution using SC and GSE103322 data. The metabolic activity of each spot was calculated using scMetabolism, and k-means clustering was used to classify all spots into hyper-, normal-, or hypometabolic regions. CD4T cell infiltration and TGF-β expression is higher in the hypermetabolic regions than in the others. Through CellPhoneDB and NicheNet cell-cell communication analysis, it was found that in the hypermetabolic region, fibroblasts can utilize the lactate produced by glycolysis of epithelial cells to transform into inflammatory cancer-associated fibroblasts (iCAFs), and the increased expression of HIF1A in iCAFs promotes the transcriptional expression of CXCL12. The secretion of CXCL12 recruits regulatory T cells (Tregs), leading to Treg infiltration and increased TGF-β secretion in the microenvironment and promotes the formation of a tumor immunosuppressive microenvironment. This study delineates the coordinate work axis of epithelial cells-iCAFs-Tregs in OSCC using SC, ST and TCGA bulk data, and highlights potential targets for therapy.
Humans ; Carcinoma, Squamous Cell/metabolism* ; Squamous Cell Carcinoma of Head and Neck ; Mouth Neoplasms/metabolism* ; Immunosuppression Therapy ; Transforming Growth Factor beta ; Head and Neck Neoplasms ; Gene Expression Profiling ; Tumor Microenvironment

Humans ; Female ; Breast Neoplasms/genetics* ; Prognosis ; Gene Expression Profiling ; Biomarkers, Tumor/genetics* ; Gene Expression Regulation, Neoplastic/genetics* ; Receptor, ErbB-2/metabolism* ; Gene Regulatory Networks/genetics*

Dental primary afferent (DPA) neurons and proprioceptive mesencephalic trigeminal nucleus (MTN) neurons, located in the trigeminal ganglion and the brainstem, respectively, are essential for controlling masticatory functions. Despite extensive transcriptomic studies on various somatosensory neurons, there is still a lack of knowledge about the molecular identities of these populations due to technical challenges in their circuit-validated isolation. Here, we employed high-depth single-cell RNA sequencing (scRNA-seq) in combination with retrograde tracing in mice to identify intrinsic transcriptional features of DPA and MTN neurons. Our transcriptome analysis revealed five major types of DPA neurons with cell type-specific gene enrichment, some of which exhibit unique mechano-nociceptive properties capable of transmitting nociception in response to innocuous mechanical stimuli in the teeth. Furthermore, we discovered cellular heterogeneity within MTN neurons that potentially contribute to their responsiveness to mechanical stretch in the masseter muscle spindles. Additionally, DPA and MTN neurons represented sensory compartments with distinct molecular profiles characterized by various ion channels, receptors, neuropeptides, and mechanoreceptors. Together, our study provides new biological insights regarding the highly specialized mechanosensory functions of DPA and MTN neurons in pain and proprioception.
Animals ; Mice ; Neurons ; Proprioception ; Gene Expression Profiling ; Pain ; Sequence Analysis, RNA

A small proportion of mononuclear diploid cardiomyocytes (MNDCMs), with regeneration potential, could persist in adult mammalian heart. However, the heterogeneity of MNDCMs and changes during development remains to be illuminated. To this end, 12 645 cardiac cells were generated from embryonic day 17.5 and postnatal days 2 and 8 mice by single-cell RNA sequencing. Three cardiac developmental paths were identified: two switching to cardiomyocytes (CM) maturation with close CM-fibroblast (FB) communications and one maintaining MNDCM status with least CM-FB communications. Proliferative MNDCMs having interactions with macrophages and non-proliferative MNDCMs (non-pMNDCMs) with minimal cell-cell communications were identified in the third path. The non-pMNDCMs possessed distinct properties: the lowest mitochondrial metabolisms, the highest glycolysis, and high expression of Myl4 and Tnni1. Single-nucleus RNA sequencing and immunohistochemical staining further proved that the Myl4⁺Tnni1⁺ MNDCMs persisted in embryonic and adult hearts. These MNDCMs were mapped to the heart by integrating the spatial and single-cell transcriptomic data. In conclusion, a novel non-pMNDCM subpopulation with minimal cell-cell communications was unveiled, highlighting the importance of microenvironment contribution to CM fate during maturation. These findings could improve the understanding of MNDCM heterogeneity and cardiac development, thus providing new clues for approaches to effective cardiac regeneration.
Animals ; Mice ; Diploidy ; Heart ; Myocytes, Cardiac/metabolism* ; Cell Communication ; Gene Expression Profiling ; Mitochondria ; Regeneration ; Mammals/genetics*

BACKGROUND: High-grade serous ovarian cancer (HGSOC) is the biggest cause of gynecological cancer-related mortality because of its extremely metastatic nature. This study aimed to explore and evaluate the characteristics of candidate factors associated with the metastasis and progression of HGSOC. METHODS: Transcriptomic data of HGSOC patients' samples collected from primary tumors and matched omental metastatic tumors were obtained from three independent studies in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were selected to evaluate the effects on the prognosis and progression of ovarian cancer using data from The Cancer Genome Atlas (TCGA) database. Hub genes' immune landscapes were estimated by the Tumor Immune Estimation Resource (TIMER) database. Finally, using 25 HGSOC patients' cancer tissues and 10 normal fallopian tube tissues, immunohistochemistry (IHC) was performed to quantify the expression levels of hub genes associated with International Federation of Gynecology and Obstetrics (FIGO) stages. RESULTS: Fourteen DEGs, ADIPOQ , ALPK2 , BARX1 , CD37 , CNR2 , COL5A3 , FABP4 , FAP , GPR68 , ITGBL1 , MOXD1 , PODNL1 , SFRP2 , and TRAF3IP3 , were upregulated in metastatic tumors in every database while CADPS , GATA4 , STAR , and TSPAN8 were downregulated. ALPK2 , FAP , SFRP2 , GATA4 , STAR , and TSPAN8 were selected as hub genes significantly associated with survival and recurrence. All hub genes were correlated with tumor microenvironment infiltration, especially cancer-associated fibroblasts and natural killer (NK) cells. Furthermore, the expression of FAP and SFRP2 was positively correlated with the International Federation of Gynecology and Obstetrics (FIGO) stage, and their increased protein expression levels in metastatic samples compared with primary tumor samples and normal tissues were confirmed by IHC ( P = 0.0002 and P = 0.0001, respectively). CONCLUSIONS This study describes screening for DEGs in HGSOC primary tumors and matched metastasis tumors using integrated bioinformatics analyses. We identified six hub genes that were correlated with the progression of HGSOC, particularly FAP and SFRP2 , which might provide effective targets to predict prognosis and provide novel insights into individual therapeutic strategies for HGSOC.
Humans ; Female ; Ovarian Neoplasms/pathology* ; Prognosis ; Gene Expression Profiling ; Transcriptome ; Tumor Microenvironment ; Receptors, G-Protein-Coupled/therapeutic use* ; Tetraspanins/genetics* ; Protein Kinases ; Integrin beta1/therapeutic use*

OBJECTIVES: To investigate possible cross-talk genes, associated pathways, and transcription factors between chronic periodontitis (CP) and chronic obstructive pulmonary disease (COPD). METHODS: The gene expression profiles of CP (GSE10334 and GSE16134) and COPD (GSE76925) were downloaded from the GEO database. Differential expression and functional clustering analyses were performed. The protein‑protein interaction (PPI) network was constructed. The core cross-talk genes were filtered using four topological analysis algorithms and modular segmentation. Then, functional clustering analysis was performed again. RESULTS: GSE10334 detected 164 differentially expressed genes (DEGs) (119 upregulated and 45 downregulated). GSE16134 identified 208 DEGs (154 upregulated and 54 downregulated). GSE76925 identified 1 408 DEGs (557 upregulated and 851 downregulated). The PPI network included 21 nodes and 20 edges. The final screening included seven cross-talk genes: CD79A, FCRLA, CD19, IRF4, CD27, SELL, and CXCL13. Relevant pathways included primary immunodeficiency, the B-cell receptor signaling pathway, and cytokine-cytokine receptor interaction. CONCLUSIONS This study indicates the probability of shared pathophysiology between CP and COPD, and their cross-talk genes, associated pathways, and transcription factors may offer novel concepts for future mechanistic investigations.
Humans ; Chronic Periodontitis/genetics* ; Gene Regulatory Networks ; Gene Expression Profiling ; Protein Interaction Maps/genetics* ; Pulmonary Disease, Chronic Obstructive/genetics*

The aldo-keto reductase super family (AKRs) has a wide range of substrate specificity. However, the systematic identification of insect AKR gene family members has not been reported. In this study, bioinformatics methods were used to predict the phylogenetic evolution, physical and chemical properties, chromosome location, conserved motifs, and gene structure of AKR family members in Bombyx mori (BmAKR). Transcriptome data or quantitative real time polymerase chain reaction (qRT-PCR) were used to analyze the expression level of BmAKR genes during different organizational periods and silkworm eggs in different developmental states. Moreover, Western blotting was used to detect the expression level of the BmAKR in silkworm eggs. The results showed that 11 BmAKR genes were identified. These genes were distributed on 4 chromosomes of the silkworm genome, all of which had the (α/β) 8-barrel motif, and their physical and chemical characteristics were relatively similar. Phylogenetic analysis showed that the BmAKR genes could be divided into 2 subgroups (AKR1 and AKR2). Transcriptome data analysis showed that the expression of BmAKR genes were quite different in different tissues and periods. Moreover, the expression analysis of BmAKR genes in silkworm eggs showed that some genes were expressed significantly higher in nondiapause eggs than in diapause eggs; but another gene, BmAKR1-1, was expressed significantly higher in diapause eggs than in nondiapause eggs. The detection of protein level found that the difference trend of BmAKR1-1 in diapause eggs and non-diapause eggs was consistent with the results of qRT-PCR. In conclusion, BmAKR1-1 was screened out as candidates through the identification and analysis of the BmAKR genes in silkworm, which may regulate silkworm egg development is worthy of further investigation.
Animals ; Bombyx/metabolism* ; Phylogeny ; Diapause ; Genes, Insect ; Gene Expression Profiling ; Insect Proteins/metabolism*

Artemisia argyi is an important medicinal and economic plant in China, with the effects of warming channels, dispersing cold, and relieving pain, inflammation, and allergy. The essential oil of this plant is rich in volatile terpenoids and widely used in moxi-bustion and healthcare products, with huge market potential. The bZIP transcription factors compose a large family in plants and are involved in the regulation of plant growth and development, stress response, and biosynthesis of secondary metabolites such as terpenoids. However, little is known about the bZIPs and their roles in A. argyi. In this study, the bZIP transcription factors in the genome of A. argyi were systematically identified, and their physicochemical properties, phylogenetic relationship, conserved motifs, and promoter-binding elements were analyzed. Candidate AarbZIP genes involved in terpenoid biosynthesis were screened out. The results showed that a total of 156 AarbZIP transcription factors were identified at the genomic level, with the lengths of 99-618 aa, the molecular weights of 11.7-67.8 kDa, and the theoretical isoelectric points of 4.56-10.16. According to the classification of bZIPs in Arabidopsis thaliana, the 156 AarbZIPs were classified into 12 subfamilies, and the members in the same subfamily had similar conserved motifs. The cis-acting elements of promoters showed that AarbZIP genes were possibly involved in light and hormonal pathways. Five AarbZIP genes that may be involved in the regulation of terpenoid biosynthesis were screened out by homologous alignment and phylogenetic analysis. The qRT-PCR results showed that the expression levels of the five AarbZIP genes varied significantly in different tissues of A. argyi. Specifically, AarbZIP29 and AarbZIP55 were highly expressed in the leaves and AarbZIP81, AarbZIP130, and AarbZIP150 in the flower buds. This study lays a foundation for the functional study of bZIP genes and their regulatory roles in the terpenoid biosynthesis in A. argyi.
Gene Expression Profiling ; Phylogeny ; Artemisia/genetics* ; Basic-Leucine Zipper Transcription Factors/metabolism* ; Terpenes ; Gene Expression Regulation, Plant

Paeonia veitchii and P. lactiflora are both original plants of the famous Chinese medicinal drug Paeoniae Radix Rubra in the Chinese Pharmacopoeia. They have important medicinal value and great potential in the flower market. The selection of stable and reliable reference genes is a necessary prerequisite for molecular research on P. veitchii. In this study, two reference genes, Actin and GAPDH, were selected as candidate genes from the transcriptome data of P. veitchii. The expression levels of the two candidate genes in different tissues(phloem, xylem, stem, leaf, petiole, and ovary) and different growth stages(bud stage, flowering stage, and dormant stage) of P. veitchii were detected using real-time fluorescence quantitative technology(qRT-PCR). Then, the stability of the expression of the two reference genes was comprehensively analyzed using geNorm, NormFinder, BestKeeper, ΔCT, and RefFinder. The results showed that the expression patterns of Actin and GAPDH were stable in different tissues and growth stages of P. veitchii. Furthermore, the expression levels of eight genes(Pv-TPS01, Pv-TPS02, Pv-CYP01, Pv-CYP02, Pv-CYP03, Pv-BAHD01, Pv-UGT01, and Pv-UGT02) in different tissues were further detected based on the transcriptome data of P. veitchii. The results showed that when Actin and GAPDH were used as reference genes, the expression trends of the eight genes in different tissues of P. veitchii were consistent, validating the reliability of Actin and GAPDH as reference genes for P. veitchii. In conclusion, this study finds that Actin and GAPDH can be used as reference genes for studying gene expression levels in different tissues and growth stages of P. veitchii.
Real-Time Polymerase Chain Reaction/methods* ; Paeonia/genetics* ; Actins/genetics* ; Reproducibility of Results ; Transcriptome ; Glyceraldehyde-3-Phosphate Dehydrogenases/genetics* ; Reference Standards ; Gene Expression Profiling/methods*