Comprehensive integration of single-cell transcriptomic data illuminates the regulatory network architecture of plant cell fate specification

Plant morphogenesis relies on precise gene expression programs at the proper time and position which is orchestrated by transcription factors(TFs)in intricate regulatory networks in a cell-type specific manner.Here we introduced a comprehensive single-cell transcriptomic atlas of Arabidopsis seedlings.This atlas is the result of meticulous integration of 63 previously published scRNA-seq datasets,addressing batch effects and conserving biological variance.This integration spans a broad spectrum of tissues,including both below-and above-ground parts.Utilizing a rigorous approach for cell type annotation,we identified 47 distinct cell types or states,largely expanding our current view of plant cell compositions.We systematically constructed cell-type specific gene regulatory networks and uncovered key regulators that act in a coordinated manner to control cell-type specific gene expression.Taken together,our study not only offers extensive plant cell atlas exploration that serves as a valuable resource,but also provides molecular insights into gene-regulatory programs that varies from different cell types.

Combining single-cell RNA-sequencing and bulk data to reveal immunity-related genes expression pattern in the systemic lupus erythematosus and target organ kidney

Background:Systemic lupus erythematosus(SLE)is a complex chronic autoimmune disease with no known cure.However,the regulatory mechanism of immunity-related genes is not fully understood in SLE.In order to explore new therapeutic targets,we used bioinformatical methods to analyze a series of data.Methods:After downloading and processing the data from Gene Expression Omnibus database,the differentially expressed genes of SLE were analyzed.CIBERSORT algorithm was used to analyze the immune infiltration of SLE.Based on single-cell RNA-sequencing data,the role of immune-related genes in SLE and its target organ(kidney)were analyzed.Key transcription factors affecting immune-related genes were identified.Cell-cell communication networks in SLE were analyzed.Results:In total,15 hub genes and 4 transcription factors were found in the bulk data.Monocytes and macrophages in GSE81622(SLE)showed more infiltration.There were four cell types were annotated in scRNA sequencing dataset(GSE135779),as follows T cells,monocyte,NK cells and B cells.Immunity-related genes were overexpressed in monocytes.Conclusion:The present study shows that immune-related genes affect SLE through monocytes and play an important role in target organ renal injury.

Comparative study of microarray and experimental data on Schwann cells in peripheral nerve degeneration and regeneration: big data analysis

A Schwann cell has regenerative capabilities and is an important cell in the peripheral nervous system.This microarray study is part of a bioinformatics study that focuses mainly on Schwann cells. Microarray data provide information on differences between microarray-based and experiment-based gene expression analyses. According to microarray data, several genes exhibit increased expression(fold change) but they are weakly expressed in experimental studies(based on morphology, protein and mRNA levels). In contrast, some genes are weakly expressed in microarray data and highly expressed in experimental studies;such genes may represent future target genes in Schwann cell studies. These studies allow us to learn about additional genes that could be used to achieve targeted results from experimental studies. In the current big data study by retrieving more than 5000 scientific articles from PubMed or NCBI, Google Scholar, and Google, 1016(up-and downregulated) genes were determined to be related to Schwann cells. However,no experiment was performed in the laboratory; rather, the present study is part of a big data analysis. Our study will contribute to our understanding of Schwann cell biology by aiding in the identification of genes.Based on a comparative analysis of all microarray data, we conclude that the microarray could be a good tool for predicting the expression and intensity of different genes of interest in actual experiments.

For robust big data analyses:a collection of 150 important pro-metastatic genes

Metastasis is the greatest contributor to cancer?related death.In the era of precision medicine,it is essential to predict and to prevent the spread of cancer cells to significantly improve patient survival.Thanks to the application of a variety of high?throughput technologies,accumulating big data enables researchers and clinicians to identify aggressive tumors as well as patients with a high risk of cancer metastasis.However,there have been few large?scale gene collection studies to enable metastasis?related analyses.In the last several years,emerging efforts have identi?fied pro?metastatic genes in a variety of cancers,providing us the ability to generate a pro?metastatic gene cluster for big data analyses.We carefully selected 285 genes with in vivo evidence of promoting metastasis reported in the literature.These genes have been investigated in different tumor types.We used two datasets downloaded from The Cancer Genome Atlas database,specifically,datasets of clear cell renal cell carcinoma and hepatocellular carcinoma,for validation tests,and excluded any genes for which elevated expression level correlated with longer overall survival in any of the datasets.Ultimately,150 pro?metastatic genes remained in our analyses.We believe this collection of pro?metastatic genes will be helpful for big data analyses,and eventually will accelerate anti?metastasis research and clinical intervention.

Identification of candidate genes controlling fiber quality traits in upland cotton through integration of meta-QTL,significant SNP and transcriptomic data

Background:Meta-analysis of quantitative trait locus(QTL)is a computational technique to identify consensus QTL and refine QTL positions on the consensus map from multiple mapping studies.The combination of meta-QTL intervals,significant SNPs and transcriptome analysis has been widely used to identify candidate genes in various plants.Results:In our study,884 QTLs associated with cotton fiber quality traits from 12 studies were used for meta-QTL analysis based on reference genome TM-1,as a result,74 meta-QTLs were identified,including 19 meta-QTLs for fiber length;18 meta-QTLs for fiber strength;11 meta-QTLs for fiber uniformity;11 meta-QTLs for fiber elongation;and 15 meta-QTLs for micronaire.Combined with 8589 significant single nucleotide polymorphisms associated with fiber quality traits collected from 15 studies,297 candidate genes were identified in the meta-QTL intervals,20 of which showed high expression levels specifically in the developing fibers.According to the function annotations,some of the 20 key candidate genes are associated with the fiber development.Conclusions:This study provides not only stable QTLs used for marker-assisted selection,but also candidate genes to uncover the molecular mechanisms for cotton fiber development.

Incorporating heterogeneous biological data sources in clustering gene expression data

In this paper, a similarity measure between genes with protein-protein interactions is pro-posed. The chip-chip data are converted into the same form of gene expression data with pear-son correlation as its similarity measure. On the basis of the similarity measures of protein- protein interaction data and chip-chip data, the combined dissimilarity measure is defined. The combined distance measure is introduced into K-means method, which can be considered as an improved K-means method. The improved K-means method and other three clustering methods are evaluated by a real dataset. Per-formance of these methods is assessed by a prediction accuracy analysis through known gene annotations. Our results show that the improved K-means method outperforms other clustering methods. The performance of the improved K-means method is also tested by varying the tuning coefficients of the combined dissimilarity measure. The results show that it is very helpful and meaningful to incorporate het-erogeneous data sources in clustering gene expression data, and those coefficients for the genome-wide or completed data sources should be given larger values when constructing the combined dissimilarity measure.

Analysis of Gene Expression Profiles of Rice Mutant SLR1 Based on Microarray Data

Gibberellins are an important class of plant hormones.They play an important regulatory role in all stages of growth and development of higher plants.The use of mutants to study gibberellin metabolism and signal transduction pathways is currently a research hotspot.This article takes the data of Affymetrix chips of rice as an example,bioinformatics method was used to study rice SLR1 mutant and mine differentially expressed wild-type genes,thus exploring the expression regulation network of gibberellin signaling pathway-related genes.

Deep Learning Enabled Microarray Gene Expression Classification for Data Science Applications

In bioinformatics applications,examination of microarray data has received significant interest to diagnose diseases.Microarray gene expression data can be defined by a massive searching space that poses a primary challenge in the appropriate selection of genes.Microarray data classification incorporates multiple disciplines such as bioinformatics,machine learning(ML),data science,and pattern classification.This paper designs an optimal deep neural network based microarray gene expression classification(ODNN-MGEC)model for bioinformatics applications.The proposed ODNN-MGEC technique performs data normalization process to normalize the data into a uniform scale.Besides,improved fruit fly optimization(IFFO)based feature selection technique is used to reduce the high dimensionality in the biomedical data.Moreover,deep neural network(DNN)model is applied for the classification of microarray gene expression data and the hyperparameter tuning of the DNN model is carried out using the Symbiotic Organisms Search(SOS)algorithm.The utilization of IFFO and SOS algorithms pave the way for accomplishing maximum gene expression classification outcomes.For examining the improved outcomes of the ODNN-MGEC technique,a wide ranging experimental analysis is made against benchmark datasets.The extensive comparison study with recent approaches demonstrates the enhanced outcomes of the ODNN-MGEC technique in terms of different measures.

A Novel Soft Clustering Approach for Gene Expression Data

Gene expression data represents a condition matrix where each rowrepresents the gene and the column shows the condition. Micro array used todetect gene expression in lab for thousands of gene at a time. Genes encode proteins which in turn will dictate the cell function. The production of messengerRNA along with processing the same are the two main stages involved in the process of gene expression. The biological networks complexity added with thevolume of data containing imprecision and outliers increases the challenges indealing with them. Clustering methods are hence essential to identify the patternspresent in massive gene data. Many techniques involve hierarchical, partitioning,grid based, density based, model based and soft clustering approaches for dealingwith the gene expression data. Understanding the gene regulation and other usefulinformation from this data can be possible only through effective clustering algorithms. Though many methods are discussed in the literature, we concentrate onproviding a soft clustering approach for analyzing the gene expression data. Thepopulation elements are grouped based on the fuzziness principle and a degree ofmembership is assigned to all the elements. An improved Fuzzy clustering byLocal Approximation of Memberships (FLAME) is proposed in this workwhich overcomes the limitations of the other approaches while dealing with thenon-linear relationships and provide better segregation of biological functions.

Modeling viscosity of methane,nitrogen,and hydrocarbon gas mixtures at ultra-high pressures and temperatures using group method of data handling and gene expression programming techniques

Accurate gas viscosity determination is an important issue in the oil and gas industries.Experimental approaches for gas viscosity measurement are timeconsuming,expensive and hardly possible at high pressures and high temperatures(HPHT).In this study,a number of correlations were developed to estimate gas viscosity by the use of group method of data handling(GMDH)type neural network and gene expression programming(GEP)techniques using a large data set containing more than 3000 experimental data points for methane,nitrogen,and hydrocarbon gas mixtures.It is worth mentioning that unlike many of viscosity correlations,the proposed ones in this study could compute gas viscosity at pressures ranging between 34 and 172 MPa and temperatures between 310 and 1300 K.Also,a comparison was performed between the results of these established models and the results of ten wellknown models reported in the literature.Average absolute relative errors of GMDH models were obtained 4.23%,0.64%,and 0.61%for hydrocarbon gas mixtures,methane,and nitrogen,respectively.In addition,graphical analyses indicate that the GMDH can predict gas viscosity with higher accuracy than GEP at HPHT conditions.Also,using leverage technique,valid,suspected and outlier data points were determined.Finally,trends of gas viscosity models at different conditions were evaluated.

Challenges Analyzing RNA-Seq Gene Expression Data

The analysis of messenger Ribonucleic acid obtained through sequencing techniques (RNA-se- quencing) data is very challenging. Once technical difficulties have been sorted, an important choice has to be made during pre-processing: Two different paths can be chosen: Transform RNA- sequencing count data to a continuous variable or continue to work with count data. For each data type, analysis tools have been developed and seem appropriate at first sight, but a deeper analysis of data distribution and structure, are a discussion worth. In this review, open questions regarding RNA-sequencing data nature are discussed and highlighted, indicating important future research topics in statistics that should be addressed for a better analysis of already available and new appearing gene expression data. Moreover, a comparative analysis of RNAseq count and transformed data is presented. This comparison indicates that transforming RNA-seq count data seems appropriate, at least for differential expression detection.

Mining and analysis of intracranial aneurysms formation and fracture-related genes based on Gene Expression Omnibus database

Objective: To explore potential genes associated with the formation and rupture of intracranial aneurysms based on the Gene Expression Omnibus (GEO) database. Methods: A total of 133 mRNA microarrays were collected from the GEO database. Differential mRNA gene analysis was performed on the data of each group in the GEO2R platform, and the common differential genes were screened and the gene ontology enrichment analysis and the Kyoto Gene and Genomic Encyclopedia pathway enrichment analysis were completed. The screened differential genes were introduced into the String online database to obtain the interaction between the proteins encoded by the differential genes. Results: Forty-two common differential genes were screened, and the main biological processes involved included the transcriptional regulation of oxidative stress, the positive regulation of chemokine production, and the positive regulation of autophagy of giant cells by RNA polymerase II promoter. Molecular functions included protein binding, RNA polymerase II transcriptional co-repressor activity, transcriptional activator activity, and protein kinase C binding. The main signal pathways covered included hypoxia-inducible factor-1 signaling pathway, glucagon signaling pathway, and metabolic pathway signaling pathway. Conclusions: The formation and rupture of the intracranial aneurysm may be initially screened with amidoxime reduction component 1, tumor necrosis factor-α-inducible protein 6, haptoglobin, mast cell membrane-expressing protein 1, zipper containing kinase, phospholipase Cβ4 and blood and nervous system expression factor-1. In addition to the previously knownintracranial aneurysms mechanisms, cellular autophagy and hypoxia inducible factor-1 pathway may also be involved in the formation of intracranial aneurysms.

A Survey on Acute Leukemia Expression Data Classification Using Ensembles

Acute leukemia is an aggressive disease that has high mortality rates worldwide.The error rate can be as high as 40%when classifying acute leukemia into its subtypes.So,there is an urgent need to support hematologists during the classification process.More than two decades ago,researchers used microarray gene expression data to classify cancer and adopted acute leukemia as a test case.The high classification accuracy they achieved confirmed that it is possible to classify cancer subtypes using microarray gene expression data.Ensemble machine learning is an effective method that combines individual classifiers to classify new samples.Ensemble classifiers are recognized as powerful algorithms with numerous advantages over traditional classifiers.Over the past few decades,researchers have focused a great deal of attention on ensemble classifiers in a wide variety of fields,including but not limited to disease diagnosis,finance,bioinformatics,healthcare,manufacturing,and geography.This paper reviews the recent ensemble classifier approaches utilized for acute leukemia gene expression data classification.Moreover,a framework for classifying acute leukemia gene expression data is proposed.The pairwise correlation gene selection method and the Rotation Forest of Bayesian Networks are both used in this framework.Experimental outcomes show that the classification accuracy achieved by the acute leukemia ensemble classifiers constructed according to the suggested framework is good compared to the classification accuracy achieved in other studies.

The expression significance of RACGAP1 gene in hepatocellular carcinoma and its prognostic effect were analyzed based on TCGA database

Objective:To explore the expression and clinical significance of RACGAP1 gene in hepatocellular carcinoma.Methods:Data about RACGAP1 gene and clinic pathological data in liver cancer were retrieved from The Cancer Genome Atlas(TCGA).The relationship between the expression of RACGAP1 gene and clinic pathological parameters,and prognosis were analyzed by R 2.15.3 software.The association between RACGAP1 gene expression and prognosis of liver cancer patients was analyzed by Kaplan-Meier survival function analysis and Cox regression analysis.Results:TCGA database was used to collect 235 cases of liver cancer with clinical pathological parameters and their corresponding RACGAP1 expression levels.After the incomplete cases and those with no detailed pathological parameters were excluded,and it was found that RACGAP1 was highly expressed in liver cancer tissues.Meanwhile,the expression of RACGAP1 in patients with liver cancer in the TCGA tumor database was further analyzed with the matching clinical data parameters.The expression level of RACGAP1 was significantly correlated with the pathological grade and T stage of liver cancer patients(all P<0.05),but was not significantly correlated with American Joint Committee on Cancer(AJCC)pathological stage and gender(P>0.05).There was a significant correlation between RACGAP1 expression level and overall survival(OS)in patients with liver cancer(P<0.05),and the overall survival time of patients with low expression was better than that of patients with high expression(P<0.05).Cox regression was used to analyze the correlation between T stage,M stage,N stage and RACGAP1 expression in patients with hepatocellular carcinoma(HCC),and RACGAP1 became an independent prognostic factor in patients with HCC(P<0.05).Conclusion:Based on the tumor-related gene information in the public database TCGA,RACGAP1 gene is highly expressed in liver cancer tissues and becomes an independent prognostic factor of liver cancer,which is expected to become an important therapeutic target of drug therapy for liver cancer.

Genome-Wide Identification of Genes Responsive to ABA and Cold/Salt Stresses in Gossypium hirsutum by Data-Mining and Expression Pattern Analysis

For making better use of nucleic acid resources of Gossypium hirsutum, a data-mining method was used to identify putative genes responsive to various abiotic stresses in G. hirsutum. Based on the compiled database including genes involved in abiotic stress response in Arabidopsis thaliana and the comprehensive analysis tool of GENEVESTIGATOR v3, 826 genes up-regulated or down-regulated significantly in roots or leaves during salt or cold treatment in Arabidopsis were identified. As compared to these 826 Arabidopsis genes annotated, 38 homologous expressed sequence tags （ESTs） from G. hirsutum were selected randomly and their expression patterns were studied using a quantitative real-time reverse transcription-polymerase chain reaction method. Among these 38 ESTs, about 55% of the genes （21 of 38） were different in response to ABA between cotton and Arabidopsis, whereas 70% of genes had similar responses to cold and salt treatments, and some of them which had not been characterized in Arabidopsis are now being investigated in gene function studies. According to these results, this approach of analyzing ESTs appears effective in large-scale identification of cotton genes involved in abiotic stress and might be adopted to determine gene functions in various biologic processes in cotton.

Identification of key genes and biological pathways in lung adenocarcinoma by integrated bioinformatics analysis

BACKGROUND The objectives of this study were to identify hub genes and biological pathways involved in lung adenocarcinoma(LUAD)via bioinformatics analysis,and investigate potential therapeutic targets.AIM To determine reliable prognostic biomarkers for early diagnosis and treatment of LUAD.METHODS To identify potential therapeutic targets for LUAD,two microarray datasets derived from the Gene Expression Omnibus(GEO)database were analyzed,GSE3116959 and GSE118370.Differentially expressed genes(DEGs)in LUAD and normal tissues were identified using the GEO2R tool.The Hiplot database was then used to generate a volcanic map of the DEGs.Weighted gene co-expression network analysis was conducted to cluster the genes in GSE116959 and GSE-118370 into different modules,and identify immune genes shared between them.A protein-protein interaction network was established using the Search Tool for the Retrieval of Interacting Genes database,then the CytoNCA and CytoHubba components of Cytoscape software were used to visualize the genes.Hub genes with high scores and co-expression were identified,and the Database for Annotation,Visualization and Integrated Discovery was used to perform enrichment analysis of these genes.The diagnostic and prognostic values of the hub genes were calculated using receiver operating characteristic curves and Kaplan-Meier survival analysis,and gene-set enrichment analysis was conducted.The University of Alabama at Birmingham Cancer data analysis portal was used to analyze relationships between the hub genes and normal specimens,as well as their expression during tumor progression.Lastly,validation of protein expression was conducted on the identified hub genes via the Human Protein Atlas database.RESULTS Three hub genes with high connectivity were identified;cellular retinoic acid binding protein 2(CRABP2),matrix metallopeptidase 12(MMP12),and DNA topoisomerase II alpha(TOP2A).High expression of these genes was associated with a poor LUAD prognosis,and the genes exhibited high diagnostic value.CONCLUSION Expression levels of CRABP2,MMP12,and TOP2A in LUAD were higher than those in normal lung tissue.This observation has diagnostic value,and is linked to poor LUAD prognosis.These genes may be biomarkers and therapeutic targets in LUAD,but further research is warranted to investigate their usefulness in these respects.

Hub genes and their key effects on prognosis of Burkitt lymphoma

BACKGROUND Burkitt lymphoma(BL)is an exceptionally aggressive malignant neoplasm that arises from either the germinal center or post-germinal center B cells.Patients with BL often present with rapid tumor growth and require high-intensity multidrug therapy combined with adequate intrathecal chemotherapy prophylaxis,however,a standard treatment program for BL has not yet been established.It is important to identify biomarkers for predicting the prognosis of BLs and discriminating patients who might benefit from the therapy.Microarray data and sequencing information from public databases could offer opportunities for the discovery of new diagnostic or therapeutic targets.AIM To identify hub genes and perform gene ontology(GO)and survival analysis in BL.METHODS Gene expression profiles and clinical traits of BL patients were collected from the Gene Expression Omnibus database.Weighted gene co-expression network analysis(WGCNA)was applied to construct gene co-expression modules,and the cytoHubba tool was used to find the hub genes.Then,the hub genes were analyzed using GO and Kyoto Encyclopedia of Genes and Genomes analysis.Additionally,a Protein-Protein Interaction network and a Genetic Interaction network were constructed.Prognostic candidate genes were identified through overall survival analysis.Finally,a nomogram was established to assess the predictive value of hub genes,and drug-gene interactions were also constructed.RESULTS In this study,we obtained 8 modules through WGCNA analysis,and there was a significant correlation between the yellow module and age.Then we identified 10 hub genes(SRC,TLR4,CD40,STAT3,SELL,CXCL10,IL2RA,IL10RA,CCR7 and FCGR2B)by cytoHubba tool.Within these hubs,two genes were found to be associated with OS(CXCL10,P=0.029 and IL2RA,P=0.0066)by survival analysis.Additionally,we combined these two hub genes and age to build a nomogram.Moreover,the drugs related to IL2RA and CXCL10 might have a potential therapeutic role in relapsed and refractory BL.CONCLUSION From WGCNA and survival analysis,we identified CXCL10 and IL2RA that might be prognostic markers for BL.

血管样本生物信息学分析鉴定烟雾病相关的潜在关键基因

目的本研究对烟雾病患者血管样本的差异表达基因(differentially expressed genes,DEGs)进行生物信息学鉴定和分析,旨在探讨烟雾病的潜在发病机制。方法本研究以烟雾病和颈内动脉瘤患者大脑血管样本为研究对象。利用R语言线性模型微阵列数据(linear models for microarray data,limma)分析包对基因表达综合数据库(gene expression omnibus,GEO)中的GSE141025数据集进行分析,该数据集涵盖4例烟雾病患者和4例颈内动脉瘤患者的大脑中动脉和颞浅动脉样本各1个,共计16个样本。选择烟雾病患者的大脑中动脉、颞浅动脉及颈内动脉瘤患者的颞浅动脉共12个样本进行DEGs筛选。通过R语言功能富集分析工具包clusterProfiler,对筛选出的DEGs进行基因本体(gene ontology,GO)富集分析和京都基因与基因组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)通路分析。利用STRING数据库构建蛋白质-蛋白质相互作用(proteinprotein interaction,PPI)网络,并使用网络可视化软件Cytoscape进行蛋白质网络的可视化和枢纽基因筛选。结果本研究在烟雾病患者的大脑中动脉与颞浅动脉样本间鉴定出138个DEGs,包括18个上调基因和120个下调基因。GO富集分析显示,以上DEGs在细胞外基质、受体配体活性和生长因子活性等方面显著富集,可能与烟雾病相关的血管病变和神经保护机制有关。KEGG通路分析提示,DEGs主要在酪氨酸代谢通路中富集。通过PPI网络分析,共筛选出9个枢纽基因,包括骨膜蛋白(periostin,POSTN)、脑源性神经营养因子(brain derived neurotrophic factor,BDNF)、血小板衍生生长因子受体α(platelet derived growth factor receptor alpha,PDGFRA)、Thy-1细胞表面抗原(Thy-1 cell surface antigen,THY1)、ⅩⅤ型胶原蛋白α1链(collagen typeⅩⅤalpha 1 chain,COL15A1)、成纤维细胞生长因子7(fibroblast growth factor 7,FGF7)、光蛋白聚糖(l umi can,LUM)、层粘连蛋白α2亚基(laminin subunit alpha 2,LAMA2)和RELN(reelin)。此外,上调基因delta样典型Notch配体4(delta like canonical Notch ligand 4,DLL4)在本研究中首次被发现可能在烟雾病中扮演重要角色,或与烟雾病的病理性血管生成有关。结论细胞外基质、生长因子及其受体的表达失调等可能参与烟雾病的发病过程。DEGs分析筛选出的枢纽基因(POSTN、BDNF、PDGFRA、THY1、COL15A1、FGF7、LUM、LAMA2、RELN)以及DLL4可能在烟雾病的病理形成过程中发挥作用。

基于生物信息学方法的子宫腺肌病关键靶点与中药筛选研究

目的通过数据挖掘结合生物信息学技术,比较子宫腺肌病患者与正常人群的基因芯片数据库,查询差异基因,并分析差异基因参与的生物过程和潜在靶点,预测治疗子宫腺肌病的潜在药物。方法利用基因表达数据库(GEO)筛选出GSE78851和GSE7307两个数据集,通过GEO2R分析芯片的共有差异基因,应用Omicshare制作差异基因热图,借助DAVID数据库对差异基因进行基因本体分析(GO)和京都基因与基因组百科全书(KEGG)分析,使用String构建蛋白质-蛋白质相互作用网络,筛选核心基因,再进一步利用Cytoscape 3.9.0的MCODE插件对差异表达基因进行模块分析。通过医学本体信息检索平台筛选治疗子宫腺肌病的中医药物。结果筛选出433个共有差异基因,差异表达基因主要富集在细胞黏附、对雌二醇的反应、细胞凋亡等生物过程和氧化磷酸化、雌激素信号通路、PI3K-Akt等信号通路。共得出p53蛋白(TP53)、细胞色素C(CYCS)、肌动蛋白(ACTB)、Ⅳ型胶原酶(MMP2)、雌激素受体(ESR1)、CD44、丝裂原活化蛋白激酶3(MAPK3)、胰岛素样生长因子-Ⅰ(IGF-Ⅰ)、细胞质动力蛋白1重链1(DYNC1H1)、血小板反应蛋白-1(THBS1)、细胞色素氧化酶4Ⅱ(COX4Ⅱ)、NDUFAB1、人原纤维蛋白-1(FBN1)、ATP5A1、NDUFB7、ATP5G1等16个核心基因。在Cytoscape 3.9.0中进行模块分析并得出两个重要模块,对其进行分析后发现模块1基因主要与氧化磷酸化信号通路相关,模块2基因主要与Rap1、Wnt信号通路相关。筛选得到丹参、当归、茯苓、三七、厚朴花、茶树根、蚕砂7味出现频次最高的中药。结论通过差异表达基因和富集分析得到的生物过程与信号通路可以为子宫腺肌病的潜在治疗靶点提供参考,通过核心基因映射得到了治疗子宫腺肌病的潜在药物,为子宫腺肌病的治疗提供新的思路。

Gene Ontology在生物数据整合中的应用

异构数据的高效整合,在生物数据呈爆炸性增长、生物数据库复杂度不断增加的今天,具有重要的理论价值和实际意义。该文基于BioDW——一个整合的生物信息学数据仓库平台,利用统一的GeneOntology语义模型,建立异构数据库之间的语义链接,在概念和联系层次上有效地解决了生物异构数据的整合问题,实现了对生物数据智能化的多重、复合和交叉检索,为生物信息的进一步研究奠定了坚实的基础。