Vamsi K. Mootha

Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles

Although genomewide RNA expression analysis has become a routine tool in biomedical research, extracting biological insight from such information remains a major challenge. Here, we describe a powerful analytical method called Gene Set Enrichment Analysis (GSEA) for interpreting gene expression data. The method derives its power by focusing on gene sets, that is, groups of genes that share common biological function, chromosomal location, or regulation. We demonstrate how GSEA yields insights into several cancer-related data sets, including leukemia and lung cancer. Notably, where single-gene analysis finds little similarity between two independent studies of patient survival in lung cancer, GSEA reveals many biological pathways in common. The GSEA method is embodied in a freely available software package, together with an initial database of 1,325 biologically defined gene sets.

PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes.

DNA microarrays can be used to identify gene expression changes characteristic of human disease. This is challenging, however, when relevant differences are subtle at the level of individual genes. We introduce an analytical strategy, Gene Set Enrichment Analysis, designed to detect modest but coordinate changes in the expression of groups of functionally related genes. Using this approach, we identify a set of genes involved in oxidative phosphorylation whose expression is coordinately decreased in human diabetic muscle. Expression of these genes is high at sites of insulin-mediated glucose disposal, activated by PGC-1α and correlated with total-body aerobic capacity. Our results associate this gene set with clinically important variation in human metabolism and illustrate the value of pathway relationships in the analysis of genomic profiling experiments.

Mechanisms controlling mitochondrial biogenesis and respiration through the thermogenic coactivator PGC-1.

Mitochondrial number and function are altered in response to external stimuli in eukaryotes. While several transcription/replication factors directly regulate mitochondrial genes, the coordination of these factors into a program responsive to the environment is not understood. We show here that PGC-1, a cold-inducible coactivator of nuclear receptors, stimulates mitochondrial biogenesis and respiration in muscle cells through an induction of uncoupling protein 2 (UCP-2) and through regulation of the nuclear respiratory factors (NRFs). PGC-1 stimulates a powerful induction of NRF-1 and NRF-2 gene expression; in addition, PGC-1 binds to and coactivates the transcriptional function of NRF-1 on the promoter for mitochondrial transcription factor A (mtTFA), a direct regulator of mitochondrial DNA replication/transcription. These data elucidate a pathway that directly links external physiological stimuli to the regulation of mitochondrial biogenesis and function.

Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals.

Comprehensive identification of all functional elements encoded in the human genome is a fundamental need in biomedical research. Here, we present a comparative analysis of the human, mouse, rat and dog genomes to create a systematic catalogue of common regulatory motifs in promoters and 3′ untranslated regions (3′ UTRs). The promoter analysis yields 174 candidate motifs, including most previously known transcription-factor binding sites and 105 new motifs. The 3′-UTR analysis yields 106 motifs likely to be involved in post-transcriptional regulation. Nearly one-half are associated with microRNAs (miRNAs), leading to the discovery of many new miRNA genes and their likely target genes. Our results suggest that previous estimates of the number of human miRNA genes were low, and that miRNAs regulate at least 20% of human genes. The overall results provide a systematic view of gene regulation in the human, which will be refined as additional mammalian genomes become available.

Metabolite profiles and the risk of developing diabetes

Emerging technologies allow the high-throughput profiling of metabolic status from a blood specimen (metabolomics). We investigated whether metabolite profiles could predict the development of diabetes. Among 2,422 normoglycemic individuals followed for 12 years, 201 developed diabetes. Amino acids, amines and other polar metabolites were profiled in baseline specimens by liquid chromatography–tandem mass spectrometry (LC-MS). Cases and controls were matched for age, body mass index and fasting glucose. Five branched-chain and aromatic amino acids had highly significant associations with future diabetes: isoleucine, leucine, valine, tyrosine and phenylalanine. A combination of three amino acids predicted future diabetes (with a more than fivefold higher risk for individuals in top quartile). The results were replicated in an independent, prospective cohort. These findings underscore the potential key role of amino acid metabolism early in the pathogenesis of diabetes and suggest that amino acid profiles could aid in diabetes risk assessment.

A mitochondrial protein compendium elucidates complex I disease biology.

Mitochondria are complex organelles whose dysfunction underlies a broad spectrum of human diseases. Identifying all of the proteins resident in this organelle and understanding how they integrate into pathways represent major challenges in cell biology. Toward this goal, we performed mass spectrometry, GFP tagging, and machine learning to create a mitochondrial compendium of 1098 genes and their protein expression across 14 mouse tissues. We link poorly characterized proteins in this inventory to known mitochondrial pathways by virtue of shared evolutionary history. Using this approach, we predict 19 proteins to be important for the function of complex I (CI) of the electron transport chain. We validate a subset of these predictions using RNAi, including C8orf38, which we further show harbors an inherited mutation in a lethal, infantile CI deficiency. Our results have important implications for understanding CI function and pathogenesis and, more generally, illustrate how our compendium can serve as a foundation for systematic investigations of mitochondria.

Defects in Adaptive Energy Metabolism with CNS-Linked Hyperactivity in PGC-1α Null Mice

PGC-1alpha is a coactivator of nuclear receptors and other transcription factors that regulates several metabolic processes, including mitochondrial biogenesis and respiration, hepatic gluconeogenesis, and muscle fiber-type switching. We show here that, while hepatocytes lacking PGC-1alpha are defective in the program of hormone-stimulated gluconeogenesis, the mice have constitutively activated gluconeogenic gene expression that is completely insensitive to normal feeding controls. C/EBPbeta is elevated in the livers of these mice and activates the gluconeogenic genes in a PGC-1alpha-independent manner. Despite having reduced mitochondrial function, PGC-1alpha null mice are paradoxically lean and resistant to diet-induced obesity. This is largely due to a profound hyperactivity displayed by the null animals and is associated with lesions in the striatal region of the brain that controls movement. These data illustrate a central role for PGC-1alpha in the control of energy metabolism but also reveal novel systemic compensatory mechanisms and pathogenic effects of impaired energy homeostasis.

Integrative genomics identifies MCU as an essential component of the mitochondrial calcium uniporter

Mitochondria from diverse organisms are capable of transporting large amounts of Ca2+ via a ruthenium-red-sensitive, membrane-potential-dependent mechanism called the uniporter. Although the uniporter’s biophysical properties have been studied extensively, its molecular composition remains elusive. We recently used comparative proteomics to identify MICU1 (also known as CBARA1), an EF-hand-containing protein that serves as a putative regulator of the uniporter. Here, we use whole-genome phylogenetic profiling, genome-wide RNA co-expression analysis and organelle-wide protein coexpression analysis to predict proteins functionally related to MICU1. All three methods converge on a novel predicted transmembrane protein, CCDC109A, that we now call ‘mitochondrial calcium uniporter’ (MCU). MCU forms oligomers in the mitochondrial inner membrane, physically interacts with MICU1, and resides within a large molecular weight complex. Silencing MCU in cultured cells or in vivo in mouse liver severely abrogates mitochondrial Ca2+ uptake, whereas mitochondrial respiration and membrane potential remain fully intact. MCU has two predicted transmembrane helices, which are separated by a highly conserved linker facing the intermembrane space. Acidic residues in this linker are required for its full activity. However, an S259A point mutation retains function but confers resistance to Ru360, the most potent inhibitor of the uniporter. Our genomic, physiological, biochemical and pharmacological data firmly establish MCU as an essential component of the mitochondrial Ca2+ uniporter.

mTOR controls mitochondrial oxidative function through a YY1-PGC-1α transcriptional complex

Transcriptional complexes that contain peroxisome-proliferator-activated receptor coactivator (PGC)-1α control mitochondrial oxidative function to maintain energy homeostasis in response to nutrient and hormonal signals. An important component in the energy and nutrient pathways is mammalian target of rapamycin (mTOR), a kinase that regulates cell growth, size and survival. However, it is unknown whether and how mTOR controls mitochondrial oxidative activities. Here we show that mTOR is necessary for the maintenance of mitochondrial oxidative function. In skeletal muscle tissues and cells, the mTOR inhibitor rapamycin decreased the gene expression of the mitochondrial transcriptional regulators PGC-1α, oestrogen-related receptor α and nuclear respiratory factors, resulting in a decrease in mitochondrial gene expression and oxygen consumption. Using computational genomics, we identified the transcription factor yin-yang 1 (YY1) as a common target of mTOR and PGC-1α. Knockdown of YY1 caused a significant decrease in mitochondrial gene expression and in respiration, and YY1 was required for rapamycin-dependent repression of those genes. Moreover, mTOR and raptor interacted with YY1, and inhibition of mTOR resulted in a failure of YY1 to interact with and be coactivated by PGC-1α. We have therefore identified a mechanism by which a nutrient sensor (mTOR) balances energy metabolism by means of the transcriptional control of mitochondrial oxidative function. These results have important implications for our understanding of how these pathways might be altered in metabolic diseases and cancer.

Energy Metabolism in Uncoupling Protein 3 Gene Knockout Mice

Uncoupling protein 3 (UCP3) is a member of the mitochondrial anion carrier superfamily. Based upon its high homology with UCP1 and its restricted tissue distribution to skeletal muscle and brown adipose tissue, UCP3 has been suggested to play important roles in regulating energy expenditure, body weight, and thermoregulation. Other postulated roles for UCP3 include regulation of fatty acid metabolism, adaptive responses to acute exercise and starvation, and prevention of reactive oxygen species (ROS) formation. To address these questions, we have generated mice lacking UCP3 (UCP3 knockout (KO) mice). Here, we provide evidence that skeletal muscle mitochondria lacking UCP3 are more coupled (i.e. increased state 3/state 4 ratio), indicating that UCP3 has uncoupling activity. In addition, production of ROS is increased in mitochondria lacking UCP3. This study demonstrates that UCP3 has uncoupling activity and that its absence may lead to increased production of ROS. Despite these effects on mitochondrial function, UCP3 does not seem to be required for body weight regulation, exercise tolerance, fatty acid oxidation, or cold-induced thermogenesis. The absence of such phenotypes in UCP3 KO mice could not be attributed to up-regulation of other UCP mRNAs. However, alternative compensatory mechanisms cannot be excluded. The consequence of increased mitochondrial coupling in UCP3 KO mice on metabolism and the possible role of yet unidentified compensatory mechanisms, remains to be determined.