Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 22 January 2016
Sec. Plant Nutrition

Fruiting Branch K+ Level Affects Cotton Fiber Elongation Through Osmoregulation

\r\nJiashuo YangJiashuo YangWei HuWei HuWenqing ZhaoWenqing ZhaoBinglin ChenBinglin ChenYouhua WangYouhua WangZhiguo Zhou*Zhiguo Zhou*Yali Meng*\r\nYali Meng*
  • Key Laboratory of Crop Physiology and Ecology, Department of Agronomy, College of Agriculture, Nanjing Agricultural University, Nanjing, China

Potassium (K) deficiency in cotton plants results in reduced fiber length. As one of the primary osmotica, K+ contributes to an increase in cell turgor pressure during fiber elongation. Therefore, it is hypothesized that fiber length is affected by K deficiency through an osmotic pathway, so in 2012 and 2013, an experiment was conducted to test this hypothesis by imposing three potassium supply regimes (0, 125, 250 kg K ha-1) on a low-K-sensitive cultivar, Siza 3, and a low-K-tolerant cultivar, Simian 3. We found that fibers were longer in the later season bolls than in the earlier ones in cotton plants grown under normal growth conditions, but later season bolls showed a greater sensitivity to low-K stress, especially the low-K sensitive genotype. We also found that the maximum velocity of fibre elongation (Vmax) is the parameter that best reflects the change in fiber elongation under K deficiency. This parameter mostly depends on cell turgor, so the content of the osmotically active solutes was analyzed accordingly. Statistical analysis showed that K+ was the major osmotic factor affecting fiber length, and malate was likely facilitating K+ accumulation into fibers, which enabled the low-K-tolerant genotype to cope with low-K stress. Moreover, the low-K-tolerant genotype tended to have greater K+ absorptive capacities in the upper fruiting branches. Based on our findings, we suggest a fertilization scheme for Gossypium hirsutum that adds extra potash fertilizer or distributes it during the development of late season bolls to mitigate K deficiency in the second half of the growth season and to enhance fiber length in late season bolls.

Introduction

Potassium (K) is a critical macronutrient for cotton production as it is vital to enzyme activation, electrochemical gradient maintenance, cell turgor generation, and osmotic regulation (Marschner, 1995; Amtmann and Blatt, 2008; Zhang et al., 2009; Römheld and Kirkby, 2010). Cotton (Gossypium hirsutum L.) is more sensitive to K deficiency compared to other crops (Cope, 1981; Cassman et al., 1990) because its developing bolls act as strong K sinks (Boquet and Moser, 2003). As a consequence of the constant decrease in soil available K caused by improved yields and insufficient K application, an increasing number of K deficiency events have been affecting cotton production worldwide (Oosterhuis et al., 1994). This issue has become a limiting factor for cotton production in China (Zhang et al., 2007).

Field tests have shown that soil K deficiency not only reduces K concentrations in cotton leaves (Reddy and Zhao, 2005), bolls (Wright, 1999), and fibers (Cassman et al., 1990), but it also affects a series of fiber quality parameters (Cassman et al., 1990; Pettigrew et al., 1996), especially fiber length (Bauer et al., 1998). In vitro research has demonstrated a specific K requirement for cotton fiber elongation (Dhindsa et al., 1975), and because K+ is one of the primary osmotically active solutes in fiber cells that powers elongation (Dhindsa et al., 1975; Ruan, 2005), fiber length is very likely to be affected by K deficiency through an osmoregulatory pathway. However, no direct evidence currently exists to prove this hypothesis. Additionally, Cassman et al. (1990) indicated that fiber length was positively related to fiber K+ concentration at maturity, whereas the K+ in fiber also plays a major role before 25 days postanthesis (DPA; Dhindsa et al., 1975). Therefore, it might be useful to uncover the quantitative relationship between fiber length and fiber K+ level during the rapid elongation stage.

Although it usually takes 40–60 days for cotton fiber to reach maturity, its final length is attained within the first 3 weeks (Ruan, 2005). Recent research has shown that adjacent elongating fibers are joined together by the adhesion of cotton fibre middle lamella (CFML) into tissue-like bundles (Avci et al., 2013), while the fiber cell elongation process follows a single-cell expansion mode (Ruan et al., 2001; Ruan, 2005; Qin and Zhu, 2011). Two inner factors have been found that can regulate fiber elongation (Braden and Smith, 2004; Zhao et al., 2012): (1) the force of cell turgor against cell walls to achieve cell expansion (Ruan et al., 2001; Fricke and Peters, 2002) and (2) the duration of cell turgor (Ruan et al., 2001; Ruan, 2005). Cell turgor is achieved through the influx of water, which is driven by osmotically active solutes (Cosgrove, 1997). Soluble sugars, K+ and malate are the three main osmotically active solutes in elongating fibers, which together account for approximately 80% of the total osmotic potential (Dhindsa et al., 1975; Ruan et al., 1997). To achieve elongation, highly expressed endo-1,4-glucanase, expansin proteins, etc., contribute to unlocking the network of cell wall polysaccharides, which permits turgor to drive fiber cell elongation (Shimizu et al., 1997; Ruan et al., 2001; Harmer et al., 2002). However, the down-regulation of these proteins during the elongation phase would hinder fiber cell lengthening (Frederic et al., 1998; Cosgrove, 2000). Secondary cell wall formation would also regulate fiber cell elongation, which usually starts at 16–20 DPA (days postanthesis) and, compared to the early stage, increases the cellulose synthesis rate by over 100-fold, which ultimately makes the cell wall too rigid to elongate (Ruan, 2007; Wang and Ruan, 2010). Additionally, plasmodesmata (PD) have been found to be gateways that control the generation of fiber cell turgor (Ruan et al., 2001). Although all of the factors mentioned above can regulate the fiber elongation process, only the accumulation of osmotic solutes can drive this process (Delmer, 1999; Martin et al., 2001; Ruan, 2005). Therefore, it is logical to study the mechanism by which K+ level affects fiber length.

Cotton bolls grown on the upper fruiting branches (late season bolls) are more prone to symptoms of K deficiency (Read et al., 2006). Previous research has shown that soil K deficiency reduced fiber length by 5.0–7.6% in late season bolls, but this reduction was only 0–1.8% in early season bolls (Read et al., 2006). More information is needed to understand fiber length on a whole plant level. In consideration of the issues mentioned above and some other key points, a 2-year test has been undertaken on a couple of severely K-deficient fields. This study aimed to determine: (a) the characteristics of the changes that occur during fiber elongation under K deficiency; (b) the effect of plant K+ status on cotton fiber elongation and the relationship between [K+] in tissues and fiber elongation properties; (c) the reason why late season bolls are sensitive to K deficiency; (d) the mechanism underlying genotypic differences in the capacity for low-K tolerance.

Materials and Methods

Cultivar and Site Description, Experimental Design and Crop Management

Two cotton (G. hirsutum L.) cultivars, Simian 3 and Siza 3, have been identified as low-K-tolerant and low-K-sensitive genotypes, respectively (Yang et al., 2014, doi: 1002-7807(2014)04-0301-09, Chinese version with English abstract, see Supplementary Figure S1). Cotton was planted in a purpose-built test field that had been potassium (K) depleted for two crop seasons, which is located at the Pailou site (yellow–brown loam) of Nanjing Agricultural University, Jiangsu Province, China. Cotton was sown on 25th Apr and transplanted on 22 May in both 2012 and 2013. The soil soluble K+ concentrations (0–40 cm) before planting in 2012 and 2013 were 92 and 86 mg kg-1, respectively. To avoid the effects of the fertilizer residues from the first year on the experiment in the second year, we separated the test field into two parts for each planting.

Three K fertilizer application levels of 0, 125, 250 kg K ha-1 per year were arranged in a randomized complete-block design with three replicates. No K applied combined with K-deficient soil created a severely K-deficient condition while the 125 and 250 kg K ha-1 treatments induced mild K deficiency and K sufficiency, respectively (Gormus, 2002; Pettigrew, 2008). Individual treatment plots were 6.6 m × 13 m size with 16 rows each, and K fertilizer (potassium sulfate, K content = 42%) was accurately hole-applied at the beginning of flowering season in each year. Nitrogen fertilization (240 kg N ha-1 per year, urea), phosphate fertilization (52 kg P ha-1 per year, calcium superphosphate), irrigation, pest control, etc., followed recommended practices. There were no indications that N supply, water stress or pest pressure negatively influenced cotton growth during the 2 years of the experiment.

Sampling Procedure and Measurement of Fiber Length

The first fruit-bearing sympodial branch was defined as the first fruiting branch (FB1). Generally, based on their growth habits and the climate in Nanjing, China, these two cultivars would have 14–16 fruiting branches before artificial topping. Therefore, the 3rd, 8th, and 13th fruiting branches of the cotton plant were tagged as representatives of the lower, middle, and upper fruiting branches, respectively. Each FB sample was cut from three replicated plants to remove the bolls, leaves, flowers, and buds before oven drying at 105°C for 0.5 h and then at 70°C to a constant weight before grinding into powder for use. All soil samples (0–40 cm) were collected from nine separate locations from each plot and mixed together and naturally air-dried to a constant weight followed by grinding into powder for use. Both the fruiting branches and soil samples were collected on 25 July, 5 August, 15 August, 25 August, and 5 September, which corresponded to 63, 74, 84, 94, and 105 days after transplantation (22 May) in both years. Cotton boll age was determined by tagging the fully opened white flower (onset of anthesis, 0 DPA).

For the measurement of immature fiber length at 5, 10, 15, 20, 24, 31, 38, and 45 DPA, bolls grown on the 3rd, 8th, and 13th fruiting branches were collected (plus an additional 17-DPA sample for the eight FB). Cotton bolls were picked and stored on ice for no more than 1 h before dissection, and bolls younger than 30 DPA were boiled with their locules in 0.1% HCl to separate the seeds and relax the fibers (Schubert et al., 1973). Thereafter, fiber length was measured by the water stream method (Thaker et al., 1989). Each seed was carefully placed on the convex surface of a watch-glass, and the fibers on seed were then streamed out with a jet of water. Fiber length was measured with a digital vernier caliper, which could be accurate to 0.01 mm, under a dissecting microscope. Measurement began from the attachment point on the epidermis of the ovule to the edge where most fibers terminate. Ten seeds were randomly selected from three bolls as replicates for the measurements.

A photoelectric stapler (Y-146 Photoelectric Stapler, Taicang Electron Apparatus Co., Ltd., Taicang, Jiangsu, China) was used to measure the length of fibers older than 30 DPA. Prior to measurement, these fiber samples were preprocessed by an oven dryer at 60°C for 0.5 h, then at 40°C for 2 days, and then placed into a standard testing room with constant temperature (20 ± 2°C) and humidity (65 ± 2%) for 48 h. Fiber samples from 30 DPA were measured by both methods to correct for differences between the two systems.

K+ Concentration Analyses

Fiber samples younger than 30 DPA were divided into those used for the fiber length measurements and the remainder, which were oven-dried through the same fiber-drying processes described above. Both the fiber and FB samples were digested by H2SO4-H2O2 (Samal et al., 2010), and soil soluble K+ was extracted by placing 2 g prepared soil into 20 ml 1 M NH4OAc at pH 7.0 and shaking the flask at 200 rpm for 5 min followed by a filtration process (Warncke and Brown, 1998). The K+ concentrations in the sample solution were measured by an atomic absorption spectrophotometer (TAS-986 Atomic Absorption Spectrophotometer, Beijing PERSEE Co., Ltd., Beijing, China), and [K+]fiber, [K+]FB and [K+]soil were all calculated based on the dry weight.

Malate and Soluble Sugar Concentration Analyses

The extraction method combined with the measurement of fiber malate content was based on Famiani et al. (2000), and the soluble sugar determination method was based on Jones et al. (1977) and Famiani et al. (2000). In summary, 50 mg of a dry fiber sample were frozen by liquid nitrogen and ground into powder and then placed in 1.5 ml of 80% ethanol containing 100 mM Hepes-KOH (pH 7.1) and 20 mM MgCl2 for 1 h at 80°C for extraction. The extract was centrifuged at 12 000 g for 5 min after cooling. The supernatant was recovered, and 150 μl of charcoal suspension (100 mg ml-1) was added with stirring, after which another centrifugation with the same settings was undertaken. The last supernatant was the extract.

Glucose, fructose, and sucrose contents were measured using the enzyme-coupled method (Jones et al., 1977; Famiani et al., 2000). First, 200 μl of extract was added into the 1 ml assay buffer containing 100 mM Hepes-KOH (pH 7.0), 5 mM MgCl2, 0.5 mM DTT, 0.02% (w/v) BSA, 100 mM ATP, and 40 mM NAD+. Then, glucose was measured by initiating the reaction by adding 3 U (in 100 μl) of hexokinase (from Saccharomyces cerevisiae, Sigma Co.) and 1 U (in 100 μl) of glucose-6-phosphate dehydrogenase (from Leuconostoc mesenteroides, Sigma Co.) into the assay mixture. After 30 min at 30°C, the absorption was measured at 340 nm, and the standard curve was made by replacing the sample extraction with a series of glucose solutions (prepared in extraction buffer). Fructose and sucrose were sequentially analyzed after glucose following the addition of 1 U (in 100 μl) of phosphoglucose isomerase (from S. cerevisiae, Sigma Co.) and 100 U (200 μl) of invertase (from S. cerevisiae, Sigma Co.), respectively, to the assay mixture. The same detection method was used for fructose and sucrose determination as for glucose. The measurement principle was based on Jones et al. (1977).

Malate (Famiani et al., 2000) was measured by adding 200 μl of extract into 1 ml of assay buffer containing 50 mM 2-amino-2-methylpropanol, 40 mM glutamate (pH 9.9), and 1 mM NAD+. The reaction was initiated by adding 10 U (in 100 μl) of glutamate-oxaloacetate transaminase (GOT, from pig heart, Sigma Co.) and 0.7 U (in 100 μl) of L-malate dehydrogenase (L-MDH, from pig heart, Sigma Co.) into the assay mixture. After 1 h at 20–25°C, the absorption was measured at 340 nm, and the standard curve was built by replacing the sample extract with a series of malate solutions (prepared in extraction buffer).

Statistical Analyses

ANOVA was performed using SPSS 20 (IBM). Figures 19 were drawn with Origin 9.0 (OriginLab), and Figure 10 was drawn with PowerPoint (Microsoft).

FIGURE 1
www.frontiersin.org

FIGURE 1. K deficiency symptoms in cotton plant (A), cotton fibers (B), and growing season temperatures and rainfall (C). Fibers shown in (B) were ginned from the 8th fruiting branch (FB8) of cv. Siza 3 and were tidied by a pair of special combs. Bar = 10 mm. K250, K125, and K0 represent 250, 125, and 0 kg K ha-1, respectively.

FIGURE 2
www.frontiersin.org

FIGURE 2. Changes in FB K+ concentration ([K+]FB, A) and soil soluble K+ concentration ([K+]soil, B) during cotton growth. Each data point in Figure 2 was determined from a sample mixture of three replicate plants (A) or nine points in one plot (B).

FIGURE 3
www.frontiersin.org

FIGURE 3. Changes in final fiber length (A) and its relationship with fiber K+ concentration ([K+]fiber) at 10 DPA (B) and FB K+ concentration ([K+]FB) at 10–15 DPA (C). CS (coefficients of stress) was calculated by the formula [Lenfibre(K250)-Lenfibre(K0))/Lenfibre(K250) × 100%, in which “Lenfibre represents fiber length. The [K+]fiber data presented above were collected at 10 DPA as [K+]fiber reached its peak value at approximately 10 DPA (see Figure 4). The [K+]FB data were collected on 25 July for FB3, on 5 August for FB8 and on 15 August for FB13, which correspond to the 10–15 DPA for the bolls grown on each FB (see arrows in Figure 2A). Bars represent ±SE (10 replicates).

FIGURE 4
www.frontiersin.org

FIGURE 4. Dynamic changes in fiber length from 5 DPA to boll opening, which grown on the 13th (A), 8th (B) and 13th (C) fruiting branches. Each data point represents the mean value from 10 seeds (ovules) with three measurements per seed. Bars represent ±SE (30 replicates).

FIGURE 5
www.frontiersin.org

FIGURE 5. Coefficients of variation (CV) of “the maximum velocity of fibre elongation” (Vmax) and “the fibre rapid elongation duration” (FRED) caused by the varying FB K+ concentrations ([K+]FB) or “FB/boll period temperature” (FB/temp.) ratios. The CV was calculated by the formula SD/Mean × 100%, and the “FB/temp.” factor represents the variation in FB that appeared with variation in boll period temperature. The mean CV value (short black bar) caused by [K+]FB (the left side) was calculated by combining the cultivars, test years and fruiting branches; the value caused by “FB/Temp.” (the right side) was calculated by combining the cultivars and test years and only accounted for K-sufficient cotton (K250).

FIGURE 6
www.frontiersin.org

FIGURE 6. Concentrations of main osmotically active solutes in cotton fiber (mg g-1 DW) from 5 to 24 DPA (A) and the results of the peak value analysis (10 DPA, B). Each data point is an average of three biological replicates. For (A), we combined all of the data from the three K treatments, three branches, two cultivars and two test years into one color with the aim of detecting their peak values. Peak values (B) were observed from 10 DPA because both the peak values and Vmax (Figure 4) consistently emerged at approximately 10 DPA. The reason for choosing soluble sugar concentration at 10 DPA as a peak value instead of 5 DPA was that the fibers did not elongate rapidly at 5 DPA. Even if there was a higher value at 5 DPA, it had no significant meaning for the elongation process. Additionally, the relatively higher value at 5 DPA might be caused by the difficulty stripping fiber at 5 DPA, so that some fiber maybe have ovule tissue residue. The three points aligned in a row (B) represent fibers on the same FB, from the same cultivar, in the same year but subject to different K treatments.

FIGURE 7
www.frontiersin.org

FIGURE 7. Coefficients of stress (CSs) of the concentrations of the main osmotically active solutes in fibers at 10 DPA when suffering from K deficiency. CS was calculated by the formula [Ci(K250) - Ci(K0))/Ci(K250) × 100%, in which “C” represents the concentration at 10 DPA, and i could be fiber K+, malate or soluble sugar. CS was calculated regardless of FB, which means it was the mean value of the data from three fruiting branches.

FIGURE 8
www.frontiersin.org

FIGURE 8. Correlations between “Vmax and fiber K+ concentration [[K+]fiber] at 10 DPA and the repeatability between the two test years (insert). The linked arrows indicate the large difference between the two cultivars in their K+ accumulation ability under K0.

FIGURE 9
www.frontiersin.org

FIGURE 9. Correlations between fiber K+ concentration ([K+]fiber) at 10 DPA and FB K+ concentration ([K+]FB) at 10–15 DPA.

FIGURE 10
www.frontiersin.org

FIGURE 10. A model of how FB K+ concentration [[K+]FB] regulates fiber elongation by osmoregulation, and a proposed K fertilization scheme based on it. FRED, fiber rapid elongation duration; Vmax, the maximum velocity of fiber elongation.

The dynamic change in fiber length was regressed by the following logistic function (1) (Naithani et al., 1982; Zhao et al., 2012), in which Len represents fiber length. Lenm is the theoretical maximal fibre length, which is a parameter in this model along with a and b.

Len=Lenm1+a×eb×DPA(1)

Eigenvalues, Vmax (the maximum velocity of fiber elongation) and FRED (fiber rapid elongation duration), were derived from formulae (2) and (3) (Zhao et al., 2012), in which DPAREI represents the initiation DPA of rapid elongation, and DPARET represents the termination DPA.

Vmax=b×Lenm4(2)
FRED=DPAREIDPARET=1b×In23a1b×In2+3a=1b×In232+32.6339×1b(3)

CV = Standard deviation (SD)/Mean × 100%.

CS = [Lenfibre(K2500)-Lenfibre(K0)]/Lenfibre(K250) × 100%, in which “Lenfibre represents fiber length. The formula can also be written as [Ci(K250)-Ci(K0)]/Ci(K250) × 100%, in which “Ci represents the concentration of fiber K+ or malate or soluble sugar at 10 DPA.

Results

K+ Levels in Soil, Cotton Fruiting Branches, and Fiber Cells

Soil K deficiency induced symptoms in leaves that corresponded with previous studies (Oosterhuis, 1995; Wright, 1999), indicating that our field experiment met the expectation that K deficiency was severe enough to affect the growth of cotton lint (Dong et al., 2004). K-deficient leaves do not produce enough carbohydrate through photosynthesis for the heavy boll load of cotton (Wright, 1999); visible differences in final fiber length were observed among the K treatments (Figure 1B).

In 2013, there were 43 days with an MDT (mean daily temperature) over 30°C, twice the 21 days in 2012 (Figure 1C). The total rainfall from 25th April to 1st November in 2013 was 720.8 mm, which exceeded the 601.2 mm in 2012. Torrential rain occurred on 9th Aug 2012 with a total rainfall of 126.3 mm, and the relatively higher MDT combined with the more regular rainfall in 2013 improved cotton growth over 2012. This was reflected in a slower reduction in [K+]FB (FB K+ concentration) in the control plants (K250 treatment, Figure 2A).

The [K+]FB continued decreasing from the flowering to the boll-opening stage, and the three fruiting branches possessed distinct [K+]FB curves (Figure 2A). Under K deficiency, significantly lower [K+]FB were found in the low-K-sensitive cultivar, Siza 3, on the 13th FB (see the red and green lines in Figure 2A). Arrows on the [K+]FB curve indicate a date within 10–15 DPA of the corresponding bolls (Figure 2A), when fibers reached peak velocity of elongation (Figure 4). In other words, it was the period when [K+]FB most likely affected fiber elongation through osmoregulation (Figure 6A).

Fiber length varied from one FB to another, but bolls on the upper fruiting branches had longer fibers (Figure 3A). Meanwhile, the fiber lengths of the bolls on the upper fruiting branches decreased more when under severe K deficiency (Figure 3A). Generally speaking, cv. Siza 3 fibers were longer than those of cv. Simian 3 under the K-sufficient condition, but they suffered a greater reduction in length when K supply became insufficient, especially on the upper fruiting branches (Figure 3A). By analyzing the correlation between final fiber length and [K+]fiber (10 DPA) or [K+]FB (10–15 DPA), both [K+]fiber and [K+]FB were found to have strong positive effects on final fiber length on all of the sampled fruiting branches (Figures 3B,C). Most of the relationships were linear except those between fiber length and the [K+]FB of the third FB for cv. Simian 3 in 2012 and 2013. However, inconspicuous relationships among fruiting branches could have been present, so further analysis of fiber elongation was needed to determine why [K+]fiber and [K+]FB regulated fiber length so efficiently.

Characteristics of Fiber Elongation Process

The elongation process reached maximum velocity at approximately 10–15 DPA, and fiber length reached a peak value after 24 DPA. These data fit a sigmoidal model very well (Figure 4), and the eigenvalues Vmax (maximum velocity of elongation), FRED (fiber rapid elongation duration), and Lenmax (maximum final fiber length) were calculated from the fitted logistic function (Table 1). In the current study, in which fiber length was directly influenced by the varying [K+]FB, it seems that both Vmax and FRED have been modified (Table 1). Through the CV analysis, we found that Vmax rather than FRED was readily influenced by [K+]FB (left in Figure 5). From the analysis of the impact of FB node on these two eigenvalues (right in Figure 5), FRED was preferentially influenced by the “FB / Temp.” factor instead of [K+]FB. We used “FB/Temp.” to represent the change in FB node and the correspondingly different boll period temperatures. In other words, the variation in fiber length from bolls from different seasons was mainly caused by their different positions on the plant and boll period temperature. Therefore, Vmax was selected as a dominant eigenvalue regulating the fiber elongation process under low-[K+]FB stress, and focusing on Vmax rather than FRED in the following analysis helped us remove the interference from the FB.

TABLE 1
www.frontiersin.org

TABLE 1. Main eigenvalues (Vmax, FRED, Lenmax) of the fiber elongation process logistic function.

Osmotically Active Solutes in Cotton Fiber

Cell turgor is achieved through the influx of water, which is driven by the accumulation of osmotically active solutes including K+, malate and soluble sugar (Ruan et al., 1997). By measuring the changes in their concentrations during rapid fiber elongation (5–24 DPA), we found that their peak values emerged at 10 DPA for K+ and malate and at 5–10 DPA for soluble sugar (Figure 6A). This coincided with the period of Vmax (Figure 4). Moreover, [K+]fiber exhibited greater variation than the other two solutes, and there was a greater change in [K+]fiber in Siza 3 when switching from K sufficiency to deficiency (Figure 6B). These results indicate that K deficiency mainly altered the K+ content in elongating fibers compared with other osmotic solutes, and this pattern was more obvious in the low-K sensitive genotype. More interestingly, the malate content in the fibers of cv. Simian 3 was consistently higher than that in cv. Siza 3 (Figure 6B) regardless of the FB or K treatment in which it occurred. This reflects a possible connection between [K+]fiber and malate, and that malate might play a role in promoting [K+]fiber accumulation in severely K-deficient soil in the low-K-tolerant cultivar, which could be meaningful for evaluating a genotype’s capacity for coping with K deficiency. Additionally, a quantitative comparison between concentrations of different osmotically active solutes using CS (CS; Figure 7) indicated that the effect of K deficiency (Figure 2A) on [K+]fiber was 2.5–3.2 times (in 2012) and 2.0–2.6 times (in 2013) greater than on malate and 2.7–3.4 times (in 2012) and 2.1–2.8 times (in 2013) greater than on soluble sugar.

Relationship Between K+ Level and Fiber Elongation

A series of correlations between Vmax and [K+]fiber (10 DPA) are shown in Figure 8. There were significant linear relationships between Vmax and [K+]fiber on each FB; interestingly, the regression lines of the two cultivars almost overlap (Figure 8). In contrast, these two cultivars also had very different Vmax values under the K0 treatment, especially for the upper fruiting branches (see the linked arrows in Figure 8), and the different amounts of K+ available for fiber elongation between cultivars might be a crucial reason for this (Figure 2A). When all of the regression lines from the three fruiting branches were combined, an approximate exponential curve appeared (Figure 8), in which the upper fruiting branches tended to have lower [K+]fiber while the lower fruiting branches tended to have higher [K+]fiber. Based on these results and a repeatability test between the two test years (insert in Figure 8), we found that Vmax had an exponential relationship with the fiber peak K+ concentration when considering all three fruiting branches. This relationship could be successfully applied to both low-K-sensitive and low-K-tolerant upland cotton genotypes. Because this relationship varied little between the two test years, we combined the data from the 2 years and produced the following function (1), which is a relationship between Vmax and fiber peak K+ concentration (for the current study, it was [K+]fiber(10 DPA)).

Vmax=3.17146.48×EXP(0.4×[K+]fiber(10DPA))R2=0.9621**,n=108(4)

To verify whether [K+]FB could efficiently regulate [K+]fiber, a correlation analysis between [K+]fiber and [K+]FB was performed (Figure 9). Although the intensity of the influence of [K+]FB on [K+]fiber changed from year to year, the difference between the two genotypes remained small (Figure 9). Due to this linear regulation, it seems very likely that [K+]FB regulated [K+]fiber through the K+-osmotic pathways.

Discussion

K+ Plays a Vital Role in Cotton Fiber Growth Through Osmoregulation

A minimum K+ level in cotton tissue is essential for plant growth (Leigh and Wyn Jones, 1984) as well as the development of cotton fiber (Dhindsa et al., 1975). Soluble sugar, malate, and K+ were considered the main active solutes in rapidly elongating fibers (Dhindsa et al., 1975; Ruan et al., 1997), so the K+ status of a cotton plant is likely to regulate fiber elongation through the osmoregulatory pathway (Cassman et al., 1990). Based on our experiment, we verified the role of K+ in the osmoregulation between the FB and elongating fibers (Figure 7). An exponential relationship was found between Vmax, and [K+]fiber (Figure 8), and a linear relationship was found between [K+]fiber and [K+]FB (Figure 9).

Vmax and FRED are important eigenvalues in the fiber elongation process due to their abilities to change fiber length (Zhao et al., 2012), but which one plays the leading role during this process is not clear. Vmax depends on cell turgor (Campas and Mahadevan, 2009), which reflects the force driving fiber elongation (Ruan et al., 2001; Fricke and Peters, 2002); FRED tends to be controlled by both the timing of secondary cell wall synthesis (Ruan, 2007; Wang and Ruan, 2010) and the structural change in the PD between the fiber cell and the epidermal cell (Ruan et al., 2001). In our study, the main effect of plant K+ status on fiber elongation was through the regulation of Vmax (Figure 5). For the purpose of building the physiological connection between fiber elongation and the K+ status of the plant, we focused our analysis on Vmax, and the results indicated that the Vmax could be osmotically adjusted by the plant K+ status (Figures 8 and 9).

Fibers of Late Season Bolls are More Sensitive to K Deficiency

Cotton (G. hirsutum) has indeterminate growth, so in China, it fruits continuously from July to October. According to the growth in body size and the increasing boll amount, cotton bolls grown at different positions along the stem suffer from different nutritional (Boquet and Moser, 2003) and environmental conditions (Davidonis et al., 2004; Zhao et al., 2012). Late season bolls are more prone to K deficiency, which is reflected in K deficiency symptoms (Pettigrew et al., 1996; Read et al., 2006). According to our findings, upper fruiting branches have a greater potential to grow longer fibers, but their fiber length decreases more dramatically under K deficiency (Figure 3A). This might be caused by their poor nutrient competitiveness for K+ in plant (Figure 2A), which promotes K deficiency disease in the upper fruiting branches.

Genotypic Differences in Fiber Elongation in Response to K Deficiency

The efficiency of K utilization varies from species to species, and great genotypic variation in K utilization can exist in a single species (Wang et al., 2012). In our study, we used two cultivars with distinct capacities for low-K tolerance (Yang et al., 2014). Cultivar Siza 3, the low-K-sensitive genotype, indeed had a more favorable fiber length trait than cv. Simian 3 under K sufficiency (Figure 3), but it exhibited a greater decline in fiber length when conditions turned to K deficiency (Figure 3). Therefore, regardless of the intrinsic differences, Siza 3 was more sensitive to K deficiency in terms of fiber length.

A correlation analysis between Vmax and [K+]fiber of the two cultivars was also performed, and the results show that cv. Siza 3 exhibited lower [K+]fiber than cv. Simian 3 under K-deficient conditions (Figure 8). This is the possible reason why cv. Siza 3 had smaller Vmax values compared with cv. Simian 3, especially under the K0 treatment. The genotypic difference in [K+]fiber between two cultivars appears to be due to their differing capacities for absorbing K+ (Figures 2 and 9). Wang et al. (2012) reported similar results in terms of cotton biomass accumulation and photosynthesis activity; in their study, a K-efficient cultivar produced 66.7% more biomass than the K-inefficient cultivar when under K deficiency but similar biomass under K sufficiency.

Interestingly, low-K-tolerant cv. Simian 3 tended to have a higher fiber malate content at 10 DPA compared with low-K-sensitive cv. Siza 3, and the pattern was constant across K treatments, fruiting branches, and years (Figure 6B). We suggest that this might explain the high capacity of cv. Simian 3 for K+ accumulation (Figure 2A). Malate has a collaborative role with K+ in osmoregulation not only in the movement of guard cells (Allaway, 1973) but also cotton fiber elongation (Dhindsa et al., 1975), and the malate accumulating in guard cells can neutralize at least half of the potassium (Raschke and Schnabl, 1978). Hence, the phenomenon of higher malate content in fibers suggests a mechanism through which the low-K-tolerant genotype has stronger capacities to accumulate K+.

Suggestions for Cotton Management Under K-Deficiency

K deficiency in soil (Figure 2B) negatively impacts K+ accumulation in cotton plants. The significantly declining [K+] in the fruiting branches affected the cotton bolls by presenting severely K-deficient growth conditions (Figure 2A), and this decline in [K+]FB led to a reduction in fiber length (Figures 3A and 4) by altering the Vmax during the fiber elongation process (Figure 10). According to our study, [K+]fiber was the factor that affected Vmax (Figure 8) and was regulated by [K+]FB (Figure 9) through K+-osmoregulation. This indicates that elevating [K+]FB above a critical value may help to grow longer fibers. However, upper fruiting branches are more prone to K deficiency disease due, basically, to a larger reduction in [K+]FB compared with lower fruiting branches (Figure 2A). Therefore, adding extra potash fertilizer or allocating it for the development of late season bolls might enhance the K+ level of the upper fruiting branches and result in grow longer fibers. This K fertilization scheme applies to cotton cultivation on K-deficient soil and to the growth of low-K-sensitive cotton cultivars.

Author Contributions

JY collected and analyzed data, and wrote and revised this manuscript; JY and WH completed the cotton plant nutrition determination; WH, WZ, BC, YW, YM, and ZZ helped analyze data and assisted in revising the manuscript; YM and ZZ designed the experiments together and final approved to publish this version. All authors have participated in this study and approved to publish this manuscript.

Funding

This study was financially supported by the National Natural Science Foundation of China (31401327), and Nanjing Agricultural University (NJAU) Youth Science and Technology Innovation Found (KJ2013001).

Conflict of Interest Statement

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.

Acknowledgments

We really appreciate Dr. Yanjiao Dai, Zhen Dai, Wenqi Zhao and all members in our lab, on their helps in the field and lab work. We also thank the technical worker Mrs. Zhang in Pailou Testing Site, Nanjing Agricultural University for her help on cotton management.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.00013

Abbreviations

CS, coefficient of stress; CV, coefficient of variation; DPA, days post anthesis; FB, fruiting branch; FRED, fibre rapid elongation duration; [K+]FB, fruiting branch K+ concentration; [K+]fibre, fibre K+ concentration; [K+]soil, soil soluble K+ concentration; Lenmax, the theoretically maximal fibre length; Lenobs, the observed final fibre length; Temp., temperature; Vmax, the maximum velocity of fibre elongation.

References

Allaway, W. G. (1973). Accumulation of malate in guard cells of Vicia faba during stomatal opening. Planta 110, 63–70. doi: 10.1007/BF00386923

PubMed Abstract | CrossRef Full Text | Google Scholar

Amtmann, A., and Blatt, M. R. (2008). Regulation of macronutrient transport. New Phytol. 181, 35–52. doi: 10.1111/j.1469-8137.2008.02666.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Avci, U., Pattathil, S., Singh, B., Brown, V. L., Hahn, M. G., and Haigler, C. H. (2013). Cotton fiber cell walls of Gossypium hirsutum and Gossypium barbadense have differences related to loosely-bound xyloglucan. PLoS ONE 8:e56315. doi: 10.1371/journal.pone.0056315

PubMed Abstract | CrossRef Full Text | Google Scholar

Bauer, P. J., May, O. L., and Camberato, J. J. (1998). Planting date and potassium fertility effects on cotton yield and fiber properties. J. Prod. Agric. 11, 415–420. doi: 10.2134/jpa1998.0415

CrossRef Full Text | Google Scholar

Boquet, D. J., and Moser, E. B. (2003). Boll retention and boll size among intrasympodial fruiting sites in cotton. Crop Sci. 43, 195–201. doi: 10.2135/cropsci2003.1950

CrossRef Full Text | Google Scholar

Braden, C. A., and Smith, C. W. (2004). Fiber length development in near-long staple upland cotton. Crop Sci. 44, 1553–1559. doi: 10.2135/cropsci2004.2032

CrossRef Full Text | Google Scholar

Campas, O., and Mahadevan, L. (2009). Shape and dynamics of tip-growing cells. Curr. Biol. 19, 2102–2107. doi: 10.1016/j.cub.2009.10.075

PubMed Abstract | CrossRef Full Text | Google Scholar

Cassman, K., Kerby, T., Roberts, B., Bryant, D., and Higashi, S. (1990). Potassium nutrition effects on lint yield and fiber quality of Acala cotton. Crop Sci. 30, 672–677. doi: 10.2135/cropsci1990.0011183X003000030039x

CrossRef Full Text | Google Scholar

Cope, J. (1981). Effects of 50 years of fertilization with phosphorus and potassium on soil test levels and yields at six locations. Soil Sci. Soc. Am. J. 45, 342–347. doi: 10.2136/sssaj1981.03615995004500020023x

CrossRef Full Text | Google Scholar

Cosgrove, D. J. (1997). Cellular mechanisms underlying growth asymmetry during stem gravitropism. Planta 203, S130–S135. doi: 10.1007/PL00008101

PubMed Abstract | CrossRef Full Text | Google Scholar

Cosgrove, D. J. (2000). Loosening of plant cell walls by expansins. Nature 407, 321–326. doi: 10.1038/35030000

PubMed Abstract | CrossRef Full Text | Google Scholar

Davidonis, G. H., Johnson, A. S., Landivar, J. A., and Fernandez, C. J. (2004). Cotton fiber quality is related to boll location and planting date. Agron. J. 96, 42–47. doi: 10.2134/agronj2004.0042

CrossRef Full Text | Google Scholar

Delmer, D. P. (1999). Cellulose biosynthesis: exciting times for a difficult field of study. Annu. Rev. Plant Biol. 50, 245–276. doi: 10.1146/annurev.arplant.50.1.245

PubMed Abstract | CrossRef Full Text | Google Scholar

Dhindsa, R. S., Beasley, C. A., and Ting, I. P. (1975). Osmoregulation in cotton fiber. Plant Physiol. 56, 394–398. doi: 10.1104/pp.56.3.394

CrossRef Full Text | Google Scholar

Dong, H., Tang, W., Li, Z., and Zhang, D. (2004). On potassium deficiency in cotton-disorder, cause and tissue diagnosis. Agric. Conspec. Sci. 69, 77–85.

Google Scholar

Famiani, F., Walker, R. P., Técsi, L., Chen, Z. H., Proietti, P., and Leegood, R. C. (2000). An immunohistochemical study of the compartmentation of metabolism during the development of grape (Vitis vinifera L.) berries. J. Exp. Bot. 51, 675–683. doi: 10.1093/jexbot/51.345.675

PubMed Abstract | CrossRef Full Text | Google Scholar

Frederic, N., Isabelle, H., Alain, J., Samantha, V., Herve, C., and Herman, H. (1998). A plasma membrane-bound putative endo-1,4-β-d-glucanase is required for normal wall assembly and cell elongation in Arabidopsis. EMBO J. 17, 5563–5576. doi: 10.1093/emboj/17.19.5563

CrossRef Full Text | Google Scholar

Fricke, W., and Peters, W. S. (2002). The biophysics of leaf growth in salt-stressed barley. A study at the cell level. Plant Physiol. 129, 374–388.

PubMed Abstract | Google Scholar

Gormus, O. (2002). Effects of rate and time of potassium application on cotton yield and quality in tuikey. J. Agron. Crop Sci. 188, 382–388. doi: 10.1046/j.1439-037X.2002.00583.x

CrossRef Full Text | Google Scholar

Harmer, S., Orford, S., and Timmis, J. (2002). Characterisation of six a-expansin genes in Gossypium hirsutum (upland cotton). Mol. Genet. Genomics 268, 1–9. doi: 10.1007/s00438-002-0721-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, M. G. K., Outlaw, W. H., and Lowry, O. H. (1977). Enzymic assay of 10-7 to 10-4 moles of sucrose in plant tissues. Plant Physiol. 60, 379–383. doi: 10.1104/pp.60.3.379

CrossRef Full Text | Google Scholar

Leigh, R. A., and Wyn Jones, R. G. (1984). A hypothesis relating critical potassium concentrations for growth to the distribution and functions of this ion in the plant cell. New Phytol. 97, 1–13. doi: 10.1111/j.1469-8137.1984.tb04103.x

CrossRef Full Text | Google Scholar

Marschner, H. (1995). Mineral Nutrition of Plants, 2nd Edn. San Diego, CA:: Academic Press.

Google Scholar

Martin, C., Bhatt, K., and Baumann, K. (2001). Shaping in plant cells. Curr. Opin. Plant Biol. 4, 540–549. doi: 10.1016/S1369-5266(00)00213-2

CrossRef Full Text | Google Scholar

Naithani, S., Rao, N. R., and Singh, Y. (1982). Physiological and biochemical changes associated with cotton fibre development. Physiol. Plant. 54, 225–229. doi: 10.1021/pr800550q

PubMed Abstract | CrossRef Full Text | Google Scholar

Oosterhuis, D. M. (1995). “Potassium nutrition of cotton in the USA, with particular reference to foliar fertilization,” in Proceedings of the World Cotton Research Conference – 1: Challenging the Future, eds G. A. Constable and N. W. Forrester (Melbourne: CSIRO), 133–146.

Google Scholar

Oosterhuis, D. M., Abaye, O., Albers, D., Baker, W., Burmester, C., Cothren, J., et al. (1994). A summary of a three-year beltwide study of soil and foliar fertilization with potassium nitrate in cotton. Beltwide Cotton Conf. 3, 1532–1533.

Google Scholar

Pettigrew, W. T. (2008). Potassium influences on yield and quality production for maize, wheat, soybean and cotton. Physiol. Plant. 133, 670–681. doi: 10.1111/j.1399-3054.2008.01073.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pettigrew, W. T., Heitholt, J. J., and Meredith, W. R. (1996). Genotypic interactions with potassium and nitrogen in cotton of varied maturity. Agron. J. 88, 89–93. doi: 10.2134/agronj1996.00021962008800010019x

CrossRef Full Text | Google Scholar

Qin, Y. M., and Zhu, Y. X. (2011). How cotton fibers elongate: a tale of linear cell-growth mode. Curr. Opin. Plant Biol. 14, 106–111. doi: 10.1016/j.pbi.2010.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Raschke, K., and Schnabl, H. (1978). Availability of chloride affects the balance between potassium chloride and potassium malate in guard cells of Vicia faba L. Plant Physiol. 62, 84–87. doi: 10.1104/pp.62.1.84

PubMed Abstract | CrossRef Full Text | Google Scholar

Read, J. J., Reddy, K. R., and Jenkins, J. N. (2006). Yield and fiber quality of upland cotton as influenced by nitrogen and potassium nutrition. Eur. J. Agron. 24, 282–290. doi: 10.1016/j.eja.2005.10.004

CrossRef Full Text | Google Scholar

Reddy, K. R., and Zhao, D. (2005). Interactive effects of elevated CO2 and potassium deficiency on photosynthesis, growth, and biomass partitioning of cotton. Field Crop Res. 94, 201–213. doi: 10.1016/j.fcr.2005.01.004

CrossRef Full Text | Google Scholar

Römheld, V., and Kirkby, E. A. (2010). Research on potassium in agriculture: needs and prospects. Plant Soil 335, 155–180. doi: 10.1007/s11104-010-0520-1

CrossRef Full Text | Google Scholar

Ruan, Y. L. (2005). Recent advances in understanding cotton fibre and seed development. Seed Sci. Res. 15, 269–280. doi: 10.1079/ssr2005217

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruan, Y. L. (2007). Rapid cell expansion and cellulose synthesis regulated by plasmodesmata and sugar: insights from the single-celled cotton fibre. Funct. Plant Biol. 34, 1–10. doi: 10.1071/FP06234

CrossRef Full Text | Google Scholar

Ruan, Y. L., Chourey, P. S., Delmer, D. P., and Perez, G. L. (1997). The differential expression of sucrose synthase in relation to diverse patterns of carbon partitioning in developing cotton seed. Plant Physiol. 115, 375–385.

PubMed Abstract | Google Scholar

Ruan, Y. L., Llewellyn, D. J., and Furbank, R. T. (2001). The control of single-celled cotton fiber elongation by developmentally reversible gating of plasmodesmata and coordinated expression of sucrose and K+ transporters and expansin. Plant Cell 13, 47–60. doi: 10.1105/tpc.13.1.47

PubMed Abstract | CrossRef Full Text | Google Scholar

Samal, D., Kovar, J. L., Steingrobe, B., Sadana, U. S., Bhadoria, P. S., and Claassen, N. (2010). Potassium uptake efficiency and dynamics in the rhizosphere of maize (Zea mays L.), wheat (Triticum aestivum L.), and sugar beet (Beta vulgaris L.) evaluated with a mechanistic model. Plant Soil 332, 105–121. doi: 10.1007/s11104-009-0277-6

CrossRef Full Text | Google Scholar

Schubert, A., Benedict, C., Berlin, J., and Kohel, R. (1973). Cotton fiber development-kinetics of cell elongation and secondary wall thickening. Crop Sci. 13, 704–709. doi: 10.2135/cropsci1973.0011183X001300060035x

CrossRef Full Text | Google Scholar

Shimizu, Y., Aotsuka, S., Hasegawa, O., Kawada, T., Sakuno, T., Sakai, F., et al. (1997). Changes in levels of mRNAs for cell wall-related enzymes in growing cotton fiber cells. Plant Cell Physiol. 38, 375–378. doi: 10.1093/oxfordjournals.pcp.a029178

PubMed Abstract | CrossRef Full Text | Google Scholar

Thaker, V. S., Saroop, S., Vaishnav, P., and Singh, Y. (1989). Genotypic variations and influence of diurnal temperature on cotton fibre development. Field Crop Res. 22, 129–141. doi: 10.1016/0378-4290(89)90063-4

CrossRef Full Text | Google Scholar

Wang, L., and Ruan, Y. L. (2010). Unraveling mechanisms of cell expansion linking solute transport, metabolism, plasmodesmtal gating and cell wall dynamics. Plant Signal. Behav. 5, 1561–1564. doi: 10.4161/psb.5.12.13568

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, N., Hua, H., Egrinya, E. A., Li, Z., Duan, L., and Tian, X. (2012). Genotypic variations in photosynthetic and physiological adjustment to potassium deficiency in cotton (Gossypium hirsutum L.). J. Photochem. Photobiol. B 110, 1–8. doi: 10.1016/j.jphotobiol.2012.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Warncke, D., and Brown, J. (1998). “Potassium and other basic cations,” in Recommended Chemical Soil Test Procedures for the North Central Regional Research Publication No. 211, eds M. Nathan and R. Gelderman (Columbia City, MO: Missouri Agricultural Experiment Station), 31–33.

Google Scholar

Wright, P. (1999). Premature senescence of cotton (Gossypium hirsutum L.) Predominantly a potassium disorder caused by an imbalance of source and sink. Plant Soil 211, 231–239. doi: 10.1023/A:1004652728420

CrossRef Full Text | Google Scholar

Yang, J., Zhao, W., Hu, W., Wang, Y., Chen, B., and Zhou, Z. (2014). Indicators of cotton (Gossypium hirsutum L.) cultivar screening for low-potassium tolerance in seedling stage and its relationship with yield and quality. Cotton Sci. 26, 283–289.

Google Scholar

Zhang, Z., Tian, X., Duan, L., Wang, B., He, Z., and Li, Z. (2007). Differential responses of conventional and Bt-transgenic cotton to potassium deficiency. J. Plant Nutr. 30, 659–670. doi: 10.1080/01904160701289206

CrossRef Full Text | Google Scholar

Zhang, Z. Y., Wang, Q. L., Li, Z. H., Duan, L. S., and Tian, X. L. (2009). Effects of potassium deficiency on root growth of cotton seedlings and its physiological mechanisms. Acta Agron. Sin. 35, 718–723. doi: 10.3724/SP.J.1006.2009.00718

CrossRef Full Text | Google Scholar

Zhao, W. Q., Wang, Y. H., Zhou, Z. G., Meng, Y. L., Chen, B. L., and Oosterhuis, D. M. (2012). Effect of nitrogen rates and flowering dates on fiber quality formation of cotton. Am. J. Exp. Agric. 2, 133–159. doi: 10.9734/AJEA/2012/954

CrossRef Full Text | Google Scholar

Keywords: cell turgor, cotton, fiber elongation, K deficiency, [K+]FB, [K+]fiber, Vmax

Citation: Yang J, Hu W, Zhao W, Chen B, Wang Y, Zhou Z and Meng Y (2016) Fruiting Branch K+ Level Affects Cotton Fiber Elongation Through Osmoregulation. Front. Plant Sci. 7:13. doi: 10.3389/fpls.2016.00013

Received: 01 July 2015; Accepted: 07 January 2016;
Published: 22 January 2016.

Edited by:

Richard William Bell, Murdoch University, Australia

Reviewed by:

Ryoung Shin, Riken Center for Sustainable Resource Science, Japan
Rumen Ivanov, Heinrich Heine University Düsseldorf, Germany

Copyright © 2016 Yang, Hu, Zhao, Chen, Wang, Zhou and Meng. 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) or licensor 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: Zhiguo Zhou, giscott@njau.edu.cn; Yali Meng, mengyl@njau.edu.cn

Disclaimer: 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.