Madiha Rashid1,2,3*†
Zubaida Yousaf1*†
Ahmad Din4
Muneeb Munawar5
Arusa Aftab1
Nadia Riaz1
Afifa Younas1Ibrahim A. Alaraidh6Mohammad K. Okla6Hamada AbdElgawad7- 1Department of Botany, Lahore College for Women University, Lahore, Pakistan
- 2Department of Botany, Division of Science and Technology, University of Education, Lahore, Pakistan
- 3Institute of Biological Sciences, University of Malaya, Kuala Lumpur, Malaysia
- 4National Institute of Food Science and Technology, Faculty of Food, Nutrition and Home Sciences, University of Agriculture, Faisalabad, Pakistan
- 5Vegetable Research Institute, Ayub Agricultural Research Institute, Faisalabad, Pakistan
- 6Department of Botany and Microbiology, College of Science, King Saud University, Riyadh, Saudi Arabia
- 7Integrated Molecular Plant Physiology Research Group, Department of Biology, University of Antwerp, Antwerp, Belgium
Leafy vegetable crops are considered as a natural source of mineral nutrients that could decrease the risk factor of many growth issues in children and adults. Spinach is globally considered as the most desirable leafy crop, due to its taste and nutrient richness along with greater nitrate contents and better nitrogen use efficiency. To evaluate the mineral nutrient efficiency of this crop, thirty genetically diverse spinach accessions were analyzed through nutritional and functional marker strategies. The accession 163,310 from Pakistan was found to be rich in minerals (sodium, calcium, potassium, zinc, and manganese) and nitrates. However, the oxalate contents were lesser in the accessions that had greater quantity of nutrients. These represented a negative correlation between mineral availability and oxalate accumulation in the leaves. To study the relationship of oxalates and minerals in the accessions, a functional marker analysis was performed, based on the genes involved in oxalate metabolism and disease resistance in spinach. High level of genetic polymorphism was observed among the accessions represented with 115 polymorphic bands out of 130 bands. Heat map clustering represented the accessions from Asian countries (Pakistan, India, China, and Iran) as the most adaptable accessions to the local environment. The correlation between nutritional and genetic analysis also revealed the nutrient richness of these accessions along with good oxalate metabolism and disease resistance. Hence, these accessions could be considered as useful genotypes in future breeding programs.
Introduction
Mineral nutrients play a vital role in proper growth and functioning of plants. A number of nutrients, including iron, magnesium, manganese, zinc, etc., act as cofactors of various enzymes, required for normal functioning of plant (Karthika et al., 2018). The deficiency of essential nutrients and vitamins, around the world, is leading toward high rate of physiological disorders, including slow growth and bones’ problem (Cousins, 2006; Abu Al-Qumboz and Abu-Naser, 2019). Therefore, a nutrient-rich food is needed as a supplement for the proper growth of both the children and adults. Leafy vegetables are a good source of nutrients that can potentially decrease the risk factor of these disorders (Odhav et al., 2007). Spinach (Spinacia oleracea L.) is a nutritionally rich leafy vegetable belonging to family Amaranthaceae. It is native to central Asia and thought to be originated from Iran (Sabaghnia et al., 2014). It is being cultivated in more than 60 countries. Moreover, green leafy vegetables are a natural source to recover the prevailing nutritional deficiency in developing countries (Arif et al., 2017). Common use of these vegetables in daily diet has been suggested as a cheapest source of micronutrients (Egbi et al., 2018). The WHO has recommended the consumption of at least 400 g of fruits and non-starchy vegetables, in order to meet the daily requirement of nutrients in human body (Arif et al., 2017).
Spinach is considered as the most desirable leafy vegetable in many countries around the globe, due to its richness in nutrients and a number of health-promoting compounds (Shi et al., 2016; Tang et al., 2019). Spinach leaves are rich with a number of vitamins, minerals, antioxidants, carotenoids, steroids, apocyanin, flavonoids, and p-coumaric acid (Avsar, 2011; Rashid et al., 2014). Its nitrates are helpful in maintaining blood pressure. Though its greater concentrations could also be toxic to human body, when nitrates are converted into nitrites (Mehri et al., 2021). It is also rich in vitamin K1, essential for blood clotting (Avsar, 2011; Rashid et al., 2014; Gunnars, 2019). Verma (2018) reported its potential uses in the remedy of urinary calculi and bowel and lung inflammation. It strongly prevents cancer cell proliferation, especially in breast cancers because of its high lutein and carotenoid contents (Avsar, 2011). The pharmacological properties of this leafy vegetable are due to its richness in a variety of phytochemicals, including its specific nutrients.
Crop variety and nutritional content of vegetable crops are particularly important for enhancing food and nutrition security on a global scale (Fresco and Baudoin, 2004). Moreover, it could be an important information for plant breeders to select improved varieties (Grant et al., 2008). However, the concentration of these nutritional components could be directly influenced by a wide range of climatic conditions (Natesh et al., 2017). The advent of molecular marker technology has transformed the study on crop’s quality (Lau et al., 2015). DNA-based markers have now been considered as a modern and more accurate tool to assess the genetic diversity irrespective to abrupt climatic influences. Among them, the functional markers (FMs) are directly linked to the allele of the trait of interest, gradually reducing the need of phenotypic validation (Varshney et al., 2005). Leafy vegetables, especially spinach, are rich in oxalates and nitrates. These important constituents have been found linked with the availability of essential nutrients in the crop (Joshi et al., 2021). Therefore, this study was conducted to assess the availability of mineral nutrients in genetically diverse spinach accessions using biochemical and functional marker strategies.
Methodology
Seeds of 30 spinach accessions (Table 1) were obtained from the gene banks of the National Agricultural Research Center (NARC), Islamabad, Pakistan and the United States Department of Agriculture (USDA) and grown in the field area of Vegetable Research Institute, Ayub Agricultural Research Institute, Faisalabad, Pakistan. After the growth of 2 months, the leaves were collected, rinsed, dried, and grounded into fine powder and subjected to nutritional analysis from Food Technology Section, Post-Harvest Research Center, Ayub Agriculture Research Institute, Faisalabad, Pakistan.
Table 1. Geographical location of thirty spinach accessions for nutritional and genetic analyses.
Proximate Analysis
Proximate analysis was performed to determine the moisture, ash, protein, fat, fiber, and carbohydrate and calorie contents of spinach accessions. Proximate contents were determined by following the methods of the Association of Official Analytical Chemists AOAC (2000).
Moisture content: Plant samples were weighed and dried in an oven at 105°C for 3 h. After drying, the sample was reweighed. The weights of sample before and after drying were used to calculate the moisture content of each sample (AOAC, 2000).
W1 = Weight of sample before drying.
W2 = Weight of sample after drying.
Ash content: About 5 g of plant sample was heated over low flame, till the fumes stop coming out. Then, it was placed in furnace overnight at 550°C and weighed once the ash turns to gray (AOAC, 2000).
Protein content: It was determined by Kjeldahl’s method, including three main steps, i.e., digestion, distillation, and titration. About 0.5–1.0 g sample was put in a digestion flask. 5 g Kjeldahl’s catalyst and 200 ml concentrated H2SO4 were added in it. Then, another tube was prepared containing all the above chemicals, except sample. The flasks were placed in inclined position and heated gently until solution clears. After cooling, 60 ml of distilled water was added in it and the flask was immediately connected to digestion bulb on condenser with tip of condenser immersed in standard acid and 5–7 drops of mix indicator in receiver. The contents were mixed thoroughly and heated till all the ammonia was distilled. The receiver solution was removed and excess standard acid distilled was titrated with NaOH solution (AOAC, 2000).
A = Volume of 0.2 N HCl used sample titration (ml).
B = Volume of 0.2 N HCl used blank titration.
N = Normality of HCl.
W is the weight of sample (g).
1.4007 is the atomic weight of nitrogen.
6.25 is the protein-nitrogen conversion factor.
Fat content: About 3–5 g sample was weighed and wrapped in filter paper. The sample was extracted with petroleum ether into soxhlet extractor. After drying of solvent, the sample was reweighed to determine fat contents (AOAC, 2000).
Fiber content: A total of 5 g of plant material was mixed with 200 ml of 1.25% sulfuric acid, heated for 30 min and filtered. The residue was washed with distilled water till it becomes acid free. Then, boiled in 200 ml of 1.25% of sodium hydroxide for 30 min, filtered, and washed with distilled water to make it alkaline free. The residue was rinsed once with 10% hydrochloric acid, twice with ethanol, and thrice with petroleum ether. Then, the residue was dried in oven at 105°C. After cooling, it was ignited in a muffle furnace for 90 min at 550°C to get the weight of fiber (AOAC, 2000).
Carbohydrate content: (Pearson, 1976; James, 1995; Otitoju, 2009).
Energy content: Energy contents were determined using bomb calorimeter. Small pellets of sample were produced on a briquette press followed by bombing with 25 atmospheric pressure (Anon, 1970).
Mineral Analysis
Mineral contents were determined by wet acid digestion method followed by atomic absorption spectroscopy (AOAC, 1990). About 2 g of powdered plant material was taken in crucible and ignited for 6 h in a muffle furnace at 550°C. The ignited residue was dissolved in 10% nitric acid (10 ml) and heated slowly for 20 min. Then, the filtrate of this material was used for the determination of mineral contents by atomic absorption spectrophotometry. The technique was used for the quantification of both the macronutrients (sodium, calcium, and potassium) and micronutrients (zinc, iron, and manganese) of spinach accessions.
Oxalate Content Determination
A total of 1 g dried plant sample was boiled for 25 min in 150 ml of water containing 27.5 ml, 6 M HCl along with 2 drops of octanol. After cooling, the mixture was diluted up to 250 ml and filtered. The filtrate of 10 ml was evaporated at 40–45°C and redissolved in 0.01 M 10 ml sulfuric acid. A solution of 0.5 mM K2Cr2O7 and 0.25 mM MnCl2 was prepared as blank solution. A total of 10 ml solution (1 ml extracted sample + 9 ml blank solution) was taken for each sample, incubated for 60 min, and then analyzed using UV-Vis spectrophotometer for the determination of oxalate contents in the sample (Popoola et al., 2014).
Nitrate Content Determination
A total of 20 g of each sample was homogenized with 100 ml of deionized water and centrifuged at 3,000 rpm for 10 min. Supernatant was collected for spectrophotometric analysis. The sample was diluted to 1/10th and 1 ml solution was added to an ion-exchange solid phase extraction (SPE) cartridge, which was preconditioned with 2 ml of deionized water. The solution was then washed out with 10 ml of deionized water, followed by 0.5 M NaCl solution to elute the nitrate components. 5 ml elute was collected for spectrometry analysis. Nitrate in aqueous solution results in an absorbance peak at wavelength of 220 nm and no absorbance at 275 nm. The net absorbance (Abs) for nitrate was calculated using the following equation (Yu et al., 2018) and converted into the nitrate concentration via a calibration curve, set by a series of standard solutions with concentrations from 0 to 25 μg/ml. The standard solutions, including a blank, were also prepared in 0.5 M sodium chloride solution.
The derived concentration was multiplied by the total weight of homogenized sample (120 g) and dilution factors and divided by the weight of plant sample (20 g) to become the total nitrate content (μg/g).
Development of Functional Markers for Genetic Analysis
A total of 110 primers were designed related to specific enzymes involved in oxalate metabolism and downy mildew disease resistance gene (fungal disease common in spinach) (Correll et al., 1994). Already reported sequences of genes in the National Center for Biotechnology Information (NCBI) database were used to synthesize primers for amplification of genes. Sequences of functionally characterized genes for ten important enzymes, involved in oxalate metabolism (i.e., glycolate oxidase, malate dehydrogenase, ascorbate peroxidase, isocitrate lyase, malate synthase, citrate synthase, ascorbate oxidase, oxalyl-CoA synthetase, formate dehydrogenase, and oxalate-CoA ligase), were used to develop 100 primer pairs to assess the genetic variability of spinach accessions. In addition, the sequence of gene responsible for regulation of downy mildew resistance in spinach (i.e., glucose-6-phosphate dehydrogenase) was also used to develop ten primer pairs for genetic characterization of accessions (Supplementary Table 1).
Deoxyribonucleic Acid Extraction
Seedlings of spinach were used to extract DNA by cetyltrimethylammonium bromide (CTAB) method (Doyle, 1990). 1 g of the plant material was put in pestle and mortar and crushed in 600 μl CTAB buffer. After fine grinding, the mixture was transferred to 1.5 μl Eppendorf tube, homogenized, and placed in water bath at 60°C for 30 min. The mixture was cooled down at room temperature for 5–7 min. After normalizing temperature, 600 μl isoamyl alcohol chloroform solution was added in the mixture and centrifuged at 13,000 rpm for 10 min. The supernatant was carefully transferred to another tube for further processing. 300 μl ice cold isopropanol was added in the tube containing supernatant solution and the mixture was kept at −20°C for 30 min. After thawing, the mixture was centrifuged at 13,000 rpm for 10 min for the formation of DNA pellet at the bottom of tube. After pellet formation, the aqueous solution was carefully removed and the pellet was washed with 200 μl of 70% ethanol solution by centrifuging at 13,000 rpm for 10 min. After washing of DNA pellet, the solution was carefully removed from the tube and the pellet was allowed to dry for about 2 h, till the smell of ethanol was completely eliminated from the tube. After complete drying, the pellet was dissolved in 100 μl of TE buffer and vortex for 1 min. After complete dissolving, 1 μl RNase solution was added in the DNA extract and vortex for 1 min to degrade all the RNA from it. DNA quality and quantity were determined by gel electrophoresis and NanoDrop spectrophotometer technique, respectively (Abdel-Latif and Osman, 2017). The DNA extract was stored at −20°C till further use.
Polymerase Chain Reaction Amplification
Reaction mixture for PCR was prepared with the volume of 20 μl by making a solution of PCR master mix, nuclease-free water, DNA sample, and forward and reverse primers. The mixture was prepared on ice. Amplification of desired DNA template was performed in PCR thermal cycler (Bio-Rad, T100 Thermal Cycler).
Gel Electrophoresis
Agarose gel of 1% concentration was prepared to run the DNA samples. Solidified gel within its casting tray was placed in gel buffer tank (Wealtec Corporation, Model: mini GES), filled with TAE buffer (Cole and Tellez, 2002). Each DNA sample (5 μl) along with negative control was loaded in the respective wells of the gel. DNA ladder of 100 base pairs (bps), with 13 discrete DNA fragments ranging from 100 to 3,000 bp, was used as the standard to determine the DNA fragment size on gel electrophoresis. After running samples, the gel was taken out from the buffer tank and observed under UV light (UV transilluminator, Cleaver Scientific Ltd.) to visualize the variation among DNA bands.
Statistical Analysis
The nutritional and genetic variabilities of spinach accessions were assessed by statistical software RStudio version 1.3.1093. The ANOVA was performed for nutritional data with 0.1, 1, and 5% levels of significance. The broad sense heritability (h2b) was calculated as σ2g/σ2p, representing low and high values as < 50 and > 80%, respectively (Riaz et al., 2021). The genetic advance was calculated as, GA = h2b × σp × K, where σp and K represented phenotypic SD and standardized selection differential constant, respectively, at 5% selection intensity (Riaz et al., 2021). The genotypic coefficient of variation (GCV) and phenotypic coefficient of variation (PCV) were calculated by the method of Ogunniyan and Olakojo (2014). Genotypic and phenotypic correlations were calculated for all the nutritional characters and represented in the form of heat map diagram. Path coefficient analysis was performed to evaluate the direct and indirect effects of the most highly correlated nutrients on oxalate contents (Riaz et al., 2021). Functional marker data of spinach accessions, based on DNA banding pattern, was analyzed through correlation coefficient to determine the genetic relationship among accessions. The genetic diversity of accessions was determined by cluster analysis of functional marker data and represented in the form of heat map diagram (Cui et al., 2020).
Results and Discussion
Proximate Analysis
Proximate composition is essential to assess the nutritional significance of plants, especially leafy vegetables (Oyugi, 2016). Significant differences were observed among the nutritional contents of accessions. Most of the variation was observed among moisture and carbohydrate contents of accessions compared to rest of the contents. Moisture content of accessions ranged from 82 in accession 163,309 (from India) to 84.69% in accession 30,526 (from Pakistan). This high level of moisture in leafy vegetables could be helpful for the activity of a number of water-soluble enzymes required for good metabolism of plants (Iheanacho and Ubebani, 2009). In these results, above 80% moisture content was observed in all of the accessions, showing its active metabolism due to aqueous environment. Similar value of moisture was reported by Hussain et al. (2009) in spinach compared to other vegetables studied. The previous findings also revealed that the abovementioned accessions from Pakistan and India, in this study, were found to be good in terms of moisture contents and involvement in active metabolism.
Fibers are important dietary constituents in terms of good impact on digestion and waste elimination. Also, these can help in prevention of a number of digestive and cardiovascular disorders (Sodipo et al., 2000). Among all the studied spinach accessions, fiber contents ranged from 2.10 in accession 648,942 (China) to 2.34% in accession 30,802 (Pakistan). The results showed the importance of Asian spinach accessions due to its richness in fiber contents. The reason might be the better adaptability of these accessions in the local environment. Compared to a previous report about Cucurbita spp. (Nwofia et al., 2012), all the accessions of spinach were found to be rich in fibers. Fat contents also varied significantly from 0.51 to 0.75% among all the spinach accessions. High quantity of fats can lead toward obesity or overweight, but its moderate amount in diet is a very good source of energy (Sodipo et al., 2000; Erukainure et al., 2011). Present results have indicated few accessions with intermediate, but slightly different levels of fats. These accessions included 227,230 (Japan), 321,020 (Taiwan), 227,383 (Iran), and 647,858 (Georgia). This slight variation might be due to different origins of spinach accessions compared to that of local origin.
Ash contents in a plant represent its richness in mineral elements (Oyugi, 2016). Significant variation was found among ash contents of spinach accessions. Maximum percentage of ash was observed in accessions of China (648,942), Pakistan (163,310, 30,802), and Germany (266,926). These values were slightly lower than that reported by Galla et al. (2017) in Indian spinach. There was also variation in protein and calorie contents within the spinach accessions. Protein contents varied from 2.41 to 2.64%, while caloric values were in the range of 11.84 to 16.76 MJ/kg. Higher level of protein and calories shows the superiority of a crop on others (Oyugi, 2016). According to Ene-Obong (1992), this diet is considered to be nutritionally fit, containing high protein and caloric values. In these results, accessions 163,310 (from Pakistan) and 648,942 (from China) were found to be rich in terms of both the protein and caloric values. Like fats and proteins, digestible carbohydrates are also an important source of energy (Oyugi, 2016). In contrast to other two energy sources (fats and proteins), more variation was observed in terms of carbohydrate content. Few accessions of Pakistan (18,654, 32,122, and 19,339), India (163,309), China (648,942), United States (184,379), and Turkey (171,860) were found to be rich in carbohydrates (i.e., above 10%). Variation in nutritional composition of plants might be due to the influence of different environmental conditions (Odhav et al., 2007).
Mineral Analysis
Minerals are the important constituents of diet for good health and proper growth. Some of the minerals are integral parts of blood, teeth, muscles, and bones. Hence, their specific amount is required for proper functioning of body systems (Soetan et al., 2010). Though in trace quantities, micronutrients play a vital role in catalytic and balancing activity of plants (Jan et al., 2011). Leafy vegetables are considered as the richest source of minerals, especially iron, calcium, sodium, potassium, zinc, and manganese (Odhav et al., 2007). Significant variation was observed among mineral composition of spinach accessions. Among micronutrients, iron (Fe) content varied greatly, followed by manganese (Mn) and zinc (Zn) contents. Variation in Fe contents ranged from 27.1 (accession 184,379, United States) to 79.9 ppm (accession 181,964, Syria). Iron, being an important element of blood hemoglobin, requires for proper functioning of circulatory and central nervous system as well (Hussain et al., 2010). Accessions 32,618 (Pakistan) and 181,964 (Syria) were found to be richest in iron. Therefore, these accessions could be selected to fulfill iron requirement through leafy vegetables. Moreover, iron content of these accessions was found to be superior than that of previous reports about spinach (Kuti and Torres, 1996; Galla et al., 2017).
Zinc is required for better growth and proper functioning of immune and reproductive systems (Black, 2003; Hambidge, 2006). The quantity of zinc varied from 3 to 8 ppm among spinach accessions. Few accessions from Pakistan (163,310, 18,649, and 18,655) and one from India (163,309) were having high zinc contents. Highest value of Mn (8.2 ppm) was found in accession 163,310 (Pakistan), while lowest value of Mn (2.9 ppm) was found in accession 227,230 (Japan). The values of these micronutrients were even greater than that of previous reports about spinach and few other vegetables from Cucurbitaceae (Hussain et al., 2010), suggesting the use of these accessions for proper functioning of body.
Quantities of macronutrients also varied greatly among spinach accessions. Significant differences were found in terms of potassium (K) contents, followed by calcium (Ca) and sodium (Na) (Table 2). All these minerals have specific roles in human body. Calcium is an important dietary constituent as being helpful in osteogenesis (Hodges et al., 1978). Potassium is helpful in controlling the risk of stroke and hypertension (National Research Council [NRC], 1989). Highest value of potassium (4,000 ppm) was observed in accession 163,310 from Pakistan, while lowest value of potassium (3,240 ppm) was observed in accession 419,218 from Hong Kong. Calcium concentration varied from 550 in accessions 30,802 (Pakistan) and 227,230 (Japan) to 1,070 ppm in accession 163,310 (Pakistan). Few Pakistani accessions (163,310, 18,649, and 18,655) showed even greater Ca and K contents than that of a previous report about chaya and spinach (Kuti and Torres, 1996). Also, the same accessions showed higher contents of Na. Sodium content among accessions ranged from 420 to 876 ppm. In contrast to one previous report (Nwofia et al., 2012) about spinach of Nigeria, quantity of Na was comparatively lesser. The variation in nutritional contents might be due to the influence of environmental factors, harvest time, or farming practices, which ultimately leads toward the overall genetic diversity (Odhav et al., 2007).
Table 2. ANOVA for nutritional contents of spinach accessions.
Oxalate and Nitrate Contents
Leafy vegetables are rich in oxalates and nitrates. However, its higher concentration greatly affects the mineral absorption. In this study, maximum oxalate contents (13.16 g/100 g dw) were found in accession 30,526 from Pakistan. The mineral composition of this accession was comparatively lower than other accessions, while another accession from Pakistan (163,310) was found to be lesser in terms of oxalate contents, but rich in minerals (Na = 876 ppm, Ca = 1,070 ppm, K = 4,000 ppm, Zn = 80 ppm, Fe = 51 ppm, and Mn = 8.2 ppm) and nitrates (2,549 ppm) as well. This also represented the effect of varied oxalate concentration on mineral and nitrate contents. Nitrates and oxalates have been considered as the antinutritional phytochemicals. However, the crops with good oxalate metabolism and better nitrogen use efficiency can reduce the harmful impacts of these antinutritional phytochemicals (Uusiku et al., 2010). In this study, the accessions with greater macro- and micronutrients represented lesser oxalates, showing its better metabolism, which were also assessed through genetic marker study. Moreover, the accessions should have better nitrogen use efficiency that represented with positive relationship of nitrates with other nutritional factors. These findings were also supported by the facts reported by Natesh et al. (2017).
Genetic Variability
The nutritional data was evaluated for variance analysis and the values of mean, variance, coefficient of variation, heritability, and genetic advance are given in Table 2 to represent the significance of various nutritional parameters. Significant variation was observed among nutritional contents of different spinach accessions. The ANOVA showed nearly equal phenotypic and genotypic variances compared to environmental variance, which was very low for all the nutritional traits. Similar observations were reported by Riaz et al. (2021) showing comparatively lesser influence of environment on the nutritional values. The phenotypic and genotypic coefficients of variations were higher for nitrates, manganese, zinc, iron, sodium, and carbohydrates. Low phenotypic and genotypic coefficients of variation were observed for moisture and oxalates. However, the environmental coefficient of variation was lower for most of the nutritional traits. Genetic heritability was higher for calories and sodium.
Correlation Analysis
Correlation analysis is a technique that is helpful in selection of positive correlated traits by giving positive or negative impacts on its correlated traits (Manggoel et al., 2012). Moreover, this analysis is helpful when highly heritable traits are linked with important nutritional contents of the leafy crops. A heat map of genotypic and phenotypic correlation coefficients was constructed among various nutritional parameters (Figure 1). The genotypic correlation coefficients were very closely correlated to the phenotypic correlation coefficients. Maximum positive correlation (97%) was observed among calcium, manganese, and zinc minerals, while negative correlation (50%) was found between moisture and carbohydrates.
Figure 1. Heat map indicating phenotypic and genotypic correlations of nutritional contents. Dark blue color is the representation of greater correlation level.
Genotypic Pathway Analysis
Besides the importance of correlation coefficients in determining the linkage of traits, the information could be a little misleading, if the observed higher correlation coefficient is due to the indirect effects of traits (Manggoel et al., 2012). Therefore, a path coefficient analysis, in addition, was performed to determine the direct and indirect effects of interrelated nutritional traits (Table 3). Among all the observed nutritional traits, seven highly linked traits with greater values of direct effects on oxalates were considered to build a path diagram, representing its relationship of direct and indirect effects with that of oxalate contents (Figure 2). The oxalate concentration in spinach accessions was found to be greatly linked with calcium concentration, as these two constituents are usually found as calcium oxalates in leafy vegetables (Ferreira et al., 2018). The least effect on oxalates was shown by proteins of spinach accessions. On the other hand, calcium and manganese had indirect effects on oxalate contents of spinach.
Table 3. Pathway analysis showing direct and indirect effects of nutritional content on oxalates of spinach accessions.
Figure 2. Path diagram representing direct and indirect effects of nutrients on oxalates.
Functional Marker Analysis
Oxalates are among the active metabolic components of leafy vegetables, especially spinach. It readily affects the absorption of mineral ions in plant tissue. The active metabolic cycle of oxalates is regulated by a number of metabolic enzymes, which maintain its concentration in plant tissue (Horner et al., 2005). Therefore, this study included the use of functional molecular markers related to the genes involved in oxalate metabolism in order to relate the oxalate contents and mineral availability in spinach accessions.
Leaves of spinach are rich in minerals, especially iron, calcium, potassium, and sodium. Potassium, sodium, and calcium minerals are usually found in association with oxalates either as soluble salts or as insoluble crystals of calcium oxalate (Ferreira et al., 2018). Oxalates are largely distributed in spinach, accounting for about 5–15% of dry weight of spinach leaves. A number of functional roles have been reported for oxalates in calcium regulation, ion balancing (Na, K), and detoxification of heavy metals, especially in case of leafy vegetables (Botstein et al., 1980; Amom and Nongdam, 2017; Dhutmal et al., 2018). Apart from positive roles, excessive accumulation of oxalates in plants could lead toward negative impacts on human health (Cai et al., 2018). Therefore, a balance should be maintained between anabolic and catabolic pathways of oxalate. A number of enzymes are involved in anabolic and catabolic pathways of oxalate. The presence or absence of its genes could be considered a criterion for selection of improved varieties. Another important consideration for breeders is the development of disease-resistant varieties. Significant variation was observed among genetic profile of accessions, especially with oxalate-related genetic markers.
Allele Number and Size
Spinach accessions showed a great deal of genetic variation based on functional molecular markers. Figure 3 displays the variation in banding pattern of oxalate-CoA ligase gene among spinach accessions. The allelic variation was represented in the form of differences in allele number and sizes (Supplementary Table 1). Number of alleles ranged from 2 to 21. Maximum number of alleles (21) was represented by the DNA product form from primer pair OxySyn-1, followed by AscPerox with 20 allele number. Minimum number of alleles was formed as a result of amplification by primer pairs GlyOx-7 and GlyOx-9. Size of amplified product ranged from 73 to 994 bp. Maximum product size (994 bp) was formed as a result of amplification by primer set AscOx-3, while minimum product length (73 bp) was shown by primer set G6PDH-2.
Figure 3. Allelic profile of oxalate-CoA ligase among spinach accessions generated by primer pair OxLig-1. L = 100 bp ladder, 1 = 321,021, 2 = 19,339, 3 = 30802, 4 = 18,649, 5 = 18,654, 6 = 30,518, 7 = 32,618, 8 = 227,230, 9 = 379,549, 10 = 18,650, 11 = 32,122, 12 = 508,504, 13 = 647,858, 14 = 604,785, 15 = 163,309, 16 = 163,310, 17 = 419218, 18 = 212,119, 19 = 647853, 20 = 30,526, 21 = 18,655, 22 = 266,926, 23 = 261,789, 24 = 181,964, 25 = 179,593, 26 = 227,383, 27 = 648,942, 28 = 648,944, 29 = 184,379, 30 = 171,860.
In these results, genetic variation was assessed in terms of DNA banding pattern. Out of 130 bands, 115 bands were found to be polymorphic representing 89% polymorphism. Most of the primers related to oxalate biosynthesis genes were found to be highly polymorphic. Functional marker AscPerox-1 was found to be highly polymorphic among all with 0.92 polymorphic information content (PIC) value and heterozygosity. It was followed by IsoLy-1 with 0.91 PIC value and heterozygosity. Minimum value of PIC (0.32) and heterozygosity (0.26) was represented by primer pair AscOx-2 (Supplementary Table 1). In a similar study, Avsar (2011) used twenty-five sequence-related amplified polymorphism (SRAP) markers for genetic characterization of ninety-five Turkish spinach accessions. The banding pattern showed 123 amplified bands, out of which 67 bands were polymorphic giving rise about 54% genetic polymorphism. In another study, Goktay (2015) characterized 176 spinach accessions through 15 simple sequence repeat (SSR) markers with 98% polymorphism. The banding pattern showed 57 polymorphic bands out of 58 bands. Hu et al. (2007) reported genetic variability of 48 spinach germplasm based on 59 target region amplification polymorphism (TRAP) markers. However, the obtained polymorphism was low (19.5%) compared to other reports. The comparison revealed the importance of functional markers as more polymorphic. Similar opinion was stated by Salgotra and Stewart (2020) about importance of functional markers due to its specificity to a particular trait.
Correlation and Cluster-Based Relationship
Table 4 represents the correlation of 30 spinach accessions based on observed genetic polymorphism. Accessions 508,504 and 184,379 were found to be highly correlated representing 91% similarity. Minimum positive correlation (1%) was shown by accessions 18,649 (with 32,618 and 181,964), 18,650 (with 379,549 and 648,942), and 163,309 (with 32,618 and 181,964). Few accessions were found to be negatively correlated in terms of genetic polymorphism. Maximum negative correlation (17%) was represented by accessions 266,926 (with 212,119 and 6,489,440) and 32,618 (with 181,964, 18,649, 18,654, and 30,526), while minimum correlation of about 1% was found in accessions 227,383, 508,504, 227,383, and 19,339 with reference to each other. Also, similar negative correlation was shown by accession 18,649 with 32,618 and 181,964.
Table 4. Correlation of thirty spinach accessions based on functional marker analysis.
Cluster analysis was performed among spinach accessions based on genetic profile and represented in the form of heat map (Figure 4). This clustering was based on homozygosity and heterozygosity of functionally characterized genes among all the accessions. Genetic marker data divided the accessions into four clusters based on allelic variation. First cluster represented accessions from Hong Kong (419,218), Pakistan (18,650, 32,122), India (163,309), Georgia (647,858), and France (261,789). Among these accessions, most of the variability was found among the genes involved in oxalate biosynthesis, including glycolate oxidase, malate dehydrogenase, ascorbate peroxidase, and isocitrate lyase. Cluster 2 represented the accessions with comparatively greater heterozygosity level in most of the genes. The accessions were belonging to Pakistan (30,802, 19,339), Belgium (179,593), Macedonia (379,549), China (648,942), Japan (227,230), and Iran (227,383). Cluster 3 was characterized with the accessions having comparatively less variability in oxalate-related genes and greater variability in downy mildew resistance gene. The cluster included the accessions from South Korea (508,504), Georgia (647,853), Afghanistan (212,119), China (648,944), Taiwan (321,020), and United States (184,379). Cluster 4 was the intermediate among all, consisting of most of the accessions from Pakistan (30,526, 18,649, 18,554, 18,655, 32,618, 163,310, and 30,518) and three from Germany (266,926), Mongolia (604,785), Turkey (171,860), and Syria (181,964). This group was at the borderline of homozygosity and heterozygosity. Most variability contributing genes included formate dehydrogenase, ascorbate peroxidase, citrate synthase, and ascorbate oxidase.
Figure 4. Heat map clustering of spinach accessions based on functionally characterized genes. X1 = Glycolate oxidase, X2 = Oxalyl CoA Synthetase, X3 = Glucose-6-Phosphate Dehydrogenase, X4 = Malate Dehydrogenase, X5 = Ascorbate Peroxidase, X6 = Iso-citrate lyase, X7 = Formate Dehydrogenase, X8 = Oxalate CoA Ligase, X9 = Malate Synthase, X10 = Citrate Synthase, X11 = Ascorbate Oxidase. Dark blue color is the representation of greater heterozygosity level.
The cluster-based relationship revealed grouping of accessions irrespective to its geographical locations. However, the variation was found to be more linked with nutritional variation of accessions. Similar results were reported by Hu et al. (2007) in its genetic marker study. Accessions from same region were clustered into the different groups, showing its independence to observed genetic variability. The reason of this variation might be due to the introduction of similar germplasm in different areas, resulting in genetic drift or natural selection. Heat map clustering represented grouping of accessions based on presence of various genes. Accessions from Pakistan, India, Georgia, Belgium, Macedonia, China, Japan, Iran, Germany, Mongolia, and Turkey were found to be more variable in terms of disease resistance gene. Besides, the accessions from India, Pakistan, Georgia, Belgium, and China were also found to be rich in minerals, especially sodium, potassium, zinc, and manganese.
Conclusion
It was concluded from the findings that the minerals are required not only for proper growth and development of plants, but also for making them resistant to various diseases. The accessions from India, Pakistan, Iran, and Belgium showed high variability in terms of oxalate metabolism-related genes. The accessions from Iran and Belgium were nutritionally more variable in terms of protein, ash, and fiber contents, along with minerals. This represented the correlation of good oxalate metabolism in leafy crops with mineral absorption. These findings were also supported by the report of Sun et al. (2019) about oxalate–mineral relationship. Globally, the accessions from Asia showed more diversity compared to the accessions from other continents. Similarly, the report of Kuwahara et al. (2014) about SSR marker diversity showed East Asian spinach germplasm as the most diverse compared to European and American germplasm. Genetically diverse accessions have the natural ability to bear a number of environmental stresses. Therefore, the abovementioned diverse accessions could be considered as useful genotypes in future breeding programs.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
MR conducted the research and wrote – original draft. ZY supervised the research. AD provided the resources for research. MM and NR interpreted data for formal analysis. AA and AY reviewed the original draft. AA, MO, and HA edited and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
We would like to extend their appreciation to the Researchers Supporting Project Number (RSP-2021/176), King Saud University, Riyadh, Saudi Arabia for providing funds.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank Lahore College for Women University, Lahore, Pakistan and Ayub Agricultural Research Institute for providing workplace for the nutritional analysis of spinach accessions. We would also like to extend their appreciation to the International Research Support Initiative Program, Higher Education Commission, Pakistan and Researchers Supporting Project Number (RSP-2021/176), King Saud University, Riyadh, Saudi Arabia for providing funds for the molecular study on spinach accessions.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.889604/full#supplementary-material
References
Abdel-Latif, A., and Osman, G. (2017). Comparison of three genomic DNA extraction methods to obtain high DNA quality from maize. Plant Methods 13, 1–9. doi: 10.1186/s13007-016-0152-4
Abu Al-Qumboz, M. N., and Abu-Naser, S. S. (2019). Spinach expert system: diseases and symptoms. Int. J. Acad. Inf. Syst. Res. 3, 16–22.
Amom, T., and Nongdam, P. (2017). The use of molecular marker methods in plants: a review. Int. J. Curr. Res. Rev. 9, 1–9.
Anon (1970). Oxygen Bomb Calorimetry and Combustion Methods. Manual No. 130. Moline, IL: Parr Instrument Co.
AOAC (1990). Official Methods of Analysis. Association of Analytical Chemists, 15th Edn. Washington, DC: AOAC.
AOAC (2000). Official Methods of Analysis. Association of Analytical Chemists, 17th Edn. Arlington, VA: AOAC.
Arif, M., Tasneem, M., Bashir, F., Yaseen, G., and Anwar, A. (2017). Evaluation of different levels of potassium and zinc fertilizer on the growth and yield of wheat. Int. J. Biosens. Bioelectron. 3, 1–5. doi: 10.15406/ijbsbe.2017.03.00057
Avsar, B. (2011). Genetic Diversity of Turkish Spinach Cultivars (Spinacia oleracea L.). Master’s thesis. Izmir: Izmir Institute of Technology.
Black, R. E. (2003). Zinc deficiency, infectious disease and mortality in developing world. J. Nutr. 133, S1485–S1489. doi: 10.1093/jn/133.5.1485S
Botstein, D., White, R. L., Skolnick, M., and Davis, R. W. (1980). Construction of a genetic linkage map in man using restriction fragment length polymorphism. Am. J. Hum. Genet. 32, 314–331.
Cai, X., Ge, C., Xu, C., Wang, X., Wang, S., and Wang, Q. (2018). Expression analysis of oxalate metabolic pathway genes reveals oxalate regulation patterns in spinach. Molecules 23:1286. doi: 10.3390/molecules23061286
Cole, K. D., and Tellez, C. M. (2002). Separation of large circular DNA by electrophoresis in agarose gels. Biotechnol. Prog. 18, 82–87. doi: 10.1021/bp010135o
Correll, J. C., Morelock, T. E., Black, M. C., Koike, S. T., Brandenberger, L. P., and Dainello, F. J. (1994). Economically important diseases of spinach. Plant Dis. 78, 653–660. doi: 10.1111/jam.12541
Cousins, R. J. (2006). “Zinc,” in Present Knowledge in Nutrition, 9th Edn, Vol. 1. eds B. A. Bowman and R. M. Russell (Washington, DC: ILSI Press), 445–457.
Cui, J., Yang, Y., Luo, S., Wang, L., Huang, R., Wen, Q., et al. (2020). Whole-genome sequencing provides insights into the genetic diversity and domestication of bitter gourd (Momordica spp.). Hortic. Res. 7:85. doi: 10.1038/s41438-020-0305-5
Dhutmal, R. R., Mundhe, A. G., and More, A. W. (2018). Molecular marker techniques: a Review. Int. J. Curr. Microbiol. Appl. Sci. 6, 816–825.
Egbi, G., Gbogbo, S., Mensah, G. E., Glover-Amengor, M., and Steiner-Asiedu, M. (2018). Effect of green leafy vegetables powder on anaemia and vitamin-A status of Ghanaian school children. BMC Nutr. 4:27. doi: 10.1186/s40795-018-0235-x
Ene-Obong, H. N. (1992). Nutritional Evaluation, Composition Pattern and Processing of Underutilized Traditional Foods with Particular Reference to the African Yam Beans (Sphenostylis sternocarpa). Ph.D. thesis. Nsukka: University of Nigeria.
Erukainure, O. L., Oke, O. V., Ajiboye, A. J., and Okafor, O. Y. (2011). Nutritional qualities and phytochemical constituents of Clerodendrum volubile, a tropical non-conventional vegetable. Int. Food Res. J. 18, 1393–1399.
Ferreira, J. F., Sandhu, D., Liu, X., and Halvorson, J. J. (2018). Spinach (Spinacia oleracea L.) response to salinity: nutritional value, physiological parameters, antioxidant capacity, and gene expression. Agriculture 8:163.
Fresco, L. O., and Baudoin, W. O. (2004). “Food and nutrition security towards human security,” in Food Security and Vegetables: A Global Perspective, eds P. Nath, P. B. Gaddagimath, and O. P. Dutta (Bangalore: Nath Agricultural Science Foundation), 7–42.
Galla, N. R., Pamidighantam, P. R., Karakala, B., Gurusiddaiah, M. R., and Akula, S. (2017). Nutritional, textural and sensory quality of biscuits supplemented with spinach (Spinacia oleracea L.). Int. J. Gastron. Food Sci. 7, 20–26. doi: 10.1016/j.ijgfs.2016.12.003
Goktay, M. (2015). Characterization of World Spinach Genetic Collection by Using Molecular Markers. Ph.D. thesis. Izmir: Izmir Institute of Technology.
Grant, C. A., Clarke, J. M., Duguid, S., and Chaney, R. L. (2008). Selecting and breeding of plant cultivars to minimize cadmium accumulation. Sci. Total Environ. 390, 301–310. doi: 10.1016/j.scitotenv.2007.10.038
Gunnars, K. (2019). Spinach 101: Nutrition Facts and Health Benefits. Health-line Newsletter. Available online at: https://www.healthline.com/nutrition/foods/spinach (accessed February 04, 2022).
Hodges, R. E., Sauberlich, H. E., Canham, J. E., Wallace, D. L., Rucker, R. B., Mejia, L. A., et al. (1978). Hematopoietic studies in vitamin A deficiency. Am. J. Clin. Nutr. 31, 876–885. doi: 10.1093/ajcn/31.5.876
Horner, H. T., Cervantes-Martinez, T., Healy, R., Reddy, M. B., Deardorff, B. L., Bailey, T. B., et al. (2005). Oxalate and phytate concentrations in seeds of soybean cultivars [Glycine max (L.) Merr.]. J. Agric. Food Chem. 53, 7870–7877. doi: 10.1021/jf051193i
Hu, J., Mou, B., and Vick, B. A. (2007). Genetic diversity of 38 spinach (Spinacia oleracea L.) germplasm accessions and 10 commercial hybrids assessed by TRAP markers. Genet. Resour. Crop Evol. 54, 1667–1674. doi: 10.1007/s10722-006-9175-4
Hussain, J., Khan, A. L., Rehman, N., Hamayun, M., Shah, T., Nisar, M., et al. (2009). Proximate and nutrient analysis of selected vegetable species: a case study of Karak region, Pakistan. Afr. J. Biotechnol. 8, 2725–2729. doi: 10.3923/pjn.2009.620.624
Hussain, J., Rehman, N. U., Khan, A. L., Hamayun, M., Hussain, S. M., and Shinwari, Z. K. (2010). Proximate and essential nutrients evaluation of selected vegetables species from Kohat region, Pakistan. Pak. J. Bot. 42, 2847–2855.
Iheanacho, K., and Ubebani, A. C. (2009). Nutritional composition of some leafy vegetable consumed in Imo-State, Nigeria. J. Appl. Sci. Environ. Manage. 13, 35–38. doi: 10.4314/jasem.v13i3.55349
Jan, G., Kahan, M., Ahmad, M., Iqbal, Z., Afzal, A., Afzal, M., et al. (2011). Nutritional analysis, micronutrients and chlorophyll contents of Cichorium intybus L. J. Med. Plants Res. 5, 2452–2456.
Joshi, V., Penalosa, A., Joshi, M., and Rodriguez, S. (2021). Regulation of Oxalate Metabolism in Spinach Revealed by RNA-Seq-Based Transcriptomic Analysis. Int. J. Mol. Sci. 22:5294. doi: 10.3390/ijms22105294
Karthika, K. S., Rashmi, I., and Parvathi, M. S. (2018). “Biological functions, uptake and transport of essential nutrients in relation to plant growth,” in Plant Nutrients and Abiotic Stress Tolerance, eds M. Hasanuzzaman, M. Fujita, H. Oku, K. Nahar, and B. Hawrylak-Nowak (Berlin: Springer), 1–49. doi: 10.1007/978-981-10-9044-8_1
Kuti, J. O., and Torres, E. S. (1996). “Potential nutritional and health benefits of tree spinach,” in Progress in New Crops, ed. J. Janick (Arlington, VA: ASHS Press), 516–520.
Kuwahara, K., Suzuki, R., Ito, Y., Mikami, T., and Onodera, Y. (2014). An analysis of genetic differentiation and geographical variation of spinach germplasm using SSR markers. Plant Genet. Resour. 12, 185–190. doi: 10.1017/S1479262113000464
Lau, W. C., Rafii, M. Y., Ismail, M. R., Puteh, A., Latif, M. A., and Ramli, A. (2015). Review of functional markers for improving cooking, eating, and the nutritional qualities of rice. Frontiers in Plant Science 6:832. doi: 10.3389/fpls.2015.00832
Manggoel, W., Uguru, M. I., Ndam, O. N., and Dasbak, M. A. (2012). Genetic variability, correlation and path coefficient analysis of some yield components of ten cowpea [Vigna unguiculata (L.) Walp] accessions. J. Plant Breed. Crop Sci. 4, 80–86. doi: 10.5897/JPBCS12.007
Mehri, F., Heshmati, A., Moradi, M., and Khaneghah, A. M. (2021). The concentration and health risk assessment of nitrate in vegetables and fruits samples of Iran. Toxin Rev. 40, 1215–1222. doi: 10.1080/15569543.2019.1673424
Natesh, H. N., Abbey, L., and Asiedu, S. K. (2017). An overview of nutritional and antinutritional factors in green leafy vegetables. Hortic. Int. J. 1:00011. doi: 10.15406/hij.2017.01.00011
Nwofia, G. E., Victoria, N. N., and Blessing, K. N. (2012). Nutritional variation in fruits and seeds of pumpkins (Cucurbita sp.) accessions from Nigeria. Pak. J. Nutr. 11, 946–956. doi: 10.3923/pjn.2012.946.956
Odhav, B., Beekrum, S., Akula, U. S., and Baijnath, H. (2007). Preliminary assessment of nutritional value of traditional leafy vegetables in KwaZulu-Natal, South Africa. J. Food Compos. Anal. 20, 430–435. doi: 10.1016/j.jfca.2006.04.015
Ogunniyan, D. J., and Olakojo, S. A. (2014). Genetic variation, heritability, genetic advance and agronomic character association of yellow elite inbred lines of maize (Zea mays L.). Niger. J. Genet. 28, 24–28. doi: 10.1016/j.nigjg.2015.06.005
Otitoju, G. T. O. (2009). Effect of dry and wet milling processing techniques on the nutrient composition and organoleptic attributes of fermented yellow maize (Zea mays). Afr. J. Food Sci. 3, 113–116.
Oyugi, J. O. (2016). Analysis of Proximate, Micronutrients and Determination of Phytochemicals in Selected Medicinal Plants in Mbita-Homabay County. Ph.D. thesis. Nairobi: Kenyatta University.
Pearson, D. (1976). The Chemical Analysis of Foods, 7th Edn. (Edinburgh: Churchill Livingstone), 488–496.
Popoola, O. E., Doherty, V. F., Odusami, J. O., and Durowoju, O. S. (2014). Oxalate content of some Nigerian tubers using titrimetric and UV-spectrophotometric methods. Acad. J. Agric. Res. 2, 054–057.
Rashid, M., Yousaf, Z., Haider, M. S., Khalid, S., Rehman, H. A., Younas, A., et al. (2014). Genetic diversity of functional food species Spinacia oleracea L. by protein markers. Nat. Prod. Res. 28, 782–787. doi: 10.1080/14786419.2014.881359
Riaz, N., Yousaf, Z., Munawar, M., Younas, A., Rashid, M., Aftab, A., et al. (2021). Development of carrot nutraceutical products as an alternative supplement for the prevention of nutritional diseases. Front. Nutr. 8:787351. doi: 10.3389/fnut.2021.787351
Sabaghnia, N., Asadi-Gharneh, H. A., and Janmohammad, M. (2014). Genetic diversity of spinach (Spinacia oleracea) landraces collected in Iran using some morphological traits. Acta Agric. Slov. 103, 101–111. doi: 10.14720/aas.2014.103.1.11
Salgotra, R. K., and Stewart, C. N. (2020). Functional markers for precision plant breeding. Int. J. Mol. Sci. 21:4792. doi: 10.3390/ijms21134792
Shi, A., Mou, B., and Correll, J. C. (2016). Association analysis for oxalate concentration in spinach. Euphytica 212, 17–28. doi: 10.1007/s10681-016-1740-0
Sodipo, O. A., Akiniyi, J. A., and Ogunbamosu, J. U. (2000). Studies on certain on certain characteristics of extracts of bark of Pansinystalia macruceras (K schemp) picrre Exbeille. Glob. J. Pure Appl. Sci. 6, 83–87. doi: 10.4314/gjpas.v6i1.16082
Soetan, K. O., Olaiya, C. O., and Oyewole, O. E. (2010). The importance of mineral elements for humans, domestic animals and plants-A review. Afr. J. Food Sci. 4, 200–222.
Sun, Q., Li, J., Finlay, R. D., and Lian, B. (2019). Oxalotrophic bacterial assemblages in the ectomycorrhizosphere of forest trees and their effects on oxalate degradation and carbon fixation potential. Chem. Geol. 514, 54–64. doi: 10.1016/j.chemgeo.2019.03.023
Tang, L., Hamid, Y., Sahito, Z. A., and Gurijala, H. (2019). Evaluation of variation in essential nutrients and hazardous materials in spinach (Spinacia oleracea L.) genotypes grown on contaminated soil for human consumption. J. Food Compos. Anal. 79, 95–106. doi: 10.1016/j.jfca.2019.03.012
Uusiku, N. P., Oelofse, A., Duodu, K. G., Bester, M. J., and Faber, M. (2010). Nutritional value of leafy vegetables of sub-Saharan Africa and their potential contribution to human health: a review. J. Food Compos. Anal. 23, 499–509. doi: 10.1016/j.jfca.2010.05.002
Varshney, R. K., Graner, A., and Sorrells, M. E. (2005). Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 23, 48–55. doi: 10.1016/j.tibtech.2004.11.005
Verma, S. (2018). A study on medicinal herb Spinacia oleraceae Linn: Amaranthaceae. J. Drug Deliv. Ther. 8, 59–61. doi: 10.22270/jddt.v8i4.1767
Keywords: accession, breeding, crop improvement, nitrate, mineral uptake, oxalate, spinach
Citation: Rashid M, Yousaf Z, Din A, Munawar M, Aftab A, Riaz N, Younas A, Alaraidh IA, Okla MK and AbdElgawad H (2022) Assessment of Mineral Nutrient Efficiency in Genetically Diverse Spinach Accessions by Biochemical and Functional Marker Strategies. Front. Plant Sci. 13:889604. doi: 10.3389/fpls.2022.889604
Received: 04 March 2022; Accepted: 26 April 2022;
Published: 30 May 2022.
Edited by:
Adnan Noor Shah, Khwaja Fareed University of Engineering and Information Technology (KFUEIT), PakistanReviewed by:
Muhammad Ayaz, Lithuanian Research Centre for Agriculture and Forestry, LithuaniaAli Raza, Fujian Agriculture and Forestry University, China
Rajesh Singhal, Indian Grassland and Fodder Research Institute (ICAR), India
Copyright © 2022 Rashid, Yousaf, Din, Munawar, Aftab, Riaz, Younas, Alaraidh, Okla and AbdElgawad. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Madiha Rashid, Ym90YW5pc3Rwa0B5YWhvby5jb20=; Zubaida Yousaf, enViYWlkYV95b3VzYWZAeWFob28uY29t; enViYWlkYS55b3VzYWZAbGN3dS5lZHUucGs=
†These authors have contributed equally to this work