InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Agricultural and Biological Sciences » "Cotton Research", book edited by Ibrokhim Y. Abdurakhmonov, ISBN 978-953-51-2737-6, Print ISBN 978-953-51-2736-9, Published: November 9, 2016 under CC BY 3.0 license. © The Author(s).

Chapter 4

Nutritional Status of Cotton Plant Assessed by Compositional Nutrient Diagnosis (CND)

By Ademar Pereira Serra, Marlene Estevão Marchetti, Manoel Carlos Gonçalves, Simone Cândido Ensinas, Bendaly Labaied Mouna, Eulene Francisco da Silva, Elaine Reis Pinheiro Lourente, Anamari Viegas de Araujo Motomiya, Alessandra Mayumi Tokura Alovisi and Flávia de Araújo Matos
DOI: 10.5772/64388

Article top


Graphic representation of the observed frequency of yields. The criteria for establishing the high-yield subpopulations were plots with (A) yield above average (4380 kg ha−1), (B) yield above average + 2/3 standard deviation (above 4650 kg ha−1), (C) yield above average + 1 standard deviation (above 4785 kg ha−1), and (D) yield above average + 4/3 standard deviation (above 4920 kg ha−1).
Figure 1. Graphic representation of the observed frequency of yields. The criteria for establishing the high-yield subpopulations were plots with (A) yield above average (4380 kg ha−1), (B) yield above average + 2/3 standard deviation (above 4650 kg ha−1), (C) yield above average + 1 standard deviation (above 4785 kg ha−1), and (D) yield above average + 4/3 standard deviation (above 4920 kg ha−1).
Distribution of principal component analysis (PCA) for row-centered log ratios. PC = principal components. Database from Serra [19].
Figure 2. Distribution of principal component analysis (PCA) for row-centered log ratios. PC = principal components. Database from Serra [19].

Nutritional Status of Cotton Plant Assessed by Compositional Nutrient Diagnosis (CND)

Ademar Pereira Serra1, Marlene Estevão Marchetti2, Manoel Carlos Gonçalves2, Simone Cândido Ensinas3, Bendaly Labaied Mouna4, Eulene Francisco da Silva5, Elaine Reis Pinheiro Lourente2, Anamari Viegas de Araujo Motomiya2, Alessandra Mayumi Tokura Alovisi2 and Flávia de Araújo Matos2
Show details


The use of compositional nutrient diagnosis (CND) to assess the nutritional status of cotton crop is quite important to improve knowledge on plant nutritional requirement and assist the fertilizer recommendation. The aim of this chapter is to introduce the possibility of using CND for cotton crop. This method has scarcely been used to assess the nutritional status of cotton plant although a few results have indicated that it can be promising. In fact, CND methodology seems to be better in the nutritional diagnosis than traditional methods such as sufficient range (SR) and critical value approach (CVA). Its efficiency has increased with the possibility of applying multivariate analysis, principal component analysis (PCA), canonical correlation, and so on. The application of PCA possibility to note some interactions among the nutrients is important for understanding the dynamics of nutrients in plants.

Keywords: multivariate analysis, nutritional balance, nutrient concentration, PCA, plant nutrition, soil science, Gossypium hirsutum L

1. Introduction

The knowledge of nutritional status of cotton is essential to understand plant nutritional requirement and to help the fertilizer recommendations to improve crop production. Without understanding cotton nutrient requirement, it is difficult to know about plant response to fertilizers. The traditional methods to assess plant nutrient status are the sufficiency range (SR) and critical value approach (CVA); however, these methods do not allow the diagnosis of the nutritional balance of the plant, and their usage is dependent on sampling time and plant development stage [1]. Additionally, in SR and CVA, it is not possible to know which nutrient is more or less required by plant, because there is no interaction among nutrients in these traditional diagnoses.

In order to improve nutritional diagnosis by conventional methods (CVA and SR), Beaufils [2] developed the diagnosis and recommendation integrated system (DRIS), which was based on “minimum law.” Through the DRIS index, it is possible to observe the plant nutritional requirement, and the most deficient nutrient is limiting the plant development and production. It is possible to consider DRIS as an early attempt to analyze the compositional space of leaf analysis results [3]. Parent and Dafir [4] developed the compositional nutrient diagnosis (CND) that is conformed to be a multivariate analysis and can combine other refined analysis to improve the nutritional diagnosis. CND uses the row-centered log ratio instead of dual ratios between nutrients [5]. CND consists in the relation between each individual nutrient concentration with the geometric mean of the other nutrients in composition. Thus, CND comprises multiple interactions among all nutrients in diagnosis. Therefore, CND method has more advantages than DRIS method due to the interaction among all nutrients instead of dual ratios mentioned above [4].

For cotton crop, Serra et al. [6, 7] have compared CVA and SR diagnosis with CND method and observed that CND shows more precise diagnosis results than CVA and SR, and the authors have reported that the use of SR may generate wide nutrient range that can disturb the accurate diagnosis.

The aim of this chapter is to introduce the CND applications to assess the nutritional status of cotton plant as a more precise method of plant diagnosis than the widely applied traditional CVA and SR methods.

2. Nutritional diagnosis in cotton plant

CVA or SR method has been defined as an interventionist approach because the production variability is explained by the variation on nutrient supply or availability of the nutrients under analysis. In CVA and SR, the other production factors are remained under unlimited levels. The understanding of the principles, considered by different diagnosis methods, and the comparison of their results are important to the application of these diagnosis tools.

In cotton and most grain crops, usually, the nutritional diagnoses are conducted using SR and CVA, which measure only individual leaf nutrient concentration of a nutrient in diagnosis without any relationship with other nutrients in the sample of leaf tissue. These traditional methods do not report the nutritional balance in plant tissue. Bates [1] noted the nutritional diagnosis through CVA was affected by the interaction of nutrients in plant tissue and other factors related to plant growth depending on the age of plant.

In the last decade, the interests for bivariate and multivariate methods of nutritional diagnosis, as DRIS and CND, increased. DRIS is based on the dual ratio between nutrients, as N/P, N/K, etc., which is considered to have a bivariate relation. DRIS method allows the assessment of the nutritional equilibrium of a plant, ranking the nutrient contents in relative order of nutritional requirement, from the most deficient nutrient to the most excessive [2].

In DRIS results, in order to express the relative nutrient balance into the plant tissue, the DRIS index is calculated through the comparison of dual nutrient ratios (N/P, P/K, K/Ca, Ca/Mg, etc.) in the sample with the DRIS norms from a reference population (N/P´, P/K´, K/Ca´, Ca/Mg´, etc.). Besides, in the DRIS, there are other bivariate methods such as modified DRIS (M-DRIS) [8] and multivariate method like the CND, which was developed by Parent and Dafir [4].

3. Steps to develop CND norms

3.1. Obtaining leaf nutrient concentration and cotton crop yield to compose the database

The first step to develop CND norms is to obtain leaf nutrient concentration and cotton yield (database). The leaf sampling may be conducted in commercial or experimental cotton field, but it is recommended that the number of samples cannot be under 30 samples of complete leaves (blade + petiole) per plot, which must be combined to form just 1 composed sample [9].

One leaf sample per plant should be taken from the fifth leaf on the main stem, during flowering of the cotton crop (stages F1–F4), according to classification of Marur and Ruano [10]. The cotton yield is assessed at the end of plant cycle, which is usually accomplished with a combine harvester in commercial field or manually in experimental plots.

In the leaf samples, the total concentrations of N, P, K, Ca, Mg, S, B, Zn, Cu, Mn, and Fe are usually determined. Nevertheless, there are cases where authors used just five nutrients (N, P, K, Ca, and Mg) [4, 11] or more, as the case of Anjaneyulu and Raghupathi [12], who used ten nutrients (N, P, K, Ca, Mg, S, Fe, Mn, Zn, and Cu).

The leaf nutrient concentration must follow normal distribution. That is why, it is necessary to transform the data of nutrient concentration from leaf tissue. However, with CND method, the calculation of row-centered log ratio corrects the non-normal distribution. As observed by Serra et al. [6], the multinutrient variable (Zi) showed 100% normal distribution in data of leaf nutrient concentration of cotton.

The size of the database (leaf nutrient concentration and yield) is not well defined in the literature. Serra et al. [6, 7] used 65 sampling from commercial plots to develop CND norms in the Western region of Bahia State, in Brazil, whereas Serra et al. [13] used 108 sampling in the database for DRIS. Khiari et al. [11] used 240 observations of commercial plots of sweet corn and five nutrient determinations, and Parent et al. [5] collected 1117 samples in potato crop to establish the database. More than quantity, the quality of the database might be more far reaching to obtain profitable CND norms. In fact, the essential to have an effective database is to acquire plots with high-yield and healthy leaves with no damage to make possible the development of the CND or DRIS norms.

3.2. Selecting the high-yielding subpopulation to develop CND norms

In order to develop the CND norms, it is necessary to define the high-yielding subpopulation into the database. The database might be divided into two subpopulations, using the mean + 0.5 standard deviation (SD) as criterion to separate the populations into a high-yielding group and low-yielding group [6].

Serra et al. [13] published a paper about the criterion of population selection for cotton crops to separate low-yielding from high-yielding population to develop DRIS norms; nevertheless, it is possible to apply these criteria for CND norms. In this criterion, they used the mean and standard deviation of the yield in the database to separate the high-yielding from low-yielding subpopulation.


Figure 1.

Graphic representation of the observed frequency of yields. The criteria for establishing the high-yield subpopulations were plots with (A) yield above average (4380 kg ha−1), (B) yield above average + 2/3 standard deviation (above 4650 kg ha−1), (C) yield above average + 1 standard deviation (above 4785 kg ha−1), and (D) yield above average + 4/3 standard deviation (above 4920 kg ha−1).

As shown in Figure 1, Serra et al. [13] used as criteria for determining the high-cotton-yield subpopulation plots with yields above the population mean (Figure 1A), yield above mean + 2/3 standard deviation (SD) (Figure 1B), yield above mean + 1 SD (Figure 1C), and yield above mean + 4/3 SD (Figure 1D). They concluded that the increasing of rigor to select the high-yielding subpopulation increases the capacity to discriminate low yielding from high-yielding subpopulation.

For other crops as well as for cotton crop, the selection of the high-yielding reference subpopulation may be carried across multiple ratios using a cumulative variance function fit to cubic [11] or Boltzmann [14] equations.

3.3. Calculation of the CND norms

The nutrient leaf tissue composition is defined by the simplex S where the sum of all components is constrained to 100%, which forms a d-dimensional nutrient arrangement, i.e., simplex (Sd) made of d + 1 nutrient proportions including d nutrients and a filling value (Rd) [4]:


where 100 is the dry matter concentration (%); N, P, K, etc., are nutrient proportions (%); and Rd is the filling value computed as follows:


The geometric mean (G) is computed as follows [4]:


After calculation of the geometric mean (G), the new expression for the multinutrient is log-transformed to generate a row-centered log ratio as follows [4]:


The sum of all row-centered log ratios must be equal to zero (VN+VP+VK++VRd=0. CND norms are the means and standard deviations (SD) of row-centered log ratios of high-yielding subpopulation from the database.

Following the procedure above, Serra et al. [6] developed the CND norms for cotton crop in the Western region of Bahia State, Brazil (Table 1).

Row-centered log ratiosMeanSD

Table 1.

Compositional nutrient diagnosis (CND) norms [mean and standard deviation (SD) of row-centered log ratio]1.

[i] - 1CND norms based on the high-yield subpopulation (>4250 kg ha−1) of cotton crop in the Western region of Bahia state in Brazil [6].

4. CND index and interpretation

4.1. Procedure to calculate CND index in cotton crop

As reported in Section 3.3, from the leaf nutrient concentration of the high-yielding subpopulation (>4250 kg ha−1), CND norms were defined (Table 1), which are the means, and standard deviations of row-centered log ratios of the nutrients concentration in leaf tissue, denoted as ´VN+´VP+´VK++´VRd and ´SDN, ´SDP, ´SDK,.... ´SDRd, respectively, were then calculated.

The CND indices, denoted as IN, IP, IK,…, IRd, were calculated from the row-centered log ratios as follows [4]:


The nutrient imbalance index (NII) of a diagnosed specimen, which is its CND r2, followed the recommendation of Parent and Dafir [4]:


The interpretation of CND is based on CND index. CND index is defined as the distance of a given nutrient Xi from its geometric mean [4], which is relative to the distance of the same nutrient from the geometric mean of the target population (reference population with high yield).

From this point of view above, it is expected that when CND index is closer to zero, Xi nutrient is less imbalanced than others in analysis. Serra et al. [6, 7] observed CND index close to zero showed higher nutritional balance. This may be related to nutrient leaf concentration and result in high relationship in cotton crop.

4.2. Interpretation of CND index by nutrient application potential response (NAPR)

Wadt [15] developed the nutrient application potential response (NAPR), which was originally used to interpret the DRIS index. However, the use of this method was extended to interpret CND index for cotton crop [6, 7]. This method proved to be efficient to accomplish the interpretation of CND index. Besides the interpretation of CND index in cotton crop, the methodology has already been applied in other crops to interpret the DRIS index, such as Eucalyptus grandis [16], coffee [17], and “cupuaçu” trees [18].

The NAPR consists in defining five groups of crop response (Table 2). The use of this methodology enables to obtain the most responsive nutrient in the case of fertilizer recommendation; however, this recommendation might be carried with the soil physical and chemical analysis to guide the decisions.

Nutritional status Criteria Type of nutrient application potential response1
Deficiency IA < 0, |IA| > NIIa2, and IA is the index of lower value Positive, with higher probability (p) 
Deficiency prone IA < 0 and |IA| > NIIa Positive, with low probability (pz) 
Sufficient |IA| ≤ NIIa Null (z) 
Excess prone IA > 0 and |IA| > NIIa Negative, with a low probability (nz) 
Excess IA > 0, |IA| > NIIa, and IA is the index of higher value Negative, with a higher probability (n) 

Table 2.

Criteria to interpret the CND index by nutrient application potential response (NAPR).

[i] - 1The NAPR was calculated based on Wadt [17].

[ii] - NIIa = nutritional imbalance index average. IA = CND index for nutrient A.

Adapted with permission from Wadt [17].

The practical interpretation for the nutritional status of the plant in Table 2 is that “deficiency” means this nutrient shows high probability to positive response if applied in soil, which expects higher cotton yields. On the other hand, “excess” means the application of the nutrient in soil may result in luxury consumption and no response ought to be shown in terms of yield or food quality [15].

NutrientsMethod of diagnosisNutrient application potential response (NAPR)1 (%)

Table 3.

Percentage of nutrient application potential response (NAPR) in cotton crop in Western region of Bahia State, Brazil [6].

[i] - 1(p) positive, with higher probability; (pz) positive, with low probability; (z) null; (nz) negative, with a low probability; (n) negative, with a higher probability.

The use of CND index interpretation for cotton crop can be a useful tool in association with soil analysis and experience of the agronomist to recommend fertilizer for cotton crop. Furthermore, the absence of nutritional diagnosis can result in excess of fertilizer as the case of P fertilizer observed by Serra et al. [6] in cotton crop in Western region of Bahia State, Brazil. In some specific plots where 180 kg P2O5 ha−1 was applied, Serra et al. [6] noted that by using NAPR the P In some specific leaf tissue was higher than plant requirement and it could be reduced by the P fertilizer doses (Table 3).

5. Nutrient interactions through principal component analysis on CND multivariable

In order to calculate the principal components analysis (PCA) of the row-centered log ratio in the high-yielding subpopulation, it was used the database (leaf nutrient content and crop yield) from Serra [19]. As reported by Parent and Dafir [4], the row-centered log ratios are compatible with PCA; based on it, this section aims at observing some interaction among the nutrients in PCA.

Row-centered log ratioHigh-yield subpopulation (n = 40)
VN 0.6288 0.4875 −0.0093
VCa−0.1868 0.7223 −0.0605
VMg0.24010.0017 0.6842
VS0.4352−0.2937 0.5691
VFe 0.8159 0.4529−0.0919
VZn 0.8078 0.4037−0.0469
VMn 0.3879 0.6686 −0.1689
VB 0.2892 −0.08710.0494
Explained variance (%)28.2518.7213.81
Accumulated variance (%)28.2546.9760.78
Selection criterion (SC)0.28360.34840.4056

Table 4.

Correlations between the row-centered log ratio and the first three principal components (PC) for the high-yield subpopulation (40 observations).

Values in boldface are the dominant in the PC loadings by setting the level of significance defined according the selection criterion (SC).

Database from Serra [19].

One of the major objectives of the PCA is to reduce the number of interdependent variables into smaller numbers of independent principal component (PC), which are linear combinations of original variables. According to Parent and Dafir [4], CND is the multivariate expansion of CVA and DRIS and fully compatible to PCA.

Based on the PCA, it was obtained three PCs with eigenvalues above 1 and accumulated variance equal 60.78% (Table 4). The definition of the loading significance was defined according to Ovalles and Collins [20]. PC loadings or correlations between row-centered log ratios and the first three principal components must have values greater than the selection criterion (SC) for significant acceptation. SC proposed by Ovalles and Collins [20] was defined as follows:


The three PCs combined explained 60.78 % of the total variance (Table 4). The first principal component was positively correlated with VN, VFe, VZn, VMn, and VB and negatively correlated to VP, S, and VCu (Table 4). The second principal component was positively correlated with VN, VCa, and VMn and negatively correlated with VP, VFe, and VZn (Table 4). Finally, the third principal component showed positively correlated with VMg and VS and negatively correlated with VK (Table 4).


Figure 2.

Distribution of principal component analysis (PCA) for row-centered log ratios. PC = principal components. Database from Serra [19].

The structures obtained by PC-1 (N+P-S-Fe+Zn+Cu-Mn+B+), PC-2 (N+P-Ca+Fe-Zn-Mn+), and PC-3 (K-Mg+S+) are supported by some rules of plant nutrition relative to interaction between nutrients, but other structures should be studied carefully for a better understanding (Figure 2).

In PC-1, it is feasible to observe positive interaction among N, Fe, and Zn; this result may be explained as reported by Marschner [21] that Fe, Zn, and protein are very highly correlated, which infer that protein is a sink for Fe and Zn. In PC-2, the inverse relation between P and Ca may be explained because of dilution and accumulation effects with plant age. The PC-3 showed the opposite direction for K and Mg inferring the antagonism effects between these nutrients (Figure 2).

N and S showed to be highly correlated with the first PC; however, S was negatively correlated with N (Figure 2A). It is important to have N availability in soil in balance with S, since this nutrient is a component of amino acids and proteins. Cysteine is the first stable compost of the assimilatory reduction of S that is the precursor of all compounds, in which S makes part in the plants [21]. Thus, S deficiency compromised the formation of amino acids. Furthermore, the N and S show essential functions in the activity of chlorophyll molecule. The S deficiency may deplete the action of chlorophyll and consequently reduction of photosynthesis [22].

The second PC was strongly influenced by positive and negative eigenvectors VN, VCa, and VMn and VP, VFe, and VZn, respectively (Figure 2C). In positive and negative quadrant, the nutrients are grouped according to the relationship among them in the plant metabolism, as example, the Zn affects the metabolism of Fe [23], or interactions among nutrients in the rhizosphere. VMg (positive) and VK (negative) affected the third PC, which showed opposite directions between these variables inferring the excessive application of K source in soil, inhibit the uptake of Mg by plant roots [21].

6. Conclusion

Being an effective diagnosis tool for the nutritional status of cotton, CND may improve the diagnosis with the application of multivariate analysis. The development of CND norms including means and standard deviation of row-centered log ratios in high-yielding subpopulation can be accomplished based on field trials or experimental database. This information in database is comprised at least by the concentration of leaf nutrients and cotton yield in well-defined plots. CND, due to its interactions with multi-nutrients, can improve the nutritional diagnosis in cotton crop better than traditional methods as CVA and SR. In literature, there are not many results published about CND in cotton crop, but researches should be encouraged worldwide to use this method because of CND importance, to know the plant nutritional requirement, and to direct the decision on fertilizer management. Moreover, the row-centered log ratio is fully compatible with principal component analysis (PCA) that can expand the results of nutrient interaction in leaf tissue of cotton crops.


1 - Bates TE. Factors affecting critical nutrient concentrations in plants and their evaluation: a review. Soil Science. 1971; 112: 116-130.
2 - Beaufils, E.R. Diagnosis and recommendation integrated system (DRIS). A general scheme of experimentation and calibration based on principles developed from research in plant nutrition. 1st ed. South Africa: University of Natal, Pietermaritzburg; 1973. 132 p.
3 - Parent LE, Natale W, Ziadi N. Compositional nutrient diagnosis of corn using the Mahalanobis distance as nutrient imbalance index. Canadian Journal of Soil Science. 2009; 89: 383-390. DOI: 10.4141/cjss08050
4 - Parent LE, Dafir MA. Theoretical concept of Compositional Nutrient Diagnosis. Journal of American Society for Horticultural Science. 1992; 117: 239-242. DOI: 117/2/239.full.pdf
5 - Parent LE, Cambouris AN, Muhawenimana A. Multivariate diagnosis of nutrient imbalance in potato crops. Soil Science Society of American Journal. 1994; 58: 1432-1438. DOI: 10.2136/sssaj1994.03615995005800050022x
6 - Serra AP, Marchetti ME, Vitorino ACT, Novelino JO, Camacho MA. Development of DRIS and CND norms and evaluation of nutritional status of cotton. Revista Brasileira de Ciência do Solo. 2010; 34: 97–104. DOI: 10.1590/ 7 S0100–06832010000100010.
7 - Serra AP, Marchetti ME, Vitorino ACT, Novelino JO, Camacho MA. Determination of normal nutrient ranges for cotton by the ChM, CND and DRIS methods. Revista Brasileira de Ciência do Solo. 2010; 34: 105–113. DOI: 10.1590/S0100–06832010000100011
8 - Hallmark WB, Walworth JL, Sumner ME, Mooy CJ, Pesek J, Shao KP. Separating limiting from non-limiting nutrients. Journal of Plant Nutrition.1987; 10: 1381–1390. DOI: 10.1080/01904168709363670
9 - Malavolta E. Mineral nutrition of plants manual. 1st ed. São Paulo: Agronômica Ceres: 2006. 638p.
10 - Marur CJ, Ruano O. A reference system for determination of developmental stages of upland cotton. Revista Brasileira de Oleaginosas e Fibrosas. 2001; 5: 313-317
11 - Khiari L, Parent L, Tremblay N. Critical compositional nutrient indexes for sweet corn at early growth stage. Agronomy Journal. 2001; 93: 809–814. DOI: 10.2134/agronj2001.934809x
12 - Anjaneyulu K, Raghupathi HB. CND and PCA approaches for multivariate diagnosis of nutrient imbalance in papaya (Carica papaya L.). Acta Horticulture. 2010; 851: 363-368. DOI: 10.17660/ActaHortic.2010.851.55
13 - Serra AP, Marchetti ME, Rojas EP, Morais HSD, Conrad VDA, Guimarães FCN. Establishing DRIS norms for cotton with different selection criteria for the reference population. Pesquisa Agropecuária Brasileira. 2013; 48: 1472- 1480. DOI: 10.1590/S0100-204X2013001100008
14 - Hernandez-Caraballo EA, Rodrıguez-Rodrıguez O, Rodrıguez-Perez V. Evaluation of the Boltzmann equation as an alternative model in the selection of the high-yield subsample within the framework of the compositional nutrient diagnosis system. Environmental and Experimental Botany. 2008; 64: 225–231. DOI: 10.1016/j.envexpbot.2008.05.010
15 - Wadt PGS. The method of mathematic chance and diagnosis and recommendation integrated system (DRIS) in nutritional evaluation of Eucalyptus planting [Thesis].Viçosa: Universidade Federal de Viçosa; 1996.
16 - Silva GGC, Neves JCL, Alvarez VH, Leite, FP. Evaluation of the universality of DRIS, M-DRIS, and CND norms. Revista Brasileira de Ciência do Solo. 2005; 29: 755–761. DOI:10.1590/S0100-06832005000500011
17 - Wadt PGS. Relationships between soil class and nutritional status of coffee crops. Revista Brasileira de Ciência do Solo. 2005; 29: 227-234. DOI: 10.1590/S0100-06832005000200008
18 - Dias JRM, Wadt PGS, Perez DV, Silva LM, Cleigiane OL. DRIS formulas for evaluation of nutritional status of cupuaçu trees. Revista Brasileira de Ciência do Solo. 2011; 35: 2083–2091. DOI: 10.1590/S0100-06832011000600023
19 - Serra AP. Evaluation of nutritional status and determination of optimal ranges for cotton [Thesis]. Dourados: Universidade Federal da Grande Dourados; 2008.
20 - Ovalles FA, Collins ME. Variability of northwest Florida soils by principal component analysis. Soil Science Society of America Journal. 1988; 52: 430–1435. DOI: 10.2136/sssaj1988.03615995005200050042x
21 - Marschner H. Mineral nutrition of higher plants. 3rd ed. London: Academic Press; 2012. 651 p.
22 - Nikiforova VJ, Kopka J, Tolstikov V, Fiehn O, Hopkins L, Hawkesford MJ, Hesse H, Hoefgen R. Systems rebalancing of metabolism in response to sulphur deprivation, as revealed by metabolome analysis of arabidopsis plants. Plant Physiology. 2005; 138: 304-318. DOI: 10.1104/pp.104.053793
23 - Rout GR, Sahoo S. Role of iron in plant growth and metabolism. Reviews in Agricultural Science. 2015; 3: 1-24. DOI: 10.7831/ras.3.1