Environmental Vibrio parahaemolyticus DNA signatures validation

Insignia is a novel DNA computational system which uses highly efficient algorithms to compare bacterial genomes and to identify specific DNA signatures to distinguish a target bacterium, or group of bacteria, from all other known bacterial species. It is currently being validated using different bacterial groups, including Vibrio spp. In this study, the genomic analysis by Insignia was conducted on Vibrio parahaemolyticus, a halophilic gram-negative bacteria which constitutes a leading cause of seafood-borne disease. Insignia was used to identify 37 V. parahaemolyticus-specific signatures and to design PCR assays to validate the representative signature sequences by TaqMan essays. The 37 assays targeted loci distributed around the genome and detected genes coding for hypothetical proteins and for proteins involved in adhesion, starvation and virulence. A panel of V. parahaemolyticus environmental strains isolated from the North Adriatic Sea (Italy) and from the Black Sea (Georgia) was used to validate the selected signatures. The signature assays revealed both sensitive and specific and the method allowed a more accurate identification of the tested bacterial strains at the species level when compared to biochemical and PCR standard methods. Using Insignia, it was possible to distinguish two different groups among the strains previously identified as V. parahaemolyticus: most of the strains were included in a "V. parahaemolyticus-like group" showing nearly all of the signatures assayed while a small group of 10 strains contained only a few of the signatures tested. By sequencing the 16S rDNA of this latter group, it was confirmed that they were not V. parahaemolyticus but in fact belonged to other Vibrio species. No significant genome-wide differences were detected between the strains isolated in Italy and in Georgia though the very different geographical origin.     

Identification of Agrobacterium spp. present within Brassica napus seed by TaqMan PCR--implications for GM screening procedures

A fluorogenic probe (fliG-P), designed within a chromosomal DNA sequence, was used in a TaqMan PCR assay to identify Agrobacterium spp. The TaqMan assay detected 58 of 59 Agrobacterium strains tested, but did not detect 13 other Rhizobiaceae strains. Seedlings were grown from seven lots of surface-sterilised Brassica napus seed. Seedlings from these samples were placed in phosphate buffer and the resulting suspensions used to inoculate broth media selective for Agrobacterium biovars 1 and 2. Lysed broths (after 48 h incubation) were used as template in the fliG TaqMan PCR to detect Agrobacterium sp. in one of the seed samples. Individual Agrobacterium strains were isolated from this sample and tested by three Ti-plasmid conventional PCR assays. None of the strains possessed a plasmid. This is the first report of Agrobacterium sp. present within the seed of B. napus, a crop routinely screened for genetically modified DNA contamination using PCR assays with Agrobacterium sequences as targets.     

Bacterial population genomics and infectious disease diagnostics

New sequencing technologies have made the production of bacterial genome sequences increasingly easy, and it can be confidently forecasted that vast genomic databases will be generated in the next few years. Here, we detail how collections of bacterial genomes from a particular species (population genomics libraries) have already been used to improve the design of several diagnostic assays for bacterial pathogens. Genome sequencing itself is also becoming more commonly used for epidemiological, forensic and clinical investigations. There is an opportunity for the further development of bioinformatic tools to bring even further value to bacterial diagnostic genomics.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

Development and comparison of two assay formats for parallel detection of four biothreat pathogens by using suspension microarrays

Microarrays provide a powerful analytical tool for the simultaneous detection of multiple pathogens. We developed diagnostic suspension microarrays for sensitive and specific detection of the biothreat pathogens Bacillus anthracis, Yersinia pestis, Francisella tularensis and Coxiella burnetii. Two assay chemistries for amplification and labeling were developed, one method using direct hybridization and the other using target-specific primer extension, combined with hybridization to universal arrays. Asymmetric PCR products for both assay chemistries were produced by using a multiplex asymmetric PCR amplifying 16 DNA signatures (16-plex). The performances of both assay chemistries were compared and their advantages and disadvantages are discussed. The developed microarrays detected multiple signature sequences and an internal control which made it possible to confidently identify the targeted pathogens and assess their virulence potential. The microarrays were highly specific and detected various strains of the targeted pathogens. Detection limits for the different pathogen signatures were similar or slightly higher compared to real-time PCR. Probit analysis showed that even a few genomic copies could be detected with 95% confidence. The microarrays detected DNA from different pathogens mixed in different ratios and from spiked or naturally contaminated samples. The assays that were developed have a potential for application in surveillance and diagnostics.     

Detection and Identification of Brenneria nigrifluens, the Causal Agent of the Shallow Bark Canker of Walnut by, PCR Amplification

A 1 kb DNA band from strains of Brenneria nigrifluens, as shown by amplification of their genomic DNA by polymerase chain reaction (PCR) using minisatellite primer designed on the minisatellite sequence of the M13 phage, was isolated, cloned and sequenced. Specific oligonucleotides (F1–C3) were selected into this 1 kb DNA sequence and used in a PCR assay to detect and identify strains of B. nigrifluens . Several strains of B. nigrifluens were assessed with F1–C3 primers producing a specific band of approximately 250 bp pairs in length. This target was successfully amplified from purified genomic DNA, from bacterial culture and directly from infected walnut bark tissue. No amplification was obtained when the PCR assay was performed on other plant-pathogenic species from the following genera Brenneria, Erwinia, Agrobacterium, Pseudomonas, Ralstonia, Pectobacterium, Xanthomonas and from walnut-associated bacteria, indicating the specificity of these primers. The PCR assay with the primers described here provides a rapid, specific and sensitive diagnostic method for B. nigrifluens and a useful tool for epidemiological studies.     

Genomics and computational molecular biology

There has been a dramatic increase in the number of completely sequenced bacterial genomes during the past two years as a result of the efforts both of public genome agencies and the pharmaceutical industry. The availability of completely sequenced genomes permits more systematic analyses of genes, evolution and genome function than was otherwise possible. Using computational methods - which are used to identify genes and their functions including statistics, sequence similarity, motifs, profiles, protein folds and probabilistic models - it is possible to develop characteristic genome signatures, assign functions to genes, identify pathogenic genes, identify metabolic pathways, develop diagnostic probes and discover potential drug-binding sites. All of these directions are critical to understanding bacterial growth, pathogenicity and host-pathogen interactions.     

炭疽芽胞杆菌染色体特异序列的筛选及实时定量检测

筛选117条炭疽芽胞杆菌(Bacillusanthracis)特异序列,经双重特异性验证后得到19条理想的特异序列(genomicsignatures),其中6条符合设计TaqMan探针建立实时定量PCR的要求,根据常规PCR检测结果选择其中C04片段与炭疽芽胞杆菌毒性质粒pX01、pX02上的pagA、capB基因建立实时定量PCR检测体系。经试验证实这一体系检测灵敏度达到每PCR反应10~100个拷贝。利用12种相关菌株评价后获得100%特异性,对10份模拟污染标本和20份对照标本检测,所有污染标本均被检出,所有对照标本均为阴性。此方法特异、灵敏、高效,在炭疽芽胞杆菌感染的诊断和环境污染的检测等领域有潜在的应用前景。     

Distribution of repetitive DNA sequences in eubacteria and application to fingerprinting of bacterial genomes.

Dispersed repetitive DNA sequences have been described recently in eubacteria. To assess the distribution and evolutionary conservation of two distinct prokaryotic repetitive elements, consensus oligonucleotides were used in polymerase chain reaction [PCR] amplification and slot blot hybridization experiments with genomic DNA from diverse eubacterial species. Oligonucleotides matching Repetitive Extragenic Palindromic [REP] elements and Enterobacterial Repetitive Intergenic Consensus [ERIC] sequences were synthesized and tested as opposing PCR primers in the amplification of eubacterial genomic DNA. REP and ERIC consensus oligonucleotides produced clearly resolvable bands by agarose gel electrophoresis following PCR amplification. These band patterns provided unambiguous DNA fingerprints of different eubacterial species and strains. Both REP and ERIC probes hybridized preferentially to genomic DNA from Gram-negative enteric bacteria and related species. Widespread distribution of these repetitive DNA elements in the genomes of various microorganisms should enable rapid identification of bacterial species and strains, and be useful for the analysis of prokaryotic genomes.     

IDENTIFICATION OF NEW TARGET SEQUENCES FOR PCR DETECTION OF VIBRIO PARAHAEMOLYTICUS BY GENOME COMPARISON

A genome comparison method was used to identify specific target sequences for the polymerase chain reaction (PCR) detection of Vibrio parahaemolyticus, and the CDS value of this bacterium was compared with that of 139 other bacterial genomes. It was found that 20 CDS of V. parahaemolyticus were relatively specific according to their E value in BLAST (a new tool for comparing protein and nucleotide sequences), and four of them were selected for the design of PCR primers. There were positive amplification products of these four pairs of primers from nine V. parahaemolyticus strains, whereas there were no amplification products from nine other Vibrionaceae strains and four non -Vibrionaceae strains. An evaluation of detection sensitivities revealed that these four pairs of primers can be used in a PCR assay for the detection of V. parahaemolyticus.

PRACTICAL APPLICATIONS

An automatic BLAST method was developed in this study, by which species-specific sequences can be screened out rapidly. In this way, new and specific genes of Vibrio parahaemolyticus were identified to be used as target sequences for PCR detection. In terms of acceptable specificity and sensitivity, the four pairs of primers were selected by screening, which can be applied in PCR assays and other molecular methods. These kinds of methods might become commercial detection products in the new future. In addition, this method for searching specific DNA sequences can also be used for the mining specific sequences in other genus and species, such as Salmonella , Staphylococcus , etc.     

A high-throughput pipeline for the design of real-time PCR signatures

Background  

Pathogen diagnostic assays based on polymerase chain reaction (PCR) technology provide high sensitivity and specificity. However, the design of these diagnostic assays is computationally intensive, requiring high-throughput methods to identify unique PCR signatures in the presence of an ever increasing availability of sequenced genomes.     

Two-dimensional DNA displays for comparisons of bacterial genomes

We have developed two whole genome-scanning techniques to aid in the discovery of polymorphisms as well as horizontally acquired genes in prokaryotic organisms.First, two-dimensional bacterial genomic display (2DBGD) was developed using restriction enzyme fragmentation to separate genomic DNA based on size, and then employing denaturing gradient gel electrophoresis (DGGE) in the second dimension to exploit differences in sequence composition. This technique was used to generate high-resolution displays that enable the direct comparison of >800 genomic fragments simultaneously and can be adapted for the high-throughput comparison of bacterial genomes. 2DBGDs are capable of detecting acquired and altered DNA, however, only in very closely related strains. If used to compare more distantly related strains (e.g. different species within a genus) numerous small changes (i.e. small deletions and point mutations) unrelated to the interesting phenotype, would encumber the comparison of 2DBGDs. For this reason asecond method, bacterial comparative genomic hybridization (BCGH), was developed to directly compare bacterial genomes to identify gain or loss of genomic DNA. BCGH relies on performing 2DBGD on a pooled sample of genomic DNA from 2 strains to be compared and subsequently hybridizing the resulting 2DBGD blot separately with DNA from each individual strain. Unique spots (hybridization signals) represent foreign DNA. The identification of novel DNA is easily achieved by excising the DNA from a dried gel followed by subsequent cloning and sequencing. 2DBGD and BCGH thus represent novel high resolution genome scanning techniques for directly identifying altered and/or acquired DNA. Published: June 15, 2003     

MRD: a microsatellite repeats database for prokaryotic and eukaryotic genomes

MRD is a database system to access the microsatellite repeats information of genomes such as archea, eubacteria, and other eukaryotic genomes whose sequence information is available in public domains. MRD stores information about simple tandemly repeated k-mer sequences where k= 1 to 6, i.e. monomer to hexamer. The web interface allows the users to search for the repeat of their interest and to know about the association of the repeat with genes and genomic regions in the specific organism. The data contains the abundance and distribution of microsatellites in the coding and non-coding regions of the genome. The exact location of repeats with respect to genomic regions of interest (such as UTR, exon, intron or intergenic regions) whichever is applicable to organism is highlighted. MRD is available on the World Wide Web at and/or . The database is designed as an open-ended system to accommodate the microsatellite repeats information of other genomes whose complete sequences will be available in future through public domain.     

Mining virulence genes using metagenomics

When a bacterial genome is compared to the metagenome of an environment it inhabits, most genes recruit at high sequence identity. In free-living bacteria (for instance marine bacteria compared against the ocean metagenome) certain genomic regions are totally absent in recruitment plots, representing therefore genes unique to individual bacterial isolates. We show that these Metagenomic Islands (MIs) are also visible in bacteria living in human hosts when their genomes are compared to sequences from the human microbiome, despite the compartmentalized structure of human-related environments such as the gut. From an applied point of view, MIs of human pathogens (e.g. those identified in enterohaemorragic Escherichia coli against the gut metagenome or in pathogenic Neisseria meningitidis against the oral metagenome) include virulence genes that appear to be absent in related strains or species present in the microbiome of healthy individuals. We propose that this strategy (i.e. recruitment analysis of pathogenic bacteria against the metagenome of healthy subjects) can be used to detect pathogenicity regions in species where the genes involved in virulence are poorly characterized. Using this approach, we detect well-known pathogenicity islands and identify new potential virulence genes in several human pathogens.     

Identification of bacterial DNA markers for the detection of human fecal pollution in water

We used genome fragment enrichment and bioinformatics to identify several microbial DNA sequences with high potential for use as markers in PCR assays for detection of human fecal contamination in water. Following competitive solution-phase hybridization of total DNA from human and pig fecal samples, 351 plasmid clones were sequenced and were determined to define 289 different genomic DNA regions. These putative human-specific fecal bacterial DNA sequences were then analyzed by dot blot hybridization, which confirmed that 98% were present in the source human fecal microbial community and absent from the original pig fecal DNA extract. Comparative sequence analyses of these sequences suggested that a large number (43.5%) were predicted to encode bacterial secreted or surface-associated proteins. Deoxyoligonucleotide primers capable of annealing to a subset of 26 of the candidate sequences predicted to encode factors involved in interactions with host cells were then used in the PCR and did not amplify markers in DNA from any additional pig fecal specimens. These 26 PCR assays exhibited a range of specificity in tests with 11 other animal sources, with more than half amplifying markers only in specimens from dogs or cats. Four assays were more specific, detecting markers only in specimens from humans, including those from 18 different human populations examined. We then demonstrated the potential utility of these assays by using them to detect human fecal contamination in several impacted watersheds.     

OrthologID: automation of genome-scale ortholog identification within a parsimony framework

MOTIVATION: The determination of gene orthology is a prerequisite for mining and utilizing the rapidly increasing amount of sequence data for genome-scale phylogenetics and comparative genomic studies. Until now, most researchers use pairwise distance comparisons algorithms, such as BLAST, COG, RBH, RSD and INPARANOID, to determine gene orthology. In contrast, orthology determination within a character-based phylogenetic framework has not been utilized on a genomic scale owing to the lack of efficiency and automation. RESULTS: We have developed OrthologID, a Web application that automates the labor-intensive procedures of gene orthology determination within a character-based phylogenetic framework, thus making character-based orthology determination on a genomic scale possible. In addition to generating gene family trees and determining orthologous gene sets for complete genomes, OrthologID can also identify diagnostic characters that define each orthologous gene set, as well as diagnostic characters that are responsible for classifying query sequences from other genomes into specific orthology groups. The OrthologID database currently includes several complete plant genomes, including Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, as well as a unicellular outgroup, Chlamydomonas reinhardtii. To improve the general utility of OrthologID beyond plant species, we plan to expand our sequence database to include the fully sequenced genomes of prokaryotes and other non-plant eukaryotes. AVAILABILITY: http://nypg.bio.nyu.edu/orthologid/     

Heteropolymeric triplex-based genomic assay to detect pathogens or single-nucleotide polymorphisms in human genomic samples

Human genomic samples are complex and are considered difficult to assay directly without denaturation or PCR amplification. We report the use of a base-specific heteropolymeric triplex, formed by native duplex genomic target and an oligonucleotide third strand probe, to assay for low copy pathogen genomes present in a sample also containing human genomic duplex DNA, or to assay human genomic duplex DNA for Single Nucleotide Polymorphisms (SNP), without PCR amplification. Wild-type and mutant probes are used to identify triplexes containing FVL G1691A, MTHFR C677T and CFTR mutations. The specific triplex structure forms rapidly at room temperature in solution and may be detected without a separation step. YOYO-1, a fluorescent bis-intercalator, promotes and signals the formation of the specific triplex. Genomic duplexes may be assayed homogeneously with single base pair resolution. The specific triple-stranded structures of the assay may approximate homologous recombination intermediates, which various models suggest may form in either the major or minor groove of the duplex. The bases of the stable duplex target are rendered specifically reactive to the bases of the probe because of the activity of intercalated YOYO-1, which is known to decondense duplex locally 1.3 fold. This may approximate the local decondensation effected by recombination proteins such as RecA in vivo. Our assay, while involving triplex formation, is sui generis, as it is not homopurine sequence-dependent, as are "canonical triplexes". Rather, the base pair-specific heteropolymeric triplex of the assay is conformation-dependent. The highly sensitive diagnostic assay we present allows for the direct detection of base sequence in genomic duplex samples, including those containing human genomic duplex DNA, thereby bypassing the inherent problems and cost associated with conventional PCR based diagnostic assays.     

Avirulence gene and insertion element-based RFLP as well as RAPD markers reveal high levels of genomic polymorphism in the rice pathogen Xanthomonas oryzae pv. oryzae

Genetic polymorphism within the genomes of bacterial pathogens determines their evolutionary potential during long-term interaction with their hosts. To investigate the level of genetic variation in Xanthomonas oryzae pv. oryzae (Xoo), the causative agent of rice bacterial blight disease, three DNA marker systems, including (i) restriction fragment length polymorphism (RFLP) of the avrBs3/PthA family genes (avrXa27), (ii) RFLP of insertion (IS) elements and (iii) random amplified polymorphic DNA (RAPD) markers, were used to detect polymorphism among 32 Xoo strains that differed in their virulence patterns. All these strains contained multiple avrXa27 homologs that were variable in copy number and genomic location. RFLP of six IS elements revealed that these mobile sequences were abundant in Xoo genomes, with 150 of the total of 165 discernable markers being variable. Thirty-eight decamer primers of RAPD amplified a total of 691 bands, with 100% of them being variable. In addition, analysis of molecular variance (AMOVA) of data from RFLP analysis of IS elements and from RAPD analysis showed that most of the genetic variation residues were within Xoo populations, rather than between populations. Although all three DNA marker systems supported that substantial variation was maintained in Xoo genomes, Mantel tests did not identify significant correlation between the similarity coefficients calculated from them. The results of the present study indicated that Xoo genomes contain a high level of genetic polymorphism, which greatly facilitates the evolution of this important pathogen during interaction with its host rice plant.     

Quantification of Bacterial Uropathogens in Preclinical Samples Using Real-Time PCR Assays

Uropathogenic Escherichia coli (UPEC) and Staphylococcus saprophyticus (S. saprophyticus) are responsible for the majority of community-acquired urinary tract infections (UTI). Agar plating, a gold standard for detection of bacterial uropathogens, is labor intensive, limited for distinguishing between environmental contaminants and pathogens, and fails to effectively detect mixed infections. A reliable method for specific and sensitive quantitative assessment of infections would allow cost-effective evaluation of large numbers of experimental samples. A methodology such as quantitative PCR (qPCR) addresses the limitations of agar plating. We developed and validated highly specific and sensitive qPCR assays to assist researchers in the evaluation of potential vaccines and interventions in preclinical models of UPEC and S. saprophyticus UTI. The developed UPEC PCR targeted a highly conserved region of the UPEC hemolysin D (hlyD) gene that reproducibly detected type strains CFT073 and J96 over a 9 log range with high precision. To quantify S. saprophyticus genomes, a separate qPCR assay targeting the Trk transport gene was developed with an 8 log range. Neither assay detected bacterial species predicted to be sample contaminants. Using our optimized workflow that includes automated steps, up to 200 urine or tissue samples can be processed in as few as 3 h. Additionally, sequence comparisons of our primers and probe to other UTI bacterial strains indicated the broad applicability of these assays. These optimized qPCR assays provide a cost-effective and time-saving method for quantification of bacterial burdens in tissues and body fluids to assess the effectiveness of candidate vaccines or interventions.     

Identification of Bacterial DNA Markers for the Detection of Human Fecal Pollution in Water

We used genome fragment enrichment and bioinformatics to identify several microbial DNA sequences with high potential for use as markers in PCR assays for detection of human fecal contamination in water. Following competitive solution-phase hybridization of total DNA from human and pig fecal samples, 351 plasmid clones were sequenced and were determined to define 289 different genomic DNA regions. These putative human-specific fecal bacterial DNA sequences were then analyzed by dot blot hybridization, which confirmed that 98% were present in the source human fecal microbial community and absent from the original pig fecal DNA extract. Comparative sequence analyses of these sequences suggested that a large number (43.5%) were predicted to encode bacterial secreted or surface-associated proteins. Deoxyoligonucleotide primers capable of annealing to a subset of 26 of the candidate sequences predicted to encode factors involved in interactions with host cells were then used in the PCR and did not amplify markers in DNA from any additional pig fecal specimens. These 26 PCR assays exhibited a range of specificity in tests with 11 other animal sources, with more than half amplifying markers only in specimens from dogs or cats. Four assays were more specific, detecting markers only in specimens from humans, including those from 18 different human populations examined. We then demonstrated the potential utility of these assays by using them to detect human fecal contamination in several impacted watersheds.