ObjectiveTo bioinformatically analyze the gene chip data of chondrocytes from osteoarthritis patients from the Gene Expression Omnibus (GEO) database, and explore the molecular mechanisms of osteoarthritis.MethodsWe searched the GEO database (up to April 23rd, 2021) for data of chondrocytes and gene expression profiling in human knee osteoarthritis via the key words of “osteoarthritis OR cartilage OR chondrocyte*”. Then, we selected the samples by our inclusion criteria. The data were normalized before analysis. After differentially expressed genes were identified, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Search Tool for the Retrival of Interacting Genes/Proteinsm, R language, Perl language, Cytoscape software, and DAVID database were used to perform differentially expressed gene analysis, functional annotation, and enrichment analysis.ResultsThe differentially expressed genes were mostly enriched in cell components and some extracellular regions, which participated in cell division, mitosis, cell proliferation and inflammatory response mainly via the regulation of protein kinase activity. The differentially expressed genes were mainly involved in the cell proliferation signaling pathway, mitogen-activated protein kinase signaling pathway, oocyte meiosis, cell cycle and so on.ConclusionsMultiple signaling pathways are involved in the changes of chondrocytes in human knee osteoarthritis, mainly about cell cycle and protein metabolism genes/pathways. Inflammatory factors and cytokines may be the most important links in the pathogenesis of osteoarthritis.
Objective To explore the pathogenesis of acute respiratory disease syndrome (ARDS) by bioinformatics analysis of neutrophil gene expression profile in order to find new therapeutic targets. Methods The gene expression chips include ARDS patients and healthy volunteers were screened from the Gene Expression Omnibus (GEO) database. The differentially expressed genes were carried out through GEO2R, OmicsBean, STRING, and Cytoscape, then enrichment analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathways was conducted to investigate the biological processes involved in ARDS via DAVID website. Results Bioinformatics analysis showed 86 differential genes achieved through the GEO2R website. Eighty-one genes were included in the STRING website for protein interaction analysis. The results of the interaction were further analyzed by Cytoscape software to obtain 11 hub genes: AHSP, ALAS2, CD177, CLEC4D, EPB42, GPR84, HBD, HVCN1, KLF1, SLC4A1, and STOM. GO analysis showed that the differential gene was enriched in the cellular component, especially the integrity of the plasma membrane. KEGG analysis showed that multiple pathways especially the cytokine receptor pathway involved in the pathogenesis of ARDS. Conclusions A variety of genes and pathways have been involved in the pathogenesis of ARDS. Eleven hub genes are screened, which may be involved in the pathogenesis of ARDS and can be used in subsequent studies.
Objective To validate the different expressions of human fxyd6 gene between normal bile duct tissues and malignant tumor tissues, and to observe the subcellular localization of human fxyd6 gene in human cholangiocarcinoma cells. MethodsThe different expressions between normal bile duct tissues and malignant tumor tissues were identified by RT-PCR. In situ polymerase chain reaction （IS-RT-PCR） was applied to detect the subcellular localization of fxyd6 gene in paraffin sections of human cholangiocarcinoma cells. Image analysis software was used to semiquantitatively determine the difference between normal and malignant tissues. ResultsHuman fxyd6 gene was highly expressed in cholangiocarcinoma tissues and lowly expressed in normal ones. There was a significant difference between the expressions of carcinoma cells and normal cells （P＜0.05）. IS-RT-PCR showed that fxyd6 gene localized in the kytoplasma of epithelial cells of human cholangiocarcinoma. ConclusionHuman fxyd6 gene may act as an essential component of the malignant transformation process in human cholangiocarcinoma.
Lung cancer is one of the malignant tumors with the greatest threat to human health, and studies have shown that some genes play an important regulatory role in the occurrence and development of lung cancer. In this paper, a LightGBM ensemble learning method is proposed to construct a prognostic model based on immune relate gene (IRG) profile data and clinical data to predict the prognostic survival rate of lung adenocarcinoma patients. First, this method used the Limma package for differential gene expression, used CoxPH regression analysis to screen the IRG to prognosis, and then used XGBoost algorithm to score the importance of the IRG features. Finally, the LASSO regression analysis was used to select IRG that could be used to construct a prognostic model, and a total of 17 IRG features were obtained that could be used to construct model. LightGBM was trained according to the IRG screened. The K-means algorithm was used to divide the patients into three groups, and the area under curve (AUC) of receiver operating characteristic (ROC) of the model output showed that the accuracy of the model in predicting the survival rates of the three groups of patients was 96%, 98% and 96%, respectively. The experimental results show that the model proposed in this paper can divide patients with lung adenocarcinoma into three groups [5-year survival rate higher than 65% (group 1), lower than 65% but higher than 30% (group 2) and lower than 30% (group 3)] and can accurately predict the 5-year survival rate of lung adenocarcinoma patients.
Objective To investigate specific changes of T cell repertoire in convalescent patients infected by influenza A (H7N9) virus. Methods Peripheral blood samples from 8 convalescent patients infected by H7N9 virus and 10 healthy donors were collected. After extracting whole DNA from these samples, arm-PCR were performed and the products were submitted to Illumina HiSeq2000 platform to produce deep sequencing data of the nucleotide sequences of complementary determining region 3 of T cell receptor β chain (TRB). Differences were compared in TRB diversity and V-D-J gene usage and similarities of sequences between the patients and the healthy donors. Results Frequency of V-D-J gene usage was different between the H7N9 patient group and the healthy group, such as TRBV30, TRBV27, and TRBV18 (Student's t test, P < 05). Main component analysis showed V-J pairing pattern was significantly different between two groups, which may have potential in identifying patients from healthy people. A considerable number of shared CDR3s were found in patient-patient pairs and normal-normal pairs, while seldom were found in patient-normal pairs. The similarity between patients was also confirmed by overlap distance analysis. Indexes for assessing diversity of immune repertoires, Shannon-Weiner index and Simpson index, were both lower in the patients (Student's t test, P < 05), suggesting that the immune system of the patients had not recovered 6 months after H7N9 infection. Compared with the healthy donors, the number of hyper-expression clones increased in the patient group, and some of them showed similarity among patients. Conclusions TRB repertoires are less diverse in patients with increased hyper-expressed clones and identifiable V-J usage pattern, which is identifiable from normal population. These results suggest that there are H7N9-specific changes in TRB repertoires of H7N9 infected patients in convalescent phase, which have potential implication in diagnosis and therapeutic T cell development.
ObjectiveTo investigate the expression of Yes-associated protein (YAP) screened by bioinformatics in rats with myocardial-ischemia reperfusion injury and establish the base for further research. MethodsThe difference of gene spectrum of rats with myocardial-ischemia reperfusion injury was analyzed by bioinformatics technique. The related signaling pathways and key genes were screened by KOBAS2.0 and KEGG. Eighteen Sprague Dawley rats were randomly divided into three groups: normal group (n=6), sham operation group (n=6) and myocardial-ischemia reperfusion injury group (n=6). The expression of target gene was detected by immunochemistry, quantitive reverse transcription polymerase chain reaction and western blotting. ResultsA total of 345 differentially expressed genes were found by bioinformatics, among which 181 were up-regulated and 164 were down-regulated. The differential genes were mainly enriched in Wnt, HIPPO, MAPK, Jak-STAT and other signaling pathways. We focused on HIPPO pathway and found that the expression of YAP increased significantly in myocardial-ischemia reperfusion injury group, compared with the normal group and sham operation group (P＜0.05). ConclusionsThe expression of YAP of HIPPO signal pathway is increased in rats with myocardial-ischemia reperfusion injury.
ObjectiveA competing endogenous RNA (ceRNA) regulatory network associated with long non-coding RNA (lncRNA) specific for lung adenocarcinoma (LUAD) was constructed based on bioinformatics methods, and the functional mechanism of actinfilament-associated protein 1-antisense RNA1 (AFAP1-AS1) in LUAD was analyzed, in order to provide a new direction for the study of LUAD therapeutic targets. MethodsThe gene chip of LUAD was downloaded from the Gene Expression Omnibus (GEO), and lncRNA and mRNA with differential expression between LUAD and normal tissues were screened using GEO2R online software, and their target genes were predicted by online databases to construct ceRNA networks and perform enrichment analysis. In cell experiments, AFAP1-AS1 was genetically knocked down and siRNA was constructed and transfected into LUAD cells A549 by cell transfection. CCK8, transwell, scratch assay and flow cytometry were used to detect the ability of cells to proliferate, invade, migrate and apoptosis. ResultsA total of 6 differentially expressed lncRNA and 494 differentially expressed mRNA were identified in the microarray of LUAD. The ceRNA network involved a total of 6 lncRNA, 22 miRNA, and 55 mRNA. Enrichment analysis revealed that mRNA was associated with cancer-related pathways. In cell assays, knockdown of AFAP1-AS1 inhibited cell proliferation, invasion, and migration, and AFAP1-AS1 promoted apoptosis. ConclusionIn this study, we construct a lncRNA-mediated ceRNA network, which may help to further investigate the mechanism of action of LUAD. In addition, through cellular experiments, AFAP1-AS1 is found to have potential as a therapeutic target for LUAD.
現已認識到免疫反應、轉錄因子核因子κB( NF-κB) 的激活、細胞因子、中性粒細胞的激活和肺泡滲入、凝血級聯反應、腎素-血管緊張素系統等多種因素構成的復雜網絡參與急性肺損傷/急性呼吸窘迫綜合征( ALI/ARDS) 的發病過程[ 1-5] 。雖然膿毒癥、創傷、肺炎等ALI/ARDS誘發因素很常見, 但僅有部分病人發生ALI/ARDS, 并且具有相似臨床特征的ALI/ARDS病人可有截然不同的結果, 這種異質性引起研究者對影響ALI/ARDS 易感性和預后的遺傳因子進行鑒別的濃厚興趣[ 6] 。由于數量龐大的表現型變異, 不完全的基因外顯率、復雜的基因-環境相互作用及高度可能的基因座不均一性而使ALI 遺傳學的研究受到挑戰[ 7] 。近年來基因組學技術被應用于ALI/ARDS 發病機制的研究, 加深了人們對ALI/ARDS的認識并有可能發展出新的治療策略以降低其發病率和病死率。
ObjectiveTo analyze the expression profile changes of osteogenic-related genes during spontaneous calcification of rat bone marrow mesenchymal stem cells (BMSCs). MethodsBMSCs were isolated from 3-day-old healthy Sprague Dawley rats;cells at the 4th generation were used to establish the spontaneous calcification model in vitro. Spontaneous calcification process was recorded by inverted phase contrast microscope observation and alizarin red staining after 7 and 14 days of culture. For gene microarray analysis, cell samples were collected at 0, 7, and 14 days after culture; the differentially expressed genes were analyzed by bioinformatics methods and validated by real-time quantitative PCR (RT-qPCR) assay. ResultsRat BMSCs calcified spontaneously in vitro. When cultured for 7 days, the cells began to aggregate and were weakly positive for alizarin red staining. After 14 days of culture, obvious cellular aggregation and typical mineralized nodules were observed, the mineralized nodules were brightly positive for alizarin red staining. A total of 576 gene probe-sets expressed differentially during spontaneous calcification, corresponding 378 rat genes. Among them, 359 gene probe-sets expressed differentially between at 0 and 7 days, while only 13 gene probe-sets expressed differentially between at 7 and 14 days. The 378 differentially expressed genes were divided into 6 modes according to their expression profiles. Moreover, according to their biological functions, differentially expressed genes related to bone cell biology could be classified into 7 major groups:angiogenesis, apoptosis, bone-related genes, cell cycle, development, cell communication, and signal pathways related to osteogenic differentiation. In cell cycle group, 12 down-regulated genes were linked with each other functionally. Matrix metalloproteinase 13 (Mmp13), secreted phosphoprotein 1 (Spp1), Cxcl12, Mmp2, Mmp3, Apoe, and Itga7 had more functional connections with other genes. The results of genes Spp1, Mgp, Mmp13, Wnt inhibitory factor 1, Cxcl12, and cyclin A2 by RT-qPCR were consistent with that of gene microarray. ConclusionThe first 7 days after rat BMSCs were seeded are a key phase determining the fate of spontaneous calcification. Multiple genes related with cell communication, bone-related genes, cell cycle, transforming growth factor-β signaling pathway, mitogen-activated protein kinase signaling pathway, and Wnt signaling pathway are involved during spontaneous calcification.
Objective To explore the mode and role of differential expression of circular RNAs (circRNAs) in myelodysplastic syndrome (MDS). Methods We preprocessed and analyzed the circRNA expression profile datasets GSE163386, GSE94591, and GSE81173 in the GEO (Gene Expression Omnibus) database. By using the circBank database and the ENCORI, miRDB, and miRWalk databases to predict microRNAs (miRNAs) that interacted with differentially expressed circRNAs and messenger RNAs (mRNAs), the circRNA-miRNA-mRNA axis was constructed. We retrieved miRNAs related to MDS in PubMed and further obtained competing endogenous RNA (ceRNA) networks related to MDS by taking intersections. Results Through analysis, 128 differentially expressed circRNAs were identified, 48 highly expressed, and 80 low expressed. Among differentially expressed circRNAs with multiple differences>10, 3 were upregulated and 11 were downregulated. Through analysis, 101 differentially expressed mRNA were identified, with 9 upregulated and 92 downregulated. Intersecting with the MDS related miRNAs retrieved by PubMed, we further obtained the MDS related ceRNA network, namely circRNA (has_circ_0061137)-miRNA (has-miR-16-5p)-mRNA (RUBCNL, TBC1D9, SLC16A6) and circRNA (has_circ_0061137)-miRNA (has-miR-125b-5p)-mRNA (CCR5, SLC16A6, IRF4), all of which were downregulated. Conclusion The ceRNA networks revealed in this study may help elucidate the circRNA mechanism in MDS.