Characterization of yellow pitahaya (Selenicereus megalanthus Haw.) genotypes under two productive systems in Colombia

Abstract Selenicereus megalanthus Haw. It is an exotic fruit tree, with productive and nutritional potential. In Colombia, there is a great phenotypic and genotypic diversity, but its genetic studies are scarce. The objective was to characterize morphologically 15 selected yellow pitahaya genotypes, under two productive systems in the open field and under cover, in the municipalities of Miraflores and Zetaquira, in Boyacá, Colombia. Quantitative characters were evaluated: plant height (PH), number of vegetative sprouts (NVS), sub-sprouts (SS), longest sprouts length (LSL), distance between areoles (DBA), width of the ribs in the apical region (WRA), width of the ribs in the middle region (WRM), width of the ribs in the basal region (WRB), height of undulations between successive areoles in a rib (HUA), number of spines per areole (NSA) and longest spine length (LSP). The results showed under the two productive systems and the evaluated localities that the variables with the highest coefficient of variation (greater than 90%) were the number of sub-sprouts, height of the undulations between successive areoles (HUA) and the longest spine length (LSP). High positive correlations were obtained between the distances areoles, the width of the ribs and the length of the spines (r>0.7). The conglomerate showed that the characteristics that define the groupings are height of the plant, the texture of the cladodes, the width of the ribs and the height of the undulations. Characters associated with the shoots and cladodes were identified, which directly influence the vegetative propagation and therefore the yield of the yellow pitahaya.

of Leovika et al. (2019) determined in Pangadaran that differences in the morphological, physiological, and anatomical structures of the plants allowed them to adapt to different agroclimatic conditions.A previous study showed that morphoagronomic characterizations could be used to determine the variation in natural pitahaya populations (Abirami et al., 2021).Sibut et al. (2023) evaluated yield capacity of six superior pitaya genotypes under edaphoclimatic conditions of the Federal district in Brazil, found genotypes with good physical characteristics, performance and adapted to local conditions.
In Colombia, knowledge of this crop mainly comes from the empirical processes of farmers who are motivated by the prices that fruits can fetch in some of the months of the year; however, increased penetration into international markets requires research on the processes of propagation, obtaining elite material, and resistance to biotic and abiotic factors, among other topics (Caetano et al., 2013).However, one of the bigger limitations is the broad morphological variation seen in the vegetative structures, which leads to confusion in identifying each species, with a lack of consensus (Gómez-Hinostrosa et al., 2014), where classification is mainly based on the number of areola ribs, the contour of the stem, the relative firmness of the stem and the size and color of the fruits; in addition, various studies on domesticated cactus species have demonstrated variations in fruit characteristics related to the domestication process, resulting in a lack of a taxonomic database (Jalgaonkar et al., 2022).Morphological and molecular characterization studies carried out on the germplasm at the Department of Boyacá have shown that there is genetic diversity in the province of Lengupá, which is the basis for the establishment of any selection process that would identify elite materials that respond to the needs of the farmer, producers, and consumers.Selenicereus megalanthus, or yellow pitahaya, a native fruit tree, constitutes a key product for the fruit and vegetable sector, with significant demand for its flavor, appearance, quality, and nutraceutical properties (Morillo et al., 2021(Morillo et al., , 2022a, b), b).In addition, it has market potential, both domestically and internationally.However, the prospects for the international market require research to raise the fruit quality and to supply elite genotypes for planting and managing the production chain with added value.
The municipality of Miraflores has the productive potential to make pitahaya an economically profitable crop, but first, it must generate strategies to promote cultivation in a technical way to position itself with the highest quality and develop research and technology for the yellow pitahaya production chain, therefore, two productive systems of yellow pitahaya have been implemented in the region in the open field and under cover, as alternatives to reduce the incidence and severity of phytosanitary problems, increase the quality of the fruit and yields, however, they do not studies have been carried out to determine the phenotypic differences of the productive systems.Therefore, the objective of this research was to morphoagronomically characterize yellow pitahaya genotypes grown in the open field and under cover in the municipality of Miraflores and Zetaquira, Boyacá in Colombia in order to identify those that present the best

Introduction
Yellow pitahaya is a member of the Cactaceae family and native to the tropical forest regions of southern Mexico, Central America, and northern South America.The crop adapts to different ecological conditions ranging from very dry regions to wet one receiving more than 3.500 mm of rainfall per year (Goenaga et al., 2020).It is an exotic fruit tree with great popularity in the national and international market due to the appearance of its fruits and its physicochemical, nutritional characteristics and bioactive components, being widely used for its organoleptic characteristics and for its added commercial value (Verona-Ruiz et al., 2020).
It is considered a colorful and attractive food, it arouses the curiosity of the palate favored by the natural color, in addition to the considerable sugar content of the fruit.Another relevant aspect is the composition of bioactives such as phenolic compounds, flavonoids and vitamin C, with quantitative concentrations, in addition to the high fiber composition (Al-Mekhlafi et al., 2021).When consumed, it can provide functional actions, supporting immunity and performing antimicrobial, antioxidant, hepatoprotective, hypoglycemic, healing and antiproliferative activities (Quispe Lupuche et al., 2021).
The dragon fruit cultivation area is expanding rapidly in many countries due to its economic potential and nutritional benefit (Hossain et al., 2021), attracting interest from countries such as the United States, Australia, Southeast Asia, Israel, and other regions (Jalgaonkar et al., 2022).Pitahaya mayor importers are the United States, Japan, European Union and Canada, while major exporters are: Israel, Mexico and Nicaragua, in addition, Colombia participate with 38% of the international exportations (AGRONET, 2022).Selenicereus production reaching 9.93 t/ ha during 2020, the principal production departments include Huila, Boyacá, Santander and Valle del Cauca (AGRONET, 2022).Boyacá is one of the departments of Colombia, where the cultivation of yellow pitahaya has the productive potential and the edaphoclimatic conditions for its production, in addition to the fact that there are producer associations that encourage its planting (Morillo et al., 2022b).Despite the productive potential, a limiting factor in the development of this crop in Colombia is the incidence of pest and diseases, low quality of the fruit, the technological level, mechanical damage during shelf life and the lack of cultivated material which generates considerable economic losses (González-Trujillo et al., 2019;Morillo et al., 2021;Álvarez-Herrera et al., 2023).
Worldwide, morphological characterization studies on yellow pitahaya have found great morphological and genetic heterogeneity in many of this fruit's characteristics, such as sweetness, size, shape, color, and number of bracts, resulting from intra-or interspecific hybridizations between different cultivated and wild materials, which make it difficult to raise quality standards for the export market, posing serious problems when it comes to improving yield and post-harvest shelf-life (Abirami et al., 2021).Studies carried out by Abirami et al. (2021) through morphological, biochemical, and molecular characterizations allowed the correct identification of four Hylocereus species.The authors Brazilian Journal of Biology, 2023, vol.83, e274152

3/12
Characterization of yellow pitahaya genotypes agronomic characteristics and can be a good production alternative under either of the two production systems.

Vegetal material
Fifteen yellow pitahaya genotypes were evaluated, which were selected in an in situ morphoagronomic characterization study carried out in the main producing municipalities of the Lengupá province: Berbeo, Zetaquira, Paéz, San Eduardo and Miraflores (Morillo et al., 2021).The 15 genotypes were established in two experimental plots in the municipalities of Miraflores and Zetaquira under two production systems in the open field (8 genotypes) and under cover (7 genotypes), the experimental design was completely randomized, with 10 repetitions per genotype and in where the treatments corresponded to the 15 selected yellow pitahaya genotypes (Table 1, Figure 1).

Morphoagronomic characterization
The quantitative descriptors were evaluated in situ: plant height (PH), number of vegetative sprouts (NVS), sub-sprouts (SS), longest sprouts length (LSL), distance between areoles (DBA), width of the ribs in the apical region (WRA), width of the ribs in the middle region (WRM), width of the ribs in the basal region (WRB), height of undulations between successive areoles in a rib (HUA), number of spines per areole (NSA) and longest spine length (LSP).These descriptors have been used successfully in the evaluation of pitahaya germplasm in Colombia (Morillo et al., 2021).

Data analysis
For the analysis of the data obtained, an analysis of descriptive statistics was initially carried out.Then the normality assumptions for the parametric analyzes were verified and the analysis of variance (ANOVA) was performed.Pearson's correlation was estimated for variables, and the significance was evaluated using a t-test (the null hypothesis was H0: r = 0; 5% significance), and plotted using R package "corrplot": Visualization of a Correlation Matrix (Version 0.84) (Taiyun and Simko, 2017).A principal component analysis was performed between the variables using the correlation matrix between the characteristics, which were plotted on a two-dimensional plane to group the accessions with R Core Team Software (2020).For the clustering analysis, the Euclidean distance and Ward's minimum distance were taken into account using the algorithms in the Factoextra package of the R program (Kassambara and Mundt, 2020).
For the joint analysis of the qualitative and quantitative descriptors, a factorial analysis of mixed data was carried out with the Factoextra package of the R program.In addition, a dendrogram was generated using a hierarchical clustering method and the Euclidean distance of the minimum variance of Ward with the FactoMineR package (Le et al., 2008).

Morphoagronomic characterization using quantitative descriptors in the system under cover in the municipalities of Miraflores and Zetaquira
The evaluation of the descriptors associated with the seed and sprouts of the 7 yellow pitahaya genotypes in the system under cover in the municipalities of Miraflores and Zetaquira showed that the variables that presented a higher percentage of variation were: sub-sprouts (SS) with a coefficient of variation greater than 94%, followed by the height of the undulations between successive areoles in a rib (HUA) in the municipality of Zetaquira with 88.12% and in Miraflores the length of the longest sprouts (LSL) with 87.41%.While the least variable descriptor was the height of the plant (PH) with approximately 30%, in both municipalities.The average height of the average plants was 137 to 138 cm, the number of sprouts was 2 sprouts per plant (NVS), the number of thorns (NSA) per areola was 1 to 2, and the length of the largest thorn (LSL) was 0.2mm (Table 2).
The Pearson correlation analysis (p≤0.05) in the system under cover in the municipalities of Miraflores and Zetaquira showed positive correlations between all the variables, in the municipality of Miraflores high, significant and positive correlations were observed between rib width and width of the rib in the apical region (WRA) and longest spine length (LSP) with r= 0.93, the number of spines per areola (NSA) and longest spine length (LSP) (r= 0.91) and  distance between areoles (DBA) and number of spines per areola (NSA) with r= 0.9 (Figure 2a).
In the municipality of Zetaquira, high, significant and positive correlations were observed between the distance between areoles (DBA) and the width rib of mid region (WRM) (r= 0.93), between the number of spines per areole (NSA) and the distance between areoles (DBA) (r= 0.92) and between number of spines per areola (NSA) and width rib in apical region (WRA) (r= 0.92) (Figure 2b).
The analysis of mixed factors in the system under cover in the municipality of Miraflores showed that 82.81% of the total variance is explained in the first two components (PC1=70.89%and PC2= 11.92%) (Figure 3a).The variables that make the greatest contribution to the variation of CP1 are plant height, number of spines per areola, and width of the rib in the apical region.While the number of sub-sprouts and the width of the rib in the basal region contribute to CP2.A cluster is observed between genotypes 3, 4 and 6 (Figure 3b).
The analysis of mixed factors in the system under cover in the municipality of Zetaquira showed that 87% of the total variance is explained in the first two components (PC1= 72.95% and CP2= 14.05%) (Figure 4a).The variables that make the greatest contribution to the variation of CP1 are distance between areoles, longest spine length, width of the rib in the basal region.While the number of sprouts and sub-sprouts and the width of the rib in the apical region contribute to CP2.A cluster is observed between genotypes 2, 3 and 6 (Figure 4b).
The cluster analysis of the 7 genotypes evaluated in the municipalities of Miraflores and Zetaquira in the system under cover grouped the genotypes into four clusters (Figure 5a, Figure 5b).Genotype 1 remained independent in a cluster in both municipalities, this genotype was characterized by presenting 1 to 2 sprouts, and without sprouts.In a cluster genotypes 2 and 3 are grouped in both municipalities, however, in Miraflores genotype 4 was added to this cluster, while genotype 6 was added in Zetaquira.In Miraflores genotype 6 remained in an independent cluster because This genotype stood out for presenting the lowest plant size with 113.60 cm, the shortest sprouts length with 22.65 cm and the shortest  height of the undulations between areolas with 0.15 cm.While in Zetaquira genotype 7 was grouped into a cluster, this genotype was characterized by presenting the highest number of sprouts with an average of 4.56 sprouts per plant, the highest number of sub-sprouts with approximately 2 sub-sprouts and a height of the undulations between successive areoles in an upper rib with 0.34 cm.A high incidence of pests and diseases, smooth cladodes texture, concave margin shape, and brown spine color were observed in both municipalities.

Morphoagronomic characterization using quantitative descriptors in the open field in the municipalities of Miraflores and Zetaquira
The evaluation of the descriptors associated with seed and sprouts of the 8 genotypes of yellow pitahaya in the open field in the municipalities of Miraflores and Zetaquira showed that the variables that presented a higher percentage of variation were: sub-sprouts (SS) with a coefficient of variation greater than 97%, followed by the height of the undulations between successive areoles in a rib (HUA) in Zetaquira with 97.48% and longest spine length with 89.33%.While the least variable descriptor was the height of the plant (PH) with approximately 35%, in both municipalities (Table 3).The average height of the average plants is 140-156 cm, the number of sprouts is 2-3 sprouts per plant (NVS), the number of spines (NSA) per areole is 1-2, and the length of the largest spine (LSP) between 0.26 to 0.52 mm.
The Pearson correlation analysis (p≤0.05) in the open field system in the municipalities of Miraflores and Zetaquira showed positive correlations between all the variables, in the municipality of Miraflores high, significant and positive correlations were observed between the distance between areoles (DBA) and number of spines per areola (NSA) with r = 0.89, between the number of spines per areola (NSA) and longest spine length (LSP) (r = 0.84) and the apical rib width (WRA) and number of spines per areola (NSA) with r= 0.83 (Figure 6a).In the municipality of Zetaquira, high, significant and positive correlations were observed between the distance between areolas (DBA) and apical rib width (WRA) (r = 0.92), between the number of spines per areola (NSA) and rib width. of basal thickness (ACGB) (r = 0.91) and between medium thickness rib width (ACGM) and basal rib width (WRB) (r = 0.9) (Figure 6b).
The analysis of mixed factors in the system under cover in the municipality of Miraflores showed that 66.19% of the total variance is explained in the first two components (PC1= 48.19% and CP2= 18.00%) (Figure 7a).The variables that make the greatest contribution to the variation of CP1 are the basal rib width, the average width of ribs and number of spines per areola.While the number of sub-sprouts, the height of the plant and the height of the undulations between areoles contribute to CP2.A cluster is observed between genotypes 1, 5 and 7 (Figure 7b).
The analysis of mixed factors in the system under cover in the municipality of Zetaquira showed that 74.66% of the total variance is explained in the first two components (PC1=61.74%and PC2= 12.92%) (Figure 8a).The variables that make the greatest contribution to the variation of the CP1 are the height of the plant, the length of the longest sprouts and the basal ribs width.While the length of the largest spine and the number of sub-sprouts contribute to CP2.A cluster is observed between genotypes 8, 2 and 3 and the cluster of genotypes 6 and 1 (Figure 8b).
The cluster analysis of the 8 genotypes evaluated in the Miraflores and Zetaquira municipalities in the open field system grouped the genotypes into four clusters (Figure 9a, Figure 9b) in the Miraflores municipality and into three clusters in the Zetaquira municipality.In Miraflores genotypes 6 and 2 each formed a cluster.Genotype 6 was characterized for being the genotype with the highest  Characterization of yellow pitahaya genotypes number of sub-sprouts with an average of 2.8 and a higher than average plant height.Genotype 2 was characterized by presenting the lowest number of sub-sprouts with 0.70 and an average plant height of 164.90 cm.
Smooth cladodes texture, concave margin shape, and brown spine color were observed in both municipalities.In the municipality of Zetaquira, genotypes 1 and 6 were grouped because they presented the lowest size with 129.00 and 107.70 cm respectively, the lowest number of sprouts and sub-sprouts with 2.0 and 0.2 sprouts per plant.The cluster formed by genotypes 2, 3 and 8 presented genotypes with plant height of 132.00 cm and 146.00 cm and width of ribs with a basal thickness between 1.69 and 1.74 cm.Finally, genotypes 4.5 and 7 presented the tallest plants with heights of 153.20, 165.30 and 160.00 cm, average number of sprouts and sub-sprouts of 3.43 and 1.30 sprouts respectively.

Discussion
Pitahaya is an exotic fruit that has been cultivated on family farms since ancient times (Tel-Zur, 2022) and currently, given the physiological characteristics of cacti, such as adaptability to drought, they have made it an economically profitable crop option in many parts of the world., in addition to its appearance, flavor and consistency of the fruit that contribute to its economic value.This has encouraged a commercial market for exotic fruits in Europe, Asia and the Americas (Jalgaonkar et al., 2022), with these also being highly valued in national markets (González-Trujillo et al., 2019).It is mainly consumed as dried fruit and can be used in fruit salads, soft drinks, jams, beer, wine, among others (El Salous et al., 2020).The pitahaya fruit, given its physical and chemical features offer a wide market of exportation as a dry fruit, or frozen pulp and its ample and diverse forms of consumption, constitute productive and viable alternatives as long as the internal market is strengthened, promoting the habit of consumption (Quispe Lupuche et al., 2021).
The importance of the cultivation of S. megalanthus for their commercialization has allowed for the development of a wide variety of studies about the quality of the fruit, commercialization, nutritional composition, cultivation, phenology, organic fertilization, physical chemistry, genetics, postharvest management, morphology and production systems (Tel-Zur, 2022).
However, even today, it is difficult to clearly identify the species and varieties of dragon fruit, due to its high rate of inter-and intra-specific hybridization and the genotype-environment interaction that has led to some confusion at the national level taxonomic worldwide (Tel-Zur, 2022).The phenotypic and genetic variability in morphological characteristics associated with the plant such as height, width ribs of the cladode, distance between areoles, number of sprouts, length of the spines and sweetness, size and color of the fruit are the product of this natural hybridization process (Tel-Zur et al., 2011).makes it difficult to increase the quality standards for the exportation market because it posed serious problems in determining their performance in handling and shelf life (Álvarez-Herrera et al., 2023).
Conventionally, morphological traits had been used to differentiate plant germplasm/ species and to elucidate their genetic relationship (Sibut et al., 2023).As in this study, high values were observed in the coefficients of variation in the characteristics associated with the cladode and the sprouts, in the evaluated locations.Several authors suggest that cladode characters such as cladode width (mm), distance between areoles (mm), number of spines, length of areoles (mm), margin ribs of cladode and its waxiness could be used for identification of Hylocereus and Selenicereus spp.(Rodrigues et al., 2021;Mejía et al., 2013).Abirami et al. (2021) in the morphological, biochemical and molecular characterization of three species of Hylocereus from India, found that cladode characters such as number of spines (Verona-Ruiz et al., 2020;Oliveira et al., 2020), length of areoles (mm) as 1-4, margin ribs of cladode (convex or concave) and its waxiness (weak or strong white waxy or light waxy) could be used for identification of three Hylocereus spp.In this study, the genotypes showed the shape of margin between areolas between concave and right, contrary to that reported in other studies on different species of Hylocereus, where the predominant forms were convex, with the presence of wax (Abirami et al., 2021).The cladode width and distance between areoles corresponded well with earlier studies from México (Castillo et al., 2005), India (Abirami et al., 2021), Puerto Rico and Colombia (Morillo et al., 2021;Betancur et al., 2020;Goenaga et al., 2020).
Multiple correspondence analyzes showed that the characteristics that most contributed to the phenotypic variation observed were the distance between areoles, the width of the ribs in the apical, middle, and basal region, and number of spines.Results that had already been reported in other studies in different countries where found high variation in characteristics of agronomic importance, even within the same species of Hylocereus and Megalanthus spp.(Goenaga et al., 2020;Abirami et al., 2021;Mejía et al., 2013), which is favorable for future breeding studies.The availability of this information would be of great assistance in developing an appropriate method for the cultivation of a particular species (Callejas-Chavero et al., 2021) Among pitahaya propagation techniques, cutting stands out, as it enables the uniformity and quality of seedlings by cloning selected genotypes from mother plants in addition to eliminating the youth period (Zem et al., 2015).Rodrigues et al. (2021) evaluated the ideal size of white-fleshed red pitahaya cladodes and the time of their collection for crop propagation by cutting in view of the physiological quality of the produced clonal plants to enable more appropriate cultural management and increase the cultivation area.Analyzing only the cladode size, a positive relationship was observed between the number of sprouts and sprouting percentage in the cladodes of the greatest studied lengths.this result corroborates the analyses carried out by Moreira et al. (2017), who studied the determination of the cutting length for the production of clonal pitahaya plants (Hylocereus costaricensis) in protected environments and observed that clonal plants formed from cladodes of greater lengths presented higher reserves and stored nutrient availabilities (Delgado and Yuyama, 2010), mainly carbohydrates, which directly contribute to the formation and growth of sprouts and roots (Nicoloso et al., 2001).
Thus, in cases where there is greater availability of material, the use of larger cladodes is recommended, since according to Balaguera et al. (2011), a greater cladode length allows a plant to reach the apex of the tutor and begin emitting secondary branches more quickly, with adult plants forming earlier, initiating the reproductive stage.The size of cuttings is related to their nutritional condition, with the amount of reserves and the number of buds varying according to the length used; normally, cuttings with greater lengths may have higher carbohydrate and endogenous auxin contents (Mayer et al., 2002).Additionally, it is known that the presence of leaves and buds in cuttings plays a key role in the formation of new root systems, and these factors are responsible for the production of assimilates and substances such as auxins, which is synthesized in these locations (Taiz et al., 2017); this analysis shows the high capacity of pitahaya wood cladodes to remain alive.
In general terms, the morphological evaluation of the yellow pitahaya genotypes under the two productive systems in the municipalities of Miraflores and Zetaquira, show differences in the expression of the phenotypic characteristics associated with the cladodes and the sprouts, the above may be due to the factors agroclimatological present in the evaluated localities.Results similar to those obtained by Morillo-Coronado et al., 2021 in the morphoagronomic characterization of different genotypes in five municipalities of the province of Lengupá.
In Colombia, the morphoagronomic characterization studies of the yellow pitahaya germplasm have shown the existence of variability which can be used within improvement programs, considering the productive systems that farmers currently have, since this in one way or another is influencing the expression of desirable characters.

Conclusions
The morphoagronomic characterization of yellow pitahaya genotypes showed morphological differences in shoots and cladodes according to the site where they were evaluated.The existence of phenotypic variation allows carrying out selection processes that lead to the identification of elite genotypes that meet the needs of producers and are an economically profitable alternative.

Figure 1 .
Figure 1.Experimental plots a. System under cover in the municipality of Zetaquira.b.Open field system in the Municipality of Zetaquira.c. System under cover in the municipality of Miraflores.d.Open field system in the Municipality of Miraflores.

Figure 2 .
Figure 2. Pearson correlation analysis (p≤0.05) between the quantitative variables evaluated in the genotypes of yellow dragon fruit under cover.a. Miraflores.b.Zetaquira.

Figure 3 .
Figure 3. Analysis of mixed factors, considering quantitative and qualitative variables evaluated in the seven genotypes of yellow pitahaya in the Miraflores municipality under cover.a. Contribution of the variables to the first two components.b.Distribution of the genotypes in the first two components.

Figure 4 .
Figure 4. Analysis of mixed factors, considering quantitative and qualitative variables evaluated in the seven genotypes of yellow pitahaya in the municipality of Zetaquira under cover.a. Contribution of the variables to the first two components.b.Distribution of the genotypes in the first two components.

Figure 5 .
Figure 5. Analysis of hierarchical clusters of the yellow pitahaya genotypes of the municipalities of Miraflores and Zetaquira in the productive system under cover.a. Miraflores.b.Zetaquira.

PH:
Plant height; NVS: Number of vegetative sprouts; SS: Sub-sprouts; LSL: Longest sprouts length; DBA: Distance between areoles; WRA: Width of the ribs in the apical region; WRM: Width of the ribs in the middle region; WRB: Width of the ribs in the basal region; HUA: Height undulations between successive areoles in the ribs; NSA: Number of spines per areole; LSP: Longest spine length; CV: Coefficient of variation.

Figure 6 .
Figure 6.Pearson correlation analysis (p≤0.05) between the quantitative variables evaluated in the open field yellow pitahaya genotypes.a. Miraflores.b.Zetaquira.

Figure 7 .
Figure 7. Analysis of mixed factors, considering quantitative and qualitative variables evaluated in the seven genotypes of yellow pitahaya in the municipality of Miraflores in the open field.a. Contribution of the variables to the first two components.b.Distribution of the genotypes in the first two components.

Figure 8 .
Figure 8. Analysis of mixed factors, taking into account quantitative and qualitative variables evaluated in the seven genotypes of yellow pitahaya in the municipality of Zetaquira in the open field.a. Contribution of the variables to the first two components.b.Distribution of the genotypes in the first two components.

Figure 9 .
Figure 9. Analysis of hierarchical clusters of the yellow pitahaya genotypes of the municipalities of Miraflores and Zetaquira in the open field production system.a. Miraflores.b.Zetaquira.

Table 1 .
Geographical location of the selected yellow pitahaya genotypes.

Table 2 .
Quantitative descriptors evaluated under cover in the municipalities of Miraflores and Zetaquira.
PH: Plant height; NVS: Number of vegetative sprouts; SS: Sub-sprouts; LSL: Longest sprouts length; DBA: Distance between areoles; WRA: Width of the ribs in the apical region; WRM: Width of the ribs in the middle region; WRB: Width of the ribs in the basal region; HUA: Height undulations between successive areoles in the ribs; NSA: Number of spines per areole; LSP: Longest spine length; CV: Coefficient of variation.

Table 3 .
Quantitative descriptors evaluated in the open field in the municipalities of Zetaquira and Miraflores.