The availability of genome sequences, annotations and knowledge of the biochemistry underlying metabolic transformations has led to the generation of metabolic network reconstructions for a wide range of organisms in bacteria, archaea, and eukaryotes. When modeled using mathematical representations, a reconstruction can simulate underlying genotype-phenotype relationships. Accordingly, genome-scale models (GEMs) can be used to predict the response of organisms to genetic and environmental variations. A bottom-up reconstruction procedure typically starts by generating a draft model from existing annotation data on a target organism. For model species, this part of the process can be straightforward, due to the abundant organism-specific biochemical data. However, the process becomes complicated for non-model less-annotated species. In this paper, we present a draft liver reconstruction, ReCodLiver0.9, of Atlantic cod (Gadus morhua), a non-model teleost fish, as a practicable guide for cases with comparably few resources. Although the reconstruction is considered a draft version, we show that it already has utility in elucidating metabolic response mechanisms to environmental toxicants by mapping gene expression data of exposure experiments to the resulting model.
Author summaryGenome-scale metabolic models (GEMs) are constructed based upon reconstructed networks that are carried out by an organism. The underlying biochemical knowledge in such networks can be transformed into mathematical models that could serve as a platform to answer biological questions. The availability of high-throughput biological data, including genomics, proteomics, and metabolomics data, supports the generation of such models for a large number of organisms. Nevertheless, challenges arise for non-model species which are typically less annotated. In this paper, we discuss these June 23, 2020 1/14 challenges and possible solutions in the context of generation of a draft liver reconstruction of Atlantic cod (Gadus morhua). We also show how experimental data, here gene expression data, can be mapped to the resulting model to understand the metabolic response of cod liver to environmental toxicants. 2 high-throughput genomic, proteomic, and metabolomic data, genome-scale metabolic 3 models (GEMs) have become fundamental tools in the systems biology of metabolism. 4 To develop a GEM, metabolic genes and their product enzymes are assembled into a 5 network of metabolic reactions carried out by an organism [42]. Such networks can be 6 used as a platform to answer relevant biological questions, by transforming biochemical 7 knowledge into a mathematical format and computing appropriate physiological 8 states [41, 46]. Using Boolean logic representation, gene-protein-reaction (GPR) 9 associations depict direct connections between genotype and metabolic capability [39]. 10 GEMs have a wide scope of applications, among which are the contextualization of 11 omics data, hypothesis-driven discovery, support for metabolic engineering, modeling 12 interact...