Zhiqiang Chang

Prediction of human protein–protein interaction by a mixed Bayesian model and its application to exploring underlying cancer-related pathway crosstalk

Protein–protein interaction (PPI) prediction method has provided an opportunity for elucidating potential biological processes and disease mechanisms. We integrated eight features involving proteomic, genomic, phenotype and functional annotation datasets by a mixed model consisting of full connected Bayesian (FCB) model and naive Bayesian model to predict human PPIs, resulting in 40 447 PPIs which contain 2740 common PPIs with the human protein reference database (HPRD) by a likelihood ratio cutoff of 512. Then we applied them to exploring underlying pathway crosstalk where pathways were derived from the pathway interaction database. Two pathway crosstalk networks (PCNs) were constructed based on PPI sets. The PPI sets were derived from two different sources. One source was strictly the HPRD database while the other source was a combination of HPRD and PPIs predicted by our mixed Bayesian method. We demonstrated that PCNs based on the mixed PPI set showed much more underlying pathway interactions than the HPRD PPI set. Furthermore, we mapped cancer-causing mutated somatic genes to PPIs between significant pathway crosstalk pairs. We extracted highly connected clusters from over-represented subnetworks of PCNs, which were enriched for mutated gene interactions that acted as crosstalk links. Most of the pathways in top ranking clusters were shown to play important roles in cancer. The clusters themselves showed coherent function categories pertaining to cancer development.

The application of gene co-expression network reconstruction based on CNVs and gene expression microarray data in breast cancer

Copy number variations (CNVs) are one type of the human genetic variations and are pervasive in the human genome. It has been confirmed that they can play a causal role in complex diseases. Previous studies of CNVs focused more on identifying the disease-specific CNV regions or candidate genes on these CNV regions, but less on the synergistic actions between genes on CNV regions and other genes. Our research combined the CNVs with related gene co-expression to reconstruct gene co-expression network by using single nucleotide polymorphism microarray datasets and gene microarray datasets of breast cancer, and then extracted the modules which connected densely inside and analyzed the functions of modules. Interestingly, all of these modules’ functions were related to breast cancer according to our enrichment analysis, and most of the genes in these modules have been reported to be involved in breast cancer. Our findings suggested that integrating CNVs and gene co-expressed relations was an available way to analyze the roles of CNV genes and their synergistic genes in breast cancer, and provided a novel insight into the pathological mechanism of breast cancer.

Functional analysis of microRNA and transcription factor synergistic regulatory network based on identifying regulatory motifs in non-small cell lung cancer

BackgroundLung cancer, especially non-small cell lung cancer, is a leading cause of malignant tumor death worldwide. Understanding the mechanisms employed by the main regulators, such as microRNAs (miRNAs) and transcription factors (TFs), still remains elusive. The patterns of their cooperation and biological functions in the synergistic regulatory network have rarely been studied.ResultsHere, we describe the first miRNA-TF synergistic regulation network in human lung cancer. We identified important regulators (MYC, NFKB1, miR-590, and miR-570) and significant miRNA-TF synergistic regulatory motifs by random simulations. The two most significant motifs were the co-regulation of miRNAs and TFs, and TF-mediated cascade regulation. We also developed an algorithm to uncover the biological functions of the human lung cancer miRNA-TF synergistic regulatory network (regulation of apoptosis, cellular protein metabolic process, and cell cycle), and the specific functions of each miRNA-TF synergistic subnetwork. We found that the miR-17 family exerted important effects in the regulation of non-small cell lung cancer, such as in proliferation and cell cycle regulation by targeting the retinoblastoma protein (RB1) and forming a feed forward loop with the E2F1 TF. We proposed a model for the miR-17 family, E2F1, and RB1 to demonstrate their potential roles in the occurrence and development of non-small cell lung cancer.ConclusionsThis work will provide a framework for constructing miRNA-TF synergistic regulatory networks, function analysis in diseases, and identification of the main regulators and regulatory motifs, which will be useful for understanding the putative regulatory motifs involving miRNAs and TFs, and for predicting new targets for cancer studies.

CNVD: text mining-based copy number variation in disease database.

Copy number variation (CNV) is a kind of chromosomal structural reorganization that has been detected, in this decade, mainly by high‐throughput biological technology. Researchers have found that CNVs are ubiquitous in many species and accumulating evidence indicates that CNVs are closely related with complex diseases. The investigation of chromosomal structural alterations has begun to reveal some important clues to the pathologic causes of diseases and to the disease process. However, many of the published studies have focused on a single disease and, so far, the experimental results have not been systematically collected or organized. Manual text mining from 6301 published papers was used to build the Copy Number Variation in Disease database (CNVD). CNVD contains CNV information for 792 diseases in 22 species from diverse types of experiments, thus, ensuring high confidence and comprehensive representation of the relationship between the CNVs and the diseases. In addition, multiple query modes and visualized results are provided in the CNVD database. With its user‐friendly interface and the integrated CNV information for different diseases, CNVD will offer a truly comprehensive platform for disease research based on chromosomal structural variations. The CNVD interface is accessible at http://bioinfo.hrbmu.edu.cn/CNVD.

Optimization of cell lines as tumour models by integrating multi-omics data

Cell lines are widely used as in vitro models of tumorigenesis. However, an increasing number of researchers have found that cell lines differ from their sourced tumour samples after long-term cell culture. The application of unsuitable cell lines in experiments will affect the experimental accuracy and the treatment of patients. Therefore, it is imperative to identify optimal cell lines for each cancer type. Here, we review the methods used to evaluate cell lines since 2005. Furthermore, gene expression, copy number and mutation profiles from The Cancer Genome Atlas and the Cancer Cell Line Encyclopedia are used to calculate similarity between tumours and cell lines. Then, the ideal cell lines to use for experiments for eight types of cancers are found by combining the results with Gene Ontology functional similarity. After verification, the optimal cell lines have the same genomic characteristics as their homologous tumour samples. The contaminated cell lines identified in previous research are also determined to be unsuitable in vitro cancer models here. Moreover, our study suggests that some of the commonly used cell lines are not suitable cancer models. In summary, we provide a reference for ideal cell lines to use in in vitro experiments and contribute to improving the accuracy of future cancer research. Furthermore, this research provides a foundation for identifying more effective treatment strategies.

Calculating phenotypic similarity between genes using hierarchical structure data based on semantic similarity.

Phenotypic similarity is correlated with a number of measures of gene function, such as relatedness at the level of direct protein-protein interaction. The phenotypic effect of a deleted or mutated gene, which is one part of gene annotation, has caught broad attention. However, there have been few measures to study phenotypic similarity with the data from Human Phenotype Ontology (HPO) database, therefore more analogous measures should be developed and investigated. We used five semantic similarity-based measures (Jiang and Conrath, Lin, Schlicker, Yu and Wu) to calculate the human phenotypic similarity between genes (PSG) with data from HPO database, and evaluated their accuracy with information of protein-protein interaction, protein complex, protein family, gene function or DNA sequence. Compared with the gene pairs that were random selected, the results of these methods were statistically significant (all P<0.001). Furthermore, we assessed the performance of these five measures by receiver operating characteristic (ROC) curve analysis, and found that most of them performed better than the previous methods. This work had proved that these measures based on semantic similarity for calculation of PSG were effective for hierarchical structure data. Our study contributes to the development and optimization of novel algorithms of PSG calculation and provides more alternative methods to researchers as well as tools and directions for PSG study.

Analysis of functional and pathway association of differential co-expressed genes: A case study in drug addiction

Drug addiction has been considered as a kind of chronic relapsing brain disease influenced by both genetic and environmental factors. At present, many causative genes and pathways related to diverse kinds of drug addiction have been discovered, while less attention has been paid to common mechanisms shared by different drugs underlying addiction. By applying a co-expression meta-analysis method to mRNA expression profiles of alcohol, cocaine, heroin addicted and normal samples, we identified significant gene co-expression pairs. As co-expression networks of drug group and control group constructed, associated function term pairs and pathway pairs reflected by co-expression pattern changes were discovered by integrating functional and pathway information respectively. The results indicated that respiratory electron transport chain, synaptic transmission, mitochondrial electron transport, signal transduction, locomotory behavior, response to amphetamine, negative regulation of cell migration, glucose regulation of insulin secretion, signaling by NGF, diabetes pathways, integration of energy metabolism, dopamine receptors may play an important role in drug addiction. In addition, the results can provide theory support for studies of addiction mechanisms.

Identification of cancer‐related potential biomarkers based on lncRNA–pseudogene–mRNA competitive networks

Accumulating evidence indicates that mRNAs and noncoding RNAs act as competitive endogenous RNAs (ceRNAs) and play a key role in tumorigenesis. However, the complex competitive relationship among genes remains unknown. In the present study, the long noncoding RNAs (lncRNAs), pseudogenes and mRNAs that compete with common microRNAs are defined as lncRNA–pseudogene–mRNA competitive triples. We find that some candidate ceRNAs, modules and triples are associated with cancers and can significantly divide patients into high‐risk and low‐risk groups; thus, they may serve as potential cancer biomarkers. In sum, the present study systematically analyzes the association between competitive triples and cancer, which provides a reference for a deeper understanding of cancer progression.

Identification of potential cancer-related pseudogenes in lung adenocarcinoma based on ceRNA hypothesis

Pseudogenes are initially regarded as non-functional genomic fossils resulted from inactivating gene mutations during evolution. Far from being silent, pseudogenes are proved to regulate the expression of protein-coding genes through function as microRNA sponge in vivo. The aim of our study was to propose an integrative systems biology approach to identify disease pseudogenes base on competitive endogenous RNA (ceRNA) hypothesis. Here, we applied our method to lung adenocarcinoma (LUAD) RNASeq data from TCGA and identified 33 candidate pseudogenes. We described the characteristics of the candidate pseudogenes and performed functional enrichment. Through analyzing neighboring genes we found these pseudogenes were surrounded by tumor genes and may involve in tumor pathway. Furthermore, the DNA methylation analysis indicated that 21 pseudogenes co-methylated with their competitive mRNAs. In the co-methylated network, we discovered 6 differentially expressed pseudogenes, which we termed potential LUAD-associated pseudogenes. We further revealed that the 3 ceRNA triples (miR-21-5p-NKAPP1-PRDM11, miR-29c-3p-MSTO2P-EZH2 and miR-29c-3p-RPLP0P2-EZH2), whose high risk groups were associated with the poor prognosis of LUAD, may be considered as potential prognostic signatures. Moreover, by integrating target information of microRNA we also provided a new perspective for the discovery of potential small molecule drugs. This work may facilitate cancer research and serve as the basis for future efforts to understand the role of pseudogenes, develop novel biomarkers and improve knowledge of tumor biology.

Identification of Biomarker and Co-Regulatory Motifs in Lung Adenocarcinoma Based on Differential Interactions

Changes in intermolecular interactions (differential interactions) may influence the progression of cancer. Specific genes and their regulatory networks may be more closely associated with cancer when taking their transcriptional and post-transcriptional levels and dynamic and static interactions into account simultaneously. In this paper, a differential interaction analysis was performed to detect lung adenocarcinoma-related genes. Furthermore, a miRNA-TF (transcription factor) synergistic regulation network was constructed to identify three kinds of co-regulated motifs, namely, triplet, crosstalk and joint. Not only were the known cancer-related miRNAs and TFs (let-7, miR-15a, miR-17, TP53, ETS1, and so on) were detected in the motifs, but also the miR-15, let-7 and miR-17 families showed a tendency to regulate the triplet, crosstalk and joint motifs, respectively. Moreover, several biological functions (i.e., cell cycle, signaling pathways and hemopoiesis) associated with the three motifs were found to be frequently targeted by the drugs for lung adenocarcinoma. Specifically, the two 4-node motifs (crosstalk and joint) based on co-expression and interaction had a closer relationship to lung adenocarcinoma, and so further research was performed on them. A 10-gene biomarker (UBC, SRC, SP1, MYC, STAT3, JUN, NR3C1, RB1, GRB2 and MAPK1) was selected from the joint motif, and a survival analysis indicated its significant association with survival. Among the ten genes, JUN, NR3C1 and GRB2 are our newly detected candidate lung adenocarcinoma-related genes. The genes, regulators and regulatory motifs detected in this work will provide potential drug targets and new strategies for individual therapy.