Abstract

The genus Capsicum provides antioxidant compounds, such as phenolics and carotenoids, into the diet. In Mexico, there is a wide diversity of species and varieties of chilli peppers, a fruit which has local cultural and gastronomic importance. In the present study, the relationship of the carotenoid and phenolic profiles with the RAPD fingerprint of three different commercial cultivars of chilli peppers of seven regions of Mexico was investigated. Through RAPD, the species of chilli were differentiated by means of different primers (OPE-18, MFG-17, MFG-18, C51, and C52). The genetic distance found with OPE 18 was in the order of 2.6. The observed differences were maintained when the chromatographic profile of carotenoids, and the molecular markers were analyzed, which suggest a close relationship between carotenoids and the genetic profile. While the chromatographic profile of phenols and the molecular markers were unable to differentiate between genotypes of chilli peppers. In addition, by using infrared spectroscopy and statistical PCA, differences explained by geographic origin were found. Thus, this method could be an alternative for identification of chilli species with respect to their geographic origin.

1. Introduction

Chilli peppers (Capsicum annuum L.) are used in a multitude of food preparations and are marketed in different regions. The greatest genetic variety of C. annuum L. can be found in Mexico, where there are chilli peppers of many different shapes, sizes, and colours [1], and they are among the most frequently consumed products. The genus Capsicum is classified into the family of Solanaceae, which is constituted by 25 wild and 5 domesticated species (C. annuum L., C. frutescens L., C. chinense Jacq, C. baccatum Jacq, and C. pubescens L.), including more than 200 varieties [2, 3]. Of these five domesticated species, C. annuum L. is the most commonly cultured [4]. Chilli peppers have been object of study mainly due to containing capsaicin, which produces different pungency levels [5], as well as carotenoids and phenolic compounds, which are used as natural pigments and antioxidant agents [6]. Chilli peppers may have different content and profiles of such compounds, depending on the genotype and variety, maturity of the fruit, and the environmental conditions of the cultures [7].

In order to preserve, manage, and improve the different chilli species, the evaluation of the extent of genetic variation within species, as this has now become a fundamental tool in biology and agriculture. Among the main criteria used for that purpose are morphological and molecular markers [8]. The RAPD (random amplified polymorphic DNA) method is a widely used technique for molecular marking and is based on the amplification of genomic DNA fragments by using primers of arbitrary nucleotide sequences; which in turn detects polymorphisms that can be employed as genetic markers without previous genetic sequences. In addition to being simple, fast, and low-cost, this method does not require radioactive markers and consumes minimum amounts of DNA [911]. However, one of the inconveniences reported about the RAPD technique, since it was published for the first time, is the low reproducibility found mainly in low-intensity bands. The most important factor affecting reproducibility in RAPD analysis is the preparation of the DNA template. Thereby, the differences in the concentration of DNA template between samples are observed as the gain or loss of some bands.

Recently, to improve the results, phytochemical markers have been combined with molecular techniques. Silva et al. [12] correlated the flavonoid content of medicinal plants from Brazil with molecular markers obtained by RAPD, whereas Ercisli et al. (2007–2008) [13, 14] obtained rapid results for evaluation of genotypic diversity and distances by using methyl esters of fatty acids as markers. Previously, the content and composition of fatty acids were tested as taxonomic markers in Hippophae rhamnoides L., showing that these substances, such as linoleic acid, can be used as biochemical markers for this species [15]. On the other hand, Füller et al. [16] analyzed the variability of morphological characteristics as well as the content of phenolic compounds and essential oils in plants from southern Brazil. Moreover, Rotini et al. [17] described the state of preservation of a meadow by comparing the content of phenolic compounds and the RAPDs.

Another essential technique for determining the biological nature of agriculture products is Fourier transform-infrared spectroscopy (FTIR), which is a technique that has been widely employed for food characterization because of several important advantages, such as its nondestructive nature and its capacity to yield structural information that constitutes a molecular fingerprint of the sample. A combination of FTIR and chemometric techniques like principal component analysis (PCA) has been employed as a direct and rapid way to discriminate properties in foods, including geographic origin, adulteration, and quality control [18].

The aim of the present study is to contribute to management and improvement of the different chilli species by determining the relationship between molecular and biochemical markers, such as the content of carotenoids and phenolic compounds among different commercial cultivars of chilli peppers from diverse geographic origins.

2. Materials and Methods

2.1. Biological Material

The samples were collected from central and north regions of Mexico (Table 1). Fresh and mature chilli peppers were used for this study. The chilli peppers were considered fresh when the samples presented bright green colour, and smooth, firm texture. The mature specimens were derived from the green samples, by exposure to sunlight and moisture-free storage in the laboratory for 25 to 35 days. The samples were considered mature when the chilli peppers had a bright red, yellow, or orange colour and less firmness and texture. The samples were washed with water and stored at −70°C.

2.2. Determination of Fruit Size

Length, width, and weight of the fruits were determined in accordance to the Mexican Official Standard [19]. Briefly, the length was taken from the base to the apex of the fruit without considering the peduncle; the width was measured at the widest part of the chilli pepper; for the weight, an analytical balance was used.

2.3. Determination of Carotenoids

The extraction of carotenoids from the chilli samples was carried out with 2 g portions of pulp and peel, which were mixed with diatomaceous earth and acetone at 10°C. Afterwards, the extract was filtered at vacuum through Whatman grade 2-filter paper and transferred to 10 mL of petroleum ether. The oil phase was extensively washed with distilled water to remove the residue of acetone, while the residual water was removed by shaking with anhydrous sodium sulphate. A 1 mL aliquot was used to determine the content of carotenoids using the molar extinction coefficient for β-carotene in hexane at 450 nm. The extracts were stored in darkness at −20°C [20].

2.4. Chromatographic Profile of Carotenoids

The extracts of carotenoids were reconstituted in 500 μL of HPLC-grade ethyl acetate and filtered through nylon membranes (0.45 μm, pore size) and then introduced into the C18 column (25 cm × 4.6 mm; 5 μm) of the HPLC equipment (Perkin Elmer Binary LC pump 250; series 200, UV/vis Detector). An isocratic elution system was established with acetonitrile-methanol-ethyl acetate (73 : 20 : 7) working at 0.6 mL/min for 70 min, and the absorbance was measured at 450 nm [21].

2.5. Determination of Free Phenolics

A total of 2 g of pulp and peel chilli peppers were mixed with 8 mL of 80% ethanol, and shaken for 10 min at 200 rpm. The mixture was centrifuged at 13000 rpm at 4°C for 10 min. The supernatant was separated by decantation and placed into an amber glass bottle. The sediment was reconstituted into 8 mL of 80% ethanol, and the procedure was repeated. Thereby, three extractions were obtained, joined, and stored under a nitrogen atmosphere until analysis [22].

Volumes of 0.75 mL of Folin-Ciocalteau reagent (1 : 10) were added to 100 μL of the extract of free phenolics, and the mixture was left to stand for 5 min in the dark. Afterwards, 0.75 mL of NaHCO3 60 g/L solution was added to neutralize the reaction. The solution was left to stand for 90 min, and the absorbance was determined at 725 nm. The result was considered as ferulic acid equivalents [23].

2.6. Chromatographic Profile of Phenolics

Aliquots of 50 μL of extract were injected in the HPLC equipment (Perkin Elmer Binary LC pump 250; series 200, UV/vis Detector) using a C18 column (25 cm × 4.6 mm; 5 μm). The mobile phase used was a binary system composed of 0.1% trifluoroacetic acid (Phase A) and 100% acetonitrile (Phase B). The absorbance was determined at 280 nm. A linear gradient elution was performed at 1 mL/min. Starting with 90% phase A and 10% phase B, the mixture was eluted to 55% phase A, 45% phase B in 50 min and to 60% phase A and 40% phase B in the last 10 min. Each compound was quantified by comparison with a calibration curve of the corresponding standard [24].

2.7. Statistical Analysis

Analyses of variance (ANOVA) were performed along with Student-Newman-Keuls test for differences between means, using SigmaStat 3.5 software. Principal component analysis (PCA) of the FTIR extract spectra was performed by using the OPUS QUANT (version 6.5) and Minitab (version 14) software.

2.8. Analysis by Infrared Spectroscopy

The dry extracts were reconstituted in petroleum ether, for carotenoids, and in ethanol, for phenolics, in order to obtain a semisolid consistency of the sample. Afterwards, the mixtures were analyzed in a FT (Fourier transform)-infrared spectrometer (Bruker Vertex 70), by using the ATR (attenuated total reflectance) sampling method previously described [25]. Briefly, a small sample amount is placed over a ZnSe crystal, where the infrared radiation is propagated and interacts with the sample to obtain the corresponding spectrum, which is averaged from several data acquisitions.

2.9. DNA Extraction

The method described by Allers and Lichten [26] was used for DNA extraction, after modifications. All samples of chilli peppers were washed with water, disinfected with sodium hypochlorite at 10%, 30%, and 70% ethanol, exhaustively rinsed with distilled water for 20 min and stored at −70°C, under sterile conditions. The pulp and peel of chilli peppers were crushed in mortars using liquid nitrogen. Portions of 100 mg per sample were mixed with 400 μL of lysis buffer (200 mM Tris-HCl pH 6.5, 70 mM EDTA pH 8, 2 mM NaCl, 2% PVP 40, 20 mM sodium metabisulphite, 1% Triton X-100, 24 mM MgCl2/MgSO4, 0.1% spermine, and 0.1% spermidine) previously heated at 60°C. The mixture was agitated and incubated at 60°C for 60 min; then it was transferred to ice, at which time 400 μL of 10% PVP was added. The mixture was agitated by inversion of the recipient and left to stand for 60 min at −20°C. Afterwards, the blend was centrifuged at 4000 rpm for 30 min, and 400 μL of each supernatant was transferred into new tubes, in order to add 5 μL of 10 mg/mL RNasa for 10 min at room temperature. Immediately, the supernatants were precipitated by addition of 200 μL of 10 M ammonium acetate and 600 μL of cold isopropanol. The mixture was agitated by inversion of the tube and incubated overnight at −20°C. The mixture was then centrifuged at 4000 rpm for 30 min at 10°C, the isopropanol eliminated, and the pellet washed with 1 mL of cold 70% ethanol, followed by centrifugation for 10 min under the same operation conditions. Finally, the ethanol was removed, and the pellet was dried for 20 min at room temperature.

The pellets were reconstituted into 200 μL of 10 mM Tris pH 8, and the DNA was precipitated again with 200 μL of PEG 8000 30%-NaCl 1.2 M at 4°C for 60 min. The mixture was then centrifuged for 30 min at 4000 rpm and 10°C. The supernatants from PEG were discarded, and the DNA pellet was washed with 1.5 mL of cold 70% ethanol, followed by centrifugation for 15 min at 4000 rpm. The ethanol was eliminated and the DNA pellet left to dry, and finally reconstituted into 200 μL of 10 mM Tris pH 8 [26].

The purity of the extracted DNA was determined by absorbance at 260/280 ratio, and values greater than 1.8 were accepted for future procedures. The DNA was quantified by its absorbance at 260 nm [27].

2.10. RAPD

For DNA amplification, a thermocycler (Bioer XP cycler) was used under the following operation settings: 1 cycle of 5 min at 94°C; 40 cycles of 1 min at 94°C; 1 min at 36°C, 2 min at 72°C; a final extension step of 10 min at 72°C. PCR reactions were performed in a volume of 25 μL containing 1X PCR buffer, 50 nM of MgCl2, 200 μM dNTP, 1.5 units of the enzyme Taq DNA polymerase (Invitrogen, Brazil), and 40–100 ng/μL DNA samples [27]. Four independent reactions were assayed using 0.4 μM oligonucleotide OPE18 (CGGCCCACGT), or 0.2 μM oligonucleotide MFG17 (CGCGTTCTTG), 0.2 μM oligonucleotide MFG18 (CGGCCCACGT), or 0.2 μM C51 (ATCAACGTACGT) and 0.2 μM C52 (GTCGACGGACGT) oligonucleotide mixture (Invitrogen) [4, 28].

The products of DNA amplification were analyzed by electrophoresis in 1.5% agarose gel, using TAE 1x buffer. The gels were stained with ethidium bromide (0.2 μL/mL) and observed with an UV-transilluminator (Stratagene Eagle Eye).

In order to increase reproducibility, optimum conditions were established for DNA extraction and amplification, without modifying Mg2+ concentration, for Taq polymerase enzyme and primers as well as for the DNA concentration in each reaction.

2.11. Data Matrix and Dendrogram

The polymorphic bands and signals of the chromatograms were considered as “present” or “absent” in a matrix of similarity. For the chromatogram, only signals above 40 milli-units of absorbance were considered “present.” The method of average grouping UPGMA was used and determined by means of the analysis Cluster by using the software Paleontological Statistics (PAST) 2.10, taking into account the Dice’s similarity coefficient and the Euclidean distance.

3. Results and Discussion

3.1. Size and Weight

The weight of the different varieties of chilli peppers ranged from 3 to 46 g. In general, samples with the best morphological characteristics (color, size, texture, and weight) were those from the northern region, except for the sample JV1 from the south-central area, the latter being the lighter weight. The length was from 4.3 to 9.5 cm, and the width was in the range of 0.1 to 3.34 cm. The three varieties of chilli showed differences in both weight and size, with remarkable greater similarities between the same species (Jalapeño and Serrano) as well as those grown in nearby areas. For instance, the Serrano chilli peppers grown in Tx, P, and V1 from the central region were quite similar. Through statistical analysis, significant differences were found in weight and width between varieties, but not in length between Serrano and Jalapeño (i.e., the same species). However, these two peppers showed significant difference with respect to Habanero (i.e., a different species) (see Table 2).

3.2. Total Content of Carotenoids

The content of carotenoids in fresh peppers varied within the interval 0.083–0.99 mg/g dry basis (DB), while mature peppers were in the range from 0.69 to 13.85 mg/g DB. The carotenoid content increased from the fresh to the mature state, as expected accordingly to the natural biosynthesis of pigments as the fruit matures [29]. In accordance with the statistical analysis, of chilli peppers in the fresh stage, there were significant differences between Jalapeño and Serrano varieties, while analysis of mature peppers, showed no differences between the three cultivars (see Table 3). This is due to the ability of fruits to synthesize carotenoids during the ripening process regardless of the species [30]. Regarding the region, fresh Jalapeño peppers from Tx and P were found to be different than those from V and S, while fresh Serrano peppers showed statistical difference in all cases except those from Jc and those from Si. Mature chilli peppers, from Si showed differences with respect to those from T, P, V1 and V2. Therefore, no association was found between the concentration of carotenoids and the geographical origin.

The results presented above could be affected also by the environmental conditions, such as temperature and relative humidity, as well as the type of soil, which of course are related to the proximity of the different geographical regions.

3.3. Identification and Quantification of Carotenoids

Among the identified carotenoids were lutein, β-cryptoxanthin, lycopene, and β-carotene, with a higher concentration of lutein and lycopene. The latter are part of the main carotenoids present in chilli [29, 31]. During the process of maturation, Lutein remained constant for Jalapeño peppers, whereas it diminished for Serrano ones. However, it has been reported that lutein disappears in mature peppers by the effect of the synthesis of pigments in chromoplasts [32]. Similarly, many of the hydroxylated carotenoids are progressively esterified with fatty acids [33]. Lycopene also tends to decrease as maturation proceeds. Contrarily β-cryptoxanthin and the β-carotene increase during the ripening process, as reported by Marín et al. [29].

According to statistical analysis, no significant differences were found for fresh peppers in concentrations of lutein, lycopene, and β-cryptoxanthin between Jalapeño and Serrano cultivars. Nonetheless, there were differences in β-carotene content, which is the major compound in several varieties of C. annuum L. [30]. On the other hand, differences were found for mature peppers only between Jalapeño and Serrano, in this case in lutein and lycopene content (Table 4).

3.4. Total Content of Phenolics

The concentration of phenolic compounds in fresh peppers ranged from 30.2 to 128.4 mg/g, whereas in the mature peppers it ranged from 38.6 to 297.8 μg/μL on dry weight basis. In most of the samples, there was a noteworthy increase in these compounds in mature chilli peppers, in agreement with the results of Deepa et al. [32]. However, other studies, such as those by Oboh et al. [5] and Zhang and Hamauzu [31] show that the content of phenolic compounds are significantly higher in the fresh peppers than in the mature ones. Given the above, the main factors involved in the variation of phenolic content are the maturation state, as well as the age of the plant, as reported by Deepa et al. [32].

When doing a comparison on the total phenolic content within the chilli varieties, significant differences were found between Serrano and Jalapeño, for both fresh and mature stages, while between Jalapeño and Serrano as a group no differences were detected with respect to Habanero. These results are inconsistent with data reported by Oboh et al. [5], where the phenolics content was significantly higher for the C. annuum L. than for the C. chinense Jacq species.

Concerning the region, differences were found in most of fresh Jalapeño peppers samples, except between P and V (V1 and V2) and between Si and Jco. For the fresh Serrano peppers, there were differences between Tx and P with respect to Jc and V2. For mature peppers, the region seems to be unimportant (see Table 5). From these results, certain patterns indicate differences between varieties and other between cultivars, in agreement with Antonious et al. [34] who reported significant differences between cultivars of C. baccatum L. and 𝐶 . chinense Jacq, and between two genotypes of the latter variety.

3.5. Identification and Quantification of Phenolic Compounds by HPLC

The phenolic compounds found in chilli peppers were gallic, protocatechuic, ferulic, o-coumaric, p-coumaric, sinapinic, trans-cinnamic and caffeic acids, quercetin, catechin, rutin, and vanillin. Many of these compounds belong to the pathway of phenylpropanoids which is characteristic of chilli species.

Among the varieties of chilli, significant differences were found between Serrano and Jalapeño in the gallic and p-coumaric acids content, indifferently of maturation stage. Regarding protocatechuic and caffeic acids as well as catechin and vanillin, there were differences only for fresh peppers. For the mature stage, the only difference was found in the content of trans-cinnamic acid between the two varieties (Table 6).

Regarding the place of origin, fresh Jalapeño peppers showed differences in the content of sinapic acid between Si samples and the rest, in the content of catechin of V samples compared to P and Tx samples (despite the fact that all three belong to the central region), and in the content of p-coumaric acid of P compared to V1, Tx, and Si. The fresh Serrano peppers showed differences in the concentration of gallic acid for samples from Si with respect to P and V. In the content of protocatechuic acid for samples from P with respect to the other specimens, and in the content of caffeic acid and rutin for V2 with respect to the other specimens. Finally, most regions showed differences in the content of p-coumaric acid, except for V2 compared to P and S, V1 compared to Jc and Tx, P compared to Si, and Jc compared to Tx, showing no association between geographical areas.

Among mature Jalapeño peppers, T showed differences with respect to Si in the content of gallic acid and with respect to V1 in protocatechuic acid. No differences were found between the mature Serrano peppers with respect to region.

By taking into account the state of maturation, for the rest of the phenolic compounds, there were not many differences as observed by the maturation stage in the content of gallic, protocatechuic, trans-cinnamic and p-coumaric acids; quercetin, rutin, and vanillin. For the rest of the phenolic compounds, there were no differences that could be explained by the maturation stage.

Within a Jalapeño variety, differences were found between fresh and mature peppers in the content of gallic acid, while no such difference existed in the content of ferulic acid. For mature specimens no differences were found only between the three varieties in the content of protocatechuic acid, caffeic acid, and catechin, while the fresh samples reported differences in the content of p-coumaric acid, catechin, and vanillin.

It is clear from the results that variations exist between the three variables and concentrations. Some interactions can be explained by the proximity of the geographical areas, as in the case of P, Tx, V, and T. Other differences may be due to the growing conditions of each geographic region or by the genotype of the chilli, as discussed earlier by Silva et al. [12]. These authors studied the content of rutin in a medicinal fruit from Brazil and found that the high variability in the concentration and the great genetic variability of the fruit may be associated not only with genetic differences, but also with environmental conditions.

3.6. Analysis of Phenolic Compounds by Infrared Spectroscopy

The typical IR spectra (600–2000 cm−1) of phenolic extracts from different geographical regions of Mexico (Figure 1) do not show a clear difference between samples. However, by applying principal component analysis (PCA) (Figure 2) to the second derivative of the spectra, four clusters were found. In the first cluster samples, P appeared with a variation axis which is parallel to the association to a second group of Tx and V, that is, with similar chemical characteristics. The third cluster, Si and T samples were able to associate and the fourth group to Jc samples, but perpendicular to the previous samples, so, the latter three were related to each other, but showed differences with the former.

We can say that the associations of geographic location could be due to the environmental conditions, that is the geographical environment determines the concentration of compounds in peppers. By joining these results with the total content of phenolics, and by chromatography, the same association for P, Tx, and V in concentration of trans-cinnamic, gallic, p-coumaric, protocatechuic acids, quercetin, and catechin can be found.

In the present study, it has been found that FTIR spectroscopy combined with the PCA method presents great potential for the verification of the geographic origin of phenolic samples. PCA gives a dimensional matrix. Regarding the principal components that were generated, that was too reduced to perform discrimination. Such discrimination was realized by using the second derivative of the spectra in the region (900–1750 cm−1). These results show that this FTIR-PCA analysis could be used as an alternative, quick, and low-cost method for identification of the geographic origin of Mexican chilli peppers (C. annuum L.) through their phenolic compounds.

For variables regarding type and maturation stage, it was not possible to find differences because the geographical variation, influenced by climatic and growing conditions, was more important than these variables.

3.7. RAPD

Among the oligonucleotides tested, OPE-18 was the most polymorphic RAPD oligonucleotide generating a total 10 bands. From the RAPD reactions carried out with a mixture of MFG-17 and MFG18 oligonucleotides and C51 to C52 oligonucleotides mixture a total of 25 and 17 bands were obtained, respectively. Independently of the oligonucleotide used, Serrano and Jalapeño chilli peppers exhibit similar RAPD profiles, whereas Habanero chilli peppers exhibited a particular profile with only 3 or 4 common bands. The majority of the bands were polymorphic, although some monomorphic bands were observed. Some examples of profiles obtained can be observed in Figure 3.

The presence/absence and similitude matrixes derived from the most informative primer (OPE-18) were obtained, and an Euclidean distance dendogram was constructed with 10 polymorphic bands ( 𝑟 = 0 . 9 4 ) (Figure 4). C. chinense Jacq. (Habanero chilli peppers) samples were clearly clustered apart from the C. annuum L. (Serrano and Jalapeño chilli peppers), confirming the clear differences (coefficient of 2.0) observed among RAPD profiles. Basically, all samples formed three clusters: the first group made by the habanero peppers and the second and third groups do not show a clear separation between cultivars. A previous RAPD analysis using 10 germplasm samples of C. annuum L from Thailand revealed a higher diversity of RAPD profiles those found in Mexican samples [4]. Moreover, other RAPD diversity analysis of C. annuum variety “Cuneo” from northwest Italy distinguished five populations [35]. Possibly, the limited profile diversity and high similitude observed in two Mexican varieties studied can be explained by selective pressure derived from culture practices and commercial interests. Evidently, RAPD confirmed the great morphological diversity observed among C. annuum L varieties around the world, but an extensive study is necessary for molecular intraspecific diversity and population genetics including many C. annuum L varieties collected from different parts of the world.

3.8. Genetic Analysis with Chromatographic Profiles

According to the dendrogram coupled to the amplification by oligonucleotide OPE-18 and the chromatographic profile of carotenoids in mature peppers ( 𝑟 = 0 . 6 9 ) (Figure 5), the separation between the Habanero peppers with respect to Jalapeño’s and Serrano’s can be observed, although it was not possible to identify between varieties within the same species or to identify an association between regions. Amorim et al. [36] found that the carotenoid content and the use of molecular markers in different genotypes of bananas, provided useful information on kin selection by crossing among different genotypes in order to develop new varieties with functional properties.

Our results could be influenced only by the RAPD because the carotenoid dendrogram did not show differences between species; however, this separation maintained between the species does not occur in the case of phenols. Therefore, the carotenoids are closely associated with molecular markers, specifically OPE 18.

4. Conclusions

The present study was able to establish a separation and differentiation between the chilli species C. annuum L. and C. chinense Jacq, through the RAPD molecular method by using different primers. Also, by chromatographic profiles of carotenoids and molecular markers in mature chilli peppers and amplification by oligonucleotide OPE 18, it was possible to distinguish between Habanero and Serrano/Jalapeño varieties, but no differences were found between C. annuum varieties. On the other hand, the PCA of phenolic compounds in peppers, showed four groups according to their place of origin. Hence this technique can be used to identify chilli varieties from different geographical areas, which might be used for the protection of rights in the variety, diversity and phylogenetic analysis, as well as for the confirmation of hybrids.

Acknowledgments

This paper was supported by the project SIP-IPN 20110329. The authors thank EDI, COFAA, and SNI. I. Troconis-Torres received a study Grant 215834 from the National Council of Science and Technology-Mexico (CONACYT).