Search

Objective To explore the gene expression profile in paraspinal muscle degeneration (PMD) and identify key ferroptosis genes. Methods RNA sequencing was performed on paraspinal muscle tissue of 3 normal and 3 PMD patients respectively to obtain differentially expressed genes. Through protein-protein interaction (PPI) and gene functional enrichment analysis, the intersection of ferroptosis genes was identified to identify key hub genes associated with ferroptosis. The diagnostic value for PMD disease was analyzed by receiver operating characteristic (ROC) curves. Results A total of 292 differentially expressed genes were identified in PMD. Among them, 125 genes were significantly downregulated and 167 genes were significantly upregulated. Bioinformatics analysis revealed that 14 differentially expressed genes were associated with ferroptosis. Among them, ferroptosis genes MUC1, ATF3 and CDKN1A were key hub genes with good specificity and sensitivity for diagnosing PMD. Functional enrichment analysis revealed that they may mediate the occurrence and progression of PMD by regulating cell apoptosis, ferroptosis and skeletal muscle tissue development and differentiation. Conclusion Ferroptosis genes MUC1, ATF3 and CDKN1A can serve as biomarkers for diagnosing PMD, providing theoretical basis for decoding the pathological mechanism of PMD and developing new drugs. [ABSTRACT FROM AUTHOR]

Objective To summarize the clinical features and treatment of genetic epilepsy with febrile seizures plus (GEFS+) caused by SCN1A gene mutations, to improve the understanding of SCN1A gene mutations, and to provide a basis for the clinical diagnosis and treatment of this disease. Methods Clinical data were collected from a child with GEFS+ and her family members who attended Department of Neurology in our hospital, and the third-generation gene sequencing technique was used to perform genetic analysis and validation of peripheral blood DNA of the child and her parents and aunt. Results A girl, aged 7 years, was admitted to the hospital due to "paroxysmal seizures lasting for more than 4 years", with loss of consciousness during seizures and more seizures during pyrexia. The seizures increased after the use of oxcarbazepine and were controlled after the admi-nistration of sodium valproate. Her father had a history of seizures, while her aunt had no clinical symptoms. Third-generation gene sequencing showed the same heterozygous mutation, c.4787G>A(p.R1596H), in the SCN1A gene in the child and her father and aunt. Conclusion The onset of GEFS+ is associated with SCN1A gene mutations, but there is heterogeneity in clinical manifestations among patients in one family, and sodium channel blockers as anti-epileptic drugs will aggravate epileptic seizures in GEFS+.

Objectives: The present study used single-cell RNA sequencing (scRNA-seq) to characterize the cellular composition of ovarian carcinosarcoma (OCS) and identify its molecular characteristics., Methods: scRNA-seq was performed in resected primary OCS for an in-depth analysis of tumor cells and the tumor microenvironment. Immunohistochemistry staining was used for validation. The scRNA-seq data of OCS were compared with those of high-grade serous ovarian carcinoma (HGSOC) tumors and other OCS tumors., Results: Both malignant epithelial and malignant mesenchymal cells were observed in the OCS patient of this study. We identified four epithelial cell subclusters with different biological roles. Among them, epithelial subcluster 4 presented high levels of breast cancer type 1 susceptibility protein homolog ( BRCA1 ) and DNA topoisomerase 2-α ( TOP2A ) expression and was related to drug resistance and cell cycle. We analyzed the interaction between epithelial and mesenchymal cells and found that fibroblast growth factor (FGF) and pleiotrophin (PTN) signalings were the main pathways contributing to communication between these cells. Moreover, we compared the malignant epithelial and mesenchymal cells of this OCS tumor with our previous published HGSOC scRNA-seq data and OCS data. All the epithelial subclusters in the OCS tumor could be found in the HGSOC samples. Notably, the mesenchymal subcluster C14 exhibited specific expression patterns in the OCS tumor, characterized by elevated expression of cytochrome P450 family 24 subfamily A member 1 ( CYP24A1 ), collagen type XXIII α1 chain ( COL23A1 ), cholecystokinin ( CCK ), bone morphogenetic protein 7 ( BMP7 ), PTN , Wnt inhibitory factor 1 ( WIF1 ), and insulin-like growth factor 2 ( IGF2 ). Moreover, this subcluster showed distinct characteristics when compared with both another previously published OCS tumor and normal ovarian tissue., Conclusions: This study provides the single-cell transcriptomics signature of human OCS, which constitutes a new resource for elucidating OCS diversity.

Objectives: RNA sequencing (RNA-seq) and bioinformatic analysis were combined and used to explore the anti-inflammatory effects and mechanisms of naringenin (Nar) in lipopolysaccharide (LPS)-stimulated human periodontal ligament stem cells (hPDLSCs)., Methods: Cell counting kit-8, quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR), and enzyme-linked immunosorbent assay (ELISA) were adopted to detect the effects of Nar on the proliferation and expression of inflammatory factors in LPS-stimulated hPDLSCs, screening for the optimal anti-inflammatory concentration of Nar. Differentially expressed genes (DEGs) were screened using |log2FC|≥1 and P ≤0.05 as criteria. Volcano plot analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, the String database, and the MCODE module of Cytoscape were utilized to select core genes and enriched pathways. The effects on the nuclear factor κB (NF-κB) signaling pathway were verified using ELISA, qRT-PCR, and Western blot., Results: Appropriate concentrations of Nar could alleviate the expression of inflammatory factors and promote the proliferation of hPDLSCs stimulated by LPS. The best anti-inflammatory effect was achieved with 20 μmol/L Nar. RNA-seq showed significant enrichment of inflammation-related signaling pathways. The anti-inflammatory effect of Nar was mediated by inhibiting the NF-κB signaling pathway, similar to the effect of the NF-κB inhibitor BAY 11-7802., Conclusions: Nar could exert its anti-inflammatory effects by inhibiting the NF-κB signaling pathway, making it a potential therapeutic option for the adjuvant treatment of periodontitis.

Over 170 chemical modifications have been discovered in various types of ribonucleic acids (RNAs), including messenger RNA (mRNA), ribosomal RNA (rRNA), transfer RNA (tRNA), and small nuclear RNA (snRNA). These RNA modifications play crucial roles in a wide range of biological processes such as gene expression regulation, RNA stability maintenance, and protein translation. RNA modifications represent a new dimension of gene expression regulation known as the "epitranscriptome". The discovery of RNA modifications and the relevant writers, erasers, and readers provides an important basis for studies on the dynamic regulation and physiological functions of RNA modifications. Owing to the development of detection technologies for RNA modifications, studies on RNA epitranscriptomes have progressed to the single-base resolution, multilayer, and full-coverage stage. Transcriptome-wide methods help discover new RNA modification sites and are of great importance for elucidating the molecular regulatory mechanisms of epitranscriptomics, exploring the disease associations of RNA modifications, and understanding their clinical applications. The existing RNA modification sequencing technologies can be categorized according to the pretreatment approach and sequencing principle as direct high-throughput sequencing, antibody-enrichment sequencing, enzyme-assisted sequencing, chemical labeling-assisted sequencing, metabolic labeling sequencing, and nanopore sequencing technologies. These methods, as well as studies on the functions of RNA modifications, have greatly expanded our understanding of epitranscriptomics. In this review, we summarize the recent progress in RNA modification detection technologies, focusing on the basic principles, advantages, and limitations of different methods. Direct high-throughput sequencing methods do not require complex RNA pretreatment and allow for the mapping of RNA modifications using conventional RNA sequencing methods. However, only a few RNA modifications can be analyzed by high-throughput sequencing. Antibody enrichment followed by high-throughput sequencing has emerged as a crucial approach for mapping RNA modifications, significantly advancing the understanding of RNA modifications and their regulatory functions in different species. However, the resolution of antibody-enrichment sequencing is limited to approximately 100-200 bp. Although chemical crosslinking techniques can achieve single-base resolution, these methods are often complex, and the specificity of the antibodies used in these methods has raised concerns. In particular, the issue of off-target binding by the antibodies requires urgent attention. Enzyme-assisted sequencing has improved the accuracy of the localization analysis of RNA modifications and enables stoichiometric detection with single-base resolution. However, the enzymes used in this technique show poor reactivity, specificity, and sequence preference. Chemical labeling sequencing has become a widely used approach for profiling RNA modifications, particularly by altering reverse transcription (RT) signatures such as RT stops, misincorporations, and deletions. Chemical-assisted sequencing provides a sequence-independent RNA modification detection strategy that enables the localization of multiple RNA modifications. Additionally, when combined with the biotin-streptavidin affinity method, low-abundance RNA modifications can be enriched and detected. Nevertheless, the specificity of many chemical reactions remains problematic, and the development of specific reaction probes for particular modifications should continue in the future to achieve the precise localization of RNA modifications. As an indirect localization method, metabolic labeling sequencing specifically localizes the sites at which modifying enzymes act, which is of great significance in the study of RNA modification functions. However, this method is limited by the intracellular labeling of RNA and cannot be applied to biological samples such as clinical tissues and blood samples. Nanopore sequencing is a direct RNA-sequencing method that does not require RT or the polymerase chain reaction (PCR). However, challenges in analyzing the data obtained from nanopore sequencing, such as the high rate of false positives, must be resolved. Discussing sequencing analysis methods for various types of RNA modifications is instructive for the future development of novel RNA modification mapping technologies, and will aid studies on the functions of RNA modifications across the entire transcriptome.

Objective: To clarify the potential correlation between biological changes of meninges in periodontitis mice and cognitive impairment by analyzing the biological changes of meninges in periodontitis mice using single-cell RNA sequencing. Methods: Thirty C57BL/6 mice were divided into two groups by using random number table method (15 mice in each group). Mice in the control group were locally administered 2% carboxyl methyl cellulose (CMC) without Porphyromonas gingivalis (Pg) on both buccal sides. A mixture of Pg W83 and 2% CMC was applied on both buccal sides in the experimental group mice three times a week, lasting for 16 weeks in total. The absorption of alveolar bone, locomotor activity and cognitive function, the activation of microglia and astrocytes in the cortex were observed and assessed. The mRNA expression levels of Occludin in meninges and brain were detected in two groups. Single-cell RNA sequencing data of meninges were processed by uniform manifold approximation and projection (UMAP). Differential genes expressions of endothelial cells were processed by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. In addition, real-time fluorescence quantitative PCR (RT-qPCR) was used to verify the expressions of transcription activating factor 3 (Atf3) and apolpoprotein L domain-containing 1 (Apold 1). Results: Methylene blue staining found the distances of buccal and palatal cement-enamel junction-alveolar bone crest in experimental mice [(185.60±17.60), (206.90±13.37) μm] increased significantly compared with the control group [(135.33±9.57), (163.05±14.98) μm] ( t= 5.02, P= 0.002; t= 4.37, P= 0.005). Open field experiment showed the total distance and average speed of mice in the experimental group [(971.88±164.57) cm, (3.25±0.55) cm/s] were not statistically significant compared with the control group [(914.24±278.81) cm, (3.05±0.93) cm/s] ( t= 0.65, P= 0.525; t= 0.65, P= 0.520). The recognition index of the experimental group [(48.02±16.92) %] was lower than the control group [(66.27±17.90) %] ( t= 2.40, P= 0.027) by novel object recognition tests. Compared with the control group [(63.56±11.88) %], the alternation of experimental group [(50.99±14.17) %] was significantly decreased in Y maze tests ( t= 2.33, P= 0.030). Immunohistochemistry results showed microglia and astrocytes were activated in the cortex of experimental mice. Compared with the control group (1.02±0.25, 1.04±0.31), the relative mRNA expressions of Occludin decreased significantly in the meninges and brain of periodontitis mice, respectively (0.61±0.10, 0.64±0.20) ( t= 3.47, P= 0.010; t= 2.66, P= 0.024). By single-cell RNA sequencing, meninges cells were divided into 11 types, such as endothelial cells, fibroblasts, immune cells and so on. Endothelial cells were the main cell types in meninges [the control group: 26.47% (1 589/6 004), the experimental group: 26.26% (807/3 073)]. Compared with the control group [5.56% (334/6 004)], the percentage of granulocytes increased in the periodontitis mice [11.65% (358/3 073)]. Using clustering analysis to further focus on endothelial cells, GO enrichment analysis revealed differential genes were mainly related to angiogenesis, cell adhesion, apoptosis and so on. KEGG enrichment analysis revealed that differential genes were related to signaling pathways of interleukin-17, relaxin and so on. The relative mRNA expressions of Atf3 and Apold1 in meninges of periodontitis mice (0.42±0.24, 0.54±0.27) were significantly lower than the control group (1.03±0.26, 1.02±0.23) ( t= 3.88, P= 0.005; t= 3.02, P= 0.017). Conclusions: The mice chronically infected with Pg W83 occurred memory impairment, neuroinflammation and changes of barrier function. In the meninges of periodontitis mice, there were infiltration of immune cells and down-regulation expressions of Atf3 and Apold1 by single-cell RNA sequencing. Meningeal immunity and changes of barrier function may play an important role in the cognitive impairment caused by periodontitis.

Objective: To analyze the differential genes and related signaling pathways of microglia subpopulations in Parkinson's disease (PD) -like mouse brains induced by paraquat (PQ) based on single-cell RNA sequencing, and provide clues to elucidate the mechanism of PQ-induced PD-like changes in the brain of animals. Methods: In September 2021, six male 6-week-old C57BL/6 mice were randomly divided into control group and experimental group (three mice in each group) . The mice were injected with saline, 10.0 mg/kg PQ intraperitoneally, once every three days, and 10 consecutive injections were used for modeling. After infection, the brains of mice were taken and 10×Genomics single-cell RNA sequencing was performed. Microglia subpopulations were screened based on gene expression characteristics, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed. The differential genes of microglia subpopulations between the experimental group and control group were further screened, and functional enrichment analysis was performed using bioinformatics tools. Mouse microglia (BV2 cells) were treated with 0, 60, 90 μmol/L PQ solution, respectively. And real-time fluorescence quantitative PCR experiments were conducted to validate the expressions of differential genes hexokinase 2 (Hk2) , ATPase H+ Transporting V0 Subunit B (Atp6v0b) and Neuregulin 1 (Nrg1) . Results: Cluster 7 and Cluster 20 were identified as microglia subpopulations based on the signature genes inositol polyphosphate-5-phosphatase d, Inpp5d (Inpp5d) and transforming growth factor beta receptor 1 (Tgfbr1) , and they reflected the microglia-activated M2 phenotype. The bioinformatics analysis showed that the characteristic genes of identified microglia subpopulations were enriched in endocytosis. In terms of molecular function, it mainly enriched in transmembrane receptor protein kinase activity and cytokine binding. The up-regulated genes of Cluster 7 were mainly enriched in lysosomal pathway, endocytosis pathway, and down-regulated genes were mainly enriched in neurodegenerative disease and other signaling pathways. The up-regulated genes of Cluster 20 were mainly enriched in signaling pathways related to PD, and down-regulated genes were mainly enriched in cyclic adenosine 3', 5'-monophosphate (cAMP) signaling pathways, neurological development, synaptic function and other signaling pathways. The results of real-time fluorescence quantitative PCR showed that the expressions of Hk2 mRNA and Atp6v0b mRNA increased and the expression of Nrg1 mRNA decreased in the 90 μmol/L PQ-treated BV2 cells compared with the 0 μmol/L, and the differences were statistically significant ( P <0.05) . Conclusion: Microglia are activated in the PQ-induced PD-like mouse model and polarized toward the M2 phenotype. And their functions are associated with lysosomal (endocytosis) , synaptic functions and the regulation of PD-related pathways.

Objectives: Abnormal programmed cell death in immune cells is associated with autoimmune diseases, but the patterns of programmed cell death in systemic lupus erythematosus (SLE) and especially lupus nephritis (LN) remain unclear. This study aims to explore the association between SLE, LN, and immune cell death patterns., Methods: Bulk RNA sequencing (bulk RNA-seq) and single-cell RNA sequencing (scRNA-seq) data were downloaded from the Gene Expression Omnibus (GEO) database. Bioinformatic analysis was conducted to explore the expression levels of genes related to 3 cell death patterns in peripheral blood mononuclear cells of SLE patients. Key cell subsets involved in the imbalance of cell death patterns were identified through scRNA-seq. Immunofluorescence was used to detect the expression levels of receptor interacting serine/threonine kinase 3 (RIPK3), mixed-lineage kinase domain-like protein (MLKL), phosphorylated MLKL (pMLKL), caspase 1 (CASP1), CD1c molecule (CD1C), C-type lectin domain containing 9A (CLEC9A), and X-C motif chemokine receptor 1 (XCR1) in dendritic cells (DC). scRNA-seq was performed on kidney tissues collected from LN patients and healthy controls (HC) at the Third Xiangya Hospital of Central South University, followed by bioinformatic analysis to identify key cell subsets involved in the imbalance of cell death patterns. Pseudotime analysis and ligand-receptor analysis were used to explore the differentiation direction and cell communication of different DC subsets. Transient transfection was used to transfect RAW264.7 cells with empty plasmid, empty plasmid+dsDNA (HSV-DNA), empty plasmid+200 μmol/L tert-butyl hydroperoxide (TBHP), stimulator of interferon genes (STING) shRNA plasmid, STING shRNA plasmid+dsDNA (HSV-DNA), and STING shRNA plasmid+200 μmol/L TBHP. Annexin V-mCherry and SYTOX Green staining were used to detect cell death in each group. Western blotting was used to detect the activation of CASP1, gasdermin D (GSDMD), RIPK3, and MLKL in each group., Results: Bioinformatic analysis showed an imbalance in 3 cell death patterns in SLE and LN patients: Pro-inflammatory pyroptosis and necroptosis were activated, while anti-inflammatory apoptosis was inhibited. The key cell subsets involved were DC subsets, particularly focusing on CLEC9A +cDC1. Immunofluorescence results showed that the expression levels of RIPK3, MLKL, and CASP1 in DCs were higher in the SLE group compared to the HC group. pMLKL and CASP1 expression levels in renal cDC1 marked by CLEC9A and XCR1 were higher in the LN group than in the HC group. Pseudotime analysis and ligand-receptor analysis suggested that the CLEC9A +cDC1 subset in LN kidney tissues originated from peripheral circulation. Annexin V-mCherry and SYTOX Green staining results showed that the number of dead cells decreased in the STING shRNA transfection group compared to the empty plasmid group in RAW264.7 cells. Western blotting results showed that the activation of CASP1, GSDMD, RIPK3, and MLKL was decreased in the STING shRNA transfection group compared to the empty plasmid group., Conclusions: This study provides novel insights into the role of CLEC9A +cDC1 in the imbalance of cell death patterns in SLE and LN.

Objective: To analyze the types and functions of CD34 + cells in full-thickness skin defect wounds of normal mice and diabetic mice by single-cell RNA sequencing. Methods: This study was an experimental study. The CD34 + cell lineage tracing mouse was produced, and the visualization of CD34 + cells under the fluorescent condition was realized. Six male CD34 + cell lineage tracing mice aged 7-8 weeks (designated as diabetic group) were intraperitoneally injected with streptozotocin to establish a diabetic model, and full-thickness skin defect wounds were prepared on their backs when they reached 13 weeks old. Another 6 male CD34 + cell lineage tracing mice aged 13 weeks (designated as control group) were also subjected to full-thickness skin defect wounds on their backs. On post-injury day (PID) 4, wound tissue was collected from 3 mice in control group and 2 mice in diabetic group, and digested to prepare single-cell suspensions. CD34 + cells were screened using fluorescence-activated cell sorting, followed by single-cell RNA sequencing. The Seurat 4.0.2 program in the R programming language was utilized for dimensionality reduction, visualization, and cell clustering analysis of CD34 + cell types, and to screen and annotate the marker genes for each CD34 + cell subpopulation. Kyoto encyclopedia of genes and genomes (KEGG) and gene ontology (GO) enrichment analysis was performed to analyze the differentially expressed genes (DEGs) of CD34 + fibroblasts (Fbs), smooth muscle cells (SMCs), keratinocytes (KCs), and chondrocyte-like cells (CLCs) in the wound tissue of two groups of mice for exploring cellular functions. Results: On PID 4, CD34 + cells in the wound tissue of both groups of mice were consisted of 7 cell types, specifically endothelial cells, Fbs, KCs, macrophages, T cells, SMCs, and CLCs. Among these, Fbs were further classified into 5 subpopulations. Compared with those in control group, the proportions of CD34 + endothelial cells, Fbs subpopulation 1, Fbs subpopulation 4, KCs, and CLCs in the wound tissue of mice were increased in diabetic group, while the proportions of CD34 + Fbs subpopulation 2, Fbs subpopulation 3, and SMCs were decreased. The marker genes for annotating CD34 + CLCs, endothelial cells, Fbs subpopulation 1, Fbs subpopulation 2, Fbs subpopulation 3, Fbs subpopulation 4, Fbs subpopulation 5, KCs, macrophages, SMCs, and T cells were respectively metastasis-associated lung adenocarcinoma transcript 1, fatty acid binding protein 4, Gremlin 1, complement component 4B, H19 imprinted maternally expressed transcript, Dickkopf Wnt signaling pathway inhibitor 2, fibromodulin, keratin 5, CD74 molecule, regulator of G protein signaling 5, and inducible T-cell co-stimulator molecule. KEGG and GO enrichment analysis revealed that, compared with those in control group, DEGs with significant differential expression (SDE) in CD34 + Fbs from the wound tissue of mice in diabetic group on PID 4 were significantly enriched in terms related to inflammatory response, extracellular matrix (ECM) organization, regulation of cell proliferation, and aging (with P values all <0.05), DEGs with SDE in CD34 + SMCs were significantly enriched in terms related to cell migration, apoptotic process, positive regulation of transcription, and phagosome (with P values all <0.05), DEGs with SDE in CD34 + KCs were significantly enriched in terms related to mitochondrial function, transcription, and neurodegenerative diseases (with P values all <0.05), and DEGs with SDE in CD34 + CLCs were significantly enriched in terms related to rhythm regulation, ECM, and viral infection (with P values all <0.05). Conclusions: CD34 + cells display high heterogeneity in the healing process of full-thickness skin defect wounds in both normal mice and diabetic mice. The significantly enriched functions of DEGs with SDE in CD34 + cell subpopulations in the wound tissue of the two mouse groups are closely related to the wound healing process.

To investigate the methods for identifying Actinotignum schaalii and to analyze its biological characteristics in order to provide a basis for the diagnosis of urinary tract infections. Methods The urine sample of a patient with urinary tract infection hospitalized in the First Teaching Hospital of Tianjin University of Traditional Chinese Medicine in January of 2020 was collected. Bacterial culture and Gram staining were used to observe the morphological characteristics. Matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) and 16S rRNA sequencing were used for the identification. The obtained sequences were compared with those in GenBank database for phylogenetic analysis. Results The case was the first reported in China for a urinary tract infection caused by Actinotignum schaalii. MALDI-TOF MS and 16S rRNA sequencing identified the strain as Actinotignum schaalii, belonging to the genera Actinotignum, which was facultative anaerobic and Gram-positive bacterium. It was considered as a new type of pathogenic bacteria causing urinary tract infection. Conclusion Actinotignum schaalii may be a new pathogen that causes urinary tract infection. [ABSTRACT FROM AUTHOR]

Objective Single-cell RNA sequencing (scRNA-Seq) and experimental verifies core genes of dendritic cells in chronic obstructive pulmonary disease (COPD). Methods scRNA-seq data GSE173896 and chip data GSE38974 were extracted from the Gene Expression Omnibus (GEO) database. GSE173896 was used to perform quality control, batch correction, dimensionality reduction clustering, cell type annotation and dendritic cell differentially expressed genes (DC-DEGs) identification. DEGs from the analysis of GSE38974 were intersected with DC-DEGs to obtain the common DC-DEGs. The diagnostic efficacy of the common DC-DEGs for COPD and their enrichment analysis were conducted. The correlation of the common DC-DEGs with activated dendritic cell (DCs), plasmacytoid dendritic cell (pDCs) and type 17 T helper(Th17) cells were analyzed. The mRNA expression level of the common DC-DEGs in the lung tissue of emphysema mice was verified. Results From GSE173896, 18 DC-DEGs were obtained between groups and from GSE38974, 646 DEGs were obtained. The intersection of the two resulted in 3 common DC-DEGs, including interleukin 1 receptor antagonist 1 (IL1RN), S100 calcicum-binding protein A8 (S100A8) and S100A9. Their respective area under curve (AUC) values were 0.841, 0.804 and 0.966. The GO and KEGG enrichment analysis mainly concentrated on chronic inflammatory response, collagen-containing extracellular matrix, receptor for advanced glycation end products (RAGE) binding, Toll-like receptor (TLR) binding and interleukin 17 (IL-17) signaling pathway. IL1RN, S100A8 and S100A9 were positively correlated with activated DCs, pDCs and Th17 cells. The results showed that the mRNA relative expression levels of IL1RN, S100A8 and S100A9 were up-regulated in the lung tissue of emphysema mice. Conclusion IL1RN, S100A8 and S100A9 may be the core genes of DCs in the pathogenesis of COPD, which potentially provide targets and a theoretical basis for subsequent COPD immunotherapy.

To explore the biological characteristics related to the pathogenesis and severity of respiratory syncytial virus (RSV) bronchiolitis by RNA sequencing of white blood cells in children with RSV bronchiolitis. This study is a case-control study. A total of 87 children diagnosed with bronchiolitis and RSV antigen positive and/or RSV nucleic acid positive in the pediatric respiratory department of the Second Affiliated Hospital of Wenzhou Medical University from October 2019 to April 2022 were selected as the case group. The case group was divided into three groups based on the condition: mild, moderate, and severe, and there were two groups according to the presence or absence of atopic symptoms: the atopic group and the non - atopic group, forty healthy children in the same period were selected as the control group. The whole blood leukocyte RNA of the children in the case group and the control group was extracted for RNA sequencing, and the data were analyzed to obtain differentially expressed genes (DEGs). Then, the immunobiological pathways and genes related to the pathogenesis, disease condition, and atopy were screened through Gene Ontology (GO) annotation, Kyoto Gene and Genome Encyclopedia (KEGG) annotation, and protein interaction network (PPI) construction methods. Construct the weighted gene co-expression network analysis (WGCNA) module to identify potential biological indicators related to disease severity.Compared with the control group, the case group had a total of 1 782 DEGs, including 1 586 upregulated genes and 196 downregulated genes. The GO pathway enrichment of DEGs is mainly enriched in molecular functions such as peroxidase activity and oxidoreductase activity. In the cytological components, it is mainly enriched in cytoplasmic vesicle lumen and secretory granule lumen. In biological processes, it is mainly enriched in processes such as neutrophil activation involved in immune responses, neutrophil degranulation, and neutrophil activation. KEGG analysis is mainly concentrated in the signal pathway of the viral protein interaction with cytokine and cytokine receptor. A PPI network was constructed to screen four genes at the core position, including CCL2, IL-10, MMP9 and JUN . The DEGs obtained by comparing different disease groups with the control group are mainly enriched in retrograde endocannabinoid signaling and cell apoptosis pathways. WGCNA analysis showed that the brown module related to oxygen saturation was most closely related to the disease, and its gene was mainly enriched in the RNA helicase retinoic acid inducible gene-I (RIG-I) like receptor signal pathway. There are 230 specific DEGs in the atopic group and 444 in the non - atopic group. KEGG enrichment analysis results show that both groups are enriched to NF-κB signaling pathway, the characteristic does not cause significant changes in immune response and transcriptome characteristics in children with RSV bronchiolitis. In conclusion, neutrophil activation, degranulation pathway and signal pathway of interaction between viral protein and cytokine and cytokine receptor are involved in the immune response of RSV bronchiolitis host. CCL2, IL-10, MMP9 and JUN genes may be associated with the pathogenesis. They might be potential biomarkers related to disease severity in RIG-I like receptors, cell apoptosis, and endogenous cannabinoid related signaling pathways.

Objective: To investigate the cell composition and the transcriptional characteristics in microenvironments of hepatic tissues in mice at late stage of Echinococcus multilocularis infection at a single-cell level., Methods: Peri-lesion and paired distal hepatic specimens were collected from two BALB/c mice (6 to 8 weeks old) infected with E. multilocularis for single-cell RNA sequencing. The Seurat package in the R software was employed for quality control of data, multi-sample integration and correction of batch effects, and uniform manifold approximation and projection (UMAP) algorithm was used for cell clustering. Cell types were annotated using classical marker genes. Differentially expressed genes were screened in each cell type through differential gene expression analysis, and the biological roles of cells were predicted using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses., Results: A total of 43 710 cells from peri-lesion and distal hepatic tissues of E. multilocularis -infected mice were analyzed, and were classified into 11 cell types, including neutrophils, T cells, macrophages, granulocyte-monocyte progenitor cells, B cells, plasma cells, basophils, hepatic stellate cells, endothelial cells, hepatocytes, and platelets. T cells were the largest population of immune cells in the microenvironment of hepatic tissues, including five CD4 + T cell subsets, two CD8 + T cell subsets and phosphoantigen-reactive γδT cells. The proportions of CD4 + helper T cells and cytotoxic CD4 + T cells decreased and the proportion of T helper 2 (Th2) cells increased in peri-lesion tissues relative to distal hepatic tissues. In addition, the differentially expressed genes in Th2 cells were associated with negative regulation of the immune system, and the highly expressed genes in cytotoxic CD4 + T cells correlated with activation of the immune system., Conclusions: Single-cell RNA sequencing deciphers the cell composition and distribution in microenvironments of hepatic tissues from mice infected with E. multilocularis , and the increased proportion of Th2 cells in peri-lesion hepatic tissues may be associated with formation of immunosuppressive microenvironments.

Objective: To depict the cell landscape and molecular biological characteristics of human intrauterine adhesion (IUA) so as to better understand its immune microenvironment and provide new inspirations for clinical treatment., Methods: Four patients with IUA who underwent hysteroscopic treatment at Dongguan Maternal and Child Health Care Hospital from February 2022 to April 2022 were selected as the study subjects. Hysteroscopy was used to collect the tissues of IUA, which were graded based on the patient's medical history, menstrual history and status of IUA. Library construction, sequencing, single cell data comparison and gene expression matrix construction were carried out in strict accordance with the single cell RNA sequencing process. Thereafter, the UMAP dimension reduction analysis of cell population and genetic analysis were carried out based on the cell types., Results: A total of 27 511 cell transcripts were obtained from four moderately graded IUA tissue samples and assigned to six cell lineages including T cells, mononuclear phagocytes, epithelial cells, fibroblasts, endothelial cells and erythrocytes. Compared with normal uterine tissue cells, the four samples showed different cell distribution, and the proportions of mononuclear phagocytes and T cells in sample IUA0202204 were significantly increased, suggesting a strong cellular immune response., Conclusion: The cell diversity and heterogeneity of moderate IUA tissues have been described. Each cell subgroup has unique molecular characteristics, which may provide new clues for further study of the pathogenesis of IUA and heterogeneity among the patients.

Objective: Immunotherapy has brought significant clinical benefits to a subset of patients, but has thus far been disappointing in the treatment of immunologically "cold" tumors. Existing biomarkers that can precisely identify these populations are insufficient. In this context, a potential cold tumor microenvironment (TME) marker FARSB was investigated to reveal its impact on TME and patients' response to immunotherapy across pan-cancer., Methods: The expression levels and mutational landscape of FARSB in pan-cancer were investigated. Kaplan-Meier and univariate Cox regression analyses were applied to analyze the prognostic significance of FARSB . Pathways affected by FARSB were investigated by gene set enrichment and variation analysis. The relationship between FARSB expression and immune infiltration was examined using the TIMER2 and R packages. Single-cell RNA sequencing (scRNA-seq) data of several cancer types from GSE72056, GSE131907, GSE132465, GSE125449 and PMID32561858 were analyzed to validate the impact of FARSB on the TME. The predictive effect of FARSB on immunotherapy efficacy was explored in 3 immune checkpoint inhibitors (ICIs)- treated cohorts (PMID32472114, GSE176307, and Riaz2017)., Results: FARSB expression was significantly higher in 25 tumor tissues than in normal tissues and was associated with poor prognosis in almost all tumor types. FARSB expression exhibited a strong association with several DNA damage repair pathways and was significantly associated with TP53 mutation in lung adenocarcinoma ( P < 0.0001, OR=2.25). FARSB characterized a typical immune desert TME and correlated with impaired expression of chemokines and chemokines receptors. Large-scale scRNA-seq analysis confirmed the immunosuppressive role of FARSB and revealed that FARSB potentially shapes the cold TME by impeding intercellular interactions. In 3 ICI-treated cohorts, FARSB demonstrated predictive value for immunotherapy., Conclusion: This study provides a pan-cancer landscape of the FARSB gene by integrated single-cell and bulk DNA sequencing analysis and elucidates its biological function to promote DNA damage repair and construct the immune desert TME, suggesting the potential value of FARSB as a novel marker for stratifying patients with poor immunotherapeutic benefits and "cold" TME.

Objectives: Uterine fluid RNA can be used as a test for endometrial receptivity, but there is still no noninvasive sampling method available. The polyvinyl alcohol (PVA) formaldehyde absorbent sponge, a medical bio-absorbent sponge with good water absorption and biophilic properties, can be used to develop a new noninvasive endometrial fluid sampler. This study aims to investigate the toxicity of PVA acetal absorbent sponges on endometrial epithelial cells and its effect on RNA sequencing (RNA-Seq)., Methods: The experimental group using PVA formaldehyde absorbent sponge was prepared into 0.005%, 0.01% and 0.02% (w/v) suspension, and 0.01%, 0.05% and 0.1% (v/v) extract groups. The control group was only the complete culture medium. Nothing was added to the blank group. In vitro cytotoxicity assay was used to evaluate the survival rate of cells. Eight patients underwent in vitro fertilization treatment in the Reproductive Center of Xiangya Hospital, Central South University from November 2019 to January 2020. The uterine fluid of each patient was aspirated. The experimental group was inhaled with sterile PVA formaldehyde absorbent sponge and then immersed RNA-later solution. The control group was directly injected into the same amount of RNA-later solution. RNA-seq and data analysis was performed later., Results: The vitro cytotoxicity assay showed that in suspension groups, there was no significance difference in cell survival between different co-culture time in 0.005% group (P=0.255). In the 0.01% and 0.02% group, there was no difference at each incubation time within 12 h (all P>0.05), but the cell survival rate was decreased at 24 h compared with 0 h (P<0.01, P<0.05). At the same co-culture time, the cell survival of the 3 concentration gradient groups were significantly lower than that of the control group (all P<0.05). The cell viability of the 0.005% concentration group was decreased less than 30% at 24 h, the 0.01% concentration group decreased more than 30% at 12 h, and the 0.02% concentration group was decreased more than 30% at 0 h. For extract groups, there was no significant difference in the survival rate within 6 h in 0.01% concentration group (all P>0.05), and the survival rate of 12 h and 24 h was lower than that of 0 h group (both P<0.01). In 0.05% group, there was no significant difference at each incubation time within 12 h (all P>0.05), but the survival rate at 24 h was lower than that at 0 h (P<0.05). There was no significant difference in survival rate at different culture time in 0.1% concentration group (P=0.082). At the same culture time, there was no significant difference in survival rate between 0.01% group and control group at 0, 3 and 24 h (all P>0.05). Except for 3 h, the survival rate of 0.05% and 0.1% groups was lower than that of control group (all P<0.05), and the decrease was all less than 30%. Uterine fluid RNA-seq showed that there was no significance difference in exonic rate, the detected genes and transcripts of RNA between the experiment groups and the control group (all P>0.05)., Conclusions: The in vitro cytotoxic of PVA formaldehyde absorbent sponge on human endometrial epithelial cell meet the national standard of the cytotoxic of medical materials. Sampling the uterine fluid with this material does not affect the RNA-Seq results. PVA formaldehyde absorbent sponge is safe and feasible when appling to the noninvasive uterine fluid sampling and RNA sequencing.

Objective: To explore the heterogeneity and growth factor regulatory network of dermal fibroblasts (dFbs) in mouse full-thickness skin defect wounds based on single-cell RNA sequencing. Methods: The experimental research methods were adopted. The normal skin tissue from 5 healthy 8-week-old male C57BL/6 mice (the same mouse age, sex, and strain below) was harvested, and the wound tissue of another 5 mice with full-thickness skin defect on the back was harvested on post injury day (PID) 7. The cell suspension was obtained by digesting the tissue with collagenase D and DNase Ⅰ, sequencing library was constructed using 10x Genomics platform, and single-cell RNA sequencing was performed by Illumina Novaseq6000 sequencer. The gene expression matrices of cells in the two kinds of tissue were obtained by analysis of Seurat 3.0 program of software R4.1.1, and two-dimensional tSNE plots classified by cell group, cell source, and gene labeling of major cells in skin were used for visual display. According to the existing literature and the CellMarker database searching, the expression of marker genes in the gene expression matrices of cells in the two kinds of tissue was analyzed, and each cell group was numbered and defined. The gene expression matrices and cell clustering information were introduced into CellChat 1.1.3 program of software R4.1.1 to analyze the intercellular communication in the two kinds of tissue and the intercellular communication involving vascular endothelial growth factor (VEGF), platelet-derived growth factor (PDGF), epidermal growth factor (EGF), and fibroblast growth factor (FGF) signal pathways in the wound tissue, the relative contribution of each pair of FGF subtypes and FGF receptor (FGFR) subtypes (hereinafter referred to as FGF ligand receptor pairs) to FGF signal network in the two kinds of tissue, and the intercellular communication in the signal pathway of FGF ligand receptor pairs with the top 2 relative contributions in the two kinds of tissue. The normal skin tissue from one healthy mouse was harvested, and the wound tissue of one mouse with full-thickness skin defect on the back was harvested on PID 7. The multiple immunofluorescence staining was performed to detect the expression and distribution of FGF7 protein and its co-localized expression with dipeptidyl peptidase 4 (DPP4), stem cell antigen 1 (SCA1), smooth muscle actin (SMA), and PDGF receptor α (PDGFRα) protein. Results: Both the normal skin tissue of healthy mice and the wound tissue of full-thickness skin defected mice on PID 7 contained 25 cell groups, but the numbers of cells in each cell group between the two kinds of tissue were different. Genes PDGFRα , platelet endothelial cell adhesion molecule 1, lymphatic endothelial hyaluronic acid receptor 1, receptor protein tyrosine phosphatase C, keratin 10, and keratin 79 all had distinct distributions on two-dimensional tSNE plots, indicating specific cell groups respectively. The 25 cell groups were numbered by C0-C24 and divided into 9 dFb subgroups and 16 non-dFb groups. dFb subgroups included C0 as interstitial progenitor cells, C5 as adipose precursor cells, and C13 as contractile muscle cells related fibroblasts, etc. Non-dFb group included C3 as neutrophils, C8 as T cells, and C18 as erythrocytes, etc. Compared with that of the normal skin tissue of healthy mice, the intercellular communication in the wound tissue of full-thickness skin defected mice on PID 7 was more and denser, and the top 3 cell groups in intercellular communication intensity were dFb subgroups C0, C1, and C2, of which all communicated with other cell groups in the wound tissue. In the wound tissue of full-thickness skin defected mice on PID 7, VEGF signals were mainly sent by the dFb subgroup C0 and received by vascular related cell groups C19 and C21, PDGF signals were mainly sent by peripheral cells C14 and received by multiple dFb subgroups, EGF signals were mainly sent by keratinocyte subgroups C9 and C11 and received by the dFb subgroup C0, and the main sender and receiver of FGF signals were the dFb subgroup C6. In the relative contribution rank of FGF ligand receptor pairs to FGF signal network in the normal skin tissue of healthy mice and the wound tissue of full-thickness skin defected mice on PID 7, FGF7-FGFR1 was the top 1, and FGF7-FGFR2 or FGF10-FGFR1 was in the second place, respectively; compared with those in the normal skin tissue, there was more intercellular communication in FGF7-FGFR1 signal pathway, while the intercellular communication in FGF7-FGFR2 and FGF10-FGFR1 signal pathways decreased slightly or did not change significantly in the wound tissue; the intercellular communication in FGF7-FGFR1 signal pathway in the wound tissue was stronger than that in FGF7-FGFR2 or FGF10-FGFR1 signal pathway; in the two kinds of tissue, FGF7 signal was mainly sent by dFb subgroups C0, C1, and C2, and received by dFb subgroups C6 and C7. Compared with that in the normal skin tissue of healthy mouse, the expression of FGF7 protein was higher in the wound tissue of full-thickness skin defected mouse on PID 7; in the normal skin tissue, FGF7 protein was mainly expressed in the skin interstitium and also expressed in the white adipose tissue near the dermis layer; in the two kinds of tissue, FGF7 protein was co-localized with DPP4 and SCA1 proteins and expressed in the skin interstitium, co-localized with PDGFRα protein and expressed in dFbs, but was not co-localized with SMA protein, with more co-localized expression of FGF7 in the wound tissue than that in the normal skin tissue. Conclusions: In the process of wound healing of mouse full-thickness skin defect wound, dFbs are highly heterogeneous, act as potential major secretory or receiving cell populations of a variety of growth factors, and have a close and complex relationship with the growth factor signal pathways. FGF7-FGFR1 signal pathway is the main FGF signal pathway in the process of wound healing, which targets and regulates multiple dFb subgroups.

The human body is colonized by densely-populated and structurally complex communities of microorganisms. The microbiota interact not only with their host cells, but also with other microbiota. Dual RNA sequencing (Dual RNA-seq) can be used to conduct simultaneous analysis of the dynamic changes of gene expression of two (or more) interactive species, and to obtain thus, through the interaction model diagram, the inter-species regulatory relationship of genes of different species, and hence the interaction mechanism between species. We herein reviewed the application status and development prospects of Dual RNA-seq in the research of intestinal, respiratory, skin and oral microbes. Since the concept of Dual RNA-seq was first introduced, the technology has been applied to a range of infection models. Direct investigation into the dynamic interactions between species at the molecular level will contribute to the better understanding of the physiological changes of pathogens and hosts during the course of infection, and thus help reveal potential new targets or biomarkers. However, the Dual RNA-seq technology is still in its early stage of development, and there are some limitations in the experimental technology. For example, due to the dynamic nature of the interaction between species, there are urgent problems awaiting solutions, such as the optimal experimental conditions, the selection of sampling sites and how to achieve real-time observation. In addition, due to the large amount of bioinformatics data of Dual RNA-seq, further research is needed to explore for ways to process the interaction information quickly and flexibly., (Copyright© by Editorial Board of Journal of Sichuan University (Medical Sciences).)

Traditional bulk RNA sequencing assesses the average expression level of genes in tissues rather than the differences in cellular responses. Accordingly, it is hard to differentiate sensitive responding cells, leading to inaccurate identification of toxicity pathways. Single-cell RNA sequencing (scRNA-seq) isolated single cells from tissue and subjected them to cell subtypes-specific transcriptome analysis. This technique in toxicological studies realizes the heterogeneous cellular responses in the tissue microenvironment upon chemical exposure. Thus it helps to identify sensitive responding cells and key molecular events, providing a powerful tool and a new perspective for exploring the mechanisms of toxicity and the modes of action. This review summarizes the development, principle, method, application and limitations of scRNA-seq in mechanistic toxicological researches, and discusses the prospect of multi-directional applications.

Objective: To clarify the effect of high-frequency repetitive transcranial magnetic stimulation (rTMS) on cognitive function in cerebral ischemic rats, and to explore its underlying mechanism by RNA sequencing. Methods: Thirty male Sprague-Dawley (SD) rats underwent transient middle cerebral artery occlusion (tMCAO). According to the Bederson score, 10 rats with a score of 1-3 were excluded, and the remaining 20 rats were then randomly divided into the tMCAO group ( n =10) and the rTMS group ( n =10). Meanwhile, 10 rats with sham operation were assigned to the sham group ( n =10). Rats in the rTMS group received 20 Hz rTMS from day 7 to day 28 after surgery. From day 28 to day 33 after the operation, Morris water maze test was performed to detect the cognitive function of rats in each group. The cortical tissues around the infarcts from the rTMS tMCAO groups were taken for RNA sequencing analysis, with 3 rats in each group. Results: The escape latency of rats in the rTMS group[ (53±4)s] and the group [(51±5)s] were significantly shorter than that of the tMCAO group[ (58±4)s, P <0.05)]. The times that the rats crossed the original platform in 60 seconds in the rTMS group[2.5 (1.5-3.0)] and sham group[3.0 (1.5-3.0)] were more than that of the tMCAO group [1.0(0.5-1.5)] ( P <0.05). RNA sequencing detected 16 significantly differentially expressed genes, including 9 up-regulated genes and 7 down-regulated genes. GO analysis showed that the functions of up-regulated genes were mainly concentrated in the processes of chemical homeostasis and cell metal ion homeostasis. While the functions of down-regulated genes mainly enriched in the inflammatory response. Conclusion: Twenty Hz rTMS can improve the cognitive function of rats with cerebral infarction, and its underlying mechanism may be related to maintaining chemical and metal ion homeostasis and regulating the polarization of microglia to reduce neuroinflammation.

The emergence of single-cell sequencing technology enables people to observe cells with unprecedented precision. However, it is difficult to capture the information on all cells and genes in one single-cell RNA sequencing (scRNA-seq) experiment. Single-cell data of a single modality cannot explain cell state and system changes in detail. The integrative analysis of single-cell data aims to address these two types of problems. Integrating multiple scRNA-seq data can collect complete cell types and provide a powerful boost for the construction of cell atlases. Integrating single-cell multimodal data can be used to study the causal relationship and gene regulation mechanism across modalities. The development and application of data integration methods helps fully explore the richness and relevance of single-cell data and discover meaningful biological changes. Based on this, this article reviews the basic principles, methods and applications of multiple scRNA-seq data integration and single-cell multimodal data integration. Moreover, the advantages and disadvantages of existing methods are discussed. Finally, the future development is prospected.

Objective: To optimize the protocol of meniscus cell extraction to enhance the efficiency of cell suspension preparation and maintain a high cell viability for single-cell RNA sequencing., Methods: We compared the efficiency of the routine cell extraction methods (short-time digestion and long-time digestion) and the optimized protocol for obtaining meniscus cell suspensions by evaluating the cell number obtained and the cell viability. Single-cell RNA sequencing datasets were analyzed to evaluate the stability of the cell suspension prepared using the optimized protocol. The reliability of the optimized protocol was assessed by comparing the single-cell RNA sequencing dataset obtained by the optimized protocol with published single-cell RNA sequencing datasets of the meniscus., Results: The optimized protocol harvested a greater number of cells (over 1×10 5 ) than the routine protocols. The cell suspension prepared with the optimized protocol showed a cell viability higher than 80%, the highest among the 3 methods. Analysis of single-cell RNA sequencing datasets showed that the ratio of the mitochondrial genes was below 20% in over 80% of the cells. CD34 + cells, MCAM + cells and COL1A1 + cells were identified in the datasets. Comparison with the publish datasets showed that the optimized protocol was capable of harvesting COL3A1 + , COL1A1 + , MYLK + , BMP2 + , CD93 + and CDK1 + cells., Conclusion: Single-cell suspension prepared from the meniscus can be stably obtained using the optimized protocol for single-cell RNA sequencing using the 10× Genomics platform.

Objective: To investigate the changes in gene expression profiles of mouse scleral fibroblasts after myopia using single-cell RNA sequencing technology and explore the mechanism of dysfunction of the scleral fibroblasts in myopia., Methods: Normal healthy C57BL/6J mice were randomly divided into negative control group and myopia model group ( n =6), and in the latter group, form deprivation myopia was induced in the right eye using translucent goggles. Single cell capture was performed in the right eye to obtain the scleral fibroblasts for RNA sequencing. The differentially expressed genes (DEGs) were screened and GO and KEGG analyses were carried out for functionally significant enrichment analysis., Results: Comparison of the gene expression profiles identified a total of 169 DEGs between the myopia model group and the negative control group ( P >0.05), including 112 up-regulated and 57 down-regulated genes. GO function analysis showed that the DEGs were involved in leukocyte aggregation, differentiation and adhesion and other inflammation-related terms; ATP metabolism and binding, redox process, oxidative stress response, oxidative phosphorylation and other GO terms related to hypoxia; protein folding, protein transport, negative regulation of protein metabolism, endoplasmic reticulum (ER), ER cavity, ER stress and other biological processes related to protein and ER stress. KEGG analysis analysis showed that the significantly enriched pathways of the DEGs involved mainly the TCA cycle, oxidative phosphorylation, PPAR signaling, oxidative phosphorylation and other pathways related to hypoxia; MAPK signaling pathway related to inflammation; leukocyte transendothelial transport; and protein-related Parkinson's disease, Huntington's disease, protein digestion and absorption pathways., Conclusion: The dysfunction of the scleral fibroblasts occurs in myopia through complex mechanisms involving inflammation, hypoxia, protein regulation, and ER stress-related gene expression and pathway regulation.

Single-cell transcriptome sequencing(scRNA-seq) can be used to analyze the expression characteristics of the transcriptome at the level of individual cell, and discover the heterogeneity of gene expression in individual cell that is &quot;diluted&quot; or averaged in study of group organization. The scRNA-seq, with the characteristics of standardization, high-throughput, and high integration, can greatly simplify the experimental operation and significantly reduce the consumption of reagents. At the same time, a variety of cells are screened and the gene expression patterns are analyzed at the single-cell level to provide a more efficient detection technique and more rich and accurate information for drug research. In the field of traditional Chinese medicine(TCM), the scRNA-seq is still a new technology, but the individual and precision concepts embodied by scRNA-seq and the theory of TCM syndrome differentiation and treatment have reached the same effect between the micro and macro aspects. This study tried to broaden the thinking for the modernization of TCM by introducing the development of scRNA-seq technology and its application in modern drug research and discussing the application prospects of scRNA-seq in TCM research.

Objective: To identify differentiation related miRNA and evaluate roles of miRNA during ATRA induced myeloid differentiation., Methods: The small RNA sequencing was used to analyze differential expressed miRNAs in ATRA induced NB4 cells. Then the several up or down-regulated miRNA were selected as the research candidates. SgRNAs targeting the genome of each miRNA were designed and NB4 cells with inducible expression of Cas9 protein were generated. After transduced sgRNA into NB4/Cas9 cells, the mutation level by PCR and surveyor assay were evaluated. The cell differentiation level was investigated by surface CD11b expression via flow cytometry., Results: A total of 410 mature miRNAs which expressed in NB4 cells were detected out after treated by ATRA, 74 miRNAs were up-regulated and 55 were down-regulated miRNAs with DNA cleavage generated by CRISPR/Cas9 was assayed directly by PCR or surveyor assay, quantitative PCR showed that the expression of miRNA was downregulated, which evaluated that gene edition successfully inhibitied the expression of mature miRNA. MiR-223 knockout showed the myeloid differentation of NB4 significantly inhibitied, while miRNA-155 knockout showed the myeloid differentation of NB4 cells significantly increased., Conclusion: CRISPR/Cas9 is a powerful tool for gene editing and can lead to miRNA knockout. Knockouts of miR-223 and miR-155 have shown a differentiation-related phenotype, and the potential mechanism is the integrative regulation of target genes.

After the first single-cell RNA sequencing of human islet in 2016, a large number of analysis of RNA sequencing at single-cell level on mouse and human islets emerged, which has brought new progress to the research of islet cell biology. It has been proved that single-cell RNA sequencing can characterize rare types of endocrine cells, find the cellular heterogeneity of typical endocrine cells and new cell subtypes, analyze the differences of endocrine cells among species, and describe different states of various types of islet cells in developmental processes and metabolic diseases more accurately. Although the detection efficiency of single-cell RNA sequencing for low abundance transcripts is relatively low at present, with the improvement of technology and the progress of analytical methods, this technology will undoubtedly become a powerful tool to explore the cellular heterogeneity, development, and biological characteristics of metabolic diseases.

Gene panel and whole exome sequencing are now commonly used to detect Mendelian disease, but the current molecular diagnostic rate of DNA sequencing is only 35%-50%. In recent years, RNA sequencing emerges as a promising diagnostic method. It can detect new pathogenic mutations, and analyze allele-specific expression. This will be helpful to understand the relationship between disease genotype and phenotype, and can complement genome sequencing in order to expand the traditional genomic diagnostic methods of Mendelian disease. RNA sequencing is expected to become a routine tool for diagnosing Mendelian diseases. This article reviews the application of RNA sequencing in the clinical diagnosis of Mendelian disease.

Steroidal saponins, which are the characteristic and main active constituents of Polygonatum, exhibit a broad range of pharmacological functions, such as regulating blood sugar, preventing cardiovascular and cerebrovascular diseases and anti-tumor. In this study, we performed RNA sequencing(RNA-Seq) analysis for the flowers, leaves, roots, and rhizomes of Polygonatum cyrtonema using the BGISEQ-500 platform to understand the biosynthesis pathway of steroidal saponins and study their key enzyme genes. The assembly of transcripts for four tissues generated 129 989 unigenes, of which 88 958 were mapped to several public databases for functional annotation, 22 813 unigenes were assigned to 53 subcategories and 64 877 unigenes were annotated to 136 pathways in KEGG database. Furthermore, 502 unigenes involved in the biosynthesis pathway of steroidal saponins were identified, of which 97 unigenes encoding 12 key enzymes. Cycloartenol synthase, the first key enzyme in the pathway of phytosterol biosynthesis, showed conserved catalytic domain and substrate binding domain based on sequence analysis and homology modeling. Differentially expressed genes(DEGs) were identified in rhizomes as compared to other tissues(flowers, leaves or roots).The 2 437 unigenes annotated by KEGG showed rhizome-specific expression, of which 35 unigenes involved in the biosynthesis of steroidal saponins. Our results greatly extend the public transcriptome dataset of Polygonatum and provide valuable information for the identification of candidate genes involved in the biosynthesis of steroidal saponins and other important secondary metabolites.

Objective: To gain insight into the transcriptional landscape including mRNA, long non-coding RNA (lncRNA), and circular RNA (circRNA) of the differentiated lens. Methods: Experiment research. The total RNAs of the differentiated lenses were extracted and purified. Total RNAs of 16-week, 23-week, and 25-week differentiated lenses were then sequenced using Illumina HiSeq 2500, and analyzed using bioinformatics tools. The top expressed and differentially expressed mRNAs and lncRNAs were screened. The expressions of overlap genes among the 16-week, 23-week, and 25-week lenses were analyzed by Venn diagram. The expression tendency of lens-specific genes was obtained and verified with real-time polymerase chain reaction. Results: A total of 67 518 311 mapped reads were obtained from differentiated lenses at 16 weeks, 99 440 160 at 23 weeks, and 67 262 320 at 25 weeks. The gene overlap expression analysis showed 740 of the top 1 000 highly expressed mRNAs, 170 of the top 300 highly expressed lncRNAs, and 69 of the top 100 highly expressed circRNAs overlapping expressed in lenses at 16, 23, and 25 weeks, respectively. Lens specific gene expression analysis revealed that the expression of crystallin (CRY) AA, CRYGA, CRYGB, CRYGC, CRYGD, CRYGEP, and CRYGS was upregulated, while the expression of gap junction (GJ) A3 and GJA8 was downregulated with the differentiation of lenses. Conclusion: The lens transcriptome profile shows that more than half of the high expressed mRNA, lncRNA and circRNA at different differentiation stages are overlapping expressed, and all of them have high expression of lens specific protein genes, such as CRY, GJ etc. (Chin J Ophthalmol, 2020, 56: 356 - 363) .

Tumors are highly complex systems. Understanding the compositions and functions of the tumor immune microenvironment is a prerequisite for effective tumor immunotherapy. Single-cell RNA sequencing can detect the transcriptome of a cell at the resolution of single-cell level,describe its functional status,and thus deepen the understanding of the composition and function of different cell clusters in tumor immune microenvironment. This article reviews the application of single-cell RNA sequencing in research on tumor immune microenvironment.

RNA is a polymeric molecule which is involved in various biological processes including the coding, decoding, regulation, and expression of genes. Whole transcriptome sequencing is the dominant method for studying RNA functions which assists researchers to deepen the exploration and analysis of gene structure and function and to reveal intrinsic links between gene expression and life phenomena. To date, extensive research has been done in animal husbandry models including swine and poultry by using high-throughput RNA sequencing technology. These studies included transcript structure, gene expression level and non-coding region function. In this review, we briefly introduce the general processes and strategies of RNA-seq and small RNA-seq technologies, and summarize the various achievements of the application of whole transcriptomics in animal husbandry.

Forensic microorganism is one of the hotspots of forensic science research. Due to its conservatism and specificity, the 16S rRNA gene is found to be an ideal marker for forensic identification. With the rapid development of high throughput sequencing technology, the research on microorganisms has been gradually applied to many fields such as environment and health care. In the field of forensic science, the results of forensic microbiology research, represented by 16S rRNA gene sequencing, are also gradually applied to forensic practice, such as biological samples identification, individual identification, postmortem interval estimation, and regional inference, which not only provide clues for the investigation of cases but also complement and assist traditional methods. This paper describes the research methods and related sequencing technologies of 16S rRNA gene sequencing, summarizes its research progress, and discusses the application value and potential of 16S rRNA in forensic science., Competing Interests: The authors of this article and the planning committee members and staff have no relevant financial relationships with commercial interests to disclose., (Copyright© by the Editorial Department of Journal of Forensic Medicine.)

The NAC family is an important transcription factor which regulate plant growth and development， signal transduction， and stress response.In this study， the protein identification， subfamily classification， the determination of physical and chemical properties， protein structure， and expression pattern of NAC family were performed using bioinformatic methods based on the RNA-seq data of ginger. The results showed that a total of 72 NAC transcription factors were identified in 271.1 Mb total nucleotides， and they could be clustered into 13 subfamilies according to the phylogenetic tree.The physical and chemical properties， structure analysis revealed that the amino acid number and isoelectric point were different among 13 NAC subfamilies; the secondary structure of NACs transcription factors mainly consist of random coil， and the tertiary structure is similar.In addition，the expression patterns of genes under different soil moisture and Ralstonia solanacearum infection showed that 23 NACs were differentially expressed， which were mainly distributed in Ⅷ，Ⅶ， and ⅩⅤ subfamilies related to plant senescence， hormone metabolism and cell wall metabolism.The results provide some valuable information for the research and development of NAC transcription factors in ginger., Competing Interests: The authors of this article and the planning committee members and staff have no relevant financial relationships with commercial interests to disclose., (Copyright© by the Chinese Pharmaceutical Association.)

Objective: To produce a comprehensive transcript dataset of Oncomelania hupensis before and after Schistosoma japonicum infection, so as to provide experimental data for perfecting genetic structural information and excavating related molecular markers of O. hupensis infected by S. japonicum ., Methods: O. hupensis snails were divided into the following 3 groups: one week after S. japonicum miracidium infection, 4 weeks after S. japonicum miracidium infection, and normal condition. Million high-quality reads were obtained from the normalized cDNA of the pooled samples, which were assembled into transcripts., Results: A total of 63 686 unigenes were identified and were classified into 4 main categories, including general functional prediction (15.36%), signal transduction mechanism (11.75%), posttranslational modification (8.89%), and functional unknown (12.20%)., Conclusions: The transcriptome information of O. hupensis snail after the invasion of S. japonicum shows that several genes are significantly up-regulated or down regulated expression, and that the availability of transcriptome information might provide a strong foundation for further understanding the schistosome-snail interaction at the molecular level.

Objective: To identify non-coding RNAs in Haloferax mediterranei through high-throughput RNA sequencing, bioinformatics analysis and molecular techniques., Methods: After H. mediterranei cells under log phase of growth were treated with different salt concentrations for 30 minutes, total RNA was extracted for the following strand-specific RNA sequencing and differential RNA sequencing. These RNA-seq data were used to identify the genome-wide ncRNAs and to predict the 5' and 3'-ends of the transcripts by bioinformatics analysis. A few selected ncRNAs were further confirmed by Northern blotting and Circularized RNA reverse transcription-PCR analysis., Results: We identified 105 highly credible ncRNAs. Expression of four ncRNAs showed difference in different salt concentrations. We confirmed the expression, length of transcripts, transcription start and termination sites of incRNA1436 and incRNA1903 by Northern blotting and CR-RT-PCR., Conclusion: We identified the ncRNAs of H. mediterranei in a genome-wide scale, including identification of a few ncRNAs involved in the responses of H. mediterranei to different salt concentrations. Our results have provided fundamental data and novel insights for future study of the function of ncRNA in haloarchaea.

Objective: Food phytochemicals as biofilm inhibitor of pathogens have been highlighted. Our study aimed to investigate the effects of resveratrol on biofilm formation of an aquatic pathogen Vibrio parahaemolyticus, and to elucidate the important regulatory genes., Methods: In the subinhibitory concentrations, the inhibition of resveratrol aganist biofilm and exopolysaccharides of V. parahaemolyticus was detected, and the differentially expressed genes were analyzed based on RNA-Seq. Four genes involved in biofilm formation was validated by qRT-PCR., Results: The minimum inhibitory concentration of resveratrol against V. parahaemolyticus was 20 μg/mL. Resveratrol at the subinhibitory concentration of 5 μg/mL and 10 μg/mL significantly decreased the biofilm development and exopolysaccharides production in V. parahaemolyticus (P < 0.05). Scanning electron microscopy micrographs showed a significant reduction of adherence and extracellular polymeric substances. RNA-Seq analysis revealed that 22.6% up-regulated and 77.4% down-regulated gene (P < 0.05) after treatment by 10 μg/mL resveratrol among 106 differential gene expressions. These differential genes in V. parahaemolyticu focused on 7 metabolic pathways, and 14 genes involved in biofilm formation were down-regulated by resveratrol, such as outer membrane protein (W, YedS, OmpK), quorum sensing (LuxS), flagellin (FlaA), fimbrial assembly protein (PilQ), hemolysin secreted protein. qRT-PCR confirmed that the expressions of luxS, trh, tlh and flaA, was significantly repressed in the presence of resveratrol, which was consistent with transcriptomics data., Conclusion: Inhibitory activity of resveratrol on biofilm formation was assicated with multiple genes and diverse cellular processes in V. parahaemolyticus. These findings suggest that resveratrol would disturb various metabolic pathways, particularly quorum sensing system, adhesion process and membrane proteins secretion, resulting in inhibition of attachment and biofilm development. The present work provided valuable information to explore molecular mechanism of resveratrol as an novel anti-biofilm compound.

Miscanthus sinensis Anderss is a perennial C4-grass. It is a promising bioenergy plant, which has been proposed as general feedstock for biomass and lignocellulosic biofuel production. In this study, the flower and leaf buds transcriptomes of Miscanthus sinensis Anderss were sequenced by the platform of Illumina HiSeq 2000. In total 98 326 Unigenes were generated by de novo assembly with an average length of 822 bp and N50 of 1 023 bp. Based on the NR, NT, Swiss-Prot, KEGG, GO and COG databases (Evalue < le-5), 74 134 (75.40%) Unigenes were annotated. A total of 45 507 Unigenes were mapped into different GO terms. In KEGG pathways identification, 36 710 sequences were assigned to 128 KEGG pathways. Sorghum bicolor (37 731, 60.86%), Zea mays (16 258, 26.22%), and Oryza sativa (3 065, 4.94%) showed high similarity to Miscanthus sinensis Anderss. And 24 photosynthesis-related enzyme genes were identified. The result provides a foundation for further characterizing the functional genes in Miscanthus sinensis Anderss.

The rhizome of Panax japonicus var. major have been used as the natural medicinal agent by Chinese traditional doctors for more than thousand years. Most of the therapeutic effects of P. japonicus var. major had been reported due to the presence of tetracyclic or pentacyclic triterpene saponins. In this study, Illumina pair-end RNA-sequencing and de novo splicing were done in order to understand the pathway of triterpenoid saponins in this species. The valid reads data of 15. 6 Gb were obtained. The 62 240 unigenes were finally obtained by de novo splicing. After annotation, we discovered 19 unigenes involved in ginsenoside backbone biosynthesis. Additionally, 69 unigenes and 18 unigenes were predicted to have potential function of cytochrome P450 and UDP-glycosyltransferase based on the annotation results, which may encode enzymes responsible for ginsenoside backbone modification. This study provides global expressed datas for P. japonicus var. major, which will contribute significantly to further genome-wide research and analysis for this species.

The transcriptome represents the whole complement of RNA transcripts in cells or tissues and reflects the expressed genes at various life stages, tissue types, physiological states, and environmental conditions. Transcriptomics study concerning medicinal plants has become the most active area in medicinal plant genome research. Transcriptome analysis provides a comprehensive understanding of gene expression and its regulation. The study of its transcriptome has great significance in solving the questions of genetic evolution, genetic breeding, ecology and so on. Here we report the application status of transcriptomics in medicinal plants based on emergence, development and methodology of transcriptomics.

RNA-Sequencing (RNA-Seq) is a newly-developed method in transcriptome research, it can afford more accurate transcription information and be more quickly by using Next-generation Sequencing (NGS) technology. RNA-Seq has been widely used in various biological fields. Genuine traditional Chinese medicines (TCM), with good quality and therapeutic effect, were always praised highly and used by famous physicians. The geo-herbalism formation of TCM is based on the product of the gene expression at specific space and time. So it has been a research hotspot to analyze the mechanism of biosynthesis through RNA-Seq in the study on the secondary metabolism of medicinal plant. This article mainly illustrates the RNA-Seq and its advantages, it also discusses the potential application in genuine TCM, and it can provide useful information for other researchers.

Cell heterogeneity is a general feature of biological tissues. Standard transcriptome analysis approaches require tens of thousands of cells to provide an average view of gene expression and ignore the information of gene expression heterogeneity. The single-cell RNA-Seq technologies profile gene expression at the single-cell level and serve as powerful tools to identify distinct phenotypic cell types within a heterogeneous population. The single-cell RNA-Seq technologies have been developed rapidly in recent years. The methodological progress includes a variety of cDNA amplification methods, the quantitative analysis of the sensitivity and noise of the technologies, and the development of the low-coverage high-throughput single-cell RNA-Seq and the in situ RNA-Seq technologies. Furthermore, the scope of application is extended from early embryonic development to tissue and organ development, immunology and oncology. In this review, we discuss recent progress in methodology and applications of the single-cell RNA-Seq technologies.

Objective: To monitor and discover medically important mosquito-borne viruses circulating in Xinjiang, China., Methods: Mosquitoes were collected from Abbey Lake wetland in Bortala, in Northern Xinjiang. Viral isolates were obtained through inoculating and serial passaging into susceptible mammalian host cells (BHK-21), identified by cytopathogenic effect (CPE) observation and plague forming assay. Genetic identification of viral isolates was conducted by RT-PCR, sequencing and phylogenetic analysis., Results: A virus strain which causing CPE on BHK-21 cells, was isolated from the predominant Culex modestus (36.6%) and tentatively designated as Abbey Lake virus. Information on molecular identification revealed that Abbey Lake virus belonged to Orthobunyavirus genus within Bunyaviridae. Partial sequences (651 bp and 980 bp) of viral genomic S and M segment showed that Abbey Lake virus was phylogenetically related to Germiston virus that uniquely found in South Africa with 90.6% nucleotides and 95.0% amino acids similarities in S segment. However, viral M segment displayed much variability with 78.6% nucleotides and 86.1% amino acid similarities, suggesting a new member of Orthobunyavirus genus was discovered in the area., Conclusion: In this study, Abbey Lake virus was isolated and characterized indicating its potential circulation nature of this newly-emerged mosquito-borne virus.

Ammopiptanthus mongolicus, a woody plant growing in the desert, plays a vital role in vegetation maintaining and restoration in the arid region in northwest China. The plant exhibits an extremely high tolerance to abiotic stress such as drought and freezing stresses, and it has been used as an ideal model for abiotic stress tolerance research in trees. MicroRNA (miRNA) is a class of approximately 21nt endogenous non-protein-coding small RNA, which plays an important role in plant growth, development and responses to environmental stresses. By now, a large number of miRNAs have been reported in many plant species, but no studies describing A.mongolicus miRNA were published. In the present study, the types, expression levels, and putative target genes of conserved miRNAs in seedlings of A. mongolicus were analyzed using small RNA deep sequencing technology and bioinformatics methods. Nineteen conserved miRNAs, which belong to 10 miRNA families, were identified, with abundance ranging from 55 to 1920269 reads. Target prediction analysis determined the target genes of 14 conserved miRNAs. The functional classification analysis indicated that the conserved miRNAs participate in the development and environmental response by regulating the biological processes including the transcription regulation, hormone signal transduction, metabolisms and stress resistance.

Objective: To learn the prevalence of infection of human and animals severe fever with thrombocytopenia syndrome bunyavirus (SFTSV) in Yantai, Shandong province, and to analyze the pathogenic features of SFTSV as well as its relationship between human and animal hosts., Methods: From April to November in 2011, 3 576 serum samples were collected from domesticated animals, including sheep, cattle, pigs, dogs, chickens, in Laizhou and Penglai areas where fever with thrombocytopenia syndrome frequently occurred among local residents. Total SFTSV antibodies and virus-specific nucleic acids of the serum were tested by ELISA and Real time RT-PCR, respectively. SFTSV infection on each animal was observed in different months. 2 590 human serum samples were also collected in Laizhou and Penglai areas, with IgG antibodies tested by ELISA. Virus was isolated with Vero cells from the serum which SFTSV viral nucleic acids were positive. S fragments were amplified by RT-PCR and sequenced, with homology analysis conducted on these sequences., Results: The overall positive rate of serum samples from animals on the total SFTSV antibodies was 40.24% (1 439/3 576) while the positive rate for specific nucleic acids was 4.56% (163/3 576). The positive rates for SFTSV antibodies were 62.78%, 52.97%, 45.56%, 28.73%, 1.45% and the positive rates for specific nucleic acids were 5.72%, 4.63%, 3.02%, 5.25% and 3.73%, in sheep, cattle, chickens, dogs, pigs, respectively. The antigens/antibodies for SFTSV in animals changed seasonally. The overall positive rate for SFTSV IgG antibody from 2 590 human samples was 5.41%. Thirteen virus strains were isolated from these serum samples (10 strains from human and 3 strains from animals). The nucleotide homology of 13S fragments' sequences ranged from 95.23% to 100.00% and the nucleotide homology with the isolates from other provinces were between 94.72% and 99.13%. The homology was considered to be high., Conclusion: High prevalence of SFTSV infections occurred both in human and domestic animals in Yantai city. The nucleotide sequences of SFTSV were highly homologous among human and domestic animals. The findings suggested that domesticated animals might serve as SFTSV proliferation and the hosts for transmission thus should be attached great importance.

In this study, the roles of hMMS2 (human methyl methanesulfonate sensitive mutant 2) gene encoding the human ubiquitin-conjugating enzyme E2 variant 2 in the drug resistance in human colon carcinoma were investigated by using a well-differentiated human colorectal carcinoma L-OHP-resistant cell line, THC8307/L-OHP. THC8307/L-OHP cells were transfected via liposome along with plasmid pcDNA6.2-GW/EmGFP-miR-MMS2 expressing both miRNA against hMMS2 and GFP, followed by real-time fluorescent quantitative PCR and immunofluorescence to select stable transfectants with significantly reduced hMMS2 expression. Compared with untransfected or pcDNA6.2-GW/EmGFP vector-transfected cells, the hMMS2-depleted cells displayed significantly (P<0.05) reduced half inhibition concentration(IC50) resistance index (RI) and colony-forming efficiency (CFE) upon treatment with oxaliplatin (L-OHP), while its relative reverse efficiency(RRE) was significantly higher (P<0.05) than the control cells, indicating compromised ability of cell proliferation. Indeed, Rho-damine 123 staining and flow cytometry analyses revealed an increased rate of apoptosis in hMMS2-depleted cells while no difference in cell proliferation or apoptosis was observed between the two control cell lines. The above observations collec-tively indicate that suppression of hMMS2 reverses L-OHP tolerance in differentiated human colorectal carcinoma cells by promoting apoptosis.

The application of next generation sequencing (NGS) technique has a great impact on epigenetic studies. Coupled with NGS, a number of sequencing-based methodologies have been developed and applied in epigenetic studies, such as Whole Genome Bisulfite Sequencing (WGBS), Reduced Representation Bisulfite Sequencing (RRBS), Methylated DNA Immunoprecipitation Sequencing (MeDIP-seq), Chromatin Immunoprecipitation-Sequencing (ChIP-seq), TAB-seq (Tet-assisted Bisulfite Sequencing), Chromosome Conformation Capture Sequencing (3C-seq) and various of 3C-seq de-rivatives, DNase1-seq/MNase-seq/FAIRE-seqand RNA Sequencing (RNA-seq). These new techniques were used to iden-tify DNA methylation patterns and a broad range of protein/nucleic acid interactions, and to analyze chromatin conforma-tion.With these new technologies, researchers have gained a broader view and better tools to investigate the distributions and dynamic changes of epigenetic markers affected by both internal and external factors. The principles and characteristics of major applications of NGS technologies on epigenetics were summarized; and the recent advances and the future direc-tions in NGS-based epigenetic studies were further discussed.