筛选差异表达基因和蛋白质的方法进展

分离和鉴定差异表达基因和蛋白质不仅有助于发现基因和蛋白质的功能,更有助于揭示某些疾病的发生机理.目前筛选差异表达基因的方法主要有差异显示PCR方法(differential display RT-PCR,DDRT-PCR)、消减杂交法(subtractive hybridization,SH)、基因芯片技术(DNA chip technique)和基因表达的系统分析(serial analysis of gene expression,SAGE)等,其中消减杂交法中又先后建立了代表性差异分析技术(representational difference analysis,RDA)、抑制消减杂交法(suppression subtractive hybridization,SSH)和获得全长基因的消减杂交法(full-length-gene-obtainable subtractive hybridization).筛选差异表达蛋白质的方法主要有双向电泳技术(two-dimentional gel electrophoresis)和噬菌体全套抗体库技术(phage display antibody repertoire library technique).这些方法各有特点,各有利弊,研究者可根据自己的需要选择适合于自己的方法.     

Avr／Cf互作诱发番茄Pto互作因子Pti编码基因的表达

植物抗病基因Pto和Cf产物具完全不同的结构域,其细胞定位也不同。二决定的抗病性产生机制的异同令人关注。采用两种过敏性反应（hypersensitive response,HR）产生系统研究了Pro互作蛋白Pti4,Pti5和Pti6编码基因在Avr/Cf互作中的时序表达：（1）通过杂交方法获取同含互补基因对Avr4/Cf-4和Avr9/Cf-9的番茄（Lycoper-sicon esculentum Mill.）种子,常温下这些种子发芽后形成的苗产生HR坏死斑。（2）先将Avr/Cf苗置于33℃下培养,此时番茄苗生长正常,不形成HR坏死斑。然后将温度降至25℃,数小时内这些苗即形成HR坏死斑。不同方法研究结果均表明,随Avr/Cf苗中HR坏死斑的形成,Pti4、Pti5和Pti6均受显诱导表达。但它们的表达水平和动态不同,这些结果表明,这些Pti在功能上互补,可能同时涉及Pto和Cf决定性抗性的调节作用。     

根据蛋白质互作网络预测前列腺癌相关蛋白质的细致功能

对于功能部分已知的前列腺癌相关蛋白质,提出了一种基于Gent Ontology的功能特异的子网构建方法来细化其功能注释。结果显示该方法能够以很高的精确率为前列腺癌相关蛋白质预测更为精细的功能。预测的相关蛋白质的功能对于指导实验研究前列腺癌的分子机制具有重要的价值。     

基于蛋白质互作知识的生物学通路扩充新方法

生物学通路被广泛应用于基因功能学研究, 但现有的生物学通路知识并不完善, 仍需进一步扩充。生物信息学预测为通路扩充提供了一种有效且经济的途径。文章提出了一种融合蛋白质-蛋白质互作知识以及Gene Ontology(GO)数据库信息进行基因通路预测的新方法。首先选取目标基因在蛋白质-蛋白质互作层面上的邻居所在的Kyoto Encyclopedia of Genes and Genomes(KEGG)通路为候选通路, 然后通过检验候选通路中的基因是否在与目标基因关联的GO节点富集来判断目标基因的通路归属。分别利用Human Protein Reference Database (HPRD)和Biological General Repository for Interaction Datasets(BioGRID)数据库中的蛋白质-蛋白质互作信息进行预测。结果表明, 在两套数据中, 随着互作邻居个数的增加, 预测的平均准确率(在所有目标基因注释的通路中被成功预测的比例)及相对准确率(在至少有一个注释通路被成功预测的基因集中, 所有注释通路均被预测正确的基因所占的比例)均呈现上升趋势。当互作邻居个数达到22时, 预测的平均准确率分别达到96.2%(HPRD)和96.3%(BioGRID), 而相对准确率分别为93.3%(HPRD)和84.1%(BioGRID)。进一步利用新版数据库对旧版数据库中被更新的89个基因进行验证, 至少有一个更新通路被预测正确的基因有50个, 其中43个基因的更新通路被完全正确预测, 相对准确率为86.0%。这些结果显示该方法是一种可靠且有效的通路扩充方法。     

基于蛋白质互作网络挖掘自闭症谱系障碍的功能模块与核心基因

自闭症谱系障碍(autism spectrum disorder, ASD)是一种具有高遗传性、临床异质性和生物复杂性的神经行为障碍类疾病。为挖掘ASD发生发展过程中的功能模块与核心基因,本文从自闭症谱系障碍疾病数据库获取ASD相关基因;利用STRING数据库构建ASD相关基因的蛋白质互作网络;通过MCODE算法对蛋白质互作网络进行模块分析并筛选核心基因;最后对各功能模块进行KEGG通路分析,根据富集到的通路类别评估功能模块之间的相互作用。结果显示, 3个疾病基因数据库筛选出182个共有基因,构建的蛋白质互作网络包含171个节点和1 041条边,其中NRXN1、GRIN2B、GRIN2A、DLG4、NLGN3、MECP2、CNTNAP2、BDNF、NLGN4X、FMR1等23个基因具有较高的连通度(degree)。从蛋白质互作网络中分析得到5个功能模块,包括68个核心基因。KEGG富集分析发现功能模块参与多个生物学通路,包括细胞黏附分子、钙离子通路、神经活性的配体-受体相互作用、多巴胺能神经突触等。分析结果提示,挖掘的ASD功能模块和核心基因大多集中在神经元活动、信号分子和信号传导等,且各模块相互作用共同影响ASD的发生发展。     

人类蛋白质结构互作网络——结构域对网络拓扑与蛋白质功能的影响

高通量的蛋白质互作数据与结构域互作数据的出现,使得在蛋白质组学领域内研究人类蛋白质结构互作网络,进一步揭示蛋白质结构与功能间的潜在关系成为可能.蛋白质上广泛分布的结构域被认为是蛋白质结构、功能以及进化的基本功能单元.然而,结合蛋白质的结构信息(例如蛋白质结构域数目、长度和覆盖率等)来研究这些表象后的内部机制仍然面临着挑战.将蛋白质分为单结构域蛋白质与多结构域蛋白质,并进一步结合蛋白质互作信息与结构域互作信息构建了人类蛋白质结构互作网络;通过与人类蛋白质互作网络进行比较,研究了人类蛋白质结构互作网络的特殊结构特征;对于单结构域蛋白质与多结构域蛋白质,分别进行了功能富集分析、功能离散度分析以及功能一致性分析等.结果发现,将结构域互作信息综合考虑进来后,人类蛋白质结构互作网络可以提供更多的单纯的蛋白质互作网络无法提供的细节信息,揭示蛋白质互作网络的复杂性.     

根据蛋白质互作网络预测乳腺癌相关蛋白质的细致功能

乳腺癌是最为常见的恶性肿瘤之一。已有的关于乳腺癌相关蛋白质的功能注释比较宽泛, 制约了乳腺癌的后续研究工作。对于已知部分功能的乳腺癌相关蛋白质, 提出了一种结合Gene Ontology功能先验知识和蛋白质互作的方法, 通过构建功能特异的局部相互作用网络来预测乳腺癌相关蛋白质的细致功能。结果显示该方法能够以很高的精确率为乳腺癌相关蛋白质预测更为精细的功能。预测的相关蛋白质的功能对于指导实验研究乳腺癌的分子机制具有重要的价值。     

互作蛋白质与复合物亚基的基因表达相关性比较分析

利用在多种应激条件下酵母的基因表达谱数据 ,分别计算互作蛋白质及复合物亚基编码基因的表达相关性。结果发现 ,相对于随机对照组 ,互作蛋白质的编码基因与蛋白质复合物的编码基因表达相关性均显著 (P <0 .0 1) ,即互作蛋白质及复合物亚基有共表达的倾向。通过比较 ,进一步发现蛋白质复合物亚基的基因表达相关性显著高于互作蛋白质的基因表达相关性 (P <0 .0 1) ,这与复合物亚基之间功能联系强于定义不甚确切的互作蛋白之间功能联系现象吻合。     

基于蛋白质互作网络的COVID-19对肺动脉高压的影响机制

通过构建肺动脉高压差异基因和冠状病毒侵入人体后免疫反应相关基因的互作网络,探索COVID-19对肺动脉高压的影响机制。首先通过Meta分析挖掘肺动脉高压相关差异表达基因;其次通过SARS-CoV侵染人体后的基因表达数据,挖掘主要功能通路;最后构建肺动脉高压差异表达基因和冠状病毒主要功能通路基因的互作网络,挖掘网络的显著功能模块。发现肺动脉高压与血管平滑肌细胞、成纤细胞、T/B细胞免疫过程、转录调节因子通路、Toll样信号通路等密切相关,互作网络发现ITGAM、HBB、VCAM1、IL1R2等基因是COVID-19感染肺动脉高压患者的重要调节基因。通过肺动脉高压与冠状病毒感染机体后蛋白质互作网络探索了COVID-19对肺动脉高压的影响机制,为肺动脉高压感染COVID-19的研究及治疗提供了新思路。     

R-Avr基因互作介导的植物抗病性

Ｒ－Ａｖｒ基因互作介导的植物抗病性李汝刚范云六伍宁丰李茂林（中国农业科学院生物技术研究中心,北京１０００８１）植物在整个生长发育周期中面临着种类繁多的病原微生物袭击,但真正能引起植物病害的病原并不多,这表明植物能抵抗大多数病原微生物袭击,表现非亲合性...     

Screening and identification of hub genes in pancreatic cancer by integrated bioinformatics analysis

Pancreatic cancer (Pa) is a malignant tumor of the digestive tract with high degree of malignancy, this study aimed to obtain the hub genes in the tumorigenesis of Pa. Microarray datasets GSE15471, GSE16515, and GSE62452 were downloaded from Gene Expression Omnibus (GEO) database, GEO2R was conducted to screen the differentially expressed genes (DEGs), and functional enrichment analyses were carried out by Database for Annotation, Visualization and Integrated Discovery (DAVID). The protein-protein interaction (PPI) network was constructed with the Search Tool for the Retrieval of Interacting Genes (STRING), and the hub genes were identified by Cytoscape. Totally 205 DEGs were identified, consisting of 51 downregulated genes and 154 upregulated genes enriched in Gene Ontology terms including extracellular matrix (ECM) organization, collagen binding, cell adhesion, and pathways associated with ECM-receptor interaction, focal adhesion, and protein digestion. Two modules in the PPI were chosen and biological process analyses showed that the module genes were mainly enriched in ECM and cell adhesion. Twenty-four hub genes were confirmed, the survival analyses from the cBioPortal online platform revealed that topoisomerase (DNA) II α (TOP2A), periostin (POSTN), plasminogen activator, urokinase (PLAU), and versican (VCAN) may be involved in the carcinogenesis and progression of Pa, and the receiver-operating characteristic curves indicated their diagnostic value for Pa. Among them, TOP2A, POSTN, and PLAU have been previously reported as biomarkers for Pa, and far too little attention has been paid to VCAN. Analysis from R2 online platform showed that Pa patients with high VCAN expression were more sensitive to gemcitabine than those with low level, suggesting that VCAN may be an indicator to guide the use of the chemotherapeutic drug. In vitro experiments also showed that the sensitivity of the VCAN siRNA group to gemcitabine was lower than that of the control group. In conclusion, this study discerned hub genes and pathways related to the development of Pa, and VCAN was identified as a novel biomarker for the diagnose and therapy of Pa.     

Identification of potentially critical genes in the development of heart failure after ST-segment elevation myocardial infarction (STEMI)

Heart failure (HF) remains a common complication after acute ST-segment elevation myocardial infarction (STEMI). Here, we aim to identify critical genes related to the developed HF in patients with STEMI using bioinformatics analysis. The microarray data of GSE59867, including peripheral blood samples from nine patients with post-infarct HF and eight patients without post-infarct HF, were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between HF and non-HF groups were screened by LIMMA package. Functional enrichment analyses of DEGs were conducted, followed by construction of a protein-protein interaction (PPI) network. The dynamic messenger RNA (mRNA) level of the hub genes during the follow-up was analyzed to further elucidate their role in HF development. A total of 58 upregulated and 75 downregulated DEGs were screen out. They were mainly enriched in biological processes about inflammatory response, extracellular matrix organization, response to cAMP, immune response, and positive regulation of cytosolic calcium ion concentration. Pathway analysis revealed that the DEGs were also involved in hematopoietic cell lineage, pathways in cancer, and extracellular matrix-receptor interaction. In the PPI network consisting of 58 nodes and 72 interactions, CXCL8 (degree = 15), THBS1 (degree = 8), FOS (degree = 7), and ITGA2B (degree = 6) were identified as the hub genes. In the comparison of patients with and without post-infarct HF, the mRNA level of these hub genes were all higher within 30 days but reached similar at 6 months after STEMI. In conclusion, CXCL8, THBS1, FOS, and ITGA2B may play important roles in the development of HF after acute STEMI.     

Integrated network analysis and machine learning approach for the identification of key genes of triple-negative breast cancer

Triple-negative breast cancer (TNBC) has attracted more attention compared with other breast cancer subtypes due to its aggressive nature, poor prognosis, and chemotherapy remains the mainstay of treatment with no other approved targeted therapy. Therefore, the study aimed to discover more promising therapeutic targets and investigating new insights of biological mechanism of TNBC. Six microarray data sets consisting of 463 non-TNBC and 405 TNBC samples were mined from Gene Expression Omnibus. The data sets were integrated by meta-analysis and identified 1075 differentially expressed genes. Protein-protein interaction network was constructed which consists of 486 nodes and 1932 edges, where 29 hub genes were obtained with high topological measures. Further, 16 features (hub genes), 12 upregulated (AURKB, CCNB2, CDC20, DDX18, EGFR, ENO1, MYC, NUP88, PLK1, PML, POLR2F, and SKP2) and four downregulated ( CCND1, GLI3, SKP1, and TGFB3) were selected through machine learning correlation based feature selection method on training data set. A naïve Bayes based classifier built using the expression profiles of 16 features (hub genes) accurately and reliably classify TNBC from non-TNBC samples in the validation test data set with a receiver operating curve of 0.93 to 0.98. Subsequently, Gene Ontology analysis revealed that the hub genes were enriched in mitotic cell cycle processes and Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that they were enriched in cell cycle pathways. Thus, the identified key hub genes and pathways highlighted in the study would enhance the understanding of molecular mechanism of TNBC which may serve as potential therapeutic target.     

Investigation of hub genes and immune status in heart transplant rejection using endomyocardial biopsies

T cell‒mediated rejection (TCMR) and antibody-mediated rejection (ABMR) are severe post-transplantation complications for heart transplantation (HTx), whose molecular and immunological pathogenesis remains unclear. In the present study, the mRNA microarray data set GSE124897 containing 645 stable, 52 TCMR and 144 ABMR endomyocardial biopsies was obtained to screen for differentially expressed genes (DEGs) between rejected and stable HTx samples and to investigate immune cell infiltration. Functional enrichment analyses indicated roles of the DEGs primarily in immune-related mechanisms. Protein-protein interaction networks were then constructed, and ICAM1, CD44, HLA-A and HLA-B were identified as hub genes using the maximal clique centrality method. Immune cell infiltration analysis revealed differences in adaptive and innate immune cell populations between TCMR, ABMR and stable HTx samples. Additionally, hub gene expression levels significantly correlated with the degree and composition of immune cell infiltration in HTx rejection samples. Furthermore, drug-gene interactions were constructed, and 12 FDA-approved drugs were predicted to target hub genes. Finally, an external GSE2596 data set was used to validate the expression of the hub genes, and ROC curves indicated all four hub genes had promising diagnostic value for HTx rejection. This study provides a comprehensive perspective of molecular and immunological regulatory mechanisms underlying HTx rejection.     

Divergence pattern of duplicate genes in protein-protein interactions follows the power law

The impact of the biological network structures on the divergence between the two copies of one duplicate gene pair involved in the networks has not been documented on a genome scale. Having analyzed the most recently updated Database of Interacting Proteins (DIP) by incorporating the information for duplicate genes of the same age in yeast, we find that there was a highly significantly positive correlation between the level of connectivity of ancient genes and the number of shared partners of their duplicates in the protein-protein interaction networks. This suggests that duplicate genes with a low ancestral connectivity tend to provide raw materials for functional novelty, whereas those duplicate genes with a high ancestral connectivity tend to create functional redundancy for a genome during the same evolutionary period. Moreover, the difference in the number of partners between two copies of a duplicate pair was found to follow a power-law distribution. This suggests that loss and gain of interacting partners for most duplicate genes with a lower level of ancestral connectivity is largely symmetrical, whereas the "hub duplicate genes" with a higher level of ancient connectivity display an asymmetrical divergence pattern in protein-protein interactions. Thus, it is clear that the protein-protein interaction network structures affect the divergence pattern of duplicate genes. Our findings also provide insights into the origin and development of biological networks.     

Understanding protein-protein interactions by genetic suppression

Protein-protein interactions influence many cellular processes and it is increasingly being felt that even a weak and remote interplay between two subunits of a protein or between two proteins in a complex may govern the fate of a particular biochemical pathway. In a bacterial system where the complete genome sequence is available, it is an arduous task to assign function to a large number of proteins. It is possible that many of them are peripherally associated with a cellular event and it is very difficult to probe such interaction. However, mutations in the genes that encode such proteins (primary mutations) are useful in these studies. Isolation of a suppressor or a second-site mutation that restores the phenotype abolished by the primary mutation could be an elegant yet simple way to follow a set of interacting proteins. Such a reversion site need not necessarily be geometrically close to the primary mutation site.     

Approach of the functional evolution of duplicated genes in Saccharomyces cerevisiae using a new classification method based on protein-protein interaction data

The concept of protein function is widely used and manipulated by biologists. However, the means of the concept and its understanding may vary depending on the level of functionality one considers (molecular, cellular, physiological, etc.). Genomic studies and new high-throughput methods of the post-genomic era provide the opportunity to shed a new light on the concept of protein function: protein-protein interactions can now be considered as pieces of incomplete but still gigantic networks and the analysis of these networks will permit the emergence of a more integrated view of protein function. In this context, we propose a new functional classification method, which, unlike usual methods based on sequence homology, allows the definition of functional classes of protein based on the identity of their interacting partners. An example of such classification will be shown and discussed for a subset of Saccharomyces cerevisiae proteins, accounting for 7% of the yeast proteome. The genome of the budding yeast contains 50% of protein-coding genes that are paralogs, including 457 pairs of duplicated genes coming probably from an ancient whole genome duplication. We will comment on the functional classification of the duplicated genes when using our method and discuss the contribution of these results to the understanding of function evolution for the duplicated genes.     

Long noncoding RNA HOXA-AS2 regulates the expression of SCN3A by sponging miR-106a in breast cancer

Breast cancer is the most commonly diagnosed cancer that affects women worldwide. This study aimed to investigate the competing endogenous RNAs (ceRNAs) mechanism in breast cancer. Microarray data were downloaded from the University of California Santa Cruz (UCSC) Xena database. The limma package was used to screen the differentially expressed messenger RNAs (DEMs) and differentially expressed long noncoding RNAs (DELs). Subsequently, functional analysis was performed using DAVID tool. After constructing the protein-protein interaction (PPI) network, we identified the major gene modules using the Cytoscape software. Univariate survival analysis in the survival package was performed. Finally, the ceRNA regulatory network was constructed to identify the critical genes. A total of 1380 DEMs and 345 DELs were identified in breast cancer samples compared with normal samples. Functional enrichment analysis showed that DEMs were mainly involved in cell division, and cell cycle. We screened four major gene modules and identified the hub nodes in these functional modules. Several DEMs (including FABP7, C4BPA, and LAMB3) and three long noncoding RNAs (lncRNAs) (LINC00092, SLC26A4.AS1, and COLCA1) exhibited significant correlation with patients' survival outcomes. In the ceRNA network, the lncRNA HOXA-AS2 regulated the expression level of SCN3A by interacting with hsa-miR-106a-5p. Thus, our study investigated the ceRNA mechanism in breast cancer. The results showed that lncRNA HOXA-AS2 might modulate the expression of SCN3A by sponging miR-106a in breast cancer.     

牵出技术在蛋白质相互作用研究中的应用

蛋白质与蛋白质的相互作用参与生命体内的许多生物学过程,关于蛋白质相互作用的研究是人们了解蛋白质功能、揭开生命奥秘的关键所在。牵出(pull-down)技术作为一种简单、经济、行之有效的一种体外验证蛋白质与蛋白质之间的相互作用的实验技术,近几年受到了科研人员的青睐。本文阐述了该技术的基本原理和技术特点,总结了近年来牵出技术在生命科学领域中的应用情况,以及由此衍生的新技术的研究进展。     

Prediction of protein-protein interactions using protein signature profiling

Protein domains are conserved and functionally independent structures that play an important role in interactions among related proteins. Domain-domain interactions have been recently used to predict protein-protein interactions (PPI). In general, the interaction probability of a pair of domains is scored using a trained scoring function. Satisfying a threshold, the protein pairs carrying those domains are regarded as "interacting". In this study, the signature contents of proteins were utilized to predict PPI pairs in Saccharomyces cerevisiae, Caenorhabditis elegans, and Homo sapiens. Similarity between protein signature patterns was scored and PPI predictions were drawn based on the binary similarity scoring function. Results show that the true positive rate of prediction by the proposed approach is approximately 32% higher than that using the maximum likelihood estimation method when compared with a test set, resulting in 22% increase in the area under the receiver operating characteristic (ROC) curve. When proteins containing one or two signatures were removed, the sensitivity of the predicted PPI pairs increased significantly. The predicted PPI pairs are on average 11 times more likely to interact than the random selection at a confidence level of 0.95, and on average 4 times better than those predicted by either phylogenetic profiling or gene expression profiling.