- Research article
- Open Access
Normal adult survival but reduced Bemisia tabaci oviposition rate on tomato lines carrying an introgression from S. habrochaites
BMC Geneticsvolume 15, Article number: 142 (2014)
Host plant resistance has been proposed as one of the most promising approaches in whitefly management. Already in 1995 two quantitative trait loci (Tv-1 and Tv-2) originating from S. habrochaites CGN1.1561 were identified that reduced the oviposition rate of the greenhouse whitefly (Trialeurodes vaporariorum). After this first study, several others identified QTLs affecting whitefly biology as well. Generally, the QTLs affecting oviposition were highly correlated with a reduction in whitefly survival and the presence of high densities of glandular trichomes type IV. The aim of our study was to further characterize Tv-1 and Tv-2, and to determine their role in resistance against Bemisia tabaci.
We selected F2 plants homozygous for the Tv-1 and Tv-2 QTL regions and did three successive backcrosses without phenotypic selection. Twenty-three F2BC3 plants were phenotyped for whitefly resistance and differences were found in oviposition rate of B. tabaci. The F2BC3 plants with the lowest oviposition rate had an introgression on Chromosome 5 in common. Further F2BC4, F2BC4S1 and F2BC4S2 families were developed, genotyped and phenotyped for adult survival, oviposition rate and trichome type and density. It was possible to confirm that an introgression on top of Chr. 5 (OR-5), between the markers rs-2009 and rs-7551, was responsible for reducing whitefly oviposition rate.
We found a region of 3.06 Mbp at the top of Chr. 5 (OR-5) associated with a reduction in the oviposition rate of B. tabaci. This reduction was independent of the presence of the QTLs Tv-1 and Tv-2 as well as of the presence of trichomes type IV. The OR-5 locus will provide new opportunities for resistance breeding against whiteflies, which is especially relevant in greenhouse cultivation.
Tomato is one of the most important vegetables worldwide. It is host for a broad range of pathogens and pests. Among the pests affecting tomato production whiteflies are the most important in terms of costs and distribution. There are more than 1500 species of whiteflies , of which Bemisia tabaci Group Mediterranean-Middle East-Asia Minor I and Trialeurodes vaporariorum (Westwood) are the biggest threats in commercial tomato production. Bemisia tabaci affects tomato production directly (i.e. phloem consumption, irregular ripening of the fruits) and indirectly (virus transmission) causing yield losses that can range from 50% to 100% of the potential production ,.
Among the possible control methods, host plant resistance has been proposed as one of the most promising for insect pest management ,. Resistance to whiteflies was found in several wild relatives of tomato (Solanum pennellii, S. habrochaites, S. lycopersicum var. cerasiforme, S. pimpinellifolium, S. galapagense) -. In these species, whitefly resistance is associated with the presence of high densities of glandular trichomes (type I, IV and VI) and with the presence of specific secondary metabolites (a.o. 7-epizingiberene, 2-tridecanone, and acyl sugars) -. The species S. habrochaites contains accessions (formerly known as Lycopersicon hirsutum fr. glabratum) that accumulate methyl ketones, of which the synthesis is located in the glandular head of type VIc trichomes -. Solanum habrochaites also contains accessions (formerly known as Lycopersicon hirsutum fr. typicum) that accumulate sesquiterpenes which are synthesised in type IV trichomes . In S. pennellii, S. pimpinellifolium and S. galapagense the synthesis of acyl sugars is associated with the presence in high densities of trichomes type IV ,,,. Although also some accessions of S. cheesmaniae accumulate high levels of acyl sugars, they lack type IV trichomes . The Mi1-2 gene, which confers resistance to several species of root-knot nematodes (Meloidogyne spp.) , plays a role in the resistance against insects, e.g. some isolates of potato aphid (Macrosiphum euphorbiae Thomas) ,, the sweet potato whitefly (B. tabaci)  and the tomato psyllid (Bactericerca cockerelli) . This resistance is independent of the presence of glandular trichomes and acyl sugar concentration .
QTL mapping studies have been carried out to identify genomic regions involved in whitefly resistance. Maliepaard et al.  focused on resistance against the greenhouse whitefly T. vaporariorum (Westwood) from S. habrochaites (CGN1.1561) and identified two QTLs reducing whitefly oviposition rate (Tv-1 on Chr. 1 and Tv-2 on Chr. 12) together with two QTLs related to trichome type IV density (TriIV-1 on Chr. 5 and TriIV-2 on Chr. 9) and one QTL for trichome type VI density (TriVI-1 on Chr. 1). After this first study, others have explored different resistance sources and more QTLs were identified. A summary of the QTLs related to whitefly resistance in tomato is given in Table 1. The use of backcross introgression lines (ILs) was also proposed as a method to identify genomic regions important for whitefly resistance. These ILs helped to identify regions and genes involved in traits related to insect resistance, like the production of monoterpenes, sesquiterpenes and acyl sugars -. However, they failed to identify regions associated to whitefly resistance in terms of adult survival or oviposition rate , supporting the observations from the QTL mapping studies that whitefly resistance is polygenic inherited and possibly epistatic interactions play a role as well. Except the QTLs described by Maliepaard et al. , all other QTLs affecting whitefly oviposition were highly correlated with a reduction in whitefly survival and/or to high densities of trichomes type IV, suggesting that the low oviposition rate is the consequence of a low survival rate ,,. To study resistance mechanisms affecting whitefly oviposition rate exclusively, we focused on the further characterization of the QTLs identified by Maliepaard et al. , and determined their role in resistance against Bemisia tabaci.
Plant materials and growing conditions
The study was based on the F2 offspring population that was created by Maliepaard et al. , it was obtained by self-pollination of a single F1 plant that was derived from a cross between S. lycopersicum (cv. Moneymaker) and S. habrochaites (CGN1.1561). We have sown again individuals of this F2 population and selected plants that were homozygous for either one or both QTLs associated to a reduction in oviposition rate using Cleaved Amplified Polymorphisms (CAPs) markers (Table 2). The selected F2, BC1 and BC2 plants were backcrossed with S. lycopersicum (cv. Moneymaker) for three generations. Plants were chosen containing at least one of the markers flanking the QTLs. The obtained F2BC3 and F2BC4 families were genotyped and phenotyped for adult survival and oviposition rate. Selected F2BC4 plants were selfed to obtain F2BC4S1 plants and F2BC4S2, which were also genotyped and phenotyped. An overview of the material development is shown in Figure 1.
The tomato plants were grown in a greenhouse in Wageningen, The Netherlands (20 ± 2°C, 70% RH, 16/8 h day/night) in 14 cm diameter pots filled with soil compost. The plants were fertilized twice a week with standard fertilizer for tomato and watered once a day. When the plants were five weeks old, they were transferred to an insect proof greenhouse. The greenhouse temperature was increased slowly from 20 to 27°C to allow plants to adapt to the higher temperature (27 ± 2°C, 70% RH, 16/8 h day/night) used during the infestation that took place one week after transfer.
A non-viruliferouswhitefly rearing (Bemisia tabaci Group Mediterranean-Middle East-Asia Minor I) was maintained on the susceptible tomato cultivar Moneymaker at Wageningen UR Plant Breeding, Wageningen, The Netherlands. The initial inoculum was obtained from a rearing at the Laboratory of Entomology, Wageningen UR, Wageningen, The Netherlands.
Whiteflies (four days old) were anesthetized using CO2. Five females were selected under a binocular and put in a clip-on cage (2.5 cm diameter and 1.0 cm high). Three cages per plant were attached to the first to third fully expanded leaf counting from the top. Five days after inoculation, the number of living and dead whiteflies was recorded and living whiteflies were removed. The number of eggs was counted, and the Oviposition rate (OR) and Adult survival (AS) were calculated according to Bas et al. . In these calculations mortality is assumed constant over time . For the analysis of AS in the F2BC3 population, a Kruskal-Wallis analysis of variance was used . A square root transformation was applied to oviposition rate (OR) prior to the data analysis and analysed by one-way ANOVA followed by a least significant difference (LSD) test.
DNA isolation and genotyping
Genomic DNA was extracted from young leaflets using the micro-prep DNA extraction protocol . The DNA concentration was adjusted to 50 ng/ul. For molecular marker analysis, three types of marker assays were used: CAPs, a custom made Infinium bead array and KASPar (KBiosciences Competitive Alelle-Specific PCR).
For CAPs the PCR reactions were carried out in a final volume of 20 μl, containing 50 ng of genomic DNA, 0.04 μl of DreamTaq polymerase (Fermentas), 2 μl 10X DreamTaq buffer (Fermentas), 0.4 μl of dNTP (5 mM) and 1 μM of each primer (20 pmol). The cycling profile was: 94°C for 3 min, followed by 30 cycles at 94°C for 30 s, 55°C for 30 s, and 72°C for 1 min, and a final extension step at 72°C for 10 min. Aliquots (5 μl) of the amplified products were digested for at least one hour at 37°C in a final volume of 15 μl with 0.5 μl of the appropriate restriction enzyme, using the buffer recommended by the supplier. Amplification and digestion products were analysed by agarose gel electrophoresis (1.5% TBE, agarose) and visualized by GelRed® staining. In Table 2 the primer sequences and the restriction enzymes used are shown.
For genome wide SNP marker analysis, an Infinium bead array was used . On this array, 5528 tomato SNPs were present. Marker analysis was carried out by Service XS Leiden, The Netherlands, according to the Illumina® Infinium HD Ultra Assay protocol (www.illumina.com). After removing missing data and monomorphic markers, 1166 SNP markers were used in the analysis. For fine mapping of the target regions, we developed KASPar assays based on SNP markers that were on the array. The chromosomal positions are according to International Tomato Annotation Group (ITAG) Release 2, official annotations on the SL2.31 version of the tomato genome  (www.solgenomics.net). The sequences flanking the SNPs can be found on http://www.plantbreeding.wur.nl/Publicaions/SNP/4072SNP-Sequences.xlsx . The KASPar assays were run by the van Haeringen lab (VHL), Wageningen, the Netherlands.
Trichomes present on the abaxial side of the leaf were classified according to type . For an estimation of trichome density, the abaxial part of three leaflets was observed under a binocular microscope and a visual scale was used to describe it. The scale used was adapted from Simmons and Gurr  and consisted of four categories: 3, Abundant (>5 per mm2); 2, sparse (5–1 per mm2); 1, very sparse (<1 per mm2), and 0, absent.
The F2BC4S2 plant were analysed for the presence and concentration of methylketones. From each plant one complete leaf (second fully expanded leaf from the top of the plant) was cut, placed immediately into an aluminium envelope and frozen in liquid nitrogen. Each sample was ground to a freeze-dried powder and stored at −80°C until processing. Tree biological replicas were used for the analysis. Each replica consisted on the mix of 5 plants per recombinant class. Per replica, 400 mg of leaf powder was put into a reaction tube with 3 ml of anhydrous dichloromethane (>99.8%, Sigma-Aldrich) as solvent and 0.75 μg per ml heptadecanoic acid methyl ester was added as internal standard. The samples were homogenized using a vortex and centrifuged at 1500 rpm for 10 min. The supernatant was filtered through a soft glass column (Pasteur capillary pipette), which contained 1 cm of silanized glass wool fibres and 2 cm sodium sulphate (Na2SO4) powder as filter. Samples were injected using a 7683 series B injector (Agilent Technologies) into a 7890 A GC (Agilent Technologies) coupled to a 5975 C MSD (Agilent Technologies). Chromatography was performed using a Zb-5MS column (Phenomenex, 30 m, 0.25 mm inner diameter, and 0.25 μm film thickness) with 5 m retention gap. Injection temperature was 250°C, and temperature of column was programmed at 45°C for 1 min, increased by 10°C min−1 to 300°C, and kept at 300°C for 7 min. Column flow was set at 1 ml min−1, using Helium as carrier. The column effluent was ionised and mass spectra was obtained from 35–400 m/z. MetAlign metabolomics software package (www.metalign.nl) was used to perform peak alignment and noise reduction, and MSClust software package (www.metalign.nl) was used for data reduction by clustering several peaks into putative metabolites. Putative metabolites were identified corresponding the obtained mass spectra to the NIST library (National Institute of Standards and Technology, Gaithersburgh, MD, USA), the Wiley online library, and the Wageningen Natural compounds spectral library. Prior to statistical analysis, the metabolites were Log transformed and auto scaled to the mean. To select metabolite compounds putatively related to whitefly preference a t-test, followed by False Discovery Rate correction .
All statistical procedures were performed using the statistical software package GenStat 16th edition. A T-test followed by a False discovery Rate  was done per marker to define the region associated to the reduced oviposition rate.
Plant material development started from F2 plants containing Tv-1, Tv-2 or both using the markers shown in Table 2. Three successive marker assisted backcrosses were carried out with selection for the presence of at least one of the markers linked to the QTL (Figure 1). Twenty-three F2BC3 plants were randomly selected for phenotyping and genotyping to confirm the presence of Tv-1 and Tv-2. As reference lines, we included S. habrochaites (CGN1.1561) and S. lycopersicum cv. Moneymaker. Accession CGN1.1561 showed low values for adult survival (AS, 0.1 ± 0.21 females/day) and oviposition rate (OR, 0.2 ± 0.30 eggs/female/day), and cv. Moneymaker showed high values for adult survival (AS, 1.0 ± 0.01 females/day) and oviposition rate (OR, 5.5 ± 0.72 eggs/female/day). Among the twenty-three F2BC3, variation was found for both parameters. For AS, only three F2BC3 plants were significantly different from cv. Moneymaker (Figure 2). Whereas, for OR a gradient was observed, with fourteen F2BC3 plants showing statistically significant lower values than cv. Moneymaker (Figure 2). To determine the position and size of the introgressions, the twenty-three F2BC3 plants were genotyped using an Infinium bead array . Several plants had an introgression of Tv-1, Tv-2 or parts thereof. None of the plants had the complete Tv-1 and Tv-2 region as defined by Maliepaard et al.  (Figure 2). The four F2BC3 plants with the lowest OR (PV101092-2, PV101088-2, PV101087-3 and PV101088-5) shared an introgression on Chr. 5, but had differences in the presence of the regions Tv-1 and/or Tv-2. One plant (PV101088-8) had the same introgression on Chr. 5, but the OR was not significantly different from cv. Moneymaker (Figure 2).
To further investigate the effect of the introgression on Chr. 5, five F2BC3 plants (PV101092-2, PV101088-2, PV101087-3, PV101093-1, PV101087-2) were selected based on OR, the presence/absence of Tv-1, Tv-2 and the presence of the introgression on Chr. 5. The plants PV101092-2, PV101088-2, PV101087-3, contained the Chr. 5 (61.27 Mbp) introgression, whereas it was smaller in PV101087-2 and not present in PV101093-1. The five plants have varying parts of Tv-1 and Tv-2 or lack these completely (PV101088-2, Figure 2). The plants were backcrossed with cv. Moneymaker to generate five F2BC4 families. All F2BC4 plants plus parental plants, CGN1.1561 and cv. Moneymaker were genotyped for Tv-1, Tv-2 and the introgression on Chr. 5, and phenotyped for adult survival, oviposition rate, trichome type and trichome density. Figure 3 shows the distribution for AS and OR and the link to the respective F2BC3 line. Clear differences were seen between cv. Moneymaker and CGN1.1561 for AS (P < 0.01) and OR (P < 0.01). In the studied F2BC4 plants, there was mainly segregation for OR with the parents on the extremes of the distribution. Genotyping showed that from the offspring of PV101088-2 (renamed to PV101392) four of the five sibling plants were heterozygous for the region on Chr. 5. These plants had an OR level comparable to CGN1.1561, the remaining plant of the five (PV101392-2), lacked the CGN1.1561 allele and had a high OR (Figure 4). To investigate a possible relation between the reduction in OR and the presence of glandular trichomes, the presence/density was determined on the parental lines and the F2BC4 plants. Accession CGN1.1561 was the only one with trichomes type IV and VIc, whereas the F2BC4 plants and the cv. Moneymaker had mainly trichomes type V and VIa. No differences were seen in the density of trichome type VIa among the F2BC4 plants and cv. Moneymaker.
As the offspring of PV101392 showed a low OR and lacked the Tv-1 and Tv-2 region (Figure 4), we focussed on the introgression on Chr. 5. To find offspring plants with a smaller introgression on Chr. 5, PV101392-1PV101392-3, PV101392-4 and PV101392-5 were selfed. Of the 275 F2BC4S1 offspring plants, 33 plants out of 61 recombinants were selected based on length differences of the introgressed region, as judged from marker analysis. The genotyping results (grouped by introgression length) and phenotyping results (OR) are shown in Figure 5. With the F2BC4S1 we could narrow down this introgression to a 3.06 Mbp region between the markers rs2009 (4.76 Mbp) and rs2071 (7.83 Mbp).
To further fine map and confirm the effect of the introgression on Chr. 5, eight F2BC4S1 plants (PV121430-4, PV121430-11, PV121433-30, PV121430-89, PV121432-26, PV121433-29, PV121433-53 and PV121434-57) with a low OR and heterozygous for parts of this region in Chr. 5 were selfed. Of the 295 F2BC4S2, 77 plants out of 154 recombinants were phenotyped based on length differences of the introgressed region, as judged from marker analysis. The results grouped by introgression length are shown in Figure 6. The F2BC4S2 plants with a S. habrochaites (CGN 1.1561) introgression on Chr. 5, between the markers rs2009 (4.76 Mbp) and rs2093 (11.8 Mbp), had an OR similar to the low levels of CGN1.1561 and in the case of plants with the cv. Moneymaker allele homozygous present, the OR was higher (Figure 6). Some of the F2BC4S2 had a reduced adult survival; however, AS and OR were not strongly correlated (R = 0.43) having plants with the introgression on Chr. 5 and with AS levels comparable to those found on cv. Moneymaker and with a significant lower OR. In the F2BC4S2, no plants were found with a further smaller introgression.
Because CGN1.1561 is member of the group of S. habrochaites accessions that accumulate methyl ketones, the F2BC4S2 families were analysed for the presence of those compounds (2-Tridecanone, 2-Undecanone, 2-Pentadecanone and 2-Dodecanone). The accession CGN1.1561 had all these methyl ketones in higher relative abundance compared to cv. Moneymaker and to the F2BC4S2 families (Additional file 1: Table S1). In addition, there were no differences on the relative abundance of these methyl ketones among the F2BC4S2 families and cv. Moneymaker.
An introgression on Chromosome 5 (OR-5)reduces whitefly oviposition rate
Using F2BC3 plants, we identified a S. habrochaites (CGN1.1561) introgression on the short arm of Chr. 5 (hereafter called OR-5), which confers a reduction in B. tabaci oviposition rate. By analysing F2BC4, F2BC4S1 and F2BC4S2 populations, we could confirm that this introgression of 3.06 Mbp is causing the reduced whitefly oviposition rate. The reduction in oviposition caused by the presence of OR-5 is independent of adult survival and the presence of trichome type IV. Plants were found on which all whiteflies were alive but a reduction in oviposition was observed (Figure 3) and none of the plants had the sticky trichomes type IV. The plants homozygous for the S. habrochaites (CGN1.1561) allele in the F2BC4S2 had a higher OR compared to plants heterozygous for this allele. This effect of over dominance might indicate an interaction between the S. habrochaites (CGN1.1561) and the S. lycopersicum allele. It would also implicate that the high level of resistance in terms of low AS and OR found in CGN1.1561 is the result of epistatic interaction between different genes. In the F2BC3 population only one plant had the OR-5 region but with OR levels similar to cv. Moneymaker. This result may be explained in several ways. First, there is the chance of a double recombination in the OR-5 region. However, no double recombination event was detected in this plant with the Infinium bead array. Secondly, there is the possibility of an epistatic effect between OR-5 and a locus different from Tv-1 or Tv-2. Finally, there is always the possibility of a phenotyping artefact.
Selection of the chromosome 5 region
For the selection of the F2BC3 plants, we used markers that are linked to the loci Tv-1 and Tv-2 loci, which are located on Chr. 1 and 12 respectively. It is therefore remarkable that we ended up with an introgression on Chr. 5, which had never actively been selected for. This may be explained by starting with F2 plants containing the introgression on Chr. 5 either homozygous or heterozygous (3 out of 4 plants have the introgression). The chance that plants in the F2BC3 still possess the introgression is 1 out of 4 or 8, which is more or less the number of Chromosome 5 containing F2BC3 plants that we found. The fact that Maliepaard et al.  did not detect the QTL for OR could be caused by the different whitefly species used (T. vaporariorum vs. B. tabaci). Different insect species or biotypes may react differently to the same host plant or odour blend, resulting in different behaviour. For example, glucosinolates can confer resistance to some insects, whereas they can be used as host and strong oviposition cues for others . In the case of whiteflies, differences were seen when compared the feeding behaviour of the Q and the B-biotype on the same host plant . Also, tomato plants carrying the Mi1-2 gene were in general more resistance to the Q-biotype than to the B-biotype .
Nature of the resistance provided by OR-5
Several QTLs related to whitefly resistance have been identified on Chr. 5 (Table 1). Maliepaard et al.  found in the region of OR-5 a QTL (TriIV-1) that increases the density of trichomes type IV. However, we did not detect any type IV trichomes on the F2BC4 plants containing the OR-5 introgression. In a backcross population of potato ((S. tuberosum × S. berthaultii) × S. berthaultii) a region on Chr. 5 was associated with a reduction in the oviposition rate and leaf consumption by the Colorado potato beetle (Leptinotarsa decemlineata) . This region also had a large effect on the density of the glandular secretory type B trichome (LOD: 19.17, explaining 35.6% of the variance), furthermore differences in the sucrose ester levels and in the presence of droplet (exudate) on the tip of the trichomes were associated with this region on Chr. 5 . For S. pennellii, two QTLs were described on Chr. 5 that are involved in acyl sugar metabolism, one (TA5) related to the total accumulation of acyl sugars and another (5) related to the proportion of 7-methyloctanonate and 9-methyldecanoncate fatty acids that are incorporated into acyl sugars ,. To check if acyl sugars were related to the reduction in oviposition rate an LC-MS chromatography analysis was done on the F2BC4S2 plants . No differences were found among the F2BC4S2 plants on the levels of acyls sugars, pointing to a different mechanism of resistance in this plant material specifically affecting whitefly oviposition rate (data not shown). As the parental accession CGN1.1561 accumulates methyl ketones, we also analysed the offspring for the presence of 2-Tridecanone, 2-Undecanone, 2-Pentadecanone and 2-Dodecanone. None of these compounds was detected at elevated levels in the offspring, excluding the option that these methyl ketones may explain the observed reduction in oviposition rate. On the 3.06 Mbp introgression of OR-5 are 258 annotated genes including R-genes, transcription factors, genes involved in acyl sugar and terpenoid metabolism which can be considered as candidate genes for reduced oviposition. To reduce the list of candidate genes and find the gene(s) responsible for the lower OR further fine mapping and functional analysis, including more detailed metabolomics is needed. However, considering the lack of recombinants found in the F2BC4S1 and F2BC4S2 populations between the markers rs-7543 (5.79 Mbp) and rs-7551 (6.15 Mbp), it might be difficult to reduce the size of the introgression.
Perspectives of OR-5 for breeding whitefly resistant varieties
Since the late nineties of the 20th century, the efforts to get whitefly resistant tomatoes have increased considerably, but so far they have been unsuccessful . The screening of genetic resources for novel whitefly resistance mechanisms has increased, going from distant wild relatives of tomato (i.e. S. pennellii, S. habrochaites) to in depth studies of several accessions of closely related species (i.e. S. galapagense, S. pimpinellifolium) -,,,. These efforts have led to the identification of specific secondary metabolites conferring resistance to whiteflies (methyl ketones, sesquiterpenes, and acyl sugars) ,,, the identification of QTLs related to resistance ,,, and in some cases to the genes involved in the synthesis of resistance related metabolites ,,,,. The identification of OR-5, affecting specifically whitefly oviposition rate and independent of the presence of trichome type IV, opens new opportunities for breeding. The OR-5 region is expected to reduce population development of B. tabaci strongly. As the reduction in oviposition is not linked to the sticky trichomes type I and IV, and the known negative effect of this type of resistance on parasitoids and predators ,-, it can be expected that this resistance will be very suitable in combination with biological control. On varieties containing the OR-5 region the B. tabaci population development will be slowed down giving the natural enemies ample opportunity to keep the population below threshold levels or even to remove developing whiteflies. Therefore, the gene will in particular be useful in protected tomato production conditions (greenhouse cultivation). For open field production, the resistances based on trichomes type I and IV will be more suitable .
We identified a region at the top of Chr. 5 (OR-5), which is associated with a reduction in the oviposition rate of B. tabaci. This reduction was independent of the presence of the QTLs Tv-1 and Tv-2 that were identified previously , as well as of the presence of trichomes type IV. The OR-5 locus will provide new opportunities for resistance breeding against whiteflies, which is especially relevant in greenhouse cultivation.
Martin JH, Mound LA: An annotated check list of the world’s whiteflies (Insecta: Hemiptera: Aleyrodidae). Zootaxa. 2007, 1492: 1-84.
Brown JK, Bird J: Whitefly-transmitted geminiviruses and associated disorders in the Americas and the Caribbean Basin. Plant Dis. 1992, 76 (3): 220-225. 10.1094/PD-76-0220.
Byrne DN, Bellows TS: Whitefly biology. Annu Rev Entomol. 1991, 36: 431-457. 10.1146/annurev.en.36.010191.002243.
Broekgaarden C, Snoeren TAL, Dicke M, Vosman B: Exploiting natural variation to identify insect-resistance genes. Plant Biotech J. 2011, 9 (8): 819-825. 10.1111/j.1467-7652.2011.00635.x.
Schoonhoven LM, van Loon JJA, Dicke M: Insect-plant Biology. 2005, Oxford University Press, Oxford
Heinz KM, Zalom FG: Variation in trichome-based resistance to Bemisia argentifolii (Homoptera: Aleyrodidae) oviposition on tomato. J Econo Entomol. 1995, 88 (5): 1494-1502. 10.1093/jee/88.5.1494.
Muigai SG, Schuster DJ, Snyder JC, Scott JW, Bassett MJ, McAuslane HJ: Mechanisms of resistance in Lycopersicon germplasm to the whitefly Bemisia argentifolii . Phytoparasitica. 2002, 30 (4): 347-360. 10.1007/BF02979682.
Freitas J, Maluf W, das Graças Cardoso M, Gomes LA, Bearzotti E: Inheritance of foliar zingiberene contents and their relationship to trichome densities and whitefly resistance in tomatoes. Euphytica. 2002, 127 (2): 275-287. 10.1023/A:1020239512598.
Muigai SG, Bassett MJ, Schuster DJ, Scott JW: Greenhouse and field screening of wild Lycopersicon germplasm for resistance to the whitefly Bemisia argentifolii . Phytoparasitica. 2003, 31 (1): 27-38. 10.1007/BF02979764.
Fernández-Muñoz R, Salinas M, Álvarez M, Cuartero J: Inheritance of resistance to two-spotted spider mite and glandular leaf trichomes in wild tomato Lycopersicon pimpinellifolium (Jusl.) Mill. J Am Soc Hortic Sci. 2003, 128 (2): 188-195.
Baldin ELL, Vendramim JD, Lourenção AL: Resistance of tomato genotypes to the whitefly Bemisia tabaci (Gennadius) biotype B (Hemiptera: Aleyrodidae). Neotrop Entomol. 2005, 34 (3): 435-441. 10.1590/S1519-566X2005000300012.
Rodriguez-Lopez MJ, Garzo E, Bonani JP, Fereres A, Fernandez-Munoz R, Moriones E: Whitefly resistance traits derived from the wild tomato Solanum pimpinellifolium affect the preference and feeding behavior of Bemisia tabaci and reduce the spread of Tomato yellow leaf curl virus. Phytopathology. 2011, 101 (10): 1191-1201. 10.1094/PHYTO-01-11-0028.
Firdaus S, Heusden A, Hidayati N, Supena E, Visser R, Vosman B: Resistance to Bemisia tabaci in tomato wild relatives. Euphytica. 2012, 187 (1): 31-45. 10.1007/s10681-012-0704-2.
Lucatti AF, Van Heusden AW, De Vos RCH, Visser RGF, Vosman B: Differences in insect resistance between tomato species endemic to the Galapagos Islands. BMC Evol Biol. 2013, 13: 175-10.1186/1471-2148-13-175.
Rodríguez-López MJ, Garzo E, Bonani JP, Fernández-Muñoz R, Moriones E, Fereres A: Acylsucrose-producing tomato plants forces Bemisia tabaci to shift its preferred settling and feeding site. PLoS ONE. 2012, 7 (3): e33064-10.1371/journal.pone.0033064.
Firdaus S, Heusden AW, Hidayati N, Supena E, Mumm R, Vos RCH, Visser RGF, Vosman B: Identification and QTL mapping of whitefly resistance components in Solanum galapagense . Theor Appl Genet. 2013, 126 (6): 1487-1501. 10.1007/s00122-013-2067-z.
Bleeker PM, Mirabella R, Diergaarde PJ, VanDoorn A, Tissier A, Kant MR, Prins M, De Vos M, Haring MA, Schuurink RC: Improved herbivore resistance in cultivated tomato with the sesquiterpene biosynthetic pathway from a wild relative. Proc Natl Acad Sci U S A. 2012, 109 (49): 20124-20129. 10.1073/pnas.1208756109.
Antonius G: Production and quantification of methyl ketones in wild tomato accessions. J Environ Sci Health. 2001, 36: 835-848. 10.1081/PFC-100107416.
Ben-Israel I, Yu G, Austin MB, Bhuiyan N, Auldridge M, Nguyen T, Schauvinhold I, Noel JP, Pichersky E, Fridman E: Multiple biochemical and morphological factors underlie the production of methylketones in tomato trichomes. Plant Physiol. 2009, 151: 1952-1964. 10.1104/pp.109.146415.
Fridman E, Wang J, Iijima Y, Froehlich JE, Gang DR, Ohlrogge J, Pichersky E: Metabolic, genomic, and biochemical analyses of glandular trichomes from the wild tomato species Lycopersicon hirsutum identify a key enzyme in the biosynthesis of methylketones. Plant Cell. 2005, 17: 1252-1267. 10.1105/tpc.104.029736.
Leckie BM, Jong DM, Mutschler MA: Quantitative trait loci increasing acylsugars in tomato breeding lines and their impacts on silverleaf whiteflies. Mol Breeding. 2012, 30: 1621-1634. 10.1007/s11032-012-9746-3.
Milligan SB, Bodeau J, Yaghoobi J, Kaloshian I, Zabel P, Williamson VM: The root-knot nematode resistance gene Mi from tomato is a member of leucine zipper, nucleotide binding, leucine-rich repeat family of plant genes. Plant Cell. 1998, 10: 1307-1319. 10.1105/tpc.10.8.1307.
Kaloshian I, Lange WH, Williamson VM: An aphid-resistance locus is tightly linked to the nematode-resistance gene, Mi, in tomato. Proc Natl Acad Sci U S A. 1995, 92: 622-625. 10.1073/pnas.92.2.622.
Rossi M, Goggin FL, Milligan SB, Kaloshian I, Ullman DE, Williamson VM: The nematode resistance gene Mi of tomato confers resistance against the potato aphid. Proc Natl Acad Sci U S A. 1998, 95: 9750-9754. 10.1073/pnas.95.17.9750.
Nombela G, Williamson VM, Muñiz M: The root-knot nematode resistance gene Mi-1.2 of tomato is responsible for resistance against the whitefly Bemisia tabaci . MPMI. 2003, 16: 645-649. 10.1094/MPMI.2003.16.7.645.
Casteel CL, Walling LL, Paine TD: Behavior and biology of the tomato psyllid, Bactericerca cockerelli, in response to the Mi-1.2 gene. Entomol Exper Appl. 2006, 121: 67-72. 10.1111/j.1570-8703.2006.00458.x.
Nombela G, Beitia F, Muñiz M: Variation in tomato host response to Bemisia tabaci (Hemiptera: Aleyrodidae) in relation to acyl sugar content and presence of the nematode and potato aphid resistance gene Mi . Bull Entomol Res. 2000, 90: 161-167. 10.1017/S0007485300000274.
Maliepaard C, Bas NJ, Van Heusden S, Kos J, Pet G, Verkerk R, Vrielink R, Zabel P, Lindhout P: Mapping of QTLs for glandular trichome densities and Trialeurodes vaporariorum (greenhouse whitefly) resistance in an F2 from Lycopersicon esculentum X Lycopersicon hirsutum f. glabratum . Heredity. 1995, 75 (4): 425-433. 10.1038/hdy.1995.155.
Schilmiller A, Shi F, Kim J, Charbonneau AL, Holmes D, Daniel Jones A, Last RL: Mass spectrometry screening reveals widespread diversity in trichome specialized metabolites of tomato chromosomal substitution lines. Plant J. 2010, 62 (3): 391-403. 10.1111/j.1365-313X.2010.04154.x.
Schilmiller AL, Charbonneau AL, Last RL: Identification of a BAHD acetyltransferase that produces protective acyl sugars in tomato trichomes. Proc Natl Acad Sci U S A. 2012, 109 (40): 16377-16382. 10.1073/pnas.1207906109.
Schilmiller AL, Schauvinhold I, Larson M, Xu R, Charbonneau AL, Schmidt A, Wilkerson C, Last RL, Pichersky E: Monoterpenes in the glandular trichomes of tomato are synthesized from a neryl diphosphate precursor rather than geranyl diphosphate. Proc Natl Acad Sci U S A. 2009, 106 (26): 10865-10870. 10.1073/pnas.0904113106.
Van der Hoeven RS, Monforte AJ, Breeden D, Tanksley SD, Steffens JC: Genetic control and evolution of sesquiterpene biosynthesis in Lycopersicon esculentum and L. hirsutum . Plant Cell. 2000, 12 (11): 2283-2294. 10.1105/tpc.12.11.2283.
Momotaz A, Scott JW, Schuster DJ: Identification of quantitative trait loci conferring resistance to Bemisia tabaci in an F2 population of Solanum lycopersicum × Solanum habrochaites accession LA1777. J Am Soc Hortic Sci. 2010, 135 (2): 134-142.
Van den Elsen F: Resistance Mechanisms against Bemisia tabaci in Wild Relatives of Tomato. 2013, Thesis Wageningen University, Wageningen, The Netherlands
Blauth SL, Churchill GA, Mutschler MA: Identification of quantitative trait loci associated with acylsugar accumulation using intraspecific populations of the wild tomato, Lycopersicon pennellii. Theor Appl Genet. 1998, 96 (3–4): 458-467. 10.1007/s001220050762.
Bas N, Mollema C, Lindhout P: Resistance in Lycopersicon hirsutum f. glabratum to the greenhouse whitefly (Trialeurodes vaporariorum) increases with plant age. Euphytica. 1992, 64 (3): 189-195. 10.1007/BF00046048.
Van Giessen WA, Mollema C, Elsey KD: Design and use of a simulation model to evaluate germoplasm for antibiotic resistance to the greenhouse whitefly (Trialeurodes vaporariorum) and the sweetpotato whitefly (Bemisia tabaci). Entomol Exp Appl. 1995, 76: 271-286. 10.1111/j.1570-7458.1995.tb01971.x.
Zar JH: Biostatistical Analysis. 2010, Pearson Education International, Upper Saddle River, NJ
Fulton T, Chunwongse J, Tanksley S: Microprep protocol for extraction of DNA from tomato and other herbaceous plants. Plant Mol Biol Rep. 1995, 13 (3): 207-209. 10.1007/BF02670897.
Viquez-Zamora M, Vosman B, van de Geest H, Bovy A, Visser R, Finkers R, van Heusden A: Tomato breeding in the genomics era: insights from a SNP array. BMC Genomics. 2013, 14 (1): 354-10.1186/1471-2164-14-354.
Consortium TG: The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012, 485 (7400): 635-641. 10.1038/nature11119.
Channarayappa C, Shivashankar G, Muniyappa V, Frist RH: Resistance of Lycopersicon species to Bemisia tabaci, a tomato leaf curl virus vector. Can J Bot. 1992, 70 (11): 2184-2192. 10.1139/b92-270.
Simmons AT, Gurr GM: Trichomes of Lycopersicon species and their hybrids: effects on pests and natural enemies. Agric For Entomol. 2005, 7 (4): 265-276. 10.1111/j.1461-9555.2005.00271.x.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc. 1995, 57 (1): 289-300.
Hopkins RJ, van Dam NM, van Loon JJA: Role of glucosinolates in insect-plant relationships and multitrophic interactions. Annu Rev Entomol. 2009, 54 (1): 57-83. 10.1146/annurev.ento.54.110807.090623.
Jiang YX, Lei H, Collar JL, Martin B, Muniz M, Fereres A: Probing and feeding behavior of two distinct biotypes of Bemisia tabaci (Homoptera: Aleyrodidae) on tomato plants. J Econo Entomol. 1999, 92 (2): 357-366. 10.1093/jee/92.2.357.
Nombela G, Beitia F, Muñiz M: A differential interaction study of Bemisia tabaci Q-biotype on commercial tomato varieties bearing or not bearing the Mi resistance gene, and comparative host responses with the B-biotype. Entomol Exp Appl. 2001, 98 (3): 339-344. 10.1046/j.1570-7458.2001.00790.x.
Bonierbale MW, Plaisted RL, Pineda O, Tanksley SD: QTL analysis of trichome-mediated insect resistance in potato. Theor Appl Genet. 1994, 87 (8): 973-987. 10.1007/BF00225792.
Mutschler MA, Doerge RW, Liu SC, Kuai JP, Liedl BE, Shapiro JA: QTL analysis of pest resistance in the wild tomato Lycopersicon pennellii: QTLs controlling acylsugar level and composition. Theor Appl Genet. 1996, 92 (6): 709-718. 10.1007/BF00226093.
Blauth SL, Steffens JC, Churchill GA, Mutschler MA: Identification of QTLs controlling acylsugar fatty acid composition in an intraspecific population of Lycopersicon pennellii (Corr.) D’Arcy. Theor Appl Genet. 1999, 99 (1–2): 373-381. 10.1007/s001220051247.
Bleeker PM, Diergaarde PJ, Ament K, Schütz S, Johne B, Dijkink J, Hiemstra H, De Gelder R, De Both MTJ, Sabelis MW, Haring MA, Schuurink RC: Tomato-produced 7-epizingiberene and R-curcumene act as repellents to whiteflies. Phytochemistry. 2011, 72 (1): 68-73. 10.1016/j.phytochem.2010.10.014.
Farrar RR, Kennedy GG, Kashyap RK: Influence of life history differences of two tachinid parasitoids of Helicoverpa zea (Boddie) (Lepidoptera: Noctuidae) on their interactions with glandular trichome/methyl ketone-based insect resistance in tomato. J Chem Ecol. 1992, 18: 499-515. 10.1007/BF00994248.
Kashyap RK, Kennedy GG, Farrar RR: Mortality and inhibition of Helicoverpa zea Egg parasitism rates by Trichogramma in relation to trichome/methyl ketone-mediated insect resistance of Lycopersicon hirsutum f. glabratum, accession PI 134417. J Chem Ecol. 1991, 17: 2381-2395. 10.1007/BF00994589.
Kennedy GG: Tomato, pests, parasitoids, and predators: Tritrophic interactions involving the genus Lycopersicon . Ann Rev Entomol. 2003, 48: 51-72. 10.1146/annurev.ento.48.091801.112733.
We would like to thank to Jarno Sinnige and Greet Steenhuis for their help during the phenotyping experiments and to Colette Broekgaarden and the two anonymous reviewers for their comments on an earlier version of this paper. AFL was financially supported by the Foundation CAPACIT-AR del NOA and by Wageningen UR Plant Breeding.
The authors declare that they have no competing interests.
The project was conceived by AFL, BV and SvH. AFL performed the experiments and data analysis. AFL and FRGMD performed the crossings and population developments. AFL and RM performed the GC-MS experiment and data analysis. AFL, BV, RM, RGFV and SvH have been involved in the writing of the manuscript. All authors read and approved the final manuscript.