A systems genetics resource and analysis of sleep regulation in the mouse.

Sleep is essential for optimal brain functioning and health, but the biological substrates through which sleep delivers these beneficial effects remain largely unknown. We used a systems genetics approach in the BXD genetic reference population (GRP) of mice and assembled a comprehensive experimental knowledge base comprising a deep "sleep-wake" phenome, central and peripheral transcriptomes, and plasma metabolome data, collected under undisturbed baseline conditions and after sleep deprivation (SD). We present analytical tools to interactively interrogate the database, visualize the molecular networks altered by sleep loss, and prioritize candidate genes. We found that a one-time, short disruption of sleep already extensively reshaped the systems genetics landscape by altering 60%-78% of the transcriptomes and the metabolome, with numerous genetic loci affecting the magnitude and direction of change. Systems genetics integrative analyses drawing on all levels of organization imply α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor trafficking and fatty acid turnover as substrates of the negative effects of insufficient sleep. Our analyses demonstrate that genetic heterogeneity and the effects of insufficient sleep itself on the transcriptome and metabolome are far more widespread than previously reported.

Systems biology of the IMIDIA biobank from organ donors and pancreatectomised patients defines a novel transcriptomic signature of islets from individuals with type 2 diabetes.

AIMS/HYPOTHESIS: Pancreatic islet beta cell failure causes type 2 diabetes in humans. To identify transcriptomic changes in type 2 diabetic islets, the Innovative Medicines Initiative for Diabetes: Improving beta-cell function and identification of diagnostic biomarkers for treatment monitoring in Diabetes (IMIDIA) consortium ( www.imidia.org ) established a comprehensive, unique multicentre biobank of human islets and pancreas tissues from organ donors and metabolically phenotyped pancreatectomised patients (PPP). METHODS: Affymetrix microarrays were used to assess the islet transcriptome of islets isolated either by enzymatic digestion from 103 organ donors (OD), including 84 non-diabetic and 19 type 2 diabetic individuals, or by laser capture microdissection (LCM) from surgical specimens of 103 PPP, including 32 non-diabetic, 36 with type 2 diabetes, 15 with impaired glucose tolerance (IGT) and 20 with recent-onset diabetes (<1 year), conceivably secondary to the pancreatic disorder leading to surgery (type 3c diabetes). Bioinformatics tools were used to (1) compare the islet transcriptome of type 2 diabetic vs non-diabetic OD and PPP as well as vs IGT and type 3c diabetes within the PPP group; and (2) identify transcription factors driving gene co-expression modules correlated with insulin secretion ex vivo and glucose tolerance in vivo. Selected genes of interest were validated for their expression and function in beta cells. RESULTS: Comparative transcriptomic analysis identified 19 genes differentially expressed (false discovery rate ≤0.05, fold change ≥1.5) in type 2 diabetic vs non-diabetic islets from OD and PPP. Nine out of these 19 dysregulated genes were not previously reported to be dysregulated in type 2 diabetic islets. Signature genes included TMEM37, which inhibited Ca2+-influx and insulin secretion in beta cells, and ARG2 and PPP1R1A, which promoted insulin secretion. Systems biology approaches identified HNF1A, PDX1 and REST as drivers of gene co-expression modules correlated with impaired insulin secretion or glucose tolerance, and 14 out of 19 differentially expressed type 2 diabetic islet signature genes were enriched in these modules. None of these signature genes was significantly dysregulated in islets of PPP with impaired glucose tolerance or type 3c diabetes. CONCLUSIONS/INTERPRETATION: These studies enabled the stringent definition of a novel transcriptomic signature of type 2 diabetic islets, regardless of islet source and isolation procedure. Lack of this signature in islets from PPP with IGT or type 3c diabetes indicates differences possibly due to peculiarities of these hyperglycaemic conditions and/or a role for duration and severity of hyperglycaemia. Alternatively, these transcriptomic changes capture, but may not precede, beta cell failure.

Avoiding the pitfalls of gene set enrichment analysis with SetRank.

BACKGROUND: The purpose of gene set enrichment analysis (GSEA) is to find general trends in the huge lists of genes or proteins generated by many functional genomics techniques and bioinformatics analyses. RESULTS: Here we present SetRank, an advanced GSEA algorithm which is able to eliminate many false positive hits. The key principle of the algorithm is that it discards gene sets that have initially been flagged as significant, if their significance is only due to the overlap with another gene set. The algorithm is explained in detail and its performance is compared to that of other methods using objective benchmarking criteria. Furthermore, we explore how sample source bias can affect the results of a GSEA analysis. CONCLUSIONS: The benchmarking results show that SetRank is a highly specific tool for GSEA. Furthermore, we show that the reliability of results can be improved by taking sample source bias into account. SetRank is available as an R package and through an online web interface.

Boolean regulatory network reconstruction using literature based knowledge with a genetic algorithm optimization method.

BACKGROUND: Prior knowledge networks (PKNs) provide a framework for the development of computational biological models, including Boolean models of regulatory networks which are the focus of this work. PKNs are created by a painstaking process of literature curation, and generally describe all relevant regulatory interactions identified using a variety of experimental conditions and systems, such as specific cell types or tissues. Certain of these regulatory interactions may not occur in all biological contexts of interest, and their presence may dramatically change the dynamical behaviour of the resulting computational model, hindering the elucidation of the underlying mechanisms and reducing the usefulness of model predictions. Methods are therefore required to generate optimized contextual network models from generic PKNs. RESULTS: We developed a new approach to generate and optimize Boolean networks, based on a given PKN. Using a genetic algorithm, a model network is built as a sub-network of the PKN and trained against experimental data to reproduce the experimentally observed behaviour in terms of attractors and the transitions that occur between them under specific perturbations. The resulting model network is therefore contextualized to the experimental conditions and constitutes a dynamical Boolean model closer to the observed biological process used to train the model than the original PKN. Such a model can then be interrogated to simulate response under perturbation, to detect stable states and their properties, to get insights into the underlying mechanisms and to generate new testable hypotheses. CONCLUSIONS: Generic PKNs attempt to synthesize knowledge of all interactions occurring in a biological process of interest, irrespective of the specific biological context. This limits their usefulness as a basis for the development of context-specific, predictive dynamical Boolean models. The optimization method presented in this article produces specific, contextualized models from generic PKNs. These contextualized models have improved utility for hypothesis generation and experimental design. The general applicability of this methodological approach makes it suitable for a variety of biological systems and of general interest for biological and medical research. Our method was implemented in the software optimusqual, available online at http://www.vital-it.ch/software/optimusqual/ .

The SwissLipids knowledgebase for lipid biology.

MOTIVATION: Lipids are a large and diverse group of biological molecules with roles in membrane formation, energy storage and signaling. Cellular lipidomes may contain tens of thousands of structures, a staggering degree of complexity whose significance is not yet fully understood. High-throughput mass spectrometry-based platforms provide a means to study this complexity, but the interpretation of lipidomic data and its integration with prior knowledge of lipid biology suffers from a lack of appropriate tools to manage the data and extract knowledge from it. RESULTS: To facilitate the description and exploration of lipidomic data and its integration with prior biological knowledge, we have developed a knowledge resource for lipids and their biology-SwissLipids. SwissLipids provides curated knowledge of lipid structures and metabolism which is used to generate an in silico library of feasible lipid structures. These are arranged in a hierarchical classification that links mass spectrometry analytical outputs to all possible lipid structures, metabolic reactions and enzymes. SwissLipids provides a reference namespace for lipidomic data publication, data exploration and hypothesis generation. The current version of SwissLipids includes over 244 000 known and theoretically possible lipid structures, over 800 proteins, and curated links to published knowledge from over 620 peer-reviewed publications. We are continually updating the SwissLipids hierarchy with new lipid categories and new expert curated knowledge. AVAILABILITY: SwissLipids is freely available at http://www.swisslipids.org/. CONTACT: alan.bridge@isb-sib.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Angiogenic activity of breast cancer patients' monocytes reverted by combined use of systems modeling and experimental approaches.

Angiogenesis plays a key role in tumor growth and cancer progression. TIE-2-expressing monocytes (TEM) have been reported to critically account for tumor vascularization and growth in mouse tumor experimental models, but the molecular basis of their pro-angiogenic activity are largely unknown. Moreover, differences in the pro-angiogenic activity between blood circulating and tumor infiltrated TEM in human patients has not been established to date, hindering the identification of specific targets for therapeutic intervention. In this work, we investigated these differences and the phenotypic reversal of breast tumor pro-angiogenic TEM to a weak pro-angiogenic phenotype by combining Boolean modelling and experimental approaches. Firstly, we show that in breast cancer patients the pro-angiogenic activity of TEM increased drastically from blood to tumor, suggesting that the tumor microenvironment shapes the highly pro-angiogenic phenotype of TEM. Secondly, we predicted in silico all minimal perturbations transitioning the highly pro-angiogenic phenotype of tumor TEM to the weak pro-angiogenic phenotype of blood TEM and vice versa. In silico predicted perturbations were validated experimentally using patient TEM. In addition, gene expression profiling of TEM transitioned to a weak pro-angiogenic phenotype confirmed that TEM are plastic cells and can be reverted to immunological potent monocytes. Finally, the relapse-free survival analysis showed a statistically significant difference between patients with tumors with high and low expression values for genes encoding transitioning proteins detected in silico and validated on patient TEM. In conclusion, the inferred TEM regulatory network accurately captured experimental TEM behavior and highlighted crosstalk between specific angiogenic and inflammatory signaling pathways of outstanding importance to control their pro-angiogenic activity. Results showed the successful in vitro reversion of such an activity by perturbation of in silico predicted target genes in tumor derived TEM, and indicated that targeting tumor TEM plasticity may constitute a novel valid therapeutic strategy in breast cancer.

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.

ExPASy: SIB bioinformatics resource portal.

ExPASy (http://www.expasy.org) has worldwide reputation as one of the main bioinformatics resources for proteomics. It has now evolved, becoming an extensible and integrative portal accessing many scientific resources, databases and software tools in different areas of life sciences. Scientists can henceforth access seamlessly a wide range of resources in many different domains, such as proteomics, genomics, phylogeny/evolution, systems biology, population genetics, transcriptomics, etc. The individual resources (databases, web-based and downloadable software tools) are hosted in a 'decentralized' way by different groups of the SIB Swiss Institute of Bioinformatics and partner institutions. Specifically, a single web portal provides a common entry point to a wide range of resources developed and operated by different SIB groups and external institutions. The portal features a search function across 'selected' resources. Additionally, the availability and usage of resources are monitored. The portal is aimed for both expert users and people who are not familiar with a specific domain in life sciences. The new web interface provides, in particular, visual guidance for newcomers to ExPASy.

Meta-Analysis of Human Cancer Single-Cell RNA-Seq Datasets Using the IMMUcan Database.

The development of single-cell RNA sequencing (scRNA-seq) technologies has greatly contributed to deciphering the tumor microenvironment (TME). An enormous amount of independent scRNA-seq studies have been published representing a valuable resource that provides opportunities for meta-analysis studies. However, the massive amount of biological information, the marked heterogeneity and variability between studies, and the technical challenges in processing heterogeneous datasets create major bottlenecks for the full exploitation of scRNA-seq data. We have developed IMMUcan scDB (https://immucanscdb.vital-it.ch), a fully integrated scRNA-seq database exclusively dedicated to human cancer and accessible to nonspecialists. IMMUcan scDB encompasses 144 datasets on 56 different cancer types, annotated in 50 fields containing precise clinical, technological, and biological information. A data processing pipeline was developed and organized in four steps: (i) data collection; (ii) data processing (quality control and sample integration); (iii) supervised cell annotation with a cell ontology classifier of the TME; and (iv) interface to analyze TME in a cancer type-specific or global manner. This framework was used to explore datasets across tumor locations in a gene-centric (CXCL13) and cell-centric (B cells) manner as well as to conduct meta-analysis studies such as ranking immune cell types and genes correlated to malignant transformation. This integrated, freely accessible, and user-friendly resource represents an unprecedented level of detailed annotation, offering vast possibilities for downstream exploitation of human cancer scRNA-seq data for discovery and validation studies. SIGNIFICANCE: The IMMUcan scDB database is an accessible supportive tool to analyze and decipher tumor-associated single-cell RNA sequencing data, allowing researchers to maximally use this data to provide new insights into cancer biology.

Transcriptomic Signature Differences Between SARS-CoV-2 and Influenza Virus Infected Patients.

The reason why most individuals with COVID-19 have relatively limited symptoms while other develop respiratory distress with life-threatening complications remains unknown. Increasing evidence suggests that COVID-19 associated adverse outcomes mainly rely on dysregulated immunity. Here, we compared transcriptomic profiles of blood cells from 103 patients with different severity levels of COVID-19 with that of 27 healthy and 22 influenza-infected individuals. Data provided a complete overview of SARS-CoV-2-induced immune signature, including a dramatic defect in IFN responses, a reduction of toxicity-related molecules in NK cells, an increased degranulation of neutrophils, a dysregulation of T cells, a dramatic increase in B cell function and immunoglobulin production, as well as an important over-expression of genes involved in metabolism and cell cycle in patients infected with SARS-CoV-2 compared to those infected with influenza viruses. These features also differed according to COVID-19 severity. Overall and specific gene expression patterns across groups can be visualized on an interactive website (https://bix.unil.ch/covid/). Collectively, these transcriptomic host responses to SARS-CoV-2 infection are discussed in the context of current studies, thereby improving our understanding of COVID-19 pathogenesis and shaping the severity level of COVID-19.

MIMAS 3.0 is a Multiomics Information Management and Annotation System.

BACKGROUND: DNA sequence integrity, mRNA concentrations and protein-DNA interactions have been subject to genome-wide analyses based on microarrays with ever increasing efficiency and reliability over the past fifteen years. However, very recently novel technologies for Ultra High-Throughput DNA Sequencing (UHTS) have been harnessed to study these phenomena with unprecedented precision. As a consequence, the extensive bioinformatics environment available for array data management, analysis, interpretation and publication must be extended to include these novel sequencing data types. DESCRIPTION: MIMAS was originally conceived as a simple, convenient and local Microarray Information Management and Annotation System focused on GeneChips for expression profiling studies. MIMAS 3.0 enables users to manage data from high-density oligonucleotide SNP Chips, expression arrays (both 3'UTR and tiling) and promoter arrays, BeadArrays as well as UHTS data using MIAME-compliant standardized vocabulary. Importantly, researchers can export data in MAGE-TAB format and upload them to the EBI's ArrayExpress certified data repository using a one-step procedure. CONCLUSION: We have vastly extended the capability of the system such that it processes the data output of six types of GeneChips (Affymetrix), two different BeadArrays for mRNA and miRNA (Illumina) and the Genome Analyzer (a popular Ultra-High Throughput DNA Sequencer, Illumina), without compromising on its flexibility and user-friendliness. MIMAS, appropriately renamed into Multiomics Information Management and Annotation System, is currently used by scientists working in approximately 50 academic laboratories and genomics platforms in Switzerland and France. MIMAS 3.0 is freely available via http://multiomics.sourceforge.net/.

DynaStI: A Dynamic Retention Time Database for Steroidomics.

: Steroidomics studies face the challenge of separating analytical compounds with very similar structures (i.e., isomers). Liquid chromatography (LC) is commonly used to this end, but the shared core structure of this family of compounds compromises effective separations among the numerous chemical analytes with comparable physico-chemical properties. Careful tuning of the mobile phase gradient and an appropriate choice of the stationary phase can be used to overcome this problem, in turn modifying the retention times in different ways for each compound. In the usual workflow, this approach is suboptimal for the annotation of features based on retention times since it requires characterizing a library of known compounds for every fine-tuned configuration. We introduce a software solution, DynaStI, that is capable of annotating liquid chromatography-mass spectrometry (LC-MS) features by dynamically generating the retention times from a database containing intrinsic properties of a library of metabolites. DynaStI uses the well-established linear solvent strength (LSS) model for reversed-phase LC. Given a list of LC-MS features and some characteristics of the LC setup, this software computes the corresponding retention times for the internal database and then annotates the features using the exact masses with predicted retention times at the working conditions. DynaStI (https://dynasti.vital-it.ch) is able to automatically calibrate its predictions to compensate for deviations in the input parameters. The database also includes identification and structural information for each annotation, such as IUPAC name, CAS number, SMILES string, metabolic pathways, and links to external metabolomic or lipidomic databases.

CATMA, a comprehensive genome-scale resource for silencing and transcript profiling of Arabidopsis genes.

BACKGROUND: The Complete Arabidopsis Transcript MicroArray (CATMA) initiative combines the efforts of laboratories in eight European countries 1 to deliver gene-specific sequence tags (GSTs) for the Arabidopsis research community. The CATMA initiative offers the power and flexibility to regularly update the GST collection according to evolving knowledge about the gene repertoire. These GST amplicons can easily be reamplified and shared, subsets can be picked at will to print dedicated arrays, and the GSTs can be cloned and used for other functional studies. This ongoing initiative has already produced approximately 24,000 GSTs that have been made publicly available for spotted microarray printing and RNA interference. RESULTS: GSTs from the CATMA version 2 repertoire (CATMAv2, created in 2002) were mapped onto the gene models from two independent Arabidopsis nuclear genome annotation efforts, TIGR5 and PSB-EuGène, to consolidate a list of genes that were targeted by previously designed CATMA tags. A total of 9,027 gene models were not tagged by any amplified CATMAv2 GST, and 2,533 amplified GSTs were no longer predicted to tag an updated gene model. To validate the efficacy of GST mapping criteria and design rules, the predicted and experimentally observed hybridization characteristics associated to GST features were correlated in transcript profiling datasets obtained with the CATMAv2 microarray, confirming the reliability of this platform. To complete the CATMA repertoire, all 9,027 gene models for which no GST had yet been designed were processed with an adjusted version of the Specific Primer and Amplicon Design Software (SPADS). A total of 5,756 novel GSTs were designed and amplified by PCR from genomic DNA. Together with the pre-existing GST collection, this new addition constitutes the CATMAv3 repertoire. It comprises 30,343 unique amplified sequences that tag 24,202 and 23,009 protein-encoding nuclear gene models in the TAIR6 and EuGène genome annotations, respectively. To cover the remaining untagged genes, we identified 543 additional GSTs using less stringent design criteria and designed 990 sequence tags matching multiple members of gene families (Gene Family Tags or GFTs) to cover any remaining untagged genes. These latter 1,533 features constitute the CATMAv4 addition. CONCLUSION: To update the CATMA GST repertoire, we designed 7,289 additional sequence tags, bringing the total number of tagged TAIR6-annotated Arabidopsis nuclear protein-coding genes to 26,173. This resource is used both for the production of spotted microarrays and the large-scale cloning of hairpin RNA silencing vectors. All information about the resulting updated CATMA repertoire is available through the CATMA database http://www.catma.org.

Arabidopsis jasmonate signaling pathway.

Jasmonates control defense gene expression and male fertility in the model plant Arabidopsis thaliana. In both cases, the involvement of the jasmonate pathway is complex, involving large-scale transcriptional reprogramming. Additionally, jasmonate signaling is hard-wired into the auxin, ethylene, and salicylate signal networks, all of which are under intense investigation in Arabidopsis. In male fertility, jasmonic acid (JA) is the essential signal intervening both at the level of anther elongation and in pollen dehiscense. A number of genes potentially involved in jasmonate-dependent anther elongation have recently been discovered. In the case of defense, at least two jasmonates, JA and its precursor 12-oxo-phytodienoic acid (OPDA), are necessary for the fine-tuning of defense gene expression in response to various microbial pathogens and arthropod herbivores. However, only OPDA is required for full resistance to some insects and fungi. Other jasmonates probably affect yet more physiological responses. A series of breakthroughs have identified the SKP/CULLIN/F-BOX (SCF), CORONATINE INSENSITIVE (COI1) complex, acting together with the CONSTITUTIVE PHOTOMORPHOGENIC 9 (COP9) signalosome, as central regulatory components of jasmonate signaling in Arabidopsis. The studies, mostly involving mutational approaches, have paved the way for suppressor screens that are expected to further extend our knowledge of jasmonate signaling. When these and other new mutants affecting jasmonate signaling are characterized, new nodes will be added to the Arabidopsis Jasmonate Signaling Pathway Connections Map, and the lists of target genes regulated by jasmonates in Arabidopsis will be expanded.

Jasmonate signaling pathway.

Jasmonates in plants are cyclic fatty acid-derived regulators structurally similar to prostaglandins in metazoans. These chemicals mediate many of plants' transcriptional responses to wounding and pathogenesis by acting as potent regulators for the expression of numerous frontline immune response genes, including those for defensins and antifungal proteins. Additionally, the pathway is critical for fertility. Ongoing genetic screens and protein-protein interaction assays are identifying components of the canonical jasmonate signaling pathway. A massive molecular machine, based on two multiprotein complexes, SCF(COI1) and the COP9 signalosome (CNS), plays a central role in jasmonate signaling. This machine functions in vivo as a ubiquitin ligase complex, probably targeting regulatory proteins, some of which are expected to be transcriptional repressors. Some defense-related mediators, notably salicylic acid, antagonize jasmonates in controlling the expression of many genes. In Arabidopsis, NONEXPRESSOR OF PR GENES (NPR1) mediates part of this interaction, with another layer of control provided further downstream by the mitogen-activated protein kinase (MAPK) homolog MPK4. Numerous other interpathway connections influence the jasmonate pathway. Insights from Arabidopsis have shown that an allele of the auxin signaling gene AXR1, for example, reduces the sensitivity of plants to jasmonate. APETALA2 (AP2)-domain transcription factors, such as ETHYLENE RESPONSE FACTOR 1 (ERF1), link the jasmonate pathway to the ethylene signaling pathway. As progress in characterizing several new mutants (some of which are hypersensitive to jasmonic acid) augments our understanding of jasmonate signaling, the Connections Map will be updated to include this new information.

Jasmonate biochemical pathway.

Plants possess an interrelated family of potent fatty acid-derived regulators-the jasmonates. These compounds, which play roles in both defense and development, are derived from tri-unsaturated fatty acids [alpha-linolenic acid (18:3) or 7Z,10Z,13Z-hexadecatrienoic acid (16:3)]. The lipoxygenase-catalyzed addition of molecular oxygen to alpha-linolenic acid initiates jasmonate synthesis by providing a 13-hydroperoxide substrate for the formation of an unstable allene oxide that is then subject to enzyme-guided cyclization to produce 12-oxo-phytodienoic acid (OPDA). OPDA has several fates, including esterification into plastid lipids or transformation into the 12-carbon co-regulator jasmonic acid (JA). JA, the best-characterized member of the family, regulates both male and female fertility (depending on the plant species) and is an important mediator of defense gene expression. JA is itself a substrate for further diverse modifications. Genetic dissection of the pathway is revealing how the different jasmonates modulate different physiological processes. Each new family member that is discovered provides another key to understanding the fine control of gene expression in immune responses, in the initiation and maintenance of long-distance signal transfer in response to wounding, in the regulation of fertility, and in the turnover, inactivation, and sequestration of jasmonates, among other processes. The Jasmonate Biochemical Pathway provides an overview of the growing jasmonate family, and new members will be included in future versions of the Connections Map.

Plasma Dihydroceramides Are Diabetes Susceptibility Biomarker Candidates in Mice and Humans.

Plasma metabolite concentrations reflect the activity of tissue metabolic pathways and their quantitative determination may be informative about pathogenic conditions. We searched for plasma lipid species whose concentrations correlate with various parameters of glucose homeostasis and susceptibility to type 2 diabetes (T2D). Shotgun lipidomic analysis of the plasma of mice from different genetic backgrounds, which develop a pre-diabetic state at different rates when metabolically stressed, led to the identification of a group of sphingolipids correlated with glucose tolerance and insulin secretion. Quantitative analysis of these and closely related lipids in the plasma of individuals from two population-based prospective cohorts revealed that specific long-chain fatty-acid-containing dihydroceramides were significantly elevated in the plasma of individuals who will progress to diabetes up to 9 years before disease onset. These lipids may serve as early biomarkers of, and help identify, metabolic deregulation in the pathogenesis of T2D.

Molecular phenotyping of multiple mouse strains under metabolic challenge uncovers a role for Elovl2 in glucose-induced insulin secretion.

OBJECTIVE: In type 2 diabetes (T2D), pancreatic ß cells become progressively dysfunctional, leading to a decline in insulin secretion over time. In this study, we aimed to identify key genes involved in pancreatic beta cell dysfunction by analyzing multiple mouse strains in parallel under metabolic stress. METHODS: Male mice from six commonly used non-diabetic mouse strains were fed a high fat or regular chow diet for three months. Pancreatic islets were extracted and phenotypic measurements were recorded at 2 days, 10 days, 30 days, and 90 days to assess diabetes progression. RNA-Seq was performed on islet tissue at each time-point and integrated with the phenotypic data in a network-based analysis. RESULTS: A module of co-expressed genes was selected for further investigation as it showed the strongest correlation to insulin secretion and oral glucose tolerance phenotypes. One of the predicted network hub genes was Elovl2, encoding Elongase of very long chain fatty acids 2. Elovl2 silencing decreased glucose-stimulated insulin secretion in mouse and human ß cell lines. CONCLUSION: Our results suggest a role for Elovl2 in ensuring normal insulin secretory responses to glucose. Moreover, the large comprehensive dataset and integrative network-based approach provides a new resource to dissect the molecular etiology of ß cell failure under metabolic stress.

Identifying biological mechanisms for favorable cancer prognosis using non-hypothesis-driven iterative survival analysis.

Survival analyses based on the Kaplan-Meier estimate have been pervasively used to support or validate the relevance of biological mechanisms in cancer research. Recently, with the appearance of gene expression high-throughput technologies, this kind of analysis has been applied to tumor transcriptomics data. In a 'bottom-up' approach, gene-expression profiles that are associated with a deregulated pathway hypothetically involved in cancer progression are first identified and then subsequently correlated with a survival effect, which statistically supports or requires the rejection of such a hypothesis. In this work, we propose a 'top-down' approach, in which the clinical outcome (survival) is the starting point that guides the identification of deregulated biological mechanisms in cancer by a non-hypothesis-driven iterative survival analysis. We show that the application of our novel method to a population of ~2,000 breast cancer patients of the METABRIC consortium allows the identification of several well-known cancer mechanisms, such as ERBB4, HNF3A and TGFB pathways, and the investigation of their paradoxical dual effect. In addition, several novel biological mechanisms are proposed as potentially involved in cancer progression. The proposed exploratory methodology can be considered both alternative and complementary to classical 'bottom-up' approaches for validation of biological hypotheses. We propose that our method may be used to better characterize cancer, and may therefore impact the future design of therapies that are truly molecularly tailored to individual patients. The method, named SURCOMED, was implemented as a web-based tool, which is publicly available at http://surcomed.vital-it.ch. R scripts are also available at http://surcomed.sourceforge.net).