A global regulatory role of gluconeogenic genes in Escherichia coli revealed by transcriptome network analysis

Efficient Fermentative Production of β-Alanine from Glucose through Multidimensional Engineering of Escherichia coli

β-Alanine, a valuable β-type amino acid, is experiencing increased demand due to its multifaceted applications in food flavoring, nutritional supplements, pharmaceuticals, and the chemical industry. Nevertheless, the sustainable biosynthesis of β-alanine currently faces challenges due to the scarcity of robust strains, attributed to the complexities of modulating multiple genes and the inherent physiological constraints. Here, systems metabolic engineering was implemented in Escherichia coli to overcome these limitations. First, an efficient l-aspartate-α-decarboxylase (ADC) was recruited for β-alanine biosynthesis. To conserve phosphoenolpyruvate flux, we subsequently modified the endogenous glucose assimilation system by inactivating the phosphotransferase system (PTS) and introducing an alternative non-PTS system, which increased β-alanine production to 1.70 g/L. The supply of key precursors, oxaloacetate and l-aspartate, was synergistically improved through comprehensive modulation, including strengthening main flux and blocking bypass metabolism, which significantly increased the β-alanine titer to 3.43 g/L. Next, the expression of ADC was optimized by promoter and untranslated region (UTR) engineering. Further transport engineering, which involved disrupting β-alanine importer CycA and heterologously expressing β-alanine exporter NCgI0580, improved β-alanine production to 8.48 g/L. Additionally, corn steep liquor was used to develop a cost-effective medium. The final strain produced 74.03 g/L β-alanine with a yield of 0.57 mol/mol glucose during fed-batch fermentation.

Conditionally unutilized proteins and their profound effects on growth and adaptation across microbial species

Protein synthesis is an important determinant of microbial growth and response that demands a high amount of metabolic and biosynthetic resources. Despite these costs, microbial species from different taxa and habitats massively synthesize proteins that are not utilized in the conditions they currently experience. Based on resource allocation models, recent studies have begun to reconcile the costs and benefits of these conditionally unutilized proteins (CUPs) in the context of varying environmental conditions. Such massive synthesis of CUPs is crucial to consider in different areas of modern microbiology, from the systematic investigation of cell physiology, via the prediction of evolution in laboratory and natural environments, to the rational design of strains in biotechnology applications.

Global Transcriptional Response of Escherichia coli Exposed In Situ to Different Low-Dose Ionizing Radiation Sources

Characterization of biological and chemical responses to ionizing radiation by various organisms is essential for potential applications in bioremediation, alternative modes of detecting nuclear material, and national security. Escherichia coli DH10β is an optimal system to study the microbial response to low-dose ionizing radiation at the transcriptional level because it is a well-characterized model bacterium and its responses to other environmental stressors, including those to higher radiation doses, have been elucidated in prior studies. In this study, RNA sequencing with downstream transcriptomic analysis (RNA-seq) was employed to characterize the global transcriptional response of stationary-phase E. coli subjected to ²³⁹Pu, ³H (tritium), and ⁵⁵Fe, at an approximate absorbed dose rate of 10 mGy day^-1 for 1 day and 15 days. Differential expression analysis identified significant changes in gene expression of E. coli for both short- and long-term exposures. Radionuclide source exposure induced differential expression in E. coli of genes involved in biosynthesis pathways of nuclear envelope components, amino acids, and siderophores, transport systems such as ABC transporters and type II secretion proteins, and initiation of stress response and regulatory systems of temperature stress, the RpoS regulon, and oxidative stress. These findings provide a basic understanding of the relationship between low-dose exposure and biological effect of a model bacterium that is critical for applications in alternative nuclear material detection and bioremediation. IMPORTANCE Escherichia coli strain DH10β, a well-characterized model bacterium, was subjected to short-term (1-day) and long-term (15-day) exposures to three different in situ radiation sources comprised of radionuclides relevant to nuclear activities to induce a measurable and identifiable genetic response. We found E. coli had both common and unique responses to the three exposures studied, suggesting both dose rate- and radionuclide-specific effects. This study is the first to provide insights into the transcriptional response of a microorganism in short- and long-term exposure to continuous low-dose ionizing radiation with multiple in situ radionuclide sources and the first to examine microbial transcriptional response in stationary phase. Moreover, this work provides a basis for the development of biosensors and informing more robust dose-response relationships to support ecological risk assessment.

Adaptation on xylose improves glucose-xylose co-utilization and ethanol production in a carbon catabolite repression (CCR) compromised ethanologenic strain

BACKGROUND

Sugar hydrolysates from lignocellulosic biomass are majorly composed of glucose and xylose that can be fermented to biofuels. Bacteria, despite having the natural ability to consume xylose are unable to consume it in presence of glucose due to a carbon catabolite repression (CCR) mechanism. This leads to overall reduced productivity as well as incomplete xylose utilization due to ethanol build-up from glucose utilization. In our effort to develop a strain for simultaneous fermentation of glucose and xylose into ethanol, we deleted ptsG in ethanologenic E. coli SSK42 to make it deficient in CCR and performed adaptive laboratory evolution to achieve accelerated growth rate, sugar consumption and ethanol production. Finally, we performed proteomics study to identify changes that might have been responsible for the observed improved phenotype of the evolved strain.

RESULTS

The parental strain of SSK42, i.e., wild-type E. coli B, did not co-utilize glucose and xylose as expected. After deleting the ptsG gene encoding the EIIBC^Glc subunit of PTS system, glucose consumption is severely affected in wild-type E. coli B. However, the ethanologenic, SSK42 strain, which was evolved in our earlier study on both glucose and xylose, didn't show such a drastic effect of EIIBC^Glc deletion, instead consumed glucose first, followed by xylose without delay for switching from one sugar to another. To improve growth on xylose and co-utilization capabilities, the ptsG deleted SSK42 was evolved on xylose. The strain evolved for 78 generations, strain SCD78, displayed significant co-utilization of glucose and xylose sugars. At the bioreactor level, the strain SCD78 produced 3-times the ethanol titer of the parent strain with significant glucose-xylose co-utilization. The rate of glucose and xylose consumption also increased 3.4-fold and 3-fold, respectively. Proteome data indicates significant upregulation of TCA cycle proteins, respiration-related proteins, and some transporters, which may have a role in increasing the total sugar consumption and co-utilization of sugars.

CONCLUSION

Through adaptive evolution, we have obtained a strain that has a significant glucose-xylose co-utilization phenotype with 3-fold higher total sugar consumption rate and ethanol production rate compared to the unevolved strain. This study also points out that adaptation on xylose is enough to impart glucose-xylose co-utilization property in CCR compromised ethanologenic strain SSK42.

C4-Dicarboxylates as Growth Substrates and Signaling Molecules for Commensal and Pathogenic Enteric Bacteria in Mammalian Intestine

The C4-dicarboxylates (C4-DC) l-aspartate and l-malate have been identified as playing an important role in the colonization of mammalian intestine by enteric bacteria, such as Escherichia coli and Salmonella enterica serovar Typhimurium, and succinate as a signaling molecule for host-enteric bacterium interaction. Thus, endogenous and exogenous fumarate respiration and related functions are required for efficient initial growth of the bacteria. l-Aspartate represents a major substrate for fumarate respiration in the intestine and a high-quality substrate for nitrogen assimilation. During nitrogen assimilation, DcuA catalyzes an l-aspartate/fumarate antiport and serves as a nitrogen shuttle for the net uptake of ammonium only, whereas DcuB acts as a redox shuttle that catalyzes the l-malate/succinate antiport during fumarate respiration. The C4-DC two-component system DcuS-DcuR is active in the intestine and responds to intestinal C4-DC levels. Moreover, in macrophages and in mice, succinate is a signal that promotes virulence and survival of S. Typhimurium and pathogenic E. coli. On the other hand, intestinal succinate is an important signaling molecule for the host and activates response and protective programs. Therefore, C4-DCs play a major role in supporting colonization of enteric bacteria and as signaling molecules for the adaptation of host physiology.

Suboptimal resource allocation in changing environments constrains response and growth in bacteria

To respond to fluctuating conditions, microbes typically need to synthesize novel proteins. As this synthesis relies on sufficient biosynthetic precursors, microbes must devise effective response strategies to manage depleting precursors. To better understand these strategies, we investigate the active response of Escherichia coli to changes in nutrient conditions, connecting transient gene expression to growth phenotypes. By synthetically modifying gene expression during changing conditions, we show how the competition by genes for the limited protein synthesis capacity constrains cellular response. Despite this constraint cells substantially express genes that are not required, trapping them in states where precursor levels are low and the genes needed to replenish the precursors are outcompeted. Contrary to common modeling assumptions, our findings highlight that cells do not optimize growth under changing environments but rather exhibit hardwired response strategies that may have evolved to promote fitness in their native environment. The constraint and the suboptimality of the cellular response uncovered provide a conceptual framework relevant for many research applications, from the prediction of evolution to the improvement of gene circuits in biotechnology.

Improved production of the non-native cofactor F420 in Escherichia coli

The deazaflavin cofactor F₄₂₀ is a low-potential, two-electron redox cofactor produced by some Archaea and Eubacteria that is involved in methanogenesis and methanotrophy, antibiotic biosynthesis, and xenobiotic metabolism. However, it is not produced by bacterial strains commonly used for industrial biocatalysis or recombinant protein production, such as Escherichia coli, limiting our ability to exploit it as an enzymatic cofactor and produce it in high yield. Here we have utilized a genome-scale metabolic model of E. coli and constraint-based metabolic modelling of cofactor F₄₂₀ biosynthesis to optimize F₄₂₀ production in E. coli. This analysis identified phospho-enol pyruvate (PEP) as a limiting precursor for F₄₂₀ biosynthesis, explaining carbon source-dependent differences in productivity. PEP availability was improved by using gluconeogenic carbon sources and overexpression of PEP synthase. By improving PEP availability, we were able to achieve a ~ 40-fold increase in the space–time yield of F₄₂₀ compared with the widely used recombinant Mycobacterium smegmatis expression system. This study establishes E. coli as an industrial F₄₂₀-production system and will allow the recombinant in vivo use of F₄₂₀-dependent enzymes for biocatalysis and protein engineering applications.

Identification of Enzymatic Bottlenecks for the Aerobic Production of Malate from Glycerol by the Systematic Gene Overexpression of Anaplerotic Enzymes in Escherichia coli

The biotechnological production of dicarboxylic acids (C4) from renewable carbon sources represents an attractive approach for the provision of these valuable compounds by green chemistry means. Glycerol has become a waste product of the biodiesel industry that serves as a highly reduced carbon source for some microorganisms. Escherichia coli is capable of consuming glycerol to produce succinate under anaerobic fermentation, but with the deletion of some tricarboxylic acid (TCA) cycle genes, it is also able to produce succinate and malate in aerobiosis. In this study, we investigate possible rate-limiting enzymes by overexpressing the C-feeding anaplerotic enzymes Ppc, MaeA, MaeB, and Pck in a mutant that lacks the succinate dehydrogenase (Sdh) enzyme. The overexpression of the TCA enzyme Mdh and the activation of the glyoxylate shunt was also examined. Using this unbiased approach, we found that phosphoenol pyruvate carboxylase (Ppc) overexpression enhances an oxidative pathway that leads to increasing succinate, while phosphoenol pyruvate carboxykinase (Pck) favors a more efficient reductive branch that produces mainly malate, at 57.5% of the theoretical maximum molar yield. The optimization of the culture medium revealed the importance of bicarbonate and pH in the production of malate. An additional mutation of the ppc gene highlights its central role in growth and C4 production.

Glycine Cleavage System and cAMP Receptor Protein Co-Regulate CRISPR/cas3 Expression to Resist Bacteriophage

The CRISPR/Cas system protects bacteria against bacteriophage and plasmids through a sophisticated mechanism where cas operon plays a crucial role consisting of cse1 and cas3. However, comprehensive studies on the regulation of cas3 operon of the Type I-E CRISPR/Cas system are scarce. Herein, we investigated the regulation of cas3 in Escherichia coli. The mutation in gcvP or crp reduced the CRISPR/Cas system interference ability and increased bacterial susceptibility to phage, when the casA operon of the CRISPR/Cas system was activated. The silence of the glycine cleavage system (GCS) encoded by gcvTHP operon reduced cas3 expression. Adding N⁵, N¹⁰-methylene tetrahydrofolate (N⁵, N¹⁰-mTHF), which is the product of GCS-catalyzed glycine, was able to activate cas3 expression. In addition, a cAMP receptor protein (CRP) encoded by crp activated cas3 expression via binding to the cas3 promoter in response to cAMP concentration. Since N⁵, N¹⁰-mTHF provides one-carbon unit for purine, we assumed GCS regulates cas3 through associating with CRP. It was evident that the mutation of gcvP failed to further reduce the cas3 expression with the crp deletion. These results illustrated a novel regulatory pathway which GCS and CRP co-regulate cas3 of the CRISPR/Cas system and contribute to the defence against invasive genetic elements, where CRP is indispensable for GCS regulation of cas3 expression.

Consistency of cybernetic variables with gene expression profiles: A more rigorous test

Diauxic growth of Escherichia coli is driven by a host of internal, complex regulatory actions. In this classic scenario of cellular control, the cell employs a rational algorithm to modulate its metabolism in a competitive fashion. Cybernetic models of metabolism, whose development now spans three decades, were first formulated to describe regulation of cells in complex, multi-substrate environments. They modeled this scenario using the hypothesis that the formation of the enzymatic machinery is regulated to maximize a return on investment. While this assumption is made on the basis of logical arguments rooted in evolutionary principles, little effort has been taken to validate if enzymes are truly synthesized in the same fashion that is predicted by cybernetic variables. This work revisits the original cybernetic models describing diauxic growth and compares their predictions of enzyme synthesis control with time series gene expression data in microarray and qRT-PCR formats. Three separate studies are made for two different strains of E. coli. The first is for the growth of E. coli BW25113 on a mixture of glucose and acetate, whose gene expression changes are metered by microarray. Another is also for the sequential consumption of glucose and acetate but involves strain MG1655 and employs qRT-PCR. The final is for E. coli MG1655 on glucose and lactose. By demonstrating how cybernetic variables for induced enzyme synthesis mimic the behavior of transcriptional data, a strong argument for using cybernetic models is made. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:858-867, 2018.

Windowed Granger causal inference strategy improves discovery of gene regulatory networks

Accurate inference of regulatory networks from experimental data facilitates the rapid characterization and understanding of biological systems. High-throughput technologies can provide a wealth of time-series data to better interrogate the complex regulatory dynamics inherent to organisms, but many network inference strategies do not effectively use temporal information. We address this limitation by introducing Sliding Window Inference for Network Generation (SWING), a generalized framework that incorporates multivariate Granger causality to infer network structure from time-series data. SWING moves beyond existing Granger methods by generating windowed models that simultaneously evaluate multiple upstream regulators at several potential time delays. We demonstrate that SWING elucidates network structure with greater accuracy in both in silico and experimentally validated in vitro systems. We estimate the apparent time delays present in each system and demonstrate that SWING infers time-delayed, gene-gene interactions that are distinct from baseline methods. By providing a temporal framework to infer the underlying directed network topology, SWING generates testable hypotheses for gene-gene influences.

Eliminating acetate formation improves citramalate production by metabolically engineered Escherichia coli

BACKGROUND

Citramalate, a chemical precursor to the industrially important methacrylic acid (MAA), can be synthesized using Escherichia coli overexpressing citramalate synthase (cimA gene). Deletion of gltA encoding citrate synthase and leuC encoding 3-isopropylmalate dehydratase were critical to achieving high citramalate yields. Acetate is an undesirable by-product potentially formed from pyruvate and acetyl-CoA, the precursors of citramalate during aerobic growth of E. coli. This study investigated strategies to minimize acetate and maximize citramalate production in E. coli mutants expressing the cimA gene.

RESULTS

Key knockouts that minimized acetate formation included acetate kinase (ackA), phosphotransacetylase (pta), and in particular pyruvate oxidase (poxB). Deletion of glucose 6-phosphate dehydrogenase (zwf) and ATP synthase (atpFH) aimed at improving glycolytic flux negatively impacted cell growth and citramalate accumulation in shake flasks. In a repetitive fed-batch process, E. coli gltA leuC ackA-pta poxB overexpressing cimA generated 54.1 g/L citramalate with a yield of 0.64 g/g glucose (78% of theoretical maximum yield), and only 1.4 g/L acetate in 87 h.

CONCLUSIONS

This study identified the gene deletions critical to reducing acetate accumulation during aerobic growth and citramalate production in metabolically engineered E. coli strains. The citramalate yield and final titer relative to acetate at the end of the fed-batch process are the highest reported to date (a mass ratio of citramalate to acetate of nearly 40) without being detrimental to citramalate productivity, significantly improving a potential process for the production of this five-carbon chemical.

Efficient anaerobic production of succinate from glycerol in engineered Escherichia coli by using dual carbon sources and limiting oxygen supply in preceding aerobic culture

Glycerol is an important resource for production of value-added bioproducts due to its large availability from the biodiesel industry as a by-product. In this study, two metabolic regulation strategies were applied in the aerobic stage of a two-stage fermentation to achieve high metabolic capacities of the pflB ldhA double mutant Escherichia coli strain overexpressing phosphoenolpyruvate carboxykinase (PCK) in the subsequent anaerobic stage: use of acetate as a co-carbon source of glycerol and restriction of oxygen supply in the PCK induction period. The succinate concentration achieved 926.7mM with a yield of 0.91mol/mol during the anaerobic stage of fermentation in a 1.5-L reactor. qRT-PCR indicated that the two strategies enhanced transcription of genes related with glycerol metabolism and succinate production. Our results showed this metabolically engineered E. coli strain has a great potential in producing succinate using glycerol as carbon source.

Loss of malic enzymes leads to metabolic imbalance and altered levels of trehalose and putrescine in the bacterium Sinorhizobium meliloti

BACKGROUND

Malic enzymes decarboxylate the tricarboxylic acid (TCA) cycle intermediate malate to the glycolytic end-product pyruvate and are well positioned to regulate metabolic flux in central carbon metabolism. Despite the wide distribution of these enzymes, their biological roles are unclear in part because the reaction catalyzed by these enzymes can be by-passed by other pathways. The N2-fixing alfalfa symbiont Sinorhizobium meliloti contains both a NAD(P)-malic enzyme (DME) and a separate NADP-malic enzyme (TME) and to help understand the role of these enzymes, we investigated growth, metabolomic, and transcriptional consequences resulting from loss of these enzymes in free-living cells.

RESULTS

Loss of DME, TME, or both enzymes had no effect on growth with the glycolytic substrate, glucose. In contrast, the dme mutants, but not tme, grew slowly on the gluconeogenic substrate succinate and this slow growth was further reduced upon the addition of glucose. The dme mutant strains incubated with succinate accumulated trehalose and hexose sugar phosphates, secreted malate, and relative to wild-type, these cells had moderately increased transcription of genes involved in gluconeogenesis and pathways that divert metabolites away from the TCA cycle. While tme mutant cells grew at the same rate as wild-type on succinate, they accumulated the compatible solute putrescine.

CONCLUSIONS

NAD(P)-malic enzyme (DME) of S. meliloti is required for efficient metabolism of succinate via the TCA cycle. In dme mutants utilizing succinate, malate accumulates and is excreted and these cells appear to increase metabolite flow via gluconeogenesis with a resulting increase in the levels of hexose-6-phosphates and trehalose. For cells utilizing succinate, TME activity alone appeared to be insufficient to produce the levels of pyruvate required for efficient TCA cycle metabolism. Putrescine was found to accumulate in tme cells growing with succinate, and whether this is related to altered levels of NADPH requires further investigation.

Metabolic design of a platform Escherichia coli strain producing various chorismate derivatives

A synthetic metabolic pathway suitable for the production of chorismate derivatives was designed in Escherichia coli. An L-phenylalanine-overproducing E. coli strain was engineered to enhance the availability of phosphoenolpyruvate (PEP), which is a key precursor in the biosynthesis of aromatic compounds in microbes. Two major reactions converting PEP to pyruvate were inactivated. Using this modified E.coli as a base strain, we tested our system by carrying out the production of salicylate, a high-demand aromatic chemical. The titer of salicylate reached 11.5 g/L in batch culture after 48 h cultivation in a 2-liter jar fermentor, and the yield from glucose as the sole carbon source exceeded 40% (mol/mol). In this test case, we found that pyruvate was synthesized primarily via salicylate formation and the reaction converting oxaloacetate to pyruvate. In order to demonstrate the generality of our designed strain, we employed this platform for the production of each of 7 different chorismate derivatives. Each of these industrially important chemicals was successfully produced to levels of 1-3g/L in test tube-scale culture.

Robust reconstruction of gene expression profiles from reporter gene data using linear inversion

MOTIVATION

Time-series observations from reporter gene experiments are commonly used for inferring and analyzing dynamical models of regulatory networks. The robust estimation of promoter activities and protein concentrations from primary data is a difficult problem due to measurement noise and the indirect relation between the measurements and quantities of biological interest.

RESULTS

We propose a general approach based on regularized linear inversion to solve a range of estimation problems in the analysis of reporter gene data, notably the inference of growth rate, promoter activity, and protein concentration profiles. We evaluate the validity of the approach using in silico simulation studies, and observe that the methods are more robust and less biased than indirect approaches usually encountered in the experimental literature based on smoothing and subsequent processing of the primary data. We apply the methods to the analysis of fluorescent reporter gene data acquired in kinetic experiments with Escherichia coli. The methods are capable of reliably reconstructing time-course profiles of growth rate, promoter activity and protein concentration from weak and noisy signals at low population volumes. Moreover, they capture critical features of those profiles, notably rapid changes in gene expression during growth transitions.

AVAILABILITY AND IMPLEMENTATION

The methods described in this article are made available as a Python package (LGPL license) and also accessible through a web interface. For more information, see https://team.inria.fr/ibis/wellinverter.

NADPH-generating systems in bacteria and archaea

Reduced nicotinamide adenine dinucleotide phosphate (NADPH) is an essential electron donor in all organisms. It provides the reducing power that drives numerous anabolic reactions, including those responsible for the biosynthesis of all major cell components and many products in biotechnology. The efficient synthesis of many of these products, however, is limited by the rate of NADPH regeneration. Hence, a thorough understanding of the reactions involved in the generation of NADPH is required to increase its turnover through rational strain improvement. Traditionally, the main engineering targets for increasing NADPH availability have included the dehydrogenase reactions of the oxidative pentose phosphate pathway and the isocitrate dehydrogenase step of the tricarboxylic acid (TCA) cycle. However, the importance of alternative NADPH-generating reactions has recently become evident. In the current review, the major canonical and non-canonical reactions involved in the production and regeneration of NADPH in prokaryotes are described, and their key enzymes are discussed. In addition, an overview of how different enzymes have been applied to increase NADPH availability and thereby enhance productivity is provided.

Acetate Exposure Determines the Diauxic Behavior of Escherichia coli during the Glucose-Acetate Transition

Growth of Escherichia coli on glucose in batch culture is accompanied by the excretion of acetate, which is consumed by the cells when glucose is exhausted. This glucose-acetate transition is classically described as a diauxie (two successive growth stages). Here, we investigated the physiological and metabolic properties of cells after glucose exhaustion through the analysis of growth parameters and gene expression. We found that E. coli cells grown on glucose in batch culture produce acetate and consume it after glucose exhaustion but do not grow on acetate. Acetate is catabolized, but key anabolic genes--such as the genes encoding enzymes of the glyoxylate shunt--are not upregulated, hence preventing growth. Both the induction of the latter anabolic genes and growth were observed only after prolonged exposure to low concentrations of acetate and could be accelerated by high acetate concentrations. We postulate that such decoupling between acetate catabolism and acetate anabolism might be an advantage for the survival of E. coli in the ever-changing environment of the intestine.

IMPORTANCE

The glucose-acetate transition is a valuable experimental model for comprehensive investigations of metabolic adaptation and a current paradigm for developing modeling approaches in systems microbiology. Yet, the work reported in our paper demonstrates that the metabolic behavior of Escherichia coli during the glucose-acetate transition is much more complex than what has been reported so far. A decoupling between acetate catabolism and acetate anabolism was observed after glucose exhaustion, which has not been reported previously. This phenomenon could represent a strategy for optimal utilization of carbon resources during colonization and persistence of E. coli in the gut and is also of significant interest for biotechnological applications.

Structural basis for methyl-donor-dependent and sequence-specific binding to tRNA substrates by knotted methyltransferase TrmD

The deep trefoil knot architecture is unique to the SpoU and tRNA methyltransferase D (TrmD) (SPOUT) family of methyltransferases (MTases) in all three domains of life. In bacteria, TrmD catalyzes the N(1)-methylguanosine (m(1)G) modification at position 37 in transfer RNAs (tRNAs) with the (36)GG(37) sequence, using S-adenosyl-l-methionine (AdoMet) as the methyl donor. The m(1)G37-modified tRNA functions properly to prevent +1 frameshift errors on the ribosome. Here we report the crystal structure of the TrmD homodimer in complex with a substrate tRNA and an AdoMet analog. Our structural analysis revealed the mechanism by which TrmD binds the substrate tRNA in an AdoMet-dependent manner. The trefoil-knot center, which is structurally conserved among SPOUT MTases, accommodates the adenosine moiety of AdoMet by loosening/retightening of the knot. The TrmD-specific regions surrounding the trefoil knot recognize the methionine moiety of AdoMet, and thereby establish the entire TrmD structure for global interactions with tRNA and sequential and specific accommodations of G37 and G36, resulting in the synthesis of m(1)G37-tRNA.

Salmonella typhimurium and Escherichia coli dissimilarity: Closely related bacteria with distinct metabolic profiles

Live attenuated strains of Salmonella typhimurium have been extensively investigated as vaccines for a number of infectious diseases. However, there is still little information available concerning aspects of their metabolism. S. typhimurium and Escherichia coli show a high degree of similarity in terms of their genome contents and metabolic networks. However, this work presents experimental evidence showing that significant differences exist in their abilities to direct carbon fluxes to biomass and energy production. It is important to study the metabolism of Salmonella to elucidate the formation of acetate and other metabolites involved in optimizing the production of biomass, essential for the development of recombinant vaccines. The metabolism of Salmonella under aerobic conditions was assessed using continuous cultures performed at dilution rates ranging from 0.1 to 0.67 h(-1), with glucose as main substrate. Acetate assimilation and glucose metabolism under anaerobic conditions were also investigated using batch cultures. Chemostat cultivations showed deviation of carbon towards acetate formation, starting at dilution rates above 0.1 h(-1). This differed from previous findings for E. coli, where acetate accumulation was only detected at dilution rates exceeding 0.4 h(-1), and was due to the lower rate of acetate assimilation by S. typhimurium under aerobic conditions. Under anaerobic conditions, both microorganisms mainly produced ethanol, acetate, and formate. A genome-scale metabolic model, reconstructed for Salmonella based on an E. coli model, provided a poor description of the mixed fermentation pattern observed during Salmonella cultures, reinforcing the different patterns of carbon utilization exhibited by these closely related bacteria.

Phenotypic bistability in Escherichia coli's central carbon metabolism

Fluctuations in intracellular molecule abundance can lead to distinct, coexisting phenotypes in isogenic populations. Although metabolism continuously adapts to unpredictable environmental changes, and although bistability was found in certain substrate-uptake pathways, central carbon metabolism is thought to operate deterministically. Here, we combine experiment and theory to demonstrate that a clonal Escherichia coli population splits into two stochastically generated phenotypic subpopulations after glucose-gluconeogenic substrate shifts. Most cells refrain from growth, entering a dormant persister state that manifests as a lag phase in the population growth curve. The subpopulation-generating mechanism resides at the metabolic core, overarches the metabolic and transcriptional networks, and only allows the growth of cells initially achieving sufficiently high gluconeogenic flux. Thus, central metabolism does not ensure the gluconeogenic growth of individual cells, but uses a population-level adaptation resulting in responsive diversification upon nutrient changes.

Metabolic evolution of two reducing equivalent-conserving pathways for high-yield succinate production in Escherichia coli

Reducing equivalents are an important cofactor for efficient synthesis of target products. During metabolic evolution to improve succinate production in Escherichia coli strains, two reducing equivalent-conserving pathways were activated to increase succinate yield. The sensitivity of pyruvate dehydrogenase to NADH inhibition was eliminated by three nucleotide mutations in the lpdA gene. Pyruvate dehydrogenase activity increased under anaerobic conditions, which provided additional NADH. The pentose phosphate pathway and transhydrogenase were activated by increased activities of transketolase and soluble transhydrogenase SthA. These data suggest that more carbon flux went through the pentose phosphate pathway, thus leading to production of more reducing equivalent in the form of NADPH, which was then converted to NADH through soluble transhydrogenase for succinate production. Reverse metabolic engineering was further performed in a parent strain, which was not metabolically evolved, to verify the effects of activating these two reducing equivalent-conserving pathways for improving succinate yield. Activating pyruvate dehydrogenase increased succinate yield from 1.12 to 1.31mol/mol, whereas activating the pentose phosphate pathway and transhydrogenase increased succinate yield from 1.12 to 1.33mol/mol. Activating these two pathways in combination led to a succinate yield of 1.5mol/mol (88% of theoretical maximum), suggesting that they exhibited a synergistic effect for improving succinate yield.

Network component analysis provides quantitative insights on an Arabidopsis transcription factor-gene regulatory network

Background

Gene regulatory networks (GRNs) are models of molecule-gene interactions instrumental in the coordination of gene expression. Transcription factor (TF)-GRNs are an important subset of GRNs that characterize gene expression as the effect of TFs acting on their target genes. Although such networks can qualitatively summarize TF-gene interactions, it is highly desirable to quantitatively determine the strengths of the interactions in a TF-GRN as well as the magnitudes of TF activities. To our knowledge, such analysis is rare in plant biology. A computational methodology developed for this purpose is network component analysis (NCA), which has been used for studying large-scale microbial TF-GRNs to obtain nontrivial, mechanistic insights. In this work, we employed NCA to quantitatively analyze a plant TF-GRN important in floral development using available regulatory information from AGRIS, by processing previously reported gene expression data from four shoot apical meristem cell types.

Results

The NCA model satisfactorily accounted for gene expression measurements in a TF-GRN of seven TFs (LFY, AG, SEPALLATA3 [SEP3], AP2, AGL15, HY5 and AP3/PI) and 55 genes. NCA found strong interactions between certain TF-gene pairs including LFY → MYB17, AG → CRC, AP2 → RD20, AGL15 → RAV2 and HY5 → HLH1, and the direction of the interaction (activation or repression) for some AGL15 targets for which this information was not previously available. The activity trends of four TFs - LFY, AG, HY5 and AP3/PI as deduced by NCA correlated well with the changes in expression levels of the genes encoding these TFs across all four cell types; such a correlation was not observed for SEP3, AP2 and AGL15.

Conclusions

For the first time, we have reported the use of NCA to quantitatively analyze a plant TF-GRN important in floral development for obtaining nontrivial information about connectivity strengths between TFs and their target genes as well as TF activity. However, since NCA relies on documented connectivity information about the underlying TF-GRN, it is currently limited in its application to larger plant networks because of the lack of documented connectivities. In the future, the identification of interactions between plant TFs and their target genes on a genome scale would allow the use of NCA to provide quantitative regulatory information about plant TF-GRNs, leading to improved insights on cellular regulatory programs.

Physiological and Molecular Timing of the Glucose to Acetate Transition in Escherichia coli

The glucose-acetate transition in Escherichia coli is a classical model of metabolic adaptation. Here, we describe the dynamics of the molecular processes involved in this metabolic transition, with a particular focus on glucose exhaustion. Although changes in the metabolome were observed before glucose exhaustion, our results point to a massive reshuffling at both the transcriptome and metabolome levels in the very first min following glucose exhaustion. A new transcriptional pattern, involving a change in genome expression in one-sixth of the E. coli genome, was established within 10 min and remained stable until the acetate was completely consumed. Changes in the metabolome took longer and stabilized 40 min after glucose exhaustion. Integration of multi-omics data revealed different modifications and timescales between the transcriptome and metabolome, but both point to a rapid adaptation of less than an hour. This work provides detailed information on the order, timing and extent of the molecular and physiological events that occur during the glucose-acetate transition and that are of particular interest for the development of dynamic models of metabolism.

Activating phosphoenolpyruvate carboxylase and phosphoenolpyruvate carboxykinase in combination for improvement of succinate production

Phosphoenolpyruvate (PEP) carboxylation is an important step in the production of succinate by Escherichia coli. Two enzymes, PEP carboxylase (PPC) and PEP carboxykinase (PCK), are responsible for PEP carboxylation. PPC has high substrate affinity and catalytic velocity but wastes the high energy of PEP. PCK has low substrate affinity and catalytic velocity but can conserve the high energy of PEP for ATP formation. In this work, the expression of both the ppc and pck genes was modulated, with multiple regulatory parts of different strengths, in order to investigate the relationship between PPC or PCK activity and succinate production. There was a positive correlation between PCK activity and succinate production. In contrast, there was a positive correlation between PPC activity and succinate production only when PPC activity was within a certain range; excessive PPC activity decreased the rates of both cell growth and succinate formation. These two enzymes were also activated in combination in order to recruit the advantages of each for the improvement of succinate production. It was demonstrated that PPC and PCK had a synergistic effect in improving succinate production.

Importance of understanding the main metabolic regulation in response to the specific pathway mutation for metabolic engineering of Escherichia coli

Recent metabolic engineering practice was briefly reviewed in particular for the useful metabolite production such as natural products and biofuel productions. With the emphasis on systems biology approach, the metabolic regulation of the main metabolic pathways in E. coli was discussed from the points of view of enzyme level (allosteric and phosphorylation/ dephosphorylation) regulation, and gene level (transcriptional) regulation. Then the effects of the specific pathway gene knockout such as pts, pgi, zwf, gnd, pyk, ppc, pckA, lpdA, pfl gene knockout on the metabolism in E. coli were overviewed from the systems biology point of view with possible application for strain improvement point.

Visualizing evolution in real time to determine the molecular mechanisms of n-butanol tolerance in Escherichia coli

Toxicity of products or feedstock components poses a challenge in the biocatalyst-based production of fuels and chemicals. The genetic determinants that are involved in increased resistance to an inhibitor form the adaptive landscape for the phenotype; so in order to engineer more robust biocatalysts, a better understanding of the adaptive landscape is required. Here, we used an adaptive laboratory evolution method called visualizing evolution in real time (VERT) to help map out part of the adaptive landscape of Escherichia coli tolerance to the biofuel n-butanol. VERT enables identification of adaptive events (population expansions triggered by adaptive mutants) via visualization of the relative proportions of different fluorescently-labeled cells. Knowledge of the occurrence of adaptive events allows for a more systematic isolation of adaptive mutants while simultaneously reducing the number of missed adaptive mutants (and the underlying adaptive mechanisms) that result from clonal interference during the course of in vitro evolution. Based on the evolutionary dynamics observed, clonal interference was found to play a significant role in shaping the population structure of E. coli during exposure to n-butanol, and VERT helped to facilitate the isolation of adaptive mutants from the population. We further combined adaptive laboratory evolution with genome shuffling to significantly enhance the desired n-butanol tolerance phenotype. Subsequent transcriptome analysis of the isolated adaptive mutants revealed different mechanisms of n-butanol resistance in different lineages. In one fluorescently-marked subpopulation, members of the Fur regulon were upregulated; which was not observed in the other subpopulation. In addition, genome sequencing of several adaptive mutants revealed the genetic basis for some of the observed transcriptome profiles. We further elucidated the potential role of the iron-related gene in n-butanol tolerance via overexpression and deletion studies and hypothesized that the upregulation of the iron-related genes indirectly led to modifications in the outer membrane, which contributed to enhanced n-butanol tolerance.

Determining PTEN functional status by network component deduced transcription factor activities

PTEN-controlled PI3K-AKT-mTOR pathway represents one of the most deregulated signaling pathways in human cancers. With many small molecule inhibitors that target PI3K-AKT-mTOR pathway being exploited clinically, sensitive and reliable ways of stratifying patients according to their PTEN functional status and determining treatment outcomes are urgently needed. Heterogeneous loss of PTEN is commonly associated with human cancers and yet PTEN can also be regulated on epigenetic, transcriptional or post-translational levels, which makes the use of simple protein or gene expression-based analyses in determining PTEN status less accurate. In this study, we used network component analysis to identify 20 transcription factors (TFs) whose activities deduced from their target gene expressions were immediately altered upon the re-expression of PTEN in a PTEN-inducible system. Interestingly, PTEN controls the activities (TFA) rather than the expression levels of majority of these TFs and these PTEN-controlled TFAs are substantially altered in prostate cancer mouse models. Importantly, the activities of these TFs can be used to predict PTEN status in human prostate, breast and brain tumor samples with enhanced reliability when compared to straightforward IHC-based or expression-based analysis. Furthermore, our analysis indicates that unique sets of PTEN-controlled TFAs significantly contribute to specific tumor types. Together, our findings reveal that TFAs may be used as “signatures” for predicting PTEN functional status and elucidate the transcriptional architectures underlying human cancers caused by PTEN loss.

Transcriptional regulation of sdiA by cAMP-receptor protein, LeuO, and environmental signals in Salmonella enterica serovar Typhimurium

The sdiA gene encodes for a LuxR-type transcription factor, which is active when bound to N-acyl homoserine lactones (AHLs). Because Salmonella enterica serovar Typhimurium does not produce AHLs, SdiA senses signals produced by other organisms. SdiA is not expressed constitutively, and response is limited to conditions in which elevated expression occurs, but little is known about the regulation of sdiA expression. Here we map the sdiA promoter and define several regulators that directly or indirectly act on the promoter. The major activator of sdiA expression is cAMP-receptor protein (CRP), and we define the CRP operator in the sdiA promoter using promoter and crp mutants. LeuO activates sdiA expression to a lesser extent than does CRP. We demonstrate that LeuO directly binds the sdiA promoter and the Rcs phosphorelay represses sdiA expression. In this study, NhaR, IlvY, and Fur affected sdiA expression indirectly and weakly. Expression in late-stationary phase depended on RpoS. AHL-dependent expression of the SdiA-regulated gene rck correlated to the observed sdiA transcriptional changes in regulator mutants. The data demonstrate that regulation of sdiA involves integration of multiple environmental and metabolic signals.

Importance of metabolic coupling for the dynamics of gene expression following a diauxic shift in Escherichia coli

Gene regulatory networks consist of direct interactions, but also include indirect interactions mediated by metabolism. We investigate to which extent these indirect interactions arising from metabolic coupling influence the dynamics of the system. To this end, we build a qualitative model of the gene regulatory network controlling carbon assimilation in Escherichia coli, and use this model to study the changes in gene expression following a diauxic shift from glucose to acetate. In particular, we compare the relative variation in the steady-state concentrations of enzymes and transcription regulators during growth on glucose and acetate, as well as the dynamic response of gene expression to the exhaustion of glucose and the subsequent assimilation of acetate. We find significant differences between the dynamics of the system in the absence and presence of metabolic coupling. This shows that interactions arising from metabolic coupling cannot be ignored when studying the dynamics of gene regulatory networks.

L-malate production by metabolically engineered Escherichia coli

Escherichia coli strains (KJ060 and KJ073) that were previously developed for succinate production have now been modified for malate production. Many unexpected changes were observed during this investigation. The initial strategy of deleting fumarase isoenzymes was ineffective, and succinate continued to accumulate. Surprisingly, a mutation in fumarate reductase alone was sufficient to redirect carbon flow into malate even in the presence of fumarase. Further deletions were needed to inactivate malic enzymes (typically gluconeogenic) and prevent conversion to pyruvate. However, deletion of these genes (sfcA and maeB) resulted in the unexpected accumulation of D-lactate despite the prior deletion of mgsA and ldhA and the absence of apparent lactate dehydrogenase activity. Although the metabolic source of this D-lactate was not identified, lactate accumulation was increased by supplementation with pyruvate and decreased by the deletion of either pyruvate kinase gene (pykA or pykF) to reduce the supply of pyruvate. Many of the gene deletions adversely affected growth and cell yield in minimal medium under anaerobic conditions, and volumetric rates of malate production remained low. The final strain (XZ658) produced 163 mM malate, with a yield of 1.0 mol (mol glucose(-1)), half of the theoretical maximum. Using a two-stage process (aerobic cell growth and anaerobic malate production), this engineered strain produced 253 mM malate (34 g liter(-1)) within 72 h, with a higher yield (1.42 mol mol(-1)) and productivity (0.47 g liter(-1) h(-1)). This malate yield and productivity are equal to or better than those of other known biocatalysts.

The carbon assimilation network in Escherichia coli is densely connected and largely sign-determined by directions of metabolic fluxes

Gene regulatory networks consist of direct interactions but also include indirect interactions mediated by metabolites and signaling molecules. We describe how these indirect interactions can be derived from a model of the underlying biochemical reaction network, using weak time-scale assumptions in combination with sensitivity criteria from metabolic control analysis. We apply this approach to a model of the carbon assimilation network in Escherichia coli. Our results show that the derived gene regulatory network is densely connected, contrary to what is usually assumed. Moreover, the network is largely sign-determined, meaning that the signs of the indirect interactions are fixed by the flux directions of biochemical reactions, independently of specific parameter values and rate laws. An inversion of the fluxes following a change in growth conditions may affect the signs of the indirect interactions though. This leads to a feedback structure that is at the same time robust to changes in the kinetic properties of enzymes and that has the flexibility to accommodate radical changes in the environment.

The regulation of gene expression is tightly interwoven with metabolism and signal transduction. A realistic view of gene regulatory networks should therefore not only include direct interactions resulting from transcription regulation, but also indirect regulatory interactions mediated by metabolic effectors and signaling molecules. Ignoring these indirect interactions during the analysis of the network dynamics may lead crucial feedback loops to be missed. We present a method for systematically deriving indirect interactions from a model of the underlying biochemical reaction network, using weak time-scale assumptions in combination with sensitivity criteria from metabolic control analysis. This approach leads to novel insights as exemplified here on the carbon assimilation network of E. coli. We show that the derived gene regulatory network is densely connected, that the signs of the indirect interactions are largely fixed by the direction of metabolic fluxes, and that a change in flux direction may invert the sign of indirect interactions. Therefore the feedback structure of the network is much more complex than usually assumed; it appears robust to changes in the kinetic properties of its components and it can be flexibly rewired when the environment changes.

Metabolic engineering for production of biorenewable fuels and chemicals: contributions of synthetic biology

Production of fuels and chemicals through microbial fermentation of plant material is a desirable alternative to petrochemical-based production. Fermentative production of biorenewable fuels and chemicals requires the engineering of biocatalysts that can quickly and efficiently convert sugars to target products at a cost that is competitive with existing petrochemical-based processes. It is also important that biocatalysts be robust to extreme fermentation conditions, biomass-derived inhibitors, and their target products. Traditional metabolic engineering has made great advances in this area, but synthetic biology has contributed and will continue to contribute to this field, particularly with next-generation biofuels. This work reviews the use of metabolic engineering and synthetic biology in biocatalyst engineering for biorenewable fuels and chemicals production, such as ethanol, butanol, acetate, lactate, succinate, alanine, and xylitol. We also examine the existing challenges in this area and discuss strategies for improving biocatalyst tolerance to chemical inhibitors.

Transcriptional patterns in both host and bacterium underlie a daily rhythm of anatomical and metabolic change in a beneficial symbiosis

Mechanisms for controlling symbiont populations are critical for maintaining the associations that exist between a host and its microbial partners. We describe here the transcriptional, metabolic, and ultrastructural characteristics of a diel rhythm that occurs in the symbiosis between the squid Euprymna scolopes and the luminous bacterium Vibrio fischeri. The rhythm is driven by the host's expulsion from its light-emitting organ of most of the symbiont population each day at dawn. The transcriptomes of both the host epithelium that supports the symbionts and the symbiont population itself were characterized and compared at four times over this daily cycle. The greatest fluctuation in gene expression of both partners occurred as the day began. Most notable was an up-regulation in the host of >50 cytoskeleton-related genes just before dawn and their subsequent down-regulation within 6 h. Examination of the epithelium by TEM revealed a corresponding restructuring, characterized by effacement and blebbing of its apical surface. After the dawn expulsion, the epithelium reestablished its polarity, and the residual symbionts began growing, repopulating the light organ. Analysis of the symbiont transcriptome suggested that the bacteria respond to the effacement by up-regulating genes associated with anaerobic respiration of glycerol; supporting this finding, lipid analysis of the symbionts' membranes indicated a direct incorporation of host-derived fatty acids. After 12 h, the metabolic signature of the symbiont population shifted to one characteristic of chitin fermentation, which continued until the following dawn. Thus, the persistent maintenance of the squid-vibrio symbiosis is tied to a dynamic diel rhythm that involves both partners.

Metabolic evolution of energy-conserving pathways for succinate production in Escherichia coli

During metabolic evolution to improve succinate production in Escherichia coli strains, significant changes in cellular metabolism were acquired that increased energy efficiency in two respects. The energy-conserving phosphoenolpyruvate (PEP) carboxykinase (pck), which normally functions in the reverse direction (gluconeogenesis; glucose repressed) during the oxidative metabolism of organic acids, evolved to become the major carboxylation pathway for succinate production. Both PCK enzyme activity and gene expression levels increased significantly in two stages because of several mutations during the metabolic evolution process. High-level expression of this enzyme-dominated CO(2) fixation and increased ATP yield (1 ATP per oxaloacetate). In addition, the native PEP-dependent phosphotransferase system for glucose uptake was inactivated by a mutation in ptsI. This glucose transport function was replaced by increased expression of the GalP permease (galP) and glucokinase (glk). Results of deleting individual transport genes confirmed that GalP served as the dominant glucose transporter in evolved strains. Using this alternative transport system would increase the pool of PEP available for redox balance. This change would also increase energy efficiency by eliminating the need to produce additional PEP from pyruvate, a reaction that requires two ATP equivalents. Together, these changes converted the wild-type E. coli fermentation pathway for succinate into a functional equivalent of the native pathway that nature evolved in succinate-producing rumen bacteria.

Reengineering Escherichia coli for Succinate Production in Mineral Salts Medium

The fermentative metabolism of glucose was redirected to succinate as the primary product without mutating any genes encoding the native mixed-acid fermentation pathway or redox reactions. Two changes in peripheral pathways were together found to increase succinate yield fivefold: (i) increased expression of phosphoenolpyruvate carboxykinase and (ii) inactivation of the glucose phosphoenolpyruvate-dependent phosphotransferase system. These two changes increased net ATP production, increased the pool of phosphoenolpyruvate available for carboxylation, and increased succinate production. Modest further improvements in succinate yield were made by inactivating the pflB gene, encoding pyruvate formate lyase, resulting in an Escherichia coli pathway that is functionally similar to the native pathway in Actinobacillus succinogenes and other succinate-producing rumen bacteria.

Glucose and glycolysis are required for the successful infection of macrophages and mice by Salmonella enterica serovar typhimurium

Salmonella is a widespread zoonotic enteropathogen that causes gastroenteritis and fatal typhoidal disease in mammals. During systemic infection of mice, Salmonella enterica serovar Typhimurium resides and replicates in macrophages within the "Salmonella-containing vacuole" (SCV). It is surprising that the substrates and metabolic pathways necessary for growth of S. Typhimurium within the SCV of macrophages have not been identified yet. To determine whether S. Typhimurium utilized sugars within the SCV, we constructed a series of S. Typhimurium mutants that lacked genes involved in sugar transport and catabolism and tested them for replication in mice and macrophages. These mutants included a mutant with a mutation in the pfkAB-encoded phosphofructokinase, which catalyzes a key committing step in glycolysis. We discovered that a pfkAB mutant is severely attenuated for replication and survival within RAW 264.7 macrophages. We also show that disruption of the phosphoenolpyruvate:carbohydrate phosphotransferase system by deletion of the ptsHI and crr genes reduces S. Typhimurium replication within RAW 264.7 macrophages. We discovered that mutants unable to catabolize glucose due to deletion of ptsHI, crr, and glk or deletion of ptsG, manXYZ, and glk showed reduced replication within RAW 264.7 macrophages. This study proves that S. Typhimurium requires glycolysis for infection of mice and macrophages and that transport of glucose is required for replication within macrophages.

Fast network component analysis (FastNCA) for gene regulatory network reconstruction from microarray data

MOTIVATION

Recently developed network component analysis (NCA) approach is promising for gene regulatory network reconstruction from microarray data. The existing NCA algorithm is an iterative method which has two potential limitations: computational instability and multiple local solutions. The subsequently developed NCA-r algorithm with Tikhonov regularization can help solve the first issue but cannot completely handle the second one. Here we develop a novel Fast Network Component Analysis (FastNCA) algorithm which has an analytical solution that is much faster and does not have the above limitations.

RESULTS

Firstly FastNCA is compared to NCA and NCA-r using synthetic data. The reconstruction of FastNCA is more accurate than that of NCA-r and comparable to that of properly converged NCA. FastNCA is not sensitive to the correlation among the input signals, while its performance does degrade a little but not as dramatically as that of NCA. Like NCA, FastNCA is not very sensitive to small inaccuracies in a priori information on the network topology. FastNCA is about several tens times faster than NCA and several hundreds times faster than NCA-r. Then, the method is applied to real yeast cell-cycle microarray data. The activities of the estimated cell-cycle regulators by FastNCA and NCA-r are compared to the semi-quantitative results obtained independently by Lee et al. (2002). It is shown here that there is a greater agreement between the results of FastNCA and Lee's, which is represented by the ratio 23/33, than that between the results of NCA-r and Lee's, which is 14/33.

AVAILABILITY

Software and supplementary materials are available from http://www.eee.hku.hk/~cqchang/FastNCA.htm