Next Article in Journal
Foliar Application of Ulva rigida Water Extracts Improves Salinity Tolerance in Wheat (Triticum durum L.)
Previous Article in Journal
Genome-Wide Identification of the NHX Gene Family in Punica granatum L. and Their Expressional Patterns under Salt Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Correlation Analysis of High-Throughput Fruit Phenomics and Biochemical Profiles in Native Peppers (Capsicum spp.) from the Primary Center of Diversification

by
Jorge González-López
,
Simón Rodríguez-Moar
and
Cristina Silvar
*
Grupo de Investigación en Bioloxía Evolutiva, Departamento de Bioloxía, Universidade da Coruña, 15071 A Coruña, Spain
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(2), 262; https://doi.org/10.3390/agronomy11020262
Submission received: 18 November 2020 / Revised: 23 January 2021 / Accepted: 27 January 2021 / Published: 30 January 2021
(This article belongs to the Section Horticultural and Floricultural Crops)

Abstract

:
The main goal of this work was to investigate the relationship between the fruit morphology and biochemical composition of peppers (Capsicum spp.). For that purpose, one hundred native varieties from the Andean region, where the genus Capsicum has its origin, were analysed for different phytochemical compounds. In addition, pepper fruits were assessed with the highly precise phenomics tool Tomato Analyzer. The collection showed a broad variability which was more evident within the C. annuum group. On average, C. frutescens accessions displayed the highest levels of solid soluble content, pH, polyphenols and antioxidant activity. The Tomato Analyzer descriptors under the categories of size, shape index, and latitudinal section, mostly contributed to the variance among Capsicum groups. C. annuum hold the larger fruits, whereas C. frutescens comprised fruits of smaller sizes. The correlation analysis revealed that biochemical traits were negatively correlated with the fruit parameters related to size, suggesting that huger fruits contain lower amounts of chemical metabolites. The multivariate approximations demonstrated that Andean peppers assorted according to morphometric and colorimetric characteristics, but independently of their species or geographical origin. Groups of valuable native varieties carrying promising traits were identified.

1. Introduction

The genus Capsicum belongs to the Solanaceae family and currently harbors more than 35 species, five of which (C. annuum L., C. chinense Jacq., C. frutescens L., C. baccatum L., and C. pubescens Ruiz et Pav.) are domesticated [1]. Capsicum is originates from South America and its genesis apparently occurred in a broad area including Peru, Ecuador and Colombia, along the Andes in nort-western South America [2]. According to Carrizo García et al. [2], the subsequent expansion of the genus took place following around the Amazon basin, towards central and south-eastern Brazil, coming back to western South America, and finally northwards to Central America. Therefore, the western Andes of South America remained as a significant center of speciation and northward species dispersal, primarily for C. annuum [2,3]. After the voyages of exploration to the New World, peppers moved worldwide, suffering from additional diversification at the secondary centers. These phenomena resulted in the vast phenotypic variability that can be observed nowadays in the form of native varieties or landraces, well adapted to specific agro-climatic conditions and consumers habits [1]. C. annuum was the most successful in this expansion, giving rise to sweet and hot pepper varieties spread all over the world. However, C. pubescens and C. baccatum remained mostly restricted to Andean South America and the Central American highlands, while C. frutescens and C. chinense are mainly cultivated as spice crops in Africa, Asia and South America [1]. Consequently, landraces exist worldwide and they likely maintain combinations of traits that were left behind when growers selected modern cultivars. In this context, it is foreseeable that native varieties from the primary centers of Capsicum speciation, i.e., Ecuador, Peru and Bolivia, preserved a higher portion of the wealth stored on primitive peppers, and therefore constitute a promising resource for the breeding of new market types with added value.
Peppers are one of the most important vegetables worldwide due to their versatility for cuisine, medicine and industry [4]. The five cultivated Capsicum species are today grown globally both as a spice and as a vegetable crop. Sweet, bell peppers have more presence in Europe and North America, while pungent varieties currently dominate the markets in South-Eastern Asia, South America, and Africa [5]. For consumption as a vegetable, most of the pepper diversity was addressed to develop fruits with a sweet taste, large size, blocky shape, thick pericarp and bright colors. In spicy chilli peppers, the predominant selected characters are hotness, fruit size and skin thickness [6]. For industrial food purposes, pepper should be dehydrated and powered, so that water content represents an essential parameter. In turn, pepper powder must possess the adequate spiciness and a specific color [7]. On the other hand, consumers are becoming increasingly interested in tastier vegetables containing healthier components and higher nutritional quality [8]. Many studies have been performed to unravel the biochemical composition of pepper fruits and this crop has demonstrated to be an important source of bioactive compounds, which possess analgesic, anti-obesity, cardio-protective, pharmacological, neurological and dietetic properties [9]. Nevertheless, the search for novel pepper genotypes with tastier fruits, more attractive colours, containing the required pungency and high concentrations of health-promoting components still remains relevant. This might help breeding programs to improve the existing commercial types, expanding the opportunities for new market segments and increasing the profitability of producers.
The organoleptic and bioactive properties of pepper fruits are largely due to the presence of phytochemical compounds. To a large extent, sugars and acids will determine the sensory quality of fruits connected to flavor, while other compounds such as carotenoids and polyphenols are responsible for the visual aspect and health-related attributes of peppers. Many of the functional features of these berries are largely attributed to the antioxidant activity of different metabolites, which play a crucial role in scavenging free radicals responsible for oxidative damage [10]. All these attributes will play a key influence on how consumers perceive the quality of the product and decide their preferences. Apart from adequate characterization of compositional profiles, an extensive and detailed assessment of fruit phenomics will also be essential since fruit size, shape and weight will largely influence the establishment of suitable market niches and the good acceptance of new pepper varieties [11]. Among the lately developed high-throughput phenotyping tools, the Tomato Analyzer (TA) software has emerged as an extremely useful system to perform accurate phenomics assessment in Solanaceae. This freeware allows semi-automatic and highly accurate scoring of a large number of quantitative traits from scanned images of fruit sections [12,13]. Although this software was initially conceived for tomato [14,15,16], it has been successfully employed to evaluate fruits of eggplant [17,18], and pepper [19,20,21,22,23]. Diversity measured by fruit TA descriptors demonstrated to be higher than that reflected by conventional descriptors. This is due to the fact that automated phenotyping covers a broader dissection of fruit architecture, allowing to measure morphometric attributes which are nearly impossible to obtain manually [16,21,24].
Several works focused on phenotypic characterization of diverse Capsicum spp. collections and their biochemical compositions have been conducted to date. However, both assessments are usually accomplished separately. As far as we know, few reports have attempted to establish some correlations between fruit morphology and chemical metabolites and this was never the primary aim of those studies [25,26,27]. Nevertheless, this knowledge is fundamental, as it will strongly determine the speed at which promising traits could be integrated into commercial pepper types. In a recent work, we employed a small collection of Ecuadorian peppers to draw relationships between the morphometrics of pepper fruits and their content in bioactive compounds [23]. Herein, we analyzed the biochemical composition of a wider panel of peppers (Capsicum spp.), and accurately determined the fruit phenomics by using a high-throughput digital assessment. Linkages between both datasets were determined. For that purpose, a broad collection of native varieties from north-western South America, which likely constitutes the origin of Capsicum, was used.

2. Materials and Methods

2.1. Plant Material

One hundred accessions depicting domesticated species and their ancestors were selected for this work: C. annuum (22), C. annuum var. glabriusculum (1), C. annuum wild (4), C. chinense (25), C. frutescens (12), C. pubescens (13), C. baccatum (4), C. baccatum var. pendulum (16), and C. baccatum var. baccatum (3). Forty-three were from Ecuador, 28 from Peru, and 18 from Bolivia. Eleven C. annuum accessions from Mexico were also included for comparative purposes (Table S1). Accessions were kindly provided by the USDA-ARS Plant Genetic Resources Conservation Unit (USA), the Institute for Conservation and Improvement of Valencian Agro-diversity (COMAV, Spain), the Centre for Genetic Resources (CGN, The Netherlands), and the Leibniz Institute of Plant Genetics and Crop Plant Research (IPK, Germany). The accession PM647 was kindly supplied by Dr. Alain Palloix (rest in peace) from the French National Institute for Agriculture, Food, and Environment (INRAE, France). When possible, they were selected to cover broad geographical regions within each country. Data on taxonomical classification and population type of each accession were extracted from genebanks databases. This collection was previously characterized for genetic diversity and molecular markers linked to disease resistance and pungency-related traits [28,29,30]. C. annuum accessions from the collection were split into pungent and non-pungent fruits according to the results of marker MAP1, linked to the Pun1 gene [31], as previously described [30].
The accessions were cultivated in a greenhouse during 2017 at the Centro de Investigaciones Agrarias de Mabegondo (Mabegondo, A Coruña, Spain) (43°15′ N, 8°18′ W). Four plants per accession were grown following a completely randomized design. Plants were drip irrigated and fertilized with a mix of nitrogen, phosphorus, and potassium before and after transplanting. Phytosanitary treatments against whiteflies, aphids, and spider mites were applied when necessary.

2.2. Conventional Descriptors Assessment

Accessions were characterized by using six conventional morphological descriptors for Capsicum following the protocol described by the International Plant Genetic Resources Institute (IPGRI) [32]. These comprised plant height (PHE), plant width (PWI), stem length (SLE), stem diameter (SDI), fruit weight (FWE) and fruit pedicel length (FPL). All parameters were meassured in centimeters (cm), except FEW, which was expressed in grams (g). Between ten and twenty-five fruits per accession, depending on their size, were harvested at maturity stage and employed for conventional characterization in the laboratory.

2.3. Fruit Characterization with Tomato Analyzer

Mature fruits, from 10 to 25 per accession, were subjected to digital phenotyping. Pepper fruits were longitudinally and transversally cut and scanned with an HP Scanjet G3110 photo scanner (Hewlett-Packard, San Jose, CA, USA) at a resolution of 300 dpi. Stored images (TIF format) were subsequently analyzed using Tomato Analyzer version 4 software [12,13]. A total of 49 fruit parameters, categorized into: basic measurements (7), fruit shape index (3), blockiness (3), homogeneity (3), proximal fruit end shape (4), distal fruit end shape (4), asymmetry (6), internal eccentricity (5), latitudinal section (5), and average color values (9) were automatically recorded (Table S2). A complete description of morphometric and colorimetric descriptors can be found elsewhere [12,13,14,33]. Default settings were used for all categories, although points were adjusted manually when the software was unable to accurately identify the outline of a trait. Individual measures of each fruit were used to obtain an average value for the corresponding accession.

2.4. Biochemical Methods

Ripe fruits were collected from different plants of the same accession and pooled into two bulks. Each group was homogenized with a blender, and filtered through a gauze. The resulting pepper juice was used for analyses. All analytical measures were performed in triplicate.

2.4.1. Soluble Solid Content, pH, and Titratable Acidity

The soluble solid content (SSC) was determined using a hand refractometer (Shibuya Optical Co., Ltd., Saitama, Japan) and expressed as equivalent °Brix. The pH was measured by direct reading using a pH meter (Crison micropH 2.000, Barcelona, Spain). Titratable acidity was calculated by titrating 25 mL pepper juice against 0.1 N NaOH, to an end point of pH 8.2, as indicated by a pH meter. Titratable acidity was expressed as percentage of acidity, using malic acid as an equivalent. Fresh pepper juice was used directly for these analyses.

2.4.2. Determination of Total Polyphenol Content

The Folin–Ciocalteu method was applied according to Singleton and Rossi [34] with slight modifications. Briefly, 1 mL on pepper juice was centrifuged at 13,000× g rpm for 5 min. Then, 50 μL supernatant was diluted with 750 μL distilled water and mixed with 50 μL Folin-Ciocalteau reagent. After an incubation time of 3 min, 150 μL of 20% Na2CO3 was added. The mixture was shaken vigorously and incubated for 2 h at room temperature in the dark. Absorbance against a blank was measured at 760 nm using a Thermo Heλios spectrophotometer. The standard curve was constructed with gallic acid. The results were expressed as gallic acid equivalents (GAE) in µg/g fresh weight (FW).

2.4.3. Trolox Equivalent Antioxidant Capacity (TEAC)

The procedure described by Miller and Rice-Evans [35], with slight modifications, was employed. Briefly, 1 mL on pepper juice was centrifuged at 13,000× g rpm for 5 min. Then, 20 μL supernatant was mixed with 980 μL of the diluted 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS)-radical solution. After an incubation time of 15 min at room temperature in darkness, the absorbance was read at 734 nm. The ABTS radical stock solution was prepared by combining 50 mL of 7 mM ABTS and 50 mL of 2.45 mM K2O8S2. The mixture was place in the dark at room temperature for 16 h to generate the radical. One milliliter of ABTS stock solution was diluted with approximately 50 mL water and the absorbance was adjusted to 0.9 at 734 nm before use. Trolox was employed for external calibration and results were expressed as µg Trollox/g FW.

2.4.4. Determination of Moisture Content and Dry Mass

Two bulks of ripe fruits per each accession comprising a variable number between 10 and 25 were independently subjected to freeze-drying in a Telstar-Cryodos-80 lyophilizer (Telstar, Barcelona, Spain) for 24–48 h. Moisture content was calculated as the difference in the sample mass before and after drying, expressed as percentage. Fruit dry weight (FDW) was determined by dividing the total dry matter by the number of fruits at each bulk.

2.5. Data Analysis

Analysis of variance (ANOVA) was performed for conventional descriptors, biochemical and digital parameters. Means, standard deviations, and coefficient of variation (expressed in percentage as the ratio between standard deviation and mean) were used for descriptive analysis of traits. Significant differences among means were detected using Tukey’s HSD test (p < 0.05). Statistical analyses were performed with Statistical Product and Service Solutions (SPSS) version 17.0 software [36]. Correlations among all evaluated parameters were estimated by using Pearson’s test at p < 0.05, after Bonferroni’s adjustment for multiple comparisons [37]. The calculation of coefficients and visualization of correlograms were performed with Rcmdr and Corrplot packages implemented in R 4.0.0 software [38]. Conventional, TA and biochemical traits were comprehensively examined through a principal component analysis (PCA). The similarity across accessions was estimated by agglomerative hierarchical cluster analysis (HCA) using the Ward coefficient. In addition, a Principal Coordinate Analysis (PCoA) for each Capsicum spp. was performed. Multivariate analyses were conducted with the computer software SPSS version 17.0 and R 4.0.0 (packages Stats, Vegan, and Ape). Graphical representation of the tree was performed with MEGA X software [39].

3. Results

3.1. Biochemical Analysis

Seven parameters were quantified on eighty-two pepper accessions (Table S3). The ANOVA found highly significant differences (p < 0.001) among them, although the greatest part of the variance was explained by fruit dry weight and antioxidant activity. The FDW ranged from 0.08 g to 6.63 g, with the lowest values being reached by CGN20808 and CAP524, and the highest by C87, all belonging to C. annuum. The average value for Trolox equivalent antioxidant capacity (TEAC) was 615.09 µg Trolox/g FW, varying from 84.86 µg Trolox/g FW in BGV005890 (C. chinense) to 4126.63 µg Trolox/g FW in the C. annuum accession CAP524. Values for the pH ranged from 4.7 in PI543184 to 6.3 in BGV014749, both C. chinense accessions. Total polyphenols content showed an average of 916.36 µg GAE/g FW, varying from 286.06 µg GAE/g FW in BGV005890 (C. chinense) to 3706.32 µg GAE/g FW in PI585257 (C. frutescens). The C. annuum accession CAP524 showed the highest values for antioxidant capacity (4126.62 µg Trolox/g FW), titratable acidity (0.57%) and SSC (11.98 °Brix), while the C. chinense accession BGV005890 exhibited the lowest values for those parameters, i.e., TEAC of 84.86 µg Trolox/g FW, acidity of 0.06% and SSC of 3.25 °Brix (Table S1). The average value for the moisture content was 82.06%, the highest being observed for the C. annuum accession C87 (90.29%), while the lowest (54.85%) for C. baccatum CGN17042.
The analysis of variance, considering the species effect, indicated significant differences (p < 0.001) among Capsicum spp. for all analytical measures, except for titratable acidity (Figure 1). C. frutescens showed the highest significant values for total phenols (3117.32 µg GAE/g FW) and antioxidant capacity (2495.44 µg Trolox/g FW). C. frutescens and C. annuum possessed the highest amounts of SSC, at 7.69 and 6.55 °Brix, respectively, while the greatest values for pH were recorded in C. frutescens (5.62) and C. chinense (5.38). The greatest percentages of moisture content were observed in C. pubescens (86.31%) and C. chinense (84.76%), and the highest FDW was recorded for C. annuum accessions (average value of 1.90 g) (Figure 1). The broadest variability was observed within the C. annuum group, which exhibited highest coefficients of variation for all traits except pH.

3.2. Morphological Assessment with Conventional and TA Descriptors

A great diversity was observed within the Andean collection after the evaluation with six quantitative conventional descriptors. The parameters which most contributed to explain the variance among genotypes were PHE, FWE, and FPL (Table S4). Thus, the PHE showed an average value of 110.36 cm, ranging from 35.67 cm in C. annuum accession C87 to 180.50 cm in C. baccatum var. pendulum PI257133. The weightiest fruits were recorded for C87 (C. annuum) with an average value of 47.89 g per fruit, while the weightless ones (0.20 g per fruit) were found in C. annuum var. glabriusculum PM647. The C. chinense accession CAP472 possessed the shortest FPL (2.17 cm) whereas the longest were observed for C. baccatum CAP1366 (9.54 cm). A coefficient of variation up to 97.92% was recorded for few, while for the remaining traits was lower than 40% (Table S4).
Belonging to a particular Capsicum species certainly determined the morphological variation among accessions. Thus, significant differences among the five species were identified for all traits (Figure 2). C. baccatum and C. frutescens accessions showed the significantly highest plant heights and stem lengths, while the shortest were observed for C. pubescens. A difference of up to 47 g was observed for FWE between the smallest C. frutescens and the biggest C. annuum fruits (Figure 2).
The Tomato Analyzer software was employed for the assessment of fruit morphology by image examination of scanned fruit sections Highly significant differences (p < 0.001) were detected among average values, except for proximal angle micro (PAMI), which exhibited differences at p < 0.05 (Table S5). Those attributes related to fruit size and shape, i.e., those under the categories of basic measurements, shape index, and latitudinal section, primarly explained the variance. Minimum values of zero were recorded for proximal indentation area (PIA), distal indentation area (DIA), obovoid (Ob), and H. Asymmetry. Ob (HAob), The coefficient of variation ranged from 1.89% (proximal eccen-tricity, PEC) to 263.39% (HAob). For each category, the descriptors with largest variability were: Area (75.30%), Fruit Shape Index Internal (48.02%), Distal End Protusion (86.43%), Tomato Pericarp Area (72.23%), Curved Fruit Shape Index (52.84%), Fruit Shape Triangle (32.72%), circular (38.81%), Proximal Indentation Area (116.14%), H. Asymmetry. Ob (263.39%), and Average a Value (AaV) (42.36%) (Table S5). The color analysis demonstrated that all color parameters displayed distinctive variations among accessions. RGB and CIELab colur spaces indicated larger amounts of redness and yellowness components. The Hue value ranged from 45.91° to 102.87° suggesting that the colors in the collection varied from red-orange to yellow-green. The majority of fruits showed moderate lightness (average luminosity value (ALV) = 50.10) and saturation of colors (average chroma (Ach) = 54.16) (Table S5).
As occurred with conventional traits, broad variability was found among the five domesticated species. Forty-four out of 49 TA descriptors exhibited highly significant differences (p < 0.001) while divergences at p < 0.01 were observed for HAob and tomato pericarp thickness ratio (TPTR). No significant differences were recorded for PAMI, DIA, and PEC, (Table 1). In general, C. frutescens and C. annuum showed the lowest and highest values, respectively, for the majority of fruit size features, including those under the basic measurements and latitudinal section categories. Considering the shape, measured as indexes fruit shape index external I (FSIEI), fruit shape index external II (FSIEII), FSI, and fruit shape index internal (FSII), the highest and significantly different values were recorded for C. annuum and C. baccatum (Table 1). This suggest that pepper fruits from those species displayed much larger height than width, whereas in C. frutescens, C. chinense and C. pubescens fruits, both attributes were not so dissimilar. Accessions from C. pubescens displayed the most circular and ellipsoid fruits, while the most rectangular were exhibited by C. annuum. The largest variability, according to the coefficient of variation, was recorded for C. annuum, which displayed CV higher than 50% for twenty-three traits. Luminosity and average luminosity values indicated that fruits from C. chinense accessions possess the brightest colors, while those from C. pubescens were the darkest. C. pubescens fruits displayed the highest amounts of red colour (higher AaV), while C. chinense fruits tended to be the most yellowish (higher Average b Value (AbV)). All species were included in the orange Hue angle range. Values for Chroma did not differ greatly, although significantly more vivid external colors were observed in C. frutescens fruits when compared to C. annuum (Table 1).

3.3. Divergences among Geographical Regions

Differences among Ecuador, Peru, Bolivia and Mexico were analyzed for C. annuum, C. baccatum, C. chinense, and C. pubescens, which have a representative number of accessions at each country. In the C. annuum group significant differences were only found for SSC, which was higher in peppers from Mexico than in those from Ecuador (Figure 3). C. chinense accessions from different countries significantly varied for their titratable acidity and moisture content, with peppers from Peru the ones with the highest acidity and lowest amounts of water. Significant differences for all biochemical traits, except for pH, were observed across C. baccatum accessions, with the Bolivian fruits possessing the highest levels of SSC, titratable acidity, polyphenols and antioxidant capacity but, the lowest values for moisture content and dry weight. The biochemical parameters quantified on C. pubescens accessions from Ecuador were not significantly distinctive from the homonym peppers growing in Peru (Figure 3).
Significant differences were observed in C. annuum for the conventional descriptors PHE, SLE, SDI, FWE, and FPL and 35 out of 49 TA descriptors. Thus, accessions from Mexico displayed the plants with the highest height and stem length, while the weightiest fruits appeared in Peruvian and Ecuadorian peppers (Table S6). Likewise, C. annuum peppers from Peru and Ecuador exhibited the largest values for the majority of parameters under the size and latitudinal section categories, although. Ecuadorian fruits presented more elongated and rectangular shapes, while peppers from Peru arise as more ellipsoid and circular, as demonstrated by the values of shape index, homogeneity, blockiness and asymmetry. Plant architecture did not significantly vary among C. baccatum accessions from the different Andean regions, although Ecuadorian peppers displayed the weightiest fruits (Table S6). Data from scanned fruits revealed significant differences for the majority of attributes, with C. baccatum peppers from Bolivia showing the lowest significant values for size, shape, and latitudinal sections. The conventional parameters PHE, SLE, FWE, and FPL significantly diverged among geographical regions in C. chinense. Thus, accessions from Ecuador exhibited taller plants than those from Bolivia. Likewise, pepper fruits from the former country had the greatest weight and largest pedicels (Table S6). The TA analysis denoted that C. chinense accessions from Ecuador possessed significantly higher values for those parameters associated to fruit size and latitudinal section, although the fruits shape indexes resulted similar across Andean countries. C. pubescens plants from Ecuador were significantly taller than those from Peru. On the contrary, Peruvian accessions showed the heaviest fruits. TA analysis revealed that C. pubescens peppers from Peru were larger, although with more roundish shapes. For all four species, the average color values denote that more reddish and darker peppers are prevalent in Ecuador, while Peruvian and Bolivian accessions contained lighter fruits with higher amounts of yellow (Table S6).

3.4. Variations between Pungent and Non-Pungent Accessions

All biochemical traits, except for pH, significantly differed among pungent and non-pungent peppers. Putative hot berries possessed greater values for SSC, titratable acidity, polyphenols and antioxidant capacity. However, sweet fruits displayed a higher moisture content and dry weight (Figure 4). Assessment of conventional attributes revealed that pungent types were significantly associated with taller plants and lighter fruits. For TA descriptors, highly significant differences (p < 0.001) were observed for 20 parameters and significant differences at p > 0.05 were recorded for 11 traits. On average, pungent peppers had smaller fruits, with thinner pericarps and more roundish shapes (Table 2). Broader variability, expressed as coefficients of variation for each trait, was found within the pungent group. Similarly, a greater variation in CIELab coordinates was evidenced in potentially hot peppers. Significant differences were found for AaV, average hue (AHue) and AChV, suggesting that sweet accessions possessed slightly more redness and brightness fruits (Table 2).

3.5. Correlation between Biochemical and Morphometric Parameters

Strong and positive significant correlations were recorded between pairs for SSC, titratable acidity, polyphenols and TEAC. These parameters negatively correlated to the moisture content and dry weight, although only the former resulted highly significant (Figure 5). Remarkably, the mean values of SSC, titratable acidity, polyphenols and TEAC showed a particularly stronger negative correlation with the TA descriptors under the categories of basic measurements and latitudinal section (tomato pericarp area (TPA) and thickness (TPT)). However, after Bonferroni correction, only width mid-height (WMH), maximum weight (MW), TPA and TPT remained significant. Conversely, the moisture content and dry weight were positively and significantly inter-related to the same TA parameters (Figure 5). Such results indicate that large fruits possess higher amounts of water, but a reduced levels of bioactive compounds. The network of correlations for TA parameters revealed that some categories were tightly linked and significantly correlated, whereas most of them remained rather independent. Hence, the parameters belonging to the size category were positive and significantly correlated with TPA, TPT, ellipsoid, circular, V. asymmetry (VAs), H. Asymmetry. Ovoid (Ov) (HAov), and the shape index attributes FSIEI, FSIEII, and curved fruit shape index (CFSI). The homogeneity of pepper fruits displayed significant associations with the internal eccentricity category and shape indexes (FSIEI, FSIEII, and CFSI), even though ellipsoid and circular correlated in a positive manner, while rectangular did it in the opposite way. Fruit weight was positive and significantly correlated to all TA parameters describing size and pericarp features, but also to moisture content and dry weight. On the contrary, a negative connection was established between that parameter and biochemical traits, although only the comparison with titratable acidity was significant. Those results suggest that weighter fruits hold thicker walls but lower amounts of bioactive metabolites. The parameters assigned to the color section were the most independent, exhibiting lack of correlation with the rest of the TA attributes, and being only significantly inter-connected among themsleves. Interestingly, the entire coefficients between pairs were positive, except for AaV, which exhibited negative correlations to the other colorimetric traits (Figure 5).

3.6. Multivariate Analysis

The principal component analysis resulted in ten principal components with eigenvalues > 1, cumulatively accounting for 89.08% of the total variance. The first component explained 28.3% of the total variance and it was positively and robustly correlated (>80%) to TA attributes related to shape, such as FSIEI, FSIEII, CFSI, circular, FSII, and lobedness degree (LD) (Figure S1). The second component, which accounted for 23.4% of the variance, was positively correlated to traits describing fruit size (perimeter (P), area (A), MW, HMW, and TPA), and weight (FEW and FDW). The third component contributed for 10.2% to the total variance, and the color parameters were clearly the principal traits responsible for the observed variability, with average green (AG), ALV, and AHue being those with a higher degree of positive correlation (Figure S1). Biochemical traits highly contribute to the fifth component, which only explained 5.1% of the variability detected across Andean peppers.
The three principal components were used to project the pepper accessions on two-dimensional plots. Pepper fruits were widely dispersed on the PCA diagrams according to their size, morphology and color (Figure 6). FSII, FSIEI and FSIEII mostly explained the distribution of the accessions on the first axis (PC1). On the other hand, Area and TPA were the main factors discriminating the pepper fruits on the second component (PC2). The arrangement of accessions on the third axis (PC3), primarily responded to the fruit colors. Belonging to a particular Capsicum spp. did not predominantly influenced the distribution of pepper accessions on the graphical space, although certain tendencies could be observed. C. annuum peppers spread all around the plot, showing the greatest diversity, which encompass fruits with variable shape and size. Hence, morphologies oscillated from highly elongated fruits with large height/width ratio (FSIEII = 7.39) to more roundish fruits (FSIEII = 1.18). Likewise, size varied from huge fruits with areas around 3000 mm2 and weights of ca. 50 g, to very small fruits (area = 82.6 mm, FEW = 0.20 g). C. baccatum and C. chinense also exhibited some level of dispersion, although the former tend to concentrate on the positive axis of PC1, while the latter mainly converged on the negative part. This indicates that C. baccatum species harbor fruits with more elongated shapes, while round and conical peppers are predominant within the C. chinense group (Figure 6). The fruit magnitude was more variable across C. baccatum accessions, whose size and weights ranged from 60.21–1985.48 mm2 to 0.24–16.75 g, respectively. Accessions from C. pubescens primarily plotted on the quadrant defined by negative PC1 and positive PC2 axes, denoting that those peppers possess medium sizes and low rates for the height/width indexes. Finally, C. frutescens comprised the smallest fruits (average Area = 300 mm2, average FWE = 1.41 g), with variable shapes from conical to elongated (average FSIEII = 1.84), which mainly converge on the negative panel delimited by the first and second components (Figure 6). The third component allowed to discriminate accessions according to their colors. Hence, yellow fruits with AHue and ALV values close to 100 and 70, respectively, were only detected for C. chinense and C. baccatum. The other species displayed red peppers with variable amounts of red color and lightness. Remarkably, the majority of accessions within the C. pubescens group exhibited the most reddish fruits with the lowest lightness (Figure 6).
Groups of promising accessions possesing interesting traits were identified with an hierarchical cluster analysis by making use of those parameters highly contributing to the PCA variance. Clustering pattern was independent of Capsicum spp. or geographical origin and was primarily based upon fruit morphology, color and biochemical content (Figure 7). Two main groups (I and II) were differentiated. The first one (I) comprised accessions with bigger and heavier fruits, possessing high moisture content but low values for the biochemical parameters. The fruit shapes varied from nearly round (pear-like shape) to triangular or slightly elongated. The sub-group I.1 consisted of C. pubescens fruits with manzano type morphology showing an intense red color with low lightness, and C. chinense fruits of the habanero type (Figure 7). The sub-group I.2 harbour those fruits with the highest sizes and weights, primarily represented by C. annuum accessions. The subdivision of group II allowed to establish two subgroups, which could be additionally partitioned. Hence, the subgroup II.1.1 clustered C. annuum and C. baccatum wild accessions with the smallest fruits and lowest amounts of water. Interestingly, that group possessed high values for the majority of biochemical components (Figure 7). Subgroup II.1.2, which included primarily C. chinense and C. frutescens, also grouped small fruits, but with higher moisture content than II.1.1 and more roundish shapes. Yellow peppers of small size and roundish shapes conformed a subgroup within II.1.2. The remaining yellow fruits, mostly C. chinense and C. baccatum, were included into the sub-group II.2.1, since they showed more elongated morphologies. Sub-group II.2.2 comprised medium-size peppers of elongated forms. It is remarkable within that group the presence of various accessions with the highest fruit shape indices, which measured the ratio of height to width. Pungent C. annuum were included into the group II.1.1, whereas the majority of non-pungent accessions were circumscribed to groups I.2 and II.2.2 (Figure 7).
The principal coordinate analysis (PCoA) carried out for each species confirmed the HCA results. Clustering of accessions mainly responded to size and shape, although some correspondences to the original geographical region could be detected (Figure S2). Hence, C. annuum peppers from Ecuador and Bolivia had bigger sizes and more elongated shapes than those from Mexico, which conversely displayed higher pungencies. C. baccatum fruits from Bolivia were smaller and more roundish than those from Ecuador. Interestingly, all C. pubescens accessions were located on the first quadrant, exhibiting fruits of medium-big sizes and roundish shapes.

4. Discussion

Chemical parameters such as SSC, pH or fruit acidity have been investigated to a lesser extent than other nutritional attributes including capsaicinoids, ascorbic acid or phenolic composition and the majority of works have been performed on the most commercially relevant C. annuum. However, these components are critical to the overall sensory quality of fruits and can be employed to determine the end-use and the preferable market for pepper commercialization [25,40]. It is well-known that SSC and acidity increase as ripening progresses, with fully ripe fruits possessing the highest levels of both components [40,41]. Significant differences among accessions and species were observed for those parameters in our Capsicum spp. panel, composed of mature fruits. The SSC and titratable acidity ranges observed for C. annuum were in agreement with those in previous reports [27,42,43]. Outstanding levels were detected in the accession CAP524, which appeared catalogued as C. frutescens by the germplasm bank. However, our earlier genetic analyses classified it as C. annuum [29,30], and recent work based on genotyping-by-sequencing suggested that this accession might be considered as a wild C. annuum (unpublished data). C. chinense showed the poorest performance for those parameters, while C. frutescens displayed similar levels to C. annuum, although such results must be considered with caution due to the low number of C. frutescens accessions assessed in this work. The highest SSC/acidity ratio, which is informative of how tasty the fruits are, was found in the C. annuum group, while C. pubescens appeared as the least tasty. The SSC mainly comprise sugars derived from degradation of polysaccharides during ripening, while acidity comes from an increasing concentration of organic acids, particularly ascorbic acid [40,44,45]. This acid, together with phenolic compounds are principal contributors to the renowned antioxidant capacity of cultivated peppers [9]. In our report, polyphenols and TEAC showed a large variation across the different accessions, but not among the five Capsicum species, which exhibited a shorter range, expect for C. frutescens. This suggested that domesticated species contained accessions with low, medium and high levels of phenols and antioxidant capacity. Such a conclusion concurred with previous works that frequently quantified the highest and lowest values for those components in accessions from the same species [46,47,48]. On the other hand, we determined those parameters on fresh fruits, which hindered the comparison to similar works, mostly performed on dry bases [44,46,49]. Indeed, Materska [50] reported that lower levels of phenolic acid derivatives are generally found in lyophilised than in fresh fruits. Nevertheless, our results are in agreement with those previously observed on ripe fresh fruits; although values recorded in our work were slightly lower [51,52,53]. The accessions of this study with the highest amounts of both parameters (ex. CAP524, CGN20808, PI585257 or PI595907) could be a promising source of antioxidants for the human diet.
A high and significant correlation between polyphenols and TEAC was obtained (R2 = 0.92, p < 0.001), in agreement with previous reports [42,46,53]. Similarly, the available literature demonstrated that antioxidant activity and polyphenols tend to concurrently increase with the amounts of ascorbic acid and capsaicin [23,45,48,54]. The latter compounds were not assessed in the present work, however the correlation found between TEAC, phenols and titratable acidity was positive, high and significant (R2 = 0.63 and 0.72, respectively at p < 0.001), suggesting that those accessions with great levels of antioxidant capacity and polyphenols might also contain high amounts of ascorbic acid. In the same way, C. annuum peppers catalogued as pungent in our study, which theoretically hold a notable quantity of capsaicinoids, showed significantly higher values for TEAC, phenols and titratable acidity. Therefore, both TEAC and phenols, which are easy to measure parameters, may act as predictors of other attributes, such as capsaicinoids and ascorbic acid, whose quantification is technically more complex.
Water is the most plentiful component in peppers. Its abundance largely influences the post-harvest quality of fruits as well as their market destination [55]. Moisture content varied among the accessions in our panel although the majority oscillated within a range of 70–90%, showing a coefficient of variation of 9.2%. Similar values were reported by other authors [56,57,58]. Minor variations were recorded among Capsicum groups, suggesting that peppers with low and high moisture levels exist across the different species. Frequently, pepper fruits must be dehydrated for commercial purposes, in particular when they are destined to powder production, in which case low levels of water are more appreciated [27]. Percentages of water content lower than 65%, i.e., levels of dry matter over 35%, were recorded for six accessions in our collection (two C. baccatum and four C. annuum). These accessions (e.g., CGN17042, CGN20808 or CGN21460) might be explored for the manufacturing of dry peppers. Interestingly, moisture content displayed strong negative correlations (p < 0.001) with chemical components, indicating that peppers with higher percentages of water will exhibit lower levels of SSC, acidity, polyphenols and antioxidant capacity. Likewise, an increment in water content would involve a reduction in the potential pungency of fruits. Comparable results were described by Tripodi et al. [27], who reported that pepper types employed for the production of powder contained higher values of dry matter and correlative amounts of capsaicin and other phytochemical traits. Small contrasts were observed among peppers from different geographical regions for their content on chemical compounds. This contradicted the conception that environmental growing conditions largely determine the biochemical composition of fruits [59,60]. However, such absence of differences might be not only be attributed to the comparable climate characteristics of the three Andean countries, but also to the similarities in consumers’ preferences, which may have conditioned the farmers’ selections and a common evolution in the fruit compositional profiles.
Fifty-five plant and fruit characters where investigated by checking traditional descriptors and highly precise data acquired with a digital tool. Conventional traits describing plant height and fruit weight mostly explained the variability among pepper accessions, as previously reported [61,62]. Regarding the TA analysis, the attributes encompassed into the categories of basic measurements, shape index, and latitudinal section, exhibited the highest contribution to the variance. Similar findings were reported not only for pepper [19,20,23] but also for tomato [16,24] and eggplant [17]. The ANOVA revealed that C. annuum statistically diverged from the other species, exhibiting the maximum values for traits related to size, while C. frutescens showed the lowest measurements for identical parameters. This is in line with results described by Tripodi and Greco [20] on a large collection of cultivated and wild peppers. Colour parameters also explained a high percentage of the variance within our Capsicum collection. This is something expected since the domestication process and subsequent selection procedures led to the generation of cultivated peppers with various ranges of colors, suitable for different uses [63]. The yellow, orange, and red colors in peppers originate from carotenoid pigments produced in the fruit during ripening, which apart from their nutritional value, can act as antioxidants preventing different human diseases [64,65]. Besides, colour will greatly determine whether peppers might be employed as colorants for food or in the cosmetics industry. The majority of Andean peppers possessed a hue angle in the orange range, although some displayed AHue values lower than 50° or higher than 80° which are indicative of intense red or pure yellow, respectively. The former group mainly comprised mainly accessions from C. pubescens, while the latter consisted of C. chinense peppers. Furthermore, a strong positive interaction was observed between AHue and Luminosity, C. chinense fruits being the lightest and C. pubescens the darkest. Similar results were reported by Meckelmann et al. [46] in a collection of native chilli peppers from Peru. Although no previous studies have undoubtedly correlated the pepper fruit color to the carotenoid accumulation, it appears that capsanthin, capsorubin and cryptocapsin could be responsible for red colors while the yellow and orange colors come from β-carotene, zeaxanthin and violaxanthin [66,67]. Therefore, those Andean pepper fruits with larger redness or yellowness components might possess notable amounts of carotenoids and they could be potentially novel sources of beneficial compounds for human health, or they might exploited for the production of natural additives. C. annuum accessions were assorted into putative pungent and non-pungent types based on the marker MAP1 [31]. Differences between sweet and hot peppers were found for all biochemical parameters (except pH), for PHE, FWE, FPL, and 31 TA descriptors. In general, pungent types had taller plants, smaller fruits in both length and width, slightly lower shape indexes, less fresh and dry weight but much higher values of SSC, titratable acidity, polyphenols and antioxidant activity. This is in line with data from Hill et al. [68], who identified 17 conserved regions across the genome of non-pungent peppers, and reported that many of those regions overlapped with QTLs for fruit size and shape, but also with genes involved in the capsaicin biosynthesis.
To date, few works have attempted to evidence the impact of fruit morphology on biochemical compositions [25,26,27]. However, to investigate the correlation between these two groups of attributes emerge as a mandatory task, since this linkage might affect the speed at which favorable attributes could be introgressed into commercial pepper cultivars. All chemical traits, except moisture content and dry weight, negatively correlated to FWE, TPA, TPT and parameters belonging to basic measurements. Conversely, positive correlations were observed to fruit indexes (FSIEI, FSIEII, CFSI), circular, DFB and color parameters except for AaV. It seems that pepper fruits with higher weights, larger sizes, thicker pericarps, low height/width ratios, and greater proportions of red color contained lower amounts of SSC, titratable acidity, polyphenols and antioxidant activity, but higher levels of water and dry matter. Negative correlations between morpho-agronomic and phytochemical traits were also evidenced by other authors [23,25,27,40]. Fruit weight was strongly and significantly (p < 0.05) correlated in a positive manner to TA parameters related to size (P, A, WMH, MW, HMW, MH, and CH), but also to PIA, TPA, and TPT. Negative correlations to fruit shape indexes, internal eccentricity and the majority of color attributes were also detected, although those were not significant. Such results suggested that heavier fruits displayed bigger sizes, thicker pericarps, and more elongated shapes, i.e., lower ratios for height to width. Vilarinho et al. [25] also reported high positive correlations between fruit weight, fruit width and pericarp thickness, supporting the data by Barchi et al. [69], who found that QTLs controlling those traits were linked on chromosome P12.
The network of associations for TA parameters revealed positive and significantly strong correlations between size-based traits and those connected to the pericarp magnitude (TPT and TPA). In addition, HMW, MH and CH showed significant positive correlations to fruit shape indexes, internal eccentricity parameters, ellipsoid and circular. Significant associations among the categories of homogeneity, fruit shape index and internal eccentricity were observed as well. All these results bring to light that smaller fruits possessed shapes that are more roundish and, have thinner walls and lower dimensions of width and height axes. Similar results were obtained in our previous work [23], and also by Tripodi and Greco [20] in a large collection of worldwide peppers. As expected, significant interactions were observed among colors parameters estimated with both RGB and CIELab color spaces implemented into TA software [33]. Unlike previous works, we could not identify strong and significant correlations among ALV, AHue and ACh [25,40,43]. However, the correlation matrix revealed that an increase in lightness significantly occurred with an increase in Hue angle, suggesting that yellow-orange peppers possessed more brilliant colors.
The most effective attributes in discriminating among accessions were determine by principal component analysis. PC1 was positively and strongly correlated to parameters defining the fruit shape, homogeneity, and internal eccentricity while PC2 displayed high coefficients with traits associated to size, latitudinal sections, FWE, and FDW. PC3 was primarily explained by those parameters related to color. Biochemical traits negatively contributed to the first two axes although with very low and non-significant coefficients. These results supported previous works and confirmed that a reduced sub-set of descriptors within the sections of basic measurements, fruit shape and latitudinal section emerged as the most informative, in such a way that they might be used alone to discriminate between Capsicum types [19,20,21,22,23]. The representation of accessions on the bi-dimensional space proved that each species was diverse for fruit morphological and colorimetric traits, as accessions of a single group presented a certain degree of dispersion, plotting at different areas of the graph. This was particularly evident for C. annuum, which displayed variable accessions from those with large fruits, thick walls, moderate height/width ratios, and triangular shapes (Bell types), to those with tiny, thin-walled fruits, high height/width ratios, and extremely elongated shapes (Cayenne type). Wild C. annuum from Mexico (PM647, CGN21526 and CGN22783) exhibited the smallest fruits, conical shapes and grouped together with C. annuum accessions CAP524 and CGN20808, which might be considered closely-related to ancestors. The fruits from C. chinense exhibited more variation in size than in shape, with fruit weight ranging from 0.48 to 19.02 g, and shapes oscillating from rectangular forms (Habanero type) to more triangular profiles. C. baccatum fruits possessed variable sizes, the wild forms (var. baccatum) being the smallest ones. Shapes were mostly elongated, with high height/width ratios. On the contrary, accessions from C. pubescens had berries of medium size with roundish shapes. Fruits with yellow and orange colours were mostly observed in C. chinense and C. baccatum, which also displayed the highest lightness. The smallest values for luminosity were attributed to C. pubescens fruits. The fact that most of the variation in fruit features occurred in C. annuum was somehow expected. This Capsicum represents the most commonly cultivated and consumed around the world and the one that has suffered from more intensive human selection pressures, who greatly modified the morphometry of fruits out of the wild ancestor C. annnuum var. glabriusculum [70].
The clustering analysis reflected similar results to those obtained in PCA, with Andean peppers being arranged independently of their species and according to their fruit morphology, color, and biochemical composition. Hence, roundish peppers of moderate size, high amounts of water and intense red appeared into group I.1, and they predominantly corresponded to C. pubescens and C. chinense accessions. Cluster I.2 grouped the peppers with the bigger sizes, higher weights, large water contents, triangular shapes and light red colors, the majority (62.5%) of them being C. annuum. Yellow peppers were distributed across clusters II.1.2 and II.2.1 and they essentially consisted of C. chinense and C. baccatum accessions of small size and weight, whose shape varied from almost round to highly elongated. The cluster II.1.1 contained the wild ancestors of C. annuum and C. baccatum, which displayed the expected tiny and slight fruits, with very low moisture contents, conical shapes and red-orange colors. Interestingly, peppers within this group showed the highest values for biochemical parameters and they were pungent according to the marker associated to the Pun1 gene. C. annuum accessions with highly elongated shapes (Cayene type) were located in cluster II.2.2 and they displayed both pungent and non-pungent types. These results confirmed our previous work [23] performed on a smaller set of Ecuadorian native landraces and agreed with those by Nankar et al. [22], who observed that a collection of Capsicum accessions from the Balkan Peninsula clustered according fruit shape, size and color. However, they differed from those reported by Tripodi and Greco [20] who noticed that a dendrogram obtained with conventional and TA descriptors split accessions according to their Capsicum species. Such differences might be attributed to the size and origin of the evaluated panels. On the other hand, these authors found a clear distinction between wild and cultivated species, what is in agreement with our work, at which C. annuum wild and C. baccatum var. baccatum accessions are assorted in separated cluster.
PCoA analyses per species considering the different geographical regions revealed certain divergences in fruit morphology according to the country of origin. Thus, C. annuum from Mexico displayed the tiniest and most pungent fruits, while Ecuador and Peru possessed the biggest non-pungent berries. Comparable results were noticed for the C. baccatum group, where the smallest fruits corresponding to C. baccatum var. baccatum accessions, were found only in Bolivia. These observations are in line with our earlier work [29], which concluded that genetic differentiations among Capsicum spp. followed a geographic pattern. Hence, higher genetic similarities were detected between Ecuadorian and Peruvian C. annuum peppers, which might be a consequence of a more recent gene flow across those regions [71]. The influence of ecogeographic distributions on the patterns of genetic differentiation within wild and domesticated C. annuum from Mexico was also detected previously and the might be related to Mexican origin of this Capsicum species [72]. Likewise, correlations between genetic variability and geographic arrangements within C. baccatum accessions from South America were established by Albrecht et al. [73], who proposed that Ecuadorian and Peruvian accessions likely derived from Bolivian wild-types. All those genetic differences might explain the variable morphology observed among Andean countries, since at each region fruits would be selected by farmers to meet consumers’ preferences and agro-climatic requirements.

5. Conclusions

The present work brings to light that native varieties from the Andean region, one of the places postulated to be the origin of the genus Capsicum, possess a substantial diversity of pepper fruits with different sizes, shapes, and colors carrying variable amounts of beneficial metabolites. These varieties might be exploited as sources of desirable traits lacking in modern cultivars. Furthermore, negative correlations between fruit geometrics and phytochemical compounds were detected, reinforcing our previous work on a small set of Capsicum spp. from Ecuador [23]. Such a contribution is essential in breeding for fruit quality, as it would determine what characters might be independently selected and the speed at which valuable traits could be transferred into commercial peppers.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4395/11/2/262/s1, Table S1: List of accessions employed in this work. Na stands for not available, Table S2: List of digital traits measured with Tomato Analyzer software [12,13], Table S3: Analysis of variance for biochemical traits (above). Mean ± SD for each accession (below). Last row indicates the average value (Coefficient Variation) for each trait, Table S4: Analysis of variance for conventional descriptors (above). Mean ± SD for each accession (below). Last row indicates the average value (Coefficient Variation) for each descriptor. PHE = Plant height, PWI = Plant width, SLE = Stem length, SDI = Stem diameter, FEW = Fruit weight, FPL = Fruit pedicel length, Table S5: Analysis of variance for TA descriptors (above). Mean ± SD for each accession (below). Last row indicates the average value (Coefficient Variation) for each descriptor. P = Perimeter, A = Area, WMH = Width mid-height, MW = Maximum width, HMW = Height mid-width, MH = Maximum height, CH = Curved height, FSIEI = Fruit shape index external I, FSIEII = Fruit shape index external II, CFSI = Curved fruit shape index, PFB = Proximal fruit blockiness, DFB = Distal fruit blockiness, FST = Fruit shape triangle, E = Ellipsoid, C = Circular, R = Rectangular, SH = Shoulder height, PAMI = Proximal angle micro, PAMA = Proximal angle macro, PIA = Proximal indentation area, DAMI = Distal angle micro, DAMA = Distal angle macro, DIA = Distal indentation area, DEP = Distal end protrusion, Ob = Obovoid, Ov = Ovoid, VAs = V. Asymmetry, HAob = H. Asymmetry. Ob, HAov = H. Asymmetry. Ov, WWP = Width widest position, EC = Eccentricity, PEC = Proximal eccentricity, DEC = Distal eccentricity, FSII = Fruit shape index internal, ECAI = Eccentricity area index, LD = Lobedness degree, TPA = Tomato pericarp area, TPAR = Tomato pericarp area ratio, TPT = Tomato pericarp thickness, TPTR = Tomato pericarp thickness ratio, AR = Average Red, AG = Average Green, AB = Average Blue, AL = Average Luminosity, ALV = Average L Value, AaV = Average a Value, AbV = Average b Value, AHue = Average Hue, Ach = Average Chroma, Table S6: Mean—Coefficient of Variation (%) for conventional and digital descriptors in C. annuum, C. baccatum, C. chinense, and C. pubescens within each geographical region. Different letters within the same column indicate significant differences at p < 0.05. Only the parameters showing significant differences among regions are represented. See Table S2 for traits’ acronyms, Figure S1: Contribution coefficients of each trait to the three Principal Components (PCs), Figure S2: Principal Coordinate Analysis (PCoA) for each Capsicum species based on twenty-six biochemical, conventional and TA descriptors. Accessions are represented according to their country of origin.

Author Contributions

C.S.: conceptualization and experiment design, J.G.-L., S.R.-M.: experiments and data analysis; C.S.: writing—original draft preparation, C.S.: writing—review and editing, C.S.: supervision, project administration, and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Spanish Institute for Agricultural and Food Research and Technology (INIA), co–financed by the European Regional Development Fund (FEDER) (Project RTA2015-00042-C02-02), and by Xunta de Galicia (Project ED431C 2018/57).

Acknowledgments

We gratefully acknowledge provision of germplasm by USDA, ARS, Plant Genetic Resources Conservation Unit (USA), Center for Genetic Resources (The Netherlands), Institute for Conservation and Improvement of Valencian Agrodiversity (Spain) and Leibniz Institute of Plant Genetics and Crop Plant Research (Germany).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bosland, P.W.; Votava, E.J. Peppers: Vegetable and spice Capsicum, 2nd ed.; CABI Publishing: Wallingford, UK, 2012. [Google Scholar]
  2. Carrizo García, C.; Barfuss, M.H.J.; Sehr, E.M.; Barboza, G.E.; Samuel, R.; Moscone, E.A.; Ehrendorfer, F. Phylogenetic relationships, diversification and expansion of chili peppers (Capsicum, Solanaceae). Ann. Bot. 2016, 118, 35–51. [Google Scholar] [CrossRef] [Green Version]
  3. Carrizo García, C.; Sterpetti, M.; Volpi, P.; Ummarino, M.; Saccardo, F. Wild Capsicums: Identification and in situ analysis of Brazilian species. In Breakthroughs in the Genetics and Breeding of Capsicum and Eggplant; Lanteri, S., Rotino, G.L., Eds.; Eucarpia: Torino, Italy, 2013; pp. 205–213. [Google Scholar]
  4. FAOSTAT. About Trade: Crops and Livestock Products. 2020. Available online: http://faostat.fao.org/site/535/default.aspx#ancor (accessed on 19 April 2020).
  5. Scossa, F.; Roda, F.; Tohge, T.; Georgiev, M.I.; Fernie, A.R. The Hot and the Colorful: Understanding the Metabolism, Genetics and Evolution of Consumer Preferred Metabolic Traits in Pepper and Related Species. Crit. Rev. Plant Sci. 2019, 38, 339–381. [Google Scholar] [CrossRef]
  6. Crosby, K.M. Pepper. In Handbook of Plant Breeding. Vegetables II; Prohens, J., Nuez, F., Eds.; Springer: New York, NY, USA, 2008; pp. 221–248. [Google Scholar]
  7. Buckenhüskes, H.J. Current requirements on paprika powder for food industry. In Capsicum: The Genus Capsicum; Krishna, A., Ed.; Taylor and Francis Ltd.: London, UK, 2003; pp. 223–230. [Google Scholar]
  8. Brugarolas, M.; Martinez-Carrasco, L.; Martinez-Poveda, A.; Ruiz, J.J. A competitive strategy for vegetable products: Traditional varieties of tomato in the local market. Span. J. Agric. Res. 2009, 7, 294–304. [Google Scholar] [CrossRef] [Green Version]
  9. Mendes, N.; de Andrade Gonçalves, É.C.B. The role of bioactive components found in peppers. Trends Food Sci. Technol. 2020, 99, 229–243. [Google Scholar] [CrossRef]
  10. Pizzino, G.; Irrera, N.; Cucinotta, M.; Pallio, G.; Mannino, F.; Arcoraci, V.; Squadrito, F.; Altavilla, D.; Bitto, A. Oxidative Stress: Harms and Benefits for Human Health. Oxid. Med. Cell. Longev. 2017, 12, 1–13. [Google Scholar] [CrossRef] [PubMed]
  11. Paran, I.; Fallik, E. Breeding for Fruit Quality in Pepper (Capsicum spp.). In Breeding for Fruit Quality; Jenks, M.A., Bebeli, B.J., Eds.; Wiley: Hoboken, NJ, USA, 2011; pp. 307–322. [Google Scholar]
  12. Brewer, M.T.; Lang, L.; Fujimura, K.; Dujmovic, N.; Gray, S.; van der Knaap, E. Development of a Controlled Vocabulary and Software Application to Analyze Fruit Shape Variation in Tomato and Other Plant Species. Plant. Physiol. 2006, 141, 15–25. [Google Scholar] [CrossRef] [Green Version]
  13. Rodríguez, G.R.; Moyseenko, J.B.; Robbins, M.D.; Morejón, N.H.; Francis, D.M.; van der Knaap, E. Tomato Analyzer: A Useful Software Application to Collect Accurate and Detailed Morphological and Colorimetric Data from Two-dimensional Objects. JoVE J. Vis. Exp. 2010, 37, e1856. [Google Scholar] [CrossRef] [Green Version]
  14. Gonzalo, M.J.; Brewer, M.T.; Anderson, C.; Sullivan, D.; Gray, S.; van der Knaap, E. Tomato Fruit Shape Analysis Using Morphometric and Morphology Attributes Implemented in Tomato Analyzer Software Program. J. Am. Soc. Hortic. Sci. 2009, 134, 77–87. [Google Scholar] [CrossRef] [Green Version]
  15. Rodríguez, G.R.; Muños, S.; Anderson, C.; Sim, S.-C.; Michel, A.; Causse, M.; Gardener, B.B.M.; Francis, D.; van der Knaap, E. Distribution of SUN, OVATE, LC, and FAS in the Tomato Germplasm and the Relationship to Fruit Shape Diversity. Plant. Physiol. 2011, 156, 275–285. [Google Scholar] [CrossRef] [Green Version]
  16. Figàs, M.R.; Prohens, J.; Raigón, M.D.; Fernández-de-Córdova, P.; Fita, A.; Soler, S. Characterization of a collection of local varieties of tomato (Solanum lycopersicum L.) using conventional descriptors and the high-throughput phenomics tool Tomato Analyzer. Genet. Resour. Crop. Evol. 2015, 62, 189–204. [Google Scholar]
  17. Hurtado, M.; Vilanova, S.; Plazas, M.; Gramazio, P.; Herraiz, F.J.; Andújar, I.; Prohens, J. Phenomics of fruit shape in eggplant (Solanum melongena L.) using Tomato Analyzer software. Sci. Hortic. 2013, 164, 625–632. [Google Scholar] [CrossRef]
  18. Plazas, M.; Andújar, I.; Vilanova, S.; Gramazio, P.; Herraiz, F.J.; Prohens, J. Conventional and phenomics characterization provides insight into the diversity and relationships of hypervariable scarlet (Solanum aethiopicum L.) and gboma (S. macrocarpon L.) eggplant complexes. Front. Plant Sci. 2014, 5, 318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Naegele, R.P.; Mitchell, J.; Hausbeck, M.K. Genetic Diversity, Population Structure, and Heritability of Fruit Traits in Capsicum annuum. PLoS ONE 2016, 11, e0156969. [Google Scholar] [CrossRef] [PubMed]
  20. Tripodi, P.; Greco, B. Large Scale Phenotyping Provides Insight into the Diversity of Vegetative and Reproductive Organs in a Wide Collection of Wild and Domesticated Peppers (Capsicum spp.). Plants 2018, 7, 103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Pereira-Dias, L.; Fita, A.; Vilanova, S.; Sánchez-López, E.; Rodríguez-Burruezo, A. Phenomics of elite heirlooms of peppers (Capsicum annuum L.) from the Spanish centre of diversity: Conventional and high-throughput digital tools towards varietal typification. Sci. Hortic. 2020, 265, 109245. [Google Scholar] [CrossRef]
  22. Nankar, A.N.; Tringovska, I.; Grozeva, S.; Todorova, V.; Kostova, D. Application of high-throughput phenotyping tool Tomato Analyzer to characterize Balkan Capsicum fruit diversity. Sci. Hortic. 2020, 260, 1–12. [Google Scholar] [CrossRef]
  23. García-González, C.A.; Silvar, C. Phytochemical Assessment of Native Ecuadorian Peppers (Capsicum spp.) and Correlation Analysis to Fruit Phenomics. Plants 2020, 9, 986–1011. [Google Scholar]
  24. Nankar, A.N.; Tringovska, I.; Grozeva, S.; Ganeva, D.; Kostova, D. Tomato Phenotypic Diversity Determined by Combined Approaches of Conventional and High-Throughput Tomato Analyzer Phenotyping. Plants 2020, 9, 197. [Google Scholar] [CrossRef] [Green Version]
  25. Vilarinho, L.B.O.; da Silva, D.J.H.; Greene, A.; Salazar, K.D.; Alves, C.; Eveleth, M.; Nichols, B.; Tehseen, S.; Khoury, J.K., Jr.; Johnson, J.V.; et al. Inheritance of fruit traits in Capsicum annuum: Heirloom cultivars as sources of quality parameters relating to pericarp shape, color, thickness, and total soluble solids. J. Am. Soc. Hortic. Sci. 2015, 140, 597–604. [Google Scholar] [CrossRef] [Green Version]
  26. Mennella, G.; D’Alessandro, A.; Francese, G.; Fontanella, D.; Parisi, M.; Tripodi, P. Occurrence of variable levels of health-promoting fruit compounds in horn-shaped Italian sweet pepper varieties assessed by a comprehensive approach. J. Sci. Food Agric. 2018, 98, 3280–3289. [Google Scholar] [CrossRef]
  27. 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]
  28. González-Pérez, S.; Garcés-Claver, A.; Mallor, C.; Sáenz de Miera, L.E.; Fayos, O.; Pomar, F.; Merino, F.; Silvar, C. New Insights into Capsicum spp. Relatedness and the Diversification Process of Capsicum annuum in Spain. PLoS ONE 2014, 9, e116276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Silvar, C.; García-González, C.A. Deciphering Genetic Diversity in the Origins of Pepper (Capsicum spp.) and Comparison with Worldwide Variability. Crop. Sci. 2016, 56, 3100–3111. [Google Scholar] [CrossRef]
  30. Silvar, C.; García-González, C.A. Screening old peppers (Capsicum spp.) for disease resistance and pungency-related traits. Sci. Hortic. 2017, 218, 249–257. [Google Scholar] [CrossRef] [Green Version]
  31. Rodríguez-Maza, M.; Garcés-Claver, A.; Park, S.W.B.; Kang, C.; Arnedo-Andrés, M.S. A versatile PCR marker for pungency in Capsicum spp. Mol. Breed. 2012, 30, 889–898. [Google Scholar] [CrossRef]
  32. IPGRI; AVRDC; CATIE. Descriptors for Capsicum (Capsicum spp.); International Plant Genetic Resources Institute: Rome, Italy, 1995. [Google Scholar]
  33. Darrigues, A.; Hall, J.; van der Knaap, E.; Francis, D.M.; Dujmovic, N.; Gray, S. Tomato analyzer-color test: A new tool for efficient digital phenotyping. J. Am. Soc. Hort. Sci. 2008, 133, 579–586. [Google Scholar] [CrossRef] [Green Version]
  34. Singleton, V.L.; Rossi, J.A. Colorimetry of Total Phenolics with Phosphomolybdic-Phosphotungstic Acid Reagents. Am. J. Enol. Vitic. 1965, 16, 144–158. [Google Scholar]
  35. Miller, N.J.; Rice-Evans, C.A. Factors influencing the antioxidant activity determined by the ABTS + radical cation assay. Free Radic. Res. 1997, 3, 195–199. [Google Scholar] [CrossRef]
  36. SPSS. SPSS Statistics for Windows; Version 17.0; SPSS Inc.: Chicago, IL, USA, 2008. [Google Scholar]
  37. Abdi, H. The Bonferonni and Šidák Corrections for Multiple Comparisons. In Encyclopedia of Measurement and Statistics; Salkind, N.J., Ed.; Sage: Thousand Oaks, CA, USA, 2007. [Google Scholar]
  38. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  39. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  40. Tadesse, T.; Hewett, E.W.; Nichols, M.A.; Fisher, K.J. Changes in physicochemical attributes of sweet pepper cv. Domino during fruit growth and development. Sci. Hortic. 2002, 93, 91–103. [Google Scholar] [CrossRef]
  41. Ghasemnezhad, M.; Sherafati, M.; Payvast, G.A. Variation in phenolic compounds, ascorbic acid and antioxidant activity of five coloured bell pepper (Capsicum annum) fruits at two different harvest times. J. Funct. Foods 2011, 3, 44–49. [Google Scholar] [CrossRef]
  42. Fox, A.J.; Pozo-Insfran, D.D.; Lee, J.; Sargent, S.; Talcott, S. Ripening-induced Chemical and Antioxidant Changes in Bell Peppers as Affected by Harvest Maturity and Postharvest Ethylene Exposure. HortScience 2005, 40, 732–736. [Google Scholar] [CrossRef] [Green Version]
  43. Corrêa, C.; Mendonça, V.Z.; Gouveia, A.M.; Carpanetti, M.G.; Tavares, A.E.; Lanna, N.D.; Evangelista, R.M.; Cardoso, A.I. Physicochemical and biochemical traits of sweet pepper hybrids as a function of harvest times. Food Chem. 2018, 257, 265–270. [Google Scholar] [CrossRef] [PubMed]
  44. Wahyuni, Y.; Ballester, A.R.; Sudarmonowati, E.; Bino, R.J.; Bovy, A.G. Metabolite biodiversity in pepper (Capsicum) fruits of thirty-two diverse accessions: Variation in health-related compounds and implications for breeding. Phytochemistry 2011, 72, 1358–1370. [Google Scholar] [CrossRef]
  45. Kantar, M.B.; Anderson, J.E.; Lucht, S.A.; Mercer, K.; Bernau, V.; Case, K.A.; Le, N.C.; Frederiksen, M.K.; DeKeyser, H.C.; Wong, Z.Z.; et al. Vitamin variation in Capsicum spp. provides opportunities to improve nutritional value of human diets. PLoS ONE 2016, 11, e0161464. [Google Scholar] [CrossRef]
  46. Meckelmann, S.W.; Riegel, D.W.; van Zonneveld, M.J.; Ríos, L.; Peña, K.; Ugas, R.; Quinonez, L.; Mueller-Seitz, E.; Petz, M. Compositional Characterization of Native Peruvian Chili Peppers (Capsicum spp.). J. Agric. Food Chem. 2013, 61, 2530–2537. [Google Scholar] [CrossRef]
  47. Loizzo, M.R.; Pugliese, A.; Bonesi, M.; Menichini, F.; Tundis, R. Evaluation of chemical profile and antioxidant activity of twenty cultivars from Capsicum annuum, Capsicum baccatum, Capsicum chacoense and Capsicum chinense: A comparison between fresh and processed peppers. LWT Food Sci. Technol. 2015, 64, 623–631. [Google Scholar] [CrossRef]
  48. Souza-Sora, G.T.; Haminiuk, C.W.I.; Vieira da Silva, M.; Ferreira Zielinski, A.A.; Almeida Gonçalves, G.; Bracht, A.; Peralta, R.M. A comparative study of the capsaicinoid and phenolic contents and in vitro antioxidant activities of the peppers of the genus Capsicum: An application of chemometrics. J. Food Sci. Technol. 2015, 52, 8086–8094. [Google Scholar] [CrossRef] [Green Version]
  49. Bogusz, S.; Libardi, S.H.; Dias, F.F.; Coutinho, J.P.; Bochi, V.C.; Rodrigues, D.; Melo, A.M.; Godoy, H.T. Brazilian Capsicum peppers: Capsaicinoid content and antioxidant activity. J. Sci. Food Agric. 2018, 98, 217–224. [Google Scholar] [CrossRef]
  50. Materska, M. Bioactive phenolics of fresh and freeze-dried sweet and semi-spicy pepper fruits (Capsicum annuum L.). J. Funct. Foods 2014, 7, 269–277. [Google Scholar] [CrossRef]
  51. Menichini, F.; Tundis, R.; Bonesi, M.; Loizzo, M.; Conforti, F.; Statti, G.; Cindio, B.; Houghton, P. The influence of fruit ripening on the phytochemical content and biological activity of Capsicum chinense Jacq. cv Habanero. Food Chem. 2009, 114, 553–560. [Google Scholar] [CrossRef]
  52. Zhuang, Y.; Chen, L.; Sun, L.; Cao, J. Bioactive characteristics and antioxidant activities of nine peppers. J. Funct. Foods 2012, 4, 331–338. [Google Scholar] [CrossRef]
  53. Carvalho, A.; Mattietto, R.A.; Rios, A.O.; Maciel, R.A.; Moresco, K.; Oliveira, T.C. Bioactive compounds and antioxidant activity of pepper (Capsicum sp.) genotypes. J. Food Sci. Technol. 2015, 52, 7457–7464. [Google Scholar] [CrossRef]
  54. Álvarez-Parrilla, E.; Rosa, L.A.; Amarowicz, R.; Shahidi, F. Antioxidant activity of fresh and processed Jalapeño and Serrano peppers. J. Agric. Food Chem. 2011, 59, 163–173. [Google Scholar] [CrossRef] [PubMed]
  55. Lownds, N.; Banaras, M.; Bosland, P. Relationships between postharvest water loss and physical properties of pepper fruit (Capsicum annuum L.). HortScience 1993, 28, 1182–1184. [Google Scholar] [CrossRef] [Green Version]
  56. Smith, D.L.; Stommel, J.; Fung, R.W.; Wang, C.Y.; Whitaker, B. Influence of cultivar and harvest method on postharvest storage quality of pepper (Capsicum annuum L.) fruit. Postharvest Biol. Technol. 2006, 42, 243–247. [Google Scholar] [CrossRef]
  57. Ziino, M.; Condurso, C.; Romeo, V.; Tripodi, G.; Verzera, A. Volatile compounds and capsaicinoid content of fresh hot peppers (Capsicum annuum L.) of different Calabrian varieties. J. Sci. Food Agric. 2009, 89, 774–780. [Google Scholar] [CrossRef]
  58. Giuffrida, D.; Dugo, P.; Torre, G.; Bignardi, C.; Cavazza, A.; Corradini, C.; Dugo, G. Characterization of 12 Capsicum varieties by evaluation of their carotenoid profile and pungency determination. Food Chem. 2013, 140, 794–802. [Google Scholar] [CrossRef]
  59. Gurung, T.; Techawongstien, S.; Suriharn, B.; Techawongstien, S. Impact of Environments on the Accumulation of Capsaicinoids in Capsicum spp. HortScience 2011, 46, 1576–1581. [Google Scholar] [CrossRef]
  60. Tripodi, P.; Cardi, T.; Bianchi, G.; Migliori, C.A.; Schiavi, M.; Rotino, G.L.; Lo Scalzo, R. Genetic and environmental factors underlying variation in yield performance and bioactive compound content of hot pepper varieties (Capsicum annuum) cultivated in two contrasting Italian locations. Eur. Food Res. Technol. 2018, 244, 1555–1567. [Google Scholar] [CrossRef]
  61. Bozokalfa, M.K.; Esiyok, D.; Turhan, K. Patterns of phenotypic variation in a germplasm collection of pepper (Capsicum annuum L.) from Turkey. Span. J. Agric. Res. 2009, 7, 83–95. [Google Scholar] [CrossRef]
  62. Rivera, A.; Monteagudo, A.B.; Igartua, E.; Taboada, A.; García-Ulloa, A.; Pomar, F.; Riveiro-Leira, M.; Silvar, C. Assessing genetic and phenotypic diversity in pepper (Capsicum annuum L.) landraces from North-West Spain. Sci. Hortic. 2016, 203, 1–11. [Google Scholar] [CrossRef] [Green Version]
  63. Paran, I.; Knaap, E.V. Genetic and molecular regulation of fruit and plant domestication traits in tomato and pepper. J. Exp. Bot. 2007, 58, 3841–3852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Hassan, N.M.; Yusof, N.A.; Yahaya, A.F.; Rozali, N.N.; Othman, R. Carotenoids of Capsicum Fruits: Pigment Profile and Health-Promoting Functional Attributes. Antioxidants 2019, 8, 469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Rodríguez-Concepcíon, M.; Avalos, J.; Bonet, M.L.; Boronat, A.; Gómez-Gómez, L.; Hornero-Méndez, D.; Limón, M.C.; Meléndez-Martínez, A.J.; Olmedilla-Alonso, B.; Palou, A.; et al. A global perspective on carotenoids: Metabolism, biotechnology, and benefits for nutrition and health. Prog. Lipid Res. 2018, 70, 62–93. [Google Scholar]
  66. Arimboor, R.; Natarajan, R.B.; Menon, K.; Chandrasekhar, L.P.; Moorkoth, V. Red pepper (Capsicum annuum) carotenoids as a source of natural food colors: Analysis and stability—A review. J. Food Sci. Technol. 2014, 52, 1258–1271. [Google Scholar] [CrossRef] [Green Version]
  67. Berry, H.M.; Rickett, D.; Baxter, C.J.; Enfissi, E.M.; Fraser, P. Carotenoid biosynthesis and sequestration in red chilli pepper fruit and its impact on colour intensity traits. J. Exp. Bot. 2019, 70, 2637–2650. [Google Scholar] [CrossRef] [Green Version]
  68. Hill, T.; Chunthawodtiporn, J.; Ashrafi, H.; Stoffel, K.; Weir, A.; Deynze, A.V. Regions Underlying Population Structure and the Genomics of Organ Size Determination in Capsicum annuum. Plant. Genome 2017, 10, 1–14. [Google Scholar] [CrossRef] [Green Version]
  69. Barchi, L.; Lefèbvre, V.; Sage-Palloix, A.; Lanteri, S.; Palloix, A. QTL analysis of plant development and fruit traits in pepper and performance of selective phenotyping. Theor. Appl. Genet. 2009, 118, 1157–1171. [Google Scholar] [CrossRef]
  70. Qin, C.; Changshui, Y.; Shen, Y.; Fang, X.; Chen, L.; Min, J.; Cheng, J.; Zhao, S.; Xu, M.; Luo, Y.; et al. Whole genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization. Proc. Natl. Acad. Sci. USA 2014, 111, 5135–5140. [Google Scholar] [CrossRef] [Green Version]
  71. Ibiza, V.P.; Blanca, J.; Cañizares, J.; Nuez, F. Taxonomy and genetic diversity of domesticated Capsicum species in the Andean region. Genet. Resour. Crop. Evol. 2012, 59, 1077–1088. [Google Scholar] [CrossRef]
  72. Kraft, K.H.; Brown, C.H.; Nabhan, G.P.; Luedeling, E.; Luna Ruiz, J.J.; d’Eeckenbrugge, G.C.; Hijmans, R.J.; Gepts, P. Multiple lines of evidence for the origin of domesticated chili pepper, Capsicum annuum, in Mexico. Proc. Natl. Acad. Sci. USA 2014, 111, 6165–6170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Albrecht, E.; Zhang, D.; Deslattes Mays, A.; Staftner, R.A.; Stommel, J.R. Genetic diversity in Capsicum baccatum is significantly influenced by its ecogeographical distribution. BMC Genet. 2012, 13, 68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Box plot of biochemical parameters. Twenty-fifth percentile, median, 75th percentile, and range minimum-maximum. Outliers (white circle) were identified as 1.5 times the interquartile range. Red dashed lines represent the mean. Different letters indicate significant differences at p < 0.05. A = C. annuum, B = C. baccatum, C = C. chinense, F = C. frutescens, P = C. pubescens. Note that C. frutescens is not represented for moisture content and fruit dry weigh (FDW) because sufficient fruits were not available for the calculation of those parameters.
Figure 1. Box plot of biochemical parameters. Twenty-fifth percentile, median, 75th percentile, and range minimum-maximum. Outliers (white circle) were identified as 1.5 times the interquartile range. Red dashed lines represent the mean. Different letters indicate significant differences at p < 0.05. A = C. annuum, B = C. baccatum, C = C. chinense, F = C. frutescens, P = C. pubescens. Note that C. frutescens is not represented for moisture content and fruit dry weigh (FDW) because sufficient fruits were not available for the calculation of those parameters.
Agronomy 11 00262 g001
Figure 2. Box plot of conventional descriptors. Twenty-fifth percentile, median, 75th percentile, and range minimum-maximum. Outliers (white circle) were identified as 1.5 times the interquartile range. Red dashed lines represent the mean. Different letters indicate significant differences at p < 0.05. A = C. annuum, B = C. baccatum, C = C. chinense, F = C. frutescens, P = C. pubescens.
Figure 2. Box plot of conventional descriptors. Twenty-fifth percentile, median, 75th percentile, and range minimum-maximum. Outliers (white circle) were identified as 1.5 times the interquartile range. Red dashed lines represent the mean. Different letters indicate significant differences at p < 0.05. A = C. annuum, B = C. baccatum, C = C. chinense, F = C. frutescens, P = C. pubescens.
Agronomy 11 00262 g002
Figure 3. Bar plots of biochemical parameters expressed per country for C. annuum, C. baccatum, and C. chinense. Different letters indicate significant differences at p < 0.05. The absence of letters indicates non-significant differences. Only those parameters showing significant differences among countries are represented. Note that C. pubescens was not considered in the figure because significant differences were not observed among countries. A = C. annuum, B = C. baccatum, C = C. chinense.
Figure 3. Bar plots of biochemical parameters expressed per country for C. annuum, C. baccatum, and C. chinense. Different letters indicate significant differences at p < 0.05. The absence of letters indicates non-significant differences. Only those parameters showing significant differences among countries are represented. Note that C. pubescens was not considered in the figure because significant differences were not observed among countries. A = C. annuum, B = C. baccatum, C = C. chinense.
Agronomy 11 00262 g003
Figure 4. Bar plots of biochemical parameters expressed according to the potential pungency of C. annuum accessions. Different letters indicate significant differences at p < 0.05. Only those parameters showing significant differences among countries are represented.
Figure 4. Bar plots of biochemical parameters expressed according to the potential pungency of C. annuum accessions. Different letters indicate significant differences at p < 0.05. Only those parameters showing significant differences among countries are represented.
Agronomy 11 00262 g004
Figure 5. Pearson′s rank correlation coefficients between pairs of fruit traits. (A) Correlation coefficients for all pairs of traits. Positive and negative correlations are indicated according to the colored scale. The size of the squares is proportional to the correlation coefficients. (B) Correlation network for most significant traits (p < 0.05). Positive correlations are displayed in grey and negative correlations in black. The size of the circles is proportional to the correlation coefficients. SSC = soluble solid content, T. acidity = titratable acidity, TEAC = trolox equivalent antioxidant capacity, Moisture C = moisture content, FWE =fruit weight. Acronyms for TA parameters can be checked at Table S2.
Figure 5. Pearson′s rank correlation coefficients between pairs of fruit traits. (A) Correlation coefficients for all pairs of traits. Positive and negative correlations are indicated according to the colored scale. The size of the squares is proportional to the correlation coefficients. (B) Correlation network for most significant traits (p < 0.05). Positive correlations are displayed in grey and negative correlations in black. The size of the circles is proportional to the correlation coefficients. SSC = soluble solid content, T. acidity = titratable acidity, TEAC = trolox equivalent antioxidant capacity, Moisture C = moisture content, FWE =fruit weight. Acronyms for TA parameters can be checked at Table S2.
Agronomy 11 00262 g005
Figure 6. Scatter plot of the PCA analysis. (A) PC1 vs. PC2; (B) PC3 vs. PC1.
Figure 6. Scatter plot of the PCA analysis. (A) PC1 vs. PC2; (B) PC3 vs. PC1.
Agronomy 11 00262 g006
Figure 7. Cluster analysis based on seven biochemical traits, one conventional fruit descriptor and eighteen TA attributes. Red branches on the dendrogram are C. annuum; grey branches, C. chinense; yellow, C. frutescens; blue, C. pubescens; green, C. baccatum. The color code matrix represents the variables standardized to z-scores. Pictures are representative accessions for each group. White lines inside de photos indicate 1 cm scale. SSC = soluble solid content, T. acidity = titratable acidity, TEAC = Trolox Equivalent Antioxidant Capacity, Moisture C = moisture content, FEW = fruit weight. Acronyms for TA parameters can be checked at Table S2.
Figure 7. Cluster analysis based on seven biochemical traits, one conventional fruit descriptor and eighteen TA attributes. Red branches on the dendrogram are C. annuum; grey branches, C. chinense; yellow, C. frutescens; blue, C. pubescens; green, C. baccatum. The color code matrix represents the variables standardized to z-scores. Pictures are representative accessions for each group. White lines inside de photos indicate 1 cm scale. SSC = soluble solid content, T. acidity = titratable acidity, TEAC = Trolox Equivalent Antioxidant Capacity, Moisture C = moisture content, FEW = fruit weight. Acronyms for TA parameters can be checked at Table S2.
Agronomy 11 00262 g007
Table 1. Analysis of variance for Tomato Analyzer (TA) descriptors. Mean and coefficient variation (in italics) for each Capsicum species. Different letters within the same row indicate significant differences at p < 0.05. See Table S2 for trait acronyms.
Table 1. Analysis of variance for Tomato Analyzer (TA) descriptors. Mean and coefficient variation (in italics) for each Capsicum species. Different letters within the same row indicate significant differences at p < 0.05. See Table S2 for trait acronyms.
TraitSum of SquaresF value C. annuumC. baccatumC. chinenseC. frutescensC. pubescens
P429,598.732.0 ***159.56a138.72ab121.49b74.05c125.77b
50.21 39.36 36.26 25.55 19.51
A66,629,009.134.8 ***1242.11a840.47b760.86b300.00c946.34ab
84.79 59.83 62.87 42.02 38.09
WMH14,766.450.6 ***19.05b16.51bc19.37b12.74c26.59ab
60.23 37.33 41.04 14.90 24.65
MW13,637.939.1 ***24.44ab19.54c22.27bc14.39d28.50ab
53.63 35.34 35.77 20.07 22.28
HMW70,333.537.7 ***56.29a52.91a43.00b25.84c40.43b
49.92 41.75 39.75 26.67 25.50
MH100,288.644.7 ***62.33a56.97a45.83b26.55c43.10b
50.32 42.25 39.89 26.29 23.96
CH96,081.542.7 ***64.07a58.18ab47.71bc28.72d45.15c
49.31 41.10 38.31 24.87 22.80
FSIEI330.4108.6 ***2.69a2.96a2.18b1.84bc1.57c
34.78 31.89 40.60 15.44 29.06
FSIEII579.779.7 ***3.36a3.32a2.53b2.01bc1.61c
51.53 35.40 52.56 19.30 35.44
CFSI789.973.1 ***3.90a3.67a2.83b2.21bc1.78c
56.95 37.12 54.88 18.62 33.60
PFB11.474.7 ***1.03a0.84b0.84b0.83b0.82b
24.42 17.14 23.31 15.50 17.49
DFB0.77.01 ***0.589b0.544b0.550b0.668a0.565b
29.31 24.83 35.05 13.61 16.75
FST33.07.6 ***1.895a1.647ab1.807a1.249b1.492ab
47.00 31.02 89.90 17.21 26.18
E0.247.4 ***0.098a0.079b0.079b0.067b0.065b
39.65 32.28 41.15 26.55 25.70
C4.7135.6 ***0.313a0.319a0.245b0.219b0.157c
27.40 27.87 45.45 19.64 49.82
R0.516.9 ***0.418c0.423c0.428bc0.467ab0.476a
24.30 19.96 20.17 17.59 13.21
SH0.012.8 ***0.018ab0.011b0.013b0.013b0.027a
166.91 228.81 197.08 121.20 101.67
PAMI30,505.51.1 ns134.39a129.38a124.94a121.95a136.91a
67.99 52.05 65.88 76.47 68.43
PAMA312,234.027.2 ***99.45c84.53c107.75bc134.32a128.90ab
62.02 45.45 49.49 27.07 48.68
PIA0.334.15 ***0.030b0.011b0.012b0.011b0.054a
230.80 258.98 220.87 121.01 127.85
Ob0.17.4 ***0.004c0.015ab0.022a0.009bc0.007c
870.35 362.11 342.52 313.75 476.11
Ov5.771.9 ***0.381a0.240b0.254b0.212b0.261b
40.97 54.41 56.73 60.64 43.68
VAs2.126.5 ***0.223a0.163ab0.147bc0.086c0.120bc
89.66 71.02 79.25 71.18 56.92
HAob0.34.7 **0.007c0.026ab0.037a0.010bc0.006c
993.76 481.46 436.67 326.04 470.86
HAov47.361.2 ***0.757a0.526b0.382bc0.181c0.281c
79.88 82.77 83.32 71.74 59.60
WWP3.958.8 ***0.277b0.392a0.391a0.370a0.356a
47.07 33.39 35.10 30.95 27.42
DAMI249,026.513.82 ***87.48ab77.38b98.19ab108.29a117.08a
76.51 89.84 69.80 41.31 51.32
DAMA450,771.957.0 ***57.69bc52.16c76.91b119.21a100.75a
71.33 84.33 65.28 52.39 36.42
DIA0.00.5 ns0.004a0.004a0.004a0.005a0.003a
293.85 131.82 266.24 137.72 216.82
DEP5.942.6 ***0.101b0.221a0.104b0.059b0.041b
187.84 97.51 179.04 143.39 227.32
EC0.214.8 ***0.733b0.748ab0.755ab0.768a0.762a
9.19 6.38 6.35 2.90 4.36
PEC0.01.4 ns0.899a0.895a0.893a0.888a0.890a
9.39 4.02 4.17 1.25 0.55
DEC0.15.2 ***0.911a0.903a0.900a0.891a0.886a
11.01 5.68 5.22 2.52 1.33
FSII557.778.56 ***3.34a3.30a2.51b1.99bc1.63c
50.51 35.74 52.99 17.90 34.42
ECA0.851.6 ***0.477a0.441b0.424b0.437b0.422b
14.60 11.77 16.52 7.31 9.19
LD94,434.3109.1 ***31.88a34.82a22.61b16.39bc11.70c
50.12 44.20 67.60 25.84 67.23
TPA6,660,826.021.1 ***474.05a359.14ab326.27b140.59c422.15ab
88.91 60.27 58.90 42.69 36.99
TPAR0.114.2 ***0.417b0.424ab0.434ab0.439a0.439a
12.84 8.98 10.41 6.88 3.59
TPT173.537.6 ***2.72b2.61b2.76b1.98c3.64a
53.88 33.27 33.65 20.00 19.47
TPTR0.04.32 **0.212a0.206a0.206a0.204a0.203a
17.80 10.30 17.10 5.87 6.72
AR104,422.899.5 ***175.39b185.74a185.45a180.04ab161.58c
8.89 9.04 9.08 7.71 9.25
AG312,572.3101.8 ***87.62b100.76a109.18a87.04b64.68c
20.08 25.49 34.29 33.48 38.45
AB45,985.140.5 ***42.26ab40.05ab47.21a36.07bc29.27c
29.45 29.66 53.24 39.36 37.53
AL57,469.9104.3 ***102.40b106.23ab109.47a101.66b89.71c
10.66 8.19 13.90 11.43 11.49
ALV28,364.8116.9 ***47.80b52.17a54.21a48.52b40.99c
11.62 14.30 18.68 15.77 16.53
AaV23,304.147.6 ***28.78ab26.05bc21.74c30.40ab33.79a
22.09 35.27 72.80 32.07 30.71
AbV10,312.744.2 ***43.36b48.29a47.36a46.92a41.41b
11.02 18.63 19.53 9.90 13.67
AHue38,146.361.9 ***56.70bc61.11ab66.18a57.55b51.14c
10.65 17.83 27.11 16.25 22.24
ACh3096.515.9 ***52.30b55.95a54.57ab56.62a54.46ab
11.35 12.15 15.72 9.94 10.00
, **, *** Significant at p < 0.01, p < 0.001, ns not significant.
Table 2. Analysis of variance for conventional and TA descriptors between pungent and non-pungent C. annuum peppers. Only those parameters showing significant differences are represented. PHE = Plant height, FEW = Fruit weight, FPL = Fruit pedicel length, See Table S2 for acronyms of TA descriptors.
Table 2. Analysis of variance for conventional and TA descriptors between pungent and non-pungent C. annuum peppers. Only those parameters showing significant differences are represented. PHE = Plant height, FEW = Fruit weight, FPL = Fruit pedicel length, See Table S2 for acronyms of TA descriptors.
TraitPungent C. annuumNon-Pungent C. annuumSig.
MeanCV %MeanCV %
PHE111.6131.1590.3237.770.05
FWE7.18141.4517.7385.64<0.001
FPL3.1226.243.6923.80<0.001
P113.9262.90208.7726.86<0.001
A722.51111.931802.4255.41<0.001
WMH14.2770.4924.1944.08<0.001
MW18.3464.3031.0235.92<0.001
HMW39.8361.4774.0526.61<0.001
MH43.9562.3882.1526.61<0.001
CH45.6660.5683.9226.41<0.001
FSIEI2.5436.662.8532.11<0.001
FSIEII3.1855.953.5646.710.019
CFSI3.6862.114.1351.470.033
PFB0.9827.641.0720.24<0.001
DFB0.6123.230.5735.060.008
FST1.6941.072.1148.21<0.001
E0.0943.080.1035.500.001
C0.3031.830.3321.69<0.001
SH0.01175.790.02156.550.019
Ov0.3545.790.4134.99<0.001
VAs0.16101.960.2973.96<0.001
HAov0.55101.430.9958.38<0.001
WWP0.3042.740.2651.250.001
DAMI98.0166.1976.1288.54<0.001
DAMA65.3762.1749.4181.29<0.001
EC0.749.410.738.840.021
FSII3.1354.203.5746.290.006
LD28.9255.1235.0844.00<0.001
TPA240.75115.09725.6455.93<0.001
TPAR0.4311.820.4012.97<0.001
TPT1.9556.293.5538.11<0.001
AaV27.9424.3429.6819.270.004
AHue57.4310.9655.9110.120.007
ACh51.7713.3452.888.730.048
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

González-López, J.; Rodríguez-Moar, S.; Silvar, C. Correlation Analysis of High-Throughput Fruit Phenomics and Biochemical Profiles in Native Peppers (Capsicum spp.) from the Primary Center of Diversification. Agronomy 2021, 11, 262. https://doi.org/10.3390/agronomy11020262

AMA Style

González-López J, Rodríguez-Moar S, Silvar C. Correlation Analysis of High-Throughput Fruit Phenomics and Biochemical Profiles in Native Peppers (Capsicum spp.) from the Primary Center of Diversification. Agronomy. 2021; 11(2):262. https://doi.org/10.3390/agronomy11020262

Chicago/Turabian Style

González-López, Jorge, Simón Rodríguez-Moar, and Cristina Silvar. 2021. "Correlation Analysis of High-Throughput Fruit Phenomics and Biochemical Profiles in Native Peppers (Capsicum spp.) from the Primary Center of Diversification" Agronomy 11, no. 2: 262. https://doi.org/10.3390/agronomy11020262

APA Style

González-López, J., Rodríguez-Moar, S., & Silvar, C. (2021). Correlation Analysis of High-Throughput Fruit Phenomics and Biochemical Profiles in Native Peppers (Capsicum spp.) from the Primary Center of Diversification. Agronomy, 11(2), 262. https://doi.org/10.3390/agronomy11020262

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop