Identifying Novel Candidate Genes Related to Apoptosis from a Protein-Protein Interaction Network

Integrated bioinformatics analysis of core regulatory elements involved in keloid formation

BACKGROUND

Keloid is a benign fibro-proliferative dermal tumor formed by an abnormal scarring response to injury and characterized by excessive collagen accumulation and invasive growth. The mechanism of keloid formation has not been fully elucidated, especially during abnormal scarring. Here, we investigated the regulatory genes, micro-RNAs (miRNAs) and transcription factors (TFs) that influence keloid development by comparing keloid and normal scar as well as keloid and normal skin.

METHODS

Gene expression profiles (GSE7890, GSE92566, GSE44270 and GSE3189) of 5 normal scar samples, 10 normal skin samples and 18 keloid samples from the Gene Expression Omnibus (GEO) database were interrogated. Differentially expressed genes (DEGs) were identified between keloid and normal skin samples as well as keloid and normal scar samples with R Project for Statistical Computing. Gene Ontology (GO) functional enrichment analysis was also performed with R software. DEG-associated protein-protein interaction (PPI) network was constructed by STRING, followed by module selection from the PPI network based on the MCODE analysis. Regulatory relationships between TF/miRNA and target genes were predicted with miRnet and cytoscape. Core regulatory genes were verified by RT-qPCR.

RESULTS

We identified 628 DEGs, of which 626 were up-regulated and 2 were down-regulated. Seven core genes [neuropeptide Y(NPY), 5-hydroxytryptamine receptor 1A(HTR1A), somatostatin (SST), adenylate cyclase 8 (ADCY8), neuromedin U receptor 1 (NMUR1), G protein subunit gamma 3 (GNG3), and G protein subunit gamma 13 (GNG13)] all belong to MCODE1 and were enriched in the "G protein coupled receptor signaling pathway" of the GO biological process category. Furthermore, nine core miRNAs (hsa-mir-124, hsa-let-7, hsa-mir-155, hsa-mir-26a, hsa-mir-941, hsa-mir-10b, hsa-mir-20, hsa-mir-31 and hsa-mir-372), and two core TFs (SP1 and TERT) were identified to play important roles in keloid formation. In the TF/miRNA-target gene network, both hsa-mir-372 and hsa-mir-20 had a regulatory effect on GNG13, ADCY8 was predicted to be target by hsa-mir-10b, and HTR1A and NPY were potentially by SP1. Furthermore, the expression of core regulatory genes (GNG13, ADCY8, HTR1A and NPY) was validated in clinical samples.

CONCLUSIONS

GNG13, ADCY8, NPY and HTR1A may act as core genes in keloid formation and these core genes establish relationship with SP1 and miRNA (hsa-mir-372, hsa-mir-20, hsa-mir-10b), which may influence multiple signaling pathways in the pathogenesis of keloid.

Identification of Therapeutic Targets for the Selective Killing of HBV-Positive Hepatocytes

The hepatitis B virus (HBV) infection is a major risk factor for cirrhosis and hepatocellular carcinoma. Most infected individuals become lifelong carriers of HBV as the drugs currently used to treat the patients can only control the disease, thereby achieving functional cure (loss of the hepatitis B surface antigen) but not complete cure (elimination of infected hepatocytes). Therefore, we aimed to identify the target genes for the selective killing of HBV-positive hepatocytes to develop a novel therapy for the treatment of HBV infection. Our strategy was to recognize the conditionally essential genes that are essential for the survival of HBV-positive hepatocytes, but non-essential for the HBV-negative hepatocytes. Using microarray gene expression data curated from the Gene Expression Omnibus database and the known essential genes from the Online GEne Essentiality database, we used two approaches, comprising the random walk with restart algorithm and the support vector machine approach, to determine the potential targets for the selective killing of HBV-positive hepatocytes. The final candidate genes list obtained using these two approaches consisted of 36 target genes, which may be conditionally essential for the cell survival of HBV-positive hepatocytes; however, this requires further experimental validation. Therefore, the genes identified in this study can be used as potential drug targets to develop novel therapeutic strategies for the treatment of HBV, and may ultimately help in achieving the elusive goal of a complete cure for hepatitis B.

New Approach for Risk Estimation Algorithms of BRCA1/2 Negativeness Detection with Modelling Supervised Machine Learning Techniques

BRCA1/2 gene testing is a difficult, expensive, and time-consuming test which requires excessive work load. The identification of the BRCA1/2 gene mutations is significantly important in the selection of treatment and the risk of secondary cancer. We aimed to develop an algorithm considering all the clinical, demographic, and genetic features of patients for identifying the BRCA1/2 negativity in the present study. An experimental dataset was created with the collection of the all clinical, demographic, and genetic features of breast cancer patients for 20 years. This dataset consisted of 125 features of 2070 high-risk breast cancer patients. All data were numeralized and normalized for detection of the BRCA1/2 negativity in the machine learning algorithm. The performance of the algorithm was identified by studying the machine learning model with the test data. k nearest neighbours (KNN) and decision tree (DT) accuracy rates of 9 features involving Dataset 2 were found to be the most effective. The removal of the unnecessary data in the dataset by reducing the number of features was shown to increase the accuracy rate of algorithm compared with the DT. BRCA1/2 negativity was identified without performing the BRCA1/2 gene test with 92.88% accuracy within minutes in high-risk breast cancer patients with this algorithm, and the test associated result waiting stress, time, and money loss were prevented. That algorithm is suggested be useful in fast performing of the treatment plans of patients and accurately in addition to speeding up the clinical practice.

Multiomic analysis of oral keratinocytes chronically exposed to shisha

BACKGROUND

Tobacco is smoked in different form including cigarettes and water pipes. One popular form of water pipe smoking especially in Middle Eastern countries is shisha smoking. Shisha has been associated with various diseases including oral cancer. However, genomic alterations and gene expression changes associated with chronic shisha exposure have not been previously investigated.

OBJECTIVES

Whole-exome sequencing and gene expression profiling of immortalized human oral keratinocytes (OKF6/TERT1) cells chronically treated with 0.5% shisha extract for a period of 8 months was undertaken to characterize molecular alterations associated with shisha exposure.

METHODS

Genomic DNA and RNA were extracted and preprocessed as per manufacturer's instruction and subjected to whole-exome and transcriptome sequencing using Illumina HiSeq2500 platform. Exome was analyzed using GATK pipeline whereas RNA-Seq data was analyzed using HiSat2 and HTSeq along with DESeq to elucidate differentially expressed genes.

RESULTS

Whole-exome sequence analysis led to identification of 521 somatic missense variants corresponding to 389 genes RNA-Seq data revealed 247 differentially expressed genes (≥2-fold, P-value<0.01) in shisha treated cells compared to parental cells. Pathway analysis of differentially expressed genes revealed that interferon-signaling pathway was significantly affected. We predict activation of MAPK1 pathway which is known to play a key role in oral cancer. We also observed allele specific expression of mutant LIMA1 based on RNA-Seq dataset.

CONCLUSION

Our findings provide insights into genomic alterations and gene expression pattern associated with oral keratinocytes chronically exposed to shisha.

Network-based method for mining novel HPV infection related genes using random walk with restart algorithm

The human papillomavirus (HPV), a common virus that infects the reproductive tract, may lead to malignant changes within the infection area in certain cases and is directly associated with such cancers as cervical cancer, anal cancer, and vaginal cancer. Identification of novel HPV infection related genes can lead to a better understanding of the specific signal pathways and cellular processes related to HPV infection, providing information for the development of more efficient therapies. In this study, several novel HPV infection related genes were predicted by a computation method based on the known genes involved in HPV infection from HPVbase. This method applied the algorithm of random walk with restart (RWR) to a protein-protein interaction (PPI) network. The candidate genes were further filtered by the permutation and association tests. These steps eliminated genes occupying special positions in the PPI network and selected key genes with strong associations to known HPV infection related genes based on the interaction confidence and functional similarity obtained from published databases, such as STRING, gene ontology (GO) terms and KEGG pathways. Our study identified 104 novel HPV infection related genes, a number of which were confirmed to relate to the infection processes and complications of HPV infection, as reported in the literature. These results demonstrate the reliability of our method in identifying HPV infection related genes. This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.

A network-based method using a random walk with restart algorithm and screening tests to identify novel genes associated with Menière's disease

As a chronic illness derived from hair cells of the inner ear, Menière’s disease (MD) negatively influences the quality of life of individuals and leads to a number of symptoms, such as dizziness, temporary hearing loss, and tinnitus. The complete identification of novel genes related to MD would help elucidate its underlying pathological mechanisms and improve its diagnosis and treatment. In this study, a network-based method was developed to identify novel MD-related genes based on known MD-related genes. A human protein-protein interaction (PPI) network was constructed using the PPI information reported in the STRING database. A classic ranking algorithm, the random walk with restart (RWR) algorithm, was employed to search for novel genes using known genes as seed nodes. To make the identified genes more reliable, a series of screening tests, including a permutation test, an interaction test and an enrichment test, were designed to select essential genes from those obtained by the RWR algorithm. As a result, several inferred genes, such as CD4, NOTCH2 and IL6, were discovered. Finally, a detailed biological analysis was performed on fifteen of the important inferred genes, which indicated their strong associations with MD.

Identifying novel fruit-related genes in Arabidopsis thaliana based on the random walk with restart algorithm

Fruit is essential for plant reproduction and is responsible for protection and dispersal of seeds. The development and maturation of fruit is tightly regulated by numerous genetic factors that respond to environmental and internal stimulation. In this study, we attempted to identify novel fruit-related genes in a model organism, Arabidopsis thaliana, using a computational method. Based on validated fruit-related genes, the random walk with restart (RWR) algorithm was applied on a protein-protein interaction (PPI) network using these genes as seeds. The identified genes with high probabilities were filtered by the permutation test and linkage tests. In the permutation test, the genes that were selected due to the structure of the PPI network were discarded. In the linkage tests, the importance of each candidate gene was measured from two aspects: (1) its functional associations with validated genes and (2) its similarity with validated genes on gene ontology (GO) terms and KEGG pathways. Finally, 255 inferred genes were obtained, subsequent extensive analysis of important genes revealed that they mainly contribute to ubiquitination (UBQ9, UBQ8, UBQ11, UBQ10), serine hydroxymethyl transfer (SHM7, SHM5, SHM6) or glycol-metabolism (HXKL2_ARATH, CSY5, GAPCP1), suggesting essential roles during the development and maturation of fruit in Arabidopsis thaliana.