Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

Differences in insect resistance between tomato species endemic to the Galapagos Islands

Alejandro F Lucatti12, Adriaan W van Heusden1, Ric CH de Vos345, Richard GF Visser1 and Ben Vosman1*

Author affiliations

1 Wageningen UR Plant Breeding, Wageningen University and Research, Centre, P.O. Box 386, Wageningen, AJ 6700, The Netherlands

2 Graduate School Experimental Plant Sciences. Wageningen Campus. Droevendaalsesteeg 1, 6708 PB, Wageningen, The Netherlands

3 Plant Research International, Business Unit Bioscience, Wageningen University and Research Centre, P.O. Box 619, Wageningen, AP 6700, The Netherlands

4 Centre for BioSystems and Genomics, P.O. Box 98, Wageningen, AB 6700, The Netherlands

5 Netherlands Metabolomics Centre, Einsteinweg 55, 2333 CC, Leiden, The Netherlands

For all author emails, please log on.

Citation and License

BMC Evolutionary Biology 2013, 13:175  doi:10.1186/1471-2148-13-175


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/13/175


Received:4 April 2013
Accepted:21 August 2013
Published:24 August 2013

© 2013 Lucatti et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The Galapagos Islands constitute a highly diverse ecosystem and a unique source of variation in the form of endemic species. There are two endemic tomato species, Solanum galapagense and S. cheesmaniae and two introduced tomato species, S. pimpinellifolium and S. lycopersicum. Morphologically the two endemic tomato species of the Galapagos Islands are clearly distinct, but molecular marker analysis showed no clear separation. Tomatoes on the Galapagos are affected by both native and exotic herbivores. Bemisia tabaci is an important introduced insect species that feeds on a wide range of plants. In this article, we address the question whether the differentiation between S. galapagense and S. cheesmaniae may be related to differences in susceptibility towards phloem-feeders and used B. tabaci as a model to evaluate this.

Results

We have characterized 12 accessions of S. galapagense, 22 of S. cheesmaniae, and one of S. lycopersicum as reference for whitefly resistance using no-choice experiments. Whitefly resistance was found in S. galapagense only and was associated with the presence of relatively high levels of acyl sugars and the presence of glandular trichomes of type I and IV. Genetic fingerprinting using 3316 SNP markers did not show a clear differentiation between the two endemic species. Acyl sugar accumulation as well as the climatic and geographical conditions at the collection sites of the accessions did not follow the morphological species boundaries.

Conclusion

Our results suggest that S. galapagense and S. cheesmaniae might be morphotypes rather than two species and that their co-existence is likely the result of selective pressure.

Keywords:
Bemisia tabaci; Solanum galapagense; Solanum cheesmaniae; Whitefly; Trichomes; Acyl sugars; Selection pressure

Background

Tomatoes are native to South America and can be found from the north of Chile/Argentina to Ecuador, including the Galapagos Islands. The Galapagos is a volcanic archipelago of 13 islands located about 1000 km from the coast of Ecuador. On this archipelago, vegetation varies among islands, altitude and cardinal direction [1,2]. The Galapagos Islands constitute a highly diverse ecosystem and a unique source of variation in the form of endemic species. There are two endemic tomato species, Solanum galapagense and S. cheesmaniae and two introduced species Solanum pimpinellifolium and S. lycopersicum. It is believed that the latter two have been introduced on the islands during the twentieth century, though it is possible that S. pimpinellifolium is present on the islands a lot longer [1]. The endemic Galapagos Island tomato species have evolved in isolation from the mainland species, resulting in clearly differing morphological features compared to the species that were introduced later. However, natural hybrids have been found [1,3]. The taxonomic status of the Galapagos’ endemic tomatoes is under debate and a historic overview is given by Darwin et al. (2003). Darwin and co-workers adopted the ‘morphological cluster’ species concept [4] to divide the two endemic tomato forms in two species S. galapagense and S. cheesmaniae. Despite the clear separation obtained on the basis of morphology, it was not possible to separate the two species with molecular markers [2,5-7].

From the total number of endemic species on the Galapagos Islands, 47% are insects [8]. However, the number of exotic insect species on the Galapagos is increasing due to human activity and now at least 463 exotic insect species can be found on the archipelago, of which 193 species are herbivores and in majority phloem feeders [9]. Of the exotic insect species 73% are naturalized or are known to feed on endemic plant species. The whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is one of the most important invasive insects on the Galapagos Islands, receiving the highest score of invasiveness due to their wide distribution, wide host range and their importance as vector of many plant viruses [9-11]. Recently it became clear that B. tabaci is not a single species but a whole complex of at least 36 cryptic species [12-15]. The first report of B. tabaci on the Galapagos dates back to 1998 [9]. Interestingly, S. galapagense was shown to be resistant to this devastating pest insect [14,16,17]. This suggests that a general mechanism conferring resistance towards insects is present in S. galapagense as the time after introduction of B. tabaci is too short for co-evolution. The resistance is associated with the presence of type IV glandular trichomes [14,17].

Selection pressure, co-evolutionary processes, population genetics, bio-geographical variables, and gene dynamics (gene flow, drift, mating systems, etc.) can affect the evolution of resistance genes and the accumulation or the presence of certain metabolites related to resistance [18,19]. For example, the evolution of genes involved in the terpenoid pathway and the relative composition of acyl sugar were suggested to be related to geographic and climatic variables under which populations of Solanum habrochaites are found [20,21]. The Galapagos Islands have proven to be a place of natural experimental conditions to answer questions related to processes like founder effect, genetic drift, divergent selection, ecological opportunity and interaction between endemic and exotic species in the speciation processes [8]. As B. tabaci has an ecological importance on the Galapagos Islands [9] and its global importance as a phloem feeding tomato herbivore, we decided to use B. tabaci as a model insect to analyse factors underlying insect resistance in relation to species delimitation between S. galapagense and S. cheesmaniae.

So far, only a limited number of accessions/populations of S. galapagense and S. cheesmaniae have been evaluated for insect resistance and therefore it is unknown if the insect resistance coincides with the species boundaries (based on the morphological differences). Neither is there any knowledge about the relation between geographical and climatic conditions today on the Galapagos and the occurrence of the two species. Recently it was shown that the whitefly resistance in an accession of S. galapagense is most likely based on the production of acyl sugars in the glandular trichomes [22], it is unknown if the relative acyl sugar concentration among the different accessions of S. galapagense and S. cheesmaniae coincides with species boundaries and insect resistance. In the present study, we address the questions raised and discuss the implication in an evolutionary context. We characterized the genetic and acyl sugar variation in 34 endemic tomato accessions and investigated if resistance and chemical variation among accessions are correlated. We demonstrate that S. galapagense is different form S. cheesmaniae in the resistance towards whiteflies and in the trichome composition. Geographic and climatic variables do not explain the distribution pattern found for the Galapagos’ endemic tomatoes. Genetic variation between the two species is almost absent and acyl sugar composition does not completely follow the morphological species boundaries. All together our results suggest that S. galapagense and S. cheesmaniae might be considered as morphotypes rather than two species and that their co-existence is likely the result of selective pressure.

Results

Resistance to whitefly

The level of whitefly resistance in accessions of S. galapagense, S. cheesmaniae and cv. Moneymaker was assessed using three parameters, namely adult survival (AS), oviposition rate (OR), and pre-adult survival (PS) (Table 1). The three parameters were highly correlated (Table 2). For AS, significant differences were found among accessions (ANOVA, p < 0.001), with survival rates ranging from 0 to 1. The lowest values for AS were found within the accessions of S. galapagense, with five accessions on which all whiteflies were dead after 5 days (AS = 0). Adult survival on all accessions of S. galapagense was statistically different from the AS on cv. Moneymaker. None of the S. cheesmaniae accessions were statistically different from cv. Moneymaker for the resistance variables (Table 1). When comparing AS at the species level, S. galapagense was the most resistant species (AS = 0.09 ± 0.12, p < 0.001).

Table 1. Adult survival, oviposition rate and pre-adult survival of the different accessions of Solanum galapagense and S. cheesmaniae

Table 2. Pearson correlation among phenotypic resistance characteristics and Acyl sugar accumulation

The oviposition rate (OR) ranged from 0 to 8.86 eggs/day/female (Table 1). All S. galapagense accessions showed a reduction in OR and were significantly different from cv. Moneymaker. In fact there were three accessions on which no eggs were deposited. All S. cheesmaniae accessions were at least as susceptible as the cv. Moneymaker. At the species level, S. galapagense showed the lowest OR (OR = 1.3 ± 1.7, p < 0.001).

For several accessions it was impossible to determine PS, as there were no eggs deposited. For two accessions of S. galapagense (LA1452 and LA1401) we observed that, although whiteflies were able to lay a few eggs, none of those eggs hatched. The S. cheesmaniae accessions were not significantly different from the cv. Moneymaker in the PS values (Table 1).

Trichome type and acyl sugar composition

A clear difference in the trichome composition was observed among the different tomato species. All accessions of S. galapagense had trichomes of both type I and IV. None of the accessions of S. cheesmaniae had trichomes type I or IV, while all of them had trichomes type V and VI in higher or lower densities (Additional file 1: Table S1). Trichome type III was absent on S. galapagense and present on 7 of the 22 accessions of S. cheesmaniae (Additional file 1: Table S1). A negative correlation between the measured resistance variables and the presence of glandular trichome types I, IV and VI was observed (Table 2).

Additional file 1: Table S1. This file includes GeneBank origin, accession numbers, collection sites and all phenotypic information, SNP array data and LC-MS data for the material used in this work. The Phylogenetic tree was uploaded to TreeBASe.org, and can be found by following this link: http://purl.org/phylo/treebase/phylows/study/TB2:S14512 webcite.

Format: XLSX Size: 559KB Download fileOpen Data

The LC-MS analysis allowed us to determine the relative abundance of 12 acyl sugars. All of them were acyl sucroses: five with 3 lateral branches (S3-acyl sucroses) and seven with 4 lateral branches (S4-acyl sucroses). A heatmap (hierarchical clustering) visualizing the results is shown in Figure 1 and the relative amounts measured can be found in Additional file 1: Table S1. In Figure 1, we see two clusters: one including 7 accessions of S. cheesmaniae and the cv. Moneymaker and the other as an intermingled cluster with all the accessions of S. galapagense and the remaining accessions of S. cheesmaniae. Total acyl sugar abundance (expressed as the total of acyl sugar peak areas per accession) was higher in plants with a higher resistance level, measured as lower values for AS, OR and PS (R2 = 0.37, p < 0.001). In general, S. galapagense accessions accumulate higher levels of acyl sugars than S. cheesmaniae accessions. Nevertheless, there were four accessions of S. cheesmaniae (LA0421, LA0746, LA0521 and LA0529) with levels of acyl sucroses that were as high as those found in S. galapagense accessions (Figure 2, Additional file 1: Table S1). Regression analysis did not indicate a specific relation between acyl sugar accumulation, neither qualitatively nor quantitatively, and the type of trichomes present on a plant. However, we observed a slight positive correlation between total acyl sugars and the presence of glandular trichomes type I and IV (Table 2).

thumbnailFigure 1. Hierarchical clustering (Pearson correlation, UPGMA), of the different acyl sucroses (columns) and the different accessions (rows).Solanum galapagense S.C. Darwin & Peralta (violet circle symbol), S. cheesmaniae L. Riley (red cross sign symbol), S. lycopersicum L cv. Moneymaker (yellow rhombus symbol). Colour key is displayed in the figure.

thumbnailFigure 2. Acyl sugar accumulation per species and accession. a) Average acyl sugar accumulation per species; values for S. galapagense and S. cheesmaniae represent means and SD of 14 and 21 plants, respectively. Cultivar Moneymaker is included as reference. Different letters indicate statistical differences according to LSD test (p < 0.05); b) PCA-biplot score of the tomato accessions based on acyl sugar accumulation. Solanum galapagense (violet circle symbol, S. cheesmaniae (red plus sign symbol) and S. lycopersicum cv. Moneymaker (yellow rhombus symbol). Lines indicate the loadings of the different acyl sucroses.

Genetic relationships and correlations among accessions

We used a SNP array [23] to determine the genetic relationships between the different accessions. The Neighbour joining tree (Figure 3), based on 3316 markers that showed polymorphisms in or among S. lycopersicum, S. pimpinellifolium, S. cheesmaniae and S. galapagense, indicated a tight cluster with all the accessions of S. galapagense and S. cheesmaniae intermingled. Two accessions of S. cheesmaniae (LA3124 and G1.1615) were clearly separated from the rest and are most likely hybrids with either S. lycopersicum or S. pimpinellifolium. After exclusion of these two deviating S. cheesmaniae accessions only 53 polymorphic markers (Additional file 2: Figure S1) were detected within S. galapagense and S. cheesmaniae. From those 53 SNPs, 44 were polymorphic among S. cheesmaniae accessions and 9 among S. galapagense accessions. Most of the polymorphic markers were randomly distributed over the two species, 7 of them were found only in one accession and 2 were unique to S. galapagense. None of the polymorphic markers were fixed in either species.

thumbnailFigure 3. Neighbour-joining tree of the accessions. The NJ tree was based on Manhattan distances using 3324 polymorphic markers. Per sample the species name is followed by the accession code (S. che. for Solanum cheesmaniae; S. gal. for S. galapagense; S. pim. for S. pimpinellifolium; S. lyc. for S. lycopersicum; S. per. for S. peruvianum; S. hab. for S. habrochaites and S. pen. for S. pennellii). Bootstrap values (higher than 50) are shown on the branches.

Additional file 2: Figure S1. Isolation by distance analysis. This figure describes the relation between the Genetic distances and the Geographic distances.

Format: PPT Size: 157KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

Correlation of species presence to geographic/climatic variables

We analysed the available geographic/climatic data for the Galapagos Islands in relation to the location where the S. galapagense and S. cheesmaniae accessions were collected, but no correlation was detected (Figure 4).

thumbnailFigure 4. PCA-biplot score of the accessions of Solanum galapagense (violet circle symbol) and S. cheesmaniae (red plus sign symbol) based on geographic/climatic variables.

Discussion

Whitefly resistance in S. galapagense, a combination of trichomes and metabolites

A difference in whitefly resistance was observed between S. galapagense (resistant) and S. cheesmaniae (susceptible). While some Solanum species are considered to be resistant, it is important to note that the level of resistance may vary between and within accessions of the same species. Next to resistance to whitefly [17], such inter- and intra-variation was also described, among others, in a screenings for white mold and late blight resistance in potato [17,24-26]. We have found resistance to whitefly in all S. galapagense accessions, which was accompanied by high densities of trichomes type IV and high acyl sugar accumulation. Based on available data and literature it is not known whether the original founder had type IV trichomes or not. Even though the accessions of S. galapagense differ in the relative amounts of the different acyl sugars present, they were all resistant. This was not the case for S. cheesmaniae. Although some accessions of S. cheesmaniae accumulated acyl sugars to levels comparable to those found in S. galapagense they were all susceptible, probably because they lack trichome type IV. Contrarily, we also observed that in some accessions of S. pimpinellifolium, although they had trichome type IV, the levels of acyl sugars and the resistance variables were not different from those found on cv. Moneymaker (data not shown). The role of glandular trichomes and acyl sugars in insect resistance has been discussed frequently [11,16,17,27,28]. In tomato, the synthesis and accumulation of acyl sugars takes place within the glandular head of the trichome [29]. Acyl sugars are non-specific resistance components providing resistance to a broad spectrum of insects of different feeding guilds (whiteflies, aphids, leaf miners, caterpillars, etc.) [30]. Based on our data it is likely that a minimum level of acyl sugars and the presence of glandular trichomes type IV are needed to achieve an effective level of resistance and a fully resistant phenotype. It also suggests that not a specific acyl sugar, but rather the total amount of acyl sugars is important for resistance. Having said so, it is possible that others metabolites, not detected by the LC-MS analysis, may be important as well.

Galapagos’ endemic tomato species cannot be differentiated by genetic analysis

The clear morphological differences between S. galapagense and S. cheesmaniae with regard to their leaf morphology, trichome composition, internode length, among others, were the reason to consider them as two distinct species [1]. However, our SNP array analysis based on 5528 markers [23], of which 3316 were polymorphic among the Lycopersicon group of Solanum sect. Lycopersicon, showed only 53 polymorphisms within the two species, with only two alleles specific to S. galapagense and no fixed alleles. As to be expected also in the NJ tree the two species could not be separated. Similar observations have been made using AFLP markers [5]. The observed low genetic variation among accessions of the two species under study is consistent with the hypothesis that these endemic species are the result of a unique founder event on the Galapagos Islands, followed by morphological divergence. The very obvious morphological differences between S. galapagense and S. cheesmaniae might be due to a limited number of genetic changes that have not been picked up by our SNP array. For instance, a single nucleotide deletion in the promoter of the PTS/TKD1 gene results in a marked change of leaf complexity as seen in S. galapagense[31]. Solanum galapagense and S. cheesmaniae are part of a monophyletic clade within Solanum sect. Lycopersicon, together with S. pimpinellifolium and S. lycopersicum[7]. Thus, based on the molecular marker analysis it might be more appropriate to consider them as morphotypes, rather than as different species. In line with this, hybrids between S. galapagense and S. cheesmaniae have been found [3]. In addition, we cannot ignore the clear morphological differences (trichomes, whitefly resistance, leaf morphology, etc.) found by us and by others [1]. All this suggest that the species status of these two tomatoes should be reconsidered carefully taking into consideration the different concepts used to define the species.

Biogeographic variables, phylogenetic studies and experimental approaches can be used to answer questions like adaptation to local environments, divergence and prediction of a phenotype based on the environment where a species was found [19,25,32-34]. We observed no correlation between the geographic location where the accession was collected and the occurrence of one of the two morphotypes. Some studies have addressed the population structure of wild tomatoes and external factors causing this. Caicedo and Schaal [35] observed a clear structure in S. pimpinellifolium populations collected from northern to southern Peru, which was consistent with a hypothesis based on genetic isolation by distance. It was also observed that genetically closely related accessions can be found far away from each other. In our analysis, we also saw a correlation between genetic distance and geographical distance (Additional file 2: Figure S1), however this has to be considered with caution as the low level of genetic variation may cause artefacts. Zuriaga et al. [36] extended the analysis of Caicedo and Schaal [35] by including accessions from Ecuador as well, and they suggested that the population structure could be better explained by ecological and climatic variables rather than by geographic distances. Similar results were obtained for populations of Solanum habrochaites[20,21] were geographical and climatic variables explain a substantial amount of the variation in terpenoids and acyl sugars. However, in our analysis, we did not see such correlation between climatically or geographical conditions at the collection sites of the accessions and the morphological species boundaries.

Can selection pressure explain morphological differences?

Starting from the generally accepted assumption [18,37] that maintenance of a constitutive resistance mechanism like the one present in S. galapagense (i.e. trichomes, acyl sugars) is energetically expensive and thus only viable when there is a high selection pressure. One possible explanation to consider is that the resistance mechanisms found in S. galapagense serves other functions as well. In other words, resistance is a secondary function of a trait evolved in response to some other biotic or abiotic pressure [18,38]. For example, a leaf surface characteristic such as high trichome density may be a direct defence against insects or pathogens like bacteria. It may also be an indirect defence mechanism against viruses although no reports mention the occurrence of whitefly transmitted viruses in S. galapagense and/or S. cheesmaniae. Trichomes might be of even greater adaptive value against other environmental stresses such as water loss. Flanders et al. [38], working with potato accessions distributed from USA to Argentina, attributed the observed resistance patterns towards potato herbivores (potato aphid, Colorado potato beetle, potato flea beetle and potato leafhopper) to variation in geographical (altitude) and climatic conditions. Specifically, they found that accessions from hot and dry areas were more resistant to Colorado potato beetle and potato flea beetle. However, considering that S. galapagense and S. cheesmaniae are found in sympatry in space and time [3,5] and that within our set of accessions there was no significant difference between S. galapagense and S. cheesmaniae for any of the climatic/geographical variables, it is unlikely that maintenance of the resistance mechanisms could be explained by this hypothesis.

Another plausible explanation is that resistance traits in tomato (trichomes, secondary metabolites) are under selective pressure. It was observed that S. galapagense is more abundant and widespread in undisturbed areas than S. cheesmaniae[5]. Although it has being proven that trichomes and their exudates also have an effect against other pathogens [39], we have focused our analysis on whiteflies. Though, B. tabaci was only recently found on the Galapagos Islands (1998) [9], numerous other herbivores sharing the same feeding guilds (i.e. Myzus persicae, Macrosiphum sp.) are present [9,40]. Recently it was shown that even the selective pressure of a single insect herbivore species (either a generalist or a specialist) can be strong enough to shift the allele frequencies of a plant population within a few generations [41-43]. Züst et al. [43], working with Arabidopsis populations and aphids (a generalist phloem-feeder), reported that the trichome densities of Arabidopsis plants in presence of herbivores remains constant over generations and decrease in an herbivore free environment. The authors also proved that aphid populations had an effect on the frequency of the plant genotypes with different aliphatic glucosinolates. It was also proven that there was an indirect effect of herbivory selection pressure in the sense that some accessions were able to compete with others only in the absence of herbivores, showing the ecological benefits of having a resistance trait only in the presence of herbivores. Agrawal et al. [41] provided evidence that herbivory by a specialist insect pest (seed predator moth) can act as a direct selective force in favour of the resistance, but also as an indirect selective force to enhance competitive ability of the plant in the presence of herbivores. All this information provides evidence to support the hypothesis that the current resistance mechanisms present in S. galapagense could be maintained by selection pressure, rather than by geographical/climatic variables. This hypothesis can explain how it is possible for the plant to maintain a high level of morphological differentiation (resistance level, trichome composition) with a relatively low genetic variation and in the presence of gene flow.

Conclusions

Our results show that whitefly resistance was found exclusively in S. galapagense accessions and that it was associated with the presence of type IV trichome and high levels of acyl sugars. Our marker and metabolomics data support the hypothesis that S. galapagense and S. cheesmaniae might be morphotypes rather than two species and that their co-existence is likely the result of a selective pressure.

Methods

Plant materials and growing conditions

In total, we evaluated 35 tomato accessions (Table 1, Additional file 1: Table S1), covering the geographical distribution as much as possible (Figure 5). These included 22 accessions of Solanum cheesmaniae L. Riley, 12 accessions of S. galapagense S.C. Darwin & Peralta, and 1 accession of S. lycopersicum L. as reference.

thumbnailFigure 5. Geographical distribution of the tomato accession on the Galapagos Islands.Solanum galapagense S.C. Darwin & Peralta (violet circle symbol); S. cheesmaniae L. Riley (red plus sign symbol). Map made with DIVA-GIS version 7.5 [44].

The accessions were grown in a greenhouse at Wageningen UR Plant Breeding, Wageningen, the Netherlands (20 ± 2°C, 70% RH, 16/8 h day/night) in 14 cm pots. The plants were fertilized twice a week and watered once a day.

When the plants were six weeks old, they were moved to an insect proof greenhouse. One week before infestation the greenhouse temperature was increased slowly (two degrees per day) from 20 till 27°C to allow plants to adapt to the higher temperature (27 ± 2°C, 70% RH, 16/8 h day/night).

Insect rearing

A non-viruliferous whitefly rearing (Bemisia tabaci Group Mediterranean-Middle East-Asia Minor I) was maintained on the tomato cultivar Moneymaker for several generations at Wageningen UR Plant Breeding, Wageningen, The Netherlands. The initial inoculum was obtained from a permanent rearing at the Laboratory of Entomology, Wageningen UR, Wageningen, The Netherlands.

No-choice test

Whiteflies were anesthetized using CO2 and four days old females were selected under a binocular microscope by the morphology of the abdomen. Five females were placed into a clip-on cage (2.5 cm in diameter and 1.0 cm in high). Three clip-on cages per plant and four plants per accession were used. The cages were placed on the first to third fully expanded leaf, counting from the top of the plant, thereby taking care not to break the plant trichomes when assembling the cages. Five days after infestation, the number of death and alive whiteflies, as well as the number of eggs was registered from that day on. The surviving whiteflies were removed from the leaves and the adult survival (AS) and oviposition rate (OR) were calculated according to Bas et al. [45]. Due to the fact that the hatching of the eggs was irregular in time, we could not assess the number of newly hatched insects per day in order to calculate the development period (DP). When almost no new whiteflies were seen, the number of empty pupae was recorded and the pre-adult survival (PS) calculated [45]. A complete randomize assay with four replicas per accession was used. Each replica consisted of the average value of three cages per plant (technical replica). The variables were analysed with a one-way ANOVA followed by a least significant difference (LSD) test [46]. An Arcsin (Sqrt) transformation was applied to the variables adult survival (AS) and pre-adult survival (PS), whereas an Sqrt (x + 1) transformation was applied to the variable oviposition rate (OR). All statistical procedures were performed using the statistical software package Infostat Professional (2010) Cordoba, Argentina.

Trichome description

Trichomes present on the abaxial side of the leaf were classified according to type [47]. For an estimation of trichome density, the abaxial part of three leaflets was observed under the binocular microscope and a visual scale was used to describe it. The scale used was adapted from Simmons and Gurr [16] 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.

Genotyping and phylogenetic analysis

Genomic DNA was extracted from one randomly selected plant per accession as described by Fulton [48]. The DNA concentration was adjusted to 50 ng/μl. For marker analysis a custom made single nucleotide polymorphism (SNP) Infinium bead array was used [23]. On this array 5528 tomato SNPs were present. The marker analysis was performed using the protocol provided by Infinium and carried out by Service XS, Leiden, The Netherlands. Marker data obtained were filtered using the following criteria: 1) monomorphic markers were removed; 2) markers were deleted when the frequency of heterozygotes per marker was equal to or higher than 5%, and 3) markers were deleted when the frequency of no-calls (NC) was equal to or higher than 50%. We choose a threshold of 50% for the NC because they can be the result of either an amplification problem (true NC) or the presence of a different, informative allele (false NC). The latter situation occurs frequently with more distantly related species [23]. After filtering a total of 3316 markers were used in the analysis. A phylogenetic tree was reconstructed by Neighbour Joining, using the Manhattan distance (Additional file 1: Table S1). The reliability of the resulting dendrogram was assessed by bootstrap analysis with 1000 replications. The analysis was carried out using the software package PAST [49]. Isolation By Distance (IBD, version 3.23) [50] was used to analyse for presence of isolation by distance.

LC-QTOF-MS analysis

Four plants per accession were used for the chemo-profiling. From each plant, one complete leaf (second fully expanded leaf from the top of the plant) was cut, placed into an aluminium envelope and immediately frozen in liquid nitrogen. Each sample was frozen with liquid nitrogen and ground to a fine powder and storage at -80°C until use. Extraction and analysis by Liquid Chromatography-Quadrupole Time of flight-Mass Spectrometry (LC-QTOF-MS) [51,52]. Four hundred mg of frozen leaf powder was put into a glass tube with 1.2 mL of methanol/formic acid solution (99.9% - 0.1%). The samples were mixed using vortex, sonicated for 15 min and centrifuged at 2500 rpm for 10 min. The supernatant was filtered using a 0.45 μm filter, injected (5 μl) using an Alliance 2795 HT instrument (Waters), separated on a Phenomenex Luna C18 (2) column (2.0 × 150 mm, 3 mm particle size) using a 5–95% gradient of acetonitrile in water (both acidified with 0.1% formic acid) in 45 min and then detected by a Water-Micromass QTOF Ultima MS with electrospray ionization in negative mode (m/z 80–1,500). Annotation of LCMS peaks corresponding to acyl sugars was done on their accurate masses as previously described in Firdaus et al. [22]. The Quanlynx tool of the Masslynx acquisition software was used to calculate the relative abundance (peak area) of the different type of acyl sugars for all samples, based on their specific mass and retention times. The three isoforms of S3-20 (I, II, III) and the four of S3-22 (III – VI), as described before in Firdaus et al. [22], were considered as single peaks for S3-20 and S3-22, respectively, due to partly overlapping chromatographic peaks of these isomeric compounds in many samples. The peak areas were LOG2 (× + 1) transformed and auto scaled to the mean. PCA-biplot was done with the software package PAST [49], and all other statistical analysis were performed with the software package GeneMaths XT (version 2.21; Applied Maths, Belgium).

Analysis of geographical distribution and climate data

The collection site information (latitude and longitude) of the accessions was obtained from the Tomato Genetic Resource Centre (TGRC). Locations and climate data were obtained from the WorldClim at 2.5 arc-min resolution [53] and were analysed as proposed by Gonzales-Vigil et al. [20], and using DIVA-GIS [44]. A regression analysis was carried out between the resistance variables (AS, OR, PS), the metabolites and the climatic/geographical variables (Altitude, Latitude, Longitude, Ann. Mean precipitation, Ann. Mean temperature).

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The project was conceived by AFL, AWvH and BV. AFL performed the experiments and data analysis. RCHdV has been responsible for the LC-MS analyses and the data processing. AFL, AWvH, BV, RCHdV and RGFV have been involved in the writing of the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We would like to thank the Tomato Genetic Resource Centre (TGRC) and the Centre for Genetic Resources (CGN) for providing plant materials. We would like to give special thanks to Jarno Sinnige, Marcela Víquez for their help during the development of this research, to Harry Jonker and Bert Schipper for their help in the LC-MS analyses. Roeland Voorrips for his advice on the statistical analysis and to Colette Broekgaarden, Sergio Rasmann and the two anonymous reviewer for critically reading and valuable comments on earlier versions of this manuscript. RCHdV was financially supported by the Netherlands Metabolomics Centre and the Centre for BioSystems Genomics, both of which are part of the Netherlands Genomics Initiative / Netherlands Organization for Scientific Research. AFL was financially supported by the Foundation CAPACIT-AR del NOA and by Wageningen UR Plant Breeding.

References

  1. Darwin SC, Knapp S, Peralta IE: Taxonomy of tomatoes in the Galápagos Islands: native and introduced species of Solanum section Lycopersicon (Solanaceae).

    Syst Biodivers 2003, 1(01):29-53. Publisher Full Text OpenURL

  2. Peralta IE, Spooner DM, Knapp S: Taxonomy of wild tomatoes and their relatives (Solanum sect. Lycopersicoides, sect. Juglandifolia, sect. Lycopersicon: Solanaceae). Ann Arbor, MI: American Society of Plant Taxonomists; 2008. OpenURL

  3. Darwin SC: The systematics and genetics of tomatoes on the Galápagos Islands (Solanum, Solanaceae). University College London: Ph.D. Thesis; 2009. PubMed Abstract | Publisher Full Text OpenURL

  4. Mallet J: A species definition for the modern synthesis.

    Trends Ecol Evol 1995, 10(7):294-299. PubMed Abstract | Publisher Full Text OpenURL

  5. Nuez F, Prohens J, Blanca JM: Relationships, origin, and diversity of Galápagos tomatoes: implications for the conservation of natural populations.

    Am J Bot 2004, 91(1):86-99. PubMed Abstract | Publisher Full Text OpenURL

  6. Peralta IE, Knapp S, Spooner DM, Lammers TG: New species of wild tomatoes (Solanum Section Lycopersicon: Solanaceae) from northern Peru.

    Syst Bot 2005, 30(2):424-434. Publisher Full Text OpenURL

  7. Spooner DM, Peralta IE, Knapp S: Comparison of AFLPs with other markers for phylogenetic inference in wild tomatoes [Solanum L. section Lycopersicon (Mill.) Wettst.].

    Taxon 2005, 54(1):43-61. Publisher Full Text OpenURL

  8. Walsh SJ, Mena CF: Science and Conservation in the Galapagos Islands : Frameworks & Perspectives. New York, NY: Springer New York; 2013. OpenURL

  9. Causton CE, Peck SB, Sinclair BJ, Roque-Albelo L, Hodgson CJ, Landry B: Alien insects: threats and implications for conservation of Galápagos Islands.

    Ann Entomol Soc Am 2006, 99(1):121-143. Publisher Full Text OpenURL

  10. Valverde RA, Sim J, Lotrakul P: Whitefly transmission of sweet potato viruses.

    Virus Res 2004, 100(1):123-128. PubMed Abstract | Publisher Full Text OpenURL

  11. 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. Publisher Full Text OpenURL

  12. Dinsdale A, Cook L, Riginos C, Buckley YM, Barro PD: Refined global analysis of Bemisia tabaci (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae) Mitochondrial Cytochrome Oxidase 1 to identify species level genetic boundaries.

    Ann Entomol Soc Am 2010, 103(2):196-208. Publisher Full Text OpenURL

  13. De Barro PJ: The Bemisia tabaci species complex: questions to guide future research.

    Journal of Integrative Agriculture 2012, 11(2):187-196. Publisher Full Text OpenURL

  14. Firdaus S, Vosman B, Hidayati N, Jaya Supena ED, Visser RG, Van Heusden AW: The Bemisia tabaci species complex: additions from different parts of the world.

    Insect Sci 2012, 00:1-11. OpenURL

  15. Liu S-s, Colvin J, De Barro PJ: Species concepts as applied to the whitefly Bemisia tabaci systematics: how many species are there?

    Journal of Integrative Agriculture 2012, 11(2):176-186. Publisher Full Text OpenURL

  16. Simmons AT, Gurr GM: Trichomes of Lycopersicon species and their hybrids: effects on pests and natural enemies.

    Agricultural and Forest Entomology 2005, 7(4):265-276. Publisher Full Text OpenURL

  17. 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. Publisher Full Text OpenURL

  18. Rausher MD: Co-evolution and plant resistance to natural enemies.

    Nature 2001, 411(6839):857-864. PubMed Abstract | Publisher Full Text OpenURL

  19. Weber MG, Agrawal AA: Phylogeny, ecology, and the coupling of comparative and experimental approaches.

    Trends Ecol Evol 2012, 27(7):394-403. PubMed Abstract | Publisher Full Text OpenURL

  20. Gonzales-Vigil E, Hufnagel DE, Kim J, Last RL, Barry CS: Evolution of TPS20-related terpene synthases influences chemical diversity in the glandular trichomes of the wild tomato relative Solanum habrochaites.

    Plant J 2012, 71(6):921-935. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Kim J, Kang K, Gonzales-Vigil E, Shi F, Daniel Jones A, Barry CS, Last RL: Striking natural diversity in glandular trichome acylsugar composition is shaped by variation at the acyltransferase2 locus in the wild tomato Solanum habrochaites.

    Plant Physiol 2012, 160(4):1854-1870. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Firdaus S, Heusden A, Hidayati N, Supena E, Mumm R, Vos RH, Visser RF, Vosman B: Identification and QTL mapping of whitefly resistance components in Solanum galapagense.

    Theor Appl Genet 2013, 126(6):1487-1501. PubMed Abstract | Publisher Full Text OpenURL

  23. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Jacobs MJ, Vosman B, Vleeshouwers VAA, Visser RF, Henken B, Berg R: A novel approach to locate Phytophthora infestans resistance genes on the potato genetic map.

    Theor Appl Genet 2010, 120(4):785-796. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Jansky SH, Simon R, Spooner DM: A test of taxonomic predictivity: Resistance to white mold in wild relatives of cultivated potato.

    Crop Sci 2006, 46(6):2561-2570. Publisher Full Text OpenURL

  26. Lokossou AA, Rietman H, Wang M, Krenek P, van der Schoot H, Henken B, Hoekstra R, Vleeshouwers VGAA, van der Vossen EAG, Visser RGF, et al.: Diversity, distribution, and evolution of Solanum bulbocastanum late blight resistance genes.

    Mol Plant Microbe Interact 2010, 23(9):1206-1216. PubMed Abstract | Publisher Full Text OpenURL

  27. 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. Publisher Full Text OpenURL

  28. Leckie B, Jong D, Mutschler M: Quantitative trait loci increasing acylsugars in tomato breeding lines and their impacts on silverleaf whiteflies.

    Mol Breeding 2012, 30(4):1621-1634. Publisher Full Text OpenURL

  29. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Glas JJ, Schimmel BCJ, Alba JM, Escobar-Bravo R, Schuurink RC, Kant MR: Plant glandular trichomes as targets for breeding or engineering of resistance to herbivores.

    Int J Mol Sci 2012, 13(12):17077-17103. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Kimura S, Koenig D, Kang J, Yoong FY, Sinha N: Natural variation in leaf morphology results from mutation of a novel KNOX gene.

    Curr Biol 2008, 18(9):672-677. PubMed Abstract | Publisher Full Text OpenURL

  32. Nakazato T, Franklin RA, Kirk BC, Housworth EA: Population structure, demographic history, and evolutionary patterns of a green-fruited tomato, Solanum peruvianum (Solanaceae), revealed by spatial genetics analyses.

    Am J Bot 2012, 99(7):1207-1216. PubMed Abstract | Publisher Full Text OpenURL

  33. Nakazato T, Housworth EA: Spatial genetics of wild tomato species reveals roles of the Andean geography on demographic history.

    Am J Bot 2011, 98(1):88-98. PubMed Abstract | Publisher Full Text OpenURL

  34. Nakazato T, Warren DL, Moyle LC: Ecological and geographic modes of species divergence in wild tomatoes.

    Am J Bot 2010, 97(4):680-693. PubMed Abstract | Publisher Full Text OpenURL

  35. Caicedo AL, Schaal BA: Heterogeneous evolutionary process affect R gene diversity in natural populations of Solanum pimpinellifolium.

    Proc Natl Acad Sci U S A 2004, 101(50):17444-17449. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Zuriaga E, Blanca J, Nuez F: Classification and phylogenetic relationships in Solanum section Lycopersicon based on AFLP and two nuclear gene sequences.

    Genet Resour Crop Evol 2009, 56(5):663-678. Publisher Full Text OpenURL

  37. Strauss SY, Rudgers JA, Lau JA, Irwin RE: Direct and ecological costs of resistance to herbivory.

    Trends Ecol Evol 2002, 17(6):278-285. Publisher Full Text OpenURL

  38. Flanders K, Radcliffe E, Hawkes J: Geographic distribution of insect resistance in potatoes.

    Euphytica 1997, 93(2):201-221. Publisher Full Text OpenURL

  39. Nonomura T, Xu L, Wada M, Kawamura S, Miyajima T, Nishitomi A, Kakutani K, Takikawa Y, Matsuda Y, Toyoda H: Trichome exudates of Lycopersicon pennellii form a chemical barrier to suppress leaf-surface germination of Oidium neolycopersici conidia.

    Plant Sci 2009, 176(1):31-37. Publisher Full Text OpenURL

  40. Wilson MR, Peck SB, Causton C: CDF Checklist of Galapagos true bugs - leafhoppers, planthoppers, aphids and scale insects. In Charles Darwin Foundation Galapagos Species Checklist. Edited by Herrera H, Jaramillo P, Tirado N, Jiménez-Uzcátegui G, Ruiz D, Guézou A, Ziemmeck F. Puerto Ayora, Galapagos: Charles Darwin Foundation/Fundación Charles Darwin; 2013. OpenURL

  41. Agrawal AA, Hastings AP, Johnson MTJ, Maron JL, Salminen JP: Insect herbivores drive real-time ecological and evolutionary change in plant populations.

    Science 2012, 338(6103):113-116. PubMed Abstract | Publisher Full Text OpenURL

  42. Hare JD: How insect herbivores drive the evolution of plants.

    Science 2012, 338(6103):50-51. PubMed Abstract | Publisher Full Text OpenURL

  43. Züst T, Heichinger C, Grossniklaus U, Harrington R, Kliebenstein DJ, Turnbull LA: Natural enemies drive geographic variation in plant defenses.

    Science 2012, 338(6103):116-119. PubMed Abstract | Publisher Full Text OpenURL

  44. Hijmans RJ, Guarino L, Mathur P: DIVA-GIS a geographic information system for the analysis of biodiversity data. Manual. Version 7.5. Lima, Peru: International Potato Center; 2012. OpenURL

  45. 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. Publisher Full Text OpenURL

  46. Zar JH: Biostatistical analysis. Upper Saddle River, NJ: Pearson Education International; 2010. OpenURL

  47. 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. Publisher Full Text OpenURL

  48. Fulton T, Chunwongse J, Tanksley S: Microprep protocol for extraction of DNA from tomato and other herbaceous plants.

    1995, 13(3):207-209.

  49. Hammer O, Harper DAT, Ryan PD: PAST: Paleontological statistics software package for education and data analysis.

    Palaeontol Electron 2001, 4(1):1-9. OpenURL

  50. Jensen J, Bohonak A, Kelley S: Isolation by distance, web service.

    BMC Genetics 2005, 6(1):13. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. De Vos RCH, Moco S, Lommen A, Keurentjes JJB, Bino RJ, Hall RD: Untargeted large-scale plant metabolomics using liquid chromatography coupled to mass spectrometry.

    Nature Protocols 2007, 2(4):778-791. PubMed Abstract | Publisher Full Text OpenURL

  52. Firdaus S, Heusden A, Hidayati N, Darmo E, Supena EDJ, Mumm R, De Vos RCH, Visser RG, Vosman B: Identification and QTL mapping of whitefly resistance components in Solanum galapagense.

    Theor Appl Genet 2013, 126:1487-1501. PubMed Abstract | Publisher Full Text OpenURL

  53. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A: Very high resolution interpolated climate surfaces for global land areas.

    International Journal of Climatology 2005, 25(15):1965-1978. Publisher Full Text OpenURL