Leif Väremo

Gene set analysis (GSA) is used to elucidate genome-wide data, in particular transcriptome data. A multitude of methods have been proposed for this step of the analysis, and many of them have been compared and evaluated. Unfortunately, there is no consolidated opinion regarding what methods should be preferred, and the variety of available GSA software and implementations pose a difficulty for the end-user who wants to try out different methods. To address this, we have developed the R package Piano that collects a range of GSA methods into the same system, for the benefit of the end-user. Further on we refine the GSA workflow by using modifications of the gene-level statistics. This enables us to divide the resulting gene set P-values into three classes, describing different aspects of gene expression directionality at gene set level. We use our fully implemented workflow to investigate the impact of the individual components of GSA by using microarray and RNA-seq data. The results show that the evaluated methods are globally similar and the major separation correlates well with our defined directionality classes. As a consequence of this, we suggest to use a consensus scoring approach, based on multiple GSA runs. In combination with the directionality classes, this constitutes a more thorough basis for an enriched biological interpretation.

Pdgfra+ oligodendrocyte precursor cells (OPCs) arise in distinct specification waves during embryogenesis in the central nervous system (CNS). It is unclear whether there is a correlation between these waves and different oligodendrocyte (OL) states at adult stages. Here, we present bulk and single-cell transcriptomics resources providing insights on how transitions between these states occur. We found that post-natal OPCs from brain and spinal cord present similar transcriptional signatures. Moreover, post-natal OPC progeny of E13.5 Pdgfra+ cells present electrophysiological and transcriptional profiles similar to OPCs derived from subsequent specification waves, indicating that Pdgfra+ pre-OPCs rewire their transcriptional network during development. Single-cell RNA-seq and lineage tracing indicates that a subset of E13.5 Pdgfra+ cells originates cells of the pericyte lineage. Thus, our results indicate that embryonic Pdgfra+ cells in the CNS give rise to distinct post-natal cell lineages, including OPCs with convergent transcriptional profiles in different CNS regions.

BACKGROUND: Succinate dehydrogenase (SDH) is a mitochondrial metabolic enzyme complex involved in both the electron transport chain and the citric acid cycle. SDH mutations resulting in enzymatic dysfunction have been found to be a predisposing factor in various hereditary cancers. Therefore, SDH has been implicated as a tumor suppressor. RESULTS: We identified that dysregulation of SDH components also occurs in serous ovarian cancer, particularly the SDH subunit SDHB. Targeted knockdown of Sdhb in mouse ovarian cancer cells resulted in enhanced proliferation and an epithelial-to-mesenchymal transition (EMT). Bioinformatics analysis revealed that decreased SDHB expression leads to a transcriptional upregulation of genes involved in metabolic networks affecting histone methylation. We confirmed that Sdhb knockdown leads to a hypermethylated epigenome that is sufficient to promote EMT. Metabolically, the loss of Sdhb resulted in reprogrammed carbon source utilization and mitochondrial dysfunction. This altered metabolic state of Sdhb knockdown cells rendered them hypersensitive to energy stress. CONCLUSIONS: These data illustrate how SDH dysfunction alters the epigenetic and metabolic landscape in ovarian cancer. By analyzing the involvement of this enzyme in transcriptional and metabolic networks, we find a metabolic Achilles' heel that can be exploited therapeutically. Analyses of this type provide an understanding how specific perturbations in cancer metabolism may lead to novel anticancer strategies.

Skeletal myocytes are metabolically active and susceptible to insulin resistance and are thus implicated in type 2 diabetes (T2D). This complex disease involves systemic metabolic changes, and their elucidation at the systems level requires genome-wide data and biological networks. Genome-scale metabolic models (GEMs) provide a network context for the integration of high-throughput data. We generated myocyte-specific RNA-sequencing data and investigated their correlation with proteome data. These data were then used to reconstruct a comprehensive myocyte GEM. Next, we performed a meta-analysis of six studies comparing muscle transcription in T2D versus healthy subjects. Transcriptional changes were mapped on the myocyte GEM, revealing extensive transcriptional regulation in T2D, particularly around pyruvate oxidation, branched-chain amino acid catabolism, and tetrahydrofolate metabolism, connected through the downregulated dihydrolipoamide dehydrogenase. Strikingly, the gene signature underlying this metabolic regulation successfully classifies the disease state of individual samples, suggesting that regulation of these pathways is a ubiquitous feature of myocytes in response to T2D.

BACKGROUND: To understand cardiac and skeletal muscle function, it is important to define and explore their molecular constituents and also to identify similarities and differences in the gene expression in these two different striated muscle tissues. Here, we have investigated the genes and proteins with elevated expression in cardiac and skeletal muscle in relation to all other major human tissues and organs using a global transcriptomics analysis complemented with antibody-based profiling to localize the corresponding proteins on a single cell level. RESULTS: Our study identified a comprehensive list of genes expressed in cardiac and skeletal muscle. The genes with elevated expression were further stratified according to their global expression pattern across the human body as well as their precise localization in the muscle tissues. The functions of the proteins encoded by the elevated genes are well in line with the physiological functions of cardiac and skeletal muscle, such as contraction, ion transport, regulation of membrane potential and actomyosin structure organization. A large fraction of the transcripts in both cardiac and skeletal muscle correspond to mitochondrial proteins involved in energy metabolism, which demonstrates the extreme specialization of these muscle tissues to provide energy for contraction. CONCLUSIONS: Our results provide a comprehensive list of genes and proteins elevated in striated muscles. A number of proteins not previously characterized in cardiac and skeletal muscle were identified and localized to specific cellular subcompartments. These proteins represent an interesting starting point for further functional analysis of their role in muscle biology and disease.