methyLImp2: faster missing value estimation for DNA methylation data.

MOTIVATION: methyLImp, a method we recently introduced for the missing value estimation of DNA methylation data, has demonstrated competitive performance in data imputation compared to the existing, general-purpose, approaches. However, imputation running time was considerably long and unfeasible in case of large datasets with numerous missing values. RESULTS: methyLImp2 made possible computations that were previously unfeasible. We achieved this by introducing two important modifications that have significantly reduced the original running time without sacrificing prediction performance. First, we implemented a chromosome-wise parallel version of methyLImp. This parallelization reduced the runtime by several 10-fold in our experiments. Then, to handle large datasets, we also introduced a mini-batch approach that uses only a subset of the samples for the imputation. Thus, it further reduces the running time from days to hours or even minutes in large datasets. AVAILABILITY AND IMPLEMENTATION: The R package methyLImp2 is under review for Bioconductor. It is currently freely available on Github https://github.com/annaplaksienko/methyLImp2.

Evaluation of different computational methods for DNA methylation-based biological age.

In recent years there has been a widespread interest in researching biomarkers of aging that could predict physiological vulnerability better than chronological age. Aging, in fact, is one of the most relevant risk factors for a wide range of maladies, and molecular surrogates of this phenotype could enable better patients stratification. Among the most promising of such biomarkers is DNA methylation-based biological age. Given the potential and variety of computational implementations (epigenetic clocks), we here present a systematic review of such clocks. Furthermore, we provide a large-scale performance comparison across different tissues and diseases in terms of age prediction accuracy and age acceleration, a measure of deviance from physiology. Our analysis offers both a state-of-the-art overview of the computational techniques developed so far and a heterogeneous picture of performances, which can be helpful in orienting future research.

Estimage: a webserver hub for the computation of methylation age.

Methylage is an epigenetic marker of biological age that exploits the correlation between the methylation state of specific CG dinucleotides (CpGs) and chronological age (in years), gestational age (in weeks), cellular age (in cell cycles or as telomere length, in kilobases). Using DNA methylation data, methylage is measurable via the so called epigenetic clocks. Importantly, alterations of the correlation between methylage and age (age acceleration or deceleration) have been stably associated with pathological states and occur long before clinical signs of diseases become overt, making epigenetic clocks a potentially disruptive tool in preventive, diagnostic and also in forensic applications. Nevertheless, methylage dependency from CpGs selection, mathematical modelling, tissue specificity and age range, still makes the potential of this biomarker limited. In order to enhance model comparisons, interchange, availability, robustness and standardization, we organized a selected set of clocks within a hub webservice, EstimAge (Estimate of methylation Age, http://estimage.iac.rm.cnr.it), which intuitively and informatively enables quick identification, computation and comparison of available clocks, with the support of standard statistics.

Vaccination in the elderly: The challenge of immune changes with aging.

The unprecedented increase of life expectancy challenges society to protect the elderly from morbidity and mortality making vaccination a crucial mean to safeguard this population. Indeed, infectious diseases, such as influenza and pneumonia, are among the top killers of elderly people in the world. Elderly individuals are more prone to severe infections and less responsive to vaccination prevention, due to immunosenescence combined with the progressive increase of a proinflammatory status characteristic of the aging process (inflammaging). These factors are responsible for most age-related diseases and correlate with poor response to vaccination. Therefore, it is of utmost interest to deepen the knowledge regarding the role of inflammaging in vaccination responsiveness to support the development of effective vaccination strategies designed for elderly. In this review we analyse the impact of age-associated factors such as inflammaging, immunosenescence and immunobiography on immune response to vaccination in the elderly, and we consider systems biology approaches as a mean for integrating a multitude of data in order to rationally design vaccination approaches specifically tailored for the elderly.

The epigenetics of inflammaging: The contribution of age-related heterochromatin loss and locus-specific remodelling and the modulation by environmental stimuli.

A growing amount of evidences indicates that inflammaging - the chronic, low grade inflammation state characteristic of the elderly - is the result of genetic as well as environmental or stochastic factors. Some of these, such as the accumulation of senescent cells that are persistent during aging or accompany its progression, seem to be sufficient to initiate the aging process and to fuel it. Others, like exposure to environmental compounds or infections, are temporary and resolve within a (relatively) short time. In both cases, however, a cellular memory of the event can be established by means of epigenetic modulation of the genome. In this review we will specifically discuss the relationship between epigenetics and inflammaging. In particular, we will show how age-associated epigenetic modifications concerned with heterochromatin loss and gene-specific remodelling, can promote inflammaging. Furthermore, we will recall how the exposure to specific nutritional, environmental and microbial stimuli can affect the rate of inflammaging through epigenetic mechanisms, touching also on the recent insight given by the concept of trained immunity.

Methylation data imputation performances under different representations and missingness patterns.

BACKGROUND: High-throughput technologies enable the cost-effective collection and analysis of DNA methylation data throughout the human genome. This naturally entails missing values management that can complicate the analysis of the data. Several general and specific imputation methods are suitable for DNA methylation data. However, there are no detailed studies of their performances under different missing data mechanisms -(completely) at random or not- and different representations of DNA methylation levels (ß and M-value). RESULTS: We make an extensive analysis of the imputation performances of seven imputation methods on simulated missing completely at random (MCAR), missing at random (MAR) and missing not at random (MNAR) methylation data. We further consider imputation performances on the popular ß- and M-value representations of methylation levels. Overall, ß-values enable better imputation performances than M-values. Imputation accuracy is lower for mid-range ß-values, while it is generally more accurate for values at the extremes of the ß-value range. The MAR values distribution is on the average more dense in the mid-range in comparison to the expected ß-value distribution. As a consequence, MAR values are on average harder to impute. CONCLUSIONS: The results of the analysis provide guidelines for the most suitable imputation approaches for DNA methylation data under different representations of DNA methylation levels and different missing data mechanisms.

Missing value estimation methods for DNA methylation data.

MOTIVATION: DNA methylation is a stable epigenetic mark with major implications in both physiological (development, aging) and pathological conditions (cancers and numerous diseases). Recent research involving methylation focuses on the development of molecular age estimation methods based on DNA methylation levels (mAge). An increasing number of studies indicate that divergences between mAge and chronological age may be associated to age-related diseases. Current advances in high-throughput technologies have allowed the characterization of DNA methylation levels throughout the human genome. However, experimental methylation profiles often contain multiple missing values that can affect the analysis of the data and also mAge estimation. Although several imputation methods exist, a major deficiency lies in the inability to cope with large datasets, such as DNA methylation chips. Specific methods for imputing missing methylation data are therefore needed. RESULTS: We present a simple and computationally efficient imputation method, metyhLImp, based on linear regression. The rationale of the approach lies in the observation that methylation levels show a high degree of inter-sample correlation. We performed a comparative study of our approach with other imputation methods on DNA methylation data of healthy and disease samples from different tissues. Performances have been assessed both in terms of imputation accuracy and in terms of the impact imputed values have on mAge estimation. In comparison to existing methods, our linear regression model proves to perform equally or better and with good computational efficiency. The results of our analysis provide recommendations for accurate estimation of missing methylation values. AVAILABILITY AND IMPLEMENTATION: The R-package methyLImp is freely available at https://github.com/pdilena/methyLImp. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Systemic evaluation of cellular reprogramming processes exploiting a novel R-tool: eegc.

MOTIVATION: Cells derived by cellular engineering, i.e. differentiation of induced pluripotent stem cells and direct lineage reprogramming, carry a tremendous potential for medical applications and in particular for regenerative therapies. These approaches consist in the definition of lineage-specific experimental protocols that, by manipulation of a limited number of biological cues-niche mimicking factors, (in)activation of transcription factors, to name a few-enforce the final expression of cell-specific (marker) molecules. To date, given the intricate complexity of biological pathways, these approaches still present imperfect reprogramming fidelity, with uncertain consequences on the functional properties of the resulting cells. RESULTS: We propose a novel tool eegc to evaluate cellular engineering processes, in a systemic rather than marker-based fashion, by integrating transcriptome profiling and functional analysis. Our method clusters genes into categories representing different states of (trans)differentiation and further performs functional and gene regulatory network analyses for each of the categories of the engineered cells, thus offering practical indications on the potential lack of the reprogramming protocol. AVAILABILITY AND IMPLEMENTATION: eegc R package is released under the GNU General Public License within the Bioconductor project, freely available at https://bioconductor.org/packages/eegc/. CONTACT: christine.nardini.rsrc@gmail.com or hongkang.k.mei@gsk.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Mechanotransduction map: simulation model, molecular pathway, gene set.

MOTIVATION: Mechanotransduction--the ability to output a biochemical signal from a mechanical input--is related to the initiation and progression of a broad spectrum of molecular events. Yet, the characterization of mechanotransduction lacks some of the most basic tools as, for instance, it can hardly be recognized by enrichment analysis tools, nor could we find any pathway representation. This greatly limits computational testing and hypothesis generation on mechanotransduction biological relevance and involvement in disease or physiological mechanisms. RESULTS: We here present a molecular map of mechanotransduction, built in CellDesigner to warrant that maximum information is embedded in a compact network format. To validate the map's necessity we tested its redundancy in comparison with existing pathways, and to estimate its sufficiency, we quantified its ability to reproduce biological events with dynamic simulations, using Signaling Petri Networks. AVAILABILITY AND IMPLEMENTATION: SMBL language map is available in the Supplementary Data: core_map.xml, basic_map.xml. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SPNConverter: a new link between static and dynamic complex network analysis.

SUMMARY: The signaling Petri net (SPN) simulator, designed to provide insights into the trends of molecules' activity levels in response to an external stimulus, contributes to the systems biology necessity of analyzing the dynamics of large-scale cellular networks. Implemented into the freely available software, BioLayout Express(3D), the simulator is publicly available and easy to use, provided the input files are prepared in the GraphML format, typically using the network editing software, yEd, and standards specific to the software. However, analysis of complex networks represented using other systems biology formatting languages (on which popular software, such as CellDesigner and Cytoscape, are based) requires manual manipulation, a step that is prone to error and limits the use of the SPN simulator in BioLayout Express(3D). To overcome this, we present a Cytoscape plug-in that enables users to automatically convert networks for analysis with the SPN simulator from the standard systems biology markup language. The automation of this step opens the SPN simulator to a far larger user group than has previously been possible. AVAILABILITY AND IMPLEMENTATION: Distributed under the GNU General Public License Version 3 at http://apps.cytoscape.org/apps/spnconverter.

Where are we in the implementation of tissue-specific epigenetic clocks?

Introduction: DNA methylation clocks presents advantageous characteristics with respect to the ambitious goal of identifying very early markers of disease, based on the concept that accelerated ageing is a reliable predictor in this sense. Methods: Such tools, being epigenomic based, are expected to be conditioned by sex and tissue specificities, and this work is about quantifying this dependency as well as that from the regression model and the size of the training set. Results: Our quantitative results indicate that elastic-net penalization is the best performing strategy, and better so when-unsurprisingly-the data set is bigger; sex does not appear to condition clocks performances and tissue specific clocks appear to perform better than generic blood clocks. Finally, when considering all trained clocks, we identified a subset of genes that, to the best of our knowledge, have not been presented yet and might deserve further investigation: CPT1A, MMP15, SHROOM3, SLIT3, and SYNGR. Conclusion: These factual starting points can be useful for the future medical translation of clocks and in particular in the debate between multi-tissue clocks, generally trained on a large majority of blood samples, and tissue-specific clocks.

MIMO: an efficient tool for molecular interaction maps overlap.

BACKGROUND: Molecular pathways represent an ensemble of interactions occurring among molecules within the cell and between cells. The identification of similarities between molecular pathways across organisms and functions has a critical role in understanding complex biological processes. For the inference of such novel information, the comparison of molecular pathways requires to account for imperfect matches (flexibility) and to efficiently handle complex network topologies. To date, these characteristics are only partially available in tools designed to compare molecular interaction maps. RESULTS: Our approach MIMO (Molecular Interaction Maps Overlap) addresses the first problem by allowing the introduction of gaps and mismatches between query and template pathways and permits -when necessary- supervised queries incorporating a priori biological information. It then addresses the second issue by relying directly on the rich graph topology described in the Systems Biology Markup Language (SBML) standard, and uses multidigraphs to efficiently handle multiple queries on biological graph databases. The algorithm has been here successfully used to highlight the contact point between various human pathways in the Reactome database. CONCLUSIONS: MIMO offers a flexible and efficient graph-matching tool for comparing complex biological pathways.

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.

Inference of gene networks--application to Bifidobacterium.

MOTIVATION: The reliable and reproducible identification of gene interaction networks represents one of the grand challenges of both modern molecular biology and computational sciences. Approaches based on careful collection of literature data and network topological analysis, applied to unicellular organisms, have proven to offer results applicable to medical therapies. However, when little a priori knowledge is available, other approaches, not relying so strongly on previous literature, must be used. We propose here a novel algorithm (based on ordinary differential equations) able to infer the interactions occurring among genes, starting from gene expression steady state data. RESULTS: The algorithm was first validated on synthetic and real benchmarks. It was then applied to the reconstruction of the core of the amino acids metabolism in Bifidobacterium longum, an essential, yet poorly known player in the human gut intestinal microbiome, known to be related to the onset of important diseases, such as metabolic syndromes. Our results show how computational approaches can offer effective tools for applications with the identification of potential new biological information. AVAILABILITY: The software is available at www.bioconductor.org and at www.picb.ac.cn/ClinicalGenomicNTW/temp2.html.

Semi-automated socio-anthropologic analysis of the medical discourse on rheumatoid arthritis: Potential impact on public health.

BACKGROUND: The debilitating effects of noncommunicable diseases (NCDs) and the accompanying chronic inflammation represent a significant obstacle for the sustainability of our development, with efforts spreading worldwide to counteract the diffusion of NCDs, as per the United Nations Sustainable Development Goals (SDG 3). In fact, despite efforts of varied intensity in numerous directions (from innovations in biotechnology to lifestyle modifications), the incidence of NCDs remains pandemic. The present work wants to contribute to addressing this major concern, with a specific focus on the fragmentation of medical approaches, via an interdisciplinary analysis of the medical discourse, i.e. the heterogenous reporting that biomedical scientific literature uses to describe the anti-inflammatory therapeutic landscape in NCDs. The aim is to better capture the roots of this compartmentalization and the power relations existing among three segregated pharmacological, experimental and unstandardized biomedical approaches to ultimately empower collaboration beyond medical specialties and possibly tap into a more ample and effective reservoir of integrated therapeutic opportunities. METHOD: Using rheumatoid arthritis (RA) as an exemplar disease, twenty-eight articles were manually translated into a nine-dimensional categorical variable of medical socio-anthropological relevance, relating in particular (but not only) to legitimacy, temporality and spatialization. This digitalized picture (9 x 28 table) of the medical discourse was further analyzed by simple automated learning approaches to identify differences and highlight commonalities among the biomedical categories. RESULTS: Interpretation of these results provides original insights, including suggestions to: empower scientific communication between unstandardized approaches and basic biology; promote the repurposing of non-pharmacological therapies to enhance robustness of experimental approaches; and align the spatial representation of diseases and therapies in pharmacology to effectively embrace the systemic approach promoted by modern personalized and preventive medicines. We hope this original work can expand and foster interdisciplinarity among public health stakeholders, ultimately contributing to the achievement of SDG3.

Early downregulation of hsa-miR-144-3p in serum from drug-naïve Parkinson's disease patients.

Advanced age represents one of the major risk factors for Parkinson's Disease. Recent biomedical studies posit a role for microRNAs, also known to be remodelled during ageing. However, the relationship between microRNA remodelling and ageing in Parkinson's Disease, has not been fully elucidated. Therefore, the aim of the present study is to unravel the relevance of microRNAs as biomarkers of Parkinson's Disease within the ageing framework. We employed Next Generation Sequencing to profile serum microRNAs from samples informative for Parkinson's Disease (recently diagnosed, drug-naïve) and healthy ageing (centenarians) plus healthy controls, age-matched with Parkinson's Disease patients. Potential microRNA candidates markers, emerging from the combination of differential expression and network analyses, were further validated in an independent cohort including both drug-naïve and advanced Parkinson's Disease patients, and healthy siblings of Parkinson's Disease patients at higher genetic risk for developing the disease. While we did not find evidences of microRNAs co-regulated in Parkinson's Disease and ageing, we report that hsa-miR-144-3p is consistently down-regulated in early Parkinson's Disease patients. Moreover, interestingly, functional analysis revealed that hsa-miR-144-3p is involved in the regulation of coagulation, a process known to be altered in Parkinson's Disease. Our results consistently show the down-regulation of hsa-mir144-3p in early Parkinson's Disease, robustly confirmed across a variety of analytical and experimental analyses. These promising results ask for further research to unveil the functional details of the involvement of hsa-mir144-3p in Parkinson's Disease.

Joint analysis of transcriptional and post- transcriptional brain tumor data: searching for emergent properties of cellular systems.

BACKGROUND: Advances in biotechnology offer a fast growing variety of high-throughput data for screening molecular activities of genomic, transcriptional, post-transcriptional and translational observations. However, to date, most computational and algorithmic efforts have been directed at mining data from each of these molecular levels (genomic, transcriptional, etc.) separately. In view of the rapid advances in technology (new generation sequencing, high-throughput proteomics) it is important to address the problem of analyzing these data as a whole, i.e. preserving the emergent properties that appear in the cellular system when all molecular levels are interacting. We analyzed one of the (currently) few datasets that provide both transcriptional and post-transcriptional data of the same samples to investigate the possibility to extract more information, using a joint analysis approach. RESULTS: We use Factor Analysis coupled with pre-established knowledge as a theoretical base to achieve this goal. Our intention is to identify structures that contain information from both mRNAs and miRNAs, and that can explain the complexity of the data. Despite the small sample available, we can show that this approach permits identification of meaningful structures, in particular two polycistronic miRNA genes related to transcriptional activity and likely to be relevant in the discrimination between gliosarcomas and other brain tumors. CONCLUSIONS: This suggests the need to develop methodologies to simultaneously mine information from different levels of biological organization, rather than linking separate analyses performed in parallel.

Identification of noninvasive imaging surrogates for brain tumor gene-expression modules.

Glioblastoma multiforme (GBM) is the most common and lethal primary brain tumor in adults. We combined neuroimaging and DNA microarray analysis to create a multidimensional map of gene-expression patterns in GBM that provided clinically relevant insights into tumor biology. Tumor contrast enhancement and mass effect predicted activation of specific hypoxia and proliferation gene-expression programs, respectively. Overexpression of EGFR, a receptor tyrosine kinase and potential therapeutic target, was also directly inferred by neuroimaging and was validated in an independent set of tumors by immunohistochemistry. Furthermore, imaging provided insights into the intratumoral distribution of gene-expression patterns within GBM. Most notably, an "infiltrative" imaging phenotype was identified that predicted patient outcome. Patients with this imaging phenotype had a greater tendency toward having multiple tumor foci and demonstrated significantly shorter survival than their counterparts. Our findings provide an in vivo portrait of genome-wide gene expression in GBM and offer a potential strategy for noninvasively selecting patients who may be candidates for individualized therapies.

The evolution of personalized healthcare and the pivotal role of European regions in its implementation.

Personalized medicine (PM) moves at the same pace of data and technology and calls for important changes in healthcare. New players are participating, providing impulse to PM. We review the conceptual foundations for PM and personalized healthcare and their evolution through scientific publications where a clear definition and the features of the different formulations are identifiable. We then examined PM policy documents of the International Consortium for Personalised Medicine and related initiatives to understand how PM stakeholders have been changing. Regional authorities and stakeholders have joined the race to deliver personalized care and are driving toward what could be termed as the next personalized healthcare. Their role as a key stakeholder in PM is expected to be pivotal.

A geroscience approach for Parkinson's disease: Conceptual framework and design of PROPAG-AGEING project.

Advanced age is the major risk factor for idiopathic Parkinson's disease (PD), but to date the biological relationship between PD and ageing remains elusive. Here we describe the rationale and the design of the H2020 funded project "PROPAG-AGEING", whose aim is to characterize the contribution of the ageing process to PD development. We summarize current evidences that support the existence of a continuum between ageing and PD and justify the use of a Geroscience approach to study PD. We focus in particular on the role of inflammaging, the chronic, low-grade inflammation characteristic of elderly physiology, which can propagate and transmit both locally and systemically. We then describe PROPAG-AGEING design, which is based on the multi-omic characterization of peripheral samples from clinically characterized drug-naïve and advanced PD, PD discordant twins, healthy controls and "super-controls", i.e. centenarians, who never showed clinical signs of motor disability, and their offspring. Omic results are then validated in a large number of samples, including in vitro models of dopaminergic neurons and healthy siblings of PD patients, who are at higher risk of developing PD, with the final aim of identifying the molecular perturbations that can deviate the trajectories of healthy ageing towards PD development.