Daniel Ziemek

An atlas of genetic influences on human blood metabolites.

Genome-wide association scans with high-throughput metabolic profiling provide unprecedented insights into how genetic variation influences metabolism and complex disease. Here we report the most comprehensive exploration of genetic loci influencing human metabolism thus far, comprising 7,824 adult individuals from 2 European population studies. We report genome-wide significant associations at 145 metabolic loci and their biochemical connectivity with more than 400 metabolites in human blood. We extensively characterize the resulting in vivo blueprint of metabolism in human blood by integrating it with information on gene expression, heritability and overlap with known loci for complex disorders, inborn errors of metabolism and pharmacological targets. We further developed a database and web-based resources for data mining and results visualization. Our findings provide new insights into the role of inherited variation in blood metabolic diversity and identify potential new opportunities for drug development and for understanding disease.

Causal reasoning on biological networks

MOTIVATION The interpretation of high-throughput datasets has remained one of the central challenges of computational biology over the past decade. Furthermore, as the amount of biological knowledge increases, it becomes more and more difficult to integrate this large body of knowledge in a meaningful manner. In this article, we propose a particular solution to both of these challenges. METHODS We integrate available biological knowledge by constructing a network of molecular interactions of a specific kind: causal interactions. The resulting causal graph can be queried to suggest molecular hypotheses that explain the variations observed in a high-throughput gene expression experiment. We show that a simple scoring function can discriminate between a large number of competing molecular hypotheses about the upstream cause of the changes observed in a gene expression profile. We then develop an analytical method for computing the statistical significance of each score. This analytical method also helps assess the effects of random or adversarial noise on the predictive power of our model. RESULTS Our results show that the causal graph we constructed from known biological literature is extremely robust to random noise and to missing or spurious information. We demonstrate the power of our causal reasoning model on two specific examples, one from a cancer dataset and the other from a cardiac hypertrophy experiment. We conclude that causal reasoning models provide a valuable addition to the biologists toolkit for the interpretation of gene expression data. AVAILABILITY AND IMPLEMENTATION R source code for the method is available upon request.

HDAC inhibitors attenuate the development of hypersensitivity in models of neuropathic pain.

Summary Intrathecal delivery of histone deacetylase inhibitors ameliorates hypersensitivity in models of neuropathic pain. This effect may be mediated at the level of the spinal cord through inhibition of HDAC1 function. ABSTRACT Histone deacetylase inhibitors (HDACIs) interfere with the epigenetic process of histone acetylation and are known to have analgesic properties in models of chronic inflammatory pain. The aim of this study was to determine whether these compounds could also affect neuropathic pain. Different class I HDACIs were delivered intrathecally into rat spinal cord in models of traumatic nerve injury and antiretroviral drug–induced peripheral neuropathy (stavudine, d4T). Mechanical and thermal hypersensitivity was attenuated by 40% to 50% as a result of HDACI treatment, but only if started before any insult. The drugs globally increased histone acetylation in the spinal cord, but appeared to have no measurable effects in relevant dorsal root ganglia in this treatment paradigm, suggesting that any potential mechanism should be sought in the central nervous system. Microarray analysis of dorsal cord RNA revealed the signature of the specific compound used (MS‐275) and suggested that its main effect was mediated through HDAC1. Taken together, these data support a role for histone acetylation in the emergence of neuropathic pain.

The Genetic Landscape of Renal Complications in Type 1 Diabetes

Diabetes is the leading cause of ESRD. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2843 subjects, we estimated that the heritability of diabetic kidney disease was 35% (P=6.4×10-3). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole-exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index (P=2.2×10-5) and the risk of type 2 diabetes (P=6.1×10-4) associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation (P=1.1×10-4). Pathway analysis implicated ascorbate and aldarate metabolism (P=9.0×10-6), and pentose and glucuronate interconversions (P=3.0×10-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.

Genome-wide association analysis of self-reported events in 6135 individuals and 252 827 controls identifies 8 loci associated with thrombosis

Thrombotic diseases are among the leading causes of morbidity and mortality in the world. To add insights into the genetic regulation of thrombotic disease, we conducted a genome-wide association study (GWAS) of 6135 self-reported blood clots events and 252 827 controls of European ancestry belonging to the 23andMe cohort of research participants. Eight loci exceeded genome-wide significance. Among the genome-wide significant results, our study replicated previously known venous thromboembolism (VTE) loci near the F5, FGA-FGG, F11, F2, PROCR and ABO genes, and the more recently discovered locus near SLC44A2 In addition, our study reports for the first time a genome-wide significant association between rs114209171, located upstream of the F8 structural gene, and thrombosis risk. Analyses of expression profiles and expression quantitative trait loci across different tissues suggested SLC44A2, ILF3 and AP1M2 as the three most plausible candidate genes for the chromosome 19 locus, our only genome-wide significant thrombosis-related locus that does not harbor likely coagulation-related genes. In addition, we present data showing that this locus also acts as a novel risk factor for stroke and coronary artery disease (CAD). In conclusion, our study reveals novel common genetic risk factors for VTE, stroke and CAD and provides evidence that self-reported data on blood clots used in a GWAS yield results that are comparable with those obtained using clinically diagnosed VTE. This observation opens up the potential for larger meta-analyses, which will enable elucidation of the genetics of thrombotic diseases, and serves as an example for the genetic study of other diseases.

Modeling the Mechanism of Action of a DGAT1 Inhibitor Using a Causal Reasoning Platform

Triglyceride accumulation is associated with obesity and type 2 diabetes. Genetic disruption of diacylglycerol acyltransferase 1 (DGAT1), which catalyzes the final reaction of triglyceride synthesis, confers dramatic resistance to high-fat diet induced obesity. Hence, DGAT1 is considered a potential therapeutic target for treating obesity and related metabolic disorders. However, the molecular events shaping the mechanism of action of DGAT1 pharmacological inhibition have not been fully explored yet. Here, we investigate the metabolic molecular mechanisms induced in response to pharmacological inhibition of DGAT1 using a recently developed computational systems biology approach, the Causal Reasoning Engine (CRE). The CRE algorithm utilizes microarray transcriptomic data and causal statements derived from the biomedical literature to infer upstream molecular events driving these transcriptional changes. The inferred upstream events (also called hypotheses) are aggregated into biological models using a set of analytical tools that allow for evaluation and integration of the hypotheses in context of their supporting evidence. In comparison to gene ontology enrichment analysis which pointed to high-level changes in metabolic processes, the CRE results provide detailed molecular hypotheses to explain the measured transcriptional changes. CRE analysis of gene expression changes in high fat habituated rats treated with a potent and selective DGAT1 inhibitor demonstrate that the majority of transcriptomic changes support a metabolic network indicative of reversal of high fat diet effects that includes a number of molecular hypotheses such as PPARG, HNF4A and SREBPs. Finally, the CRE-generated molecular hypotheses from DGAT1 inhibitor treated rats were found to capture the major molecular characteristics of DGAT1 deficient mice, supporting a phenotype of decreased lipid and increased insulin sensitivity.

Causal reasoning on biological networks: interpreting transcriptional changes

Over the past decade gene expression data sets have been generated at an increasing pace. In addition to ever increasing data generation, the biomedical literature is growing exponentially. The PubMed database (Sayers et al., 2010) comprises more than 20 million citations as of October 2010. The goal of our method is the prediction of putative upstream regulators of observed expression changes based on a set of over 400,000 causal relationships. The resulting putative regulators constitute directly testable hypotheses for follow-up.

Brown Fat Determination and Development from Muscle Precursor Cells by Novel Action of Bone Morphogenetic Protein 6

Brown adipose tissue (BAT) plays a pivotal role in promoting energy expenditure by the virtue of uncoupling protein-1 (UCP-1) that differentiates BAT from its energy storing white adipose tissue (WAT) counterpart. The clinical implication of “classical” BAT (originates from Myf5 positive myoblastic lineage) or the “beige” fat (originates through trans-differentiation of WAT) activation in improving metabolic parameters is now becoming apparent. However, the inducers and endogenous molecular determinants that govern the lineage commitment and differentiation of classical BAT remain obscure. We report here that in the absence of any forced gene expression, stimulation with bone morphogenetic protein 6 (BMP6) induces brown fat differentiation from skeletal muscle precursor cells of murine and human origins. Through a comprehensive transcriptional profiling approach, we have discovered that two days of BMP6 stimulation in C2C12 myoblast cells is sufficient to induce genes characteristic of brown preadipocytes. This developmental switch is modulated in part by newly identified regulators, Optineurin (Optn) and Cyclooxygenase-2 (Cox2). Furthermore, pathway analyses using the Causal Reasoning Engine (CRE) identified additional potential causal drivers of this BMP6 induced commitment switch. Subsequent analyses to decipher key pathway that facilitates terminal differentiation of these BMP6 primed cells identified a key role for Insulin Like Growth Factor-1 Receptor (IGF-1R). Collectively these data highlight a therapeutically innovative role for BMP6 by providing a means to enhance the amount of myogenic lineage derived brown fat.

Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease

Recent advances in highly multiplexed immunoassays have allowed systematic large-scale measurement of hundreds of plasma proteins in large cohort studies. In combination with genotyping, such studies offer the prospect to 1) identify mechanisms involved with regulation of protein expression in plasma, and 2) determine whether the plasma proteins are likely to be causally implicated in disease. We report here the results of genome-wide association (GWA) studies of 83 proteins considered relevant to cardiovascular disease (CVD), measured in 3,394 individuals with multiple CVD risk factors. We identified 79 genome-wide significant (p<5e-8) association signals, 55 of which replicated at P<0.0007 in separate validation studies (n = 2,639 individuals). Using automated text mining, manual curation, and network-based methods incorporating information on expression quantitative trait loci (eQTL), we propose plausible causal mechanisms for 25 trans-acting loci, including a potential post-translational regulation of stem cell factor by matrix metalloproteinase 9 and receptor-ligand pairs such as RANK-RANK ligand. Using public GWA study data, we further evaluate all 79 loci for their causal effect on coronary artery disease, and highlight several potentially causal associations. Overall, a majority of the plasma proteins studied showed evidence of regulation at the genetic level. Our results enable future studies of the causal architecture of human disease, which in turn should aid discovery of new drug targets.

Silencing of enzymes involved in ceramide biosynthesis causes distinct global alterations of lipid homeostasis and gene expression

Dysregulation of ceramide synthesis has been associated with metabolic disorders such as atherosclerosis and diabetes. We examined the changes in lipid homeostasis and gene expression in Huh7 hepatocytes when the synthesis of ceramide is perturbed by knocking down serine pal mitoyltransferase subunits 1, 2, and 3 (SPTLC123) or dihydroceramide desaturase 1 (DEGS1). Although knocking down all SPTLC subunits is necessary to reduce total ceramides significantly, depleting DEGS1 is sufficient to produce a similar outcome. Lipidomic analysis of distribution and speciation of multiple lipid classes indicates an increase in phospholipids in SPTLC123-silenced cells, whereas DEGS1 depletion leads to the accumulation of sphingolipid intermediates, free fatty acids, and diacylglycerol. When cer amide synthesis is disrupted, the transcriptional profiles indicate inhibition in biosynthetic processes, downregulation of genes involved in general endomembrane traffi cking, and upregulation of endocytosis and endosomal recycling. SPTLC123 silencing strongly affects the expression of genes involved with lipid metabolism. Changes in amino acid, sugar, and nucleotide metabolism, as well as vesicle trafficking between organelles, are more prominent in DEGS1-silenced cells. These studies are the first to provide a direct and comprehensive understanding at the lipidomic and transcriptomic levels of how Huh7 hepatocytes respond to changes in the inhibition of ceramide synthesis.