Transcriptional and epigenetic dynamics during specification of human embryonic stem cells.

Differentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

CEBS update: curated toxicology database with enhanced tools for data integration.

The Chemical Effects in Biological Systems database (CEBS) contains extensive toxicology study results and metadata from the Division of the National Toxicology Program (NTP) and other studies of environmental health interest. This resource grants public access to search and collate data from over 10 250 studies for 12 750 test articles (chemicals, environmental agents). CEBS has made considerable strides over the last 5 years to integrate growing internal data repositories into data warehouses and data marts to better serve the public with high quality curated datasets. This effort includes harmonizing legacy terms and metadata to current standards, mapping test articles to external identifiers, and aligning terms to OBO (Open Biological and Biomedical Ontology) Foundry ontologies. The data are made available through the CEBS Homepage (https://cebs.niehs.nih.gov/cebs/), guided search applications, flat files on FTP (file transfer protocol), and APIs (application programming interface) for user access and to provide a bridge for computational tools. The user interface is intuitive with a single search bar to query keywords related to study metadata, publications, and data availability. Results are consolidated to single pages for each test article with NTP conclusions, publications, individual studies, data collections, and links to related test articles and projects available together.

CEBS: a comprehensive annotated database of toxicological data.

The Chemical Effects in Biological Systems database (CEBS) is a comprehensive and unique toxicology resource that compiles individual and summary animal data from the National Toxicology Program (NTP) testing program and other depositors into a single electronic repository. CEBS has undergone significant updates in recent years and currently contains over 11 000 test articles (exposure agents) and over 8000 studies including all available NTP carcinogenicity, short-term toxicity and genetic toxicity studies. Study data provided to CEBS are manually curated, accessioned and subject to quality assurance review prior to release to ensure high quality. The CEBS database has two main components: data collection and data delivery. To accommodate the breadth of data produced by NTP, the CEBS data collection component is an integrated relational design that allows the flexibility to capture any type of electronic data (to date). The data delivery component of the database comprises a series of dedicated user interface tables containing pre-processed data that support each component of the user interface. The user interface has been updated to include a series of nine Guided Search tools that allow access to NTP summary and conclusion data and larger non-NTP datasets. The CEBS database can be accessed online at http://www.niehs.nih.gov/research/resources/databases/cebs/.

Discovery and characterization of artifactual mutations in deep coverage targeted capture sequencing data due to oxidative DNA damage during sample preparation.

As researchers begin probing deep coverage sequencing data for increasingly rare mutations and subclonal events, the fidelity of next generation sequencing (NGS) laboratory methods will become increasingly critical. Although error rates for sequencing and polymerase chain reaction (PCR) are well documented, the effects that DNA extraction and other library preparation steps could have on downstream sequence integrity have not been thoroughly evaluated. Here, we describe the discovery of novel C > A/G > T transversion artifacts found at low allelic fractions in targeted capture data. Characteristics such as sequencer read orientation and presence in both tumor and normal samples strongly indicated a non-biological mechanism. We identified the source as oxidation of DNA during acoustic shearing in samples containing reactive contaminants from the extraction process. We show generation of 8-oxoguanine (8-oxoG) lesions during DNA shearing, present analysis tools to detect oxidation in sequencing data and suggest methods to reduce DNA oxidation through the introduction of antioxidants. Further, informatics methods are presented to confidently filter these artifacts from sequencing data sets. Though only seen in a low percentage of reads in affected samples, such artifacts could have profoundly deleterious effects on the ability to confidently call rare mutations, and eliminating other possible sources of artifacts should become a priority for the research community.

ToxicR: A computational platform in R for computational toxicology and dose-response analyses.

The need to analyze the complex relationships observed in high-throughput toxicogenomic and other omic platforms has resulted in an explosion of methodological advances in computational toxicology. However, advancements in the literature often outpace the development of software researchers can implement in their pipelines, and existing software is frequently based on pre-specified workflows built from well-vetted assumptions that may not be optimal for novel research questions. Accordingly, there is a need for a stable platform and open-source codebase attached to a programming language that allows users to program new algorithms. To fill this gap, the Biostatistics and Computational Biology Branch of the National Institute of Environmental Health Sciences, in cooperation with the National Toxicology Program (NTP) and US Environmental Protection Agency (EPA), developed ToxicR, an open-source R programming package. The ToxicR platform implements many of the standard analyses used by the NTP and EPA, including dose-response analyses for continuous and dichotomous data that employ Bayesian, maximum likelihood, and model averaging methods, as well as many standard tests the NTP uses in rodent toxicology and carcinogenicity studies, such as the poly-K and Jonckheere trend tests. ToxicR is built on the same codebase as current versions of the EPA's Benchmark Dose software and NTP's BMDExpress software but has increased flexibility because it directly accesses this software. To demonstrate ToxicR, we developed a custom workflow to illustrate its capabilities for analyzing toxicogenomic data. The unique features of ToxicR will allow researchers in other fields to add modules, increasing its functionality in the future.

Sources of variance in baseline gene expression in the rodent liver.

The use of gene expression profiling in both clinical and laboratory settings would be enhanced by better characterization of variation due to individual, environmental, and technical factors. Analysis of microarray data from untreated or vehicle-treated animals within the control arm of toxicogenomics studies has yielded useful information on baseline fluctuations in liver gene expression in the rodent. Here, studies which highlight contributions of different factors to gene expression variability in the rodent liver are discussed including a large meta-analysis of rat liver, which identified genes that vary in control animals in the absence of chemical treatment. Genes and their pathways that are the most and least variable were identified in a number of these studies. Life stage, fasting, sex, diet, circadian rhythm and liver lobe source can profoundly influence gene expression in the liver. Recognition of biological and technical factors that contribute to variability of background gene expression can help the investigator in the design of an experiment that maximizes sensitivity and reduces the influence of confounders that may lead to misinterpretation of genomic changes. The factors that contribute to variability in liver gene expression in rodents are likely analogous to those contributing to human interindividual variability in drug response and chemical toxicity. Identification of batteries of genes that are altered in a variety of background conditions could be used to predict responses to drugs and chemicals in appropriate models of the human liver.

Maximizing biomarker discovery by minimizing gene signatures.

BACKGROUND: The use of gene signatures can potentially be of considerable value in the field of clinical diagnosis. However, gene signatures defined with different methods can be quite various even when applied the same disease and the same endpoint. Previous studies have shown that the correct selection of subsets of genes from microarray data is key for the accurate classification of disease phenotypes, and a number of methods have been proposed for the purpose. However, these methods refine the subsets by only considering each single feature, and they do not confirm the association between the genes identified in each gene signature and the phenotype of the disease. We proposed an innovative new method termed Minimize Feature's Size (MFS) based on multiple level similarity analyses and association between the genes and disease for breast cancer endpoints by comparing classifier models generated from the second phase of MicroArray Quality Control (MAQC-II), trying to develop effective meta-analysis strategies to transform the MAQC-II signatures into a robust and reliable set of biomarker for clinical applications. RESULTS: We analyzed the similarity of the multiple gene signatures in an endpoint and between the two endpoints of breast cancer at probe and gene levels, the results indicate that disease-related genes can be preferably selected as the components of gene signature, and that the gene signatures for the two endpoints could be interchangeable. The minimized signatures were built at probe level by using MFS for each endpoint. By applying the approach, we generated a much smaller set of gene signature with the similar predictive power compared with those gene signatures from MAQC-II. CONCLUSIONS: Our results indicate that gene signatures of both large and small sizes could perform equally well in clinical applications. Besides, consistency and biological significances can be detected among different gene signatures, reflecting the studying endpoints. New classifiers built with MFS exhibit improved performance with both internal and external validation, suggesting that MFS method generally reduces redundancies for features within gene signatures and improves the performance of the model. Consequently, our strategy will be beneficial for the microarray-based clinical applications.

CEBS--Chemical Effects in Biological Systems: a public data repository integrating study design and toxicity data with microarray and proteomics data.

CEBS (Chemical Effects in Biological Systems) is an integrated public repository for toxicogenomics data, including the study design and timeline, clinical chemistry and histopathology findings and microarray and proteomics data. CEBS contains data derived from studies of chemicals and of genetic alterations, and is compatible with clinical and environmental studies. CEBS is designed to permit the user to query the data using the study conditions, the subject responses and then, having identified an appropriate set of subjects, to move to the microarray module of CEBS to carry out gene signature and pathway analysis. Scope of CEBS: CEBS currently holds 22 studies of rats, four studies of mice and one study of Caenorhabditis elegans. CEBS can also accommodate data from studies of human subjects. Toxicogenomics studies currently in CEBS comprise over 4000 microarray hybridizations, and 75 2D gel images annotated with protein identification performed by MALDI and MS/MS. CEBS contains raw microarray data collected in accordance with MIAME guidelines and provides tools for data selection, pre-processing and analysis resulting in annotated lists of genes of interest. Additionally, clinical chemistry and histopathology findings from over 1500 animals are included in CEBS. CEBS/BID: The BID (Biomedical Investigation Database) is another component of the CEBS system. BID is a relational database used to load and curate study data prior to export to CEBS, in addition to capturing and displaying novel data types such as PCR data, or additional fields of interest, including those defined by the HESI Toxicogenomics Committee (in preparation). BID has been shared with Health Canada and the US Environmental Protection Agency. CEBS is available at http://cebs.niehs.nih.gov. BID can be accessed via the user interface from https://dir-apps.niehs.nih.gov/arc/. Requests for a copy of BID and for depositing data into CEBS or BID are available at http://www.niehs.nih.gov/cebs-df/.

Transcriptomic analysis of pathways regulated by toll-like receptor 4 in a murine model of chronic pulmonary inflammation and carcinogenesis.

BACKGROUND: Therapeutic strategies exist for human pulmonary neoplasia, however due to the heterogeneity of the disease, most are not very effective. The innate immunity gene, toll-like receptor 4 (TLR4), protects against chronic pulmonary inflammation and tumorigenesis in mice, but the mechanism is unclear. This study was designed to identify TLR4-mediated gene expression pathways that may be used as prognostic indicators of susceptibility to lung tumorigenesis in mice and provide insight into the mechanism. METHODS: Whole lung mRNA was isolated from C.C3H-Tlr4(Lps-d) (BALB(Lps-d); Tlr4 mutant) and BALB/c (Tlr4 normal) mice following butylated hydroxytoluene (BHT)-treatment (four weekly ip. injections; 150-200 mg/kg/each; "promotion"). mRNA from micro-dissected tumors (adenomas) and adjacent uninvolved tissue from both strains were also compared 27 wks after a single carcinogen injection (3-methylcholanthrene (MCA), 10 microg/g; "control") or followed by BHT (6 weekly ip. injections; 125-200 mg/kg/each; "progression"). Bronchoalveolar lavage fluid was analyzed for inflammatory cell content and total protein determination, a marker of lung hyperpermeability; inflammation was also assessed using immunohistochemical staining for macrophages (F4/80) and lymphocytes (CD3) in mice bearing tumors (progression). RESULTS: During promotion, the majority of genes identified in the BALB(Lps-d) compared to BALB/c mice (P < 0.05) were involved in epithelial growth factor receptor (EGFR) signaling (e.g. epiregulin (Ereg)), secreted phosphoprotein 1(Spp1)), which can lead to cell growth and eventual tumor development. Inflammation was significantly higher in BALB(Lps-d) compared to BALB/c mice during progression, similar to the observed response during tumor promotion in these strains. Increases in genes involved in signaling through the EGFR pathway (e.g. Ereg, Spp1) were also observed during progression in addition to continued inflammation, chemotactic, and immune response gene expression in the BALB(Lps-d) versus BALB/c mice (P < 0.05), which appears to provide more favorable conditions for cell growth and tumor development. In support of these findings, the BALB/c mice also had significantly reduced expression of many immune response and inflammatory genes in both the tumors and uninvolved tissue. CONCLUSION: This transcriptomic study determined the protective effect of TLR4 in lung carcinogenesis inhibition of multiple pathways including EGFR (e.g. Ereg), inflammatory response genes (e.g. Cxcl5), chemotaxis (e.g. Ccr1) and other cell proliferation genes (e.g. Arg1, Pthlh). Future studies will determine the utility of these pathways as indicators of immune system deficiencies and tumorigenesis.

Sources of variation in baseline gene expression levels from toxicogenomics study control animals across multiple laboratories.

BACKGROUND: The use of gene expression profiling in both clinical and laboratory settings would be enhanced by better characterization of variance due to individual, environmental, and technical factors. Meta-analysis of microarray data from untreated or vehicle-treated animals within the control arm of toxicogenomics studies could yield useful information on baseline fluctuations in gene expression, although control animal data has not been available on a scale and in a form best served for data-mining. RESULTS: A dataset of control animal microarray expression data was assembled by a working group of the Health and Environmental Sciences Institute's Technical Committee on the Application of Genomics in Mechanism Based Risk Assessment in order to provide a public resource for assessments of variability in baseline gene expression. Data from over 500 Affymetrix microarrays from control rat liver and kidney were collected from 16 different institutions. Thirty-five biological and technical factors were obtained for each animal, describing a wide range of study characteristics, and a subset were evaluated in detail for their contribution to total variability using multivariate statistical and graphical techniques. CONCLUSION: The study factors that emerged as key sources of variability included gender, organ section, strain, and fasting state. These and other study factors were identified as key descriptors that should be included in the minimal information about a toxicogenomics study needed for interpretation of results by an independent source. Genes that are the most and least variable, gender-selective, or altered by fasting were also identified and functionally categorized. Better characterization of gene expression variability in control animals will aid in the design of toxicogenomics studies and in the interpretation of their results.

The first RSBI (ISA-TAB) workshop: "can a simple format work for complex studies?".

This article summarizes the motivation for, and the proceedings of, the first ISA-TAB workshop held December 6-8, 2007, at the EBI, Cambridge, UK. This exploratory workshop, organized by members of the Microarray Gene Expression Data (MGED) Society's Reporting Structure for Biological Investigations (RSBI) working group, brought together a group of developers of a range of collaborative systems to discuss the use of a common format to address the pressing need of reporting and communicating data and metadata from biological, biomedical, and environmental studies employing combinations of genomics, transcriptomics, proteomics, and metabolomics technologies along with more conventional methodologies. The expertise of the participants comprised database development, data management, and hands-on experience in the development of data communication standards. The workshop's outcomes are set to help formalize the proposed Investigation, Study, Assay (ISA)-TAB tab-delimited format for representing and communicating experimental metadata. This article is part of the special issue of OMICS on the activities of the Genomics Standards Consortium (GSC).

Compound cytotoxicity profiling using quantitative high-throughput screening.

BACKGROUND: The propensity of compounds to produce adverse health effects in humans is generally evaluated using animal-based test methods. Such methods can be relatively expensive, low-throughput, and associated with pain suffered by the treated animals. In addition, differences in species biology may confound extrapolation to human health effects. OBJECTIVE: The National Toxicology Program and the National Institutes of Health Chemical Genomics Center are collaborating to identify a battery of cell-based screens to prioritize compounds for further toxicologic evaluation. METHODS: A collection of 1,408 compounds previously tested in one or more traditional toxicologic assays were profiled for cytotoxicity using quantitative high-throughput screening (qHTS) in 13 human and rodent cell types derived from six common targets of xenobiotic toxicity (liver, blood, kidney, nerve, lung, skin). Selected cytotoxicants were further tested to define response kinetics. RESULTS: qHTS of these compounds produced robust and reproducible results, which allowed cross-compound, cross-cell type, and cross-species comparisons. Some compounds were cytotoxic to all cell types at similar concentrations, whereas others exhibited species- or cell type-specific cytotoxicity. Closely related cell types and analogous cell types in human and rodent frequently showed different patterns of cytotoxicity. Some compounds inducing similar levels of cytotoxicity showed distinct time dependence in kinetic studies, consistent with known mechanisms of toxicity. CONCLUSIONS: The generation of high-quality cytotoxicity data on this large library of known compounds using qHTS demonstrates the potential of this methodology to profile a much broader array of assays and compounds, which, in aggregate, may be valuable for prioritizing compounds for further toxicologic evaluation, identifying compounds with particular mechanisms of action, and potentially predicting in vivo biological response.

Towards standards for data exchange and integration and their impact on a public database such as CEBS (Chemical Effects in Biological Systems).

Integration, re-use and meta-analysis of high content study data, typical of DNA microarray studies, can increase its scientific utility. Access to study data and design parameters would enhance the mining of data integrated across studies. However, without standards for which data to include in exchange, and common exchange formats, publication of high content data is time-consuming and often prohibitive. The MGED Society (www.mged.org) was formed in response to the widespread publication of microarray data, and the recognition of the utility of data re-use for meta-analysis. The NIEHS has developed the Chemical Effects in Biological Systems (CEBS) database, which can manage and integrate study data and design from biological and biomedical studies. As community standards are developed for study data and metadata it will become increasingly straightforward to publish high content data in CEBS, where they will be available for meta-analysis. Different exchange formats for study data are being developed: Standard for Exchange of Nonclinical Data (SEND; www.cdisc.org); Tox-ML (www.Leadscope.com) and Simple Investigation Formatted Text (SIFT) from the NIEHS. Data integration can be done at the level of conclusions about responsive genes and phenotypes, and this workflow is supported by CEBS. CEBS also integrates raw and preprocessed data within a given platform. The utility and a method for integrating data within and across DNA microarray studies is shown in an example analysis using DrugMatrix data deposited in CEBS by Iconix Pharmaceuticals.

Toward a checklist for exchange and interpretation of data from a toxicology study.

Data from toxicology and toxicogenomics studies are valuable, and can be combined for meta-analysis using public data repositories such as Chemical Effects in Biological Systems Knowledgebase, ArrayExpress, and Gene Expression Omnibus. In order to fully utilize the data for secondary analysis, it is necessary to have a description of the study and good annotation of the accompanying data. This study annotation permits sophisticated cross-study comparison and analysis, and allows data from comparable subjects to be identified and fully understood. The Minimal Information About a Microarray Experiment Standard was proposed to permit deposition and sharing of microarray data. We propose the first step toward an analogous standard for a toxicogenomics/toxicology study, by describing a checklist of information that best practices would suggest be included with the study data. When the information in this checklist is deposited together with the study data, the checklist information helps the public explore the study data in context of time, or identify data from similarly treated subjects, and also explore/identify potential sources of experimental variability. The proposed checklist summarizes useful information to include when sharing study data for publication, deposition into a database, or electronic exchange with collaborators. It is not a description of how to carry out an experiment, but a definition of how to describe an experiment. It is anticipated that once a toxicology checklist is accepted and put into use, then toxicology databases can be configured to require and output these fields, making it straightforward to annotate data for interpretation by others.

The TGx-28.65 biomarker online application for analysis of transcriptomics data to identify DNA damage-inducing chemicals in human cell cultures.

The TGx-28.65 biomarker is a 65-gene expression profile generated from testing 28 model chemicals (13 that cause DNA damage and 15 that do not) in human TK6 cells. It is used to predict whether a chemical induces DNA damage or not. We expanded availability to the biomarker by developing the online TGx-28.65 biomarker application for predicting the DNA damage inducing (DDI) potential of suspect toxicants tested in p53-proficient human cells and assessing putative mode(s) of action (MOA). Applications like this that analyse gene expression data to predict the hazard potential of test chemicals hold great promise for risk assessment paradigms. The TGx-28.65 biomarker interfaces with an analytical tool to predict the probability that a test chemical can directly or indirectly induce DNA damage. User submitted in vitro microarray data are compared to the 28-chemical x 65-gene signature profile and the probability that the data fit the profile for a DDI or a non-DDI (NDDI) chemical is calculated. The results are displayed in the Results Table, which includes the classification probability and hyperlinks to view heatmaps, hierarchical clustering, and principal component analyses of user-input data in the context of the reference profile. The heatmaps and cluster plots, along with the corresponding text data files of fold changes in gene expression and Euclidean distances can be downloaded. Review of the test chemical data in relationship to the biomarker allows rapid identification of key gene alterations associated with DNA damage as well as chemicals in the reference set that produced a similar response. Environ. Mol. Mutagen. 58:529-535, 2017. © 2017 Wiley Periodicals, Inc.

A strategy capitalizing on synergies: the Reporting Structure for Biological Investigation (RSBI) working group.

In this article we present the Reporting Structure for Biological Investigation (RSBI), a working group under the Microarray Gene Expression Data (MGED) Society umbrella. RSBI brings together several communities to tackle the challenges associated with integrating data and representing complex biological investigations, employing multiple OMICS technologies. Currently, RSBI includes environmental genomics, nutrigenomics and toxicogenomics communities, where independent activities are underway to develop databases and establish data communication standards within their respective domains. The RSBI working group has been conceived as a "single point of focus" for these communities, conforming to general accepted view that duplication and incompatibility should be avoided where possible. This endeavour has aimed to synergize insular solutions into one common terminology between biologically driven standardisation efforts and has also resulted in strong collaborations and shared understanding between those in the technological domain. Through extensive liaisons with many standards efforts, several threads have been woven with the hope that ultimately technology-centered standards and their specific extensions into biological domains of interest will not only stand alone, but will also be able to function together, as interchangeable modules.

Development of FuGO: an ontology for functional genomics investigations.

The development of the Functional Genomics Investigation Ontology (FuGO) is a collaborative, international effort that will provide a resource for annotating functional genomics investigations, including the study design, protocols and instrumentation used, the data generated and the types of analysis performed on the data. FuGO will contain both terms that are universal to all functional genomics investigations and those that are domain specific. In this way, the ontology will serve as the "semantic glue" to provide a common understanding of data from across these disparate data sources. In addition, FuGO will reference out to existing mature ontologies to avoid the need to duplicate these resources, and will do so in such a way as to enable their ease of use in annotation. This project is in the early stages of development; the paper will describe efforts to initiate the project, the scope and organization of the project, the work accomplished to date, and the challenges encountered, as well as future plans.

Exploration of the gene expression correlates of chronic unexplained fatigue using factor analysis.

OBJECTIVE: To identify biomarkers of chronic fatigue syndrome (CFS) and related disorders through analysis of microarray data, pathology test results and self-report symptom profiles. METHOD: To empirically derive the symptom domains of the illnesses, factor analysis was performed on responses to self-report questionnaires (multidimensional fatigue inventory, Centers for Disease Control and Prevention (CDC) symptom inventory and Zung depression scale) before validation with independent datasets. Gene expression patterns that distinguished subjects across each factor dimension were then sought. RESULTS: A four-factor solution was favored, featuring 'fatigue' and 'mood disturbance' factors. Scores on these factors correlated with measures of disability on the Short Form (SF)-36. A total of 57 genes that distinguished subjects along each factor dimension were identified, although the separation was significant only for subjects beyond the extreme (15th and 85th) percentiles of severity. Clustering of laboratory parameters with expression of these genes revealed associations with serum measurements of pH, electrolytes, glucose, urea, creatinine, and liver enzymes (aspartate amino transferase [AST] and alanine amino transferase [AST]); as well as hematocrit and white cell count. CONCLUSION: CFS is a complex syndrome that cannot simply be associated with changes in individual laboratory tests or expression levels of individual genes. No clear association with gene expression and individual symptom domains was found. However, analysis of such multifacetted datasets is likely to be an important means to elucidate the pathogenesis of CFS.

Gene expression profile exploration of a large dataset on chronic fatigue syndrome.

OBJECTIVE: To gain understanding of the molecular basis of chronic fatigue syndrome (CFS) through gene expression analysis using a large microarray data set in conjunction with clinically administrated questionnaires. METHOD: Data from the Wichita (KS, USA) CFS Surveillance Study was used, comprising 167 participants with two self-report questionnaires (multidimensional fatigue inventory [MFI] and Zung depression scale [Zung]), microarray data, empiric classification, and others. Microarray data was analyzed using bioinformatics tools from ArrayTrack. RESULTS: Correspondence analysis was applied to the MFI questionnaire to select the 23 samples having either the most or the least fatigue, and to the Zung questionnaire to select the 26 samples having either the most or least depression; ten samples were common, resulting in a total of 39 samples. The MFI and Zung-based CFS/non-CFS (NF) classifications on the 39 samples were consistent with the empiric classification. Two differentially-expressed gene lists were determined, 188 fatigue-related genes and 164 depression-related genes, which shared 24 common genes and involved 11 common pathways. Principal component analysis based on 24 genes clearly separates 39 samples with respect to their likelihood to be CFS. Most of the 24 genes are not previously reported for CFS, yet their functions are consistent with the prevailing model of CFS, such as immune response, apoptosis, ion channel activity, signal transduction, cell-cell signaling, regulation of cell growth and neuronal activity. Hierarchical cluster analysis was performed based on 24 genes to classify 128 (=167-39) unassigned samples. Several of the 11 identified common pathways are supported by earlier findings for CFS, such as cytokine-cytokine receptor interaction and neuroactive ligand-receptor interaction. Importantly, most of the 11 common pathways are interrelated, suggesting complex biological mechanisms associated with CFS. CONCLUSION: Bioinformatics is critical in this study to select definitive sample groups, analyze gene expression data and gain insight into biological mechanisms. The 24 identified common genes and 11 common pathways could be important in future studies of CFS at the molecular level.

The Ontology for Biomedical Investigations.

The Ontology for Biomedical Investigations (OBI) is an ontology that provides terms with precisely defined meanings to describe all aspects of how investigations in the biological and medical domains are conducted. OBI re-uses ontologies that provide a representation of biomedical knowledge from the Open Biological and Biomedical Ontologies (OBO) project and adds the ability to describe how this knowledge was derived. We here describe the state of OBI and several applications that are using it, such as adding semantic expressivity to existing databases, building data entry forms, and enabling interoperability between knowledge resources. OBI covers all phases of the investigation process, such as planning, execution and reporting. It represents information and material entities that participate in these processes, as well as roles and functions. Prior to OBI, it was not possible to use a single internally consistent resource that could be applied to multiple types of experiments for these applications. OBI has made this possible by creating terms for entities involved in biological and medical investigations and by importing parts of other biomedical ontologies such as GO, Chemical Entities of Biological Interest (ChEBI) and Phenotype Attribute and Trait Ontology (PATO) without altering their meaning. OBI is being used in a wide range of projects covering genomics, multi-omics, immunology, and catalogs of services. OBI has also spawned other ontologies (Information Artifact Ontology) and methods for importing parts of ontologies (Minimum information to reference an external ontology term (MIREOT)). The OBI project is an open cross-disciplinary collaborative effort, encompassing multiple research communities from around the globe. To date, OBI has created 2366 classes and 40 relations along with textual and formal definitions. The OBI Consortium maintains a web resource (http://obi-ontology.org) providing details on the people, policies, and issues being addressed in association with OBI. The current release of OBI is available at http://purl.obolibrary.org/obo/obi.owl.