ObjectiveTo reveal the potential mechanism of cisplatin resistance in non-small cell lung cancer A549 cells by comparing the expression profiles of wild-type A549 cells and cisplatin-resistant A549 cells (A549/DPP) through whole transcriptome sequencing analysis.MethodsThe cisplatin resistant A549 (A549/DDP) cell line was first established. Then, the whole-transcriptome analysis was conducted both on A549 and A549/DDP cells. Next, the differentially expressed RNAs of lncRNA-seq, circRNA-seq, and miRNA-seq data were identified, respectively, followed by functional enrichment analysis. Finally, a comprehensive analysis based on the whole transcriptome data was performed and the construction of the ceRNA network was carried out.ResultsA total of 4 517 lncRNA, 123 circRNA, and 145 miRNA were differentially expressed in A549/DDP cells compared with the A549 cell line. These different RNAs were significantly enriched in cancer-related pathways. The ceRNA network contained 12 miRNAs, 4 circRNAs, 23 lncRNAs, and 9 mRNA nodes, of which hsa-miR-125a-5p and hsa-miR-125b-5p were important miRNAs based on the topological analysis.ConclusionTumor necrosis factor signaling pathway and p53 signaling pathway are involved in A549/DPP resistance. Hsa-miR-125a-5p and hsa-miR-125b-5p may be potential targets for reversing cisplatin resistance.
RNA can be labeled by more than 170 chemical modifications after transcription, and these chemical modifications are collectively referred to as RNA modifications. It opened a new chapter of epigenetic research and became a major research hotspot in recent years. RNA modification regulates the expression of genes from the transcriptome level by regulating the fate of RNA, thus participating in many biological processes and disease occurrence and development. With the deepening of research, the diversity and complexity of RNA modification, as well as its physiological significance and potential as a therapeutic target, can not be ignored.
ObjectiveTo study the abnormal biological pathways of intrahepatic cholangiocarcinoma (ICC) from the transcriptomics level and identify genes associated with the prognosis of ICC.MethodsThe differentially expressed genes were screened by t test and fold change method, then KEGG functional enrichment analysis was performed on related genes. The STRING database was applied to construct protein interaction network and find the hub nodes of the network by calculating the degree, betweenness, and closeness of each node. Kaplan-Meier survival analysis was performed using log-rank test to identify prognostic genes related to ICC.ResultsAll of 1 134 differentially expressed genes were overlapped in 3 datasets, which were mainly involved in 15 pathways, including DNA replication, cell cycle, drug metabolism, RNA transport, etc. signaling pathways and amino acid synthesis. According to protein interaction network analysis, TAF1, GRB2, E2F4, HNF4A, MYC, and TP53 genes were hub nodes. As GRB2 and TP53 genes were also the death related genes of ICC, it was found that patients with lower GRB2 gene expression had a better overall survival than those with higher GRB2 gene expression (P=0.040 9), while patients with lower TP53 had a worse overall survival than those with higher TP53 gene expression (P=0.027 3), which were also verified in the TCGA database.ConclusionsThe abnormal cell metabolism is notably related to the tumorigenesis of ICC. TAF1, GRB2, E2F4, HNF4A, MYC, and TP53 are the key genes in the carcinogenesis and progression of ICC. Expressions of GRB2 and TP53 genes are associated with the prognosis of ICC.
Objective To study the differential expression profiling of the transcripts modified by m5C methylation in a rat model of N-methyl-D-aspartate (NMDA)-induced retinal excitotoxicity. MethodsA total of 65 Sprague Dawley male rats aged 7-8 weeks were randomly divided into two groups: normal control group and NMDA group. The right eye (model eye) of rats in the NMDA group were injected with 50.0 mmol/L of NMDA 3 μl in the vitreous cavity, while in the normal control group, equal volume of normal saline was injected into the vitreous cavity. After 1 week of the injection, the optic nerve conduction function of rats was detected by visual evoked potential. The whole structure of rat retina was observed by hematoxylin-eosin staining, and the thickness of each retinal layer and the number of retinal ganglion cell layer were detected. The number of β3 tubulin immunofluorescence positive cells was detected by immunofluorescence staining on retinal stretched preparation. Total RNA was extracted from the retinas of normal control group and NMDA group, and high-throughput m5C modified RNA was sequenced, and bioinformatics analysis was performed. The relative expression levels of SLFN3, PLXNB3, CD36 and HIC2 mRNA in retina were detected by real-time quantitative polymerase chain reaction. The comparison between the two groups was performed using an unpaired t test. ResultsThe P1 latency of control group and NMDA group were (117.86±6.48) and (148.46±3.78) ms, and the amplitudes were (42.57±2.41) and (8.68±0.63) μV, respectively. Compared with the normal control group, the latency period was prolonged and the amplitude was significantly decreased in the NMDA group, with statistical significance (P＜0.001). In normal control group, retinal ganglion cells (RGC) were uniformly arranged with large round nuclei. In NMDA group, the volume of retinal RGC was atrophied and the number of RGC was reduced. The total retinal thickness in the control group and NMDA group was (207.51±12.76) μm and (187.51±12.54) μm, respectively. The number of β3 tubulin positive cells was 79.86±6.56 and 29.36±2.16, respectively. Compared with normal control group, the total retinal thickness and the number of β3 tubulin positive cells in NMDA group were decreased, with statistical significance (P＜0.001). Compared with the control group, 576 differentially expressed m5C mRNA were screened in the NMDA group, among which 230 up-regulated and 346 down-regulated genes were detected, respectively. The results of biological information analysis showed that compared with the control group, the upregulated m5C mRNA in the NMDA group was mainly involved in biological processes such as perception and cell-cell adhesion, and was mainly concentrated in the cytokine-cytokine receptor interaction and neural active ligand-receptor interaction pathway. The biological processes in which down-regulated m5C mRNA was mainly involved in biological processes such as G-protein-coupled receptor signaling pathway and cell communication, which were mainly concentrated in primary immune deficiency pathway and neural active ligand-receptor interaction pathway. Real-time quantitative polymerase chain reaction detection results showed that compared with the normal control group, the relative expression levels of SLFN3 and PLXNB3 mRNA in the retina of rats in NMDA group were significantly increased, while the relative expression levels of CD36 and HIC2 mRNA were significantly decreased, with statistical significance (P＜0.05). ConclusionIn NMDA induced retinal excitatory toxicity rat models, m5C modified retinal transcriptome showed abnormal expression.
Objective To explore the role and clinical significance of cell-cycle dependent kinase 1 (CDK1) and its upstream and downstream molecules in the development of malignant peripheral nerve sheath tumor (MPNST) through the analysis of clinical tissue samples. Methods A total of 56 tumor samples from MPNST patients (“Tianjin” dataset) who underwent surgical resection, confirmed by histology and pathology between September 2011 and March 2020, along with 17 normal tissue samples, were selected as the research subjects. MPNST-related hub genes were identified through transcriptome sequencing, bioinformatics analysis, immunohistochemistry staining, and survival analysis, and their expression levels and prognostic associations were analyzed. Results Transcriptome sequencing and bioinformatics analysis revealed that upregulated genes in MPNST were predominantly enriched in cell cycle-related pathways, with CDK1 occupying a central position among all differentially expressed genes. Further differential analysis demonstrated that CDK1 mRNA expression in sarcoma tissues was significantly higher than in normal tissues [based on searching the cancer genome atlas (TCGA) dataset, P<0.05]. In MPNST tissues, CDK1 mRNA expression was not only significantly higher than in normal tissues (based on Tianjin, GSE141438 datasets, P<0.05), but also significantly higher than in neurofibromatosis (NF) and plexiform neurofibromas (PNF) (based on GSE66743 and GSE145064 datasets, P<0.05). Immunohistochemical staining results indicated that the expression rate of CDK1 protein in MPNST tissues was 40.31%. Survival analysis results demonstrated that CDK1 expression was associated with poor prognosis. The survival time of MPNST patients with high CDK1 mRNA expression was significantly lower than that of the low expression group (P<0.05), and the overall survival trend of patients with positive CDK1 protein expression was worse than that of patients with negative CDK1 expression. Additionally, differential analysis of CDK family genes (CDK1-8) revealed that only CDK1 was significantly upregulated in MPNST, NF, and PNF. Conclusion Increased expression of CDK1 is associated with poor prognosis in MPNST patients. Compared to other CDK family members, CDK1 exhibits a unique expression pattern, suggesting its potential as a therapeutic target for MPNST.
Single-cell transcriptome sequencing (scRNA-seq) can resolve the expression characteristics of cells in tissues with single-cell precision, enabling researchers to quantify cellular heterogeneity within populations with higher resolution, revealing potentially heterogeneous cell populations and the dynamics of complex tissues. However, the presence of a large number of technical zeros in scRNA-seq data will have an impact on downstream analysis of cell clustering, differential genes, cell annotation, and pseudotime, hindering the discovery of meaningful biological signals. The main idea to solve this problem is to make use of the potential correlation between cells and genes, and to impute the technical zeros through the observed data. Based on this, this paper reviewed the basic methods of imputing technical zeros in the scRNA-seq data and discussed the advantages and disadvantages of the existing methods. Finally, recommendations and perspectives on the use and development of the method were provided.
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.
目的 研究糖尿病、阳痿、慢性肾炎患者中筛选出的典型肾阳虚证病例的转录组学特征，揭示出肾阳虚证“同证异治”的生物学基础。 方法 分别对9例肾阳虚证患者进行Agilent人444k表达谱芯片实验，对差异表达基因进行基因本休论（GO）、Pathway分析。 结果 找出332条共同差异表达基因，其中有注释的基因为181条。通过GO分析发现肾阳虚证在免疫系统、氨基酸分解和合成、脂类代谢、生殖、能量代谢及肿瘤相关的基因有密切联系，通过Pathway分析发现与肾阳虚证相关的信号通路有39个。 结论 肾阳虚证可能导致免疫系统低下，氨基酸分解和合成、脂类代谢、生殖、能量代谢功能降低，以及与肿瘤形成相关。
Objective To study the effects of connective tissue growth factor (CTGF) on retinal Müller cells based on transcriptome analysis of RNA-seq technology.MethodsRetinal Müller cells were divided into the control group and the CTGF treatment group which was continuously cultured with 10 ng/ml of CTGF for 24 h. The influence of CTGF on the migration of Müller cells were tested by scratching experiments. The RNA transcriptome analysis was applied to complete transcriptome sequencing, differentially expressed genes and functional enrichment analysis of the two groups of cells. HiSeq sequencing technology was used to sequence the whole transcriptome of the two groups of cells to obtain biological big data, and analyze the differentially expressed miRNAs on this basis. The functions and signal pathways of differential miRNAs were analyzed through gene annotation (GO) functional significance enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway significant enrichment analysis. Based on transcriptome data, genes with differential expression multiples in the top ten between the two groups were screened out, and the expression of bone morphogenetic protein 4 (BMP4) gene was verified by real time fluorescence quantification PCR (qRT-PCR), immunofluorescence and Western blot.ResultsAfter CTGF stimulation of Müller cells, cell viability and mobility which compared with the control group were significantly increased, with statistically significant differences (t=3.453, P＜0.05). The differential gene expression profile of CTGF induced Müller cells was obtained by RNA transcriptome analysis. Comparing the sequencing results of the two groups, it was found that 325 differentially expressed genes included 152 up-regulated genes and 173 down-regulated genes. The results of GO functional significance enrichment analysis showed that the functions of differential miRNA were mainly divided into three categories: biological processes, cellular components, and molecular functions. These differentially expressed genes were involved in signaling between nervous systems, adhesion between cells, and the interaction between cytokines and their receptors. These differentially expressed genes were involved in different metabolic pathways and biological processes such as tissue inflammation and fibrosis. BMP4 gene was seected for verification through immunofluorescence, qRT-PCR and western blot. The results showed that the expression of BMP4 was significantly higher than that in the control group, and the difference was statistically significant (t=39.490, 10.110, 5.470; P=0.004, 0.001, 0.006).ConclusionCTGF promotes cell proliferation and migration by up-regulating the expression of BMP4 in Müller cells, leading to tissue fibrosis and inducing inflammation.
Non-small cell lung cancer (NSCLC) is one of the most common types of cancer in the world and is an important cause for cancer death. Although the application of immunotherapy in recent years has greatly improved the prognosis of NSCLC, there are still huge challenges in the treatment of NSCLC. The immune microenvironment plays an important role in the process of NSCLC development, infiltration and metastasis, and they can interact and influence each other, forming a vicious circle. Notably, single-cell RNA sequencing enables high-resolution analysis of individual cells and is of great value in revealing cell types, cell evolution trajectories, molecular mechanisms of cell differentiation, and intercellular regulation within the immune microenvironment. Single-cell RNA sequencing is expected to uncover more promising immunotherapies. This article reviews the important researches and latest achievements of single-cell RNA sequencing in the immune microenvironment of NSCLC, and aims to explore the significance of applying single-cell RNA sequencing to analyze the immune microenvironment of NSCLC.