Novel genotyping algorithms for rare variants significantly improve the accuracy of Applied Biosystems™ Axiom™ array genotyping calls: Retrospective evaluation of UK Biobank array data.

The UK Biobank genotyped about 500k participants using Applied Biosystems Axiom microarrays. Participants were subsequently sequenced by the UK Biobank Exome Sequencing Consortium. Axiom genotyping was highly accurate in comparison to sequencing results, for almost 100,000 variants both directly genotyped on the UK Biobank Axiom array and via whole exome sequencing. However, in a study using the exome sequencing results of the first 50k individuals as reference (truth), it was observed that the positive predictive value (PPV) decreased along with the number of heterozygous array calls per variant. We developed a novel addition to the genotyping algorithm, Rare Heterozygous Adjusted (RHA), to significantly improve PPV in variants with minor allele frequency below 0.01%. The improvement in PPV was roughly equal when comparing to the exome sequencing of 50k individuals, or to the more recent ~200k individuals. Sensitivity was higher in the 200k data. The improved calling algorithm, along with enhanced quality control of array probesets, significantly improved the positive predictive value and the sensitivity of array data, making it suitable for the detection of ultra-rare variants.

The liver pharmacological and xenobiotic gene response repertoire.

We have used a supervised classification approach to systematically mine a large microarray database derived from livers of compound-treated rats. Thirty-four distinct signatures (classifiers) for pharmacological and toxicological end points can be identified. Just 200 genes are sufficient to classify these end points. Signatures were enriched in xenobiotic and immune response genes and contain un-annotated genes, indicating that not all key genes in the liver xenobiotic responses have been characterized. Many signatures with equal classification capabilities but with no gene in common can be derived for the same phenotypic end point. The analysis of the union of all genes present in these signatures can reveal the underlying biology of that end point as illustrated here using liver fibrosis signatures. Our approach using the whole genome and a diverse set of compounds allows a comprehensive view of most pharmacological and toxicological questions and is applicable to other situations such as disease and development.

Large-scale integrated databases supporting drug discovery.

Over the past 15 years, genomics, combinatorial chemistry and high-throughput automation have transformed the setting for drug discovery, from an information-poor to a data-rich environment. The next challenge for informatics scientists is to convert the large amount of disparate data produced into useful, integrated information. Consolidation of the different types of information related to drug discovery requires a good working knowledge of database technology, the existence of accepted data standards for achieving uniformity and a complete understanding of the different data systems that are already available. Chemogenomic databases represent the first example of truly integrated systems that make 'omic' technologies directly relevant to small-molecule drug discovery. Researchers within drug discovery programs now have an opportunity to take advantage of new information domains, through the advance and adoption of integrated chemogenomic databases.

Development of a large-scale chemogenomics database to improve drug candidate selection and to understand mechanisms of chemical toxicity and action.

Successful drug discovery requires accurate decision making in order to advance the best candidates from initial lead identification to final approval. Chemogenomics, the use of genomic tools in pharmacology and toxicology, offers a promising enhancement to traditional methods of target identification/validation, lead identification, efficacy evaluation, and toxicity assessment. To realize the value of chemogenomics information, a contextual database is needed to relate the physiological outcomes induced by diverse compounds to the gene expression patterns measured in the same animals. Massively parallel gene expression characterization coupled with traditional assessments of drug candidates provides additional, important mechanistic information, and therefore a means to increase the accuracy of critical decisions. A large-scale chemogenomics database developed from in vivo treated rats provides the context and supporting data to enhance and accelerate accurate interpretation of mechanisms of toxicity and pharmacology of chemicals and drugs. To date, approximately 600 different compounds, including more than 400 FDA approved drugs, 60 drugs approved in Europe and Japan, 25 withdrawn drugs, and 100 toxicants, have been profiled in up to 7 different tissues of rats (representing over 3200 different drug-dose-time-tissue combinations). Accomplishing this task required evaluating and improving a number of in vivo and microarray protocols, including over 80 rigorous quality control steps. The utility of pairing clinical pathology assessments with gene expression data is illustrated using three anti-neoplastic drugs: carmustine, methotrexate, and thioguanine, which had similar effects on the blood compartment, but diverse effects on hepatotoxicity. We will demonstrate that gene expression events monitored in the liver can be used to predict pathological events occurring in that tissue as well as in hematopoietic tissues.