Signalling pathway database usability: lessons learned

Testing the Significance of Ranked Gene Sets in Genome-wide Transcriptome

Background

Popular gene set enrichment analysis approaches assumed that genes in the gene set contributed to the statistics equally. However, the genes in the transcription factors (TFs) derived gene sets, or gene sets constructed by TF targets identified by the ChIP-Seq experiment, have a rank attribute, as each of these genes have been assigned with a p-value which indicates the true or false possibilities of the ownerships of the genes belong to the gene sets.

Objectives

Ignoring the rank information during the enrichment analysis will lead to improper statistical inference. We address this issue by developing of new method to test the significance of ranked gene sets in genome-wide transcriptome profiling data.

Methods

A method was proposed by first creating ranked gene sets and gene lists and then applying weighted Kendall's tau rank correlation statistics to the test. After introducing top-down weights to the genes in the gene set, a new software called "Flaver" was developed.

Results

Theoretical properties of the proposed method were established, and its differences over the GSEA approach were demonstrated when analyzing the transcriptome profiling data across 55 human tissues and 176 human cell-lines. The results indicated that the TFs identified by our method have higher tendency to be differentially expressed across the tissues analyzed than its competitors. It significantly outperforms the well-known gene set enrichment analyzing tools, GOStats (9%) and GSEA (17%), in analyzing well-documented human RNA transcriptome datasets.

Conclusions

The method is outstanding in detecting gene sets of which the gene ranks were correlated with the expression levels of the genes in the transcriptome data.

Explainable Drug Repurposing Approach From Biased Random Walks

Drug repurposing is a highly active research area, aiming at finding novel uses for drugs that have been previously developed for other therapeutic purposes. Despite the flourishing of methodologies, success is still partial, and different approaches offer, each, peculiar advantages. In this composite landscape, we present a novel methodology focusing on an efficient mathematical procedure based on gene similarity scores and biased random walks which rely on robust drug-gene-disease association data sets. The recommendation mechanism is further unveiled by means of the Markov chain underlying the random walk process, hence providing explainability about how findings are suggested. Performances evaluation and the analysis of a case study on rheumatoid arthritis show that our approach is accurate in providing useful recommendations and is computationally efficient, compared to the state of the art of drug repurposing approaches.

Systematic assessment of pathway databases, based on a diverse collection of user-submitted experiments

A knowledge-based grouping of genes into pathways or functional units is essential for describing and understanding cellular complexity. However, it is not always clear a priori how and at what level of specificity functionally interconnected genes should be partitioned into pathways, for a given application. Here, we assess and compare nine existing and two conceptually novel functional classification systems, with respect to their discovery power and generality in gene set enrichment testing. We base our assessment on a collection of nearly 2000 functional genomics datasets provided by users of the STRING database. With these real-life and diverse queries, we assess which systems typically provide the most specific and complete enrichment results. We find many structural and performance differences between classification systems. Overall, the well-established, hierarchically organized pathway annotation systems yield the best enrichment performance, despite covering substantial parts of the human genome in general terms only. On the other hand, the more recent unsupervised annotation systems perform strongest in understudied areas and organisms, and in detecting more specific pathways, albeit with less informative labels.

Application of text mining to develop AOP-based mucus hypersecretion genesets and confirmation with in vitro and clinical samples

Mucus hypersecretion contributes to lung function impairment observed in COPD (chronic obstructive pulmonary disease), a tobacco smoking-related disease. A detailed mucus hypersecretion adverse outcome pathway (AOP) has been constructed from literature reviews, experimental and clinical data, mapping key events (KEs) across biological organisational hierarchy leading to an adverse outcome. AOPs can guide the development of biomarkers that are potentially predictive of diseases and support the assessment frameworks of nicotine products including electronic cigarettes. Here, we describe a method employing manual literature curation supported by a focused automated text mining approach to identify genes involved in 5 KEs contributing to decreased lung function observed in tobacco-related COPD. KE genesets were subsequently confirmed by unsupervised clustering against 3 different transcriptomic datasets including (1) in vitro acute cigarette smoke and e-cigarette aerosol exposure, (2) in vitro repeated incubation with IL-13, and (3) lung biopsies from COPD and healthy patients. The 5 KE genesets were demonstrated to be predictive of cigarette smoke exposure and mucus hypersecretion in vitro, and less conclusively predict the COPD status of lung biopsies. In conclusion, using a focused automated text mining and curation approach with experimental and clinical data supports the development of risk assessment strategies utilising AOPs.

The transcriptional response to oxidative stress is part of, but not sufficient for, insulin resistance in adipocytes

Insulin resistance is a major risk factor for metabolic diseases such as Type 2 diabetes. Although the underlying mechanisms of insulin resistance remain elusive, oxidative stress is a unifying driver by which numerous extrinsic signals and cellular stresses trigger insulin resistance. Consequently, we sought to understand the cellular response to oxidative stress and its role in insulin resistance. Using cultured 3T3-L1 adipocytes, we established a model of physiologically-derived oxidative stress by inhibiting the cycling of glutathione and thioredoxin, which induced insulin resistance as measured by impaired insulin-stimulated 2-deoxyglucose uptake. Using time-resolved transcriptomics, we found > 2000 genes differentially-expressed over 24 hours, with specific metabolic and signalling pathways enriched at different times. We explored this coordination using a knowledge-based hierarchical-clustering approach to generate a temporal transcriptional cascade and identify key transcription factors responding to oxidative stress. This response shared many similarities with changes observed in distinct insulin resistance models. However, an anti-oxidant reversed insulin resistance phenotypically but not transcriptionally, implying that the transcriptional response to oxidative stress is insufficient for insulin resistance. This suggests that the primary site by which oxidative stress impairs insulin action occurs post-transcriptionally, warranting a multi-level ‘trans-omic’ approach when studying time-resolved responses to cellular perturbations.

Spatio-Temporal Gene Expression Profiling during In Vivo Early Ovarian Folliculogenesis: Integrated Transcriptomic Study and Molecular Signature of Early Follicular Growth

Background

The successful achievement of early ovarian folliculogenesis is important for fertility and reproductive life span. This complex biological process requires the appropriate expression of numerous genes at each developmental stage, in each follicular compartment. Relatively little is known at present about the molecular mechanisms that drive this process, and most gene expression studies have been performed in rodents and without considering the different follicular compartments.

Results

We used RNA-seq technology to explore the sheep transcriptome during early ovarian follicular development in the two main compartments: oocytes and granulosa cells. We documented the differential expression of 3,015 genes during this phase and described the gene expression dynamic specific to these compartments. We showed that important steps occurred during primary/secondary transition in sheep. We also described the in vivo molecular course of a number of pathways. In oocytes, these pathways documented the chronology of the acquisition of meiotic competence, migration and cellular organization, while in granulosa cells they concerned adhesion, the formation of cytoplasmic projections and steroid synthesis. This study proposes the involvement in this process of several members of the integrin and BMP families. The expression of genes such as Kruppel-like factor 9 (KLF9) and BMP binding endothelial regulator (BMPER) was highlighted for the first time during early follicular development, and their proteins were also predicted to be involved in gene regulation. Finally, we selected a data set of 24 biomarkers that enabled the discrimination of early follicular stages and thus offer a molecular signature of early follicular growth. This set of biomarkers includes known genes such as SPO11 meiotic protein covalently bound to DSB (SPO11), bone morphogenetic protein 15 (BMP15) and WEE1 homolog 2 (S. pombe)(WEE2) which play critical roles in follicular development but other biomarkers are also likely to play significant roles in this process.

Conclusions

To our knowledge, this is the first in vivo spatio-temporal exploration of transcriptomes derived from early follicles in sheep.

Interlog protein network: an evolutionary benchmark of protein interaction networks for the evaluation of clustering algorithms

BACKGROUND

In the field of network science, exploring principal and crucial modules or communities is critical in the deduction of relationships and organization of complex networks. This approach expands an arena, and thus allows further study of biological functions in the field of network biology. As the clustering algorithms that are currently employed in finding modules have innate uncertainties, external and internal validations are necessary.

METHODS

Sequence and network structure alignment, has been used to define the Interlog Protein Network (IPN). This network is an evolutionarily conserved network with communal nodes and less false-positive links. In the current study, the IPN is employed as an evolution-based benchmark in the validation of the module finding methods. The clustering results of five algorithms; Markov Clustering (MCL), Restricted Neighborhood Search Clustering (RNSC), Cartographic Representation (CR), Laplacian Dynamics (LD) and Genetic Algorithm; to find communities in Protein-Protein Interaction networks (GAPPI) are assessed by IPN in four distinct Protein-Protein Interaction Networks (PPINs).

RESULTS

The MCL shows a more accurate algorithm based on this evolutionary benchmarking approach. Also, the biological relevance of proteins in the IPN modules generated by MCL is compatible with biological standard databases such as Gene Ontology, KEGG and Reactome.

CONCLUSION

In this study, the IPN shows its potential for validation of clustering algorithms due to its biological logic and straightforward implementation.

Comparison of human cell signaling pathway databases--evolution, drawbacks and challenges

Elucidating the complexities of cell signaling pathways is of immense importance to gain understanding about various biological phenomenon, such as dynamics of gene/protein expression regulation, cell fate determination, embryogenesis and disease progression. The successful completion of human genome project has also helped experimental and theoretical biologists to analyze various important pathways. To advance this study, during the past two decades, systematic collections of pathway data from experimental studies have been compiled and distributed freely by several databases, which also integrate various computational tools for further analysis. Despite significant advancements, there exist several drawbacks and challenges, such as pathway data heterogeneity, annotation, regular update and automated image reconstructions, which motivated us to perform a thorough review on popular and actively functioning 24 cell signaling databases. Based on two major characteristics, pathway information and technical details, freely accessible data from commercial and academic databases are examined to understand their evolution and enrichment. This review not only helps to identify some novel and useful features, which are not yet included in any of the databases but also highlights their current limitations and subsequently propose the reasonable solutions for future database development, which could be useful to the whole scientific community.

Multi-omic landscape of rheumatoid arthritis: re-evaluation of drug adverse effects

Objective: To provide a frame to estimate the systemic impact (side/adverse events) of (novel) therapeutic targets by taking into consideration drugs potential on the numerous districts involved in rheumatoid arthritis (RA) from the inflammatory and immune response to the gut-intestinal (GI) microbiome.

Methods: We curated the collection of molecules from high-throughput screens of diverse (multi-omic) biochemical origin, experimentally associated to RA. Starting from such collection we generated RA-related protein-protein interaction (PPI) networks (interactomes) based on experimental PPI data. Pharmacological treatment simulation, topological and functional analyses were further run to gain insight into the proteins most affected by therapy and by multi-omic modeling.

Results: Simulation on the administration of MTX results in the activation of expected (apoptosis) and adverse (nitrogenous metabolism alteration) effects. Growth factor receptor-bound protein 2 (GRB2) and Interleukin-1 Receptor Associated Kinase-4 (IRAK4, already an RA target) emerge as relevant nodes. The former controls the activation of inflammatory, proliferative and degenerative pathways in host and pathogens. The latter controls immune alterations and blocks innate response to pathogens.

Conclusions: This multi-omic map properly recollects in a single analytical picture known, yet complex, information like the adverse/side effects of MTX, and provides a reliable platform for in silico hypothesis testing or recommendation on novel therapies. These results can support the development of RA translational research in the design of validation experiments and clinical trials, as such we identify GRB2 as a robust potential new target for RA for its ability to control both synovial degeneracy and dysbiosis, and, conversely, warn on the usage of IRAK4-inhibitors recently promoted, as this involves potential adverse effects in the form of impaired innate response to pathogens.