- 1Key Laboratory of Crop Genetics and Physiology of Jiangsu Province/Co-Innovation Center for Modern Production Technology of Grain Crops of Jiangsu Province, Yangzhou University, Yangzhou, Jiangsu, China
- 2Research Department for Agro-ecological Security, Jiangsu Coastal Area Institute of Agricultural Sciences, Yancheng, Jiangsu, China
- 3Faculty of Horticultural Science and Technology, Suzhou Polytechnic Institute of Agriculture, Suzhou, Jiangsu, China
- 4College of Environmental Science and Engineering, Yangzhou University, Yangzhou, Jiangsu, China
- 5State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing, Jiangsu, China
Introduction: This study explores how elevated CO2 concentration may alter the source-sink dynamics in rice by providing additional carbon for photosynthesis, thereby affecting nutrient absorption and distribution.
Methods: A free-air CO2 enrichment experiment was conducted on a japonica cultivar Wuyunjing 27 in 2017 and 2018 growing seasons. The plants were exposed to ambient and elevated CO2 level (increased by 200 μmol·mol-1) and two source-sink manipulation treatments (control with no leaf cutting and cutting off the top three leaves at heading).
Results: The elevated CO2 significantly increased the above-ground biomass and the straw non-structural carbohydrate concentration by an average of 19.3% and 12.5%, respectively. Significant changes in the concentrations of N, S, Fe, and Zn in straw were noted under elevated CO2, with average decreases by 7.1, 7.2, 11.6, and 10.1%, respectively. The exposure to elevated CO2 significantly enhanced the elements accumulation, yet it had minimal impact on their distribution across different organs. When compared to intact rice, removing the top three leaves at heading reduced the above-ground biomass by 36.8% and the straw non-structural carbohydrate content by 44.8%. Leaf-cutting generally increased the concentration of elements in stem, leaf, and grain, likely due to a concentration effect from reduced biomass and carbohydrate accumulation. Leaf-cutting reduced element accumulation and shifted element allocation in rice organs. It increased the proportion of elements in stems while reduced their presence in leaves and grains.
Discussion: Our study suggests that a dilution effect may cause a decrease in mineral elements concentrations under elevated CO2 because of the increase in biomass and carbohydrates.
Introduction
Atmospheric carbon dioxide (CO2) levels have risen significantly from approximately 280 μmol·mol-1 during the Industrial Revolution to about 420 μmol·mol-1 today (NOAA, 2024), with projections suggesting a rise to 550 μmol·mol-1 by 2050 (IPCC, 2021). As a critical substrate for photosynthesis, elevated CO2 (eCO2) has profound implications for rice, one of the main stable food crops globally (Fan et al., 2016). Although eCO2 can enhance rice yield (Hu et al., 2021, Hu et al., 2022a; Ainsworth and Long, 2021), it often reduces the concentration of crucial mineral elements such as N, S, Mg, Fe, Cu, and Zn (Hu et al., 2022b; Hu et al., 2024). This reduction could pose ‘invisible hunger’ risks (Loladze, 2014) and impact the biochemical interactions between rice and soil, thus influencing overall soil health and sustainability. Given these factors, investigating the dynamics of mineral element absorption and distribution in rice under eCO2 conditions is vital for ensuring food security and improving the quality of arable land in the face of climate change.
Leaves under eCO2 conditions produce increased non-structural carbohydrates (NSC) (Zhu et al., 2016), which are transferred to the grains during the filling stage, resulting in a ‘dilution effect’ that associated with the reduction of concentrations of elements such as N, Fe, S and Zn (Zhang et al., 2013; Myers et al., 2014; Chaturvedi et al., 2017; Zhu et al., 2018; Ziska, 2022; Gojon et al., 2023). However, findings from free air CO2 enrichment (FACE) experiments have shown mixed outcomes. Some studies found no significant decrease in mineral element concentration in grains, while others reported increased concentrations of certain elements (Lieffering et al, 2004; Wang et al., 2014). This variability highlights the need for ongoing research in this domain, especially considering the potential influence of climatic variations across different years (Yang et al., 2007; Li et al., 2023).
Several studies conducted under FACE conditions have shown a general trend of increased accumulation of elements such as Fe, Mn, Ca, Mg, and Zn in rice stems and panicles but not in leaves (Lieffering et al, 2004; Zhou et al., 2021; Cui et al., 2023; Wang et al., 2023). Furthermore, studies suggest an enhanced distribution of these elements to rice panicle, while their presence in stem and leaf diminishes (Li et al., 2020; Jiang et al., 2020a). A study found that potassium (K) distribution in leaf decreased significantly by 16.3% under eCO2 conditions, while no substantial change was observed in stem and panicle (Pang et al., 2005). This indicates differential responses among elements to eCO2. The stem and leaf, which exposed to high CO2 concentrations more directly and longer than grain enclosed in a glume shell, may be more affected by eCO2. According to a study by Ujiie et al (2019), the concentrations of S, Mg, and Mn in straw decreased by 1.5-4 times as much as in grain under eCO2 conditions. Nutrient elements in grains are absorbed from the soil by roots and initially accumulate in vegetative organs such as stems and leaves before being transferred to grain. Therefore, to understand the changes in grain element composition by eCO2, it is essential to study the element accumulation and distribution in vegetative organs in the response to high CO2 concentrations. However, research in this area is limited and has mainly focused on nitrogen (Li et al., 2015), with less attention given to other elements.
The increase in atmospheric CO2 concentration has significant implications for rice photosynthesis (Dingkuhn et al., 2020; Fabre et al., 2020; Yin et al., 2022) and the absorption and transport of elements within the plant (Wei et al., 2018; Guo et al., 2015). Prior FACE studies showed that the concentrations of elements (such as N and Cu) in grain decreased under eCO2 while increased under partial defoliation treatment (Gao et al., 2021a), demonstrating the influence of alterations in the source-sink relationship on grain element concentrations. The accumulation of nutrients in grain during the grain-filling stage is mainly accomplished through mineral remobilization from vegetative organs and/or uptake by roots (Wu et al., 2010; Yuan et al., 2012). According to Sperotto et al. (2013), reducing source-sink ratio by flag leaf removal at anthesis does not change Fe and Zn concentration or content in mature seeds. It is possible that flag leaves are preferential but are not essential as a source of metals, which can probably be compensated by other leaves remobilization and/or continuous uptake by roots. However, it is not clear whether this assumption remains viable in the case of a significant increase in carbohydrate release. Moreover, it remains unclear how the source-sink relationship influences element absorption and distribution throughout the rice plant under eCO2 concentrations. This study employed a FACE setting with two distinct CO2 levels, ambient CO2 (aCO2) and elevated CO2 (eCO2, aCO2 + 200 μmol·mol-1), and established two source-sink ratio treatments: a conventional source-sink ratio (CK) and a reduced source-sink ratio by cutting off the top three leaves (LC). The investigation, conducted over two consecutive years, aims to elucidate the impacts of elevated CO2 and source-sink relationships on mineral elements uptake and distribution in rice.
Materials and methods
Experimental site and FACE system
This study was conducted at the FACE system located in Yangzhou (119°42′0″E, 32°35′5″N), Jiangsu Province, China. The site featured a soil composition of organic carbon (24.8 g·kg-1), total nitrogen (1.13 g·kg-1), total phosphorus (0.54 g·kg-1), and total potassium (9.7 g·kg-1), with available nitrogen, phosphorous, and potassium concentrations of 122.4 mg·kg-1, 15.1 mg·kg-1, and 56.5 mg·kg-1 respectively, and a pH of 6.9. The FACE system comprised six plots distributed across paddies with similar soil characteristics and agronomic history. Three plots were randomly assigned to elevated CO2 treatment and three to ambient conditions, each measuring approximately 80 m2. To prevent CO2 cross-contamination, the centers of the FACE and ambient plots were spaced 90 meters apart. Pure CO2 was delivered into the center of the FACE plots through perimeter pipelines. CO2 concentrations were controlled by a computer system adjusting for atmospheric levels, wind direction, and speed, ensuring a consistent 200 μmol·mol-1 above ambient CO2 levels at crop canopy height throughout the main growth period of rice. CO2 fumigation was conducted daily from seedling transplantation to maturity, from sunrise to sunset. The average daily temperatures from June 1 to October 31 were 24.9°C in 2017 and 25.5 °C in 2018. The average monthly precipitation was 149.1 mm in 2017 and 120.1 mm in 2018, with monthly sunshine durations of 176.5 and 200.2 hours, respectively, as detailed in Figure 1.
Figure 1. Daily mean temperature (A), monthly mean precipitation (B) and monthly mean sunshine duration (C) at the experiment site during 2017-2018 growing seasons.
Crop cultivation and defoliation treatment
In each experimental season, the japonica cultivar Wuyunjing 27, a popular cultivar in the study region, was used. Seeds were sown in a nursing paddy, and seedlings were grown for one month under ambient conditions. Seedlings were manually transplanted to the experimental plots at a density of two per hill around June 20 each season. Fertilization included urea (46% N) and a balanced NPK compound (N: P2O5: K2O = 15:15:15) applied at a rate of 22.5 g N·m-2. In each growing season, nitrogen was administered into three applications: 40% as a basal dressing one day before transplanting, 30% as a top dressing at the early tillering stage, and 30% as a top dressing at panicle ignition. Phosphorus and potassium were also applied as compound fertilizers P2O5 and K2O at a rate of 9 g·m-2 one day before transplanting. The water regime was meticulously managed with paddy fields submerged to a depth of 5 cm from about June 17 to July 20, followed by a period of wet-dry cycles from about July 21 to August 10 through natural drainage and intermittent irrigation. Routine pest and disease management was routinely conducted throughout the growing seasons. Fertilizer application and water regime specifics are available in Gao et al. (2021a).
At the heading stage, 30 plants per plot were selected for the study. To investigate source-sink dynamics, half of these plants underwent a leaf-cutting treatment there the top three leaves were removed (LC), and the other half were left intact as controls (CK). This method ensured uniformity in tiller number and developmental stage across treatments, critical for evaluating the impact of altered source-sink relationships on growth and yield.
Plant sampling and parameter measurements
At maturity, five hills, each representing the average number of tillers per hill, were selected for destructive sampling. These samples were systematically divided into leaves, stems, and panicles. To determine the accurate dry weights, stems and leaves were oven-dried at 105 °C for 30 min, followed by 80 °C for 72 h. Panicles were air-dried to a constant weight. After drying, panicles were carefully threshed and dehulled to produce brown rice. Above-ground biomass was calculated from the dry weights of each organ. Stems and leaves were ground to less than 0.15 mm using a stainless-steel mill, and brown rice was ground into powder using a vibration disk mill (TS1000, Germany) with a 100-mesh sieve, preparing them for subsequent analyses.
The non-structural carbohydrate (NSC) content was determined as the sum of total soluble sugars and starch. The anthrone H2SO4 method, as described by Yoshida et al. (1976), was employed to measure soluble sugars and starch in the stems and leaves. Nitrogen (N) concentrations were determined using a Kjeltec 8400 autoanalyzer (FOSS Analytical AB, Sweden) after hydrolysis by the Kjeldahl method. The concentrations of other elements were determined as follows: 0.50 g of flour was weighed into the lining tube of a microwave digestion apparatus (MARS5, CEM Corporation, Matthews, USA), followed by the addition of 5 mL of 65% nitric acid, 3 mL of ultrapure water, and three drops of hydrogen peroxide; the mixture was then placed into a microwave digestion-meter for high-temperature digestion. Post-digestion, the volume of the digestion liquid was increased to 50 mL and filtered using quantitative filter paper. The filtrate’s mineral content was analyzed using an inductively coupled plasma emission spectrometer (iCAP6300, Thermo Fisher Scientific, USA). This analysis included the macroelements phosphorus (P), sulfur (S), potassium (K), calcium (Ca), and magnesium (Mg) and the microelements boron (B), copper (Cu), iron (Fe), manganese (Mn), and zinc (Zn). Element accumulations were calculated by multiplying the biomass of each organ by its elemental concentration. Elemental allocations were then determined based on the total accumulation across all sampled organs.
Statistical analyses
The field experiment was designed as a completely randomized setup with a split-plot arrangement. The CO2 levels were treated as the main plot variable, while source-sink treatments acted as subplots, each with three replications. Analysis of variance (ANOVA) was conducted using the SPSS statistical software (Version 20.0, SPSS Inc., Chicago, USA) to assess the effects of these treatments. Post-hoc comparisons among treatment means were performed using the Least Significant Difference (LSD) test. Statistical significance of the effects was indicated as follows: highly significant (** P < 0.01), significant (* P < 0.05), and marginally significant (+ P < 0.1). Additionally, Pearson’s correlation coefficients were calculated to explore the relationships among the different measured parameters.
Results
Above-ground biomass response
In the present two-year study, the above-ground biomass (AGB) of rice was higher in 2018 than in 2017 (Figure 2 and Table 1). As compared to the ambient, on average, eCO2 increased the AGB by19.3% significantly, and the increase in 2017 and 2018 were close to each other. Among the different plant organs, eCO2 increased the biomass of stem, leaf, and grain by 23.0%, 10.3%, and 18.3%, respectively. The magnitude of AGB increase by eCO2 was larger for LC crops (+23.7%) compared to CK crops (+16.6%), though no interaction of CO2 with LC was detected (Table 1). Compared with CK crops, the LC treatment significantly reduced AGB by 36.8% on average, with reductions of 32.7% and 40.3% observed in 2017 and 2018, respectively (Table 1).
Figure 2. Response of rice biomass at maturity to elevated CO2 under CK (no leaf cutting) and LC (cutting off top three leaves at heading) in 2017 and 2018 growing seasons. Each bar in the figure represents the mean values across three plots for ambient CO2 (AC) or elevated CO2 (EC, AC + 200 μmol·mol-1); vertical bars represent standard error (n = 3). Statistically significant effects are indicated as ** P < 0.01; * P < 0.05; + P < 0.1.
Table 1. Analysis of variance (ANOVA) results of effects of elevated CO2 and LC (leaf-cutting treatment at heading) on rice biomass at maturity in 2017 and 2018 growing seasons.
Element concentrations
The concentrations of eleven mineral elements including N, Ca, K, Mg, P, S, B, Cu, Fe, Mn, and Zn in rice straw and grain were determined in the present study (Figures 3, 4; Supplementary Tables S1, S2). Apart from Mg, B, and Zn, the concentrations of other elements in straw were lower in 2018 than in 2017. Relative to the concentrations of grain elements in 2017, higher K, Mg, P, and B concentrations but lower N, Ca Fe, Mn, and Zn concentration were found in 2018. Among these elements, N, S, Fe, and Zn concentrations in straw were significantly altered by eCO2 when averaged across different source-sink crops in two years, with average decreases by 7.1, 7.2, 11.6, and 10.1%, respectively. Response of these elements in grain was not significant over the growing seasons of 2017 and 2018 (Supplementary Table S1). Compared to CK crops, the concentrations of N, K, Mg, P, S, and Zn in grain showed a significant increase in response to LC treatment, ranging from 3.6% to 15.1% on average, similar responses of stem and leaf were observed in two seasons.
Figure 3. Response of macroelements concentrations in rice straw, stem, and leaf at maturity to elevated CO2 under CK (no leaf cutting) and LC (cutting off top three leaves) in 2017 and 2018 growing seasons. Each bar in the figure represents the mean values across three plots for ambient CO2 (AC) or elevated CO2 (EC, AC + 200 μmol·mol-1); vertical bars represent standard error (n = 3). Statistically significant effects are indicated as ** P < 0.01; * P < 0.05; + P < 0.1.
Figure 4. Response of microelements concentrations in rice straw, stem, and leaf at maturity to elevated CO2 under CK (no leaf cutting) and LC (cutting off top three leaves) in 2017 and 2018 growing seasons. Each bar in the figure represents the mean values across three plots for ambient CO2 (AC) or elevated CO2 (EC, AC+200 μmol·mol-1); vertical bars represent standard error (n = 3). Statistically significant effects are indicated as ** P < 0.01; * P < 0.05; + P < 0.1.
Non-structural carbohydrates
On average, eCO2 increased straw non-structural carbohydrates (NSC) content by 12.5% (Figure 5), with a more pronounced CO2 effect occurring in 2018 (+16.2%, P < 0.01) than in 2017 (+7.3%, P > 0.1), which resulted in a significant CO2 × year interaction (Supplementary Table S2). In addition, this parameter showed a significant CO2 × organ interaction: it was increased in stem by 14.2% (P < 0.01) by eCO2, but had no effects in leaf. A strong LC effect was observed on straw NSC content, with a large decrease by 46.8% and 42.7% in 2017 and 2018, respectively. Under the LC treatment, the NSC decrease was greater in stem (-49.0%) than in leaf (-35.7%), with a notable significant LC × organ interaction (Supplementary Table S2).
Figure 5. Response of non-structural carbohydrate (NSC) concentrations in rice straw (A), stem (B), and leaf (C) at maturity to elevated CO2 under CK (no leaf cutting) and LC (cutting off top three leaves) in 2017 and 2018 growing seasons. Each bar in the figure represents the mean values across three plots for ambient CO2 (AC) or elevated CO2 (EC, AC + 200 μmol·mol-1); vertical bars represent standard error (n = 3). Statistically significant effects are indicated as * P < 0.05; + P < 0.1.
Correlation analysis
Pearson correlation analyses were conducted to investigate the relationship between the concentrations of elements and the biomass or NSC in different organs (Figure 6). The concentrations of N, Ca, K, P, S, Fe, and Mn in stems and leaves were negatively correlated with biomass and NSC, whereas concentrations of N, S, Fe, Mn, and Zn in grain were negatively correlated with grain yield. In contrast, a significant positive linear correlation between biomass and NSC concentration was observed in stems and leaves.
Figure 6. Correlation analysis of elements concentration in different parts of rice with biomass, non-structure carbohydrate (NSC) concentrations of stem and leaf under CO2 and leaf-cutting treatments. Red and blue showed significant positive and negative correlation, while white showed no significant correlation. S-Biomass, L-Biomass, and G-Biomass represent stem, leaf, and grain biomass respectively. S-NSC and L-NSC represent NSC concentrations in stem and leaves, respectively. S-N, S-Ca, S-K, S-Mg, S-P, S-S, S-B, S-Cu, S-Fe, S-Mn, and S-Zn represent concentrations of N, Ca, K, Mg, P, S, B, Cu, Fe, Mn, and Zn in stem, respectively. L-N, L-Ca, L-K,L-Mg, L-P, L-S, L-B, L-Cu, L-Fe, L-Mn, and L-Zn represent concentrations of N, Ca, K, Mg, P, S, B, Cu, Fe, Mn, and Zn in leaves, respectively. G-N, G-Ca, G-K, G-Mg, G-P, G-S, G-B, G-Cu, G-Fe, G-Mn, and G-Zn represent concentrations of N, Ca, K, Mg, P, S, B, Cu, Fe, Mn, and Zn in grain, respectively. * and ** indicate significance at P < 0.05, and P < 0.01, respectively. (n = 24).
Element accumulation and allocation
We multiplied the elemental concentration by biomass of each organ to calculate the total amounts of the various elements removed by harvesting (Table 2). Compared with the ambient, eCO2 significantly increased accumulation of elements in both years, with the increase ranging from 5.8% to 29.9% on average. In contrast, the LC treatment significantly decreased the accumulation of mineral elements, with similar magnitude in 2017 and 2018.
Table 2. Response of mineral elements accumulations of rice aboveground to elevated CO2 and LC (cutting off top three leaves at heading) in 2017 and 2018 growing seasons.
The effects of all main factors (CO2, year, and LC) on the elements distributions at maturity differed among the three organs (Figure 7 and Supplementary Table S3). Averaged across two years and different source-sink ratio crops, eCO2 significantly increased Zn allocation to grain by 7.2%, and Cu and Fe allocation to leaf by 22.1% and 16.1%, while respectively decreased Fe and Mn allocations to stem by 3.1% and 3.7%. Significant year effects were observed on the elements percentages of each organs at maturity (Supplementary Table S3). Relative to these elements fractions in 2017, higher fractions of grain K, P, and B, and stem K, Mn, and Zn, and leaf K and Mn were found in 2018. As compared to CK crops, LC treatment significantly altered the elements distributions. Overall, LC treatment significant increased elements allocation to stem but decreased those to grain and leaf, and similar trends were seen in 2017 and 2018, though interactions of LC with year were detected on some elements allocations (Supplementary Table S3).
Figure 7. The mineral elements allocations in rice stem, leaf, and grain as affected by elevated CO2 and LC (cutting off top three leaves) in 2017 and 2018 growing seasons. AC and EC refer to ambient CO2 and elevated CO2, respectively. CK and LC refer to no leaf cutting and cutting off top three leaves, respectively. Each bar in the figure represents the mean values across three plots for ambient CO2 or elevated CO2.
Discussion
Above-ground biomass
Increasing atmospheric CO2 concentration has a fertilization effect that contributes to rice growth, as CO2 is essential for photosynthesis. Numerous studies consistently demonstrate this effect, including those in artificial climate chambers, open-top chambers (OTCs), and free air CO2 enrichment (FACE) platforms (Ainsworth and Long, 2021). Among these platforms, FACE is generally believed that FACE is closer to the natural ecological environment, thus truly reflecting the growth of rice under future elevated CO2 conditions. A meta-analysis of 20 years FACE studies suggested that, on average, elevated CO2 (eCO2) increased above-ground biomass of rice by 19%, and the increase is generally more pronounced in panicles than in stems or leaves (Hu et al., 2022a). From the perspective of grain yield alone, an average increase of 16% under eCO2 conditions was observed by Hu et al. (2021). The yield enhancement by elevated CO2 varied across different rice varieties (Ainsworth and Long, 2021; Yang et al., 2007; Kim et al., 2003), with hybrid rice (24.7%) was significantly greater than conventional rice (14.2%) (Hu et al., 2021). The leaf-cutting (LC) treatment led to a notable reduction in biomass across all above-ground plant parts due to decreased photosynthetic productivity and assimilate availability in the stem and panicle, culminating in reduced biomass.
While the ‘CO2 fertilization’ effect is well-documented to increase yields, the actual gains in field conditions were often lower than expected (Long et al., 2006; Ainsworth and Long, 2021; Sonnewald and Fernie, 2018). This discrepancy was increasingly attributed to the source-sink ratio of rice (Lv et al., 2020; Yin et al., 2022). Previous FACE studies have demonstrated that modifying the source-sink ratio by thinning flowers or cutting off leaves significantly affected rice yield response to eCO2 (Gao et al., 2021b). Specifically, reducing the source-sink ratio appeared more effective in enhancing biomass under eCO2, suggesting potential strategies for maximising yield benefits in future high CO2 conditions (Dingkuhn et al., 2020; Lawlor and Paul, 2014; Paul et al., 2020).
Element concentration dynamics
Nitrogen (N) is an essential element of rice protein, an important element affecting leaf photosynthesis (Cai et al., 2020), and a nutrient element strongly influenced by elevated CO2 concentration (Gojon et al., 2023). Research has indicated that elevated CO2 enhanced the synergistic cycles of nitrogen and carbon (Cui et al., 2023), boosting the translocation of N from vegetative parts to grains (Xu et al., 2016; Jiang et al., 2020b), which consequently reduces its concentration in stems and leaves. Additionally, eCO2 often results in reduced stomatal conductance, diminishing leaf transpiration and potentially lowering N uptake due to decreased soil nutrient demand (Pang et al., 2006; Pozo et al., 2007). Besides N, in this study, we observed consistent decreases in S, Cu, Fe, Mn, and Zn under eCO2, which corroborates findings from other FACE studies (Uddling et al., 2018; Myers et al., 2014; Zhu et al., 2018). These declines are typically attributed to the ‘dilution effect’ (Pang et al., 2005; Chaturvedi et al., 2017; Loladze, 2014; Li et al., 2015; Ujiie et al., 2019; Jiang et al., 2020b), where increased plant biomass leads to lower nutrient concentrations. Correlation analysis confirmed significantly negative relationships between the biomass and the concentrations of N, S, Fe, and Mn in plant tissue.
Conversely, the concentrations of K and P increased, which does not align with the typical dilution effect narrative (Ren et al., 2007; Wang et al., 2023). K is crucial for carbohydrate synthesis and transport, and P, a key player in energy transfer and photosynthetic regulation (Ellsworth et al., 2022), has shown higher demand and uptake under eCO2 (Shimono et al., 2009; Wang et al., 2023). This suggests that eCO2 affects growth and modifies nutrient uptake priorities based on physiological needs, indicating a complex interaction between CO2 levels and nutrient dynamics that transcends simple biomass accumulation. Interestingly, the leaf-cutting (LC) treatment increased nearly all measured mineral element concentrations across different plant organs, contrasting with the broader biomass reduction and demonstrating a concentration effect due to reduced tissue mass. However, reductions in K and Mn in leaves under LC treatments suggest a decline in photosynthetic efficiency and sugar transport, which affects the demand for these specific nutrients.
The variability observed in element concentrations year-over-year also highlights the significant influence of annual climatic conditions on nutrient dynamics. For instance, the lower element content in stems in 2018 compared to 2017 could be linked to higher biomass production driven by less rainfall and more sunlight (Li et al., 2015; Wang et al., 2024), illustrating the intricate connections between environmental factors, plant physiology, and nutrient uptake under changing climatic conditions. This nuanced understanding of how elevated CO2 and defoliation affect nutrient concentrations provides valuable insights into managing crop nutrition in the face of global climate change, emphasizing the need for adaptive strategies that consider both environmental influences and intrinsic plant responses.
Non-structural carbohydrate dynamics
Non-structural carbohydrates (NSCs) are crucial for plant growth, serving as primary metabolites that regulate energy and structural functions within plant systems. Consistent with Zhu et al. (2016) findings, an increase in NSC concentrations in rice straw was observed in the present study when exposed to eCO2, with increments ranging between 22-70%. Notably, our results indicated a lower increase in NSC concentration than previous studies, likely due to varietal differences. Specifically, the japonica rice exhibited a lesser response to eCO2 compared to the indica rice. From the perspective of source-sink balance, compared with CK crops, the NSC concentration of straw was significantly reduced by LC treatment. This variation could be attributed to the source-sink dynamics within the plants, where the top three functional leaves, crucial for photosynthetic activity, significantly influence NSC production and distribution. Under eCO2 conditions, enhanced CO2 levels act as a vital carbon source for photosynthesis, potentially leading to an increased synthesis of assimilates. However, if the sink capacity, primarily the grain, does not adequately accommodate this increased production, excess carbohydrates may accumulate in the straw, thereby limiting potential yield improvements.
Elemental accumulation and allocation
Increased carbon capture under eCO2 conditions leads to greater carbohydrate accumulation in rice, while will the accumulation of mineral elements increase as a result? This study showed that elemental accumulation in rice enhanced with an increased carbon source provided by eCO2 yet diminishes when the source is decreased by LC. This trend aligns with biomass response patterns under both conditions. Notably, eCO2 has been shown to increase the absorption of N and P from the soil, altering the nutrient cycling and potentially impacting soil quality adversely (Cui et al., 2023; Zhang et al., 2022; Gojon et al., 2023). Over time, enhanced CO2 levels could significantly reduce soil available phosphorus, as demonstrated in a 15-year FACE study (Wang et al., 2023). The adventitious root length and adventitious root number could be increased under eCO2 by 25-37% (Yang et al., 2008), which influences the rhizosphere microbial environment and promotes soil organic matter decomposition, underscores rice dynamic changes in nutrient uptake (Drissner et al, 2007; Jin et al., 2014). These findings suggest that future strategies in rice cultivation should focus on improving soil organic matter and fertilizer utilization to mitigate nutrient elements loss under climate change conditions.
The distribution of mineral elements within rice plants is notably influenced by eCO2, with varying effects across different organs and varieties. This study indicates a shift in the allocation of elements like Fe, Mn, and Zn, with a reduction in stem and an increase in grain. These changes could be attributed to the interaction between CO2 levels and environmental factors across different years. Interestingly, while the proportion of elements like N, Ca, and Zn in grains fluctuated annually, eCO2 generally promoted the retention of Fe, Mn, Cu, Zn, and Ca in rice husks, reducing their availability in brown rice (Jiang et al., 2020a), suggesting a need for further investigation into the effects of eCO2 on element distribution in ‘flow’ organs such as branches and husks. Additionally, the observed inconsistencies in element allocation under LC treatment result from the complex interactions between photosynthetic activity, nutrient uptake, and transport processes, which can lead to a concentration effect in grains despite unchanged dry matter distribution in panicles. This phenomenon underscores the intricate balance between source and sink dynamics, necessitating further exploration to fully understand the impacts of climatic variations on crop quality and yield.
Conclusion
Results from FACE experiments conducted over two consecutive rice seasons showed that an increase in atmospheric CO2 concentration resulted in a significant increase in rice biomass and NSC concentration of straw. However, the concentrations of mineral elements showed a decreasing trend under eCO2, and the response of the stem was more pronounced than that of the leaf and grain. When compared to naturally growing rice, treatments that reduced the source-sink ratio by cutting off leaves significantly decreased rice biomass and straw NSC concentration but increased the element concentrations. From the perspective of source-sink balance, the increased source (eCO2) treatment may dilute some mineral elements, while the reduced source (LC) treatment may concentrate them. Increased CO2 concentration causes rice to accumulate more mineral elements from the soil, exacerbating soil mineral deficiency risk. Integrating these findings, we consider that future strategies in rice cultivation should focus on improving soil organic matter and fertilizer utilization, in order to mitigate nutrient elements loss under climate change conditions.
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.
Author contributions
BG: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing. SH: Data curation, Writing – review & editing. MZ: Data curation, Writing – review & editing. LJ: Data curation, Funding acquisition, Writing – review & editing. YXW: Funding acquisition, Supervision, Validation, Writing – review & editing. JZ: Supervision, Validation, Writing – review & editing. XS: Data curation, Writing – review & editing. KW: Data curation, Writing – review & editing. YLW: Funding acquisition, Supervision, Writing – review & editing. LY: Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (No.32372216, 32172102, 31671618, and 31571597) and the Priority Academic Program Development of Jiangsu Higher Education Institutions.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1450893/full#supplementary-material
References
Ainsworth, E. A., Long, S. P. (2021). 30 years of free-air carbon dioxide enrichment (FACE): What have we learned about future crop productivity and its potential for adaptation? Global Change Biol. 27, 27–49. doi: 10.1111/gcb.15375
Cai, C., Li, G., Di, L. J., Ding, Y. J., Fu, L., Guo, X. H., et al. (2020). The acclimation of leaf photosynthesis of wheat and rice to seasonal temperature changes in T-FACE environments. Global Change Biol. 26, 539–556. doi: 10.1111/gcb.14830
Chaturvedi, A. K., Bahuguna, R. N., Pal, M., Shah, D., Maurya, S., Jagadish, K. S. V. (2017). Elevated CO2 and heat stress interactions affect grain yield, quality and mineral nutrient composition in rice under field conditions. Field Crops Res. 206, 149–157. doi: 10.1016/j.fcr.2017.02.018
Cui, J. L., Zhang, X. M., Reis, S., Wang, C., Wang, S. T., He, P. Y., et al. (2023). Nitrogen cycles in global croplands altered by elevated CO2. Nat. Sustainability 6, 1166–1176. doi: 10.1038/s41893-023-01154-0
Dingkuhn, M., Luquet, D., Fabre, F., Muller, B., Yin, X. Y., Paul, M. J. (2020). The case for improving crop carbon sink strength or plasticity for a CO2-rich future. Curr. Opin. Plant Biol. 56, 259–272. doi: 10.1016/j.pbi.2020.05.012
Drissner, D., Blum, H., Tscherko, D., Kandeler, E. (2007). Nine years of enriched CO2 changes the function and structural diversity of soil microorganisms in a grassland. Eur. J. Soil Sci. 58, 260–269. doi: 10.1111/j.1365-2389.2006.00838.x
Ellsworth, D. S., Crous, K. Y., Kauwe, M. G. D., Verryckt, L. T., Goll, D., Zaehle, S., et al. (2022). Convergence in phosphorus constraints to photosynthesis in forests around the world. Nat. Commun. 13, 1–12. doi: 10.1038/s41467-022-32545-0
Fabre, D., Dingkuhn, M., Yin, X. Y., Clément-Vidal, A., Roques, S., Soutiras, A., et al. (2020). Genotypic variation in source and sink traits affects the response of photosynthesis and growth to elevated atmospheric CO2. Plant Cell Environ. 43, 579–593. doi: 10.1111/pce.13693
Fan, X. R., Tang, Z., Tan, Y. W., Zhang, Y., Luo, B. B., Yang, M., et al. (2016). Overexpression of a pH-sensitive nitrate transporter in rice increases crop yields. Proc. Natl. Acad. Sci. United States America 113, 7118–7123. doi: 10.1073/pnas.1525184113
Gao, B., Hu, S. W., Jing, L. Q., Niu, X. C., Wang, Y. X., Zhu, J. G., et al. (2021b). Alterations in source-sink relations affect rice yield response to elevated CO2: A free-air CO2 enrichment study. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.700159
Gao, B., Hu, S. W., Jing, L. Q., Wang, Y. X., Zhu, J. G., Li, H. Y., et al. (2021a). Impact of elevated CO2 and reducing the source-sink ratio by partial defoliation on rice grain quality – A 3-year free-air CO2 enrichment study. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.788104
Gojon, A., Cassan, O., Bach, L., Lejay, L., Martin, A. (2023). The decline of plant mineral nutrition under rising CO2: physiological and molecular aspects of a bad deal. Trends Plant Sci. 28, 185–198. doi: 10.1016/j.tplants.2022.09.002
Guo, J., Zhang, M. Q., Wang, X. W., Zhang, W. J. (2015). A possible mechanism of mineral responses to elevated atmospheric CO2 in rice grains. J. Integr. Agric. 14, 50–57. doi: 10.1016/S2095-3119(14)60846-7
Hu, S. W., Chen, W., Tong, K. C., Wang, Y. X., Jing, L. Q., Wang, Y. L., et al. (2022a). Response of rice growth and leaf physiology to elevated CO2 concentrations: A meta-analysis of 20-year FACE studies. Sci. Total Environ. 807, 151017. doi: 10.1016/j.scitotenv.2021.151017
Hu, S. W., Li, T., Wang, Y. X., Gao, B., Jing, L. Q., Zhu, J. G., et al. (2024). Effects of free air CO2 enrichment (FACE) on grain yield and quality of hybrid rice. Field Crops Res. 306, 109237. doi: 10.1016/j.fcr.2023.109237
Hu, S. W., Tong, K. C., Chen, W., Wang, Y. X., Wang, Y. L., Yang, L. X. (2022b). Response of rice grain quality to elevated atmospheric CO2 concentration: A meta-analysis of 20-year FACE studies. Field Crops Res. 284, 108562. doi: 10.1016/j.fcr.2022.108562
Hu, S. W., Wang, Y. X., Yang, L. X. (2021). Response of rice yield traits to elevated atmospheric CO2 concentration and its interaction with cultivar, nitrogen application rate and temperature: A meta-analysis of 20 years FACE studies. Sci. Total Environ. 764, 142797. doi: 10.1016/j.scitotenv.2020.142797
IPCC. (2021). Climate change 2021: the physical science basis. Contribution of working group I to the sixth assessment report of the intergovernmental panel on climate change. Eds. Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S. (Cambridge, United Kingdom & New York, NY, USA: Cambridge University Press), 81. doi: 10.1017/9781009157896
Jiang, Q., Zhang, J. S., Xu, X., Liu, G., Zhu, J. G. (2020b). Effects of free-air CO2 enrichment (FACE) and nitrogen (N) supply on N uptake and utilization of indica and japonica cultivars (Oryza sativa L.). Ecol. Processes 9, 35. doi: 10.1186/s13717-020-00238-5
Jiang, Q., Zhu, J. G., Zhu, C. W., Liu, G., Zhang, J. S., Xu, X. (2020a). Effects of free-air CO2 enrichment(FACE) on mineral element partitioning and rhizosphere availability of rice (Oryza sativa L.). Soils 52, 552–560. doi: 10.13758/j.cnki.tr.2020.03.019
Jin, J., Tang, C. X., Robertson, A., Franks, A. E., Armstrong, R., Sale, P. (2014). Increased microbial activity contributes to phosphorus immobilization in the rhizosphere of wheat under elevated CO2. Soil Biol. Biochem. 75, 292–299. doi: 10.1016/j.soilbio.2014.04.019
Kim, H. Y., Lieffering, M., Kobayashi, K., Okada, M., Miura, S. (2003). Seasonal changes in the effects of elevated CO2 on rice at three levels of nitrogen supply: A free air CO2 enrichment (FACE) experiment. Global Change Biol. 9, 826–837. doi: 10.1046/j.1365-2486.2003.00641.x
Lawlor, D. W., Paul, M. J. (2014). Source/sink interactions underpin crop yield: the case for trehalose 6-phosphate/SnRK1 in improvement of wheat. Front. Plant Sci. 4. doi: 10.3389/fpls.2014.00418
Li, C. H., Zhu, J. G., Zeng, Q., Liu, G. (2020). Changes in microelement availability in a paddy field exposed to long-term atmospheric CO2 enrichment. J. Soils Sediments 20, 2439–2445. doi: 10.1007/s11368-020-02601-7
Li, M., Jia, Y. X., Xia, X. Y., Wang, B., Song, C. Y., Zhu, B., et al. (2023). Effects of elevated atmospheric CO2 concentration and warming on the content of iron, zinc, and phytic acid in double cropping rice grains. J. Agro-Environment Sci. 42, 1195–1207. doi: 10.11654/jaes.2022-0758
Li, P., Han, X., Zong, Y. Z., Li, H. Y., Lin, E., Han, Y. H., et al. (2015). Effects of free-air CO2 enrichment (FACE) on the uptake and utilization of N, P and K in Vigna radiata. Agriculture Ecosyst. Environ. 202, 120–125. doi: 10.1016/j.agee.2015.01.004
Lieffering, M., Kim, H. Y., Kobayashi, K., Okada, M. (2004). The impact of elevated CO2 on the elemental concentrations of field-grown rice grains. Field Crops Res. 88, 279–286. doi: 10.1016/j.fcr.2004.01.004
Loladze, I. (2014). Hidden shift of the ionome of plants exposed to elevated CO2 depletes minerals at the base of human nutrition. eLife 3, e02245. doi: 10.7554/eLife.02245.001
Long, S. P., Ainsworth, E. A., Leakey, A. D. B., Nösberger, J., Ort, D. R. (2006). Food for thought: lower-than-expected crop yield stimulation with rising CO2 concentrations. Science 312, 1918–1921. doi: 10.1126/science.1114722
Lv, C. H., Huang, Y., Sun, W. J., Zhu, J. G. (2020). Response of rice yield and yield components to elevated [CO2]: A synthesis of updated data from FACE experiments. Eur. J. Agron. 112, 125961. doi: 10.1016/j.eja.2019.125961
Myers, S. S., Zanobetti, A., Kloog, I., Huybers, P., Leakey, A. D. B., Bloom, A. J., et al. (2014). Increasing CO2 threatens human nutrition. Nature 510, 139–142. doi: 10.1038/nature13179
NOAA. (2024). Trends in CO2. Available online at: https://www.esrl.noaa.gov/gmd/ccgg/trends/ (Accessed June 05, 2024).
Pang, J., Zhu, J. G., Xie, Z. B., Chen, G. P., Liu, G., Zhang, Y. L. (2005). Effects of elevated pCO2 on nutrient uptake by rice and nutrient contents in rice grain. Chin. J. Rice Sci. 19, 350–354. doi: 10.16819/j.1001-7216.2005.04.012
Pang, J., Zhu, J. G., Xie, Z. B., Liu, G., Zhang, Y. L., Chen, G. P., et al. (2006). A new explanation of the N concentration decrease in tissues of rice (Oryza sativa L.) exposed to elevated atmospheric pCO2. Environ. Exp. Bot. 57, 98–105. doi: 10.1016/j.envexpbot.2005.04.004
Paul, M. J., Watson, A., Griffiths, C. A. (2020). Linking fundamental science to crop improvement through understanding source and sink traits and their integration for yield enhancement. J. Exp. Bot. 71, 2270–2280. doi: 10.1093/jxb/erz480
Pozo, A. D., Pérez, P., Gutiérrez, D., Alonso, A., Morcuende, R., Martínez-Carrasco, R. (2007). Gas exchange acclimation to elevated CO2 in upper-sunlit and lower-shaded canopy leaves in relation to nitrogen acquisition and partitioning in wheat grown in field chambers. Environ. Exp. Bot. 59, 371–380. doi: 10.1016/j.envexpbot.2006.04.009
Ren, S. R., Zhu, J. G., LI, H. X., Xie, Z. B., Zeng, Q., Wang, X. Z. (2007). Effect of free-air CO2 enrichment (FACE) on miner al elements of xylem sap of rice. J. Agro-Environment Sci. 26, 1849–1853. doi: 10.3321/j.issn:1672-2043.2007.05.050
Shimono, H., Okada, M., Yamakawa, Y., Nakamura, H., Kobayashi, K., Hasegawa, T. (2009). Genotypic variation in rice yield enhancement by elevated CO2 relates to growth before heading, and not to maturity group. J. Exp. Bot. 60, 523–532. doi: 10.1093/jxb/ern288
Sonnewald, U., Fernie, A. R. (2018). Next-generation strategies for understanding and influencing source-sink relations in crop plants. Curr. Opin. Plant Biol. 43, 63–70. doi: 10.1016/j.pbi.2018.01.004
Sperotto, R. A., Ricachenevsky, F. K., Waldow, V., de, A., Müller, A. L. H., Dressler, V. L., et al. (2013). Rice grain Fe, Mn and Zn accumulation: How important are flag leaves and seed number? Plant Soil Environ. 59, 262–266. doi: 10.17221/841/2012-PSE
Uddling, J., Broberg, M. C., Feng, Z. Z., Pleijel, H. (2018). Crop quality under rising atmospheric CO2. Curr. Opin. Plant Biol. 45, 262–267. doi: 10.1016/j.pbi.2018.06.001
Ujiie, K., Ishimaru, K., Hirotsu, N., Nagasaka, S., Miyakoshi, Y., Ota, M., et al. (2019). How elevated CO2 affects our nutrition in rice, and how we can deal with it. PloS One 14, e0212840. doi: 10.1371/journal.pone.0212840
Wang, X. X., Cai, C., Song, L., Zhou, W., Yang, X., Gu, X. Y., et al. (2024). Responses of rice grain yield and quality to factorial combinations of ambient and elevated CO2 and temperature in T-FACE environments. Field Crops Res. 309, 109328. doi: 10.1016/j.fcr.2024.109328
Wang, Y., Huang, Y. Y., Song, L., Yuan, J. H., Li, W., Zhu, Y. G., et al. (2023). Reduced phosphorus availability in paddy soils under atmospheric CO2 enrichment. Nat. Geosci. 16, 162–168. doi: 10.1038/s41561-022-01105-y
Wang, Y. X., Song, Q. L., Frei, M., Shao, Z. S., Yang, L. X. (2014). Effects of elevated ozone,carbon dioxide,and the combination of bothon the grain quality of Chinese hybrid rice. Environ. pollut. 189, 9–17. doi: 10.1016/j.envpol.2014.02.016
Wei, H. H., Meng, T. Y., Li, X. Y., Dai, Q. G., Zhang, H. C., Yin, X. Y. (2018). Sink-source relationship during rice grainfilling is associated with grain nitrogen concentration. Field Crops Res. 215, 23–38. doi: 10.1016/j.fcr.2017.09.029
Wu, C. Y., Lu, L. L., Yang, X. E., Feng, Y., Wei, Y. Y., Hao, H. L., et al. (2010). Uptake, translocation, and remobilization of zinc absorbed at different growth stages by rice genotypes of different Zn densities. J. Agric. Food Chem. 58, 6767–6773. doi: 10.1021/jf100017e
Xu, B. B., You, C. C., Ding, Y. F., Wang, S. H. (2016). Effect of source-sink manipulation on translocation of carbohydrate and nitrogen, phosphors, potassium in vegetative organs of conventional Japonica rice after anthesis. Scientia Agricultura Sin. 49, 643–656. doi: 10.3864/j.issn.0578-1752.2016.04.004
Yang, L. X., Wang, Y. L., Huang, J. Y., Zhu, J. G., Yang, H. J., Liu, G., et al. (2007). Seasonal changes in the effects of free-air CO2 enrichment (FACE) on phosphorus uptake and utilization of rice at three levels of nitrogen fertilization. Field Crops Res. 102, 141–150. doi: 10.1016/j.fcr.2007.03.004
Yang, L. X., Wang, Y. L., Kobayashi, K., Zhu, J. G., Huang, J. Y., Yang, H. J., et al. (2008). Seasonal changes in the effects of free-air CO2 enrichment (FACE) on growth, morphology and physiology of rice root at three levels of nitrogen fertilization. Global Change Biol. 14, 1844–1853. doi: 10.1111/j.1365-2486.2008.01624.x
Yin, X. Y., Gu, J. F., Dingkuhn, M., Struik, P. C. (2022). A model-guided holistic review of exploiting natural variation of photosynthesis traits in crop improvement. J. Exp. Bot. 73, 3173–3188. doi: 10.1093/jxb/erac109
Yoshida, S., Forno, D. A., Cock, J. H., Gomez, K. A. (1976). Laboratory manual for physiological studies of rice. 3rd ed.; (Laguna, Philippines: The International Rice Research Institute), 46–49. Available at: http://books.irri.org/9711040352_content.pdf.
Yuan, L., Wu, L. H., Yang, C. L., Lv, Q. (2012). Effects of iron and zinc foliar applications on rice plants and their grain accumulation and grain nutritional quality. J. Science Food Agric. 93, 254–261. doi: 10.1002/jsfa.5749
Zhang, G. Y., Sakai, H., Tokida, T., Usui, Y., Zhu, C. W., Nakamura, H., et al. (2013). The effects of free-air CO2 enrichment (FACE) on carbon and nitrogen accumulation in grains of rice (Oryza sativa L.). J. Exp. Bot. 64, 3179–3188. doi: 10.1093/jxb/ert154
Zhang, J. Y., Li, Y. S., Yu, Z. H., Adams, J., Tang, C. X., Wang, G. H., et al. (2022). Elevated atmospheric CO2 and warming enhance the acquisition of soil-derived nitrogen rather than urea fertilizer by rice cultivars. Agric. For. Meteorology 324, 109117. doi: 10.1016/j.agrformet.2022.109117
Zhou, J., Shu, X. W., Xu, G. P., Lai, S. K., Yang, L. X., Dong, G. C., et al. (2021). Effect of elevated CO2 on potassium absorption and utilization in different rice varieties. J. Agro-Environment Sci. 40, 736–746. doi: 10.11654/jaes.2020-0980
Zhu, C. W., Kobayashi, K., Loladze, I., Zhu, J. G., Jiang, Q., Xu, X., et al. (2018). Carbon dioxide (CO2) levels this century will alter the protein, micronutrients, and vitamin content of rice grains with potential health consequences for the poorest rice-dependent countries. Sci. Adv. 4, eaaq1012. doi: 10.1126/sciadv.aaq1012
Zhu, C. W., Xu, X., Wang, D., Zhu, J. G., Liu, G., Seneweera, S. (2016). Elevated atmospheric CO2 stimulates sugar accumulation and cellulose degradation rates of rice straw. GCB Bioenergy 8, 579–587. doi: 10.1111/gcbb.12277
Keywords: free-air CO2 enrichment, rice, source and sink, mineral elements, absorption and distribution
Citation: Gao B, Hu S, Zhou M, Jing L, Wang Y, Zhu J, Sun X, Wang K, Wang Y and Yang L (2024) Impacts of elevated CO2 and partial defoliation on mineral element composition in rice. Front. Plant Sci. 15:1450893. doi: 10.3389/fpls.2024.1450893
Received: 18 June 2024; Accepted: 30 October 2024;
Published: 18 November 2024.
Edited by:
M. J. I. Shohag, University of Florida, United StatesReviewed by:
Tingting Chen, Chinese Academy of Agricultural Sciences, ChinaYansheng Li, Chinese Academy of Sciences (CAS), China
Copyright © 2024 Gao, Hu, Zhou, Jing, Wang, Zhu, Sun, Wang, Wang and Yang. 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: Lianxin Yang, lxyang@yzu.edu.cn