2020 Volume 70 Issue 4 Pages 462-473
Bacterial wilt, caused by the Ralstonia pseudosolanacearum species complex, is an important vascular disease that limits tomato production in tropical and subtropical regions. Two major quantitative trait loci (QTL) of bacterial wilt resistance on chromosome 6 (Bwr-6) and 12 (Bwr-12) were previously identified in Solanum lycopersicum ‘Hawaii 7996’; however, marker-assisted breeding for bacterial wilt resistance is not well established. To dissect the QTL, six cleaved amplified polymorphic sites (CAPS) and derived CAPS (dCAPS) markers within the Bwr-6 region and one dCAPS marker near Bwr-12 were developed, and resistance levels in 117 tomato cultivars were evaluated. Two markers, RsR6-5 on chromosome 6 and RsR12-1 on chromosome 12, were selected based on the genotypic and phenotypic analysis. The combination of RsR6-5 and RsR12-1 effectively distinguishes resistant and susceptible cultivars. Furthermore, the efficiency of the two markers was validated in the F3 generation derived from the F2 population between E6203 (susceptible) and Hawaii 7998 (resistant). Resistant alleles at both loci led to the resistance to bacterial wilt. These markers will facilitate marker-assisted breeding of tomato resistant to bacterial wilt.
Bacterial wilt, caused by Ralstonia solanacearum, is one of the most destructive diseases that affects many plant species. The R. solanacearum species complex (RSSC) has been classified into races, biovars, and phylotypes based on host range, carbon source usage, 16S-23S rRNA gene sequence of the strains, respectively (Cho et al. 2018, Hayward 1991, Jeong et al. 2007). Genomic analysis and proteomic profiling of various strains of the pathogen collected from different countries classified RSSC into three species: R. solanacearum (Phylotype II), R. pseudosolanacearum (Phylotypes I and III), and R. syzygii (Phylotype IV). Further genomic analysis of the species classified R.syzygii into three subspecies named syzygii, indonesiensis, and celebesensis (Prior et al. 2016).
The pathogen has a wide host range. Tomato and other Solanaceae plants are major hosts. The disease threatens the cultivation of these crops in tropical and subtropical regions and heated greenhouses in temperate regions because high temperatures are better suited for the pathogen and disease development. As a result, the pathogen causes substantial economic losses (Hayward 1991, Lopes and Rossato 2018). The pathogen moves into plant roots via natural openings, such as hydathodes, or damaged areas and proliferates in the xylem tissues. It then damages the xylem tissues and blocks the water flow, leading to the total collapse and death of susceptible plants (Bae et al. 2015, Lowe-Power et al. 2018). Xylem colonization and spread are necessary for bacterial wilt disease progress because mutations in xylem colonization rendered pathogen strains incapable of causing wilting in tomato plants (Schell 2000). Evaluation of core collections of the three fruit vegetables of Solanaceae crops (tomato, eggplant, and pepper) against different strains of the pathogen showed that resistance of tomato collections is low compared with eggplant and pepper (Lebeau et al. 2011).
Different control strategies, such as chemical, biological, and cultural practices, can reduce bacterial wilt severity, but none of them are effective (Yuliar et al. 2015). Development and use of resistant cultivars is the most effective approach to control bacterial wilt (Abebe et al. 2016, Huet 2014, Scott et al. 2005, Wang et al. 2018). Breeding for bacterial wilt-resistant cultivars has been challenging due to the polygenic nature of resistance, broad host range, variability of the pathogen strains, and effect of environmental factors that directly influence the phenotypic expression of the disease (Danesh et al. 1994, Fegan and Prior 2005, Hayward 1991, Lee et al. 2011, Thoquet et al. 1996b, Tran and Kim 2010). Moreover, the resistance locus is linked to undesirable horticultural traits (Scott et al. 2005). Solanum lycopersicum ‘Hawaii 7996’ (hereafter Hawaii 7996) is a stable resistant resource against bacterial wilt across various geographic locations and different bacterial strains with the highest average survival rate of 97% (Wang et al. 1998). Analysis of bacterial wilt resistance using F2, F3, and recombinant inbred line (RIL) populations derived from a cross between Hawaii 7996 (resistant) and West Virginia 700 (susceptible) identified QTL on chromosomes 4, 6, and 11 (Thoquet et al. 1996a); 3, 4, 6, 8, and 10 (Thoquet et al. 1996b); 6 (Mangin et al. 1999); 6, 8, and 12 (Wang et al. 2000); 3, 4, and 6 (Carmeille et al. 2006); and 3, 6, and 12 (Wang et al. 2013). Among them, a QTL on chromosome 6 (Bwr-6) was stable when measured with different phenotyping criteria (area under disease progress curve and bacterial colonization), in different bacterial strains (race 3-phylotype II and race 1-phylotype I), and as measured in two different seasons (hot and cold) (Carmeille et al. 2006). Bwr-12 was an active QTL specifically against Pss4 (race 1, biovar 4) (Wang et al. 2000). QTL on chromosome 6 (Bwr-6) and 12 (Bwr-12) are thought to be responsible for the stable resistance of Hawaii 7996. Bwr-12 covers 2.8 cM between the SSR markers SLM12-12 and SLM12-2, controlling more than 50% of the phenotypic variation in some trials. Bwr-6 was localized between SLM6-124 and SLM6-110 covering 15.5 cM and controlling up to 22.2% of the phenotypic variation (Wang et al. 2013). Although two QTL, Bwr-6 and Bwr-12, were repeatedly confirmed as major contributors to bacterial wilt resistance in Hawaii 7996, the genetic nature of these critical QTL remained unidentified. We have previously conducted whole-genome resequencing of two susceptible cultivars (Heinz 1706 and BWS-3) and seven resistant cultivars (Hawaii 7996, Hawaii 7998, 10-BA-3-33, 10-BA-4-24, BWR-1, BWR-22, and BWR-23) of tomato and identified genome-wide single nucleotide polymorphisms (SNPs) in resistant and susceptible groups of cultivars (Kim et al. 2018). The highest number of polymorphic SNPs in coding regions were found on chromosome 12 (168 SNPs) followed by chromosome 6 (53 SNPs). These SNPs might be associated with resistance to bacterial wilt. Based on the SNP information generated by re-sequencing, an HRM marker (KHU-1) that is tightly linked to Bwr-12 was developed; however, no tightly linked SNP-based molecular marker was developed to trace Bwr-6 resistance due to the large interval (~12.7 Mbp) (Kim et al. 2018).
In this study, we further analyzed the Bwr-6 region to dissect and develop a diagnostic molecular marker for this important QTL. Cleaved amplified polymorphic site (CAPS) and derived CAPS (dCAPS) markers were developed within Bwr-6 and screened for their diagnostic potential using 117 tomato genotypes. Phenotypic and genotypic analysis using a wide range of germplasms enabled us to select RsR6-5 as a diagnostic marker for Bwr-6 among the newly developed markers in the region. Consequently, this marker, in combination with RsR12-1, effectively distinguished bacterial wilt-resistant and wilt-susceptible tomato cultivars. The newly developed marker RsR6-5 together with RsR12-1 will promote marker-assisted breeding of tomato by targeting two major resistance QTL against bacterial wilt.
In total, 117 tomato cultivars were collected (either seed or genomic DNA sample) from different sources, including the Tomato Genetics Resource Center (TGRC), UC Davis; Kyung Hee University, Korea; National Agrobiodiversity Center (RDA-Genebank), Korea; and various commercial seed companies in Korea. The phenotype of 27 cultivars was confirmed by inoculation test in this study, 12 cultivars were inferred from previous reports and the phenotype of 78 cultivars was received from the respective company/supplier along with the genomic DNA sample (Table 1). Seeds were first sown in Petri dishes for germination, and germinated seeds were transferred to 128 cell seedling trays “(28 × 28 × 40, bottom 15 mm)” filled with bio mix (JM bio, Korea). The seedlings were grown in the Agricultural Experiment Station of Kyungpook National University in a glasshouse at an average temperature of 25–28°C and 16–8 h light-dark cycles. The seedlings were moved from the glasshouse to the growth chamber 3–4 days before inoculation for acclimatization to growth chamber conditions where they were kept post-inoculation. Four-week-old seedlings were used for inoculation.
Phenotyping | Number of cultivars | Bacterial wilt phenotype | |
---|---|---|---|
Resistant | Susceptible | ||
This study | 27 | 11 | 16 |
Previous report | 12 | 4 | 8 |
Company/supplier | 78 | 2 | 76 |
Total | 117 | 17 | 100 |
R. pseudosolanacearum strain SL882, classified as race 1, biovar 4, and phylotype I (Lee et al. 2011), was cultured on casamino acid-peptone-glucose (CPG) medium (casamino acid, 1 g; peptone, 10 g; glucose, 5 g; and agar, 15 g per liter of distilled water) and incubated at 28°C for 48 h (Kelman 1954). The bacterial culture from the Petri dish (90 × 15 mm) was rinsed with distilled water and washed using a cotton swab to make the inoculum suspension and its concentration was adjusted to approximately 108 CFU/ml (OD600~0.1) using a NanoDrop 2000/UV-Vis spectrophotometer (ThermoFisher Scientific, Waltham, MA, USA).
Seedlings at one month after sowing were inoculated by dipping the roots in the bacterial suspension (Caldwell et al. 2017). Each seedling was grown in a 128-cell seedling tray (28 × 28 × 40, bottom 15 mm) was pulled out, and its roots were dipped in the bacterial suspension. The inoculated seedlings were transplanted into 50-cell seedling trays (45 × 45 × 50, bottom 32 mm) and kept in a growth chamber (temperature = 28°C, relative humidity = 70%, and 16–8 h light-dark cycles). The disease severity was evaluated based on a disease scale of 0 to 4, where 0 = no visible symptoms; 1 = 25% of leaves wilting; 2 = 50% of leaves wilting; 3 = 75% leaves wilting; 4 = all foliage is wilted, and the plant dies (Morel et al. 2018). The disease scale was determined based on visual observation of the degree of wilting. The average value of disease severity for ten plants was calculated per each line. Cultivars with mean disease severity scores of <2 were classified as resistant, while those with scores >2 were classified as susceptible to bacterial wilt.
Genomic DNA extractionGenomic DNA of 27 cultivars which were phenotyped in this study was extracted from young leaf tissues with a modified cetyltrimethylammonium bromide (CTAB) method (Murray and Thompson 1980). The concentration and quality of genomic DNA were measured using NanoDrop 2000/UV–Vis spectrophotometer (ThermoFisher Scientific, Waltham, MA, USA).
Sequence analysis of Bwr-6 region and marker developmentNon-synonymous SNPs between resistant and susceptible groups of tomato cultivars near Bwr-6 within 15 candidate genes (18 SNPs) were previously identified by whole-genome resequencing (Kim et al. 2018). Each candidate gene contained one non-synonymous SNP except two genes, Solyc06g051110.1 and Solyc06g051140.2, which have two and three SNPs, respectively. The nucleotide changes between bacterial wilt-resistant and wilt-susceptible groups of tomato varieties for these SNPs, along with the respective amino acid changes, are presented in Table 2. To dissect Bwr-6 and develop diagnostic markers, selected SNPs were converted to CAPS/dCAPS markers. Based on six non-synonymous SNPs located at 24669159, 34389374, 34399541, 35950028, 37049726, and 37186202, respectively) in Bwr-6 (Kim et al. 2018), five CAPS (RsR6-1~RsR6-5) and one dCAPS (RsR6-6) markers were developed. An HRM molecular marker (KHU-1), tightly linked to bacterial wilt resistance on chromosome 12, was developed based on an SNP (A/G) located at 2,941,301 bp (Kim et al. 2018). We converted this HRM marker, KHU-1, to a dCAPS marker (RsR12-1) to trace Bwr-12 in this study (Table 3). The sequence of the target genes was retrieved from the Sol Genomics Network (https://solgenomics.net/). dCAPS finder 2.0 (http://helix.wustl.edu/dcaps/dcaps.html) NEBcutter (http://nc2.neb.com/NEBcutter2/) were used to find the appropriate restriction enzymes for the respective SNP site.
Candidate gene | SNP position (bp) | Nucleotide change | Amino acid change | Gene annotation (ITAG2.4) | ||
---|---|---|---|---|---|---|
Susceptible | Resistant | Susceptible | Resistant | |||
Solyc06g035530.2 | 24,482,686 | T | C | F | L | Gibberellin 20-oxidase-2 |
Solyc06g035620.2 | 24,667,815 | A | G | Y | C | Scarecrow-like 1 transcription factor |
Solyc06g035630.1 | 24,669,159 | T | C | L | P | GRAS family transcription factor |
Solyc06g036060.2 | 25,438,944 | A | G | I | V | Zinc finger family protein |
Solyc06g048580.1 | 31,287,788 | T | C | L | S | Unknown protein |
Solyc06g051110.1 | 34,213,688 | G | T | L | F | Unknown protein |
Solyc06g051110.1 | 34,214,085 | C | T | A | V | Unknown protein |
Solyc06g051140.2 | 34,285,452 | A | C | E | A | Ubiquitin-conjugating enzyme 22 |
Solyc06g051140.2 | 34,285,469 | G | A | M | I | Ubiquitin-conjugating enzyme 22 |
Solyc06g051140.2 | 34,285,734 | C | A | A | D | Ubiquitin-conjugating enzyme 22 |
Solyc06g051150.1 | 34,291,655 | C | A | L | I | Pentatricopeptide repeat-containing protein |
Solyc06g051190.1 | 34,389,374 | A | G | T | A | RNA-dependent RNA polymerase family protein |
Solyc06g051210.2 | 34,399,541 | C | T | P | S | Bromodomain-containing protein |
Solyc06g053210.2 | 35,950,028 | T | C | X | Q | Ubiquitin | Homeobox leucine zipper protein |
Solyc06g054000.1 | 36,868,039 | G | A | V | I | Unknown protein |
Solyc06g054200.1 | 37,021,520 | A | G | R | S | Calmodulin protein kinase |
Solyc06g054230.2 | 37,049,726 | T | G | D | E | Calmodulin protein kinase |
Solyc06g054400.2 | 37,186,202 | T | C | I | T | Translation initiation factor |
Solyc12g009690.1 | 2,941,301 | A | G | H | R | LRR receptor-like serine/threonine-protein kinase |
Marker name | SNP position (bp) | Marker type | Primer sequence (5ʹ→3ʹ) | Tm (°C) | Restriction enzyme | Expected size (bp) | ||
---|---|---|---|---|---|---|---|---|
Susceptible | Resistant | |||||||
RsR6-1 | 24,669,159 | CAPS | F: | GGAAATATTGGTTACAATCCAGTG | 57.5 | MnlI | 227 | 173, 54 |
R: | GAATACAACAAATCACTACCGGTC | 59.3 | ||||||
RsR6-2 | 34,389,374 | CAPS | F: | CTTCTTGATAGGACGACGTGATAT | 59.3 | RsaI | 87, 116 | 203 |
R: | CAATCAACGGATCACCCATTTTTC | 59.3 | ||||||
RsR6-3 | 34,399,541 | CAPS | F: | CTCTTTTTGCCAGATCTTGAATAG | 57.5 | MnlI | 214 | 116, 98 |
R: | CCATAGGTCAGCATCAAATTTCAA | 57.5 | ||||||
RsR6-4 | 35,950,028 | CAPS | F: | GTTTTCCTTGCAAATCATTTTGGC | 57.5 | MseI | 116, 97 | 213 |
R: | GTATATGTTGAGTTCACAATTCCC | 57.5 | ||||||
RsR6-5 | 37,049,726 | CAPS | F: | CTCAGAAACTGGATAAACTCGAAG | 59.3 | HinfI | 204 | 129, 75 |
R: | GGAGAAAGCAGCCAGCCATTTTT | 60.6 | ||||||
RsR6-6 | 37,186,202 | dCAPS | F: | CGGTGATGAGCAGGATTGATAAAA | 59.3 | HpyCH4III | 234 | 200, 34 |
R: | AGTCTTGGCCTTTGACGTGAAAGTGACACAAGAAG | 60.6 | ||||||
RsR12-1 | 2,941,301 | dCAPS | F: | GTTACACGAACAAGCTTAAATTTCTAGATTTATCCC | 58.8 | AciI | 203 | 168, 35 |
R | GTAATCAATTCGAAGGACCTGTC | 64.9 |
PCR reactions were carried out according to the manufacturer’s instructions (SolGent Co., Ltd., Daejeon, Korea) in a total volume of 25 μl containing 1 μl genomic DNA, 2.5 μl 10X e-Tag reaction buffer, 0.5 μl of 10 mM dNTP mix, 1 μl of each forward and reverse primers, 0.125 μl Solg e-Taq DNA polymerase, and 18.875 μl of ddH2O. PCR amplification was carried out using a Bio-Rad T100 thermocycler (Bio-Rad Laboratories, Inc.) with the following conditions: denaturing for 3 min at 95°C, followed by 34 cycles of 30 s at 95°C denaturation, 30 s at annealing temperature (which varied for different primer sets (Table 3)), 1 min at 72°C extension, and a final elongation step at 72°C for 5 min. PCR products were digested with the respective restriction enzymes. The reaction mixture consisted of 5 μl template PCR product, 1 μl reaction buffer, 0.1 μl of restriction enzyme, and 3.9 μl ddH2O. The mixture was incubated at 37°C for 16 hrs and was carried out using Bio-Rad T100 thermal cycler (Bio-Rad Laboratories, Inc.). The digested product was mixed with 2 μl 6X DNA loading buffer and subjected to gel electrophoresis on a 3% agarose gel to visualize the polymorphic DNA bands. The details of CAPS/dCAPS marker information, including primer sequences, restriction enzymes, and the expected band sizes for resistant and susceptible tomato groups, are presented in Table 3.
Selection of F3 and marker validationTo confirm the efficiency of the two QTLs, F2 generation developed from a cross between E6203 (susceptible) and Hawaii 7998 (resistant) were screened with two markers RsR6-5 (Bwr-6) and RsR12-1 (Bwr-12). Ten F2 plants homozygous for both markers (five resistant and five susceptible), four F2 plants harboring only RsR6-5, and four F2 plants harboring only RsR12-1 were selected for harvesting the F3 seeds. The F3 generation were inoculated and evaluated for disease resistance. Ten plants of each F3 progeny were inoculated. The resistant and susceptible parents were included as controls during disease evaluation. The mean disease severity of ten plants was used to designate the resistance level of each progeny. Furthermore, the markers were evaluated for their diagnostic value using tomato cultivars (Yang et al. 2015).
The resistance level of 27 cultivars against bacterial wilt disease was inoculated and confirmed in this study. Hawaii 7996, B-Blocking, Shincheonggang, BWR-20, Spider, High Power, 10-BA-3-33, 10-BA-4-24, IT 201664, Hawaii 7998, and Fighting were resistant with low mean disease severity scores (<2) while UC-134, LA1589, Purple Calabash, Florida8516, Heinz 1706, A-1, E6203, Moneymaker, Super Dotaerang, Anahu, Red Strong, Bluck Plum, Gold Nugget, VF-36, M82, and Dotaerang Red were susceptible with a mean disease severity score of >2.0 (Fig. 1). Previously reported susceptible cultivars, such as Moneymaker, Heinz 1706, and Super Dotaerang (Han et al. 2009, Kim et al. 2018), showed the highest mean disease severity indicating the presence of adequate disease pressure on inoculated plants. In addition, phenotypic information of 12 cultivars was inferred from previous reports and 78 cultivars were received from respective companies. In total, the collection includes 17 resistant and 100 susceptible cultivars (Table 1).
Mean disease severity of tomato cultivars infected with R. pseudosolanacearum strain SL882, one month after inoculation. The disease severity was scored based on the disease scale of 0 to 4, where 0 = no visible symptoms; 1 = 25% of leaves wilting; 2 = 50% of leaves wilting; 3 = 75% of leaves wilting; and 4 = all foliage is wilted, and the plant dies.
The primer sets were tested for polymorphism using six susceptible and five resistant cultivars. All primer sets resulted in a clear polymorphic band between susceptible and resistant tomato cultivars after enzyme digestion (Fig. 2). For selecting an accurate and reliable marker for Bwr-6, six markers (RsR6-1~6-6) were screened using 117 tomato cultivars (Table 4). The Bwr-12 genotype of these cultivars was determined using RsR12-1. Hawaii 7996, Hawaii 7998, 10-BA-3-33, 10-BA-4-24, BWR-20, BWR-1, BWR-22, and BWR-23 had homozygous resistant genotypes with all six markers near Bwr-6. IT201664, High Power, Spider, SVTX6258, Super High Power, B-Blocking, Shincheonggang, Fighting, and Geumgang are either susceptible or heterozygous to RsR6-1, RsR6-2, and RsR6-3 while they are resistant to RsR6-4, RsR6-5, and RsR6-6 except for High Power, which is heterozygous to RsR6-4 and RsR6-5. Comparing the six markers based on the genotype of resistant cultivars, RsR6-1, RsR6-2, and RsR6-3 did not seem better candidate markers for Bwr-6 because IT 201664, High Power, Spider, and Super High Power are susceptible to these markers.
Marker analysis of the newly developed CAPS/dCAPS in 11 tomato genotypes near Bwr-6 (RsR6-1~RsR6-6) and Bwr-12 (RsR12-1). Lanes 1 to 6 represent the bacterial wilt-susceptible group (M82, E6203, UC-134, VF-36, Gold Nugget, and Moneymaker) and lanes 7 to 11 represent the bacterial wilt-resistant group (Hawaii 7996, Hawaii 7998, 10-BA-3-33, 10-BA-4-24, and BWR-20). M, 1 kb DNA size marker. The arrows on the left indicate a 200 bp standard fragment.
No. | Tomato cultivar/line | Type | Company/Suppliera | SNP marker genotypeb | Referencec | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|
RsR6-1 | RsR6-2 | RsR6-3 | RsR6-4 | RsR6-5 | RsR6-6 | RsR12-1 | |||||
Bacterial wilt resistant | |||||||||||
1 | Hawaii 7996 | Inbred line | KHU | R | R | R | R | R | R | R | This study |
2 | Hawaii 7998 | Inbred line | KHU | R | R | R | R | R | R | R | This study |
3 | 10-BA-3-33 | Inbred line | KHU | R | R | R | R | R | R | R | This study |
4 | 10-BA-4-24 | Inbred line | KHU | R | R | R | R | R | R | R | This study |
5 | BWR-20 | Inbred line | KHU | R | R | R | R | R | R | R | This study |
6 | BWR-1 | Inbred line | KHU | R | R | R | R | R | R | R | Kim et al. (2018) |
7 | BWR-22 | Inbred line | KHU | R | R | R | R | R | R | R | Kim et al. (2018) |
8 | BWR-23 | Inbred line | KHU | R | R | R | R | R | R | R | Kim et al. (2018) |
9 | IT 201664 | Inbred line | RDA | H | S | S | R | R | R | R | This study |
10 | High Power | F1 hybrid | Dae Yeon seed co. | H | S | S | H | H | R | R | This study |
11 | Super High Power | F1 hybrid | Dae Yeon seed co. | S | S | S | R | R | R | R | Kim et al. (2018) |
12 | Spider | F1 hybrid | Takii Korea | S | S | S | R | R | R | R | This study |
13 | B-Blocking | F1 hybrid | Takii Korea | H | H | H | R | R | R | H | This study |
14 | Shincheonggang | F1 hybrid | Farm Hannong | H | H | H | R | R | R | H | This study |
15 | Fighting | F1 hybrid | Takii Korea | H | H | H | R | R | R | H | This study |
16 | SVTX6258 | F1 hybrid | Monsanto Korea | H | H | H | R | R | R | R | Supplier |
17 | Geumgang | F1 hybrid | Monsanto Korea | H | H | H | R | R | R | H | Supplier |
Bacterial wilt susceptible | |||||||||||
18 | M82 | Inbred line | TGRC | S | S | S | S | S | S | S | This study |
19 | E6203 | Inbred line | TGRC | S | S | S | S | S | S | S | This study |
20 | VF36 | Inbred line | TGRC | R | S | S | S | S | S | S | This study |
21 | Moneymaker | Inbred line | TGRC | R | S | S | S | S | S | S | This study |
22 | A-1 | Inbred line | TGRC | S | S | S | S | S | S | S | This study |
23 | Anahu | Inbred line | TGRC | S | S | S | S | S | R | S | This study |
24 | Black Plum | Inbred line | TGRC | S | S | S | R | S | R | S | This study |
25 | Florida 8516 | Inbred line | TGRC | S | S | S | S | S | S | S | This study |
26 | Gold Nugget | Inbred line | TGRC | R | S | S | S | S | R | S | This study |
27 | Purple Calabash | Inbred line | TGRC | H | H | S | H | S | R | S | This study |
28 | UC-134 | Inbred line | TGRC | S | S | S | S | S | S | S | This study |
29 | Heinz 1706 | Inbred line | TGRC | S | S | S | S | S | S | S | This study |
30 | New Yorker | Inbred line | TGRC | S | S | S | S | S | R | S | Jung et al. (2014) |
31 | Miniheuksu | F1 hybrid | Asia Seed Co., Ltd. | R | S | S | S | S | R | S | Kim et al. (2018) |
32 | TY Unique | F1 hybrid | Asia Seed Co., Ltd. | R | S | S | S | H | R | H | Supplier |
33 | Shinsugar Yellow | F1 hybrid | Asia Seed Co., Ltd. | H | H | H | H | H | R | H | Supplier |
34 | Red Strong | F1 hybrid | Bunong Seed | R | S | R | R | S | R | R | This study |
35 | Sun Star | F1 hybrid | Bunong Seed | H | S | S | S | S | H | S | Supplier |
36 | Black Eagle | F1 hybrid | Bunong Seed | R | S | H | H | S | S | H | Supplier |
37 | Tamla | F1 hybrid | Bunong Seed | H | S | H | H | S | H | S | Supplier |
38 | Bntoskna | F1 hybrid | Bunong Seed | H | S | S | H | S | H | H | Supplier |
39 | TY Izzang | F1 hybrid | Bunong Seed | H | S | H | H | S | R | S | Supplier |
40 | Candy Plus | F1 hybrid | Bunong Seed | H | S | H | H | S | H | S | Supplier |
41 | Super Star | F1 hybrid | Bunong Seed | H | S | S | S | S | H | H | Supplier |
42 | Red Zenith | F1 hybrid | Bunong Seed | H | S | H | H | S | R | H | Supplier |
43 | TY Hunter | F1 hybrid | Bunong Seed | R | S | H | H | S | H | H | Supplier |
44 | TY One Top | F1 hybrid | Bunong Seed | S | S | S | H | S | H | H | Supplier |
45 | Yureka | F1 hybrid | Bunong Seed | S | S | H | H | S | H | H | Supplier |
46 | Black Ace | F1 hybrid | Bunong Seed | R | S | S | S | S | S | H | Supplier |
47 | Oasis | F1 hybrid | Bunong Seed | S | S | H | S | S | H | S | Supplier |
48 | TY Megaton | F1 hybrid | Bunong Seed | S | S | S | S | S | R | S | Supplier |
49 | Gold Sugar | F1 hybrid | Bunong Seed | R | S | H | H | R | R | S | Supplier |
50 | Dotaerang Red | F1 hybrid | Dong seo seed | H | H | S | S | S | R | S | This study |
51 | Black Ball | F1 hybrid | Dongoh Seed | R | S | H | H | S | H | H | Supplier |
52 | Kolmi | F1 hybrid | Dongoh Seed | S | S | H | R | H | R | S | Supplier |
53 | Sky Ball | F1 hybrid | Dongoh Seed | R | S | H | H | S | R | R | Supplier |
54 | Starbuck | F1 hybrid | Farm Hannong | R | S | S | S | S | S | S | Supplier |
55 | Olleh TY | F1 hybrid | Farm Hannong | S | S | H | H | S | R | S | Supplier |
56 | Rafito | F1 hybrid | Farm Hannong | H | S | S | S | S | H | R | Supplier |
57 | Big Wonderful | F1 hybrid | Gonong Seed | R | S | H | H | S | H | S | Supplier |
58 | TY Carnival | F1 hybrid | Gyeongwon | R | S | H | H | S | H | H | Supplier |
59 | Legend Summer | F1 hybrid | Haesung Seed Plus | S | S | S | S | S | R | S | Supplier |
60 | Daewang | F1 hybrid | Jeil Seed Bio | H | S | S | S | S | S | S | Supplier |
61 | Hongboseok | F1 hybrid | Jeil Seed Bio | S | S | S | S | S | S | S | Supplier |
62 | Jeilheukjinju | F1 hybrid | Jeil Seed Bio | S | S | S | R | S | R | S | Supplier |
63 | Dotaerang Myeongpum | F1 hybrid | Jeil Seed Bio | S | S | S | S | S | S | S | Supplier |
64 | Heukryong | F1 hybrid | Jeil Seed Bio | S | S | S | S | S | S | H | Supplier |
65 | Minijaok | F1 hybrid | Jeil Seed Bio | S | S | S | R | S | R | S | Supplier |
66 | Sinheukjinju | F1 hybrid | Jeil Seed Bio | H | R | S | R | R | R | S | Supplier |
67 | Super Dotaerang | F1 hybrid | Koregon seed | H | H | S | S | S | R | S | This study |
68 | Lezaforta F1 | F1 hybrid | Mifko seed | S | S | S | S | S | H | S | Supplier |
69 | Unicorn | F1 hybrid | Monsanto Korea | H | S | S | S | S | R | S | Supplier |
70 | SV 7160 TC | F1 hybrid | Monsanto Korea | R | S | H | R | R | R | S | Supplier |
71 | SV02444 TG | F1 hybrid | Monsanto Korea | S | S | H | S | S | R | H | Supplier |
72 | Bacchus | F1 hybrid | Monsanto Korea | R | S | H | H | S | R | S | Supplier |
73 | SV4224 TH | F1 hybrid | Monsanto Korea | S | S | H | S | S | R | H | Supplier |
74 | SV0339TG | F1 hybrid | Monsanto Korea | S | S | H | S | S | R | H | Supplier |
75 | Tiara | F1 hybrid | Nongwoo Bio Co., Ltd. | H | S | S | S | S | R | S | Kim et al. (2018) |
76 | TY SenseQ | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | H | H | S | H | S | Supplier |
77 | Redpang | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | S | S | S | R | S | Supplier |
78 | Titichal | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | H | H | S | R | S | Supplier |
79 | TY Altorang | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | H | H | S | R | S | Supplier |
80 | Beta Tiny | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | S | H | H | H | S | Supplier |
81 | TY Tiny | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | H | H | S | H | S | Supplier |
82 | Cupirang | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | S | S | S | R | S | Supplier |
83 | Minichal | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | S | H | H | H | S | Supplier |
84 | TY Sispen | F1 hybrid | Nongwoo Bio Co., Ltd. | S | S | H | H | S | R | S | Supplier |
85 | Black Change | F1 hybrid | Nongwoo Bio Co., Ltd. | H | S | H | H | S | H | H | Supplier |
86 | Mulya | F1 hybrid | RDA | S | S | S | S | S | R | S | Kim et al. (2018) |
87 | Sigyo 1 ho | F1 hybrid | RDA | S | S | S | S | S | R | S | Kim et al. (2018) |
88 | Broadley | F1 hybrid | RDA | R | S | S | S | S | R | S | Kim et al. (2018) |
89 | Yulwon | F1 hybrid | RDA | S | H | H | R | H | R | S | Kim et al. (2018) |
90 | Hoyong | F1 hybrid | Sakata Korea | H | S | S | S | S | R | H | Kim et al. (2018) |
91 | Tosama | F1 hybrid | Sakata Korea | S | S | S | H | S | H | S | Supplier |
92 | Super Sun Road | F1 hybrid | Sakata Korea | H | S | S | H | S | H | H | Supplier |
93 | Super Top | F1 hybrid | Sakata Korea | S | S | S | S | S | H | H | Supplier |
94 | Lokousan Maru | F1 hybrid | Sakata Korea | S | S | S | S | S | H | S | Supplier |
95 | Taiyau | F1 hybrid | Sakata Korea | S | S | S | S | S | H | H | Supplier |
96 | Taihu | F1 hybrid | Sakata Korea | R | S | H | S | S | H | H | Supplier |
97 | Super Top | F1 hybrid | Sakata Korea | S | S | S | S | S | H | H | Supplier |
98 | Tiger | F1 hybrid | Samsung Seeds | H | S | S | S | S | R | S | Supplier |
99 | Chalstone TY | F1 hybrid | Sky seed | S | S | H | H | S | R | H | Supplier |
100 | TY Marathon | F1 hybrid | Sky seed | S | S | S | S | S | R | H | Supplier |
101 | Rapsody | F1 hybrid | Syngenta Korea | H | S | S | H | S | S | S | Supplier |
102 | Madison | F1 hybrid | Syngenta Korea | H | S | S | S | S | S | S | Supplier |
103 | Ricophin-9 | F1 hybrid | Syngenta Korea | H | S | S | H | H | R | S | Supplier |
104 | Duine | F1 hybrid | Syngenta Korea | H | S | S | S | S | H | S | Supplier |
105 | Dafnis | F1 hybrid | Syngenta Korea | S | S | H | H | S | H | H | Supplier |
106 | Komodo | F1 hybrid | Syngenta Korea | S | S | H | H | S | H | H | Supplier |
107 | Tory | F1 hybrid | Syngenta Korea | S | S | H | H | S | H | H | Supplier |
108 | Mamirio | F1 hybrid | Syngenta Korea | R | S | H | H | S | H | H | Supplier |
109 | European Rapsodie | F1 hybrid | Syngenta Korea | H | S | S | H | S | S | S | Supplier |
110 | Trio Plus | F1 hybrid | Taeyang seed | R | S | H | H | S | R | H | Supplier |
111 | Kang Jeok | F1 hybrid | Taeyang seed | R | S | H | H | S | H | S | Supplier |
112 | Dotaerang TY Winner | F1 hybrid | Takii Korea | S | S | H | S | S | R | S | Supplier |
113 | Doterang Plus | F1 hybrid | Takii Korea | S | S | S | S | S | R | S | Supplier |
114 | Dotaerang Solar | F1 hybrid | Takii Korea | S | S | S | S | S | R | H | Supplier |
115 | Cuty | F1 hybrid | Takii Korea | S | H | H | H | H | R | S | Supplier |
116 | Dotaerang Diamond | F1 hybrid | Takii Korea | S | S | S | S | S | R | S | Supplier |
117 | LA1589 | Wild species | TGRC | R | S | R | R | R | R | S | This study |
a Company/supplier of seed or DNA sample: KHU = Kyung Hee University, RDA = National Agrobiodiversity Center (RDA-Genebank), TGRC = Tomato Genetics Resource Center. All tomato genotypes belong to Solanum lycopersicum species except LA1589 (Solanum pimpinellifolium).
b SNP marker genotype: R = resistant, S = susceptible, H = heterozygous.
c Reference for the phenotypic information of the tomato cultivars used in the study. Supplier’s phenotypic information was obtained via personal communication from the company.
Therefore, we considered RsR6-4, RsR6-5, and RsR6-6 for further analysis using the susceptible set of cultivars. Almost all bacterial wilt-susceptible cultivars had susceptible or heterozygous genotypes to RsR6-5 except Gold Sugar, Sinheukjinju, SV7160TC, and LA1589. Red Strong shows the resistant genotype to RsR6-4 and the susceptible genotype to RsR6-5. Cultivars resistant to RsR6-6 and heterozygous with RsR12-1, such as SV02444 TG, SV4224 TH, and SV0339TG, are expected to be resistant to bacterial wilt, but all were susceptible. In addition, Red Strong and SkyBall have homozygous resistant genotypes with RsR6-6 and RsR12-1, although both exhibit susceptible phenotype (Table 4).
The presence of resistant alleles in both Bwr-6 and Bwr-12 resulted in resistant phenotype while the absence of either of the two resulted in susceptible phenotype. This marker analysis indicated that cultivars homozygous resistant to RsR6-5 and either a homozygous resistant or heterozygous genotype to RsR12-1 exhibited a resistant phenotype. However, homozygous susceptible or heterozygous genotypes with RsR6-5 exhibited a susceptible phenotype regardless of the RsR12-1 genotype. One exception to this was High Power, which exhibits a resistant phenotype although it has a heterozygous genotype with RsR6-5. The genotype of RsR6-5 is highly correlated with the bacterial wilt phenotype of tomato cultivars and selected for tracing Bwr-6. The diagnostic accuracy of the markers was evaluated and RsR6-5 and RsR12-1 combination was resulted in 94.1% true positive rate and 100% true negative rate (Supplemental Table 1). Taken together, RsR6-5 and RsR12-1 should be used for effective marker-assisted selection of bacterial wilt resistance in tomatoes
Validation of the RsR6-5 and RsR12-1 markers using F3 populationsTo validate the efficiency of the two markers for selecting resistant lines in the segregating population, 10 F3 generations homozygous resistant and susceptible (each five lines), and eight F3 generation carrying only Bwr-6 or Bwr-12 (each four lines) were developed from E6203 (susceptible) and Hawaii 7998 (resistant) were selected based on RsR6-5 and RsR12-1 genotype. Hawaii 7998 is one of the entries in international set of bacterial wilt resistant lines evaluated in twelve fields and showed an average of 90% survival rate (Wang et al. 1998). The resistance of Hawaii 7998 and Hawaii 7996 was derived from the same origin, PI 127805A (S. pimpinellifolium) (Daunay et al. 2010, Scott et al. 2005). Hence, genetic resistance to bacterial wilt in these two lines might be governed by same gene.
Ten plants of each F3 generation, along with the two parental lines, were inoculated with R. pseudosolanacearum strain SL882 for disease evaluation. The two parental lines exhibited distinct differences in bacterial wilt resistance as expected, with mean disease scores of 4.0 ± 0.0 and 1.2 ± 0.61, respectively. The mean disease severity between homozygous resistant and susceptible genotypes in the F3 generation was significantly different (Fig. 3A). Lines with homozygous resistant genotypes exhibited highly resistant phenotypes with mean disease severity scores ranging from 0.4 ± 0.40 to 1.6 ± 0.65, which were not significantly different from that of the resistant parent. On the other hand, homozygous susceptible F3 exhibited highly susceptible phenotypes with mean disease severity scores of 4.0 ± 0.0, which were similar to the susceptible parent (Fig. 3A, 3B). F3 progenies carrying resistant allele only in RsR6-5 or RsR12-1 were susceptible (Supplemental Fig. 1). These results suggest that the combination of RsR6-5 and RsR12-1, which are associated with Bwr-6 and Bwr-12, respectively, are predictive for bacterial wilt resistance.
Evaluation of F3 generation selected by RsR6-5 and RsR12-1 for resistance to R. pseudosolanacearum. (A) Mean disease severity of F3 generation developed from a crossing between E6203 and Hawaii 7998 one month after inoculation with R. pseudosolanacearum strain SL882. The genotypes of the F3 lines are either homozygous resistant or susceptible at RsR6-5 and RsR12-1. The F3 lines are followed by RR (carrying resistant alleles of both Bwr6 and Bwr-12) or SS (carrying susceptible alleles at both loci). Different letters on the bars indicate a significant difference (P < 0.05) of mean disease severity. (B) Photographs of homozygous F3 generation and parental lines one month after inoculation with SL882.
Genetic resistance is the most effective control strategy for bacterial wilt of tomato, and multiple breeding programs have been engaged in developing resistant lines by incorporating resistance from different resistant sources (Daunay et al. 2010, Wang et al. 1998). Genomic regions linked to bacterial wilt resistance in the well-known resistant cultivar, Hawaii 7996, were identified on different chromosomes, and some regions were detected against specific strains (Thoquet et al. 1996a, 1996b, Wang et al. 2000). Bwr-3 and Bwr-4 were associated with resistance against phylotype II strains, while Bwr-12 was specific to Phylotype-I (Carmeille et al. 2006). In contrast, Bwr-6 was associated with various strains from different phylotypes (I and II) (Thoquet et al. 1996a) and consistently detected under various conditions (Carmeille et al. 2006, Geethanjali et al. 2010, Mangin et al. 1999, Wang et al. 2013). Genetic analysis of bacterial wilt resistance in S. lycopersicum var. cerasiforme ‘L285’ identified three QTL on chromosomes 6, 7, and 10 (Danesh et al. 1994). The genomic regions associated with bacterial wilt resistance on chromosome 6 in both Hawaii 7996 and L285 were co-localized. Furthermore, the most stable QTL for bacterial wilt resistance in eggplant was also identified in chromosome 6 and syntenic with Bwr-6 of tomato (Salgon et al. 2018). All these suggest the significance of Bwr-6 in bacterial wilt resistance in tomato and likely in other Solanaceae crops. The development of functional markers for such broad-spectrum QTL is essential to facilitate marker-assisted breeding.
Genetic analysis using only segregating populations is time-consuming, cost-inefficient, and the developed markers may be specific to certain resistant lines (Pascual et al. 2016). In this regard, we used a complementary approach to dissect the Bwr-6 region, using germplasm collections and F3 generation to validate our result. SNP-based CAPS/dCAPS markers near Bwr-6 were developed and validated. A total of 117 tomato germplasms were screened with newly developed markers for Bwr-6 genotypes, and the corresponding phenotypic information was used to explore the efficiency of each marker.
Among 17 resistant cultivars used to screen the markers in the Bwr-6 region, four cultivars (IT 201664, High Power, Super High Power, and Spider) had susceptible genotypes with RsR6-1, RsR6-2, and RsR6-3. On the other hand, all resistant cultivars had homozygous resistant genotypes with RsR6-4, RsR6-5, and RsR6-6 except High Power, which has a heterozygous genotype with RsR6-4 and RsR6-5. Similarly, all cultivars exhibiting resistant phenotypes have homozygous resistant genotypes with RsR12-1, except B-Blocking, Shincheonggang, Fighting, and Geumgang which are heterozygous to this marker. The true positive rate of RsR6-1, RsR6-2, and RsR6-3 combined with RsR12-1 was 47.1% whereas that of RsR6-4 and RsR6-5 was 94.1% and that of RsR6-6 was 100%. Based on these observations, we hypothesized that the three markers (RsR6-4, RsR6-5, and RsR6-6) are better predictors of resistance conferred by Bwr-6 than RsR6-1, RsR6-2, and RsR6-3.
The markers were further compared using the susceptible panel of germplasms. Some cultivars exhibiting susceptible phenotypes, such as Red Strong, are resistant to RsR6-4, RsR6-6, and RsR12-1. However, all cultivars having homozygous resistant genotypes with RsR6-5 and RsR12-1 exhibited resistant phenotypes. In summary, 99% and 87% true negative rate was obtained for RsR6-4 and RsR6-6 while RsR6-5 resulted in 100% true negative rate genotypes in combination with RsR12-1 (Supplemental Table 1). These results suggest that RsR6-5 is the best diagnostic marker to trace Bwr-6 associated with bacterial wilt resistance. Diagnostic markers developed based only on a segregating population may not fully correlate with the trait when tested in diverse germplasms, hindering their utilization for marker-assisted selection in a broad set of breeding germplasms (Niewohner et al. 1995). Therefore, utilization of a wide range of germplasms, including inbred lines, commercial F1 hybrids, and wild species, for validating developed markers is essential before deployment to end-users, including breeders and farmers.
The diagnostic potential of RsR6-5 coupled with RsR12-1 for bacterial wilt resistance in tomato was tested using a broader set of germplasms (Bartkiewicz et al. 2018, Yang et al. 2015) and can be used for marker-assisted selection in commercial breeding programs. Without determining the genotyping results for another major resistance QTL, Bwr-12, determination of Bwr-6 genotype with the RsR6-5 marker alone is able to predict resistant and susceptible phenotypes in the tomato cultivars used in this study with 94.1% and 96% accuracy, respectively. High Power, SV7160TC, Gold Sugar, Sinheukjinju, and LA1589 showed non-matching genotypes. On the other hand, Bwr-12 genotyping determined by RsR12-1 alone was able to predict resistant and susceptible phenotypes with 100% and 66% accuracy, respectively. Resistance conferred by the Bwr-12 genotype shows a dominant inheritance pattern (Kim et al. 2018). RsR6-5 and RsR12-1 combination was resulted in 94.1% of true positive rate and 100% true negative rate showing the highest diagnostic accuracy compared to other marker combinations.
Heterozygous to RsR6-5 and either heterozygous or homozygous resistant to RsR12-1 yielded susceptible phenotypes except for High Power, which suggests that the Bwr-6 resistance allele might be recessive; however, this should be further validated using a segregating population. Recessive gene resistance to bacterial disease has been reported in Arabidopsis and rice. R. solanacearum resistance in Arabidopsis thaliana is governed by a recessively inherited gene (RRS1-R) (Deslandes et al. 2002). Similarly, among more than 43 resistance genes identified so far for bacterial blight in rice caused by Xanthomonas oryzae pv. oryzae, 16 are inherited as recessive traits (Kim 2018, Vikal and Bhatia 2017).
Extensive QTL mapping studies have been conducted to identify genomic regions associated with bacterial wilt resistance (Carmeille et al. 2006, Geethanjali et al. 2010, Mangin et al. 1999, Thoquet et al. 1996a, 1996b, Wang et al. 2000, 2013). The release of the tomato reference genome (Tomato Genome Consortium 2012) and availability of whole-genome resequencing data for various tomato cultivars (Lin et al. 2014, The 100 Tomato Genome Sequencing Consortium 2014) facilitated comparisons among tomato genotypes using genome-wide SNP markers for different traits. The adequate number of SNPs among tomato genotypes can be used to saturate markers nearby previously identified QTL regions. Accordingly, whole-genome resequencing of bacterial wilt resistant and susceptible tomato cultivars revealed genome-wide SNPs that were candidates for distinguishing the two groups of tomato, with the highest number of non-synonymous SNPs identified on chromosomes 12 and 6 (Kim et al. 2018). Analysis of SNPs near Bwr-12 in the same study discovered molecular marker (KHU-1) tightly linked to bacterial wilt resistance; this marker was used to discriminate resistant and susceptible tomato cultivars. The analysis of SNPs near Bwr-6 and the development of diagnostic markers in this study will pave the way toward identifying candidate genes and facilitating resistance gene pyramiding.
RsR6-5 is located in the coding region of Solyc06g054230.2, which encodes a putative calmodulin protein kinase, suggesting that this gene could be a possible candidate gene for bacterial wilt resistance in tomato. The nucleotide substitution of guanine to thymine at 128 bp changes the amino acid aspartate to glutamate. The majority of plant disease resistance genes encode nucleotide-binding site-leucine-rich repeat (NBS-LRR) proteins (McHale et al. 2006). However, the role of calmodulin proteins in the response of plants to biotic and abiotic stresses has also been reported (Cheval et al. 2013, Zeng et al. 2015). Calmodulin is a calcium-binding protein and regulates downstream calcium signal-related responses. The expression of tomato calmodulin genes is significantly altered upon pathogen infection. Functional analysis revealed that the silencing of SlCaM2 (Solyc10g081170.1.1) and SlCaM6 (Solyc03g098050.2.1) in tomato reduced its resistance to tomato rattle virus and Pythium aphanidermatum and decreased the expression of downstream signaling and defense-related genes (Zhao et al. 2013). Transcriptome analysis of bacterial wilt-resistant (LS-89) and susceptible (Ponderosa) cultivars indicated an approximately 30-fold increase of a putative calmodulin-binding family protein in response to R. solanacearum infection in resistant cultivars, while the analogous response in susceptible cultivars was very limited (Ishihara et al. 2012). These findings suggest that Solyc06g054230.2 may play an important role in bacterial wilt resistance in tomato; although, further analysis is required to elucidate the gene function.
In conclusion, SNPs near Bwr-6 were analyzed to search for markers tightly linked to this QTL. A total of 117 tomato germplasms were used to validate newly developed markers near this QTL. Among the analyzed markers, RsR6-5 is tightly linked to bacterial wilt resistance derived from Bwr-6. Consequently, this marker, in combination with RsR12-1, effectively predicted bacterial wilt-resistant and susceptible cultivars. The significance of these markers was further validated using an F3 generation developed from a crossing between resistant and susceptible parents. F3 lines that had resistant genotypes with RsR6-5 and RsR12-1 exhibited resistant phenotypes, while susceptible to the same markers exhibited susceptible phenotypes. The SNP-based diagnostic marker to Bwr-6 was not identified in previous studies, and the newly developed marker in this study (RsR6-5) will help to trace this locus in marker-assisted breeding of tomato cultivars that are resistant to the devastating effects of bacterial wilt.
JML conceived and designed the experiments. AMA, JC and YK performed the experiments. CSO, IY, and ISN provided experimental materials. AMA and JML wrote and revised the manuscript, and all co-authors contributed and approved the final draft of the manuscript.
This work was supported by the Golden Seed Project (Center for Horticultural Seed Development) funded by the Ministry of Agriculture, Food and Rural Affairs of Korea.