From raw materials to validated system: the construction of a genomic library and microarray to interpret systemic perturbations in Northern bobwhite.

The limited availability of genomic tools and data for nonmodel species impedes computational and systems biology approaches in nonmodel organisms. Here we describe the development, functional annotation, and utilization of genomic tools for the avian wildlife species Northern bobwhite (Colinus virginianus) to determine the molecular impacts of exposure to 2,6-dinitrotoluene (2,6-DNT), a field contaminant of military concern. Massively parallel pyrosequencing of a normalized multitissue library of Northern bobwhite cDNAs yielded 71,384 unique transcripts that were annotated with gene ontology (GO), pathway information, and protein domain analysis. Comparative genome analyses with model organisms revealed functional homologies in 8,825 unique Northern bobwhite genes that are orthologous to 48% of Gallus gallus protein-coding genes. Pathway analysis and GO enrichment of genes differentially expressed in livers of birds exposed for 60 days (d) to 10 and 60 mg/kg/d 2,6-DNT revealed several impacts validated by RT-qPCR including: prostaglandin pathway-mediated inflammation, increased expression of a heme synthesis pathway in response to anemia, and a shift in energy metabolism toward protein catabolism via inhibition of control points for glucose and lipid metabolic pathways, PCK1 and PPARGC1, respectively. This research effort provides the first comprehensive annotated gene library for Northern bobwhite. Transcript expression analysis provided insights into the metabolic perturbations underlying several observed toxicological phenotypes in a 2,6-DNT exposure case study. Furthermore, the systemic impact of dinitrotoluenes on liver function appears conserved across species as PPAR signaling is similarly affected in fathead minnow liver tissue after exposure to 2,4-DNT.

[1]  J. Andrew Royle,et al.  USING THE NORTH AMERICAN BREEDING BIRD SURVEY AS A TOOL FOR CONSERVATION: A CRITIQUE OF BART ET AL. (2004) , 2005 .

[2]  Tao Cai,et al.  Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary , 2005, Bioinform..

[3]  Jianmin Wu,et al.  KOBAS server: a web-based platform for automated annotation and pathway identification , 2006, Nucleic Acids Res..

[4]  H. Ellegren Sequencing goes 454 and takes large‐scale genomics into the wild , 2008, Molecular ecology.

[5]  R. Bartrons,et al.  Pck1 Gene Silencing in the Liver Improves Glycemia Control, Insulin Sensitivity, and Dyslipidemia in db/db Mice , 2008, Diabetes.

[6]  Jade Buchanan-Carter,et al.  Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx , 2009, BMC Genomics.

[7]  Brad T. Sherman,et al.  DAVID: Database for Annotation, Visualization, and Integrated Discovery , 2003, Genome Biology.

[8]  Maciej Szymanski,et al.  Noncoding regulatory RNAs database , 2003, Nucleic Acids Res..

[9]  Matthew A. Bazar,et al.  Sublethal effects of subacute exposure to RDX (1,3,5‐trinitro‐1,3,5‐triazine) in the northern bobwhite (Colinus virginianus) , 2009, Environmental toxicology and chemistry.

[10]  Huiyun Liang,et al.  Staying Current PGC-1 : a key regulator of energy metabolism , 2006 .

[11]  J. Burnside,et al.  Chicken genomics resource: sequencing and annotation of 35,407 ESTs from single and multiple tissue cDNA libraries and CAP3 assembly of a chicken gene index. , 2006, Physiological genomics.

[12]  Gabriel Moreno-Hagelsieb,et al.  Choosing BLAST options for better detection of orthologs as reciprocal best hits , 2008, Bioinform..

[13]  Robert J Kavlock,et al.  Computational toxicology--a state of the science mini review. , 2008, Toxicological sciences : an official journal of the Society of Toxicology.

[14]  H. Leis,et al.  Prostaglandin endoperoxide synthase-2 contributes to the endothelin/sarafotoxin-induced prostaglandin E2 synthesis in mouse osteoblastic cells (MC3T3-E1): evidence for a protein tyrosine kinase-signaling pathway and involvement of protein kinase C. , 1998, Endocrinology.

[15]  M. T. Fernández –Perea,et al.  Influence of avian reproduction ecotoxicological endpoints in the assessment of plant protection products , 2009, Journal of environmental science and health. Part. B, Pesticides, food contaminants, and agricultural wastes.

[16]  Richard D Beger,et al.  Gene expression profiles in fathead minnow exposed to 2,4-DNT: correlation with toxicity in mammals. , 2006, Toxicological sciences : an official journal of the Society of Toxicology.

[17]  P. Greaves Comprar Histopathology of Preclinical Toxicity Studies, 3rd Edition. Interpretation and Relevance in Drug Safety Evaluation | Peter Greaves | 9780444527714 | Elsevier Science , 2007 .

[18]  Huaiyu Mi,et al.  Ontology annotation: mapping genomic regions to biological function. , 2007, Current opinion in chemical biology.

[19]  T. Werner Bioinformatics applications for pathway analysis of microarray data. , 2008, Current opinion in biotechnology.

[20]  Johan Robbens,et al.  Ecotoxicological risk assessment using DNA chips and cellular reporters. , 2007, Trends in biotechnology.

[21]  Feng Liu,et al.  Evolutionary and biomedical implications of a Schistosoma japonicum complementary DNA resource , 2003, Nature Genetics.

[22]  James R. Knight,et al.  Genome sequencing in microfabricated high-density picolitre reactors , 2005, Nature.

[23]  F. M. Cretella,et al.  Nitroaromatic munition compounds: environmental effects and screening values. , 1999, Reviews of environmental contamination and toxicology.

[24]  Robert J. Moore,et al.  Application of chicken microarrays for gene expression analysis in other avian species , 2009, BMC Genomics.

[25]  Youping Deng,et al.  Neurotoxicogenomic investigations to assess mechanisms of action of the munitions constituents RDX and 2,6-DNT in Northern bobwhite (Colinus virginianus). , 2009, Toxicological sciences : an official journal of the Society of Toxicology.

[26]  Yan P. Yuan,et al.  Predicting function: from genes to genomes and back. , 1998, Journal of molecular biology.

[27]  C. V. Jongeneel,et al.  ESTScan: A Program for Detecting, Evaluating, and Reconstructing Potential Coding Regions in EST Sequences , 1999, ISMB.

[28]  Michael I. Jordan,et al.  Protein Molecular Function Prediction by Bayesian Phylogenomics , 2005, PLoS Comput. Biol..

[29]  Christophe Perin,et al.  Phylogenomics of plant genomes: a methodology for genome-wide searches for orthologs in plants , 2008, BMC Genomics.

[30]  O. Ilkayeva,et al.  Metabolic profiling of PPARα−/− mice reveals defects in carnitine and amino acid homeostasis that are partially reversed by oral carnitine supplementation , 2009, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[31]  Ryan A. Rapp,et al.  A Global Assembly of Cotton ESTs , 2016 .

[32]  G. Reddy,et al.  Toxicity of 2,4,6-TrinitrotoIuene (TNT) in Hispid Cotton Rats (Sigmodon hispidus): Hematological, Biochemical, and Pathological Effects , 2000 .

[33]  Alexander R. Pico,et al.  GenMAPP 2: new features and resources for pathway analysis , 2007, BMC Bioinformatics.

[34]  Rolf Apweiler,et al.  InterProScan: protein domains identifier , 2005, Nucleic Acids Res..

[35]  Lin Fang,et al.  WEGO: a web tool for plotting GO annotations , 2006, Nucleic Acids Res..

[36]  Matthew A. Bazar,et al.  Effects of subchronic exposure to 2,6‐dinitrotoluene in the northern bobwhite (Colinus Virginianus) , 2007, Environmental toxicology and chemistry.

[37]  Li Liu,et al.  Construction of a robust microarray from a non‐model species largemouth bass, Micropterus salmoides (Lacèpede), using pyrosequencing technology , 2008 .

[38]  J. Marden,et al.  Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing , 2008, Molecular ecology.

[39]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[40]  W Walter Lorenz,et al.  SAGE Analysis of Transcriptome Responses in Arabidopsis Roots Exposed to 2,4,6-Trinitrotoluene1 , 2003, Plant Physiology.

[41]  J. Steevens,et al.  The role of metabolism in the toxicity of 2,4,6-trinitrotoluene and its degradation products to the aquatic amphipod Hyalella azteca. , 2008, Ecotoxicology and environmental safety.

[42]  Paolo Uva,et al.  An annotated cDNA library and microarray for large-scale gene-expression studies in the ant Solenopsis invicta , 2007, Genome Biology.

[43]  S. Holladay,et al.  Fate and the biochemical effects of 2,4,6-trinitrotoluene exposure to tiger salamanders (Ambystoma tigrinum). , 2000, Ecotoxicology and environmental safety.

[44]  Michael R. Walsh,et al.  Distribution and Fate of Energetics on DoD Test and Training Ranges: Interim Report 5 , 2006 .

[45]  Shoba Ranganathan,et al.  ESTExplorer: an expressed sequence tag (EST) assembly and annotation platform , 2007, Nucleic Acids Res..

[46]  Kiyoko F. Aoki-Kinoshita,et al.  From genomics to chemical genomics: new developments in KEGG , 2005, Nucleic Acids Res..

[47]  Arun Rawat,et al.  Novel implementation of conditional co-regulation by graph theory to derive co-expressed genes from microarray data , 2008, BMC Bioinformatics.

[48]  J. Levine,et al.  Annexin I is a local mediator in neural-endocrine feedback control of inflammation. , 1998, Journal of neurophysiology.

[49]  E. Wright,et al.  Structure of the human Na+/glucose cotransporter gene SGLT1. , 1994, The Journal of biological chemistry.

[50]  Thomas H Hutchinson,et al.  Systems biology meets stress ecology: linking molecular and organismal stress responses in Daphnia magna , 2008, Genome Biology.

[51]  J. Minna,et al.  DNA methylation in health, disease, and cancer. , 2007, Current molecular medicine.

[52]  Rolf Apweiler,et al.  InterProScan - an integration platform for the signature-recognition methods in InterPro , 2001, Bioinform..

[53]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[54]  P. Greaves Liver and Pancreas , 2012 .

[55]  Pierre Baldi,et al.  A Bayesian framework for the analysis of microarray expression data: regularized t -test and statistical inferences of gene changes , 2001, Bioinform..

[56]  G. Elder Genetic defects in the porphyrias: types and significance. , 1998, Clinics in dermatology.

[57]  M. Timmermans,et al.  Collembase: a repository for springtail genomics and soil quality assessment , 2007, BMC Genomics.

[58]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[59]  X. Huang,et al.  CAP3: A DNA sequence assembly program. , 1999, Genome research.

[60]  Gene Ontology Consortium,et al.  The Gene Ontology (GO) project in 2006 , 2005, Nucleic Acids Res..

[61]  P. Greaves 9 – Liver and Pancreas , 2007 .

[62]  Mark S. Johnson,et al.  Subacute toxicity of oral 2,6‐dinitrotoluene and 1,3,5‐trinitro‐1,3,5‐triazine (RDX) exposure to the northern bobwhite (Colinus virginianus) , 2007, Environmental toxicology and chemistry.

[63]  Paul A Clemons,et al.  The Connectivity Map: Using Gene-Expression Signatures to Connect Small Molecules, Genes, and Disease , 2006, Science.

[64]  Steven C. Lawlor,et al.  GenMAPP, a new tool for viewing and analyzing microarray data on biological pathways , 2002, Nature Genetics.

[65]  Matthew A. Bazar,et al.  Influence of Oral 2,4-Dinitrotoluene Exposure to the Northern Bobwhite (Colinus virginianus) , 2005, International journal of toxicology.

[66]  Wu-chun Feng,et al.  The design, implementation, and evaluation of mpiBLAST , 2003 .