Large-scale data-driven integrative framework for extracting essential targets and processes from disease-associated gene data sets

Network-driven analysis of human-Plasmodium falciparum interactome: processes for malaria drug discovery and extracting in silico targets

BACKGROUND

The emergence and spread of malaria drug resistance have resulted in the need to understand disease mechanisms and importantly identify essential targets and potential drug candidates. Malaria infection involves the complex interaction between the host and pathogen, thus, functional interactions between human and Plasmodium falciparum is essential to obtain a holistic view of the genetic architecture of malaria. Several functional interaction studies have extended the understanding of malaria disease and integrating such datasets would provide further insights towards understanding drug resistance and/or genetic resistance/susceptibility, disease pathogenesis, and drug discovery.

METHODS

This study curated and analysed data including pathogen and host selective genes, host and pathogen protein sequence data, protein-protein interaction datasets, and drug data from literature and databases to perform human-host and P. falciparum network-based analysis. An integrative computational framework is presented that was developed and found to be reasonably accurate based on various evaluations, applications, and experimental evidence of outputs produced, from data-driven analysis.

RESULTS

This approach revealed 8 hub protein targets essential for parasite and human host-directed malaria drug therapy. In a semantic similarity approach, 26 potential repurposable drugs involved in regulating host immune response to inflammatory-driven disorders and/or inhibiting residual malaria infection that can be appropriated for malaria treatment. Further analysis of host-pathogen network shortest paths enabled the prediction of immune-related biological processes and pathways subverted by P. falciparum to increase its within-host survival.

CONCLUSIONS

Host-pathogen network analysis reveals potential drug targets and biological processes and pathways subverted by P. falciparum to enhance its within malaria host survival. The results presented have implications for drug discovery and will inform experimental studies.

Investigations of Kidney Dysfunction-Related Gene Variants in Sickle Cell Disease Patients in Cameroon (Sub-Saharan Africa)

Background

Renal dysfunctions are associated with increased morbidity and mortality in sickle cell disease (SCD). Early detection and subsequent management of SCD patients at risk for renal failure and dysfunctions are essential, however, predictors that can identify patients at risk of developing renal dysfunction are not fully understood.

Methods

In this study, we have investigated the association of 31 known kidney dysfunctions-related variants detected in African Americans from multi-ethnic genome wide studies (GWAS) meta-analysis, to kidney-dysfunctions in a group of 413 Cameroonian patients with SCD. Systems level bioinformatics analyses were performed, employing protein-protein interaction networks to further interrogate the putative associations.

Results

Up to 61% of these patients had micro-albuminuria, 2.4% proteinuria, 71% glomerular hyperfiltration, and 5.9% had renal failure. Six variants are significantly associated with the two quantifiable phenotypes of kidney dysfunction (eGFR and crude-albuminuria): A1CF-rs10994860 (P = 0.02020), SYPL2-rs12136063 (P = 0.04208), and APOL1 (G1)-rs73885319 (P = 0.04610) are associated with eGFR; and WNT7A-rs6795744 (P = 0.03730), TMEM60-rs6465825 (P = 0.02340), and APOL1 (G2)-rs71785313 (P = 0.03803) observed to be protective against micro-albuminuria. We identified a protein-protein interaction sub-network containing three of these gene variants: APOL1, SYPL2, and WNT7A, connected to the Nuclear factor NF-kappa-B p105 subunit (NFKB1), revealed to be essential and might indirectly influence extreme phenotypes. Interestingly, clinical variables, including body mass index (BMI), systolic blood pressure, vaso-occlusive crisis (VOC), and haemoglobin (Hb), explain better the kidney phenotypic variations in this SCD population.

Conclusion

This study highlights a strong contribution of haematological indices (Hb level), anthropometric variables (BMI, blood pressure), and clinical events (i.e., vaso-occlusive crisis) to kidney dysfunctions in SCD, rather than known genetic factors. Only 6/31 characterised gene-variants are associated with kidney dysfunction phenotypes in SCD samples from Cameroon. The data reveal and emphasise the urgent need to extend GWAS studies in populations of African ancestries living in Africa, and particularly for kidney dysfunctions in SCD.

IHP-PING-generating integrated human protein-protein interaction networks on-the-fly

Advances in high-throughput sequencing technologies have resulted in an exponential growth of publicly accessible biological datasets. In the 'big data' driven 'post-genomic' context, much work is being done to explore human protein-protein interactions (PPIs) for a systems level based analysis to uncover useful signals and gain more insights to advance current knowledge and answer specific biological and health questions. These PPIs are experimentally or computationally predicted, stored in different online databases and some of PPI resources are updated regularly. As with many biological datasets, such regular updates continuously render older PPI datasets potentially outdated. Moreover, while many of these interactions are shared between these online resources, each resource includes its own identified PPIs and none of these databases exhaustively contains all existing human PPI maps. In this context, it is essential to enable the integration of or combining interaction datasets from different resources, to generate a PPI map with increased coverage and confidence. To allow researchers to produce an integrated human PPI datasets in real-time, we introduce the integrated human protein-protein interaction network generator (IHP-PING) tool. IHP-PING is a flexible python package which generates a human PPI network from freely available online resources. This tool extracts and integrates heterogeneous PPI datasets to generate a unified PPI network, which is stored locally for further applications.

Computational/in silico methods in drug target and lead prediction

Drug-like compounds are most of the time denied approval and use owing to the unexpected clinical side effects and cross-reactivity observed during clinical trials. These unexpected outcomes resulting in significant increase in attrition rate centralizes on the selected drug targets. These targets may be disease candidate proteins or genes, biological pathways, disease-associated microRNAs, disease-related biomarkers, abnormal molecular phenotypes, crucial nodes of biological network or molecular functions. This is generally linked to several factors, including incomplete knowledge on the drug targets and unpredicted pharmacokinetic expressions upon target interaction or off-target effects. A method used to identify targets, especially for polygenic diseases, is essential and constitutes a major bottleneck in drug development with the fundamental stage being the identification and validation of drug targets of interest for further downstream processes. Thus, various computational methods have been developed to complement experimental approaches in drug discovery. Here, we present an overview of various computational methods and tools applied in predicting or validating drug targets and drug-like molecules. We provide an overview on their advantages and compare these methods to identify effective methods which likely lead to optimal results. We also explore major sources of drug failure considering the challenges and opportunities involved. This review might guide researchers on selecting the most efficient approach or technique during the computational drug discovery process.

An NMF-Based Methodology for Selecting Biomarkers in the Landscape of Genes of Heterogeneous Cancer-Associated Fibroblast Populations

The rapid development of high-performance technologies has greatly promoted studies of molecular oncology producing large amounts of data. Even if these data are publicly available, they need to be processed and studied to extract information useful to better understand mechanisms of pathogenesis of complex diseases, such as tumors. In this article, we illustrated a procedure for mining biologically meaningful biomarkers from microarray datasets of different tumor histotypes. The proposed methodology allows to automatically identify a subset of potentially informative genes from microarray data matrices, which differs either in the number of rows (genes) and of columns (patients). The methodology integrates nonnegative matrix factorization method, a functional enrichment analysis web tool with a properly designed gene extraction procedure to allow the analysis of omics input data with different row size. The proposed methodology has been used to mine microarray of solid tumors of different embryonic origin to verify the presence of common genes characterizing the heterogeneity of cancer-associated fibroblasts. These automatically extracted biomarkers could be used to suggest appropriate therapies to inactivate the state of active fibroblasts, thus avoiding their action on tumor progression.

Common transcriptional programs and the role of chemokine (C-C motif) ligand 20 (CCL20) in cell migration of cholangiocarcinoma

The incidence of cholangiocarcinoma (CCA) has risen in many countries, but there is still no appropriate screening and treatment available. The growing number of microarray data published todays can be a powerful resource for the discovery of biomarkers to tackle challenges in the management of CCA. This study analyzed multiple microarray datasets to identify the common transcriptional networks in CCA and select a possible biomarker for functional study in CCA cell lines. A systematic searching identified 4 microarray datasets from Gene Expression Omnibus (GEO) repository and PubMed articles. Differential expression analysis between tumor and normal tissues was performed in each dataset. In order to characterize the common expression pattern, differentially expressed genes (DEGs) from all datasets were combined and visualized by hierarchical clustering and heatmap. Gene enrichment analysis performed in each cluster revealed that over-expressed DEGs were enriched in cell cycle, cell migration and response to cytokines while under-expressed DEGs were enriched in metabolic processes such as oxidation-reduction, lipid, and drug. To explain tumor characteristics, genes enriched in cell migration and response to cytokines were further investigated. Among these genes, CCL20 was selected for functional study because its role has never been studied in CCA. Moreover, its signaling may be regulated by disrupting its only receptor, CCR6. Treatment with recombinant CCL20 induced higher cell migration and increased expression of N-cad. In contrast, knockdown of CCR6 by siRNA reduced cell migration ability and decreased N-cadherin level. Altogether, these results suggested the contribution of CCL20/CCR6 signaling in cell migration through epithelial-mesenchymal transition process. Thus, CCL20/CCR6 signaling might be a target for the management of CCA.