Using genomic sequencing for classical genetics in E. coli K12

Cell Geometry and Membrane Protein Crowding Constrain Growth Rate, Overflow Metabolism, Respiration, and Maintenance Energy

A metabolic theory is presented for predicting maximum growth rate, overflow metabolism, respiration efficiency, and maintenance energy flux based on the intersection of cell geometry, membrane protein crowding, and metabolism. The importance of cytosolic macromolecular crowding on phenotype has been established in the literature but the importance of surface area has been largely overlooked due to incomplete knowledge of membrane properties. We demonstrate that the capacity of the membrane to host proteins increases with growth rate offsetting decreases in surface area-to-volume ratios (SA:V). This increase in membrane protein is hypothesized to be essential to competitive Escherichia coli phenotypes. The presented membrane-centric theory uses biophysical properties and metabolic systems analysis to successfully predict the phenotypes of E. coli K-12 strains, MG1655 and NCM3722, which are genetically similar but have SA:V ratios that differ up to 30%, maximum growth rates on glucose media that differ by 40%, and overflow phenotypes that start at growth rates that differ by 80%. These analyses did not consider cytosolic macromolecular crowding, highlighting the distinct properties of the presented theory. Cell geometry and membrane protein crowding are significant biophysical constraints on phenotype and provide a theoretical framework for improved understanding and control of cell biology.

Direct measurement of dynamic attractant gradients reveals breakdown of the Patlak-Keller-Segel chemotaxis model

Chemotactic bacteria not only navigate chemical gradients, but also shape their environments by consuming and secreting attractants. Investigating how these processes influence the dynamics of bacterial populations has been challenging because of a lack of experimental methods for measuring spatial profiles of chemoattractants in real time. Here, we use a fluorescent sensor for aspartate to directly measure bacterially generated chemoattractant gradients during collective migration. Our measurements show that the standard Patlak-Keller-Segel model for collective chemotactic bacterial migration breaks down at high cell densities. To address this, we propose modifications to the model that consider the impact of cell density on bacterial chemotaxis and attractant consumption. With these changes, the model explains our experimental data across all cell densities, offering insight into chemotactic dynamics. Our findings highlight the significance of considering cell density effects on bacterial behavior, and the potential for fluorescent metabolite sensors to shed light on the complex emergent dynamics of bacterial communities.

High quality genomes produced from single MinION flow cells clarify polyploid and demographic histories of critically endangered Fraxinus (ash) species

With populations of threatened and endangered species declining worldwide, efforts are being made to generate high quality genomic records of these species before they are lost forever. Here, we demonstrate that data from single Oxford Nanopore Technologies (ONT) MinION flow cells can, even in the absence of highly accurate short DNA-read polishing, produce high quality de novo plant genome assemblies adequate for downstream analyses, such as synteny and ploidy evaluations, paleodemographic analyses, and phylogenomics. This study focuses on three North American ash tree species in the genus Fraxinus (Oleaceae) that were recently added to the International Union for Conservation of Nature (IUCN) Red List as critically endangered. Our results support a hexaploidy event at the base of the Oleaceae as well as a subsequent whole genome duplication shared by Syringa, Osmanthus, Olea, and Fraxinus. Finally, we demonstrate the use of ONT long-read sequencing data to reveal patterns in demographic history.

Functional decomposition of metabolism allows a system-level quantification of fluxes and protein allocation towards specific metabolic functions

Quantifying the contribution of individual molecular components to complex cellular processes is a grand challenge in systems biology. Here we establish a general theoretical framework (Functional Decomposition of Metabolism, FDM) to quantify the contribution of every metabolic reaction to metabolic functions, e.g. the synthesis of biomass building blocks. FDM allowed for a detailed quantification of the energy and biosynthesis budget for growing Escherichia coli cells. Surprisingly, the ATP generated during the biosynthesis of building blocks from glucose almost balances the demand from protein synthesis, the largest energy expenditure known for growing cells. This leaves the bulk of the energy generated by fermentation and respiration unaccounted for, thus challenging the common notion that energy is a key growth-limiting resource. Moreover, FDM together with proteomics enables the quantification of enzymes contributing towards each metabolic function, allowing for a first-principle formulation of a coarse-grained model of global protein allocation based on the structure of the metabolic network.

Direct measurement of dynamic attractant gradients reveals breakdown of the Patlak-Keller-Segel chemotaxis model

Chemotactic bacteria not only navigate chemical gradients, but also shape their environments by consuming and secreting attractants. Investigating how these processes influence the dynamics of bacterial populations has been challenging because of a lack of experimental methods for measuring spatial profiles of chemoattractants in real time. Here, we use a fluorescent sensor for aspartate to directly measure bacterially generated chemoattractant gradients during collective migration. Our measurements show that the standard Patlak-Keller-Segel model for collective chemotactic bacterial migration breaks down at high cell densities. To address this, we propose modifications to the model that consider the impact of cell density on bacterial chemotaxis and attractant consumption. With these changes, the model explains our experimental data across all cell densities, offering new insight into chemotactic dynamics. Our findings highlight the significance of considering cell density effects on bacterial behavior, and the potential for fluorescent metabolite sensors to shed light on the complex emergent dynamics of bacterial communities.

Revisiting Fur Regulon Leads to a Comprehensive Understanding of Iron and Fur Regulation

Iron is an essential element because it functions as a cofactor of many enzymes, but excess iron causes cell damage. Iron hemostasis in Escherichia coli was transcriptionally maintained by the ferric uptake regulator (Fur). Despite having been studied extensively, the comprehensive physiological roles and mechanisms of Fur-coordinated iron metabolism still remain obscure. In this work, by integrating a high-resolution transcriptomic study of the Fur wild-type and knockout Escherichia coli K-12 strains in the presence or absence of iron with high-throughput ChIP-seq assay and physiological studies, we revisited the regulatory roles of iron and Fur systematically and discovered several intriguing features of Fur regulation. The size of the Fur regulon was expanded greatly, and significant discrepancies were observed to exist between the regulations of Fur on the genes under its direct repression and activation. Fur showed stronger binding strength to the genes under its repression, and genes that were repressed by Fur were more sensitive to Fur and iron regulation as compared to the genes that were activated by Fur. Finally, we found that Fur linked iron metabolism to many essential processes, and the systemic regulations of Fur on carbon metabolism, respiration, and motility were further validated or discussed. These results highlight how Fur and Fur-controlled iron metabolism affect many cellular processes in a systematic way.

Dissociation kinetics of small-molecule inhibitors in Escherichia coli is coupled to physiological state of cells

Bioactive small-molecule inhibitors represent a treasure chest for future drugs. In vitro high-throughput screening is a common approach to identify the small-molecule inhibitors that bind tightly to purified targets. Here, we investigate the inhibitor-target binding/unbinding kinetics in E. coli cells using a benzimidazole-derivative DNA inhibitor as a model system. We find that its unbinding rate is not constant but depends on cell growth rate. This dependence is mediated by the cellular activity, forming a feedback loop with the inhibitor's activity. In accordance with this feedback, we find cell-to-cell heterogeneity in inhibitor-target interaction, leading to co-existence of two distinct subpopulations: actively growing cells that dissociate the inhibitors from the targets and non-growing cells that do not. We find similar heterogeneity for other clinical DNA inhibitors. Our studies reveal a mechanism that couples inhibitor-target kinetics to cell physiology and demonstrate the significant effect of this coupling on drug efficacy.

Navigating the CoGe Online Software Suite for Polyploidy Research

The CoGe software suite at genomevolution.org hosts a number of tools that facilitate genomic research on plant and animal whole-genome multiplication-polyploidy. SynMap permits analysis and visualization of two-way syntenic dotplot alignments of genomes, includes many options and data/graphics download possibilities, and even permits three-genome synteny maps and interactive views. FractBias is a tool that operates within SynMap that permits calculation and graphic display of genome fragments (such as chromosomes) of one species mapped to another, displaying both blockwise homology depths and the extent of syntenic gene (syntelog) loss following polyploidy events. SynMap macrosynteny results can segue into the microsynteny tool GEvo, which provides genome-browser-like views of homologous genome blocks. CoGe FeatView allows call-up of given gene features already stored in the CoGe resource, and CoGeBlast permits searches for additional features that can be analyzed or downloaded further. Links from these tools can be fed into SynFind, which can find syntenic blocks surrounding a feature across multiple specified genomes while also simultaneously providing overall genome-wide syntenic depth calculations that can be interpreted to reflect polyploidy levels. Here, we describe basic use of these tools on the CoGe software suite.

Tracking of Diversity and Evolution in the Brown Rot Fungi Monilinia fructicola, Monilinia fructigena, and Monilinia laxa

Monilinia species are among the most devastating fungi worldwide as they cause brown rot and blossom blight on fruit trees. To understand the molecular bases of their pathogenic lifestyles, we compared the newly assembled genomes of single strains of Monilinia fructicola, M. fructigena and M. laxa, with those of Botrytis cinerea and Sclerotinia sclerotiorum, as the closest species within Sclerotiniaceae. Phylogenomic analysis of orthologous proteins and syntenic investigation suggest that M. laxa is closer to M. fructigena than M. fructicola, and is closest to the other investigated Sclerotiniaceae species. This indicates that M. laxa was the earliest result of the speciation process. Distinct evolutionary profiles were observed for transposable elements (TEs). M. fructicola and M. laxa showed older bursts of TE insertions, which were affected (mainly in M. fructicola) by repeat-induced point (RIP) mutation gene silencing mechanisms. These suggested frequent occurrence of the sexual process in M. fructicola. More recent TE expansion linked with low RIP action was observed in M. fructigena, with very little in S. sclerotiorum and B. cinerea. The detection of active non-syntenic TEs is indicative of horizontal gene transfer and has resulted in alterations in specific gene functions. Analysis of candidate effectors, biosynthetic gene clusters for secondary metabolites and carbohydrate-active enzymes, indicated that Monilinia genus has multiple virulence mechanisms to infect host plants, including toxins, cell-death elicitor, putative virulence factors and cell-wall-degrading enzymes. Some species-specific pathogenic factors might explain differences in terms of host plant and organ preferences between M. fructigena and the other two Monilinia species.

Gain of Spontaneous clpX Mutations Boosting Motility via Adaption to Environments in Escherichia coli

Motility is finely regulated and is crucial to bacterial processes including colonization and biofilm formation. There is a trade-off between motility and growth in bacteria with molecular mechanisms not fully understood. Hypermotile Escherichia coli could be isolated by evolving non-motile cells on soft agar plates. Most of the isolates carried mutations located upstream of the flhDC promoter region, which upregulate the transcriptional expression of the master regulator of the flagellum biosynthesis, FlhDC. Here, we identified that spontaneous mutations in clpX boosted the motility of E. coli largely, inducing several folds of changes in swimming speed. Among the mutations identified, we further elucidated the molecular mechanism underlying the ClpX^V78F mutation on the regulation of E. coli motility. We found that the V78F mutation affected ATP binding to ClpX, resulting in the inability of the mutated ClpXP protease to degrade FlhD as indicated by both structure modeling and in vitro protein degradation assays. Moreover, our proteomic data indicated that the ClpX^V78F mutation elevated the stability of known ClpXP targets to various degrees with FlhD as one of the most affected. In addition, the specific tag at the C-terminus of FlhD being recognized for ClpXP degradation was identified. Finally, our transcriptome data characterized that the enhanced expression of the motility genes in the ClpX^V78F mutations was intrinsically accompanied by the reduced expression of stress resistance genes relating to the reduced fitness of the hypermotile strains. A similar pattern was observed for previously isolated hypermotile E. coli strains showing high expression of flhDC at the transcriptional level. Hence, clpX appears to be a hot locus comparable to the upstream of the flhDC promoter region evolved to boost bacterial motility, and our finding provides insight into the reduced fitness of the hypermotile bacteria.

Small RNA GcvB Regulates Oxidative Stress Response of Escherichia coli

Small non-translated regulatory RNAs control plenty of bacterial vital activities. The small RNA GcvB has been extensively studied, indicating the multifaceted roles of GcvB beyond amino acid metabolism. However, few reported GcvB-dependent regulation in minimal medium. Here, by applying a high-resolution RNA-seq assay, we compared the transcriptomes of a wild-type Escherichia coli K-12 strain and its gcvB deletion derivative grown in minimal medium and identified putative targets responding to GcvB, including flu, a determinant gene of auto-aggregation. The following molecular studies and the enhanced auto-aggregation ability of the gcvB knockout strain further substantiated the induced expression of these genes. Intriguingly, the reduced expression of OxyR (the oxidative stress regulator) in the gcvB knockout strain was identified to account for the increased expression of flu. Additionally, GcvB was characterized to up-regulate the expression of OxyR at the translational level. Accordingly, compared to the wild type, the GcvB deletion strain was more sensitive to oxidative stress and lost some its ability to eliminate endogenous reactive oxygen species. Taken together, we reveal that GcvB regulates oxidative stress response by up-regulating OxyR expression. Our findings provide an insight into the diversity of GcvB regulation and add an additional layer to the regulation of OxyR.

Active Efflux Leads to Heterogeneous Dissipation of Proton Motive Force by Protonophores in Bacteria

Various toxic compounds disrupt bacterial physiology. While bacteria harbor defense mechanisms to mitigate the toxicity, these mechanisms are often coupled to the physiological state of the cells and become ineffective when the physiology is severely disrupted. Here, we characterized such feedback by exposing Escherichia coli to protonophores. Protonophores dissipate the proton motive force (PMF), a fundamental force that drives physiological functions. We found that E. coli cells responded to protonophores heterogeneously, resulting in bimodal distributions of cell growth, substrate transport, and motility. Furthermore, we showed that this heterogeneous response required active efflux systems. The analysis of underlying interactions indicated the heterogeneous response results from efflux-mediated positive feedback between PMF and protonophores' action. Our studies have broad implications for bacterial adaptation to stress, including antibiotics. IMPORTANCE An electrochemical proton gradient across the cytoplasmic membrane, alternatively known as proton motive force, energizes vital cellular processes in bacteria, including ATP synthesis, nutrient uptake, and cell division. Therefore, a wide range of organisms produce the agents that collapse the proton motive force, protonophores, to gain a competitive advantage. Studies have shown that protonophores have significant effects on microbial competition, host-pathogen interaction, and antibiotic action and resistance. Furthermore, protonophores are extensively used in various laboratory studies to perturb bacterial physiology. Here, we have characterized cell growth, substrate transport, and motility of Escherichia coli cells exposed to protonophores. Our findings demonstrate heterogeneous effects of protonophores on cell physiology and the underlying mechanism.

From coarse to fine: the absolute Escherichia coli proteome under diverse growth conditions

Accurate measurements of cellular protein concentrations are invaluable to quantitative studies of gene expression and physiology in living cells. Here, we developed a versatile mass spectrometric workflow based on data-independent acquisition proteomics (DIA/SWATH) together with a novel protein inference algorithm (xTop). We used this workflow to accurately quantify absolute protein abundances in Escherichia coli for > 2,000 proteins over > 60 growth conditions, including nutrient limitations, non-metabolic stresses, and non-planktonic states. The resulting high-quality dataset of protein mass fractions allowed us to characterize proteome responses from a coarse (groups of related proteins) to a fine (individual) protein level. Hereby, a plethora of novel biological findings could be elucidated, including the generic upregulation of low-abundant proteins under various metabolic limitations, the non-specificity of catabolic enzymes upregulated under carbon limitation, the lack of large-scale proteome reallocation under stress compared to nutrient limitations, as well as surprising strain-dependent effects important for biofilm formation. These results present valuable resources for the systems biology community and can be used for future multi-omics studies of gene regulation and metabolic control in E. coli.

Global coordination of metabolic pathways in Escherichia coli by active and passive regulation

Microorganisms adjust metabolic activity to cope with diverse environments. While many studies have provided insights into how individual pathways are regulated, the mechanisms that give rise to coordinated metabolic responses are poorly understood. Here, we identify the regulatory mechanisms that coordinate catabolism and anabolism in Escherichia coli. Integrating protein, metabolite, and flux changes in genetically implemented catabolic or anabolic limitations, we show that combined global and local mechanisms coordinate the response to metabolic limitations. To allocate proteomic resources between catabolism and anabolism, E. coli uses a simple global gene regulatory program. Surprisingly, this program is largely implemented by a single transcription factor, Crp, which directly activates the expression of catabolic enzymes and indirectly reduces the expression of anabolic enzymes by passively sequestering cellular resources needed for their synthesis. However, metabolic fluxes are not controlled by this regulatory program alone; instead, fluxes are adjusted mostly through passive changes in the local metabolite concentrations. These mechanisms constitute a simple but effective global regulatory program that coarsely partitions resources between different parts of metabolism while ensuring robust coordination of individual metabolic reactions.

Metabolic engineering of Escherichia coli for the utilization of ethanol

BACKGROUND

The fuel ethanol industry has made tremendous progress in the last decades. Ethanol can be obtained by fermentation using a variety of biomass materials as the feedstocks. However, few studies have been conducted on ethanol utilization by microorganisms. The price of petroleum-derived ethanol, easily made by the hydrolysis of ethylene, is even lower than that of bioethanol. If ethanol can be metabolized by microorganisms to produce value-added chemicals, it will open a new door for the utilization of inexpensive ethanol resources.

RESULTS

We constructed an engineered Escherichia coli strain which could utilize ethanol as the sole carbon source. The alcohol dehydrogenase and aldehyde dehydrogenase from Aspergillus nidulans was introduced into E. coli and the recombinant strain acquired the ability to grow on ethanol. Cell growth continued when ethanol was supplied after glucose starvation and 2.24 g L-1 of ethanol was further consumed during the shake-flasks fermentation process. Then ethanol was further used for the production of mevalonic acid by heterologously expressing its biosynthetic pathway. Deuterium-labeled ethanol-D6 as the feedstock confirmed that mevalonic acid was synthesized from ethanol.

CONCLUSIONS

This study demonstrated the possibility of using ethanol as the carbon source by engineered E. coli strains. It can serve as the basis for the construction of more robust strains in the future though the catabolic capacity of ethanol should be further improved.

A universal trade-off between growth and lag in fluctuating environments

The rate of cell growth is crucial for bacterial fitness and drives the allocation of bacterial resources, affecting, for example, the expression levels of proteins dedicated to metabolism and biosynthesis1,2. It is unclear, however, what ultimately determines growth rates in different environmental conditions. Moreover, increasing evidence suggests that other objectives are also important3-7, such as the rate of physiological adaptation to changing environments8,9. A common challenge for cells is that these objectives cannot be independently optimized, and maximizing one often reduces another. Many such trade-offs have indeed been hypothesized on the basis of qualitative correlative studies8-11. Here we report a trade-off between steady-state growth rate and physiological adaptability in Escherichia coli, observed when a growing culture is abruptly shifted from a preferred carbon source such as glucose to fermentation products such as acetate. These metabolic transitions, common for enteric bacteria, are often accompanied by multi-hour lags before growth resumes. Metabolomic analysis reveals that long lags result from the depletion of key metabolites that follows the sudden reversal in the central carbon flux owing to the imposed nutrient shifts. A model of sequential flux limitation not only explains the observed trade-off between growth and adaptability, but also allows quantitative predictions regarding the universal occurrence of such tradeoffs, based on the opposing enzyme requirements of glycolysis versus gluconeogenesis. We validate these predictions experimentally for many different nutrient shifts in E. coli, as well as for other respiro-fermentative microorganisms, including Bacillus subtilis and Saccharomyces cerevisiae.

Chemotaxis as a navigation strategy to boost range expansion

Bacterial chemotaxis, the directed movement of cells along gradients of chemoattractants, is among the best-characterized subjects in molecular biology^1-10, but much less is known about its physiological roles¹¹. It is commonly seen as a starvation response when nutrients run out, or as an escape response from harmful situations^12-16. Here we identify an alternative role of chemotaxis by systematically examining the spatiotemporal dynamics of Escherichia coli in soft agar^12,17,18. Chemotaxis in nutrient-replete conditions promotes the expansion of bacterial populations into unoccupied territories well before nutrients run out in the current environment. Low levels of chemoattractants act as aroma-like cues in this process, establishing the direction and enhancing the speed of population movement along the self-generated attractant gradients. This process of navigated range expansion spreads faster and yields larger population gains than unguided expansion following the canonical Fisher-Kolmogorov dynamics^19,20 and is therefore a general strategy to promote population growth in spatially extended, nutrient-replete environments.

Conserved Pheromone Production, Response and Degradation by Streptococcus mutans

Streptococcus mutans, a bacterium with high cariogenic potential, coordinates competence for natural transformation and bacteriocin production via the XIP and CSP pheromones. CSP is effective in inducing bacteriocin responses but not competence in chemically defined media (CDM). This is in contrast to XIP, which is a strong inducer of competence in CDM but can also stimulate bacteriocin genes as a late response. Interconnections between the pathways activated by the two pheromones have been characterized in certain detail in S. mutans UA159, but it is mostly unknown whether such findings are representative for the species. In this study, we used bioassays based on luciferase reporters for the bacteriocin gene cipB and the alternative sigma factor sigX to investigate various S. mutans isolates for production and response to CSP and XIP pheromones in CDM. Similar to S. mutans UA159, endogenous CSP was undetectable in the culture supernatants of all tested strains. During optimization of the bioassay using the cipB reporter, we discovered that the activity of exogenous CSP used as a standard was reduced over time during S. mutans growth. Using a FRET-CSP reporter peptide, we found that S. mutans UA159 was able to degrade CSP, and that such activity was not significantly different in isogenic mutants with deletion of the protease gene htrA or the competence genes sigX, oppD, and comR. CSP cleavage was also detected in all the wild type strains, indicating that this is a conserved feature in S. mutans. For the XIP pheromone, endogenous production was observed in the supernatants of all 34 tested strains at peak concentrations in culture supernatants that varied between 200 and 26000 nM. Transformation in the presence of exogenous XIP was detected in all but one of the isolates. The efficiency of transformation varied, however, among the different strains, and for those with the highest transformation rates, endogenous XIP peak concentrations in the supernatants were above 2000 nM XIP. We conclude that XIP production and inducing effect on transformation, as well as the ability to degrade CSP, are conserved functions among different S. mutans isolates. Understanding the functionality and conservation of pheromone systems in S. mutans may lead to novel strategies to prevent or treat unbalances in oral microbiomes that may favor diseases.

Maintenance of translational elongation rate underlies the survival of Escherichia coli during oxidative stress

To cope with harsh circumstances, bacterial cells must initiate cellular stress response programs, which demands the de novo synthesis of many stress defense proteins. Reactive oxygen species (ROS) is a universal environmental stressor for both prokaryotic cells and eukaryotic cells. However, the physiological burden that limits the survival of bacterial cells during oxidative stress remains elusive. Here we quantitatively characterize the cell growth and translational elongation rate of Escherichia coli cells treated with different doses of hydrogen peroxide. Cell growth is immediately arrested by low to moderate levels of hydrogen peroxide, but completely recovers after a certain lag time. The lag time depends positively on the dose of hydrogen peroxide. During the lag time, translational elongation rate drops by as much as ∼90% at initial stage and recovers to its normal state later, a phenomenon resulting from the dramatic alteration in cellular tRNA pools during oxidative stress. However, translational elongation is completely stalled at a certain threshold-level of hydrogen peroxide, at which cells ultimately fail to resume growth. Although the mRNA transcription of oxidative defense genes in oxyR regulon is dramatically induced upon hydrogen peroxide treatment, the extreme slow-down of translational elongation during high levels of hydrogen peroxide has severely compromised the timely synthesis of those oxidative defense proteins. Our study demonstrates that the tRNA-limited translational elongation is a key physiological bottleneck that the bacteria must overcome to counteract ROS, and the maintenance of translational elongation rate for timely synthesis of stress defense proteins is crucial for cells to smoothly get over the oxidative stress.

Growth suppression by altered (p)ppGpp levels results from non-optimal resource allocation in Escherichia coli

Understanding how bacteria coordinate gene expression with biomass growth to adapt to various stress conditions remains a grand challenge in biology. Stress response is often associated with dramatic accumulation of cellular guanosine tetra- or penta-phosphate (p)ppGpp (also known as 'magic spot'), which is a key second messenger participating in regulating various biochemical and physiological processes of bacteria. Despite of the extensive studies on the mechanism of gene regulation by (p)ppGpp during stringent response, the connection between (p)ppGpp and bacterial steady-state exponential growth remains elusive. Here, we establish a versatile genetic approach to systematically perturb the (p)ppGpp level of Escherichia coli through titrating either the single-function (p)ppGpp synthetase or the singe-function (p)ppGpp hydrolase and quantitatively characterize cell growth and gene expression. Strikingly, increased and decreased (p)ppGpp levels both cause remarkable growth suppression of E. coli. From a coarse-grained insight, we demonstrate that increased (p)ppGpp levels limit ribosome synthesis while decreased (p)ppGpp levels limit the expression of metabolic proteins, both resulting in non-optimal resource allocation. Our study reveals a profound role of (p)ppGpp in regulating bacterial growth through governing global resource allocation. Moreover, we highlight the Mesh1 (p)ppGpp hydrolase from Drosophila melanogaster as a powerful genetic tool for interrogating bacterial (p)ppGpp physiology.

Mechanistic Origin of Cell-Size Control and Homeostasis in Bacteria

Evolutionarily divergent bacteria share a common phenomenological strategy for cell-size homeostasis under steady-state conditions. In the presence of inherent physiological stochasticity, cells following this "adder" principle gradually return to their steady-state size by adding a constant volume between birth and division, regardless of their size at birth. However, the mechanism of the adder has been unknown despite intense efforts. In this work, we show that the adder is a direct consequence of two general processes in biology: (1) threshold-accumulation of initiators and precursors required for cell division to a respective fixed number-and (2) balanced biosynthesis-maintenance of their production proportional to volume growth. This mechanism is naturally robust to static growth inhibition but also allows us to "reprogram" cell-size homeostasis in a quantitatively predictive manner in both Gram-negative Escherichia coli and Gram-positive Bacillus subtilis. By generating dynamic oscillations in the concentration of the division protein FtsZ, we were able to oscillate cell size at division and systematically break the adder. In contrast, periodic induction of replication initiator protein DnaA caused oscillations in cell size at initiation but did not alter division size or the adder. Finally, we were able to restore the adder phenotype in slow-growing E. coli, the only known steady-state growth condition wherein E. coli significantly deviates from the adder, by repressing active degradation of division proteins. Together, these results show that cell division and replication initiation are independently controlled at the gene-expression level and that division processes exclusively drive cell-size homeostasis in bacteria. VIDEO ABSTRACT.

Spatiotemporal establishment of dense bacterial colonies growing on hard agar

The physical interactions of growing bacterial cells with each other and with their surroundings significantly affect the structure and dynamics of biofilms. Here a 3D agent-based model is formulated to describe the establishment of simple bacterial colonies expanding by the physical force of their growth. With a single set of parameters, the model captures key dynamical features of colony growth by non-motile, non EPS-producing E. coli cells on hard agar. The model, supported by experiment on colony growth in different types and concentrations of nutrients, suggests that radial colony expansion is not limited by nutrients as commonly believed, but by mechanical forces. Nutrient penetration instead governs vertical colony growth, through thin layers of vertically oriented cells lifting up their ancestors from the bottom. Overall, the model provides a versatile platform to investigate the influences of metabolic and environmental factors on the growth and morphology of bacterial colonies.

Distinct mechanisms coordinate transcription and translation under carbon and nitrogen starvation in Escherichia coli

Bacteria adapt to environmental stress by producing proteins that provide stress protection. However, stress can severely perturb the kinetics of gene expression, disrupting protein production. Here, we characterized how Escherichia coli mitigates such perturbations under nutrient stress through the kinetic coordination of transcription and translation. We observed that, when translation became limiting under nitrogen starvation, transcription elongation slowed accordingly. This slowdown was mediated by (p)ppGpp, the alarmone whose primary role is thought to be promoter regulation. This kinetic coordination by (p)ppGpp was critical for the robust synthesis of gene products. Surprisingly, under carbon starvation, (p)ppGpp was dispensable for robust synthesis. Characterization of the underlying kinetics revealed that under carbon starvation, transcription became limiting, and translation aided transcription elongation. This mechanism naturally coordinated transcription with translation, alleviating the need for (p)ppGpp as a mediator. These contrasting mechanisms for coordination resulted in the condition-dependent effects of (p)ppGpp on global protein synthesis and starvation survival. Our findings reveal a kinetic aspect of gene expression plasticity, establishing (p)ppGpp as a condition-dependent global effector of gene expression.

Antibiotic-induced population fluctuations and stochastic clearance of bacteria

Effective antibiotic use that minimizes treatment failures remains a challenge. A better understanding of how bacterial populations respond to antibiotics is necessary. Previous studies of large bacterial populations established the deterministic framework of pharmacodynamics. Here, characterizing the dynamics of population extinction, we demonstrated the stochastic nature of eradicating bacteria with antibiotics. Antibiotics known to kill bacteria (bactericidal) induced population fluctuations. Thus, at high antibiotic concentrations, the dynamics of bacterial clearance were heterogeneous. At low concentrations, clearance still occurred with a non-zero probability. These striking outcomes of population fluctuations were well captured by our probabilistic model. Our model further suggested a strategy to facilitate eradication by increasing extinction probability. We experimentally tested this prediction for antibiotic-susceptible and clinically-isolated resistant bacteria. This new knowledge exposes fundamental limits in our ability to predict bacterial eradication. Additionally, it demonstrates the potential of using antibiotic concentrations that were previously deemed inefficacious to eradicate bacteria.

Phylogenomics of the olive tree (Olea europaea) reveals the relative contribution of ancient allo- and autopolyploidization events

BACKGROUND

Polyploidization is one of the major evolutionary processes that shape eukaryotic genomes, being particularly common in plants. Polyploids can arise through direct genome doubling within a species (autopolyploidization) or through the merging of genomes from distinct species after hybridization (allopolyploidization). The relative contribution of both mechanisms in plant evolution is debated. Here we used phylogenomics to dissect the tempo and mode of duplications in the genome of the olive tree (Olea europaea), one of the first domesticated Mediterranean fruit trees.

RESULTS

Our results depict a complex scenario involving at least three past polyploidization events, of which two-at the bases of the family Oleaceae and the tribe Oleeae, respectively-are likely to be the result of ancient allopolyploidization. A more recent polyploidization involves specifically the olive tree and relatives.

CONCLUSION

Our results show the power of phylogenomics to distinguish between allo- and auto polyploidization events and clarify the contributions of duplications in the evolutionary history of the olive tree.

The emergence of metabolic heterogeneity and diverse growth responses in isogenic bacterial cells

Microorganisms adapt to frequent environmental changes through population diversification. Previous studies demonstrated phenotypic diversity in a clonal population and its important effects on microbial ecology. However, the dynamic changes of phenotypic composition have rarely been characterized. Also, cellular variations and environmental factors responsible for phenotypic diversity remain poorly understood. Here, we studied phenotypic diversity driven by metabolic heterogeneity. We characterized metabolic activities and growth kinetics of starved Escherichia coli cells subject to nutrient upshift at single-cell resolution. We observed three subpopulations with distinct metabolic activities and growth phenotypes. One subpopulation was metabolically active and immediately grew upon nutrient upshift. One subpopulation was metabolically inactive and non-viable. The other subpopulation was metabolically partially active, and did not grow upon nutrient upshift. The ratio of these subpopulations changed dynamically during starvation. A long-term observation of cells with partial metabolic activities indicated that their metabolism was later spontaneously restored, leading to growth recovery. Further investigations showed that oxidative stress can induce the emergence of a subpopulation with partial metabolic activities. Our findings reveal the emergence of metabolic heterogeneity and associated dynamic changes in phenotypic composition. In addition, the results shed new light on microbial dormancy, which has important implications in microbial ecology and biomedicine.

Genome sequence and analysis of Escherichia coli production strain LS5218

Escherichia coli strain LS5218 is a useful host for the production of fatty acid derived products, but the genetics underlying this utility have not been fully investigated. Here, we report the genome sequence of LS5218 and a list of large mutations and single nucleotide permutations (SNPs) relative to E. coli K-12 strain MG1655. We discuss how genetic differences may affect the physiological differences between LS5218 and MG1655. We find that LS5218 is more closely related to E. coli strain NCM3722 and suspect that small genetic differences between K-12 derived strains may have a significant impact on metabolic engineering efforts.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination

Background

Root and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected “orphan” crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop.

Results

To accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage.

Conclusions

Guinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0419-x) contains supplementary material, which is available to authorized users.

Reduction of translating ribosomes enables Escherichia coli to maintain elongation rates during slow growth

Bacteria growing under different conditions experience a broad range of demand on the rate of protein synthesis, which profoundly affects cellular resource allocation. During fast growth, protein synthesis has long been known to be modulated by adjusting the ribosome content, with the vast majority of ribosomes engaged at a near-maximal rate of elongation. Here, we systematically characterize protein synthesis by Escherichia coli, focusing on slow-growth conditions. We establish that the translational elongation rate decreases as growth slows, exhibiting a Michaelis-Menten dependence on the abundance of the cellular translational apparatus. However, an appreciable elongation rate is maintained even towards zero growth, including the stationary phase. This maintenance, critical for timely protein synthesis in harsh environments, is accompanied by a drastic reduction in the fraction of active ribosomes. Interestingly, well-known antibiotics such as chloramphenicol also cause a substantial reduction in the pool of active ribosomes, instead of slowing down translational elongation as commonly thought.

The draft genome of the C3 panicoid grass species Dichanthelium oligosanthes

BACKGROUND

Comparisons between C3 and C4 grasses often utilize C3 species from the subfamilies Ehrhartoideae or Pooideae and C4 species from the subfamily Panicoideae, two clades that diverged over 50 million years ago. The divergence of the C3 panicoid grass Dichanthelium oligosanthes from the independent C4 lineages represented by Setaria viridis and Sorghum bicolor occurred approximately 15 million years ago, which is significantly more recent than members of the Bambusoideae, Ehrhartoideae, and Pooideae subfamilies. D. oligosanthes is ideally placed within the panicoid clade for comparative studies of C3 and C4 grasses.

RESULTS

We report the assembly of the nuclear and chloroplast genomes of D. oligosanthes, from high-throughput short read sequencing data and a comparative transcriptomics analysis of the developing leaf of D. oligosanthes, S. viridis, and S. bicolor. Physiological and anatomical characterizations verified that D. oligosanthes utilizes the C3 pathway for carbon fixation and lacks Kranz anatomy. Expression profiles of transcription factors along developing leaves of D. oligosanthes and S. viridis were compared with previously published data from S. bicolor, Zea mays, and Oryza sativa to identify a small suite of transcription factors that likely acquired functions specifically related to C4 photosynthesis.

CONCLUSIONS

The phylogenetic location of D. oligosanthes makes it an ideal C3 plant for comparative analysis of C4 evolution in the panicoid grasses. This genome will not only provide a better C3 species for comparisons with C4 panicoid grasses, but also highlights the power of using high-throughput sequencing to address questions in evolutionary biology.

Deciphering the Principles of Bacterial Nitrogen Dietary Preferences: a Strategy for Nutrient Containment

A fundamental question in microbial physiology concerns why organisms prefer certain nutrients to others. For example, among different nitrogen sources, ammonium is the preferred nitrogen source, supporting fast growth, whereas alternative nitrogen sources, such as certain amino acids, are considered to be poor nitrogen sources, supporting much slower exponential growth. However, the physiological/regulatory logic behind such nitrogen dietary choices remains elusive. In this study, by engineering Escherichia coli, we switched the dietary preferences toward amino acids, with growth rates equivalent to that of the wild-type strain grown on ammonia. However, when the engineered strain was cultured together with wild-type E. coli, this growth advantage was diminished as a consequence of ammonium leakage from the transport-and-catabolism (TC)-enhanced (TCE) cells, which are preferentially utilized by wild-type bacteria. Our results reveal that the nitrogen regulatory (Ntr) system fine tunes the expression of amino acid transport and catabolism components to match the flux through the ammonia assimilation pathway such that essential nutrients are retained, but, as a consequence, the fast growth rate on amino acids is sacrificed.

Bacteria exhibit different growth rates under various nutrient conditions. These environmentally related behaviors reflect the coordination between metabolism and the underlying regulatory networks. In the present study, we investigated the intertwined nitrogen metabolic and nitrogen regulatory systems to understand the growth differences between rich and poor nitrogen sources. Although maximal growth rate is considered to be evolutionarily advantageous for bacteria (as remarked by François Jacob, who said that the “dream” of every cell is to become two cells), we showed that negative-feedback loops in the regulatory system inhibit growth rates on amino acids. We demonstrated that in the absence of regulatory feedback, amino acids are capable of supporting fast growth rates, but this results in ammonia leaking out from cells as “waste,” benefiting the growth of competitors. These findings provide important insights into the regulatory logic that controls metabolic flux and ensures nutrient containment but consequently sacrifices growth rate.

Insertion sequence-caused large-scale rearrangements in the genome of Escherichia coli

A majority of large-scale bacterial genome rearrangements involve mobile genetic elements such as insertion sequence (IS) elements. Here we report novel insertions and excisions of IS elements and recombination between homologous IS elements identified in a large collection of Escherichia coli mutation accumulation lines by analysis of whole genome shotgun sequencing data. Based on 857 identified events (758 IS insertions, 98 recombinations and 1 excision), we estimate that the rate of IS insertion is 3.5 × 10(-4) insertions per genome per generation and the rate of IS homologous recombination is 4.5 × 10(-5) recombinations per genome per generation. These events are mostly contributed by the IS elements IS1, IS2, IS5 and IS186 Spatial analysis of new insertions suggest that transposition is biased to proximal insertions, and the length spectrum of IS-caused deletions is largely explained by local hopping. For any of the ISs studied there is no region of the circular genome that is favored or disfavored for new insertions but there are notable hotspots for deletions. Some elements have preferences for non-coding sequence or for the beginning and end of coding regions, largely explained by target site motifs. Interestingly, transposition and deletion rates remain constant across the wild-type and 12 mutant E. coli lines, each deficient in a distinct DNA repair pathway. Finally, we characterized the target sites of four IS families, confirming previous results and characterizing a highly specific pattern at IS186 target-sites, 5'-GGGG(N6/N7)CCCC-3'. We also detected 48 long deletions not involving IS elements.

Overflow metabolism in Escherichia coli results from efficient proteome allocation

Overflow metabolism refers to the seemingly wasteful strategy in which cells use fermentation instead of the more efficient respiration to generate energy, despite the availability of oxygen. Known as Warburg effect in the context of cancer growth, this phenomenon occurs ubiquitously for fast growing cells, including bacteria, fungi, and mammalian cells, but its origin has remained mysterious despite decades of research. Here we study metabolic overflow in E. coli and show that it is a global physiological response used to cope with changing proteomic demands of energy biogenesis and biomass synthesis under different growth conditions. A simple model of proteomic resource allocation can quantitatively account for all of the observed behaviors and accurately predict responses to novel perturbations. The key hypothesis of the model, that the proteome cost of energy biogenesis by respiration exceeds that by fermentation, is quantitatively confirmed by direct measurement of protein abundances via quantitative mass spectrometry.

Inflating bacterial cells by increased protein synthesis

Understanding how the homeostasis of cellular size and composition is accomplished by different organisms is an outstanding challenge in biology. For exponentially growing Escherichia coli cells, it is long known that the size of cells exhibits a strong positive relation with their growth rates in different nutrient conditions. Here, we characterized cell sizes in a set of orthogonal growth limitations. We report that cell size and mass exhibit positive or negative dependences with growth rate depending on the growth limitation applied. In particular, synthesizing large amounts of “useless” proteins led to an inversion of the canonical, positive relation, with slow growing cells enlarged 7- to 8-fold compared to cells growing at similar rates under nutrient limitation. Strikingly, this increase in cell size was accompanied by a 3- to 4-fold increase in cellular DNA content at slow growth, reaching up to an amount equivalent to ∼8 chromosomes per cell. Despite drastic changes in cell mass and macromolecular composition, cellular dry mass density remained constant. Our findings reveal an important role of protein synthesis in cell division control.

Complete Genome Sequence of Escherichia coli NCM3722

Escherichia coli NCM3722 is a prototrophic K-12 strain with robust physiologic phenotypes. We report the complete 4,678,045-bp chromosome and 67,545-bp F-like plasmid of this unique model organism.

Survival kinetics of starving bacteria is biphasic and density-dependent

In the lifecycle of microorganisms, prolonged starvation is prevalent and sustaining life during starvation periods is a vital task. In the literature, it is commonly assumed that survival kinetics of starving microbes follows exponential decay. This assumption, however, has not been rigorously tested. Currently, it is not clear under what circumstances this assumption is true. Also, it is not known when such survival kinetics deviates from exponential decay and if it deviates, what underlying mechanisms for the deviation are. Here, to address these issues, we quantitatively characterized dynamics of survival and death of starving E. coli cells. The results show that the assumption--starving cells die exponentially--is true only at high cell density. At low density, starving cells persevere for extended periods of time, before dying rapidly exponentially. Detailed analyses show intriguing quantitative characteristics of the density-dependent and biphasic survival kinetics, including that the period of the perseverance is inversely proportional to cell density. These characteristics further lead us to identification of key underlying processes relevant for the perseverance of starving cells. Then, using mathematical modeling, we show how these processes contribute to the density-dependent and biphasic survival kinetics observed. Importantly, our model reveals a thrifty strategy employed by bacteria, by which upon sensing impending depletion of a substrate, the limiting substrate is conserved and utilized later during starvation to delay cell death. These findings advance quantitative understanding of survival of microbes in oligotrophic environments and facilitate quantitative analysis and prediction of microbial dynamics in nature. Furthermore, they prompt revision of previous models used to analyze and predict population dynamics of microbes.

Polyploidy-associated genome modifications during land plant evolution

The occurrence of polyploidy in land plant evolution has led to an acceleration of genome modifications relative to other crown eukaryotes and is correlated with key innovations in plant evolution. Extensive genome resources provide for relating genomic changes to the origins of novel morphological and physiological features of plants. Ancestral gene contents for key nodes of the plant family tree are inferred. Pervasive polyploidy in angiosperms appears likely to be the major factor generating novel angiosperm genes and expanding some gene families. However, most gene families lose most duplicated copies in a quasi-neutral process, and a few families are actively selected for single-copy status. One of the great challenges of evolutionary genomics is to link genome modifications to speciation, diversification and the morphological and/or physiological innovations that collectively compose biodiversity. Rapid accumulation of genomic data and its ongoing investigation may greatly improve the resolution at which evolutionary approaches can contribute to the identification of specific genes responsible for particular innovations. The resulting, more 'particulate' understanding of plant evolution, may elevate to a new level fundamental knowledge of botanical diversity, including economically important traits in the crop plants that sustain humanity.

Mechanism for nitrogen isotope fractionation during ammonium assimilation by Escherichia coli K12

Organisms that use ammonium as the sole nitrogen source discriminate between [(15)N] and [(14)N] ammonium. This selectivity leaves an isotopic signature in their biomass that depends on the external concentration of ammonium. To dissect how differences in discrimination arise molecularly, we examined a wild-type (WT) strain of Escherichia coli K12 and mutant strains with lesions affecting ammonium-assimilatory proteins. We used isotope ratio mass spectrometry (MS) to assess the nitrogen isotopic composition of cell material when the strains were grown in batch culture at either high or low external concentrations of NH3 (achieved by controlling total NH4Cl and pH of the medium). At high NH3 (≥ 0.89 µM), discrimination against the heavy isotope by the WT strain (-19.2‰) can be accounted for by the equilibrium isotope effect for dissociation of NH4(+) to NH3 + H(+). NH3 equilibrates across the cytoplasmic membrane, and glutamine synthetase does not manifest an isotope effect in vivo. At low NH3 (≤ 0.18 µM), discrimination reflects an isotope effect for the NH4(+) channel AmtB (-14.1‰). By making E. coli dependent on the low-affinity ammonium-assimilatory pathway, we determined that biosynthetic glutamate dehydrogenase has an inverse isotope effect in vivo (+8.8‰). Likewise, by making unmediated diffusion of NH3 across the cytoplasmic membrane rate-limiting for cell growth in a mutant strain lacking AmtB, we could deduce an in vivo isotope effect for transport of NH3 across the membrane (-10.9‰). The paper presents the raw data from which our conclusions were drawn and discusses the assumptions underlying them.

A cost-effective and universal strategy for complete prokaryotic genomic sequencing proposed by computer simulation

BACKGROUND

Pyrosequencing techniques allow scientists to perform prokaryotic genome sequencing to achieve the draft genomic sequences within a few days. However, the assemblies with shotgun sequencing are usually composed of hundreds of contigs. A further multiplex PCR procedure is needed to fill all the gaps and link contigs into complete chromosomal sequence, which is the basis for prokaryotic comparative genomic studies. In this article, we study various pyrosequencing strategies by simulated assembling from 100 prokaryotic genomes.

FINDINGS

Simulation study shows that a single end 454 Jr. run combined with a paired end 454 Jr. run (8 kb library) can produce: 1) ~90% of 100 assemblies with < 10 scaffolds and ~95% of 100 assemblies with < 150 contigs; 2) average contig N50 size is over 331 kb; 3) average single base accuracy is > 99.99%; 4) average false gene duplication rate is < 0.7%; 5) average false gene loss rate is < 0.4%.

CONCLUSIONS

A single end 454 Jr. run combined with a paired end 454 Jr. run (8 kb library) is a cost-effective way for prokaryotic whole genome sequencing. This strategy provides solution to produce high quality draft assemblies for most of prokaryotic organisms within days. Due to the small number of assembled scaffolds, the following multiplex PCR procedure (for gap filling) would be easy. As a result, large scale prokaryotic whole genome sequencing projects may be finished within weeks.

Need-based activation of ammonium uptake in Escherichia coli

The efficient sequestration of nutrients is vital for the growth and survival of microorganisms. Some nutrients, such as CO2 and NH3, are readily diffusible across the cell membrane. The large membrane permeability of these nutrients obviates the need of transporters when the ambient level is high. When the ambient level is low, however, maintaining a high intracellular nutrient level against passive back diffusion is both challenging and costly. Here, we study the delicate management of ammonium (NH4+/NH3) sequestration by E. coli cells using microfluidic chemostats. We find that as the ambient ammonium concentration is reduced, E. coli cells first maximize their ability to assimilate the gaseous NH3 diffusing into the cytoplasm and then abruptly activate ammonium transport. The onset of transport varies under different growth conditions, but always occurring just as needed to maintain growth. Quantitative modeling of known interactions reveals an integral feedback mechanism by which this need-based uptake strategy is implemented. This novel strategy ensures that the expensive cost of upholding the internal ammonium concentration against back diffusion is kept at a minimum.

The iPlant Collaborative: Cyberinfrastructure for Plant Biology

The iPlant Collaborative (iPlant) is a United States National Science Foundation (NSF) funded project that aims to create an innovative, comprehensive, and foundational cyberinfrastructure in support of plant biology research (PSCIC, 2006). iPlant is developing cyberinfrastructure that uniquely enables scientists throughout the diverse fields that comprise plant biology to address Grand Challenges in new ways, to stimulate and facilitate cross-disciplinary research, to promote biology and computer science research interactions, and to train the next generation of scientists on the use of cyberinfrastructure in research and education. Meeting humanity's projected demands for agricultural and forest products and the expectation that natural ecosystems be managed sustainably will require synergies from the application of information technologies. The iPlant cyberinfrastructure design is based on an unprecedented period of research community input, and leverages developments in high-performance computing, data storage, and cyberinfrastructure for the physical sciences. iPlant is an open-source project with application programming interfaces that allow the community to extend the infrastructure to meet its needs. iPlant is sponsoring community-driven workshops addressing specific scientific questions via analysis tool integration and hypothesis testing. These workshops teach researchers how to add bioinformatics tools and/or datasets into the iPlant cyberinfrastructure enabling plant scientists to perform complex analyses on large datasets without the need to master the command-line or high-performance computational services.