Physiological Traits of Thirty-Five Tomato Accessions in Response to Low Temperature

Tomato is often exposed to diverse abiotic stresses and cold stress is one of harsh environmental stresses. Abnormal low temperature affects tomato growth and development, leading to, e.g., physiological disorders, flower drops, and abnormal fruit morphology, and causing a decrease in tomato yield and fruit quality. It is important to identify low temperature-(LT) tolerant tomato (Solanum lycopersicum L.) cultivars relying on different fruit types. In this study, our focus was to analyze the physiological traits of 35 tomato accessions with three different fruit types (cherry, medium, and large sizes) under night temperature set-points of 15 ◦C for control temperature (CT) and 10 ◦C for LT, respectively. Plant heights (PH) of most tomato accessions in LT were remarkably decreased compared to those in CT. The leaf length (LL) and leaf width (LW) were reduced depending on the genotypes under LT. In addition, the number of fruits (NFR), fruit set (FS), fruit yield (FY), and marketable yield (MY) were negatively affected in LT. The variation was further investigated by the correlation, the principal component (PCA), and the cluster analysis. Interestingly, positive correlations between different vegetative and reproductive traits were uncovered. Multivariate analysis including the PCA and hierarchical clustering classified the LT-treated 35 tomato accessions into four major groups. The identified accessions were associated with vegetative and reproductive parameters on positive directions. The results might be utilized for establishing breeding programs on selecting LT-tolerant tomato cultivars with different selection indices relying on fruit types during vegetative and/or reproductive stages.


Introduction
Tomato (Solanum lycopersicum L.) is the most essential horticultural vegetable in the world, providing vitamin A and C, rich minerals, low calories, lycopene, and β-carotene [1]. The importance of tomato crops has been gradually emerging due to increasing commercial, marketable, and dietary values [2]. The cultivated area of tomatoes is widely expanded among agricultural crops. According to Food and Agriculture Organization (FAO, http:// www.fao.org/faostat/ (accessed on 7 May 2021)) in 2019 and Korean Statistical Information Service (KOSIS, https://kosis.kr/eng/ (accessed on 7 May 2021)) in 2021, the cultivated area and tomato production reached approximately 4.8 million hectares and 182 million metric tons in the world and around 6000 hectares and 40,020,000 metric tons in Korea, respectively. However, due to global warming and climate changes, the unpredictable agriculture weather, such as low and high temperatures, has critically limited the yields and the area of agricultural cultivation of tomato plants [3][4][5].
Tomato plant is one of sessile organisms, which experiences multiple abiotic stresses including cold stress, heat stress, high salinity stress, and drought stress during the pe-riods of vegetative and reproductive growth [6][7][8]. Low temperature (LT) is a critical factor for maintaining and improving the crop yield of tomato plants during the periods of growth and development stages [3,4]. LT (0-20 • C) above the freezing temperature (below 0 • C) plays an important role in the leaf morphology [9], the truss appearance and growth [9][10][11], and the fruit development [12][13][14][15] during vegetative and reproductive stages. Recent studies have demonstrated that LT significantly influenced plant height (PH), plant diameter (PD), leaf length (LL), and leaf width (LW) in vegetative parameters [9,10,16,17] and flowering time (FT) [16,18], the number of flowers (NFL) and fruits (NFR) [18,19], fruit set (FS) [19], and fruit yield (FY) [19,20] in reproductive parameters. Moreover, the relationships of the same traits during either vegetative or reproductive stages have been investigated, but the correlation of vegetative traits with reproductive traits remains unexplored under LT condition.
In Wanju, Korea, the maximum and minimum average temperatures have been ranging from 2.8-8.7 • C and from −6.3-0.2 • C in the winter season since 1970, respectively. Moreover, the average temperature has been ranging from −1.7-3.9 • C (https: //www.weather.go.kr/w/index.do (accessed on 9 August 2021)). Due to such extreme and fluctuated low temperatures of the area in the winter season, the heating demand is dramatically increased at night time [21,22] and the heating cost occupied approximately 20-30% of the managing cost of winter tomato cultivation from 2017-2019 [23]. The temperature control is one of the major considerable factors for the tomato cultivation in the winter greenhouse. Approximately 15 • C in the winter is maintained for the optimal temperature set-points, which provide tomatoes with conditions to grow healthy without severe cold stress [9,24,25]. The studies on the optimal temperature set-point have reported that the reduction of temperature by around 2 • C in greenhouse allowed to reduce around 16% of the winter heating cost for tomato cultivation [16,26], implying that a reduction of temperature from 15 • C to 10 • C in winter greenhouse would lead to a significant decrease in the heating cost of tomato cultivation in agriculture. However, a few studies have dissected the relationship of winter tomato cultivation and night low temperature (NLT) (10 • C) in greenhouse [4,16,21,22]. It is reasonable that practical breeding programs for low temperature (LT)-tolerant tomato cultivars economically consider keeping low temperature (10 • C) during the night.
It is essential to utilize the large number of genotypes with various fruit types to provide proper indices to breeders for selecting LT-tolerant tomatoes. Although studies have been performed on the effect of LT on tomato plants with the large number of tomato genotypes [3,16,22], most have determined the impact of LT with limited genotypes [17,25]. In addition to this, only several studies have reported the selection criteria for LT-tolerant cultivars depending on fruit types [16,27]. Thus, it is still required to study the physiological aspects in response to LT with large-scale accessions and different fruit types of tomatoes.
In this work, we investigated the physiological traits of 35 tomato genotypes with different fruit types, which were grown in two different greenhouse conditions with night temperature set-points at 10 • C for low temperature (LT) and 15 • C for control treatment (CT), respectively, and analyzed the vegetative parameters of PH, SD LL, and LW and the reproductive parameters of NFL, NFR, FS, FY, and MY with different fruit types through the correlation coefficient, PCA, and cluster analysis. This research would be helpful for understanding the relationships of multiple variables with vegetative and reproductive parameters and could be further utilized for breeding programs on selecting LT-tolerant tomato cultivars with different selection criteria relying on fruit types in greenhouse condition.

Plant Material and Growth Conditions
The plant material and growth conditions were followed as previously described in [18]. A total of 35 tomato breeding lines from the National Institute of Horticultural and Herbal Science (NIHHS, Wanju, Korea, 35 • 83 N, 127 • 03 E) were used in this research (Table S1). All accessions were classified into two wild (<10 g), twenty cherry (10-30 g), eleven medium (31-80 g), and two large (>81 g) depending on fruit sizes [18]. The seeds of 35 accessions were sown on 31 August 2020 in plastic trays (52 × 26 cm in size, 6 × 6 cm cells with pot volume 5 L) containing 1:1 sand and commercial bed soil (Bio Sangto, Seoul, Korea) containing coco peat (47.2%), peat moss (35%), zeolite (7%), vermiculite (10.0%), dolomite (0.6%), humectant (0.006%), and fertilizers (0.194%) containing 270 mg kg −1 of N, P, and K, respectively. A liter of water was provided to each tray daily, and the trays were placed in a glasshouse (26/18 • C in day/night with relative humidity within 65-70%). Seedlings with 20-25 cm height and first truss were transplanted on 28 October 2020. The seedlings were transferred into two plastic film greenhouses, where night temperature set-point was maintained at 15 • C for 14 days in both greenhouses, adapting the seedlings to new environment conditions. Subsequently, night temperature set-point of each greenhouse was controlled for low temperature (LT) at 10 • C and control temperature (CT) at 15 • C, respectively. Tomato seedlings of five plants per accession were planted with a plant distance of 140 cm by 40 cm between plants in both LT and CT greenhouses. All tomato accessions were randomly selected and planted with keeping the same arrangement of the accessions between LT and CT greenhouses.
The temperature was monitored in both LT and CT greenhouses during the periods of whole growth and development using a data logger ( Figure S1A) (WatchDog 1450, Spectrum Technologies Inc., Aurora, CO, USA) and the climate data of the region was presented during the periods of tomato growth and development ( Figure S1B). Night time temperature was maintained by a heating machine (Model TKP-800, Tae Kwang Machine Co. LTD., Daegu, Korea) when the temperature lowered below 10 • C and 15 • C and overall the relative humidity (RH) was approximately within 40% to 60% in both greenhouses, respectively.
The appearance of insects and diseases was observed daily and the chemicals were periodically applied to control them once every two weeks following the manufacturer's instructions. We used Jijon (Farmhannong, Korea, Spiromesifen 20%) for the whitefly and the aphid with the dilution of 1:2000, Rampage (Hankooksamgong, Korea, Chlorofenapyr 5%,) for the thrip with the dilution of 1:1000, Butina (Farmhannong, Korea, Iminoctadine tris 30%) for the leaf mold and the gray mold, with the dilution of 1:1000, and Deomani (Farmhannong, Korea, Polyoxin B 50%) for the powderly mildew and the fungi with the dilution of 1:5000.

Data collection on Vegetative and Reproductive Growth
The vegetative parameters of plant height (PH), leaf length (LL), leaf width (LW), and stem diameter (SD) were measured using 70-d-old plants after transplanting (DAT) from five plants per accession in both greenhouses. The reproductive parameters of the number of flowers (NFL), the number of fruits (NFR), fruit set (FS), fruit yield (FY), marketable yield (MY), and output of marketable yield (OMY) were evaluated by calculating from the third to six trusses of each plant. Differences in FS and FY parameters between plants grown in LT and CT greenhouses were calculated by subtracting the index of FS and FY of CT from LT, respectively [16]. Fruit set (FS, %) with diameter ≥ 0.5 cm was calculated as follows [28]: Fruit set (%) = (The number of fruits / The number of flowers) × 100. Fruit yield (FY) was determined by the sum of fresh weight (FW) in kg of all fruits harvested from the third to sixth trusses from five individual plants.

Data Analysis
The significance of difference in vegetative parameters of PH, SD, LL, and LW, and reproductive parameters of NFL, NFR, FS, TY, MY, and OMY under LT and CT was assessed as described in the figure legends with Student's t-test and the analysis of correlation coefficients was performed among the total population (n = 35) for correlation coefficient with EXCEL 2016 (Microsoft, WA, USA). Principal components analysis (PCA) was implemented using SPSS (IBM SPSS v27.0., Chicago, IL, USA) to study the sample patterns and cluster analysis. The adequacy of the samples was tested by The Kasier-Meyer-Olkin (KMO) and the Bartlett's test of Sphericity (BTS) was conducted to estimate the relationship between variables.

The Analysis of the Vegetative Traits with Different Fruit Type
To study the vegetative traits including plant height (PH), stem diameter (SD), leaf length (LL), and leaf width (LW) of tomato plants under LT condition, we analyzed 35 tomato accessions with different fruit types classified into wild, cherry, medium, and large fruit size (Table S1). The PH of most tomato accessions in LT were remarkably reduced at 70 days after transplanting (DAT) compared to those in CT except for T32 accession in medium size ( Figure 1A), while SD in LT were not significantly different from those in CT except for T21 accession in cherry size ( Figure 1B). Next, the LL and LW were investigated at 70 DAT and 25 tomato accessions decreased under LT, whereas 10 tomato accessions including T04, T07, T09, T12, T13, T19, and T22 in cherry size, T28 and T31 in medium size, and T35 in large size were not significantly influenced by LT (Figure 2A,B).

The Analysis of the Reproductive Traits with Different Fruit Type
The reproductive traits including the number of flowers (NFL), the number of fruits (NFR), fruit set (FS), fruit yield (FY), and marketable yield (MY) were evaluated among 35 tomato accessions under LT and CT conditions. The NFL of T04 and T20 accessions in LT positively increased, whereas the NFL of T15, T29, and T31 accessions in LT negatively decreased compared to that in the CT condition ( Figure 3A). The NFR was also assessed among 35 tomato populations. Only the T11 accession in cherry type increased more in LT than that in CT, whereas the NFR of 2 genotypes (T01 and T02) in wild, 13 genotypes (T04,  T08, T09, T10, T14, T15, T16, T17, T18, T19, T20, T21, and T22) in cherry size, 4 genotypes (T27, T30, T31, and T33) in medium size, and 1 genotype (T35) in large size under LT were remarkably reduced compared to that in CT ( Figure 3B). On the other hand, the NFR of 7 genotypes (T03, T05, T06, T07, T12, and T13) in cherry type, 7 genotypes (T23, T24, T25, T26, T28, T29, and T32) in medium, and 1 genotype (T34) in large type were identified with no significant difference between LT and CT conditions ( Figure 3B). In addition, the fruit set (FS) ratio was calculated with all tomato accessions. Interestingly, the FS of only the T11 accession in cherry type was significantly higher in LT compared with that in CT ( Figure 3C). In contrast to this, the FS of 2 genotypes (T01 and T02) in wild, 15 genotypes (T04, T07, T08, T09, T10, T12, T14, T15, T16, T17, T18, T19, T20, and T21, T22) in cherry type, 3 genotypes (T27, T30, and T31) in medium type, and 2 genotypes (T34 and T35) in large type under LT dramatically decreased compared to those in CT ( Figure 3C). However, the FS of 4 genotypes (T03, T05, T06, and T13) in cherry type and 7 genotypes (T23, T24, T25, T26, T28, T29, and T32) in medium type were not observed with significant difference in both LT and CT conditions.

The Analysis of the Reproductive Traits with Different Fruit Type
The reproductive traits including the number of flowers (NFL), the number of fruits (NFR), fruit set (FS), fruit yield (FY), and marketable yield (MY) were evaluated among 35 tomato accessions under LT and CT conditions. The NFL of T04 and T20 accessions in LT positively increased, whereas the NFL of T15, T29, and T31 accessions in LT negatively decreased compared to that in the CT condition ( Figure 3A). The NFR was also assessed among 35 tomato populations. Only the T11 accession in cherry type increased more in LT than that in CT, whereas the NFR of 2 genotypes (T01 and T02) in wild, 13 genotypes (T04, T08, T09, T10, T14, T15, T16, T17, T18, T19, T20, T21, and T22) in cherry size, 4 genotypes (T27, T30, T31, and T33) in medium size, and 1 genotype (T35) in large size under LT were remarkably reduced compared to that in CT ( Figure 3B). On the other hand, the NFR of 7 genotypes (T03, T05, T06, T07, T12, and T13) in cherry type, 7 genotypes (T23, T24, T25, T26, T28, T29, and T32) in medium, and 1 genotype (T34) in large type were Leaf length and width were measured at 70 days. Significant differences were evaluated with Student's t-test (p ≤ 0.05, p ≤ 0.01, and p ≤ 0.001) and are denoted by *, **, and ***, respectively. NS indicates not significant and bars indicate ± standard deviation (n = 5).

The Principal Component Analysis (PCA) of Physiological Traits
In order to approach the multivariate analysis, the principal component analysis (PCA) was performed together with the correlation matrix of nine variables measured under LT conditions ( Figure S3). The Kaisere-Meyere-Olkin (KMO) for sampling adequacy on vegetative and reproductive score data was 0.643 in LT. Bartlett's Test of Sphericity (BTS) was significantly lower than 0.001. The first three principal components showed that the eigenvalues were greater than 1 (Table S2). Total variance of the data was explained with the 67.65% which were composed of 41.02% from component factor 1 (PC1) and 26.63% from component 2 (PC2), respectively, and the traits (scores > 0.30) were loaded onto PC1 and PC2 (Table S2). PC1 and PC2 were clearly separate from vegetative and reproductive parameters with positive trends ( Figure 5). PC1 was primarily combined with vegetative parameters (LL, LW, and SD) and PC2 was associated with reproductive parameters (FS, NFR, and FY).

Clustering Analysis
A score plot was obtained using the component factor 1 and factor 2 in LT, displaying 35 tomato accessions with different fruit types ( Figure 5). The factors were further subjected to hierarchical cluster analysis using the Euclidean distance matrix via Ward's method of agglomeration ( Figure 6). Based on the dendrogram result, four major groups were observed in LT. The cluster 2 consisted of vegetative parameters related to SD, LL, and LW and the cluster 3 was associated with reproductive parameters related to NFR, FS, and FY.

The Principal Component Analysis (PCA) of Physiological Traits
In order to approach the multivariate analysis, the principal component analysis (PCA) was performed together with the correlation matrix of nine variables measured under LT conditions ( Figure S3). The Kaisere-Meyere-Olkin (KMO) for sampling adequacy on vegetative and reproductive score data was 0.643 in LT. Bartlett's Test of Sphericity (BTS) was significantly lower than 0.001. The first three principal components showed that the eigenvalues were greater than 1 (Table S2). Total variance of the data was explained with the 67.65% which were composed of 41.02% from component factor 1 (PC1) and 26.63% from component 2 (PC2), respectively, and the traits (scores > 0.30) were  (Table S2). PC1 and PC2 were clearly separate from vegetative and reproductive parameters with positive trends ( Figure 5). PC1 was primarily combined with vegetative parameters (LL, LW, and SD) and PC2 was associated with reproductive parameters (FS, NFR, and FY).

Clustering Analysis
A score plot was obtained using the component factor 1 and factor 2 in LT, displaying 35 tomato accessions with different fruit types ( Figure 5). The factors were further subjected to hierarchical cluster analysis using the Euclidean distance matrix via Ward's method of agglomeration ( Figure 6). Based on the dendrogram result, four major groups were observed in LT. The cluster 2 consisted of vegetative parameters related to SD, LL, and LW and the cluster 3 was associated with reproductive parameters related to NFR, FS, and FY.

Discussion
Tomato plants have evolved against adverse environment factors and current tomato plants harbor several mechanisms to overcome cold stress [29][30][31][32][33]. Previous studies have mainly focused on the response to low temperature stress with limited genotypes and fruit types as well as a short period of the treatment during specific growth stages

Discussion
Tomato plants have evolved against adverse environment factors and current tomato plants harbor several mechanisms to overcome cold stress [29][30][31][32][33]. Previous studies have mainly focused on the response to low temperature stress with limited genotypes and fruit types as well as a short period of the treatment during specific growth stages [11,17,25,34]. Moreover, recent studies have not been reported in regards to night low temperature (NLT) which can save the heating costs of the tomato cultivation in winter greenhouse. In this study, we assessed the vegetative and reproductive traits of 35 tomato accessions with different fruit types during whole plant growth stages under NLT.
LT noticeably influenced vegetative traits including PH, LL, and LW in tomato plants at 70 DAT compared to CT. In particular, the effect of LT on the PH in LT condition resulted in a significant decrease in most accessions except for one accession, T32, for which there was no significant difference between the LT and CT conditions ( Figure 1A). Previous studies have shown that the growth of tomato plants is significantly retarded in LT conditions as shown in this study [3,21,22,25]. Moreover, the researchers showed that the inhibited growth of PH results from a decrease in the number of leaves [35,36]. It is of interest to further investigate the correlation of PH and the number of tomato leaves between LT-tolerant and-sensitive cultivars. The effect of LT on LL and LW in LT exhibits that the growth inhibited 100% and 50% of accession in wild, 60% and 40% accession in cherry types, and 50% accession and 0% accession in large types, respectively ( Figure 2). In line with our results, previous studies have also demonstrated that LT-treated plants lead to a reduction in leaf area and the decreased leaf area ratio which influences the decrease in relative growth rate by lowering photosynthesis [21,22,37]. Next, our important endeavor is to determine how photosynthetic parameters including chlorophyll contents and gas exchange rates are associated with LL and LW in LT-tolerant and -sensitive cultivars within the same fruit types or different fruit types. Intriguingly, our result shows that LT did not affect the SD of most accessions ( Figure 1B). The result is inconsistent with a previous study which has determined the difference in the SD of tomato plants under LT and CT conditions. On the other hand, a report supported that a minimum temperature around −2 • C to −4 • C results in no stem diameter difference in woody plants, whereas a minimum temperature around −5 • C to −6 • C causes the reduction of stem diameter [38]. The minimum temperature in our previous research was lower by around 3 • C until 10 weeks than that in this study [16]. Studies on SD of tomato plants have not been reported under LT conditions. Further studies are necessary to determine low temperature set points affecting the stem diameter in tomato plants.
We evaluated the reproductive parameters including NFL, NFR, and FS in LT and CT. Previous reports have determined the increased NFL in declining average temperature during 24 h [20], day temperature [19], and night temperature [39]. However, a study revealed that NFL in the first truss was not influenced by air temperature, whereas the NFL was increased by low root temperature [40]. Moreover, our previous study exhibits that the increased NFL and decreased NFL result from a genotype-specific interaction with LT [16]. Our current result shows that LT did not affect NFL in most accessions ( Figure 3A), except that the NFL of two and three accessions was increased and decreased under LT, respectively ( Figure 3A). Although we do not know the exact mechanistic function of the effect of LT on NFL, our finding with previous results imply that the effect of low air temperature on NFL was not notably influenced in our study. Further studies are required to determine how air and root low temperature impact on NFL with fruit types. Previous studies have shown that NFR and FS play an important role in determining LT tolerant-tomato cultivars, because NFR and FS are highly related to FY [16,35]. Similarly, our results show that the effect of LT on NFR led to a reduction in 70.0% of cherry type and 36.36% of medium type, respectively ( Figure 3B). As well, FS also led to a reduction in 80% of cherry type and 36.36% of medium type in LT ( Figure 3C), showing that NFR is closely correlated with the FS. Moreover, our findings reveal that both NFR and FS are associated with FY (Supplementary Figure S3). Considering that the NFLs were not different in most accessions regardless of fruit types, the poor quality of flower pollens in LT might influence the decreased FS among the accessions [12,18,20], but the reduced FS might be not influenced by ovule viability and/or stigma level [41]. Furthermore, the effect of LT resulted in the significantly declined FY and MY, regardless of fruit type (Figure 4). A previous research has demonstrated that the effect of LT reduces FY due to the reduced number of pollens, the reduced pollen germination rate, aborted pollens, and the malformation of pollen tube which affect FS later on [18,20]. Although we cannot completely rule out the involvement of ovule and stigma development in FS and FY, next, it is worth determining the mechanistic role of pollen germination and pollen grain viability involved in fruit set and fruit development in response to LT. Collectively, our aforementioned results suggest that the impact of LT on vegetative and reproductive traits may be associated with the effect of genotype, not with fruit types.
Previous researches have been determined the traits during vegetative or reproductive stage in LT and the question of how the vegetative growth is associated with fruit yield remains unanswered. We applied the correlation matrix of vegetative and reproductive traits to the PCA analysis [42]. The first two PCA analyses explained 41.02% and 26.63% of the total physiological variables among 9 traits determined in LT. Six traits (with score > 0.30) were loaded onto PC1 and PC2 (Table S2). Our finding showed that the angle between the vectors of traits including LL, LW, and SD in PC1 and NFR, FY, and FS in PC2 are lower than 90 • . The vegetative and reproductive traits are positively correlated to some extent and the FY is strongly correlated with NFR and FS. In line with our results, previous researches have shown that the effect of LT on reproductive traits led to the reduction in NFR and FS, together with FY, indicating that these traits are positively associated with LT [20]. Moreover, a study reported that seed germination and vegetative growth of tomato are not correlated under cold stress, which suggests that each growth and development stage should be selected and assessed for plant breeding [11]. It is therefore likely that FY in LT is not strongly connected with vegetative traits. Subsequently, the biplot analysis was further drawn to assess the multivariate relationship among 35 tomato accessions and it displayed that 35 tomato genotypes are not grouped with fruit types ( Figure 5). Rather, the accessions are sporadically located in a plot regardless of fruit types. This might result from the tomato accessions that possess the genetic diversity within the same fruit type [43]. Moreover, the cluster analysis exhibited four major clusters ( Figure 6) and the cluster 1 and 4 tended to be negatively involved in vegetative and reproductive traits and/or be less associated with them on positive trends. The cluster 2 and 3 showed vegetative (SD, LL, and LW) and reproductive parameters (NFR, FS, and FY) with positive directions, respectively. Interestingly, 5 out of 6 were cherry type in cluster 3, but different fruit types were grouped in cluster 2, implying that cherry types in cluster 3, with the consideration of reproductive parameters, are more tolerant to LT than other accessions. Although we cannot completely exclude a possibility that the results from the PCA and clustering analysis might be influenced by some secondary symptoms of LT stress response, the selected genotypes in cluster 2 and 3 could be further utilized for the breeding program to select LT-tolerant tomatoes with the parameters in the greenhouse.

Conclusions
The current study has determined the physiological traits of thirty-five tomato accessions with different fruit types in the response to night low temperature, which is economically important for the tomato cultivation in winter greenhouse. Based on correlation coefficient, PCA, and cluster analysis, some accessions were closely involved in vegetative and/or reproductive parameters depending on genotypes. Future researches will be required to evaluate more accessions of large fruit type in LT and the selected accessions will have to be used for the determination of physiological and molecular functions, combined with DNA-and RNA-seq.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/10 .3390/agriculture11080792/s1, Figure S1: Air temperature was measured in LT and CT greenhouse during the period of tomato growth and development, respectively. Figure S2: The analysis of difference in FS and FY among 35 tomato accessions in LT and CT. Figure S3: The correlations coefficients between vegetative and reproductive traits in total population of tomatoes in LT, Table  S1: Tomato accessions for the evaluation of physiological traits against low temperature in winter 2020-2021. Table S2: Loading matrix associated with the principal components analysis (PCA) for 9 physiological traits.