Jurka Lesjak
Daniel F. Calderini- 1Graduate School, Faculty of Agricultural Sciences, Universidad Austral de Chile, Valdivia, Chile
- 2Institute of Plant Production and Protection, Universidad Austral de Chile, Valdivia, Chile
Quinoa high nutritive value increases interest worldwide, especially as a crop that could potentially feature in different cropping systems, however, climate change, particularly rising temperatures, challenges this and other crop species. Currently, only limited knowledge exists regarding the grain yield and other key traits response to higher temperatures of this crop, especially to increased night temperatures. In this context, the main objective of this study was to evaluate the effect of increased night temperature on quinoa yield, grain number, individual grain weight and processes involved in crop growth under the environmental conditions (control treatment) and night thermal increase at two phases: flowering (T1) and grain filling (T2) in southern Chile. A commercial genotype, Regalona, and a quinoa accession (Cod. BO5, N°191, grain bank from Semillas Baer, hereby referred to as Accession) were used, due to their adaptability to Southern Chilean conditions and contrasting grain yield potential, grain weight and size of plants. Temperature was increased ≈4°C above the ambient from 8 pm until 9 am the next morning. Control treatments reached a high grain yield (600 and 397 g m-2, i.e., Regalona and Accession). Temperature increase reduced grain yield by 31% under T1 treatment and 12% when under T2 in Regalona and 23 and 26% in Accession, respectively. Aboveground biomass was negatively affected by the thermal treatments and a positive linear association was found between grain yield and aboveground biomass across treatments. By contrast, the harvest index was unaffected either by genotype, or by thermal treatments. Grain number was significantly affected between treatments and this key trait was linearly associated with grain yield. On the other hand, grain weight showed a narrow range of variation across treatments. Additionally, leaf area index was not affected, but significant differences were found in SPAD values at the end of T1 treatment, compared to control. Little change was found in the harvest index, individual grain weight, grain protein content or water soluble carbohydrates in response to the increased night temperature in this crop.
Introduction
Quinoa (Chenopodium quinoa Willd.), called “the mother grain” by the natives of the Andean highlands of Bolivia, Peru, Chile, and Ecuador (Cusack, 1984), has reached global recognition as a “superfood” for the future (Shearman, 2014). Quinoa production has almost tripled in the last two decades due largely to prices rising in response to the growing demand (FAOSTAT, 2014). This pseudo-cereal gained global attention due to its high nutritive value, associated with high protein quality (Ruales and Nair, 1992). Quinoa contains all the essential amino acids, trace elements and vitamins and is gluten-free (FAO, 2013) which makes it an excellent crop to face the challenge of increasing quantity and quality of food production globally (FAO, 2013).
Global warming is a central challenge for grain crops such as quinoa. It is expected that world average temperatures will increase between 1 and 3.7°C by the end of this century. Within these processes daily minimum temperatures are projected to increase faster than daily maximum ones (IPCC, 2014). It is noteworthy that significant increases in annual night average temperatures were recorded between 1979 and 2003 (Alexander et al., 2006), particularly in South America (IPCC, 2007; Marengo et al., 2010). Expected temperature increases in the Southern Hemisphere are foreseen for the end of spring and during the summer, from November to February (Parry et al., 2007), with the highest daily maximum temperature increase of 1–4°C (IPCC, 2014). Given that temperature is an important environmental factor affecting plant development, growth and yield (Hunt et al., 2001), it is expected that warmer climate conditions will accelerate plant development and consequently reduce the length of crop phenophases, decreasing yield (Fuhrer, 2003; Ortiz, 2008; Hatfield and Prueger, 2015) as it was demonstrated in wheat (García et al., 2015, 2016). Projections, on how and to what extent temperature increases will affect staple food crops worldwide, are crucial to future food security and cropping system adaptations. The impact of higher temperature on grain yield and associated traits has so far been assessed by simulation and field experiments for variety of crops: e.g., wheat (Asseng et al., 2009; Lizana and Calderini, 2013; García et al., 2015); sunflower (Rondanini et al., 2006; Van der Merwe, 2010); rice (Peng et al., 2004) and canola (Gan et al., 2004); however, the response of quinoa to augmented night temperatures is still unknown.
In Chile, quinoa is associated with diverse production areas as it is cultivated in a wide diversity of ecological zones. One of the two highest quinoa-yielding regions is southern Chile (Bazile, 2013), recording an average yield of 1.9 ton ha-1 (von Baer et al., 2009). This area (from 37 to 41°S) is also one of the most promising in the context of high yields of various crops, mainly due to favorable temperature conditions (Bustos et al., 2013; Mera et al., 2014). In the context of expected temperature increases in the Southern Hemisphere presented above, a temperature increase might therefore affect the flowering and grain filling periods of quinoa in the southern area of Chile, decreasing yield potential.
Wide genetic variability and plasticity which allow adaptation to different environments are the reasons why quinoa can grow under extreme climate conditions and poor soils (Jacobsen, 2003; Danielsen and Ames, 2004; Fuentes et al., 2009). The ability of this crop to cope with heat stress, however, has not yet been extensively researched. Mujica et al. (2001) reported that quinoa is cultivated between 15 and 20°C in the Andean environment and GxE interaction was demonstrated by Bertero et al. (2004) for quinoa grain yield and grain size through a temperature range between 9 and 22°C across three continents. It has also been suggested that Chilean sea level cultivars are less sensitive to the combination of higher temperatures and long day photoperiod, which may explain its wide adaptation (Bertero et al., 1999; Jacobsen, 2003). Nevertheless, Fuentes et al. (2012) suggested that the grain yield of quinoa in the Atacama Desert is low due to the negative effects of high-temperature stress around flowering (ranging from 31 to 35°C). To the best of our knowledge, no experimental evaluation has been performed to assess these assumptions, highlighting the need to evaluate the grain yield response of quinoa to projected temperatures.
Grain yield integrates two main components, grain number per m2 and average grain weight, where yield variations in grain crop species are usually associated with grain number (Peltonen-Sainio et al., 2007; Sadras, 2007; Araus et al., 2008; Gambín and Borrás, 2010). Crop resources are distributed between these two components during the grain set period resulting in a trade-off between them (Gambín and Borrás, 2010), with differences in assimilate availability per grain, i.e., a measure of the source size for the growing grains (Gambín et al., 2006). A strong positive relationship between yield and grain number has also been found in quinoa, whereas grain weight has shown a conservative behavior (Mignone and Bertero, 2007). Recently, Bertero and Ruiz (2008) reported that grain number in quinoa is highly sensitive to environmental factors from first anthesis to the end of flowering, followed by the grain filling period. Differences in grain number can also be attributed to the partitioning of biomass into the reproductive structures during the critical phase (D’Andrea et al., 2006). An improved understanding of yield responses to higher temperature during different phenological phases could uncover important clues about the sensitivity of the grain yield of quinoa to global warming.
As discussed above, yield potential of quinoa is very promising in southern Chile. Taking into account the scarcity of data on the response of this crop to environmental constrains as global warming, the objective of this study was to evaluate the impact of increased temperature on yield, grain number, grain weight and associated traits of two Chilean quinoas under the high-yielding environment of southern Chile.
Materials and Methods
Site Location, Experimental Design, and Plant Material
An experiment was carried out at the Experimental Agricultural Station of the Universidad Austral de Chile, in Valdivia (39° 47′S, 73° 14′W, 19 m asl) during the 2012/13 growing season. Treatments consisted of two genotypes and three thermal conditions; i.e., control, heating around anthesis (T1) and heating around grain filling (T2). These treatments were arranged in a 2 × 3 factorial combination in a Completely Randomized Design with three replications (i.e., 18 experimental units).
In order to achieve the proposed objective, a commercial variety, Regalona (Semillas Baer Temuco, Chile), was sown, as that this variety is the topmost sown quinoa in the south of Chile (virtually the only one in this area of the country). Also, an accession (Cod. BO5, N°191, grain bank from Semillas Baer), was also included in the experiment. This was done with the aim to increase our understanding of the response of the Chilean sea level genotypes to higher temperature, since this genotype is a source of germplasm for the quinoa breeding program. In addition, Regalona has presented a higher grain weight than Accession but similar phenology in previous evaluations. Thermal treatments consisted of: control at ambient temperature, and two independent increased temperature conditions, i.e., at 12 days after the first anthesis (T1) and at 12 days after the end of flowering, i.e., during the grain filling period (T2) (Supplementary Figure A1). Both treatments were applied for a 12-day period as in previous study of wheat in southern Chile (Lizana and Calderini, 2013). Temperature was increased at night, by approximately 4°C above the ambient.
Sowing and Management Conditions
The experiment was sown on October 30, 2012 in a Typic Hapludand soil, developed from volcanic ashes, pH 5.8 and organic matter by 15% (for more details see Valle et al., 2009; Dörner et al., 2015). Seedling emergence was recorded on November 13, 2012. The experimental units (plots), oriented north–south, were 2.8 m long and 4.2 m wide, 35 cm between rows (11 rows per plot) and 20 cm between plants, with 13 plants per row at a plant rate of 12.1 plants m-2. To avoid aluminum and biotic constrains, 4kg m-2 of CaO and 50g m-2 of fumigant (DAZOMET) were applied in the soil before sowing. Plots were kept free of fungal diseases and insect pests through the entire growing season as recommended by Semillas Baer Co. (Chile). Weeds were removed by hand. To avoid water stress, plots were regularly watered throughout the growing season to complement seasonal rainfall. Experimental units were fertilized before sowing with 20 kg P ha-1 and 18 kg N ha-1 and two fertilizations with urea totaling 200 kg N ha-1, were applied at 30 and 60 days after sowing.
The thermal treatment was performed using transparent polyethylene (100 μm thick) chambers, of 2.8 × 2.8 and 1.5 m (length, width, and height) equipped with electric heaters controlled by thermostats as in prior experiments (Savin et al., 1996; Ugarte et al., 2007; Lizana and Calderini, 2013). The tops of the chambers were set at night (8 pm) and removed in the morning (9 am). This was done to increase the night temperature only, avoiding a daily maximum temperature above 30°C (heat shock) in the chambers and to prevent reductions in solar radiation by the polyethylene film. Heating was controlled by thermal sensors placed at panicle height, connected to a temperature regulator as in Lizana and Calderini (2013). Temperature was recorded, both in the chambers and outside, by data loggers (Cavadevices, Buenos Aires, Argentina). Mean air temperature in the field was recorded by the weather station set (Davis Vantage Pro, USA) positioned 50 m from the experiment. Average day temperature during the treatments was calculated by using start temperature of the treatments (from 8 pm until 9 am the next morning) and ambient temperatures.
Measurements
Development and Biomass
Developmental stages were recorded as: emergence, visible flower bud, first anthesis (at least one flower open on 50% of the plants in the main inflorescence from each plot), end of anthesis (no more flowers open on at least 50% of the plants from each plot), and physiological maturity.
These stages also defined four developmental phases: vegetative, reproductive, flowering and grain filling (Bertero and Ruiz, 2008). Grain filling duration was calculated by subtracting the days up to anthesis from the time to physiological maturity.
Biomass was sampled at each plot by harvesting five contiguous plants in central rows, not including border plants. In treatment T1, samplings were made at: (i) 14 days after visible floral bud, (ii) the beginning of the heating treatment (12 days after first anthesis), (iii) 6 days after starting the treatment, (iv) the end of the treatment, (v) 23 days after the treatment and (vi) physiological maturity for both genotypes. Therefore, samples for the T1 treatment were obtained at 59 (the beginning of the treatment), 65, 68, 71, 83, 94, 100, and 128 (physiological maturity) days after emergence (DAE). In the T2 treatment, samplings were made at (i) 14 days after visible floral bud, (ii) the end of anthesis (beginning of the grain filling period), (iii) 12 days after the beginning of the grain filling period (beginning of the treatment), (iv) 6 days after sampling iii, (v) the end of the treatment, and (vi) physiological maturity. Thus in this treatment samples were harvested at 76 (the beginning of grain filling), 88 (beginning of the treatment), 91, 94, 97, 100, 115, and 128 DAE. Plants from the corresponding control plots were similarly sampled. At harvest maturity, 10 contiguous plants were sampled per plot.
Plant biomass was divided into: green leaves, senescent leaves, stem (main stem and branches) and inflorescences (main stem and branches) when present. To calculate leaf area index (LAI), leaf area of green leaves was measured by an electronic leaf area meter (LI 3100, Licor Inc., Lincoln, NE, USA) and accordingly calculated. Grain weight was measured in each plant; from which 3 × 50 grains were taken from the middle position of the panicle. Grains were sampled with petal leaves to weigh fresh biomass. Once dried in an air-forced drying oven at 60°C for 72 h, grains were weighed to record dry biomass. Grain number per m-2 was calculated as the ratio between grain yield and average grain weight at harvest of each plot. Aboveground biomass was calculated as the sum of all organs. Grain yield of the main stem and branches were determined at the end of the growth cycle separately, whereas total yield was considered their sum. The harvest index (HI) of each experimental unit was calculated as the ratio between grain yield and aboveground biomass at maturity. Crop growth rate (CGR) was calculated as the slope of the linear regression between biomass and time (in days).
Chlorophyll Content
A SPAD portable chlorophyll meter (SPAD-502 leaf chlorophyll meter, Minolta Corporation, Osaka, Japan) was used to measure the light transmission spectrum of chlorophyll a and b at two leaf positions: (i) the top third and (ii) from the middle third of the main stem, at each evaluation (five plants per plot) in thermal treatment. The measured values (arbitrary units) showed a linear relationship with total chlorophyll content in quinoa (Bertero, 2001).
Water-Soluble Carbohydrates (WSC)
Water-soluble carbohydrates (glucose, fructose, and sucrose) were extracted from the main stem and branches at the beginning and end of the thermal treatments. Five plants per plot were sampled for each evaluation. WSC was measured according to the Analysis of Agricultural Materials (A.O.A.C., 1996. Official Method; Methods 14, Carbohydrates, Soluble, in Herbage).
Grain Protein Content
Grain protein concentration (%) at harvest maturity was measured from the bulk grain at each replicate, using the Kjeldhal method and a conversion factor of 6.25, i.e., equivalent to 0.16 g nitrogen per gram of protein (Merrill and Watt, 1973).
Data Analysis
The effects of increased temperature treatments were analyzed using a factorial analysis of variance (ANOVA) with two genotypes and three temperature regimes. An ANOVA was performed using the STATGRAPHICS® Centurion XV software, (StatPoint Inc., 2005). Multiple means were compared with Fisher’s Least Significant Difference (LSD) test, where the probability level of p ≤ 0.05 was considered statistically significant. Regression analyses were used to measure the association between variables.
Results
Environmental Conditions and Crop Phenology
Environmental conditions, other than temperature during the thermal treatments, were similar for both genotypes across the growing season. Thermal treatments increased the daily average temperature in the chambers by 2.1°C above the ambient temperature for the 12 days of treatment across genotypes and phenophases (Table 1 and Supplementary Figure A1). With regard to night temperatures (between 8 pm and 9 am), heating increased temperature in the chambers by 4.3°C (T1) and 4.2°C (T2), across the 12 days treatment in Regalona, and by 3.8°C (T1) and 3.7°C (T2) in Accession (Table 1 and Supplementary Figure A1).
TABLE 1. Mean temperature of control (from emergence to physiological maturity at ambient temperature) and increased temperature (T1 and T2) treatments in Regalona and Accession.
The average temperature for cv. Regalona during the crop cycle (from emergence to physiological maturity) was 15.9°C at the control treatment, and 16.1°C under thermal increase (under both T1 and T2). The average temperature for Accession between emergence and maturity was 15.6°C in the control plants, 15.9°C under T1 and 15.8°C under the T2 treatments (Supplementary Figure A1). Average incident PAR, between emergence and physiological maturity, was 11 MJ m-2 day-1 for Regalona and 10 MJ m-2 day-1 for Accession (Supplementary Figure A1). The crop cycle was divided into four phases: vegetative, reproductive, flowering and grain filling, following Bertero and Ruiz (2008). Regalona reached visible flower at 25 days after emergence, while Accession showed a longer vegetative phase, i.e., 35 days across treatments (Figure 1). Increased temperature negligibly shortened the growing cycle in both genotypes with no significant difference (P > 0.05). Regalona growth cycle was shortened by 3 days under both temperature treatments (T1, T2), while Accession accelerated the growth cycle by 4 days under T1 and by 2 days under T2 (Figure 1).
FIGURE 1. Phenological phases duration from emergence in control and increased temperature treatments (T1 and T2) of Regalona and Accession. DAE, days after emergence.
Effect of Increased Temperature Treatments on Grain Yield, Aboveground Biomass, and Yield Components at Harvest
Grain yield was affected (P < 0.001) by both genotype and thermal treatments. In the controls, Regalona exceeded Accession yield by 51% (Table 2). Increased temperature at flowering reduced grain yield by 31% in Regalona and 23% in Accession. When temperature was increased during grain filling, grain yield of Regalona was not significantly affected by heating, though reduced by 13%, while Accession showed similar sensitivity (-26%) to that at flowering (Table 2). Like grain yield, aboveground biomass was also affected (P < 0.001) by genotype and thermal treatments. The sensitivity of biomass to increased temperature at flowering (T1) was higher in Regalona (29%) than in Accession (17%), while during grain filling (T2) biomass decreased by 15 and 21% in Regalona and Accession, respectively (Table 2). A positive linear association was found between grain yield and aboveground biomass across treatments when they were plotted together (Figure 2). By contrast, the HI was not affected by either genotype or increased temperature treatments (Table 2), and no association was found between grain yield and HI. Plant height differed between genotype and treatment (p < 0.05), where the shorter plants were found in the more productive cultivar, Regalona (Table 2).
TABLE 2. Yield and yield components at harvest of Regalona and Accession under temperature treatments (Control: C; Anthesis: T1 and Grain Filling: T2).
FIGURE 2. Relationship between grain yield and above-ground biomass (A), grain number (B), and plant height (C). Bars show the standard error (SE) of the means.
With regard to yield components, significant differences were found in grain number per square meter between genotypes and among thermal treatments (Table 2). In Regalona, grain number was significantly (P < 0.05) affected by high temperatures only under T1 (27% reduction); however, in Accession grain number was sensitive to heating under both T1 (20%) and T2 (22%). Interestingly, a positive association (P = 0.01) was found between grain yield and grain number across the treatments (Figure 2).
As expected, grain weight significantly differed between the genotypes (Table 2). Despite grain weight being reduced by 5 and 4% due to the thermal treatments in flowering and grain filling, these reductions were not significant (Table 2). In addition, no interaction between genotype x thermal treatments was observed for this trait. A linear relationship was found between yield and grain weight (Figure 2), but this association was mainly due to differences between genotypes. Therefore, changes in grain yield were explained by grain number.
Biomass Partitioning in Main Stems and Branches
A complementary analysis was made by assessing the contributions of main stems and branches to grain yield and aboveground biomass across the treatments. In Regalona, panicles from the main stems and branches of the control plants contributed almost the same to grain yield, i.e., 294.5 and 305.5 g m-2, respectively (Figure 3). Main stems and branches of Accession showed a greater contrast than in Regalona for both grain yield and biomass partitioning in the control plants. Main stems of Accession over-yielded branches in both grain yield and aboveground biomass, but particularly in the latter (Figure 3).
FIGURE 3. Grain yield (lines) and biomass (bars) of main stem (black) and branches (gray) of control and heating treatments (T1 and T2) in Regalona (A), Accession (B). Bars show the SE of the means.
Increased temperature treatments significantly (P < 0.05) affected the grain yield of main stems, but not biomass (P > 0.05; Figure 3) in Regalona. Both traits were decreased by 23% under T1 and by 16 and 7% under T2, respectively. Grain yield and biomass of branches were both significantly affected (P < 0.05) by the thermal treatment, showing a linear relationship between them (R2 = 0.50; P = 0.03). Grain yield of branches was decreased by 38% under T1, and 17% under T2, while biomass of branches decreased by 30% under T1 and 23% under T2, compared to the controls. Grain yield and biomass of the main stems were not affected in Accession and no association was found between these traits (Figure 3). The difference in main stem grain yield was -8% than control under T1 and -17% under T2, while biomass decreased by -20% under T1 and by -8% under T2. No difference was found in the biomass of branches, while grain yield was affected (P < 0.05) and a positive association (R2= 0.85; P < 0.001) was observed between them. Additionally, the biomass of branches decreased by 30 and 32% under T1 and T2, respectively, while grain yield was more affected and decreased by 42% under T1 and 37% under T2.
Aboveground Biomass Accumulation under Thermal Treatments
Taking into account the association between grain yield and biomass shown in Figure 4A, and to obtain data on the sensitivity of biomass production at higher temperatures in quinoa, we also evaluated the time-course of aboveground biomass was evaluated in the present study. At the beginning of T1, the aboveground biomass was similar (P > 0.05) between the control plots and thermal treatments in Regalona and Accession (Figure 4), however, at the end of heating, i.e., 12 days after the beginning of T1, lower (P < 0.05) biomass was recorded in both genotypes (-17 and -13% for Regalona and Accession). Differences in biomass between T1 and the control treatments were increasing until harvest (Figure 4), when the control treatment showed significantly higher biomass (29%; P < 0.01) than T1 in Regalona and Accession (17%; P < 0.05) supporting the key hypothesis on the sensitivity of the biomass to higher temperature (Table 2).
FIGURE 4. Time-course of above-ground biomass offor control and T1 (A) and T2 (B) treatments in cv. Regalona (A,B) and Accession (C,D) days after emergence (DAE). Arrows indicate the beginning and the end of heat treatments. Bars show the SE of the means.
As under T1, aboveground biomass under T2 was similar to the control treatment at the beginning of the treatment in Regalona and Accession (Figure 4) and significantly lower biomass was found at the end of heating in Regalona (-16%) and Accession (-20%; P < 0.05) (Figure 4).
The CGR calculated from 50 DAE to maturity was significantly (P = 0.007) affected by the thermal treatments (Figure 5), whereas the effect of heating on the whole crop cycle duration, i.e., from emergence to maturity, was negligible (Figure 1). Heating decreased CGR by 23% under T1 and 6% under T2 for Regalona during the period from emergence to maturity, while this trait was not affected in Accession. Calculating CGR only during the heating treatment, i.e., 12 day-period, this rate decreased by 25% (P = 0.001) under T1 and by 63% (P = 0.03) under T2 for Regalona (Figure 6). In Accession, CGR was significantly decreased during T1, i.e., 45% (P = 0.03) and similarly decreased (43%) under T2, though not significant.
FIGURE 5. Crop growth rate (CGR) of cv. Regalona and Accession (a) under control, temperature treatments (T1 and T2) calculated for entire crop cycle.
FIGURE 6. Crop growth rate of cv. Regalona and Accession under control, temperature treatments (T1 and T2) calculated for the 12 days thermal treatment period. Bars indicate SE of the means.
Leaf area index and chlorophyll concentration (estimated by the SPAD values) were also measured. LAI was not affected under the T1 treatment, neither in Regalona or Accession (Figure 7). However, SPAD values were lower in both genotypes (P = 0.01 for Regalona and P = 0.04 for Accession), compared to control. LAI values were clearly lower under T2 (P = 0.01). At the end of T2, LAI reached 56% reduction compared to the control, mainly due to leaf senescence, which begins at the end of flowering (beginning of the grain filling period).
FIGURE 7. Time-course of LAI (leaf area index) and SPAD values under T1 treatment in cv. Regalona (A,B) and Accession (C,D) in days after emergence. Bars show the SE of the means.
Water-Soluble Carbohydrates (WSC)
As heating might affect WSC, which contributes to the size of the source for growing grains, this trait was measured by individual stem category. Similar WSC concentration was found at the beginning control and T1 treatment plots of each genotype (data not shown). Immediately after heating, the evaluation of WSC showed no differences in the WSC concentration in either main stems (P = 0.55) or branches (P = 0.47) between the control treatment (1.3 and 0.7 g/plant in main stem and branches, respectively) and T1 (1.0 and 0.6, respectively) treatment in Regalona. Similar results were found in Accession, where higher temperatures had no effect on either main stems (P = 0.72; 1.7 and 1.3 g/plant in main stems and branches, respectively) or branches (P = 0.20; 0.2 and 0.16 g/plant). Similar to T1, the WSC concentration was not affected (P > 0.05) by T2 in the main stems or branches between the control and thermal treatments in either Regalona or Accession (data no shown).
Grain Protein Concentration
Genotypes showed differences in grain protein concentration, where Accession reached higher values (Table 2). The higher temperature did not affect the grain protein concentration (Table 2). A positive association was found between the protein concentration and grain weight (P = 0.03), whereas no association (P = 0.07) between protein concentration and grain yield was found. However, it is important to highlight that the protein content-grain weight relationship was mostly caused by genotype differences in the protein concentration.
Discussion
This study evaluated the effect of increased night temperature on yield and associated traits of quinoa and addressed the gap in literature on the response of quinoa to elevated temperatures, particularly for Chilean sea level cultivars.
Increased temperature treatments were chosen in this study considering to the future climate change scenarios forecasted for southern Chile (Department of Geophysics, University of Chile, 2006; IPCC, 2014). The heated treatments consisted on ≈4°C above the ambient temperature at night for a period of 12 days, which accounted for around 2°C increase above the 24 h ambient air temperatures (Table 1). The main findings show that grain yield of quinoa is highly sensitive to increased night temperatures, especially during flowering. Taking into account this original data and its potential relevance for future research, we also compared the results were compared to evaluations carried out on other crops.
Increased temperature had a negative effect on grain yield which was explained by the adverse impact of high temperature on biomass production, what is in line with previous studies evaluating wheat (Calderini et al., 1999; Peltonen-Sainio et al., 2007; Gambín and Borrás, 2010; Sadras and Slafer, 2012). It has also been shown that chickpea is highly susceptible to heat stress during flowering and grain filling due to the negative effect of temperature on biomass production (Wang et al., 2006; Devasirvatham et al., 2012). Interestingly, the HI was insensitive (p > 0.05) to the thermal treatments in our study, even when plants were heated during grain filling.
The negative impact of high temperature on the crop biomass of quinoa could be explained by the effect on either the crop cycle duration or CGR. A great deal of evidence has shown that the rise of temperature increases the development rate which, in turn, reduces biomass production and grain yield (e.g., Fischer and Maurer, 1976). For example, previous experiments of wheat have shown that higher temperatures increase the rate of development, reducing crop radiation interception with concomitant reductions in dry matter accumulation and wheat grain yield (Farooq et al., 2011; García et al., 2015). This response is different from our study on quinoa, where only CGR was affected by temperature increase and there was no impact on the crop cycle. However, given the genetic variability shown by different traits in this specie (De Santis et al., 2016), we cannot discard genetic variability of CGR and the development rate to increased temperature in quinoa. The biggest impact found in this study was in the biomass of branches; which was considerably smaller in the main stems. As LAI was not significantly affected at the end of the thermal treatment under T1, it may be safe to suggest that heating affected the photosynthesis capacity of the plants. This would also be true for T2, as the senescence process was accelerated by the thermal treatment.
We also found that grain yield was more closely associated with grain number (R2 = 0.92; P = 0.001) than grain weight in this experiment, which is consistent with previous studies on quinoa (Bertero and Ruiz, 2008; Gómez et al., 2011) and other grain crops (Peltonen-Sainio et al., 2007; Sadras, 2007; Gambín and Borrás, 2010; Sadras and Slafer, 2012). Several studies have reported that grain yield is reduced due to higher temperatures at flowering by decreasing grain set in crop species like chickpea (Wang et al., 2006), canola (Morrison and Stewart, 2002), buckwheat (Michiyama and Sakurai, 1999; Michiyama et al., 2001), soybean (Jiang and Egli, 1995), and temperate cereals (Calderini et al., 1999; Wheeler et al., 2000; Ugarte et al., 2007; Lizana and Calderini, 2013; García et al., 2015). Regalona and Accession were highly sensitive to elevated temperatures during flowering, adjusting grain set (from main stem and branch panicles). As suggested by Bertero and Ruiz (2008), the efficiency of quinoa grain setting can be affected during grain filling, altering the final grain number by grain abortion. In soybean, abortion and abscission can occur even in immature pods, with the latter also being the most vulnerable (McBlain and Hume, 1981); by contrast, immature grains are more likely to be aborted than larger grains (Van Steveninck, 1959; Abernathy et al., 1977; Westgate and Peterson, 1993). This can be caused by different environmental stresses such as high temperatures (Duthion and Pigeaire, 1991; Nilsen and Orcutt, 1996). In addition, the quinoa commercial supply chain in Europe and North America recommends sowing quinoa in a mild climate, such as in the Loire Valley in France and the northern part of Minnesota, where summer temperatures do not exceed 23°C on average (Abbottagra, 2014; SFA, 2014). Furthermore, average day temperatures in Valdivia from November to March (almost the entire growing season) are exceedingly lower than in the Loire Valley from May to September. In the present study, the average day temperature was below 23°C even in thermal treatments, but both biomass and grain yield were negatively impacted. This suggests that higher night temperatures (by approximately ≈4°C above the ambient temperature) for a short period will negatively affect these traits. On the other hand, the temperature increase treatments had no effect on grain protein concentration.
The grain weight of quinoa was highly conservative compared to grain number in this study, which is in line with other studies (Mignone and Bertero, 2007; Bertero and Ruiz, 2008; Lopez, 2008; Spehar and Rocha, 2009; Gómez et al., 2011). Therefore, quinoa shows sensitivity to elevated temperatures only through an adjustment of the grain number, perhaps as an avoidance mechanism or as a strategy for tolerating this abiotic stress (Jacobsen et al., 2009). Nevertheless, significant changes were observed in quinoa grain weight by Munir (2011) and Gómez et al. (2013) in experiments evaluating sowing dates in contrasting environments (Faisalabad, Pakistan and Buenos Aires, Argentina).
The sensitivity of quinoa grain yield to temperature during flowering found in the present study (a reduction of 31.3% in Regalona and 23.5% in Accession) is consistent with previous manipulations of the source-sink ratio around this phenological stage (Mignone and Bertero, 2007; Bertero and Ruiz, 2008). Therefore, the importance of biomass accumulation during flowering is a key for grain number determination in quinoa.
Conclusion
Grain yield of the most sown cultivar of quinoa in southern Chile (Regalona) and a source of germplasm like Accession were highly sensitive to increased night temperature around flowering; however, this constraint had a different impact at grain filling. Thermal treatments (by approximately ≈4°C above the ambient temperature) affected biomass production and grain number only when applied during flowering in Regalona and when applied either at flowering and grain filling in Accession. In addition, grain yield was closely associated with aboveground biomass and CGR, and there was little change in the HI across treatments. Individual grain weight and grain protein content were not sensitive to either the T1 or T2 treatments.
Author Contributions
JL proposed the experiment, carried out the measurements and analyzed the data. DC proposed the hypothesis, collaborated with JL in planning the experiment and measurement and collaborated in the writing of the manuscript.
Funding
The experiment and publication of this research were partially supported by Project FONDECYT 1141048 (Chilean Technical and Scientific Research Council, CONICYT) competitive grant.
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 would like to acknowledge the assistance of Eusebio Miranda Chihuaicura from the experimental station of the Universidad Austral de Chile for his technical support, and Ingrid von Baer for providing the seeds. We also thank Helen Lowry and Jasna Magić for revising the English version of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.00352/full#supplementary-material
References
Abbottagra (2014). Available at: www.abbottagra.com/index.html
Abernathy, R. H., Palmer, R. G., Shibles, R., and Anderson, I. C. (1977). Histological observations on abscising and retaining soybean flowers. Can. J. Plant Sci. 57, 713–716. doi: 10.4141/cjps77-104
Alexander, L. V., Zhang, X., and Peterson, T. C. (2006). Global observed changes in daily climate extremes of temperature and precipitation. J. Geophys. Res. 111:D05109. doi: 10.1038/nature09763
Araus, J. L., Slafer, G. A., Royo, C., and Serret, M. D. (2008). Breeding for yield potential and stress adaptation in cereals. Crit. Rev. Plant Sci. 27, 377–412. doi: 10.1080/07352680802467736
Asseng, S., Cao, W., Zhang, W., and Ludwig, F. (2009). “Crop physiology, modelling and climate change: impact and adaptation strategies,” in Crop Physiology: Applications for Genetic Improvement and Agronomy, eds V. O. Sadras and D. F. Calderini (London: Elsevier), 511–543. doi: 10.1016/b978-0-12-374431-9.00020-7
Bazile, D. (2013). “La amplia diversidad genética de la quinua Chenopodium quinoa Willd y su expansión a nivel mundial,” in Proceedings of the IV. Congreso mundial de la quinua, 8-12 July 2013, Ibarra.
Bertero, H. D. (2001). “Variabilidad intraespecífica en variables asociadas a la generación de biomasa. Asociación con el origen de los cultivares,” in ‘Memorias Primer Taller Internacional en Quinua: Recursos Genéticos y Sistemas de Producción’, eds S. E. Jacobsen, A. Mujica, and Z. Portillo (Lima: Universidad Nacional del Altiplano), 265–272.
Bertero, H. D., De la Vega, A. J., Correa, G., Jacobsen, S. E., and Mujica, A. (2004). Genotype and genotype-by-environment interaction effects for grain yield and grain size of quinoa (Chenopodium quinoa Willd.) as revealed by pattern analysis of international multi-environment trials. Field Crops Res. 89, 299–318. doi: 10.1016/j.fcr.2004.02.006
Bertero, H. D., King, R. W., and Hall, A. J. (1999). Modelling photoperiod and temperature responses of flowering in quinoa (Chenopodium quinoa Willd.). Field Crops Res. 63, 19–34. doi: 10.1016/S0378-4290(99)00024-6
Bertero, H. D., and Ruiz, R. A. (2008). Determination of grain number in sea level quinoa (Chenopodium quinoa Willd.) cultivars. Eur. J. Agron. 28, 186–194. doi: 10.1016/j.eja.2007.07.002
Bustos, D. V., Hasan, A. K., Reynolds, M. P., and Calderini, D. F. (2013). Combining high grain number and weight through a DH-population to improve grain yield potential of wheat in high-yielding environments. Field Crops Res. 145, 106–115. doi: 10.1016/j.fcr.2013.01.015
Calderini, D. F., Abeledo, L. G., Savin, R., and Slafer, G. A. (1999). Carpel size and temperature in preanthesis modify potential grain weight in wheat. J. Agric. Sci. 132, 453–460. doi: 10.1017/S0021859699006504
D’Andrea, K. E., Otegui, M. E., Cirilo, A. G., and Eyerhabide, G. (2006). Genotypic variability in morphological and physiological traits among maize inbred lines-nitrogen responses. Crop Sci. 46, 1266–1276. doi: 10.2135/cropsci2005.07-0195
Danielsen, S., and Ames, T. (2004). Mildew (Peronospora farinosa) of Quinua (Chenopodioum quinoa) in the Andean Region. Provo, UT: Brigham Young University, 29.
De Santis, G., D’Ambrosio, T., Rinaldi, M., and Rascio, A. (2016). Heritabilities of morphological and quality traits and interrelationships with yield in quinoa (Chenopodium quinoa Willd.) genotypes in the Mediterranean environment. J. Cereal Sci. 70, 177–185. doi: 10.1016/j.jcs.2016.06.003
Devasirvatham, V., Tan, D. K. Y., Gaur, P. M., Raju, T. N., and Trethowan, R. M. (2012). High temperature tolerance in chickpea and its implications for plant improvement. Crop Pasture Sci. 63, 419–428. doi: 10.1071/CP11218
Dörner, J., Dec, D., Peng, X., and Horn, R. (2015). Effect of land use change on the dynamic behaviour of structural properties of an Andisol in southern Chile under saturated and unsaturated hydraulic conditions. Geoderma 159, 189–197. doi: 10.1016/j.geoderma.2010.07.011
Duthion, C., and Pigeaire, A. (1991). Grain lengths corresponding to the final stage in grain abortion of three grain legumes. Crop Sci. 31, 1579–1583. doi: 10.2135/cropsci1991.0011183X003100060040x
FAO (2013). Food and Agriculture Organization of the United Nations. Available at: http://www.fao.org/quinoa/es/
FAOSTAT (2014). Food and Agricultural Organization of the United Nations. The Statistics Division of the FAO. Available at: http://www.fao.org
Farooq, M., Bramley, H., Palta, J. A., and Siddique, K. H. M. (2011). Heat stress in wheat during reproductive and grain-filling phases. Crit. Rev. Plant Sci. 30, 1–17. doi: 10.1080/07352689.2011.615687
Fischer, R. A., and Maurer, R. O. (1976). Crop temperature modification and yield potential in a dwarf spring wheat. Crop Sci. 16, 855–859. doi: 10.2135/cropsci1976.0011183X001600060031x
Fuentes, F. F., Bazile, D., Bhargava, A., and Martínez, E. A. (2012). Implications of farmers’ seed exchanges for on-farm conservation of quinoa, as revealed by its genetic diversity in Chile. J. Agric. Sci. 150, 702–716. doi: 10.1017/S0021859612000056
Fuentes, F. F., Maughan, P. J., and Jellen, E. R. (2009). Diversidad genética y recursos genéticos para el mejoramiento de la quinoa (Chenopodium quinoa Willd). Revista Geográfica Valparaíso 42, 20–33.
Fuhrer, J. (2003). Agroecosystem response to combinations of elevated CO2, ozone, and global climate change. Agric. Ecosyst. Environ. 97, 1–20. doi: 10.1016/S0167-8809(03)00125-7
Gambín, B. L., and Borrás, L. (2010). Resource distribution and the trade-off between grain number and grain weight: a comparison across crop species. Anna. Appl. Biol. 156, 91–102. doi: 10.1111/j.1744-7348.2009.00367.x
Gambín, B. L., Borrás, L., and Otegui, M. E. (2006). Source-sink relations and kernel weight differences in maize temperate hybrids. Field Crops Res. 95, 316–326. doi: 10.1016/j.fcr.2005.04.002
Gan, Y., Angadi, S. V., Cutforth, H. W., Potts, D., Angadi, V. V., and Mc Donald, C. L. (2004). Canola and mustard response to short period of high temperature and water stress at different developmental stages. Can. J. Plant Sci. 84, 697–704. doi: 10.4141/P03-109
García, G. A., Dreccer, M. F., Miralles, D. J., and Serrago, R. A. (2015). High night temperatures during grain number determination reduce wheat and barley grain yield: a field study. Glob. Change Biol. 21, 4153–4164. doi: 10.1111/gcb.13009
García, G. A., Serrago, R. A., Dreccer, F. M., and Miralles, D. J. (2016). Post-anthesis warm nights reduce grain weight in field-grown wheat and barley. Field Crops Res. 195, 50–59. doi: 10.1016/j.fcr.2016.06.002
Gómez, M. B., Aguirre-Castro, P., Mignone, C., and Bertero, H. D. (2011). Can yield potential be increased by manipulation of reproductive partitioning in quinoa (Chenopodium quinoa)? Evidence from gibberellic acid synthesis inhibition using Paclobutrazol. Funct. Plant Biol. 38, 420–430. doi: 10.1071/FP10168
Gómez, M. B., Cresta, A., and Bertero, H. D. (2013). “Efecto de fechas de siembra y genotipos sobre el peso final de semillas en cultivares de quínoa adaptados a ambientes templados,” in Proceedings of the IV. Congreso Mundial de la Quinua, 8-12 Julio 2013, Ibarra.
Hatfield, J. L., and Prueger, J. H. (2015). Temperature extremes: effect on plant growth and development. Weather Clim. Extremes 10, 4–10. doi: 10.1016/j.wace.2015.08.001
Hunt, L. A., Yan, W., and McMaster, G. S. (2001). “Simulating response to temperature,” in Proceedings of the Maize and Wheat Systems Modeling Workshop CIMMYT, Symposium Proceedings, April 23-25, ed. J. W. White (El Batán: Natural Resources Group, Geographic Information Systems, Series 03-01), 23–29.
IPCC (2007). “Climate Change 2007: impacts, adaptation and vulnerability,” in Proceedings of the Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, eds M. L. Parry, O. F. Canziani, J. P. Palutikof, P. J. van der Linden, and C. E. Hanson (Cambridge: Cambridge University Press), 976.
IPCC (2014). “Central and South America,” in Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part B: Regional Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. [eds V. R. Barros, C. B. Field, D. J. Dokken, M. D. Mastrandrea, K. J. Mach, T. E. Bilir, et al. ]. Cambridge: Cambridge University Press, 1499–1566.
Jacobsen, S. E. (2003). The worldwide potential for quinoa. Food Rev. Int. 19, 167–177. doi: 10.1081/FRI-120018883
Jacobsen, S. E., Liu, F., and Jensen, C. R. (2009). Does root-sourced ABA play a role for regulation of stomata under drought in quinoa (Chenopodium quinoa Willd.). Sci. Hortic. 122, 281–287. doi: 10.1016/j.scienta.2009.05.019
Jiang, H., and Egli, D. B. (1995). Soybean grain number and crop growth rate during flowering. Agron. J. 87, 264–267. doi: 10.2134/agronj1995.00021962008700020020x
Lizana, X. C., and Calderini, D. F. (2013). Yield and grain quality of wheat in response to increased temperatures at key periods for grain number and grain weight determination: considerations for the climatic change scenarios of Chile. J. Agric. Sci. 151, 209–221. doi: 10.1017/S0021859612000639
Lopez, C. A. (2008). Efecto de la Manipulación de las Relaciones Fuente Destino al Comienzo del Llenado Sobre el Peso de Las Semillas y Variables Asociadas en Quinoa (Chenopodium quinoa Willd). Tesis der grado, Universidad de Buenos Aires, Buenos Aires.
Marengo, J. A., Rusticucci, M., Penalba, O., and Renom, M. (2010). An intercomparison of observed and simulated extreme rainfall and temperature events during the last half of the twentieth century: Part 2: historical trends. Clim. Change 98, 509–529. doi: 10.1007/s10584-009-9743-7
McBlain, B. A., and Hume, D. J. (1981). Reproductive abortion, yield components and nitrogen content in three early soybean cultivars. Can. J. Plant Sci. 61, 449–505. doi: 10.4141/cjps81-072
Mera, M., Lizana, X. C., and Calderini, D. F. (2014). “Cropping systems in environments with high yield potential of Southern Chile,” in Crop Physiology: Applications for Genetic Improvement and Agronomy, 2nd Edn, eds V. O. Sadras and D. F. Calderini (Amsterdam: Elsevier Science Publishing Co Inc.), 564.
Merrill, A. L., and Watt, B. K. (1973). Energy Value of Foods: Basis and Derivation. Agriculture Handbook No. 74. Washington, DC: ARS, US Department of Agriculture.
Michiyama, H., Arikuni, M., and Hirano, T. (2001). “Effect of air temperature on the growth, flowering and ripening in common buckwheat,” in Proceedings of the 8th International Symposium on Buckwheat (ISB) (Chunchon: Kangwon National University), 138–142.
Michiyama, H., and Sakurai, S. (1999). Effect of day and night temperatures on the growth and development of common buckwheat (Fagopyrum esculentum Moench). Jpn. J. Crop Sci. 68, 401–407. doi: 10.1626/jcs.68.401
Mignone, C. M., and Bertero, H. D. (2007). “Identificación del período crítico de determinación del rendimiento en quínoas de nivel del mar,” in Proceedings of the Congreso Internacional de la Quinua, 23-26 October 2007, Iquique.
Morrison, M. J., and Stewart, D. W. (2002). Heat stress during flowering in summer Brassica. Crop Sci. 42, 797–803. doi: 10.2135/cropsci2002.0797
Mujica, A., Jacobsen, S. E., Izquierdo, J., and Marathee, J. P. (2001). Resultados de la Prueba Americana y Europea de la Quinua. Santiago: FAO, 51.
Munir, H. (2011). Introduction and Assessment of Quinoa (Chenopodium quinoa Willd.) as a Potential Climate Proof Grain Crop. Ph.D. thesis, University of Agriculture Faisalabad, Faisalabad. 223.
Nilsen, E. T., and Orcutt, D. M. (1996). The Physiology of Plants Under Stress: Abiotic Factors. New York, NY: John Wiley and Sons, Inc., 689.
Parry, M. L., Canziani, O. F., Palutikof, J. P., van der Linden, P. J., and Hanson, C. E. (2007). “Technical Summary. Climate Change 2007: impacts, adaptation and vulnerability,” in Proceedings of the Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, eds M. L. Parry, O. F. Canziani, J. P. Palutikof, P. J. van der Linden, and C. E. Hanson (Cambridge: Cambridge University Press), 23–78.
Peltonen-Sainio, P., Kangas, A., Salo, Y., and Jauhiainen, L. (2007). Grain number dominates grain weight in temperate cereal yield determination: evidence based on 30 years of multi-location trials. Field Crops Res. 100, 179–188. doi: 10.1016/j.fcr.2006.07.002
Peng, S. B., Huang, J. L., and Sheehy, J. E. (2004). Rice yields decline with higher night temperature from global warming. Proc. Natl. Acad. Sci. U.S.A. 101, 9971–9975. doi: 10.1073/pnas.0403720101
Rondanini, D. P., Mantese, A., Savin, R., and Hall, A. J. (2006). Responses of sunflower yield and grain quality to alternating day/night high temperature regimes during grain filling: effects of timing, duration and intensity of exposure to stress. Field Crops Res. 96, 48–62. doi: 10.1016/j.fcr.2005.05.006
Ruales, J. B., and Nair, M. (1992). Nutritional quality of the protein in quinoa (Chenopodium quinoa Willd.) seeds. Plant Foods Hum. Nutr. 42, 1–11. doi: 10.1007/BF02196067
Sadras, V. O. (2007). Evolutionary aspects of the trade-off between grain size and number in crops. Field Crops Res. 100, 125–138. doi: 10.1016/j.fcr.2006.07.004
Sadras, V. O., and Slafer, G. A. (2012). Environmental modulation of yield components in cereals: heritabilities reveal a hierarchy of phenotypic plasticities. Field Crops Res. 127, 215–224. doi: 10.1016/j.fcr.2011.11.014
Savin, R., Stone, P. J., and Nicolas, M. E. (1996). Responses of grain growth and malting quality of barley to short periods of high temperature in field studies using portable chambers. Aust. J. Agric. Res. 47, 465–477. doi: 10.1071/AR9960465
SFA (2014). Sustainable Farming Association. Available at: www.sfa-mn.org/
Shearman, S. (2014). Quinoa, Chia Seeds and Kale: Superfoods or Supermarketing? The Guardian. Available at: www.theguardian.com
Spehar, C. R., and Rocha, J. E. S. (2009). Effect of sowing density on plant growth and development of quinoa, genotype 4.5, in the Brazilian savannah highlands. Biosci. J25, 53–58.
StatPoint Inc. (2005). Statgraphics Centurion XV Software, Version 15.2.06. Herndon, VA: StatPoint, Inc.
Ugarte, C., Calderini, D. F., and Slafer, G. A. (2007). Grain weight and grain number responsiveness to pre-anthesis temperature in wheat, barley and triticale. Field Crop Res. 100, 240–248. doi: 10.1016/j.fcr.2006.07.010
Valle, S., Carrasco, J., Pinochet, D., and Calderini, D. F. (2009). Grain yield above-ground and root biomass of Al-tolerant and Al-sensitive wheat cultivars under different soil aluminium concentrations at field conditions. Plant Soil 318, 299–310. doi: 10.1007/s11104-008-9841-8
Van der Merwe, R. (2010). Genotype by Environment Interaction for Oil Quality in High Oleic Acid Sunflower Lines. Doctoral thesis, Department of Plant Sciences (Plant Breeding) University of the Free State, Bloemfontein, 260.
Van Steveninck, R. F. M. (1959). Factors affecting abscission of reproductive organs in yellow lupins (Lupinus luteus L.). III. Endogenous growth substances in virus-infected and healthy plants and their effect on abscission. J. Exp. Bot. 10, 367–376. doi: 10.1093/jxb/10.3.367
von Baer, I., Bazile, D., and Martinez, E. (2009). Cuarenta años de mejoramiento de quínoa (Chenopodium quinoa Willd.) en la Araucanía: origen de “La Regalona-B”. Revista Geográfica Valparaíso 42, 34–44.
Wang, J., Gan, Y. T., Clarke, F., and McDonald, C. L. (2006). Response of chickpea yield to high temperature stress during reproductive development. Crop Sci. 46, 2171–2178. doi: 10.2135/cropsci2006.02.0092
Westgate, M. E., and Peterson, C. M. (1993). Flower and pod development in water-deficient soybeans (Glycine max L. Merr.). J. Exp. Bot. 44, 109–117. doi: 10.1093/jxb/44.1.109
Keywords: Chenopodium quinoa Willd, grain yield, grain number, temperature increase, climate change, protein concentration
Citation: Lesjak J and Calderini DF (2017) Increased Night Temperature Negatively Affects Grain Yield, Biomass and Grain Number in Chilean Quinoa. Front. Plant Sci. 8:352. doi: 10.3389/fpls.2017.00352
Received: 04 November 2016; Accepted: 28 February 2017;
Published: 23 March 2017.
Edited by:
Marcello Mastrorilli, Consiglio per la Ricerca in Agricoltura e l’Analisi dell’Economia Agraria (CREA), ItalyReviewed by:
Sergio Saia, Consiglio per la Ricerca in Agricoltura e l’Analisi dell’Economia Agraria (CREA), ItalyRobert John French, Department of Agriculture and Food, Australia
Copyright © 2017 Lesjak and Calderini. 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: Jurka Lesjak, jurka.lesjak@gmail.com