Rapid 16S rRNA next-generation sequencing of polymicrobial clinical samples for diagnosis of complex bacterial infections.

Classifying individual bacterial species comprising complex, polymicrobial patient specimens remains a challenge for culture-based and molecular microbiology techniques in common clinical use. We therefore adapted practices from metagenomics research to rapidly catalog the bacterial composition of clinical specimens directly from patients, without need for prior culture. We have combined a semiconductor deep sequencing protocol that produces reads spanning 16S ribosomal RNA gene variable regions 1 and 2 (∼360 bp) with a de-noising pipeline that significantly improves the fraction of error-free sequences. The resulting sequences can be used to perform accurate genus- or species-level taxonomic assignment. We explore the microbial composition of challenging, heterogeneous clinical specimens by deep sequencing, culture-based strain typing, and Sanger sequencing of bulk PCR product. We report that deep sequencing can catalog bacterial species in mixed specimens from which usable data cannot be obtained by conventional clinical methods. Deep sequencing a collection of sputum samples from cystic fibrosis (CF) patients reveals well-described CF pathogens in specimens where they were not detected by standard clinical culture methods, especially for low-prevalence or fastidious bacteria. We also found that sputa submitted for CF diagnostic workup can be divided into a limited number of groups based on the phylogenetic composition of the airway microbiota, suggesting that metagenomic profiling may prove useful as a clinical diagnostic strategy in the future. The described method is sufficiently rapid (theoretically compatible with same-day turnaround times) and inexpensive for routine clinical use.

Bactobolin resistance is conferred by mutations in the L2 ribosomal protein.

UNLABELLED: Burkholderia thailandensis produces a family of polyketide-peptide molecules called bactobolins, some of which are potent antibiotics. We found that growth of B. thailandensis at 30°C versus that at 37°C resulted in increased production of bactobolins. We purified the three most abundant bactobolins and determined their activities against a battery of bacteria and mouse fibroblasts. Two of the three compounds showed strong activities against both bacteria and fibroblasts. The third analog was much less potent in both assays. These results suggested that the target of bactobolins might be conserved across bacteria and mammalian cells. To learn about the mechanism of bactobolin activity, we isolated four spontaneous bactobolin-resistant Bacillus subtilis mutants. We used genomic sequencing technology to show that each of the four resistant variants had mutations in rplB, which codes for the 50S ribosome-associated L2 protein. Ectopic expression of a mutant rplB gene in wild-type B. subtilis conferred bactobolin resistance. Finally, the L2 mutations did not confer resistance to other antibiotics known to interfere with ribosome function. Our data indicate that bactobolins target the L2 protein or a nearby site and that this is not the target of other antibiotics. We presume that the mammalian target of bactobolins involves the eukaryotic homolog of L2 (L8e). IMPORTANCE: Currently available antibiotics target surprisingly few cellular functions, and there is a need to identify novel antibiotic targets. We have been interested in the Burkholderia thailandensis bactobolins, and we sought to learn about the target of bactobolin activity by mapping spontaneous resistance mutations in the bactobolin-sensitive Bacillus subtilis. Our results indicate that the bactobolin target is the 50S ribosome-associated L2 protein or a region of the ribosome affected by L2. Bactobolin-resistant mutants are not resistant to other known ribosome inhibitors. Our evidence indicates that bactobolins interact with a novel antibiotic target.

Evolution of Burkholderia pseudomallei in recurrent melioidosis.

Burkholderia pseudomallei, the etiologic agent of human melioidosis, is capable of causing severe acute infection with overwhelming septicemia leading to death. A high rate of recurrent disease occurs in adult patients, most often due to recrudescence of the initial infecting strain. Pathogen persistence and evolution during such relapsing infections are not well understood. Bacterial cells present in the primary inoculum and in late infections may differ greatly, as has been observed in chronic disease, or they may be genetically similar. To test these alternative models, we conducted whole-genome comparisons of clonal primary and relapse B. pseudomallei isolates recovered six months to six years apart from four adult Thai patients. We found differences within each of the four pairs, and some, including a 330 Kb deletion, affected substantial portions of the genome. Many of the changes were associated with increased antibiotic resistance. We also found evidence of positive selection for deleterious mutations in a TetR family transcriptional regulator from a set of 107 additional B. pseudomallei strains. As part of the study, we sequenced to base-pair accuracy the genome of B. pseudomallei strain 1026b, the model used for genetic studies of B. pseudomallei pathogenesis and antibiotic resistance. Our findings provide new insights into pathogen evolution during long-term infections and have important implications for the development of intervention strategies to combat recurrent melioidosis.

Evidence for diversifying selection at the pyoverdine locus of Pseudomonas aeruginosa.

Pyoverdine is the primary siderophore of the gram-negative bacterium Pseudomonas aeruginosa. The pyoverdine region was recently identified as the most divergent locus alignable between strains in the P. aeruginosa genome. Here we report the nucleotide sequence and analysis of more than 50 kb in the pyoverdine region from nine strains of P. aeruginosa. There are three divergent sequence types in the pyoverdine region, which correspond to the three structural types of pyoverdine. The pyoverdine outer membrane receptor fpvA may be driving diversity at the locus: it is the most divergent alignable gene in the region, is the only gene that showed substantial intratype variation that did not appear to be generated by recombination, and shows evidence of positive selection. The hypothetical membrane protein PA2403 also shows evidence of positive selection; residues on one side of the membrane after protein folding are under positive selection. R', previously identified as a type IV strain, is clearly derived from a type III strain via a 3.4-kb deletion which removes one amino acid from the pyoverdine side chain peptide. This deletion represents a natural modification of the product of a nonribosomal peptide synthetase enzyme, whose consequences are predictive from the DNA sequence. There is also linkage disequilibrium between the pyoverdine region and pvdY, a pyoverdine gene separated by 30 kb from the pyoverdine region. The pyoverdine region shows evidence of horizontal transfer; we propose that some alleles in the region were introduced from other soil bacteria and have been subsequently maintained by diversifying selection.

Whole-genome sequence variation among multiple isolates of Pseudomonas aeruginosa.

Whole-genome shotgun sequencing was used to study the sequence variation of three Pseudomonas aeruginosa isolates, two from clonal infections of cystic fibrosis patients and one from an aquatic environment, relative to the genomic sequence of reference strain PAO1. The majority of the PAO1 genome is represented in these strains; however, at least three prominent islands of PAO1-specific sequence are apparent. Conversely, approximately 10% of the sequencing reads derived from each isolate fail to align with the PAO1 backbone. While average sequence variation among all strains is roughly 0.5%, regions of pronounced differences were evident in whole-genome scans of nucleotide diversity. We analyzed two such divergent loci, the pyoverdine and O-antigen biosynthesis regions, by complete resequencing. A thorough analysis of isolates collected over time from one of the cystic fibrosis patients revealed independent mutations resulting in the loss of O-antigen synthesis alternating with a mucoid phenotype. Overall, we conclude that most of the PAO1 genome represents a core P. aeruginosa backbone sequence while the strains addressed in this study possess additional genetic material that accounts for at least 10% of their genomes. Approximately half of these additional sequences are novel.

Genetic variation at the O-antigen biosynthetic locus in Pseudomonas aeruginosa.

The outer carbohydrate layer, or O antigen, of Pseudomonas aeruginosa varies markedly in different isolates of these bacteria, and at least 20 distinct O-antigen serotypes have been described. Previous studies have indicated that the major enzymes responsible for O-antigen synthesis are encoded in a cluster of genes that occupy a common genetic locus. We used targeted yeast recombinational cloning to isolate this locus from the 20 internationally recognized serotype strains. DNA sequencing of these isolated segments revealed that at least 11 highly divergent gene clusters occupy this region. Homology searches of the encoded protein products indicated that these gene clusters are likely to direct O-antigen biosynthesis. The O15 serotype strains lack functional gene clusters in the region analyzed, suggesting that O-antigen biosynthesis genes for this serotype are harbored in a different portion of the genome. The overall pattern underscores the plasticity of the P. aeruginosa genome, in which a specific site in a well-conserved genomic region can be occupied by any of numerous islands of functionally related DNA with diverse sequences.

Linker-mediated recombinational subcloning of large DNA fragments using yeast.

The homologous recombination pathway in yeast is an ideal tool for the sequence-specific assembly of plasmids. Complementary 80-nucleotide oligonucleotides that overlap a vector and a target fragment were found to serve as efficient recombination linkers for fragment subcloning. Using electroporation, single-stranded 80-mers were adequate for routine plasmid construction. A cycloheximide-based counterselection was introduced to increase the specificity of cloning by homologous recombination relative to nonspecific vector background. Reconstruction experiments suggest this counterselection increased cloning specificity by 100-fold. Cycloheximide counterselection was used in conjunction with 80-bp linkers to subclone targeted regions from bacterial artificial chromosomes. This technology may find broad application in the final stages of completing the Human Genome Sequencing Project and in applications of BAC clones to the functional analysis of complex genomes.