Conventional agriculture applies large inputs of agrochemical products, especially inorganic fertilizers, to maintain global food production, affecting the health of natural and agricultural ecosystems, the human population, and jeopardizing the economy of resource-poor farmers  . Organic agriculture has been proposed as a safe and healthy alternative but the crop yield under this agricultural system is variable due to the influence and interaction of crop type, cultivation system, soil physicochemical characteristics and environmental variation    . Of these factors, soil quality constitutes a key factor in the results obtained with organic agriculture  . Physical soil and biological characteristics are improved with organic fertilizers, like vermicompost, which is the process of bio-oxidation and stabilization of organic matter in controlled conditions of the joint action of specific species of earthworm and microorganisms to accelerate the nutrients transformation   .
This kind of fertilization also increases the colonization of plant roots by arbuscular mycorrhizal (AM) fungi, their abundance and diversity  . This plant-fungi interaction plays a key role in the soil ecosystem services  , because it maintains soil fertility and plant health. But, when intensive farming techniques are used in conventional agriculture systems these soil ecosystem services are degraded   .
The present study quantifies and compares the effects of vermicompost fertilization in the organic system against those obtained using chemical fertilizer in the conventional system on the changes in AM fungi diversity and intra root colonization, as well as the growth, yield and production of green tomato. We hypothesized that organic agricultural system, with vermicompost and leachate, increases AM fungi diversity and tomato plant growth, with similar yield as with agrochemicals.
2. Materials and Methods
2.1. Site Description
The study was conducted in El Chichonal community (Figure 1), 18˚30'52''N and 89˚31'24"W, within the buffer zone of the Calakmul Biosphere Reserve (CBR) and in the municipality of Calakmul, Campeche, Mexico. There are more than 100 ejido communities in the area surrounding the CBR  . The land exhibits a plain topography with an altitude of 250 meters above sea level. The parent
Figure 1. Site description, experimental design and treatments of organic and chemical agriculture systems of green tomato in Calakmul, Mexico. Red: NF = no fertilization, Blue: OTH = organic treatment with vermicompost and leachate, Green: OTV = organic treatment with vermicompost, and Yellow: CST = inorganic fertilizer treatment.
material of the soil is limestone rock, which is very poor in iron (Fe), silica, aluminum (Al), phosphorus (P) and micronutrients (zinc and copper). The main soil type is leptosol rendzico, with a high content of organic matter and is alkaline  . It contains deciduous tropical forest vegetation with different structural appearances that reflect the variation in soil and other biophysical conditions  . A sub humid warm climate predominates (AW1) according to Köppen classification scheme  , and on the experimental period the mean temperature registered was 33˚C and the mean rainfall precipitation was 45 mm  .
2.2. Site Selection
The criteria used to select the crop planting site included the use of organic techniques by the farmer and a homogenous terrain with less than 5% of slope  . An agricultural plot with an area of 2400 m2 was chosen at a distance of 2 km from El Chichonal settlement. A year before our study, organic fertilization with vermicompost and tillage was applied to the soil for corn cultivation (Zea mays, L.), and the experimental site had been converted ten years before from conventional system to organic system production. Experimental units comprised plots of 25 m2. A two-meter-wide alley on each side separated each plot and ditches with a depth of 20 cm were dug around each plot to avoid the possibility of any effect caused by the contiguous treatment. All the experimental units were irrigated with plastic hoses and the plants were watered every 7 days. At the beginning the study, the physicochemical characteristics of the soil on the experimental site was determined and the initial values were: 11.05% of organic matter (OM), 0.55% of nitrogen (N), 7.47 pH, 19.80 mg∙kg−1 of phosphorus (P), 0.45 cmol∙kg−1 of potassium (K), 0.14 dS∙m−1 electrical conductivity (EC), 48.78 cmol∙kg−1 cation exchange capacity (CEC), and texture was 40.30% Sand, 31.70% Clay, 28.00% Silt, and 0.98 g∙mL−1 bulk density (BD).
2.3. Preparation of Organic Fertilizer
The organic fertilizer used was obtained from compost prepared using 40% corn stubble (Z. mays), 40% Taiwan grass (Pennisetum purpureum, Schumacher), and 20% Arnica leaves (Thitonia diversifolia, Hemsl). These materials were obtained on the farm and crushed with a Koltler hammer mill (2500 rpm and 14 HP) and mixed in a 1:1 proportion with dry sheep dung (Ovis aries, L.). The mixture was then composted for one month on a concrete bed and subsequently inoculated with 3 kg of red earthworms (Lumbricus rubellus, Hoffmeister), approximately 1500 red earthworms, in the substrate. The substrate was maintained around 70% to 80% of humidity with water irrigation, until vermicompost was formed through the transformation of the organic matter by the action of the microorganisms present in the digestive tract of worms. Simultaneously, during vermicompost formation, leachates were collected by runoff from the floor of the composter, in which a plastic tube was connected to a tank of 500 liters. Finally, samples of vermicompost and leachate were collected for chemical analysis (Table 1).
2.4. Experimental Design and Treatments
The tomato crop was sown in January 2015 and harvested in April 2015. A Westar seed of green tomato was used (variety “rendidora”). The seeds were sown in 5-cm-deep holes spaced 50 cm between them in each row. The distance between rows was 100 cm; ten seeds per hole were initially sown following the farmer’s agricultural technique. Subsequently, the best plants were selected, leaving 5 plants per hole when they had their first 5 true leaves in average and a mean height of 15 cm. Twenty experimental plots (5 × 5 m each one) were established. Each plot had 50 holes. A fully randomized block design was used, utilizing the slight (<5%) slope of the land as a blocking criterion  . The following treatments with five repetitions were tested: a) control, without fertilizer application (NF); b) application of 13,000 kg∙ha−1 of organic compost (vermicompost) equivalent to 160 kg∙N∙ha−1 (OTV); c) application of vermicompost and leachate
Table 1. Chemical characteristics of the organic fertilizer used on the tomato cultivation.
EC = Electrical conductivity (dS∙m−1), OM = Organic Matter (%), N = Total nitrogen (%), P = Available Phosphorus (mg∙kg−1), K = Potassium (cmol∙kg−1); -- Not available.
at 20 L∙ha−1 (OTH) and d) application of 600 kg∙ha−1 of chemical (inorganic) fertilizer containing N, equivalent to 160 kg∙N∙ha−1 (CST) (Figure 1).
Fertilizer application and crop maintenance were carried out following Gûemes-Guillen et al.,  recommendations, included in the green tomato growing guide prepared by National Research Institute of Forestry, Agriculture and Livestock (Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias―INIFAP). Vermicompost was applied manually in each plant at 20 days of growth, and leachate was prepared as liquid fertilizer by dilution of 1 liter of leachate in 20 liters of water, and then was sprayed with a backpack spray pump at the same time that vermicompost was applied. Pest control was achieved by applying Bacillum® at 125 mL∙ha−1 to the parcels under the organic treatments, and Falvus® was applied at 30 g∙ha−1 to the inorganic treatment. The control treatment (NF) had no experienced pest management to avoid a confused effect due to the application of agrochemicals.
2.5. Arbuscular Mycorrhizal (AM) Fungi Analysis
AM fungi identification and morphologic analysis was conducted on soil samples that were randomly collected from the experimental plots. The quantification of spores and species diversity were analyzed in soil samples that were collected at a depth of 15 cm. The fine roots were taken, at the end of the crop cycle, at 15 cm depth and the soil was washed in a sieve with running water to obtain only fine root samples. The fine root samples were dried at shade and placed in 100-mL jars containing 5% acetic acid and stored at 4˚C until analyzed in the laboratory.
For the laboratory analysis, spores present in 100-g soil samples were extracted and counted by sieving method, decanting wet and sucrose gradient centrifugation, as suggested by Moreira et al.,  . The spores were detached using a dissection needle under a Nikon SMZ800 stereoscopic microscope and were mounted in microscope slides with concentration polyvinyl alcohol solutions. Where necessary, the spores were dyed using Melzer’s reagent. Taxonomic determination was performed based on observations of the form, size and coloration of the spores under a transmitted light microscope (100×) (Nikon Eclipse E600) and were collated with the guidelines proposed by Moreira et al.,  and the International Culture Collection of Arbuscular Mycorrhizal Fungi (2014) (http://invam.caf.wvu.edu).
To assess the similarity in AM fungal communities amongst treatments, a Jaccard’s index was calculated and defined following formula:
Jaccard’s index = a/(a + b + c)
where a = number of AM fungal species common for two treatments, b = number of AM fungal species present in the first treatments and absent in the second, and c = number of AM fungal species absent in the first treatments and present in the second treatments   .
To estimate the percentage of AM fungi colonization fine roots, segments were examined under a transmitted light microscope (100×) (Nikon Eclipse E600). The presence or absence of structures such as hyphae, arbuscules and spores were determined in three equidistant visual fields for each segment  . Additionally, the length of hyphae (µm/mm root) was measured in each segment following the procedures described by Vega-Frutis and Guevara  .
2.6. Crop Yield
Yield was recorded after fruit maturation on the plant. Fruits from 10 plants corresponding to each treatment were counted and weighed using a Torrey MFQ-40 electronic scale. Representative samples of each treatment were obtained, and the dry matter content (DM) was determined using a forced-air drying oven. Simultaneously, the equatorial and polar diameters of the fruits were measured using Vernier calipers. Fruit size was estimated based on the following formula:
V = 4/3 (π*A*2B)
where A corresponds to the polar radius (height), B corresponds to the equatorial radius (width and depth), and π = 3.1416  .
2.7. Crop Survival, Growth and Production
Plant survival was determined based on the percentage of plants that died throughout the growing period. Plant height was measured with a tape measure from the base of each plant to the highest branch. Measurements were made upon application of the experimental treatments and after the growing period. Relative growth rate (RGR) was estimated using the following formula:
where A2 and A1 represent plant height at times 2 (end of the growing period) and 1 (when were applied the experimental treatments); and t2 and t1 represent the days at the end (t2) of the growing period and the moment of the treatment application (t1)   .
Above-ground and root biomass were also recorded throughout the total harvest of the plants. Above-ground biomass was quantified as stem, leaf, flowers and fruits production. Root biomass was obtained from the extraction of 20 cm long × 20 cm width × 15 cm depth soil cubes around each plant, helped with a shovel flat. The soil was removed manually with a brush and washed in a sieve with running water to obtain fresh roots. Fresh roots were dried at shade for 24 hours and weighted using a Torrey MFQ-40 digital scale. Later, to obtain dry matter content, roots were dried using a forced-air drying oven during 3 days at 50˚C, until constant weight were obtained.
2.8. Soil Analysis
During soil preparation to sown and at 50 days after starting the experiment, three soil samples per treatment were collected for physicochemical analysis. We used a representative random sampling with 95% confidence interval (CI) to collect the soil samples, it were taken at 20 cm from the tomato plant stem and a depth of 15 cm. The collected soil was sieved through a 5-mm sieve and then dried in the shade at room temperature. The samples were then sent to the laboratory to determine available phosphorus (P) by Olsen’s method, organic matter (OM) was determined by Walkley-Black method, soil pH was measured in water (2:1 ratio) soil:water suspension, for total nitrogen (N) determination samples was digested and quantified by titration according to method of micro-Kjeldahl, potassium (K) was determined by ammonium acetate (1 N, pH 7), electrical conductivity (EC) was measured in a 1:5 (soil:deionized water) soil solution with a conductivity meter, cation exchange capacity (CEC) was determined as the total exchangeable base and ammonium acetate (1 N, pH 7), soil texture (clay, sand and silt percentages) were determined using a 10-g subsample and the Bouyoucos hydrometer method, and bulk density (BD) was determined by test tube of the volume of the dry soil mass displaced, according to the Mexican Official Standard (Norma Oficial Mexicana) NOM-021-RECNAT-2000   .
2.9. Statistical Analysis
Data on soil physicochemical variables (OM, pH, P, N, K, EC, CEC and BD), were analyzed through one-way analysis of variance. The data of crop yield (fruit weight and size) and production (plant survival, relative growth rate, above-ground and root biomass) were analyzed through a two way analysis of variance in a fully randomized block design and the means of treatments were tested with Tukey’s range test  using the statistical software R 3.1.2  . Data on AM fungi colonization (hyphae, length hyphae, arbuscules, spores in root, species diversity, and spores in soil) were analyzed using the Kruskal-Wallis test, and treatment medians were tested using the Wilcoxon test  .
The spores data counts for each species of AM fungi and soil physicochemical variables (P, OM, pH, K, EC, BD, and silt), were correlated using a canonical correspondence analysis (CCA) performed with R software, using “vegan” package  . Soil physicochemical variables were tested for outliers (Cleveland dotplot), normality (histogram), homogeneity (conditional boxplot), relationships (pair plots using “psych” package  ) and collinearity (Variance inflation factor using “car” package  ) according to Zuur et al.  .
3.1. Soil Arbuscular Mycorrhizal (AM) Fungi Diversity and Root Colonization
Twenty-one species and 16 morphotypes of AM fungi were found. OTH treatment had the largest number of species and morphotypes (17 and 11, respectively), whereas CST registered fewer species and morphotypes (7 and 4, respectively) (Figure 2).
The genera Glomus and Acaulospora were best represented in all treatments. With inorganic fertilization a reduced diversity of AM fungi in the soil was found; whereas organic fertilizer increased the AM fungi diversity, mainly by the presence of Claroideoglomus sp., Septoglomus sp., Ambispora sp., Gigaspora sp., Archaeospora sp., and Racocetra sp. (Figure 3). Also, with organic fertilization the following OTV-associated species and morphotypes were found: Acaulospora sp. 1, Ambispora sp., Claroideoglomus sp., and Glomus sp. 6. The species and morphospecies associated with OTH treatment were: Claroideoglomus claroideum, Claroideoglomus luteum, Funneliformis coronatum, Glomus sp. 5, Racocetra sp., Acaulospora denticulata and Archaeospora sp. Only Acaulospora rehmii was uniquely found in the NF treatment, and no unique species was found with CST treatment (Figure 2). Six species common to all treatments included the following: Funneliformis geosporum, Glomus glomerulatum, Rhizophagus intraradices, Glomus macrocarpum, and two Glomus sp. (Figure 2).
The effect of organic agricultural system on the presence of mycorrhizal interaction with the fine roots of tomato plants differed significantly (p < 0.05) between treatments. Plants subjected to organic agricultural system exhibited a higher percentage of total colonization (p = 0.03), vesicles (p = 0.02), arbuscules (p = 0.03) (Figure 4(a)) and length of hyphae (µm/mm root) (p = 0.02) (Figure 4(b)); compared with unfertilized plants. And the percentage of colonization by vesicles and arbuscules was significantly higher (p < 0.05) in the OTV and OTH treatments compared with conventional agricultural system (CST).
No significant differences were observed between the treatments in respect to number of species (p = 0.50) and number of spores (in 100 g) present in the soil (p = 0.65) (Figure 4(b)). But, the Jaccard index of similarity (Table 2) indicated that organic treatments had more AM fungal species in common than no fertilization or inorganic fertilization (CST). Nevertheless, there were significant differences in the diversity of AM fungal species.
3.2. Physicochemical Soil Characteristics and AM Fungi Species Relation
The relationship between soil physicochemical variables and the axis that resulted from the canonical correspondence analysis (CCA) are shown in Table 3 and Figure 5. The first component CCA1 was mainly explained by K (−0.740)
Figure 2. Species and morphotypes of Arbuscular Mycorrhizal fungi present on the soil in organic and conventional agricultural systems of green tomato (P. ixocarpa) cultivation. Red: NF = no fertilization, Blue: OTH = organic treatment with vermicompost and leachate, Green: OTV = organic treatment with vermicompost, and Yellow: CST = inorganic fertilizer treatment.
Figure 3. Number of species for genera of arbuscular mycorrhizal fungi present in the soil with four treatments for green tomato (P. ixocarpa) cultivation. NF = no fertilization, OTH = organic treatment with vermicompost and leachate, OTV = organic treatment with vermicompost, CST = inorganic fertilizer treatment.
meaning that lowers concentration of K promoted greater AM fungi diversity. On the other hand, variables that explained the maximum variance in CCA2 component were EC (0.792), BD (0.547), P (0.542) and they were positively related (Table 3).
Figure 4. Effect of organic and conventional agricultural system on the AM fungi root colonization on green tomato (P. ixocarpa) plants, and their effect on the number of fungi species and spores present in the soil. NF = no fertilization, OTH = organic treatment with vermicompost and leachate, OTV = organic treatment with vermicompost and CST = inorganic fertilizer treatment. (a) Root colonization; white: total (%); light grey: vesicles (%); grey: arbuscules (%); and dark grey: spores (%). (b) White: species (Spp, number) present in soil; light grey: spores (S, number) present in soil; and grey: length (l) of hyphae (µm/mm) in the roots. Kruskal-Wallis test; abc different letters on the bar of the same color indicate significant (p < 0.05) differences; ?: median, ┬ :1st and 3rd quartiles; outlayers.
Table 2. Jaccard’s index values of AM fungi diversity associated of the soil on organic and conventional agricultural systems of green tomato plants.
NF = no fertilization, OTH = organic treatment with vermicompost and leachate, OTV = organic treatment with vermicompost and CST = inorganic fertilizer treatment.
Silt was significantly correlated with the presence of AM fungi in the soil. K and EC showed a negative relation whit AM fungi diversity, in contrast with the silt (Table 4).
CCA showed that several AM fungi species were related negatively with the P, silt, pH, OM and BD, with independence of the treatments used (Figure 5). The
Table 3. Canonical correspondence analysis of the soil physicochemical characteristic effects in AM fungi diversity considering all tomato cultivation treatments.
P = available phosphorus, OM = organic matter, K = potassium, EC = electrical conductivity and BD = bulk density.
Figure 5. Canonical correlation analysis (CCA) showed the relationship between AM fungi species and physicochemical soil characteristics, under organic and conventional agricultural system of green tomato cultivation. Blue: P (mg∙kg−1), K (cmol∙kg−1), OM (%), pH (H2O), EC (dS∙m−1), BD (g∙mL−1), silt (%). Gray: AM fungi species.
total variation value of the data was 3.86, the sum of all canonical eigenvalues was 2.09, all seven explanatory variables explained 54% of the total variation in the data and the first two (canonical) eigenvalues were CCA1: 0.50 and CCA2: 0.43. The first two canonical axes explained 24% of the variation that was related with the soil variables.
Table 4. Canonical correspondence coefficients of the first two ordination axes showing the main contribution of the physicochemical characteristics on the tomato cultivation.
P = available phosphorus, OM = organic Mater, K = potassium, EC = electrical conductivity and BD = bulk density. *pα < 0.05.
3.3. Development and Yield of Tomato Fruits
No significant differences were found between treatments with respect to total number of fruits (p = 0.47), number of mature fruits (p = 0.11), weight of mature fruits (p = 0.28), equatorial diameter (p = 0.33), polar diameter (p = 0.40) or fruit size (p = 0.39) (Table 5).
3.4. Mortality, Growth and Plant Biomass
Plant mortality was 18.8% to 22.6% and did not differ significantly between treatments (p = 0.95). Although final plant growth was not significantly influenced by treatment (p = 0.95), differences in the relative growth rate (RGR) (p = 0.03) were observed according to treatment. OTV presented a higher RGR (10.47 cm) with respect to OTH and NF and was similar to CST treatment (8.18 cm) (Table 6).
No significant differences in total above-ground and root biomass were found (P = 0.07 and P = 0.77, respectively) between treatments at the end of the crop cycle.
3.5. Soil Physicochemical Characteristics
The soil predominantly comprised clay and loam. The OM, pH, P, N, and K soil content, and the EC, CEC and BD were similar (P > 0.05) in all experimental plots (Table 7).
4.1. Soil Arbuscular Mycorrhizal (AM) Fungi Diversity and Root Colonization
AM fungi are present in the soil and colonize plant roots of cultivated plants. They increase the physicochemical quality of the soil, stabilize soil micro aggregates and improve the mobilization and transformation of nutrients  . In addition, AM fungi promote biota interactions in the rhizosphere  ; the type of
Table 5. Effect of organic and conventional agricultural system on mean yield and development of green tomato (P. ixocarpa) crop.
NF = no fertilization, OTH = organic treatment with vermicompost and leachate, OTV = organic treatment with vermicompost and CST = inorganic fertilizer treatment. F: degrees of freedom of the origins of variance. S.E.: Standard error of the mean differences, (pα < 0.05), 1Dry weight.
Table 6. Effect of organic and conventional agricultural system in plant mortality, plant height, relative growth rate and biomass production of green tomato (P. ixocarpa) crop.
NF = no fertilization, OTH = organic treatment with vermicompost and leachate, OTV = organic treatment with vermicompost, CST = inorganic fertilizer treatment. S.E.: Standard error of the mean differences (Pα < 0.05). RGR = Relative Growth Rate of 43 days. 1mean dry weight; abcValues in rows followed by the same letter are not significantly different (P < 0.05).
Table 7. The physicochemical characteristics of the soil on organic and conventional agricultural systems of green tomato plants.
NF = no fertilization, OTH = organic treatment with vermicompost and leachate, OTV = organic treatment with vermicompost, and CST = inorganic fertilization treatment. S.E.: Standard error of the mean differences. (Pα < 0.05). OM = organic matter, P = available phosphorus, N = total nitrogen, K = potassium, EC = electrical conductivity, CEC = cation exchange capacity and BD = bulk density.
organic fertilizer used, influences the diversity of AM taxa in soils under organic agriculture. The use of compost in organic fertilization promotes a greater AM fungi species    .
The favorable impact of organic fertilizer application on the soil’s microorganism community, especially on the AM fungi community represents a potential short-term additional benefit   . The number of species of AM fungi have been reported for other agroecosystems under organic farming and is similar to the findings in this study  . However, as already mentioned, the application of organic fertilizers favors a greater diversity and richness of AM fungi species in the soil, whereas the application of inorganic fertilizers negatively affects diversity      , as the conventional agriculture affect the viability and growth of some taxa in the rhizosphere     . These changes are also observed in short-cycle crops, as in this study with green tomato, where the application of inorganic chemicals (CST) reduced the number of species of AM fungi.
The greater AM fungi abundance of the genera Glomus (Glomus glomerulatum and Glomus macrocarpum) and Acaulospora (Acaulospora delicata) in the treatments is due to differences in fungal strategies regarding to the use of host carbon sources under different physicochemical soil conditions  . Consequently, the genus Acaulospora is commonly observed in agricultural fields; similarly, the genus Glomus has been recorded in disturbed conditions because it can colonize roots and produce fungal structures in less time when compared with other genera   . These genera have been found in various crops under both organic and inorganic systems  , and in this study was not the exception (Figure 3). However, the abundance of Funneliformis geosporum, which was associated with CST, might reflect a greater ability of this species to tolerate changes in the soil nutrients under conventional agriculture; furthermore, this species might have a cost-benefit advantage in favoring crops under conditions of low soil P concentrations   . However, when organic fertilization is applied (OTH, OTV), the organic matter and P increase together with AM fungi diversity,  leading to the presence of the genus Claroideoglomus spp., as well as the other species that were uniquely found in these treatments, such as Acaulospora mella, Archaeospora trappei, and Glomus multicaule and the genera Ambispora spp. and Racocetra spp. (Figure 2). This finding might indicate that these taxa are more sensitive to chemical management and respond favorably to the greater soil nutrient availability mostly P and MO, resulting from the use of organic agricultural system , such as vermicompost   ; that could help to explain the similarity in species composition on organic treatments with Jaccard index showed in AM fungal community. In addition, fertilization with vermicompost promotes greater colonization and the formation of fungal structures (vesicles and arbuscules) in plant roots (Figure 4(a))    , and stimulates the growth of fine roots  , thus allowing the plant to explore a greater area and to have greater interaction and contact surface with AM fungi   .
4.2. Physicochemical Soil Characteristics and AM Fungi Diversity Relation
The CCA showed a relation between physicochemical soil characteristics and AM fungi spores’ diversity. K, CE and silt were the environmental factors related to AM fungi community composition. Most species were not positively associated with elevated P, pH, OM and silt concentrations. However, soil results showed these variables had higher values in organic treatments which suggested that diversity of AM fungi increased with the application of vermicompost and leachate and may benefit AM fungi development on soil and crop yield  . Avio, et al.  showed indirectly evidenced of AM fungi community structure may be affected by the conventional agriculture especially with the use of chemical fertilizers altering the concentration of soil mineral nutrients and physical characteristics. The organic management can promote higher AMF diversity and abundance than conventional agriculture, because it improves the quality of the soil   . The CCA showed too that no particular AM fungi genera had an association with some soil variable but the rare species may be more sensitive to conventional management  .
4.3. Development and Yield of Tomato Fruits
In the tropics and subtropics the lack of differences on tomato yield between organic and conventional agriculture is highly influenced by the precipitation, irrigation management, soil texture, the characteristic of organic input and the elapsed time after conversion from conventional systems to organic systems   . The use of vermicompost favors the growth, quality and yield of leaves, fruits and grains in crops   because it provides organic matter and micro- and macronutrients, and promotes soil fungi and bacterial activity 
In countries or geographical regions like Calakmul, where predominates 1000 - 1500 mm rain precipitation it is hardly to find differences in yield between organic and conventional systems  . Also, it has been observed that the irrigation management, as we supplied to tomato crop in this study, promotes lower yields on organic conditions. This is due to lower N use efficiency to produce fruits resulting in similar yield on organic treatments as respect to conventional system  . Other key factor that could influence these results and the quality of the fruit was the soil texture. In our study the soil was mainly composed of clay and loam, this fine textured soil decrease crop yield on organic system and it is consistently reported on other countries  . In addition, the conversion time from conventional to organic system could influence the lack of differences in fruit production in our study. As previously stated, the experimental site was converted from conventional system to organic system ten years before the present. In this respect, it has been observed that at least five years are required to equate the production of organic systems with a conventional system, and yet, after ten years organic farming only showed a trend to increase productivity  . Also, AM fungi community composition in soils of organic system could influence the crop production. It has been reported that some species of the genus Glomus  associated with organic treatments affect fruit production and maturation, resulting in a trend of lower total fruit production, number of mature fruits and fruit weight. On the other hand, our results differ from the findings of Carrubba et al.  , they found a favorable effect of organic system compared to conventional systems regarding the performance of coriander (Coriandrum sativum, L.). Also, our data regarding fruit size and development differed from the findings reported by Singh et al.  , who obtained better size and quality of strawberries (Fragaria × ananassa, Weston) using vermicompost compared to inorganic fertilization. These differences might be mediated by the environmental variations and soil physicochemical characteristics where was carried out the comparing studies   . These studies were carried out on semi-arid environments and soils with coarse texture; in this conditions crops yield boost when received organic fertilization    . Although organic management did not result in greater fruit yield and development, it is important to note that the CST increases the dependency on external inputs, increase the cost of the cultivation and jeopardizes the development of sustainable agriculture  . In contrast, it is possible to maintain production similar to that achieved using inorganic fertilizers by using vermicompost. This organic fertilizer reduces external dependency, increases the integration of diverse agroecological components, and improve the health of the soil   ; promoting a sustainable agriculture development and maintaining soil fertility, one of the key ecosystem service of agroecological systems.
4.4. Mortality, Growth and Plant Biomass
Fertilization with vermicompost and leachate provides humic acids, which favor the absorption of micro- and macronutrients, control viral infections and dampen the effects of water stress, thereby reducing plant mortality   . In addition, the application of organic fertilizers (OTV, OTH) increases the rate of plant growth, possibly due to the action of soil microorganisms that facilitate nutrient availability and AM fungi symbiosis, which in turn favor the absorption of essential minerals for plant growth   . However, the use of organic fertilizers can result in low carbon assimilation in the plants due to some mycorrhizae of the genus Glomus; this results decrease in the reduced production of above-ground biomass compared with inorganic fertilizers  , as observed in this study. The findings described here are consistent with findings reported for other crops, where greater production of above-ground biomass was observed for coriander (Coriandrum sativum L.), wheat (Triticum spp.) and medicinal plants in response to the application of inorganic fertilizers compared with organic fertilizers    . Despite lower above-ground biomass production, organic management with vermicompost favors the development and amount of root biomass in plants, perhaps due to a greater availability of soil nutrients   mediated by AM fungi, through their mycelia network that improve P uptake on plant growth   .
4.5. Soil Physicochemical Characteristics
The management of arable land influences the mineralization of organic matter and the availability of macro and micronutrients   . When soil is managed organically, the chemical characteristics of the soil and its nutrient quality and quantity remain unchanged  . But under medium and long-term organic management has been found an improved in chemical, physical and biological soil quality, increasing soil organic carbon, total N content and enhancing the soil organic matter, water holding capacity, aggregate stability and field-saturated hydraulic conductivity   . However, management with agrochemicals reduces mineralization processes and lowers soil quality, thus affecting OM, N, P, and K contents   . This effect was not clearly observed in the present study and the lack of the expected effects was most likely masked by the organic management that was applied to soil in the year before this study. In this regard, studies of longer duration are warranted to determine the effect of organic fertilization regime on the edaphic nutrient dynamics  , physical quality of the soil  and the activity of microorganisms that affect nutrient availability  .
On the organic agricultural systems, the use of vermicompost and lechate promotes a greater soil arbuscular mycorrhizal diversity, which interacts with the plant in the rhizosphere and results in higher plant growth rates with similar crop development and production that are obtained with conventional agriculture. Therefore, the use of vermicompost on organic agricultural systems can positively impact the production of green tomato in buffer zones of protected natural areas and this has implications for the future development of sustainable agriculture involving agroecological production systems.
The authors would like to thank the El Chichonal community, Calakmul, Campeche and the Institute of Ecology (Instituto de Ecología A.C); and the College of the Southern Border (El Colegio de la Frontera Sur) for facilitating the research and for technical support. Also thank the National Council for Science and Technology (Consejo Nacional de Ciencia y Tecnología, CONACYT) for a scholarship granted (588020) to the first author. We greatly appreciate the aid of Guadalupe Pérez Escobar, Miguel A. López Anaya and Jesús Carmona de la Torre for the chemical analyses, and Dra. Nelly Tremblay for the statistical analysis.