Growing mixtures of annual arable crop species or genotypes is a promising way to improve crop production without increasing agricultural inputs. To design optimal crop mixtures, choices of species, genotypes, sowing proportion, plant arrangement, and sowing date need to be made but field experiments alone are not sufficient to explore such a large range of factors. Crop modeling allows to study, understand, and ultimately design cropping systems and is an established method for sole crops. Recently, modeling started to be applied to annual crop mixtures as well. Here, we review to what extent crop simulation models and individual-based models are suitable to capture and predict the specificities of annual crop mixtures. We argued that (1) the crop mixture spatio-temporal heterogeneity (influencing the occurrence of ecological processes) determines the choice of the modeling approach (plant or crop centered). (2) Only few crop models (adapted from sole crop models) and individual-based models currently exist to simulate annual crop mixtures. Crop models are mainly used to address issues related to both crop mixtures management and the integration of crop mixtures into larger scales such as the rotation. In contrast, individual-based models are mainly used to identify plant traits involved in crop mixture performance and to quantify the relative contribution of the different ecological processes (niche complementarity, facilitation, competition, plasticity) to crop mixture functioning. This review highlights that modeling of annual crop mixtures is in its infancy and gives to model users some important keys to choose the model based on the questions they want to answer, with awareness of the strengths and weaknesses of each of the modeling approaches.

Background and Aims Predicting light partitioning in crop mixtures is a critical step in improving the productivity of such complex systems, and light interception has been shown to be closely linked to plant architecture. The aim of the present work was to analyse the relationships between plant architecture and light partitioning within wheat–pea (Triticum aestivum–Pisum sativum) mixtures. An existing model for wheat was utilized and a new model for pea morphogenesis was developed. Both models were then used to assess the effects of architectural variations in light partitioning. Methods First, a deterministic model (L-Pea) was developed in order to obtain dynamic reconstructions of pea architecture. The L-Pea model is based on L-systems formalism and consists of modules for ‘vegetative development’ and ‘organ extension’. A tripartite simulator was then built up from pea and wheat models interfaced with a radiative transfer model. Architectural parameters from both plant models, selected on the basis of their contribution to leaf area index (LAI), height and leaf geometry, were then modified in order to generate contrasting architectures of wheat and pea. Key results By scaling down the analysis to the organ level, it could be shown that the number of branches/tillers and length of internodes significantly determined the partitioning of light within mixtures. Temporal relationships between light partitioning and the LAI and height of the different species showed that light capture was mainly related to the architectural traits involved in plant LAI during the early stages of development, and in plant height during the onset of interspecific competition. Conclusions In silico experiments enabled the study of the intrinsic effects of architectural parameters on the partitioning of light in crop mixtures of wheat and pea. The findings show that plant architecture is an important criterion for the identification/breeding of plant ideotypes, particularly with respect to light partitioning.

Backgrounds and Aims: A major challenge when supporting the development of intercropping systems remains the design of efficient species mixtures. The ecological processes that sustain overyielding of legume-based mixtures compared to pure crops are well known, but their links to plant traits remain to be unravelled. A common assumption is that enhancing trait divergence among species for resource acquisition when assembling plant mixtures should increase species complementarity and improve community performance. Methods: The Virtual Grassland model was used to assess how divergence in trait values between species on four physiological functions (namely light and mineral N acquisition, temporal development, and C–N use efficiency) affected overyielding and mixture stability in legume-based binary mixtures. A first step allowed us to identify the model parameters that were most important to interspecies competition. A second step involved testing the impact of convergent and divergent parameter (or trait) values between species on virtual mixture performance. Results: Maximal overyielding was achieved in cases where trait values were divergent for the physiological functions controlling N acquisition and temporal development but convergent for light interception. It was also found that trait divergence should not affect competitive abilities of legume and non-legumes at random. Indeed, random trait combinations frequently led to reduced mixture yields when compared to a perfectly convergent neutral model. Combinations with the highest overyielding also tended to be associated with mixture instability and decreasing legume biomass proportion. Achieving both high overyielding and mixture stability was only found to be possible under low or moderate N levels, using combinations of traits adapted to each environment. Conclusions: No simple assembly rule based on trait divergence could be confirmed. Plant models able to infer plant–plant interactions can be helpful for the identification of major interaction traits and the definition of ideotypes adapted to a targeted intercropping system.

Background and Aims Improving crops requires better linking of traits and metabolic processes to whole plant performance. In this paper, we present CN-Wheat, a comprehensive and mechanistic model of carbon (C) and nitrogen (N) metabolism within wheat culms after anthesis. Methods The culm is described by modules that represent the roots, photosynthetic organs and grains. Each of them includes structural, storage and mobile materials. Fluxes of C and N among modules occur through a common pool and through transpiration flow. Metabolite variations are represented by differential equations that depend on the physiological processes occurring in each module. A challenging aspect of CN-Wheat lies in the regulation of these processes by metabolite concentrations and the environment perceived by organs. Key Results CN-Wheat simulates the distribution of C and N into wheat culms in relation to photosynthesis, N uptake, metabolite turnover, root exudation and tissue death. Regulation of physiological activities by local concentrations of metabolites appears to be a valuable feature for understanding how the behaviour of the whole plant can emerge from local rules. Conclusions The originality of CN-Wheat is that it proposes an integrated view of plant functioning based on a mechanistic approach. The formalization of each process can be further refined in the future as knowledge progresses. This approach is expected to strengthen our capacity to understand plant responses to their environment and investigate plant traits adapted to changes in agronomical practices or environmental conditions. A companion paper will evaluate the model.

Most studies dealing with light partitioning in intercropping systems have used statistical models based on the turbid medium approach, thus assuming homogeneous canopies. However, these models could not be directly validated although spatial heterogeneities could arise in such canopies. The aim of the present study was to assess the ability of the turbid medium approach to accurately estimate light partitioning within grass–legume mixed canopies. Three contrasted mixtures of wheat–pea, tall fescue–alfalfa and tall fescue–clover were sown according to various patterns and densities. Three-dimensional plant mock-ups were derived from magnetic digitizations carried out at different stages of development. The benchmarks for light interception efficiency (LIE) estimates were provided by the combination of a light projective model and plant mock-ups, which also provided the inputs of a turbid medium model (SIRASCA), i.e. leaf area index and inclination. SIRASCA was set to gradually account for vertical heterogeneity of the foliage, i.e. the canopy was described as one, two or ten horizontal layers of leaves. Mixtures exhibited various and heterogeneous profiles of foliar distribution, leaf inclination and component species height. Nevertheless, most of the LIE was satisfactorily predicted by SIRASCA. Biased estimations were, however, observed for (1) grass species and (2) tall fescue–alfalfa mixtures grown at high density. Most of the discrepancies were due to vertical heterogeneities and were corrected by increasing the vertical description of canopies although, in practice, this would require time-consuming measurements. The turbid medium analogy could be successfully used in a wide range of canopies. However, a more detailed description of the canopy is required for mixtures exhibiting vertical stratifications and inter-/intra-species foliage overlapping. Architectural models remain a relevant tool for studying light partitioning in intercropping systems that exhibit strong vertical heterogeneities. Moreover, these models offer the possibility to integrate the effects of microclimate variations on plant growth.

There is currently a strong societal demand for sustainability, quality, and safety in bread wheat production. To address these challenges, new and innovative knowledge, resources, tools, and methods to facilitate breeding are needed. This starts with the development of high throughput genomic tools including single nucleotide polymorphism (SNP) arrays, high density molecular marker maps, and full genome sequences. Such powerful tools are essential to perform genome-wide association studies (GWAS), to implement genomic and phenomic selection, and to characterize the worldwide diversity. This is also useful to breeders to broaden the genetic basis of elite varieties through the introduction of novel sources of genetic diversity. Improvement in varieties particularly relies on the detection of genomic regions involved in agronomical traits including tolerance to biotic (diseases and pests) and abiotic (drought, nutrient deficiency, high temperature) stresses. When enough resolution is achieved, this can result in the identification of candidate genes that could further be characterized to identify relevant alleles. Breeding must also now be approached through in silico modeling to simulate plant development, investigate genotype × environment interactions, and introduce marker–trait linkage information in the models to better implement genomic selection. Breeders must be aware of new developments and the information must be made available to the world wheat community to develop new high-yielding varieties that can meet the challenge of higher wheat production in a sustainable and fluctuating agricultural context. In this review, we compiled all knowledge and tools produced during the BREEDWHEAT project to show how they may contribute to face this challenge in the coming years.