PhiSiGns: an online tool to identify signature genes in phages and design PCR primers for examining phage diversity

Phages from Ganges River curtail in vitro biofilms and planktonic growth of drug resistant Klebsiella pneumoniae in a zebrafish infection model

Bacteriophages are a promising alternative for curtailing infections caused by multi drug resistant (MDR) bacteria. The objective of the present study is to evaluate phage populations from water bodies to inhibit planktonic and biofilm mode of growth of drug resistant Klebsiella pneumoniae in vitro and curtail planktonic growth in vivo in a zebrafish model. Phage specific to K. pneumoniae (MTCC 432) was isolated from Ganges River (designated as KpG). One-step growth curve, in vitro time kill curve study and in vivo infection model were performed to evaluate the ability of phage to curtail planktonic growth. Crystal violet assay and colony biofilm assay were performed to determine the action of phages on biofilms. KpG phages had a greater burst size, better bactericidal potential and enhanced inhibitory effect against biofilms formed at liquid air and solid air interfaces. In vitro time kill assay showed a 3 log decline and a 6 log decline in K. pneumoniae colony counts, when phages were administered individually and in combination with streptomycin, respectively. In vivo injection of KpG phages revealed that it did not pose any toxicity to zebrafish as evidenced by liver/brain enzyme profiles and by histopathological analysis. The muscle tissue of zebrafish, infected with K. pneumoniae and treated with KpG phages alone and in combination with streptomycin showed a significant 77.7% and 97.2% decline in CFU/ml, respectively, relative to untreated control. Our study reveals that KpG phages has the potential to curtail plantonic and biofilm mode of growth in higher animal models.

Bioinformatic Tools and Guidelines for the Design of Fluorescence In Situ Hybridization Probes

Fluorescence in situ hybridization (FISH) is a well-established technique that allows the detection of microorganisms in diverse types of samples (e.g., clinical, food, environmental samples, and biofilm communities). The FISH probe design is an essential step in this technique. For this, two strategies can be used, the manual form based on multiple sequence alignment to identify conserved regions and programs/software specifically developed for the selection of the sequence of the probe. Additionally, databases/software for the theoretical evaluation of the probes in terms of specificity, sensitivity, and thermodynamic parameters (melting temperature and Gibbs free energy change) are used. The purpose of this chapter is to describe the essential steps and guidelines for the design of FISH probes (e.g., DNA and Nucleic Acid Mimic (NAM) probes), and its theoretical evaluation through the application of diverse bioinformatic tools.

Development and evaluation of taxon-specific primers for the selected Caudovirales taxa

The phage taxonomy is primarily based on the morphology derived from Transmission Electron Microscopic (TEM) studies. TEM based characterization is authentic and accepted by scientific community. However, TEM based identification is expensive and time consuming. After the phage isolation, before analysis TEM, a DNA based rapid method could be introduced. The DNA based method could dramatically reduce the number of samples analyzed by TEM and thereby increase the speed and reduce the cost of identification. In the present work, four environmental phage isolates were identified based on TEM studies and genome size. The identification of these four phages was validated using DNA based method. The taxon-specific DNA markers were identified through multiple sequence alignments. The primers were designed at conserved genes (DNA polymerase or integrase) of 4 different phage taxa viz. family Ackermannviridae, genus Jerseyvirus, genus T4virus, and genus P22virus. These primers were evaluated using both in vitro and in silico approach for the amplification of the target taxons. Majority of the primer sets were found to amplify member species of the targeted taxa in vitro. In In silico analysis, six primer sets intended for identification of family Ackermannviridae showed positive amplification of ≥86.7% classified species. Further, the primers targeting the genus Jerseyvirus and T4virus showed the amplification of 53.8% and ≥84.6% species, respectively. The present work is a case study performed to explore the possibility of use of taxon-specific primers for identification and taxonomic studies of newly isolated phages to supplement the TEM.

MRPrimerW: a tool for rapid design of valid high-quality primers for multiple target qPCR experiments

Design of high-quality primers for multiple target sequences is essential for qPCR experiments, but is challenging due to the need to consider both homology tests on off-target sequences and the same stringent filtering constraints on the primers. Existing web servers for primer design have major drawbacks, including requiring the use of BLAST-like tools for homology tests, lack of support for ranking of primers, TaqMan probes and simultaneous design of primers against multiple targets. Due to the large-scale computational overhead, the few web servers supporting homology tests use heuristic approaches or perform homology tests within a limited scope. Here, we describe the MRPrimerW, which performs complete homology testing, supports batch design of primers for multi-target qPCR experiments, supports design of TaqMan probes and ranks the resulting primers to return the top-1 best primers to the user. To ensure high accuracy, we adopted the core algorithm of a previously reported MapReduce-based method, MRPrimer, but completely redesigned it to allow users to receive query results quickly in a web interface, without requiring a MapReduce cluster or a long computation. MRPrimerW provides primer design services and a complete set of 341 963 135 in silico validated primers covering 99% of human and mouse genes. Free access: http://MRPrimerW.com.

MG-RAST, a Metagenomics Service for Analysis of Microbial Community Structure and Function

Approaches in molecular biology, particularly those that deal with high-throughput sequencing of entire microbial communities (the field of metagenomics), are rapidly advancing our understanding of the composition and functional content of microbial communities involved in climate change, environmental pollution, human health, biotechnology, etc. Metagenomics provides researchers with the most complete picture of the taxonomic (i.e., what organisms are there) and functional (i.e., what are those organisms doing) composition of natively sampled microbial communities, making it possible to perform investigations that include organisms that were previously intractable to laboratory-controlled culturing; currently, these constitute the vast majority of all microbes on the planet. All organisms contained in environmental samples are sequenced in a culture-independent manner, most often with 16S ribosomal amplicon methods to investigate the taxonomic or whole-genome shotgun-based methods to investigate the functional content of sampled communities. Metagenomics allows researchers to characterize the community composition and functional content of microbial communities, but it cannot show which functional processes are active; however, near parallel developments in transcriptomics promise a dramatic increase in our knowledge in this area as well. Since 2008, MG-RAST (Meyer et al., BMC Bioinformatics 9:386, 2008) has served as a public resource for annotation and analysis of metagenomic sequence data, providing a repository that currently houses more than 150,000 data sets (containing 60+ tera-base-pairs) with more than 23,000 publically available. MG-RAST, or the metagenomics RAST (rapid annotation using subsystems technology) server makes it possible for users to upload raw metagenomic sequence data in (preferably) fastq or fasta format. Assessments of sequence quality, annotation with respect to multiple reference databases, are performed automatically with minimal input from the user (see Subheading 4 at the end of this chapter for more details). Post-annotation analysis and visualization are also possible, directly through the web interface, or with tools like matR (metagenomic analysis tools for R, covered later in this chapter) that utilize the MG-RAST API ( http://api.metagenomics.anl.gov/api.html ) to easily download data from any stage in the MG-RAST processing pipeline. Over the years, MG-RAST has undergone substantial revisions to keep pace with the dramatic growth in the number, size, and types of sequence data that accompany constantly evolving developments in metagenomics and related -omic sciences (e.g., metatranscriptomics).

Computational approaches to predict bacteriophage-host relationships

Metagenomics has changed the face of virus discovery by enabling the accurate identification of viral genome sequences without requiring isolation of the viruses. As a result, metagenomic virus discovery leaves the first and most fundamental question about any novel virus unanswered: What host does the virus infect? The diversity of the global virosphere and the volumes of data obtained in metagenomic sequencing projects demand computational tools for virus–host prediction. We focus on bacteriophages (phages, viruses that infect bacteria), the most abundant and diverse group of viruses found in environmental metagenomes. By analyzing 820 phages with annotated hosts, we review and assess the predictive power of in silico phage–host signals. Sequence homology approaches are the most effective at identifying known phage–host pairs. Compositional and abundance-based methods contain significant signal for phage–host classification, providing opportunities for analyzing the unknowns in viral metagenomes. Together, these computational approaches further our knowledge of the interactions between phages and their hosts. Importantly, we find that all reviewed signals significantly link phages to their hosts, illustrating how current knowledge and insights about the interaction mechanisms and ecology of coevolving phages and bacteria can be exploited to predict phage–host relationships, with potential relevance for medical and industrial applications.

New viruses infecting bacteria are increasingly being discovered in many environments through sequence-based explorations. To understand their role in microbial ecosystems, computational tools are indispensable to prioritize and guide experimental efforts. This review assesses and discusses a range of bioinformatic approaches to predict bacteriophage–host relationships when all that is known is their genome sequence.

Using signature genes as tools to assess environmental viral ecology and diversity

Viruses (including bacteriophages) are the most abundant biological entities on the planet. As such, they are thought to have a major impact on all aspects of microbial community structure and function. Despite this critical role in ecosystem processes, the study of virus/phage diversity has lagged far behind parallel studies of the bacterial and eukaryotic kingdoms, largely due to the absence of any universal phylogenetic marker. Here we review the development and use of signature genes to investigate viral diversity, as a viable strategy for data sets of specific virus groups. Genes that have been used include those encoding structural proteins, such as portal protein, major capsid protein, and tail sheath protein, auxiliary metabolism genes, such as psbA, psbB,and phoH, and several polymerase genes. These marker genes have been used in combination with PCR-based fingerprinting and/or sequencing strategies to investigate spatial, temporal, and seasonal variations and diversity in a wide range of habitats.

MRPrimer: a MapReduce-based method for the thorough design of valid and ranked primers for PCR

Primer design is a fundamental technique that is widely used for polymerase chain reaction (PCR). Although many methods have been proposed for primer design, they require a great deal of manual effort to generate feasible and valid primers, including homology tests on off-target sequences using BLAST-like tools. That approach is inconvenient for many target sequences of quantitative PCR (qPCR) due to considering the same stringent and allele-invariant constraints. To address this issue, we propose an entirely new method called MRPrimer that can design all feasible and valid primer pairs existing in a DNA database at once, while simultaneously checking a multitude of filtering constraints and validating primer specificity. Furthermore, MRPrimer suggests the best primer pair for each target sequence, based on a ranking method. Through qPCR analysis using 343 primer pairs and the corresponding sequencing and comparative analyses, we showed that the primer pairs designed by MRPrimer are very stable and effective for qPCR. In addition, MRPrimer is computationally efficient and scalable and therefore useful for quickly constructing an entire collection of feasible and valid primers for frequently updated databases like RefSeq. Furthermore, we suggest that MRPrimer can be utilized conveniently for experiments requiring primer design, especially real-time qPCR.

Deep sequencing of the viral phoH gene reveals temporal variation, depth-specific composition, and persistent dominance of the same viral phoH genes in the Sargasso Sea

Deep sequencing of the viral phoH gene, a host-derived auxiliary metabolic gene, was used to track viral diversity throughout the water column at the Bermuda Atlantic Time-series Study (BATS) site in the summer (September) and winter (March) of three years. Viral phoH sequences reveal differences in the viral communities throughout a depth profile and between seasons in the same year. Variation was also detected between the same seasons in subsequent years, though these differences were not as great as the summer/winter distinctions. Over 3,600 phoH operational taxonomic units (OTUs; 97% sequence identity) were identified. Despite high richness, most phoH sequences belong to a few large, common OTUs whereas the majority of the OTUs are small and rare. While many OTUs make sporadic appearances at just a few times or depths, a small number of OTUs dominate the community throughout the seasons, depths, and years.

Diversity of environmental single-stranded DNA phages revealed by PCR amplification of the partial major capsid protein

The small single-stranded DNA (ssDNA) bacteriophages of the subfamily Gokushovirinae were traditionally perceived as narrowly targeted, niche-specific viruses infecting obligate parasitic bacteria, such as Chlamydia. The advent of metagenomics revealed gokushoviruses to be widespread in global environmental samples. This study expands knowledge of gokushovirus diversity in the environment by developing a degenerate PCR assay to amplify a portion of the major capsid protein (MCP) gene of gokushoviruses. Over 500 amplicons were sequenced from 10 environmental samples (sediments, sewage, seawater and freshwater), revealing the ubiquity and high diversity of this understudied phage group. Residue-level conservation data generated from multiple alignments was combined with a predicted 3D structure, revealing a tendency for structurally internal residues to be more highly conserved than surface-presenting protein-protein or viral-host interaction domains. Aggregating this data set into a phylogenetic framework, many gokushovirus MCP clades contained samples from multiple environments, although distinct clades dominated the different samples. Antarctic sediment samples contained the most diverse gokushovirus communities, whereas freshwater springs from Florida were the least diverse. Whether the observed diversity is being driven by environmental factors or host-binding interactions remains an open question. The high environmental diversity of this previously overlooked ssDNA viral group necessitates further research elucidating their natural hosts and exploring their ecological roles.

Viral communities associated with human pericardial fluids in idiopathic pericarditis

Pericarditis is a common human disease defined by inflammation of the pericardium. Currently, 40% to 85% of pericarditis cases have no identified etiology. Most of these cases are thought to be caused by an infection of undetected, unsuspected or unknown viruses. In this work, we used a culture- and sequence-independent approach to investigate the viral DNA communities present in human pericardial fluids. Seven viral metagenomes were generated from the pericardial fluid of patients affected by pericarditis of unknown etiology and one metagenome was generated from the pericardial fluid of a sudden infant death case. As a positive control we generated one metagenome from the pericardial fluid of a patient affected by pericarditis caused by herpesvirus type 3. Furthermore, we used as negative controls a total of 6 pericardial fluids from 6 different individuals affected by pericarditis of non-infectious origin: 5 of them were sequenced as a unique pool and the remaining one was sequenced separately. The results showed a significant presence of torque teno viruses especially in one patient, while herpesviruses and papillomaviruses were present in the positive control. Co-infections by different genotypes of the same viral type (torque teno viruses) or different viruses (herpesviruses and papillomaviruses) were observed. Sequences related to bacteriophages infecting Staphylococcus, Enterobacteria, Streptococcus, Burkholderia and Pseudomonas were also detected in three patients. This study detected torque teno viruses and papillomaviruses, for the first time, in human pericardial fluids.

Complete genome sequence of the naphthalene-degrading bacterium Pseudomonas stutzeri AN10 (CCUG 29243)

Pseudomonas stutzeri AN10 (CCUG 29243) can be considered a model strain for aerobic naphthalene degradation. We report the complete genome sequence of this bacterium. Its 4.71-Mb chromosome provides insights into other biodegradative capabilities of strain AN10 (i.e., benzoate catabolism) and suggests a high number of horizontal gene transfer events.

Gemi: PCR primers prediction from multiple alignments

Designing primers and probes for polymerase chain reaction (PCR) is a preliminary and critical step that requires the identification of highly conserved regions in a given set of sequences. This task can be challenging if the targeted sequences display a high level of diversity, as frequently encountered in microbiologic studies. We developed Gemi, an automated, fast, and easy-to-use bioinformatics tool with a user-friendly interface to design primers and probes based on multiple aligned sequences. This tool can be used for the purpose of real-time and conventional PCR and can deal efficiently with large sets of sequences of a large size.

De-MetaST-BLAST: a tool for the validation of degenerate primer sets and data mining of publicly available metagenomes

Development and use of primer sets to amplify nucleic acid sequences of interest is fundamental to studies spanning many life science disciplines. As such, the validation of primer sets is essential. Several computer programs have been created to aid in the initial selection of primer sequences that may or may not require multiple nucleotide combinations (i.e., degeneracies). Conversely, validation of primer specificity has remained largely unchanged for several decades, and there are currently few available programs that allows for an evaluation of primers containing degenerate nucleotide bases. To alleviate this gap, we developed the program De-MetaST that performs an in silico amplification using user defined nucleotide sequence dataset(s) and primer sequences that may contain degenerate bases. The program returns an output file that contains the in silico amplicons. When De-MetaST is paired with NCBI’s BLAST (De-MetaST-BLAST), the program also returns the top 10 nr NCBI database hits for each recovered in silico amplicon. While the original motivation for development of this search tool was degenerate primer validation using the wealth of nucleotide sequences available in environmental metagenome and metatranscriptome databases, this search tool has potential utility in many data mining applications.