A database for the provisional identification of species using only genotypes: web-based genome profiling
© Watanabe et al., licensee BioMed Central Ltd 2002
Received: 28 August 2001
Accepted: 5 December 2001
Published: 28 January 2002
For a long time one could not imagine being able to identify species on the basis of genotype only as there were no technological means to do so. But conventional phenotype-based identification requires much effort and a high level of skill, making it almost impossible to analyze a huge number of organisms, as, for example, in microbe-related biological disciplines. Comparative analysis of 16S rRNA has been changing the situation, however. We report here an approach that will allow rapid and accurate phylogenetic comparison of any unknown strain to all known type strains, enabling tentative assignments of strains to species. The approach is based on two main technologies: genome profiling and Internet-based databases.
A complete procedure for provisional identification of species using only their genomes is presented, using random polymerase chain reaction, temperature-gradient gel electrophoresis, image processing to generate 'species-identification dots' (spiddos) and data processing. A database website for this purpose was also constructed and operated successfully. The protocol was standardized to make the system reproducible and reliable. The overall methodology thus established has remarkable aspects in that it enables non-experts to obtain an initial species identification without a lot of effort and is self-developing; that is, species can be determined more definitively as the database is used more and accumulates more genome profiles.
We have devised a methodology that enables provisional identification of species on the basis of their genotypes only. It is most useful for microbe-related disciplines as they face the most serious difficulties in species identification.
A biological species is usually defined in principle as a set of actually or potentially interbreeding organisms, but as interbreeding is very difficult to measure, species have in practice been identified by their phenotypic traits. Until recently, progress in most microbe-related disciplines has been hampered by the enormous effort needed to identify less prominent traits. We are now in an age when we can identify species based on the genome (genotype) , although this does not change the principle that taxonomy is defined by phenotypes : according to the generally accepted rules of taxonomy, a strain belongs to a species if it falls within the range of phenotypes that define that species. This situation has been brought about by the success of the ribosomal RNA approach to phylogenetics [3,4,5]. Well-conserved molecules, such as 16S rRNA in particular, have been used to give a species a molecular identifier and to draw phylogenetic relationships. The 16S rRNA-based approach has been widely accepted and has proved successful in phylogenetic tree-making and even in identifying species. In this context, the Ribosomal Database Project has been established . There are other similar approaches, such as one based on the gyrase gene  and multilocus sequence typing . Nonetheless, it has been impossible in practice to analyze all the constituents of a microbial population, not only because of the huge size of such populations (more than 108 cells per ml) but also because of lack of suitable methodology. Although there are methods other than gene and genome sequencing for analyzing genomes, such as restriction-fragment length polymorphism (RFLP), amplified fragment-length polymorphism (AFLP), Octamer-based genome scanning (OBGS), random polymerase chain reaction (PCR) and others [9,10,11,12], most cannot be used to identify species without a knowledge of phenotypic traits. In reality, there is no general methodology that enables us to identify species by genotype only, although many approaches use genotypic information (DNA sequences) to complement phenotypic information.
We have recently demonstrated the possibility of species identification by genotype using genome profiling , which is a temperature-gradient gel electrophoresis (TGGE) analysis of random PCR products . In particular, the use of 'species-identification dots' (spiddos), which are feature points in genome profiles, is very useful for objective and reproducible data processing [15,16]. We present here a universal method for provisional genotype-based species identification based on these technological advances and using the Internet environment, which enables us to identify species in general. This paper also presents the important concepts of genome distance and genome sequence space, which are essential for species identification based on genotype.
Results and discussion
Key concepts of the on-web genome profiling
As the extent of these changes is roughly proportional to the evolutionary time since the species diverged, we can expect that the summation of the displacement of each spiddo is approximately proportional to the time since divergence, and thus to the genome-to-genome distance. By using a sufficient number of spiddos, we can obtain statistically reliable results. Empirically, we know that 8-10 spiddos, which can be obtained from a single genome profile, can be significant. However, since the more spiddos the better the result, we tentatively made it a rule to adopt four genome profiles (≈ 32-40 spiddos) - that is, four random-PCR products - as a current standard of initial species identification. Therefore, PaSS has the theoretical and empirical basis to be used as a measure of similarity between genomes, although the extent of its effectiveness remains to be shown experimentally as data accumulates. We have introduced a measure of distance, d', obtained from PaSS as formulated in Equation 2 (see Materials and methods), for the sake of convenience .
We call d' a genome sub-distance because it is based not on the whole but a part of the genome sequence. Thus, we introduce (true) genome distance, d, as in Equation 3 (see Materials and methods). Genome distance must have a close relationship with genetic distance, as defined by Nei and others [17,18,19], although there is a difference in the definition. The genetic distance based on sequences is basically the Hamming distance (the number of different letters at each corresponding position of two sequences of letters that are optimally aligned) between two nucleotide (or amino-acid) sequences. In aligning sequences arbitrariness is introduced, depending on the algorithm and parameters used . Another constraint on genetic distance is that it is usually obtained from a limited number of genes, although that is also the case for genome distance. As genome distance is easier to obtain in practice using our method, it should be easier to obtain a lot of data on it compared with genetic distance. On the basis of d (in practice d'), we can construct phylogenetic trees and genome sequence space (an imaginary spherical space in which all the genomes (individuals) can be uniquely located in a finite manner based on the distance between genomes, providing clusters of species (K.N., unpublished observations)).
Although the applicability and effectiveness of genome distance for such purposes needs to be further investigated, it is obvious that an organism that has near-zero genome distance from a certain standard species, as an average over four or more genome sub-distances obtained from as many genome profiles, can be easily assigned to that same species with a high level of confidence. We are not claiming, however, to be able to give the correct taxonomical name to any species using this method. The greater the number of strains registered in the database, the more easily will a species be assigned. Basically, no special efforts, except expanding the database and using sophisticated algorithms, are necessary to raise the proportion of correct assignments. This is the self-developing nature of the database. Therefore, this methodology has two potential great advantages for tentative species identification: first, expertise is not always necessary; and second, database building can be carried out in a self-developing manner (that is, by acquiring more and more accurate data on species) with no waste of information.
Materials and methods
The principle that a species can be identified on the basis of its similarity to a standard species remains unchanged in the shift from phenotype-based to genotype-based methodology. Therefore, the essence of our methodology resides in finding a sufficiently closely related species by way of a measure of similarity - a pattern similarity score (PaSS). Note that this genotype-based methodology cannot define species under the current taxonomy regime, in which phenotype is used as the defining characteristic of species .
General protocol for on-web genome profiling
Although genome profiling is the basic technology for our current purpose, provisional species identification based on genotype can be fulfilled only by using computer-aided database technology, which is most effectively constructed in the Internet environment. As this methodology is based on use by a large number of scientists, the protocol must be designed to be reproducible and easy to carry out. The processes have been deliberately designed with this in view, and are presented on our website .
Genome profiling consists of two basic technologies: random PCR and TGGE, which have been well established [22,23,24]. However, if it is to be used for the purpose of general and universal applications, well-defined standardization is absolutely required to obtain significant results. We have carried out such standardization for genome profiling. The main topics included in the protocol are: preparation of genome DNAs; the set of primers used for random PCR and the internal reference DNAs used for TGGE; experimental conditions for random PCR; and the experimental conditions for TGGE. The protocol also includes the related procedures (extraction of spiddos, calculation of PaSS, and others).
Preparation of genomic DNA
Briefly, the alkaline extraction method was selected for simplicity as follows: 10 mg of cells or tissue are placed in an eppendorf tube and heated for 1 min at 100°C. The cells are mixed with 10 μg 0.5 M NaOH and stirred for 1 min (or 5 min or so for stiffer cells such as yeast) using a microhomogenizer, if necessary, with added quartz sand. Immediately, a 5 μl aliquot of the lysate is mixed with 495 μl 100 mM Tris-HCl (pH 8.0). Usually, a 3 μl aliquot of the mixture thus obtained is used as a template for 100-μl-scale PCR. In some cases, such as Escherichia coli, which does not have a strong cell envelope, these cell-breakdown processes can even be omitted and the cells can be directly used in PCR. In other cases, such as fungi, thorough mechanical treatment (grinding with quartz sand) is needed. Thus, minimal and common procedures are preferred as much as possible for simplicity and generality in so far as they are consistent with the purity and integrity of the DNA samples. DNA samples thus prepared were shown to be identical with those DNAs prepared by the more elaborate conventional method of Thomas  as a PCR template . This seems quite natural, as PCR can be carried out successfully in the presence of contaminating proteins or polysaccharides, irrespective of the DNA cleavages introduced, unless the regions of DNA to be amplified are completely cleaved. Nonspecific binding of proteins, which gives footprint effects, will change the yield but not the molecular ratio of random PCR products as long as the binding is totally stochastic. We also adopt a universal, convenient definition for genome DNA - that it is composed of all DNAs thus prepared, including dynamic elements such as satellite and organelle DNAs, and is irrespective of haploid or diploid status of the cells. Therefore, the DNA samples for genome profiling can be prepared in a common, technically well-defined method for all organisms.
Set of primers for random PCR
Technically important restrictions are introduced by selecting a standard set of primers for random PCR (T.W., A.S., M.N. and K.N., unpublished observations). It is important to carry out random PCR with all kinds of organisms using the same primers so that all species can be compared on the same platform. We have initially selected four oligonucleotides (pfM12: dAGAACGCGCCTG; pfM19: dCAGGGCGCGTAC; d(TGC)3; d(T3G3)2) as a standard set of random PCR primers. The primers pfM12 and pfM19 were selected on the basis of the abundant experimental background on them, whereas d(TGC)3 and d(T3G3)2 were rather theoretically favored (K.N. and A.S., unpublished observations). 'Oligonucleotide-stickiness analysis', which monitors oligonucleotide-binding sites along the template DNA (K.N. and A.S., unpublished observations), was exploited to determine the universal primers and moderately sticky oligonucleotides were selected. These four primers can be fluorescently labeled for convenience. More primers can be used to obtain more detailed information or to supplement insufficient information provided by the four primers about particular pairs of organisms. The information provided by such extra primers can explore in a more detailed manner the local landscape in genome sequence space. In contrast, the standard primers give us rough relationship between any pair of organisms.
Internal reference DNAs
Internal reference bands, which are provided by DNAs of a known melting pattern, are used to calibrate each genome profile, giving highly reproducible results .
Conditions for random PCR
Random PCR is usually carried out under standard conditions: 10 ng template DNA, 50 pmol primer DNA, 250 μM of each dNTP, 50 mM Tris-HCl (pH 8.8), 15 mM (NH4)2SO4, 10 mM MgCl2, 0.45% Triton X-100, 200 μg/ml bovine serum albumin and 2 units of Taq DNA polymerase (Biotech International). PCR was carried out in 30 cycles of 30 sec at 94°C, 2 min at 28°C and 2 min at 47°C, using a thermal cycler PTC-100TM (MJ Research, MA). Annealing temperature can be attenuated depending on the size of the template DNA (in general, the larger the template, the greater the number of DNA fragments generated by random PCR).
Experimental conditions for TGGE
TGGE analysis of random PCR products is carried out with co-migrating internal reference DNAs. TGGE can be either the conventional type or a micronized type . At least two feature points are extracted from the band pattern of the internal reference DNA(s), and then used for calibration of genome profiles or species identification dots (spiddos)  as described below. After calibration, sufficiently high reproducibility of the pattern of spiddos is guaranteed .
Extraction of spiddos
Calculation of PaSS and genome distance
A set of spiddos (around ten), assigned to a genome profile on a computer display, is processed to calculate the normalized mobility and temperature of each point. A measure of similarity of two genomes - the PaSS - is introduced as follows.
d' = (1 - PaSS)/PaSS (2)
Where d'(i) is the ith genome sub-distance obtained with the ith primer used for random PCR.
Computer-aided data acquisition
This study was supported in part by a Grant-in-Aid (09272203) from the Ministry of Education, Science, Sports and Culture of Japan. M.N. was supported by the Japan Society for Promotion of Science (13001147).
- Woese CR, Kandler O, Wheelis ML: Towards a natural system of organisms: Proposal for the domains Archaea, Bacteria, and Eukarya. Proc Natl Acad Sci USA. 1990, 87: 4576-4579.PubMedPubMed CentralView ArticleGoogle Scholar
- Wayne LG, Brenner DJ, Colwell RR, Grimont PAD, Kandler O, Krichevsky MI, Moore LH, Moore WEC, Murray RGE, Stackebrandt E, et al: Report of the ad hoc committee on reconciliation of approaches to bacterial systematics. Int J Syst Bacteriol. 1987, 37: 463-464.View ArticleGoogle Scholar
- Olsen GJ, Overbeek R, Larsen N, Marsh TL, McCaughey MJ, Maciukenas MA, Kuan WM, Macke TJ, Xing Y, Woese CR: The Ribosomal Database Project. Nucleic Acids Res. 1992, Suppl 20: 199-200.Google Scholar
- Preparata RM, Meyer EB, Preparata FP, Simon EM, Vossbrinck CR, Nanney DL: Ciliate evolution: the ribosomal phylogenies of the tetrahymenine ciliates. J Mol Evol. 1989, 28: 427-441.PubMedView ArticleGoogle Scholar
- Martinez JG, Bescos I, Sala JJR, Valera FR: RISSC: a novel database for ribosomal 16S-23S RNA genes spacer regions. Nucleic Acids Res. 2001, 29: 178-180. 10.1093/nar/29.1.178.View ArticleGoogle Scholar
- Maidak BL, Cole JR, Lilburn TG, Parker CT, Saxman PR, Farris RJ, Garrity GM, Olsen GJ, Schmidt TM, Tiedje JM: The RDP-II (Ribosomal Database Project). Nucleic Acids Res. 2001, 29: 173-174. 10.1093/nar/29.1.173.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamamoto S, Harayama S: PCR amplification and direct sequencing of gyrB genes with universal primers and their application to the detection and taxonomic analysis of Pseudomonas putida strains. Appl Environ Microbiol. 1995, 61: 1104-1109.PubMedPubMed CentralGoogle Scholar
- Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA, et al: Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998, 95: 3140-3145. 10.1073/pnas.95.6.3140.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsipouras P: Restriction fragment length polymorphisms. Methods Enzymol. 1987, 145: 205-213.PubMedView ArticleGoogle Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, van de Lee T, Hornes M, Fri-jters A, Pot J, Peleman J, Kuiper M, Zabeau M: AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res. 1995, 23: 4407-4414.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim J, Nietfeldt J, Benson AK: Octamer-based genome scanning distinguishes a unique subpopulation of Escherichia coli O157:H7 strains in cattle. Proc Natl Acad Sci USA. 1999, 96: 13288-13293. 10.1073/pnas.96.23.13288.PubMedPubMed CentralView ArticleGoogle Scholar
- Williams JG, Kubelik AR, Livak KJ, Rafalski JA, Tingey SV: DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucleic Acids Res. 1990, 18: 6531-6535.PubMedPubMed CentralView ArticleGoogle Scholar
- Nishigaki K, Naimuddin M, Hamano K: Genome profiling: a realistic solution for genotype-based identification of species. J Biochem. 2000, 128: 107-112.PubMedView ArticleGoogle Scholar
- Nishigaki K, Amano N, Takasawa T: DNA profiling: an approach of systematic characterization, classification, and comparison of genomic DNAs. Chem Lett. 1991, 1991: 1097-1100.View ArticleGoogle Scholar
- Naimuddin M, Kurazono T, Zhang Y, Watanabe T, Yamaguchi M, Nishigaki K: Species-identification dots: a potent tool for developing genome microbiology. Gene. 2000, 261: 243-250. 10.1016/S0378-1119(00)00502-3.PubMedView ArticleGoogle Scholar
- Biyani M, Nishigaki K: Hundredfold productivity of genome analysis by introduction of microtemperature-gradient gel electrophoresis. Electrophoresis. 2001, 22: 23-28. 10.1002/1522-2683(200101)22:1<23::AID-ELPS23>3.0.CO;2-Z.PubMedView ArticleGoogle Scholar
- Barnabas J, Goodman M, Moore GW: Descent of mammalian alpha globin chain sequences investigated by the maximum parsimony method. J Mol Biol. 1972, 69: 249-278.PubMedView ArticleGoogle Scholar
- Nei M, Chakraborty R: Genetic distance and electrophoretic identity of proteins between taxa. J Mol Evol. 1973, 2: 323-328.PubMedView ArticleGoogle Scholar
- Tateno Y, Nei M, Tajima F: Accuracy of estimated phylogenetic trees from molecular data. I. Distantly related species. J Mol Evol. 1982, 18: 387-404.PubMedView ArticleGoogle Scholar
- Argos P: Sensitive methods for determining the relatedness of proteins with limited sequence homology. Curr Opin Biotechnol. 1994, 5: 361-371.PubMedView ArticleGoogle Scholar
- On-Web GP. [http://gp.fms.saitama-u.ac.jp]
- Nishigaki K, Tsubota M, Miura T, Chonan Y, Husimi Y: Structural analysis of nucleic acids by precise denaturing gradient gel electrophoresis: I. Methodology. J Biochem. 1992, 111: 144-150.PubMedGoogle Scholar
- Henco K, Harders J, Wiese U, Riesner D: Temperature gradient gel electrophoresis (TGGE) for the detection of polymorphic DNA and RNA. Methods Mol Biol. 1994, 31: 211-228. 10.1385/0-89603-258-2:211.PubMedGoogle Scholar
- Wartell RM, Hosseini S, Powell S, Zhu J: Detecting single base substitutions, mismatches and bulges in DNA by temperature gradient gel electrophoresis and related methods. J Chromatogr A. 1998, 806: 169-185. 10.1016/S0021-9673(98)00149-6.PubMedView ArticleGoogle Scholar
- Berns KI, Thomas CA: Isolation of higher molecular weight DNA from Hemophilus influenzae. J Mol Biol. 1965, 11: 476-490.PubMedView ArticleGoogle Scholar
- Hamano K, Takasawa T, Kurazono T, Okuyama Y, Nishigaki K: Genome profiling- establishment and practical evaluation of its methodology. Nikkashi. 1996, 1996: 54-61.Google Scholar
- Wada A, Yabuki S, Husimi Y: Fine structure in the thermal denaturation of DNA: high temperature-resolution spectrophotometric studies. CRC Crit Rev Biochem. 1980, 9: 87-144.PubMedView ArticleGoogle Scholar