Integrative network analysis reveals time-dependent molecular events underlying left ventricular remodeling in post-myocardial infarction patients.

To elucidate the time-resolved molecular events underlying the LV remodeling (LVR) process, we developed a large-scale network model that integrates the 24 molecular variables (plasma proteins and non-coding RNAs) collected in the REVE-2 study at four time points (baseline, 1month, 3months and 1year) after MI. The REVE-2 network model was built by extending the set of REVE-2 variables with their mechanistic context based on known molecular interactions (1310 nodes and 8639 edges). Changes in the molecular variables between the group of patients with high LVR (>20%) and low LVR (<20%) were used to identify active network modules within the clusters associated with progression of LVR, enabling assessment of time-resolved molecular changes. Although the majority of molecular changes occur at the baseline, two network modules specifically show an increasing number of active molecules throughout the post-MI follow up: one involved in muscle filament sliding, containing the major troponin forms and tropomyosin proteins, and the other associated with extracellular matrix disassembly, including matrix metalloproteinases, tissue inhibitors of metalloproteinases and laminin proteins. For the first time, integrative network analysis of molecular variables collected in REVE-2 patients with known molecular interactions allows insight into time-dependent mechanisms associated with LVR following MI, linking specific processes with LV structure alteration. In addition, the REVE-2 network model provides a shortlist of prioritized putative novel biomarker candidates for detection of LVR after MI event associated with a high risk of heart failure and is a valuable resource for further hypothesis generation.

The Network Library: a framework to rapidly integrate network biology resources.

MOTIVATION: Much of the biological knowledge accumulated over the last decades is stored in different databases governed by various organizations and institutes. Integrating and connecting these vast knowledge repositories is an extremely useful method to support life sciences research and help formulate novel hypotheses. RESULTS: We developed the Network Library (NL), a framework and toolset to rapidly integrate different knowledge sources to build a network biology resource that matches a specific research question. As a use-case we explore the interactions of genes related to heart failure with miRNAs and diseases through the integration of 6 databases. AVAILABILITY AND IMPLEMENTATION: The NL is open-source, developed in Java and available on Github (https://github.com/gsummer). CONTACT: georg.summer@gmail.com.

cyNeo4j: connecting Neo4j and Cytoscape.

UNLABELLED: We developed cyNeo4j, a Cytoscape App to link Cytoscape and Neo4j databases to utilize the performance and storage capacities Neo4j offers. We implemented a Neo4j NetworkAnalyzer, ForceAtlas2 layout and Cypher component to demonstrate the possibilities a distributed setup of Cytoscape and Neo4j have. AVAILABILITY AND IMPLEMENTATION: The app is available from the Cytoscape App Store at http://apps.cytoscape.org/apps/cyneo4j, the Neo4j plugins at www.github.com/gsummer/cyneo4j-parent and the community and commercial editions of Neo4j can be found at http://www.neo4j.com. CONTACT: georg.summer@gmail.com.

Drugs that reverse disease transcriptomic signatures are more effective in a mouse model of dyslipidemia.

High-throughput omics have proven invaluable in studying human disease, and yet day-to-day clinical practice still relies on physiological, non-omic markers. The metabolic syndrome, for example, is diagnosed and monitored by blood and urine indices such as blood cholesterol levels. Nevertheless, the association between the molecular and the physiological manifestations of the disease, especially in response to treatment, has not been investigated in a systematic manner. To this end, we studied a mouse model of diet-induced dyslipidemia and atherosclerosis that was subject to various drug treatments relevant to the disease in question. Both physiological data and gene expression data (from the liver and white adipose) were analyzed and compared. We find that treatments that restore gene expression patterns to their norm are associated with the successful restoration of physiological markers to their baselines. This holds in a tissue-specific manner­treatments that reverse the transcriptomic signatures of the disease in a particular tissue are associated with positive physiological effects in that tissue. Further, treatments that introduce large non-restorative gene expression alterations are associated with unfavorable physiological outcomes. These results provide a sound basis to in silico methods that rely on omic metrics for drug repurposing and drug discovery by searching for compounds that reverse a disease's omic signatures. Moreover, they highlight the need to develop drugs that restore the global cellular state to its healthy norm rather than rectify particular disease phenotypes.

A user-friendly workflow for analysis of Illumina gene expression bead array data available at the arrayanalysis.org portal.

BACKGROUND: Illumina whole-genome expression bead arrays are a widely used platform for transcriptomics. Most of the tools available for the analysis of the resulting data are not easily applicable by less experienced users. ArrayAnalysis.org provides researchers with an easy-to-use and comprehensive interface to the functionality of R and Bioconductor packages for microarray data analysis. As a modular open source project, it allows developers to contribute modules that provide support for additional types of data or extend workflows. RESULTS: To enable data analysis of Illumina bead arrays for a broad user community, we have developed a module for ArrayAnalysis.org that provides a free and user-friendly web interface for quality control and pre-processing for these arrays. This module can be used together with existing modules for statistical and pathway analysis to provide a full workflow for Illumina gene expression data analysis. The module accepts data exported from Illumina's GenomeStudio, and provides the user with quality control plots and normalized data. The outputs are directly linked to the existing statistics module of ArrayAnalysis.org, but can also be downloaded for further downstream analysis in third-party tools. CONCLUSIONS: The Illumina bead arrays analysis module is available at http://www.arrayanalysis.org . A user guide, a tutorial demonstrating the analysis of an example dataset, and R scripts are available. The module can be used as a starting point for statistical evaluation and pathway analysis provided on the website or to generate processed input data for a broad range of applications in life sciences research.

Identification of prognostic and diagnostic biomarkers of glucose intolerance in ApoE3Leiden mice.

The prevalence of diabetes mellitus Type 2 could be significantly reduced by early identification of subjects at risk, allowing for better prevention and earlier treatment. Glucose intolerance (GI) is a hallmark of the prediabetic stage. This study aims at identifying 1) prognostic biomarkers predicting the risk of developing GI later in life and 2) diagnostic biomarkers reflecting the degree of already manifest GI. To this end, disease development was followed over time in mice, and biomarkers were identified using lipidomics and transcriptomics. Young adult ApoE3Leiden mice were treated a high-fat diet for 12 wk to induce GI. Blood was collected before and during disease development. The individual extent of GI was determined with a glucose tolerance test and the area under the curve (AUC) was calculated for each animal. Subject-specific AUC values were correlated to the plasma lipidome (t = 0) and the white blood cell (WBC) transcriptome (t = 0, 6, and 12 wk) to identify prognostic and diagnostic biomarkers, respectively. The plasma ratio of specific free fatty acids prior to high-fat feeding (C16:1/C16:0, C18:1/C18:0 and C18:2/C22:6) was significantly correlated with the AUC and predictive for future GI. Subsequently, the expression level of specific WBC genes (Acss2, Arfgap1, Tfrc, Cox6b2, Barhl2, Abcb4, Cyp4b1, Sars2, Fgf16, and Tceal8) reflected the individual degree of GI during disease progression. Specific plasma free fatty acids as well as their ratio can be used to predict future GI. The expression levels of specific WBC genes can serve as easy accessible markers to diagnose and monitor already existing GI.

Systems biology analysis unravels the complementary action of combined rosuvastatin and ezetimibe therapy.

AIMS: Combination-drug therapy takes advantage of the complementary action of their individual components, thereby potentiating its therapeutic effect. Potential disadvantages include side effects that are not foreseen on basis of the data available from drug monotherapy. Here, we used a systems biology approach to understand both the efficacy and the side effects of a cholesterol-lowering drug-combination therapy on the basis of the biological pathways and molecular processes affected by each drug alone or in combination. METHODS AND RESULTS: ApoE*3Leiden transgenic mice, a mouse model with human-like cholesterol-lowering drug responses, were treated with rosuvastatin and ezetimibe, alone and in combination. Analyses included functional responses, viz. effects on cardiovascular risk factors, inflammation, and atherosclerosis, and measurement of global gene expression, and identification of enriched biological pathways and molecular processes. Combination therapy reduced plasma cholesterol, plasma inflammation markers, and atherosclerosis stronger than the single drugs did. Systems biology analysis at the level of biological processes shows that the therapeutic benefit of combined therapy is largely the result of additivity of the complementary mechanisms of action of the two single drugs. Importantly, combination therapy also exerted a significant effect on 16 additional and mostly NF-κB-linked signaling processes, 11 of which tended to be regulated in a similar direction with monotherapy. CONCLUSION: This study shows that gene expression analysis together with bioinformatics pathway analysis has the potential to help predict and identify drug combination-specific complementary and side effects.

Improved genome-wide localization by ChIP-chip using double-round T7 RNA polymerase-based amplification.

Chromatin immunoprecipitation combined with DNA microarrays (ChIP-chip) is a powerful technique to detect in vivo protein-DNA interactions. Due to low yields, ChIP assays of transcription factors generally require amplification of immunoprecipitated genomic DNA. Here, we present an adapted linear amplification method that involves two rounds of T7 RNA polymerase amplification (double-T7). Using this we could successfully amplify as little as 0.4 ng of ChIP DNA to sufficient amounts for microarray analysis. In addition, we compared the double-T7 method to the ligation-mediated polymerase chain reaction (LM-PCR) method in a ChIP-chip of the yeast transcription factor Gsm1p. The double-T7 protocol showed lower noise levels and stronger binding signals compared to LM-PCR. Both LM-PCR and double-T7 identified strongly bound genomic regions, but the double-T7 method increased sensitivity and specificity to allow detection of weaker binding sites.

Filling gaps in PPAR-alpha signaling through comparative nutrigenomics analysis.

BACKGROUND: The application of high-throughput genomic tools in nutrition research is a widespread practice. However, it is becoming increasingly clear that the outcome of individual expression studies is insufficient for the comprehensive understanding of such a complex field. Currently, the availability of the large amounts of expression data in public repositories has opened up new challenges on microarray data analyses. We have focused on PPARalpha, a ligand-activated transcription factor functioning as fatty acid sensor controlling the gene expression regulation of a large set of genes in various metabolic organs such as liver, small intestine or heart. The function of PPARalpha is strictly connected to the function of its target genes and, although many of these have already been identified, major elements of its physiological function remain to be uncovered. To further investigate the function of PPARalpha, we have applied a cross-species meta-analysis approach to integrate sixteen microarray datasets studying high fat diet and PPARalpha signal perturbations in different organisms. RESULTS: We identified 164 genes (MDEGs) that were differentially expressed in a constant way in response to a high fat diet or to perturbations in PPARs signalling. In particular, we found five genes in yeast which were highly conserved and homologous of PPARalpha targets in mammals, potential candidates to be used as models for the equivalent mammalian genes. Moreover, a screening of the MDEGs for all known transcription factor binding sites and the comparison with a human genome-wide screening of Peroxisome Proliferating Response Elements (PPRE), enabled us to identify, 20 new potential candidate genes that show, both binding site, both change in expression in the condition studied. Lastly, we found a non random localization of the differentially expressed genes in the genome. CONCLUSION: The results presented are potentially of great interest to resume the currently available expression data, exploiting the power of in silico analysis filtered by evolutionary conservation. The analysis enabled us to indicate potential gene candidates that could fill in the gaps with regards to the signalling of PPARalpha and, moreover, the non-random localization of the differentially expressed genes in the genome, suggest that epigenetic mechanisms are of importance in the regulation of the transcription operated by PPARalpha.

Molecular characterization of trimellitic anhydride-induced respiratory allergy in Brown Norway rats.

To contribute to the hazard identification of low molecular weight (LMW) respiratory allergens, respiratory allergy induced by trimellitic anhydride (TMA) was characterized by whole genome analysis of lung tissue and blood proteomics in Brown Norway rats. Dermal sensitization (50% and 25% w/v) with TMA and an inhalation challenge of 15 mg/m(3) TMA-induced apneas, laryngeal inflammation, increased numbers of eosinophils, neutrophils and macrophages in bronchoalveolar lavage (BAL), and increased immunoglobulin E levels in serum and lung tissue. Whole genome analysis of lung, sampled 24 hours after challenge, showed expression changes of not only genes belonging to several Gene Ontology groups with up-regulation of inflammatory-associated genes and those associated with lung remodeling but also genes involved in downsizing these processes. Blood proteomics reflected activation of inflammation-inhibiting pathways. Unsensitized animals challenged with TMA exhibited also an increased number of macrophages in BAL, but gene expression in the above-mentioned gene pathways was unchanged or down-regulated. The authors conclude that parameters for lung remodeling can be a valuable tool in hazard identification of LMW respiratory allergens.

Propelling the paradigm shift from reductionism to systems nutrition.

The complex physiology of living organisms represents a challenge for mechanistic understanding of the action of dietary bioactives in the human body and of their possible role in health and disease. Animal, cell, and microbial models have been extensively used to address questions that could not be pursued experimentally in humans, posing an additional level of complexity in translation of the results to healthy and diseased metabolism. The past few decades have witnessed a surge in development of increasingly sensitive molecular techniques and bioinformatic tools for storing, managing, and analyzing increasingly large datasets. Application of such powerful means to molecular nutrition research led to a major leap in study designs and experimental approaches yielding experimental data connecting dietary components to human health. Scientific journals bear major responsibilities in the advancement of science. As primary actors of dissemination to the scientific community, journals can impose rigid criteria for publishing only sound, reliable, and reproducible data. Journal policies are meant to guide potential authors to adopt the most updated standardization guidelines and shared best practices. Such policies evolve in parallel with the evolution of novel approaches and emerging challenges and therefore require constant updating. We highlight in this manuscript the major scientific issues that led to formulating new, updated journal policies for Genes & Nutrition, a journal which targets the growing field of nutritional systems biology interfacing personalized nutrition and preventive medicine, with the ultimate goal of promoting health and preventing or treating disease. We focus here on relevant issues requiring standardization in nutrition research. We also introduce new sections on human genetic variation and nutritional bioinformatics which follow the evolution of nutritional science into the twenty-first century.

Apolipoprotein E genotype status affects habitual human blood mononuclear cell gene expression and its response to fish oil intervention.

SCOPE: People who carry the apolipoprotein E4 (APOE4) single nucleotide polymorphism have an increased risk of cardiovascular disease (CVD). Fish-oil supplementation may help in the prevention of CVD, though interindividual differences in the response to n-3 PUFAs have been observed. We aimed to assess the impact of APOE genotype on peripheral blood mononuclear cell whole genome gene expression at baseline and following a fish-oil intervention. METHODS AND RESULTS: Participants received 6 months of fish-oil supplementation containing 1800 mg of eicosapentaenoic acid and docosahexaenoic acid per day. APOE genotype and peripheral blood mononuclear cell whole genome gene expression before and after supplementation were measured. We characterized the differences in gene expression profiles in carriers of APOE4 (N = 8) compared to noncarriers (N = 15). At baseline, 1320 genes were differentially expressed and the fish-oil supplementation differentially regulated 866 genes between APOE4 carriers and noncarriers. Gene set enrichment analysis showed that carriers had a higher gene expression of cholesterol biosynthesis and IFN signaling pathways. Fish-oil supplementation reduced expression of IFN-related genes in carriers only. CONCLUSION: The increased expression of IFN signaling and cholesterol biosynthesis pathways might explain part of the association between APOE4 and CVD. Fish-oil supplementation may particularly benefit APOE4 carriers by decreasing expression of IFN-related genes.

Network-based integration of molecular and physiological data elucidates regulatory mechanisms underlying adaptation to high-fat diet.

Health is influenced by interplay of molecular, physiological and environmental factors. To effectively maintain health and prevent disease, health-relevant relations need to be understood at multiple levels of biological complexity. Network-based methods provide a powerful platform for integration and mining of data and knowledge characterizing different aspects of health. Previously, we have reported physiological and gene expression changes associated with adaptation of murine epididymal white adipose tissue (eWAT) to 5 days and 12 weeks of high-fat diet (HFD) and low-fat diet feeding (Voigt et al. in Mol Nutr Food Res 57:1423-1434, 2013. doi: 10.1002/mnfr.201200671 ). In the current study, we apply network analysis on this dataset to comprehensively characterize mechanisms driving the short- and long-term adaptation of eWAT to HFD across multiple levels of complexity. We built a three-layered interaction network comprising enriched biological processes, their transcriptional regulators and associated changes in physiological parameters. The multi-layered network model reveals that early eWAT adaptation to HFD feeding involves major changes at a molecular level, including activation of TGF-ß signalling pathway, immune and stress response and downregulation of mitochondrial functioning. Upon prolonged HFD intake, initial transcriptional response tails off, mitochondrial functioning is even further diminished, and in turn the relation between eWAT gene expression and physiological changes becomes more prominent. In particular, eWAT weight and total energy intake negatively correlate with cellular respiration process, revealing mitochondrial dysfunction as a hallmark of late eWAT adaptation to HFD. Apart from global understanding of the time-resolved adaptation to HFD, the multi-layered network model allows several novel mechanistic hypotheses to emerge: (1) early activation of TGF-ß signalling as a trigger for structural and morphological changes in mitochondrial organization in eWAT, (2) modulation of cellular respiration as an intervention strategy to effectively deal with excess dietary fat and (3) discovery of putative intervention targets, such those in pathways related to appetite control. In conclusion, the generated network model comprehensively characterizes eWAT adaptation to high-fat diet, spanning from global aspects to mechanistic details. Being open to further exploration by the research community, it provides a resource of health-relevant interactions ready to be used in a broad range of research applications.

White adipose tissue reference network: a knowledge resource for exploring health-relevant relations.

Optimal health is maintained by interaction of multiple intrinsic and environmental factors at different levels of complexity-from molecular, to physiological, to social. Understanding and quantification of these interactions will aid design of successful health interventions. We introduce the reference network concept as a platform for multi-level exploration of biological relations relevant for metabolic health, by integration and mining of biological interactions derived from public resources and context-specific experimental data. A White Adipose Tissue Health Reference Network (WATRefNet) was constructed as a resource for discovery and prioritization of mechanism-based biomarkers for white adipose tissue (WAT) health status and the effect of food and drug compounds on WAT health status. The WATRefNet (6,797 nodes and 32,171 edges) is based on (1) experimental data obtained from 10 studies addressing different adiposity states, (2) seven public knowledge bases of molecular interactions, (3) expert's definitions of five physiologically relevant processes key to WAT health, namely WAT expandability, Oxidative capacity, Metabolic state, Oxidative stress and Tissue inflammation, and (4) a collection of relevant biomarkers of these processes identified by BIOCLAIMS ( http://bioclaims.uib.es ). The WATRefNet comprehends multiple layers of biological complexity as it contains various types of nodes and edges that represent different biological levels and interactions. We have validated the reference network by showing overrepresentation with anti-obesity drug targets, pathology-associated genes and differentially expressed genes from an external disease model dataset. The resulting network has been used to extract subnetworks specific to the above-mentioned expert-defined physiological processes. Each of these process-specific signatures represents a mechanistically supported composite biomarker for assessing and quantifying the effect of interventions on a physiological aspect that determines WAT health status. Following this principle, five anti-diabetic drug interventions and one diet intervention were scored for the match of their expression signature to the five biomarker signatures derived from the WATRefNet. This confirmed previous observations of successful intervention by dietary lifestyle and revealed WAT-specific effects of drug interventions. The WATRefNet represents a sustainable knowledge resource for extraction of relevant relationships such as mechanisms of action, nutrient intervention targets and biomarkers and for assessment of health effects for support of health claims made on food products.

Network signatures link hepatic effects of anti-diabetic interventions with systemic disease parameters.

BACKGROUND: Multifactorial diseases such as type 2 diabetes mellitus (T2DM), are driven by a complex network of interconnected mechanisms that translate to a diverse range of complications at the physiological level. To optimally treat T2DM, pharmacological interventions should, ideally, target key nodes in this network that act as determinants of disease progression. RESULTS: We set out to discover key nodes in molecular networks based on the hepatic transcriptome dataset from a preclinical study in obese LDLR-/- mice recently published by Radonjic et al. Here, we focus on comparing efficacy of anti-diabetic dietary (DLI) and two drug treatments, namely PPARA agonist fenofibrate and LXR agonist T0901317. By combining knowledge-based and data-driven networks with a random walks based algorithm, we extracted network signatures that link the DLI and two drug interventions to dyslipidemia-related disease parameters. CONCLUSIONS: This study identified specific and prioritized sets of key nodes in hepatic molecular networks underlying T2DM, uncovering pathways that are to be modulated by targeted T2DM drug interventions in order to modulate the complex disease phenotype.

A systems biology approach to understand the pathophysiological mechanisms of cardiac pathological hypertrophy associated with rosiglitazone.

BACKGROUND: Cardiac pathological hypertrophy is associated with a significantly increased risk of coronary heart disease and has been observed in diabetic patients treated with rosiglitazone whereas most published studies do not suggest a similar increase in risk of cardiovascular events in pioglitazone-treated diabetic subjects. This study sought to understand the pathophysiological and molecular mechanisms underlying the disparate cardiovascular effects of rosiglitazone and pioglitazone and yield knowledge as to the causative nature of rosiglitazone-associated cardiac hypertrophy. METHODS: We used a high-fat diet-induced pre-diabetic mouse model to allow bioinformatics analysis of the transcriptome of the heart of mice treated with rosiglitazone or pioglitazone. RESULTS: Our data show that rosiglitazone and pioglitazone both markedly improved systemic markers for glucose homeostasis, fasting plasma glucose and insulin, and the urinary excretion of albumin. Only rosiglitazone, but not pioglitazone, tended to increase atherosclerosis and induced pathological cardiac hypertrophy, based on a significant increase in heart weight and increased expression of the validated markers, ANP and BNP. Functional enrichment analysis of the rosiglitazone-specific cardiac gene expression suggests that a shift in cardiac energy metabolism, in particular decreased fatty acid oxidation toward increased glucose utilization as indicated by down regulation of relevant PPARα and PGC1α target genes. This underlies the rosiglitazone-associated pathological hypertrophic cardiac phenotype in the current study. CONCLUSION: Application of a systems biology approach uncovered a shift in energy metabolism by rosiglitazone that may impact cardiac pathological hypertrophy.

Consensus statement understanding health and malnutrition through a systems approach: the ENOUGH program for early life.

Nutrition research, like most biomedical disciplines, adopted and often uses experimental approaches based on Beadle and Tatum's one gene-one polypeptide hypothesis, thereby reducing biological processes to single reactions or pathways. Systems thinking is needed to understand the complexity of health and disease processes requiring measurements of physiological processes, as well as environmental and social factors, which may alter the expression of genetic information. Analysis of physiological processes with omics technologies to assess systems' responses has only become available over the past decade and remains costly. Studies of environmental and social conditions known to alter health are often not connected to biomedical research. While these facts are widely accepted, developing and conducting comprehensive research programs for health are often beyond financial and human resources of single research groups. We propose a new research program on essential nutrients for optimal underpinning of growth and health (ENOUGH) that will use systems approaches with more comprehensive measurements and biostatistical analysis of the many biological and environmental factors that influence undernutrition. Creating a knowledge base for nutrition and health is a necessary first step toward developing solutions targeted to different populations in diverse social and physical environments for the two billion undernourished people in developed and developing economies.

Delay and impairment in brain development and function in rat offspring after maternal exposure to methylmercury.

Maternal exposure to the neurotoxin methylmercury (MeHg) has been shown to have adverse effects on neural development of the offspring in man. Little is known about the underlying mechanisms by which MeHg affects the developing brain. To explore the neurodevelopmental defects and the underlying mechanism associated with MeHg exposure, the cerebellum and cerebrum of Wistar rat pups were analyzed by [(18)F]FDG PET functional imaging, field potential analysis, and microarray gene expression profiling. Female rat pups were exposed to MeHg via maternal diet during intrauterinal and lactational period (from gestational day 6 to postnatal day (PND)10), and their brain tissues were sampled for the analysis at weaning (PND18-21) and adulthood (PND61-70). The [(18)F]FDG PET imaging and field potential analysis suggested a delay in brain activity and impaired neural function by MeHg. Genome-wide transcriptome analysis substantiated these findings by showing (1) a delay in the onset of gene expression related to neural development, and (2) alterations in pathways related to both structural and functional aspects of nervous system development. The latter included changes in gene expression of developmental regulators, developmental phase-associated genes, small GTPase signaling molecules, and representatives of all processes required for synaptic transmission. These findings were observed at dose levels at which only marginal changes in conventional developmental toxicity endpoints were detected. Therefore, the approaches applied in this study are promising in terms of yielding increased sensitivity compared with classical developmental toxicity tests.

Differential effects of drug interventions and dietary lifestyle in developing type 2 diabetes and complications: a systems biology analysis in LDLr-/- mice.

Excess caloric intake leads to metabolic overload and is associated with development of type 2 diabetes (T2DM). Current disease management concentrates on risk factors of the disease such as blood glucose, however with limited success. We hypothesize that normalizing blood glucose levels by itself is insufficient to reduce the development of T2DM and complications, and that removal of the metabolic overload with dietary interventions may be more efficacious. We explored the efficacy and systems effects of pharmaceutical interventions versus dietary lifestyle intervention (DLI) in developing T2DM and complications. To mimic the situation in humans, high fat diet (HFD)-fed LDLr-/- mice with already established disease phenotype were treated with ten different drugs mixed into HFD or subjected to DLI (switch to low-fat chow), for 7 weeks. Interventions were compared to untreated reference mice kept on HFD or chow only. Although most of the drugs improved HFD-induced hyperglycemia, drugs only partially affected other risk factors and also had limited effect on disease progression towards microalbuminuria, hepatosteatosis and atherosclerosis. By contrast, DLI normalized T2DM risk factors, fully reversed hepatosteatosis and microalbuminuria, and tended to attenuate atherogenesis. The comprehensive beneficial effect of DLI was reflected by normalized metabolite profiles in plasma and liver. Analysis of disease pathways in liver confirmed reversion of the metabolic distortions with DLI. This study demonstrates that the pathogenesis of T2DM towards complications is reversible with DLI and highlights the differential effects of current pharmacotherapies and their limitation to resolve the disease.

Evaluation of multiple variate selection methods from a biological perspective: a nutrigenomics case study.

Genomics-based technologies produce large amounts of data. To interpret the results and identify the most important variates related to phenotypes of interest, various multivariate regression and variate selection methods are used. Although inspected for statistical performance, the relevance of multivariate models in interpreting biological data sets often remains elusive. We compare various multivariate regression and variate selection methods applied to a nutrigenomics data set in terms of performance, utility and biological interpretability. The studied data set comprised hepatic transcriptome (10,072 predictor variates) and plasma protein concentrations [2 dependent variates: Leptin (LEP) and Tissue inhibitor of metalloproteinase 1 (TIMP-1)] collected during a high-fat diet study in ApoE3Leiden mice. The multivariate regression methods used were: partial least squares "PLS"; a genetic algorithm-based multiple linear regression, "GA-MLR"; two least-angle shrinkage methods, "LASSO" and "ELASTIC NET"; and a variant of PLS that uses covariance-based variate selection, "CovProc." Two methods of ranking the genes for Gene Set Enrichment Analysis (GSEA) were also investigated: either by their correlation with the protein data or by the stability of the PLS regression coefficients. The regression methods performed similarly, with CovProc and GA performing the best and worst, respectively (R-squared values based on "double cross-validation" predictions of 0.762 and 0.451 for LEP; and 0.701 and 0.482 for TIMP-1). CovProc, LASSO and ELASTIC NET all produced parsimonious regression models and consistently identified small subsets of variates, with high commonality between the methods. Comparison of the gene ranking approaches found a high degree of agreement, with PLS-based ranking finding fewer significant gene sets. We recommend the use of CovProc for variate selection, in tandem with univariate methods, and the use of correlation-based ranking for GSEA-like pathway analysis methods.