Systematic analysis and prediction of longevity genes in Caenorhabditis elegans

Unique Patterns in Amino Acid Sequences of Aging-Related Proteins

Aging has strong genetic components and the list of genes that may regulate the aging process is collected in the GenAge database. There may be characteristic patterns in the amino acid sequences of aging-related proteins that distinguish them from other proteins and this information will lead to a better understanding of the aging process. To test this hypothesis, human protein sequences are extracted from the UniProt database and the relative frequency of every amino acid residue in aging-related proteins and the remaining proteins is calculated. The main observation is that the mean relative frequency of aspartic acid (D) is consistently higher, while the mean relative frequencies of tryptophan (W) and leucine (L) are consistently lower in aging-related proteins compared to the non-aging-related proteins for the human and four examined model organisms. It is also observed that the mean relative frequency of aspartic acid is higher, while the mean relative frequency of tryptophan is lower in pro-longevity proteins compared to anti-longevity proteins in model organisms. Finally, it is found that aging-related proteins tend to be longer than non-aging-related proteins. It is hoped that this analysis initiates further computational and experimental research to explore the underlying mechanisms of these findings.

A meta-analysis of RNA-Seq studies to identify novel genes that regulate aging

Aging is a ubiquitous biological process that limits the maximal lifespan of most organisms. Significant efforts by many groups have identified mechanisms that, when triggered by natural or artificial stimuli, are sufficient to either enhance or decrease maximal lifespan. Previous aging studies using the nematode Caenorhabditis elegans (C. elegans) generated a wealth of publicly available transcriptomics datasets linking changes in gene expression to lifespan regulation. However, a comprehensive comparison of these datasets across studies in the context of aging biology is missing. Here, we carry out a systematic meta-analysis of over 1200 bulk RNA sequencing (RNASeq) samples obtained from 74 peer-reviewed publications on aging-related transcriptomic changes in C. elegans. Using both differential expression analyses and machine learning approaches, we mine the pooled data for novel pro-longevity genes. We find that both approaches identify known and propose novel pro-longevity genes. Further, we find that inter-lab experimental variance complicates the application of machine learning algorithms, a limitation that was not solved using bulk RNA-Seq batch correction and normalization techniques. Taken as a whole, our results indicate that machine learning approaches may hold promise for the identification of genes that regulate aging but will require more sophisticated batch correction strategies or standardized input data to reliably identify novel pro-longevity genes.

Investigating the role of Simpson's paradox in the analysis of top-ranked features in high-dimensional bioinformatics datasets

An important problem in bioinformatics consists of identifying the most important features (or predictors), among a large number of features in a given classification dataset. This problem is often addressed by using a machine learning-based feature ranking method to identify a small set of top-ranked predictors (i.e. the most relevant features for classification). The large number of studies in this area has, however, an important limitation: they ignore the possibility that the top-ranked predictors occur in an instance of Simpson's paradox, where the positive or negative association between a predictor and a class variable reverses sign upon conditional on each of the values of a third (confounder) variable. In this work, we review and investigate the role of Simpson's paradox in the analysis of top-ranked predictors in high-dimensional bioinformatics datasets, in order to avoid the potential danger of misinterpreting an association between a predictor and the class variable. We perform computational experiments using four well-known feature ranking methods from the machine learning field and five high-dimensional datasets of ageing-related genes, where the predictors are Gene Ontology terms. The results show that occurrences of Simpson's paradox involving top-ranked predictors are much more common for one of the feature ranking methods.

ORS-Pred: An optimized reduced scheme-based identifier for antioxidant proteins

Antioxidant proteins can terminate a chain of reactions caused by free radicals and protect cells from damage. To identify antioxidant proteins rapidly, a computational model was proposed based on the optimized recoding scheme, sequence information and machine learning methods. First, over 600 recoding schemes were collected to build a scheme set. Then, the original sequence was recoded as a reduced expression whose g-gap dipeptides (g = 0, 1, 2) were used as the features of proteins. Furthermore, a random forest method was used to evaluate the classification ability of the obtained dipeptide features. After going through all schemes, the best predictive performance scheme was chosen as the optimized reduction scheme. Finally, for the RF method, a grid search strategy was used to select a better parameter combination to identify antioxidant proteins. In the experiment, the present method correctly recognized 90.13-99.87% of the antioxidant samples. Other experimental results also proved that the present method was efficient to identify antioxidant proteins. Finally, we also developed a web server that was freely accessible to researchers.

Network analysis in aged C. elegans reveals candidate regulatory genes of ageing

Ageing is a biological process guided by genetic and environmental factors that ultimately lead to adverse outcomes for organismal lifespan and healthspan. Determination of molecular pathways that are affected with age and increase disease susceptibility is crucial. The gene expression profile of the ideal ageing model, namely the nematode Caenorhabditis elegans mapped with the microarray technology initially led to the identification of age-dependent gene expression alterations that characterize the nematode's ageing process. The list of differentially expressed genes was then utilized to construct a network of molecular interactions with their first neighbors/interactors using the interactions listed in the WormBase database. The subsequent network analysis resulted in the unbiased selection of 110 candidate genes, among which well-known ageing regulators appeared. More importantly, our approach revealed candidates that have never been linked to ageing before, thus suggesting promising potential targets/ageing regulators.

Computational tools for geroscience

The rapid progress of the past three decades has led the geroscience field near a point where human interventions in aging are plausible. Advances across scientific areas, such as high throughput "-omics" approaches, have led to an exponentially increasing quantity of data available for biogerontologists. To best translate the lifespan and healthspan extending interventions discovered by basic scientists into preventative medicine, it is imperative that the current data are comprehensively utilized to generate testable hypotheses about translational interventions. Building a translational pipeline for geroscience will require both systematic efforts to identify interventions that extend healthspan across taxa and diagnostics that can identify patients who may benefit from interventions prior to the onset of an age-related morbidity. Databases and computational tools that organize and analyze both the wealth of information available on basic biogerontology research and clinical data on aging populations will be critical in developing such a pipeline. Here, we review the current landscape of databases and computational resources available for translational aging research. We discuss key platforms and tools available for aging research, with a focus on how each tool can be used in concert with hypothesis driven experiments to move closer to human interventions in aging.

Computational characterization and identification of human polycystic ovary syndrome genes

Human polycystic ovary syndrome (PCOS) is a highly heritable disease regulated by genetic and environmental factors. Identifying PCOS genes is time consuming and costly in wet-lab. Developing an algorithm to predict PCOS candidates will be helpful. In this study, for the first time, we systematically analyzed properties of human PCOS genes. Compared with genes not yet known to be involved in PCOS regulation, known PCOS genes display distinguishing characteristics: (i) they tend to be located at network center; (ii) they tend to interact with each other; (iii) they tend to enrich in certain biological processes. Based on these features, we developed a machine-learning algorithm to predict new PCOS genes. 233 PCOS candidates were predicted with a posterior probability >0.9. Evidence supporting 7 of the top 10 predictions has been found.

Prediction and characterization of human ageing-related proteins by using machine learning

Ageing has a huge impact on human health and economy, but its molecular basis - regulation and mechanism - is still poorly understood. By today, more than three hundred genes (almost all of them function as protein-coding genes) have been related to human ageing. Although individual ageing-related genes or some small subsets of these genes have been intensively studied, their analysis as a whole has been highly limited. To fill this gap, for each human protein we extracted 21000 protein features from various databases, and using these data as an input to state-of-the-art machine learning methods, we classified human proteins as ageing-related or non-ageing-related. We found a simple classification model based on only 36 protein features, such as the "number of ageing-related interaction partners", "response to oxidative stress", "damaged DNA binding", "rhythmic process" and "extracellular region". Predicted values of the model quantify the relevance of a given protein in the regulation or mechanisms of the human ageing process. Furthermore, we identified new candidate proteins having strong computational evidence of their important role in ageing. Some of them, like Cytochrome b-245 light chain (CY24A) and Endoribonuclease ZC3H12A (ZC12A) have no previous ageing-associated annotations.

Towards understanding the lifespan extension by reduced insulin signaling: bioinformatics analysis of DAF-16/FOXO direct targets in Caenorhabditis elegans

DAF-16, the C. elegans FOXO transcription factor, is an important determinant in aging and longevity. In this work, we manually curated FOXODB http://lyh.pkmu.cn/foxodb/, a database of FOXO direct targets. It now covers 208 genes. Bioinformatics analysis on 109 DAF-16 direct targets in C. elegans found interesting results. (i) DAF-16 and transcription factor PQM-1 co-regulate some targets. (ii) Seventeen targets directly regulate lifespan. (iii) Four targets are involved in lifespan extension induced by dietary restriction. And (iv) DAF-16 direct targets might play global roles in lifespan regulation.

Network-based characterization and prediction of human DNA repair genes and pathways

Network biology is a useful strategy to understand cell’s functional organization. In this study, for the first time, we successfully introduced network approaches to study properties of human DNA repair genes. Compared with non-DNA repair genes, we found distinguishing features for DNA repair genes: (i) they tend to have higher degrees; (ii) they tend to be located at global network center; (iii) they tend to interact directly with each other. Based on these features, we developed the first algorithm to predict new DNA repair genes. We tested several machine-learning models and found that support vector machine with kernel function of radial basis function (RBF) achieve the best performance, with precision = 0.74 and area under curve (AUC) = 0.96. In the end, we applied the algorithm to predict new DNA repair genes and got 32 new candidates. Literature supporting four of the predictions was found. We believe the network approaches introduced here might open a new avenue to understand DNA repair genes and pathways. The suggested algorithm and the predicted genes might be helpful for scientists in the field.

A review of supervised machine learning applied to ageing research

Broadly speaking, supervised machine learning is the computational task of learning correlations between variables in annotated data (the training set), and using this information to create a predictive model capable of inferring annotations for new data, whose annotations are not known. Ageing is a complex process that affects nearly all animal species. This process can be studied at several levels of abstraction, in different organisms and with different objectives in mind. Not surprisingly, the diversity of the supervised machine learning algorithms applied to answer biological questions reflects the complexities of the underlying ageing processes being studied. Many works using supervised machine learning to study the ageing process have been recently published, so it is timely to review these works, to discuss their main findings and weaknesses. In summary, the main findings of the reviewed papers are: the link between specific types of DNA repair and ageing; ageing-related proteins tend to be highly connected and seem to play a central role in molecular pathways; ageing/longevity is linked with autophagy and apoptosis, nutrient receptor genes, and copper and iron ion transport. Additionally, several biomarkers of ageing were found by machine learning. Despite some interesting machine learning results, we also identified a weakness of current works on this topic: only one of the reviewed papers has corroborated the computational results of machine learning algorithms through wet-lab experiments. In conclusion, supervised machine learning has contributed to advance our knowledge and has provided novel insights on ageing, yet future work should have a greater emphasis in validating the predictions.

Systematic Characterization and Prediction of Human Hypertension Genes

Hypertension is a major cardiovascular risk factor and accounts for a large part of cardiovascular mortality. In this work, we analyzed the properties of hypertension genes and found that when compared with genes not yet known to be involved in hypertension regulation, known hypertension genes display distinguishing features: (1) hypertension genes tend to be located at network center; (2) hypertension genes tend to interact with each other; and (3) hypertension genes tend to enrich in certain biological processes and show certain phenotypes. Based on these features, we developed a machine-learning algorithm to predict new hypertension genes. One hundred and seventy-seven candidates were predicted with a posterior probability >0.9. Evidence supporting 17 of the predictions has been found.

PPUS: a web server to predict PUS-specific pseudouridine sites

MOTIVATION

Pseudouridine (Ψ), catalyzed by pseudouridine synthase (PUS), is the most abundant RNA modification and has important cellular functions. Developing an algorithm to identify Ψ sites is an important work. And it is better if the algorithm could assign which PUS modifies the Ψ sites. Here, we developed PPUS (http://lyh.pkmu.cn/ppus/), the first web server to predict PUS-specific Ψ sites. PPUS: employed support vector machine as the classifier and used nucleotides around Ψ sites as the features. Currently, PPUS: could accurately predict new Ψ sites for PUS1, PUS4 and PUS7 in yeast and PUS4 in human. PPUS: is well designed and friendly to user.

AVAILABILITY AND IMPLEMENTATION

Our web server is available freely for non-commercial purposes at: http://lyh.pkmu.cn/ppus/

CONTACT

liyanhui@bjmu.edu.cn or cuiqinghua@hsc.pku.edu.cn.

Predicting the Pro-Longevity or Anti-Longevity Effect of Model Organism Genes with New Hierarchical Feature Selection Methods

Ageing is a highly complex biological process that is still poorly understood. With the growing amount of ageing-related data available on the web, in particular concerning the genetics of ageing, it is timely to apply data mining methods to that data, in order to try to discover novel patterns that may assist ageing research. In this work, we introduce new hierarchical feature selection methods for the classification task of data mining and apply them to ageing-related data from four model organisms: Caenorhabditis elegans (worm), Saccharomyces cerevisiae (yeast), Drosophila melanogaster (fly), and Mus musculus (mouse). The main novel aspect of the proposed feature selection methods is that they exploit hierarchical relationships in the set of features (Gene Ontology terms) in order to improve the predictive accuracy of the Naïve Bayes and 1-Nearest Neighbour (1-NN) classifiers, which are used to classify model organisms' genes into pro-longevity or anti-longevity genes. The results show that our hierarchical feature selection methods, when used together with Naïve Bayes and 1-NN classifiers, obtain higher predictive accuracy than the standard (without feature selection) Naïve Bayes and 1-NN classifiers, respectively. We also discuss the biological relevance of a number of Gene Ontology terms very frequently selected by our algorithms in our datasets.

An evidence-based approach to identify aging-related genes in Caenorhabditis elegans

Background

Extensive studies have been carried out on Caenorhabditis elegans as a model organism to elucidate mechanisms of aging and the effects of perturbing known aging-related genes on lifespan and behavior. This research has generated large amounts of experimental data that is increasingly difficult to integrate and analyze with existing databases and domain knowledge. To address this challenge, we demonstrate a scalable and effective approach for automatic evidence gathering and evaluation that leverages existing experimental data and literature-curated facts to identify genes involved in aging and lifespan regulation in C. elegans.

Results

We developed a semantic knowledge base for aging by integrating data about C. elegans genes from WormBase with data about 2005 human and model organism genes from GenAge and 149 genes from GenDR, and with the Bio2RDF network of linked data for the life sciences. Using HyQue (a Semantic Web tool for hypothesis-based querying and evaluation) to interrogate this knowledge base, we examined 48,231 C. elegans genes for their role in modulating lifespan and aging. HyQue identified 24 novel but well-supported candidate aging-related genes for further experimental validation.

Conclusions

We use semantic technologies to discover candidate aging genes whose effects on lifespan are not yet well understood. Our customized HyQue system, the aging research knowledge base it operates over, and HyQue evaluations of all C. elegans genes are freely available at http://hyque.semanticscience.org.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0469-4) contains supplementary material, which is available to authorized users.

Deciphering the effects of gene deletion on yeast longevity using network and machine learning approaches

Longevity is one of the most basic and one of the most essential properties of all living organisms. Identification of genes that regulate longevity would increase understanding of the mechanisms of aging, so as to help facilitate anti-aging intervention and extend the life span. In this study, based on the network features and the biochemical/physicochemical features of the deletion network and deletion genes, as well as their functional features, a two-layer model was developed for predicting the deletion effects on yeast longevity. The first stage of our prediction approach was to identify whether the deletion of one gene would change the life span of yeast; if it did, the second stage of our procedure would automatically proceed to predict whether the deletion of one gene would increase or decrease the life span. It was observed by analyzing the predicted results that the functional features (such as mitochondrial function and chromatin silencing), the network features (such as the edge density and edge weight density of the deletion network), and the local centrality of deletion gene, would have important impact for predicting the deletion effects on longevity. It is anticipated that our model may become a useful tool for studying longevity from the angle of genes and networks. Moreover, it has not escaped our notice that, after some modification, the current model can also be used to study many other phenotype prediction problems from the angle of systems biology.