Search Google Scholar | Search PubMed
84.
Willerth, S M; Pedro, H A; Pachter, L; Humeau, L M; Arkin, A P; Schaffer, D V
Đevelopment of a low bias method for characterizing viral populations using next generation sequencing technology Journal Article
In: PLoS ONE, vol. 5, no. 10, pp. e13564, 2010.
@article{pmid21042592,
title = {Đevelopment of a low bias method for characterizing viral populations using next generation sequencing technology},
author = {S M Willerth and H A Pedro and L Pachter and L M Humeau and A P Arkin and D V Schaffer},
year = {2010},
date = {2010-00-01},
journal = {PLoS ONE},
volume = {5},
number = {10},
pages = {e13564},
abstract = {With an estimated 38 million people worldwide currently infected with human immunodeficiency virus (HIV), and an additional 4.1 million people becoming infected each year, it is important to understand how this virus mutates and develops resistance in order to design successful therapies. We report a novel experimental method for amplifying full-length HIV genomes without the use of sequence-specific primers for high throughput DNA sequencing, followed by assembly of full length viral genome sequences from the resulting large dataset. Illumina was chosen for sequencing due to its ability to provide greater coverage of the HIV genome compared to prior methods, allowing for more comprehensive characterization of the heterogeneity present in the HIV samples analyzed. Our novel amplification method in combination with Illumina sequencing was used to analyze two HIV populations: a homogenous HIV population based on the canonical NL4-3 strain and a heterogeneous viral population obtained from a HIV patient's infected T cells. In addition, the resulting sequence was analyzed using a new computational approach to obtain a consensus sequence and several metrics of diversity. This study demonstrates how a lower bias amplification method in combination with next generation DNA sequencing provides in-depth, complete coverage of the HIV genome, enabling a stronger characterization of the quasispecies present in a clinically relevant HIV population as well as future study of how HIV mutates in response to a selective pressure.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
With an estimated 38 million people worldwide currently infected with human immunodeficiency virus (HIV), and an additional 4.1 million people becoming infected each year, it is important to understand how this virus mutates and develops resistance in order to design successful therapies. We report a novel experimental method for amplifying full-length HIV genomes without the use of sequence-specific primers for high throughput DNA sequencing, followed by assembly of full length viral genome sequences from the resulting large dataset. Illumina was chosen for sequencing due to its ability to provide greater coverage of the HIV genome compared to prior methods, allowing for more comprehensive characterization of the heterogeneity present in the HIV samples analyzed. Our novel amplification method in combination with Illumina sequencing was used to analyze two HIV populations: a homogenous HIV population based on the canonical NL4-3 strain and a heterogeneous viral population obtained from a HIV patient's infected T cells. In addition, the resulting sequence was analyzed using a new computational approach to obtain a consensus sequence and several metrics of diversity. This study demonstrates how a lower bias amplification method in combination with next generation DNA sequencing provides in-depth, complete coverage of the HIV genome, enabling a stronger characterization of the quasispecies present in a clinically relevant HIV population as well as future study of how HIV mutates in response to a selective pressure.
83.
Computational models of ĦIV-1 resistance to gene therapy elucidate therapy design principles Journal Article
In: PLoS Comput. Biol., vol. 6, no. 8, 2010.
@article{pmid20711350,
title = {Computational models of ĦIV-1 resistance to gene therapy elucidate therapy design principles},
year = {2010},
date = {2010-00-01},
journal = {PLoS Comput. Biol.},
volume = {6},
number = {8},
abstract = {Gene therapy is an emerging alternative to conventional anti-HIV-1 drugs, and can potentially control the virus while alleviating major limitations of current approaches. Yet, HIV-1's ability to rapidly acquire mutations and escape therapy presents a critical challenge to any novel treatment paradigm. Viral escape is thus a key consideration in the design of any gene-based technique. We develop a computational model of HIV's evolutionary dynamics in vivo in the presence of a genetic therapy to explore the impact of therapy parameters and strategies on the development of resistance. Our model is generic and captures the properties of a broad class of gene-based agents that inhibit early stages of the viral life cycle. We highlight the differences in viral resistance dynamics between gene and standard antiretroviral therapies, and identify key factors that impact long-term viral suppression. In particular, we underscore the importance of mutationally-induced viral fitness losses in cells that are not genetically modified, as these can severely constrain the replication of resistant virus. We also propose and investigate a novel treatment strategy that leverages upon gene therapy's unique capacity to deliver different genes to distinct cell populations, and we find that such a strategy can dramatically improve efficacy when used judiciously within a certain parametric regime. Finally, we revisit a previously-suggested idea of improving clinical outcomes by boosting the proliferation of the genetically-modified cells, but we find that such an approach has mixed effects on resistance dynamics. Our results provide insights into the short- and long-term effects of gene therapy and the role of its key properties in the evolution of resistance, which can serve as guidelines for the choice and optimization of effective therapeutic agents.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Gene therapy is an emerging alternative to conventional anti-HIV-1 drugs, and can potentially control the virus while alleviating major limitations of current approaches. Yet, HIV-1's ability to rapidly acquire mutations and escape therapy presents a critical challenge to any novel treatment paradigm. Viral escape is thus a key consideration in the design of any gene-based technique. We develop a computational model of HIV's evolutionary dynamics in vivo in the presence of a genetic therapy to explore the impact of therapy parameters and strategies on the development of resistance. Our model is generic and captures the properties of a broad class of gene-based agents that inhibit early stages of the viral life cycle. We highlight the differences in viral resistance dynamics between gene and standard antiretroviral therapies, and identify key factors that impact long-term viral suppression. In particular, we underscore the importance of mutationally-induced viral fitness losses in cells that are not genetically modified, as these can severely constrain the replication of resistant virus. We also propose and investigate a novel treatment strategy that leverages upon gene therapy's unique capacity to deliver different genes to distinct cell populations, and we find that such a strategy can dramatically improve efficacy when used judiciously within a certain parametric regime. Finally, we revisit a previously-suggested idea of improving clinical outcomes by boosting the proliferation of the genetically-modified cells, but we find that such an approach has mixed effects on resistance dynamics. Our results provide insights into the short- and long-term effects of gene therapy and the role of its key properties in the evolution of resistance, which can serve as guidelines for the choice and optimization of effective therapeutic agents.
82.
A universal ŦagModule collection for parallel genetic analysis of microorganisms Journal Article
In: Nucleic Acids Res., vol. 38, no. 14, pp. e146, 2010.
@article{pmid20494978,
title = {A universal ŦagModule collection for parallel genetic analysis of microorganisms},
year = {2010},
date = {2010-00-01},
journal = {Nucleic Acids Res.},
volume = {38},
number = {14},
pages = {e146},
abstract = {Systems-level analyses of non-model microorganisms are limited by the existence of numerous uncharacterized genes and a corresponding over-reliance on automated computational annotations. One solution to this challenge is to disrupt gene function using DNA tag technology, which has been highly successful in parallelizing reverse genetics in Saccharomyces cerevisiae and has led to discoveries in gene function, genetic interactions and drug mechanism of action. To extend the yeast DNA tag methodology to a wide variety of microorganisms and applications, we have created a universal, sequence-verified TagModule collection. A hallmark of the 4280 TagModules is that they are cloned into a Gateway entry vector, thus facilitating rapid transfer to any compatible genetic system. Here, we describe the application of the TagModules to rapidly generate tagged mutants by transposon mutagenesis in the metal-reducing bacterium Shewanella oneidensis MR-1 and the pathogenic yeast Candida albicans. Our results demonstrate the optimal hybridization properties of the TagModule collection, the flexibility in applying the strategy to diverse microorganisms and the biological insights that can be gained from fitness profiling tagged mutant collections. The publicly available TagModule collection is a platform-independent resource for the functional genomics of a wide range of microbial systems in the post-genome era.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Systems-level analyses of non-model microorganisms are limited by the existence of numerous uncharacterized genes and a corresponding over-reliance on automated computational annotations. One solution to this challenge is to disrupt gene function using DNA tag technology, which has been highly successful in parallelizing reverse genetics in Saccharomyces cerevisiae and has led to discoveries in gene function, genetic interactions and drug mechanism of action. To extend the yeast DNA tag methodology to a wide variety of microorganisms and applications, we have created a universal, sequence-verified TagModule collection. A hallmark of the 4280 TagModules is that they are cloned into a Gateway entry vector, thus facilitating rapid transfer to any compatible genetic system. Here, we describe the application of the TagModules to rapidly generate tagged mutants by transposon mutagenesis in the metal-reducing bacterium Shewanella oneidensis MR-1 and the pathogenic yeast Candida albicans. Our results demonstrate the optimal hybridization properties of the TagModule collection, the flexibility in applying the strategy to diverse microorganisms and the biological insights that can be gained from fitness profiling tagged mutant collections. The publicly available TagModule collection is a platform-independent resource for the functional genomics of a wide range of microbial systems in the post-genome era.
81.
Zhou, A; He, Z; Redding-Johanson, A M; Mukhopadhyay, A; Hemme, C L; Joachimiak, M P; Luo, F; Deng, Y; Bender, K S; He, Q; Keasling, J D; Stahl, D A; Fields, M W; Hazen, T C; Arkin, A P; Wall, J D; Zhou, J
Ħydrogen peroxide-induced oxidative stress responses in Đesulfovibrio vulgaris Ħildenborough Journal Article
In: Environ. Microbiol., vol. 12, no. 10, pp. 2645–2657, 2010.
@article{pmid20482586,
title = {Ħydrogen peroxide-induced oxidative stress responses in Đesulfovibrio vulgaris Ħildenborough},
author = {A Zhou and Z He and A M Redding-Johanson and A Mukhopadhyay and C L Hemme and M P Joachimiak and F Luo and Y Deng and K S Bender and Q He and J D Keasling and D A Stahl and M W Fields and T C Hazen and A P Arkin and J D Wall and J Zhou},
year = {2010},
date = {2010-00-01},
journal = {Environ. Microbiol.},
volume = {12},
number = {10},
pages = {2645--2657},
abstract = {To understand how sulphate-reducing bacteria respond to oxidative stresses, the responses of Desulfovibrio vulgaris Hildenborough to H(2)O(2)-induced stresses were investigated with transcriptomic, proteomic and genetic approaches. H(2)O(2) and induced chemical species (e.g. polysulfide, ROS) and redox potential shift increased the expressions of the genes involved in detoxification, thioredoxin-dependent reduction system, protein and DNA repair, and decreased those involved in sulfate reduction, lactate oxidation and protein synthesis. A gene coexpression network analysis revealed complicated network interactions among differentially expressed genes, and suggested possible importance of several hypothetical genes in H(2)O(2) stress. Also, most of the genes in PerR and Fur regulons were highly induced, and the abundance of a Fur regulon protein increased. Mutant analysis suggested that PerR and Fur are functionally overlapped in response to stresses induced by H(2)O(2) and reaction products, and the upregulation of thioredoxin-dependent reduction genes was independent of PerR or Fur. It appears that induction of those stress response genes could contribute to the increased resistance of deletion mutants to H(2)O(2)-induced stresses. In addition, a conceptual cellular model of D. vulgaris responses to H(2)O(2) stress was constructed to illustrate that this bacterium may employ a complicated molecular mechanism to defend against the H(2)O(2)-induced stresses.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
To understand how sulphate-reducing bacteria respond to oxidative stresses, the responses of Desulfovibrio vulgaris Hildenborough to H(2)O(2)-induced stresses were investigated with transcriptomic, proteomic and genetic approaches. H(2)O(2) and induced chemical species (e.g. polysulfide, ROS) and redox potential shift increased the expressions of the genes involved in detoxification, thioredoxin-dependent reduction system, protein and DNA repair, and decreased those involved in sulfate reduction, lactate oxidation and protein synthesis. A gene coexpression network analysis revealed complicated network interactions among differentially expressed genes, and suggested possible importance of several hypothetical genes in H(2)O(2) stress. Also, most of the genes in PerR and Fur regulons were highly induced, and the abundance of a Fur regulon protein increased. Mutant analysis suggested that PerR and Fur are functionally overlapped in response to stresses induced by H(2)O(2) and reaction products, and the upregulation of thioredoxin-dependent reduction genes was independent of PerR or Fur. It appears that induction of those stress response genes could contribute to the increased resistance of deletion mutants to H(2)O(2)-induced stresses. In addition, a conceptual cellular model of D. vulgaris responses to H(2)O(2) stress was constructed to illustrate that this bacterium may employ a complicated molecular mechanism to defend against the H(2)O(2)-induced stresses.
80.
Walker, C B; Stolyar, S; Chivian, D; Pinel, N; Gabster, J A; Dehal, P S; He, Z; Yang, Z K; Yen, H C; Zhou, J; Wall, J D; Hazen, T C; Arkin, A P; Stahl, D A
Contribution of mobile genetic elements to Đesulfovibrio vulgaris genome plasticity Journal Article
In: Environ. Microbiol., vol. 11, no. 9, pp. 2244–2252, 2009.
@article{pmid19737303,
title = {Contribution of mobile genetic elements to Đesulfovibrio vulgaris genome plasticity},
author = {C B Walker and S Stolyar and D Chivian and N Pinel and J A Gabster and P S Dehal and Z He and Z K Yang and H C Yen and J Zhou and J D Wall and T C Hazen and A P Arkin and D A Stahl},
year = {2009},
date = {2009-09-01},
journal = {Environ. Microbiol.},
volume = {11},
number = {9},
pages = {2244--2252},
abstract = {The genome of Desulfovibrio vulgaris strain DePue, a sulfate-reducing Deltaproteobacterium isolated from heavy metal-impacted lake sediment, was completely sequenced and compared with the type strain D. vulgaris Hildenborough. The two genomes share a high degree of relatedness and synteny, but harbour distinct prophage and signatures of past phage encounters. In addition to a highly variable phage contribution, the genome of strain DePue contains a cluster of open-reading frames not found in strain Hildenborough coding for the production and export of a capsule exopolysaccharide, possibly of relevance to heavy metal resistance. Comparative whole-genome microarray analysis on four additional D. vulgaris strains established greater interstrain variation within regions associated with phage insertion and exopolysaccharide biosynthesis.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
The genome of Desulfovibrio vulgaris strain DePue, a sulfate-reducing Deltaproteobacterium isolated from heavy metal-impacted lake sediment, was completely sequenced and compared with the type strain D. vulgaris Hildenborough. The two genomes share a high degree of relatedness and synteny, but harbour distinct prophage and signatures of past phage encounters. In addition to a highly variable phage contribution, the genome of strain DePue contains a cluster of open-reading frames not found in strain Hildenborough coding for the production and export of a capsule exopolysaccharide, possibly of relevance to heavy metal resistance. Comparative whole-genome microarray analysis on four additional D. vulgaris strains established greater interstrain variation within regions associated with phage insertion and exopolysaccharide biosynthesis.
79.
Walker, C B; He, Z; Yang, Z K; Ringbauer, J A; He, Q; Zhou, J; Voordouw, G; Wall, J D; Arkin, A P; Hazen, T C; Stolyar, S; Stahl, D A
Ŧhe electron transfer system of syntrophically grown Đesulfovibrio vulgaris Journal Article
In: J. Bacteriol., vol. 191, no. 18, pp. 5793–5801, 2009.
@article{pmid19581361,
title = {Ŧhe electron transfer system of syntrophically grown Đesulfovibrio vulgaris},
author = {C B Walker and Z He and Z K Yang and J A Ringbauer and Q He and J Zhou and G Voordouw and J D Wall and A P Arkin and T C Hazen and S Stolyar and D A Stahl},
year = {2009},
date = {2009-09-01},
journal = {J. Bacteriol.},
volume = {191},
number = {18},
pages = {5793--5801},
abstract = {Interspecies hydrogen transfer between organisms producing and consuming hydrogen promotes the decomposition of organic matter in most anoxic environments. Although syntrophic coupling between hydrogen producers and consumers is a major feature of the carbon cycle, mechanisms for energy recovery at the extremely low free energies of reactions typical of these anaerobic communities have not been established. In this study, comparative transcriptional analysis of a model sulfate-reducing microbe, Desulfovibrio vulgaris Hildenborough, suggested the use of alternative electron transfer systems dependent on growth modality. During syntrophic growth on lactate with a hydrogenotrophic methanogen, numerous genes involved in electron transfer and energy generation were upregulated in D. vulgaris compared with their expression in sulfate-limited monocultures. In particular, genes coding for the putative membrane-bound Coo hydrogenase, two periplasmic hydrogenases (Hyd and Hyn), and the well-characterized high-molecular-weight cytochrome (Hmc) were among the most highly expressed and upregulated genes. Additionally, a predicted operon containing genes involved in lactate transport and oxidation exhibited upregulation, further suggesting an alternative pathway for electrons derived from lactate oxidation during syntrophic growth. Mutations in a subset of genes coding for Coo, Hmc, Hyd, and Hyn impaired or severely limited syntrophic growth but had little effect on growth via sulfate respiration. These results demonstrate that syntrophic growth and sulfate respiration use largely independent energy generation pathways and imply that to understand microbial processes that sustain nutrient cycling, lifestyles not captured in pure culture must be considered.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Interspecies hydrogen transfer between organisms producing and consuming hydrogen promotes the decomposition of organic matter in most anoxic environments. Although syntrophic coupling between hydrogen producers and consumers is a major feature of the carbon cycle, mechanisms for energy recovery at the extremely low free energies of reactions typical of these anaerobic communities have not been established. In this study, comparative transcriptional analysis of a model sulfate-reducing microbe, Desulfovibrio vulgaris Hildenborough, suggested the use of alternative electron transfer systems dependent on growth modality. During syntrophic growth on lactate with a hydrogenotrophic methanogen, numerous genes involved in electron transfer and energy generation were upregulated in D. vulgaris compared with their expression in sulfate-limited monocultures. In particular, genes coding for the putative membrane-bound Coo hydrogenase, two periplasmic hydrogenases (Hyd and Hyn), and the well-characterized high-molecular-weight cytochrome (Hmc) were among the most highly expressed and upregulated genes. Additionally, a predicted operon containing genes involved in lactate transport and oxidation exhibited upregulation, further suggesting an alternative pathway for electrons derived from lactate oxidation during syntrophic growth. Mutations in a subset of genes coding for Coo, Hmc, Hyd, and Hyn impaired or severely limited syntrophic growth but had little effect on growth via sulfate respiration. These results demonstrate that syntrophic growth and sulfate respiration use largely independent energy generation pathways and imply that to understand microbial processes that sustain nutrient cycling, lifestyles not captured in pure culture must be considered.
78.
Price, M N; Dehal, P S; Arkin, A P
FastŦree: computing large minimum evolution trees with profiles instead of a distance matrix Journal Article
In: Mol. Biol. Evol., vol. 26, no. 7, pp. 1641–1650, 2009.
@article{pmid19377059,
title = {FastŦree: computing large minimum evolution trees with profiles instead of a distance matrix},
author = {M N Price and P S Dehal and A P Arkin},
year = {2009},
date = {2009-07-01},
journal = {Mol. Biol. Evol.},
volume = {26},
number = {7},
pages = {1641--1650},
abstract = {Gene families are growing rapidly, but standard methods for inferring phylogenies do not scale to alignments with over 10,000 sequences. We present FastTree, a method for constructing large phylogenies and for estimating their reliability. Instead of storing a distance matrix, FastTree stores sequence profiles of internal nodes in the tree. FastTree uses these profiles to implement Neighbor-Joining and uses heuristics to quickly identify candidate joins. FastTree then uses nearest neighbor interchanges to reduce the length of the tree. For an alignment with N sequences, L sites, and a different characters, a distance matrix requires O(N(2)) space and O(N(2)L) time, but FastTree requires just O(NLa + N ) memory and O(N log (N)La) time. To estimate the tree's reliability, FastTree uses local bootstrapping, which gives another 100-fold speedup over a distance matrix. For example, FastTree computed a tree and support values for 158,022 distinct 16S ribosomal RNAs in 17 h and 2.4 GB of memory. Just computing pairwise Jukes-Cantor distances and storing them, without inferring a tree or bootstrapping, would require 17 h and 50 GB of memory. In simulations, FastTree was slightly more accurate than Neighbor-Joining, BIONJ, or FastME; on genuine alignments, FastTree's topologies had higher likelihoods. FastTree is available at http://microbesonline.org/fasttree.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Gene families are growing rapidly, but standard methods for inferring phylogenies do not scale to alignments with over 10,000 sequences. We present FastTree, a method for constructing large phylogenies and for estimating their reliability. Instead of storing a distance matrix, FastTree stores sequence profiles of internal nodes in the tree. FastTree uses these profiles to implement Neighbor-Joining and uses heuristics to quickly identify candidate joins. FastTree then uses nearest neighbor interchanges to reduce the length of the tree. For an alignment with N sequences, L sites, and a different characters, a distance matrix requires O(N(2)) space and O(N(2)L) time, but FastTree requires just O(NLa + N ) memory and O(N log (N)La) time. To estimate the tree's reliability, FastTree uses local bootstrapping, which gives another 100-fold speedup over a distance matrix. For example, FastTree computed a tree and support values for 158,022 distinct 16S ribosomal RNAs in 17 h and 2.4 GB of memory. Just computing pairwise Jukes-Cantor distances and storing them, without inferring a tree or bootstrapping, would require 17 h and 50 GB of memory. In simulations, FastTree was slightly more accurate than Neighbor-Joining, BIONJ, or FastME; on genuine alignments, FastTree's topologies had higher likelihoods. FastTree is available at http://microbesonline.org/fasttree.
77.
Elias, D A; Mukhopadhyay, A; Joachimiak, M P; Drury, E C; Redding, A M; Yen, H C; Fields, M W; Hazen, T C; Arkin, A P; Keasling, J D; Wall, J D
Expression profiling of hypothetical genes in Đesulfovibrio vulgaris leads to improved functional annotation Journal Article
In: Nucleic Acids Res., vol. 37, no. 9, pp. 2926–2939, 2009.
@article{pmid19293273,
title = {Expression profiling of hypothetical genes in Đesulfovibrio vulgaris leads to improved functional annotation},
author = {D A Elias and A Mukhopadhyay and M P Joachimiak and E C Drury and A M Redding and H C Yen and M W Fields and T C Hazen and A P Arkin and J D Keasling and J D Wall},
year = {2009},
date = {2009-05-01},
journal = {Nucleic Acids Res.},
volume = {37},
number = {9},
pages = {2926--2939},
abstract = {Hypothetical (HyP) and conserved HyP genes account for >30% of sequenced bacterial genomes. For the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough, 347 of the 3634 genes were annotated as conserved HyP (9.5%) along with 887 HyP genes (24.4%). Given the large fraction of the genome, it is plausible that some of these genes serve critical cellular roles. The study goals were to determine which genes were expressed and provide a more functionally based annotation. To accomplish this, expression profiles of 1234 HyP and conserved genes were used from transcriptomic datasets of 11 environmental stresses, complemented with shotgun LC-MS/MS and AMT tag proteomic data. Genes were divided into putatively polycistronic operons and those predicted to be monocistronic, then classified by basal expression levels and grouped according to changes in expression for one or multiple stresses. One thousand two hundred and twelve of these genes were transcribed with 786 producing detectable proteins. There was no evidence for expression of 17 predicted genes. Except for the latter, monocistronic gene annotation was expanded using the above criteria along with matching Clusters of Orthologous Groups. Polycistronic genes were annotated in the same manner with inferences from their proximity to more confidently annotated genes. Two targeted deletion mutants were used as test cases to determine the relevance of the inferred functional annotations.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Hypothetical (HyP) and conserved HyP genes account for >30% of sequenced bacterial genomes. For the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough, 347 of the 3634 genes were annotated as conserved HyP (9.5%) along with 887 HyP genes (24.4%). Given the large fraction of the genome, it is plausible that some of these genes serve critical cellular roles. The study goals were to determine which genes were expressed and provide a more functionally based annotation. To accomplish this, expression profiles of 1234 HyP and conserved genes were used from transcriptomic datasets of 11 environmental stresses, complemented with shotgun LC-MS/MS and AMT tag proteomic data. Genes were divided into putatively polycistronic operons and those predicted to be monocistronic, then classified by basal expression levels and grouped according to changes in expression for one or multiple stresses. One thousand two hundred and twelve of these genes were transcribed with 786 producing detectable proteins. There was no evidence for expression of 17 predicted genes. Except for the latter, monocistronic gene annotation was expanded using the above criteria along with matching Clusters of Orthologous Groups. Polycistronic genes were annotated in the same manner with inferences from their proximity to more confidently annotated genes. Two targeted deletion mutants were used as test cases to determine the relevance of the inferred functional annotations.
76.
Yang, Y; Harris, D P; Luo, F; Xiong, W; Joachimiak, M; Wu, L; Dehal, P; Jacobsen, J; Yang, Z; Palumbo, A V; Arkin, A P; Zhou, J
Snapshot of iron response in Shewanella oneidensis by gene network reconstruction Journal Article
In: BMC Genomics, vol. 10, pp. 131, 2009.
@article{pmid19321007,
title = {Snapshot of iron response in Shewanella oneidensis by gene network reconstruction},
author = {Y Yang and D P Harris and F Luo and W Xiong and M Joachimiak and L Wu and P Dehal and J Jacobsen and Z Yang and A V Palumbo and A P Arkin and J Zhou},
year = {2009},
date = {2009-03-01},
journal = {BMC Genomics},
volume = {10},
pages = {131},
abstract = {Iron homeostasis of Shewanella oneidensis, a gamma-proteobacterium possessing high iron content, is regulated by a global transcription factor Fur. However, knowledge is incomplete about other biological pathways that respond to changes in iron concentration, as well as details of the responses. In this work, we integrate physiological, transcriptomics and genetic approaches to delineate the iron response of S. oneidensis. We show that the iron response in S. oneidensis is a rapid process. Temporal gene expression profiles were examined for iron depletion and repletion, and a gene co-expression network was reconstructed. Modules of iron acquisition systems, anaerobic energy metabolism and protein degradation were the most noteworthy in the gene network. Bioinformatics analyses suggested that genes in each of the modules might be regulated by DNA-binding proteins Fur, CRP and RpoH, respectively. Closer inspection of these modules revealed a transcriptional regulator (SO2426) involved in iron acquisition and ten transcriptional factors involved in anaerobic energy metabolism. Selected genes in the network were analyzed by genetic studies. Disruption of genes encoding a putative alcaligin biosynthesis protein (SO3032) and a gene previously implicated in protein degradation (SO2017) led to severe growth deficiency under iron depletion conditions. Disruption of a novel transcriptional factor (SO1415) caused deficiency in both anaerobic iron reduction and growth with thiosulfate or TMAO as an electronic acceptor, suggesting that SO1415 is required for specific branches of anaerobic energy metabolism pathways. Using a reconstructed gene network, we identified major biological pathways that were differentially expressed during iron depletion and repletion. Genetic studies not only demonstrated the importance of iron acquisition and protein degradation for iron depletion, but also characterized a novel transcriptional factor (SO1415) with a role in anaerobic energy metabolism.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Iron homeostasis of Shewanella oneidensis, a gamma-proteobacterium possessing high iron content, is regulated by a global transcription factor Fur. However, knowledge is incomplete about other biological pathways that respond to changes in iron concentration, as well as details of the responses. In this work, we integrate physiological, transcriptomics and genetic approaches to delineate the iron response of S. oneidensis. We show that the iron response in S. oneidensis is a rapid process. Temporal gene expression profiles were examined for iron depletion and repletion, and a gene co-expression network was reconstructed. Modules of iron acquisition systems, anaerobic energy metabolism and protein degradation were the most noteworthy in the gene network. Bioinformatics analyses suggested that genes in each of the modules might be regulated by DNA-binding proteins Fur, CRP and RpoH, respectively. Closer inspection of these modules revealed a transcriptional regulator (SO2426) involved in iron acquisition and ten transcriptional factors involved in anaerobic energy metabolism. Selected genes in the network were analyzed by genetic studies. Disruption of genes encoding a putative alcaligin biosynthesis protein (SO3032) and a gene previously implicated in protein degradation (SO2017) led to severe growth deficiency under iron depletion conditions. Disruption of a novel transcriptional factor (SO1415) caused deficiency in both anaerobic iron reduction and growth with thiosulfate or TMAO as an electronic acceptor, suggesting that SO1415 is required for specific branches of anaerobic energy metabolism pathways. Using a reconstructed gene network, we identified major biological pathways that were differentially expressed during iron depletion and repletion. Genetic studies not only demonstrated the importance of iron acquisition and protein degradation for iron depletion, but also characterized a novel transcriptional factor (SO1415) with a role in anaerobic energy metabolism.
75.
Jo, W J; Kim, J H; Oh, E; Jaramillo, D; Holman, P; Loguinov, A V; Arkin, A P; Nislow, C; Giaever, G; Vulpe, C D
Novel insights into iron metabolism by integrating deletome and transcriptome analysis in an iron deficiency model of the yeast Saccharomyces cerevisiae Journal Article
In: BMC Genomics, vol. 10, pp. 130, 2009.
@article{pmid19321002,
title = {Novel insights into iron metabolism by integrating deletome and transcriptome analysis in an iron deficiency model of the yeast Saccharomyces cerevisiae},
author = {W J Jo and J H Kim and E Oh and D Jaramillo and P Holman and A V Loguinov and A P Arkin and C Nislow and G Giaever and C D Vulpe},
year = {2009},
date = {2009-03-01},
journal = {BMC Genomics},
volume = {10},
pages = {130},
abstract = {Iron-deficiency anemia is the most prevalent form of anemia world-wide. The yeast Saccharomyces cerevisiae has been used as a model of cellular iron deficiency, in part because many of its cellular pathways are conserved. To better understand how cells respond to changes in iron availability, we profiled the yeast genome with a parallel analysis of homozygous deletion mutants to identify essential components and cellular processes required for optimal growth under iron-limited conditions. To complement this analysis, we compared those genes identified as important for fitness to those that were differentially-expressed in the same conditions. The resulting analysis provides a global perspective on the cellular processes involved in iron metabolism. Using functional profiling, we identified several genes known to be involved in high affinity iron uptake, in addition to novel genes that may play a role in iron metabolism. Our results provide support for the primary involvement in iron homeostasis of vacuolar and endosomal compartments, as well as vesicular transport to and from these compartments. We also observed an unexpected importance of the peroxisome for growth in iron-limited media. Although these components were essential for growth in low-iron conditions, most of them were not differentially-expressed. Genes with altered expression in iron deficiency were mainly associated with iron uptake and transport mechanisms, with little overlap with those that were functionally required. To better understand this relationship, we used expression-profiling of selected mutants that exhibited slow growth in iron-deficient conditions, and as a result, obtained additional insight into the roles of CTI6, DAP1, MRS4 and YHR045W in iron metabolism. Comparison between functional and gene expression data in iron deficiency highlighted the complementary utility of these two approaches to identify important functional components. This should be taken into consideration when designing and analyzing data from these type of studies. We used this and other published data to develop a molecular interaction network of iron metabolism in yeast.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Iron-deficiency anemia is the most prevalent form of anemia world-wide. The yeast Saccharomyces cerevisiae has been used as a model of cellular iron deficiency, in part because many of its cellular pathways are conserved. To better understand how cells respond to changes in iron availability, we profiled the yeast genome with a parallel analysis of homozygous deletion mutants to identify essential components and cellular processes required for optimal growth under iron-limited conditions. To complement this analysis, we compared those genes identified as important for fitness to those that were differentially-expressed in the same conditions. The resulting analysis provides a global perspective on the cellular processes involved in iron metabolism. Using functional profiling, we identified several genes known to be involved in high affinity iron uptake, in addition to novel genes that may play a role in iron metabolism. Our results provide support for the primary involvement in iron homeostasis of vacuolar and endosomal compartments, as well as vesicular transport to and from these compartments. We also observed an unexpected importance of the peroxisome for growth in iron-limited media. Although these components were essential for growth in low-iron conditions, most of them were not differentially-expressed. Genes with altered expression in iron deficiency were mainly associated with iron uptake and transport mechanisms, with little overlap with those that were functionally required. To better understand this relationship, we used expression-profiling of selected mutants that exhibited slow growth in iron-deficient conditions, and as a result, obtained additional insight into the roles of CTI6, DAP1, MRS4 and YHR045W in iron metabolism. Comparison between functional and gene expression data in iron deficiency highlighted the complementary utility of these two approaches to identify important functional components. This should be taken into consideration when designing and analyzing data from these type of studies. We used this and other published data to develop a molecular interaction network of iron metabolism in yeast.
74.
Skerker, J M; Lucks, J B; Arkin, A P
Evolution, ecology and the engineered organism: lessons for synthetic biology Journal Article
In: Genome Biol., vol. 10, no. 11, pp. 114, 2009.
@article{pmid19941672,
title = {Evolution, ecology and the engineered organism: lessons for synthetic biology},
author = {J M Skerker and J B Lucks and A P Arkin},
year = {2009},
date = {2009-01-01},
journal = {Genome Biol.},
volume = {10},
number = {11},
pages = {114},
abstract = {As the scope and complexity of synthetic biology grows, an understanding of evolution and ecology will be critical to its success.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
As the scope and complexity of synthetic biology grows, an understanding of evolution and ecology will be critical to its success.
73.
Burnett, J C; Miller-Jensen, K; Shah, P S; Arkin, A P; Schaffer, D V
Control of stochastic gene expression by host factors at the ĦIV promoter Journal Article
In: PLoS Pathog., vol. 5, no. 1, pp. e1000260, 2009.
@article{pmid19132086,
title = {Control of stochastic gene expression by host factors at the ĦIV promoter},
author = {J C Burnett and K Miller-Jensen and P S Shah and A P Arkin and D V Schaffer},
year = {2009},
date = {2009-01-01},
journal = {PLoS Pathog.},
volume = {5},
number = {1},
pages = {e1000260},
abstract = {The HIV promoter within the viral long terminal repeat (LTR) orchestrates many aspects of the viral life cycle, from the dynamics of viral gene expression and replication to the establishment of a latent state. In particular, after viral integration into the host genome, stochastic fluctuations in viral gene expression amplified by the Tat positive feedback loop can contribute to the formation of either a productive, transactivated state or an inactive state. In a significant fraction of cells harboring an integrated copy of the HIV-1 model provirus (LTR-GFP-IRES-Tat), this bimodal gene expression profile is dynamic, as cells spontaneously and continuously flip between active (Bright) and inactive (Off) expression modes. Furthermore, these switching dynamics may contribute to the establishment and maintenance of proviral latency, because after viral integration long delays in gene expression can occur before viral transactivation. The HIV-1 promoter contains cis-acting Sp1 and NF-kappaB elements that regulate gene expression via the recruitment of both activating and repressing complexes. We hypothesized that interplay in the recruitment of such positive and negative factors could modulate the stability of the Bright and Off modes and thereby alter the sensitivity of viral gene expression to stochastic fluctuations in the Tat feedback loop. Using model lentivirus variants with mutations introduced in the Sp1 and NF-kappaB elements, we employed flow cytometry, mRNA quantification, pharmacological perturbations, and chromatin immunoprecipitation to reveal significant functional differences in contributions of each site to viral gene regulation. Specifically, the Sp1 sites apparently stabilize both the Bright and the Off states, such that their mutation promotes noisy gene expression and reduction in the regulation of histone acetylation and deacetylation. Furthermore, the NF-kappaB sites exhibit distinct properties, with kappaB site I serving a stronger activating role than kappaB site II. Moreover, Sp1 site III plays a particularly important role in the recruitment of both p300 and RelA to the promoter. Finally, analysis of 362 clonal cell populations infected with the viral variants revealed that mutations in any of the Sp1 sites yield a 6-fold higher frequency of clonal bifurcation compared to that of the wild-type promoter. Thus, each Sp1 and NF-kappaB site differentially contributes to the regulation of viral gene expression, and Sp1 sites functionally "dampen" transcriptional noise and thereby modulate the frequency and maintenance of this model of viral latency. These results may have biomedical implications for the treatment of HIV latency.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
The HIV promoter within the viral long terminal repeat (LTR) orchestrates many aspects of the viral life cycle, from the dynamics of viral gene expression and replication to the establishment of a latent state. In particular, after viral integration into the host genome, stochastic fluctuations in viral gene expression amplified by the Tat positive feedback loop can contribute to the formation of either a productive, transactivated state or an inactive state. In a significant fraction of cells harboring an integrated copy of the HIV-1 model provirus (LTR-GFP-IRES-Tat), this bimodal gene expression profile is dynamic, as cells spontaneously and continuously flip between active (Bright) and inactive (Off) expression modes. Furthermore, these switching dynamics may contribute to the establishment and maintenance of proviral latency, because after viral integration long delays in gene expression can occur before viral transactivation. The HIV-1 promoter contains cis-acting Sp1 and NF-kappaB elements that regulate gene expression via the recruitment of both activating and repressing complexes. We hypothesized that interplay in the recruitment of such positive and negative factors could modulate the stability of the Bright and Off modes and thereby alter the sensitivity of viral gene expression to stochastic fluctuations in the Tat feedback loop. Using model lentivirus variants with mutations introduced in the Sp1 and NF-kappaB elements, we employed flow cytometry, mRNA quantification, pharmacological perturbations, and chromatin immunoprecipitation to reveal significant functional differences in contributions of each site to viral gene regulation. Specifically, the Sp1 sites apparently stabilize both the Bright and the Off states, such that their mutation promotes noisy gene expression and reduction in the regulation of histone acetylation and deacetylation. Furthermore, the NF-kappaB sites exhibit distinct properties, with kappaB site I serving a stronger activating role than kappaB site II. Moreover, Sp1 site III plays a particularly important role in the recruitment of both p300 and RelA to the promoter. Finally, analysis of 362 clonal cell populations infected with the viral variants revealed that mutations in any of the Sp1 sites yield a 6-fold higher frequency of clonal bifurcation compared to that of the wild-type promoter. Thus, each Sp1 and NF-kappaB site differentially contributes to the regulation of viral gene expression, and Sp1 sites functionally "dampen" transcriptional noise and thereby modulate the frequency and maintenance of this model of viral latency. These results may have biomedical implications for the treatment of HIV latency.
72.
Kazakov, A E; Rodionov, D A; Alm, E; Arkin, A P; Dubchak, I; Gelfand, M S
Comparative genomics of regulation of fatty acid and branched-chain amino acid utilization in proteobacteria Journal Article
In: J. Bacteriol., vol. 191, no. 1, pp. 52–64, 2009.
@article{pmid18820024,
title = {Comparative genomics of regulation of fatty acid and branched-chain amino acid utilization in proteobacteria},
author = {A E Kazakov and D A Rodionov and E Alm and A P Arkin and I Dubchak and M S Gelfand},
year = {2009},
date = {2009-01-01},
journal = {J. Bacteriol.},
volume = {191},
number = {1},
pages = {52--64},
abstract = {Bacteria can use branched-chain amino acids (ILV, i.e., isoleucine, leucine, valine) and fatty acids (FAs) as sole carbon and energy sources converting ILV into acetyl-coenzyme A (CoA), propanoyl-CoA, and propionyl-CoA, respectively. In this work, we used the comparative genomic approach to identify candidate transcriptional factors and DNA motifs that control ILV and FA utilization pathways in proteobacteria. The metabolic regulons were characterized based on the identification and comparison of candidate transcription factor binding sites in groups of phylogenetically related genomes. The reconstructed ILV/FA regulatory network demonstrates considerable variability and involves six transcriptional factors from the MerR, TetR, and GntR families binding to 11 distinct DNA motifs. The ILV degradation genes in gamma- and betaproteobacteria are regulated mainly by a novel regulator from the MerR family (e.g., LiuR in Pseudomonas aeruginosa) (40 species); in addition, the TetR-type regulator LiuQ was identified in some betaproteobacteria (eight species). Besides the core set of ILV utilization genes, the LiuR regulon in some lineages is expanded to include genes from other metabolic pathways, such as the glyoxylate shunt and glutamate synthase in Shewanella species. The FA degradation genes are controlled by four regulators including FadR in gammaproteobacteria (34 species), PsrA in gamma- and betaproteobacteria (45 species), FadP in betaproteobacteria (14 species), and LiuR orthologs in alphaproteobacteria (22 species). The remarkable variability of the regulatory systems associated with the FA degradation pathway is discussed from functional and evolutionary points of view.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Bacteria can use branched-chain amino acids (ILV, i.e., isoleucine, leucine, valine) and fatty acids (FAs) as sole carbon and energy sources converting ILV into acetyl-coenzyme A (CoA), propanoyl-CoA, and propionyl-CoA, respectively. In this work, we used the comparative genomic approach to identify candidate transcriptional factors and DNA motifs that control ILV and FA utilization pathways in proteobacteria. The metabolic regulons were characterized based on the identification and comparison of candidate transcription factor binding sites in groups of phylogenetically related genomes. The reconstructed ILV/FA regulatory network demonstrates considerable variability and involves six transcriptional factors from the MerR, TetR, and GntR families binding to 11 distinct DNA motifs. The ILV degradation genes in gamma- and betaproteobacteria are regulated mainly by a novel regulator from the MerR family (e.g., LiuR in Pseudomonas aeruginosa) (40 species); in addition, the TetR-type regulator LiuQ was identified in some betaproteobacteria (eight species). Besides the core set of ILV utilization genes, the LiuR regulon in some lineages is expanded to include genes from other metabolic pathways, such as the glyoxylate shunt and glutamate synthase in Shewanella species. The FA degradation genes are controlled by four regulators including FadR in gammaproteobacteria (34 species), PsrA in gamma- and betaproteobacteria (45 species), FadP in betaproteobacteria (14 species), and LiuR orthologs in alphaproteobacteria (22 species). The remarkable variability of the regulatory systems associated with the FA degradation pathway is discussed from functional and evolutionary points of view.
71.
Bischofs, I B; Hug, J A; Liu, A W; Wolf, D M; Arkin, A P
Complexity in bacterial cell-cell communication: quorum signal integration and subpopulation signaling in the Bacillus subtilis phosphorelay Journal Article
In: Proc. Natl. Acad. Sci. U.S.A., vol. 106, no. 16, pp. 6459–6464, 2009.
@article{pmid19380751,
title = {Complexity in bacterial cell-cell communication: quorum signal integration and subpopulation signaling in the Bacillus subtilis phosphorelay},
author = {I B Bischofs and J A Hug and A W Liu and D M Wolf and A P Arkin},
year = {2009},
date = {2009-00-01},
journal = {Proc. Natl. Acad. Sci. U.S.A.},
volume = {106},
number = {16},
pages = {6459--6464},
abstract = {A common form of quorum sensing in gram-positive bacteria is mediated by peptides that act as phosphatase regulators (Phr) of receptor aspartyl phosphatases (Raps). In Bacillus subtilis, several Phr signals are integrated in sporulation phosphorelay signal transduction. We theoretically demonstrate that the phosphorelay can act as a computational machine performing a sensitive division operation of kinase-encoded signals by quorum-modulated Rap signals, indicative of cells computing a "food per cell" estimate to decide whether to enter sporulation. We predict expression from the rapA-phrA operon to bifurcate as relative environmental signals change in a developing population. We experimentally observe that the rapA-phrA operon is heterogeneously induced in sporulating microcolonies. Uninduced cells sporulate rather synchronously early on, whereas the RapA/PhrA subpopulation sporulates less synchronously throughout later stationary phase. Moreover, we show that cells sustain PhrA expression during periods of active growth. Together with the model, these findings suggest that the phosphorelay may normalize environmental signals by the size of the (sub)population actively competing for nutrients (as signaled by PhrA). Generalizing this concept, the various Phrs could facilitate subpopulation communication in dense isogenic communities to control the physiological strategies followed by differentiated subpopulations by interpreting (environmental) signals based on the spatiotemporal community structure.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
A common form of quorum sensing in gram-positive bacteria is mediated by peptides that act as phosphatase regulators (Phr) of receptor aspartyl phosphatases (Raps). In Bacillus subtilis, several Phr signals are integrated in sporulation phosphorelay signal transduction. We theoretically demonstrate that the phosphorelay can act as a computational machine performing a sensitive division operation of kinase-encoded signals by quorum-modulated Rap signals, indicative of cells computing a "food per cell" estimate to decide whether to enter sporulation. We predict expression from the rapA-phrA operon to bifurcate as relative environmental signals change in a developing population. We experimentally observe that the rapA-phrA operon is heterogeneously induced in sporulating microcolonies. Uninduced cells sporulate rather synchronously early on, whereas the RapA/PhrA subpopulation sporulates less synchronously throughout later stationary phase. Moreover, we show that cells sustain PhrA expression during periods of active growth. Together with the model, these findings suggest that the phosphorelay may normalize environmental signals by the size of the (sub)population actively competing for nutrients (as signaled by PhrA). Generalizing this concept, the various Phrs could facilitate subpopulation communication in dense isogenic communities to control the physiological strategies followed by differentiated subpopulations by interpreting (environmental) signals based on the spatiotemporal community structure.
70.
Ross, J; Arkin, A P
Complex systems: from chemistry to systems biology Journal Article
In: Proc. Natl. Acad. Sci. U.S.A., vol. 106, no. 16, pp. 6433–6434, 2009.
BibTeX | Tags:
@article{pmid19380716,
title = {Complex systems: from chemistry to systems biology},
author = {J Ross and A P Arkin},
year = {2009},
date = {2009-00-01},
journal = {Proc. Natl. Acad. Sci. U.S.A.},
volume = {106},
number = {16},
pages = {6433--6434},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
69.
Lucks, J B; Qi, L; Whitaker, W R; Arkin, A P
Ŧoward scalable parts families for predictable design of biological circuits Journal Article
In: Curr. Opin. Microbiol., vol. 11, no. 6, pp. 567–573, 2008.
@article{pmid18983935,
title = {Ŧoward scalable parts families for predictable design of biological circuits},
author = {J B Lucks and L Qi and W R Whitaker and A P Arkin},
year = {2008},
date = {2008-12-01},
journal = {Curr. Opin. Microbiol.},
volume = {11},
number = {6},
pages = {567--573},
abstract = {Our current ability to engineer biological circuits is hindered by design cycles that are costly in terms of time and money, with constructs failing to operate as desired, or evolving away from the desired function once deployed. Synthetic biologists seek to understand biological design principles and use them to create technologies that increase the efficiency of the genetic engineering design cycle. Central to the approach is the creation of biological parts--encapsulated functions that can be composited together to create new pathways with predictable behaviors. We define five desirable characteristics of biological parts--independence, reliability, tunability, orthogonality and composability, and review studies of small natural and synthetic biological circuits that provide insights into each of these characteristics. We propose that the creation of appropriate sets of families of parts with these properties is a prerequisite for efficient, predictable engineering of new function in cells and will enable a large increase in the sophistication of genetic engineering applications.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Our current ability to engineer biological circuits is hindered by design cycles that are costly in terms of time and money, with constructs failing to operate as desired, or evolving away from the desired function once deployed. Synthetic biologists seek to understand biological design principles and use them to create technologies that increase the efficiency of the genetic engineering design cycle. Central to the approach is the creation of biological parts--encapsulated functions that can be composited together to create new pathways with predictable behaviors. We define five desirable characteristics of biological parts--independence, reliability, tunability, orthogonality and composability, and review studies of small natural and synthetic biological circuits that provide insights into each of these characteristics. We propose that the creation of appropriate sets of families of parts with these properties is a prerequisite for efficient, predictable engineering of new function in cells and will enable a large increase in the sophistication of genetic engineering applications.
68.
Flaherty, P; Radhakrishnan, M L; Dinh, T; Rebres, R A; Roach, T I; Jordan, M I; Arkin, A P
A dual receptor crosstalk model of G-protein-coupled signal transduction Journal Article
In: PLoS Comput. Biol., vol. 4, no. 9, pp. e1000185, 2008.
@article{pmid18818727,
title = {A dual receptor crosstalk model of G-protein-coupled signal transduction},
author = {P Flaherty and M L Radhakrishnan and T Dinh and R A Rebres and T I Roach and M I Jordan and A P Arkin},
year = {2008},
date = {2008-09-01},
journal = {PLoS Comput. Biol.},
volume = {4},
number = {9},
pages = {e1000185},
abstract = {Macrophage cells that are stimulated by two different ligands that bind to G-protein-coupled receptors (GPCRs) usually respond as if the stimulus effects are additive, but for a minority of ligand combinations the response is synergistic. The G-protein-coupled receptor system integrates signaling cues from the environment to actuate cell morphology, gene expression, ion homeostasis, and other physiological states. We analyze the effects of the two signaling molecules complement factors 5a (C5a) and uridine diphosphate (UDP) on the intracellular second messenger calcium to elucidate the principles that govern the processing of multiple signals by GPCRs. We have developed a formal hypothesis, in the form of a kinetic model, for the mechanism of action of this GPCR signal transduction system using data obtained from RAW264.7 macrophage cells. Bayesian statistical methods are employed to represent uncertainty in both data and model parameters and formally tie the model to experimental data. When the model is also used as a tool in the design of experiments, it predicts a synergistic region in the calcium peak height dose response that results when cells are simultaneously stimulated by C5a and UDP. An analysis of the model reveals a potential mechanism for crosstalk between the Galphai-coupled C5a receptor and the Galphaq-coupled UDP receptor signaling systems that results in synergistic calcium release.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Macrophage cells that are stimulated by two different ligands that bind to G-protein-coupled receptors (GPCRs) usually respond as if the stimulus effects are additive, but for a minority of ligand combinations the response is synergistic. The G-protein-coupled receptor system integrates signaling cues from the environment to actuate cell morphology, gene expression, ion homeostasis, and other physiological states. We analyze the effects of the two signaling molecules complement factors 5a (C5a) and uridine diphosphate (UDP) on the intracellular second messenger calcium to elucidate the principles that govern the processing of multiple signals by GPCRs. We have developed a formal hypothesis, in the form of a kinetic model, for the mechanism of action of this GPCR signal transduction system using data obtained from RAW264.7 macrophage cells. Bayesian statistical methods are employed to represent uncertainty in both data and model parameters and formally tie the model to experimental data. When the model is also used as a tool in the design of experiments, it predicts a synergistic region in the calcium peak height dose response that results when cells are simultaneously stimulated by C5a and UDP. An analysis of the model reveals a potential mechanism for crosstalk between the Galphai-coupled C5a receptor and the Galphaq-coupled UDP receptor signaling systems that results in synergistic calcium release.
67.
Ham, T S; Lee, S K; Keasling, J D; Arkin, A P
Đesign and construction of a double inversion recombination switch for heritable sequential genetic memory Journal Article
In: PLoS ONE, vol. 3, no. 7, pp. e2815, 2008.
@article{pmid18665232,
title = {Đesign and construction of a double inversion recombination switch for heritable sequential genetic memory},
author = {T S Ham and S K Lee and J D Keasling and A P Arkin},
year = {2008},
date = {2008-07-01},
journal = {PLoS ONE},
volume = {3},
number = {7},
pages = {e2815},
abstract = {Inversion recombination elements present unique opportunities for computing and information encoding in biological systems. They provide distinct binary states that are encoded into the DNA sequence itself, allowing us to overcome limitations posed by other biological memory or logic gate systems. Further, it is in theory possible to create complex sequential logics by careful positioning of recombinase recognition sites in the sequence. In this work, we describe the design and synthesis of an inversion switch using the fim and hin inversion recombination systems to create a heritable sequential memory switch. We have integrated the two inversion systems in an overlapping manner, creating a switch that can have multiple states. The switch is capable of transitioning from state to state in a manner analogous to a finite state machine, while encoding the state information into DNA. This switch does not require protein expression to maintain its state, and "remembers" its state even upon cell death. We were able to demonstrate transition into three out of the five possible states showing the feasibility of such a switch. We demonstrate that a heritable memory system that encodes its state into DNA is possible, and that inversion recombination system could be a starting point for more complex memory circuits. Although the circuit did not fully behave as expected, we showed that a multi-state, temporal memory is achievable.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Inversion recombination elements present unique opportunities for computing and information encoding in biological systems. They provide distinct binary states that are encoded into the DNA sequence itself, allowing us to overcome limitations posed by other biological memory or logic gate systems. Further, it is in theory possible to create complex sequential logics by careful positioning of recombinase recognition sites in the sequence. In this work, we describe the design and synthesis of an inversion switch using the fim and hin inversion recombination systems to create a heritable sequential memory switch. We have integrated the two inversion systems in an overlapping manner, creating a switch that can have multiple states. The switch is capable of transitioning from state to state in a manner analogous to a finite state machine, while encoding the state information into DNA. This switch does not require protein expression to maintain its state, and "remembers" its state even upon cell death. We were able to demonstrate transition into three out of the five possible states showing the feasibility of such a switch. We demonstrate that a heritable memory system that encodes its state into DNA is possible, and that inversion recombination system could be a starting point for more complex memory circuits. Although the circuit did not fully behave as expected, we showed that a multi-state, temporal memory is achievable.
66.
Singh, A H; Wolf, D M; Wang, P; Arkin, A P
Modularity of stress response evolution Journal Article
In: Proc. Natl. Acad. Sci. U.S.A., vol. 105, no. 21, pp. 7500–7505, 2008.
@article{pmid18495925,
title = {Modularity of stress response evolution},
author = {A H Singh and D M Wolf and P Wang and A P Arkin},
year = {2008},
date = {2008-05-01},
journal = {Proc. Natl. Acad. Sci. U.S.A.},
volume = {105},
number = {21},
pages = {7500--7505},
abstract = {Responses to extracellular stress directly confer survival fitness by means of complex regulatory networks. Despite their complexity, the networks must be evolvable because of changing ecological and environmental pressures. Although the regulatory networks underlying stress responses are characterized extensively, their mechanism of evolution remains poorly understood. Here, we examine the evolution of three candidate stress response networks (chemotaxis, competence for DNA uptake, and endospore formation) by analyzing their phylogenetic distribution across several hundred diverse bacterial and archaeal lineages. We report that genes in the chemotaxis and sporulation networks group into well defined evolutionary modules with distinct functions, phenotypes, and substitution rates as compared with control sets of randomly chosen genes. The evolutionary modules vary in both number and cohesiveness among the three pathways. Chemotaxis has five coherent modules whose distribution among species shows a clear pattern of interdependence and rewiring. Sporulation, by contrast, is nearly monolithic and seems to be inherited vertically, with three weak modules constituting early and late stages of the pathway. Competence does not seem to exhibit well defined modules either at or below the pathway level. Many of the detected modules are better understood in engineering terms than in protein functional terms, as we demonstrate using a control-based ontology that classifies gene function according to roles such as "sensor," "regulator," and äctuator." Moreover, we show that combinations of the modules predict phenotype, yet surprisingly do not necessarily correlate with phylogenetic inheritance. The architectures of these three pathways are therefore emblematic of different modes and constraints on evolution.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Responses to extracellular stress directly confer survival fitness by means of complex regulatory networks. Despite their complexity, the networks must be evolvable because of changing ecological and environmental pressures. Although the regulatory networks underlying stress responses are characterized extensively, their mechanism of evolution remains poorly understood. Here, we examine the evolution of three candidate stress response networks (chemotaxis, competence for DNA uptake, and endospore formation) by analyzing their phylogenetic distribution across several hundred diverse bacterial and archaeal lineages. We report that genes in the chemotaxis and sporulation networks group into well defined evolutionary modules with distinct functions, phenotypes, and substitution rates as compared with control sets of randomly chosen genes. The evolutionary modules vary in both number and cohesiveness among the three pathways. Chemotaxis has five coherent modules whose distribution among species shows a clear pattern of interdependence and rewiring. Sporulation, by contrast, is nearly monolithic and seems to be inherited vertically, with three weak modules constituting early and late stages of the pathway. Competence does not seem to exhibit well defined modules either at or below the pathway level. Many of the detected modules are better understood in engineering terms than in protein functional terms, as we demonstrate using a control-based ontology that classifies gene function according to roles such as "sensor," "regulator," and äctuator." Moreover, we show that combinations of the modules predict phenotype, yet surprisingly do not necessarily correlate with phylogenetic inheritance. The architectures of these three pathways are therefore emblematic of different modes and constraints on evolution.
65.
Wolf, D M; Fontaine-Bodin, L; Bischofs, I; Price, G; Keasling, J; Arkin, A P
Memory in microbes: quantifying history-dependent behavior in a bacterium Journal Article
In: PLoS ONE, vol. 3, no. 2, pp. e1700, 2008.
@article{pmid18324309,
title = {Memory in microbes: quantifying history-dependent behavior in a bacterium},
author = {D M Wolf and L Fontaine-Bodin and I Bischofs and G Price and J Keasling and A P Arkin},
year = {2008},
date = {2008-02-01},
journal = {PLoS ONE},
volume = {3},
number = {2},
pages = {e1700},
abstract = {Memory is usually associated with higher organisms rather than bacteria. However, evidence is mounting that many regulatory networks within bacteria are capable of complex dynamics and multi-stable behaviors that have been linked to memory in other systems. Moreover, it is recognized that bacteria that have experienced different environmental histories may respond differently to current conditions. These "memory" effects may be more than incidental to the regulatory mechanisms controlling acclimation or to the status of the metabolic stores. Rather, they may be regulated by the cell and confer fitness to the organism in the evolutionary game it participates in. Here, we propose that history-dependent behavior is a potentially important manifestation of memory, worth classifying and quantifying. To this end, we develop an information-theory based conceptual framework for measuring both the persistence of memory in microbes and the amount of information about the past encoded in history-dependent dynamics. This method produces a phenomenological measure of cellular memory without regard to the specific cellular mechanisms encoding it. We then apply this framework to a strain of Bacillus subtilis engineered to report on commitment to sporulation and degradative enzyme (AprE) synthesis and estimate the capacity of these systems and growth dynamics to 'remember' 10 distinct cell histories prior to application of a common stressor. The analysis suggests that B. subtilis remembers, both in short and long term, aspects of its cell history, and that this memory is distributed differently among the observables. While this study does not examine the mechanistic bases for memory, it presents a framework for quantifying memory in cellular behaviors and is thus a starting point for studying new questions about cellular regulation and evolutionary strategy.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Memory is usually associated with higher organisms rather than bacteria. However, evidence is mounting that many regulatory networks within bacteria are capable of complex dynamics and multi-stable behaviors that have been linked to memory in other systems. Moreover, it is recognized that bacteria that have experienced different environmental histories may respond differently to current conditions. These "memory" effects may be more than incidental to the regulatory mechanisms controlling acclimation or to the status of the metabolic stores. Rather, they may be regulated by the cell and confer fitness to the organism in the evolutionary game it participates in. Here, we propose that history-dependent behavior is a potentially important manifestation of memory, worth classifying and quantifying. To this end, we develop an information-theory based conceptual framework for measuring both the persistence of memory in microbes and the amount of information about the past encoded in history-dependent dynamics. This method produces a phenomenological measure of cellular memory without regard to the specific cellular mechanisms encoding it. We then apply this framework to a strain of Bacillus subtilis engineered to report on commitment to sporulation and degradative enzyme (AprE) synthesis and estimate the capacity of these systems and growth dynamics to 'remember' 10 distinct cell histories prior to application of a common stressor. The analysis suggests that B. subtilis remembers, both in short and long term, aspects of its cell history, and that this memory is distributed differently among the observables. While this study does not examine the mechanistic bases for memory, it presents a framework for quantifying memory in cellular behaviors and is thus a starting point for studying new questions about cellular regulation and evolutionary strategy.
