Sex-specific differences in myocardial glucose metabolic rate in non-diabetic, pre-diabetic and type 2 diabetic subjects

BACKGROUND

Evidence has shown that women with type 2 diabetes (T2DM) have a higher excess risk for cardiovascular disease (CVD) than men with T2DM. Subjects with either T2DM or prediabetes exhibit myocardial insulin resistance, but it is still unsettled whether sex-related differences in myocardial insulin resistance occur in diabetic and prediabetic subjects.

METHODS

We aimed to evaluate sex-related differences in myocardial glucose metabolic rate (MRGlu), assessed using dynamic PET with 18F-FDG combined with euglycemic-hyperinsulinemic clamp, in subjects with normal glucose tolerance (NGT; n = 20), prediabetes (n = 11), and T2DM (n = 26).

RESULTS

Women with prediabetes or T2DM exhibited greater relative differences in myocardial MRGlu than men with prediabetes or T2DM when compared with their NGT counterparts. As compared with women with NGT, those with prediabetes exhibited an age-adjusted 35% lower myocardial MRGlu value (P = 0.04) and women with T2DM a 74% lower value (P = 0.006), respectively. Conversely, as compared with men with NGT, men with T2DM exhibited a 40% lower myocardial MRGlu value (P = 0.004), while no significant difference was observed between men with NGT and prediabetes. The statistical test for interaction between sex and glucose tolerance on myocardial MRGlu (P < 0.0001) was significant suggesting a sex-specific association.

CONCLUSIONS

Our data suggest that deterioration of glucose homeostasis in women is associated with a greater impairment in myocardial glucose metabolism as compared with men. The sex-specific myocardial insulin resistance could be an important factor responsible for the greater effect of T2DM on the excess risk of cardiovascular disease in women than in men.

Machine learning pipeline to analyze clinical and proteomics data: experiences on a prostate cancer case

Proteomic-based analysis is used to identify biomarkers in blood samples and tissues. Data produced by devices such as mass spectrometry requires platforms to identify and quantify proteins (or peptides). Clinical information can be related to mass spectrometry data to identify diseases at an early stage. Machine learning techniques can be used to support physicians and biologists in studying and classifying pathologies. We present the application of machine learning techniques to define a pipeline aimed at studying and classifying proteomics data enriched using clinical information. The pipeline allows users to relate established blood biomarkers with clinical parameters and proteomics data. The proposed pipeline entails three main phases: (i) feature selection, (ii) models training, and (iii) models ensembling. We report the experience of applying such a pipeline to prostate-related diseases. Models have been trained on several biological datasets. We report experimental results about two datasets that result from the integration of clinical and mass spectrometry-based data in the contexts of serum and urine analysis. The pipeline receives input data from blood analytes, tissue samples, proteomic analysis, and urine biomarkers. It then trains different models for feature selection, classification and voting. The presented pipeline has been applied on two datasets obtained in a 2 years research project which aimed to extract hidden information from mass spectrometry, serum, and urine samples from hundreds of patients. We report results on analyzing prostate datasets serum with 143 samples, including 79 PCa and 84 BPH patients, and an urine dataset with 121 samples, including 67 PCa and 54 BPH patients. As results pipeline allowed to identify interesting peptides in the two datasets, 6 for the first one and 2 for the second one. The best model for both serum (AUC=0.87, Accuracy=0.83, F1=0.81, Sensitivity=0.84, Specificity=0.81) and urine (AUC=0.88, Accuracy=0.83, F1=0.83, Sensitivity=0.85, Specificity=0.80) datasets showed good predictive performances. We made the pipeline code available on GitHub and we are confident that it will be successfully adopted in similar clinical setups.

The Omicron XBB.1 Variant and Its Descendants: Genomic Mutations, Rapid Dissemination and Notable Characteristics

The SARS-CoV-2 virus, which is a major threat to human health, has undergone many mutations during the replication process due to errors in the replication steps and modifications in the structure of viral proteins. The XBB variant was identified for the first time in Singapore in the fall of 2022. It was then detected in other countries, including the United States, Canada, and the United Kingdom. We study the impact of sequence changes on spike protein structure on the subvariants of XBB, with particular attention to the velocity of variant diffusion and virus activity with respect to its diffusion. We examine the structural and functional distinctions of the variants in three different conformations: (i) spike glycoprotein in complex with ACE2 (1-up state), (ii) spike glycoprotein (closed-1 state), and (iii) S protein (open-1 state). We also estimate the affinity binding between the spike protein and ACE2. The market binding affinity observed in specific variants raises questions about the efficacy of current vaccines in preparing the immune system for virus variant recognition. This work may be useful in devising strategies to manage the ongoing COVID-19 pandemic. To stay ahead of the virus evolution, further research and surveillance should be carried out to adjust public health measures accordingly.

REHABS: An Innovative and User-Friendly Device for Rehabilitation

Rehabilitation is a complex set of interventions involving the assessment, management, and treatment of injuries. It aims to support and facilitate an individual's recovery process by restoring a physiological function, e.g., limb movement, compromised by physical impairments, injuries or diseases to a condition as close to normal as possible. Innovative devices and solutions make the rehabilitation process of patients easier during their daily activities. Devices support physicians and physiotherapists in monitoring and measuring patients' physical improvements during rehabilitation. In this context, we report the design and implementation of a low-cost rehabilitation system, which is a programmable device designed to support tele-rehabilitation of the upper limbs. The proposed system includes a mechanism to acquire and analyze data and signals related to rehabilitation processes.

Development of a predictive model to distinguish prostate cancer from benign prostatic hyperplasia by integrating serum glycoproteomics and clinical variables

BACKGROUND

Prostate Cancer (PCa) represents the second leading cause of cancer-related death in men. Prostate-specific antigen (PSA) serum testing, currently used for PCa screening, lacks the necessary sensitivity and specificity. New non-invasive diagnostic tools able to discriminate tumoral from benign conditions and aggressive (AG-PCa) from indolent forms of PCa (NAG-PCa) are required to avoid unnecessary biopsies.

METHODS

In this work, 32 formerly N-glycosylated peptides were quantified by PRM (parallel reaction monitoring) in 163 serum samples (79 from PCa patients and 84 from individuals affected by benign prostatic hyperplasia (BPH)) in two technical replicates. These potential biomarker candidates were prioritized through a multi-stage biomarker discovery pipeline articulated in: discovery, LC-PRM assay development and verification phases. Because of the well-established involvement of glycoproteins in cancer development and progression, the proteomic analysis was focused on glycoproteins enriched by TiO2 (titanium dioxide) strategy.

RESULTS

Machine learning algorithms have been applied to the combined matrix comprising proteomic and clinical variables, resulting in a predictive model based on six proteomic variables (RNASE1, LAMP2, LUM, MASP1, NCAM1, GPLD1) and five clinical variables (prostate dimension, proPSA, free-PSA, total-PSA, free/total-PSA) able to distinguish PCa from BPH with an area under the Receiver Operating Characteristic (ROC) curve of 0.93. This model outperformed PSA alone which, on the same sample set, was able to discriminate PCa from BPH with an AUC of 0.79. To improve the clinical managing of PCa patients, an explorative small-scale analysis (79 samples) aimed at distinguishing AG-PCa from NAG-PCa was conducted. A predictor of PCa aggressiveness based on the combination of 7 proteomic variables (FCN3, LGALS3BP, AZU1, C6, LAMB1, CHL1, POSTN) and proPSA was developed (AUC of 0.69).

CONCLUSIONS

To address the impelling need of more sensitive and specific serum diagnostic tests, a predictive model combining proteomic and clinical variables was developed. A preliminary evaluation to build a new tool able to discriminate aggressive presentations of PCa from tumors with benign behavior was exploited. This predictor displayed moderate performances, but no conclusions can be drawn due to the limited number of the sample cohort. Data are available via ProteomeXchange with identifier PXD035935.

Genome mining conformance to metabolite profile of Bacillus strains to control potato pathogens

Biocontrol agents are safe and effective methods for controlling plant disease pathogens, such as Fusarium solani, which causes dry wilt, and Pectobacterium spp., responsible for potato soft rot disease. Discovering agents that can effectively control both fungal and bacterial pathogens in potatoes has always presented a challenge. Biological controls were investigated using 500 bacterial strains isolated from rhizospheric microbial communities, along with two promising biocontrol strains: Pseudomonas (T17-4 and VUPf5). Bacillus velezensis (Q12 and US1) and Pseudomonas chlororaphis VUPf5 exhibited the highest inhibition of fungal growth and pathogenicity in both laboratory (48%, 48%, 38%) and greenhouse (100%, 85%, 90%) settings. Q12 demonstrated better control against bacterial pathogens in vivo (approximately 50%). Whole-genome sequencing of Q12 and US1 revealed a genome size of approximately 4.1 Mb. Q12 had 4413 gene IDs and 4300 coding sequences, while US1 had 4369 gene IDs and 4255 coding sequences. Q12 exhibited a higher number of genes classified under functional subcategories related to stress response, cell wall, capsule, levansucrase synthesis, and polysaccharide metabolism. Both Q12 and US1 contained eleven secondary metabolite gene clusters as identified by the antiSMASH and RAST servers. Notably, Q12 possessed the antibacterial locillomycin and iturin A gene clusters, which were absent in US1. This genetic information suggests that Q12 may have a more pronounced control over bacterial pathogens compared to US1. Metabolic profiling of the superior strains, as determined by LC/MS/MS, validated our genetic findings. The investigated strains produced compounds such as iturin A, bacillomycin D, surfactin, fengycin, phenazine derivatives, etc. These compounds reduced spore production and caused deformation of the hyphae in F. solani. In contrast, B. velezensis UR1, which lacked the production of surfactin, fengycin, and iturin, did not affect these structures and failed to inhibit the growth of any pathogens. Our findings suggest that locillomycin and iturin A may contribute to the enhanced control of bacterial pectolytic rot by Q12.

Strategies and Trends in COVID-19 Vaccination Delivery: What We Learn and What We May Use for the Future

Vaccination has been the most effective way to control the outbreak of the COVID-19 pandemic. The numbers and types of vaccines have reached considerable proportions, even if the question of vaccine procedures and frequency still needs to be resolved. We have come to learn the necessity of defining vaccination distribution strategies with regard to COVID-19 that could be used for any future pandemics of similar gravity. In fact, vaccine monitoring implies the existence of a strategy that should be measurable in terms of input and output, based on a mathematical model, including death rates, the spread of infections, symptoms, hospitalization, and so on. This paper addresses the issue of vaccine diffusion and strategies for monitoring the pandemic. It provides a description of the importance and take up of vaccines and the links between procedures and the containment of COVID-19 variants, as well as the long-term effects. Finally, the paper focuses on the global scenario in a world undergoing profound social and political change, with particular attention on current and future health provision. This contribution would represent an example of vaccination experiences, which can be useful in other pandemic or epidemiological contexts.

A systems biology approach to pathogenesis of gastric cancer: gene network modeling and pathway analysis

BACKGROUND

Gastric cancer (GC) ranks among the most common malignancies worldwide. This study aimed to find critical genes/pathways in GC pathogenesis.

METHODS

Gene interactions were analyzed, and the protein-protein interaction network was drawn. Then enrichment analysis of the hub genes was performed and network cluster analysis and promoter analysis of the hub genes were done. Age/sex analysis was done on the identified genes.

RESULTS

Eleven hub genes in GC were identified in the current study (ATP5A1, ATP5B, ATP5D, MT-ATP8, COX7A2, COX6C, ND4, ND6, NDUFS3, RPL8, and RPS16), mostly involved in mitochondrial functions. There was no report on the ATP5D, ND6, NDUFS3, RPL8, and RPS16 in GC. Our results showed that the most affected processes in GC are the metabolic processes, and the oxidative phosphorylation pathway was considerably enriched which showed the significance of mitochondria in GC pathogenesis. Most of the affected pathways in GC were also involved in neurodegenerative diseases. Promoter analysis showed that negative regulation of signal transduction might play an important role in GC pathogenesis. In the analysis of the basal expression pattern of the selected genes whose basal expression presented a change during the age, we found that a change in age may be an indicator of changes in disease insurgence and/or progression at different ages.

CONCLUSIONS

These results might open up new insights into GC pathogenesis. The identified genes might be novel diagnostic/prognostic biomarkers or potential therapeutic targets for GC. This work, being based on bioinformatics analysis act as a hypothesis generator that requires further clinical validation.

SARS-CoV-2 protein structure and sequence mutations: Evolutionary analysis and effects on virus variants

The structure and sequence of proteins strongly influence their biological functions. New models and algorithms can help researchers in understanding how the evolution of sequences and structures is related to changes in functions. Recently, studies of SARS-CoV-2 Spike (S) protein structures have been performed to predict binding receptors and infection activity in COVID-19, hence the scientific interest in the effects of virus mutations due to sequence, structure and vaccination arises. However, there is the need for models and tools to study the links between the evolution of S protein sequence, structure and functions, and virus transmissibility and the effects of vaccination. As studies on S protein have been generated a large amount of relevant information, we propose in this work to use Protein Contact Networks (PCNs) to relate protein structures with biological properties by means of network topology properties. Topological properties are used to compare the structural changes with sequence changes. We find that both node centrality and community extraction analysis can be used to relate protein stability and functionality with sequence mutations. Starting from this we compare structural evolution to sequence changes and study mutations from a temporal perspective focusing on virus variants. Finally by applying our model to the Omicron variant we report a timeline correlation between Omicron and the vaccination campaign.

Aligning Cross-Species Interactomes for Studying Complex and Chronic Diseases

Neurodegenerative diseases (NDs) are a group of complex disorders characterized by the progressive degeneration and dysfunction of neurons in the central nervous system. NDs encompass many conditions, including Alzheimer's disease and Parkinson's disease. Alzheimer's disease (AD) is a complex disease affecting almost forty million people worldwide. AD is characterized by a progressive decline of cognitive functions related to the loss of connections between nerve cells caused by the prevalence of extracellular Aβ plaques and intracellular neurofibrillary tangles plaques. Parkinson's disease (PD) is a neurodegenerative disorder that primarily affects the movement of an individual. The exact cause of Parkinson's disease is not fully understood, but it is believed to involve a combination of genetic and environmental factors. Some cases of PD are linked to mutations in the LRRK2, PARKIN and other genes, which are associated with familial forms of the disease. Different research studies have applied the Protein Protein Interaction (PPI) networks to understand different aspects of disease progression. For instance, Caenorhabditis elegans is widely used as a model organism for the study of AD due to roughly 38% of its genes having a human ortholog. This study's goal consists of comparing PPI network of C. elegans and human by applying computational techniques, widely used for the analysis of PPI networks between species, such as Local Network Alignment (LNA). For this aim, we used L-HetNetAligner algorithm to build a local alignment among two PPI networks, i.e., C. elegans and human PPI networks associated with AD and PD built-in silicon. The results show that L-HetNetAligner can find local alignments representing functionally related subregions. In conclusion, since local alignment enables the extraction of functionally related modules, the method can be used to study complex disease progression.

Analysis of age-dependent gene-expression in human tissues for studying diabetes comorbidities

The study of the relationship between type 2 diabetes mellitus (T2DM) disease and other pathologies (comorbidities), together with patient age variation, poses a challenge for medical research. There is evidence that patients affected by T2DM are more likely to develop comorbidities as they grow older. Variation of gene expression can be correlated to changes in T2DM comorbidities insurgence and progression. Understanding gene expression changes requires the analysis of large heterogeneous data at different scales as well as the integration of different data sources into network medicine models. Hence, we designed a framework to shed light on uncertainties related to age effects and comorbidity by integrating existing data sources with novel algorithms. The framework is based on integrating and analysing existing data sources under the hypothesis that changes in the basal expression of genes may be responsible for the higher prevalence of comorbidities in older patients. Using the proposed framework, we selected genes related to comorbidities from existing databases, and then analysed their expression with age at the tissues level. We found a set of genes that changes significantly in certain specific tissues over time. We also reconstructed the associated protein interaction networks and the related pathways for each tissue. Using this mechanistic framework, we detected interesting pathways related to T2DM whose genes change their expression with age. We also found many pathways related to insulin regulation and brain activities, which can be used to develop specific therapies. To the best of our knowledge, this is the first study that analyses such genes at the tissue level together with age variations.

Computational analysis of the sequence-structure relation in SARS-CoV-2 spike protein using protein contact networks

The structure of proteins impacts directly on the function they perform. Mutations in the primary sequence can provoke structural changes with consequent modification of functional properties. SARS-CoV-2 proteins have been extensively studied during the pandemic. This wide dataset, related to sequence and structure, has enabled joint sequence-structure analysis. In this work, we focus on the SARS-CoV-2 S (Spike) protein and the relations between sequence mutations and structure variations, in order to shed light on the structural changes stemming from the position of mutated amino acid residues in three different SARS-CoV-2 strains. We propose the use of protein contact network (PCN) formalism to: (i) obtain a global metric space and compare various molecular entities, (ii) give a structural explanation of the observed phenotype, and (iii) provide context dependent descriptors of single mutations. PCNs have been used to compare sequence and structure of the Alpha, Delta, and Omicron SARS-CoV-2 variants, and we found that omicron has a unique mutational pattern leading to different structural consequences from mutations of other strains. The non-random distribution of changes in network centrality along the chain has allowed to shed light on the structural (and functional) consequences of mutations.

How Well Did the Healthcare System Respond to the Healthcare Needs of Older People with and without Dementia during the COVID-19 Pandemic? The Perception of Healthcare Providers and Older People from the SI4CARE Project in the ADRION Region

One major challenge during the COVID-19 pandemic was the limited accessibility to healthcare facilities, especially for the older population. The aim of the current study was the exploration of the extent to which the healthcare systems responded to the healthcare needs of the older people with or without cognitive impairment and their caregivers in the Adrion/Ionian region. Data were collected through e-questionnaires regarding the adequacy of the healthcare system and were anonymously administered to older individuals and stakeholder providers in the following countries: Slovenia, Italy (Calabria), Croatia, Bosnia and Herzegovina, Greece, Montenegro, and Serbia. Overall, 722 older people and 267 healthcare stakeholders participated in the study. During the COVID-19 pandemic, both healthcare stakeholders and the older population claimed that the healthcare needs of the older people and their caregivers increased dramatically in all countries, especially in Italy (Calabria), Croatia and BiH. According to our results, countries from the Adrion/Ionian regions faced significant challenges to adjust to the special needs of the older people during the COVID-19 pandemic, which was possibly due to limited accessibility opportunities to healthcare facilities. These results highlight the need for the development of alternative ways of providing medical assistance and supervision when in-person care is not possible.

Impaired insulin-stimulated myocardial glucose metabolic rate is associated with reduced estimated myocardial energetic efficiency in subjects with different degrees of glucose tolerance

BACKGROUND

Alterations in myocardial mechano-energetic efficiency (MEEi), which represents the capability of the left ventricles to convert the chemical energy obtained by oxidative metabolism into mechanical work, have been associated with cardiovascular disease. Although whole-body insulin resistance has been related to impaired myocardial MEEi, it is unknown the relationship between cardiac insulin resistance and MEEi. Aim of this study was to evaluate the relationship between insulin-stimulated myocardial glucose metabolic rate (MrGlu) and myocardial MEEi in subjects having different degrees of glucose tolerance.

METHODS

We evaluated insulin-stimulated myocardial MrGlu using cardiac dynamic positron emission tomography (PET) with 18F-Fluorodeoxyglucose (18F-FDG) combined with euglycemic-hyperinsulinemic clamp, and myocardial MEEi in 57 individuals without history of coronary heart disease having different degrees of glucose tolerance. The subjects were stratified into tertiles according to their myocardial MrGlu values.

RESULTS

After adjusting for age, gender and BMI, subjects in I tertile showed a decrease in myocardial MEEi (0.31 ± 0.05 vs 0.42 ± 0.14 ml/s*g, P = 0.02), and an increase in myocardial oxygen consumption (MVO2) (10,153 ± 1375 vs 7816 ± 1229 mmHg*bpm, P < 0.0001) as compared with subjects in III tertile. Univariate correlations showed that insulin-stimulated myocardial MrGlu was positively correlated with MEEi and whole-body glucose disposal, and negatively correlated with waist circumference, fasting plasma glucose, HbA1c and MVO2. In a multivariate regression analysis running a model including several CV risk factors, the only variable that remained significantly associated with MEEi was myocardial MrGlu (β 0.346; P = 0.01).

CONCLUSIONS

These data suggest that an impairment in insulin-stimulated myocardial glucose metabolism is an independent contributor of depressed myocardial MEEi in subjects without history of CHD.

Temporal networks in biology and medicine: a survey on models, algorithms, and tools

The use of static graphs for modelling and analysis of biological and biomedical data plays a key role in biomedical research. However, many real-world scenarios present dynamic behaviours resulting in both node and edges modification as well as feature evolution. Consequently, ad-hoc models for capturing these evolutions along the time have been introduced, also referred to as dynamic, temporal, time-varying graphs. Here, we focus on temporal graphs, i.e., graphs whose evolution is represented by a sequence of time-ordered snapshots. Each snapshot represents a graph active in a particular timestamp. We survey temporal graph models and related algorithms, presenting fundamentals aspects and the recent advances. We formally define temporal graphs, focusing on the problem setting and we present their main applications in biology and medicine. We also present temporal graph embedding and the application to recent problems such as epidemic modelling. Finally, we further state some promising research directions in the area. Main results of this study include a systematic review of fundamental temporal network problems and their algorithmic solutions considered in the literature, in particular those having application in computational biology and medicine. We also include the main software developed in this context.

Design and Implementation of a New Local Alignment Algorithm for Multilayer Networks

Network alignment (NA) is a popular research field that aims to develop algorithms for comparing networks. Applications of network alignment span many fields, from biology to social network analysis. NA comes in two forms: global network alignment (GNA), which aims to find a global similarity, and LNA, which aims to find local regions of similarity. Recently, there has been an increasing interest in introducing complex network models such as multilayer networks. Multilayer networks are common in many application scenarios, such as modelling of relations among people in a social network or representing the interplay of different molecules in a cell or different cells in the brain. Consequently, the need to introduce algorithms for the comparison of such multilayer networks, i.e., local network alignment, arises. Existing algorithms for LNA do not perform well on multilayer networks since they cannot consider inter-layer edges. Thus, we propose local alignment of multilayer networks (MultiLoAl), a novel algorithm for the local alignment of multilayer networks. We define the local alignment of multilayer networks and propose a heuristic for solving it. We present an extensive assessment indicating the strength of the algorithm. Furthermore, we implemented a synthetic multilayer network generator to build the data for the algorithm's evaluation.

PCN-Miner: an open-source extensible tool for the analysis of Protein Contact Networks

MOTIVATION

Protein Contact Network (PCN) is a powerful method for analysing the structure and function of proteins, with a specific focus on disclosing the molecular features of allosteric regulation through the discovery of modular substructures. The importance of PCN analysis has been shown in many contexts, such as the analysis of SARS-CoV-2 Spike protein and its complexes with the Angiotensin Converting Enzyme 2 (ACE2) human receptors. Even if there exist many software tools implementing such methods, there is a growing need for the introduction of tools integrating existing approaches.

RESULTS

We present PCN-Miner, a software tool implemented in the Python programming language, able to (i) import protein structures from the Protein Data Bank; (ii) generate the corresponding PCN; (iii) model, analyse and visualize PCNs and related protein structures by using a set of known algorithms and metrics. The PCN-Miner can cover a large set of applications: from clustering to embedding and subsequent analysis.

AVAILABILITY AND IMPLEMENTATION

The PCN-Miner tool is freely available at the following GitHub repository: https://github.com/hguzzi/ProteinContactNetworks. It is also available in the Python Package Index (PyPI) repository.

Structural genetics of circulating variants affecting the SARS-CoV-2 spike/human ACE2 complex

SARS-CoV-2 entry in human cells is mediated by the interaction between the viral Spike protein and the human ACE2 receptor. This mechanism evolved from the ancestor bat coronavirus and is currently one of the main targets for antiviral strategies. However, there currently exist several Spike protein variants in the SARS-CoV-2 population as the result of mutations, and it is unclear if these variants may exert a specific effect on the affinity with ACE2 which, in turn, is also characterized by multiple alleles in the human population. In the current study, the GBPM analysis, originally developed for highlighting host-guest interaction features, has been applied to define the key amino acids responsible for the Spike/ACE2 molecular recognition, using four different crystallographic structures. Then, we intersected these structural results with the current mutational status, based on more than 295,000 sequenced cases, in the SARS-CoV-2 population. We identified several Spike mutations interacting with ACE2 and mutated in at least 20 distinct patients: S477N, N439K, N501Y, Y453F, E484K, K417N, S477I and G476S. Among these, mutation N501Y in particular is one of the events characterizing SARS-CoV-2 lineage B.1.1.7, which has recently risen in frequency in Europe. We also identified five ACE2 rare variants that may affect interaction with Spike and susceptibility to infection: S19P, E37K, M82I, E329G and G352V.Communicated by Ramaswamy H. Sarma.

Metabolic Syndrome Is Associated With Impaired Insulin-Stimulated Myocardial Glucose Metabolic Rate in Individuals With Type 2 Diabetes: A Cardiac Dynamic 18F-FDG-PET Study

Metabolic syndrome is a condition characterized by a clustering of metabolic abnormalities associated with an increased risk of type 2 diabetes and cardiovascular disease. An impaired insulin-stimulated myocardial glucose metabolism has been shown to be a risk factor for the development of cardiovascular disease in patients with type 2 diabetes. Whether cardiac insulin resistance occurs in subjects with metabolic syndrome remains uncertain. To investigate this issue, we evaluated myocardial glucose metabolic rate using cardiac dynamic 18F-FDG-PET combined with euglycemic-hyperinsulinemic clamp in three groups: a group of normal glucose tolerant individuals without metabolic syndrome (n = 10), a group of individuals with type 2 diabetes and metabolic syndrome (n = 19), and a group of subjects with type 2 diabetes without metabolic syndrome (n = 6). After adjusting for age and gender, individuals with type 2 diabetes and metabolic syndrome exhibited a significant reduction in insulin-stimulated myocardial glucose metabolic rate (10.5 ± 9.04 μmol/min/100 g) as compared with both control subjects (32.9 ± 9.7 μmol/min/100 g; P &lt; 0.0001) and subjects with type 2 diabetes without metabolic syndrome (25.15 ± 4.92 μmol/min/100 g; P = 0.01). Conversely, as compared with control subjects (13.01 ± 8.53 mg/min x Kg FFM), both diabetic individuals with metabolic syndrome (3.06 ± 1.7 mg/min × Kg FFM, P = 0.008) and those without metabolic syndrome (2.91 ± 1.54 mg/min × Kg FFM, P = 0.01) exhibited a significant reduction in whole-body insulin-stimulated glucose disposal, while no difference was observed between the 2 groups of subjects with type 2 diabetes with or without metabolic syndrome. Univariate correlations showed that myocardial glucose metabolism was positively correlated with insulin-stimulated glucose disposal (r = 0.488, P = 0.003), and negatively correlated with the presence of metabolic syndrome (r = −0.743, P &lt; 0.0001) and with its individual components. In conclusion, our data suggest that an impaired myocardial glucose metabolism may represent an early cardio-metabolic defect in individuals with the coexistence of type 2 diabetes and metabolic syndrome, regardless of whole-body insulin resistance.

Disease spreading modeling and analysis: a survey

MOTIVATION

The control of the diffusion of diseases is a critical subject of a broad research area, which involves both clinical and political aspects. It makes wide use of computational tools, such as ordinary differential equations, stochastic simulation frameworks and graph theory, and interaction data, from molecular to social granularity levels, to model the ways diseases arise and spread. The coronavirus disease 2019 (COVID-19) is a perfect testbench example to show how these models may help avoid severe lockdown by suggesting, for instance, the best strategies of vaccine prioritization.

RESULTS

Here, we focus on and discuss some graph-based epidemiological models and show how their use may significantly improve the disease spreading control. We offer some examples related to the recent COVID-19 pandemic and discuss how to generalize them to other diseases.

Beyond COVID-19 Pandemic: Topology-aware optimization of vaccination strategy for minimizing virus spreading

•

Paper discusses the relevance of the adoption of ad-hoc vaccination strategies.

•

Paper shows how to evaluate the impact of different vaccination strategy by considering network-based models.

•

Tailored interventions, e.g., vaccination, applied on central nodes of these networks may efficiently stop the propagation of an infection.

•

The way node "centrality" is defined is the key to curb infection spreading.

The mitigation of an infectious disease spreading has recently gained considerable attention from the research community. It may be obtained by adopting sanitary measurements (e.g., vaccination, wearing masks), social rules (e.g., social distancing), together with an extensive vaccination campaign. Vaccination is currently the primary way for mitigating the Coronavirus Disease (COVID-19) outbreak without severe lockdown. Its effectiveness also depends on the number and timeliness of administrations and thus demands strict prioritization criteria. Almost all countries have prioritized similar classes of exposed workers: healthcare professionals and the elderly, obtaining to maximize the survival of patients and years of life saved. Nevertheless, the virus is currently spreading at high rates, and any prioritization criterion so far adopted did not account for the structural organization of the contact networks.

We reckon that a network where nodes are people while the edges represent their social contacts may efficiently model the virus’s spreading. It is known that tailored interventions (e.g., vaccination) on central nodes may efficiently stop the propagation, thereby eliminating the “bridge edges.” We then introduce such a model and consider both synthetic and real datasets. We present the benefits of a topology-aware versus an age-based vaccination strategy to mitigate the spreading of the virus. The code is available at https://github.com/mazzalab/playgrounds.

Modeling multi-scale data via a network of networks

MOTIVATION

Prediction of node and graph labels are prominent network science tasks. Data analyzed in these tasks are sometimes related: entities represented by nodes in a higher-level (higher scale) network can themselves be modeled as networks at a lower level. We argue that systems involving such entities should be integrated with a 'network of networks' (NoNs) representation. Then, we ask whether entity label prediction using multi-level NoN data via our proposed approaches is more accurate than using each of single-level node and graph data alone, i.e. than traditional node label prediction on the higher-level network and graph label prediction on the lower-level networks. To obtain data, we develop the first synthetic NoN generator and construct a real biological NoN. We evaluate accuracy of considered approaches when predicting artificial labels from the synthetic NoNs and proteins' functions from the biological NoN.

RESULTS

For the synthetic NoNs, our NoN approaches outperform or are as good as node- and network-level ones depending on the NoN properties. For the biological NoN, our NoN approaches outperform the single-level approaches for just under half of the protein functions, and for 30% of the functions, only our NoN approaches make meaningful predictions, while node- and network-level ones achieve random accuracy. So, NoN-based data integration is important.

AVAILABILITY AND IMPLEMENTATION

The software and data are available at https://nd.edu/~cone/NoNs.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Using dual-network-analyser for communities detecting in dual networks

BACKGROUND

Representations of the relationships among data using networks are widely used in several research fields such as computational biology, medical informatics and social network mining. Recently, complex networks have been introduced to better capture the insights of the modelled scenarios. Among others, dual networks (DNs) consist of mapping information as pairs of networks containing the same set of nodes but with different edges: one, called physical network, has unweighted edges, while the other, called conceptual network, has weighted edges.

RESULTS

We focus on DNs and we propose a tool to find common subgraphs (aka communities) in DNs with particular properties. The tool, called Dual-Network-Analyser, is based on the identification of communities that induce optimal modular subgraphs in the conceptual network and connected subgraphs in the physical one. It includes the Louvain algorithm applied to the considered case. The Dual-Network-Analyser can be used to study DNs, to find common modular communities. We report results on using the tool to identify communities on synthetic DNs as well as real cases in social networks and biological data.

CONCLUSION

The proposed method has been tested by using synthetic and biological networks. Results demonstrate that it is well able to detect meaningful information from DNs.

Analyzing host-viral interactome of SARS-CoV-2 for identifying vulnerable host proteins during COVID-19 pathogenesis

The development of therapeutic targets for COVID-19 relies on understanding the molecular mechanism of pathogenesis. Identifying genes or proteins involved in the infection mechanism is the key to shedding light on the complex molecular mechanisms. The combined effort of many laboratories distributed throughout the world has produced protein and genetic interactions. We integrated available results and obtained a host protein-protein interaction network composed of 1432 human proteins. Next, we performed network centrality analysis to identify critical proteins in the derived network. Finally, we performed a functional enrichment analysis of central proteins. We observed that the identified proteins are primarily associated with several crucial pathways, including cellular process, signaling transduction, neurodegenerative diseases. We focused on the proteins that are involved in human respiratory tract diseases. We highlighted many potential therapeutic targets, including RBX1, HSPA5, ITCH, RAB7A, RAB5A, RAB8A, PSMC5, CAPZB, CANX, IGF2R, and HSPA1A, which are central and also associated with multiple diseases.

Exploiting the molecular basis of age and gender differences in outcomes of SARS-CoV-2 infections

Motivation: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (coronavirus disease, 2019; COVID-19) is associated with adverse outcomes in patients. It has been observed that lethality seems to be related to the age of patients. While ageing has been extensively demonstrated to be accompanied by some modifications at the gene expression level, a possible link with COVID-19 manifestation still need to be investigated at the molecular level. Objectives: This study aims to shed out light on a possible link between the increased COVID-19 lethality and the molecular changes that occur in elderly people. Methods: We considered public datasets of ageing-related genes and their expression at the tissue level. We selected human proteins interacting with viral ones that are known to be related to the ageing process. Finally, we investigated changes in the expression level of coding genes at the tissue, gender and age level. Results: We observed a significant intersection between some SARS-CoV-2 interactors and ageing-related genes, suggesting that those genes are particularly affected by COVID-19 infection. Our analysis evidenced that virus infection particularly involves ageing molecular mechanisms centred around proteins EEF2, NPM1, HMGA1, HMGA2, APEX1, CHEK1, PRKDC, and GPX4. We found that HMGA1 and NPM1 have different expressions in the lung of males, while HMGA1, APEX1, CHEK1, EEF2, and NPM1 present changes in expression in males due to ageing effects. Conclusion: Our study generated a mechanistic framework to clarify the correlation between COVID-19 incidence in elderly patients and molecular mechanisms of ageing. We also provide testable hypotheses for future investigation and pharmacological solutions tailored to specific age ranges.

BioPAX-Parser: parsing and enrichment analysis of BioPAX pathways

SUMMARY

Biological pathways are fundamental for learning about healthy and disease states. Many existing formats support automatic software analysis of biological pathways, e.g. BioPAX (Biological Pathway Exchange). Although some algorithms are available as web application or stand-alone tools, no general graphical application for the parsing of BioPAX pathway data exists. Also, very few tools can perform pathway enrichment analysis (PEA) using pathway encoded in the BioPAX format. To fill this gap, we introduce BiP (BioPAX-Parser), an automatic and graphical software tool aimed at performing the parsing and accessing of BioPAX pathway data, along with PEA by using information coming from pathways encoded in BioPAX.

AVAILABILITY AND IMPLEMENTATION

BiP is freely available for academic and non-profit organizations at https://gitlab.com/giuseppeagapito/bip under the LGPL 2.1, the GNU Lesser General Public License.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Non-coding RNAs in cancer: platforms and strategies for investigating the genomic "dark matter"

The discovery of the role of non-coding RNAs (ncRNAs) in the onset and progression of malignancies is a promising frontier of cancer genetics. It is clear that ncRNAs are candidates for therapeutic intervention, since they may act as biomarkers or key regulators of cancer gene network. Recently, profiling and sequencing of ncRNAs disclosed deep deregulation in human cancers mostly due to aberrant mechanisms of ncRNAs biogenesis, such as amplification, deletion, abnormal epigenetic or transcriptional regulation. Although dysregulated ncRNAs may promote hallmarks of cancer as oncogenes or antagonize them as tumor suppressors, the mechanisms behind these events remain to be clarified. The development of new bioinformatic tools as well as novel molecular technologies is a challenging opportunity to disclose the role of the "dark matter" of the genome. In this review, we focus on currently available platforms, computational analyses and experimental strategies to investigate ncRNAs in cancer. We highlight the differences among experimental approaches aimed to dissect miRNAs and lncRNAs, which are the most studied ncRNAs. These two classes indeed need different investigation taking into account their intrinsic characteristics, such as length, structures and also the interacting molecules. Finally, we discuss the relevance of ncRNAs in clinical practice by considering promises and challenges behind the bench to bedside translation.

DMETTM Genotyping: Tools for Biomarkers Discovery in the Era of Precision Medicine

The knowledge of genetic variants in genes involved in drug metabolism may be translated into reduction of adverse drug reactions, increase of efficacy, healthcare outcomes improvement and economic benefits. Many high-throughput tools are available for the genotyping of Single Nucleotide Polymorphisms (SNPs) known to be related to drugs and xenobiotics metabolism. DMET^TM platform represents an example of SNPs panel to discover biomarkers correlated to efficacy or toxicity in common and rare diseases. The difficulty in analyzing the mole of information generated by DMET^TM platform led to the development and implementation of algorithms and tools for statistical and data mining analysis. These softwares allow efficient handling of the omics data to validate the explorative SNPs identified by DMET assay and to correlate them with drug efficacy, toxicity and/or cancer susceptibility. In this review we present a suite of bioinformatic frameworks for the preprocessing and analysis of DMET-SNPs data. In particular, we introduce a workflow that uses the GenoMetric Query Language, a high-level query language specifically designed for genomics, able to query public datasets (such as ENCODE, TCGA, GENCODE annotation dataset, etc.) as well as to combine them with private datasets (e.g., output from Affymetrix® DMET^TM Platform).

L-HetNetAligner: A novel algorithm for Local Alignment of Heterogeneous Biological Networks

Networks are largely used for modelling and analysing a wide range of biological data. As a consequence, many different research efforts have resulted in the introduction of a large number of algorithms for analysis and comparison of networks. Many of these algorithms can deal with networks with a single class of nodes and edges, also referred to as homogeneous networks. Recently, many different approaches tried to integrate into a single model the interplay of different molecules. A possible formalism to model such a scenario comes from node/edge coloured networks (also known as heterogeneous networks) implemented as node/ edge-coloured graphs. Therefore, the need for the introduction of algorithms able to compare heterogeneous networks arises. We here focus on the local comparison of heterogeneous networks, and we formulate it as a network alignment problem. To the best of our knowledge, the local alignment of heterogeneous networks has not been explored in the past. We here propose L-HetNetAligner a novel algorithm that receives as input two heterogeneous networks (node-coloured graphs) and builds a local alignment of them. We also implemented and tested our algorithm. Our results confirm that our method builds high-quality alignments. The following website *contains Supplementary File 1 material and the code.

Survey of local and global biological network alignment: the need to reconcile the two sides of the same coin

Analogous to genomic sequence alignment that allows for across-species transfer of biological knowledge between conserved sequence regions, biological network alignment can be used to guide the knowledge transfer between conserved regions of molecular networks of different species. Hence, biological network alignment can be used to redefine the traditional notion of a sequence-based homology to a new notion of network-based homology. Analogous to genomic sequence alignment, there exist local and global biological network alignments. Here, we survey prominent and recent computational approaches of each network alignment type and discuss their (dis)advantages. Then, as it was recently shown that the two approach types are complementary, in the sense that they capture different slices of cellular functioning, we discuss the need to reconcile the two network alignment types and present a recent first step in this direction. We conclude with some open research problems on this topic and comment on the usefulness of network alignment in other domains besides computational biology.

From Single Level Analysis to Multi-Omics Integrative Approaches: A Powerful Strategy towards the Precision Oncology

Integration of multi-omics data from different molecular levels with clinical data, as well as epidemiologic risk factors, represents an accurate and promising methodology to understand the complexity of biological systems of human diseases, including cancer. By the extensive use of novel technologic platforms, a large number of multidimensional data can be derived from analysis of health and disease systems. Comprehensive analysis of multi-omics data in an integrated framework, which includes cumulative effects in the context of biological pathways, is therefore eagerly awaited. This strategy could allow the identification of pathway-addiction of cancer cells that may be amenable to therapeutic intervention. However, translation into clinical settings requires an optimized integration of omics data with clinical vision to fully exploit precision cancer medicine. We will discuss the available technical approach and more recent developments in the specific field.

An Innovative Framework for Bioimage Annotation and Studies

The collection and analysis of clinical data are needed to investigate diseases and to define medical protocols and treatments. Bioimages, medical annotations and patient history are clinical data acquired and studied to perform a correct diagnosis and to propose an appropriate therapy. Currently, hospital departments manage these data using legacy systems which do not often allow data integration among different departments or health structures. Thus, in many cases clinical information sharing and exchange are difficult to implement. This is also the case for biomedical images for which data integration or data overlapping is usually not available. Image annotations and comparison can be crucial for physicians in many case studies. In this paper, a general purpose framework for bioimage management and annotations is proposed. Moreover, a simple-to-use information system has been developed to integrate clinical and diagnosis codes. The framework allows physicians (1) to integrate DICOM images from different platforms and (2) to report notes and highlights directly on images, thus offering, among the others, to query and compare similar clinical cases. This contribution is the result of a framework aimed to support oncologists in managing DICOM images and clinical data from different departments. Data integration is performed using a here-proposed XML-based module also utilized to trace temporal changes in image annotations.

A Parallel Software Pipeline for DMET Microarray Genotyping Data Analysis

Personalized medicine is an aspect of the P4 medicine (predictive, preventive, personalized and participatory) based precisely on the customization of all medical characters of each subject. In personalized medicine, the development of medical treatments and drugs is tailored to the individual characteristics and needs of each subject, according to the study of diseases at different scales from genotype to phenotype scale. To make concrete the goal of personalized medicine, it is necessary to employ high-throughput methodologies such as Next Generation Sequencing (NGS), Genome-Wide Association Studies (GWAS), Mass Spectrometry or Microarrays, that are able to investigate a single disease from a broader perspective. A side effect of high-throughput methodologies is the massive amount of data produced for each single experiment, that poses several challenges (e.g., high execution time and required memory) to bioinformatic software. Thus a main requirement of modern bioinformatic softwares, is the use of good software engineering methods and efficient programming techniques, able to face those challenges, that include the use of parallel programming and efficient and compact data structures. This paper presents the design and the experimentation of a comprehensive software pipeline, named microPipe, for the preprocessing, annotation and analysis of microarray-based Single Nucleotide Polymorphism (SNP) genotyping data. A use case in pharmacogenomics is presented. The main advantages of using microPipe are: the reduction of errors that may happen when trying to make data compatible among different tools; the possibility to analyze in parallel huge datasets; the easy annotation and integration of data. microPipe is available under Creative Commons license, and is freely downloadable for academic and not-for-profit institutions.

GLAlign: A Novel Algorithm for Local Network Alignment

Networks are successfully used as a modelling framework in many application domains. For instance, Protein-Protein Interaction Networks (PPINs) model the set of interactions among proteins in a cell. A critical application of network analysis is the comparison among PPINs of different organisms to reveal similarities among the underlying biological processes. Algorithms for comparing networks (also referred to as network aligners) fall into two main classes: global aligners, which aim to compare two networks on a global scale, and local aligners that evidence single sub-regions of similarity among networks. The possibility to improve the performance of the aligners by mixing the two approaches is a growing research area. In our previous work, we started to explore the possibility to use global alignment to improve the local one. We here examine further this possibility by using topological information extracted from global alignment to guide the steps of the local alignment. Therefore, we present GLAlign (Global Local Aligner), a methodology that improves the performances of local network aligners by exploiting a preliminary global alignment. Furthermore, we provide an implementation of GLAlign. As a proof-of-principle, we evaluated the performance of the GLAlign prototype. Results show that GLAlign methodology outperforms the state-of-the-art local alignment algorithms. GLAlign is publicly available for academic use and can be downloaded here: https://sites.google.com/site/globallocalalignment/.

Genetic variants associated with Fabry disease progression despite enzyme replacement therapy

Enzyme replacement therapy (ERT) has been widely used for the treatment of Fabry disease, a rare X-linked recessive disorder due to absent or reduced activity of lysosomal enzyme α-galactosidase A. It is still unclear why some patients under ERT show disease progression typically with renal, cardiovascular and cerebrovascular dysfunctions. Here, we investigated the involvement of drug absorption, distribution, metabolism, and excretion gene variants in response variability to ERT, genotyping 37 patients with the Affymetrix Drug Metabolizing Enzyme and Transporters (DMET) Plus microarray. We found three single nucleotide polymorphisms in human alcohol dehydrogenase (ADH)4 gene (rs1126670, rs1126671, rs2032349) and one in ADH5 gene (rs2602836) associated with disease progression (p < 0.05). Our data provide a basic tool for identification of patient with ERT non-response risk that may represent a framework for personalized treatment of this rare disease.

An extensive assessment of network alignment algorithms for comparison of brain connectomes

Background

Recently the study of the complex system of connections in neural systems, i.e. the connectome, has gained a central role in neurosciences. The modeling and analysis of connectomes are therefore a growing area. Here we focus on the representation of connectomes by using graph theory formalisms. Macroscopic human brain connectomes are usually derived from neuroimages; the analyzed brains are co-registered in the image domain and brought to a common anatomical space. An atlas is then applied in order to define anatomically meaningful regions that will serve as the nodes of the network - this process is referred to as parcellation. The atlas-based parcellations present some known limitations in cases of early brain development and abnormal anatomy. Consequently, it has been recently proposed to perform atlas-free random brain parcellation into nodes and align brains in the network space instead of the anatomical image space, as a way to deal with the unknown correspondences of the parcels. Such process requires modeling of the brain using graph theory and the subsequent comparison of the structure of graphs. The latter step may be modeled as a network alignment (NA) problem.

Results

In this work, we first define the problem formally, then we test six existing state of the art of network aligners on diffusion MRI-derived brain networks. We compare the performances of algorithms by assessing six topological measures. We also evaluated the robustness of algorithms to alterations of the dataset.

Conclusion

The results confirm that NA algorithms may be applied in cases of atlas-free parcellation for a fully network-driven comparison of connectomes. The analysis shows MAGNA++ is the best global alignment algorithm. The paper presented a new analysis methodology that uses network alignment for validating atlas-free parcellation brain connectomes. The methodology has been experimented on several brain datasets.

Genetic variants associated with gastrointestinal symptoms in Fabry disease

Gastrointestinal symptoms (GIS) are often among the earliest presenting events in Fabry disease (FD), an X-linked lysosomal disorder caused by the deficiency of α-galactosidase A. Despite recent advances in clinical and molecular characterization of FD, the pathophysiology of the GIS is still poorly understood. To shed light either on differential clinical presentation or on intervariability of GIS in FD, we genotyped 1936 genetic markers across 231 genes that encode for drug-metabolizing enzymes and drug transport proteins in 49 FD patients, using the DMET Plus platform. All nine single nucleotide polymorphisms (SNPs) mapped within four genes showed statistically significant differences in genotype frequencies between FD patients who experienced GIS and patients without GIS: ABCB11 (odd ratio (OR) = 18.07, P = 0,0019; OR = 8.21, P = 0,0083; OR=8.21, P = 0,0083; OR = 8.21, P = 0,0083),SLCO1B1 (OR = 9.23, P = 0,0065; OR = 5.08, P = 0,0289; OR = 8.21, P = 0,0083), NR1I3 (OR = 5.40, P = 0,0191) and ABCC5 (OR = 14.44, P = 0,0060). This is the first study that investigates the relationships between genetic heterogeneity in drug absorption, distribution, metabolism and excretion (ADME) related genes and GIS in FD. Our findings provide a novel genetic variant framework which warrants further investigation for precision medicine in FD.

OSAnalyzer: A Bioinformatics Tool for the Analysis of Gene Polymorphisms Enriched with Clinical Outcomes

BACKGROUND

The identification of biomarkers for the estimation of cancer patients' survival is a crucial problem in modern oncology. Recently, the Affymetrix DMET (Drug Metabolizing Enzymes and Transporters) microarray platform has offered the possibility to determine the ADME (absorption, distribution, metabolism, and excretion) gene variants of a patient and to correlate them with drug-dependent adverse events. Therefore, the analysis of survival distribution of patients starting from their profile obtained using DMET data may reveal important information to clinicians about possible correlations among drug response, survival rate, and gene variants.

METHODS

In order to provide support to this analysis we developed OSAnalyzer, a software tool able to compute the overall survival (OS) and progression-free survival (PFS) of cancer patients and evaluate their association with ADME gene variants.

RESULTS

The tool is able to perform an automatic analysis of DMET data enriched with survival events. Moreover, results are ranked according to statistical significance obtained by comparing the area under the curves that is computed by using the log-rank test, allowing a quick and easy analysis and visualization of high-throughput data.

CONCLUSIONS

Finally, we present a case study to highlight the usefulness of OSAnalyzer when analyzing a large cohort of patients.

Extracting Cross-Ontology Weighted Association Rules from Gene Ontology Annotations

Gene Ontology (GO) is a structured repository of concepts (GO Terms) that are associated to one or more gene products through a process referred to as annotation. The analysis of annotated data is an important opportunity for bioinformatics. There are different approaches of analysis, among those, the use of association rules (AR) which provides useful knowledge, discovering biologically relevant associations between terms of GO, not previously known. In a previous work, we introduced GO-WAR (Gene Ontology-based Weighted Association Rules), a methodology for extracting weighted association rules from ontology-based annotated datasets. We here adapt the GO-WAR algorithm to mine cross-ontology association rules, i.e., rules that involve GO terms present in the three sub-ontologies of GO. We conduct a deep performance evaluation of GO-WAR by mining publicly available GO annotated datasets, showing how GO-WAR outperforms current state of the art approaches.

Biological Network Inference from Microarray Data, Current Solutions, and Assessments

Currently in bioinformatics and systems biology there is a growing interest for the analysis of associations among biological molecules at a network level. A main research in this area is represented by the inference of biological networks from experimental data. Biological network inference aims to reconstruct network of interactions (or associations) among biological molecules (e.g., genes or proteins) starting from experimental observations. The current scenario is characterized by a growing number of algorithms for the inference, while few attention has been posed on the determination of fair assessments and comparisons. Current assessments are usually based on the comparison of the algorithms using reference networks or gold standard datasets. Here we survey some selected inference algorithms and we compare current assessments. We also present a systematic listing of freely available inference and assessment tools for easy reference. Finally we outline some possible future directions of research, such as the use of a prior knowledge into the assessment process.

Non-alcoholic fatty liver disease severity, central fat mass and adinopectin: a close relationship

AIM

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in the general population. Overweight is a common condition in patients with NAFLD, and body composition (BC) assessment is useful to evaluate nutritional status and the efficacy of nutritional strategies. A valid tool for assessing BC is dual-energy X-ray absorptiometry (DXA). Adiponectin has been shown to be relevant to the pathogenesis of NAFLD. The aim of this observational study is to define the relationship between the severity of NAFLD, the central fat mass evaluated by DXA, and the circulating levels of adiponectin.

METHODS

The study was carried out in 31 overweight patients. The degree of liver steatosis was evaluated by ultrasound (US) examination. Anthropometric parameters were measured according to standard methods. Fasting glucose and insulin level were used also to calculate insulin resistance (IR), according to the homeostasis model assessment-insulin resistance (HOMA-IR). The enzyme-linked immunosorbent assay technique was performed to dose fasting serum levels of adiponectin.

RESULTS

NAFLD progression was significantly associated with increased central fat (p<0.05). Using DXA, we quantified the regional distribution of adipose tissue and found the expected association between central fat and the US severity of NAFLD. Serum levels of adiponectin, were inversely related to NAFLD progression (p<0.05).

CONCLUSION

BC evaluated by anthropometry and DXA, may be used as indicator of NAFLD severity in overweight patients. The evaluation of BC in clinical practice, can improve the nutritional strategies and follow-up. In the clinical setting adiponectin may represent a potential marker for the staging of NAFLD.

Using GO-WAR for mining cross-ontology weighted association rules

The Gene Ontology (GO) is a structured repository of concepts (GO terms) that are associated to one or more gene products. The process of association is referred to as annotation. The relevance and the specificity of both GO terms and annotations are evaluated by a measure defined as information content (IC). The analysis of annotated data is thus an important challenge for bioinformatics. There exist different approaches of analysis. From those, the use of association rules (AR) may provide useful knowledge, and it has been used in some applications, e.g. improving the quality of annotations. Nevertheless classical association rules algorithms do not take into account the source of annotation nor the importance yielding to the generation of candidate rules with low IC. This paper presents GO-WAR (Gene Ontology-based Weighted Association Rules) a methodology for extracting weighted association rules. GO-WAR can extract association rules with a high level of IC without loss of support and confidence from a dataset of annotated data. A case study on using of GO-WAR on publicly available GO annotation datasets is used to demonstrate that our method outperforms current state of the art approaches.

Vestibular disorders in euthyroid patients with Hashimoto's thyroiditis: role of thyroid autoimmunity

INTRODUCTION

A relationship between vestibular disorders and thyroid autoimmunity independently from thyroid function has been postulated.

AIM

To shed more light on the actual relationship between vestibular lesions and Hashimoto's thyroiditis (HT) regardless of thyroid function.

METHODS

Forty-seven patients with HT (89·4% F; aged 48·3 ± 12·7 years), 21 with multinodular goitre (MNG; 57·1% F; 54·1 ± 9·8 years) and 30 healthy volunteers (56·7% F; 50·7 ± 13·9 years) were enrolled. Inclusion criteria were the presence of normal thyroid function tests and no clinical history of vestibular dysfunction. Each subject was submitted to complete vestibular evaluation [Caloric Test, Vestibular evoked myogenic potentials (VEMPs), Head Shaking Test (HST)].

RESULTS

52·2% of HT patients showed an alteration of VEMPs and 44·7% of caloric test (P < 0·0001 for both). None of the MNG patients showed any vestibular alteration, while one healthy control showed an altered caloric test. A correlation was found between vestibular alterations of HT patients and the degree of serum TPOAb level, not affected by age and serum TSH value. By logistic regression analysis, the absence of thyroid autoimmunity significantly reduced the risk of vestibular alterations: HR 0.19 (95%CI: 0·003-0.25, P = 0·0004) for caloric test; HR 0·07 (95%CI: 0·02-0·425, P < 0·0001) for VEMPs; and HR 0·22 (95%CI: 0·06-0·7, P = 0·01) for HST.

CONCLUSION

In euthyroid HT patients, a significant relationship between subclinical vestibular damage and the degree of TPOAb titre was documented. This finding suggests that circulating antithyroid autoantibodies may represent a risk factor for developing vestibular dysfunction. An accurate vestibular evaluation of HT patients with or without symptoms is therefore warranted.

Serum adipokine levels in overweight patients and their relationship with non-alcoholic fatty liver disease

AIM

Non-alcoholic fatty liver disease (NAFLD) is a relevant public health matter in Western countries. The pathogenetic link between visceral fat, insulin resistance (IR) and NAFLD has been reported in literature. However, there are contradictions on the changes of adipokine levels in serum related to the presence of NAFLD. The aim of the present study was to evaluate the serum concentrations of a selected set of adipokines, that is, adiponectin, leptin, resistin and the pro-inflammatory cytokine interleukin-6 (IL-6) in overweight patients, and to clarify their relationship with NAFLD.

METHODS

Fasting serum levels of adipokines were determined in 42 consecutive overweight patients and in 25 lean controls. The degree of ultrasound (US) liver steatosis was graded according to the Hamaguchi score.

RESULTS

Liver steatosis was detected in 33 patients (78%) by US examination. Twelve patients with elevated transaminases levels showed significantly higher values of IR, leptin and resistin levels (P<0.05). Patients with steatosis presented a significantly higher leptin and a lower adiponectin levels (P<0.05) than controls. A significant inverse correlation was found between US steatosis progression and adiponectin and resistin levels (p<0.05). Considering the multiple logistic regression, adiponectin and leptin were good predictors to detect the presence of steatosis (p<0.05).

CONCLUSION

Our data support the concept that adipokine level changes are closely linked with IR. In addition, serum adiponectin and leptin levels may be used as diagnostic markers to determine the presence of NAFLD in overweight patients.

Improving the Robustness of Local Network Alignment: Design and Extensive Assessment of a Markov Clustering-Based Approach

The analysis of protein behavior at the network level had been applied to elucidate the mechanisms of protein interaction that are similar in different species. Published network alignment algorithms proved to be able to recapitulate known conserved modules and protein complexes, and infer new conserved interactions confirmed by wet lab experiments. In the meantime, however, a plethora of continuously evolving protein-protein interaction (PPI) data sets have been developed, each featuring different levels of completeness and reliability. For instance, algorithms performance may vary significantly when changing the data set used in their assessment. Moreover, existing papers did not deeply investigate the robustness of alignment algorithms. For instance, some algorithms performances vary significantly when changing the data set used in their assessment. In this work, we design an extensive assessment of current algorithms discussing the robustness of the results on the basis of input networks. We also present AlignMCL, a local network alignment algorithm based on an improved model of alignment graph and Markov Clustering. AlignMCL performs better than other state-of-the-art local alignment algorithms over different updated data sets. In addition, AlignMCL features high levels of robustness, producing similar results regardless the selected data set.