Graduate Program in Bioinformatics, Boston University, Boston, Massachusetts, United States of America.
PLoS Comput Biol. 2013;9(7):e1003126. doi: 10.1371/journal.pcbi.1003126. Epub 2013 Jul 18.
The filamentous fungus Neurospora crassa played a central role in the development of twentieth-century genetics, biochemistry and molecular biology, and continues to serve as a model organism for eukaryotic biology. Here, we have reconstructed a genome-scale model of its metabolism. This model consists of 836 metabolic genes, 257 pathways, 6 cellular compartments, and is supported by extensive manual curation of 491 literature citations. To aid our reconstruction, we developed three optimization-based algorithms, which together comprise Fast Automated Reconstruction of Metabolism (FARM). These algorithms are: LInear MEtabolite Dilution Flux Balance Analysis (limed-FBA), which predicts flux while linearly accounting for metabolite dilution; One-step functional Pruning (OnePrune), which removes blocked reactions with a single compact linear program; and Consistent Reproduction Of growth/no-growth Phenotype (CROP), which reconciles differences between in silico and experimental gene essentiality faster than previous approaches. Against an independent test set of more than 300 essential/non-essential genes that were not used to train the model, the model displays 93% sensitivity and specificity. We also used the model to simulate the biochemical genetics experiments originally performed on Neurospora by comprehensively predicting nutrient rescue of essential genes and synthetic lethal interactions, and we provide detailed pathway-based mechanistic explanations of our predictions. Our model provides a reliable computational framework for the integration and interpretation of ongoing experimental efforts in Neurospora, and we anticipate that our methods will substantially reduce the manual effort required to develop high-quality genome-scale metabolic models for other organisms.
丝状真菌粗糙脉孢菌(Neurospora crassa)在 20 世纪的遗传学、生物化学和分子生物学的发展中发挥了核心作用,并且仍然是真核生物学的模式生物。在这里,我们重建了它的代谢的基因组规模模型。该模型包含 836 个代谢基因、257 条途径、6 个细胞区室,并且通过对 491 篇文献引用的广泛手动编目得到支持。为了帮助我们进行重建,我们开发了三种基于优化的算法,它们共同构成了快速自动代谢重建(Fast Automated Reconstruction of Metabolism,FARM)。这些算法是:线性代谢物稀释通量平衡分析(LInear MEtabolite Dilution Flux Balance Analysis,limed-FBA),它在线性考虑代谢物稀释的同时预测通量;一步功能修剪(One-step functional Pruning,OnePrune),它使用单个紧凑的线性规划去除被阻塞的反应;以及生长/无生长表型的一致性再现(Consistent Reproduction Of growth/no-growth Phenotype,CROP),它比以前的方法更快地协调了计算机模拟和实验基因必需性之间的差异。针对一组独立的、超过 300 个未用于模型训练的必需/非必需基因的测试集,该模型显示出 93%的灵敏度和特异性。我们还使用该模型模拟了最初在粗糙脉孢菌上进行的生物化学遗传学实验,全面预测了必需基因的营养补救和合成致死相互作用,并提供了我们预测的详细基于途径的机制解释。我们的模型为粗糙脉孢菌正在进行的实验工作的整合和解释提供了一个可靠的计算框架,我们预计我们的方法将大大减少为其他生物开发高质量基因组规模代谢模型所需的人工努力。