Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2005, p. 2092-2103, Vol. 43, No. 5
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.5.2092-2103.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Steven F. Hurst,
Timothy J. Lott,
David W. Warnock, and
Christine J. Morrison*
Mycotic Diseases Branch, Division of Bacterial and Mycotic Diseases, National Center for Infectious Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia
Received 25 October 2004/ Accepted 27 December 2004
|
|
|---|
1-nucleotide divergence in the D1-D2 region from closely related but different species. In contrast, only 5 of the species examined exhibited a
1-nucleotide divergence from sibling species in their ITS1 or ITS2 sequences. Although the GenBank database currently lacks ITS sequence entries for some species, and major improvement in the quality and accuracy of GenBank entries is needed, current identification of medically important Aspergillus species using GenBank reference data seems more reliable using ITS query sequences than D1-D2 sequences, especially for the identification of closely related species. |
|
|---|
The current laboratory identification of Aspergillus species is based on macroscopic colonial and microscopic morphological characteristics (7, 20, 45). Over 180 different species in at least 16 subgeneric groups or sections can be distinguished (35, 37, 38), including approximately 30 species which are recognized as opportunistic pathogens of humans (7). Many clinical laboratories use traditional phenotypic methods of identification and can differentiate only the more common Aspergillus species; the delineation of less common species must be referred to specialist laboratories. In addition, species identification by traditional phenotypic methods may require several weeks to accomplish, and in a small number of cases, isolates may not conidiate, obstructing species identification (20, 38). Therefore, rapid molecular approaches are now being developed to replace the need for culture by detecting and identifying Aspergillus species DNA directly from clinical materials.
A number of targets for the molecular identification of aspergilli have been investigated including the mitochondrial cytochrome b gene (46, 57, 62), a putative aflatoxin pathway regulatory gene (aflR) (4), the DNA topoisomerase II gene (TOP2) (19), the ß-tubulin gene (11), and various rRNA gene regions (18). The most promising targets to date have been the 5' end of the large-subunit rRNA gene (D1-D2 region) (35) and the internal transcribed spacer 1 and 2 (ITS1 and ITS2) regions between the small- and large-subunit rRNA genes (18, 60). The use of DNA sequence diversity in the ribosomal regions as an aid to species identification has been exploited using PCR amplification of targets followed by either fragment length analysis (1, 16, 22, 39, 52, 56), DNA probe hybridization (15, 30, 32, 61), or DNA sequence analysis (3, 17, 35, 43, 55, 61). In general, DNA-based approaches were found to provide more reliable and faster species identifications than culture-based methods. Nonetheless, some diagnostic problems were encountered, resulting from similar amplicon or fragment lengths for different taxa (16, 39, 52), unexpected cross-hybridization results (9, 16, 30, 32, 42, 61), and ambiguous sequencing-based identifications because of significant intraspecies heterogeneity (44, 53, 54) or virtually identical sequences for apparently distinct organisms (14, 17, 35, 43). The validation of test specificity was further compromised by the absence of a comprehensive test strain panel; most studies did not include type strains or authenticated culture collection reference strains of the target species or of closely related species in their examinations. Finally, although the D1-D2 and the ITS1 and ITS2 regions have been analyzed separately for some species of Aspergillus, a systematic evaluation and comparison of sequences from all three ribosomal regions for their usefulness in the identification and differentiation of the most medically important Aspergillus species has not been published to date.
Therefore, the present investigation assessed the utility of the ribosomal D1-D2, ITS1, and ITS2 regions as targets for the molecular identification of 13 potentially invasive Aspergillus species (Aspergillus candidus, Aspergillus chevalieri, Aspergillus flavipes, Aspergillus flavus, Aspergillus fumigatus, Aspergillus granulosus, Aspergillus nidulans, Aspergillus niger, Aspergillus restrictus, Aspergillus sydowii, Aspergillus terreus, Aspergillus ustus, and Aspergillus versicolor). We acquired DNA sequence information concerning all three ribosomal regions for each of the above species and used this information to conduct DNA sequence alignments, pairwise nucleotide sequence analyses, comparative GenBank database searches, and an examination for sequence length polymorphisms among Aspergillus species.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Ribosomal Aspergillus sequences determined in this study
|
PCR amplification of the D1-D2 region of the large-subunit (28S) rRNA gene. Seminested PCR using broad-range primer pairs ITS1 (5' TCC GTA GGT GAA CCT GCG G) and D2R (5' TTG GTC CGT GTT TCA AGA CG) followed by D1 (5' GCA TAT CAA TAA GCG GAG GA) and D2R generated the D1-D2 region amplicons for sequencing (36, 60). PCR amplification was conducted in a GeneAmp model 9700 thermal cycler (Perkin Elmer Applied Biosystems, Foster City, Calif.), and PCR components were obtained from Roche Molecular Biochemicals, Indianapolis, Ind. The PCR mix contained 0.2 mM concentrations of each deoxynucleoside triphosphate, 1.25 U of Taq DNA polymerase, 0.2 µM concentrations of each primer, 2 µl of a 1:50 dilution of the original DNA extract (>5 ng of DNA for the initial PCR or 2 µl of undiluted amplicon from the initial PCR when products were reamplified in the seminested PCR), and the appropriate amount of PCR buffer (final concentration, 10 mM Tris-HCl [pH 8.3], 1.5 mM MgCl2, 50 mM KCl) to bring the final volume to 50 µl. One sample without DNA template was always included as a negative control. Amplification with the primer pair ITS1 and D2R was performed using a denaturation step of 95°C for 5 min, followed by 30 cycles of 95°C for 30 s, 51°C for 30 s, and 72°C for 150 s. A final extension step at 72°C for 10 min was then used. Reamplification with primer pair D1 and D2R for the seminested PCR amplification included denaturation at 95°C for 5 min, followed by 30 cycles of 95°C for 30 s, 60°C for 45 s, and 72°C for 90 s. A final extension step at 72°C for 7 min was then employed.
PCR amplification of the ribosomal ITS1 and ITS2 regions. The broad-range primer pairs ITS5 (5' GGA AGT AAA AGT CGT AAC AAG G) and ITS2 (5' GCT GCG TTC TTC ATC GAT GC) or ITS3 (5' GCA TCG ATG AAG AAC GCA GC) and ITS4 (5' TCC TCC GCT TAT TGA TAT GC) were used to amplify ribosomal internal transcribed spacer regions 1 and 2, respectively (60). Some samples required seminested PCR amplification with primer pair ITS5 and ITS4 followed by ITS5 and ITS2 or ITS3 and ITS4 to generate sufficient quantities of PCR amplicon for DNA sequencing. Amplifications and reamplifications were performed as described above, except for the primer pair used and minor differences in the thermal cycling conditions employed: initial denaturation at 95°C for 5 min, followed by 30 cycles of 95°C for 30 s, 58°C for 30 s, and 72°C for 1 min. A final extension step at 72°C for 5 min was then conducted.
DNA sequencing. All PCR products were purified before DNA sequence analysis using a QIAquick PCR purification kit (Qiagen) according to the manufacturer's instructions. Purified amplicons were then sequenced on both strands using the same primers as described above. BigDye terminator cycle sequencing Ready Reaction kits (Perkin Elmer Applied Biosystems) were employed as recommended by the manufacturer. All cycle sequencing reactions were performed on a GeneAmp model 9700 thermal cycler using an initial denaturation at 96°C for 5 s, followed by 30 cycles of 96°C for 10 s, 50°C for 5 s, and 60°C for 4 min. Products were purified using a Dye-Ex spin kit (Qiagen), dried in a vacuum centrifuge, and resuspended in either template suppression reagent (Perkin Elmer Applied Biosystems; D1/D2R products) or formamide (Perkin Elmer Applied Biosystems; ITS products). Products were then analyzed on an automated capillary DNA sequencer (ABI Prism 310 genetic analyzer; Perkin Elmer Applied Biosystems) according to the manufacturer's directions. Comparative sequence analysis and GenBank searches were assisted by the Genetics Computer Group software package (FASTA, BESTFIT, STRETCHER, and PILEUP algorithms; University of Wisconsin, Madison), the Clustal W alignment program (51), and the nucleotide-nucleotide Basic Local Alignment Search Tool (BLAST) algorithm (blastn) (2, 63). Calculation of pairwise sequence identity was restricted to full-length reference data of Aspergillus species currently sanctioned by the International Commission of Penicillium and Aspergillus (37). Authoritative Aspergillus reference sequences were taken from the work of Peterson (35), Haugland et al. (15), and Pazoutova et al. (34).
Nucleotide sequence accession numbers. The Aspergillus large-subunit D1-D2, ITS1, and ITS2 ribosomal regions determined in this study were deposited in GenBank and assigned the accession numbers listed in Table 1. The species designations are in compliance with currently sanctioned names published by the International Commission of Penicillium and Aspergillus (37).
|
|
|---|
Use of the D1-D2, ITS1, and ITS2 DNA sequence length to discriminate among Aspergillus species investigated in this study. The D1-D2 sequencing results for all Aspergillus species studied showed identical overall lengths for all species and strains investigated (i.e., 542 nucleotides [nt]) with the exception of A. nidulans (total length, 543 nt [Table 1]). In contrast, the ITS1 region ranged in overall length from 142 nt (A. chevalieri) to 186 nt (A. terreus), and sequence length was more divergent among Aspergillus species in the ITS1 region than in the D1-D2 region; i.e., in the ITS1 region, 5 of the 13 species studied demonstrated a unique sequence length (A. candidus, A. chevalieri, A. granulosus, A. nidulans, and A. terreus) compared to only one species (A. nidulans) in the D1-D2 region. In addition, the overall sequence length for the type strain of A. ustus (NRRL 275) not only was unique compared to other Aspergillus species but differed by 6 nt from its own sequevars (i.e., 161 nt versus 155 nt); these data suggest substantial divergence of these sequevars from the type strain. However, all of the 7 remaining species and the two divergent A. ustus sequevars (I and II) shared identical overall sequence length with at least one strain from each of the 13 species examined (e.g., A. flavipes and A. niger, 185 nt [Table 1]). Identical or similar overall sequence lengths were observed for A. candidus, A. flavus, and A. restrictus (180 to 181 nt), A. flavipes, A. fumigatus, A. niger, and A. terreus (184 to 186 nt), and A. granulosus, A. nidulans, A. sydowii, A. ustus (sequevars I and II), and A. versicolor (153 to 156 nt). In some instances, overall sequence length differences varied more within a given species than between species and ranged from 1 nt (A. flavipes, A. nidulans) to 6 nt (A. ustus) (Table 1).
In comparison, ITS2 regions varied in overall length from 160 nt (A. flavipes) to 177 nt (A. terreus) (Table 1). Overall sequence length was slightly less divergent among the Aspergillus species in the ITS2 than in the ITS1 region but was more divergent than in the D1-D2 region (Table 1). Three of the 13 species examined demonstrated a unique overall sequence length in the ITS2 region (A. chevalieri, A. granulosus, and A. terreus), and 5 shared identical overall sequence length with at least 1 strain from each of the 13 species examined (e.g., A. chevalieri and A. versicolor, 167 nt). Identical or similar sequence lengths were observed for A. candidus, A. chevalieri, A. flavus, A. fumigatus, A. nidulans, A. niger, A. sydowii, and A. versicolor (167 to 169 nt) and among A. granulosus, A. restrictus, and A. ustus (170 to 172 nt). In some instances, similar to the results obtained in the ITS1 analysis, overall sequence length differences in the ITS2 region varied more within a given species than between species and ranged from 1 nt (e.g., A. niger) to 3 nt (A. flavipes) (Table 1).
Therefore, as noted for data generated using the D1-D2 and ITS1 regions, the overall sequence length of the ITS2 region alone could not be used reliably to differentiate among all 13 Aspergillus species examined. However, when data was combined for both the ITS1 and ITS2 regions, including or excluding the intervening, conserved, 5.8S rRNA gene, all 13 Aspergillus type strain sequences could be differentiated from one another (ITS1 plus ITS2 overall sequence length [Table 1]). The only overlaps in overall sequence length were between A. candidus (505 nt) and one non-type strain sequevar of A. flavipes (A. flavipes I; 505 nt) and between a non-type strain sequevar of A. nidulans (A. nidulans II, 480 nt) and A. sydowii.
Pairwise nucleotide sequence comparison of the D1-D2, ITS1, and ITS2 ribosomal regions among Aspergillus species investigated in this study. Table 2 shows the results of pairwise nucleotide sequence analysis of the D1-D2 region, and Tables 3 and 4 show the same analysis but for the ITS1 and ITS2 ribosomal regions. D1-D2 sequences were found to be highly conserved among the aspergilli (91.9 to 99.6% identity [Table 2]). The two species with the greatest similarity in nucleotide sequence were A. ustus and A. granulosus (99.6% identity); A. sydowii and A. versicolor sequences were also very similar and shared 99.1% sequence identity. Those species that were the most dissimilar in D1-D2 sequence were A. sydowii and A. restrictus (91.9% sequence identity). Phylogenetically related species, i.e., representatives of the Aspergillus subgenus Nidulantes (35), demonstrated more than 97% sequence identity: A. granulosus, 98.3 to 99.6% identity with A. ustus, A. sydowii, A. versicolor, and A. nidulans; A. ustus, 98.0 to 98.7% identity with A. sydowii, A. versicolor, and A. nidulans; and A. sydowii, 99.1% identity with A. versicolor. In addition, A. sydowii and A. versicolor shared 97.4 and 97.6% sequence identity, respectively, with A. nidulans.
|
View this table: [in a new window] |
TABLE 2. Pairwise sequence comparison in D1-D2 regions between medically important Aspergillus species investigated in this study
|
|
View this table: [in a new window] |
TABLE 3. Pairwise sequence comparison in ITS1 regions between medically important Aspergillus species investigated in this study
|
|
View this table: [in a new window] |
TABLE 4. Pairwise sequence comparison in ITS2 regions between medically important Aspergillus species investigated in this study
|
On the other hand, pairwise nucleotide sequence analysis of the ITS2 region demonstrated intermediate variation among Aspergillus species compared to the D1-D2 and ITS1 regions (Table 4). Sequence identities among the 13 Aspergillus species ranged from 75.6 to 98.3% identity. The two species with the greatest similarity in nucleotide sequence were A. granulosus and A. ustus (98.3% identity), agreeing with the D1-D2 similarity ranking; those that were the most dissimilar were A. terreus and A. restrictus (75.6% identity). Eight pairs of species, all belonging to the subgeneric Nidulantes group, demonstrated sequence identities between 93% and 98.3%: A. granulosus (93.6 to 98.3% identity with A. versicolor, A. sydowii, A. nidulans, and A. ustus), A. ustus (94.8% identity with A. nidulans), A. sydowii (95.9 and 97% identity with A. nidulans and A. versicolor, respectively), and A. nidulans (95.3% identity with A. versicolor).
Alignment of DNA sequences of the D1-D2, ITS1, and ITS2 ribosomal regions among Aspergillus species investigated in this study.
DNA sequence alignments of the D1 and D2 regions and of the ITS1 and ITS2 regions were conducted to identify areas within each region which displayed the greatest diversity and which might best discriminate among the Aspergillus species examined. Although some interspecies sequence divergence was observed in the more conserved D1 region, the most significant sequence divergence among Aspergillus species occurred in the more variable D2 region (i.e., in the D1 region, 27 [19.7%] divergent sites occurred over 137 aligned positions; in the D2 region, 55 [26.3%] divergent sites occurred over 209 aligned positions) (Fig. 1). Closely related taxa differed by
3 nt in the D2 region, e.g., A. granulosus versus A. ustus (1 nt) and A. sydowii versus A. versicolor (3 nt). Intraspecies sequence variability was
2 nt over the entire D1-D2 region (A. flavipes, A. sydowii, and A. nidulans differed by 1 to 2 nt from their respective type strains; remaining species differed by 0 nt).
![]() View larger version (46K): [in a new window] |
FIG. 1. Alignment of complete Aspergillus D1-D2 regions illustrating the sequence divergence among medically important species. The highly conserved region of the 28S rRNA gene intervening between the D1 and D2 regions has been omitted. Dots indicate identical nucleotides compared to the leader sequence (A. fumigatus strain ATCC 1022); dashes indicate alignment gaps. Sequence data were derived from this study (for GenBank accession numbers, see Table 1).
|
3 nt (A. flavipes: 3 nt, <2%; A. niger: 3 nt, <2%; A. nidulans: 2 nt, <2%) except for A. ustus strains that differed from the type strain sequence by 11 to 14 aligned positions in the ITS1 region.
![]() View larger version (62K): [in a new window] |
FIG. 2. Alignment of complete Aspergillus ITS1-ITS2 regions illustrating the sequence divergence among medically important species. The highly conserved intervening 5.8S rRNA gene has been omitted. Dots symbolize identical nucleotides compared to the leader sequence (A. fumigatus strain ATCC 1022); dashes indicate alignment gaps. Sequence data were derived from this study (for GenBank accession numbers, see Table 1).
|
5 nt (A. ustus, 5 nt, <3%; A. nidulans and A. niger, 1 nt, <1%; remaining species, 0 nt), except for A. flavipes strains that differed from the type strain sequence by 4 to 8 aligned nucleotides in the ITS2 region.
Comparative GenBank analysis of D1-D2, ITS1, and ITS2 ribosomal sequences.
To determine the usefulness of a public database to serve as a guide to species identification by comparative DNA sequence analysis, D1-D2, ITS1, and ITS2 sequences for each strain of the 13 Aspergillus species studied were used to conduct BLAST searches of the GenBank database (Table 5). All BLAST search results, using our own D1-D2 data as the query sequences, gave top-ranking scores (100% sequence identity) with corresponding reference sequences (35). GenBank searches revealed that sequences for strains of A. flavipes, A. sydowii, and A. versicolor were clearly distinct from those of other Aspergillus species or related taxa (i.e.,
4 nt or
99.3% similarity with any other GenBank entry for a different species); however, sequences for the remaining species (A. candidus, A. chevalieri, A. flavus, A. fumigatus, A. granulosus, A. nidulans, A. niger, A. restrictus, A. terreus, and A. ustus) were found to be identical with, or very similar to, reference sequences from closely related but different taxa (molecular siblings) (Table 5). For instance, A. flavus D1-D2 sequences were not distinguished from sequences representing Aspergillus oryzae, Aspergillus parasiticus, Aspergillus sojae, Aspergillus subolivaceus, or Aspergillus terricola (100% identity in the D1-D2 region, Table 5). Only 1 nt difference (99.8% identity) was observed between our D1-D2 sequences for strains of A. flavus and GenBank entries for Aspergillus tamari and Aspergillus flavofurcatus (Table 5). In total, 8 of the 13 Aspergillus species examined demonstrated identical D1-D2 sequences with at least one molecular sibling and the remainder of the species gave sequence identities of 99.3 to 99.8% with molecular siblings (Table 5).
|
View this table: [in a new window] |
TABLE 5. Molecular siblings of medically important Aspergillus species investigated in this study
|
Unlike GenBank search results using query sequences from the D1-D2 region, where 8 of the 13 species gave 100% identical sequence data for at least one molecular sibling and all remaining species gave sequence similarities of 99.3% to 99.8% with their respective molecular siblings (Table 5), ITS1 sequences gave 100% identity with molecular siblings for only 4 of the 13 Aspergillus species studied (A. chevalieri, A. flavus, A. nidulans, and A. niger); the remaining species gave sequence similarities ranging from 84.3% (A. granulosus) to 99.4% (A. restrictus) with sibling species (Table 5). ITS2 sequences gave 100% identity with molecular siblings for only 5 of the 13 Aspergillus species studied (A. chevalieri, A. flavus, A. nidulans, A. niger and A. ustus); the remaining species gave sequence similarities ranging from 91.5% (A. terreus) to 99.4% (A. niger) with sibling species (Table 5). However, it must be noted that ITS reference data were not available in GenBank for all sibling species, especially for Eurotium, Emericella, and Neosartorya species, at the time of our study.
Assessment of the quality of D1-D2, ITS1 and ITS2 sequence entries in GenBank. Inspection of BLAST alignments, and of multiple alignments generated with Aspergillus species D1-D2, ITS1, and ITS2 data from GenBank, revealed that many Aspergillus species sequences in this database had truncated ends and/or heterogeneities at positions found to be conserved at the generic or subgeneric level among the reference sequences of type strains and of authenticated culture collection strains (15, 35). These data indicate that many of the current GenBank sequences are not complete and/or contain errors in sequence data for the ribosomal regions studied. Further, BLAST search analysis of some query sequences gave best-match results with reference sequences of corresponding type strains but also resulted in excellent scores for sequences assigned to (i) unrelated or distantly related taxa, i.e., Cordyceps sinensis AJ488259, AJ488260, and AJ488268 (A. chevalieri query), Aureobasidium mansonii AF121288 (A. flavus query), and A. wentii U03522 and U03523, Arthrobotrys species U72602, Gliocladium cibotii AF021264, AF048739, and AF048738, and Verticillium bulbillosum AF048741 (A. niger query), or (ii) more closely related but different taxa, i.e., A. nidulans AF455499 (A. sydowii query), A. nidulans AF455424, AF455505, and AY452983 (A. versicolor query), and A. sydowii AJ312221 (A. nidulans query). When reference alignments were used for type strains, all of these GenBank entries were shown to exhibit significant nucleotide differences compared to respective reference sequences. These results indicate the existence of mislabeled ribosomal Aspergillus species data in GenBank ("false" molecular siblings). We could not determine, however, whether these ambiguous sequence annotations were due to misnamed or misidentified isolates or to manipulation or editing errors occurring during sequence analysis or GenBank submission.
|
|
|---|
Primary structural analysis of our Aspergillus species D1-D2 sequences revealed negligible intraspecies variability and recognizable interspecies divergence within the D2 region as described for other fungi (12, 24, 25, 36, 42). Although comparison of the D1-D2 sequences determined in this study with those in GenBank showed complete identity with corresponding reference data (35), most BLAST search results were inconclusive at the species level due to identical or very similar (>99% identical) reference data which were assigned to other organisms than the query sequences (molecular siblings). Nonetheless, these sibling organisms most frequently belonged to the same Aspergillus group as the query species (35) and supported the identification to at least the subgeneric group level. Consequently, results obtained by sequencing the D1-D2 region should be interpreted accordingly. The same test specificity applies to currently available large-subunit rRNA gene probes (42) which target D2 sequences conserved among medically important species and their molecular siblings. In addition, other targets with conserved interspecies sequences are expected to exhibit analogous diagnostic limitations, such as those reported previously using the mitochondrial cytochrome b gene (46, 58).
Unlike the D1-D2 region, the two ITS regions investigated were found to exhibit localized hypervariable regions containing most of the intraspecies sequence diversity. Nonetheless, most ITS amplicons showed a consistent sequence length among strains of the same species in one or both spacer regions; exceptions included sequevars of A. flavipes, A. nidulans, A. niger, and A. ustus. These hypervariable regions also contained most of the interspecies sequence diversity; however, similar sequence lengths in one or both spacer regions often occurred, making the overall size of ITS amplicons an imperfect tool to differentiate medically important Aspergillus species. This was especially evident when data from only one spacer region were analyzed. It is also not known if sibling species or other potential pathogens possess similar amplicon lengths that would confound a specific identification because, at present, GenBank data are incomplete in the ITS region. Our work therefore supports a previous report suggesting that A. fumigatus and A. flavus could not be differentiated by ITS amplicon length alone and that application of a specific DNA capture probe may be required for species identification (16). Alternatively, ITS amplicons have been characterized using single-strand conformation polymorphism (SSCP) analysis exploiting both size and sequence differences (23, 56); however, clinical diagnostic applications may be compromised by intraspecies SSCP pattern variability such as that described recently for strains of A. fumigatus and A. flavus (39). Multiple SSCP patterns are also predicted to occur within species of A. flavipes, A. nidulans, A. niger, A. ustus, and other Aspergillus species that exhibit sequence differences among strains of the same species (this study; 15, 34, 44, 53, 54).
The presence of intraspecies ITS variability did not hamper Aspergillus species identification by comparative sequence analysis, as indicated by similar BLAST search results for all sequevars of a particular species. All ITS sequences determined in this study yielded top-ranking BLAST scores with corresponding reference data (15, 34) regardless of the spacer region analyzed (ITS1 or ITS2) except for A. granulosus because its ITS sequence was not available in GenBank at the time of this study. Reference sequences of organisms most closely related to, but different from, the query species usually exhibited less than 99% sequence identity in at least one spacer region. Nonetheless, some ITS analyses were inconclusive at the species level because similar GenBank reference sequences existed for different organisms (molecular siblings of A. chevalieri, A. flavus, A. nidulans, and A. niger). Additional ambiguity arose from the existence of probably misclassified or misnamed GenBank sequences exhibiting excellent BLAST scores with some of our query sequences (A. chevalieri, A. flavus, A. nidulans, A. niger, A. sydowii, and A. versicolor) and yet assigned to different taxa (including Arthrobotrys, Aspergillus, Aureobasidium, Cordyceps, Gliocladium, and Verticillium species"false" molecular siblings [see Results]). Truncated and low-quality GenBank sequences were disregarded in the present study and therefore did not affect identification; however, truncated and low-quality GenBank entries clearly limit the utility of this database to assess inter- and intraspecies sequence similarity. Overall, ITS comparative sequence analysis, which included both spacer regions, showed better species differentiation than use of a single spacer region alone. However, differentiation between some more closely related Aspergillus species, especially among Eurotium and Emericella species, may require analysis of less conserved targets such as the ribosomal external transcribed spacer regions or intergenic spacer regions (18).
Considering the results presented above, the development of ribosomal approaches targeting Aspergillus groups rather than particular species would seem more practical. Ribosomal sequence analysis of the 13 medically important species investigated in this study allowed unambiguous identification of 11 distinct groups (A. flavipes and remaining groups as listed in the Results section) using either the D1-D2 or the ITS regions. Further, these molecular groups largely reflected relationships established previously by phenotypic classification (38, 45). Indeed, A. fumigatus group-specific primers (64), various sets of species-specific primers and group-specific probes (15), or panfungal primers in combination with group-specific probes (6, 42) may provide sufficient differentiation for diagnostic applications in a conventional clinical setting.
Nevertheless, ribosomal approaches for Aspergillus species group identification still requires important improvement in the GenBank database, including submission of additional sequence information for the ITS1 and ITS2 regions for sibling species, implementation of quality control measures for sequence entries using comparisons to known reference sequences, authoritative nomenclatural updates, and expert-based reclassification of mislabeled sequences (correct classification of currently "false" molecular siblings). Accurate sequence databases would not only benefit clinical diagnostic applications but also facilitate the ongoing effort to improve the current taxonomy of the genus Aspergillus, which still lacks complete definition by molecular criteria. Indeed, traditional classification and identification of Aspergillus species is based on morphological criteria usually expressed in a continuum rather than as discrete entities as discussed in detail by Raper and Fennell (38). To some extent, this may explain the intraspecies ribosomal divergence found in A. flavipes, A. niger, and A. ustus as well as within other species (this study; 15, 33, 34, 35, 44, 53, 54) and the potential diagnostic problems encountered with species-specific probes (30, 32). A true species definition for molecular taxonomic purposes is lacking at present for the aspergilli. Because DNA-DNA hybridization studies are excessively laborious and difficult to perform and interpret, and because many of the medically important Aspergillus species studied have no known sexual state (precluding the establishment of a biological species concept based on genetic crosses), other means will need to be employed to define species within the genus Aspergillus.
In conclusion, the medically important Aspergillus species investigated in this study can largely be identified at the traditional group level using ribosomal D1-D2 or ITS regions as diagnostic targets. Upon major improvement of present sequence databases, identification to the species level should be feasible through ITS sequence analysis; some clinical isolates, however, may require additional analyses. Such improved databases should enable a more accurate identification of Aspergillus species, thus promoting further advancements in the molecular diagnosis, epidemiology, and clinical management of invasive aspergillosis.
We thank Stephen W. Peterson and Arvind A. Padhye for kindly providing reference strains and helpful advice concerning the taxonomy of Aspergillus species. Brandi Reddick is acknowledged for excellent technical assistance.
Present address: Institut de Microbiologie, Centre Hospitalier Universitaire Vaudois, 1011 Lausanne, Switzerland. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»