Bacterial Species Determination from DNA-DNA Hybridization by Using Genome Fragments and DNA Microarrays

Whole genomic DNA-DNA hybridization has been a cornerstone of bacterial species determination but is not widely used because it is not easily implemented. We have developed a method based on random genome fragments and DNA microarray technology that overcomes the disadvantages of whole-genome DNA-DNA hybridization. Reference genomes of four fluorescent Pseudomonas species were fragmented, and 60 to 96 genome fragments of approximately 1 kb from each strain were spotted on microarrays. Genomes from 12 well-characterized fluorescent Pseudomonas strains were labeled with Cy dyes and hybridized to the arrays. Cluster analysis of the hybridization profiles revealed taxonomic relationships between bacterial strains tested at species to strain level resolution, suggesting that this approach is useful for the identification of bacteria as well as determining the genetic distance among bacteria. Since arrays can contain thousands of DNA spots, a single array has the potential for broad identification capacity. In addition, the method does not require laborious cross-hybridizations and can provide an open database of hybridization profiles, avoiding the limitations of traditional DNA-DNA hybridization.     

Bacterial species determination from DNA-DNA hybridization by using genome fragments and DNA microarrays

Whole genomic DNA-DNA hybridization has been a cornerstone of bacterial species determination but is not widely used because it is not easily implemented. We have developed a method based on random genome fragments and DNA microarray technology that overcomes the disadvantages of whole-genome DNA-DNA hybridization. Reference genomes of four fluorescent Pseudomonas species were fragmented, and 60 to 96 genome fragments of approximately 1 kb from each strain were spotted on microarrays. Genomes from 12 well-characterized fluorescent Pseudomonas strains were labeled with Cy dyes and hybridized to the arrays. Cluster analysis of the hybridization profiles revealed taxonomic relationships between bacterial strains tested at species to strain level resolution, suggesting that this approach is useful for the identification of bacteria as well as determining the genetic distance among bacteria. Since arrays can contain thousands of DNA spots, a single array has the potential for broad identification capacity. In addition, the method does not require laborious cross-hybridizations and can provide an open database of hybridization profiles, avoiding the limitations of traditional DNA-DNA hybridization.     

Development and evaluation of functional gene arrays for detection of selected genes in the environment.

To determine the potential of DNA array technology for assessing functional gene diversity and distribution, a prototype microarray was constructed with genes involved in nitrogen cycling: nitrite reductase (nirS and nirK) genes, ammonia mono-oxygenase (amoA) genes, and methane mono-oxygenase (pmoA) genes from pure cultures and those cloned from marine sediments. In experiments using glass slide microarrays, genes possessing less than 80 to 85% sequence identity were differentiated under hybridization conditions of high stringency (65 degrees C). The detection limit for nirS genes was approximately 1 ng of pure genomic DNA and 25 ng of soil community DNA using our optimized protocol. A linear quantitative relationship (r(2) = 0.89 to 0.94) was observed between signal intensity and target DNA concentration over a range of 1 to 100 ng for genomic DNA (or genomic DNA equivalent) from both pure cultures and mixed communities. However, the quantitative capacity of microarrays for measuring the relative abundance of targeted genes in complex environmental samples is less clear due to divergent target sequences. Sequence divergence and probe length affected hybridization signal intensity within a certain range of sequence identity and size, respectively. This prototype functional gene array did reveal differences in the apparent distribution of nir and amoA and pmoA gene families in sediment and soil samples. Our results indicate that glass-based microarray hybridization has potential as a tool for revealing functional gene composition in natural microbial communities; however, more work is needed to improve sensitivity and quantitation and to understand the associated issue of specificity.     

Specificity of randomly generated genomic DNA fragment probes on a DNA array

The use of randomly generated DNA fragment sequences as probes on DNA arrays offers a unique potential for exploring unsequenced microorganisms. In this study, the detection specificity was evaluated with respect to probe-target sequence similarity using genomic DNAs of four Pseudomonas strains. Genome fragments averaging 2000?bp were found to be specific enough to discriminate 85-90% similarity under highly stringent hybridization conditions. Such stringent conditions compromised signal intensities; however, specific signals remained detectable at the highest stringency (at 75?°C hybridization) with negligible false negatives. These results suggest that, without any probe design or selection, genomic fragments can provide a reasonable specificity for microbial diagnostics or species delineation by DNA-DNA similarities.     

Bordetella species are distinguished by patterns of substantial gene loss and host adaptation

Pathogens of the bacterial genus Bordetella cause respiratory disease in humans and animals. Although virulence and host specificity vary across the genus, the genetic determinants of this diversity remain unidentified. To identify genes that may underlie key phenotypic differences between these species and clarify their evolutionary relationships, we performed a comparative analysis of genome content in 42 Bordetella strains by hybridization of genomic DNA to a microarray representing the genomes of three Bordetella species and by subtractive hybridization. Here we show that B. pertussis and B. parapertussis are predominantly differentiated from B. bronchiseptica by large, species-specific regions of difference, many of which encode or direct synthesis of surface structures, including lipopolysaccharide O antigen, which may be important determinants of host specificity. The species also exhibit sequence diversity at a number of surface protein-encoding loci, including the fimbrial major subunit gene, fim2. Gene loss, rather than gene acquisition, accompanied by the proliferation of transposons, has played a fundamental role in the evolution of the pathogenic bordetellae and may represent a conserved evolutionary mechanism among other groups of microbial pathogens.     

Development and Evaluation of Functional Gene Arrays for Detection of Selected Genes in the Environment

To determine the potential of DNA array technology for assessing functional gene diversity and distribution, a prototype microarray was constructed with genes involved in nitrogen cycling: nitrite reductase (nirS and nirK) genes, ammonia mono-oxygenase (amoA) genes, and methane mono-oxygenase (pmoA) genes from pure cultures and those cloned from marine sediments. In experiments using glass slide microarrays, genes possessing less than 80 to 85% sequence identity were differentiated under hybridization conditions of high stringency (65°C). The detection limit for nirS genes was approximately 1 ng of pure genomic DNA and 25 ng of soil community DNA using our optimized protocol. A linear quantitative relationship (r² = 0.89 to 0.94) was observed between signal intensity and target DNA concentration over a range of 1 to 100 ng for genomic DNA (or genomic DNA equivalent) from both pure cultures and mixed communities. However, the quantitative capacity of microarrays for measuring the relative abundance of targeted genes in complex environmental samples is less clear due to divergent target sequences. Sequence divergence and probe length affected hybridization signal intensity within a certain range of sequence identity and size, respectively. This prototype functional gene array did reveal differences in the apparent distribution of nir and amoA and pmoA gene families in sediment and soil samples. Our results indicate that glass-based microarray hybridization has potential as a tool for revealing functional gene composition in natural microbial communities; however, more work is needed to improve sensitivity and quantitation and to understand the associated issue of specificity.     

中国吐鲁番两株产木聚糖酶的极端耐碱Bacillus halodurans的分类鉴定

通过生理生化实验、16S rDNA序列分析和同源性杂交,将分离到的XJU-1和XJU-80菌种进行了分类鉴定。XJU-1和XJU-80具有较宽的pH生长范围(分别是pH4·5 ~12·6和pH3·8 ~12·6) ,其G C mol%含量分别是40·5mol%和42·2mol%。16S rDNA序列分析和DNA-DNA同源杂交结果表明,XJU-1和XJU-80与Bacillushalodurans C-125和Bacillus halodurans DSM497T具有较高的同源性(99 %) ;两者之间也具有85 %的相关性,但其与Bacillus halodurans C-125和Bacillus halodurans DSM497T分别具有81·3 %和71·5 %的相关性。基于以上结果,将两株分离菌株分类为Bacillus halodurans的两个新品系。     

Genetic diversity of Rhodopirellula strains

Rhodopirellula baltica SH1^T is a marine planctomycete with 7,325 genes in its genome. Ten strains of the genus Rhodopirellula were studied in whole genome microarray experiments to assess the extent of their genetic relatedness to R. baltica SH1^T. DNA of strains which were previously affiliated with the species R. baltica (OTU A) hybridized with 3,645–5,728 genes of the type strain on the microarray. Strains SH398 and 6C (OTU B), representing a closely related species with an average nucleotide identity of 88 %, showed less hybridization signals: 1,816 and 3,302 genes gave a hybridization signal, respectively. Comparative genomics of eight permanent draft genomes revealed the presence of over 4,000 proteins common in R. baltica SH1^T and strains of OTU A or B. The genus Rhodopirellula is characterized by large genomes, with over 7,000 genes per genome and a core genome of around 3000 genes. Individual Rhodopirellula strains have a large portion of strain-specific genes.     

Phylogenetic analysis of proteolytic Acinetobacter strains based on the sequence of genes encoding aminoglycoside 6'-N-acetyltransferases

The sequence of seven aac(6')-I genes encoding aminoglycoside 6'-N-acetyltransferases from proteolytic Acinetobacter strains including genomic species 14, 15, 16, and 17 and from ungrouped proteolytic strains 631, 640, and BM2722 was determined. Pulsed-field gel electrophoresis of genomic DNA of these strains and of Acinetobacter sp. 6 CIP A165 digested with SfiI followed by hybridization with rRNA and aac(6')-I specific probes indicated that these genes were located in the chromosome. Phylogenetic analysis of the genes indicated that aac(6')-I of A. baumannii, Acinetobacter ungrouped strain 631, and Acinetobacter sp. 16 formed a cluster (91.5 to 92.3% identity) whereas aac(6')-I of Acinetobacter sp. 15, sp. 17, and Acinetobacter ungrouped strain BM2722 formed another cluster (90.7 to 94.6% identity). A third cluster was constituted by A. haemolyticus and Acinetobacter sp. 6 (83.6% identity). The phylogeny drawn from aac(6')-I sequences was consistent with that based on DNA-DNA hybridization and phenotype comparison. The aac(6')-I genes were all species specific except for aac(6')-Ih located in a 13.7-kb non conjugative plasmid from A. baumannii BM2686. We conclude that aac(6')-I genes may be suitable for identification at the species level and for analysis of the phylogenetic relationships of Acinetobacter.     

Quantitative microarray-based DNA-DNA hybridization assay for measuring genetic distances among bacterial species and its application to the identification of family Enterobacteriaceae

Quantitative DNA-DNA hybridization to measure the genetic distances among bacterial species is indispensable for taxonomical determination. In the current studies, we developed a method to determine bacterial DNA relatedness on a glass microarray. Reference DNAs representing a total 93 species of Enterobacteriaceae were arrayed on a glass microplate, and signal intensities were measured after 2 hr of hybridization with Cy3-labeled bacterial DNAs. All immobilized DNAs from members of the family Enterobacteriaceae were identified by this method except for DNAs from Yersinia pseudotuberculosis and Y. pestis. These results suggest that quantitative microarray hybridization could be an alternative to conventional DNA-DNA hybridization for measuring chromosome relatedness among bacterial species.     

Similarity and differences in the Lactobacillus acidophilus group identified by polyphasic analysis and comparative genomics

A set of lactobacilli were investigated by polyphasic analysis. Multilocus sequence analysis, DNA typing, microarray analysis, and in silico whole-genome alignments provided a remarkably consistent pattern of similarity within the Lactobacillus acidophilus complex. On microarray analysis, 17 and 5% of the genes from Lactobacillus johnsonii strain NCC533 represented variable and strain-specific genes, respectively, when tested against four independent isolates of L. johnsonii. When projected on the NCC533 genome map, about 10 large clusters of variable genes were identified, and they were enriched around the terminus of replication. A quarter of the variable genes and two-thirds of the strain-specific genes were associated with mobile DNA. Signatures for horizontal gene transfer and modular evolution were found in prophages and in DNA from the exopolysaccharide biosynthesis cluster. On microarray hybridizations, Lactobacillus gasseri strains showed a shift to significantly lower fluorescence intensities than the L. johnsonii test strains, and only genes encoding very conserved cellular functions from L. acidophilus hybridized to the L. johnsonii array. In-silico comparative genomics showed extensive protein sequence similarity and genome synteny of L. johnsonii with L. gasseri, L. acidophilus, and Lactobacillus delbrueckii; moderate synteny with Lactobacillus casei; and scattered X-type sharing of protein sequence identity with the other sequenced lactobacilli. The observation of a stepwise decrease in similarity between the members of the L. acidophilus group suggests a strong element of vertical evolution in a natural phylogenetic group. Modern whole-genome-based techniques are thus a useful adjunct to the clarification of taxonomical relationships in problematic bacterial groups.     

Diversity of Micromonospora strains isolated from nitrogen fixing nodules and rhizosphere of Pisum sativum analyzed by multilocus sequence analysis

It was recently reported that Micromonospora inhabits the intracellular tissues of nitrogen fixing nodules of the wild legume Lupinus angustifolius. To determine if Micromonospora populations are also present in nitrogen fixing nodules of cultivated legumes such as Pisum sativum, we carried out the isolation of this actinobacterium from P. sativum plants collected in two man-managed fields in the region of Castilla and León (Spain). In this work, we describe the isolation of 93 Micromonospora strains recovered from nitrogen fixing nodules and the rhizosphere of P. sativum. The genomic diversity of the strains was analyzed by amplified ribosomal DNA restriction analysis (ARDRA). Forty-six isolates and 34 reference strains were further analyzed using a multilocus sequence analysis scheme developed to address the phylogeny of the genus Micromonospora and to evaluate the species distribution in the two studied habitats. The MLSA results were evaluated by DNA-DNA hybridization to determine their usefulness for the delineation of Micromonospora at the species level. In most cases, DDH values below 70% were obtained with strains that shared a sequence similarity of 98.5% or less. Thus, MLSA studies clearly supported the established taxonomy of the genus Micromonospora and indicated that genomic species could be delineated as groups of strains that share > 98.5% sequence similarity based on the 5 genes selected. The species diversity of the strains isolated from both the rhizosphere and nodules was very high and in many cases the new strains could not be related to any of the currently described species.     

Isolation and characterization of a DNA probe for Staphylococcus aureus subspecies aureus biovar

The gene encoding for polynucleotide phosphorylase (pnp) of a new biovar of Staphylococcus aureus subsp. aureus (NBSA) has been isolated from a genomic library of strain M280(0). The coding region consisted of a 1094-bp HindIII-HindIII DNA fragment encoding for a protein of 277 amino acids with a calculated molecular mass of 29.5 kDa. The nucleotide sequence of the structural gene, contained a continuous open reading frame of 836 bp, showed significant homology with the genes of bacterial polynucleotide phosphorylase from Bacillus subtilis (67.7% identity), from Haemophilus influenzae (62.4% identity), from Pseudomonas luminescens (61.6% identity), and from Escherichia coli (59.7% identity). DNA-DNA and DNA-colony slot-blot hybridizations demonstrated that the pnp gene, employed as a molecular probe, is specific for the identification of NBSA strains.     

Genetic diversity of endophytic diazotrophs of the wild rice, Oryza alta and identification of the new diazotroph, Acinetobacter oryzae sp. nov

Thirty-three endophytic diazotrophs were isolated from surface-sterilized leaves, stem, and roots of wild rice Oryza alta. The SDS-PAGE profile of total protein and insertion sequence-based polymerase chain reaction (IS-PCR) fingerprinting grouped the isolates into four clusters (I-IV). The 16S rRNA gene sequence homology of the representative strains B21, B31, B1, and B23 of clusters I, II, III, and IV were assigned to Pseudomonas oleovorans (99.2% similarity), Burkholderia fungorum (99.4% similarity), Enterobacter cloacae (98.9% similarity), and Acinetobacter johnsonii (98.4% similarity), respectively. The results showed wide genetic diversity of the putative diazotrophic strains of the wild rice, O. alta, and the strains of cluster IV are the first report of nitrogen-fixing Acinetobacter species. The cell size, phenotypic characters, total protein profile, genomic DNA fingerprinting, DNA-DNA hybridization, and antibiotic resistance differentiated strain B23(T) from its closest relatives A. johnsonii LMG999(T) and Acinetobacter haemolyticus LMG996(T). The DNA-DNA hybridization also distinguished the strain B23(T) from the closely related Acinetobacter species. Based on these data, a novel species, Acinetobacter oryzae sp. nov., and strain B23(T) (=LMG25575(T)?=?CGMCC1.10689(T)) as the type strain were proposed.     

Cloning and characterization of resistance gene analogs from Avena species.

Sequences analogous to plant resistance genes of the NBS-LRR class were cloned from the genomic DNA of 11 Avena species with different genomes and levels of ploidy. Three pairs of degenerate primers were used, based on conserved DNA sequence motifs belonging to the NBS domain, and 33 sequences were identified. These were subdivided into 7 classes depending on nucleotide sequence identity. Despite the high level of degeneracy, the primers behaved in a highly selective way; the majority of sequences from the different species obtained with every primer combination showed strong identity and were considered homologous. For most species, only one sequence of each class was identified in each genome, suggesting that duplicated sequences are fairly divergent. The strong identity among specific NBS sequences precludes any conclusions being made on the evolution of these species. The genomic organization of the RGA sequences was explored using those of A. strigosa as probes in Southern blots involving digested DNA from 15 Avena species. The hybridization patterns showed wide diversity both among sequences within a species and among species for each sequence. However, the dendrogram generated using the RFLPs showed relationships among species to be in good agreement with those previously established using other molecular markers.     

Identification,analysis, and utilization of conserved ortholog set markers for comparative genomics in higher plants

We have screened a large tomato EST database against the Arabidopsis genomic sequence and report here the identification of a set of 1025 genes (referred to as a conserved ortholog set, or COS markers) that are single or low copy in both genomes (as determined by computational screens and DNA gel blot hybridization) and that have remained relatively stable in sequence since the early radiation of dicotyledonous plants. These genes were annotated, and a large portion could be assigned to putative functional categories associated with basic metabolic processes, such as energy-generating processes and the biosynthesis and degradation of cellular building blocks. We further demonstrate, through computational screens (e.g., against a Medicago truncatula database) and direct hybridization on genomic DNA of diverse plant species, that these COS markers also are conserved in the genomes of other plant families. Finally, we show that this gene set can be used for comparative mapping studies between highly divergent genomes such as those of tomato and Arabidopsis. This set of COS markers, identified computationally and experimentally, may further studies on comparative genomes and phylogenetics and elucidate the nature of genes conserved throughout plant evolution.     

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     

Evolution of paralogous genes: Reconstruction of genome rearrangements through comparison of multiple genomes within Staphylococcus aureus

Analysis of evolution of paralogous genes in a genome is central to our understanding of genome evolution. Comparison of closely related bacterial genomes, which has provided clues as to how genome sequences evolve under natural conditions, would help in such an analysis. With species Staphylococcus aureus, whole-genome sequences have been decoded for seven strains. We compared their DNA sequences to detect large genome polymorphisms and to deduce mechanisms of genome rearrangements that have formed each of them. We first compared strains N315 and Mu50, which make one of the most closely related strain pairs, at the single-nucleotide resolution to catalogue all the middle-sized (more than 10 bp) to large genome polymorphisms such as indels and substitutions. These polymorphisms include two paralogous gene sets, one in a tandem paralogue gene cluster for toxins in a genomic island and the other in a ribosomal RNA operon. We also focused on two other tandem paralogue gene clusters and type I restriction-modification (RM) genes on the genomic islands. Then we reconstructed rearrangement events responsible for these polymorphisms, in the paralogous genes and the others, with reference to the other five genomes. For the tandem paralogue gene clusters, we were able to infer sequences for homologous recombination generating the change in the repeat number. These sequences were conserved among the repeated paralogous units likely because of their functional importance. The sequence specificity (S) subunit of type I RM systems showed recombination, likely at the homology of a conserved region, between the two variable regions for sequence specificity. We also noticed novel alleles in the ribosomal RNA operons and suggested a role for illegitimate recombination in their formation. These results revealed importance of recombination involving long conserved sequence in the evolution of paralogous genes in the genome.