Sequence Listing is submitted in triplicate on CD-ROM and is herein incorporated by reference in its entirety. Five tables (14, 15, 19, 20, and 21) submitted on CD-ROM are also incorporated into the specification by reference in their entirety. The files bmx—2003_seq_list.txt, Table 14.txt, Table 15.txt, Table 19.txt, Table 20.txt, and Table 21.txt were saved on Apr. 23, 2004, and are respectively 5138, 3150, 4003, 2049, 1346, and 3517 kilobytes.
Microorganisms are classically identified by their ability to utilize different substrates as a source of carbon and nitrogen through the use of biochemical tests such as the API20E™ system (bioMérieux). For susceptibility testing, clinical microbiology laboratories use methods including disk diffusion, agar dilution and broth microdilution. Although identifications based on biochemical testing and antibacterial susceptibility tests are cost-effective, generally two days are required to obtain preliminary results due to the necessity of two successive overnight incubations to identify the bacteria from clinical specimens as well as to determine their susceptibility to antimicrobial agents. There are some commercially available automated systems (i.e., the MicroScan™ system from Dade Behring and the Vitek™ system from bioMérieux) which use sophisticated and expensive apparatus for faster microbial identification and susceptibility testing (Stager et al., 1992, Clin. Microbiol. Rev. 5: 302-327). These systems require shorter incubations periods, thereby allowing most bacterial identifications and susceptibility testing to be performed in less than 6 hours. Nevertheless, these faster systems always require the primary isolation of the bacteria or fungi as a pure culture, a process which takes at least 18 hours for a pure culture or 2 days for a mixed culture.
Most clinical samples are in the form of blood or urine samples. The remaining samples are in the form of such biological fluids as sputum, pus, cerebrospinal fluid, synovial fluid and the like. The biochemical and susceptibility testing for a urine sample typically requires 18-24 hours of incubation and for blood upwards of 6 to 7 days.
Thus there exists an obvious need for rapid and accurate diagnostic tests for the detection and identification of pathogens. DNA tests are preferred because these tests can be performed more rapidly and accurately than the standard biochemical and susceptibility tests. Thus new DNA tests capable of discriminating between microorganisms are needed.
Bacterial ribosomes contain at least three distinct RNA molecules: 5S, 16S and 23S rRNAs. Historically, these names were chosen with reference to their sedimentation rate, which is reflective of the size of the molecule. However, the true size of the ribosomes from one organism to another varies substantially. Nevertheless, the terminology of 5S, 16S and 23S rRNA is used to describe the ribosomes of all bacteria.
A genetic comparison of the 16S subunits of various bacterial species has shown that there are highly conserved regions intercalated with regions of average and low homology, even in cases of related species. In fact, 16S RNA genes have been used for analyzing the evolutionary relationship between microorganisms. Research groups have used differentially hybridizing DNA probes in order to identify unknown microorganisms based on the hybridization patterns of ribosomal RNA. However, nucleic acid hybridization is an imprecise technique and is ill suited for distinguishing between closely related species and strains of organisms.
Methods of distinguishing between genera and strains for purposes of identification and classification differ and there is no set method. Classification at the genus or species level may be based on DNA/DNA hybridization, whereas identification of the subject organism may be based on a phenotypic character of the organism. Serological reactions, which have only limited value in classification, have enormous value for identification of a particular organism. Serological methods include slide agglutination tests, fluorescent antibody techniques and other serological methods. Although these methods can be performed simply and rapidly, their specificity is frequently not absolute and additional confirmation by physiological or biochemical tests is usually required.
Recently, an effort has been made to identify probes which will differentiate genera, groups, species and strains based on the genetic make-up of the organism. The ability to find probes that distinguish between related species and strains is further complicated by the fact that public databases, such as GenBank, possess accuracy and completeness problems. These problems arise at least because of DNA sequencing errors, and because many bacteria have two or more 16S ribosomal RNA loci in their genomes. Sequence variations may occur between the different copies of the gene present in the same genome, i.e. polymorphisms.
In view of these issues, a great need remains for methods and reagents which can be used to differentiate bacteria based on genus, species and strain for diagnostic, prognostic, environmental, agricultural, and research purposes. Herein are provided reagents and methods for systematically using such reagents to identify bacteria based on genus, species and strain.
One aspect contemplates a plurality of 16S polynucleotides immobilized to a solid support, wherein the plurality of 16S polynucleotides are subsequences of 16S rDNA and each 16S polynucleotide individually comprises at least one distinguishing moiety, which differentiates between microorganisms by genus, group, species, strain and/or isolate. The polynucleotide is preferably an oligomer of about 11 to about 45 nucleotides, and more preferably between 15-30 nucleotides. The plurality can include 10-100 or 5 to 1×106 or more polynucleotides, and any number inbetween.
Another aspect contemplates a method of detecting the presence of a microorganism and determining an isolate, a strain, a species, a group, or a genus of a microorganism in a sample suspected of containing the microorganism comprising the steps of: (A) selecting at least one primer pair to amplify at least a portion of a 16S rDNA of the sample; (B) amplifying the 16S rDNA of the sample with the at least one primer pair; (C) contacting the amplified rDNA with at least one isolated nucleic acid comprising at least one distinguishing moiety; (D) incubating the amplified rDNA and the isolated nucleic acid under hybridizing conditions which allow hybridization in a sequence-specific manner between the sample and the at least one isolated nucleic acid to form a hybridization product; (E) detecting presence of the hybridization product and thereby one or more distinguishing moieties of the microorganism; and (F) determining the isolate, strain, species, group, and/or genus of the microorganism by the presence of the one or more distinguishing moieties.
In yet another aspect, a kit is contemplated. The kit is for the detection and identification of at least one microorganism by genus, group, species, strain and/or isolate in a sample and comprises: (A) at least one primer pair for amplification of at least a portion of a 16S rRNA of the microorganism; (B) two or more nucleic acids comprising at least two critical residues of a 16S rDNA which distinguish the microorganism by genus, group, species, strain or isolate; (C) a hybridization buffer to allow sequence-specific hybridization between the probes and the nucleic acids present in the sample, or to allow sequence-specific hybridization between the probes and the nucleic acids of amplified products of the sample; and (D) a detection moiety.
Another embodiment is a composition comprising a plurality of probes of Table 14 and/or Table 15, (or Table 20 and/or Table 21), wherein each probe comprises at least one distinguishing moiety and wherein the plurality of probes are immobilized on a substrate. The substrate can be a bead, plate, slide, microtube, in the form of an affinity column, and the like. Preferably, the plurality of probes comprises probes that are about 15 to about 45 and more preferably 20 to about 30 nucleotides in length.
Another aspect of the invention contemplates a method of diagnosing a subject and determining the microorganism causing an infection in the subject comprising the steps of: (A) obtaining a sample from the subject; (B) screening the sample for the microorganism using a kit described herein.
Yet a further aspect of the invention contemplates a method of identifying distinguishing moieties in a 16S bacterial rRNA or rDNA comprising the steps of: (A) obtaining a nucleotide sequence of a genetic locus shared by two or more different bacterial strains, species, or genera; (B) dividing the nucleotide sequence into a set of oligomers of length “n” which overlap by “x” nucleotides, wherein “x” is at least one nucleotide less than “n” and wherein said overlapping oligomers span the length of the sequence of the genetic locus; (C) comparing an oligomer using a comparative algorithm against at least one database of nucleotide sequences for that locus from a plurality of bacterial strains, species, or genera, wherein the nucleotide sequences are stored in at least one database; and (D) determining whether the oligomer has a nucleotide sequence which matches, or has no more than one mismatch with, a portion of all available nucleotide sequences for the locus of the strain, species, or genus of origin, or whether the nucleotide sequence has at least two mismatches when aligned with any other strain, species, or genus, wherein the at least two mismatches when aligned correspond to distinguishing moieties which differentiate between strain, species or genus.
This invention provides new methods for identifying genera, group, species, strain, and isolate specific markers based on bacterial rRNA and rDNA sequences. Preferred bacterial rRNA or rDNA is 16S rRNA or rDNA, however in all instances when discussing 16S rRNA or rDNA, unless otherwise noted, 23S rRNA or rDNA is also contemplated.
It is a further object of the invention to provide a method a method of identifying a distinguishing moiety in a 16S or a 23S bacterial rRNA or rDNA comprising:
In a preferred embodiment of the invention “n” is 30 nucleotides and “x” is 15 nucleotides or alternatively “n” is 20 nucleotides and “x” is 19 nucleotides. However “x” can be 9 to 19 nucleotides and “n” is 10 to 30 nucleotides.
The method further comprises analyzing the short oligomers by identifying short oligomers that require the largest number of nucleotide changes to match a different strain, species, or genus than that from which the oligomer derived. The mismatch is preferably one or more. The sequences thus identified can be genus, group, species, strain or isolate specific.
Another embodiment of the invention includes nucleic acids identified by the above method which are isolate, strain, species, group or genus specific. Some of these nucleic acids may be in the forms of probes of sufficient length to bind in a sequence-specific manner to a polynucleotide in a sample or to a polynucleotide amplified from a sample. Optimally such probes will comprise a detectable label and range from 15 to about 60 nucleotides or any number in between.
Another aspect of the invention comprises a method of detecting the presence of and determining the strain, species, group, or genus identity of a microorganism in a sample suspected of containing said microorganism. Such a method can comprise the following steps: (i) optionally selecting and using at least one primer pair to amplify at least a portion of a 16S rDNA of the sample; (ii) contacting the amplified or unamplified sample DNA with at least one of the isolated nucleic acid discussed above; (iii) incubating the amplified or unamplified sample DNA and the isolated nucleic acid under hybridizing conditions which allow hybridization in a sequence-specific manner between the sample and said at least one isolated nucleic acid to form a hybridization product; and (iv) detecting presence of the hybridization product as an indication of the identity of said microorganism. The microorganisms are preferentially from Table 1, but can be any human bacterial pathogen. Alternatively, the sample can be obtained from food or is a biological sample taken from a subject, an environmental sample, or a plant.
Another aspect of the invention contemplates a kit for the detection and identification of at least one microorganism from a set of microorganisms of Table 1 in a sample. The kit comprises comprising: (i) at least one primer pair for amplification of at least a portion of a 16S or a 23S rRNA of a microorganism of Table 1; (ii) a composition comprising at least one probe having at least one distinguishing moiety; (iii) a hybridization buffer to allow sequence-specific hybridization between the probes in said composition and the nucleic acid present in the sample or amplified products of the sample nucleic acid; and (iv) a detection moiety. In a preferred example, the microorganisms are from Table 1.
A further aspect of the invention contemplates a method of using a probe of Tables 14, 15, 20, or 21 (all of which are attached to the specification as separate documents) comprising at least one 16S rRNA distinguishing moiety to detect the presence of a bacterium in a sample. The sample can be an environmental sample, food sample, or a biological sample obtained from an animal. Alternatively, the probes can be used in a method to determine a treatment protocol in a subject believed to suffer from a bacterial infection.
Yet another aspect of the invention includes a database comprising nucleic acid sequences with distinguishing moieties which distinguish a microorganism based on genus, species, or strain and wherein the database comprises at least one sequence from Tables 2, 14, 15, 20, or 21 (all of which are attached to the specification as separate documents). Another database contemplated is a reference database which comprises Table 14 and/or Table 15 or Table 20 and/or Table 21 (in a relational form with a means for querying said reference database).
In another embodiment of the invention, a computerized storage and retrieval system of biological information, comprising: a data entry means; a display means; a programmable central processing unit; and a data storage means having oligomers comprising 16S rRNA sequences with distinguishing moieties and annotated information on attributes of the rRNA sequences electronically stored in a relational database.
Another object of the invention provides for a computerized storage and retrieval system of biological information, comprising: a data entry means; a display means; a programmable central processing unit; and a data storage means having oligomers comprising 16S rRNA sequences with distinguishing moieties and annotated information on attributes of the rRNA sequences electronically stored in a relational database.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
The methods and reagents disclosed herein are for use in differentiating bacteria by genus, group, species and strain. Bacteria are described as single cells or simple associations of similar cells forming a group defined by cellular and not organismal properties. The nucleoplasm (genophore) of a bacterium is never separated from the cytoplasm by a unit-membrane system (nuclear membrane). The super kingdom of bacteria is classified in descending order by phylum, class, family, genus, group, species and strain.
There are four major identification categories of bacteria: (1) Gram-negative eubacteria that have cell walls, (2) Gram-positive eubacteria that have cell walls, (3) eubacteria lacking cell walls such as mycoplasmas, and (4) archaeobacteria. The preferred bacterial targets of the methods and reagents of this invention are 1 through 3.
Of the preferred categories, these can then be subdivided again into groups, any of which can be subjected to the methods and reagents derived therefrom using the methods of the application. For Category 1, the groups include (a) the spirochetes, (b) aerobic/microaerophilic, motile helical/vibroid gram-negative bacteria, (c) non-motile (or rarely motile), gram-negative curved bacteria, (d) gram-negative aerobic/microaerophilic rods and cocci, (e) facultatively anaerobic Gram-negative rods, (f) gram-negative anaerobic straight, curved and helical rods, (g) dissimilatory sulfate- or sulfur-reducing bacteria, (h) anaerobic gram-negative cocci, (i) the rickettsias and chlamydias, (j) anoxygenic phototrophic bacteria, (k) oxygenic phototrophic bacteria, (l) aerobic chemolithotrophic bacteria and associated organisms, (m) budding and/or appendaged bacteria, (n) sheathed bacteria, (o) nonphotosynthetic, non-fruiting gliding bacteria, and (p) fruiting gliding bacteria: the myxobacteria. Of these groups, the preferred groups that include, among others, bacteria of medical importance, are: (a), (b), (c), (d), (e), (f), and (h).
Category 2 can be divided into 6 groups as follows: (a) gram-positive cocci, (b) endospore-forming gram-positive rods and cocci, (c) regular, non-sporing gram-positive rods, (d) irregular, non-sporing gram-positive rods, (e) the mycobacteria, and (f) actinomycetes. Of these groups, the preferred groups that include, among others, bacteria of medical importance are (a), (b), (c), and (d). Category 3 consists of mycoplasmas. For a breakdown and further description of the groups and the members of the groups and descriptions of each, see Bergey's Manual of Determinative Bacteriology (9th ed., John G. Holt et al., eds. Philadelphia, 1994). More preferred organisms are presented infra.
1. Definitions
The “sample” may be any biological material taken either directly from the infected human being (or animal), or after culturing (enrichment), a sample taken from food or feed, an environmental sample comprising a mixture of bacteria, or a plant sample. Biological material may be, e.g., expectorations of any kind, broncheolavages, blood, skin tissue, biopsies, lymphocyte blood culture material, urine, fecal samples, sputum, and the like. Said samples may be prepared or extracted according to any of the techniques known in the art. Sample is also meant to include samples from food, environmental samples or plant samples.
By “subject” is meant to include any vertebrate or invertebrate capable of being infected by a bacterium. Preferred subjects include agricultural animals (e.g., birds, pigs, cattle, sheep, goats, bison, horses and the like), and mammals (e.g., dogs, cats, etc.) including primates (humans).
By “distinguishing moiety” is meant a nucleic acid(s) sequence difference in the sequence of mature 16S rRNA, or in the portion of rDNA which encodes it, that differentiates a bacterium's 16S rRNA or rDNA from another bacterium's 16S rRNA. The distinguishing moiety may be family or genus specific. More preferably the distinguishing moiety is group, species, subspecies or strain specific.
The probes of the invention preferably distinguish between genus, species and strains of bacteria. However, the probes may also be used to distinguish between families and classes of bacteria for classification purposes. By “family specific” is meant a characteristic, preferably at the nucleic acid level, which distinguishes families of bacteria based on their ribosomal DNA or RNA.
By “genus-specific” is meant to include a distinguishing characteristic of a genus which allows differentiation between genera of bacterium based on their ribosomal DNA or RNA. It would be understood that such a genus-specific moiety would also identify the family of a bacterium of that particular genus. Such genus specific moieties may or may not be capable of identifying all members in that family.
By “group specific” is meant an oligonucleotide sequence which is specific towards two or more members of a species. “Group specific” may include species from different genera. By “microorganism group” is meant at least two members of a genus (e.g., Bacteroides) or even members of at least two different genera (e.g., Bacteroides and Prevotella), provided that the species involved are clinically and therapeutically similar. For example, probes, which identify members of a group consisting of Streptococcus pneumoniae, Haemophilus influenzae, and Moraxella catarrhalis, would be useful, because these members of the group all exhibit a similar clinical outcome (namely, pneumonia), which is treated similarly. Thus, knowledge that a bacterium is a member of this group would provide clinicians with useful information. Probes, which identify a microorganism group, will typically hybridize to rRNA or rDNA of a plurality of microorganism species tested.
By “species-specific” is meant to include a distinguishing characteristic of a genus, which allows differentiation between species of bacterium based on their ribosomal DNA or RNA (e.g., Cedecea davisaes, Cedecea lapagei and Cedecea neteri). It would be understood that such a species-specific moiety would also identify the genus of a bacterium of that particular species, but may not necessarily be capable of identifying all members of that genus.
By “strain specific” is meant to include distinguishing characteristics of a bacterium which allows differentiation between strains of bacterium. It would be understood that such a strain-specific moieties would also identify the species of that particular strain, but may not necessarily be capable of identifying all members of that species.
By “rRNA” and “ribosomal RNA” is meant the structural RNA components of the ribosome. Prokaryotes such as bacteria have 5S rRNA and 23S rRNA species in the large subunit and a 16S rRNA species in the small subunit. The large subunit in, for example, E. coli is known to sediment at 50S due to the presence of 23S rRNA, 5S rRNA and 31 proteins. The small subunit, in E. coli, is 30S due to the presence of the 16S rRNA and 21 ribosomal proteins. In E. coli the 16S rRNA is 1,542 nucleotides. This rRNA is synthesized like mRNA using a RNA polymerase that uses DNA as a template.
By “rDNA” is meant typically to refer to the DNA coding for ribosomal RNA. However, in certain instances “rDNA” may refer to recombinant DNA.
By “r-proteins” and “ribosomal proteins” is meant the proteins that associate with the large and small ribosomal subunits.
As used throughout the application and claims, the term “probe” will refer to synthetic or biologically produced nucleic acids, between 10 and 250 bases in length, which by design or selection, contain specific nucleotide sequences that allow specific and preferential hybridization under predetermined conditions to target nucleic acid sequences, and optionally contain a moiety for detection or for enhancing assay performance. A minimum of ten nucleotides is generally necessary in order to statistically obtain specificity and form stable hybridization products, and a maximum of 250 nucleotides generally represents an upper limit for sequences in which reaction parameters can be adjusted to determine mismatched sequences and preferential hybridization. Therefore, in general, a preferred length of a probe will be between 10 and 250 nucleotides. A more preferred length of probe is between 15 and 60 oligonucleotides (i.e., 15-mers and 60-mers and any range in between). Probes may optionally contain certain constituents that pertain to their proper or optimal functioning under certain assay conditions. For example, probes may be modified to improve their resistance to nuclease degradation (e.g., such as by end-capping), to carry detection ligands (e.g., such as fluorescein, 32P, biotin, and like labeling) or to facilitate their capture onto a solid support (e.g., poly-deoxyadenosine “tails”).
“Preferential hybridization” or “hybridizing preferentially” means that hybridization with the intended target nucleic acid results in a hybridization reaction product which is more stable than any hybridization reaction products resulting from hybridization with a non-target nucleic acid under identical conditions. It is well within the skill of the ordinary artisan to compare stability of hybridization reaction products and evaluate which one is more stable, i.e., determine which one has bound “preferentially”. By “specific hybridization” is meant that a nucleotide sequence will hybridize to a predetermined target sequence and will not substantially hybridize to a non-target sequence.
“Specifically discriminate” means that a probe will substantially hybridize to a predetermined target sequence and will not substantially hybridize to a non-target sequence.
“Hybridization” is a process by which, under predetermined reaction conditions, two partially or completely complementary strands of nucleic acids are allowed to come together in an antiparallel fashion to form a double stranded nucleic acid with specific and stable hydrogen bonds, following explicit rules pertaining to which nucleic acids bases may pair with one another. “Substantial hybridization” means that the amount of hybridization observed will be such that one observing the results would consider the result positive in a clinical setting. Data which is considered “background noise” is not substantial hybridization.
“Stringent hybridization conditions” means approximately 35° C. to 65° C. in a salt solution of approximately 0.9 M NaCl. Stringency may also be governed by such reaction parameters as the concentration and type of ionic species present in the hybridization solution, the types and concentrations of denaturing agents present, and the temperature of hybridization. Generally as hybridization conditions become more stringent, probes of greater length and/or containing fewer mismatched nucleotides are preferred if stable hybrids are to be formed. As a rule, the stringency of the conditions under which a hybridization is to take place will dictate certain characteristics of the preferred probes to be employed. Such relationships are well understood and can be readily manipulated by those skilled in the art.
As used herein, the terms “homology” and “homologous to” are meant to refer to the degree of similarity between two or more nucleic acid sequences, and is not meant to imply any taxonomic relatedness between organisms. The degree of similarity is expressed as a percentage, i.e., 90% homology between two sequences will mean that 90% of the bases of the first sequence are identically matched to the bases of the second sequence.
“Target” or “target molecule” as used herein in the diagnostic sense, refers to a molecule of interest, i.e., the molecule whose presence one wishes to determine. In a therapeutic sense, the term “target” or “target molecule” refers to a molecule associated with a disease or with an organism causing a disease.
“Biological binding pair” as used herein refers to any pair of molecules, which exhibit mutual affinity or binding capacity. A biological binding pair is capable of forming a complex under binding conditions. For the purposes of the present application, the term “ligand” will refer to one molecule of the biological binding pair; and the terms “antiligand” or “receptor” will refer to the opposite molecule of the biological binding pair. For example, without limitation, embodiments of the present invention have application in nucleic acid hybridization assays where the biological binding pair includes two complementary nucleic acids. One of the nucleic acids is designated the ligand and the other nucleic acid is designated the antiligand or receptor. One of the nucleic acids may also be a target molecule. The designation of ligand or antiligand is a matter of arbitrary convenience. The biological binding pair may include antigens and antibodies, drugs and drug receptor sites, and enzymes and enzyme substrates, to name a few.
The term “label” refers to a chemical moiety which is capable of detection including, by way of example, without limitation, radioactive isotopes (e.g., 32P), enzymes (e.g., horseradish peroxidase), luminescent agents (e.g., luciferase), precipitating agents, and dyes (e.g., FAM which is derived from fluoresceine and TAMRA which is derived from rhodamine). The term “agent” is used in a broad sense, including any chemical moiety, which participates in reactions that lead to a detectable response. The term “cofactor” is used broadly to include any chemical moiety, which participates in reactions with the label.
The term “amplify” is used in the broad sense to mean creating an amplification product, which may include by way of example, additional target molecules, or target-like molecules, capable of functioning in a manner like the target molecule, or a molecule subject to detection steps in place of the target molecule, which molecules are created by virtue of the presence of the target molecule in the sample. In the situation where the target is a polynucleotide, additional target, or target-like molecules, or molecules subject to detection can be made enzymatically with DNA or RNA polymerase.
The term “support” when used alone, includes conventional supports such as filters, dipsticks, microarrays, beads and membranes as well as retrievable supports.
2. Method of Identifying Group, Genus, Species and Strain Specific Polymorphisms
2.1 In Silico Methods. A systematic method is described herein for identifying regions of a 16S sequence, for any bacterial organism, as being specific to particular genera or species, and in some instance, specific to a particular strain or isolate. As illustrated in
As illustrated in
Oligomers exhibiting 0-1 nucleotide mismatches when compared to nucleotide sequences associated with 0-z number of species are categorized as either species-specific oligomers or group-specific oligomers, as shown in
More specific examples of the aforementioned method are described below (see Example 2 and Example 3). Example 2 specifically provides an example where n equals 30 and x equals 15. Example 3 specifically provides an example where n equals 20 and x equals 19.
2.2 Experimental verification. Experimental verification can then be performed using, for example, real time PCR (TaqMan®), spot format assays (e.g., a single hybridization test on a membrane) or low-density DNA microarray format (e.g., multiple probes spotted on a chip). The TaqMan® technique was developed by Perkin-Elmer Applied-Biosystems and relies on fluorescence resonance energy transfer (FRET). It requires the use of two primers and a fluorescent probe. The TaqMan® probes possess a fluorescent reporter at the 5′ end and a fluorescent quencher at the 3′ end. When irradiated, the excited fluorescent dye transfers energy to the nearby quenching dye molecule rather than fluorescing, resulting in a non-fluorescent substrate. When the probe hybridizes to its complementary sequence on the target DNA, the Taq DNA polymerase starts to digest the probe separating the reporter and the quencher so that no FRET occurs and the fluorescent signal of the reporter can be detected and measured. The more PCR product formed, the higher the fluorescent signal observed.
To accurately test the specificity of a TaqMan® probe, the PCR primers used in the TaqMan® assay must represent sequences that are present in all species to be tested. In order to verify the presence of conserved sequences upstream and downstream from the 30-mer oligonucleotide sequence from which to select PCR primers, each 30-mer was extended by 60 bases on each side to generate a 150-mer oligonucleotide. A multiple sequence alignment using CLUSTALW was generated using the particular 150-mer oligonucleotide to be tested and the 20-40 closest sequences in the selected database(s). CLUSTALW is the more recent version of CLUSTAL, with the W standing for weighting to represent the ability of the program to provide weights to the sequence and program parameters (Higgins et al., 1996, Meth. Enz. 266: 383-402). Another program that can be used is for multiple sequence alignment is PILEUP, and the like.
3. Organisms
The methods described infra can be utilized with any bacteria, such as those discussed in John G. Holt et al., Bergey's Manual of Determinative Bacteriology (Philadelphia, Pa., 9th ed., 1994). The approximately 600 species and about 80 genera, whose rDNA sequences were analyzed to discover the markers of this invention, are provided in Table 7. Preferred genera include Actinobacillus spp., Actinomyces spp., Bacteroides spp., Campylobacter spp., Cardiobacterium spp., Clostridium spp., Eikenella spp., Enterobacter spp., Enterococcus spp., Escherichia spp., Fusobacterium spp., Haemophilus spp., Kingella spp., Klebsiella spp., Moraxella spp., Neisseria spp., Oligella spp., Prevotella spp., Propionibacterium spp., Pseudomonas spp., Staphylococcus spp., Streptococcus and Wolinella spp. Any of the species within these individual groups of species can be analyzed using the methods, reagents, and the like described herein. Preferred species, include the following:
Actinomyces species occur mainly in the oral cavity and on mucous membranes of warm-blooded vertebrates. Actinomyces spp. commonly cause pyogenic infections in association with other concomitant bacteria. Differentiation of the genus Actinomyces is difficult because of variable test reactions and because some species are best differentiated by protein gel electrophoresis, which requires culture and thus slows analysis.
Bacteroides are isolated from a wide range of anaerobic habitats including gingival crevices, the intestinal tract, sewage sludge and infective and purulent conditions in humans and animals.
Campylobacter is found in the reproductive organs, intestinal tract and oral cavity of humans and animals.
Cardiobacterium are associated with endocarditis in humans.
Cedecea were first thought to be intermediate between typical Serratia species and S. fonticola. They also are similar to Ewingella and thus can be hard to distinguish by traditional phenotypic methods.
Clostridum species are widespread in the environment. Many species produce potent endotoxins and some are pathogenic for animals, because of either wound infections or the absorption of toxins. There are over 100 species in this genus and differentiation between the species needs to be carried out only with well trained personnel because of the exacting growth conditions and test procedures required for distinguishing the species.
Eikenella can be opportunistic pathogens that cause infections in the human mouth and intestine.
Enterobacter species are widely distributed in nature, occurring in fresh water, soil, sewage, plants, vegetables and animal and human feces. Several species, most notably E. cloacae, E. sakazakii, E. aerogenes, E. agglomerans, and E. gergoviae are opportunistic pathogens causing urinary tract infections and occasionally septicemia and meningitis.
Enterococcus species occur widely in the environment, particularly in feces of vertebrates. Enterococcus species are also sometimes the cause of pyogenic infections.
Although Escherichia species are part of the normal intestinal flora of warm-blooded animals and in the case of E. blattae, of cockroaches, Escherichia are also responsible for diarrheal disease, and are major contributors to urinary tract infections and nosocomial infections, including septicemia and meningitis.
Fusobacterium is found principally in the gingival sulcus and in the intestinal and genital tracts. Species have also been isolated from blood cultures and various purulent lesions in humans and animals and from tropical ulcers. Because the typical fusiform or spindle-shaped appearance of these organisms is not shown by all fusobacteria, it causes difficulty in distinguishing some species from strains of Bacteroides, Clostridium and Eubacterium where the cells are easily decolorized. Differentiation of the species belonging to the Fusobacterium genus is also difficult because of the general lack of reactivity in conventional tests. This is especially true for subspecies of F. nucleatum: subsp. nucleatum, subsp. polymorphym and subsp. vincentii.
The genus Haemophilus includes H. influenzae, which is the leading cause of meningitis in children. It also is known to cause septicemic conditions, otitis media, sinusitis, and chronic bronchitis. H. influenzae biovar aegyptius is mainly responsible for conjunctivitis, a highly transmissible eye infection. H. ducreyi is the causative agent of the venereal disease soft chancre or chancroid.
Kingella species are Gram-negative cells, but because they have tendency to resist decolorization, they are difficult to characterize clinically. These organisms occur in human mucous membranes of the upper respiratory tract and are known to be susceptible to penicillin. Consequently, rapid identification of Kingella spp. would allow clinicians to prescribe drugs such as penicillin over more pharmacologically toxic antibiotics.
Klebsiella species can be found in human feces, clinical specimens, soil, water, grain, fruits and vegetables. K. pneumoniae and K. oxytoca are opportunistic pathogens that can cause bacteremia, pneumonia and urinary tract and other human infections. Klebsiella species are also known for causing nosocomial infections in urological, neonatal, intensive care and geriatric patients.
Species of Moraxella are Gram-negative, but as with other microorganisms often have a tendency to resist decolorization of the dye, thus making a clinical diagnosis difficult. Diagnosis of a Moraxella spp. infection is useful, as they are usually highly sensitive to penicillin. These organisms are responsible for parasitic infections of mucous membranes of humans and other warm-blooded animals.
Species of Neisseria are also Gram-negative, and also have a tendency to resist decolorization of the dye, thus making a clinical diagnosis difficult. These organisms typically are inhabitants of the mucous membranes of mammals. Some species are primary pathogens for humans.
The Oligella genus was not created until 1987 (Rossaue et al., 1987, Int. J. Syst. Bacteriol. 37: 198-210) and contains two species: O. urethralis and O. ureolytica. These bacterium are isolated mainly from the genitourinary tract of humans.
Pantoea is a genus that was not created until 1989 by Gavini et al. (Int. J. Syst. Bacteriol. 39: 337-45). It included two species at that time, Pantoea agglomerans (also known as Enterobacter agglomerans, Erwinia herbicola, and Erwinia milletiae) and Pantoea dispersa. Since its creation, the Pantoa genus has grown to include at least seven species. The species of Pantoea are isolated from plant surfaces, seeds, soil, and water, as well as from animals and human wounds, blood and urine. These microorganisms are also considered opportunistic human pathogens.
The genus Prevotella was created in 1990 by Shah and Collins (Int. J. Syst. Bacteriol. 40: 205-208) with the species P. melaninogenica (formerly Bacteroides melaninogenicus). At least fifteen other species have been added to the Prevotella genus from the genus Bacteroides. Consequently, discrimination between at least these two genera is difficult. Additional classification based on species is required and would be aid by species specific probes.
Propionibacterium species are found mainly in cheese and dairy food products and on the human skin. Propionibacterium are readily confused with some species of Cornebacterium or Clostridium.
Proteus is a Gram-negative genus of organisms that occur in intestines of humans and a wide variety of animals. Species are also found in manure, soil and polluted waters. P. myxofaciens has been isolated from gypsy moth larvae. Species, which are human pathogens, cause urinary tract infections. Certain species also act as secondary invaders causing septic lesions often in burn patients. The species, Proteus penneri, was first created in 1982 by Hickman et al., J. Clin. Microbiol. 15: 1097-1102 and Int. J. Syst. Bacteriol. 33: 438-440.
The genus, Providencia, consists of microorganisms that have been isolated from human diarrhetic stools, urinary tract infections, wounds, burns and bacteremias, and from penguins. It is thus considered to include human pathogen species. The species Providencia heimbachae was created in 1986 by Muller et al., (Int. J. Syst. Bacteriol. 36: 252-6) and Providencia rustigianii (formerly known as biogroup 3 of Providencia alcalifaciens) was created in 1983 by Hickman-Brenner et al., (J. Clin. Microbiol. 17: 1057-60; and Int. J. Syst. Bacteriol. 33: 672-4).
Pseudomonas spp. are widely distributed in nature. Some species are pathogenic for humans, animals or plants.
Staphylococcus species are mainly associated with the skin and mucous membranes of warm-blooded vertebrates but are often isolated from food products, dust and water. Some species are opportunistic pathogens in humans and animals or are known to produce extracellular toxins. Consequently, subspecies and strain determination is important for clinical diagnosis, prognosis and treatment of Staphylococcus infections.
Streptococcus is a complex genus. The genus in Bergey's Manual of Systematic Bacteriology included microorganisms that were members of the Enterococcus, Lactococcus and Streptococcus genera. These genera broadly encompass enterococci, the lactic streptococci and the pyrogenic, oral and anaerobic streptococci. Additionally, the genus Melissococcus (q.v.) contains an organism previously known as Streptococcus pluton. Additionally, most of the anaerobic streptococci (whose speciation is currently confused) will be transferred to genera such as Peptostreptococcus, but they are currently included in the Streptococcus genus. Streptococci are parasites of vertebrates that mainly inhabit the mouth and upper respiratory tract. Some species are pathogenic for humans and animals. Various antigens associated with Lancefield serological groups are characteristic of some of the species and are required for accurate identification of Streptococcus species. There is difficulty in differentiating the strains that belong to the pyrogenic, oral and anaerobic groups of Streptococcus especially. For example, many species especially in the oral group are undergoing active study, with consequent rearrangement of taxonomy and continuing emendation of description so that a number of areas are still not very clear. See John G. Holt et al., Bergey's Manual of Determinative Bacteriology (9th ed., Lippincott Williams & Wilkins 1994).
Wolinella species are isolated from the bovine rumen, human gingival sulcus, dental root canal infections and other clinical material.
The genus Yokenella was not included in Bergey's Manual of Systematic Bacteriology. The genus was created in 1984 by Kosako et al., (Japan. J. Med. Sci. Biol. 37: 117-24; Int. J. Med. Sci. Biol. 35: 223-5, 1985). It included one species, Y. regensburgei. The genus Koserella with its single species, K. trabulsii, was created in 1985 by Hickman-Brenner et al., (J. Clin. Microbiol. 21: 39-42; Int. J. Syst. Bacteriol. 35: 223-5, 1985). Y. regensburgei and K. trabulsii were shown to be subject synonyms in 1987 by Kosako et al., (Int. J. Syst. Bacteriol. 37: 127-9). Other synonyms for this organism are Hafnia hybridization group 3, Enteric Group 45 (CDC), and NIH biogroup 9 (Japan). Today, the predominant name is J. regensburgei. The organism is isolated from human wounds, urine, sputum and stool and insect intestine, although its clinical significance is unknown. However, there is difficulty characterizing the microorganism. Both human and insect strains of Yokenella were first thought to be Hafnia alvei or a Hafnia-like species. Yokenella is also somewhat similar to species in the genera Citrobacter and Escherichia. Finally, it is further difficult to characterize Yokenella because several strains frequently give delayed (i.e. 3-7 days) positive reactions for several biochemical tests. Thus, a method of positively identifying this organism will be helpful clinically.
Novel sequences have been identified for the following species using the methods described herein.
4. Amplification
Amplification of genes (or DNA and RNA) can be performed using polymerase chain reaction (PCR) and other rapid amplification procedures known in the art, such as ligase chain reaction (LCR), transcription-mediated amplification (TMA), self-sustained sequence replication transcription, self-sustained sequence replication (3SR), nucleic acid sequence-based amplification (NASBA), strand displacement amplification (SDA), branched DNA (bDNA), cycling probe technology (CPR), solid phase amplification (SPA), rolling circle amplification technology (RCA), solid phase RCA, anchored SDA and nuclease dependent signal amplification (DNSA) (Lee et al., 1997, Nucleic Acid Amplification Technologies: Application to Disease Diagnosis, Eaton Publishing, Boston, Mass.; Persing et al., 1993 Diagnostic Molecular Microbiology: Principles and Applications, Amer. Soc. Microbiol., Washington, D.C.; Westin et al., 2000, Nat. Biotechnol. 18:199-204). In the instance of amplifying rRNA, the transcripts can first be reverse-transcribed using RT-PCR and then any of the above cited amplification methods can be used to amplify the DNA transcripts obtained from RT-PCR. Any of these methods of rapid amplification can be used according to the basic concept of identifying sequence specific changes which discriminate microorganisms bases on genus, species, group, family and the like.
For DNA amplification by the widely used PCR (polymerase chain reaction) method, primer pairs are derived either from the DNA sequences surrounding the probe (e.g., TaqMan® primer design) or from data bank sequences. Prior to synthesis, the potential primer pairs can be analyzed by using the program Oligo™ 4.0 (National Biosciences, Plymouth, Minn.) to verify that they are likely candidates for PCR amplifications.
During DNA amplification by PCR, two oligonucleotide primers binding respectively to each strand of the denatured double-stranded target DNA from the bacterial genome are used to amplify exponentially in vitro the target DNA by successive thermal cycles allowing denaturation of the DNA, annealing of the primers and synthesis of two new DNA strands of the target DNA at each cycle (Persing et al, 1993, Diagnostic Molecular Microbiology: Principles and Applications, American Society for Microbiology, Washington, D.C.). Briefly, the PCR protocols may be as follows. Clinical specimens or bacterial colonies or extracted genomic DNA were added directly to the 50 uL PCR reaction mixtures containing 50 mM KCl, 10 mM Tris-HCl pH 8.3, 1.5 mM MgCl2, 0.4 uM of each of the two primers, 200 uM of each of the four dNTPs and 1.25 Units of Taq DNA polymerase (Perkin Elmer).
PCR reactions can then subjected to thermal cycling (e.g., 3 min at 95° C. followed by 30 cycles of 30 sec at 95° C., 30 sec at 55° C. and 45 sec at 72° C.) using, for example, a Perkin Elmer Gene Amp System 2400 thermal cycler. Amplified DNA products can then be analyzed by standard ethidium bromide-stained agarose gel electrophoresis. It is clear that other methods for the detection of specific amplification products, which may be faster and more practical for routine diagnosis, may be used. Such methods may be based on the detection of fluorescence signal after amplification (e.g., TaqMan® system from Perkin Elmer or Amplisensor™ from Biotronics Technology Corp., Lowell, Mass.) or liquid hybridization with an oligonucleotide probe binding to internal sequences of the specific amplification product. These novel probes can be generated preferably from genus-specific or species-specific fragment probes. Methods based on the detection of fluorescence are particularly useful for utilization in routine diagnosis as they are, very rapid and quantitative and can be automated.
To assure PCR efficiency, glycerol or dimethyl sulfoxide (DMSO) or other related solvents, can be used to increase the sensitivity of the PCR and to overcome problems associated with the amplification of target with a high GC content or with strong secondary structures. The concentration ranges for glycerol and DMSO are about 5-10% (v/v) and about 3-10% (v/v), respectively. For the PCR reaction mixture, the concentration ranges for the amplification primers and the MgCl2 are about 0.1-1.0 μM and about 1.5-3.5 mM, respectively. Modifications of the standard PCR protocol using external and nested primers (i.e., nested PCR) or using more than one primer pair (i.e., multiplex PCR) may also be used (see Persing et al., 1993, Diagnostic Molecular Microbiology: Principles and Applications, American Society for Microbiology, Washington, D.C.).
The person skilled in the art of DNA amplification knows the existence of other rapid amplification procedures such as ligase chain reaction (LCR), transcription-based amplification systems (TAS), self-sustained sequence replication (3SR), nucleic acid sequence-based amplification (NASBA), strand displacement amplification (SDA) and branched DNA (bDNA) (Persing et al., 1993). The scope of this invention is not limited to the use of amplification by PCR, but rather includes the use of any rapid nucleic acid amplification methods or any other procedures, which may be used to increase rapidity and sensitivity of the tests. Any oligonucleotides suitable for nucleic acid amplification can also use approaches other than PCR and are contemplated for use herein.
Amplification products are classically detected by standard ethidium bromide-stained agarose gel electrophoresis. It is clear that other methods for the detection of specific amplification products, which may be faster and more practical, can be used. Such methods may be based on detection of fluorescence after or during amplification. One simple method for monitoring amplified DNA is to measure its rate of formation by measuring the increase in fluorescence of intercalating agents such as ethidium bromide or SYBR® Green I (Molecular Probes). If more specific detection is required, fluorescence-based technologies can monitor the appearance of a specific product during the reaction. The use of dual-labeled fluorogenic probes such as in the TaqMan® system (Applied Biosystems), which utilizes the 5′ to 3′ exonuclease activity of the Taq polymerase, is a good example (Livak et al., 1995, PCR Methods Appl. 4: 357-62). TaqMan® can be performed during amplification and this “real time” detection can be done in a single closed tube, which eliminates post-PCR sample handling and consequently prevents the risk of amplicon carryover. Several other fluorescence-based detection methods can be used in real-time.
Fluorescence resonance energy transfer (FRET) is the principle behind the use of adjacent hybridization probes (Wittwer et al., 1997 BioTechniques 22: 130-1, 138-8) and molecular beacons (Tyagi et al., 1996, Nature Biotech. 14: 303-8). Adjacent hybridization probes are designed to be internal to the amplification primers. The 3′ end of one probe is labeled with a donor fluorophore while the 5′ end of an adjacent probe is labeled with an acceptor fluorophore. When the two probes are specifically hybridized in close proximity (spaced by 1 to 5 nucleotides) the donor fluorophores, which has been excited by an external light source emits light that is absorbed by a second acceptor that emit more fluorescence and yields a FRET signal.
Molecular beacons possess a stem-and-loop structure where the loop is the probe and at the bottom of the stem a fluorescent moiety is at one end while a quenching moiety is at the other end. The beacons undergo a fluorogenic conformation change when they hybridize their targets hence separating the fluorochrome from its quencher. The FRET principle is also used in an air thermal cycler with a built-in fluorometer (Wittwer et al., 1997). The amplification and detection are extremely rapid as reactions are performed in capillaries; it takes only 18 min. to complete 45 cycles. Those techniques are suitable especially in the case where few pathogens are searched for. Boehringer-Roche Inc. sell the LightCycler™, and Cepheid makes the SmartCycler. These two apparatus are capable of rapid cycle PCR combined with fluorescent SYBR® Green I or FRET detection.
Microbial pathogen detection and identification may also be performed by solid support or liquid hybridization using, e.g., genus-, species-, and/or strain-specific internal DNA probes, which hybridize to an amplification product. Such probes may be generated from any sequence discussed herein and designed to specifically hybridize to DNA amplification products that are objects of the present invention. The oligonucleotide probes may be labeled with biotin or with digoxigenin or with any other suitable reporter molecule. Hybridization on a solid support is amenable to miniaturization.
Preferred methods include oligonucleotide microarray technology. Currently, low to medium density arrays are available and can be used to capture fluorescent labeled amplicons (Heller et al., “An integrated microelectronics hybridization system for genomic research and diagnostic applications,” in Harrison et al., 1998, Micrototal Analysis Systems '98 (Kluwer Academic Publisher, Dordrecht). Detection methods of hybridization are not limited to fluorescence. Potentiometry, colorimetry and plasmon resonance are some examples of alternative detection methods. In addition to detection by hybridization, nucleic acid microarrays could be used to perform rapid sequencing by hybridization.
4.1 Amplification Primers. Amplification primers need to be designed in order to amplify the DNA regions that contain probes (i.e., species-specific oligonucleotide sequences) of interest. Multiple sequence alignments are made using the 16S (or 23S) rDNA sequence of the target species and 16S sequences of other species phylogenetically related to the target species, as well as species phylogenetically distant from the target species.
A typical PCR amplification requires the use of two primers, one called a Forward primer (typically designated by “F”) and the other one called a Reverse primer (typically designated by “R”). The Forward primer is located upstream of the probe sequence, the Reverse primer is located downstream of the probe sequence.
Suitable primers strongly depend on the type of applications/technologies used for probe testing. If the species-specificity exclusively relies on the probe sequence, regions suitable for amplification primers should be as conserved as possible among all species included in the multiple sequence alignment. If, on the other hand, the user wishes to use the locus amplification process to achieve some of the specificity of the diagnostic test, then the regions suitable for amplification primers should be conserved for the target species and different for all the other species from which the target species should be differentiated.
Once the regions for amplification primers are selected, several characteristics need to be checked to ensure the formation of a good PCR product. These characteristics include one or more of the following:
Probes are classified as primers and/or probes for use in amplifying and/or isolating genus-, species- and strain-specific moieties. Then there are the probes used for diagnostic, prognostic, and research purposes that allow the categorization of the bacterium based on phylum, family, genus, species and strain. The latter probes will contain a distinguishing moiety that will allow the probe to distinguish between families, genera, species and/or strains of bacteria. However, regardless of the category of probe, most will share similar characteristics and will follow the similar guidelines for their preparation and selection.
In general, “primers” as used herein are oligonucleotides that are used for amplifying the regions containing the differentiation moieties. “Probes” are oligonucleotides that are used to distinguish between genera, groups, species, and strains of bacteria. For the instant invention, primers can range in size from 15 nucleotides to 50 nucleotides, and more preferably from 17 nucleotides to 30 nucleotides, and most preferably from 20 nucleotides to 25 nucleotides. Probe size can range in size from 10 to 100 nucleotides, and more preferably from 15 nucleotides to 30 nucleotides, and most preferably from 20 nucleotides to 25 nucleotides.
Because the extent and specificity of hybridization reactions such as those described herein are affected by a number of factors, manipulation of one or more of those factors will determine the exact sensitivity and specificity of a particular probe, whether perfectly complementary (i.e., an oligonucleotide chain in which all of the bases are able to form base pairs with a sequence of bases in another polynucleotide chain) to its target or not. The importance and effect of various assay conditions, explained further herein, are known to those skilled in the art.
First, the stability of the [probe:target] nucleic acid hybrid should be chosen to be compatible with the assay conditions. This may be accomplished by avoiding long A and T rich sequences, by terminating the hybrids with G:C base pairs, and by designing the probe with an appropriate Tm. The beginning and end points of the probe should be chosen so that the length and % GC result in a Tm about 2-10° C. higher than the temperature at which the final assay will be performed. The base composition of the probe is significant because G-C base pairs exhibit greater thermal stability as compared to A-T base pairs due to additional hydrogen bonding. Thus, hybridization involving complementary nucleic acids of higher G-C content will be stable at higher temperatures. Probe selection for use in the in silico method is dependent on the overlapping nature of the oligomers (e.g., 30-mer) used. Once the distinguishing moiety is identified, then suitable diagnostic, prognostic and research probes can be prepared, according to these guidelines and what is generally known in the art, which comprise the distinguishing moiety (e.g., one or more nucleotides which distinguish, for example, different species of bacteria).
Conditions such as ionic strength and incubation temperature under which a probe will be used should also be taken into account in constructing a probe. It is known that hybridization will increase as the ionic strength of the reaction mixture increases, and that the thermal stability of the hybrids will increase with increasing ionic strength. On the other hand, chemical reagents, such as formamide, urea, DMSO and alcohols, which disrupt hydrogen bonds, will increase the stringency of hybridization. Destabilization of the hydrogen bonds by such reagents can greatly reduce the Tm. In general, optimal hybridization for synthetic oligonucleotide probes of about 10-50 bases in length occurs approximately 5° C. below the melting temperature (Tm) for a given duplex. Incubation at temperatures below the optimum may allow mismatched base sequences to hybridize and can therefore result in reduced specificity.
It is desirable to have probes that hybridize only under conditions of high stringency. Under high stringency conditions, only highly complementary nucleic acid hybrids will form; hybrids without a sufficient degree of complementarity will not form. Accordingly, the stringency of the assay conditions determines the amount of complementarity needed between two nucleic acid strands forming a hybrid. Stringency is chosen to maximize the difference in stability between the hybrid formed with the target and the non-target nucleic acid. In some examples of the current invention, it may be necessary to detect single base pair changes. In those instances, conditions of very high stringency are needed.
Second, probes should be positioned so as to minimize the stability of the [probe:non-target] nucleic acid hybrid. This may be accomplished by minimizing the length of perfect complementarity to non-target organisms, avoiding GC rich regions of homology to non-target sequences, and by positioning the probe to span as many destabilizing mismatches as possible. Whether a probe sequence is useful to detect only a specific type of organism depends largely on the thermal stability difference between [probe:target] hybrids and [probe:non-target] hybrids. In designing probes, the differences in these Tm values should be as large as possible (e.g., at least about 2° C. and preferably about 5° C.).
The length of the target nucleic acid sequence and, accordingly, the length of the probe sequence can also be important. In some cases, there may be several sequences from a particular region, varying in location and length, which will yield probes with the desired hybridization characteristics. In other cases, one sequence may be significantly better than another which differs merely by a single base. While it is possible for nucleic acids that are not perfectly complementary to hybridize, the longest stretch of perfectly complementary base sequence will normally primarily determine hybrid stability. While oligonucleotide probes of different lengths and base composition may be used, oligonucleotide probes preferred in this invention are between about 10 to about 50 bases in length, and most preferred 20-30 bases in length, and are sufficiently homologous to the target nucleic acid.
Third, regions in the target DNA or RNA that are known to form strong internal structures inhibitory to hybridization are less preferred. In the preferred embodiment, the probe will hybridize to an rDNA. However, hybridization between probes and rRNA, DNA prepared from rDNA, and cDNA prepared from rRNA are also contemplated. Likewise, probes with extensive self-complementarity should be avoided. As explained above, hybridization is the association of two single strands of complementary nucleic acids to form a hydrogen bonded double strand. It is implicit that if one of the two strands is wholly or partially involved in a hybrid that it will be less able to participate in formation of a new hybrid. There can be intramolecular and intermolecular hybrids formed within the molecules of one type of probe if there is sufficient self-complementarity. Such structures can be avoided through careful probe design. By designing a probe so that a substantial portion of the sequence of interest is single stranded, the rate and extent of hybridization may be greatly increased. Computer programs are available to search for this type of interaction. However, in certain instances, it may not be possible to avoid this type of interaction.
The probes of the present invention are designed for attaining optimal performance under the same hybridization conditions so that they can be used in sets for simultaneous hybridization; this highly increases the usability of these probes and results in a significant gain in time and labor. Evidently, when other hybridization conditions should be preferred, all probes should be adapted accordingly by adding or deleting a number of nucleotides at their extremities. It should be understood that these concomitant adaptations should give rise to essentially the same result, namely that the respective probes still hybridize specifically with the defined target. Such adaptations might also be necessary if the amplified material should be RNA (e.g., rRNA) in nature and not DNA.
The hybridization conditions can be monitored relying upon several parameters, such as the nature and concentration of the components of the media, and the temperatures under which the hybrids are formed and washed.
The hybridization and wash temperature is limited in upper value depending on the sequence of the probe (i.e., its nucleic acid composition, kind and length). The maximum hybridization or wash temperature of the probes described in the present invention ranges from about 40° C. to about 60° C., more preferably from about 45° C. to about 55° C., in the specific hybridization and wash media as described in the Examples section. At higher temperatures, duplexing (i.e., formation of the hybrids) competes with the dissociation (i.e., or denaturation) of the hybrid formed between the probe and the target.
In a preferred hybridization medium of the invention, containing 3×SSC and 20% formamide, hybridization temperatures can range from about 45° C. to about 55° C., with a preferred hybridization temperature of 50° C. A preferred wash medium contains 3×SSC and 20% formamide, and preferred wash temperatures are the same as the preferred hybridization temperatures, i.e., preferably between 45° C. and 55° C., and most preferably 50° C.
However, when modifications are introduced, be it either in the probes or in the media, the temperatures at which the probes can be used to obtain the required specificity should be changed according to known relationships, such as those described in the following reference: B. Hames and S. Higgins (eds.), Nucleic Acid Hybridization—A Practical Approach, IRL Press, Oxford, U.K., 1985.
The selected nucleic acid probes derived from the 16S rRNA using the in silico methods and provided herein include SEQ ID NOS: 1-48664. As described in the examples section, some of these probes show a better sensitivity and/or specificity than others, as reflected in the overall quality score for each probe (see Tables 14 and Table 20 which are attached to the specification as a separate document) which incorporates the estimated robustness of the probe and its predicted ability to detect most strains of a species. The better probes are therefore preferentially used in methods to detect the organism of interest in a biological sample. However, it is possible that for certain applications (e.g., epidemiology, substrain typing, and groups which comprise two or more species and potentially even spanning genera) a set of probes including the less specific and/or less sensitive probes may be very informative.
Preferably said probes share at least 70% or 80% identity with the complementary sequence. More preferably, the probes share at least 90% identity (e.g., 3 mismatches in a 30 base probe). More preferably the probes are at least about 95%, 96%, 97%, 98%, or 99% identital to the exact complement of the target sequence to be detected. Most preferred, the probe is homologous to the target sequence (i.e. has 100% sequence identity). Preferably, the target sequences are either ribosomal RNA or ribosomal DNA or amplified versions thereof. A probe with an apparent mismatch at one nucleotide position in its matching species rDNA sequence is still considered specific for the purposes of this invention because (1) available sequence databases contain occasional sequence errors (i.e., the apparent mismatch may not be an actual mismatch), and (2) in practice the hybridization conditions and the method determine the actual specificity of a probe. Therefore, it is important to maintain some flexibility in selecting probes for potential use, and to qualify them practically by testing them under actual hybridization conditions.
Preferably, these probes are about 5 to 60 nucleotides long, more preferably from about 10 to 30 nucleotides. The nucleotides as used in the present invention may be ribonucleotides, deoxyribonucleotides and modified nucleotides such as inosine or nucleotides containing modified groups, which do not essentially alter their hybridization characteristics. Also contemplated is the use of peptide nucleic acid probes (PNA). Moreover, it is obvious to one skilled in the art that any of the below-specified probes can be used as such, or in their complementary form, or in their RNA form (wherein T is replaced by U).
The probes according to the invention can be formed by cloning of recombinant plasmids containing inserts including the corresponding nucleotide sequences, if need be by cleaving the latter out from the cloned plasmids upon using the requisite nucleases and recovering them, e.g., by fractionation according to molecular weight. The probes according to the present invention can also be synthesized chemically, for instance by the conventional phosphotriester method.
The term “polynucleic acid” as used herein corresponds to either double-stranded or single-stranded DNA, rDNA, cDNA or genomic DNA or RNA, containing at least 5, 10, 15, 20, 25, 30, 40, 45 or 50 contiguous nucleotides or any length in between. A polynucleic acid that is smaller than 100 nucleotides in length is often also referred to as an oligonucleotide. Thus, in some instances polynucleotides greater than 100 bases are contemplated. Additionally the terms oligonucleotide and polynucleotide may be used interchangeably. Single stranded polynucleic acid sequences are always represented in the current invention from the 5′ end to the 3′ end.
The term “sensitivity” refers to the number of false negatives: i.e. if 1 of 100 strains to be detected is not detected, the test shows a sensitivity of [(100−1)/100]%=99%.
The term “specificity” refers to the number of false positives: i.e. if out of 100 strains detected, 2 actually belong to species for which the test is not designed, the specificity of the test is [(100−2)/100]%=98%.
The probes selected as being “preferential” show a sensitivity and specificity of more than 80%, preferably more than 90% and most preferably more than 95%. The higher the robustness number (i.e., the number of nucleotide changes required for a species-specific oligonucleotide sequence to match a second species), the more likely the probe will have a high specificity score. The precise relationship between robustness number and specificity depends on the method used to apply the species-specific oligonucleotide sequences of this invention. In general, probes with a robustness number of 2 or more are useful in hybridization-based methods that can distinguish 1-2 nucleotide differences such as the Affymetrix Gene Chip System. However, use of probes with higher robustness numbers increases the likelihood of achieving high specificity in a diagnostic test.
The term “primer” refers to a single stranded DNA oligonucleotide sequence capable of acting as a point of initiation for synthesis of a primer extension product which is complementary to the nucleic acid strand to be copied. The length and the sequence of the primer must be such that they allow to prime the synthesis of the extension products. Preferably the primer is about 5-50 nucleotides long. Specific length and sequence will depend on the complexity of the required DNA or RNA targets, as well as on the conditions of primer use such as temperature and ionic strength. The fact that amplification primers do not have to match exactly with the corresponding template sequence to warrant proper amplification is amply documented in the literature. Methods of amplifying products using primers are known in the art.
The oligonucleotides used as primers or probes may also comprise nucleotide analogues such as phosphorothioates, alkylphosphorothioates or peptide nucleic acids or may contain intercalating agents.
As most other variations or modifications introduced into the nucleic acid sequences of the invention, these variations may necessitate adaptions with respect to the conditions under which the oligonucleotide should be used to obtain the required specificity and sensitivity. However, the eventual results of hybridization will be essentially the same as those obtained with the unmodified oligonucleotides.
The introduction of these modifications may be advantageous in order to positively influence characteristics such as hybridization kinetics, reversibility of the hybrid-formation, biological stability of the oligonucleotide molecules, etc.
The term “solid support” can refer to any substrate to which an oligonucleotide probe can be coupled, provided that it retains its hybridization characteristics and provided that the background level of hybridization remains low. Usually the solid substrate will be a microtiter plate, a membrane (e.g. nylon or nitrocellulose), slide or microarray, or a microsphere (bead). Prior to application to the membrane or fixation it may be convenient to modify the nucleic acid probe in order to facilitate fixation or improve the hybridization efficiency. Such modifications may encompass homopolymer tailing, coupling with different reactive groups such as aliphatic groups, —NH2 groups, —SH groups, carboxylic groups, or coupling with biotin, haptens or proteins.
The term “labeled” refers to the use of labeled nucleic acids. Labeling may be carried out by the use of labeled nucleotides incorporated during the polymerase step of the amplification or by the use of labeled primers, or by any other method known to the person skilled in the art. The nature of the label may be isotopic (e.g., 32P, 35S, etc.) or non-isotopic (e.g., biotin, digoxigenin, etc.).
6. Preparation of Samples
In other embodiments of the invention, the sequences identified by these methods can be used in the form of probes and in kits for the identification of a specific family, genus, group, species or strain of a bacterium in a sample.
Typically, detection of the microorganism is performed by first isolating DNA from the sample. The sample can be, for example, from food, a mixed environmental sample or a biological sample from a subject (e.g., urine or blood). The cells from a sample are typically pelleted and resuspended in a lysis buffer (e.g., 10 mM Tris-HCl, 10 mM NaCl, 50 mM EDTA at pH 8.0). The cells are subjected to an enzymatic digest which can comprise, e.g., 25 U Lysostaphin, 30 μg N-Acetylmuramidase, 400 μg Achromopeptidase, and 600 μg lysozyme in a final volume of 306 μL of lysis buffer. The digestion is carried out for 30-45 minutes at 37° C. Lysates are extracted with chloroform/phenol, then ethanol precipitated. The DNA is redissolved in a Tris/EDTA buffer and the final concentration is determined by spectrophotometric measurement. Variations of this method would by known to the skilled artisan.
6.1 Amplification and Labeling of DNA. An amplification reaction is generally described in PCR P
Alternatively, rRNA from a sample can be used by first processing it using RT-PCR to produce DNA, or probes can be allowed to directly hybridize to the rRNA once it is isolated from the sample. Different labeling techniques for hybridization probe have been developed. In the past, labeling of nucleic acids by the enzymatic incorporation of radioactive isotopes (32P, 33P, 35S, 3H, and 125I) was the method of choice. However, in an effort to move away from the drawbacks associated to the use of radioactive products, DNA labeling techniques with non-radioactive products have been successfully developed. These non-isotopic labeling alternatives include among others the biotin-avidin/streptavidin system, the digoxigenin (DIG)/enzyme-labeled anti-DIG antibodies, the use of fluorescent dyes such as fluorescein, rhodamine, and cyanine dyes (Cy3 and Cy5). Nucleic acid labeling can be done by nick translation, random prime labeling, end-labeling reactions, and PCR. Several direct labeling kits are commercially available from different companies (Amersham Pharmacia Biotech Inc, Piscataway, N.J.; Molecular Probes Inc, Eugene, Oreg.; Vysis Inc, Downers Grove, Ill.) and make DNA labeling relatively easy. Labeled nucleic acid probes have been used in a wide range of applications in medicine, food industry, and environmental studies. For instance, probes labeled with fluorophore tags have been used in Fluorescent In Situ Hybridization (FISH) assay to detect bacteria in hemocultures (Oliveira et al., J Clin Microbiol. 2002, 40: 247-51).
To facilitate identification of a strain of a microorganism from a pattern of amplification products generated from the genome of that organism, the pattern should contain elements common to all members of the species in conjunction with elements which differentiate the unique strains within that species This only needs to be performed if it is necessary to distinguish strains. In the event that it is not necessary to distinguish between the strain of pathogen, this step can be omitted. This step is not necessary for the identification of many types of pathogens. In most cases, it will be preferable to include all strains in the species identification test.
This attribute makes it possible to identify a new strain of a species, which is not contained within the reference pattern database, providing that the common elements of that species pattern are contained within the reference database. If different strains of the same species each generate completely unique patterns with no common species elements, then the patterns themselves will not identify the species of an unknown strain. However, again typically there will not be a need to identify the strain. Additional methods of amplifying and labeling DNA are described infra and would be known to the artisan of ordinary skill.
6.2 Assays for Use with Plants Bacteria also cause problems in agriculture. For example, the bacterium, Acidovorax avenae, is a seed borne pathogen of several hosts including oats, corn, millet, wheat, sugarcane, rice and melons. It causes bacterial stripe of rice, leaf blight of oats, red stripe disease of sugarcane and millet and brown stripe of Setaria italica. It is an example of an organism which is hard to test for as it can be overgrown by other plant pathogens, such as the rice pathogens of Pantoea (Erwinia) herbicola, Burkholdera (Pseudomonas) glumae, B. fuscovaginae, or P. Syringae pv. syringae. Although these pathogens are known to produce distinct disease symptoms, field diagnosis remains difficult. Isolation of a pathogen is difficult in the presence of coexisting epiphytes. Thus, more sensitive testing methods are needed. Differentiation of microorganisms based on phyla, genus and species is important as well as testing of seed lots to determine whether seed batches are contaminated. Preferred bacteria include, but are not limited to, Pseudomonas fluorescens, Pseudomonas syringae species of the Burkholderia cepacia complex (Bcc), Erwinia species, Xanthomonas Species, Agrobacterium species and Clavibacter species.
DNA amplification and biological sample testing can be carried out generally as follows. In brief, the DNA amplification process is carried out by (a) providing a biological sample comprising bacterial cells or extracted DNA for standard PCR or cells amplified by growing on an agar medium for BIO-PCR; (b) amplifying a target sequence of the DNA to provide DNA amplification products carrying the selected target DNA sequence; and (c) detecting the presence of the DNA amplification products as an indication of the presence of, e.g., A. avenae.
The biological sample may either be bacteria cells or extracted genomic DNA. The biological sample may be a test sample suspected of containing bacterial cells, and thus the DNA of the bacterial cells, or a test sample containing extracted DNA.
The BIO-PCR method combines biological pre-amplification of the PCR target organism with enzymatic amplification of the PCR target. Briefly, the advantages of the BIO-PCR method, over those of the standard PCR assay, include the detection of live cells only, a 100-1000 fold increase in sensitivity, and elimination of PCR inhibitors associated with plant samples thereby eliminating false negatives. Sample processing can be further simplified by directly processing the samples comprising the expanded cells without further DNA extraction. However, even if a DNA extraction step is included, an advantage of the BIO-PCR methodology is that the DNA is extracted from a growing, viable population of cells or microorganisms. The enhanced sensitivity of the BIO-PCR method is particularly valuable, for example, in those screening situations where the monetary value of a particular seed type is high, and thus it is desirable to test the smallest quantity of seeds possible.
The preamplification step involves a plating step on an agar growth medium (or a liquid medium) prior to PCR analysis. A single cell per 0.1 ml can be detected because the single cell multiplies into a colony containing over 1000 cells on the agar medium.
Bacteria are recovered from suspect seeds of rice and watermelon by first soaking 1000-2000 seeds in 0.02% TWEEN 20 solution (ratio of 5 ml/g) for 4 hours at 40° C., pipetting aliquots of 0.1 ml of the bacterial extracts either onto plates of general purpose agar media, such as, KB or nutrient agar (NA), or onto a semi-selective medium such as EBB, described below, and then cloning. For BIO-PCR, each of five plates is washed three times with 1.0 ml of water and the resulting 14-15 ml of wash solution can be used for PCR-amplification with or without further DNA extraction or sample processing. Similarly, for standard PCR, either the DNA can be extracted or intact cells can be used. Since only pinpoint-size colonies are needed, incubation time ranges from only 10-15 hr for fast growing bacteria to 24-48 hrs for most plant pathogenic bacteria, depending on the media. Since the incubation time is short, few other bacterial colonies are present.
EBB medium is a preferred, semi-selective medium to use for the preamplification step. The protocol for preparing the EBB Medium is as follows:
In the preferred method, the enzymatic amplification of the DNA sequence is by polymerase chain reaction (PCR), as described herein.
For the binding and amplification, the biological sample (bacterial cells or extracted DNA) is provided in an aqueous buffer formulated with an effective amount of a divalent cation, which is preferable MgCl2, preferably at a concentration of about 0.05-5 mM; an effective amount of DNA polymerase with Taq DNA polymerase being preferred in the form of native purified enzyme or a synthesized form such as AMPLITAQ® (Perkin-Elmer), an effective amount of dNTPs as a nucleotide source, including, dATP, dCTP, dGTP, and dTTP, preferably in a saturating concentration, preferably about 200 FM per dNTP; and an effective amount of one or a pair of oligonucleotide primers. The reaction mixture containing the annealed primer(s) is reacted with a DNA polymerase in a thermocycler. Each PCR cycle begins with a DNA denaturation step of 94° C. for 30 s, followed by a primer annealing step at 60° C. for 30 s, and a DNA elongation step at 72° C. for 45 s.
If designed properly, a single product results. This product is preferably about 450-550 kb in size, whose termini are defined by the oligonucleotide primer(s), and whose length is defined by the distance between the two primers or the length of time of the amplification reaction. The gene sequence then serves as a template for the next amplification cycle.
The amplified DNA product is optionally separated from the reaction mixture and then analyzed. The amplified gene sequence may be visualized, for example, by electrophoresis in an agarose or polyacrylamide gel or by other like techniques, known and used in the art.
The amplified gene sequence may be directly or indirectly labeled by incorporation of an appropriate visualizing label, as for example, a radioactive, calorimetric, fluorometric or luminescent signal, or the like. In addition, the gel may be stained during or after electrophoresis with a visualizing dye such as ethidium bromide or silver stain, wherein the resulting bands by be visualized under ultraviolet light. The amplification of the gene sequence can be performed using PCR as described herein, or as would be known in the art.
To prove the identity of the amplified DNA product, a Southern blot assay should be conducted. The amplified products are separated by electrophoresis on a polyacrylamide or agarose gel, transferred to a membrane such as a nitrocellulose or nylon membrane, reacted with an oligonucleotide probe, and stained as above. The amplified products may also be detected by reverse blotting hybridization (dot blot) in which an oligonucleotide probe specific to the gene sequence is adhered to a nitrocellulose or polyvinylchloride (PVC) support, such as a multi-well plate, and then the sample containing labeled amplified product is added, reacted, washed to remove unbound substance, and a labeled amplified product attached to the probe or the gene sequence imaged by standard methods. A major advantage of TaqMan® PCR is that the technology is based on hybridization; therefore, a Southern blot assay is not needed.
The detection of amplified gene product in the sample is evidence of the presence of the bacterium, e.g., A. avenae subspecies, in the biological sample. When combined with BIO-PCR, the method is useful in diagnosing presence of viable cells of, e.g., A. avenae.
The primers and amplification method can further be useful for evaluating and monitoring the efficacy of any treatments utilized to eliminate the pathogen. In this method, biological samples are obtained from seeds, or other biological samples, prior to treatment and from seeds, or other biological samples that have undergone treatment with a treatment protocol. In addition, biological samples can be obtained from seeds at several time points during treatment. DNA amplification products of a target sequence of pathogen from all samples are analyzed for the presence of the pathogen. Results from samples obtained prior, during and after treatment are compared in order to determine efficacy of the treatment protocol.
7. Diagnostic, Prognostic and Classification Assays
It is also contemplated that the distinguishing family-, genus-, species- and/or strain-specific moieties be used in kits, such as a diagnostic kit. Additionally these probes can be used in kits helpful identifying a strain, which can be useful for classification or determining the appropriateness of a particular drug regimen to be administered to a subject. These kits may comprise any one or more of the following components:
Herein is provided a list of the genus-, group-, species- and strain-specific rRNA moieties thus identified (see Tables 14, 15, 20, and 21 which are attached to the specification as separate documents).
8. 16S Nucleic Acid Microarrays
8.1. Nucleic Acid Microarrays. The present invention provides for compositions comprising a plurality of polynucleotide probes which contain preferably two or more ribosomal moieties, e.g., SEQ ID NOS: 1-48664. The microarrays can be prepared for entire phylum, families, genera, species, groups, strains, or isolates of bacteria. Alternatively, diagnostic microarrays can be prepared based on infection phenotypes, e.g., all bacteria that produce skin lesions or particular types of skin lesions.
In one embodiment, the arrays can be prepared based on the origin of the biological sample, which would indicate only a certain grouping of organisms. Thus, microarrays can be created for screening body fluids, bronchial aspirates, cerebrospinal fluid, genital swabs, nares swabs, wounds, sputum, stool (e.g., for enteric pathogens, Yersinia, Aeromonas and Pleisomonas for example), throat samples, urine, miscellaneous sterile sites (e.g., surgical specimens), and blood culture (e.g., for aerobes, anaerobes, yeast, fungus, Mycobacteria). Arrays can also be created to screen cultures obtained from such biological samples. Preparations of such cultures and for obtaining such biological samples are known and routine in the art.
Biological samples can be obtained and screened based on clinical symptoms as determined by a health worker. For example, specific microarrays may be prepared based on inflammatory response, immune response or complement response, such as described in Table 3 below, which differentiate infectious diseases from conditions not linked to infection.
Microarrays can also be prepared based on means by which an organism is introduced into a host (e.g., bites, scratches, burns and environmental organisms), associated with surgery (e.g., prosthetic devices, infectious complications of solid organ transplantation) or result from suppressed immunity (e.g. AIDS).
Additional microarrays can be prepared for organisms common for causing sepsis, which results in 300,000 to 500,000 incidences a year, and about 100,000 deaths per year in the United States alone. Such microarrays would typically comprise nucleic acids as discussed herein from Gram-positive cocci (10-20% of the cases), Gram-negative, bacteria and fungi. Treatment of sepsis must be rapid and specific in order to eradicate the organisms from the blood stream and thus there is a continued need to develop faster and more accurate diagnostic screening methods.
Arrays can be prepared which are specific towards infectious diseases of the upper respiratory tract. These arrays would comprise nucleic acids identified by the methods described herein of organisms which cause nasal infections (commonly caused by Mycobacterium tuberculosis, M. leprae, Pseudomonus malleli, K. rhinoscleromatis, Rhinospordium seeberi, Actinomyces israelii, Cryptococcus neoformans, Blastomyces hominis), paranasal sinus infections (e.g., commonly caused by S. pneumoniae, Haemophilus influenazae, streptococci, and Moraxella), ear infections (commonly caused by S. pneumoniae, H. influenzae, Branhamella catarrhalis, Streptococcus and S. aureus). Similar arrays can be prepared to identify the microorganism causing an intraabdominal infection or abscess, acute infectious diarrheal disease and bacterial food poisoning (e.g., common infectious causing agents include Vibrio cholerae, E. coli, Clostridium perfringens, Bacillus cereus, Staphylococcus aureus, Aeromonas hydrophila, Plesiomonas shigelloides, Giardia lamblia, Cryptosporidum, Shigella spp., Salmonella enteritidis, Campylobacter jejuni, Vibrio parahemolyticus, Clostridium difficile, Entaemoeba histolytica, Salmonella typhi, Yersinia enterocolitica). Additional arrays may be prepared that are directed towards screening for microorganisms responsible for a sexually transmitted infection (e.g., commonly caused by Neisseria gonorrhoeae, Chlamydia trachomatis, Treponema palidum, Calymmatobacterium granulomatis, Ureaplasma urealyticum, Mycoplasma hominis, Gardnerella vaginalis, Shigella spp., and Campylobacter spp.), pelvic inflammatory disease, urinary tract infections (e.g., commonly caused by Gram-negative bacilli) and pyelonephritis. Additionally, arrays may be prepared to screen for the microorganism responsible for infectious arthritis in a host (such as Gram-postive cocci induced arthritis, gonococcal arthritis, chronic monarticular arthritis, viral arthritis, and spirochetal arthritis).
Alternatively, microarrays can be prepared based on, for example, bacterial classification. For example, arrays can be prepared for diseases that are caused by Gram-postive bacteria. Such infectious diseases include pneumococcal infections (e.g., caused by Streptococcus pneumoniae, and include such conditions as pneumococcal pneumonia, pneumococcal meningitis, pneumococcal peritonitis, and pneumococcal endocarditis), staphylococcal infections (e.g., commonly caused by S. saprophyticus, S. aureaus, S. epidermidis and include conditions such as staphylococcal scalded-skin syndrome, toxic shock syndrome, bacteremia, endocarditis, osteomyelitis, pneumonia and urinary tract infections), streptococcal infections (e.g., caused commonly by S. pyogenes, S. agalactiae, S. faecalis, S. equi, S. bovis, S. canis, S. mutans, S. anguis, S. milleri and include conditions such as pharyngitis, cellulitis, urinary tract infections, bacteremia, endocarditis, sinusitis, pneumonia, and meningitis), corynebacterial infections (C diphtheriae which causes diphtheria), anthrax (Bacillus anthracis), Listeria monocytogenes infections (e.g., contributes to conditions including sepsis, central nervous system infections, and endocarditis).
Also contemplated are microarrays, which distinguish diseases caused by Gram-negative bacteria. Gram-negative bacterial infections include, but are not limited to, meningococcal infections (e.g., Neisseria meningitides and N. lactamica), gonococcal infections (e.g., Neisseria gonorrhoeae), Moraxella (also known as Branhamella) infections (e.g., M. catarrhalis, M. osloensis, M. nonliquefaciens, M. osloensis and M. lacunata), Haemophilus infections (e.g., including H. influenzae, and Hacek group infections which include infections by H. phrophilis, H. paraphrophilis, H. parainfluenzae, Actinobacillus actinomycetemcomitans, Cardiobacterium hominis, Eikenella corrodens and Kingella kingae), and Legionella infections (e.g., Legionnaire's disease). Another embodiment contemplates microarrays that distinguish between diseases caused by Gram-negative enteric bacilli (e.g., E. coli, Klebsiella, Enterobacter, Serratia, Proteus, Morganella, Providencia and Acinetobacter).
Microarrays for distinguishing Pseudomonas species may also be prepared. The common Pseudomonas organisms that infect humans include P. aeruginosa, P. cepacia, Xanthoomonas maltophilia, P. pseudomallei, and P. mallei. These organisms are responsible for bacteremia, endocarditis, central nervous system infections, ear infections, eye infections, bone and joint infections, gastrointestinal infections, urinary tract infections, and skin/tissue infections.
Arrays for distinguishing Salmonella organisms are also contemplated, as these organisms are responsible for a large variety of human infections including typhoid (or enteric) fevers, focal systemic infections, septicemias, and gastroenteritis. S. typhi, S. paratyphi A, S. paratyphi B, S. typhimurium, and S. enteritidis are the most common species responsible for infection.
Also contemplated are arrays for distinguishing Shigella, which is responsible to Shigellosis and which is closely related to E. coli such that they cannot normally be distinguished by DNA hybridization methods. Common species involved in human disease include S. flexneri type 2A, S. dysenteriae, and S. sonnei.
Another embodiment contemplates an array that screens for mycobacterial diseases, which have become a problem, for example in patiesnts who are immune suppressed (e.g., AIDS patients). Mycobacteria are also responsible for leprosy (also known as Hansen's disease; M. leprae) and tuberculosis (caused by M. tuberculosis). Other Mycobacteria responsible for diseases include M. haemophilum, M. kansasii, M. marinum, M. scrofulaceum, M. szulgai, M. ulcerans, M. xenopi, M. asiaticum, and M. simiae. Also contemplated are microarrays that distinguish between the various species and strains of mycobacteria.
Microarrays can also be prepared which can distinguish between different spirochetes. Spirochetes are responsible for a wide range of infectious diseases including syphilis (caused by Treponema pallidum), leptospirosis (caused by Leptospira and which encompass Weil's disease and canicola fever), relapsing fever (caused by Borrelia), and Lyme borreliosis. In the instance of leptospirosis, Leptospira has only one species, L. interrogans, but has two complexes, interrogans and biflexa. The interrogans complex contains the pathogenic strains while the biflexa complex includes saphrophytic strains. Thus microarrays are preferred which can distinguish between the two complexes and potentially also the strains found in each of those complexes.
For Borrelia induced infections, there are several organisms responsible for the infections that are introduced either by louses or ticks. The most common Borrelia organisms include B. recurrentis, B. duttoni, B. hermsii, B. parkeri, and B. tuicatae. Borrelia are also responsible for Lyme disease (specifically B. burgdorferi). Thus preferred microarrays may contain nucleic acids for both Lyme disease and relapsing fever inducing microorganisms.
Microarrays can be prepared for use in distinguishing members of the genus of Chlamydia. Although classifed as bacteria, Chlamydia has its own order, Chlamydiales. The Chlamydia genus includes three species: C. psittaci, C. trachmatis and C. pneumoniae. C. psittaci can be found in numerous avian and mammalian species, but only the avian strains have been shown to infect humans causing a condition known as psittacosis. C. trachmatis is exclusively a human pathogen and is recognized to cause trachoma (a contagious form of conjunctivitis). It is also known to be one of the most common bacterial sexually transmitted diseases in the United States with an estimated 3 to 4 million cases annually. Such genital infections are responsible for urethritis, proctits, epididymitis in men, and mucopurulent cervicitis, acute salpingitis, bartholinitis and Fitz-Hugh-Curtis syndrome in women. C. pneumoniae is a fastidious species that appears to be a frequent cause of upper respiratory tract infections and pneumonia primarily in children and young adults.
Other organisms contemplated for use in such microarrays are provided in Tables 2 and 7. Microarrays using any of these organisms would preferentially be prepared based on any of the above-described compilations (e.g., distinguishing between strains and species, diagnostic screening based on origin of introduction, similar clinical presentation between infectious organisms, or similarity based on relative proximity on the phylogenic tree).
Preferably, the plurality of polynucleotide probes comprise at least a portion of one or more of the sequences of Tables 14, 15, 20, or 21 (attached to the specification as separate documents) or a fragment thereof.
A microarray can be used for large scale genetic or gene expression analysis of a large number of target polynucleotides. These microarrays can also be used in the diagnosis and/or prognosis of diseases and in the monitoring of treatments.
When the composition of the invention is employed as hybridizable array elements in a microarray, the array elements are organized in an ordered fashion so that each element is present at a specified location on the substrate. Because the array elements are at specified locations on the substrate, the hybridization patterns and intensities can be interpreted in terms of the presence or absence of particular nucleic acid sequences and can be correlated with, for example, a particular bacterial infection or for classification, a particular group or strain of bacterium.
The composition comprising a plurality of polynucleotide probes can also be used to purify a subpopulation of rRNAs or rDNAs, DNAs, genomic fragments and the like, in a sample. Typically, samples will include target polynucleotides of interest and other nucleic acids which may enhance the hybridization background. Therefore, it may be advantageous to remove these nucleic acids from the sample. One method for removing the additional nucleic acids is by hybridizing the sample containing target polynucleotides with immobilized polynucleotide probes under hybridizing conditions. Those nucleic acids that do not hybridize to the polynucleotide probes are washed away. At a later point, the immobilized target polynucleotide probes can be released in the form of purified target polynucleotides.
8.1.1. Method for Selecting Polynucleotide Probes. This section describes the selection of probe sequences for the plurality of polynucleotide probes. In one embodiment, the probe sequences are selected based on robustness and likelihood of being species-specific.
The resulting composition can comprise polynucleotide probes that are not redundant, i.e., there is no more than one polynucleotide probe to represent a particular distinguishing moiety for any particular rRNA and/or rDNA. Alternatively, and preferably, the composition can contain polynucleotide probes that are redundant, i.e., a rDNA is represented by more than one polynucleotide probe, because there are multiple combinations of distinguishing moieties known for that bacterium's rDNA.
The selected polynucleotide probes may be manipulated further to optimize their performance as hybridization probes. Probes that may not hybridize effectively under hybridization conditions due to secondary structure are avoided. To optimize probe selection, the sequences are examined using a computer algorithm to identify portions of genes without potential secondary structure. Such computer algorithms are well known in the art, such as Oligo® 4.06 software (National Biosciences) LASERGENE® software (DNASTAR®), and more preferably Primer 3 (Whitehead Institute and Howard Hughes Institute). These programs can search nucleotide sequences to identify stem-loop structures and tandem repeats and to analyze G+C content of the sequence; those sequences with a G+C content greater than 60% are preferably excluded from use in microarray compositions. Alternatively, the probes can be optimized by trial and error. Experiments can be performed to determine whether probes and complementary target polynucleotides hybridize optimally under experimental conditions.
8.1.2. Polynucleotide Probes. This section describes the polynucleotide probes for use in a nucleic acid microarray. The polynucleotide probes can be genomic DNA or amplified fragments of genomic DNA, which includes portions the bacterial genome responsible for producing rRNA, preferably rDNA, or rRNA, or DNA or cDNAs derived from rRNA, or complements of any thereof. The probes can also include peptide nucleic acids, branched DNAs and the like. The polynucleotide probes can be sense or antisense polynucleotide probes to the 16S rRNA or rDNA sequences of the subject bacterium. Where target polynucleotides are double stranded, the probes may be either sense or antisense strands. Where the target polynucleotides are single stranded, the nucleotide probes are complementary single strands. The preferred lengths of the probes will range between 12 and 60 nucleotides, more preferably between 12 and 50 nucleotides, more preferably between 15 and 35 nucleotides and most preferably between 20 and 30 nucleotides, and any range in between.
In another embodiment, the polynucleotide probes are plasmids. In this case, the size of the DNA sequence of interest, i.e., the insert sequence excluding the vector DNA and its regulatory sequences, may vary from about 15 to 2,000 nucleotides, more preferably from about 15 to 150 nucleotides.
The polynucleotide probes can be prepared by a variety of synthetic or enzymatic schemes that are well known in the art. The probes can be synthesized, in whole or in part, using chemical methods well known in the art. Alternatively, the probes can be generated, in whole or in part, enzymatically.
Nucleotide analogues can be incorporated into the polynucleotide probes by methods well known in the art. The only requirement is that the incorporated nucleotide analogues must serve to base pair with target polynucleotide sequences. For example, certain guanine nucleotides can be substituted with hypoxanthine which base pairs with cytosine residues. However, these base pairs are less stable than those between guanine and cytosine. Alternatively, adenine nucleotides can be substituted with 2,6-diaminopurine which can form stronger base pairs than those between adenine and thymidine.
Additionally, the polynucleotide probes can include nucleotides that have been derivatized chemically or enzymatically. Typical chemical modifications include derivatization with acyl, alkyl, aryl or amino groups.
The polynucleotide probes can be immobilized on a substrate. Preferred substrates are any suitable rigid or semi-rigid support including membranes, filters, chips, slides, wafers, fibers, magnetic or nonmagnetic beads, glass or plastic beads, gels, tubing, microtubes, plates, polymers, microparticles and capillaries. The substrate can have a variety of surface forms, such as wells, trenches, pins, channels and pores, to which the polynucleotide probes are bound. Preferably, the substrates are optically transparent.
Probes can be synthesized, in whole or in part, on the surface of a substrate using a chemical coupling procedure and a piezoelectric printing apparatus, such as that described in PCT publication WO95/251116 (Baldeschweiler et al.). Mechanical coupling is also an option. See U.S. Pat. No. 5,143,854 and PCT Application Nos. WO 90/15070 and WO 92/10092. Alternatively, the probe can be synthesized on a substrate surface using a self-addressable electronic device that controls when reagents are added (Heller et al. U.S. Pat. No. 5,605,662).
The arrays can have any density of oligonucleotide probes. Arrays can have probes of greater than about 100, preferably greater than about 1,000, more preferably greater than about 15,000, and most preferably greater than about 25,000 different oligonucleotide probes. Arrays with 50,000, 65,000, 200,000 or even 1,000,000 or more different probes are also contemplated. Such arrays generally comprise a probe density of typically greater than about 60, more generally greater than about 100, most generally greater than about 600, often greater than about 1,000, more oftren greater than about 5,000, preferably more than about 10,000 probes per cm. Some arrays may be preferred wherein the density is about 100,000 to 400,000 probes per cm2.
DNAs corresponding to rRNA or rDNA can be arranged and then immobilized on a substrate. The probes can be immobilized by covalent means such as by chemical bonding procedures or using ultraviolet (UV) light. In one such method, a DNA is bound to a glass surface, that has been modified to contain epoxide or aldehyde groups. In another case, a cDNA probe is placed on a polylysine coated surface and then UV cross-linked (Shalon et al., PCT publication WO95/35505). In yet another method, a DNA is actively transported from a solution to a given position on a substrate by electrical means (Heller et al. U.S. Pat. No. 5,605,662). Alternatively, individual DNA clones can be gridded on a filter. Cells are lysed, proteins and cellular components degraded, and the DNA coupled to the filter by UV cross-linking or other method known in the art.
Furthermore, the probes do not have to be directly bound to the substrate, but rather can be bound to the substrate through a linker group. The linker groups are typically about 6 to 50 atoms long to provide exposure to the attached polynucleotide probe. Preferred linker groups include ethylene glycol oligomers, diamines, diacids and the like. Reactive groups on the substrate surface react with one of the terminal portions of the linker to bind the linker to the substrate. The other terminal portion of the linker is then functionalized for binding the polynucleotide probe.
The polynucleotide probes can be attached to a substrate by dispensing reagents for probe synthesis on the substrate surface or by dispensing preformed DNA fragments or clones on the substrate surface. Typical dispensers include a micropipette delivering solution to the substrate with a robotic system to control the position of the micropipette with respect to the substrate. There can be a multiplicity of dispensers so that reagents can be delivered to the reaction regions simultaneously.
8.2 Hybridization Array Design
One of skill in the art will appreciate the enormous number of array options that are suitable for use with the disclosed methods and compositions. An array, including high density arrays, will typically include a number of probes that specifically hybridize to the target nucleic acid of interest. Such arrays should also preferably include one or more control probes. Such probes include but are not limited to test probes, normalization controls, mismatch controls, and sample preparation/amplification controls.
Test probes are oligonucleotide probes having sequences complementary to particular microorganism sequences whose expression they are designed to detect. Such test probes are typically from 5 to about 50 nucleotides, more preferably from about 10 to about 40 nucleotides and most preferably from about 15 to about 40 nucleotides in length.
Normalization controls are oligonucleotide probes that are perfectly complementary to labeled reference oligonucleotides that are added to the nucleic acid sample. The signals obtained from the normalization controls after hybridization provide a control for variations in hybrization conditions, label intensity, reading efficiency and other factors that may cause the signal of the hybridization to vary between arrays. In one embodiement, signals (e.g., fluorescence intensity) read from all other probes in the array are divided by the signal from the control probes thereby normalizing the measurements. Although any probe may be used as a normalization control, it is recognized that hybridization efficiencies vary with base composition and probe length. Thus, preferred normalization probes are selected to reflect the average length of the other probes present in the array. Normalization probes can be localized at any position in the array or at multiple positions throught the array to control for spatial variation in hybridization efficiency.
Mismatch controls can also be used in the array. Mismatch controls are oligonucleotide probes identical to their corresponding test or control probes except for the presence of one or more mismatched bases. A mismatched base is a base selected so that it is not complementary to the corresponding base in the target sequence to which the probe would otherwise specifically hybridize. One or more mismatches are selected such that under appropriate hybridization conditions the control probe would be expected to hybridize with its target sequence, but the mismatch probe would not hybridize (or would hybridize to statistically significant lesser extent). Preferred mismatch probes contain a central mismatch. For example, with a 20-mer probe, a corresponding mismatch probe will have the identical sequence except for a single base mismatch at any of positions 6 to 14.
Sample preparation and amplification controls are probes that are complementary to sequences of control genes or gene subsequences because they do not normally occur in the nucleic acids of the particular biological sample being analyzed. For example, the use of a eukaryotic gene in a sample believed to contain a microorganism. The sample is then spiked with a known amount of the nucleic acid to which the sample preparation/amplification control probe is directed before processing.
8.3 Labeling Nucleic Acids
Hybridized nucleic acids can be detected using one or more labels attached to the sample nucleic acids. The labels may be incorporated by any of a number of means known in the art. Preferred embodiments include incorporating the label during the amplification step in preparing the sample nucleic acids.
Alternatively, a label may be added directly to the original nucleic acid sample (e.g., mRNA, polyA mRNA, cDNA, etc.) or to the amplification product after the amplification is completed. Means of attaching labels are well known and include, but are not limited to, nick attachment (ligation) of a nucleic acid linker joining the sample nucleic acid to a label (e.g, a fluorophore).
Detection labels suitable for use include any composition detectable by means of spectroscopic, photochemical, biochemical, immunochemical, electrical, optical, or chemical means. Commonly preferred labels include biotin for staining with labeled strepavidin conjugates, magnetic beads (e.g., Dynabeads™), fluorescent dyes (e.g., fluorescein, Texas red, rhodamine, green fluorescent protein and the like), radiolabels (e.g., 3H, 125I, 32P, or 14C), enzymes (e.g., horseradish peroxidase, alkaline phosphatase and the like), and calorimetric labels such as colloidal gold.
8.4 Hybridization and Detection
Hybridization of the probes is preferably at low stringency. Hybridization stringency includes for example, 20° C. to about 50° C., more preferably 30° C. to about 40° C., and most preferably about 37° C. and 6×SSPE-T buffer (0.9 M NaCl, 60 mM Na2HPO4, 6 mM EDTA, 0.005% Triton X-100, pH 7.6) or less. Optimally washes can be performed thereafter at high stringency or progressively increasing stringency until a desired level of hybridization specificity is reached.
As the number of probes on a chip can be quite large, preferred chips may include only thoe probes that are needed, e.g, those probes relating to 16S or 23S. This will reduce the total number probes and the necessary size of the chip. This may also result in a different tiling strategy used for a particular chip design. For conceptual simplicity, probes are often arranged in order of the sequence in a lane across the chip. However, this tiling strategy is not required. Probes can be randomly placed on the chip. For discussion on tiling, see U.S. Pat. No. 6,228,575.
Probe length can vary. Probes usually have a single complementary segment having a length of at least 3 nucleotides, and more usually at least 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, or 30 bases exhibiting complimentarity.
In some chips, it may be optimal to have the probes all the same length, while in others it may be optimal to have probe sets of varying lenth. The probe length may vary by group or may vary by individual probe. For example, some chips may contain gour groups of probes having sizes of 20-mers, 22-mers, 25-mers and 30-mers. Other chips may have different size probes within the same group of four probes. Additional methods of designing chips comprising the distinguishing moieties or of the present invention would be known in the art. See e.g., U.S. Pat. No. 6,228,575.
Additionally, the arrays are also contemplated as having other components than probes, for example linkers attaching the probes to a support.
8.5 Microarrays and Kits for Their Use
Many of these microarrays can be used as research tools, as diagnostics to determine the presence of a microorganism in a sample (e.g., food sample or patient sample), as a means of following a patient's response to therapy and so forth. Thus these microarrays can be used in the form of kits. Kits can include an array of immobilized oligonucleotide probes complementary to subsequences of 16S and/or 23S. The kit may also include instructions describing the use of the array for the detection and/or quantification of expression levels of these sequences. The kit may optimally contain one or more of the following: buffers, hybridization mix, wash and read solutions, labels, labeling reagents (e.g., enzymes), control nucleic acids, and software for probe selection, array reading or data analysis.
9. Spot Tests and Other High Through-Put Testing
Potential high through put testing methods include but are not limited to (a) The Nanosphere Spot Assay (Nanosphere Inc.) that uses proprietary gold nanoparticle probe technology for a colorimetric detection of amplified DNA sequences; (b) a new fluorescence in situ hybridization (FISH) method with peptide nucleic acid (PNA) probes for identification of bacteria, which is based on a fluorescein-labeled PNA probe that targets a species-specific sequence of the 16S rRNA of the target species (e.g., S. aureus) (Oliveira et al., 2002 J. Clin. Microbiol. 40(1): 247-51); and c) a low-density microarray, in which a set of up to 1,000 probes are spotted on an array and tested with PCR products. Modifications of these would be readily apparent to the artisan of ordinary skill.
10. 16S Database Management System
Other aspects of the present invention involves a database management system for storing, maintaining, accessing and processing information relating to 16S rDNA sequences, where the information may be used to associate or distinguish moieties or residues from other moieties or residues based, for example, on nucleic acid sequence, location, relative position or any of a number of other characteristics.
Databases have been employed for storing, organizing and accessing biological data for quite some time. A database is a collection of information which has been organized in such a way that retrieval of the information stored therein is relatively quick and easy. In order to retrieve the information, however, a database management system is needed. A database management system is actually a collection of programs (i.e., computer programs) that enables one to enter, organize, select and access the information. For the purposes of the present invention, the term database and database management system are used interchangeably. Also, for purposes of the present invention, the information entered, organized, selected and accessed, as stated above, relates to 16S rDNA sequences, and in particular see Tables 14, 15, 19, 20, and 21.
Databases are traditionally organized into fields, records and files. A field is a space, e.g., a location in a memory device associated with one or more computer systems which make up the database management system, that has been allocated for a particular data item or piece of information. In general, a field is the smallest unit of information of data that one can access. A field may contain, for example, text information, such as a particular nucleic acid sequence; graphical information, such as data defining a three dimensional molecular representation of the nucleic acid sequence; or numeric information, such as the molecular weight of the nucleic acid sequence.
A collection of fields is generally referred to as a record. In the present invention, a single record may pertain to a particular full-length 16S rRNA/rDNA sequence or fragments thereof. The record may, in turn, contain several fields. As suggested above, one field may contain the nucleic acid sequence itself. Another field may contain the molecular weight of the sequence. Still another may contain data that defines a three dimensional molecular representation of the sequence. Yet another may define the classification of the microorganism (e.g., the bacterium) from which the 16S rRNA/rDNA sequence was isolated (i.e., the family, genus, species and/or strain of the microorganism). Still another may contain data that identifies the source of the sequence. In no way, however, is this list of examples intended to be exhaustive. In a relational database management system, records are often called tuples.
Generally speaking, a file is a collection of data or information that has been given a particular name, called a filename. In a database management system, a file may contain a collection of records. The records contained in that file may have a common attribute. For example, each record may contain 16S rRNA/rDNA data pertaining to or associated with a particular strain of bacteria.
The database management system of the present invention may be a relational database management system or, what is referred to as a flat-file database. A relational database is one that comprises multiple tables of data and/or information, where the word table simply refers to data arranged in rows and columns, and where the records contained in each table may have a different format. In a flat-file database, the data is contained in a single table. In accordance with exemplary embodiments of the present invention, the 16S rRNA/rDNA sequence data and/or information may be contained in a single table or spread across multiple tables. In one specific instance, one table may contain newly isolated 16S rRNA/rDNA sequences, or fragments thereof, where the sequence data associated with this table may be referred to as an internal database. Another table may contain 16S rRNA/rDNA sequences retrieved from one or more publically available or external databases, such as GenBank or the Ribosomal Database Project (RDP) database.
A more detailed description will now follow, with respect to generating, analyzing, annotating, storing, organizing and utilizing the 16S rRNA/rDNA sequence data and/or information in the database management system of the present invention.
10.1 Generating Raw Sequence Data. Raw sequence data refers to unedited nucleic acid sequence information. Raw sequence data may be obtained by sequencing isolated or amplified 16S rDNA or rRNA that has been subjected to PCR or RT-PCR and made into DNA or cDNA respectively. Raw sequence data may also be obtained through a variety of other methods, including the acquisition of sequence data from external sources. The preferred method for generating raw sequence data from a biological sample includes the steps of: isolating genomic DNA, preparing a template by PCR and, therefrom, generating a corresponding nucleic acid sequence.
10.2 Automated Bioanalysis. Once raw sequence data is generated for a given sample, as described in section 10.1 above, the raw sequence reads are preferably edited and annotated before the information is entered (i.e., stored) in the database. More specifically, the process maybe divided into two levels of processing: 1) editing raw sequence data and 2) annotating and organizing the edited sequences data. These processing steps are collective referred to herein as automated bioanalysis.
10.2.1 Editing. Automated bioanalysis processing of the 16S rRNA/rDNA sequence data prior to annotating and storing the data in the database according to exemplary embodiments of the present invention may involve, for example, sequence editing, sequence masking, clipping portions of sequences, and removing artifacts associated with cloning and sequencing, as would be apparent to one skilled in the art. Furthermore, the editing process may include the functional arrangement of sequences, for instance, through clustering and master clustering. Still further, the editing process may involve the use of existing sequences, that have been previously edited, to extend a sequence and identify other existing sequences related thereto.
10.2.2 Annotating and Organizing. After edited the raw sequences, the additional information related to a given sequence may be introduced by way of annotations. Annotations may include, as suggested above, information relating to the classification of the microorganism (e.g., the bacterium) from which the 16S rRNA/rDNA sequence was isolated (i.e., the family, genus, species and/or strain of the microorganism), the molecular weight of the sequence and the name(s) associated with the group responsible for isolating the sequence. Additionally, annotations may include information that relates a given sequence to other sequences within a given classification or classifications, disease/phenotype or disease relationship, pathology, histology or epidemiology relationships and methods used in editing the sequence.
The sequence information and the associated annotations pertaining to a given sequence are stored in the database. The annotations, assigned through the automated bioanalysis, may contain information on the bacterium from which the genes encoding the rRNA are expressed, the relationship of the bacterium to other members of the genus, species and/or strain, and preparation techniques. The sequences from rDNA libraries are preferably organized by bacterium category.
10.3 Database Organization. The database of the present system utilizes the capabilities of modern computers by storing genetic information in association with a large amount of related information. In a preferred embodiment, the information on essentially all the steps of isolating and amplifying DNA, and identifying rRNA or rDNA sequences comprising distinguishing moieties are stored in various relational tables. The database can also allow a user to access information pertinent to the sequences.
Both sequences and information annotating the sequences are stored in a database such as a relational database. Data is stored in the database in a functional arrangement that allows the user to store, track, and manipulate the rRNA or rDNA sequences and annotated information comprising information regarding distinguishing moieties. Users can access one or more databases via an integrated network, e.g. an Ethernet network. The workstations are typically computers, preferably personal computers, that include data entry means, output devices, display, CPU, memory (i.e., RAM and ROM) and interfaces to the network.
In the preferred embodiment of the present invention the relational database is stored at a file server connected to network. Computers are linked, via an integrated network, to a computer that grants access to the storage unit of the internal database of the present invention. The access computer preferably includes CPU, a memory means, interfaces to the network, and input and output devices. Reference databases illustrate sources of data that, for example, may be searched during use of the database. See for example, Table 7, which can be used as a reference database. Such a database is preferably in a relational form with a means, typically computerized, for querying the reference database. Tables 12, 14, 15, 20, and 21 comprise data that can be related within the database
10.4 Access to the Curated Database. The curated database preferably has a user-friendly interface, which is preferably created in HTML for access with Web browsers (e.g., Netscape® or Internet Explorer®).
10.5 Exemplary Full-length Sequences Stored in the Database. The sequences stored within the database provide information useful for designing probes to reveal the identity of various bacteria. Information of this nature is extremely powerful, as it can be utilized in clinical diagnostics, prognostics, patient treatment, etc. For example, certain types of bacterial species and/or strains may be resistant to some antibiotics, but may respond to other antibiotics. Determining this information would be clinically useful in prescribing the correct drug to treat the animal subject.
10.6 Use of the Internal Database. The structure and methods of data entry of the database allow many different types of analysis to be performed, both within the internal database and between sequences in the internal database and sequences in publicly available databases. The automated bioanalysis of the sequences enhances this analysis by masking or removing sequence elements that may hinder meaningful comparisons. The organization of the database facilitates analysis by providing mechanisms by which queries may be done quickly and efficiently, both within the internal database and with other external databases. The relational nature of the internal database thus provides a more comprehensive analysis, without the need to reformulate each search for each separate database.
10.6.1 Query Sequence Comparison. DNA sequence comparisons can involve comparing sequences within the internal database or comparing sequences with those in external databases.
Data relating to sequence comparison is organized and stored in the sequence comparison portion of the database. This storage area includes tables containing information about the quality of the sequence matches in sequence match logs, as well as tables containing information about other features of compared sequences. The sequence comparison portion also contains information found during accession of external databases (e.g., GenBank, RDP, PathoGenome™ (Genome Therapeutics Corp.), GenSeq (Derwent)). These databases may provide information on similarity, or rRNA domains, of the compared sequences that may be predictive of activity.
A sequence comparison results in a text file with details and summaries with respect to the possible relationships between the query sequence and any database sequences identified through the comparison. Preferably, comparisons are done with the aid of algorithms and software such as BLAST, FASTA, Boyer-Moore, or Smith-Waterman. Most preferably, parameters are used which permit the identification of highly related sequences with 0, 1, or more mismatches with the query sequence as discussed herein.
10.6.2 Species and Group Specific Oligonucleotide Generation. The full-length gene sequence can be split into various segment sizes with a given overlap length. These segments are then compared by, for example, BLAST against sequence databases. The databases that can be used for the pairwise comparisons of the oligonucleotides include GenBank, Ribosomal Database Project (RDP), PathoGenome™ and the internal reference database as described herein and in the Examples. The reports produced by the BLAST comparisons are then parsed to determine which oligonucleotides are either species or group-specific. The criterion for a species-specific oligonucleotide is that it exhibit only 0 or 1 mismatch when its sequence is aligned with at least one sequence from a single species (or no species if no additional sequences are available for the particular species under study) and exhibit two or more mismatches when its sequence is aligned with those from any other species. In other words, to be recognized as a “species-specific oligonucleotide” for a first species, an oligonucleotide hit to a second species must contain at least two mismatches. Likewise, a group-specific oligonucleotide is defined using the same mismatch criteria, but permitting 0 or 1 mismatch for five or fewer species. To better determine the validity of the species or groups-specific oligonucleotides, one can use multiple database comparisons because the content of different databases may not overlap entirely. The validity of the group- or species-specific oligonucleotide sequence is improved by comparison to a larger number of sequences from different databases.
Although the present invention has been described in detail with reference to examples below, it is understood that various modifications can be made without departing from the spirit of the invention, and would be readily known to the skilled artisan.
Direct sequencing of uncloned PCR generated template for variant sequence discovery is described by, for example, Wrischnik et al., 1987 Nuc. Acids Res. 15(2): 529-42; Gibbs et al., 1989 Proc. Nat'l Acad. Sci. USA 86(6): 1919-23; and Rogall et al., 1990 J. Gen. Microbiol. 136(Pt 9): 1915-20. This method has been applied to both prokaryotic and eukaryotic systems. Software to support this process is widely available (Nickerson et al., 1997 Nuc. Acids Res. 25(14): 2745-51). We have adapted this strategy for discovery of 16S RNA variation in numerous bacterial species.
PCR primers for amplification of the 16S rRNA gene were designed using the E. coli (ATCC11775) and S. aureus (ATCC12066) 16S rDNA sequences. Three tiers of amplicons were designed for complete coverage of the 16S gene as demonstrated in
Forward PCR primers are tailed with the M13-21 sequence and reverse primers with the M13-28 sequence to facilitate production sequencing. In addition, fourteen walking primers are available for use on the tier three amplicon. These primers are used to generate sequence reads directly from the 1,400 nucleotides (nt) fragment. All sequence reads are generated with BigDye Terminator chemistry (Applied Biosystems, Foster City, Calif.) run on MegaBACE 1000 instruments (Amersham, Piscataway, N.J.). Sequence reads representing the 16S sequence of individual bacterial isolates and passing Phred (Ewing et al., 1998, Genome Res. 8: 175-185; and Ewing et al., 1998 Genome Res. 8: 186-194) quality criteria of 175 total 30-quality bases, or 75 contiguous 30-quality bases, are assembled using polyphred. The assemblies are considered complete when a consensus sequence of at least 1,380 nt with an average of at least 3.5-fold Phred 20-quality data coverage is achieved.
A summary of all sequences that meet the above criteria was generated, along with information on the origin of the isolate from which each sequence was generated (the “source species”), and is represented in Table 7. Each strain, from which 16S rDNA sequence was obtained, was either a clinical isolate from the bioMérieux bacterial collection or was obtained from a reference collection center such as the ATCC. For all “type strains”, a “T” has been added after the strain collection number (e.g., ATCC 25238T).
DNA sequences of the 16S ribosomal loci of 1,214 bacterial samples representing 545 different species were generated and stored in an internal database. As outlined in the schematic diagram of
In order to identify conserved sequences upstream and downstream from the 30-mer species-specific oligonucleotide sequences, each 30-mer oligonucleotide was extended by 60 bases on each side to generate a 150-mer oligonucleotide. A multiple sequence alignment using CLUSTALW was generated using the particular 150-mer oligonucleotide to be tested and the 20 and 40 closest sequences in the GenBank database and in the internal database described above, respectively. This information was used to design primers for TaqMan®.
Superimposition of the species-specific regions on the structure of 16S revealed that they fall in non-base-paired loops with considerable tolerance for nucleotide differences. These hypervariable regions are known, but no systematic examination of 16S segments to identify those that are species specific has been reported.
DNA sequences of the 16S ribosomal loci of 1,324 bacterial samples representing 585 different species were processed in silico to yield 1,859,805 fragments of 20 nucleotide (nt) each, with an overlap of 19 nt (n=19). A pair-wise comparison was conducted using BLAST on the fragments against the Ribosomal Database Project (RDP) database. The BLAST reports generated from the comparisons were parsed using PERL programming language to identify oligonucleotides that matched a portion of the 16S rDNA sequence of a bacterial species with a criterion of permitting no more than 1 mismatch. Fragments that met this criterion and only matched one or zero species (or unidentified/uncultured/unknown entries) within that criterion were called “species-specific oligos”. The analysis discovered 90,079 fragments that met this criterion in RDP database. These 90,079 fragments were then blasted against the GenBank database to further qualify these oligonucleotides as being species-specific, i.e., to identify those which matched one or zero species (or unidentified/uncultured/unknown entries) by this criterion. Of these, only 37,072 remained species-specific after the comparison with GenBank. Of these species-specific oligonucleotides, 25,346 were unique and specific for 648 different species, including 68 out of the 81 priority species denoted in Table 1.
In order to identify conserved sequences upstream and downstream from the 20-mer species-specific oligonucleotide sequences, each 20-mer oligonucleotide was extended by 60 bases on each side to generate a 140-mer oligonucleotide. A multiple sequence alignment using CLUSTALW was generated using the particular 140-mer oligonucleotide to be tested and the 20 and 40 closest sequences in the GenBank database and the internal database, respectively.
Overlapping species specific oligonucleotides (SSO's) were converted into species specific regions (SSR's) as follows. Apparently species-specific oligonucleotide (20-mers) that were consecutively overlapping (e.g., nucleotides 1-20 and nucleotides 2-21, etc. of a 16S rDNA sequence) were combined into a large segment (i.e., 1-21 nucleotides) and termed a “species-specific region” (“SSR”). Any 20-mer oligonucleotide sequence selected from such an SSR will itself be species-specific for the particular species.
In order to distinguish strain- and species-specific oligonucleotides, a curated database was constructed from Ribosomal Database Project (RDP) 16S rDNA sequences. The RDP database consisted of aprroximatedly 30,000 16S rDNA sequences, but the majority of these sequences (i.e., ˜20,000) were not full length. Accordingly the RDP sequences were mapped to the same E. coli rDNA sequence as used to assign nucleotide numbers to the SSOs and the GSOs. Once the sequences were mapped, a comparison was performed to determine if each species-specific oligonucleotide was present at the same nucleotide coordinate positions in all entries for the species. The results indicated that 9,879 SSOs representing 232 species were species-specific using the prioritization ranking of 1-25. The other 18,191 SSOs were in the last ranking of 25-30, therefore being strain specific. The resulting curated RDP database (cRDP) was used to further predict the degree of species-specificity of probes.
A count of the number of entries for each species was conducted and the number stored in a table. Oligonucleotides that had been identified as being species-specific using the methods disclosed herein were compared using BLASTN2 (Washington University, St. Louis, Mo.) (default parameters except that “expectation value”=1000) against the sequences in this database. The results obtained using the BLAST algorithm were then parsed and examined for the number of species hit and compared with the entries for a given species represented. If a putative species-specific oligonucleotide hit all entries for a given species, then it was designated a species-specific oligonucleotide (SSO); but if the oligonucleotide did not hit all of the entries, then it was designated as a putative strain-specific oligonucleotide and placed on lower priority for further analysis. Species-specific oligonucleotides that are identified can then be used in, for example, diagnostic testing of bacterial infection or bacterial contamination. For example, if the number of E. coli entries in the database were five, then to be a species-specific oligonucleotide for the species E. coli, an oligonucleotide would have “hit” all 5 E. coli entries in the database by using BLAST or using a similar algorithm. Alternatively, if the oligonucleotide hit only 4 of the 5 entries, then the oligonucleotide would be a putative strain-specific oligonucleotide. Further analysis would be required to verify whether or not the oligonucleotide might be species-specific. For example, sequence errors in the database sequences in some cases could result in artificial designation of an oligonucleotide as a putative strain-specific oligonucleotide. 16S rDNA sequences would have to be obtained from additional strains of the species in question and then compared by BLAST or similar algorithm to the putative species-specific oligonucleotide in order to determine whether it is truly species-specific.
Species Specific Oligonucleotide “Robustness”. To prioritize the species-specific oligonucleotides according to their expected utility in a diagnostic test, BLAST reports from comparisons of the species-specific oligonucleotides with the database sequences were parsed, and a table was constructed indicating the number of mismatched nucleotides (i.e., distinguishing moieties) observed (up to a total of 5) in the BLAST alignment with 16S rDNA sequences from other bacterial species. The number of nucleotide changes required for a species-specific oligonucleotide to perfectly match the 16S rDNA sequence of another species is defined as its “robustness” in a diagnostic test. The larger the number of nucleotide changes required for a probe to perfectly align with a different species, the more robust the resulting diagnostic test. In other words, there is a greater margin of error or more tolerance in such a test. For example, hybridization conditions need not be as rigidly defined for an SSO with a high robustness number as for a species-specific oligonucleotide whose sequence is only 1 nucleotide change away from a perfect match with a second species.
Species Specific Oligonucleotide Quality Score. To further prioritize the species specific-oligonucleotides, a quality score was generated using information from both the “robustness” number and the predicted species-specificity. A scoring rating system was created based on 5 different robustness scores combined with 6 levels of predicted species-specificity. This generated a scoring rating system of scores ranging from 1 to 30, with 1 being the most preferred species-specific oligonucleotide and 30 being the least preferred. The 6 levels of predicted species-specificity are calculated based on the total number of entries and the total hits percentage, which is the percentage of the total full length cRDP database entries for a species which are “hit” (i.e., align, permitting no more than 1 mismatch) by the probe in question. The first specificity level represents oligonucleotides with greater than two full-length entries and a 100% total hits percentage. The second level represents oligonucleotides with only one full-length entry and a 100% total hits percentage. The third level contains oligonucleotides lacking any full-length entries; therefore, their species-specificity cannot be predicted at this time. Level four represents oligonucleotides with more than 10 entries and a total hits percentage of at least 85%. The fifth level contains entries from 4 to 9 and a total hits percentage of at least 75%. The sixth level represents oligonucleotides that are most likely strain-specific rather than species-specific and do not meet the above criteria. The “robustness” score indicates the number of nucleotide changes required for the probe to match another species sequence in either RDP or GenBank (i.e. MM2 is a mismatch of two nucleotides between the probe and the rDNA sequence of the next best species match, i.e., two nucleotide changes would be required for a perfect match to the next best species; MM>5 indicates that more than 5 nucleotide changes are required for the probe to match another species). In general, a probe with a robustness score of MM>5 is more preferred than a probe with a score of MM2. Therefore, a species-specific olignucleotide with a quality score of 1 would exhibit a “robustness” score of MM>5 and be at level one in predicted species-specificity. Table 8 summarizes the different quality scores according to “robustness” and predicted species-specificity, and Table 9 shows specific examples for each quality score. Table 14 provides the quality score information for all of the species-specific markers (attached to the specification as a separate document). Table 15 (attaches to the specification as a separate document) provides the number and identity of the species that each group-specific oligonucleotide marker identifies. This data was further refined a discussed in Example 17 below. All group-specific markers were chosen to match no more than five different species in either the GenBank or RDP databases, with the criterion that each marker exhibits no more than one mismatch when its sequence is aligned with the 16S rDNA sequence of each matching species. Table 15 lists the species that each group-specific oligonucleotide marker identifies. All group-specific markers were chosen to match no more than five different species in GenBank or RDP databases, with the criterion that each marker exhibits no more than one mismatch when its sequence is aligned with the 16S rDNA sequence of each matching species. Shown in Table 15 are the names of the species, which each group-specific marker matches as compared to 16S rDNA sequences in GenBank and RDP.
Enterococcus
Enterococcus
Enterococcus faecalis
faecalis
faecalis
Haemophilus
Haemophilus
Haemophilus ducreyi
ducreyi
ducreyi
Streptococcus
Streptococcus
Streptococcus
dysgalactiae
dysgalactiae
dysgalactiae
Streptococcus
Streptococcus
Streptococcus
thoraltensis
thoraltensis
thoraltensis
Prevotella disiens
Prevotella
Prevotella disiens
disiens
Atopobium
Atopobium
Atopobium parvulum
parvulum
parvulum
Corynebacterium
renale
Plesiomonas
Plesiomonas
Plesiomonas
shigelloides
shigelloides
shigelloides
Bacteroides
Bacteroides
Bacteroides stercoris
stercoris
stercoris
Streptococcus
Streptococcus
Streptococcus
dysgalactiae
dysgalactiae
dysgalactiae
Streptococcus
Streptococcus
Streptococcus suis
suis
suis
Streptococcus
Streptococcus
Streptococcus
dysgalactiae
dysgalactiae
dysgalactiae
Enterococcus
Enterococcus
Enterococcus
saccharolyticus
saccharolyticus
saccharolyticus
Streptococcus
Streptococcus
Streptococcus suis
suis
suis
Aerococcus
Aerococcus
Aerococcus urinae
urinae
urinae
Paenibacillus
Paenibacillus
Paenibacillus
azotofixans
azotofixans
azotofixans
Streptococcus
Streptococcus
Streptococcus
dysgalactiae
dysgalactiae
dysgalactiae
Campylobacter
Campylobacter
Campylobacter
sputorum
sputorum
sputorum
The system of naming species represented by species-specific oligonucleotides was done as follows. To further evaluate the species for which the oligonucleotide probe was specific, BLAST searches were run against GenBank and RDP for each oligonucleotide probe, and the name of the species exhibiting the best BLAST hit was compared to the species name of the isolate from which the oligonucleotide probe was derived, (i.e., the species name of this isolate in the internal database, SNID). In most cases, the same species name was found in the internal database, the RDP database, and the GenBank database. However, the species name for an rDNA sequence may vary in various databases for several reasons including, for example, systematic nomenclature changes as well as identification errors. In case of discrepancies in the species names between the three databases, the following criteria were used to select the species name for the purposes of this invention. The following definitions are helpful for describing the set of criteria: (1) let “SNIDseq” and “SNIDprobe” be the species name in the internal database (i.e., that assigned by diagnostic tests to the isolate from which the 16S rDNA sequence was determined and the probe was derived) for the 16S rDNA sequence and the probe, respectively, (2) let “BLASTseq” and “BLASTprobe” be the names of the species among the top 10 GenBank and RDP BLAST hits for the 16S rDNA sequence and the probe, respectively, and (3) let “bestBLASTseq” and “bestBLASTprobe” be the name of the species which is the best GenBank and RDP BLAST hits for the 16S rDNA sequence and the probe, respectively; and (4) let a single isolate be designated as “A” and multiple probes derived from the 16S rDNA sequence of isolate “A” be designated as “A1”, “A2”, etc. The following criteria serve as useful guides for assigning a species name to each probe, but one skilled in the art will recognize that confirmation must be obtained experimentally by use of the probes for hybridization-based testing. This naming criteria was also utilized in preparation of the attached Sequence Listing for the Source Identifier
TaqMan® probes and primers were designed from the 150-mer oligonucleotide sequence using PrimerExpress software (version 1.5, Applied Biosystems, Foster City, Calif.). Probes and primers were designed with the following characteristics (as recommended by Applied Biosystems):
The choice of TaqMan® probes can be done using PrimerExpress software or manually by the user; in which case, the user indicates to the software the TaqMan® probe sequence location, and the software searches for primers fitting the characteristics described above. If no primers corresponding to the criteria described above are found, the PrimerExpress software parameters are adjusted. For instance, the appropriate Tm may not be found surrounding some TaqMan® probe unless the requirements for Tm are adjusted (preferably, using a Tm between 53° C. and 60° C. and more than 2° C. difference between the Tms of the forward and reverse primers).
Using this method, twenty sets of probes and primers were designed for the following species: Bacteroides merdae, Bacteroides stercoris, Bacteroides thetaiotaomicron, Bacteroides vulgatus, Cardiobacterium hominis, Clostridium botulinum, Clostridium difficile, Clostridium septicum, Clostridium tetani, Eikenella corrodens, Enterococcus faecalis, Haemophilus ducreyi, Kingella denitrificans, Morganella morganii, Neisseria gonorrhoeae, Oligella urethralis, Proteus mirabilis, Providencia stuartii, and Streptococcus agalactiae.
Using sets of GSOs, rDNA sequences were identified to the species level as described further below.
Number of species covered by SSOs, GSOs or sets of GSOs. For species that were lacking a SSO (i.e., 76 species), the 20,594 GSOs were analyzed to determine if sets of GSOs could identify a particular species. The following table represents the 46 species in which species-specificity could be determined by using sets of GSOs.
Of the 567 species represented by the 1,324 completed 16S rDNA, 491 of these species had at least one SSO. Another 46 species were identified at the species level using a combination of two or more GSOs, as reflected in the third column of Table 10. Twenty species could only be characterized as belonging to a subset containing multiple species. Only 10 species were not able to be determined using the above analysis: Cedecea neteri, Citrobacter gillenii, Citrobacter murliniae, Enterobacter asburiae, Enterobacter intermedius, Kluyvera ascorbata, Kluyvera cryocrescens, Moraxella (Branhamella) catarrhalis, Pantoea dispersa, and Pseudomonas alcaligenes. Thus the probes utilized were able to discriminate 557 out of 567 species tested.
The species-specificity and robustness of probes was evaluated by bioinformatic analysis. Each 30-mer oligonucleotide was compared by BLASTN2 (Washington University, St. Louis, Mo.) (default parameters except that “expectation value”=1000) to sequences present in proprietary or public databases. The number of mismatches was derived from the BLAST report (see Example 4). As discussed, the robustness of the 30-mer increases as the number of mismatches increases. Table 11 below illustrates the number of mismatches needed for the probes to hit the next closest species in RDP database.
The species-specificity of probes was also verified at the experimental level using the TaqMan® method.
For each probe, one or several isolates from the target species and several isolates of related species were tested. The isolates tested in TaqMan® assays were obtained from the internal database described above. These isolates represented clinical strains obtained from different clinical laboratories and collection strains obtained from the American Type Culture Collection (ATCC) or other reference institutes. Identification of these strains was done by phenotypic methods, which included classical biochemical reference tests and/or commercial tests. The TaqMan® results, number of isolates and species tested for each probe are shown in Table 12.
*Explanation for isolates with late or no amplification:
Some of the discrepancies between TaqMan® results and the expected identification (i.e. the species the isolate belongs to) Were further investigated by sequencing of 16S rDNA. For instance, all isolates labeled as Bacteroides vulgatus were positive in a TaqMan® assay using a probe specific to B. vulgatus (see
TaqMan® reactions were conducted in 50 μL volumes in 96-well optical-grade plate (PE Applied Biosystems) using a ABI7700 sequence detector (PE Applied Biosystems). Each 50 μL reaction mix contained 900 nM forward primer (Seq ID 48716), 900 nM reverse primer (Seq ID 48741), 250 nM TaqMan® probe (Seq ID 273) labeled with FAM in the 5′ direction and TAMRA in the 3′ direction, and 25 μL of TaqMan® universal PCR master mix at a concentration of 2× (PE Applied Biosystems) that contained MgCl2, AmpErase® uracil-N-glycosylase (UNG), AmpliTaq Gold® DNA polymerase, reference dye ROX, and a mix of dNTP including dUTP, and 10 ng of DNA sample.
A typical TaqMan® reaction consisted of three steps: (1) 2 min. at 50° C. to activate the UNG enzyme that prevents the reamplification of carryover-PCR products by removing any uracyl present in double-stranded DNA; (2) 10 min at 95° C. to activate the AmpliTaq Gold® enzyme; and (3) 40 cycles of 15 sec at 95° C. (denaturation) and 1 min at 60° C. (annealing and extension).
Data analysis was conducted using the ABI 7900 Prism software (Applied Biosystems). The baseline fluorescence was determined by the software between the 3rd and 15th PCR cycles. The baseline can be recalculated between different cycles (3 and 12 for example) if the first positive samples have a Ct of less than 15. Once the baseline is set correctly, the software automatically sets the threshold at 10 standard deviations above the mean baseline fluorescence. The threshold can also be adjusted manually by the user. The cycle threshold (Ct) represents the cycle number at which there is a significant increase of fluorescence above threshold.
The software examines the fluorescence intensity of reporter and quencher dyes and calculates the increase in normalized reporter emission intensity over the course of the amplification. The ΔRn is calculated as follows:
A non-template control was used as a negative control.
TaqMan® reactions were conducted in 50 μL volumes in 96-well optical-grade plate (PE Applied Biosystems) using an ABI7700 sequence detector (PE Applied Biosystems). Each 50 μL reaction mix contained 900 nM forward primer (Seq ID 48720), 900 nM reverse primer (Seq ID 48745), 250 nM TaqMan® probe (Seq ID 5711) labeled with FAM in the 5′ direction and TAMRA in the 3′ direction, 25 μL of TaqMan® universal PCR master mix at a concentration of 2× (PE Applied Biosystems) that contained MgCl2, AmpErase® uracil-N-glycosylase (UNG), AmpliTaq Gold® DNA polymerase, reference dye ROX, and a mix of dNTP including dUTP, and 10 ng of DNA sample. The results expressed as ΔRn (indicator of the magnitude of the signal generated by the PCR) were then plotted versus time, represented by cycle number, to produce a continuous measure of PCR amplification.
A non-template (no DNA) control was used as a negative control. As shown in
The TaqMan® assay for detecting DNA sequence differences is based on real-time PCR which can be also be used to quantify the amount of target DNA originally present in a sample to test. Therefore, it is important to determine the optimal DNA concentration to use in each test, and different isolates should be tested at the same DNA concentration in order to make an accurate comparison of their threshold cycles (Ct). The threshold cycle (Ct) is defined as the PCR cycle at which a fluorescence signal passes a preset value (threshold). The concentration of DNA stock solutions was first determined using a spectrophotomer, reading DNA absorbance at 260 nm. Serial dilutions of genomic DNA from one Bacteroides stercoris isolate were used as templates for real-time PCR.
Each of 1,241 16S sequences has been subjected to a search for homologous sequences in the GenBank database using the BLAST algorithm. In order to determine if the phenotypic identification was supported by the 16S rDNA sequence information, the species represented by the first ten BLAST hits were compared to the species name assigned to the samples as a result of phenotypic characterization. The following categories: CORRECT TYPE, PROBABLY CORRECT, CLOSE MATCH, MISTYPED, POSSIBLE MISTYPE, NOVEL SEQUENCE, and VAGUE SAMPLE NAME were defined based on the presence of 16S rDNA sequences for the expected species in GenBank and/or RDP and the Blast hit position obtained for this given species (Table 16A).
BLAST result hits 4 to 10 were further investigated to determine whether or not a clear assignment can be made from the additional information obtained.
Only 69% of the isolates had a 16S sequence that gave the best sequence homology with the 16S sequence of the expected species. In other words, the identification of the isolate by 16S rDNA sequencing agreed with the identification results obtained by phenotyping for about two-thirds of all the isolates tested. For as many as 31% of the isolates, the species name defined by phenotypic methods did not match the information obtained using a BLAST analysis, at least at the level of the top three hits obtained.
There are several possible explanations for these findings including (1) the species not represented in public databases, (2) sequencing errors, (3) some 16S sequences are too short in the public databases, (4) nomenclature changes not integrated, (5) closely related species not differentiated by 16S rDNA, (6) species not claimed as identifiable with the phenotypic identification products used, (7) unclear identification based on phenotypic analysis, and/or (8) wrong phenotypic identification.
These results underscore the need to check a reference database for accuracy. These results also emphasize the need to develop rapid diagnostic tests based on molecular method in complement and/or in replacement of some phenotypic tests in order to improve the diagnostic of bacterial infections and orient towards an appropriate and efficient antibiotherapy.
The 16S rDNA probes described in this patent can serve as the basis for the design of research and commercial kits for bacterial identification. The probes are suitable to be tested in low-density microarray systems, such as those developed by Apibio (Grenoble, France). Several probes (up to 1,000) can be spotted on a solid substrate. The substrate or support actually contains several wells and several probes (up to 16) can be spotted in each well of the support. Depending on the number of probes specific to a particular species, one well can be dedicated to a particular bacterial species. Biotinylated PCR products will be prepared for each of the bacterial strains to be identified. The resulting PCR products will then be hybridized to the species-specific oligonucleotide (SSO) probes in each well of the solid support. The probe-target hybridization will be evaluated by a streptavidin-coupled detection with colorimetry or fluorescence. A positive signal is expected when the DNA target tested and the probes spotted on the support belong to the same species. Negative and positive controls preferably are both adhered or spotted on each DNA chip or chip substrate.
The 16S rDNA probes described in this application can be used on a high-density DNA chip system, such as GeneChip® (Affymetrix). Since the length of probes, the Tm and secondary structures of probes are important for preparing a GeneChip® array, the probes can be derived from 20-mer, 30-mer and species-specific regions (SSR), such as those disclosed herein. The species-specific oligonucleotides (SSOs) and group-specific oligonucleotides (GSOs) covering clinically relevant species were used to build a GeneChip® (Affymetrix) chip. One of these species was Streptococcus agalactiae, which is responsible for severe infections in neonates. For a given species-specific oligonucleotide or group-specific oligonucleotide, the original 20-mer sequence and four 17-mer probes derived from the 20-mer were included onto the chip. Each probe called a perfect match (PM) cell was associated to another probe called a mismatch (MM) cell that differs from the PM by just one nucleotide (probe substitution position). The pair of PM and MM cells defines an atom on the chip. For Streptococcus agalactiae species, twenty 20-mer species-specific oligonucleotides (SEQ ID NOS: 274, 393, 665, 703, 1076, 1173, 1409, 1508, 1576, 1614, 1838, 2199, 2213, 2541, 2839, 2857, 2914, 3078, 3157, and 3232) and thirty-two 20-mer group-specific oligonucleotides (SEQ ID NOS: 28304, 28306, 29325, 29448, 29632, 30338, 30633, 32212, 33065, 33324, 33370, 34524, 34870, 35688, 35949, 36303, 36583, 36647, 37293, 37551, 39226, 39998, 41902, 42485, 42884, 43386, 43488, 43620, 44569, 45263, 45776, 46945, 47112, and 47975) and their 17-mer derivatives were used on the chip.
The assay was conducted as follows. A Streptococcus agalactiae isolate was cultured on a blood agar plate overnight. After culture, cells (600 μl of a 0.5 McFarland suspension) were mechanically lysed with beads. Cell lysate (2 μL) was used to amplify 16S rDNA gene by PCR. Each 50 μL reaction mix contained 300 nM forward primer (5′-AGA GTT TGA TCA TGG CTC AG-3′), 300 nM reverse primer (5′-GAA GGA GGT GAT CCA ACC GCA-3′), 1×PCR buffer (Expand High Fidelity Buffer Roche) that contained 1.5 mM MgCl2, a mix of 200 nM dNTP, 1.75 unit of Taq DNA polymerase (Expand High Fidelity Taq DNA Polymerase, Roche), and 2 μL of DNA lysate. The 16S rDNA amplified fragment with a size of about 1,500 bp was then labeled with a biotinylated compound (bioPMDAM, bioMerieux), and cleaved with 0.1N HCl in order to obtain small fragments (i.e., less than 100 bp) to hybridize with the probes on the chip.
Labeled and cleaved PCR samples were then hybridized on the Affymetrix chip at 45° C. for 30 min in a hybridization buffer containing 7× saline-sodium phosphate-EDTA (SSPE), 3 M betaine, 0.01 M dodecyltri-methylammonium bromide (DTAB), and antifoam. Hybridized chips were then washed in 4×SSPE buffer containing 0.01% Triton, followed by washes in 6×SSPE buffer with 0.01% Triton. The hybridized chip was then stained with Streptavidine-phycoerythrin (DAKO) and laser scanned at 570 nm on a scanner (Hewlet Packard). The scan was recorded as a pixel image.
Data were analyzed using a bioMerieux software similar to GeneChip Analysis Suite, version 3.3 (Affymetrix). Average background noise was substracted from all the cells. A probe was scored as positive when 1) the intensity of the hybridization signal of the perfectly matched probe cell (PM) was greater than or equal to 1.2 times the intensity of the mismatch control cell (MM); and 2) the intensity of the PM was 8 times higher than the noise value. Data analysis was based on the comparaison of the hybridization signals. The full length 16S rDNA sequences of the species represented on the chip are stored into a reference database. Since each probe on the chip is derived from one 16S rDNA sequence of this database, it can easily be matched to the relevant 16S rDNA sequence and species.
In the present example, the reference database used contained 40 full-length 16S rDNAs sequences from 20 species (i.e., Burkholderia cepacia, Burkholderia gladioli, Burkholderia multivorans, Cardiobacterium hominis, Enterobacter sakazakii, Enterococcus faecalis, Gemella morbillorum, Granulicatella adiacens, Kocuria rosea, Mannheimia haemolytica, Pantoea agglomerans, Proteus penneri, Providencia stuartii, Pseudomonas aeruginosa, Ralstonia pickettii, Sphingobacterium multivorum, Staphylococcus chromogenes, Staphylococcus simulans, Streptococcus agalactiae, and Vibrio cholerae). All the positive hybridization signals of the test samples are then correlated to species that contain the probes on the chip with the highest sequence homology. These sequences are full length 16S rDNA sequences without SEQ ID NO and are included in Table 16B. Each species is represented on the chip by its different species-specific and group-specific oligonucleotides. Since some species have a higher number of species-specific oligonucleotides than others and some species lack group-specific oligonucleotides, the number of nucleotide positions tested on the chip varies from one species to species. Below is the percent homology ranking obtained for the type strain of Streptococcus agalactiae (ATCC 13813T).
*“SSO”: species-specific oligonucleotide; “GSO”: group-specific oligonucleotide
Species-specific probes and group-specific probes were analyzed separately as well as in combination. In all cases, the best homology score was obtained for Streptococcus agalactiae. These results confirm that the 52 probes cited above and specific to Streptococcus agalactiae or to a group of species including S. agalactiae can be used on a microarray for the detection of Streptococcus agalactiae in a sample. This methodology can be applied to a large range of species. Such microarrays can be used to diagnose the cause of a bacterial infection.
The LightCycler method is a real-time PCR method. It uses two different species-specific oligonucleotide (SSO) probes that are labeled with different dyes. These two hybridization probes hybridize to their complementary DNA template in a head-to-tail arrangement; i.e., the two probes are just separated by one or two nucleotides. Because of this close proximity and after excitation by a blue light source, the donor dye (fluorescein) located on the first probe transfers its energy to the acceptor dye (Red640) located on the second probe. The transfer of energy causes fluorescent light to be emitted. When the probes cannot find their complementary sequence on the DNA template, no hybridization occurs such that no fluorescent signal is produced. The intensity of the fluorescent signal is proportional to the amount of PCR product generated during the cycle. Similar to the TaqMan® assay, a cycle threshold (Ct) is obtained.
In Example 13, the 16S rDNA sequence of Staphylococus simulans strains was used to select two species-specific hybridiztion probes spaced by 1 nucleotide. At the same time, two non-species-specific PCR primers were designed as controls and used in a LightCycler assay with the two hybridization probes. Sequences of these probes and primers are presented in Table 17.
Real-time PCR data obtained with this combination of probes and primers on different Staphylococcal species are presented in
The pyrosequencing method is a simple and easy-to-use sequencing method by synthesis. This method uses a mixture of four enzymes and specific substrates. At first, a sequencing primer is hybridized to a single-stranded PCR product (e.g., DNA template). With an enzyme cascade system, light is produced whenever a nucleotide is incorporated and forms a base pair with a complementary nucleotide located on the DNA template (as discussed in Example 13 above). Each peak represents the incorporation of the nucleotide mentioned on the X-axis during the sequencing reaction. The size of the peak determines the number of nucleotides incorporated. The software also gives a schematic representation of the pyrogram for each possibility of the single nucleotide polymorphism (SNP).
In this example, six species of Enterococcus (i.e., Enterococccus faecalis, E. durans, E. faecium, E. avium, E. gallinarum, and E. saccharolyticus) were selected to determine if pyrosequencing was a suitable method to validate (at the experimental level) the species specificity of some SSOs (signature sequences). A multiple sequence alignment of 16S rDNA sequences of the six different species of Enteroccci was performed using ClustaIX and GeneDoc softwares. Variable regions that contained SSOs were identified. Variation at one nucleotide position between the different species defines a distinguishing moiety also known as a single nucleotide polymorphism (SNP). Four SNP positions were selected to test the pyrosequencing method.
PCR primers and sequencing primers were designed to amplify and sequence the 16S rDNA areas that contain a SNP. The size of PCR product should not exceed 500 bp in order to minimize the risk of mispriming for the sequencing primer. The pyrosequencing primer needs to be localized just upstream of the SNP site as the maximum length of the sequence obtained by pyrosequencing is about 40-50 bases. Pyrosequencing was done on a Pyrosequencer using the SNP detection software.
Two isolates per species for the six different species of Enterococci were tested. An example of pyrogram is presented in
Interestingly, by testing only 4 SNPs, it was possible to differentiate six species of Enterococci. These results underscore the high potential of pyrosequencing for bacterial typing and rapid evaluation of signature sequences, provided that adequate PCR and sequencing primers are designed and used for the group of species to screen.
To identify nucleotides that are critical for distinguishing species (i.e., distinguishing moieties), comparisons with the GenBank sequence database were made using BLAST for each of the 28,070 Species Specific Oligonucleotides. For each SSO, the top GenBank BLAST hits were compared with the SSO sequence, and the number of variations at each nucleotide position counted. The higher the count at any position of the 30-mer (e.g., Position 1 being represented by the column entitled “pos—01”; Postion 2 being represented in the column entitled “pos—02” and so forth), the more critical th distinguishing moiety is in discriminating other species. These counts are represented in Table 19, which is attached to the specification as a separate document. For certain sequences, only a 20-mer was used and analyzed, consequently no information is provided for certain sequences at postions 21-30.
The results acquired through this method may change over time due to the entry of new sequence information into the public domain. Thus, one skilled in the art will recognize that the assignment of robustness of distinguishing moieties is a process that should be updated periodically by methods such as the one described here. Even if the status of some of the distinguishing moieties described here changes from species-specific to genus-specific or group-specific, a sufficient number of probes containing distinguishing moieties are provided to maintain the diagnostic utility of the invention.
The data from Example 4 was analyzed further in view of additional sequence information released to the public. Sequences from a total of 567 different bacterial species were examined. No sequence could be obtained for 7 isolates representing 4 species: Corynebacterium minutissimum (2 isolates), Kocuria kristinae (1 isolate), Micrococcus lylae (3 isolates), and Sphingobacterium spiritivorum (1 isolate). The 16 S loci of these 7 isolates could not be amplified using the 16S universial primers as discussed above.
From the collection of 1,325 16S rDNA sequences, 2,724 30-mer SSOs, 25,346 20-mer SSOs, and 20,954 20-mer GSOs were identified in silico using the methods discussed above. The 30-mer and 20-mer SSos were specific for 325 and 648 different species respectively. At least one SSO was identified for each of 491 of the 567 total species examined. Another 46 species were identified as the species level using a combination of two or more GSOs. Tweny species could only be characterized by using GSOs as belonging to a subset containing multiple sepces. Only 10 species were not identifiable at all using any combination of the SSOs and GSOs (i.e., Cedecea neteri, Citrobacter gillenii, Citrobacter murliniae, Enterobacter asburiae, Enterobacter intermedius, Kluvera ascorbata, Kluyvera cryocrescens, Moraxella (Branbamella) catarrhalis, Pantoea dispersa, and Pseudomonas alcaligenes.). The SSO and GSO probes developed were capable of discriminating 557 out of the 567 species.
To ensure the elimination of SSOs that were merely strain specific, further comparisons using all sequences for a given species in both public and properiteary 16S databases was performed. A total of 15,127 SSOs were found to be species-specific for 271 species after comparison to the proprietary sequence data and to publicaly available 16S rDNA sequences in GenBank. When the analysis was extended to include partial 16S rDNA sequences in RDP, a total of 9,879 SSOs representing 232 species were species specific. The information based on this further analysis is provided in Tables 20 and 21, which are similar to Tables 14 and 15 respectively.
All references cited herein as well as U.S. Ser. No. 60/464,955 filed Apr. 24, 2003 are herein incorporated in their entirety for all purposes.
Number | Date | Country | |
---|---|---|---|
60464955 | Apr 2003 | US |