1
|
Harris EB, Ewool KKK, Bowden LC, Fierro J, Johnson D, Meinzer M, Tayler S, Grose JH. Genomic and Proteomic Analysis of Six Vi01-like Phages Reveals Wide Host Range and Multiple Tail Spike Proteins. Viruses 2024; 16:289. [PMID: 38400064 PMCID: PMC10892097 DOI: 10.3390/v16020289] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 02/05/2024] [Accepted: 02/07/2024] [Indexed: 02/25/2024] Open
Abstract
Enterobacteriaceae is a large family of Gram-negative bacteria composed of many pathogens, including Salmonella and Shigella. Here, we characterize six bacteriophages that infect Enterobacteriaceae, which were isolated from wastewater plants in the Wasatch front (Utah, United States). These phages are highly similar to the Kuttervirus vB_SenM_Vi01 (Vi01), which was isolated using wastewater from Kiel, Germany. The phages vary little in genome size and are between 157 kb and 164 kb, which is consistent with the sizes of other phages in the Vi01-like phage family. These six phages were characterized through genomic and proteomic comparison, mass spectrometry, and both laboratory and clinical host range studies. While their proteomes are largely unstudied, mass spectrometry analysis confirmed the production of five hypothetical proteins, several of which unveiled a potential operon that suggests a ferritin-mediated entry system on the Vi01-like phage family tail. However, no dependence on this pathway was observed for the single host tested herein. While unable to infect every genus of Enterobacteriaceae tested, these phages are extraordinarily broad ranged, with several demonstrating the ability to infect Salmonella enterica and Citrobacter freundii strains with generally high efficiency, as well as several clinical Salmonella enterica isolates, most likely due to their multiple tail fibers.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | - Julianne H. Grose
- Department of Microbiology and Molecular Biology, Brigham Young University, Provo, UT 84604, USA; (E.B.H.); (K.K.K.E.)
| |
Collapse
|
2
|
Assaf R, Xia F, Stevens R. Detecting operons in bacterial genomes via visual representation learning. Sci Rep 2021; 11:2124. [PMID: 33483546 PMCID: PMC7822928 DOI: 10.1038/s41598-021-81169-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2020] [Accepted: 12/30/2020] [Indexed: 12/05/2022] Open
Abstract
Contiguous genes in prokaryotes are often arranged into operons. Detecting operons plays a critical role in inferring gene functionality and regulatory networks. Human experts annotate operons by visually inspecting gene neighborhoods across pileups of related genomes. These visual representations capture the inter-genic distance, strand direction, gene size, functional relatedness, and gene neighborhood conservation, which are the most prominent operon features mentioned in the literature. By studying these features, an expert can then decide whether a genomic region is part of an operon. We propose a deep learning based method named Operon Hunter that uses visual representations of genomic fragments to make operon predictions. Using transfer learning and data augmentation techniques facilitates leveraging the powerful neural networks trained on image datasets by re-training them on a more limited dataset of extensively validated operons. Our method outperforms the previously reported state-of-the-art tools, especially when it comes to predicting full operons and their boundaries accurately. Furthermore, our approach makes it possible to visually identify the features influencing the network’s decisions to be subsequently cross-checked by human experts.
Collapse
Affiliation(s)
- Rida Assaf
- Department of Computer Science, University of Chicago, Chicago, 60637, USA.
| | - Fangfang Xia
- Computing Environment and Life Sciences Division, Argonne National Laboratory, Lemont, 60439, USA.,Data Science and Learning Division, Argonne National Laboratory, Lemont, 60439, USA
| | - Rick Stevens
- The University of Chicago Consortium for Advanced Science and Engineering, University of Chicago, Chicago, 60637, USA.,Computing Environment and Life Sciences Division, Argonne National Laboratory, Lemont, 60439, USA
| |
Collapse
|
3
|
Shmakov SA, Faure G, Makarova KS, Wolf YI, Severinov KV, Koonin EV. Systematic prediction of functionally linked genes in bacterial and archaeal genomes. Nat Protoc 2019; 14:3013-3031. [PMID: 31520072 DOI: 10.1038/s41596-019-0211-1] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Accepted: 06/13/2019] [Indexed: 11/09/2022]
Abstract
Functionally linked genes in bacterial and archaeal genomes are often organized into operons. However, the composition and architecture of operons are highly variable and frequently differ even among closely related genomes. Therefore, to efficiently extract reliable functional predictions for uncharacterized genes from comparative analyses of the rapidly growing genomic databases, dedicated computational approaches are required. We developed a protocol to systematically and automatically identify genes that are likely to be functionally associated with a 'bait' gene or locus by using relevance metrics. Given a set of bait loci and a genomic database defined by the user, this protocol compares the genomic neighborhoods of the baits to identify genes that are likely to be functionally linked to the baits by calculating the abundance of a given gene within and outside the bait neighborhoods and the distance to the bait. We exemplify the performance of the protocol with three test cases, namely, genes linked to CRISPR-Cas systems using the 'CRISPRicity' metric, genes associated with archaeal proviruses and genes linked to Argonaute genes in halobacteria. The protocol can be run by users with basic computational skills. The computational cost depends on the sizes of the genomic dataset and the list of reference loci and can vary from one CPU-hour to hundreds of hours on a supercomputer.
Collapse
Affiliation(s)
- Sergey A Shmakov
- National Center for Biotechnology Information, National Library of Medicine, Bethesda, MD, USA.,Skolkovo Institute of Science and Technology, Skolkovo, Russia
| | - Guilhem Faure
- National Center for Biotechnology Information, National Library of Medicine, Bethesda, MD, USA.,Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Kira S Makarova
- National Center for Biotechnology Information, National Library of Medicine, Bethesda, MD, USA
| | - Yuri I Wolf
- National Center for Biotechnology Information, National Library of Medicine, Bethesda, MD, USA
| | - Konstantin V Severinov
- Skolkovo Institute of Science and Technology, Skolkovo, Russia.,Waksman Institute of Microbiology, Rutgers, The State University of New Jersey, Piscataway, NJ, USA.,Institute of Molecular Genetics, Russian Academy of Sciences, Moscow, Russia
| | - Eugene V Koonin
- National Center for Biotechnology Information, National Library of Medicine, Bethesda, MD, USA.
| |
Collapse
|
4
|
Chuang LY, Yang CH, Tsai JH, Yang CH. Operon prediction using chaos embedded particle swarm optimization. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2013; 10:1299-1309. [PMID: 24384714 DOI: 10.1109/tcbb.2013.63] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Operons contain valuable information for drug design and determining protein functions. Genes within an operon are co-transcribed to a single-strand mRNA and must be coregulated. The identification of operons is, thus, critical for a detailed understanding of the gene regulations. However, currently used experimental methods for operon detection are generally difficult to implement and time consuming. In this paper, we propose a chaotic binary particle swarm optimization (CBPSO) to predict operons in bacterial genomes. The intergenic distance, participation in the same metabolic pathway and the cluster of orthologous groups (COG) properties of the Escherichia coli genome are used to design a fitness function. Furthermore, the Bacillus subtilis, Pseudomonas aeruginosa PA01, Staphylococcus aureus and Mycobacterium tuberculosis genomes are tested and evaluated for accuracy, sensitivity, and specificity. The computational results indicate that the proposed method works effectively in terms of enhancing the performance of the operon prediction. The proposed method also achieved a good balance between sensitivity and specificity when compared to methods from the literature.
Collapse
Affiliation(s)
| | - Cheng-Huei Yang
- National Kaohsiung Institute of Marine Technology, Kaohsiung
| | - Jui-Hung Tsai
- National Kaohsiung University of Applied Sciences, Kaohsiung
| | - Cheng-Hong Yang
- National Kaohsiung University of Applied Sciences, Kaohsiung
| |
Collapse
|
5
|
Abstract
An operon is a fundamental unit of transcription and contains specific functional genes for the construction and regulation of networks at the entire genome level. The correct prediction of operons is vital for understanding gene regulations and functions in newly sequenced genomes. As experimental methods for operon detection tend to be nontrivial and time consuming, various methods for operon prediction have been proposed in the literature. In this study, a binary particle swarm optimization is used for operon prediction in bacterial genomes. The intergenic distance, participation in the same metabolic pathway, the cluster of orthologous groups, the gene length ratio and the operon length are used to design a fitness function. We trained the proper values on the Escherichia coli genome, and used the above five properties to implement feature selection. Finally, our study used the intergenic distance, metabolic pathway and the gene length ratio property to predict operons. Experimental results show that the prediction accuracy of this method reached 92.1%, 93.3% and 95.9% on the Bacillus subtilis genome, the Pseudomonas aeruginosa PA01 genome and the Staphylococcus aureus genome, respectively. This method has enabled us to predict operons with high accuracy for these three genomes, for which only limited data on the properties of the operon structure exists.
Collapse
Affiliation(s)
- Li-Yeh Chuang
- Department of Chemical Engineering, I-Shou University, Kaohsiung, Taiwan
| | | | | |
Collapse
|
6
|
Taboada B, Verde C, Merino E. High accuracy operon prediction method based on STRING database scores. Nucleic Acids Res 2010; 38:e130. [PMID: 20385580 PMCID: PMC2896540 DOI: 10.1093/nar/gkq254] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We present a simple and highly accurate computational method for operon prediction, based on intergenic distances and functional relationships between the protein products of contiguous genes, as defined by STRING database (Jensen,L.J., Kuhn,M., Stark,M., Chaffron,S., Creevey,C., Muller,J., Doerks,T., Julien,P., Roth,A., Simonovic,M. et al. (2009) STRING 8–a global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res., 37, D412–D416). These two parameters were used to train a neural network on a subset of experimentally characterized Escherichia coli and Bacillus subtilis operons. Our predictive model was successfully tested on the set of experimentally defined operons in E. coli and B. subtilis, with accuracies of 94.6 and 93.3%, respectively. As far as we know, these are the highest accuracies ever obtained for predicting bacterial operons. Furthermore, in order to evaluate the predictable accuracy of our model when using an organism's data set for the training procedure, and a different organism's data set for testing, we repeated the E. coli operon prediction analysis using a neural network trained with B. subtilis data, and a B. subtilis analysis using a neural network trained with E. coli data. Even for these cases, the accuracies reached with our method were outstandingly high, 91.5 and 93%, respectively. These results show the potential use of our method for accurately predicting the operons of any other organism. Our operon predictions for fully-sequenced genomes are available at http://operons.ibt.unam.mx/OperonPredictor/.
Collapse
Affiliation(s)
- Blanca Taboada
- Centro de Ciencias Aplicadas y Desarrollo Tecnológico, Universidad Nacional Autónoma de México, México, D.F., México
| | | | | |
Collapse
|
7
|
Hamer R, Chen PY, Armitage JP, Reinert G, Deane CM. Deciphering chemotaxis pathways using cross species comparisons. BMC SYSTEMS BIOLOGY 2010; 4:3. [PMID: 20064255 PMCID: PMC2829493 DOI: 10.1186/1752-0509-4-3] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/10/2009] [Accepted: 01/11/2010] [Indexed: 12/29/2022]
Abstract
Background Chemotaxis is the process by which motile bacteria sense their chemical environment and move towards more favourable conditions. Escherichia coli utilises a single sensory pathway, but little is known about signalling pathways in species with more complex systems. Results To investigate whether chemotaxis pathways in other bacteria follow the E. coli paradigm, we analysed 206 species encoding at least 1 homologue of each of the 5 core chemotaxis proteins (CheA, CheB, CheR, CheW and CheY). 61 species encode more than one of all of these 5 proteins, suggesting they have multiple chemotaxis pathways. Operon information is not available for most bacteria, so we developed a novel statistical approach to cluster che genes into putative operons. Using operon-based models, we reconstructed putative chemotaxis pathways for all 206 species. We show that cheA-cheW and cheR-cheB have strong preferences to occur in the same operon as two-gene blocks, which may reflect a functional requirement for co-transcription. However, other che genes, most notably cheY, are more dispersed on the genome. Comparison of our operons with shuffled equivalents demonstrates that specific patterns of genomic location may be a determining factor for the observed in vivo chemotaxis pathways. We then examined the chemotaxis pathways of Rhodobacter sphaeroides. Here, the PpfA protein is known to be critical for correct partitioning of proteins in the cytoplasmically-localised pathway. We found ppfA in che operons of many species, suggesting that partitioning of cytoplasmic Che protein clusters is common. We also examined the apparently non-typical chemotaxis components, CheA3, CheA4 and CheY6. We found that though variants of CheA proteins are rare, the CheY6 variant may be a common type of CheY, with a significantly disordered C-terminal region which may be functionally significant. Conclusions We find that many bacterial species potentially have multiple chemotaxis pathways, with grouping of che genes into operons likely to be a major factor in keeping signalling pathways distinct. Gene order is highly conserved with cheA-cheW and cheR-cheB blocks, perhaps reflecting functional linkage. CheY behaves differently to other Che proteins, both in its genomic location and its putative protein interactions, which should be considered when modelling chemotaxis pathways.
Collapse
Affiliation(s)
- Rebecca Hamer
- Department of Statistics, University of Oxford, Oxford, UK
| | | | | | | | | |
Collapse
|
8
|
Proudfoot M, Sanders SA, Singer A, Zhang R, Brown G, Binkowski A, Xu L, Lukin JA, Murzin AG, Joachimiak A, Arrowsmith CH, Edwards AM, Savchenko AV, Yakunin AF. Biochemical and structural characterization of a novel family of cystathionine beta-synthase domain proteins fused to a Zn ribbon-like domain. J Mol Biol 2007; 375:301-15. [PMID: 18021800 DOI: 10.1016/j.jmb.2007.10.060] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2007] [Revised: 10/23/2007] [Accepted: 10/24/2007] [Indexed: 12/22/2022]
Abstract
We have identified a novel family of proteins, in which the N-terminal cystathionine beta-synthase (CBS) domain is fused to the C-terminal Zn ribbon domain. Four proteins were overexpressed in Escherichia coli and purified: TA0289 from Thermoplasma acidophilum, TV1335 from Thermoplasma volcanium, PF1953 from Pyrococcus furiosus, and PH0267 from Pyrococcus horikoshii. The purified proteins had a red/purple color in solution and an absorption spectrum typical of rubredoxins (Rds). Metal analysis of purified proteins revealed the presence of several metals, with iron and zinc being the most abundant metals (2-67% of iron and 12-74% of zinc). Crystal structures of both mercury- and iron-bound TA0289 (1.5-2.0 A resolution) revealed a dimeric protein whose intersubunit contacts are formed exclusively by the alpha-helices of two cystathionine beta-synthase subdomains, whereas the C-terminal domain has a classical Zn ribbon planar architecture. All proteins were reversibly reduced by chemical reductants (ascorbate or dithionite) or by the general Rd reductase NorW from E. coli in the presence of NADH. Reduced TA0289 was found to be capable of transferring electrons to cytochrome C from horse heart. Likewise, the purified Zn ribbon protein KTI11 from Saccharomyces cerevisiae had a purple color in solution and an Rd-like absorption spectrum, contained both iron and zinc, and was reduced by the Rd reductase NorW from E. coli. Thus, recombinant Zn ribbon domains from archaea and yeast demonstrate an Rd-like electron carrier activity in vitro. We suggest that, in vivo, some Zn ribbon domains might also bind iron and therefore possess an electron carrier activity, adding another physiological role to this large family of important proteins.
Collapse
Affiliation(s)
- Michael Proudfoot
- Banting and Best Department of Medical Research, University of Toronto, 112 College Street, Room 72, Toronto, ON, Canada
| | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
9
|
Harrington ED, Singh AH, Doerks T, Letunic I, von Mering C, Jensen LJ, Raes J, Bork P. Quantitative assessment of protein function prediction from metagenomics shotgun sequences. Proc Natl Acad Sci U S A 2007; 104:13913-8. [PMID: 17717083 PMCID: PMC1955820 DOI: 10.1073/pnas.0702636104] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
To assess the potential of protein function prediction in environmental genomics data, we analyzed shotgun sequences from four diverse and complex habitats. Using homology searches as well as customized gene neighborhood methods that incorporate intergenic and evolutionary distances, we inferred specific functions for 76% of the 1.4 million predicted ORFs in these samples (83% when nonspecific functions are considered). Surprisingly, these fractions are only slightly smaller than the corresponding ones in completely sequenced genomes (83% and 86%, respectively, by using the same methodology) and considerably higher than previously thought. For as many as 75,448 ORFs (5% of the total), only neighborhood methods can assign functions, illustrated here by a previously undescribed gene associated with the well characterized heme biosynthesis operon and a potential transcription factor that might regulate a coupling between fatty acid biosynthesis and degradation. Our results further suggest that, although functions can be inferred for most proteins on earth, many functions remain to be discovered in numerous small, rare protein families.
Collapse
Affiliation(s)
- E. D. Harrington
- *Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany; and
| | - A. H. Singh
- *Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany; and
| | - T. Doerks
- *Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany; and
| | - I. Letunic
- *Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany; and
| | - C. von Mering
- *Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany; and
| | - L. J. Jensen
- *Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany; and
| | - J. Raes
- *Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany; and
| | - P. Bork
- *Structural and Computational Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany; and
- Max Delbrück Centre for Molecular Medicine, D-13092 Berlin, Germany
- To whom correspondence should be addressed. E-mail:
| |
Collapse
|
10
|
Miller LD, Yost CK, Hynes MF, Alexandre G. The major chemotaxis gene cluster of Rhizobium leguminosarum bv. viciae is essential for competitive nodulation. Mol Microbiol 2007; 63:348-62. [PMID: 17163982 DOI: 10.1111/j.1365-2958.2006.05515.x] [Citation(s) in RCA: 84] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Rhizobium leguminosarum biovar viciae strain 3841 is a motile alpha-proteobacterium that can establish a nitrogen-fixing symbiosis within the roots of pea plants. In order to determine the contribution of chemotaxis to the lifestyle of R. leguminosarum, we have characterized the function of two chemotaxis gene clusters (che1 and che2) in controlling motility behaviour. We have found that both chemotaxis gene clusters modulate the motility swimming bias of R. leguminosarum cells and that the che1 cluster is the major pathway controlling swimming bias and chemotaxis. The che2 cluster also contributes to swimming bias, but has a minor effect on chemotaxis. Using competitive nodulation assays, we have demonstrated that a functional che1 cluster, but not the che2 cluster, promotes competitive nodulation of the peas. This finding implies that the environmental cue(s) triggering chemotaxis of R. leguminosarum bv. viciae cells towards the roots of pea and facilitating colonization are likely to be processed through the che1 cluster despite the contribution of both che clusters to swimming behaviour. A phylogenetic analysis of the distribution of che1 and che2 orthologues in the alpha-proteobacteria together with our results allow us to propose that che1 homologues are major controllers of chemotaxis and host association in the Rhizobiaceae.
Collapse
Affiliation(s)
- Lance D Miller
- Department of Biochemistry, Cellular and Molecular Biology and Department of Microbiology, The University of Tennessee, Knoxville, TN 37996, USA
| | | | | | | |
Collapse
|