Background: Understanding ethanol tolerance in microorganisms is important for the improvement of bioethanol production. Hence, we performed parallel-evolution experiments using Escherichia coli cells under ethanol stress to determine the phenotypic changes necessary for ethanol tolerance. Results: After cultivation of 1,000 generations under 5% ethanol stress, we obtained 6 ethanol-tolerant strains that showed an approximately 2-fold increase in their specific growth rate in comparison with their ancestor. Expression analysis using microarrays revealed that common expression changes occurred during the adaptive evolution to the ethanol stress environment. Biosynthetic pathways of amino acids, including tryptophan, histidine, and branched-chain amino acids, were commonly up-regulated in the tolerant strains, suggesting that activating these pathways is involved in the development of ethanol tolerance. In support of this hypothesis, supplementation of isoleucine, tryptophan, and histidine to the culture medium increased the specific growth rate under ethanol stress. Furthermore, genes related to iron ion metabolism were commonly up-regulated in the tolerant strains, which suggests the change in intracellular redox state during adaptive evolution. Conclusions: The common phenotypic changes in the ethanol-tolerant strains we identified could provide a fundamental basis for designing ethanol-tolerant strains for industrial purposes.

Databases (DBs) are required by various omics fields because the volume of molecular biology data is increasing rapidly. In this study, we provide instructions for users and describe the current status of our metabolite activity DB. To facilitate a comprehensive understanding of the interactions between the metabolites of organisms and the chemical-level contribution of metabolites to human health, we constructed a metabolite activity DB known as the KNApSAcK Metabolite Activity DB. It comprises 9,584 triplet relationships (metabolite–biological activity–target species), including 2,356 metabolites, 140 activity categories, 2,963 specific descriptions of biological activities and 778 target species. Approximately 46% of the activities described in the DB are related to chemical ecology, most of which are attributed to antimicrobial agents and plant growth regulators. The majority of the metabolites with antimicrobial activities are flavonoids and phenylpropanoids. The metabolites with plant growth regulatory effects include plant hormones. Over half of the DB contents are related to human health care and medicine. The five largest groups are toxins, anticancer agents, nervous system agents, cardiovascular agents and non-therapeutic agents, such as flavors and fragrances. The KNApSAcK Metabolite Activity DB is integrated within the KNApSAcK Family DBs to facilitate further systematized research in various omics fields, especially metabolomics, nutrigenomics and foodomics. The KNApSAcK Metabolite Activity DB could also be utilized for developing novel drugs and materials, as well as for identifying viable drug resources and other useful compounds.

Studies on plant metabolites have attracted significant attention in recent years. Over the past 8 years, we have constructed a unique metabolite database, called KNApSAcK, that contains information on the relationships between metabolites and their expressing organism(s). In the present paper, we introduce KNApSAcK-3D, which contains the three-dimensional (3D) structures of all of the metabolic compounds included in the original KNApSAcK database. The 3D structure for each compound was optimized using the Merck Molecular Force Field (MMFF94), and a multiobjective genetic algorithm was used to search extensively for possible conformations and locate the global minimum. The resulting set of structures may be used for docking studies to identify new and potentially unexpected binding sites for target proteins. The 3D structures may also be utilized for more qualitative studies, such as the estimation of biological activities using 3D-QSAR. The database can be accessed via a link from the KNApSAcK Family website (http://kanaya.naist.jp/KNApSAcK_Family/) or directory at http://kanaya.naist.jp/knapsack3d/.

It remains to be determined experimentally whether increasing fitness is related to positive selection, while stationary fitness is related to neutral evolution. Long-term laboratory evolution in Escherichia coli was performed under conditions of thermal stress under defined laboratory conditions. The complete cell growth data showed common continuous fitness recovery to every 2°C or 4°C stepwise temperature upshift, finally resulting in an evolved E. coli strain with an improved upper temperature limit as high as 45.9°C after 523 days of serial transfer, equivalent to 7,560 generations, in minimal medium. Two-phase fitness dynamics, a rapid growth recovery phase followed by a gradual increasing growth phase, was clearly observed at diverse temperatures throughout the entire evolutionary process. Whole-genome sequence analysis revealed the transition from positive to neutral in mutation fixation, accompanied with a considerable escalation of spontaneous substitution rate in the late fitness recovery phase. It suggested that continually increasing fitness not always resulted in the reduction of genetic diversity due to the sequential takeovers by fit mutants, but caused the accumulation of a considerable number of mutations that facilitated the neutral evolution. The detailed results of a two-year in vitro thermal adaptive evolution experiment are described. A laboratory-evolved E. coli strain with an improved upper temperature limit, as high as 45.9°C, was acquired after 523 days of serial transfer, equivalent to 7,560 generations, in nutrient-limited medium. The complete daily records of cell growth exhibited universal two-phase fitness dynamics, a rapid growth recovery phase followed by a gradual increasing growth phase, throughout the entire evolutionary process. Genome-sequence analysis not only showed considerable escalation of the spontaneous substitution rate, but also revealed the transition from positive to nearly neutral in mutation fixation. Particularly, even with the rising fitness of bacterial cells, neutrality was observed in molecular evolution. These observations suggested that a discrete evolutionary mode occurred in the continuous evolutionary route, linking Darwinian adaptive selection with Kimura's neutral evolution. Such transition from beneficial to neutral path may be adopted as an evolutionary strategy robust to rigorous environmental changes.

Next-generation sequencing technologies enable the rapid cost-effective production of sequence data. To evaluate the performance of these sequencing technologies, investigation of the quality of sequence reads obtained from these methods is important. In this study, we analyzed the quality of sequence reads and SNP detection performance using three commercially available next-generation sequencers, i.e., Roche Genome Sequencer FLX System (FLX), Illumina Genome Analyzer (GA), and Applied Biosystems SOLiD system (SOLiD). A common genomic DNA sample obtained from Escherichia coli strain DH1 was applied to these sequencers. The obtained sequence reads were aligned to the complete genome sequence of E. coli DH1, to evaluate the accuracy and sequence bias of these sequence methods. We found that the fraction of “junk” data, which could not be aligned to the reference genome, was largest in the data set of SOLiD, in which about half of reads could not be aligned. Among data sets after alignment to the reference, sequence accuracy was poorest in GA data sets, suggesting relatively low fidelity of the elongation reaction in the GA method. Furthermore, by aligning the sequence reads to the E. coli strain W3110, we screened sequence differences between two E. coli strains using data sets of three different next-generation platforms. The results revealed that the detected sequence differences were similar among these three methods, while the sequence coverage required for the detection was significantly small in the FLX data set. These results provided valuable information on the quality of short sequence reads and the performance of SNP detection in three next-generation sequencing platforms.

To understand how two organisms that have not previously been in contact can establish mutualism, it is first necessary to examine temporal changes in their phenotypes during the establishment of mutualism. Instead of tracing back the history of known, well-established, natural mutualisms, we experimentally simulated the development of mutualism using two genetically-engineered auxotrophic strains of Escherichia coli, which mimic two organisms that have never met before but later establish mutualism. In the development of this synthetic mutualism, one strain, approximately 10 hours after meeting the partner strain, started oversupplying a metabolite essential for the partner's growth, eventually leading to the successive growth of both strains. This cooperative phenotype adaptively appeared only after encountering the partner strain but before the growth of the strain itself. By transcriptome analysis, we found that the cooperative phenotype of the strain was not accompanied by the local activation of the biosynthesis and transport of the oversupplied metabolite but rather by the global activation of anabolic metabolism. This study demonstrates that an organism has the potential to adapt its phenotype after the first encounter with another organism to establish mutualism before its extinction. As diverse organisms inevitably encounter each other in nature, this potential would play an important role in the establishment of a nascent mutualism in nature.