-
PDF
- Split View
-
Views
-
Cite
Cite
Bart Lievens, Petra M. Houterman, Martijn Rep, Effector gene screening allows unambiguous identification of Fusarium oxysporum f. sp. lycopersici races and discrimination from other formae speciales, FEMS Microbiology Letters, Volume 300, Issue 2, November 2009, Pages 201–215, https://doi.org/10.1111/j.1574-6968.2009.01783.x
Close - Share Icon Share
Abstract
During infection of tomato, the fungus Fusarium oxysporum f. sp. lycopersici secretes several unique proteins, called ‘secreted in xylem’ (Six) proteins, into the xylem sap. At least some of these proteins promote virulence towards tomato and among them, all predicted avirulence proteins that can trigger disease resistance in tomato have been found. In this study, a large, worldwide collection of F. oxysporum isolates was screened for the presence of seven SIX genes (SIX1–SIX7). The results convincingly show that identification of F. oxysporum formae speciales and races based on host-specific virulence genes can be very robust. SIX1, SIX2, SIX3 and SIX5 can be used for unambiguous identification of the forma specialis lycopersici. In addition, SIX4 can be used for the identification of race 1 strains, while polymorphisms in SIX3 can be exploited to differentiate race 2 from race 3 strains. For SIX6 and SIX7, close homologs were found in a few other formae speciales, suggesting that these genes may play a more general role in pathogenicity. Host specificity may be determined by the unique SIX genes, possibly in combination with the absence of genes that trigger resistance in the host.
Introduction
Fusarium oxysporum Schlechtend:Fr is an asexual fungus that is common in soils worldwide. Collectively, F. oxysporum strains can cause wilt or root, bulb or foot rot in a wide variety of plant species, among which are several economically important crops (Gordon & Martyn, 1997). Individual strains of F. oxysporum, however, usually infect only one or a few host species. Pathogenic strains have therefore been grouped into host-specific forms called formae speciales (f. spp.), which are sometimes divided further into races based on cultivar specificity (Armstrong & Armstrong, 1981; Di Pietro et al., 2003; Michielse & Rep, 2009). Fusarium oxysporum strains have also been assigned to vegetative compatibility groups (VCGs) (Puhalla, 1985), which correspond to clonal lineages of the fungus (Correll, 1991; Koenig et al., 1997; Kistler et al., 1998; Katan & Katan, 1999). While a particular forma specialis may cause disease in a certain plant species, strains belonging to other formae speciales may have a harmless or even a beneficial relation to the same species and vice versa (Recorbet et al., 2003). Therefore, discrimination between strains pathogenic and nonpathogenic towards a specific crop is essential in order to prevent unnecessary disease control efforts.
Classically, identification of pathogenic F. oxysporum isolates is based on pathogenicity testing (Recorbet et al., 2003), which is time consuming and laborious. In addition, as presently over 70 formae speciales have been described, an enormous number of plant species and cultivars should be used for correct strain identification (Fravel et al., 2003). Therefore, attempts are increasingly being made to replace these methods with molecular identification techniques (Lievens et al., 2008). Unfortunately, molecular discrimination of F. oxysporum isolates is seriously complicated by the polyphyletic nature of many formae speciales, such that isolates belonging to different formae speciales may be more related than isolates belonging to the same forma specialis (Kistler et al., 1997; Lievens et al., 2008). Ideally, molecular identification of F. oxysporum strains is based on DNA sequences directly related to (host-specific) pathogenicity or nonpathogenicity (Recorbet et al., 2003; Lievens et al., 2008).
In many cases, the ability of a fungus to infect particular plant species depends on specific genes encoding host-determining ‘virulence factors’ that distinguish virulent from avirulent strains. These include small secreted proteins, called effectors, and enzymes involved in the synthesis of host-specific toxins (van der Does & Rep, 2007). Recently, several in planta secreted proteins have been identified in F. oxysporum f. sp. lycopersici (Sacc.) W.C. Snyder and H.N. Hans. Fusarium oxysporum f. sp. lycopersici is the causal agent of Fusarium wilt in tomato and has been reported in at least 32 countries (Jones et al., 1991). Several polymorphic resistance genes have been identified in tomato that each confers resistance against a subset of F. oxysporum f. sp. lycopersici strains. These resistance genes include I (for immunity), I-1, I-2 and I-3 (Huang & Lindhout, 1997). Races are named historically according to the resistance gene that is effective against them: the I gene and the (unlinked) I-1 gene are effective only against race 1; I-2 confers resistance to race 2 (which overcomes I and I-1); and I-3 confers resistance to race 3 (which overcomes I, I-1 and I-2) (Rep et al., 2005). Race 1 was initially described in 1886 (Booth, 1971). Race 2 was first reported in 1945 in Ohio (Alexander & Tucker, 1945) and race 3 was originally observed in Australia in 1978 (Grattidge & O'Brien, 1982). Subsequently, the different races have been reported in tomato crops worldwide. The first in planta secreted protein that was identified in F. oxysporum f. sp. lycopersici, called ‘secreted in xylem 1’ (Six1), is a small cysteine-rich protein required for full virulence on tomato (Rep et al., 2005). In addition, recognition of the protein by tomato plants carrying the resistance gene I-3 leads to disease resistance (Rep et al., 2004). Therefore, Six1 is also called Avr3 to indicate its gene-for-gene relationship with the I-3 resistance gene. More recently, additional fungal proteins were identified from xylem sap of infected plants, encompassing the small secreted proteins Six2, Six3, Six4, Six5, Six6 and Six7, an arabinanase, an oxidoreductase and a serine protease (Houterman et al., 2007; van der Does, unpublished data). While the function of most of these small proteins is unknown so far, for two of them, an avirulence function has been established using gene knockout experiments. These are Six4 (Avr1), which is required for I and I-1-mediated resistance (Houterman et al., 2008), and Six3 (Avr2), which is required for I-2-mediated resistance (Houterman et al., 2009). In addition, Six4/Avr1 was found to suppress I-2- and I-3-mediated disease resistance (Houterman et al., 2008). Race 2 strains are thought to have arisen through loss of AVR1 (SIX4) from race 1 strains, while race 3 strains appear to have evolved from race 2 through point mutations in AVR2 (SIX3) (Houterman et al., 2009). The strong link between these SIX genes and pathogenicity towards tomato makes them excellent markers for host- and cultivar-specific pathogenicity.
In this study, we show the usefulness of effector genes for reliable identification of host-specific fungal pathogens. Particularly, we show that F. oxysporum f. sp. lycopersici and its races can be unambiguously identified based on the above-mentioned SIX genes. The robustness of this approach is underscored by the fact that, like many other formae speciales, F. oxysporum f. sp. lycopersici strains do not have a common ancestor within the F. oxysporum species complex (O'Donnell et al., 1998; van der Does et al., 2008).
Materials and methods
Fungal isolates
A worldwide collection of 270 F. oxysporum strains, obtained from diverse geographic origins, was assembled in this study (Table 1). One hundred and sixty-four strains isolated from tomato were used, encompassing 15 avirulent isolates, 75 F. oxysporum f. sp. lycopersici strains and 74 isolates belonging to F. oxysporum f. sp. radicis-lycopersici, representing most known VCGs of these formae speciales (Katan & Katan, 1999). While both formae speciales share tomato as the same host, they cause different symptoms: F. oxysporum f. sp. lycopersici causes wilt and F. oxysporum f. sp. radicis-lycopersici causes root and foot rot (Menzies et al., 1990). In addition, 106 F. oxysporum isolates of 14 other formae speciales were included in our study as well as seven isolates of three other Fusarium species (Table 1). For most of these isolates, pathogenicity, vegetative compatibility and genetic diversity have been assessed in previous studies (e.g. Katan et al., 1991; Marlatt et al., 1996; O'Donnell et al., 1998; Katan & Katan, 1999; Vakalounakis & Fragkiadakis, 1999; Baayen et al., 2000; Cai et al., 2003; Vakalounakis et al., 2004; Balmas et al., 2005; Kawabe et al., 2005; Lievens et al., 2007; van der Does et al., 2008). Isolates were grown on potato dextrose agar containing 0.1 mg mL−1 streptomycin sulfate in the dark at 22 °C.
Fusarium oxysporum strains used in this study
ATCC, American Type Culture Collection, Manassas, VA; CABI, Centre for Agriculture and Bioscience International, Surrey, UK; CBS, Centraalbureau voor Schimmelcultures, Utrecht, the Netherlands; DSM, Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH, Braunschweig, Germany; MAFF, Genetic Resources Management Section, GenBank, National Institute of Agrobiological Sciences, Tsukuba, Ibaraki, Japan; MUCL, Mycotèque de l'Université Catholique de Louvain, Louvain-la-Neuve, Belgium; NRRL, Agriculture Research Service Culture Collection, National Center for Agricultural Utilization Research, Peoria, IL.
Kawabe et al. (2005) showed that the VCG 0030 tester isolate NRRL 26037 and some other isolates of VCG 0030 were weakly compatible with VCG 0032 tester isolate FOL-MM66. In addition, FOL-MM59 (VCG 0032) was found to be a bridging isolate that can also form heterokaryons with some members of VCG 0030 (Cai et al., 2003; Kawabe et al., 2005). sc, self-compatible; si, self-incompatible.
If known, race is mentioned for F. oxysporum f. sp. lycopersici. A pathogenicity test on a tomato line carrying the gene revealed that isolate CBS 645.78 was originally misidentified as race 1, while isolates CBS 646.78 and MX295 should be designated as race 1 isolates instead of race 2 and race 3, respectively. In addition, nine isolates of unknown race were tested and turned out to be race 1. These strains included CBS 164.85, CBS 165.85, CBS 758.68, DSM 62338, FOL 00/60309/1, FOL-295A, FRC-0-1113A, FRC-0-1113N and NRRL 26034.
A.C.M.C., A.C.M. Clerckx, Plant Research International, Wageningen, the Netherlands; B.L., B. Lievens, Scientia Terrae Research Institute, Sint-Katelijne-Waver, Belgium; C.H.B., C.H. Beckman, University of Rhode Island, Kingston, RI; D.A.L., D.A. Lawn, Asgrow Seed Company, Sinaloa, Mexico; D.J.M., D.J. McGrath, The University of Queensland, St. Lucia, Qld, Australia; D.J.V., D.J. Vakalounakis, National Agricultural Research Foundation (N.AG.RE.F.), Plant Protection Institute, Heraklion, Crete, Greece; D.M.E., D.M. Elgersma, Willie Commelin Scholten Phytopathological Laboratory, the Netherlands; F.V., F. Vey, Laboratoire National de la Protection des Végétaux, Fleury-les-Aubrais Cédex, France; H.C.K., H.C. Kistler, University of Minnesota, St. Paul, MN; IPO-DL, Research Institute for Plant Protection (currently Plant Research International), Wageningen, the Netherlands; J.P.J., J.P. Jones, University of Florida, Gainesville, FL; J.W.S., J.W. Scott, University of Florida, Gainesville, FL; K.S.E., K.S. Elias, Louisiana State University Agricultural Center, Baton Rouge, LA; L.D., L. Davidse; M.I.G.R., M.I.G. Roncero, Universidad de Cordoba, Cordoba, Spain; N.A.M.V.S., N.A.M. van Stekelenburg, the Netherlands; P.S., P. Suleman, University of Maryland, College Park, MD; Q.M., Q. Mighelli, University of Sassari, Sassari, Italy; R.J.S., R.J. Scheffer, S&G Seeds BV, Enkhuizen, the Netherlands; T.K., T. Katan, Hebrew University of Jerusalem, Jerusalem, Israel; Y.T.W., Y.-T. Wang, Soochow University, Taipei, Taiwan; Z.K.P., Z.K. Punja, Simon Fraser University, Burnaby, BC, Canada.
+, amplicon detected; −, no amplicon detected. A fragment of 992, 749, 608, 967, 667, 793 and 862 bp is generated for SIX1, SIX2, SIX3, SIX4, SIX5, SIX6 and SIX7, respectively. As a check for DNA quality, all DNA samples were successfully subjected to PCR analysis using the universal primers ITS5 and ITS4, targeting the rRNA gene (White et al., 1990).
Sequence determined.
Previously misidentified as F. oxysporum f. sp. radicis-lycopersici, recently classified as F. oxysporum f. sp. lycopersici (Kawabe et al., 2005; van der Does et al., 2008).
Previously misidentified as F. oxysporum f. sp. lycopersici, recently determined as nonpathogenic to tomato (van der Does et al., 2008).
Fusarium oxysporum strains used in this study
ATCC, American Type Culture Collection, Manassas, VA; CABI, Centre for Agriculture and Bioscience International, Surrey, UK; CBS, Centraalbureau voor Schimmelcultures, Utrecht, the Netherlands; DSM, Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH, Braunschweig, Germany; MAFF, Genetic Resources Management Section, GenBank, National Institute of Agrobiological Sciences, Tsukuba, Ibaraki, Japan; MUCL, Mycotèque de l'Université Catholique de Louvain, Louvain-la-Neuve, Belgium; NRRL, Agriculture Research Service Culture Collection, National Center for Agricultural Utilization Research, Peoria, IL.
Kawabe et al. (2005) showed that the VCG 0030 tester isolate NRRL 26037 and some other isolates of VCG 0030 were weakly compatible with VCG 0032 tester isolate FOL-MM66. In addition, FOL-MM59 (VCG 0032) was found to be a bridging isolate that can also form heterokaryons with some members of VCG 0030 (Cai et al., 2003; Kawabe et al., 2005). sc, self-compatible; si, self-incompatible.
If known, race is mentioned for F. oxysporum f. sp. lycopersici. A pathogenicity test on a tomato line carrying the gene revealed that isolate CBS 645.78 was originally misidentified as race 1, while isolates CBS 646.78 and MX295 should be designated as race 1 isolates instead of race 2 and race 3, respectively. In addition, nine isolates of unknown race were tested and turned out to be race 1. These strains included CBS 164.85, CBS 165.85, CBS 758.68, DSM 62338, FOL 00/60309/1, FOL-295A, FRC-0-1113A, FRC-0-1113N and NRRL 26034.
A.C.M.C., A.C.M. Clerckx, Plant Research International, Wageningen, the Netherlands; B.L., B. Lievens, Scientia Terrae Research Institute, Sint-Katelijne-Waver, Belgium; C.H.B., C.H. Beckman, University of Rhode Island, Kingston, RI; D.A.L., D.A. Lawn, Asgrow Seed Company, Sinaloa, Mexico; D.J.M., D.J. McGrath, The University of Queensland, St. Lucia, Qld, Australia; D.J.V., D.J. Vakalounakis, National Agricultural Research Foundation (N.AG.RE.F.), Plant Protection Institute, Heraklion, Crete, Greece; D.M.E., D.M. Elgersma, Willie Commelin Scholten Phytopathological Laboratory, the Netherlands; F.V., F. Vey, Laboratoire National de la Protection des Végétaux, Fleury-les-Aubrais Cédex, France; H.C.K., H.C. Kistler, University of Minnesota, St. Paul, MN; IPO-DL, Research Institute for Plant Protection (currently Plant Research International), Wageningen, the Netherlands; J.P.J., J.P. Jones, University of Florida, Gainesville, FL; J.W.S., J.W. Scott, University of Florida, Gainesville, FL; K.S.E., K.S. Elias, Louisiana State University Agricultural Center, Baton Rouge, LA; L.D., L. Davidse; M.I.G.R., M.I.G. Roncero, Universidad de Cordoba, Cordoba, Spain; N.A.M.V.S., N.A.M. van Stekelenburg, the Netherlands; P.S., P. Suleman, University of Maryland, College Park, MD; Q.M., Q. Mighelli, University of Sassari, Sassari, Italy; R.J.S., R.J. Scheffer, S&G Seeds BV, Enkhuizen, the Netherlands; T.K., T. Katan, Hebrew University of Jerusalem, Jerusalem, Israel; Y.T.W., Y.-T. Wang, Soochow University, Taipei, Taiwan; Z.K.P., Z.K. Punja, Simon Fraser University, Burnaby, BC, Canada.
+, amplicon detected; −, no amplicon detected. A fragment of 992, 749, 608, 967, 667, 793 and 862 bp is generated for SIX1, SIX2, SIX3, SIX4, SIX5, SIX6 and SIX7, respectively. As a check for DNA quality, all DNA samples were successfully subjected to PCR analysis using the universal primers ITS5 and ITS4, targeting the rRNA gene (White et al., 1990).
Sequence determined.
Previously misidentified as F. oxysporum f. sp. radicis-lycopersici, recently classified as F. oxysporum f. sp. lycopersici (Kawabe et al., 2005; van der Does et al., 2008).
Previously misidentified as F. oxysporum f. sp. lycopersici, recently determined as nonpathogenic to tomato (van der Does et al., 2008).
PCR analysis and sequencing
All isolates listed in Table 1 were subjected to PCR analysis using primers that were designed previously, amplifying SIX1, SIX2 and SIX3 (Rep et al., 2004; van der Does et al., 2008;Table 2). In addition, all isolates were screened using primers that were located just outside the respective ORFs of SIX4, SIX5, SIX6 and SIX7 (Table 2). In order to check DNA quality, PCR amplification was performed using the universal primers ITS5 and ITS4, which anneal to conserved regions of the 18S and 28S rRNA genes, respectively (White et al., 1990). Genomic DNA was extracted using the phenol–chloroform extraction method as described previously (Lievens et al., 2003) and the yield was determined spectrophotometrically. PCR amplification was carried out in a reaction volume of 20 μL, containing 0.15 mM of each dNTP, 0.5 μM of each primer, 1 × Titanium Taq DNA polymerase, 1 × Titanium Taq PCR buffer (Clontech Laboratories, Palo Alto, CA) and 5 ng genomic DNA. Thermal cycling conditions consisted of 2 min at 94 °C, followed by 35 cycles of 45 s at 94 °C, 45 s at 59 °C and 45 s at 72 °C, with a final elongation step at 72 °C for 10 min. Amplified products were resolved electrophoretically in a 1.5% agarose gel. All reactions were performed twice. Following PCR amplification, some amplicons were cloned into the pCR2.1-TOPO vector using the TOPO TA cloning kit (Invitrogen, San Diego, CA) according to the manufacturer's instructions. Subsequently, both amplicon strands were sequenced using the M13 forward and reverse primer. DNA sequences were aligned using the clustalw algorithm to assess the potential intraspecific sequence variety. Sequences for each gene variant obtained is this study were deposited in GenBank under the accession numbers GQ268948–GQ268960.
Primers used in this study
| Code | Sequence (5′–3′) | Target | Amplicon length (bp) | Reference |
| P12-F2B | GTATCCCTCCGGATTTTGAGC | SIX1 | 992 | van der Does et al. (2008) |
| P12-R1 | AATAGAGCCTGCAAAGCATG | SIX1 | Rep et al. (2004) | |
| SIX2-F2 | CAACGCCGTTTGAATAAGCA | SIX2 | 749 | van der Does et al. (2008) |
| SIX2-R2 | TCTATCCGCTTTCTTCTCTC | SIX2 | van der Does et al. (2008) | |
| SIX3-F1 | CCAGCCAGAAGGCCAGTTT | SIX3 | 608 | van der Does et al. (2008) |
| SIX3-R2 | GGCAATTAACCACTCTGCC | SIX3 | van der Does et al. (2008) | |
| SIX4-F1 | TCAGGCTTCACTTAGCATAC | SIX4 | 967 | – |
| SIX4-R1 | GCCGACCGAAAAACCCTAA | SIX4 | – | |
| SIX5-F1 | ACACGCTCTACTACTCTTCA | SIX5 | 667 | – |
| SIX5-R1 | GAAAACCTCAACGCGGCAAA | SIX5 | – | |
| SIX6-F1 | CTCTCCTGAACCATCAACTT | SIX6 | 793 | – |
| SIX6-R1 | CAAGACCAGGTGTAGGCATT | SIX6 | – | |
| SIX7-F1 | CATCTTTTCGCCGACTTGGT | SIX7 | 862 | – |
| SIX7-R1 | CTTAGCACCCTTGAGTAACT | SIX7 | – |
| Code | Sequence (5′–3′) | Target | Amplicon length (bp) | Reference |
| P12-F2B | GTATCCCTCCGGATTTTGAGC | SIX1 | 992 | van der Does et al. (2008) |
| P12-R1 | AATAGAGCCTGCAAAGCATG | SIX1 | Rep et al. (2004) | |
| SIX2-F2 | CAACGCCGTTTGAATAAGCA | SIX2 | 749 | van der Does et al. (2008) |
| SIX2-R2 | TCTATCCGCTTTCTTCTCTC | SIX2 | van der Does et al. (2008) | |
| SIX3-F1 | CCAGCCAGAAGGCCAGTTT | SIX3 | 608 | van der Does et al. (2008) |
| SIX3-R2 | GGCAATTAACCACTCTGCC | SIX3 | van der Does et al. (2008) | |
| SIX4-F1 | TCAGGCTTCACTTAGCATAC | SIX4 | 967 | – |
| SIX4-R1 | GCCGACCGAAAAACCCTAA | SIX4 | – | |
| SIX5-F1 | ACACGCTCTACTACTCTTCA | SIX5 | 667 | – |
| SIX5-R1 | GAAAACCTCAACGCGGCAAA | SIX5 | – | |
| SIX6-F1 | CTCTCCTGAACCATCAACTT | SIX6 | 793 | – |
| SIX6-R1 | CAAGACCAGGTGTAGGCATT | SIX6 | – | |
| SIX7-F1 | CATCTTTTCGCCGACTTGGT | SIX7 | 862 | – |
| SIX7-R1 | CTTAGCACCCTTGAGTAACT | SIX7 | – |
Primers used in this study
| Code | Sequence (5′–3′) | Target | Amplicon length (bp) | Reference |
| P12-F2B | GTATCCCTCCGGATTTTGAGC | SIX1 | 992 | van der Does et al. (2008) |
| P12-R1 | AATAGAGCCTGCAAAGCATG | SIX1 | Rep et al. (2004) | |
| SIX2-F2 | CAACGCCGTTTGAATAAGCA | SIX2 | 749 | van der Does et al. (2008) |
| SIX2-R2 | TCTATCCGCTTTCTTCTCTC | SIX2 | van der Does et al. (2008) | |
| SIX3-F1 | CCAGCCAGAAGGCCAGTTT | SIX3 | 608 | van der Does et al. (2008) |
| SIX3-R2 | GGCAATTAACCACTCTGCC | SIX3 | van der Does et al. (2008) | |
| SIX4-F1 | TCAGGCTTCACTTAGCATAC | SIX4 | 967 | – |
| SIX4-R1 | GCCGACCGAAAAACCCTAA | SIX4 | – | |
| SIX5-F1 | ACACGCTCTACTACTCTTCA | SIX5 | 667 | – |
| SIX5-R1 | GAAAACCTCAACGCGGCAAA | SIX5 | – | |
| SIX6-F1 | CTCTCCTGAACCATCAACTT | SIX6 | 793 | – |
| SIX6-R1 | CAAGACCAGGTGTAGGCATT | SIX6 | – | |
| SIX7-F1 | CATCTTTTCGCCGACTTGGT | SIX7 | 862 | – |
| SIX7-R1 | CTTAGCACCCTTGAGTAACT | SIX7 | – |
| Code | Sequence (5′–3′) | Target | Amplicon length (bp) | Reference |
| P12-F2B | GTATCCCTCCGGATTTTGAGC | SIX1 | 992 | van der Does et al. (2008) |
| P12-R1 | AATAGAGCCTGCAAAGCATG | SIX1 | Rep et al. (2004) | |
| SIX2-F2 | CAACGCCGTTTGAATAAGCA | SIX2 | 749 | van der Does et al. (2008) |
| SIX2-R2 | TCTATCCGCTTTCTTCTCTC | SIX2 | van der Does et al. (2008) | |
| SIX3-F1 | CCAGCCAGAAGGCCAGTTT | SIX3 | 608 | van der Does et al. (2008) |
| SIX3-R2 | GGCAATTAACCACTCTGCC | SIX3 | van der Does et al. (2008) | |
| SIX4-F1 | TCAGGCTTCACTTAGCATAC | SIX4 | 967 | – |
| SIX4-R1 | GCCGACCGAAAAACCCTAA | SIX4 | – | |
| SIX5-F1 | ACACGCTCTACTACTCTTCA | SIX5 | 667 | – |
| SIX5-R1 | GAAAACCTCAACGCGGCAAA | SIX5 | – | |
| SIX6-F1 | CTCTCCTGAACCATCAACTT | SIX6 | 793 | – |
| SIX6-R1 | CAAGACCAGGTGTAGGCATT | SIX6 | – | |
| SIX7-F1 | CATCTTTTCGCCGACTTGGT | SIX7 | 862 | – |
| SIX7-R1 | CTTAGCACCCTTGAGTAACT | SIX7 | – |
Race determination/verification
Race 1 verification of a selection of F. oxysporum f. sp. lycopersici isolates was performed using the tomato lines GCR161, which contains the I gene and is resistant to race 1 isolates, and C32, which has general susceptibility to F. oxysporum f. sp. lycopersici (Kroon & Elgersma, 1993). Pathogenicity was tested using the root dip method (Wellman, 1939). Briefly, spores were collected from 5-day-old cultures in Czapek Dox broth (Difco) and used for root inoculation of 10-day-old plants at a spore density of 107 mL−1. The seedlings were then potted individually and grown at 25 °C in a greenhouse. Three weeks after inoculation, disease was scored in two ways, including average plant weight above the cotyledons and phenotype scoring according to a disease index ranging from 0 (no symptoms) to 4 (heavily diseased or dead) (Rep et al., 2004).
Results and discussion
During infection, F. oxysporum f. sp. lycopersici secretes several unique, small proteins into the xylem sap that promote virulence of the fungus towards tomato (Rep et al., 2004; Houterman et al., 2007). Among them, all three predicted avirulence genes were found, including AVR1=SIX4, AVR2=SIX3 and AVR3=SIX1 (Rep et al., 2004; Houterman et al., 2008, 2009). The linkage between this group of genes and pathogenicity makes them potential host- and cultivar-specific pathogenicity markers. This hypothesis was tested by assessing the presence of the different SIX genes in a large, worldwide collection of fungal isolates (Table 1). All isolates listed in Table 1 were tested for the presence of the different SIX genes by PCR assays, using primers that anneal just outside the ORF. Most SIX genes were found to be present in the forma specialis lycopersici, but not in other formae speciales or nonpathogenic isolates (Table 1). While SIX1–SIX5 are exclusively present in F. oxysporum f. sp. lycopersici, SIX6 and SIX7 amplicons were also generated from isolates of a few other formae speciales, including lilii (SIX7; 862 bp), melonis (SIX6; 793 bp) and radicis-cucumerinum (SIX6; 793 bp) (Table 1). Sequencing of these homologs revealed that they are identical within each forma specialis, but different between formae speciales, except for the SIX6 homologs of the formae speciales melonis and radicis-cucumerinum, which are identical. This high conservation is in line with the high degree of conservation that was also observed for these and other SIX genes between isolates of F. oxysporum f. sp. lycopersici (Houterman et al., 2008, 2009; van der Does et al., 2008; this study). In comparison with F. oxysporum f. sp. lycopersici, 91% DNA sequence identity (ORF) and 84% amino acid identity was found for the SIX7 homolog in F. oxysporum f. sp. lilii. For the SIX6 homologs, 95% DNA sequence identity (ORF) and 90% amino acid identity was found between the formae speciales melonis/radicis-cucumerinum and lycopersici. Additional screening revealed that isolates of F. oxysporum f. sp. vasinfectum, pathogenic on cotton, also carry a homolog of SIX6 (J. Ellis, pers. commun.). In accordance with Houterman et al. (2008), SIX4 (AVR1) was generally found to be present in race 1 isolates only (Table 1), supporting its gene-for-gene relationship with the I resistance gene (Houterman et al., 2009). Remarkably, SIX4 was absent in isolate CBS 645.78, which was previously classified as F. oxysporum f. sp. lycopersici race 1 (Table 1). On the other hand, a positive PCR was obtained for isolate CBS 646.78 and isolate MX395, which were previously determined as race 2 and race 3, respectively (Table 1). To check the racial identity of these isolates, a pathogenicty test on a tomato line carrying the I gene was performed, revealing that isolate CBS 645.78 was able to overcome the I gene and was therefore originally misidentified as race 1 isolate, while the last two isolates should be designated as race 1 isolates because they are avirulent on the I tomato line. In order to further assess the link between the unique presence of SIX4 and race 1 designation (i.e. avirulence on an I tomato line), nine additional isolates of unknown race that were shown to have SIX4 by PCR (i.e. CBS 164.85, CBS 165.85, CBS 758.68, DSM 62338, FOL 00/60309/1, FOL-295A, FRC-0-1113A, FRC-0-113N and NRRL 26034) were tested for their ability to cause disease on tomato seedlings containing the I resistance gene. All isolates turned out to be race 1, further demonstrating a perfect correlation between race 1 and the presence of SIX4 as well as the ability to detect misidentifications with our PCR assay. The results of two of these pathogenicity assays are presented in Fig. 1. In order to differentiate race 2 from race 3 isolates, SIX3 (AVR2) can be exploited. Indeed, race 2 strains were found to contain identical SIX3 sequences (also identical to those of race 1 strains), which differ from race 3 strains by single point mutations (Houterman et al., 2009; Fig. 2). Three variants of this gene were found within race 3 isolates differing in a single nucleotide from race 1 and race 2 isolates (within the ORF: G121>A, G134>A and G137>C; Fig. 2). Based on these nucleotide differences, three forward primers, SIX3-G121A-F2 (5′-ACGGGGTAACCCATATTGCA-3′), SIX3-G134A-F2 (5′-TTGCGTGTTTCCCGGCCA-3′) and SIX3-G137C-F1 (5′-GCGTGTTTCCCGGCCGCCC-3′), were developed and combined with the reverse primer SIX3-R2 (Table 2), enabling unambiguous PCR differentiation of race 2 and race 3 isolates when using stringent PCR conditions (i.e. using the PCR program described above, with the exception of an annealing temperature of 67 °C, an elongation time of 2 s and 30 instead of 35 cycles) (Table 3).
Race determination of Fusarium oxysporum f. sp. lycopersici isolates by pathogenicity tests on tomato. To verify or determine the race of individual F. oxysporum f. sp. lycopersici isolates, tomato lines with single resistance genes were inoculated at the seedling stage and disease was scored 3 weeks after inoculation. Disease severity was measured by average plant weight above the cotyledons as compared with mock inoculation (mock). In the two experiments shown here (a and b), the plant line GCR161 was used to determine race 1 identity – this line contains only the I gene and is therefore resistant only to race 1. As controls, previously verified F. oxysporum f. sp. lycopersici isolates were included: IPO1530/B1 (race 1; alternative designation: Fol004), D2 (race 2; alternative designation: Fol007) and 5397 (race 3; alternative designation: Fol029). To assess the pathogenicity of each isolate, they were also inoculated on the generally susceptible tomato line C32 (i.e. without resistance genes against F. oxysporum f. sp. lycopersici). Apart from the race 2 and race 3 control strains, all the other strains tested in these two experiments were race 1, except for CBS 645.78, which is race 2 (avirulent on an I-2 line, results not shown). Error bars indicate SD (n=20).
Sequence alignment of part of the different versions of the SIX3 gene found in Fusarium oxysporum f. sp. lycopersici. The LSU-7 sequence represents race 1 and race 2 sequences, whereas isolates FOL-MM10, 14844 and IPO3 represent race 3 strains. Variable nucleotides are underlined, marked in bold and with asterisks.
Differentiation of Fusarium oxysporum f. sp. lycopersici race 3 from race 1 or race 2 isolates by PCR
+, amplicon detected; −, no amplicon detected. A positive signal for one of the three PCR reactions is expected for race 3 isolates. A fragment of 429, 414 or 412 bp is generated when using SIX3-G121A-F2/SIX3-R2, SIX3-G134A-F2/SIX3-R2 or SIX3-G137C-F1/SIX3-R2, respectively.
Differentiation of Fusarium oxysporum f. sp. lycopersici race 3 from race 1 or race 2 isolates by PCR
+, amplicon detected; −, no amplicon detected. A positive signal for one of the three PCR reactions is expected for race 3 isolates. A fragment of 429, 414 or 412 bp is generated when using SIX3-G121A-F2/SIX3-R2, SIX3-G134A-F2/SIX3-R2 or SIX3-G137C-F1/SIX3-R2, respectively.
Rapid detection and reliable identification of potential plant pathogens is required for taking appropriate and timely disease management measures. As DNA sequence variation in the commonly used housekeeping genes such as those for rRNA, β-tubulin or translation elongation factor-1α is not sufficient for the unambiguous identification of most formae speciales (O'Donnell et al., 1998; Lievens et al., 2007), other approaches are increasingly being used. These include, for example, methods based on the use of molecular markers identified by genotyping or amplification of transposon insertions (Lievens et al., 2008). Recently, another PCR method, based on polymorphisms in two genes encoding cell wall-degrading enzymes (PG1 and PGX4), has been developed that distinguishes Japanese F. oxysporum f. sp. lycopersici strains (and its races) from Japanese F. oxysporum f. sp. radicis-lycopersici strains (Hirano & Arie, 2006). However, cross-reaction with other formae speciales was observed (Hirano & Arie, 2006), limiting the use of this method in practice. In contrast, our PCR screen showed a 100% success rate in discriminating F. oxysporum f. sp. lycopersici from other formae speciales (targeting SIX1, SIX2, SIX3 and/or SIX5) and identification of the different races (targeting SIX3 and SIX4). In addition, an experiment in which the presence of the pathogen was assessed in plant tissue showed that our PCR assays have the potential to detect and identify F. oxysporum f. sp. lycopersici in environmental samples, even before the plants have developed disease symptoms (data not shown). Consequently, our results pose potential benefits for tomato breeders and growers. These include, for example, accurate pathogen identification ensuring that effective control measures can be adopted or, conversely, unnecessary efforts can be avoided to control populations of F. oxysporum that are harmless (Lievens et al., 2008). In addition, our results can be used in breeding programs to evaluate or monitor disease resistance against specific pathogenic forms or races, for example by the development of quantitative real-time PCR assays enabling detection and quantification of pathogen biomass in planta (Brouwer et al., 2003). Apart from these practical implications, our results also shed new light on the possible origin of formae speciales of F. oxysporum. All SIX genes, except SIX4, are located on a single, relatively small chromosome (chromosome 14 in the sequenced F. oxysporum f. sp. lycopersici isolate 4287; http://www.broad.mit.edu/annotation/genome/fusarium_group/) (van der Does et al., unpublished data), and most of them are exclusively present in F. oxysporum f. sp. lycopersici (Table 1). Nevertheless, a few SIX genes (SIX6 and SIX7) have close homologs in a few other formae speciales (Table 1). The ORFs of these homologs were found to be intact, indicating that these genes may have a function in pathogenicity and may be part of (a) different pathogenicity chromosome(s), with a gene content partly overlapping with chromosome 14. The SIX genes shared between different formae speciales may play a more general role in pathogenicity, while host specificity may be determined by a combination of unique genes. There are at least 122 transposons and 227 predicted genes, including the SIX genes, on chromosome 14 (M. Rep, unpublished data). The high transposon content of ‘pathogenicity’ chromosomes may explain why in several studies specific transposon-based markers could be linked to pathogenicity (Fernandez et al., 1998; Chiocchetti et al., 1999; Pasquali et al., 2004, 2007). Evolution of an ancestral pathogenicity chromosome within the F. oxysporum species complex and transfer between different VCGs (van der Does et al., 2008) may have given rise to a variety of host-specific pathogenicity chromosomes. Genome sequencing of isolates belonging to different formae speciales should reveal whether or not this is a likely scenario.
Acknowledgements
The authors thank the ‘Vlaams Instituut voor de bevordering van het Wetenschappelijk-Technologisch Onderzoek-Vlaanderen’ (IWT-040169) and De Ceuster Corp. (Sint-Katelijne-Waver, Belgium) for financial support. In addition, we are grateful to the donors of the Fusarium strains.
References
Editor: Bernard Paul

