Dissecting the Genotypic and Environmental Factors Underpinning the Quantitative Trait Variation in a Set of Wild Tomato (Solanum habrochaites LA1777) Introgression Lines
by
, and
Pasquale Tripodi
1,*
,
Antonella Vitiello
2, 
Bruno D’Onofrio
1,
Mario Parisi
1
and
Maria Cammareri
2 
1
CREA Research Centre for Vegetable and Ornamental Crops, 84098 Pontecagnano (SA), Italy
2
Institute of Biosciences and BioResources (CNR-IBBR), National Research Council of Italy, Via Universita’ 133, 80055 Portici (NA), Italy
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(1), 38; https://doi.org/10.3390/agronomy11010038
Submission received: 20 November 2020
/
Revised: 22 December 2020
/
Accepted: 24 December 2020
/
Published: 27 December 2020
(This article belongs to the Special Issue Application of Wild-Species Germplasm Genetic Variation in Crop Improvement)
Abstract
:Exotic libraries have proven a powerful tool for the exploitation of wild relatives and quantitative trait loci (QTLs) detection in crop species. In early 2000, an introgression line (IL) population of the wild tomato Solanum habrochaites (SH) acc. LA1777 was developed and made publicly available. Despite the potentiality of the donor parent, so far, these lines have been poorly explored for their agronomic performance and for the identification of genomic regions underlying the variation of quantitative traits (QTLs). Here, we report the evaluation of 19 morpho-agronomic and chemical traits on a set of 39 ILs grown in three consecutive field seasons with the aim to: (a) Determine the overall phenotypic performances of the studied collection, (b) estimate the influence of the genotype (G) and the year of cultivation (Y) and their interaction on the traits analyzed, (c) investigate the plasticity of the traits, and (d) identify whole-genome QTLs in the wild SH background. The ILs showed lower productivity compared to the control genotype, while no major effects were found for the morphological fruit-related traits. Instead, a general increase in the soluble solids content was observed. The combined analysis of G × Y highlighted a major effect of the genotype on trait variation, although yield-related traits were more influenced by environmental factors. In total, 75 associations for 17 traits were detected. Major QTLs increasing soluble solids, pericarp thickness, and trichome density were respectively found on chromosomes 1, 5, and 11 with a percentage variation (PV) of 24.01%, 32.49%, and 200%. Furthermore, different QTLs increasing the color intensity and fruit shape were detected. These results suggest that SH could be a potential source of favorable alleles for qualitative traits despite its inferior phenotype compared to the cultivated parent. The evaluated set of SH LA1777 ILs is a potential for novel allele discovery in wild tomatoes and for breeding purposes towards the exploitation of the available introgressions and for the pyramiding of traits.
1. Introduction
Tomatoes (Solanum lycopersicum, Solanaceae; 2n = 2x = 24) are considered one of the world’s most important vegetable crops widely spread with a production estimation of ~180 million tons and particularly appreciated for their nutritional properties [1]. As a model in plant research, the tomato has been pioneering for exploiting exotic germplasms toward the dissection of the genetic factors underpinning the variations of complex traits. Since the early nineties, the interspecific hybridization of S. lycopersicum with wild relatives has led to the development of various recombinant mapping populations able to enhance the identification of loci involved in agronomic and quality traits, as well as in biotic and abiotic stresses resistance [2,3]. Among the different types of experimental populations, introgression lines (ILs) were the most successful for allele transfer in precision breeding, providing, furthermore, an excellent material for quantitative trait loci (QTLs) mapping. ILs harbor a single marker-defined chromosome segment from the donor wild parent in the cultivated background [4]. The accuracy for mapping polygenes relies on the direct association of the phenotypic variation to the introgressed segment, minimizing the epistatic effects from other wild genome regions. Ideally, an IL population is supposed to cover the whole wild genome and is a permanent resource, being developed after subsequent backcross and self-fertilization cycles. This allows tests across multiple locations and seasons, providing a powerful tool to estimate the genotype (G) and the environmental (E) effects. In fact, dissection of the G and E factors and their interaction is essential for determining those traits highly inherited in the offspring and for the estimation of phenotypic plasticity as the capacity of a genotype to express different phenotypes according to different environmental conditions [5]. Therefore, ILs play a key role for successful genetic improvement programs, facilitating, furthermore, the pyramid breeding of favorable genes [6]. In this context, the genetic variability enclosed in wild relatives can be exploited to address the future challenges of agriculture, such as the growth of the global population, the adaptation to climate changes, and the demand for quality and sustainability for food supply [7]. In recent years, ILs and near-isogenic lines (NILs) sets were developed in tomatoes through hybridization with different wild species, including S. pennellii [8], S. habrochaites [9,10,11,12], S. lycopersicoides [13,14], S. pimpinellifolium [15], and, more recently, S. sitiens [16].
Solanum habrochaites S. Knapp and D.M. Spooner (previously, Lycopersicum hirsutum Dunal) is a potential source of beneficial alleles in tomato breeding, including quality traits [17], chilling tolerance [18,19], and disease resistances [3,20]. There are over 90 accessions distributed in the native origin area for the species [21]; among these, acc. LA1777 is a green-fruited accession widely studied for resistance to the main tomato and potato pathogens, including begomoviruses and fungi [22,23,24]. It is also reported as resistant to insect pests, thanks to the presence of trichomes and of secondary metabolites produced by their glandules [25,26,27]. Finally, this wild species is a valuable source of beneficial alleles involved in fruit ripening [28] and an increase of soluble solids [10,29,30]. Monforte and Tanksley [9] developed a set of 99 near-isogenic and backcross recombinant inbred lines from a cross between a processing tomato cultivar E6203 and S. habrochaites LA1777. The population provides coverage of more than 85% of the wild genome based on genotyping with ~93 RFLP (Restriction Fragment Length Polymorphism) markers. Despite that the population is publicly available, the enclosed potentialities for mapping plant and fruit traits have not yet been broadly investigated. Indeed, the identification of QTLs and genes affecting yield-related components have been, so far, mainly focused on the S. pennellii ILs, which have been extensively exploited for the identification of over 3000 QTLs [31,32] and for map-based cloning of the first two QTLs in tomatoes: fw2.2 for fruit weight [33] and brix9-2-5 for the soluble solids content [34]. Attempts toward the assessment of the agronomic performance of a subset of LA1777 ILs in hybrid combinations have been reported [35], although no QTL analysis has been performed. Instead, specific ILs were used for the fine mapping of QTLs for agronomic traits on the distal parts of chromosomes 1 and 4 [10,11]. Therefore, a lack of information at the whole-genome level still occurs for the LA1777 introgression library.
In the present study, we report the evaluation of 39 ILs for nineteen agronomic, morphological, and quality-related traits over three consecutive field seasons with the following objectives: (i) Determine the phenotypic performances of the selected set, (ii) estimate the genotypic and environmental effects and the plasticity of the traits, and (iii) identify wild introgressions responsible for the quantitative variations of the traits studied. The set of ILs was chosen based on the minimal number of introgressions able to cover, as much as possible, the whole genome, in agreement with previously reported data [9,36].
2. Materials and Methods
2.1. Plant Material
Plant material consisted of 39 S. habrochaites (acc. LA1777) (hereafter, SH) introgression lines and by the recurrent parent S. lycopersicum cv E6203 (hereafter, SL). The population was developed through an advanced backcross scheme derived from the cross between SL and SH. From the generation of BC2S3 families and further marker-assisted selection, lines containing a single or few segments of the wild relative in the cultivated background were obtained [9].
2.2. Field Trials and Phenotypic Characterization
Plants were grown in independent trials across three consecutive seasons at the experimental farm of the Research Centre for Vegetable Crops located at Battipaglia (Salerno, Italy). Geographical and pedological characteristics of the site and climatic conditions of each season are reported in Table 2.
Field tests were conducted using a completely randomized experimental design with 10 plants from each IL-SH and 40 plants from the control SL. Each year, the sowing was done in March, and the seedlings were transplanted in April in single rows, adopting distances of 100 cm between the rows and 30 cm along the rows. For plant fertilization, 200 kg ha−1 of N, 200 kg ha−1 of P2O5, and 120 kg ha−1 of K2O were applied according to what was established for the same site [37]. Irrigation scheduling was based on crop evapotranspiration. A total of 270, 260, and 290 mm of water were applied in 2014, 2015, and 2016, respectively, by drip irrigation, restoring 40% of the total available water depleted. The effective control of diseases and harmful insects was pursued, with the aim of having healthy plants until the end of the cycle. Other cultivation techniques included stakes as support and galvanized wires [38].
Agronomical and morphological traits were assessed when 95% of the plants of the control genotype (SL) had ripened fruits. Yield-related traits were performed on each plant and included: total yield (TY) (in kg), assessing the total weight of fruits (green and red), red yield (RY) (in kg) as the weight of whole, ripened fruits at commercial status, and green yield (GY) (in kg), including all unripe fruits at harvest time (in kg). Uniformity of fruit maturity (UM) was then measured as follows:
GY = TY − RY,
UM = (RY/TY) × 100,
The weight of the vegetative part (PW) (in kg) was obtained weighing the epigeal portion of the plants after harvesting all fruits; biomass (BM) (in kg) was measured as:
BM = TY + PW,
The fruit traits assessment included the average weight (FW) (in gr) obtained by dividing the total yield by the number of fruits harvested, the length (FL) and the width (FD) (in cm) measured on 10 fruits/plant by using a ruler, and the fruit shape (FS) index as follows:
where values < 1 indicate fruits more flattened, and values > 1 indicate fruits more elongated. External (EC) and internal (IC) colors were measured by scale (1 = yellow fruits and 5 dark red fruits), as well as pericarp thickness (PER) and puffiness (PUF) (1 = thin pericarp, low puffiness and 3 = thick pericarp, high puffiness).
FS = FL/FD,
Chemical traits were measured on a bulk of 10 well-ripened fruits/plants. Fruits were washed and dried and then sliced and homogenized in a Waring blender (2-L capacity; Model HGB140, Waring Commercial, Stamford, CT, USA) for 1 min. Soluble solids content (SSC) was measured using a digital refractometer (Refracto 30PX, Mettler-Toledo, Novate Milanese, Italy), and the results were expressed as °Brix per 100 g of fresh weight (fw). The pH was determined using a pH-Matic 23 titroprocessor equipped with a pH electrode with a temperature sensor (model 5011T) (Crison Instruments, Barcelona, Spain).
The estimation of the yield of soluble solids produced per plant (BY) (in gr) was estimated as follows:
BY = RY × SSC,
Trichome density (TRIC) for each plant was estimated at the base of the stem and on the pedicel of flowers using a visual scale (1 = no or low density of trichomes and 3 = high density of trichomes).
2.3. Data Analyses
All phenotypic traits were subjected to an analysis of variance (ANOVA) test. Trait means for the control and the set of ILs were compared by using Tukey’s HSD (honest significant difference) test. Results with p less than 0.05 were considered statistically significant. Broad-sense heritability (H2) was estimated as:
where “VG” is the genotypic component of variance, and “VP” is the total phenotypic variance. The coefficient of variation (CV) in percentage was calculated as:
in which “SD” is the standard deviation, and “m” is the mean for a trait.
H2 = VG/VP,
CV = (SD/m) × 100,
A two-factorial linear model was used to determine the main effects of the genotype (G), season of cultivation (Y), and their interactions (G × Y) for each trait. Mean square values (MS) were used to estimate the magnitude of the observed effect, while the total sum of squares in percentage (TSS%) was calculated by dividing the TSS of the variable by the total TSS.
The coefficient of relative phenotypic plasticity (CRP) for each trait was calculated according to the formula of Dingemanse et al. [39] as follows:
where “Vi” is the variance of the individual (i), and “VPp” is the overall phenotypic variance of the population. The CRP index close to 0 indicates no plasticity (high stability); on the contrary, higher values indicate an increase of plasticity (decrease of stability).
CRP = Vi/VPp,
Correlations across the genotypes for phenotypic traits were calculated using the Pearson’s test at p < 0.01 after Bonferroni’s correction for multiple comparisons. The correlogram was constructed and visualized using the Corrplot package implemented in R version 3.0.2 (R Development Core Team). Principal component analysis (PCA) was carried out to determine what are most effective descriptors in discriminating among accessions using the computer package XLSTAT 2012.1 and visualizing the similarities among accessions. Experimental data were statistically elaborated using R. The relationships between pairs of phenotypic and molecular data matrices were computed by the Mantel test using Pearson’s r-value [40].
Chromosomal regions responsible for the quantitative variation for traits measured were detected comparing the mean phenotypic values of the tested ILs in each experiment to the recurrent parent SL using Dunnett’s test with Type I error α ≤ 0.05 [41]. When an IL showed significantly different effects compared to the control in one or more seasons, putative quantitative trait loci (QTLs) were assumed in the chromosomal region covered by the IL. If a QTL was detected in multiple ILs, the likely position was inferred on overlapping chromosome segments. For each IL, the percentage variation (PV) compared to the control, and underlying a QTL region was calculated as:
where, for each trait, “m” is the mean calculated across all replications.
PV = [m(IL) − m(E6203)/m(E6203)] × 100,
3. Results
3.1. Phenotypic Variation and Trait Performances of SL and SH ILs
In Table 3, the ANOVA results and the parameters mean, minimum, maximum, coefficient of variation (CV), and broad sense of heritability (H2) specify the performance of the recurrent parent SL “E6203” and the 39 SH ILs set for each trait. Considering the average of three years of cultivation, no significant differences between the sets of ILs and SL were found for GY, UM, PW, IC, EC, and TRIC, while PUF and SSC, as well as FS, differed at p < 0.05 and p < 0.01, respectively. In average, the wild SH alleles conferred less productivity when introgressed in the SL background, while no major effect was noticed for the morphological fruit-related traits. Instead, a general increase in the soluble solids content was detected in the tested ILs, although no effect was found for the production of the soluble solid. Across all years, the CV% was lower than 36% for fruit traits in both SH ILs and the recurrent parent. A higher variation was shown by yield-related traits, reaching peaks for GY. Overall, the lowest coefficient of variation was found for FS, highlighting very low or absent morphological segregation within the ILs. The broad-sense heritability ranged from 0.84 (GY and CI) to 0.96 (FS). The fruits and chemical traits showed, on average, a greater heritability (Table 3).
The mean, ranges, CV, and results of the post hoc Tukey’s considering each independent season in both SL and the set of SH ILs are reported in Table S1. Highly significant differences (p < 0.001) between the ILs and the recurrent parent were found for 17 out of the 19 traits analyzed. Only the fruit shape and trichomes differed at p < 0.01 and p < 0.05, respectively (Table S1).
The average phenotypic values and the percentage variation compared to the control revealed a large variability of the studied ILs set for the considered traits (Tables S2 and S3).
In all trials, the introgression lines were less productive than the control. Overall, a decrease of the yield was observed in both E6203 and the ILs from the first to the third season. The same trend was observed for plant weight and biomass, with the lowest values exhibited during Y3. On the contrary, the weight of the fruits was higher during the second year of cultivation with respect to the other seasons, showing greater values in the control genotype across all trials (Table S2). The soluble solids content was found to be higher than the control in almost all lines during the three-year evaluation period (Figure 1). Only for Y2, variable levels were observed. Overall, the soluble solids values reached contents up to +55% with respect to E6203.
3.2. Genotypic and Environmental Components Underlying Trait Variation
The results of the combined analysis of variance for the traits evaluated across three years of cultivation are given in Table 4. A high level of significance (p < 0.01) was found for the genotypic component (G), while the variation due to the cultivation season (Y) was lower for trichomes and fruit characteristics showing p < 0.05 for FS and CE, and there was no significance for the CI. As for the G × Y interaction, highly significant differences (p < 0.01) were found in all traits, except for FW and pH, which showed a p threshold < 0.05 and no significance, respectively.
For most traits, the main source of variability was due to G, which accounted, on average, for 61.90% of the total variation, expressed by TSS%, ranging from 12.35% for GY to 90.08% for TRIC. The traits under a strong genotypic effect were FL, CE, and CI, exhibiting TSS values above 80%. On average, G accounted for 43.01%, 75.28%, and 60.94% for yield-related traits, fruit characteristics, and chemical features, respectively. The partitioning of the TSS% in the other components affecting the variation indicated that Y and G × Y accounted for, on average, 16.47% and 21.77%, respectively, resulting in, generally, less influence of Y and G × Y than G on the analyzed parameters.
Among the assessed traits, GY and RY were the traits most affected by seasonal fluctuations, showing a TSS of 65.21% and 37.72%, respectively, while for all fruit traits, except puffiness, the effect was very little, with a variation of 1.96% on average. The G × Y interaction was, on average, higher for yield traits (TSS equal to 25.91%) with respect to fruit characteristics (TSS equal to 20.99%) and their main chemical features (TSS equal to 18.46%). Among these, HI and PER were those for which the variation was highly influenced by the G × Y interaction, with values of 33.30% and 36.02%, respectively.
The index of phenotypic plasticity (CRP) (Table S4) was highly variable for the traits tested in the 39 ILs and the control, ranging from 0.05 (UM) to 5.40 (FW). On average, the fruit shape and color traits were more stable, with average plasticity values below 0.62. The yield traits, instead, exhibited higher plasticity, with average indexes above 0.79. In the genotypes tested, most of the values were lower than 1.00. Several ILs had CRP less than 1.00 for all the assessed traits, such as TA1128, TA1536, TA1537, and TA1312. We instead found values above 2.00 for TA1105 on chromosome 2 and for several lines on chromosomes 2, 4, and 5. Only TA1316 (chromosome 8) had CRP values < 0.93 for all traits. The yield traits exhibited high plasticity for TA1276 (HI = 4.68), TA1542 (PW = 4.16), and TA1562 (UM = 4.88), while the maximum value (5.40) was found for FW in TA1548. Instead, the fruit shape and color traits had lower values of plasticity, with an average index below 0.61.
3.3. Multivariate Analysis and Correlations between Traits
The principal component analysis (PCA) in the first two dimensions explained 50.94% of the total variance (Figure 2).
The genotypes were evenly distributed on the two axes of the biplot. The first component accounted for 27.95% of the total variance and was positively correlated with all the traits except SSC and BY. The second component accounted for 22.98% of the variance, being positively correlated with all yield-related traits and negatively correlated with the remaining ones. For the entire set analyzed, the total phenotypic variation was explained by the first 19 components (data not shown). The yield-related traits explained the largest part of the variation in the first two components, showing a contribution of the variables with values above 12% (Table S5). The projection of genotypes on the two-dimensional PCA graph showed the wide range of phenotypic variations enclosed in the SH IL set, highlighting how wild alleles can improve diverse traits. For each growing season, correlograms between pairs of variables were generated using Bonferroni correction (p < 0.01) (Figure 3). Correlations mainly occurred among the same trait categories. In all cultivation trials, strong negative correlations were found for UM and GY, PW and HI, and PW and UM, while strong positive correlations were found between FD and FL, FW and FD, CE and CI, and RY and HY, as well as among BY, BM, TY, and RY. Only in the first and second years of cultivation, TRIC was correlated negatively to UM and positively to GY (Figure 3a,b), while, during the third year, it was negatively correlated with BM (Figure 3c). Furthermore, a significant positive correlation between FW and PER was detected in the third season. In all trials, the soluble solids were negatively correlated with the yield traits and harvest index, although a significance was only observed during the last year. The degree of relationship between the elements of the phenotypic matrices revealed significant correlations between the data of the independent growing seasons, highlighting the strongest correlations between the second and third seasons of the evaluation: Y1/Y2 (r = 0.2, p < 0.0001), Y1/Y3 (r = 0.361, p < 0.0001), and Y2/Y3 (r = 0.607, p < 0.0001).
3.4. Genetic Regions Underlying the Phenotypic Variations
The overlaps between the introgressed fragments of the studied ILs determined the chromosomal regions (or BIN) carrying the wild alleles into the cultivated background. The size of each BIN was determined considering the overlapping regions, plus the half-intervals between the flanking markers and nearest mapped loci. The markers’ positions were updated based on the Tomato-EXPEN reference map [32]. In total, 61 BINs were detected, ranging from a minimum of three on chromosomes 2, 5, 7, 8, 9, and 10 to a maximum of seven on chromosome 1 (Figure 4). Based on the mean comparisons of each ILs with the recurrent parent (Table S2), we detected 75 associations for all traits except for puffiness and green yield. Forty-two associated regions were found across two independent seasons, whereas 33 were across all trials, covering 17 genomic regions on nine chromosomes and including sixteen out of the nineteen evaluated traits. Except for chromosome 8, wild alleles were responsible for the quantitative variations for the studied traits in two or three trials (Table 5).
3.4.1. Yield Traits
For the total yield and red yield, in all cases, the wild alleles were responsible for the reduction of productivity, with the PV ranging from −52.35% to −33.51% for TY and from −52.30% to −41.28% for RY. In total, six QTLs were found to be significant in two growing seasons, while three robust associations were found across the three years, being located in both the short and long arm of chromosome 1 (ty1.2 and ry1.3) and on chromosome 10 (ry10.1).
For the uniformity of maturity, two regions on chromosomes 1 and 4 were responsible for a decrease of the traits, with a PV −15% in both instances. The wild alleles exerted the highest effect on chromosome 1, reducing the trait across all the seasons.
As concerns the biomass and harvest index, a total of nine QTLs were detected, with the larger effect of the wild alleles on the former trait leading to a reduction ranging from −52.03% to −40.24%, whereas, for the latter, the PV ranged from −25.52% to −12.14%. Two QTLs had significant effects during the three consecutive seasons, being responsible for the reduction of the traits for −40.24% (bm1.2) and −25.52% (hi3.4). On the contrary, only for the plant weights, the wild alleles located in a 15-cM chromosomal region on chromosome 3 were responsible for an increase (PW3.6), with a PV of 73.92%.
3.4.2. Fruit Traits
Numerous alleles were found to be responsible for the variations of the fruit traits. The fruit weight was affected by eight significant QTLs; in all cases, the SH alleles reduced the fruit sizes. The genomic regions underlying the FW were identified on seven chromosomes. Five QTLs were robust across the three years, with the strongest effect found for fw11.1, which accounted for a percentage of the phenotypic variance lower than −43.14%. The lowest effect was found for fw3.5, with a PV of −27.91%. Three QTLs were identified in two seasons on chromosomes 9, 10, and 11, with PV% ranging from −35.63% (fw10.3) to −32.01% (fw11.3).
For fruit morphological traits, a total of nine and seven QTLs were found to be significantly associated with the fruit lengths and fruit widths, respectively. The SH alleles decreased the two dimensions in all cases, leading to fruits with small sizes. Robust QTLs were found across the three years. A robust cluster was found on chromosome 7 (fl7.3 and fd7.3), with a PV of −19% in both instances. Closely located QTLs were found on chromosome 3 for both traits (fl3.4 and fd3.5) and on the bottom part of chromosome 12 for the fruit widths (fl12.3 and fd12.4). Overall, the reduction of the lengths and the widths of the fruits were lower than 20% for all the detected genomic regions.
In total, five QTLs for the fruit shape were detected; for three identified across all the trials, the SH alleles caused the fruit to be more flattened (decreasing of the FS). Among these, fs3.1 and fs11.1 had the highest and the lowest effects, respectively, with PV% ranging from −12.73% to −5.96%. On the bottom of chromosome 1, we found alleles leading to more elongated fruits, with increased shapes of 7.42% (FS1.5) and 7.28% (FS1.4) during seasons two and three, respectively.
The wild SH alleles were responsible for both decreasing and increasing the intensity of the fruit colors. Four QTLs for the internal and external colors of the fruits were found to colocalize at the bottom of chromosomes 1 and 2. In all instances, the wild alleles led to a reduction of color intensity, from red to orange/yellow. The strongest effects were observed on chromosome 2, with a PV of −53.95% for the internal color across the three years of evaluation. Minor effects were instead found on chromosome 1, with variations on growing seasons two and three for the internal (ic1.4) and external (ec1.4) colors, respectively. The QTLs involved in the increase of the fruit color intensity were detected across two seasons on the long arms of chromosomes 4, 9, and 11, with a PV ranging from 18.97 (EC11.3) to 25.05 (IC4.5). Although not colocalizing, neighboring QTLs were found at the base of chromosome 4 in a region spanning 96.75–119.0 cM.
Finally, the QTLs involved in pericarp thickness were identified on chromosomes 5 and 11, respectively. The former was detected across all seasons, increasing the pericarp thickness up to 32.49%, while the latter was detected only in two trials, leading to a reduction of −38.69%.
3.4.3. Chemical Traits
For the chemical traits, we found 11 associated regions in two or three growing seasons. In all cases, the SH alleles increased the soluble solids with a PV in the range of 24–27%. Most of the significant QTLs were detected during two seasons. Only SSC1.5, on the long arm of chromosome 1, had the greatest overall effect, being detected across all seasons. Despite the increase of the soluble solids, the product brix × red yield decreased in percentage, ranging from −57.73% to −39.33%. In tomatoes, this trait provides the estimation of the paste amount during processing. We also found four QTLs associated to a reduction of pH on the long arms of chromosomes 6 and 9 in two seasons (average PV ~4.5%) and on chromosomes 1 and 11 across three seasons, with an average PV of 8%.
3.4.4. Trichomes
The amount of trichomes is of interest, as it provides an estimate of the level of tolerance/resistance to virus-carrying insects such as aphids and thrips [42]. This protection mechanism is exercised both as a mechanical opposition and through the production of volatile repellent substances contained in the trichromatic glands [26]. Although there are very advanced methods of investigation based on microscopic techniques, the scale we used serves to give an estimate of the putative QTLs underlying the trait. Three highly significant QTLs were found on chromosomes 4, 7, and 10 and were related to TA1280, TA1304, and TA1551, respectively. In all cases, the wild alleles provided an increase of the trichomes number in comparison to the control in a percentage range of 100–200%.
4. Discussion
4.1. Performance of S. habrochaites Introgression Lines
Changes of the agricultural scenario in the coming decades will require efforts from breeders in finding novel sources of variations for the improvement of crops [7]. Wild tomato species have evolved in a wide range of habitats and, as a result of natural selection, represent a useful reservoir of genes for various agro-morphological, physio-biochemical, and resistances traits [2]. This exotic genome, once introgressed in the cultivated background, can be successfully exploited, breaking up the association with undesirable traits (linkage drag) and overcoming any occurring breeding barriers [43]. In the present study, a comprehensive approach aimed at assessing agronomic qualitative performances was carried out in multi-season trials on a set of 39 S. habrochaites LA1777 ILs. We firstly examined the performance of the ILs by focusing on the genetic and environmental factors affecting trait variations. Subsequently, we investigated those genomic regions underlying the studied traits and identified putative QTLs. A wide range of diversity was found among ILs for the considered traits with a medium-high level of heritability, which was maximized for fruit shape features. Despite the presence of genetically diverse lines and of additional introgressions for 26% of the lines, the fruit traits were minimally affected by environmental factors. In the set analyzed, we found trait plasticity depending on the chromosomal regions and type of traits, due likely to the different types of gene regions interacting with the environment.
Overall, the whole set of introgression lines resulted in being very promising for soluble solid contents rather than yield traits. The level of soluble solids measured as brix degrees is a measure of the total content of sugars and, in a minor extent, of the organic acids, together with small amounts of dissolved vitamins, fructans, proteins, pigments, phenolics, and minerals [44]. All these compounds influence the flavor and sweetness, which are important precursors of sensory quality [45].
Four lines at the bottom of chromosomes 1 (TA523 and TA1223) and 4 (TA1473 and TA1475) were previously used as a basis to fine map the agronomic traits [10,11]. Our results agree with earlier observations reporting an increase of soluble solids and a decrease of yield traits at the bottom of chromosome 1 in the overlapping region between TA523 and TA1223 [10]. Moreover, we confirmed a general decrease in fruit weight and an increase in the soluble solids content and color intensity at the bottom of chromosome 4 [11]. Regarding productivity, the differences in performances of the control genotype and the IL set in each cultivation year were almost stable. We observed a decrease of the yield and uniformity of maturity from the first to the last season. As expected, similar trends of correlations were found across the years, particularly within the same category of traits. Interestingly, significative negative correlations were found between the yield traits and soluble solids in the last growing season characterized by less rainfall, humidity, and average temperatures. It is likely that the combination of these effects led to a general decrease of TY and RY compared to the previous years in both E6203 and the ILs. Instead, the first season was characterized by the highest average temperatures and humidity, leading to the best yield values.
The soluble solids tended to be stable across the three seasons, although slightly higher in the less productive ones. Within the ILs, the contrasting values observed in the second year mainly concerned lines with multiple introgressions and/or low brix values in the other seasons, suggesting a major effect of the environment in these cases. Previous evidence in tomatoes reported how increasing the water stress led to the increase of the soluble solids content [46]. Moreover, the rainfall and temperature changes could significantly impact the overall performance of the tomatoes [47].
Although the used drip irrigation restored the entire water requirements at the root level, it is likely that the whole soil moisture was affected by the conditions that occurred in the last season, leading to the observed correlations. Furthermore, it is known that S. habrochaites is a cold-tolerant species; therefore, it is possible to suggest that lower temperatures may have a role in the increase of the soluble solids. Since no previous evidence was reported in S. habrochaites acc. LA1777, further investigation is required.
4.2. QTL Mapping
Despite S. habrochaites being recognized as an important source of alleles of agronomic interest, the potential of the existing exotic libraries has not yet been extensively exploited. We identified 75 regions likely with underlying trait variations, and the multi-evaluation trials allowed the identification of 33 robust QTLs, removing any possibility of false discoveries.
Our results agree with previous studies involving diverse types of exotic libraries and experimental mapping populations developed from wild relatives [8,29,30,48,49,50]. Many QTLs reducing the yield performances and fruit weights confirm the preponderant effect of the wild Solanum alleles on these traits. We found specific regions associated with an increase in plant weights, soluble solids contents, internal and external color intensities, pericarp thickness, and trichome density. The results of this study confirmed that S. habrochaites could be a potential source of favorable alleles for horticultural traits, despite the inferior phenotype compared to the cultivated parent. Some genomic regions seem to be conserved across wild tomato relatives; as an example, the increase of soluble solids at the bottom of chromosome 1 fall in a region previously reported carrying alleles improving the trait in other wild species, including S. pennellii, S. neorickii, and S. peruvianum [8,48,49,50]. Furthermore, the additional QTLs detected on chromosomes 3 (SSC3.6) and 10 (SSC10.1) were also reported in S. pennellii introgression lines [8] and, for the one on chromosome 10, in advanced backcross lines of S. peruvianum [48,50]. These results confirm the existence of shared regions in various tomato wild species carrying alleles for improving the quality of cultivated ones. Clusters of QTLs involved in fruit weight and morphology were found at the long arms of chromosomes 3, 7, 11, and 12, suggesting a gene linkage. For all, the underlying alleles led to a reduction of fruit weights and shapes in agreement with the wild phenotype, which has an average weight fifty-fold less than the recurrent parent and rounder fruits. Fruit color is one of the main attributes of quality and an indicator of the contents of health-related compounds. Diverse orange, yellow, pink, and black tomatoes are becoming increasingly popular by consumers who appreciate a rich and varied diet. The decrease of internal and external fruit color at the bottom of chromosome 2 was associated with fw2.2, a major gene modulating fruit sizes [33]. The same gene was reported being associated with increased color in S. pimpinellifolium due to a pleiotropic effect between the decreasing fruit size and pigment concentration [51]. Pleiotropy between the fruit size and color occurs also in SH, although with the opposite effect. In S. pennellii, a decrease in color was shown by the mutation of yellow flesh (r) on chromosome 3, responsible for the loss of function of the phytoene synthase, a key enzyme of the carotenoid biosynthetic pathway [52]. The mutation was carried by IL3-2, which defines the same region of TA1539 and TA1541. However, these two lines had additional introgressions, which likely masked a possible effect of any occurring mutation.
Therefore, the additional alleles part of the SH genomic background could be involved in fruit pigmentation. A strong QTL responsible for the increase of trichome density was found on chromosome 10. Recently, Barrantes et al. [15] identified a major QTL in the same chromosomal region in S. pimpinellifolium ILs. In agreement with our study, the trait was scored by the same scale. Trichomes are associated with insect resistance, although no candidate gene has been reported on chromosome 10. These results suggest that further studies are needed for a better investigation of the underlying region, although it cannot be excluded that this increase in the density of trichomes does not correspond to a resistance.
4.3. Future Prospects for Using S. habrochaites ILs in Tomato Breeding
Beyond QTL identification, the advantage of using introgression lines relies on the possibility of their direct use in breeding programs as parent lines. This offers an invaluable opportunity for transferring new traits for the improvement of modern varieties. The set of SH ILs analyzed in this study resulted in being promising for the diverse phenotyped traits. Furthermore, the number of QTLs detected in specific lines facilitates gene pyramiding; as an example, the improvement of soluble solids and pericarp thickness could be obtained by crossing TA1287 either with TA1223 or TA523, which, in addition, showed a slight increase of pH. The synergy of these QTLs could improve the processing yields, with less acidic fruits with greater pulp consistency. For a successful introgression program, it is very important to choose phenotypes stable across a wide range of environmental conditions. To that end, QTLs with large genetic effects are preferred. However, it must be considered, the occurrence of constraints such as pleiotropy, linkage-drag, epistatic effects, and relationships with other beneficial traits. For example, introgressing soluble solids is often challenging due to inverse correlations with the yield [53]. Moreover, pleiotropic effects may occur between the fruit size and color [51].
Despite one of the advantages in holding only a fraction of the entire wild genome being the reduction of fertility problems [8], it must be recognized that several S. habrochaites ILs still carry multiple wild segments, as well as large introgressed fragments, covering, in some instances, more than half the chromosome. Therefore, the linkage between the beneficial genes introduced during the backcrossing with the deleterious ones (linkage drag) could reduce the potentiality of ILs for both breeding and QTL mapping purposes. In the current set, we demonstrated how the availability of lines with reduced introgressions led to a major accuracy in the identification and positioning of QTLs with positive effects (e.g., SSC1.5 and PER5.1). While linkage drag occurs for large introgression lines such as TA1304, with both the increasing of the trichome density and reduction of the fruit weight and size.
Our results suggest, in some instances, the necessity of generating new sub-ILs able to break up the linkage occurring with undesired traits and facilitate the fine mapping of QTLs. To that end, genomic platforms such as SNP (single nucleotide polymorphism) arrays or reduced-representation sequencing (e.g., ddRAD-seq: Double digest restriction-site associated DNA and GBS: Genotyping by sequencing) could serve as a key step for the generation of high-density marker data, speeding up the deep exploitation of the current S. habrochaites exotic library. These tools have been successfully applied in tomatoes and other related species, providing large-scale and cost-effective genome scans and generating novel knowledge for genomic-assisted breeding [54,55].
5. Conclusions
The present study aims to provide a broad overview of S. habrochaites introgression lines in relation to the morpho-agronomic performance and chemical properties. Through the investigation of phenotypic performances, the G × Y interaction, and the QTL analysis, we provided novel information for genetic and breeding purposes. This study represents a step toward the exploitation of the hidden properties of S. habrochaites, to be implemented through the development of near-isogenic lines assisted by genomics approaches.
Supplementary Materials
The following are available online at https://www.mdpi.com/2073-4395/11/1/38/s1: Table S1: Parameters describing the trait performances of the recurrent parent “E6203” and the 39 SH ILs for each growing season. For each trait, the significant differences between the ILs and the recurrent parent is reported. Means in the same column with different letters are significantly different, according to Tukey’s honest significant difference test; Table S2: Mean values for the introgression lines and the control genotype in each growing season (Y1, Y2, and Y3). In bold and with asterisks are indicated those means that significantly differ compared with the control, according to Dunnett’s test; Table S3: Percent of phenotypic variation with respect to the control for the introgression lines in each growing season (Y1, Y2, and Y3). In bold are indicated those values corresponding to significant differences according to Dunnett’s test; Table S4: The phenotypic plasticity index (PPI) for the traits analyzed in 39 S. habrochaites ILs and the control S. lycopersicum E6203. PPI values > 1.0 are in bold and gray highlighted. For each IL is indicated the chromosomes holding introgressions; Table S5: Variable contribution (VarPC), correlation coefficient (CorrPC), and Eigenvectors (Egv) for fruit descriptors in the two first principal components. Highlighted traits explained a variation >12% in the first or second PC.
Author Contributions
Conceptualization, P.T.; methodology, P.T., M.P., and B.D.; validation, P.T.; formal analysis, P.T.; data curation, P.T., M.C., and M.P.; writing—original draft preparation, P.T.; writing—review and editing, P.T., M.P., A.V., and M.C.; and supervision, P.T. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by funded by the Italian Ministry of Agriculture, Food and Forestry, grant name “RGV-FAO”. Part of this research is carried out in the context of a PhD funded by CREA—Council for Agricultural Research and Economics.
Acknowledgments
The authors wish to acknowledge Silvana Grandillo of CNR-IBBR, Luigi Frusciante of the University of Naples, and Massimo Zaccardelli of CREA-OF. We furthermore thank Armida Del Galdo, Rosaria Macellaro, Luigi Aliberti, and Giovanni Ragosta for the valuable support.
Conflicts of Interest
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.
References
- FAOSTAT. Available online: http://www.fao.org/faostat/ (accessed on 15 November 2020).
- Grandillo, S.; Chetelat, R.; Knapp, S.; Spooner, D.; Peralta, I.; Cammareri, M.; Perez, O.; Termolino, P.; Tripodi, P.; Chiusano, M.; et al. Solanum sect. Lycopersicon. In Wild Crop Relatives and Breeding Resources: Vegetables; Kole, C., Ed.; Springer: Berlin/Heidelberg, Germany, 2011; pp. 129–215. [Google Scholar]
- Foolad, M.R.; Panthee, D.R. Marker-assisted selection in tomato breeding. Crit. Rev. Plant Sci. 2012, 31, 93–123. [Google Scholar] [CrossRef]
- Zamir, D. Improving plant breeding with exotic genetic libraries. Nat. Rev. Genet. 2001, 2, 983–989. [Google Scholar] [CrossRef] [PubMed]
- Tripodi, P.; Ficcadenti, N.; Rotino, G.L.; Festa, G.; Bertone, A.; Pepe, A.; Caramanico, R.; Migliori, C.A.; Spadafora, D.; Schiavi, M.; et al. Genotypic and environmental effects on the agronomic, health-related compounds and antioxidant properties of chilli peppers for diverse market destinations. J. Sci. Food Agric. 2019, 99, 4550–4560. [Google Scholar] [CrossRef] [PubMed]
- Gur, A.; Zamir, D. Unused natural variation can lift yield barriers in plant breeding. PLoS Biol. 2004, 2, 1610–1615. [Google Scholar] [CrossRef]
- D’Agostino, N.; Tripodi, P. NGS-based genotyping, high throughput phenotyping and genome-wide association studies laid the foundations for next-generation breeding in horticultural crops. Diversity 2017, 9, 38. [Google Scholar] [CrossRef] [Green Version]
- Eshed, Y.; Zamir, D. An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield-associated QTL. Genetics 1995, 141, 1147–1162. [Google Scholar]
- Monforte, A.J.; Tanksley, S.D. Development of a set of near isogenic and backcross recombinant inbred lines containing most of the Lycopersicon hirsutum genome in a L. esculentum genetic background: A tool for gene mapping and gene discovery. Genome 2000, 43, 803–813. [Google Scholar] [CrossRef]
- Monforte, A.J.; Tanksley, S.D. Fine mapping of a quantitative trait locus (QTL) from Lycopersicon hirsutum chromosome 1 affecting fruit characteristics and agronomic traits: Breaking linkage among QTLs affecting different traits and dissection of heterosis for yield. Theor. Appl. Genet. 2000, 100, 471–479. [Google Scholar] [CrossRef]
- Monforte, A.J.; Friedman, E.; Zamir, D.; Tanksley, S.D. Comparison of a set of allelic QTL-NILs for chromosome 4 of tomato: Deductions about natural variation and implications for germplasm utilization. Theor. Appl. Genet. 2001, 102, 572–590. [Google Scholar] [CrossRef]
- Finkers, R.; van Heusden, A.W.; Meijer-Dekens, F.; van Kan, J.A.L.; Maris, P.; Lindhout, P. The construction of a Solanum habrochaites LYC4 introgression line population and the identification of QTLs for resistance to Botrytis cinerea. Theor. Appl. Genet. 2007, 114, 1071–1080. [Google Scholar] [CrossRef] [Green Version]
- Chetelat, R.T.; Meglic, V. Molecular mapping of chromosome segments introgressed from Solanum lycopersicoides into cultivated tomato (Lycopersicon esculentum). Theor. Appl. Genet. 2000, 100, 232–241. [Google Scholar] [CrossRef]
- Canady, M.A.; Meglic, V.; Chetelat, R.T. A library of Solanum lycopersicoides introgression lines in cultivated tomato. Genome 2005, 48, 685–697. [Google Scholar] [CrossRef] [PubMed]
- Barrantes, W.; López-Casado, G.; García-Martínez, S.; Alonso, A.; Rubio, F.; Ruiz, J.J.; Fernández-Muñoz, R.; Granell, A.; Monforte, A.J. Exploring New Alleles Involved in Tomato Fruit Quality in an Introgression Line Library of Solanum pimpinellifolium. Front. Plant Sci. 2016, 7, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Chetelat, R.T.; Qin, X.; Tan, M.; Burkart-Waco, D.; Moritama, Y.; Huo, X.; Wills, T.; Pertuzé, R. Introgression lines of Solanum sitiens; a wild nightshade of the Atacama Desert; in the genome of cultivated tomato. Plant J. 2019, 100, 836–850. [Google Scholar] [CrossRef]
- Stevens, R.; Buret, M.; Duffé, P.; Garchery, C.; Baldet, P.; Rothan, C.; Causse, M. Candidate genes and Quantitative Trait Loci affecting fruit ascorbic acid content in three tomato populations. Plant Physiol. 2007, 143, 1943–1953. [Google Scholar] [CrossRef] [Green Version]
- Goodstal, J.F.; Kohler, G.; Randall, L.; Bloom, A.; Clair, D.S. A major QTL introgressed from wild Lycopersicon hirsutum confers chilling tolerance to cultivated tomato (Lycopersicon esculentum). Theor. Appl. Genet. 2005, 111, 898–905. [Google Scholar] [CrossRef]
- Easlon, H.M.; Asensio, J.S.; St Clair, D.A.; Bloom, A.J. Chilling-induced water stress: Variation in shoot turgor maintenance among wild tomato species from diverse habitats. Am. J. Bot. 2013, 100, 1991–1999. [Google Scholar] [CrossRef]
- Vidavski, F.; Czosnek, H.; Gazit, S.; Levy, D.; Lapidot, M. Pyramiding of genes conferring resistance to Tomato yellow leaf curl virus from different wild tomato species. Plant Breed. 2008, 127, 625–631. [Google Scholar] [CrossRef]
- Sifres, A.; Blanca, J.; Nuez, F. Pattern of genetic variability of Solanum habrochaites in its natural area of distribution. Genet. Resour. Crop Evol. 2011, 58, 347–360. [Google Scholar] [CrossRef]
- Momotaz, A.; Scott, J.W.; Schuster, D.J. Searching for silverleaf whitefly and begomovirus resistance genes from Lycopersicon hirsutum accession LA1777. Acta Hortic. (ISHS) 2004, 695, 417–422. [Google Scholar] [CrossRef]
- Thapa, S.P.; Miyao, E.M.; Davis, R.M.; Coaker, G. Identification of QTLs controlling resistance to Pseudomonas syringae pv. tomato race 1 strains from the wild tomato; Solanum habrochaites LA1777. Theor. Appl. Genet. 2015, 128, 681–692. [Google Scholar] [CrossRef] [PubMed]
- Bao, Z.; Meng, F.; Strickler, S.R.; Dunham, D.M.; Munkvold, K.R.; Martin, G.B. Identification of a candidate gene in Solanum habrochaites for resistance to a race 1 strain of Pseudomonas syringae pv. tomato. Plant Genome 2015, 8, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Avila, C.A.; Marconi, T.G.; Viloria, Z.; Kurpis, J.; Del Rio, S.Y. Bactericera cockerelli resistance in the wild tomato Solanum habrochaites is polygenic and influenced by the presence of Candidatus Liberibacter solanacearum. Sci. Rep. 2019, 9, 14031. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Kim, J.; Kang, K.; Gonzales-Vigil, E.; Shi, F.; Jones, A.D.; Barry, C.S.; Last, R.L. 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, 1854–1870. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Zabel, S.; Brandt, W.; Porzel, A.; Athmer, B.; Kortbeek, R.; Bleeker, P.; Tissier, A. Two novel 7-epi-zingiberene derivatives with biological activity from Solanum habrochaites are produced by a single cytochrome P450 monooxygenase. Biorxiv 2020. [Google Scholar] [CrossRef]
- Dal Cin, V.; Kevany, B.; Fei, Z.; Klee, H.J. Identification of Solanum habrochaites loci that quantitatively influence tomato fruit ripening-associated ethylene emissions. Theor. Appl. Genet. 2009, 119, 1183–1192. [Google Scholar] [CrossRef]
- Bernacchi, D.; Beck-Bunn, T.; Emmatty, D.; Eshed, Y.; Inai, S.; Lopez, J.; Petiard, V.; Sayama, H.; Uhlig, J.; Zamir, D.; et al. Advanced backcross QTL analysis of tomato. II. Evaluation of near-isogenic lines carrying single-donor introgressions for desirable wild QTL-alleles derived from Lycopersicon hirsutum and L. pimpinellifolium. Theor. Appl. Genet. 1998, 97, 170–180. [Google Scholar] [CrossRef]
- Bernacchi, D.; Beck-Bunn, T.; Eshed, Y.; Lopez, J.; Petiard, V.; Uhlig, J.; Zamir, D.; Tanksley, S.D. Advanced backcross QTL analysis in tomato I. Identification of QTLs for traits of agronomic importance from Lycopersicon hirsutum. Theor. Appl. Genet. 1998, 97, 381–397. [Google Scholar] [CrossRef]
- Alseekh, S.; Ofner, I.; Pleban, T.; Tripodi, P.; Di Dato, F.; Cammareri, M.; Mohammad, A.; Grandillo, S.; Fernie, A.R.; Zamir, D. Resolution by recombination: Breaking up Solanum pennellii introgressions. Trends Plant Sci. 2013, 18, 536–538. [Google Scholar] [CrossRef]
- SOLGENOMICS. Available online: http://www.solegenomics.net (accessed on 15 November 2020).
- Frary, A.; Nesbitt, T.C.; Grandillo, S.; Knaap, E.; Cong, B.; Liu, J.; Meller, J.; Elber, R.; Alpert, K.B.; Tanksley, S.D. fw2.2: A quantitative trait locus key to the evolution of tomato fruit size. Science 2000, 289, 85–88. [Google Scholar] [CrossRef] [Green Version]
- Fridman, E.; Carrari, F.; Liu, Y.S.; Fernie, A.R.; Zamir, D. Zooming in on a quantitative trait for tomato yield using interspecific introgressions. Science 2004, 305, 1786–1789. [Google Scholar] [CrossRef] [PubMed]
- Hanson, P.M.; Sitathanim, K.; Sadashivam, A.T.; Yang, R.Y.; Grahamm, E.; Ledesma, D. Performance of Solanum habrochaites LA1777 introgression line hybrids for marketable tomato fruit yield in Asia. Euphytica 2007, 158, 167–178. [Google Scholar] [CrossRef]
- Tomato Genetic Resource Center. Available online: https//tgrc.ucdavis.edu (accessed on 15 November 2020).
- Ronga, D.; Parisi, M.; Pentangelo, A.; Mori, M.; Di Mola, I. Effects of nitrogen management on biomass production and dry matter distribution of processing tomato cropped in Southern Italy. Agronomy 2019, 9, 855. [Google Scholar] [CrossRef] [Green Version]
- Parisi, M.; Di Dato, F.; Ricci, S.; Mennella, G.; Cardi, T.; Tripodi, P. A multi-trait characterization of the ‘Friariello’ landrace: A Mediterranean resource for sweet pepper breeding. Plant Genet. Resour. 2017, 15, 165–176. [Google Scholar] [CrossRef]
- Dingemanse, N.J.; Kazem, A.J.M.; Réalem, D.; Wright, J. Behavioural reaction norms, animal personality meets individual plasticity. Trends Ecol. Evol. 2010, 25, 81–89. [Google Scholar] [CrossRef] [PubMed]
- Mantel, N.A. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar] [PubMed]
- Dunnett, C.W. A multiple comparison procedure for comparing several treatments with a control. J. Am. Stat. Assoc. 1955, 50, 1096–1121. [Google Scholar] [CrossRef]
- Simmons, A.T.; Gurr, G.M. Trichomes of Lycopersicon species and their hybrids, effects on pests and natural enemies. Agric. For. Entomol. 2005, 7, 265–276. [Google Scholar] [CrossRef]
- Prohens, J.; Gramazio, P.; Plazas, M.; Dempewolf, H.; Kilian, B.; Díez, M.J.; Fita, A.; Herraiz, F.J.; Rodríguez-Burruezo, B.; Soler, S.; et al. Introgressiomics, a new approach for using crop wild relatives in breeding for adaptation to climate change. Euphytica 2017, 213, 158. [Google Scholar] [CrossRef]
- Magwaza, L.S.; Opara, U.L. Analytical methods for determination of sugars and sweetness of horticultural products—A review. Sci. Hortic. 2015, 184, 179–192. [Google Scholar] [CrossRef]
- Gautier, H.; Diakou-Verdin, V.; Bénard, C.; Reich, M.; Buret, M.; Bourgaud, F.; Poëssel, J.L.; Caris-Veyrat, C.; Génard, M. How does tomato quality (sugar.; acid.; and nutritional quality) vary with ripening stage.; temperature.; and irradiance? J. Agric. Food Chem. 2008, 56, 1241–1250. [Google Scholar] [CrossRef] [PubMed]
- Zhang, H.; Xiong, Y.; Huang, G.; Xu, X.; Huang, Q. Effects of water stress on processing tomatoes yield, quality and water use efficiency with plastic mulched drip irrigation in sandy soil of the Hetao Irrigation District. Agric. Water Manag. 2017, 179, 205–214. [Google Scholar] [CrossRef]
- Cammarano, D.; Ronga, D.; Di Mola, I.; Mori, M.; Parisi, M. Impact of climate change on water and nitrogen use efficiencies of processing tomato cultivated in Italy. Agric. Water Manag. 2020, 241, 106336. [Google Scholar] [CrossRef]
- Fulton, T.M.; Beck-Bunn, T.; Emmatty, D.; Eshed, Y.; Lopez, J.; Petiard, V.; Uhlig, J.; Zamir, D.; Tanksley, S.D. QTL analysis of an advanced backcross of Lycopersicon peruvianum to the cultivated tomato and comparisons with QTLs found in other wild species. Theor. Appl. Genet. 1997, 95, 881–894. [Google Scholar] [CrossRef]
- Fulton, T.M.; Grandillo, S.; Beck-Bunn, T.; Fridman, E.; Frampton, A.; Lopez, J.; Petiard, V.; Uhlig, J.; Zamir, D.; Tanksley, S.D. Advanced backcross QTL analysis of a Lycopersicon esculentum × Lycopersicon parviflorum cross. Theor. Appl. Genet. 2000, 100, 1025–1042. [Google Scholar] [CrossRef]
- Fulton, T.M.; Buchelim, E.; Voirolm, E.; Lopez, J.; Pétiard, V.; Tanksley, S.D. Quantitative trait loci (QTL) affecting sugars.; organic acids and other biochemical properties possibly contributing to flavor.; identified in four advanced backcross populations of tomato. Euphytica 2002, 127, 163–177. [Google Scholar] [CrossRef]
- Tanksley, S.D.; Grandillo, S.; Fulton, T.M.; Zamir, D.; Eshed, Y.; Petiard, V.; López, J.; Beck-Bunn, T. Advanced backcross QTL analysis in a cross between an elite processing line of tomato and its wild relative L. pimpinellifolium. Theor. Appl. Genet. 1996, 92, 213–224. [Google Scholar] [CrossRef]
- Liu, J.; Gur, A.; Ronen, G.; Causse, M.; Damidaux, R.; Buret, M.; Hirschberg, J.; Zamir, D. There is more to tomato fruit colour than candidate carotenoid genes. Plant Biotechnol. J. 2003, 1, 195–207. [Google Scholar] [CrossRef] [Green Version]
- Liabeuf, D.; Francis, D.M. The use of historical datasets to develop multi-trait selection models in processing tomato. Euphytica 2017, 213, 100. [Google Scholar] [CrossRef] [Green Version]
- Esposito, S.; Cardi, T.; Campanelli, G.; Sestili, S.; Díez, M.J.; Soler, S.; Prohens, J.; Tripodi, P. ddRAD sequencing-based genotyping for population structure analysis in cultivated tomato provides new insights into the genomic diversity of Mediterranean ‘da serbo’ type long shelf-life germplasm. Hortic. Res. 2020, 7, 134. [Google Scholar] [CrossRef]
- Colonna, V.; D’Agostino, N.; Garrison, E.; Albrechtsen, A.; Meisner, J.M.; Facchiano, A.; Cardi, T.; Tripodi, P. Genomic diversity and novel genome-wide association with fruit morphology in Capsicum, from 746k polymorphic sites. Sci. Rep. 2019, 9, 10067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1.
Phenotypic differences (in percentage with respect to the SL recurrent parent as the control) for the soluble solids content in 39 Solanum habrochaites introgression lines tested in three growing seasons. Blue and orange-colored bars with black edges indicate the positive and negative significant differences at p < 0.05, respectively.
Figure 1.
Phenotypic differences (in percentage with respect to the SL recurrent parent as the control) for the soluble solids content in 39 Solanum habrochaites introgression lines tested in three growing seasons. Blue and orange-colored bars with black edges indicate the positive and negative significant differences at p < 0.05, respectively.
Figure 2.
Phenotypic variability in the studied introgression lines and the recurrent parent. Loading plot of the first (PC1) and second (PC2) principal components showed the variations for the 19 morphological and agronomic traits in the set of 39 introgressions lines and the control (cv E6203) Introgressions lines (ILs) carrying single introgressions are indicated in green font with a triangle symbol, whereas ILs carrying multiple fragments are indicated in blue font with a circular symbol. The control is indicated with a square in red font. The direction and distance from the center of the biplot indicate how each trait contributes to the first two components. The explanation for trait acronyms is reported in Table 3.
Figure 2.
Phenotypic variability in the studied introgression lines and the recurrent parent. Loading plot of the first (PC1) and second (PC2) principal components showed the variations for the 19 morphological and agronomic traits in the set of 39 introgressions lines and the control (cv E6203) Introgressions lines (ILs) carrying single introgressions are indicated in green font with a triangle symbol, whereas ILs carrying multiple fragments are indicated in blue font with a circular symbol. The control is indicated with a square in red font. The direction and distance from the center of the biplot indicate how each trait contributes to the first two components. The explanation for trait acronyms is reported in Table 3.
Figure 3.
Pearson’s rank correlation coefficients between the pairs of phenotypes. The correlation coefficients are indicated in each cell. Colored correlations are those with p-values < 0.01 after Bonferroni correction. The color intensity is directly proportional to the coefficients. On the right side of the correlogram, the color legend shows the correlation coefficients and the corresponding colors. (a) The correlogram for the traits scored in year 1, (b) year 2, and (c) year 3. The trait acronyms are in Table 3.
Figure 3.
Pearson’s rank correlation coefficients between the pairs of phenotypes. The correlation coefficients are indicated in each cell. Colored correlations are those with p-values < 0.01 after Bonferroni correction. The color intensity is directly proportional to the coefficients. On the right side of the correlogram, the color legend shows the correlation coefficients and the corresponding colors. (a) The correlogram for the traits scored in year 1, (b) year 2, and (c) year 3. The trait acronyms are in Table 3.
Figure 4.
Genetic map showing the position of the markers delimiting the ILs and used for the quantitative trait loci (QTL) analysis. Markers are indicated on the right side of each chromosome, whereas map distances in centiMorgan (cM) are indicated on the left side and are based on the Tomato-EXPEN 2000 and Tomato-EXPEN 1992 (with asterisks) mapping populations [32]. At the bottom of each chromosome, the total genetic length is mentioned. Black bars indicated ILs carrying a single homozygous introgression. Dashed lines indicate ILs with multiple introgressions. The empty part of each chromosome indicates the regions not covered by wild introgressions. Centromeres are indicated with gray dots. Light blue bars on the left of each chromosome indicate the genomic positions of the markers based on tomato reference annotation SL4.0. The total length in megabases is reported at the bottom. Only robust QTLs detected across the three consecutive seasons are reported. The QTLs are in italics and defined by the trait abbreviation, followed by the chromosomal BIN (region delimiting single or overlapping introgressions). In uppercase and blue font are the QTLs increasing the trait performances. In red and in lowercase are the indicated QTLs decreasing the traits.
Figure 4.
Genetic map showing the position of the markers delimiting the ILs and used for the quantitative trait loci (QTL) analysis. Markers are indicated on the right side of each chromosome, whereas map distances in centiMorgan (cM) are indicated on the left side and are based on the Tomato-EXPEN 2000 and Tomato-EXPEN 1992 (with asterisks) mapping populations [32]. At the bottom of each chromosome, the total genetic length is mentioned. Black bars indicated ILs carrying a single homozygous introgression. Dashed lines indicate ILs with multiple introgressions. The empty part of each chromosome indicates the regions not covered by wild introgressions. Centromeres are indicated with gray dots. Light blue bars on the left of each chromosome indicate the genomic positions of the markers based on tomato reference annotation SL4.0. The total length in megabases is reported at the bottom. Only robust QTLs detected across the three consecutive seasons are reported. The QTLs are in italics and defined by the trait abbreviation, followed by the chromosomal BIN (region delimiting single or overlapping introgressions). In uppercase and blue font are the QTLs increasing the trait performances. In red and in lowercase are the indicated QTLs decreasing the traits.
Table 1.
Chromosomal location, number, size, and percentage of the coverage of the wild fragments in the Solanum habrochaites LA1777 introgression lines phenotyped in the present study.
Table 1.
Chromosomal location, number, size, and percentage of the coverage of the wild fragments in the Solanum habrochaites LA1777 introgression lines phenotyped in the present study.
| LA Pedigree TGRC a | TA IL Name b | Chromosomes Bearing Wild Segments c | N° of Donor Segments c | Introgression Size (cM) d | Total Percentage of Coverage e |
|---|---|---|---|---|---|
| LA3921 | TA1105 | 2 | 1 | 118.50 | 8.11 |
| LA3981 | TA1116 | 5 | 1 | 47.50 | 3.25 |
| LA3969 | TA1121 | 12 | 1 | 87.25 | 5.97 |
| LA3919 | TA1128 | 1 | 1 | 41.00 | 2.81 |
| LA3916 | TA1223 | 1 | 1 | 59.50 | 4.07 |
| LA3913 | TA1258 | 1 | 1 | 103.00 | 7.05 |
| LA3922 | TA1266 | 2 | 1 | 61.45 | 4.21 |
| LA3926 | TA1276 | 3 | 1 | 88.00 | 6.02 |
| LA3931 | TA1280 | 4 | 1 | 67.30 | 4.61 |
| LA3938 | TA1287 | 5 | 1 | 20.00 | 1.37 |
| LA3939 | TA1293 | 5 | 1 | 47.50 | 3.25 |
| LA3948 | TA1303 | 7 | 1 | 62.90 | 4.31 |
| LA3949 | TA1304 | 7 | 1 | 113.00 | 7.73 |
| LA3951 | TA1312 | 7 | 1 | 22.90 | 1.57 |
| LA3953 | TA1316 | 8 | 1 | 46.25 | 3.17 |
| LA3955 | TA1320 | 8 | 1 | 9.50 | 0.65 |
| LA3956 | TA1324 | 9 | 1 | 61.10 | 4.18 |
| LA3991 | TA1326 | 9 | 1 | 46.85 | 3.21 |
| LA3958 | TA1330 | 9, 11 | 1 | 45.90, 57.50 | 7.08 |
| LA3936 | TA1475 | 4 | 1 | 34.35 | 2.35 |
| LA3937 | TA1473 | 4 | 1 | 33.50 | 2.29 |
| LA3917 | TA1535 | 1, 2, 12 | 4 | 17.25, 21.00, 28.25 | 4.55 |
| LA3920 | TA1536 | 1 | 1 | 10.00 | 0.68 |
| LA3923 | TA1537 | 2 | 1 | 21.00 | 1.44 |
| LA3944 | TA1539 | 3, 6 | 2 | 39.20, 21.30 | 4.14 |
| LA3929 | TA1541 | 3, 8 | 2 | 42.00, 9.50 | 3.52 |
| LA3933 | TA1542 | 4 | 1 | 13.75 | 0.94 |
| LA3941 | TA1543 | 5 | 1 | 19.50 | 1.33 |
| LA3945 | TA1545 | 6, 10 | 2 | 50.95, 7.50 | 4.00 |
| LA3954 | TA1548 | 8 | 1 | 66.50 | 4.55 |
| LA3961 | TA1551 | 10 | 1 | 78.50 | 5.37 |
| LA3995 | TA1553 | 1, 11, 12 | 3 | 103.00, 24.50, 12.00 | 9.55 |
| LA3966 | TA1554 | 10, 11, 12 | 3 | 29.00, 31.00, 41.50 | 6.95 |
| LA3965 | TA1555 | 2, 10, 11 | 3 | 21.50, 29.00, 9.50 | 4.11 |
| LA3947 | TA1559 | 6 | 1 | 21.75 | 1.49 |
| LA4005 | TA1562 | 4 | 1 | 26.25 | 1.80 |
| LA4002 | TA1645 | 1, 8, 12 | 3 | 11.50, 9.50, 30.75 | 3.54 |
| LA4004 | TA1649 | 2, 3, 6 | 3 | 21.00, 56.55, 5.00 | 5.65 |
| LA3914 | TA523 | 1 | 1 | 50.25 | 3.44 |
a TGRC = Tomato Genetic Resource Center; LA is the accession prefix used by the TGRC germplasm bank. b TA is the name of each introgression line (IL) as reported by Monforte and Tanksley [9]. c Information reported in Monforte et al. (2000) and the Tomato Genetic Resource Center. d determined considering the chromosomal region covered by each introgression, plus the half-intervals between the flanking markers and nearest mapped loci. e Considering the tomato genome size equal to 1461 cM.
Table 2.
Geographic coordinates, pedological characteristics, and climatic conditions of the cultivation environments of the present study. For temperature, precipitation, and humidity, the mean values of the grown season of each year (April–August) are reported.
Table 2.
Geographic coordinates, pedological characteristics, and climatic conditions of the cultivation environments of the present study. For temperature, precipitation, and humidity, the mean values of the grown season of each year (April–August) are reported.
| Location | Battipaglia | |||
| Geographical coordinates | 40°37′ N; 14°58′ E | |||
| M.a.s.l. * | 65 | |||
| Soil regions + | Cambisol, Regosol, Calcisol, Phaeozem, Luvisol | |||
| Year | Y1 (2014) | Y2 (2015) | Y3 (2016) | Average 50 Years # |
| Maximum temperature (°C) | 28.10 | 27.67 | 26.66 | 21.10 |
| Minimum temperature (°C) | 17.50 | 17.14 | 16.02 | 11.50 |
| Precipitation (mm) | 12.60 | 14.17 | 8.53 | 16.50 |
| Humidity (%) | 67.30 | 65.76 | 63.98 | 70.64 |
* Meters above sea level. + http://esdac.jrc.ec.europa.eu/content/carta-dei-suoli-ditalia-soil-map-italy (accessed on 26 December 2020); in brackets is the soil type number. # Average values of 1966–2016.
Table 3.
Descriptive statistics of the traits analyzed in the 39 introgression lines (IL) and the control (SL) during three consecutive growing seasons.
Table 3.
Descriptive statistics of the traits analyzed in the 39 introgression lines (IL) and the control (SL) during three consecutive growing seasons.
| Traits a | # | MS b | R b | F Value c | H2 b | Mean SL | Range SL | CV b | Mean IL | Range IL | CV b | Prob > F |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Yield | ||||||||||||
| Total yield | TY | 12.13 | 0.22 | 8.46 ** | 0.89 | 4.11 | 4.7–3.3 | 31.69 | 3.41 | 4.16–2.76 | 38.58 | *** |
| Red yield | RY | 13.53 | 0.22 | 8.49 ** | 0.89 | 3.77 | 4.67–2.96 | 36.53 | 3.16 | 4.06–2.32 | 44.19 | *** |
| Green yield | GY | 1.47 | 0.18 | 5.36 ** | 0.84 | 0.50 | 0.82–0.03 | 95.77 | 0.56 | 0.99–0.24 | 102.67 | ns |
| Uniformity maturity | UM | 974.95 | 0.21 | 8.33 ** | 0.89 | 91.01 | 99.26–80.45 | 11.55 | 90.98 | 96.54–83.43 | 13.38 | ns |
| Plant weight | PW | 0.77 | 0.19 | 7.01 ** | 0.88 | 0.67 | 0.86–0.45 | 68.56 | 0.68 | 0.75–0.55 | 51.43 | ns |
| Biomass | BM | 13.68 | 0.19 | 7.37 ** | 0.88 | 4.79 | 5.45–3.75 | 33.43 | 4.08 | 4.89–3.31 | 35.88 | *** |
| Harvest index | HI | 0.06 | 0.27 | 11.50 ** | 0.92 | 0.86 | 0.88–0.84 | 6.93 | 0.83 | 0.85–0.81 | 10.66 | *** |
| Fruit | ||||||||||||
| Fruit weight | FW | 2029.75 | 0.25 | 10.28 ** | 0.91 | 74.27 | 83.52–67.85 | 20.16 | 60.06 | 62.96–57.08 | 25.76 | *** |
| Fruit length | FL | 3.05 | 0.07 | 21.02 ** | 0.95 | 5.68 | 6.01–5.52 | 7.17 | 5.09 | 5.13–5.04 | 9.04 | *** |
| Fruit diameter | FD | 1.52 | 0.24 | 15.92 ** | 0.94 | 4.94 | 5.23–4.79 | 7.37 | 4.52 | 4.57–4.45 | 7.83 | *** |
| Fruit shape | FS | 0.07 | 0.41 | 21.51 ** | 0.96 | 1.15 | 1.15–1.15 | 4.29 | 1.13 | 1.13–1.12 | 6.62 | ** |
| External color | CE | 2.57 | 0.49 | 18.77 ** | 0.94 | 2.92 | 3.75–2.44 | 20.61 | 2.90 | 3.27–2.67 | 19.23 | ns |
| Internal color | CI | 1.69 | 0.60 | 13.73 ** | 0.84 | 3.11 | 3.81–2.73 | 17.04 | 3.09 | 3.36–2.93 | 19.18 | ns |
| Pericarp | PER | 1.51 | 0.22 | 8.41 ** | 0.89 | 2.07 | 2.22–1.67 | 22.76 | 1.90 | 1.93–1.87 | 24.56 | *** |
| Puffiness | PUF | 2.52 | 0.27 | 11.20 ** | 0.92 | 1.68 | 2.16–1.24 | 35.18 | 1.57 | 1.76–1.33 | 34.48 | * |
| Chemical | ||||||||||||
| Soluble solids | SSC | 6.86 | 0.31 | 13.53 ** | 0.93 | 4.99 | 5.73–4.36 | 16.52 | 5.44 | 5.67–5.05 | 15.30 | *** |
| pH | pH | 1.44 | 0.36 | 36.05 ** | 0.94 | 4.52 | 4.6–4.41 | 4.21 | 4.40 | 4.47–4.34 | 4.57 | *** |
| SSC per plant | BY | 248.62 | 0.22 | 5.01 * | 0.89 | 18.23 | 19.9–15.09 | 36.54 | 16.73 | 19.96–12.87 | 41.99 | * |
| Stem | ||||||||||||
| Trichomes | TRIC | 2.99 | 0.17 | 6.37 ** | 0.86 | 1.00 | 1.00–1.00 | 0.00 | 1.30 | 3.00–1.00 | 31.89 | ns |
# Acronyms. a TY, RY, GY, PW, and BM are expressed in Kg; FW in g; FL and FD in cm; SSC is expressed in °Brix; CE, CI, PER, PUF, and TRIC are in scale; and FS, HI, and BY are derived indices. b MS: mean squares, R: R-square, CV: coefficient of variation, and H2: heritability. c * Indicate significance at p < 0.05, ** indicate significance at p < 0.01, and *** indicate significance at p < 0.001; ns = not significant.
Table 4.
Analysis of variance and significant levels for the genotypic (G) and environmental effects due to growing season (Y) for Solanum habrochaites (SH) ILs and the SL control, and the combined effects for the traits evaluated.
Table 4.
Analysis of variance and significant levels for the genotypic (G) and environmental effects due to growing season (Y) for Solanum habrochaites (SH) ILs and the SL control, and the combined effects for the traits evaluated.
| Trait | Genotype (G) df = 39 | Year (Y) df = 2 | G × Y df = 78 | Total Error df = 1108 | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| MS | TSS% | F | MS | TSS% | F | MS | TSS% | F | ||
| Yield | ||||||||||
| TY | 11.60 | 45.49 | 11.57 ** | 151.98 | 30.57 | 151.72 ** | 3.05 | 23.94 | 3.04 ** | 1109.89 |
| RY | 12.34 | 38.77 | 13.07 ** | 234.09 | 37.72 | 248.08 ** | 3.74 | 23.51 | 3.97 ** | 1045.53 |
| GY | 0.56 | 12.35 | 3.92 ** | 94.32 | 65.21 | 663.08 ** | 0.46 | 22.44 | 3.21 ** | 125.60 |
| UM | 946.65 | 39.85 | 13.1 ** | 13,639.01 | 29.45 | 188.68 ** | 364.62 | 30.70 | 5.04 ** | 80,094.20 |
| PW | 0.76 | 58.29 | 7.97 ** | 4.28 | 16.74 | 44.61 ** | 0.16 | 24.97 | 1.71 ** | 106.19 |
| BM | 12.69 | 43.66 | 9.43 ** | 191.67 | 33.82 | 142.5 ** | 3.27 | 22.53 | 2.43 ** | 1490.35 |
| HI | 0.08 | 62.64 | 17.76 ** | 0.10 | 4.07 | 22.48 ** | 0.02 | 33.30 | 4.72 ** | 4.85 |
| Fruit | ||||||||||
| FW | 2110.45 | 76.39 | 11.4 ** | 3122.66 | 5.80 | 16.87 ** | 245.99 | 17.81 | 1.33 * | 205,111.16 |
| FL | 3.11 | 81.24 | 23.99 ** | 0.95 | 1.27 | 7.33 ** | 0.33 | 17.49 | 2.58 ** | 22,962.97 |
| FD | 1.59 | 78.27 | 18.65 ** | 1.60 | 4.03 | 18.59 ** | 0.19 | 17.70 | 2.10 ** | 15,829.19 |
| FS | 0.07 | 73.70 | 25.90 ** | 0.01 | 0.62 | 4.23 * | 0.01 | 25.68 | 4.51 ** | 2.87 |
| CE | 2.08 | 85.50 | 19.98 ** | 0.65 | 0.68 | 6.29 * | 0.33 | 13.80 | 3.23 ** | 110.28 |
| CI | 2.92 | 89.75 | 23.70 ** | 0.05 | 0.04 | 0.38 NS | 0.33 | 12.96 | 2.69 ** | 90.69 |
| PER | 1.49 | 62.66 | 9.39 ** | 0.61 | 1.31 | 3.84 * | 0.43 | 36.02 | 2.70 ** | 176.04 |
| PUF | 2.48 | 54.72 | 15.84 ** | 16.58 | 18.79 | 106.03 ** | 0.60 | 26.50 | 3.83 ** | 173.06 |
| Chemical | ||||||||||
| SSC | 6.43 | 58.74 | 17.97 ** | 37.24 | 17.43 | 104.04 ** | 1.30 | 23.83 | 3.65 ** | 395.49 |
| pH | 0.44 | 77.47 | 18.29 ** | 3.09 | 14.33 | 128.6 ** | 0.02 | 8.20 | 0.96 NS | 25.51 |
| BY | 307.53 | 46.61 | 10.40 ** | 3863.56 | 30.03 | 130.6 ** | 77.08 | 23.36 | 2.61 ** | 32,776.91 |
| Stem | ||||||||||
| TRIC | 3.13 | 90.08 | 12.68 ** | 32.15 | 1.01 | 3.25 NS | 2.49 | 8.91 | 10.08 * | 27.32 |
df: degrees of freedom, * significant at p < 0.05, ** significant at p < 0.01, and NS: not significant. MS: mean squares, Err: error, TSS: total sum of squares, and F = F ratio.
Table 5.
List of putative chromosome regions underlying the variations of the traits (QTL) identified across two or three independent trials. PV = percent of the phenotypic variance estimated and - = not available. In bold are indicated the positive PV (%). In capital letters are indicated those QTLs that increased the performance of the traits compared to the control. BIN sizes are determined considering the overlapping regions plus the half-intervals between flanking markers and the nearest mapped loci (see Figure 4).
Table 5.
List of putative chromosome regions underlying the variations of the traits (QTL) identified across two or three independent trials. PV = percent of the phenotypic variance estimated and - = not available. In bold are indicated the positive PV (%). In capital letters are indicated those QTLs that increased the performance of the traits compared to the control. BIN sizes are determined considering the overlapping regions plus the half-intervals between flanking markers and the nearest mapped loci (see Figure 4).
| Trait | ILs | Marker(s) | Chr | BIN n° | BIN Position (cM) | BIN Size (cM) | QTL Name * | PV (%) 2 Years | PV (%) 3 Years |
|---|---|---|---|---|---|---|---|---|---|
| TY | TA1128 | TG184 | 1 | 1.2 | 10–40.5 | 30.5 | ty1.2 | −39.15 | |
| TA1535 | TG59 | 1 | 1.3 | 61–79.25 | 18.3 | ty1.3 | −50.62 | ||
| TA1276-TA1649 | TG42, TG359, CT243 | 3 | 3.5 | 68.0–156.0 | 88.0 | ty3.5 | −52.35 | ||
| TA1551 | TG313, CT234, TG408 | 10 | 10.1 | 0.00–57.0 | 57 | ty10.1 | −45.21 | ||
| TA1554 | TG651 | 11 | 11.1 | 11.0–36.0 | 25 | ty11.1 | −44.51 | ||
| TA1121 | CT211, TG283, CT287a | 12 | 12.3 | 32.25–78.5 | 43.25 | ty12.3 | −33.51 | ||
| RY | TA1535 | TG59 | 1 | 1.3 | 61–79.25 | 18.3 | ry1.3 | −51.38 | |
| TA1276-TA1649 | TG42, TG359, CT243 | 3 | 3.5 | 68.0–156.0 | 88.0 | ry3.5 | −52.30 | ||
| TA1551 | TG313, CT234, TG408 | 10 | 10.1 | 0.00–57.0 | 57 | ry10.1 | −41.28 | ||
| UM | TA1535 | TG59 | 1 | 1.3 | 61–79.25 | 18.3 | um1.3 | −14.25 | |
| TA1280 | TG15, TG370, TG287 | 4 | 4.1 | 3.25–70.5 | 67.25 | um4.1 | −9.95 | ||
| PW | TA1649 | TG244 | 3 | 3.6 | 156.0–171 | 15.0 | PW3.6 | 73.92 | |
| BM | TA1128 | TG184 | 1 | 1.2 | 10–40.5 | 30.5 | bm1.2 | −40.24 | |
| TA1535 | TG59 | 1 | 1.3 | 61–79.25 | 18.3 | bm1.3 | −44.70 | ||
| TA1276 | TG42, TG359 | 3 | 3.4 | 68.0–114.45 | 46.45 | bm3.4 | −52.03 | ||
| TA1551 | TG313, CT234, TG408 | 10 | 10.1 | 0.00–57.0 | 57 | bm10.1 | −40.79 | ||
| TA1554 | TG651 | 11 | 11.1 | 11.0–36.0 | 25 | bm11.1 | −47.15 | ||
| HI | TA1535 | TG59 | 1 | 1.3 | 61–79.25 | 18.3 | hi1.3 | −12.14 | |
| TA1276 | TG42, TG359 | 3 | 3.4 | 68.0–114.45 | 46.45 | hi3.4 | −25.52 | ||
| TA1542-TA1562 | TG264 | 4 | 4.3 | 70.5–84.25 | 13.75 | hi4.3 | −16.43 | ||
| TA1649 | TG244 | 3 | 3.6 | 156.0–171 | 15.0 | hi3.6 | −17.14 | ||
| FW | TA1266-TA1537-TA1649-TA1535 | TG620 | 2 | 2.4 | 110.5–131.5 | 21.0 | fw2.4a | −30.49 | |
| TA1276-TA1649 | TG42, TG359, CT243 | 3 | 3.5 | 68.0–156.0 | 88.0 | fw3.5 | −27.91 | ||
| TA1303-TA1304 | TG216, CT195 | 7 | 7.3 | 50.4–113 | 62.6 | fw7.3 | −42.33 | ||
| TA1330 | CT198, CT112 | 9 | 9.3 | 63.1–109.0 | 45.9 | fw9.3 | −32.29 | ||
| TA1545- TA1554-TA1555 | TG233 | 10 | 10.3 | 78.5–86.0 | 7.5 | fw10.3 | −35.63 | ||
| TA1554 | TG651 | 11 | 11.1 | 11.0–36.0 | 25 | fw11.1 | −43.14 | ||
| TA1330-TA1553 | TG36 | 11 | 11.3 | 93.5–69.0 | 24.5 | fw11.3 | −32.01 | ||
| TA1554-TA1221 | TG296 | 12 | 12.4 | 78.5–108.0 | 29.5 | fw12.4 | −36.44 | ||
| FL | TA1541 | TG479 | 3 | 3.1 | 0.0–16.0 | 16.0 | fl3.1 | −13.82 | |
| TA1276 | TG42, TG359 | 3 | 3.4 | 68.0–114.45 | 46.45 | fl3.4 | −16.85 | ||
| TA1280 | TG15, TG370, TG287 | 4 | 4.1 | 3.25–70.5 | 67.25 | fl4.1 | −9.95 | ||
| TA1543 | TG60 | 5 | 5.3 | 99.5–119.0 | 19.5 | fl5.3 | −9.86 | ||
| TA1545 | TG352, TG164 | 6 | 6.3 | 21.75–50.95 | 29.2 | fl6.3 | −11.38 | ||
| TA1303-TA1304 | TG216, CT195 | 7 | 7.3 | 50.4–113 | 62.6 | fl7.3 | −19.05 | ||
| TA1545-TA1554-TA1555 | TG233 | 10 | 10.3 | 78.5–86.0 | 7.5 | fl10.3 | −18.7 | ||
| TA1330-TA1553 | TG36 | 11 | 11.3 | 93.5–69.0 | 24.5 | fl11.3 | −15.92 | ||
| TA1121 | CT211, TG283, CT287a | 12 | 12.3 | 32.25–78.5 | 43.25 | fl12.3 | −14.83 | ||
| FD | TA1105 | CT140, TG554, TG553 | 2 | 2.1 | 8.0–70.15 | 62.2 | fd2.1 | −17.32 | |
| TA1276-TA1649 | TG42, TG359, CT243 | 3 | 3.5 | 68.0–156.0 | 88.0 | fd3.5 | −14.71 | ||
| TA1473 | TG464, TG163 | 4 | 4.7 | 119–135 | 16.0 | fd4.7 | −11.21 | ||
| TA1303-TA1304 | TG216, CT195 | 7 | 7.3 | 50.4–113 | 62.6 | fd7.3 | −19.03 | ||
| TA1554 | TG651 | 11 | 11.1 | 11.0–36.0 | 25 | fd11.1 | −16.82 | ||
| TA1555 | TG393 | 11 | 11.4 | 93.5–104.0 | 10.5 | fd11.4 | −15.53 | ||
| TA1554- TA1221 | TG296 | 12 | 12.4 | 78.5–108.0 | 29.5 | fd12.4 | −16.14 | ||
| FS | TA1258 | TG607 | 1 | 1.4 | 79.25–114–75 | 35.5 | FS1.4 | 7.28 | |
| TA1258-TA523-TA1223 | TG245, TG17 | 1 | 1.5 | 114.75–153.5 | 38.8 | FS1.5 | 7.42 | ||
| TA1266-TA1537-TA1649-TA1535 * | TG620 | 2 | 2.4 | 110.5–131.5 | 21.0 | fs2.4b | −11.89 | ||
| TA1541 | TG479 | 3 | 3.1 | 0.0–16.0 | 16.0 | fs3.1 | −5.96 | ||
| TA1554 | TG651 | 11 | 11.1 | 11.0–36.0 | 25 | fs11.1 | −12.73 | ||
| CE | TA1258 | TG607 | 1 | 1.4 | 79.25–114–75 | 35.5 | ec1.4 | −15.11 | |
| TA1266-TA1537-TA1649-TA1535 * | TG620 | 2 | 2.4 | 110.5–131.5 | 21.0 | ec2.4c | −38.56 | ||
| TA1473-TA1475 | CD39 | 4 | 4.6 | 101.5–119 | 17.5 | EC4.6 | 24.27 | ||
| TA1330 | CT198, CT112 | 9 | 9.3 | 63.1–109.0 | 45.9 | EC9.3 | 21.6 | ||
| TA1330-TA1553 | TG36 | 11 | 11.3 | 69.0–93.5 | 24.5 | EC11.3 | 18.97 | ||
| CI | TA1535 | TG59 | 1 | 1.3 | 61–79.25 | 18.3 | ic1.3 | −25.85 | |
| TA1258 | TG607 | 1 | 1.4 | 79.25–114–75 | 35.5 | ic1.4 | −25.19 | ||
| TA1266-TA1537-TA1649-TA1535 * | TG620 | 2 | 2.4 | 110.5–131.5 | 21.0 | ic2.4c | −53.95 | ||
| TA1475 | TG345 | 4 | 4.5 | 96.75–101.5 | 4.75 | IC4.5 | 25.08 | ||
| PER | TA1287 | CT101 | 5 | 5.1 | 0.0–20.0 | 20.0 | PER5.1 | 32.49 | |
| TA1554 | TG651 | 11 | 11.1 | 11.0–36.0 | 25 | per11.1 | −38.69 | ||
| SSC | TA1258-TA523-TA1223 | TG245, TG17 | 1 | 1.5 | 114.75–153.5 | 38.8 | SSC1.5 | 24.01 | |
| TA1649 | TG244 | 3 | 3.6 | 156.0–171 | 15.0 | SSC3.6 | 26.36 | ||
| TA1543 | TG60 | 5 | 5.3 | 99.5–119.0 | 19.5 | SSC5.3 | 25.37 | ||
| TA1551 | TG313, CT234, TG408 | 10 | 10.1 | 0.00–57.0 | 57 | SSC10.1 | 26.47 | ||
| PH | TA1535 | TG59 | 1 | 1.3 | 61–79.25 | 18.3 | ph1.3 | −8.91 | |
| TA1473 | TG464, TG163 | 4 | 4.7 | 119–135 | 16.0 | ph4.7 | −4.73 | ||
| TA1554 | TG651 | 11 | 11.1 | 11.0–36.0 | 25 | ph11.1 | −8.50 | ||
| TA1330 | CT198, CT112 | 9 | 9.3 | 63.1–109.0 | 45.9 | ph9.3 | −4.24 | ||
| BY | TA1128 | TG184 | 1 | 1.2 | 10–40.5 | 30.5 | by1.2 | −39.33 | |
| TA1535 | TG59 | 1 | 1.3 | 61–79.25 | 18.3 | by1.3 | −49.15 | ||
| TA1276 | TG42, TG359 | 3 | 3.4 | 68.0–114.45 | 46.45 | by3.4 | −57.73 | ||
| TRIC | TA1280 | TG15, TG370, TG287 | 4 | 4.1 | 3.25–70.5 | 67.25 | TRIC4.1 | 100 | |
| TA1304 | TG202 | 7 | 7.2 | 23.9–50.4 | 26.5 | TRIC7.2 | 125 | ||
| TA1551 | TG313, CT234, TG408 | 10 | 10.1 | 0.00–57.0 | 57 | TRIC10.1 | 200 |
* QTL name = trait abbreviation + BIN n°. a Not considering the additional wild introgression from TA1535. b Not considering the additional wild introgression from TA1535–TA1649. c Not considering the additional wild introgression from TA1649.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Tripodi, P.; Vitiello, A.; D’Onofrio, B.; Parisi, M.; Cammareri, M. Dissecting the Genotypic and Environmental Factors Underpinning the Quantitative Trait Variation in a Set of Wild Tomato (Solanum habrochaites LA1777) Introgression Lines. Agronomy 2021, 11, 38. https://doi.org/10.3390/agronomy11010038
AMA Style
Tripodi P, Vitiello A, D’Onofrio B, Parisi M, Cammareri M. Dissecting the Genotypic and Environmental Factors Underpinning the Quantitative Trait Variation in a Set of Wild Tomato (Solanum habrochaites LA1777) Introgression Lines. Agronomy. 2021; 11(1):38. https://doi.org/10.3390/agronomy11010038
Chicago/Turabian StyleTripodi, Pasquale, Antonella Vitiello, Bruno D’Onofrio, Mario Parisi, and Maria Cammareri. 2021. "Dissecting the Genotypic and Environmental Factors Underpinning the Quantitative Trait Variation in a Set of Wild Tomato (Solanum habrochaites LA1777) Introgression Lines" Agronomy 11, no. 1: 38. https://doi.org/10.3390/agronomy11010038
APA StyleTripodi, P., Vitiello, A., D’Onofrio, B., Parisi, M., & Cammareri, M. (2021). Dissecting the Genotypic and Environmental Factors Underpinning the Quantitative Trait Variation in a Set of Wild Tomato (Solanum habrochaites LA1777) Introgression Lines. Agronomy, 11(1), 38. https://doi.org/10.3390/agronomy11010038
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
