We present the RAVEN (Reconstruction, Analysis and Visualization of Metabolic Networks) Toolbox: a software suite that allows for semi-automated reconstruction of genome-scale models. It makes use of published models and/or the KEGG database, coupled with extensive gap-filling and quality control features. The software suite also contains methods for visualizing simulation results and omics data, as well as a range of methods for performing simulations and analyzing the results. The software is a useful tool for system-wide data analysis in a metabolic context and for streamlined reconstruction of metabolic networks based on protein homology. The RAVEN Toolbox workflow was applied in order to reconstruct a genome-scale metabolic model for the important microbial cell factory Penicillium chrysogenum Wisconsin54-1255. The model was validated in a bibliomic study of in total 440 references, and it comprises 1471 unique biochemical reactions and 1006 ORFs. It was then used to study the roles of ATP and NADPH in the biosynthesis of penicillin, and to identify potential metabolic engineering targets for maximization of penicillin production.

Saccharomyces cerevisiae CEN.PK 113-7D is widely used for metabolic engineering and systems biology research in industry and academia. We sequenced, assembled, annotated and analyzed its genome. Single-nucleotide variations (SNV), insertions/deletions (indels) and differences in genome organization compared to the reference strain S. cerevisiae S288C were analyzed. In addition to a few large deletions and duplications, nearly 3000 indels were identified in the CEN.PK113-7D genome relative to S288C. These differences were overrepresented in genes whose functions are related to transcriptional regulation and chromatin remodelling. Some of these variations were caused by unstable tandem repeats, suggesting an innate evolvability of the corresponding genes. Besides a previously characterized mutation in adenylate cyclase, the CEN.PK113-7D genome sequence revealed a significant enrichment of non-synonymous mutations in genes encoding for components of the cAMP signalling pathway. Some phenotypic characteristics of the CEN.PK113-7D strains were explained by the presence of additional specific metabolic genes relative to S288C. In particular, the presence of the BIO1 and BIO6 genes correlated with a biotin prototrophy of CEN.PK113-7D. Furthermore, the copy number, chromosomal location and sequences of the MAL loci were resolved. The assembled sequence reveals that CEN.PK113-7D has a mosaic genome that combines characteristics of laboratory strains and wild-industrial strains.

Identification of the underlying molecular mechanisms for a derived phenotype by adaptive evolution is difficult. Here, we performed a systems-level inquiry into the metabolic changes occurring in the yeast Saccharomyces cerevisiae as a result of its adaptive evolution to increase its specific growth rate on galactose and related these changes to the acquired phenotypic properties. Three evolved mutants (62A, 62B, and 62C) with higher specific growth rates and faster specific galactose uptake were isolated. The evolved mutants were compared with a reference strain and two engineered strains, SO16 and PGM2, which also showed higher galactose uptake rate in previous studies. The profile of intermediates in galactose metabolism was similar in evolved and engineered mutants, whereas reserve carbohydrates metabolism was specifically elevated in the evolved mutants and one evolved strain showed changes in ergosterol biosynthesis. Mutations were identified in proteins involved in the global carbon sensing Ras/PKA pathway, which is known to regulate the reserve carbohydrates metabolism. We evaluated one of the identified mutations, RAS2 Tyr112, and this mutation resulted in an increased specific growth rate on galactose. These results show that adaptive evolution results in the utilization of unpredicted routes to accommodate increased galactose flux in contrast to rationally engineered strains. Our study demonstrates that adaptive evolution represents a valuable alternative to rational design in bioengineering of improved strains and, that through systems biology, it is possible to identify mutations in evolved strain that can serve as unforeseen metabolic engineering targets for improving microbial strains for production of biofuels and chemicals. I n the field of industrial biotechnology, there is a need to develop efficient cell factories for the production of fuels and chemicals. Even though the concept of metabolic engineering (1) is frequently used in both academia and industry for the development of unique cell factories, evolutionary engineering methods are still widely performed (2). The power of adaptive evolution, sometimes in combination with metabolic engineering, is well illustrated in several recent examples (3, 4). Despite its advantages, conventional random mutagenesis and screening are hampered by the difficulties associated with finding the underlying molecular mechanisms for a derived phenotype and, hence, the combination of adaptive evolution with more rational approaches like metabolic engineering is attractive. Tools from systems biology and the ability to perform deep sequencing of several strains have offered new opportunities for establishing links between genotype and phenotype and, hereby, allow for combinations of random and rational approaches to strain improvement (5, 6).Understanding the evolutionary strategies of a cell to metabolize nonfavored carbon sources is an integral part of strain development in cost efficient bioprocesses. Galactose is an abundant sugar in nonfood crops (7), an...

With the development of next-generation sequencing (NGS) techniques, many software tools have emerged for the discovery of novel microRNAs (miRNAs) and for analyzing the miRNAs expression profiles. An overall evaluation of these diverse software tools is lacking. In this study, we evaluated eight software tools based on their common feature and key algorithms. Three deep-sequencing data sets were collected from different species and used to assess the computational time, sensitivity and accuracy of detecting known miRNAs as well as their capacity for predicting novel miRNAs. Our results provide useful information for researchers to facilitate their selection of the optimal software tools for miRNA analysis depending on their specific requirements, i.e. novel miRNAs discovery or miRNA expression profile analysis of sequencing data sets.