Automated detection of acute respiratory distress syndrome from chest X-Rays using Directionality Measure and deep learning features.

Acute respiratory distress syndrome (ARDS) is a life-threatening lung injury with global prevalence and high mortality. Chest x-rays (CXR) are critical in the early diagnosis and treatment of ARDS. However, imaging findings may not result in proper identification of ARDS due to a number of reasons, including nonspecific appearance of radiological features, ambiguity in a patient's case due to the pathological stage of the disease, and poor inter-rater reliability from interpretations of CXRs by multiple clinical experts. This study demonstrates the potential capability of methodologies in artificial intelligence, machine learning, and image processing to overcome these challenges and quantitatively assess CXRs for presence of ARDS. We propose and describe Directionality Measure, a novel feature engineering technique used to capture the "cloud-like" appearance of diffuse alveolar damage as a mathematical concept. This study also examines the effectiveness of using an off-the-shelf, pretrained deep learning model as a feature extractor in addition to standard features extracted from the histogram and gray-level co-occurrence matrix (GLCM). Data was collected from hospitalized patients at Michigan Medicine's intensive care unit and the cohort's inclusion criteria was specifically designed to be representative of patients at risk of developing ARDS. Multiple machine learning models were used to evaluate these features with 5-fold cross-validation and the final performance was reported on a hold-out, temporally distinct test set. With AdaBoost, Directionality Measure achieved an accuracy of 78% and AUC of 74% - outperforming classification results using features from the histogram (75% accuracy and 73% AUC), GLCM (76% accuracy and 73% AUC), and ResNet-50 (77% accuracy and 73% AUC). Further experimental results demonstrated that using all feature sets in combination achieved the best overall performance, yielding an accuracy of 83% and AUC of 79% with AdaBoost. These results demonstrate the potential capability of using the proposed methodologies to complement current clinical analysis for detection of ARDS from CXRs.

Valproic acid-induced changes of 4D nuclear morphology in astrocyte cells.

Histone deacetylase inhibitors, such as valproic acid (VPA), have important clinical therapeutic and cellular reprogramming applications. They induce chromatin reorganization that is associated with altered cellular morphology. However, there is a lack of comprehensive characterization of VPA-induced changes of nuclear size and shape. Here, we quantify 3D nuclear morphology of primary human astrocyte cells treated with VPA over time (hence, 4D). We compared volumetric and surface-based representations and identified seven features that jointly discriminate between normal and treated cells with 85% accuracy on day 7. From day 3, treated nuclei were more elongated and flattened and then continued to morphologically diverge from controls over time, becoming larger and more irregular. On day 7, most of the size and shape descriptors demonstrated significant differences between treated and untreated cells, including a 24% increase in volume and 6% reduction in extent (shape regularity) for treated nuclei. Overall, we show that 4D morphometry can capture how chromatin reorganization modulates the size and shape of the nucleus over time. These nuclear structural alterations may serve as a biomarker for histone (de-)acetylation events and provide insights into mechanisms of astrocytes-to-neurons reprogramming.

Robust segmentation of lung in chest x-ray: applications in analysis of acute respiratory distress syndrome.

BACKGROUND: This study outlines an image processing algorithm for accurate and consistent lung segmentation in chest radiographs of critically ill adults and children typically obscured by medical equipment. In particular, this work focuses on applications in analysis of acute respiratory distress syndrome - a critical illness with a mortality rate of 40% that affects 200,000 patients in the United States and 3 million globally each year. METHODS: Chest radiographs were obtained from critically ill adults (n = 100), adults diagnosed with acute respiratory distress syndrome (ARDS) (n = 25), and children (n = 100) hospitalized at Michigan Medicine. Physicians annotated the lung field of each radiograph to establish the ground truth. A Total Variation-based Active Contour (TVAC) lung segmentation algorithm was developed and compared to multiple state-of-the-art methods including a deep learning model (U-Net), a random walker algorithm, and an active spline model, using the Sørensen-Dice coefficient to measure segmentation accuracy. RESULTS: The TVAC algorithm accurately segmented lung fields in all patients in the study. For the adult cohort, an averaged Dice coefficient of 0.86 ±0.04 (min: 0.76) was reported for TVAC, 0.89 ±0.12 (min: 0.01) for U-Net, 0.74 ±0.19 (min: 0.15) for the random walker algorithm, and 0.64 ±0.17 (min: 0.20) for the active spline model. For the pediatric cohort, a Dice coefficient of 0.85 ±0.04 (min: 0.75) was reported for TVAC, 0.87 ±0.09 (min: 0.56) for U-Net, 0.67 ±0.18 (min: 0.18) for the random walker algorithm, and 0.61 ±0.18 (min: 0.18) for the active spline model. CONCLUSION: The proposed algorithm demonstrates the most consistent performance of all segmentation methods tested. These results suggest that TVAC can accurately identify lung fields in chest radiographs in critically ill adults and children.

Druggable Transcriptional Networks in the Human Neurogenic Epigenome.

Chromosome conformation capture methods have revealed the dynamics of genome architecture which is spatially organized into topologically associated domains, with gene regulation mediated by enhancer-promoter pairs in chromatin space. New evidence shows that endogenous hormones and several xenobiotics act within circumscribed topological domains of the spatial genome, impacting subsets of the chromatin contacts of enhancer-gene promoter pairs in cis and trans Results from the National Institutes of Health-funded PsychENCODE project and the study of chromatin remodeling complexes have converged to provide a clearer understanding of the organization of the neurogenic epigenome in humans. Neuropsychiatric diseases, including schizophrenia, bipolar spectrum disorder, autism spectrum disorder, attention deficit hyperactivity disorder, and other neuropsychiatric disorders are significantly associated with mutations in neurogenic transcriptional networks. In this review, we have reanalyzed the results from publications of the PsychENCODE Consortium using pharmacoinformatics network analysis to better understand druggable targets that control neurogenic transcriptional networks. We found that valproic acid and other psychotropic drugs directly alter these networks, including chromatin remodeling complexes, transcription factors, and other epigenetic modifiers. We envision a new generation of CNS therapeutics targeted at neurogenic transcriptional control networks, including druggable parts of chromatin remodeling complexes and master transcription factor-controlled pharmacogenomic networks. This may provide a route to the modification of interconnected gene pathways impacted by disease in patients with neuropsychiatric and neurodegenerative disorders. Direct and indirect therapeutic strategies to modify the master regulators of neurogenic transcriptional control networks may ultimately help extend the life span of CNS neurons impacted by disease.

Publisher Correction: 3D Shape Modeling for Cell Nuclear Morphological Analysis and Classification.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

3D Shape Modeling for Cell Nuclear Morphological Analysis and Classification.

Quantitative analysis of morphological changes in a cell nucleus is important for the understanding of nuclear architecture and its relationship with pathological conditions such as cancer. However, dimensionality of imaging data, together with a great variability of nuclear shapes, presents challenges for 3D morphological analysis. Thus, there is a compelling need for robust 3D nuclear morphometric techniques to carry out population-wide analysis. We propose a new approach that combines modeling, analysis, and interpretation of morphometric characteristics of cell nuclei and nucleoli in 3D. We used robust surface reconstruction that allows accurate approximation of 3D object boundary. Then, we computed geometric morphological measures characterizing the form of cell nuclei and nucleoli. Using these features, we compared over 450 nuclei with about 1,000 nucleoli of epithelial and mesenchymal prostate cancer cells, as well as 1,000 nuclei with over 2,000 nucleoli from serum-starved and proliferating fibroblast cells. Classification of sets of 9 and 15 cells achieved accuracy of 95.4% and 98%, respectively, for prostate cancer cells, and 95% and 98% for fibroblast cells. To our knowledge, this is the first attempt to combine these methods for 3D nuclear shape modeling and morphometry into a highly parallel pipeline workflow for morphometric analysis of thousands of nuclei and nucleoli in 3D.

Deep learning in pharmacogenomics: from gene regulation to patient stratification.

This Perspective provides examples of current and future applications of deep learning in pharmacogenomics, including: identification of novel regulatory variants located in noncoding domains of the genome and their function as applied to pharmacoepigenomics; patient stratification from medical records; and the mechanistic prediction of drug response, targets and their interactions. Deep learning encapsulates a family of machine learning algorithms that has transformed many important subfields of artificial intelligence over the last decade, and has demonstrated breakthrough performance improvements on a wide range of tasks in biomedicine. We anticipate that in the future, deep learning will be widely used to predict personalized drug response and optimize medication selection and dosing, using knowledge extracted from large and complex molecular, epidemiological, clinical and demographic datasets.

The pharmacoepigenomics informatics pipeline defines a pathway of novel and known warfarin pharmacogenomics variants.

AIM: 'Pharmacoepigenomics' methods informed by omics datasets and pre-existing knowledge have yielded discoveries in neuropsychiatric pharmacogenomics. Now we evaluate the generality of these methods by discovering an extended warfarin pharmacogenomics pathway. MATERIALS & METHODS: We developed the pharmacoepigenomics informatics pipeline, a scalable multi-omics variant screening pipeline for pharmacogenomics, and conducted an experiment in the genomics of warfarin. RESULTS: We discovered known and novel pharmacogenomics variants and genes, both coding and regulatory, for warfarin response, including adverse events. Such genes and variants cluster in a warfarin response pathway consolidating known and novel warfarin response variants and genes. CONCLUSION: These results can inform a new warfarin test. The pharmacoepigenomics informatics pipeline may be able to discover new pharmacogenomics markers in other drug-disease systems.

Mining the topography and dynamics of the 4D Nucleome to identify novel CNS drug pathways.

The pharmacoepigenome can be defined as the active, noncoding province of the genome including canonical spatial and temporal regulatory mechanisms of gene regulation that respond to xenobiotic stimuli. Many psychotropic drugs that have been in clinical use for decades have ill-defined mechanisms of action that are beginning to be resolved as we understand the transcriptional hierarchy and dynamics of the nucleus. In this review, we describe spatial, temporal and biomechanical mechanisms mediated by psychotropic medications. Focus is placed on a bioinformatics pipeline that can be used both for detection of pharmacoepigenomic variants that discretize drug response and adverse events to improve pharmacogenomic testing, and for the discovery of novel CNS therapeutics. This approach integrates the functional topology and dynamics of the transcriptional hierarchy of the pharmacoepigenome, gene variant-driven identification of pharmacogenomic regulatory domains, and mesoscale mapping for the discovery of novel CNS pharmacodynamic pathways in human brain. Examples of the application of this pipeline are provided, including the discovery of valproic acid (VPA) mediated transcriptional reprogramming of neuronal cell fate following injury, and mapping of a CNS pathway glutamatergic pathway for the mood stabilizer lithium. These examples in regulatory pharmacoepigenomics illustrate how ongoing research using the 4D nucleome provides a foundation to further insight into previously unrecognized psychotropic drug pharmacodynamic pathways in the human CNS.

Network Reconstruction Reveals that Valproic Acid Activates Neurogenic Transcriptional Programs in Adult Brain Following Traumatic Injury.

OBJECTIVES: To determine the mechanism of action of valproic acid (VPA) in the adult central nervous system (CNS) following traumatic brain injury (TBI) and hemorrhagic shock (HS). METHODS: Data were analyzed from different sources, including experiments in a porcine model, data from postmortem human brain, published studies, public and commercial databases. RESULTS: The transcriptional program in the CNS following TBI, HS, and VPA treatment includes activation of regulatory pathways that enhance neurogenesis and suppress gliogenesis. Genes which encode the transcription factors (TFs) that specify neuronal cell fate, including MEF2D, MYT1L, NEUROD1, PAX6 and TBR1, and their target genes, are induced by VPA. VPA represses genes responsible for oligodendrogenesis, maintenance of white matter, T-cell activation, angiogenesis, and endothelial cell proliferation, adhesion and chemotaxis. NEUROD1 has regulatory interactions with 38% of the genes regulated by VPA in a swine model of TBI and HS in adult brain. Hi-C spatial mapping of a VPA pharmacogenomic SNP in the GRIN2B gene shows it is part of a transcriptional hub that contacts 12 genes that mediate chromatin-mediated neurogenesis and neuroplasticity. CONCLUSIONS: Following TBI and HS, this study shows that VPA administration acts in the adult brain through differential activation of TFs responsible for neurogenesis, genes responsible for neuroplasticity, and repression of TFs that specify oligodendrocyte cell fate, endothelial cell chemotaxis and angiogenesis. Short title: Mechanism of action of valproic acid in traumatic brain injury.

Linking MedDRA(®)-Coded Clinical Phenotypes to Biological Mechanisms by the Ontology of Adverse Events: A Pilot Study on Tyrosine Kinase Inhibitors.

INTRODUCTION: A translational bioinformatics challenge exists in connecting population and individual clinical phenotypes in various formats to biological mechanisms. The Medical Dictionary for Regulatory Activities (MedDRA(®)) is the default dictionary for adverse event (AE) reporting in the US Food and Drug Administration Adverse Event Reporting System (FAERS). The ontology of adverse events (OAE) represents AEs as pathological processes occurring after drug exposures. OBJECTIVES: The aim of this work was to establish a semantic framework to link biological mechanisms to phenotypes of AEs by combining OAE with MedDRA(®) in FAERS data analysis. We investigated the AEs associated with tyrosine kinase inhibitors (TKIs) and monoclonal antibodies (mAbs) targeting tyrosine kinases. The five selected TKIs/mAbs (i.e., dasatinib, imatinib, lapatinib, cetuximab, and trastuzumab) are known to induce impaired ventricular function (non-QT) cardiotoxicity. RESULTS: Statistical analysis of FAERS data identified 1053 distinct MedDRA(®) terms significantly associated with TKIs/mAbs, where 884 did not have corresponding OAE terms. We manually annotated these terms, added them to OAE by the standard OAE development strategy, and mapped them to MedDRA(®). The data integration to provide insights into molecular mechanisms of drug-associated AEs was performed by including linkages in OAE for all related AE terms to MedDRA(®) and the existing ontologies, including the human phenotype ontology (HP), Uber anatomy ontology (UBERON), and gene ontology (GO). Sixteen AEs were shared by all five TKIs/mAbs, and each of 17 cardiotoxicity AEs was associated with at least one TKI/mAb. As an example, we analyzed "cardiac failure" using the relations established in OAE with other ontologies and demonstrated that one of the biological processes associated with cardiac failure maps to the genes associated with heart contraction. CONCLUSION: By expanding the existing OAE ontological design, our TKI use case demonstrated that the combination of OAE and MedDRA(®) provides a semantic framework to link clinical phenotypes of adverse drug events to biological mechanisms.

Epigenomic mapping and effect sizes of noncoding variants associated with psychotropic drug response.

AIM: To provide insight into potential regulatory mechanisms of gene expression underlying addiction, analgesia, psychotropic drug response and adverse drug events, genome-wide association studies searching for variants associated with these phenotypes has been undertaken with limited success. We undertook analysis of these results with the aim of applying epigenetic knowledge to aid variant discovery and interpretation. METHODS: We applied conditional imputation to results from 26 genome-wide association studies and three candidate gene-association studies. The analysis workflow included data from chromatin conformation capture, chromatin state annotation, DNase I hypersensitivity, hypomethylation, anatomical localization and biochronicity. We also made use of chromatin state data from the epigenome roadmap, transcription factor-binding data, spatial maps from published Hi-C datasets and 'guilt by association' methods. RESULTS: We identified 31 pharmacoepigenomic SNPs from a total of 2024 variants in linkage disequilibrium with lead SNPs, of which only 6% were coding variants. Interrogation of chromatin state using our workflow and the epigenome roadmap showed agreement on 34 of 35 tissue assignments to regulatory elements including enhancers and promoters. Loop boundary domains were inferred by association with CTCF (CCCTC-binding factor) and cohesin, suggesting proximity to topologically associating domain boundaries and enhancer clusters. Spatial interactions between enhancer-promoter pairs detected both known and previously unknown mechanisms. Addiction and analgesia SNPs were common in relevant populations and exhibited large effect sizes, whereas a SNP located in the promoter of the SLC1A2 gene exhibited a moderate effect size for lithium response in bipolar disorder in patients of European ancestry. SNPs associated with drug-induced organ injury were rare but exhibited the largest effect sizes, consistent with the published literature. CONCLUSION: This work demonstrates that an in silico bioinformatics-based approach using integrative analysis of a diversity of molecular and morphological data types can discover pharmacoepigenomic variants that are suitable candidates for further validation in cell lines, animal models and human clinical trials.

A glutamatergic network mediates lithium response in bipolar disorder as defined by epigenome pathway analysis.

AIM: A regulatory network in the human brain mediating lithium response in bipolar patients was revealed by analysis of functional SNPs from genome-wide association studies (GWAS) and published gene association studies, followed by epigenome mapping. METHODS: An initial set of 23,312 SNPs in linkage disequilibrium with lead SNPs, and sub-threshold GWAS SNPs rescued by pathway analysis, were studied in the same populations. These were assessed using our workflow and annotation by the epigenome roadmap consortium. RESULTS: Twenty-seven percent of 802 SNPs that were associated with lithium response (13 published studies gene association studies and two GWAS) were shared in common with 1281 SNPs from 18 GWAS examining psychiatric disorders and adverse events associated with lithium treatment. Nineteen SNPs were annotated as active regulatory elements such as enhancers and promoters in a tissue-specific manner. They were located within noncoding regions of ten genes: ANK3, ARNTL, CACNA1C, CACNG2, CDKN1A, CREB1, GRIA2, GSK3B, NR1D1 and SLC1A2. Following gene set enrichment and pathway analysis, these genes were found to be significantly associated (p = 10(-27); Fisher exact test) with an AMPA2 glutamate receptor network in human brain. Our workflow results showed concordance with annotation of regulatory elements from the epigenome roadmap. Analysis of cognate mRNA and enhancer RNA exhibited patterns consistent with an integrated pathway in human brain. CONCLUSION: This pharmacoepigenomic regulatory pathway is located in the same brain regions that exhibit tissue volume loss in bipolar disorder. Although in silico analysis requires biological validation, the approach provides value for identification of candidate variants that may be used in pharmacogenomic testing to identify bipolar patients likely to respond to lithium.

The epigenome, 4D nucleome and next-generation neuropsychiatric pharmacogenomics.

The 4D nucleome has the potential to render challenges in neuropsychiatric pharmacogenomics more tractable. The epigenome roadmap consortium has demonstrated the critical role that noncoding regions of the human genome play in determination of human phenotype. Chromosome conformation capture methods have revealed the 4D organization of the nucleus, bringing interactions between distant regulatory elements into close spatial proximity in a periodic manner. These functional interactions have the potential to elucidate mechanisms of CNS drug response and side effects that previously have been unrecognized. This perspective assesses recent advances likely to reveal novel pharmacodynamic regulatory pathways in human brain, charting a future new avenue of pharmacogenomics research, using the spatial and temporal architecture of the human epigenome as its foundation.

CLO: The cell line ontology.

BACKGROUND: Cell lines have been widely used in biomedical research. The community-based Cell Line Ontology (CLO) is a member of the OBO Foundry library that covers the domain of cell lines. Since its publication two years ago, significant updates have been made, including new groups joining the CLO consortium, new cell line cells, upper level alignment with the Cell Ontology (CL) and the Ontology for Biomedical Investigation, and logical extensions. CONSTRUCTION AND CONTENT: Collaboration among the CLO, CL, and OBI has established consensus definitions of cell line-specific terms such as 'cell line', 'cell line cell', 'cell line culturing', and 'mortal' vs. 'immortal cell line cell'. A cell line is a genetically stable cultured cell population that contains individual cell line cells. The hierarchical structure of the CLO is built based on the hierarchy of the in vivo cell types defined in CL and tissue types (from which cell line cells are derived) defined in the UBERON cross-species anatomy ontology. The new hierarchical structure makes it easier to browse, query, and perform automated classification. We have recently added classes representing more than 2,000 cell line cells from the RIKEN BRC Cell Bank to CLO. Overall, the CLO now contains ~38,000 classes of specific cell line cells derived from over 200 in vivo cell types from various organisms. UTILITY AND DISCUSSION: The CLO has been applied to different biomedical research studies. Example case studies include annotation and analysis of EBI ArrayExpress data, bioassays, and host-vaccine/pathogen interaction. CLO's utility goes beyond a catalogue of cell line types. The alignment of the CLO with related ontologies combined with the use of ontological reasoners will support sophisticated inferencing to advance translational informatics development.

tranSMART: An Open Source and Community-Driven Informatics and Data Sharing Platform for Clinical and Translational Research.

tranSMART is an emerging global open source public private partnership community developing a comprehensive informatics-based analysis and data-sharing cloud platform for clinical and translational research. The tranSMART consortium includes pharmaceutical and other companies, not-for-profits, academic entities, patient advocacy groups, and government stakeholders. The tranSMART value proposition relies on the concept that the global community of users, developers, and stakeholders are the best source of innovation for applications and for useful data. Continued development and use of the tranSMART platform will create a means to enable "pre-competitive" data sharing broadly, saving money and, potentially accelerating research translation to cures. Significant transformative effects of tranSMART includes 1) allowing for all its user community to benefit from experts globally, 2) capturing the best of innovation in analytic tools, 3) a growing 'big data' resource, 4) convergent standards, and 5) new informatics-enabled translational science in the pharma, academic, and not-for-profit sectors.

Ontology-based combinatorial comparative analysis of adverse events associated with killed and live influenza vaccines.

Vaccine adverse events (VAEs) are adverse bodily changes occurring after vaccination. Understanding the adverse event (AE) profiles is a crucial step to identify serious AEs. Two different types of seasonal influenza vaccines have been used on the market: trivalent (killed) inactivated influenza vaccine (TIV) and trivalent live attenuated influenza vaccine (LAIV). Different adverse event profiles induced by these two groups of seasonal influenza vaccines were studied based on the data drawn from the CDC Vaccine Adverse Event Report System (VAERS). Extracted from VAERS were 37,621 AE reports for four TIVs (Afluria, Fluarix, Fluvirin, and Fluzone) and 3,707 AE reports for the only LAIV (FluMist). The AE report data were analyzed by a novel combinatorial, ontology-based detection of AE method (CODAE). CODAE detects AEs using Proportional Reporting Ratio (PRR), Chi-square significance test, and base level filtration, and groups identified AEs by ontology-based hierarchical classification. In total, 48 TIV-enriched and 68 LAIV-enriched AEs were identified (PRR>2, Chi-square score >4, and the number of cases >0.2% of total reports). These AE terms were classified using the Ontology of Adverse Events (OAE), MedDRA, and SNOMED-CT. The OAE method provided better classification results than the two other methods. Thirteen out of 48 TIV-enriched AEs were related to neurological and muscular processing such as paralysis, movement disorders, and muscular weakness. In contrast, 15 out of 68 LAIV-enriched AEs were associated with inflammatory response and respiratory system disorders. There were evidences of two severe adverse events (Guillain-Barre Syndrome and paralysis) present in TIV. Although these severe adverse events were at low incidence rate, they were found to be more significantly enriched in TIV-vaccinated patients than LAIV-vaccinated patients. Therefore, our novel combinatorial bioinformatics analysis discovered that LAIV had lower chance of inducing these two severe adverse events than TIV. In addition, our meta-analysis found that all previously reported positive correlation between GBS and influenza vaccine immunization were based on trivalent influenza vaccines instead of monovalent influenza vaccines.

Researchers' needs for resource discovery and collaboration tools: a qualitative investigation of translational scientists.

BACKGROUND: A critical aspect of clinical and translational science (CTS) is interdisciplinary and collaborative research, which increasingly requires a wide range of computational and human resources. However, few studies have systematically analyzed such resource needs of CTS researchers. OBJECTIVE: To improve our understanding of CTS researchers' needs for computational and human resources in order to build useful and useable supporting informatics tools. METHODS: We conducted semistructured interviews of 30 CTS researchers from the University of Michigan, followed by qualitative analysis of the interview transcripts. RESULTS: The analysis identified three recurring themes: the need for the federation of information, the need to address information overload, and the need to humanize computing, including strong and well-informed views about the use of social networking tools for research collaboration. These findings helped us to narrow down the available design choices for assisting CTS researchers, and helped to identify potential deficiencies of well-known theoretical frameworks used to guide our study, with suggestions for future remedies. CONCLUSIONS: The user needs identified through the study, along with concrete design suggestions, provided key design, methodological, and theoretical insights, which are being used to guide the design and development of a CTS resource portal. The results and interview instrument should be useful to other institutions with Clinical and Translational Science Awards that face similar challenges related to helping CTS researchers make more effective use of computational and human resources.

The NIH National Center for Integrative Biomedical Informatics (NCIBI).

The National Center for Integrative and Biomedical Informatics (NCIBI) is one of the eight NCBCs. NCIBI supports information access and data analysis for biomedical researchers, enabling them to build computational and knowledge models of biological systems to address the Driving Biological Problems (DBPs). The NCIBI DBPs have included prostate cancer progression, organ-specific complications of type 1 and 2 diabetes, bipolar disorder, and metabolic analysis of obesity syndrome. Collaborating with these and other partners, NCIBI has developed a series of software tools for exploratory analysis, concept visualization, and literature searches, as well as core database and web services resources. Many of our training and outreach initiatives have been in collaboration with the Research Centers at Minority Institutions (RCMI), integrating NCIBI and RCMI faculty and students, culminating each year in an annual workshop. Our future directions include focusing on the TranSMART data sharing and analysis initiative.

Integrated metabolome and transcriptome analysis of the NCI60 dataset.

BACKGROUND: Metabolite profiles can be used for identifying molecular signatures and mechanisms underlying diseases since they reflect the outcome of complex upstream genomic, transcriptomic, proteomic and environmental events. The scarcity of publicly accessible large scale metabolome datasets related to human disease has been a major obstacle for assessing the potential of metabolites as biomarkers as well as understanding the molecular events underlying disease-related metabolic changes. The availability of metabolite and gene expression profiles for the NCI-60 cell lines offers the possibility of identifying significant metabolome and transcriptome features and discovering unique molecular processes related to different cancer types. METHODS: We utilized a combination of analytical methods in the R statistical package to evaluate metabolic features associated with cancer cell lines from different tissue origins, identify metabolite-gene correlations and detect outliers cell lines based on metabolome and transcriptome data. Statistical analysis results are integrated with metabolic pathway annotations as well as COSMIC and Tumorscape databases to explore associated molecular mechanisms. RESULTS: Our analysis reveals that although the NCI-60 metabolome dataset is quite noisy comparing with microarray-based transcriptome data, it does contain tissue origin specific signatures. We also identified biologically meaningful gene-metabolite associations. Most remarkably, several abnormal gene-metabolite relationships identified by our approach can be directly linked to known gene mutations and copy number variations in the corresponding cell lines. CONCLUSIONS: Our results suggest that integrative metabolome and transcriptome analysis is a powerful method for understanding molecular machinery underlying various pathophysiological processes. We expect the availability of large scale metabolome data in the coming years will significantly promote the discovery of novel biomarkers, which will in turn improve the understanding of molecular mechanism underlying diseases.