Tom Slezak

MvirDB—a microbial database of protein toxins, virulence factors and antibiotic resistance genes for bio-defence applications

Knowledge of toxins, virulence factors and antibiotic resistance genes is essential for bio-defense applications aimed at identifying ‘functional’ signatures for characterizing emerging or engineered pathogens. Whereas genetic signatures identify a pathogen, functional signatures identify what a pathogen is capable of. To facilitate rapid identification of sequences and characterization of genes for signature discovery, we have collected all publicly available (as of this writing), organized sequences representing known toxins, virulence factors, and antibiotic resistance genes in one convenient database, which we believe will be of use to the bio-defense research community. MvirDB integrates DNA and protein sequence information from Tox-Prot, SCORPION, the PRINTS virulence factors, VFDB, TVFac, Islander, ARGO and a subset of VIDA. Entries in MvirDB are hyperlinked back to their original sources. A blast tool allows the user to blast against all DNA or protein sequences in MvirDB, and a browser tool allows the user to search the database to retrieve virulence factor descriptions, sequences, and classifications, and to download sequences of interest. MvirDB has an automated weekly update mechanism. Each protein sequence in MvirDB is annotated using our fully automated protein annotation system and is linked to that systems browser tool. MvirDB can be accessed at .

kSNP3.0: SNP detection and phylogenetic analysis of genomes without genome alignment or reference genome

UNLABELLED We announce the release of kSNP3.0, a program for SNP identification and phylogenetic analysis without genome alignment or the requirement for reference genomes. kSNP3.0 is a significantly improved version of kSNP v2. AVAILABILITY AND IMPLEMENTATION kSNP3.0 is implemented as a package of stand-alone executables for Linux and Mac OS X under the open-source BSD license. The executable packages, source code and a full User Guide are freely available at https://sourceforge.net/projects/ksnp/files/ CONTACT barryghall@gmail.com.

A microbial detection array (MDA) for viral and bacterial detection

BackgroundIdentifying the bacteria and viruses present in a complex sample is useful in disease diagnostics, product safety, environmental characterization, and research. Array-based methods have proven utility to detect in a single assay at a reasonable cost any microbe from the thousands that have been sequenced.MethodsWe designed a pan-Microbial Detection Array (MDA) to detect all known viruses (including phages), bacteria and plasmids and developed a novel statistical analysis method to identify mixtures of organisms from complex samples hybridized to the array. The array has broader coverage of bacterial and viral targets and is based on more recent sequence data and more probes per target than other microbial detection/discovery arrays in the literature. Family-specific probes were selected for all sequenced viral and bacterial complete genomes, segments, and plasmids. Probes were designed to tolerate some sequence variation to enable detection of divergent species with homology to sequenced organisms, and to have no significant matches to the human genome sequence.ResultsIn blinded testing on spiked samples with single or multiple viruses, the MDA was able to correctly identify species or strains. In clinical fecal, serum, and respiratory samples, the MDA was able to detect and characterize multiple viruses, phage, and bacteria in a sample to the family and species level, as confirmed by PCR.ConclusionsThe MDA can be used to identify the suite of viruses and bacteria present in complex samples.

Limitations of TaqMan PCR for Detecting Divergent Viral Pathogens Illustrated by Hepatitis A, B, C, and E Viruses and Human Immunodeficiency Virus

ABSTRACT Recent events illustrate the imperative to rapidly and accurately detect and identify pathogens during disease outbreaks, whether they are natural or engineered. Particularly for our primary goal of detecting bioterrorist releases, detection techniques must be both species-wide (capable of detecting all known strains of a given species) and species specific. Due to classification restrictions on the publication of data for species that may pose a bioterror threat, we illustrate the challenges of finding such assays using five nonthreat organisms that are nevertheless of public health concern: human immunodeficiency virus (HIV) and four species of hepatitis viruses. Fluorogenic probe-based PCR assays (TaqMan; Perkin-Elmer Corp., Applied Biosystems, Foster City, Calif.) may be sensitive, fast methods for the identification of species in which the genome is conserved among strains, such as hepatitis A virus. For species such as HIV, however, the strains are highly divergent. We use computational methods to show that nine TaqMan primer and probe sequences, or signatures, are needed to ensure that all strains will be detected, but this is an unfeasible number, considering the cost of TaqMan probes. Strains of hepatitis B, C, and E viruses show intermediate divergence, so that two to three TaqMan signatures are required to detect all strains of each virus. We conclude that for species such as hepatitis A virus with high levels of sequence conservation among strains, signatures can be found computationally for detection by the TaqMan assay, which is a sensitive, rapid, and cost-effective method. However, for species such as HIV with substantial genetic divergence among strains, the TaqMan assay becomes unfeasible and alternative detection methods may be required. We compare the TaqMan assay with some of the alternative nucleic acid-based detection techniques of microarray, chip, and bead technologies in terms of sensitivity, speed, and cost.

Cynomolgus macaque as an animal model for severe acute respiratory syndrome

Background The emergence of severe acute respiratory syndrome (SARS) in 2002 and 2003 affected global health and caused major economic disruption. Adequate animal models are required to study the underlying pathogenesis of SARS-associated coronavirus (SARS-CoV) infection and to develop effective vaccines and therapeutics. We report the first findings of measurable clinical disease in nonhuman primates (NHPs) infected with SARS-CoV. Methods and Findings In order to characterize clinically relevant parameters of SARS-CoV infection in NHPs, we infected cynomolgus macaques with SARS-CoV in three groups: Group I was infected in the nares and bronchus, group II in the nares and conjunctiva, and group III intravenously. Nonhuman primates in groups I and II developed mild to moderate symptomatic illness. All NHPs demonstrated evidence of viral replication and developed neutralizing antibodies. Chest radiographs from several animals in groups I and II revealed unifocal or multifocal pneumonia that peaked between days 8 and 10 postinfection. Clinical laboratory tests were not significantly changed. Overall, inoculation by a mucosal route produced more prominent disease than did intravenous inoculation. Half of the group I animals were infected with a recombinant infectious clone SARS-CoV derived from the SARS-CoV Urbani strain. This infectious clone produced disease indistinguishable from wild-type Urbani strain. Conclusions SARS-CoV infection of cynomolgus macaques did not reproduce the severe illness seen in the majority of adult human cases of SARS; however, our results suggest similarities to the milder syndrome of SARS-CoV infection characteristically seen in young children.

Multiniche crowding in genetic algorithms and its application to the assembly of dna restriction-fragments

The determination of the sequence of all nucleotide base-pairs in a DNA molecule, from restriction-fragment data, is a complex task and can be posed as the problem of finding the optima of a multimodal function. A genetic algorithm that uses multiniche crowding permits us to do this. Performance of this algorithm is first tested using a standard suite of test functions. The algorithm is next tested using two data sets obtained from the Human Genome Project at the Lawrence Livermore National Laboratory. The new method holds promise in automating the sequencing computations.

Cosmid assembly and anchoring to human chromosome 21

A human chromosome 21-specific cosmid library from the Lawrence Livermore National Laboratory has been analyzed by two complementary methods, fingerprinting and hybridization; 40% coverage of the entire chromosome 21 has been achieved. To prepare a contig pool, approximately 9300 cosmid clones randomly selected from the library were fingerprinted and automatically assembled into 467 overlapping sets by the fluorescence-tagged restriction fragment method. The average size of the overlapping sets was 9.5 cosmids with minimal tiling paths consisting of 5.4 cosmids with a 10-kb extension each. However, as many as 10% of overlaps within members were estimated to be false. For regional localization, we hybridized gridded arrays of cosmids with inter-Alu-PCR probes obtained from YAC clones and somatic cell hybrids and assigned 592 cosmids to 26 subregions of 21q. Of these, 371 clones were incorporated into 139 contigs, anchoring the total 1864 cosmids to the subregion. The remaining 221 clones were mapped as orphans. To correlate the cytogenetic, YAC, and cosmid maps on 21q, the translocation breakpoints of the chromosomes contained in the somatic cell hybrids were mapped with respect to the STS content of the YACs. From the gene cluster regions, 176 ribosomal and 25 alphoid clones were isolated by hybridization. Together, these sets of anchored contigs and cosmids will provide a valuable resource for construction of a high-resolution map and for isolation of genes of interest from chromosome 21.

AS2TS system for protein structure modeling and analysis

We present a set of programs and a website designed to facilitate protein structure comparison and protein structure modeling efforts. Our protein structure analysis and comparison services use the LGA (local-global alignment) program to search for regions of local similarity and to evaluate the level of structural similarity between compared protein structures. To facilitate the homology-based protein structure modeling process, our AL2TS service translates given sequence–structure alignment data into the standard Protein Data Bank (PDB) atom records (coordinates). For a given sequence of amino acids, the AS2TS (amino acid sequence to tertiary structure) system calculates (e.g. using PSI-BLAST PDB analysis) a list of the closest proteins from the PDB, and then a set of draft 3D models is automatically created. Web services are available at .

Conserved amino acid markers from past influenza pandemic strains

BackgroundFinding the amino acid mutations that affect the severity of influenza infections remains an open and challenging problem. Of special interest is better understanding how current circulating influenza strains could evolve into a new pandemic strain. Influenza proteomes from distinct viral phenotype classes were searched for class specific amino acid mutations conserved in past pandemics, using reverse engineered linear classifiers.ResultsThirty-four amino acid markers associated with host specificity and high mortality rate were found. Some markers had little impact on distinguishing the functional classes by themselves, however in combination with other mutations they improved class prediction. Pairwise combinations of influenza genomes were checked for reassortment and mutation events needed to acquire the pandemic conserved markers. Evolutionary pathways involving H1N1 human and swine strains mixed with avian strains show the potential to acquire the pandemic markers with a double reassortment and one or two amino acid mutations.ConclusionThe small mutation combinations found at multiple protein positions associated with viral phenotype indicate that surveillance tools could monitor genetic variation beyond single point mutations to track influenza strains. Finding that certain strain combinations have the potential to acquire pandemic conserved markers through a limited number of reassortment and mutation events illustrates the potential for reassortment and mutation events to lead to new circulating influenza strains.

A Functional Gene Array for Detection of Bacterial Virulence Elements

Emerging known and unknown pathogens create profound threats to public health. Platforms for rapid detection and characterization of microbial agents are critically needed to prevent and respond to disease outbreaks. Available detection technologies cannot provide broad functional information about known or novel organisms. As a step toward developing such a system, we have produced and tested a series of high-density functional gene arrays to detect elements of virulence and antibiotic resistance mechanisms. Our first generation array targets genes from Escherichia coli strains K12 and CFT073, Enterococcus faecalis and Staphylococcus aureus. We determined optimal probe design parameters for gene family detection and discrimination. When tested with organisms at varying phylogenetic distances from the four target strains, the array detected orthologs for the majority of targeted gene families present in bacteria belonging to the same taxonomic family. In combination with whole-genome amplification, the array detects femtogram concentrations of purified DNA, either spiked in to an aerosol sample background, or in combinations from one or more of the four target organisms. This is the first report of a high density NimbleGen microarray system targeting microbial antibiotic resistance and virulence mechanisms. By targeting virulence gene families as well as genes unique to specific biothreat agents, these arrays will provide important data about the pathogenic potential and drug resistance profiles of unknown organisms in environmental samples.