WALTer: a three-dimensional wheat model to study competition for light through the prediction of tillering dynamics

Efficient Bayesian automatic calibration of a functional-structural wheat model using an adaptive design and a metamodelling approach

Functional-structural plant models are increasingly being used by plant scientists to address a wide variety of questions. However, the calibration of these complex models is often challenging, mainly because of their high computational cost, and, as a result, error propagation is usually ignored. Here we applied an automatic method to the calibration of WALTer: a functional-structural wheat model that simulates the plasticity of tillering in response to competition for light. We used a Bayesian calibration method to jointly estimate the values of five parameters and quantify their uncertainty by fitting the model outputs to tillering dynamics data. We made recourse to Gaussian process metamodels in order to alleviate the computational cost of WALTer. These metamodels are built from an adaptive design that consists of successive runs of WALTer chosen by an efficient global optimization algorithm specifically adapted to this particular calibration task. The method presented here performed well on both synthetic and experimental data. It is an efficient approach for the calibration of WALTer and should be of interest for the calibration of other functional-structural plant models.

Optimizing plant spatial competition can change phytohormone content and promote tillering, thereby improving wheat yield

As an important type of interplant competition, line-spacing shrinkage and row-spacing expansion (LSRE) can increase the number of tillers and improve resource utilization efficiency in wheat. Wheat tillering is closely related to various phytohormones. However, it is unclear whether LSRE regulates phytohormones and their relationship to tillering and wheat yield. This study evaluated tillering characteristics, phytohormone content in tiller nodes at the pre-winter stage, and grain yield factors for the winter wheat variety Malan1. We used a two-factor randomized block trial design with two sowing spacings of 15 cm (15RS, conventional treatment) and 7.5 cm (7.5RS, LSRE treatment) at the same density and three sowing-date groups (SD1, SD2, and SD3). LSRE significantly promoted wheat tillering and biomass at the pre-winter stage (average increases of 14.5% and 20.9% in the three sowing-date groups, respectively) and shortened the accumulated temperature required for a single tiller. Changes in the levels of phytohormones, including decreased gibberellin and indole acetic acid and increased zeatin riboside and strigolactones, were determined by high-performance liquid chromatography and were shown to be responsible for the tillering process under LSRE treatment in winter wheat. LSRE treatment can improve crop yield by increasing the number of spikes per unit area and grain weight. Our results clarified the changes in tillering and phytohormones content of winter wheat under LSRE treatment and their correlation with grain yield. This study also provides insights into the physiological mechanisms of alleviating inter-plant competition to improve crop yield.

Interactions Between Environment and Genetic Diversity in Perennial Grass Phenology: A Review of Processes at Plant Scale and Modeling

In perennial grasses, the reproductive development consists of major phenological stages which highly determine the seasonal variations of grassland biomass production in terms of quantity and quality. The reproductive development is regulated by climatic conditions through complex interactions subjected to high genetic diversity. Understanding these interactions and their impact on plant development and growth is essential to optimize grassland management and identify the potential consequences of climate change. Here, we review the main stages of reproductive development, from floral induction to heading, i.e., spike emergence, considering the effect of the environmental conditions and the genetic diversity observed in perennial grasses. We first describe the determinants and consequences of reproductive development at individual tiller scale before examining the interactions between plant tillers and their impact on grassland perenniality. Then, we review the available grassland models through their ability to account for the complexity of reproductive development and genetic × environmental interactions. This review shows that (1) The reproductive development of perennial grasses is characterized by a large intraspecific diversity which has the same order of magnitude as the diversity observed between species or environmental conditions. (2) The reproductive development is determined by complex interactions between the processes of floral induction and morphogenesis of the tiller. (3) The perenniality of a plant is dependent on the reproductive behavior of each tiller. (4) Published models only partly explain the complex interactions between morphogenesis and climate on reproductive development. (5) Introducing more explicitly the underlying processes involved in reproductive development in models would improve our ability to anticipate grassland behavior in future growth conditions.

Functional-Structural Plant Modeling Highlights How Diversity in Leaf Dimensions and Tillering Capability Could Promote the Efficiency of Wheat Cultivar Mixtures

Increasing the cultivated diversity has been identified as a major leverage for the agroecological transition as it can help improve the resilience of low input cropping systems. For wheat, which is the most cultivated crop worldwide in terms of harvested area, the use of cultivar mixtures is spreading in several countries, but studies have seldom focused on establishing mixing rules based on plant architecture. Yet, the aerial architecture of plants and the overall canopy structure are critical for field performance as they greatly influence light interception, plant interactions and yield. The very high number of trait combinations in wheat mixtures makes it difficult to conduct experimentations on this issue, which is why a modeling approach appears to be an appropriate solution. In this study, we used WALTer, a functional structural plant model (FSPM), to simulate wheat cultivar mixtures and try to better understand how differences between cultivars in key traits of the aerial architecture influence mixture performance. We simulated balanced binary mixtures of cultivars differing for different critical plant traits: final height, leaf dimensions, leaf insertion angle and tillering capability. Our study highlights the impact of the leaf dimensions and the tillering capability on the performance of the simulated mixtures, which suggests that traits impacting the plants' leaf area index (LAI) have more influence on the performance of the stand than traits impacting the arrangement of the leaves. Our results show that the performance of mixtures is very variable depending on the values of the explored architectural traits. In particular, the best performances were achieved by mixing cultivars with different leaf dimensions and different tillering capability, which is in agreement with numerous studies linking the diversity of functional traits in plant communities to their productivity. However, some of the worst performances were also achieved by mixing varieties differing in their aerial architecture, which suggests that diversity is not a sufficient criterion to design efficient mixtures. Overall, these results highlight the importance of simulation-based explorations for establishing assembly rules to design efficient mixtures.

Regulation of 2,4-D Isooctyl Ester on Triticum aestivum and Aegilops tauschii Tillering and Endogenous Phytohormonal Responses

Tillering is an important agronomic trait essential for the yield of Triticum aestivum and the propagation of Aegilops tauschii. However, the effect of phytohormones on T. aestivum and Ae. tauschii tillering and the underlying regulatory mechanisms remain poorly understood. In the study, we found that T. aestivum and Ae. tauschii exhibited different tillering sensitivities to the auxin herbicide 2,4-D isooctyl ester. At 3 days post-application, tiller bud growth was inhibited by 77.50% in T. aestivum, corresponding to 2.0-fold greater inhibition than that in Ae. tauschii (38.71%). Transcriptome analysis showed that differentially expressed genes (DEGs) in the T. aestivum response to 2,4-D isooctyl ester were mainly enriched in plant hormone metabolism and signal transduction pathways, but similar changes were not observed in Ae. tauschii. Among that, the auxin biosynthesis and signaling induced by 2,4-D isooctyl ester was quite different between the two species. A total of nine candidate genes involved in varied tillering responses were selected from the DEGs and validated by quantitative real-time PCR. Endogenous hormone levels were assayed to further verify the RNA-seq results. After 2,4-D isooctyl ester treatment, a significant increase in abscisic acid (ABA) levels was observed in T. aestivum, whereas ABA levels were relatively stable in Ae. tauschii. The herbicide induced more cytokinin (CTK) accumulation in Ae. tauschii than in T. aestivum. External ABA clearly restricted tiller bud growth in both T. aestivum and Ae. tauschii, while 6-benzyl aminopurine had no significant effect. These results indicate that ABA and CTK may be related with 2,4-D isooctyl ester-regulated tillering differences between the two species, which will help to further understand the mechanism of the auxin-mediated regulation of tillering.

Effects of plant density on tillering in the weed grass Aegilops tauschii Coss. and its phytohormonal regulation

Aegilops tauschii Coss, a notorious wheat field weed, poses a serious threat to wheat in China. Tillers are an important agronomic tool for yield. In this study, a total of 12 Ae. tauschii populations were collected from China to investigate the effect of plant density on tiller occurrence and its phytohormonal regulation. We assayed the growth parameters of Ae. tauschii and the levels of endogenous hormones at different plant densities. The results showed that plant density had a significant effect on the quantity and quality of Ae. tauschii seeds produced per plant. In particular, the tiller and spike numbers per plant were negatively affected by plant density (P < 0.0001). The contents of 13 endogenous hormones in the tiller nodes changed in response to plant density. Among them, indole-3-acetic acid (IAA) and gibberellin (GA) positively responded to plant density. However, the reverse result was found for cytokinin (CTK). Interestingly, phylogenetic tree analysis of auxin (AeYUCCA), CK (AeIPT) and GA (AeCPS) biosynthesis related genes found that phylogenies in the Gramineae for the three different genes were various, hinting at horizontal gene transfer. Moreover, the dynamics of the expression of AeYUCCA, AeIPT and AeCPS were roughly consistent with their phytohormone contents during tillering stage. When externally sprayed on plants of Ae. tauschii, 2,4-D isooctyl ester and GA₃ markedly reduced its tillering while 6-BA had no significant effect.

Regulation of Aegilops tauschii Coss Tiller Bud Growth by Plant Density: Transcriptomic, Physiological and Phytohormonal Responses

Aegilops tauschii Coss is one of the most hazardous weeds that severely infests wheat fields in China. The tillering ability of Ae. tauschii strongly affects the occurrence and spread by influencing its seed output. In this study, Ae. tauschii was sown at low plant density (LPD) and high plant density (HPD) to investigate the effect of plant density on tiller bud outgrowth and its potential regulators using RNA-Seq. Additionally, the chlorophyll content and photosynthesis, soluble sugar and phytohormone levels were also determined at different plant densities. The results showed that an increased plant density significantly inhibited the elongation of tiller buds in the axil of the first leaf at 15 days after planting, with 7.69 mm at LPD and 1.69 mm at HPD. A total of seven putative tiller-related genes were selected and validated using quantitative real-time PCR. Furthermore, chlorophyll levels, photosynthetic efficiency, and soluble sugar contents were distinctly inhibited by HPD in Ae. tauschii, which may be responsible for the restriction of tiller bud growth. In addition, differentially expressed genes (DEGs) were markedly enriched in indole-3-acetic acid (IAA), abscisic acid (ABA), and gibberellin metabolism and signaling. Accordingly, the levels of ABA and gibberellin A3 in Ae. tauschii were strikingly higher at HPD compared with those at LPD, yet the reverse tendency was observed for IAA. Undoubtedly, such results will be highly beneficial for illuminating the underlying regulators of the Ae. tauschii tillering response to plant density and may provide new ideas for the control of this weed in the future.