Revista Mexicana Ciencias Agrícolas   volume 12   number 2   February 15 - March 31, 2021

DOI: https://doi.org/10.29312/remexca.v12i2.2848

Article

Planting density in verdolaga growth

Erika Lagunes-Fortiz

Clemente Villanueva-Verduzco1

Edgar Ricardo Lagunes-Fortiz1

Erika Janet Zamora-Macorra1

Norma Ávila-Alistac1

Evert Villanueva-Sánchez2

1Chapingo Autonomous University. Highway México-Texcoco km 38.5, Chapingo, Texcoco, State of Mexico. ZC. 56230. Tel. 595 9521500, ext. 6133. 2CONACYT-Chapingo Autonomous University. National Laboratory for Research and Food and Forestry Service. Highway México-Texcoco km 38.5, Chapingo, State of Mexico, Mexico. ZC. 56230.

§Corresponding author: elagunesf@chapingo.mx.

Abstract

In Mexico there are more than 358 species of plants, of which only tender leaves and stems are consumed, known as ‘quelites’. Given the productive, economic and nutritional potential of these species, including verdolaga, factors that affect their growth and development need to be determined. Verdolaga is an alternative as a vegetable due to its nutritional, productive and economic potential; its yield will be determined mainly by the development of stems and leaves, so the knowledge generated about the factors that affect its growth is of great agronomic importance. This research assessed the effect of three planting densities on verdolaga growth. Data for three phenological stages of verdolaga development (vegetative development, vegetative maturity and flowering) were considered to determine whether there were significant differences between the mean treatments evaluated by an Anova and to determine whether variance homogeneity existed, the Levene method was carried out. Plants planted at a high planting density (295 plants m-1 linear) had higher height, smaller stem diameter, fewer leaves (large and small) and lower branch development (branch length) compared to those planted at medium and low planting densities (250 and 134 plants m-1 linear); however, the number of leaves on the branches depends on the level of the canopy where they are located. Planting density is a contributing factor in determining the growth of the ‘Chapingo’ verdolaga.

Keywords: phenological stages, quelites, verdolaga.

Reception date: January 2021

Acceptance date: March 2021

Introduction

In Mexico there are more than 500 species of wild edible plants, whose use dates back to pre-Hispanic times given its great cultural, culinary and nutritional relevance, including quelites, nahuatl ‘quiltil’, of which only tender leaves, stems and flowers are consumed (Bye and Linares, 2002;  Balcázar-Quiñones et al., 2020), as they are an important source of nutrition (Ca, K, Mg, Fe, Mn, Zn, P), protein, dietary fiber, vitamins A, C and bioactive compounds, such as phenols, betalains, flavonoids, alkaloids, omega 3, omega 6, carotenes, ascorbic acid, α-tocopherol, glutathione, etc., which give these plants a high antioxidant capacity, in addition to providing dishes with different fragrance, colors and flavors (Bye and Linares, 2002;  Santiago-Saenz et al., 2019;  Balcázar-Quiñones et al., 2020).

Given their high productive potential and the existence of an ancestral culture of consumption of some of these crops (verdolagas, quintoniles, papalos, epazotes, etc.) in the Mexican population of the center of the country, where there are 25 million inhabitants in The City and The State of Mexico alone (INEGI, 2015), the Educational and Research Institutions must develop technology, varieties , intensive production systems and agronomic practices that provide a comprehensive knowledge of the different species of quelites and that help to promote the optimal use of their economic potential and as a food plant.

In Mexico, commercial verdolaga production is concentrated in the states of Morelos (171 ha), Mexico City (158 ha), State of Mexico (45 ha) and Baja California (46 ha) with yields ranging from 12 to 14.86 udm ha-1 and a production of 5 474.15 t, with a production value of just over 26 361.30 thousand pesos (SIAP, 2019).

In our country, there are two well-defined producing regions (contrasting in climate) with local creole varieties domesticated under cultivation: the region of Xochimilco-San Gregorio-Mixquic (temperate climate) and that of Cuautla, Morelos (warm climate) (Villanueva and Ramírez, 2003). The active and continuous domestication of verdolaga in the different producing regions has allowed the obtaining of genotypes with thicker leaves, habits of erect growth, greater precocity and resistance to diseases, such as the crop ‘Chapingo’ of verdolaga, whose bromatological analyses have shown that it also contains high levels of iron, magnesium, omega 3 and omega 6 (Solís et al., 2016).

Studies conducted under greenhouse (without heating) during the winter cycle indicate that it is possible to produce up to 20 kg m-2 of vegetables in 50 days (200 t ha-1), with a population density of 500 m-2 (Villanueva and Ramirez, 2003). The value of the winter verdolaga harvest can reach up to $8.00 per kilogram to the producer, although it drops dramatically at times of the year of its natural production (Villanueva and Ramirez, 2003; Mera et al., 2010).

The study of the growth of verdolaga and the factors that affect it is of great importance, since its yield will be determined mostly by the development of the vegetative part of the plant, which in turn will depend on multiple factors, some of them typical of the genotype, others of the environment and various conditions of crop management.

Recently Montoya-García et al. (2017) assessed the effect of N, P and K fertilization on verdolaga production, found that the height and yield of verdolaga increase in response to the application of applied nitrogen and the harvest date. Montoya-García et al. (2018) report that it is possible to maximize the nutraceutical content, as well as the concentration of fatty acids and antioxidants with fertilizations with nitrogen, phosphorus and potassium, so factors such as fertilization not only influence the growth of the plant, but also on the potentiation of its properties for the benefit of health.

Knowledge about crop growth forms the basis for the optimal planning and execution of cultural work to be carried out during the production process, in order to obtain the best yields. The genotype, planting date, planting density, locality and climate in which plants develop are factors that can influence the amount of foliar tissue, growth habit, part of the plant and development of other structures.

The most important factors affecting the duration of the stages of verdolaga development include the genotype (inside and between populations), the environment (climate and soil) and its interactions, which can vary the growth habit and precocity of the plant. Some other characteristics of the plant such as the habit of growth, flavor, phyllotaxis and stem, may vary depending on whether the plant develops as weeds or cultivated (Villanueva and Ramírez, 2003). Plant growth involves at the physiological level a series of changes and irreversible reactions in the plant, on which the final agronomic characteristics (phenotype) and the potential yield of the different genotypes will depend.

The evaluation of the growth of a plant can be carried out through direct measures (plant height, stem diameter, number of leaves, number of branches, foliar area, fresh or dry mass, etc.) and indirect (net assimilation rate, crop growth rate, relative growth rate, etc.) (Rodríguez, 2000).

Planting density is a factor that can influence crop yield such as verdolaga, as it interferes with plant physiological events and directly with accumulation of dry matter from the different organs of the plant, and also determines competition between plants by light, water and nutrition (Rodríguez, 2000; Ayala et al., 2004). Based on the above, the objective of this research was to evaluate the effect of three planting densities on the growth of the crop 'Chapingo' of verdolaga, as well as characterize its dynamics and habit of growth.  

Materials and methods

The experiment was established in a batch of the Experimental Agricultural Field of the Chapingo Autonomous University, Mexico at 19° 29’ 27” north latitude and 98° 52’ 23” west longitude, at an altitude of 2 240 meters above sea level (msnm). The modifications made by Garcia (1988) to the Köppen climate classification system describe the site of the establishment of the experiment as a C(W)(w)b(i)g, which means that it is the driest of temperate subhumid with rains in summer, an average annual temperature of 15.6 °C  and -3° to 18 °C in the coldest month (January), a cool summer, the temperature of the warmest month is less than 22 °C (may), with little annual thermal oscillation ranging from 5° to 7 °C. The soil is classified as Mollic Ustifluvent, with a depth of 150 cm and a slimy clay loam texture on the Horizon Ap (Survey Soil Staff, 2010).

40 furrows ±1 of 30 m long and a distance between furrows 0.6 m with verdolaga of the crop ‘Chapingo’ were planting. The sampling sites were established at random, so that each experimental unit had a different planting density, it is delimited 1 m for each experimental unit and it is counted the number of plants in each experimental unit, leaving three different planting densities (treatments): high (295 plants. linear meter), medium (250 plants m-1 linear) and low (134 plants m-1 linear). In each experimental unit, data were taken to 20 plants (repeats).

The following variables were evaluated every 4 days from planting and up to 48 days after planting (dds). Variables studied were total plant height (AT), was measured in centimeters (cm), from ground level to the end of the main stem; number of branches (NR), the branches that were inserted directly into the main stem were counted; stem diameter (DT), was measured in millimeters (mm) with the help of a Vernier, in the second internode of the plant, just below the third knot; number of large leaves and small leaves (NHG and NHCH, respectively), the large and small leaves of the plant were counted, considering as large leaves those between 1.5 to 2 cm wide and 2.5 to 3 cm long. Average length of the branches of the n knot (LPRn), was taken the average, in centimeters, of the length of the two branches that were inserted into the n knot of the main stem, numbering the knots from the base to the tip of the plant. The average number of leaves of the n knot branches (NHPn), refers to the average number of leaves of each of the branches that were inserted into the n knot of the main stem, numbering the knots of the base to the tip of the plant.

For statistical analysis, data were taken for three phenological stages of verdolaga development: vegetative development (19 dds; 52.1 Units Heat, UC), vegetative maturity (31 dds, 100.65 UC) and flowering (39 dds, 39 142.55 UC) of verdolaga of the crop ‘Chapingo’. The data was analyzed with the SAS System for Windows 9.1 (SAS Institute, 2002) package, to determine whether there were statistically significant differences between the averages of the three planting densities evaluated an Anova was performed and to determine whether there was equality of its variances the Levene method was used, both with a significance level of p≤ 0.05.

Results and discussion

Figure 1 shows the growth of the verdolaga in plant height (AT), number of branches (NR) and stem diameter (DT); in Figure 2, the number of large leaves (NHG) and number of small leaves (NHCH), while in Figure 3, the growth of the branches of the main stem at the height of the 1st and 4th knot, counting from the base to the tip of the plant. The growth dynamics represented in these figures show that it is between 20 and 40 dds where the growth of the verdolaga is accelerated so, although the variables were evaluated every 4 days, the Anova and the Levene test of variance comparison was performed with the data obtained in the phenological stages covering that period (vegetative development, vegetative maturity and flowering).

Plant height

Throughout the life cycle of the plant, the height increases constantly and exponentially, particularly between 20 and 40 dds, which match with the beginning of the vegetative development stage and until the beginning of flowering (Figure 1a). The AT mean comparison test shows that only during the vegetative maturity stage there were statistically significant differences between the means of the three planting densities (p< 0.05, Table 1). The highest AT was obtained in the high planting density (295 plants m linear) relative to the mean and low densities (250 and 134 plants m linear, respectively).

Figure 1. Growth dynamics of a) total plant height (cm), b) number of branches and c) stem diameter (in cm) of the verdolaga in three planting densities (high: 295; medium: 250 and low: 134 plants m linear, respectively).

Table 1. Anova of total height of the plant (AT, in cm), number of branches (NR), stem diameter (DT, in cm), number of large leaves (NHG) and number of small leaves (NHCH) at 19 dds (vegetative development), 31 dds (vegetative maturity) and 39 dds (start of flowering) of verdolaga of the variety ‘Chapingo’ sown in three planting densities (high: 295; medium: 250 and low: 134 plants m linear, respectively).

Planting density

AT

NR

DT

NHG

NHCH

M

Anova (F/Pr value)

M

Anova (F/Pr value)

M

Anova (F/Pr value)

M

Anova (F/Pr value)

M

Anova (F/Pr value)

Vegetative development

High

12.29

0.20/0.8197

9.4

0.45/0.6419

0.45

0.64/0.5335

10.75

2.02/0.1416

41

4.39/0.0169

Medium

11.6

8.65

0.42

10.55

38.4

Low

11.4

9

0.48

12.25

54.3

Vegetative maturity

High

28.35

3.68/0.0313

11.55

2.68/0.077

0.81

2.62/0.0813

11.5

1.31/0.2768

107.15

0.45/0.6376

Medium

25.26

11.6

0.84

12.15

111.45

Low

28.6

12.45

0.9

13.75

114.3

Flowering

High

34.51

2.54/0.0877

12.05

0.66/0.5205

1

0.64/0.5335

10.15

2.79/0.07

227.1

0.04/0.9633

Medium

31.69

12

1.02

12.35

224.8

Low

32.48

11.5

1.05

11.5

222.5

M= average. Bold data group indicate statistically significant differences between them (p< 0.05).

A preliminary analysis of the data revealed that the dispersion of observations at the three planting densities might be different, so a Levene test was performed. It was concluded that the variances are statistically different so planting density is a contributing factor in determining the plant height of the verdolaga (Table 2).

In this sense, Rodríguez (2000) reports that the increase in population density generates the production of longer plants, while Favaro and Pilatti (1997) indicate that tomato plants planted in high planting density, originate a higher proportion of far-red light (730 nm) in their lower strata, affecting phytochrome activity, which is expressed by an increase in the length of internodes and consequently in plant height.

Number of branches


Noda and Gildardo (2008)
 when studying  the growth of Morus alba in different planting densities, report that, plants develop more branches in low planting densities, since little competition by light and the arrangement of direct solar radiation favors the development of lateral branches; however, in this research, no statistically significant differences were found in the number of branches (NR, Table 1) determined by planting density, probably because verdolaga ‘Chapingo’ being a previously selected material, has well-defined structural characteristics, such as erect growth habit, opposite phyllotaxis, long stem, thick stem, succulent green stem and arrangement of its opposite lateral branches, which it preserves very stable.

Between 20 and 30 dds the development of branches is accelerated (period that coincide with the vegetative development stage) and subsequently increases, but to a lesser extent, remaining stable from 30 dds until the end of its life cycle. The verdolaga ‘Chapingo’ develops 8 to 9 branches in the vegetative development stage, 12 branches at the beginning of flowering (Table 2) and up to 14 branches at the end of its cycle (48 dds, Figure 1).

Table 2. Standard deviation of total plant height (AT, in cm), number of branches (NR), stem diameter (DT, in cm), number of large leaves (NHG) and number of small leaves (NHCH) at 19 dds (vegetative development), 31 dds (vegetative maturity) and 39 dds (flowering) of verdolaga of the variety ‘Chapingo’ sown in three planting densities (high: 295; medium: 250 and low: 134 plants m linear, respectively).

Planting density

AT

NR

DT

NHG

NHCH

DS

Levene

(F/Pr value)

DS

Levene(F/Pr value)

DS

Levene(F/Pr value)

DS

Levene (F/Pr value)

DS

Levene (F/Pr value)

Vegetative development

High

4.94

0.06/0.9373

2.6

0.05/0.9536

0.2

0.09/0.9142

2.27

1.64/0.2024

16.62

0.37/0.693

Medium

4.77

2.46

0.17

2.61

18.64

Low

4.42

2.47

0.18

3.7

19.27

Vegetative maturity

High

4.28

2.94/0.061

1.5

0.53/0.5896

0.19

0.77/0.465

2.37

1.8/0.1741

25.4

0.27/0.768

Medium

3.17

1.43

0.1

6.96

25.26

Low

5.28

1.19

0.1

2.68

20.75

Flowering

High

5.35

3.48/0.0377

1.68

0.46/0.6364

0.2

0.09/0.9142

2.77

1.67/0.1978

34.13

0.99/0.3786

Medium

3.49

1.84

0.17

2.91

62.67

Low

3.79

1.47

0.18

3.65

60.59

M= average. Bold data group indicate group of variances with statistically significant differences between them (p<  0.05).

Authors mention that one of the most sensitive factors to density variations is the number of branches per plant, an increase in the population of plants produces a decrease in the number of branches (Ayala et al., 2004; Kakiuchi and Kobata, 2004) which coincides with the results obtained in this research as the number of branches was lower in verdolaga plants sown in high planting density compared to medium and low planting densities.

Stem diameter

From 10 dds the diameter of the stem increases exponentially, from 20 and 40 dds there is a similar accelerated growth in the three planting densities, the growth stops from the 43 dds (capsule development stage) in the high planting density (295 plants m-1 linear) (Figure 1c). No highly significant differences in stem diameter were found between the verdolagas sown in the three planting densities (Table 1).

According to García and Watson (2003), the thickness of the stem depends on the genotype, the environmental and nutritional conditions of the soil, states that in maize the resistance that the plant presents to the acame, depends to a large extent on the stem diameter and that in addition the stem diameter tends to decrease when the planting density is increased, due to competition between the plants.

High planting densities and light competition with weeds cause elongation of stems, longer internodes and taller plants, reducing the thickness of the stems and increasing the chances of acame of the plants, this coincides by observing the results obtained in this research at 48 dds of verdolaga sown in three different planting densities. The diameter of the stem is a very important parameter since in addition to containing the ducts that channel the transport of water and nutrients to the different organs of the plant, it will also influence the resistance to bending of the plant when affected by external factors such as wind or rain. The verdolaga crop ‘Chapingo’, has an erect growth, so an optimal thickening of the stem will favor the development of all the structural organs of the plant.

Number of leaves

In the three planting densities evaluated, from 10 dds to the end of the verdolaga life cycle, the number of large leaves (NHG) is maintained between 10 and 18 large leaves, being at 27 dds where there are the most large leaves (18 leaves) which coincides with the end of the vegetative development stage, then some leaves fall (mainly those located at the base of the lower branches), an event that coincides with the insinuation stage (development of visible floral buttons on the terminal part of the main stem) and flowering, finally towards the end of the plant’s life cycle, the number of large leaves is maintained between 10 and 14 leaves according to planting density (Figure 2a).

The dynamics of the development of the small leaves (NCH) is different (Figure 2b) because from the beginning of the vegetative development stage (18 dds) and until the insinuation stage (35 dds) the number of small leaves is constant in the three planting densities, however from 35 to 48 dds (seed maturation stage) the number of small leaves increases, perhaps due to the demand that represents the production of photo-assimilated for the development of plant reproductive structures during the hint and maturation stages of the seed.

Figure 2. Growth dynamics of a) number of large leaves and c) number of small leaves in three planting densities (high: 295; medium: 250 and low: 134 plants m-1 linear, respectively).

The verdolaga plants planted in the low planting density (134 plants m linear) are the ones that reach the greatest number of large and small leaves, followed by those planted in medium and high planting densities, behavior similar to what was reported by Díaz et al. (1999), which point out that biomass production decreases as density increases, due to increased competition between plants by light, CO2, water and minerals.

Length and number of branch leaves

During the development of the verdolaga a total of 8 knots were developed on the main stem, on which two lateral branches were developed, as well as new leaves on these branches. In this study only those branches located in the 1st and 4th knot were measured (numbering from the base to the tip of the plant) because they are the ones that develop in greater magnitude compared to the branches of the upper knots.

The growth of branches and leaves in the 1st knot (LPR1 and NHP1, respectively) begins to be noticeable from 10 dds (Figure 3a and 3b). On the other hand, the growth of the branches and leaves of the 4th knot (LPR4 and NHP4, respectively) is presented between 15 and 20 dds (Figure 3c and d); that is, during the vegetative development stage.

Figure 3. Growth dynamics of the average length of the two branches of knot a) 1 (LPR1) and c) 4 (LPR4) and number of average leaves of the two branches of knot a) 1 (NHP1) and d) 4 (NHP4) of the verdolaga in three planting densities (high: 295; medium: 250 and low: 134 plants m linear, respectively).

In the LPR1variable, statistically significant differences were found between the three planting densities (p< 0.05), in the vegetative development, vegetative maturity and flowering phases, while for the NHP1 statistically significant differences were found (p< 0.05) between the sowing densities evaluated in the vegetative development stages, vegetative maturity, but not in the flowering stage.  For the LPR4 and NHP4 variables, no statistically significant differences were found between the mean of the three planting densities evaluated with the exception of the LPR4 variable in the vegetative maturity stage (p< 0.05, Table 3).

Table 3. Length and number of leaves of the branches at 19 dds (vegetative development), 31 dds (vegetative maturity) and 39 dds (flowering) of verdolaga of the 'Chapingo' variety sown in three planting densities (high: 295; medium: 250 and low: 134 plants m linear, respectively).

Planting density

LPR1

NHP1

LPR4

NHP4

M

Anova (F/Pr value)

M

Anova (F/Pr value)

M

Anova (F/Pr value)

M

Anova (F/Pr value)

Vegetative development

High

1.13

8.17/0.0008

5.35

14.187/<0.0001

0.3

0.44/0.6488

3.85

0.1/0.9083

Medium

1.28

5.9

0.33

4.1

Low

3.1

10.8

0.26

3.8

Vegetative maturity

High

6.6

17.34/<.0001

12.1

7.86/0.001

5.06

3.29/0.0444

9.9

3.06/0.0546

Medium

8.54

12.35

7.15

10.75

Low

16.1

23.35

5.63

13

Flowering

High

17.22

4.4/0.0169

19.75

2.49/0.0916

12.23

3.04/0.0558

28

2.59/0.0836

Medium

12.49

16.1

10.1

20.4

Low

19.46

24

12.6

21.8

LPRn= average length of the two branches of the n knot; NHPn= number of average leaves of the two branches of the n knot. Bold data group indicate statistically significant differences between them (p< 0.05).

According to the variance comparison test for the variables LPR1, NHPR1, LPR4 and NHPR4, only in a few cases there were statistically significant differences, so the planting density influences the number of branches and leaves that grow on the plant (Table 4). Poulain (1984) notes that high planting densities generate higher plants, but with fewer branches, as happened in the growth of the branches in the 1st and 4th knot (LPR1 and LPR4, respectively) where plants sown in high planting density reached lower length compared to those sown in medium and low densities.

The branches of the first knot of those plants planted in low density contained the largest number of leaves (NHP1) compared to those of the medium and low densities, possibly because their length and availability of space allowed them to produce more leaves. On the other hand, the number of leaves of the branches of the 4th knot (NHPR4) was higher in the verdolaga plants sown in the high planting density, compared to the medium and low densities, probably because the level of the canopy where they were located influenced the development of the leaves to supply the demand for photo-assimilated necessary for the formation of structures typical of the subsequent phenological stages.

Table 4. Standard deviation of branch length (in cm) and number of branch leaves at 19 dds (vegetative development), 31 dds (vegetative maturity) and 39 dds (flowering) of verdolaga of the ‘Chapingo’ variety sown in three planting densities (high: 295; medium: 250 and low: 134 plants m linear, respectively).

Planting density

LPR1

NHP1

LPR4

NHP4

DS

Levene

(F/Pr value)

DS

Levene

(F/Pr value)

DS

Levene (F/Pr value)

DS

Levene (F/Pr value)

Vegetative development

High

1.01

1.62/0.2072

2.6

2.76/0.0721

0.2

0.21/0.8091

2.3

0.41/0.6638

Medium

1.86

3.6

0.21

2.51

Low

2.09

4.29

0.25

2.12

Vegetative maturity

High

5.92

0.61/0.5485

8.96

6.74/0.0023

2.03

3.61/0.0334

3.31

2.69/0.0767

Medium

5.82

14.15

1.96

5.6

Low

4.28

5.86

3.66

2.83

Flowering

High

4.67

5.24/0.0082

8.41

2.83/0.0675

3.06

0.3/0.7406

13.01

0.54/0.5857

Medium

7.18

9.1

3.41

10.76

Low

10.01

14.83

3.57

11.02

LPRn= average length of the two branches of the n knot; NHPn= number of average leaves of the two branches of the n knot. Bold data group indicate group of variances with statistically significant differences between them (p< 0.05).

Conclusions

There are significant differences between planting densities for total height (AT), number of small leaves (NHCH), average length of the two branches of knot 1 and 4 (LPR1 and LPR1, respectively) and number of average leaves of the two branches of knot 1 (NHP1). High planting densities generate taller plants, smaller stem diameter, fewer leaves (large and small) and less branch development (branch length); however, the number of leaves developed in the branches depends on the level of the canopy where they are located.

Cited literature

Ayala, D. I. M.; Rey, B. L. y Durán, Q. C. 2004. Efecto de la densidad de siembra sobre el crecimiento, desarrollo y productividad de dos materiales de palma de aceite Elaeis guineensis Jacq. Rev. Palmas. 25:66-73. https://publicaciones.fedepalma.org/ index.php/palmas/article/view/1067.

Balcázar-Quiñones, A; White-Olascoaga, L.; Chávez-Mejía, C. y Zepeda-Gómez, C. 2020. Los quelites: riqueza de especies y conocimiento tradicional en la comunidad otomí de San Pedro Arriba, Temoaya, Estado de México. Polibotánica. 49:219-242. https://doi.org/10.18387/polibotanica.49.14.

Bye, R. y Linares, E. 2002. Los quelites, plantas comestibles de México. Una reflexión sobre intercambio cultural. Biodiversitas. Comisión Nacional para el Conocimiento y Uso de la Biodiversidad (CONABIO). 31:11-14.

Díaz, L.; Viloria, Z. A. y Arteaga, R. L. 1999. Crecimiento vegetativo del pimentón en función de la densidad de plantas y edad del cultivo. Bioagro. 11(2):69-73.

Favaro, J. C. y Pilati R. A. 1997. Cultivo de tomate. In: cultivos bajo invernaderos. Centro de Publicaciones Universidad Nacional del Litoral y Editorial Hemisferio Sur. 2da (Ed.). Buenos Aires, Argentina. 7-33 pp.

García, E. 1988. Los climas de México. Sistemas de Información Geográfica, SA. México, DF. 1-16 pp.

García, M. y Watson, C. E. 2003. Herencia de la resistencia al acame de raíces en maíz dulce (Zea mays L). Revista Científica UDO Agrícola. 3:24-33.

INEGI. 2015. Instituto Nacional de Estadística y Geografía. Encuesta Intercensal. Dirección General de Estadísticas Sociodemográficas. Encuesta Intercensal 2015. Instituto Nacional de Estadística y Geografía. México, DF. https://www.inegi.org.mx/programas/ intercensal/2015/.

Kakiuchi, J. and Kobata, R. 2004. Shading and thinning effects on seed and shoot dry matter increase in determinate soybean during the seed filling period. Agron. J. 96:398-405.

Mera, O. L. M.; Castro, D. L.; Bye, R. A. y Villanueva, V. C. 2010. Importancia de la verdolaga en México. Universidad Autónoma Chapingo. Chapingo, Estado de México. 28 p.

Montoya-García, C. O.; Volke-Haller, V. H.; Trinidad-Santos, A. and Villanueva-Verduzco, C. 2018. Change in the contents of fatty acids and antioxidant capacity of purslane in relation to fertilization. Scientia Horticulturae. 234:152-159.

Montoya-García, C. O.; Volke-Haller, V. H.; Trinidad-Santos, A.; Villanueva-Verduzco, C. y Sánchez-Escudero, J. 2017. Respuesta de la verdolaga (Portulaca oleracea L.) a la fertilización con NPK. Rev. Fitotec. Mex. 40(3):325-332.

Noda, Y. y Giraldo, M. 2008. Efecto de la densidad de siembra en el establecimiento de morera para su inclusión en sistemas ganaderos. Zootecnia Tropical. 26(3):339-341. http://ve.scielo.org/scielo.php?script=sci_arttext&pid=S079872692008000300039&lng=es&tlng=es.

Poulain, D. 1984. Influence of density on the growth and development of winter field bean (Vicia faba L.), 159-167 pp. In: Hebblethwaite, P. D.; Dawkins, T. C. K. Heath, M. C. and Lockwood. G. (Ed). Vicia faba: agronomy, phisiology and breeding. Martinus Nijhoff/ Dr. Junk, W. Publishers, Netherlands. 316 p.

Rodríguez, L. 2000. Densidad de población vegetal y producción de materia seca. Revista COMALFI. 27(1-2):31-38.

Santiago-Sáenz, Y. O.; Hernández-Fuentes, A. D.; López-Palestina, C. U.; Garrido-Cauich, J. H.; Alatorre-Cruz, J. M. y Monroy-Torres, R. 2019. Importancia nutricional y actividad biológica de los compuestos bioactivos de quelites consumidos en México. Revista Chilena de Nutrición. 46(5):593-605. https://dx.doi.org/10.4067/S0717-75182019000500593.

SAS Institute. 2002. SAS versión 9.1.3 help and documentation. Cary, NC, SAS Institute Inc 633 p.

SIAP. 2019. Servicio de Información Agroalimentaria y Pesquera. Producción Agrícola. Servicio de Información Agroalimentaria y Pesquera, SAGARPA. Ciudad de México. https://nube.siap.gob.mx/cierreagricola.

Soil Survey Staff. 2010. Claves para la taxonomía de suelos. 11a (Ed.). Ortiz-Solorio. C. A. and Gutiérrez-Castorena, M. C. (Trads.). USDA Natural Resources Conservation Service, Washington, DC. USA. 365 p.

Solís, A. G.; Mera, O. M.; Castro, L. D.; Pinzón, L. L.; Noguez, H. R.; Casanova, G. E. y Méndez, V. R. 2016. Conservación y aprovechamiento sostenible de especies vegetales tradicionales de México. Servicio Nacional de Inspección y Certificación de Semillas y Universidad Autónoma Chapingo. Chapingo, Estado de México. 108 p.

Villanueva, V. C. y Ramírez, A, M. 2003. Verdolaga hidropónica y tecnología de producción superintensiva de germinados en invernadero. In: Memorias del XIX Presentación de trabajos de investigación, producción y servicio. Departamento de Fitotecnia. Chapingo, Estado de México. 68-72 pp.