Search

Abstract Lynch syndrome is rarely associated with rectal cancer (RC) and thus, metachronous RC has been scarcely investigated. This study aimed to analyze the mucosal immune microenvironment in sporadic and metachronous RC. We analyzed the mucosal immune microenvironment in the 25 metachronous RCs present in the IMMUNOREACT 1 and 2 multicentre observational studies (624 patients). A panel of immune markers was retrospectively investigated at immunohistochemistry: CD3, CD4, CD8, CD8b, Tbet, FoxP3, PD-L1, MSH6, and PMS2 and CD80. Single-cell suspensions were subjected to flow-cytometry to determine the proportion of epithelial cells (pan-cytokeratin) acting as antigen-presenting cells (expressing CD80, CD86, HLA-ABC) and the proportion of activated CD8 + T cells (CD8 + positive for CD28, CD38), inhibitory T cells (CD3 + CTLA-4+) of activated CD4 + T helper cells (CD4 + CD25+) and activated T regulatory cells (CD4 + CD25 + FoxP3+). No mismatch repair gene deficiencies were observed in the patients. The previous history of colorectal adenoma was significantly more frequent in metachronous RC. In healthy epithelial cells, HLA-ABC expression was significantly higher in patients with metachronous RC. In therapy-naïve metachronous RC patients, a significantly lower level of circulating lymphocytes and CD3 + T-cell infiltration in the healthy mucosa surrounding the RC was observed compared to patients with non-metachronous cancer. Our study supports the hypothesis that metachronous RC can occur in a cancerization field in patients with weak systemic and local immune systems. The peculiar site of RC makes the mismatch-repair genes deficiency in metachronous cancer onset less relevant.

Abstract Background Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease that results in death within a short time span (3-5 years). One of the major challenges in treating ALS is its highly heterogeneous disease progression and the lack of effective prognostic tools to forecast it. The main aim of this study was, then, to test the feasibility of predicting relevant clinical outcomes that characterize the progression of ALS with a two-year prediction horizon via artificial intelligence techniques using routine visits data. Methods Three classification problems were considered: predicting death (binary problem), predicting death or percutaneous endoscopic gastrostomy (PEG) (multiclass problem), and predicting death or non-invasive ventilation (NIV) (multiclass problem). Two supervised learning models, a logistic regression (LR) and a deep learning multilayer perceptron (MLP), were trained ensuring technical robustness and reproducibility. Moreover, to provide insights into model explainability and result interpretability, model coefficients for LR and Shapley values for both LR and MLP were considered to characterize the relationship between each variable and the outcome. Results On the one hand, predicting death was successful as both models yielded F1 scores and accuracy well above 0.7. The model explainability analysis performed for this outcome allowed for the understanding of how different methodological approaches consider the input variables when performing the prediction. On the other hand, predicting death alongside PEG or NIV proved to be much more challenging (F1 scores and accuracy in the 0.4-0.6 interval). Conclusions In conclusion, predicting death due to ALS proved to be feasible. However, predicting PEG or NIV in a multiclass fashion proved to be unfeasible with these data, regardless of the complexity of the methodological approach. The observed results suggest a potential ceiling on the amount of information extractable from the database, e.g., due to the intrinsic difficulty of the prediction tasks at hand, or to the absence of crucial predictors that are, however, not currently collected during routine practice.

Abstract Comorbidity is widespread in the ageing population, implying multiple and complex medical needs for individuals and a public health burden. Determining risk factors and predicting comorbidity development can help identify at-risk subjects and design prevention strategies. Using socio-demographic and clinical data from approximately 11,000 subjects monitored over 11 years in the English Longitudinal Study of Ageing, we develop a dynamic Bayesian network (DBN) to model the onset and interaction of three cardio-metabolic comorbidities, namely type 2 diabetes (T2D), hypertension, and heart problems. The DBN allows us to identify risk factors for developing each morbidity, simulate ageing progression over time, and stratify the population based on the risk of outcome occurrence. By applying hierarchical agglomerative clustering to the simulated, dynamic risk of experiencing morbidities, we identified patients with similar risk patterns and the variables contributing to their discrimination. The network reveals a direct joint effect of biomarkers and lifestyle on outcomes over time, such as the impact of fasting glucose, HbA1c, and BMI on T2D development. Mediated cross-relationships between comorbidities also emerge, showcasing the interconnected nature of these health issues. The model presents good calibration and discrimination ability, particularly in predicting the onset of T2D (iAUC-ROC = 0.828, iAUC-PR = 0.294) and survival (iAUC-ROC = 0.827, iAUC-PR = 0.311). Stratification analysis unveils two distinct clusters for all comorbidities, effectively discriminated by variables like HbA1c for T2D and age at baseline for heart problems. The developed DBN constitutes an effective, highly-explainable predictive risk tool for simulating and stratifying the dynamic risk of developing cardio-metabolic comorbidities. Its use could help identify the effects of risk factors and develop health policies that prevent the occurrence of comorbidities.

Abstract Writing notes is the most widespread method to report clinical events. Therefore, most of the information about the disease history of a patient remains locked behind free-form text. Natural language processing (NLP) provides a solution to automatically transform free-form text into structured data. In the present work, electronic healthcare records data of patients with diabetes were used to develop deep-learning based NLP models to automatically identify, within free-form text describing routine visits, the occurrence of hospitalisations related to cardiovascular disease (CVDs), an outcome of diabetes. Four possible time windows of increasing level of expected difficulty were considered: infinite, 24 months, 12 months, and 6 months. Model performance was evaluated by means of the area under the precision recall curve, as well as precision, recall, and F1-score after thresholding. Results showed that the proposed NLP approach was successful for both the infinite and 24-month windows, while, as expected, performance deteriorated with shorter time windows. Possible clinical applications of tools based on the proposed NLP approach include the retrospective filling of medical records with respect to a patient’s CVD history for epidemiological and research purposes as well as for clinical decision making.

Abstract Background The ageing global population presents significant public health challenges, especially in relation to the subjective wellbeing of the elderly. In this study, our aim was to investigate the potential for developing a model to forecast the two-year variation of the perceived wellbeing of individuals aged over 50. We also aimed to identify the variables that predict changes in subjective wellbeing, as measured by the CASP-12 scale, over a two-year period. Methods Data from the European SHARE project were used, specifically the demographic, health, social and financial variables of 9422 subjects. The subjective wellbeing was measured through the CASP-12 scale. The study outcome was defined as binary, i.e., worsening/not worsening of the variation of CASP-12 in 2 years. Logistic regression, logistic regression with LASSO regularisation, and random forest were considered candidate models. Performance was assessed in terms of accuracy in correctly predicting the outcome, Area Under the Curve (AUC), and F1 score. Results The best-performing model was the random forest, achieving an accuracy of 65%, AUC = 0.659, and F1 = 0.710. All models proved to be able to generalise both across subjects and over time. The most predictive variables were the CASP-12 score at baseline, the presence of depression and financial difficulties. Conclusions While we identify the random forest model as the more suitable, given the similarity of performance, the models based on logistic regression or on logistic regression with LASSO regularisation are also possible options.

Abstract The outbreak of Coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, forced us to face a pandemic with unprecedented social, economic, and public health consequences. Several nations have launched campaigns to immunize millions of people using various vaccines to prevent infections. Meanwhile, therapeutic approaches and discoveries continuously arise; however, identifying infected patients that are going to experience the more severe outcomes of COVID-19 is still a major need, to focus therapeutic efforts, reducing hospitalization and mitigating drug adverse effects. Microbial communities colonizing the respiratory tract exert significant effects on host immune responses, influencing the susceptibility to infectious agents. Through 16S rDNAseq we characterized the upper airways’ microbiota of 192 subjects with nasopharyngeal swab positive for SARS-CoV-2. Patients were divided into groups based on the presence of symptoms, pneumonia severity, and need for oxygen therapy or intubation. Indeed, unlike most of the literature, our study focuses on identifying microbial signatures predictive of disease progression rather than on the probability of infection itself, for which a consensus is lacking. Diversity, differential abundance, and network analysis at different taxonomic levels were synergistically adopted, in a robust bioinformatic pipeline, highlighting novel possible taxa correlated with patients’ disease progression to intubation.

Abstract Human cellular reprogramming to induced pluripotency is still an inefficient process, which has hindered studying the role of critical intermediate stages. Here we take advantage of high efficiency reprogramming in microfluidics and temporal multi-omics to identify and resolve distinct sub-populations and their interactions. We perform secretome analysis and single-cell transcriptomics to show functional extrinsic pathways of protein communication between reprogramming sub-populations and the re-shaping of a permissive extracellular environment. We pinpoint the HGF/MET/STAT3 axis as a potent enhancer of reprogramming, which acts via HGF accumulation within the confined system of microfluidics, and in conventional dishes needs to be supplied exogenously to enhance efficiency. Our data suggest that human cellular reprogramming is a transcription factor-driven process that it is deeply dependent on extracellular context and cell population determinants.

Abstract A growing body of evidence links gut microbiota changes with inflammatory bowel disease (IBD), raising the potential benefit of exploiting metagenomics data for non-invasive IBD diagnostics. The sbv IMPROVER metagenomics diagnosis for inflammatory bowel disease challenge investigated computational metagenomics methods for discriminating IBD and nonIBD subjects. Participants in this challenge were given independent training and test metagenomics data from IBD and nonIBD subjects, which could be wither either raw read data (sub-challenge 1, SC1) or processed Taxonomy- and Function-based profiles (sub-challenge 2, SC2). A total of 81 anonymized submissions were received between September 2019 and March 2020. Most participants’ predictions performed better than random predictions in classifying IBD versus nonIBD, Ulcerative Colitis (UC) versus nonIBD, and Crohn’s Disease (CD) versus nonIBD. However, discrimination between UC and CD remains challenging, with the classification quality similar to the set of random predictions. We analyzed the class prediction accuracy, the metagenomics features by the teams, and computational methods used. These results will be openly shared with the scientific community to help advance IBD research and illustrate the application of a range of computational methodologies for effective metagenomic classification.

Abstract Background Amyotrophic Lateral Sclerosis (ALS) is a neurodegenerative disease whose spreading and progression mechanisms are still unclear. The ability to predict ALS prognosis would improve the patients’ quality of life and support clinicians in planning treatments. In this paper, we investigate ALS evolution trajectories using Process Mining (PM) techniques enriched to both easily mine processes and automatically reveal how the pathways differentiate according to patients’ characteristics. Methods We consider data collected in two distinct data sources, namely the Pooled Resource Open-Access ALS Clinical Trials (PRO-ACT) dataset and a real-world clinical register (ALS–BS) including data of patients followed up in two tertiary clinical centers of Brescia (Italy). With a focus on the functional abilities progressively impaired as the disease progresses, we use two Process Discovery methods, namely the Directly-Follows Graph and the CareFlow Miner, to mine the population disease trajectories on the PRO-ACT dataset. We characterize the impairment trajectories in terms of patterns, timing, and probabilities, and investigate the effect of some patients’ characteristics at onset on the followed paths. Finally, we perform a comparative study of the impairment trajectories mined in PRO-ACT versus ALS–BS. Results We delineate the progression pathways on PRO-ACT, identifying the predominant disabilities at different stages of the disease: for instance, 85% of patients enter the trials without disabilities, and 48% of them experience the impairment of Walking/Self-care abilities first. We then test how a spinal onset increases the risk of experiencing the loss of Walking/Self-care ability as first impairment (52% vs. 27% of patients develop it as the first impairment in the spinal vs. the bulbar cohorts, respectively), as well as how an older age at onset corresponds to a more rapid progression to death. When compared, the PRO-ACT and the ALS–BS patient populations present some similarities in terms of natural progression of the disease, as well as some differences in terms of observed trajectories plausibly due to the trial scheduling and recruitment criteria. Conclusions We exploited PM to provide an overview of the evolution scenarios of an ALS trial population and to preliminary compare it to the progression observed in a clinical cohort. Future work will focus on further improving the understanding of the disease progression mechanisms, by including additional real-world subjects as well as by extending the set of events considered in the impairment trajectories.

Abstract Background Results of cardiovascular outcome trials enabled a shift from “treat-to-target” to “treat-to-benefit” paradigm in the management of type 2 diabetes (T2D). However, studies validating such approach are limited. Here, we examined whether treatment according to international recommendations for the pharmacological management of T2D had an impact on long-term outcomes. Methods This was an observational study conducted on outpatient data collected in 2008–2018 (i.e. prior to the “treat-to-benefit” shift). We defined 6 domains of treatment based on the ADA/EASD consensus covering all disease stages: first- and second-line treatment, intensification, use of insulin, cardioprotective, and weight-affecting drugs. At each visit, patients were included in Group 1 if at least one domain deviated from recommendation or in Group 2 if aligned with recommendations. We used Cox proportional hazard models with time-dependent co-variates or Cox marginal structural models (with inverse-probability of treatment weighing evaluated at each visit) to adjust for confounding factors and evaluate three outcomes: major adverse cardiovascular events (MACE), hospitalization for heart failure or cardiovascular mortality (HF-CVM), and all-cause mortality. Results We included 5419 patients, on average 66-year old, 41% women, with a baseline diabetes duration of 7.6 years. Only 11.7% had pre-existing cardiovascular disease. During a median follow-up of 7.3 years, patients were seen 12 times at the clinic, and we recorded 1325 MACE, 1593 HF-CVM, and 917 deaths. By the end of the study, each patient spent on average 63.6% of time in Group 1. In the fully adjusted model, being always in Group 2 was associated with a 45% lower risk of MACE (HR 0.55; 95% C.I. 0.46–0.66; p

Abstract Amyotrophic lateral sclerosis (ALS) is a highly complex and heterogeneous neurodegenerative disease that affects motor neurons. Since life expectancy is relatively low, it is essential to promptly understand the course of the disease to better target the patient’s treatment. Predictive models for disease progression are thus of great interest. One of the most extensive and well-studied open-access data resources for ALS is the Pooled Resource Open-Access ALS Clinical Trials (PRO-ACT) repository. In 2015, the DREAM-Phil Bowen ALS Prediction Prize4Life Challenge was held on PRO-ACT data, where competitors were asked to develop machine learning algorithms to predict disease progression measured through the slope of the ALSFRS score between 3 and 12 months. However, although it has already been successfully applied in several studies on ALS patients, to the best of our knowledge deep learning approaches still remain unexplored on the ALSFRS slope prediction in PRO-ACT cohort. Here, we investigate how deep learning models perform in predicting ALS progression using the PRO-ACT data. We developed three models based on different architectures that showed comparable or better performance with respect to the state-of-the-art models, thus representing a valid alternative to predict ALS disease progression.

Abstract Aim Treatment algorithms define lines of glucose lowering medications (GLM) for the management of type 2 diabetes (T2D), but whether therapeutic trajectories are associated with major adverse cardiovascular events (MACE) is unclear. We explored whether the temporal resolution of GLM usage discriminates patients who experienced a 4P-MACE (heart failure, myocardial infarction, stroke, death for all causes). Methods We used an administrative database (Veneto region, North-East Italy, 2011–2018) and implemented recurrent neural networks (RNN) with outcome-specific attention maps. The model input included age, sex, diabetes duration, and a matrix of GLM pattern before the 4P-MACE or censoring. Model output was discrimination, reported as area under receiver characteristic curve (AUROC). Attention maps were produced to show medications whose time-resolved trajectories were the most important for discrimination. Results The analysis was conducted on 147,135 patients for training and model selection and on 10,000 patients for validation. Collected data spanned a period of ~ 6 years. The RNN model efficiently discriminated temporal patterns of GLM ending in a 4P-MACE vs. those ending in an event-free censoring with an AUROC of 0.911 (95% C.I. 0.904–0.919). This excellent performance was significantly better than that of other models not incorporating time-resolved GLM trajectories: (i) a logistic regression on the bag-of-words encoding all GLM ever taken by the patient (AUROC 0.754; 95% C.I. 0.743–0.765); (ii) a model including the sequence of GLM without temporal relationships (AUROC 0.749; 95% C.I. 0.737–0.761); (iii) a RNN model with the same construction rules but including a time-inverted or randomised order of GLM. Attention maps identified the time-resolved pattern of most common first-line (metformin), second-line (sulphonylureas) GLM, and insulin (glargine) as those determining discrimination capacity. Conclusions The time-resolved pattern of GLM use identified patients with subsequent cardiovascular events better than the mere list or sequence of prescribed GLM. Thus, a patient’s therapeutic trajectory could determine disease outcomes.

Abstract Predicting the risk of cardiovascular complications, in particular heart failure hospitalisation (HHF), can improve the management of type 2 diabetes (T2D). Most predictive models proposed so far rely on clinical data not available at the higher Institutional level. Therefore, it is of interest to assess the risk of HHF in people with T2D using administrative claims data only, which are more easily obtainable and could allow public health systems to identify high-risk individuals. In this paper, the administrative claims of > 175,000 patients with T2D were used to develop a new risk score for HHF based on Cox regression. Internal validation on the administrative data cohort yielded satisfactory results in terms of discrimination (max AUROC = 0.792, C-index = 0.786) and calibration (Hosmer–Lemeshow test p value

Abstract Background 16S rRNA-gene sequencing is a valuable approach to characterize the taxonomic content of the whole bacterial population inhabiting a metabolic and spatial niche, providing an important opportunity to study bacteria and their role in many health and environmental mechanisms. The analysis of data produced by amplicon sequencing, however, brings very specific methodological issues that need to be properly addressed to obtain reliable biological conclusions. Among these, 16S count data tend to be very sparse, with many null values reflecting species that are present but got unobserved due to the multiplexing constraints. However, current data workflows do not consider a step in which the information about unobserved species is recovered. Results In this work, we evaluate for the first time the effects of introducing in the 16S data workflow a new preprocessing step, zero-imputation, to recover this lost information. Due to the lack of published zero-imputation methods specifically designed for 16S count data, we considered a set of zero-imputation strategies available for other frameworks, and benchmarked them using in silico 16S count data reflecting different experimental designs. Additionally, we assessed the effect of combining zero-imputation and normalization, i.e. the only preprocessing step in current 16S workflow. Overall, we benchmarked 35 16S preprocessing pipelines assessing their ability to handle data sparsity, identify species presence/absence, recovery sample proportional abundance distributions, and improve typical downstream analyses such as computation of alpha and beta diversity indices and differential abundance analysis. Conclusions The results clearly show that 16S data analysis greatly benefits from a properly-performed zero-imputation step, despite the choice of the right zero-imputation method having a pivotal role. In addition, we identify a set of best-performing pipelines that could be a valuable indication for data analysts.

Abstract Aim We aimed to compare cardiovascular outcomes of patients with type 2 diabetes (T2D) who initiated GLP-1 receptor agonists (GLP-1RA) or basal insulin (BI) under routine care. Methods We accessed the administrative claims database of the Veneto Region (Italy) to identify new users of GLP-1RA or BI in 2014–2018. Propensity score matching (PSM) was implemented to obtain two cohorts of patients with superimposable characteristics. The primary endpoint was the 3-point major adverse cardiovascular events (3P-MACE). Secondary endpoints included 3P-MACE components, hospitalization for heart failure, revascularizations, and adverse events. Results From a background population of 5,242,201 citizens, 330,193 were identified as having diabetes. PSM produced two very well matched cohorts of 4063 patients each, who initiated GLP-1RA or BI after an average of 2.5 other diabetes drug classes. Patients were 63-year-old and only 15% had a baseline history of cardiovascular disease. During a median follow-up of 24 months in the intention-to-treat analysis, 3P-MACE occurred less frequently in the GLP-1RA cohort (HR versus BI 0.59; 95% CI 0.50–0.71; p

Abstract Background The ability to rapidly adapt to adverse environmental conditions represents the key of success of many pathogens and, in particular, of Mycobacterium tuberculosis. Upon exposition to heat shock, antibiotics or other sources of stress, appropriate responses in terms of genes transcription and proteins activity are activated leading part of a genetically identical bacterial population to express a different phenotype, namely to develop persistence. When the stress response network is mathematically described by an ordinary differential equations model, development of persistence in the bacterial population is associated with bistability of the model, since different emerging phenotypes are represented by different stable steady states. Results In this work, we develop a mathematical model of SigE stress response network that incorporates interactions not considered in mathematical models currently available in the literature. We provide, through involved analytical computations, accurate approximations of the system’s nullclines, and exploit the obtained expressions to determine, in a reliable though computationally efficient way, the number of equilibrium points of the system. Conclusions Theoretical analysis and perturbation experiments point out the crucial role played by the degradation pathway involving RseA, the anti-sigma factor of SigE, for coexistence of two stable equilibria and the emergence of bistability. Our results also indicate that a fine control on RseA concentration is a necessary requirement in order for the system to exhibit bistability.

Background and Objectives: This study aimed to investigate the causes of continuous deep fluctuations in the absolute lymphocyte count (ALC) in an untreated patient with Chronic Lymphocytic Leukemia (CLL), who has had a favorable prognosis since the time of diagnosis. Up until now, the patient has voluntarily chosen to adopt a predominantly vegetarian and fruitarian diet, along with prolonged periods of total fasting (ranging from 4 to 39 days) each year. Materials and Methods: For this purpose, we decided to analyze the whole transcriptome profiling of peripheral blood (PB) CD19+ cells from the patient (#1) at different time-points vs. the same cells of five other untreated CLL patients who followed a varied diet. Consequently, the CLL patients were categorized as follows: the 1st group comprised patient #1 at 20 different time-points (16 time-points during nutrition and 4 time-points during fasting), whereas the 2nd group included only one time point for each of the patients (#2, #3, #4, #5, and #6) as they followed a varied diet. We performed microarray experiments using a powerful tool, the Affymetrix Human Clariom™ D Pico Assay, to generate high-fidelity biomarker signatures. Statistical analysis was employed to identify differentially expressed genes and to perform sample clustering. Results: The lymphocytosis trend in patient #1 showed recurring fluctuations since the time of diagnosis. Interestingly, we observed that approximately 4–6 weeks after the conclusion of fasting periods, the absolute lymphocyte count was reduced by about half. The gene expression profiling analysis revealed that nine genes were statistically differently expressed between the 1st group and the 2nd group. Specifically, IGLC3, RPS26, CHPT1, and PCDH9 were under expressed in the 1st group compared to the 2nd group of CLL patients. Conversely, IGHV3-43, IGKV3D-20, PLEKHA1, CYBB, and GABRB2 were over-expressed in the 1st group when compared to the 2nd group of CLL patients. Furthermore, clustering analysis validated that all the samples from patient #1 clustered together, showing clear separation from the samples of the other CLL patients. Conclusions: This study unveiled a small gene expression signature consisting of nine genes that distinguished an untreated CLL patient who followed prolonged periods of total fasting, maintaining a gradual growth trend of lymphocytosis, compared to five untreated CLL patients with a varied diet. Future investigations focusing on patient #1 could potentially shed light on the role of prolonged periodic fasting and the implication of this specific gene signature in sustaining the lymphocytosis trend and the favorable course of the disease.

The development of increasingly efficient and cost-effective high throughput DNA sequencing techniques has enhanced the possibility of studying complex microbial systems. Recently, researchers have shown great interest in studying the microorganisms that characterise different ecological niches. Differential abundance analysis aims to find the differences in the abundance of each taxa between two classes of subjects or samples, assigning a significance value to each comparison. Several bioinformatic methods have been specifically developed, taking into account the challenges of microbiome data, such as sparsity, the different sequencing depth constraint between samples and compositionality. Differential abundance analysis has led to important conclusions in different fields, from health to the environment. However, the lack of a known biological truth makes it difficult to validate the results obtained. In this work we exploit metaSPARSim, a microbial sequencing count data simulator, to simulate data with differential abundance features between experimental groups. We perform a complete comparison of recently developed and established methods on a common benchmark with great effort to the reliability of both the simulated scenarios and the evaluation metrics. The performance overview includes the investigation of numerous scenarios, studying the effect on methods' results on the main covariates such as sample size, percentage of differentially abundant features, sequencing depth, feature variability, normalisation approach and ecological niches. Mainly, we find that methods show a good control of the type I error and, generally, also of the false discovery rate at high sample size, while recall seem to depend on the dataset and sample size.

Abstract Background Clinical registers constitute an invaluable resource in the medical data-driven decision making context. Accurate machine learning and data mining approaches on these data can lead to faster diagnosis, definition of tailored interventions, and improved outcome prediction. A typical issue when implementing such approaches is the almost unavoidable presence of missing values in the collected data. In this work, we propose an imputation algorithm based on a mutual information-weighted k-nearest neighbours approach, able to handle the simultaneous presence of missing information in different types of variables. We developed and validated the method on a clinical register, constituted by the information collected over subsequent screening visits of a cohort of patients affected by amyotrophic lateral sclerosis. Methods For each subject with missing data to be imputed, we create a feature vector constituted by the information collected over his/her first three months of visits. This vector is used as sample in a k-nearest neighbours procedure, in order to select, among the other patients, the ones with the most similar temporal evolution of the disease over time. An ad hoc similarity metric was implemented for the sample comparison, capable of handling the different nature of the data, the presence of multiple missing values and include the cross-information among features captured by the mutual information statistic. Results We validated the proposed imputation method on an independent test set, comparing its performance with those of three state-of-the-art competitors, resulting in better performance. We further assessed the validity of our algorithm by comparing the performance of a survival classifier built on the data imputed with our method versus the one built on the data imputed with the best-performing competitor. Conclusions Imputation of missing data is a crucial –and often mandatory– step when working with real-world datasets. The algorithm proposed in this work could effectively impute an amyotrophic lateral sclerosis clinical dataset, by handling the temporal and the mixed-type nature of the data and by exploiting the cross-information among features. We also showed how the imputation quality can affect a machine learning task.

Abstract Background Cardiovascular outcome trials in high-risk patients showed that some GLP-1 receptor agonists (GLP-1RA), but not dipeptidyl-peptidase-4 inhibitors (DPP-4i), can prevent cardiovascular events in type 2 diabetes (T2D). Since no trial has directly compared these two classes of drugs, we performed a comparative outcome analysis using real-world data. Methods From a database of ~ 5 million people from North-East Italy, we retrospectively identified initiators of GLP-1RA or DPP-4i from 2011 to 2018. We obtained two balanced cohorts by 1:1 propensity score matching. The primary outcome was the 3-point major adverse cardiovascular events (3P-MACE; a composite of death, myocardial infarction, or stroke). 3P-MACE components and hospitalization for heart failure were secondary outcomes. Results From 330,193 individuals with T2D, we extracted two matched cohorts of 2807 GLP-1RA and 2807 DPP-4i initiators, followed for a median of 18 months. On average, patients were 63 years old, 60% male; 15% had pre-existing cardiovascular disease. The rate of 3P-MACE was lower in patients treated with GLP-1RA compared to DPP4i (23.5 vs. 34.9 events per 1000 person-years; HR: 0.67; 95% C.I. 0.53–0.86; p = 0.002). Rates of myocardial infarction (HR 0.67; 95% C.I. 0.50–0.91; p = 0.011) and all-cause death (HR 0.58; 95% C.I. 0.35–0.96; p = 0.034) were lower among GLP-1RA initiators. The as-treated and intention-to-treat approaches yielded similar results. Conclusions Patients initiating a GLP-1RA in clinical practice had better cardiovascular outcomes than similar patients who initiated a DPP-4i. These data strongly confirm findings from cardiovascular outcome trials in a lower risk population.

Interferon-α (IFN-α) comprises a family of 13 cytokines involved in the modulation of antiviral, immune, and anticancer responses by orchestrating a complex transcriptional network. The activation of IFN-α signaling pathway in endothelial cells results in decreased proliferation and migration, ultimately leading to suppression of angiogenesis. In this study, we knocked-down the expression of seven established or candidate modulators of IFN-α response in endothelial cells to reconstruct a gene regulatory network and to investigate the antiangiogenic activity of IFN-α. This genetic perturbation approach, along with the analysis of interferon-induced gene expression dynamics, highlighted a complex and highly interconnected network, in which the angiostatic chemokine C-X-C Motif Chemokine Ligand 10 (CXCL10) was a central node targeted by multiple modulators. IFN-α-induced secretion of CXCL10 protein by endothelial cells was blunted by the silencing of Signal Transducer and Activator of Transcription 1 (STAT1) and of Interferon Regulatory Factor 1 (IRF1) and it was exacerbated by the silencing of Ubiquitin Specific Peptidase 18 (USP18). In vitro sprouting assay, which mimics in vivo angiogenesis, confirmed STAT1 as a positive modulator and USP18 as a negative modulator of IFN-α-mediated sprouting suppression. Our data reveal an unprecedented physiological regulation of angiogenesis in endothelial cells through a tonic IFN-α signaling, whose enhancement could represent a viable strategy to suppress tumor neoangiogenesis.

Abstract Background In the last few years, 16S rRNA gene sequencing (16S rDNA-seq) has seen a surprisingly rapid increase in election rate as a methodology to perform microbial community studies. Despite the considerable popularity of this technique, an exiguous number of specific tools are currently available for proper 16S rDNA-seq count data preprocessing and simulation. Indeed, the great majority of tools have been developed adapting methodologies previously used for bulk RNA-seq data, with poor assessment of their applicability in the metagenomics field. For such tools and the few ones specifically developed for 16S rDNA-seq data, performance assessment is challenging, mainly due to the complex nature of the data and the lack of realistic simulation models. In fact, to the best of our knowledge, no software thought for data simulation are available to directly obtain synthetic 16S rDNA-seq count tables that properly model heavy sparsity and compositionality typical of these data. Results In this paper we present metaSPARSim, a sparse count matrix simulator intended for usage in development of 16S rDNA-seq metagenomic data processing pipelines. metaSPARSim implements a new generative process that models the sequencing process with a Multivariate Hypergeometric distribution in order to realistically simulate 16S rDNA-seq count table, resembling real experimental data compositionality and sparsity. It provides ready-to-use count matrices and comes with the possibility to reproduce different pre-coded scenarios and to estimate simulation parameters from real experimental data. The tool is made available at http://sysbiobig.dei.unipd.it/?q=Software#metaSPARSimand https://gitlab.com/sysbiobig/metasparsim. Conclusion metaSPARSim is able to generate count matrices resembling real 16S rDNA-seq data. The availability of count data simulators is extremely valuable both for methods developers, for which a ground truth for tools validation is needed, and for users who want to assess state of the art analysis tools for choosing the most accurate one. Thus, we believe that metaSPARSim is a valuable tool for researchers involved in developing, testing and using robust and reliable data analysis methods in the context of 16S rRNA gene sequencing.

Abstract Background Complication screening is recommended for patients with type 2 diabetes (T2D), but the optimal screening intensity and schedules are unknown. In this study, we evaluated whether intensive versus standard complication screening affects long-term cardiovascular outcomes. Methods In this observational study, we included 368 T2D patients referred for intensive screening provided as a 1-day session of clinical–instrumental evaluation of diabetic complications, followed by dedicated counseling. From a total of 4906 patients, we selected control T2D patients who underwent standard complication screening at different visits, by 2:1 propensity score matching. The primary endpoint was the 4p-MACE, defined as cardiovascular mortality, or non-fatal myocardial infarction, stroke, or heart failure. The Cox proportional regression analyses was used to compare outcome occurrence in the two groups, adjusted for residual confounders. Results 357 patients from the intensive screening group (out of 368) were matched with 683 patients in the standard screening group. Clinical characteristics were well balanced between the two groups, except for a slightly higher prevalence of microangiopathy in the intensive group (56% vs 50%; standardized mean difference 0.11, p = 0.1). Median follow-up was 5.6 years. The adjusted incidence of 4p-MACE was significantly lower in the intensive versus standard screening group (HR 0.70; 95% CI 0.52–0.95; p = 0.02). All components of the primary endpoint had nominally lower rates in the intensive versus standard screening group, which was particularly significant for heart failure (HR 0.43; 95% CI 0.22–0.83; p = 0.01). Conclusion Among T2D patients attending a specialist outpatient clinic, intensive complication screening is followed by better long-term cardiovascular outcomes. No significant effect was noted for cardiovascular and all-cause mortality and the benefit was mainly driven by a reduced rate of hospitalization for heart failure.

Abstract Background Amyotrophic lateral sclerosis (ALS) is an adult-onset neurodegenerative disease progressively affecting upper and lower motor neurons in the brain and spinal cord. Mean life expectancy is three to five years, with paralysis of muscles, respiratory failure and loss of vital functions being the common causes of death. Clinical manifestations of ALS are heterogeneous due to the mix of anatomic regions involvement and the variability in disease course; consequently, diagnosis and prognosis at the level of individual patient is really challenging. Prediction of ALS progression and stratification of patients into meaningful subgroups have been long-standing interests to clinical practice, research and drug development. Methods We developed a Dynamic Bayesian Network (DBN) model on more than 4500 ALS patients included in the Pooled Resource Open-Access ALS Clinical Trials Database (PRO-ACT), in order to detect probabilistic relationships among clinical variables and identify risk factors related to survival and loss of vital functions. Furthermore, the DBN was used to simulate the temporal evolution of an ALS cohort predicting survival and the time to impairment of vital functions (communication, swallowing, gait and respiration). A first attempt to stratify patients by risk factors and simulate the progression of ALS subgroups was also implemented. Results The DBN model provided the prediction of ALS most probable trajectories over time in terms of important clinical outcomes, including survival and loss of autonomy in functional domains. Furthermore, it allowed the identification of biomarkers related to patients’ clinical status as well as vital functions, and unrevealed their probabilistic relationships. For instance, DBN found that bicarbonate and calcium levels influence survival time; moreover, the model evidenced dependencies over time among phosphorus level, movement impairment and creatinine. Finally, our model provided a tool to stratify patients into subgroups of different prognosis studying the effect of specific variables, or combinations of them, on either survival time or time to loss of autonomy in specific functional domains. Conclusions The analysis of the risk factors and the simulation allowed by our DBN model might enable better support for ALS prognosis as well as a deeper insight into disease manifestations, in a context of a personalized medicine approach.

Fecal microbiota transplantation (FMT) represents a very promising approach to decreasing disease activity in canine chronic enteropathies (CE). However, the relationship between remission mechanisms and microbiome changes has not been elucidated yet. The main objective of this study was to report the clinical effects of oral freeze-dried FMT in CE dogs, comparing the fecal microbiomes of three groups: pre-FMT CE-affected dogs, post-FMT dogs, and healthy dogs. Diversity analysis, differential abundance analysis, and machine learning algorithms were applied to investigate the differences in microbiome composition between healthy and pre-FMT samples, while Canine Chronic Enteropathy Clinical Activity Index (CCECAI) changes and microbial diversity metrics were used to evaluate FMT effects. In the healthy/pre-FMT comparison, significant differences were noted in alpha and beta diversity and a list of differentially abundant taxa was identified, while machine learning algorithms predicted sample categories with 0.97 (random forest) and 0.87 (sPLS-DA) accuracy. Clinical signs of improvement were observed in 74% (20/27) of CE-affected dogs, together with a statistically significant decrease in CCECAI (median value from 5 to 2 median). Alpha and beta diversity variations between pre- and post-FMT were observed for each receiver, with a high heterogeneity in the response. This highlighted the necessity for further research on a larger dataset that could identify different healing patterns of microbiome changes.

Newly differentiated pancreatic β cells lack proper insulin secretion profiles of mature functional β cells. The global gene expression differences between paired immature and mature β cells have been studied, but the dynamics of transcriptional events, correlating with temporal development of glucose-stimulated insulin secretion (GSIS), remain to be fully defined. This aspect is important to identify which genes and pathways are necessary for β-cell development or for maturation, as defective insulin secretion is linked with diseases such as diabetes. In this study, we assayed through RNA sequencing the global gene expression across six β-cell developmental stages in mice, spanning from β-cell progenitor to mature β cells. A computational pipeline then selected genes differentially expressed with respect to progenitors and clustered them into groups with distinct temporal patterns associated with biological functions and pathways. These patterns were finally correlated with experimental GSIS, calcium influx, and insulin granule formation data. Gene expression temporal profiling revealed the timing of important biological processes across β-cell maturation, such as the deregulation of β-cell developmental pathways and the activation of molecular machineries for vesicle biosynthesis and transport, signal transduction of transmembrane receptors, and glucose-induced Ca2+ influx, which were established over a week before β-cell maturation completes. In particular, β cells developed robust insulin secretion at high glucose several days after birth, coincident with the establishment of glucose-induced calcium influx. Yet the neonatal β cells displayed high basal insulin secretion, which decreased to the low levels found in mature β cells only a week later. Different genes associated with calcium-mediated processes, whose alterations are linked with insulin resistance and deregulation of glucose homeostasis, showed increased expression across β-cell stages, in accordance with the temporal acquisition of proper GSIS. Our temporal gene expression pattern analysis provided a comprehensive database of the underlying molecular components and biological mechanisms driving β-cell maturation at different temporal stages, which are fundamental for better control of the in vitro production of functional β cells from human embryonic stem/induced pluripotent cell for transplantation-based type 1 diabetes therapy.

Introduction Many predictive models for incident type 2 diabetes (T2D) exist, but these models are not used frequently for public health management. Barriers to their application include (1) the problem of model choice (some models are applicable only to certain ethnic groups), (2) missing input variables, and (3) the lack of calibration. While (1) and (2) drives to missing predictions, (3) causes inaccurate incidence predictions. In this paper, a combined T2D risk model for public health management that addresses these three issues is developed.Research design and methods The combined T2D risk model combines eight existing predictive models by weighted average to overcome the problem of missing incidence predictions. Moreover, the combined model implements a simple recalibration strategy in which the risk scores are rescaled based on the T2D incidence in the target population. The performance of the combined model was compared with that of the eight existing models using data from two test datasets extracted from the Multi-Ethnic Study of Atherosclerosis (MESA; n=1031) and the English Longitudinal Study of Ageing (ELSA; n=4820). Metrics of discrimination, calibration, and missing incidence predictions were used for the assessment.Results The combined T2D model performed well in terms of both discrimination (concordance index: 0.83 on MESA; 0.77 on ELSA) and calibration (expected to observed event ratio: 1.00 on MESA; 1.17 on ELSA), similarly to the best-performing existing models. However, while the existing models yielded a large percentage of missing predictions (17%–45% on MESA; 63%–64% on ELSA), this was negligible with the combined model (0% on MESA, 4% on ELSA).Conclusions Leveraging on existing literature T2D predictive models, a simple approach based on risk score rescaling and averaging was shown to provide accurate and robust incidence predictions, overcoming the problem of recalibration and missing predictions in practical application of predictive models.

Introduction Sodium glucose cotransporter-2 inhibitors (SGLT2i) and glucagon-like peptide-1 receptor agonists (GLP-1RA) protect type 2 diabetic (T2D) patients from cardiovascular events, but no trial has directly compared their cardiovascular effects. We aimed to address this gap using real-world data.Research design and methods We performed a retrospective real-world study on a population of ~5 million inhabitants from North-East Italy. We identified T2D patients who received new prescription of SGLT2i or GLP-1RA from 2014 to 2018. SGLT2i and GLP-1RA initiators were matched 1:1 by propensity scores. The primary outcome was a composite of all-cause death, myocardial infarction, and stroke (three-point major adverse cardiovascular events (3P-MACE)). Secondary endpoints were each component of the primary endpoint, hospitalization for heart failure (HF), revascularization, hospitalization for cardiovascular causes, and adverse events.Results From a population of 330 193 diabetic patients, we followed 8596 SGLT2i and GLP-1RA matched initiators for a median of 13 months. Patients in both groups were on average 63 years old, 63% men, and 18% had pre-existing cardiovascular disease. T2D patients treated with SGLT2i versus GLP-1RA, experienced a lower rate of 3P-MACE (HR 0.68; 95% CI 0.61 to 0.99; p=0.043), myocardial infarction (HR 0.72; 95% CI 0.53 to 0.98; p=0.035), hospitalization for HF (HR 0.59; 95% CI 0.35 to 0.99; p=0.048), and hospitalization for cardiovascular causes (HR 0.82; 95% CI 0.69 to 0.99; p=0.037). Adverse events were not significantly different between the two groups.Conclusions In the absence of dedicated trials, this observational study suggests that SGLT2i may be more effective than GLP-1RA in improving cardiovascular outcomes of T2D.Trial registration number NCT04184947.

Abstract Background Targeted amplicon sequencing of the 16S ribosomal RNA gene is one of the key tools for studying microbial diversity. The accuracy of this approach strongly depends on the choice of primer pairs and, in particular, on the balance between efficiency, specificity and sensitivity in the amplification of the different bacterial 16S sequences contained in a sample. There is thus the need for computational methods to design optimal bacterial 16S primers able to take into account the knowledge provided by the new sequencing technologies. Results We propose here a computational method for optimizing the choice of primer sets, based on multi-objective optimization, which simultaneously: 1) maximizes efficiency and specificity of target amplification; 2) maximizes the number of different bacterial 16S sequences matched by at least one primer; 3) minimizes the differences in the number of primers matching each bacterial 16S sequence. Our algorithm can be applied to any desired amplicon length without affecting computational performance. The source code of the developed algorithm is released as the mopo16S software tool (Multi-Objective Primer Optimization for 16S experiments) under the GNU General Public License and is available at http://sysbiobig.dei.unipd.it/?q=Software#mopo16S. Conclusions Results show that our strategy is able to find better primer pairs than the ones available in the literature according to all three optimization criteria. We also experimentally validated three of the primer pairs identified by our method on multiple bacterial species, belonging to different genera and phyla. Results confirm the predicted efficiency and the ability to maximize the number of different bacterial 16S sequences matched by primers.

When building a predictive model for predicting a clinical outcome using machine learning techniques, the model developers are often interested in ranking the features according to their predictive ability. A commonly used approach to obtain a robust variable ranking is to apply recursive feature elimination (RFE) on multiple resamplings of the training set and then to aggregate the ranking results using the Borda count method. However, the presence of highly correlated features in the training set can deteriorate the ranking performance. In this work, we propose a variant of the method based on RFE and Borda count that takes into account the correlation between variables during the ranking procedure in order to improve the ranking performance in the presence of highly correlated features. The proposed algorithm is tested on simulated datasets in which the true variable importance is known and compared to the standard RFE-Borda count method. According to the root mean square error between the estimated rank and the true (i.e., simulated) feature importance, the proposed algorithm overcomes the standard RFE-Borda count method. Finally, the proposed algorithm is applied to a case study related to the development of a predictive model of type 2 diabetes onset.

Chronic myeloid leukemia (CML) is characterized by the constitutive tyrosine kinase activity of the oncoprotein BCR-ABL1 in myeloid progenitor cells that activates multiple signal transduction pathways leading to the leukemic phenotype. The tyrosine-kinase inhibitor (TKI) nilotinib inhibits the tyrosine kinase activity of BCR-ABL1 in CML patients. Despite the success of nilotinib treatment in patients with chronic-phase (CP) CML, a population of Philadelphia-positive (Ph+) quiescent stem cells escapes the drug activity and can lead to drug resistance. The molecular mechanism by which these quiescent cells remain insensitive is poorly understood. The aim of this study was to compare the gene expression profiling (GEP) of bone marrow (BM) CD34+/lin- cells from CP-CML patients at diagnosis and after 12 months of nilotinib treatment by microarray, in order to identify gene expression changes and the dysregulation of pathways due to nilotinib action. We selected BM CD34+/lin- cells from 78 CP-CML patients at diagnosis and after 12 months of first-line nilotinib therapy and microarray analysis was performed. GEP bioinformatic analyses identified 2,959 differently expressed probes and functional clustering determined some significantly enriched pathways between diagnosis and 12 months of nilotinib treatment. Among these pathways, we observed the under expression of 26 genes encoding proteins belonging to the cell cycle after 12 months of nilotinib treatment which led to the up-regulation of chromosome replication, cell proliferation, DNA replication, and DNA damage checkpoint at diagnosis. We demonstrated the under expression of the ATP-binding cassette (ABC) transporters ABCC4, ABCC5, and ABCD3 encoding proteins which pumped drugs out of the cells after 12 months of nilotinib. Moreover, GEP data demonstrated the deregulation of genes involved in the JAK-STAT signaling pathway. The down-regulation of JAK2, IL7, STAM, PIK3CA, PTPN11, RAF1, and SOS1 key genes after 12 months of nilotinib could demonstrate the up-regulation of cell cycle, proliferation and differentiation via MAPK and PI3K-AKT signaling pathways at diagnosis.

Chronic enteropathies (CE) are gastrointestinal diseases that afflict about one in five dogs in Europe. Conventional therapeutic approaches include dietary intervention, pharmacological treatment and probiotic supplements. The patient response can be highly variable and the interventions are often not resolutive. Moreover, the therapeutic strategy is usually planned (and gradually corrected) based on the patient’s response to empirical treatment, with few indirect gut health indicators useful to drive clinicians’ decisions. The ever-diminishing cost of high-throughput sequencing (HTS) allows clinicians to directly follow and characterise the evolution of the whole gut microbial community in order to highlight possible weaknesses. In this framework, faecal microbiome transplantation (FMT) is emerging as a feasible solution to CE, based on the implant of a balanced, eubiotic microbial community from a healthy donor to a dysbiotic patient. In this study, we report the promising results of FMT carried out in a 9-year-old dog suffering from CE for the last 3 years. The patient underwent a two-cycle oral treatment of FMT and the microbiota evolution was monitored by 16S rRNA gene sequencing both prior to FMT and after the two administrations. We evaluated the variation of microbial composition by calculating three different alpha diversity indices and compared the patient and donor data to a healthy control population of 94 dogs. After FMT, the patient’s microbiome and clinical parameters gradually shifted to values similar to those observed in healthy dogs. Symptoms disappeared during a follow-up period of six months after the second FMT. We believe that this study opens the door for potential applications of FMT in clinical veterinary practice and highlights the need to improve our knowledge on this relevant topic.

Autophagy is a primordial eukaryotic pathway, which provides the immune system with multiple mechanisms for the elimination of invading pathogens including Mycobacterium tuberculosis (Mtb). As a consequence, Mtb has evolved different strategies to hijack the autophagy process. Given the crucial role of human primary dendritic cells (DC) in host immunity control, we characterized Mtb-DC interplay by studying the contribution of cellular microRNAs (miRNAs) in the post-transcriptional regulation of autophagy related genes. From the expression profile of de-regulated miRNAs obtained in Mtb-infected human DC, we identified 7 miRNAs whose expression was previously found to be altered in specimens of TB patients. Among them, gene ontology analysis showed that miR-155, miR-155* and miR-146a target mRNAs with a significant enrichment in biological processes linked to autophagy. Interestingly, miR-155 was significantly stimulated by live and virulent Mtb and enriched in polysome-associated RNA fraction, where actively translated mRNAs reside. The putative pair interaction among the E2 conjugating enzyme involved in LC3-lipidation and autophagosome formation-ATG3-and miR-155 arose by target prediction analysis, was confirmed by both luciferase reporter assay and Atg3 immunoblotting analysis of miR-155-transfected DC, which showed also a consistent Atg3 protein and LC3 lipidated form reduction. Late in infection, when miR-155 expression peaked, both the level of Atg3 and the number of LC3 puncta per cell (autophagosomes) decreased dramatically. In accordance, miR-155 silencing rescued autophagosome number in Mtb infected DC and enhanced autolysosome fusion, thereby supporting a previously unidentified role of the miR-155 as inhibitor of ATG3 expression. Taken together, our findings suggest how Mtb can manipulate cellular miRNA expression to regulate Atg3 for its own survival, and highlight the importance to develop novel therapeutic strategies against tuberculosis that would boost autophagy.

The sequencing of the transcriptomes of single-cells, or single-cell RNA-sequencing, has now become the dominant technology for the identification of novel cell types and for the study of stochastic gene expression. In recent years, various tools for analyzing single-cell RNA-sequencing data have been proposed, many of them with the purpose of performing differentially expression analysis. In this work, we compare four different tools for single-cell RNA-sequencing differential expression, together with two popular methods originally developed for the analysis of bulk RNA-sequencing data, but largely applied to single-cell data. We discuss results obtained on two real and one synthetic dataset, along with considerations about the perspectives of single-cell differential expression analysis. In particular, we explore the methods performance in four different scenarios, mimicking different unimodal or bimodal distributions of the data, as characteristic of single-cell transcriptomics. We observed marked differences between the selected methods in terms of precision and recall, the number of detected differentially expressed genes and the overall performance. Globally, the results obtained in our study suggest that is difficult to identify a best performing tool and that efforts are needed to improve the methodologies for single-cell RNA-sequencing data analysis and gain better accuracy of results.

The task of the DREAM4 (Dialogue for Reverse Engineering Assessments and Methods) "Predictive signaling network modeling" challenge was to develop a method that, from single-stimulus/inhibitor data, reconstructs a cause-effect network to be used to predict the protein activity level in multi-stimulus/inhibitor experimental conditions. The method presented in this paper, one of the best performing in this challenge, consists of 3 steps: 1. Boolean tables are inferred from single-stimulus/inhibitor data to classify whether a particular combination of stimulus and inhibitor is affecting the protein. 2. A cause-effect network is reconstructed starting from these tables. 3. Training data are linearly combined according to rules inferred from the reconstructed network. This method, although simple, permits one to achieve a good performance providing reasonable predictions based on a reconstructed network compatible with knowledge from the literature. It can be potentially used to predict how signaling pathways are affected by different ligands and how this response is altered by diseases.

MotivationThe identification of robust lists of molecular biomarkers related to a disease is a fundamental step for early diagnosis and treatment. However, methodologies for the discovery of biomarkers using microarray data often provide results with limited overlap. These differences are imputable to 1) dataset size (few subjects with respect to the number of features); 2) heterogeneity of the disease; 3) heterogeneity of experimental protocols and computational pipelines employed in the analysis. In this paper, we focus on the first two issues and assess, both on simulated (through an in silico regulation network model) and real clinical datasets, the consistency of candidate biomarkers provided by a number of different methods.MethodsWe extensively simulated the effect of heterogeneity characteristic of complex diseases on different sets of microarray data. Heterogeneity was reproduced by simulating both intrinsic variability of the population and the alteration of regulatory mechanisms. Population variability was simulated by modeling evolution of a pool of subjects; then, a subset of them underwent alterations in regulatory mechanisms so as to mimic the disease state.ResultsThe simulated data allowed us to outline advantages and drawbacks of different methods across multiple studies and varying number of samples and to evaluate precision of feature selection on a benchmark with known biomarkers. Although comparable classification accuracy was reached by different methods, the use of external cross-validation loops is helpful in finding features with a higher degree of precision and stability. Application to real data confirmed these results.

Insulin action on protein synthesis (translation of transcripts) and post-translational modifications, especially of those involving the reversible modifications such as phosphorylation of various signaling proteins, are extensively studied but insulin effect on transcription of genes, especially of transcriptional temporal patterns remains to be fully defined.To identify significant transcriptional temporal patterns we utilized primary differentiated rat skeletal muscle myotubes which were treated with insulin and samples were collected every 20 min for 8 hours. Pooled samples at every hour were analyzed by gene array approach to measure transcript levels. The patterns of transcript levels were analyzed based on a novel method that integrates selection, clustering, and functional annotation to find the main temporal patterns associated to functional groups of differentially expressed genes. 326 genes were found to be differentially expressed in response to in vitro insulin administration in skeletal muscle myotubes. Approximately 20% of the genes that were differentially expressed were identified as belonging to the insulin signaling pathway. Characteristic transcriptional temporal patterns include: (a) a slow and gradual decrease in gene expression, (b) a gradual increase in gene expression reaching a peak at about 5 hours and then reaching a plateau or an initial decrease and other different variable pattern of increase in gene expression over time.The new method allows identifying characteristic dynamic responses to insulin stimulus, common to a number of genes and associated to the same functional group. The results demonstrate that insulin treatment elicited different clusters of gene transcript profile supporting a temporal regulation of gene expression by insulin in skeletal muscle cells.

BACKGROUND: In diabetes chronic hyperinsulinemia contributes to the instability of the atherosclerotic plaque and stimulates cellular proliferation through the activation of the MAP kinases, which in turn regulate cellular proliferation. However, it is not known whether insulin itself could increase the transcription of specific genes for cellular proliferation in the endothelium. Hence, the characterization of transcriptional modifications in endothelium is an important step for a better understanding of the mechanism of insulin action and the relationship between endothelial cell dysfunction and insulin resistance. METHODOLOGY AND PRINCIPAL FINDINGS: The transcriptional response of endothelial cells in the 440 minutes following insulin stimulation was monitored using microarrays and compared to a control condition. About 1700 genes were selected as differentially expressed based on their treated minus control profile, thus allowing the detection of even small but systematic changes in gene expression. Genes were clustered in 7 groups according to their time expression profile and classified into 15 functional categories that can support the biological effects of insulin, based on Gene Ontology enrichment analysis. In terms of endothelial function, the most prominent processes affected were NADH dehydrogenase activity, N-terminal myristoylation domain binding, nitric-oxide synthase regulator activity and growth factor binding. Pathway-based enrichment analysis revealed "Electron Transport Chain" significantly enriched. Results were validated on genes belonging to "Electron Transport Chain" pathway, using quantitative RT-PCR. CONCLUSIONS: As far as we know, this is the first systematic study in the literature monitoring transcriptional response to insulin in endothelial cells, in a time series microarray experiment. Since chronic hyperinsulinemia contributes to the instability of the atherosclerotic plaque and stimulates cellular proliferation, some of the genes identified in the present work are potential novel candidates in diabetes complications related to endothelial dysfunction.