| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Centro de Investigação das Ferrugens do Cafeeiro, Instituto de Investigação Científica Tropical, Quinta do Marquês, 2784-505 Oeiras, Portugal, and Departamento de Biologia, Centro de Estudos da Macaronésia, Universidade da Madeira, 9000-390 Funchal, Portugal
Ana Ribeiro
Vítor M.P. Várzea
Carlos J. Rodrigues, Jr.
Centro de Investigação das Ferrugens do Cafeeiro, Instituto de Investigação Científica Tropical, Quinta do Marquês, 2784-505 Oeiras, Portugal
| ABSTRACT |
|---|
|
|
|---|
Random amplified polymorphic DNA (RAPD) was used to assess the genetic structure of Hemileia vastatrix populations. Forty-five rust isolates with different virulence spectra and from different hosts and geographical regions were analyzed. Out of 45 bands, generated with three RAPD primers, 35 (78%) were polymorphic and scored as molecular markers. Cluster analysis exhibits unstructured variability of this pathogen with regard to physiological race, geographical origin or host. The genotypic diversity (H') inferred from Shannons index was higher than gene diversity (Ht), suggesting that diversity is distributed among clonal lineages. Estimates of gene diversity in Africa and Asia populations were higher in total (Ht) as compared to within population diversity (Hs). Genetic differentiation was considerable among coffee rust isolates from Africa (Gst = 0.865) and Asia (Gst = 0.768) but not among isolates from South America (Gst = 0.266). We concluded that genetic diversity in H. vastatrix was moderately low and that the genetic differentiation among populations shows that asexual reproduction is likely to play an important role in the population biology of this fungus. This should be taken into account for the development of breeding programs.
Key words: coffee rust, DNA polymorphism, genetic variability, molecular markers, population structure
| INTRODUCTION |
|---|
|
|
|---|
H. vastatrix life cycle is understood incompletely. The fungus survives as urediniospores, uredia and mycelia, whereas the majority of rusts usually have five spore stages and two hosts. Despite the occasional production of teliospores and basidiospores, under cool dry conditions, no alternate host has been discovered. The basidiospores germinate in vitro but do not infect coffee leaves. Because basidiospores do not have a known function, it has been thought that physiological races arise as a result of mutation rather than genetic recombination (Rodrigues Jr. et al 2001
).
Disease control with fungicides is possible but inconvenient, marginally effective and environmentally hazardous (Thurston 1998
). Therefore efforts have been directed toward breeding for resistance in coffee plants. Despite breakthroughs achieved with the production of varieties derived from Híbrio de Timor (HDT) like Catimor, Sarchimor, Colombia and derivatives of these varieties, new virulent races of the pathogen continue to arise (Rodrigues Jr. et al 2001
, Várzea et al 2001
).
Evidence for variation in pathogenicity among H. vastatrix isolates was reported in the early 1930s (Rodrigues Jr. et al 1975
). Further evidence of pathogenic variability has been recognized by distinct phenotypes among Coffea sp. plants used as differentials (Bettencourt and Rodrigues Jr. 1988
, Várzea et al 2001
). The differentials comprise a collection of coffee germplasm on which different rust isolates are classified into physiological races, based on individual strains virulence. These cumbersome screening methods prompted us to carry out molecular studies because direct analysis of DNA polymorphisms is a more reliable approach to detect genetic variation in organisms, especially those with limited morphological characters such as H. vastatrix, and little is known about the diversity of this fungus.
Standard molecular methods such as random amplified polymorphic DNA (RAPD), amplified fragment length polymorphism (AFLP), restriction fragment length polymorphism (RFLP) and microsatellite analysis have been used to estimate genetic diversity and give insights into the population structure of important rusts such as Puccinia recondita (Kolmer et al 1995
), P. striiformis (Steele et al 2001
, Justesen et al 2002
), Melampsora epitea (Pei et al 1997
, Samils et al 2001
), Cronartium ribicola (Hamelin et al 1998
, Kinloch et al 1998
), C. flaccidum (Moricca and Ragazzi 1998
) and Peridermium pini (Hantula et al 1998
, Moricca and Ragazzi 1998
).
In the present study we used RAPD markers to examine molecular polymorphisms among isolates of H. vastatrix from different coffee growing regions and assess broader levels of genetic diversity. This is the first time that RAPDs have been applied to population genetic studies of this species. The purpose was to increase our knowledge about the pathogen variability and differentiation that might be used to investigate pathogen evolution and to design strategies for disease management as well as for developing new cultivars.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
RAPD analysis.
Standard arbitrary 10-mer oligonucleotides (Operon Technologies Inc., USA) were tested for RAPD analysis using the method of Williams et al (1990)
. Reactions were carried out in a final volume of 25 µL containing 20 mM Tris-HCl pH 8.4, 50 mM KCl, 3 mM MgCl2, 0.2 mM dNTPs, 0.3 µM primer, 50 ng of genomic DNA and 1 U of Taq DNA polymerase (Invitrogen). Amplifications were performed in a PTC-100 thermocycler (MJ Research Inc.) programmed for an initial denaturation cycle (2 min at 94 C) followed by 40 cycles of denaturation (1 min at 94 C), annealing (1 min at 35 C) and extension (2 min at 72 C) with a final extension of 5 min at 72 C. Each reaction was performed at least twice to assess the consistency of the band profiles. The reproducibility of the technique was calculated Si2= n/2t (Sneath and Johnson 1972
), where Si2 represents the variances of a primer set between replicate isolates, n the number of isolates with discrepancies and t the total number of duplicate isolates. Reaction products (10 µL) were resolved by electrophoresis in 1.2% agarose gels stained with ethidium bromide in 1x TAE buffer at 90 V for 140 min. All reactions included negative controls in which DNA or primer was omitted in the amplification reaction mixture. Because the urediniospores in this study were obtained from greenhouse-grown infected Coffea differentials, and in such conditions, contamination with the hyperparasitic fungus Verticillium may occur, a positive control with Verticillium sp. DNA also was included.
Data analysis.
Differences in fingerprinting patterns between isolates were assessed visually. Polymorphisms including faint bands that could be scored unequivocally were included in the analyses. Presumed homologous bands were scored as present (1) or absent (0) to create a binary matrix. The binary matrix was analyzed with the computer program NTSYSpc version 2.02h (Applied Biostatistics Inc., USA). An unweighted pair group arithmetic mean method (UPGMA) cluster analysis was performed based on the DICEs similarity coefficient SD= [2a/(2a + b + c)], where a is the number of bands common for isolate x and y, b is the number of bands present only in isolate x, and c is the number of bands present only in isolate y. Dendrograms were generated with the tree display option (TREE). Goodness of fit of the cluster analysis was determined by computing a cophenetic value matrix. A cophenetic correlation of r > 0.9 was considered a good fit. Clustering and statistical analyses were performed on data generated from each primer set separately and on the combined data primer sets. The ability of the most informative primers to differentiate between isolates was assessed by calculating their Resolving power (Rp) (Prevost and Wilkinson 1999
). Rp was calculated with the formula Rp =
Ib, where Ib (band informativeness) takes the value of 1[2x (0.5p)], and p is the proportion of isolates containing the band. The genetic variability was quantified in two ways: (i) the percentage of polymorphic loci (representing polymorphic bands) and (ii) Shannon information index (H'denotes the diversity of RAPD markers) (Lewontin 1972
). Calculations of these two diversity estimates were performed with Popgene version 1.31 (Yeh et al 1997
). Correlation coefficients were determined with the XLSTAT program version 6.1.8. The occurrence of any population structure was studied by calculating total gene diversity (Ht) of the polymorphic loci, mean gene diversity for each RAPD marker within populations (Hs) and genetic differentiation (Gst) across populations as a proportion of total diversity (Nei 1973
), using the Popgene program.
| RESULTS |
|---|
|
|
|---|
|
|
The patterns obtained with primer OPD-05 for isolates of the same race suggested that this primer has the ability to produce fingerprints specific to particular physiological races. For example most bands generated from isolates 37, 395 and 1821, identified as race III, were similar to other H. vastatrix isolates but yielded a distinct band at 0.85 kb (FIG. 1B
). Representatives of race XV, isolates 70 and 772 (FIG. 1B
), yielded a prominent band about 1.8 kb. In contrast isolates belonging to different physiological races, such as 2531 (race XLII) and 2542 (race XLI), lacked a 1.2 kb band (FIG. 1B
).
Multiple distinct DNA fragments were found among H. vastatrix isolates examined with primer OPA-09. Of the 19 bands analyzed 17 (89.5%) were polymorphic across all the isolates. Some DNA fragments were present exclusively in isolates from the same physiological race, which could be considered race specific. For example isolate 70 (FIG. 1C
), from race XV revealed a polymorphic band of 1.1 kb. Three bands (2.6 kb, 1.5 kb, 0.7 kb) similarly were specific for race II (isolate 1065 and isolate 1126, data not shown), while one band of about 1.2 kb was present in isolates 166, 178 and 2571 from race VIII (data not shown).
Cluster analysis.
Cluster analyses performed with the RAPD data matrix generated by each primer revealed that only primer OPA-09 grouped some isolates according to physiological races, although several races could not be distinguished (data not shown). The dendrogram obtained with all RAPD markers (FIG. 2
) showed that genetic relatedness of the 45 coffee rust isolates ranged from 0.7 to 1.0. The cophenetic correlation coefficient was 0.93, indicating the dendrogram was an excellent fit representation of the original data. Reproducibility of the RAPD technique averaged 92% and we used this value to set the level for isolate differentiation, so 12 (C1 to C12) main clusters were delineated. Judging from cluster analysis, isolate 2307 (cluster C12) is genetically unique, whereas extensive clonality was found in the larger collection from India (cluster C2). Isolates classified in different physiological races, such as 292 (XVII) and 292a (XXIII), were identical at the molecular level and were closely associated within the same cluster (C6, FIG. 2
). Isolates 2537(XL) and 2191 (XXXVII) similarly were indistinguishable at the DNA level (FIG. 2
, cluster C5). In several cases there was high genetic similarity between isolates from different geographical regions (e.g., cluster C1) and hosts (e.g., cluster C7). Few rust isolates were exceptions and could be grouped according to the indicated physiological race. However overall results of cluster analysis showed no correlation with physiological races (Spermans rank correlation coefficient: r = 0.127, P = 0.202), geographical origin (r = 0.159, P = 0.148) or host (r = 0.043, P = 0.390).
|
|
In the C. arabica population the Shannon diversity index was superior to the HDT population (TABLE III
). This variation was reflected by the higher polymorphism detected among isolates from C. arabica (64.4%) when compared with those derived from HDT (42.2%). A positive correlation was detected between the Shannon index and the percentage of polymorphism for the HDT population but not for the C. arabica population (Pearson rank correlation: r = 0.999, P = 0.011; r = 0.969, P = 0.080, respectively). We suggest that RAPD markers were unevenly distributed in the former population. In addition the coefficient of gene differentiation (Gst) was similar in rusts from C. arabica and HDT, indicating that these populations were highly different but isolates were relatively similar.
Genotypic diversity varied among geographical regions, based on Shannon indices (TABLE III
). The diversity within regions ranged from 0.140 in African to 0.241 in South American populations. A significant association (P < 0.05) was found between the Shannon index and the percentage of polymorphism for Africa and South America populations but not for Asian population (Pearson rank correlation: r = 0.801, P = 0.204). This absence of correlation probably was due to the larger number of physiological races and isolates in Asian population. Compared to gene diversity (Ht) the higher genotypic diversity (H') observed suggested that genetic diversity is distributed mainly among clonal lineages. Genetic subdivision analysis revealed the lowest level of genetic differentiation within the South American region (Gst = 0.266). This value may reflect that the isolates originated from two countries and mainly from HDT derivatives or that low genetic differentiation possibly was due to a more recent origin. In addition total gene diversity (Ht = 0.170) was attributable mostly to diversity within population (Hs = 0.125), indicating that isolates were likely to differ genetically.
| DISCUSSION |
|---|
|
|
|---|
Our results suggest that RAPD markers can be used to evaluate genetic diversity within the species H. vastatrix through estimates of variation at multiple loci (35 in this study) across the genome. The number of alleles per locus, a potentially more sensitive measure of genetic diversity, cannot be determined from our dominant marker data. Nevertheless dominant markers predictably can underestimate genetic diversity (Krutovskii et al 1999
, Wu et al 1999
), and if this is true for H. vastatrix then the true genetic diversity might be higher than found in this study.
The molecular variation yielded up to 78% polymorphic bands, which represents a high degree of polymorphism uncommon in rust fungi (Kinloch et al 1998
, Park et al 2000
, Li et al 2001
, Steele et al 2001
). The fingerprints obtained with primer OPF-06 identified polymorphic differences among isolates but generated only a limited pattern of one or two additional bands. Primers OPA-09 and OPD-05 conversely showed a higher resolution. Primer OPA-09 can be used to unequivocally discriminate races XV and VIII, whereas primer OPD-05 discriminates races III and XV. However, to be used as race-specific markers, further surveying is required of a larger number of isolates.
Comparison of the H. vastatrix isolates by means of RAPD (FIG. 2
) showed that rust isolates classified in the same physiological race by means of conventional methods (e.g., isolates from races XXIII and XXXVII) seemed to differ substantially at the molecular level. This variability could indicate the occurrence of genetic changes independent of pathogenicity. Cluster analysis suggests a widespread occurrence of homogeneous isolates of H. vastatrix, probably resulting from the introduction of isolates between countries as a consequence of germplasm exchange or long-distance dispersal (Brown and Hovmoller 2002
). Large dispersal with no apparent grouping, regardless of host or location, also was observed with isolates of Melampsora epitea (Samils et al 2001
), Botrytis elliptica (Huang et al 2001
) and Plasmodiophora brassicae (Manzanares-Dauleux et al 2001
).
The level of uniformity and the partitioning of genetic diversity suggest a clonal reproduction in the H. vastatrix population studied. The fungus has a population structure consistent with an asexual reproduction. This agrees with the absence of biological evidence for sexual behavior under natural and greenhouse conditions (Rodrigues Jr. et al 2001
). Considering the ability of long-distance dispersal (Brown and Hovmoller 2002
) clones can be spread worldwide, undergoing differentiation via mutation, asexual recombination or other sources of genetic variation (mechanisms occurring in asexually reproducing populations and favoring genetic exchange). Important phytopathogenic fungi, such as Cronartium ribicola (Hamelin et al 2000
), Melampsora spp. (Pei et al 2000
), Puccinia recondita (Park et al 2000
), P. striiformis (Hovmoller et al 2002
, Justesen et al 2002
), are known to have clonal lineages and little genetic diversity worldwide.
Substantial genetic differentiation was observed between the African and Asian populations, contrary to populations within the South American region. High genetic differentiation was observed in asexually propagated populations (Wyand and Brown 2003
). The differences observed in our study could be related to population age and level of gene flow among populations. H. vastatrix urediniospores are thought to have spread by air currents from Angola, where coffee leaf rust first was reported in 1966, to Brazil in 1970 (Bowden et al 1971
). In this study populations from South American region mainly were formed by isolates collected from HDT derivatives (78%). These isolates were introduced to Latin America about 30 years ago (Rodrigues Jr. et al 2001
). Since 1990 Colombian plants, derived from HDT, started showing susceptibility to rust (Rodrigues Jr. et al 2001
). The moderate degree of differentiation in South American populations is consistent with a relatively recent introduction of the pathogen and with recent observations concerning the susceptibility of cultivars derived from HDT to coffee rust (Rodrigues Jr. et al 2001
, Várzea et al 2001
). Moreover, considerable variation within South America populations suggests a high evolutionary rate of H. vastatrix. Such pathogen populations are more likely to overcome genetic resistance (McDonald and Linde 2002
) and might explain the frequent loss of resistance to the rust among Coffea genotypes.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Corresponding author. E-mail: mgouveia{at}uma.pt
| LITERATURE CITED |
|---|
|
|
|---|
Bowden J, Gregory PH, Johnson CG. 1971. Possible wind transport of coffee leaf rust across the Atlantic Ocean. Nature 224:500501.
Brown J, Whan J, Kenny M, Merriman P. 1995. The effect of coffee leaf rust on foliation and yield of coffee in Papua New Guinea. Crop Prot 14:589592.[CrossRef]
Brown JKM, Hovmoller MS. 2002. Aerial dispersal of pathogens on the global and continental scales and its impact on plant disease. Science
297:537541.
Chen X, Line RF, Leung H. 1995. Virulence and polymorphic DNA relationships of Puccinia striiformis f.sp. hordei to other rusts. Phytopathology 85:13351342.[CrossRef]
Hamelin RC, Hunt RS, Geils BW, Jensen GD, Jacobi V, Lecours N. 2000. Barrier to gene flow between eastern and western populations of Cronartium ribicola in North America. Phytopathology 90:10731078.[Medline]
, Dusabenyagasani M, Ettouil K. 1998. Fine-level genetic structure of white pine blister rust populations. Phytopathology 88:11871191.[Medline]
Hantula J, Niemi EM, Kaitera J, Jalkanen R, Kurkela T. 1998. Genetic variation of the resin top fungus in Finland as determined by random amplified microsatellites (RAMS). Eur J Forest Pathol 28:361372.[CrossRef]
Hovmoller MS, Justesen AF, Brown JKM. 2002. Clonality and long-distance migration of Puccinia striiformis f.sp. tritici in northwest Europe. Plant Pathol 51:2432.[CrossRef]
Huang J, Hsieh TF, Chastagner GA, Hsiang T. 2001. Clonal and sexual propagation in Botrytis elliptica. Mycol Res 105:833842.[CrossRef]
Justesen AF, Ridout CJ, Hovmoller MS. 2002. The recent history of Puccinia striiformis f. sp. tritici in Denmark as revealed by disease incidence and AFLP markers. Plant Pathol 51:1323.[CrossRef]
Kinloch BB, Westfall RD, White EE, Gitzedanner MA, Dupper GE, Foord BM, Hodgskiss PD. 1998. Genetics of Cronartium ribicola. IV. Population structure in western North America. Can J Bot 76:9198.[CrossRef]
Kolmer JA, Liu JQ. 2000. Virulence and molecular polymorphism in international collections of the wheat leaf rust fungus Puccinia triticina. Phytopathology 90:427436.[Medline]
, , Sies M. 1995. Virulence and molecular polymorphism in Puccinia recondita f. sp. tritici in Canada. Phytopathology 85:276285.[CrossRef]
Krutovskii KV, Erofeeva SY, Aagaard JE, Strauss SH. 1999. Simulation of effects of dominance on estimates of population genetic diversity and differentiation. J Hered 90:449502.
Lewontin RC. 1972. The apportionment of human diversity. Evol Biol 6:381394.
Li CX, Yeh FC, Hiratsuka Y. 2001. Random amplified polymorphic DNA variability among geographic isolates of western gall rust fungus in Canada. Can J Forest Res 31:13041311.[CrossRef]
Manzanares-Dauleux MJ, Divaret I, Baron F, Thomas G. 2001. Assessment of biological and molecular variability between and within field isolates of Plasmodiophora brassicae. Plant Pathol 50:165173.[CrossRef]
McDonald BA, Linde C. 2002. Pathogen population genetics, evolutionary potential, and durable resistance. Ann Rev Phytopathol 40:349379.[CrossRef][Medline]
Moricca S, Ragazzi A. 1998. Use of RFLP and SSCP analysis to differentiate the pine rusts Cronartium flaccidum and Peridermium pini. Mycol Res 102:666670.[CrossRef]
Nei M. 1973. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci
70:33213323.
Park RF, Jahoor A, Felsenstein FG. 2000. Population structure of Puccinia recondita in western Europe during 1995, as assessed by variability in pathogenicity and molecular markers. J Phytopathol 148:169179.[CrossRef]
Pei MH, Whelan MJ, Halford NG, Royle DJ. 1997. Distinction between stem- and leaf infecting forms of Melampsora rust on Salix viminalis using RAPD markers. Mycol Res 101:710.[CrossRef]
, Yuan ZW, Hunter T, Ruiz C. 2000. Heterogeneous nature of a new pathotype of Melampsora rust on Salix revealed by AFLP. Eur J Plant Pathol 106:771779.[CrossRef]
Prevost A, Wilkinson MJ. 1999. A new system of comparing PCR primers applied to ISSR fingerprinting of potato cultivars. Theor Appl Genet 98:107112.[CrossRef]
Raeder U, Broda P. 1985. Rapid preparation of DNA from filamentous fungi. Lett Appl Microbiol 1:1720.
Rodrigues Jr CJ, Bettencourt AJ, Rijo L. 1975. Races of the pathogen and resistance to coffee rust. Ann Rev Phytopathol 13:4970.[CrossRef]
, Várzea VMP, Silva MC, Guerra-Guimarães L, Rocheta M, Marques DV. 2001. Recent advances on coffee leaf rust. In: Prakash NS, Raghuramulu Y, Devasia J, eds. Proceedings of the International Scientific Symposium on Coffee. Bangalore Central Coffee Research Institute, Coffee Board. p 179193.
Samils B, Lagercrantz U, Lascoux M, Gullberg U. 2001. Genetic structure of Melampsora epitea populations in Swedish Salix viminalis plantations. Eur J Plant Pathol 107:399409.[CrossRef]
Sneath PHA, Johnson R. 1972. The influence of numerical taxonomic similarities of errors in microbiological tests. J Gene Microbiol 72:377392.
Steele KA, Humphreys E, Wellings CR, Dickinson MJ. 2001. Support for a stepwise mutation model for pathogen evolution in Australasian Puccinia striiformis f. sp. tritici by use of molecular markers. Plant Pathol 50:174180.[CrossRef]
Thurston HD. 1998. Coffee rust. In: Tropical Plant Diseases. 2nd ed. Minnesota: APS Press. p 123128.
Várzea VMP, Rodrigues Jr CJ, Silva MC, Gouveia M, Marques DV, Guerra-Guimarães L, Ribeiro A. 2002. Resistência do cafeeiro à H. vastatrix. In: Zambolin L, ed. O estado da arte de tecnologias na produção de café. Universidade Federal de Viçosa, Minas Gerais. p 297313.
, Gouveia M, Rodrigues Jr C, Ribeiro A, Marques D. 2001. Isolamentos de H. vastatrix virulentos em derivados do Híbrido de Timor. Fitopatologia Brasileira 26(suppl.):780 (abstract).
Williams JGK, Kubelik AR, Livak KJ, Rafalski JA, Tingey SV. 1990. DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucl Acids Res
18: 65316535.
Wu J, Krutovskii KV, Strauss SH. 1999. Nuclear DNA diversity, population differentiation, and phylogenetic relationship in the California closed-cone pines based on RAPD and allozyme markers. Genome 41:893908.
Wyand RA, Brown JKM. 2003. Genetic and forma specialis diversity in Blumeria graminis of cereals and its implication for host-pathogen co-evolution. Mol Plant Pathol 4:187198.[CrossRef]
Yeh FC, Yang RC, Boyle TBJ, Ye ZH, Mao JX. 1997. POPGENE, The User Friendly Software for Population Genetic Analysis. Molecular Biology and Biotechnology Center, University of Alberta, Canada.
This article has been cited by other articles:
![]() |
M. Gally, A. M. Ramos, D. Dokmetzian, and S. E. Lopez Genetic variability of Phytophthora sojae isolates from Argentina. Mycologia, November 1, 2007; 99(6): 877 - 883. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |