| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Botany Department, University of Wisconsin-Madison, Madison, Wisconsin 53706
Andrea Gargas
Symbiology, LLC, Middleton, Wisconsin 53562
| ABSTRACT |
|---|
|
|
|---|
Oligonucleotide microarray based on ITS2 rDNA sequences would be extremely useful in identifying fungi within soil samples. However ITS2 contains phylogenetic information and duplication of sequences among taxa make false positive detections likely unless a way could be found to identify taxon-specific portions of the ITS2 sequence a priori. Examination of component ITS2 sequences suggested one method of identifying species-specific probes. Analysis of 168 fungal ITS2 sequences showed that all 168 ITS2 rRNA sequences could be folded to produce similar secondary structures of 3–4 loops. Unique probes occurred most often in the second loop. While the loop 2 sequence was unique in all taxa, there were partial congeneric and intergeneric duplicates. Evidence for a decrease in duplicates with increasing phylogenetic distance was mixed. From the evidence, 2 or 3 disjunct oligonucleotide probes from the loop 2 sequence might be sufficient to identify most fungal species. This combination appears minimally susceptible to false positives and conceivably could be extended to design probes to identify any eukaryotic species.
Key words: ITS2 secondary structure, microarray design, oligonucleotide probes, soil fungi
| INTRODUCTION |
|---|
|
|
|---|
An obvious approach is to create an oligonucleotide microarray that contains taxon-specific probes. Such an array would be simple to use; DNA could be extracted from a soil sample, hybridized to the array and the fungi present could be read from the array (although admittedly only known fungi would be found). Because the sequences for the ITS2 region of many fungi and oomycetes are readily available and the sequences are easy to extract from samples, this region appears to be a good target for finding unique probes. However ITS2 sequences do carry some phylogenetic information (Coleman 2003
, Schultz et al 2005
) and fungi that share identical regions of their ITS2 sequences will share probes. The probes ideally should contain only autapomorphies for species or strains of interest, not synapomorphies at higher phylogenetic levels such as genera, families and phyla.
Finding appropriate autapomorphies within the ITS2 region is even more difficult than it first appears. Given that most species of soil fungi are unknown (Hawksworth et al 1995
) the array designer has to create probes that not only will register the taxon of interest but will not respond to some unknown fungus, especially if the array is designed to sample wild soils. The only apparent solution is to find a region of ITS2 that contains many autapomorphies, probes based on that region would have a higher a priori probability of being taxon specific.
One possible solution is to use ITS2 loop structure to find regions containing many autapomorphies. The ITS2 region appears to fold in comparable ways across a number of species, including Saccharomyces cerevisiae (Joseph et al 1999
), plants (Mai and Coleman 1997
, Coleman 2003
), green algae (Coleman and Mai 1997
), Drosophila (Young and Coleman 2004
) and recently in a broad survey of 5000 sequences across the eukaryotes (Schultz et al 2005
). Numbers of autapomorphies appear to vary by position within the secondary structure (Coleman and Mai 1997
, Mai and Coleman 1997
, Coleman 2003
, Young and Coleman 2004
) and secondary structure has been used as a species-level character within the fungal genus Polyporus (Krüger and Gargas 2004
).
Following their lead we tested whether folded structures could be used to find regions with high numbers of autapomorphies likely to generate unique probes among fungal taxa. Our analysis of the ITS2 secondary structures, sequences and probes focused on three questions: (i) Do fungal ITS2 sequences have a common folding pattern? (ii) Does the number of duplicates per microarray probe depend on its position in its parent ITS2 sequence and structure? and (iii) Is there a correlation between number of duplicates and phylogenetic distance? In other words, if two fungi share a probe, is it more likely that they are near relatives? This last feature might be useful in designing an array because it would indicate that any fungus generating a false positive would be more likely to have a near relative on the array, rather than being a random organism that happened to match a particular probe.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Probes.— Numbers of probe duplicates within the dataset were calculated with Microsoft Access, both across the entire ITS2 sequence and by secondary structure. The number of duplicates was combined with probe numbering and mapping to determine the number of duplicates per nucleotide position along each ITS2 sequence. Mean number and standard deviation of duplicates per probe position were calculated for all sequences with Microsoft Excel. Because these ITS2 sequences were 130–334 nucleotides long their secondary structures also varied in length, and comparing duplicate numbers base pair by base pair was not practical. Thus we also compared the number of duplicates within each secondary structural feature with ANOVA with Type III sums of squares (to compensate for unequal sequence lengths within each structure) with differences compared by Tukeys HSD test. These tests were run with S-plus version 6.0 (2001, Insightful Corp., Seattle, Washington).
Probes and phylogenetic distance.— The ITS2 dataset could not be used to generate a phylogeny to test the correlation between the number of duplicates and phylogenetic distance because such a test would be circular. Phylogenetic distance therefore was determined in two ways.
The first approach was to use another gene to create the phylogeny and measure distances. Thirty-five representative species for the distance measurements and 35 proxy species (near relatives) were used to calculate phylogenetic distances (TABLE II
). A distance matrix was obtained from an alignment of these 35 species with heuristic search of PAUP* 4.0b 10 (Swofford 2002) using default settings with the distance optimality criterion set to minimum evolution. The number of duplicates was regressed against phylogenetic distance with two methods, a linear regression and a correlation whose probabilities were calculated from 999 permutations of the two datasets. This second method was implemented by Dr Bret Larget (University of Wisconsin at Madison).
|
| RESULTS |
|---|
|
|
|---|
|
|
|
|
Overall in no case was the entire L2 region shared between truly different genera. Although some parts of the L2 region might be shared, duplication events inevitably were clustered at either the beginning or end of the ITS2 sequence. In a few cases up to half of the sequence was duplicated among species, but in those cases sequences at the other end of the loop were unique.
| DISCUSSION |
|---|
|
|
|---|
Interspecific and intergeneric duplications did occur within L2, but this presents no major problems. A single probe cannot unambiguously identify a fungal species. Multiple, separated probes, for instance from the beginning, middle, and end of L2, likely will be sufficient to identify a species because we found no evidence of species that share entire loop 2 sequences. It is possible that groups of species share entire loop 2 sequences, but they did not appear in this sample. However any such groups identified in the future likely will appear in the literature and future researchers will be able to avoid problematic taxa. Searches of GenBank and other databases for related sequences should be a preliminary part of any probes designed with this technique.
The regression results do not unequivocally support the hypothesis that the number of duplicate probes decreases with increasing phylogenetic distance. The results of the duplicate counts in the entire dataset suggested otherwise, but that dataset could not be used to calculate phylogenetic distances in our research. More thorough sampling of sequences outside the ITS2 region likely will provide the accessions necessary for future re-examination of this hypothesis.
The fundamental conclusion is that ITS2 sequences can be used to design probes that are, if not species-specific, close to that level. The second loop of ITS2 has been shown to contain a high number of autapomorphies. By using a number of probes based on the loop 2 sequence, it should be possible to detect any known fungal species. Moreover a set of 2–3 loop 2 probes is resistant to false positives, a benefit given that unknown fungal species likely will be present in most environmental samples. While this analysis focused on probes for microarray design, these results are adaptable to any technology that uses oligonucleotide probes to sample communities. Finally, because eukaryotes in a number of phyla have been shown to have ITS2 sequences with similar secondary structures, ITS2 loop 2 sequences could be suitable for creating probes for any eukaryote.
|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Corresponding author. Tel: (310) 427-0624. E-mail: fclandis{at}wisc.edu
| LITERATURE CITED |
|---|
|
|
|---|
Bodrossy L, Sessitsch A. 2004. Oligonucleotide microarrays in microbial diagnostics. Curr Op Microbiol 7:245–254.[CrossRef][Medline]
Brasier CM. 1991. Ophiostoma novo-ulmi sp. nov., causative agent of the current wave of Dutch elm disease. Mycopathologia 115:151–161.[CrossRef]
Buisman CJ. 1932. Ceratostomella ulmi, de geslachtelijke vorm van Graphium ulmi Shwarz. Tijdshcr Plantenziek 38:1–5.
Coleman AW. 2003. ITS2 is a double-edged tool for eukaryote evolutionary comparisons. Trend Genet 19: 370–375.[CrossRef][Medline]
———, Mai JC. 1997. Ribosomal DNA ITS-1 and ITS-2 sequence comparisons as a tool for predicting genetic relatedness. J Mol Evol 45:168–177.[CrossRef][Medline]
Grime JP, Mackey JML, Hillier SH, Read DJ. 1987. Floristic diversity in a model system using experimental microcosms. Nature 328:420–422.[CrossRef]
Hartnett DC, Wilson GWT. 2002. The role of mycorrhizas in plant community structure and dynamics: lessons from grasslands. Plant Soil 244:319–331.[CrossRef]
Hawksworth DL, Kirk PM, Sutton BC, Pegler DN. 1995. Ainsworth & Bisbys Dictionary of the Fungi. Wallingford, UK: CAB International.
Ingrouille M. 1995. Historical Ecology of the British Flora. London: Chapman and Hall.
Joseph N, Krauskopf E, Vera MI, Michot B. 1999. Ribosomal internal transcribed spacer 2 (ITS2) exhibits a common core of secondary structure in vertebrates and yeast. Nucl Acid Res 27:4533–4540.
Krüger D, Gargas A. 2004. The basidiomycete genus Polyporus— an emendation based on phylogeny and putative secondary structure of ribosomal RNA molecules. Feddes Repert 115:530–546.[CrossRef]
Landis FC, ———, Givnish TJ. 2004. Relationships among arbuscular mycorrhizal fungi, vascular plants, and environmental conditions in oak savannas. New Phytologist 164:493–504.[CrossRef]
Mai JC, Coleman AW. 1997. The internal transcribed spacer 2 exhibits a common secondary structure in green algae and flowering plants. J Mol Evol 44:258–271.[CrossRef][Medline]
Mathews DH, Sabina J, Zuker M, Turner DH. 1999. Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J Mol Biol 288:911–940.[CrossRef][Medline]
Rizzo DM, Garbelotto M. 2003. Sudden oak death: endangering California and Oregon forest ecosystems. Front Ecol Environ 1:197–204.
———, ———, Davidson JM, Slaughter GW, Koike ST. 2002. Phytophthora ramorum as the cause of extensive mortality of Quercus spp. and Lithocarpus densiflorus in California. Plant Dis 86:205–214.[CrossRef]
Schleifer KH. 2004. Microbial diversity: facts, problems and prospects. Sys Appl Microbiol 27:3–9.[CrossRef]
Schultz J, Maisel S, Gerlach D, Muller T, Wolf M. 2005. A common core of secondary structure of the internal transcribed spacer 2 (ITS2) throughout the Eukaryota. RNA 11:361–364.
———, Klironomos JN, Ursic M, Moutoglis P, Streitwolf-Engel R, Boller T, Wiemken A, Sanders IR. 1998. Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396: 69–72.[CrossRef]
van der Heijden MGA. 2002. Arbuscular mycorrhizal fungi as a determinant of plant diversity: in search of underlying mechanisms and general principles. In: van der Heijden MGA, Sanders IR, eds. Mycorrhizal ecology. Ecological studies 157. Berlin: Springer-Verlag. p 243–265.
Weste G, Brown K, Kennedy J, Walshe T. 2002. Phytophthora cinnamomi infestation—a 24-year study of vegetation change in forests and woodlands of the Grampians, Western Victoria. Austral J Bot 50:247–274.[CrossRef]
Wills RT. 1993. The ecological impact of Phytophthora cinnamomi in the Stirling Range National Park, Western Australia. Austral J Ecol 18:145–159.[CrossRef]
Young I, Coleman AW. 2004. The advantages of the ITS2 region of the nuclear rDNA cistron for analysis of phylogenetic relationships of insects: a Drosophila example. Mol Phylogenet Evol 30:236–242.[CrossRef][Medline]
Zhou JH. 2003. Microarrays for bacterial detection and microbial community analysis. Curr Op Microbiol 6: 288–294.[CrossRef][Medline]
Zuker M. 2003. Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acid Res 31:3406–3415.
This article has been cited by other articles:
![]() |
C. Koetschan, F. Forster, A. Keller, T. Schleicher, B. Ruderisch, R. Schwarz, T. Muller, M. Wolf, and J. Schultz The ITS2 Database III--sequences and structures for phylogeny Nucleic Acids Res., January 1, 2010; 38(suppl_1): D275 - D279. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |