Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2004, p. 255-261, Vol. 70, No. 1
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.1.255-261.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Biotechnology Research Unit, Centro Internacional de Agricultura Tropical (CIAT) and Institut de la Recherche pour le Développement (IRD), Cali, Colombia
Received 14 April 2003/ Accepted 22 October 2003
|
|
|---|
|
|
|---|
In previous studies, X. axonopodis pv. manihotis populations have been characterized by restriction fragment length polymorphism (RFLP) (24, 27, 28) and virulence variation (20). The RFLP technique, using the pathogenicity gene pthB as a probe, is precise and discriminatory in differentiating among X. axonopodis pv. manihotis strains (28). It permitted the description of basic population genetic parameters and processes like geographical differentiation of the X. axonopodis pv. manihotis subpopulations in ECZs and the migration of strains between ECZs (22, 24). The role of host selection in X. axonopodis pv. manihotis population structure was also reported (18). The study of the virulence variation led to the conclusions that the cassava-X. axonopodis pv. manihotis interaction is pathotype-cultivar specific and that different pathotypes are present in Colombia and Venezuela (20, 22, 29). However, changes in the pathogen population structure over time have not been studied, and until now, little was known regarding the population dynamics of X. axonopodis pv. manihotis.
Understanding the mechanisms by which bacterial populations evolve to overcome host resistance is critical for developing strategies to improve the durability of resistance. This study examines the genetic change in X. axonopodis pv. manihotis populations over time from 1995 to 1999 by using RFLP and virulence assays. Temporal dynamics were studied within and between years. Since changes in pathogen pathotypes have long-term implications in disease management and varietal improvement, the pathotype composition among X. axonopodis pv. manihotis populations was also addressed.
|
|
|---|
![]() View larger version (10K): [in a new window] |
FIG. 1. Geographical distribution of the six cassava fields sampled in Colombia. In each ECZ, strains were collected at different sites. A, Cajibio; B, Mondomo; C, Carimagua; D, Villavicencio; E, Pivijay; F, Santander de Quilichao.
|
|
View this table: [in a new window] |
TABLE 1. Number of isolates, genetic diversities, disease incidence, number of haplotypes, and chi-square analysis of haplotypic frequencies for X. axonopodis pv. manihotis populations collected at six sites in four ECZs in Colombia
|
Characterization of isolates. (i) RFLP with pthB as probe.
Genomic DNA extractions, digestion, electrophoresis, and Southern blotting were performed as previously described (24). Membranes were hybridized with a plasmid probe, designated pthB, corresponding to a 5.6-kb EcoRI fragment. DNA sequencing of this fragment showed that the pthB gene is a member of the Xanthomonas avr-pth gene family harboring 12.5 internal 102-bp repeats and flanked by near-identical 980-bp direct repeats. The nucleotide sequence of the pthB gene is available in the GenBank database (accession number AF012325). Conditions of hybridization and washes were as previously described (24).
(ii) Virulence patterns.
Virulence patterns of 96 strains representing the 45 different RFLP haplotypes described in ECZs 1, 2, 2-4, and 5 from 1995 to 1999 were determined. Strains were inoculated as previously described on a set of differential cultivars previously selected for ECZs 1 and 2 (20). Briefly, plants were arranged in the greenhouse according to a randomized complete block design with five replicates. Noninoculated plants were included in each experiment as controls. Cassava stems were inoculated by inserting a toothpick contaminated by passage through a 24-h-old culture of the appropriate isolate. Symptoms were monitored 7, 14, and 30 days after inoculation, and disease was rated according to the following scale: 0, healthy plant, no reaction observed; 1, dark area or necrosis around the inoculation point; 2, gum exudates on the stem; 3, wilting of one or two leaves and exudates; 4, more than two leaves wilted; and 5, complete wilting and dieback. Cultivars with at least one rating of 4 or higher in a replicate were considered susceptible (20).
The cassava differential cultivars were the same as previously used for the characterization of pathotypes in ECZs 1 and 2 (20). For ECZ 5, we used three varieties previously studied (20) and included a larger set of cassava cultivars which had been selected for their adaptation to ECZ 5 (CM7436-7, CM7438-1, MBRA383, MCOL2061, SM1846-12, SM524-1, SM616-22, and CM7595-1). For ECZ 2-4, nine cultivars were selected based on their adaptation to the ECZ conditions (CM523-7, CM2177-2, CM6740-3, SM1219-9, CM7951-5, CM3306-4, CM7514-7, and CM849-1).
Statistical analysis.
Each distinct RFLP banding pattern, obtained after hybridization with probe pthB, was regarded as a haplotype. The haplotype diversity was estimated using the Nei and Tajima index (16), H = [n/(n - 1)](1 -
Xi2), where Xi is the proportion of the ith distinct pthB haplotype and n is the number of strains. Significance testing for differences in haplotype frequencies over years in each site or ECZ was performed by a chi-square analysis.
Three analyses of molecular variance (AMOVAs) (Arlequin version 1.1; S. Schneider, J. M. Kruffer, D. Roessli, and L. Excoffier, Genetics and Biometry Laboratory, University of Geneva, Geneva, Switzerland [http://lgb.unige.ch/arlequin/]) were performed. First, AMOVA was used to estimate variance components for RFLP-pthB haplotypes, partitioning the variation among ECZs and fields within ECZs and among fields. A second AMOVA was performed to study changes of X. axonopodis pv. manihotis population structure over a cropping cycle. The variation was partitioned into two populations, Mondomo and Pivijay, each including the two samplings performed at the beginning and at the end of the crop cycle. A third AMOVA was performed in order to assess if the genetic structure of X. axonopodis pv. manihotis had changed over the 5 years under study; strains collected within a year were considered a single population. The analysis was repeated for each ECZ. Genetic distance, calculated as Euclidean metric distance, was computed between all strains and between locations. The Euclidean metric distance between two individuals is equivalent to their total number of observed band differences for RFLP-pthB analysis.
|
|
|---|
The analysis of X. axonopodis pv. manihotis strains collected from 1995 to 1999 revealed extensive genetic variability among subpopulations of X. axonopodis pv. manihotis in different ECZs. Genetic diversities calculated for each site were highest for Carimagua and Villavicencio (ECZ 2) and lowest for Cajibio (ECZ 5) and Santander de Quilichao (ECZ 2-4) (Table 1). Except in Mondomo, genetic diversities decreased over time (Table 1). The AMOVA showed that 55.4% of the haplotype diversity was within locations, 34.8% was among ECZs, and 9.8% was among locations within ECZs. In order to determine if the population structure changed over time, an AMOVA was performed taking all isolates collected in 1 year in an ECZ as a population. More variation was ascribed to differences within the years than to differences among the years for all the ECZs, indicating that genetic diversity was maintained from year to year. The percentages of variation due to differences among and within years of collection were 13.6 and 86.4 for ECZ 2, 17.7 and 82.3 for ECZ 1, 41.6 and 58.4 for ECZ 2-4, and 42.2 and 57.8 for ECZ 5, respectively.
In some cases, disease incidence index was the average of data collected during various samplings in the same field. Disease severity varied from year to year in each ECZ. The most dramatic variation was observed in ECZ 1. In 1997 disease symptoms were rarely observed and disease severity was very low (disease severity of 1) compared to 1996 (disease severity of 3) and 1999 (disease severity of 4,5) (Table 1). The origin of this variation is not known and is difficult to estimate. However, this variation in the disease severity index in ECZ 1 can be correlated with the climatic conditions observed in ECZ 1 in 1997. Differences between minimum and maximum temperatures were low at the beginning of the year, and the amount of rainfall in July was significantly lower, almost nonexistent, compared to 1996 and 1998 (Fig. 2).
![]() View larger version (26K): [in a new window] |
FIG. 2. Rainfall and difference between minimum and maximum temperatures at Pivijay (ECZ 1) in 1996, 1997, and 1998.
|
|
View this table: [in a new window] |
TABLE 2. Haplotypes and genetic diversities (H) of the strains of X. axonopodis pv. manihotis populations collected in two samplings during a cycle of production
|
|
View this table: [in a new window] |
TABLE 3. X. axonopodis pv. manihotis pathotypes defined in ECZ 1
|
|
View this table: [in a new window] |
TABLE 6. X. axonopodis pv. manihotis pathotypes defined in ECZ 2-4
|
|
View this table: [in a new window] |
TABLE 4. X. axonopodis pv. manihotis pathotypes defined in ECZ 2
|
|
View this table: [in a new window] |
TABLE 5. X. axonopodis pv. manihotis pathotypes defined in ECZ 5
|
|
View this table: [in a new window] |
TABLE 7. Distribution of pathotypes of X. axonopodis pv. manihotis in Colombia during 1995 to 1999
|
|
|
|---|
In a previous study we identified migrations as a major factor influencing the pathogen population structure and maintaining the high genetic diversity (24). In this study RFLP variation in the X. axonopodis pv. manihotis populations indicated that a high level of variation at the DNA level has remained constant over time and we confirmed that migrations continue to shape the X. axonopodis pv. manihotis population structure. The presence of identical RFLP patterns in several fields of the same or different ECZs indicated pathogen migration. AMOVA confirmed the existence of migrations of strains between locations of the same ECZ. However, in comparing results obtained in 1995 and 1996 (24) with results from this study, it becomes evident that, over time, bacterial populations in sites of the same ECZ become more divergent and populations from different ECZs become less divergent. This change is a consequence of the increasing rate of migration of the pathogen among sites located in different ECZs. If this situation persists, cultivars to be deployed will have to be tested not only with the strains representing the pathotypes found in one ECZ but also with pathotypes from various ECZs.
Genetic changes can be subtle, resulting in simple numerical changes of haplotypes or pathotypes, or be dramatic and affect the pathogen population structure. Contingency chi-square tests showed that differences in haplotype frequencies over the 4 years were statistically significant (Table 1), indicating that X. axonopodis pv. manihotis populations at the site level were unstable. These changes in haplotype frequencies occur rapidly, in less than a year, and during a crop cycle a new population structure is established. In addition, AMOVA showed that genetic variation was higher within a year than among years, confirming how rapidly X. axonopodis pv. manihotis populations are able to evolve. We therefore recommend monitoring the pathogen population each year at the end of the crop cycle.
The influence of host selection or environmental factors on the X. axonopodis pv. manihotis population structure is difficult to estimate. However, as previously mentioned some correlations can be drawn between changes in the population structure and changes in the host population or in the climatic conditions. In fact, in a previous study (20) a correlation was established between host diversity and X. axonopodis pv. manihotis variation. The management and deployment of new cassava resistant genotypes had a significant impact on the population structure of X. axonopodis pv. manihotis over time. In Mondomo (ECZ 5), changes in haplotype frequencies were probably influenced by the change in the host population and the introduction of new cassava genotypes. During our first study in 1995 to 1996, only two varieties were cultivated in farmers' fields. The introduction of new genotypes was followed by the increase in haplotype diversity (Table 1) (20, 24). The same situation occurred in Villavicencio (ECZ 2) (Table 1), where populations changed after the introduction in 1997 of a pair-cross population of 150 individuals derived from an intraspecific cross between two noninbred cassava (M. esculenta Crantz) lines (11). In ECZs 2 and 5, weather conditions were very similar over time (data not shown), and changes in X. axonopodis pv. manihotis population structure may not be explained by changes in the environmental conditions alone.
Environmental factors may also have a major influence on disease development and population fluctuations, causing a genetic drift of the population. Temperature and rainfall are the most important factors that affect the severity of CBB (14, 25). The epiphytic phase of X. axonopodis pv. manihotis has important implications in the epidemiology of the disease, and the growth of epiphytical bacterial populations is favored by high humidity and large differences in day-night temperatures (7). In 1997 weather conditions were unfavorable for the development of the disease because of low rainfall as a consequence of the El Niño phenomenon. Consequently a very low percentage of plants showed symptoms at Pivijay (21) and the disease severity index was particularly low. We observed that some haplotypes disappeared while new haplotypes and pathotypes were detected. The predominant haplotype also changed. In Pivijay (ECZ 1) genetic drift had a significant impact on pathogen population structure.
The molecular analysis of diversity showed the potential for X. axonopodis pv. manihotis populations to change. If the haplotypic composition of the population changed during the crop cycles, the pathotypic composition could also change. We showed that the pathotypic composition really changed and the majority of new pathotypes appeared in ECZs 1 and 2 after 1997. As previously stated, a monitoring of the population at the end of the crop cycle and the resistance screening of cassava varieties with the new pathotypes will be the recommended strategy. Instability is the general rule in X. axonopodis pv. manihotis populations, as shown in this study, and genetics had and probably will continue to have a significant impact on the X. axonopodis pv. manihotis population structure. After 1997 in ECZ 1, some haplotype C6 strains that were dominant before 1997 persisted, exhibiting rare pathotypes. Because of this situation we believe that more than one resistance genotype should be deployed at the same time to prevent these rare pathotypes that are potential sources of resistance breakdown from becoming dominant.
In ECZs 5 and 2-4, the high level of virulence showed by the ECZ 5 and 2-4 strains suggests a need for characterizing more suitable sources of resistance. For ECZ 2 more differential cultivars are needed, first because X. axonopodis pv. manihotis populations in ECZ 2 are highly variable and second because, with five differential cultivars, only a limited number of pathotypes can be described. It was the case with ECZ 5, where in a previous study only one pathotype was described (20). By increasing the number of differential cultivars, we characterize new pathotypes that have been present since 1995. In general, results did not show a clear correlation between haplotypes and pathotypes. It could be interesting to know whether the molecular diversity in the avr gene structure of X. axonopodis pv. manihotis may account for the phenotypic diversity observed. The RFLP-EcoRI analysis performed in this study may not detect changes in the copies of the avr genes that are causing a variation in the virulence patterns of the strains. Changes in the avr structure that cause a variation in the virulence pattern of the strain can be due to single mutations at the 3' terminus of the avrXa7 alleles (4, 26). One strategy could be to sequence the copies of the plasmid-borne avr gene copies in the new pathotypes. It could also be interesting to determine if some of these avr copies in X. axonopodis pv. manihotis have a measurable fitness function.
This study is a part of our continuing effort to understand the population structure of the X. axonopodis pv. manihotis. Knowledge of pathogen population structure and dynamics is useful in deciding the way in which genotypes are deployed in time and space (13). As shown above, a molecular monitoring of the haplotypic composition of the population reflects a change of the pathotypes in each ECZ. If new pathotypes appear in the populations, these must be characterized and included in the set of strains used to screen cultivars to be deployed.
Monitoring the pathogen population requires extensive and repeated surveys. The development and recent application of geographical information systems, geostatistics, and the Kriging technique to plant pathology offer an alternative approach to predict and monitor changes in plant pathogen populations (12, 17). For the future it is important to develop a spatial simulation model to understand the influence of various ecological and genetic factors on the distribution and potential change of X. axonopodis pv. manihotis population structure. Recently we developed a database containing information on pathogen diversity and CBB incidence, host diversity and distribution, and climatic conditions for all the ECZs. This is the first step towards a better understanding of the biology and management of CBB.
This research was supported by grants from IRD and CIAT and by a doctoral fellowship awarded to S. Restrepo by the Institut de la Recherche pour le Développement (IRD).
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»