ABSTRACT: To be considered drought-tolerant, a tomato cultivar is required to present some level of tolerance at all developmental stages of plant growth. Since drought tolerance is a stage-specific phenomenon, genotype assessment must be performed separately at all developmental stages. In this study, we used a multi-trait index based on factor analysis and genotype-ideotype distance (FAI-BLUP index) to properly rank 49 tomato genotypes according to their tolerance to drought stress conditions at germination and seedling stages. Seeds of 47 introgression lines, which cultivar M82 is considered drought-sensitive, and the Solanum pennellii wild accession LA 716, which is considered drought-tolerant, were subjected to a control condition, where seeds were soaked in distilled water, and a drought condition, where seeds were soaked in a polyethylene glycol (PEG) solution (-0.3 MPa). Drought stress, induced by PEG, had a significant impact on all nine germination and growth performance-related traits; there was a reduction in shoot length (SL), total length (TL), initial germination percentage (IGP), final germination percentage (FGP), germination velocity index (GVI), and germination rate (GR). In contrast, the root-to-shoot ratio (R/S) and time to reach 50% germination (T50) increased under drought stress. Root length (RL) was less affected by drought, and in some genotypes, it was even increased. As expected, LA 716 ranked closest to the drought tolerance ideotype. IL 1-4-18, IL 2-3, IL 1-2, IL 9-2, and IL 10-1 were the most drought-tolerant at the germination stage. These results will serve as guidance for breeders who are aiming at developing drought-resistant tomato cultivars.
Keywords: Solanum lycopersicum L., Solanum pennellii Corr., genotype-ideotype distance, drought tolerance, water stress.
PLANT BREEDING
Multi-trait selection of tomato introgression lines under drought-induced conditions at germination and seedling stages
Received: 17 September 2020
Accepted: 27 January 2021
Drought is one of the most significant environmental constraints that limits crop growth and yield in agriculture worldwide (Thabet, Moursi, Karam, Graner, & Alqudah, 2018). It is already known that a growing world population requires the production of more food by using less area and inputs such as water and fertilizers. The development of cultivars with tolerance to water deficit is, therefore, a sustainable and economically viable approach to enhance yield and ensure food security for the human population (Merchuk-Ovnat et al., 2016). Thus, plant breeding for drought tolerance is a promising strategy since it allows crop yield in water-limited environments and also makes it possible to grow crops using less water (Cattivelli et al., 2008).
Tomato (Solanum lycopersicum L.) is one of the most important vegetables in the world owing to its socioeconomic and nutritional importance (Rigano et al., 2016). Cultivated tomato is sensitive to drought stress from seed germination to harvest (Foolad, Zhang, & Subbiah, 2003). Tomato plants require approximately 0.89 L to 2.3 L of water daily during their growing cycle (Ayankojo, Morgan, Ozores-Hampton, & Migliaccio, 2018; Goyal & Sharma, 2018), with qualitative and quantitative losses during a drought period (George, Ahmad Jatoi, & Uddin Siddiqui, 2013; Zhou et al., 2019).
The development of drought-tolerant tomato plants faces several challenges. First, the trait is controlled by many genes with different effects and is highly influenced by environmental variation (Florido Bacallao & Bao Fundora, 2014; Foolad & Foolad, 2004). Thus, the precision and repeatability of field trials to assess drought tolerance are reduced because variations in environmental conditions such as rainfall, temperature, salinity, and nutrient availability cannot be fully controlled (Mohammadi, 2016). Another problem is that drought tolerance is a developmentally regulated stage-specific phenomenon. Therefore, specific stages throughout the ontogeny of plants must be evaluated separately for the assessment of drought tolerance and identification of its genetic components (Foolad, 2007).
Most of the recent tomato cultivars are drought-sensitive at all stages of plant development; the stage of seed germination and early seedling growth are the most sensitive (Mishra, Rai, Kumar, Singh, & Pandey, 2016). Genetic variability among the modern cultivated tomato is limited owing to natural and artificial selection during domestication and evolution (Steinhauser et al., 2011), which makes tomato wild relatives a good source of drought tolerance genes. Solanum pennellii Corr. is one of the crossable wild relatives of the cultivated tomato that displays drought tolerance (Bolger et al., 2014). Eshed and Zamir (1995) generated a collection of introgression lines (ILs) in which defined genomic segments of wild species replaced homologous regions in the background of the cultivated tomato variety M82 (S. lycopersicum). The IL population provides complete coverage of the wild species genome, aiding the investigation of the wild gene reservoir. Moreover, any phenotypic difference found between an IL and the recurrent parent M82 subjected to a drought condition can be attributed solely to the introgressed chromosomal segment of the drought-tolerant S. pennellii (Lippman, Semel, & Zamir, 2007). Thus, they are great sources for identifying genomic regions associated with drought tolerance and interactions between genomic regions (Schauer et al., 2006).
The screening of drought-tolerant genotypes during germination and seedling stages is often performed using polyethylene glycol (PEG) to simulate drought stress conditions (Ali et al., 2018; Feng, Wang, Shao, Zhang, & Tang, 2018; Ghebremariam, Liang, Li, Li, & Qin, 2013; Toscano, Romano, Tribulato, & Patanè, 2017). Drought tolerance in these trials was assessed by evaluating the genotype performance for several traits by using data from germination and seedling growth tests on a daily basis (Foolad et al., 2003), and genotype ranking for all the evaluated traits is performed simultaneously using selection indexes (Baker, 1986). Most of the available selection indexes require orthogonality among traits. Germination and seedling traits, however, are usually multicollinear. The FAI-BLUP, a multi-trait index based on factor analysis and ideotype design, proposed by Rocha, Machado, and Carneiro (2018), is used in this study; it combines factor analysis and genotype-ideotype design by using a multi-trait approach free from multicollinearity and has been successfully used for multi-trait selection of genotypes for several purposes (Oliveira et al., 2019; Rocha et al., 2019; Silva et al., 2018; Woyann et al., 2019).
Therefore, the purpose of the present study was to use the FAI-BLUP index to select drought-tolerant ILs and then identify the tomato genome regions associated with drought tolerance at germination and seedling stages. To be considered a drought-resistant, a tomato genotype must exhibit a fair level of drought resistance at all developmental stages of plant growth. Our findings will serve as guidance for breeders aiming at developing drought-resistant tomato cultivars.
Forty-seven ILs from the collection developed by Eshed and Zamir (1995) covering twelve tomato chromosomes (Figure 1) were chosen and tested for drought tolerance at germination and seedling stages. The cultivar M82, considered to be drought-sensitive, and S. pennellii accession LA716, considered to be drought-tolerant, were included in the experiment as testers.
Seeds of all genotypes were produced under uniform and optimum environmental conditions at the Research and Extension Farm Unit Horta Velha in Universidade Federal de Viçosa (UFV), Minas Gerais State, Brazil (20°45’14’ S; 42°52’53’ W). Seeds were collected from healthy mature fruits and subjected to natural fermentation for two days to remove the mucilaginous material surrounding the seed coat. After fermentation, the solution was passed through a fine-mesh sieve and washed with tap water to remove the mucilage. The cleaned seeds were then dried on filter paper at ambient temperature and stored in a storage room (13°C and 30% RH) inside paper bags until further use.

The study was performed in the Seed Laboratory at UFV. To assess drought tolerance at the germination and seedling stages, seeds of all genotypes were subjected to two different conditions: a control condition, where seeds were soaked in distilled water, and a drought condition, where seeds were soaked in a PEG solution (-0.3 MPa) to induce drought stress. The osmotic potential of -0.3 MPa was chosen here because it was reported to be the most suitable for selecting drought-tolerant tomato genotypes during germination and seedling stages (Borba, Maciel, Marquez, Junior, & Nogueira, 2017).
The experiment was arranged in a completely randomized design with four replicates. Each replicate consisted of an acrylic square box (11 x 11 x 3 cm) containing 25 seeds. The sterilized acrylic boxes were lined up with two layers of paper towels, moistened with distilled water in the control condition and a PEG-6,000 solution (-0.3 MPa) in the drought condition. The boxes were then placed inside plastic bags to avoid water loss and incubated in a growth chamber at a constant temperature of 25°C. Distilled water was sprayed on the paper to replace the evaporated water whenever it was noticed that they were dry. After incubation, the boxes were examined daily for 14 days to determine the number of germinated seeds (normal seedlings). On the 14th day, seedling length was measured. Both daily germination and seedling length data were analyzed using the SeedCalc R package (Silva, Medeiros, & Oliveira, 2019) to estimate nine germination and growth performance-related traits, as described in Table 1.

To assess genotype response to drought, the means of each parameter of drought-stressed seeds were divided by the means of their respective controls, as described by Thabet et al. (2018). Data were subjected to analysis of variance (ANOVA) to confirm significant differences between genotypes for each trait by using the GENES software (Cruz, 2013).
The univariate statistical model used to predict the genotypic values was denoted by:
where: y is the data vector, m is the vector of the effects of the measurement- combined replicates (assumed as fixed) added to the overall mean; g = vector of genotype effects (assumed as random); ɛ = vector of residue (random); X and Z represent the incidence matrices for these effects.
The transformed means were then used to rank the ILs according to their level of tolerance to drought using a multi-trait index based on factor analysis and genotype-ideotype distance (FAI-BLUP index), using the R script provided by Rocha et al. (2018) and R software (R Core Team, 2017). This methodology was adopted because it deals very well with multicollinearity problems, considers the correlation structure observed, and drives genotype ranking closer to the one hypothesized by researchers as the ideotype.
Principal component analysis was used to extract factorial loads of the phenotypic correlation matrix, obtained from the predicted genotypic values. Varimax criterion described by Kaiser (1958) was used for the analytic rotation, and the calculation of the factor scores of the weighted least squares was performed using the method proposed by Bartlett (1938).
The number of ideotypes was defined based on the combination of desirable and undesirable factors according to the selection purpose. The following algorithm gives the number of ideotypes: NI = 2n, where: NI = number of ideotypes and n = number of factors.
The drought-tolerant ideotype was built considering the maximum transformed mean values for IGP, GVI FGP, SL, RL, TL, and R/S and minimum mean values for T50 and GR; it is expected that a drought-tolerant genotype will show this performance. In the present study, besides the drought-tolerant ideotype, we further analyzed the ideotype for each of the four factors individually to better understand the performance of the genotype under drought.
After ideotype determination, genotype-ideotype distances were estimated and converted into spatial probability, enabling genotype ranking. The following algorithm was used:
where: Pij is the probability of the ith genotype (i = 1, 2, ... , n) to be similar to the jth ideotype (j = 1, 2,. . ., m); dij = Genotype-ideotype distance from the ith genotype to the jth ideotype, based on standardized mean Euclidean distance.
The ANOVA detected significant genotype effect (p < 0.05) for all the evaluated traits in the two conditions. Therefore, all the nine traits were used in the following analyses.
Drought stress induced by the PEG 6,000 solution (-0.3 MPa) had a considerable impact on all evaluated traits (Figure 2). Under drought conditions, SL was reduced for most genotypes, which probably led to reductions in TL. RL was less affected by drought, and for some genotypes, it was even increased. The R/S increased in almost all the genotypes (Figure 2), but decrease in SL was more evident. IGP, FGP, GVI, and GR were also severely reduced under drought-induced conditions for all genotypes, as shown in Figure 2. IGP showed the highest reduction under drought condition compared to the other traits (Figure 2). Consistent with this result, the T50 of all genotypes increased under drought condition (Figure 2). Even though different performances were observed between ILs, in general, a decrease in germination percentage and germination velocity, followed by a reduction in seedling growth occurred under drought stress (Figure 2).

As for the principal component analysis, we found that from the eigenvalues and cumulative frequency for the nine principal components obtained via the phenotypic correlation matrix, the first four principal components had eigenvalues higher than one (Table 2). Thus, according to Kaiser’s criterion (Kaiser, 1958), the data may be condensed (dimensional reduction) into four factors. The cumulative frequency for the first four principal components or communality mean (common variance) was higher than 84%, indicating that they are sufficient to represent 84% of all data variability (Table 2).

After varimax rotation (Table 3), a high genetic correlation for the first factor was observed among the traits IGP, FGP, and GVI, and we named this factor germination. For the second factor, a high genetic correlation was observed among SL, RL, and TL, and we named it size. The third factor was named velocity of germination, to which T50 and GR were strongly correlated. The fourth factor was named root-to-shoot ratio and was formed only by R/S. The genetic correlations among traits within a factor can be given in the same and/or opposite directions.

Table 4 shows genotype ranking according to the FAI-BLUP index and the probability associated with the distance from a given genotype to the ideotype. Besides the drought-tolerant ideotype, which shows desirable traits for all the five factors simultaneously, we decided to analyze the ideotype for each factor individually in order to better understand genotype performance under drought.
The accession LA716 (S. pennellii) was ranked closest to the drought tolerance ideotype during germination and seedling stages, followed by the introgression lines IL 1-4-18, IL 2-3, IL 1-2, IL 9-2, and IL 10-1. For the germination ideotype, the closest genotype was IL 4-2, whereas M82 was ranked 26th and LA716 44th. For velocity of germination, IL 11-4-1 was ranked closest to the ideotype, cultivar M82 was ranked 24th, and LA716 was the furthest genotype from the germination ideotype. The IL 10-2-2 was ranked closest to the length ideotype, and the cultivar M82 ranked 18th and LA716 47th. IL 1-4-18 was ranked closest to the root-to-shoot ratio ideotype.

Water scarcity for agricultural production is a rising concern, especially for high water-demanding crops such as tomato, and plant breeding appears to be an economically viable approach to ensure food production in such arid areas. Since genetic tolerance to drought is a stage-specific phenomenon, the first step of a breeding program for drought tolerance is to identify drought-tolerant genotypes during the initial stages of plant growth. In this study, we used a multicollinearity free multi-trait index to properly rank 50 tomato genotypes according to their level of drought tolerance during germination and early seedling stages.
As expected, we found a significant genotype effect for all germination and initial seedling-growth traits, revealing genetic variability within the IL population. This genetic variability among ILs is associated with the small genomic segments of S. pennellii that replaced the homologous regions in the M82 background (Eshed & Zamir, 1995).
The low water potential of the external environment, promoted by PEG, made water less available to the seeds; hence, germination became either delayed or occurred at a reduced rate (Partheeban, Chandrasekhar, Jeyakumar, Ravikesavan, & Gnanam, 2017), which explains why drought imposed by PEG had a considerable impact on all evaluated traits (Figure 2). These drought effects on germination parameters were also recurrent for several other crops (Khodarahmpour, 2011). Low water potentials also result in reduced seedling growth owing to restricted cell enlargement and cell division (Borba et al., 2017). Consistent with other studies, shoot and root growths were inhibited by drought, and because shoot growth is more sensitive than root growth, the root-to-shoot ratio typically increases (Li, Mu, & Lin, 2014; Naveed, Mitter, Reichenauer, Wieczorek, & Sessitsch, 2014). This means that under drought stress, plants allocate more resources for root growth than for shoot growth in order to enhance water acquisition and limit evaporation (Lynch, 2007; Sharp, 2004).
Since genotype performance under drought conditions was very diverse, as shown in Figure 2, a given genotype can show an increase for some traits and a decrease for other traits. Therefore, we decided to utilize a multi-trait approach to properly rank and select the most drought-resistant genotypes. The FAI-BLUP index was chosen over other multivariate approaches because it performs the selection without multicollinearity problems. Moreover, it allows breeders to design an ideotype based on purpose and according to their dataset. In this way, the genotype ranking, associated simultaneously with spatial probability for multiple traits, is performed without the need to assign weights to each of the evaluated traits (Rocha et al., 2018).
This method was first proposed and applied for multi-trait selection of elephant grass (Pennisetum purpureum Schum.) for bioenergy purposes (Rocha et al., 2018). It was also recently used for multi-trait selection in several other studies. Woyann et al. (2019) used the FAI-BLUP index to select soybean (Glycine max [L.] Merrill) genotypes that best fit biodiesel production. Silva et al. (2018) evaluated lines and hybrids of biomass sorghum (Sorghum bicolor (L.) Moench) by using the same index. Oliveira et al. (2019) also used this index for the selection of photosensitive sorghum hybrids for energy cogeneration. Rocha et al. (2019) used the FAI-BLUP index for the selection of superior inbred progenies of common beans (Phaseolus vulgaris L.).
Multicollinearity is a common problem when handling multiple traits, even more, when traits are not orthogonal, as in the case of our study. It provides inflated errors, leading to inaccurate results regarding significance tests of predictors, and hence affects the selection process (Dormann et al., 2013; Prunier, Colyn, Legendre, Nimon, & Flamand, 2015). The use of factor analysis, as in the FAIBLUP index, produces uncorrelated or orthogonal axes among the final factors, generating scores without multicollinearity (Dormann et al., 2013). In addition, the approach, based on genotype-to-ideotype distance, focuses on multiple traits simultaneously, allowing breeders to consider all the evaluated traits and to underly plant morphological and physiological nuances (Van Oijen & Höglind, 2016). In fact, by analyzing drought tolerance ranking, germination, velocity of germination, size, and root-to-shoot ideotypes, we could better understand the behavior of ILs under drought conditions.
IL 4-2 was ranked closest to the germination ideotype, whereas M82 was ranked 26th and LA716 44th. Introgressed fragments of S. pennellii in the M82 background promoted an increase in the germination of 25 out of the 47 ILs, and 23 out of 47 ILs in the velocity of germination. However, these fragments alone did not result in the desired germination traits, for example, LA716 was one of the worst genotypes concerning germination and velocity of germination (Table 4). One possible hypothesis is that these traits are controlled by more than one gene and, hence, the increased germination and velocity of germination observed is probably due to the interaction between S. pennellii introgressed fragments and the rest of the M82 genome.
IL 10-2-2 was ranked closest to the length ideotype, M82 was ranked 18th and LA716 47th (Table 4). Therefore, we may infer that the S. pennellii introgressed fragments in M82 promoted an increase in seedling length-related traits for 17 out of the 47 evaluated ILs. However, as previously discussed, the M82 genome seems to be essential for its expression. Another explanation for why LA716 was found between the worst genotypes for seedling growth is that LA 716 seeds are smaller than the seeds from the other evaluated genotypes. Larger seeds tend to have a higher amount of nutritional reserves, which enable greater initial seedling development even under stress conditions such as drought (Pereira, Pereira, & Dias, 2013).
IL 1-4-18 was ranked closest to the root-to-shoot ratio ideotype (Table 4). The same IL was also closest to the drought tolerance ideotype. Some studies have shown that tolerant genotypes tend to sustain root-to-shoot ratios, whereas sensitive genotypes show a decrease in root-to-shoot ratios (Shamim, Saqlan, Habib-Ur-Rehman, & Waheed, 2014; Rahman, Nawata, & Sakuratani, 1999). In other crops, plant growth performance was positively associated with well-developed root systems and early seedling vigor, both of which can help to improve stress tolerance (Chloupek, Dostál, Středa, Psota, & Dvořáčková, 2010). Thus, focusing on such traits would be an efficient approach for phenotyping a large number of genotypes under drought conditions in a short time. Hence, selection based on shoot and root traits would be useful to identify genotypes for better performance under drought stress conditions. An important finding from this study is that the relationship between shoot and root, which have already been successfully implemented in existing high-throughput screening platforms, appears to be the most effective way to screen drought-tolerant genotypes prior to more laborious field trials. Similar findings were also observed for cotton and maize (Avramova et al., 2016; Riaz et al., 2013).
Although LA 716 did not perform well in terms of seedling germination and growth, it was ranked closest to the drought tolerance ideotype during germination (Table 4). As previously discussed, LA 716 seeds are smaller than those of the other genotypes. Due to the degree of contact between seed and water, sensitivity to drought is inversely related to seed size. This means that small seeds tend to be more efficient at absorbing water owing to an increased contact surface area between seed and substrate (Jiang & Su, 2018). Lima, Santiago, Araújo, and Teixeira (2005) and Pereira et al. (2013) observed that smaller seeds were less affected by drought because in smaller seeds, the amount of water required to reach a water content is smaller than that required by larger seeds.
Considering a 15% selection intensity, we selected ILs 1-4-18, IL 2-3, IL 1-2, IL 9-2, and IL 10-1 as drought-tolerant genotypes. The name of the ILs is related to the position and chromosome where the introgression was made. For example, IL1-1 means that one fragment of chromosome 1 of S. pennellii was introduced into the background genotype of M82 (Eshed & Zamir, 1995). Thus, in our findings, fragments of different S. pennellii chromosomes (1, 2, 9, and 10) and different regions within a chromosome are involved in determining drought tolerance during germination stage. Our findings are consistent because it is known that plants respond to drought by expressing quantitative traits, which involve multiple genes and physiological mechanisms. In addition, according to Foolad and Foolad (2004), drought tolerance during germination is genetically controlled, and additivity is the main genetic component. As M82 was ranked far from the drought tolerance ideotype, the improved performance of the selected ILs is related to the introgressed chromosomal segments of S. pennellii. The introgressions contain genes that contribute to the increased drought of these ILs.
In summary, the genotypes IL 1-4-18, IL 2-3, IL 1-2, IL 9-2, and IL 10-1 were selected as drought-tolerant during germination and seedling stages. They will be evaluated at other developmental stages of plant growth and be used as a source of drought tolerance genes in tomato breeding programs. The fragment 4-18 on chromosome 1, the fragment 3 on chromosome 2, the fragment 2 on chromosomes 1 and 9, and the fragment 1 on chromosome 10 can be used to identify more accurately the portions of S. pennellii genome that enhance tolerance to drought in tomato plants. Given the information of our study and the genome sequencing of S. lycopersicum and S. pennellii, the identification of candidate genes that control morphological and physiological traits related to drought tolerance at seed level will be made possible by exploring the map of the introgressed region. The identification of polymorphisms in the key genes controlling these traits can also be used to guide tomato breeding programs for drought tolerance.
The authors would like to acknowledge the Department of Plant Science at Universidade Federal de Viçosa (UFV). This work was supported by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) (Finance Code 001), Fundação de Amparo a Pesquisa do Estado de Minas Gerais (FAPEMIG) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).
* Author for Correspondence. E-mail: herikapaulapessoa@gmail.com





