Crops In Silico: Generating Virtual Crops Using an Integrative and Multi-scale Modeling Platform

Regulation of root growth and elongation in wheat

Currently, the control of rhizosphere selection on farms has been applied to achieve enhancements in phenotype, extending from improvements in single root characteristics to the dynamic nature of entire crop systems. Several specific signals, regulatory elements, and mechanisms that regulate the initiation, morphogenesis, and growth of new lateral or adventitious root species have been identified, but much more work remains. Today, phenotyping technology drives the development of root traits. Available models for simulation can support all phenotyping decisions (root trait improvement). The detection and use of markers for quantitative trait loci (QTLs) are effective for enhancing selection efficiency and increasing reproductive genetic gains. Furthermore, QTLs may help wheat breeders select the appropriate roots for efficient nutrient acquisition. Single-nucleotide polymorphisms (SNPs) or alignment of sequences can only be helpful when they are associated with phenotypic variation for root development and elongation. Here, we focus on major root development processes and detail important new insights recently generated regarding the wheat genome. The first part of this review paper discusses the root morphology, apical meristem, transcriptional control, auxin distribution, phenotyping of the root system, and simulation models. In the second part, the molecular genetics of the wheat root system, SNPs, TFs, and QTLs related to root development as well as genome editing (GE) techniques for the improvement of root traits in wheat are discussed. Finally, we address the effect of omics strategies on root biomass production and summarize existing knowledge of the main molecular mechanisms involved in wheat root development and elongation.

Predicting the unpredictable: the regulatory nature and promiscuity of herbicide cross resistance

The emergence of herbicide-resistant weeds is a significant threat to modern agriculture. Cross resistance, a phenomenon where resistance to one herbicide confers resistance to another, is a particular concern owing to its unpredictability. Nontarget-site (NTS) cross resistance is especially challenging to predict, as it arises from genes that encode enzymes that do not directly involve the herbicide target site and can affect multiple herbicides. Recent advancements in genomic and structural biology techniques could provide new venues for predicting NTS resistance in weed species. In this review, we present an overview of the latest approaches that could be used. We discuss the use of genomic and epigenomics techniques such as ATAC-seq and DAP-seq to identify transcription factors and cis-regulatory elements associated with resistance traits. Enzyme/protein structure prediction and docking analysis are discussed as an initial step for predicting herbicide binding affinities with key enzymes to identify candidates for subsequent in vitro validation. We also provide example analyses that can be deployed toward elucidating cross resistance and its regulatory patterns. Ultimately, our review provides important insights into the latest scientific advancements and potential directions for predicting and managing herbicide cross resistance in weeds. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Model-assisted analysis for tuning anthocyanin composition in grape berries

Anthocyanin composition is responsible for the red colour of grape berries and wines, and contributes to their organoleptic quality. However, anthocyanin biosynthesis is under genetic, developmental and environmental regulation, making its targeted fine-tuning challenging. We constructed a mechanistic model to simulate the dynamics of anthocyanin composition throughout grape ripening in Vitis vinifera, employing a consensus anthocyanin biosynthesis pathway. The model was calibrated and validated using six datasets from eight cultivars and 37 growth conditions. Tuning the transformation and degradation parameters allowed us to accurately simulate the accumulation process of each individual anthocyanin under different environmental conditions. The model parameters were robust across environments for each genotype. The coefficients of determination (R2) for the simulated versus observed values for the six datasets ranged from 0.92 to 0.99, while the relative root mean square errors (RRMSEs) were between 16.8 and 42.1 %. The leave-one-out cross-validation for three datasets showed R2 values of 0.99, 0.96 and 0.91, and RRMSE values of 28.8, 32.9 and 26.4 %, respectively, suggesting a high prediction quality of the model. Model analysis showed that the anthocyanin profiles of diverse genotypes are relatively stable in response to parameter perturbations. Virtual experiments further suggested that targeted anthocyanin profiles may be reached by manipulating a minimum of three parameters, in a genotype-dependent manner. This model presents a promising methodology for characterizing the temporal progression of anthocyanin composition, while also offering a logical foundation for bioengineering endeavours focused on precisely adjusting the anthocyanin composition of grapes.

How to engineer the unknown: Advancing a quantitative and predictive understanding of plant and soil biology to address climate change

Our basic understanding of carbon cycling in the biosphere remains qualitative and incomplete, precluding our ability to effectively engineer novel solutions to climate change. How can we attempt to engineer the unknown? This challenge has been faced before in plant biology, providing a roadmap to guide future efforts. We use examples from over a century of photosynthesis research to illustrate the key principles that will set future plant engineering on a solid footing, namely, an effort to identify the key control variables, quantify the effects of systematically tuning these variables, and use theory to account for these observations. The main contributions of plant synthetic biology will stem not from delivering desired genotypes but from enabling the kind of predictive understanding necessary to rationally design these genotypes in the first place. Only then will synthetic plant biology be able to live up to its promise.

Cell-type-specific metabolism in plants

Every plant organ contains tens of different cell types, each with a specialized function. These functions are intrinsically associated with specific metabolic flux distributions that permit the synthesis of the ATP, reducing equivalents and biosynthetic precursors demanded by the cell. Investigating such cell-type-specific metabolism is complicated by the mosaic of different cells within each tissue combined with the relative scarcity of certain types. However, techniques for the isolation of specific cells, their analysis in situ by microscopy, or modeling of their function in silico have permitted insight into cell-type-specific metabolism. In this review we present some of the methods used in the analysis of cell-type-specific metabolism before describing what we know about metabolism in several cell types that have been studied in depth; (i) leaf source and sink cells; (ii) glandular trichomes that are capable of rapid synthesis of specialized metabolites; (iii) guard cells that must accumulate large quantities of the osmolytes needed for stomatal opening; (iv) cells of seeds involved in storage of reserves; and (v) the mesophyll and bundle sheath cells of C4 plants that participate in a CO₂ concentrating cycle. Metabolism is discussed in terms of its principal features, connection to cell function and what factors affect the flux distribution. Demand for precursors and energy, availability of substrates and suppression of deleterious processes are identified as key factors in shaping cell-type-specific metabolism.

Modelling plants across scales of biological organisation for guiding crop improvement

Grain yield improvement in globally important staple crops is critical in the coming decades if production is to keep pace with growing demand; so there is increasing interest in understanding and manipulating plant growth and developmental traits for better crop productivity. However, this is confounded by complex cross-scale feedback regulations and a limited ability to evaluate the consequences of manipulation on crop production. Plant/crop modelling could hold the key to deepening our understanding of dynamic trait-crop-environment interactions and predictive capabilities for supporting genetic manipulation. Using photosynthesis and crop growth as an example, this review summarises past and present experimental and modelling work, bringing about a model-guided crop improvement thrust, encompassing research into: (1) advancing cross-scale plant/crop modelling that connects across biological scales of organisation using a trait dissection-integration modelling principle; (2) improving the reliability of predicted molecular-trait-crop-environment system dynamics with experimental validation; and (3) innovative model application in synergy with cross-scale experimentation to evaluate G×M×E and predict yield outcomes of genetic intervention (or lack of it) for strategising further molecular and breeding efforts. The possible future roles of cross-scale plant/crop modelling in maximising crop improvement are discussed.

Lettuce Production in Intelligent Greenhouses-3D Imaging and Computer Vision for Plant Spacing Decisions

Recent studies indicate that food demand will increase by 35-56% over the period 2010-2050 due to population increase, economic development, and urbanization. Greenhouse systems allow for the sustainable intensification of food production with demonstrated high crop production per cultivation area. Breakthroughs in resource-efficient fresh food production merging horticultural and AI expertise take place with the international competition "Autonomous Greenhouse Challenge". This paper describes and analyzes the results of the third edition of this competition. The competition's goal is the realization of the highest net profit in fully autonomous lettuce production. Two cultivation cycles were conducted in six high-tech greenhouse compartments with operational greenhouse decision-making realized at a distance and individually by algorithms of international participating teams. Algorithms were developed based on time series sensor data of the greenhouse climate and crop images. High crop yield and quality, short growing cycles, and low use of resources such as energy for heating, electricity for artificial light, and CO2 were decisive in realizing the competition's goal. The results highlight the importance of plant spacing and the moment of harvest decisions in promoting high crop growth rates while optimizing greenhouse occupation and resource use. In this paper, images taken with depth cameras (RealSense) for each greenhouse were used by computer vision algorithms (Deepabv3+ implemented in detectron2 v0.6) in deciding optimum plant spacing and the moment of harvest. The resulting plant height and coverage could be accurately estimated with an R2 of 0.976, and a mIoU of 98.2, respectively. These two traits were used to develop a light loss and harvest indicator to support remote decision-making. The light loss indicator could be used as a decision tool for timely spacing. Several traits were combined for the harvest indicator, ultimately resulting in a fresh weight estimation with a mean absolute error of 22 g. The proposed non-invasively estimated indicators presented in this article are promising traits to be used towards full autonomation of a dynamic commercial lettuce growing environment. Computer vision algorithms act as a catalyst in remote and non-invasive sensing of crop parameters, decisive for automated, objective, standardized, and data-driven decision making. However, spectral indexes describing lettuces growth and larger datasets than the currently accessible are crucial to address existing shortcomings between academic and industrial production systems that have been encountered in this work.

Breeding crops for drought-affected environments and improved climate resilience

Breeding climate-resilient crops with improved levels of abiotic and biotic stress resistance as a response to climate change presents both opportunities and challenges. Applying the framework of the "breeder's equation," which is used to predict the response to selection for a breeding program cycle, we review methodologies and strategies that have been used to successfully breed crops with improved levels of drought resistance, where the target population of environments (TPEs) is a spatially and temporally heterogeneous mixture of drought-affected and favorable (water-sufficient) environments. Long-term improvement of temperate maize for the US corn belt is used as a case study and compared with progress for other crops and geographies. Integration of trait information across scales, from genomes to ecosystems, is needed to accurately predict yield outcomes for genotypes within the current and future TPEs. This will require transdisciplinary teams to explore, identify, and exploit novel opportunities to accelerate breeding program outcomes; both improved germplasm resources and improved products (cultivars, hybrids, clones, and populations) that outperform and replace the products in use by farmers, in combination with modified agronomic management strategies suited to their local environments.

A cross-scale analysis to understand and quantify the effects of photosynthetic enhancement on crop growth and yield across environments

Photosynthetic manipulation provides new opportunities for enhancing crop yield. However, understanding and quantifying the importance of individual and multiple manipulations on the seasonal biomass growth and yield performance of target crops across variable production environments is limited. Using a state-of-the-art cross-scale model in the APSIM platform we predicted the impact of altering photosynthesis on the enzyme-limited (A_c ) and electron transport-limited (A_j ) rates, seasonal dynamics in canopy photosynthesis, biomass growth, and yield formation via large multiyear-by-location crop growth simulations. A broad list of promising strategies to improve photosynthesis for C₃ wheat and C₄ sorghum were simulated. In the top decile of seasonal outcomes, yield gains were predicted to be modest, ranging between 0% and 8%, depending on the manipulation and crop type. We report how photosynthetic enhancement can affect the timing and severity of water and nitrogen stress on the growing crop, resulting in nonintuitive seasonal crop dynamics and yield outcomes. We predicted that strategies enhancing A_c alone generate more consistent but smaller yield gains across all water and nitrogen environments, A_j enhancement alone generates larger gains but is undesirable in more marginal environments. Large increases in both A_c and A_j generate the highest gains across all environments. Yield outcomes of the tested manipulation strategies were predicted and compared for realistic Australian wheat and sorghum production. This study uniquely unpacks complex cross-scale interactions between photosynthesis and seasonal crop dynamics and improves understanding and quantification of the potential impact of photosynthesis traits (or lack of it) for crop improvement research.

Enhancing control systems of higher plant culture chambers via multilevel structural mechanistic modelling

Modelling higher plant growth is of strategic interest for modern agriculture as well as for the development of bioregenerative life support systems for space applications, where crop growth is expected to play an essential role. The capability of constraint-based metabolic models to cope the diel dynamics of plants growth is integrated into a multilevel modelling approach including mass and energy transfer and enzyme kinetics. Lactuca sativa is used as an exemplary crop to validate, with experimental data, the approach presented as well as to design a novel model-based predictive control strategy embedding metabolic information. The proposed modelling strategy predicts with high accuracy the dynamics of gas exchange and the distribution of fluxes in the metabolic network whereas the control architecture presented can be useful to manage higher plants chambers and open new ways of merging metabolome and control algorithms.

Plasticity in stomatal behaviour across a gradient of water supply is consistent among field-grown maize inbred lines with varying stomatal patterning

Stomata regulate leaf CO₂ assimilation (A) and water loss. The Ball-Berry and Medlyn models predict stomatal conductance (g_s ) with a slope parameter (m or g₁ ) that reflects the sensitivity of g_s to A, atmospheric CO₂  and humidity, and is inversely related to water use efficiency (WUE). This study addressed knowledge gaps about what the values of m and g₁ are in C₄ crops under field conditions, as well as how they vary among genotypes and with drought stress. Four inbred maize genotypes were unexpectedly consistent in how m and g₁ decreased as water supply decreased. This was despite genotypic variation in stomatal patterning, A and g_s . m and g₁ were strongly correlated with soil water content, moderately correlated with predawn leaf water potential (Ψ_pd ), but not correlated with midday leaf water potential (Ψ_md ). This implied that m and g₁ respond to long-term water supply more than short-term drought stress. The conserved nature of m and g₁ across anatomically diverse genotypes and water supplies suggests there is flexibility in structure-function relationships underpinning WUE. This evidence can guide the simulation of maize g_s across a range of water supply in the primary maize growing region and inform efforts to improve WUE.

3dCAP-Wheat: An Open-Source Comprehensive Computational Framework Precisely Quantifies Wheat Foliar, Nonfoliar, and Canopy Photosynthesis

Canopy photosynthesis is the sum of photosynthesis of all above-ground photosynthetic tissues. Quantitative roles of nonfoliar tissues in canopy photosynthesis remain elusive due to methodology limitations. Here, we develop the first complete canopy photosynthesis model incorporating all above-ground photosynthetic tissues and validate this model on wheat with state-of-the-art gas exchange measurement facilities. The new model precisely predicts wheat canopy gas exchange rates at different growth stages, weather conditions, and canopy architectural perturbations. Using the model, we systematically study (1) the contribution of both foliar and nonfoliar tissues to wheat canopy photosynthesis and (2) the responses of wheat canopy photosynthesis to plant physiological and architectural changes. We found that (1) at tillering, heading, and milking stages, nonfoliar tissues can contribute ~4, ~32, and ~50% of daily gross canopy photosynthesis (A_cgross; ~2, ~15, and ~-13% of daily net canopy photosynthesis, A_cnet) and absorb ~6, ~42, and ~60% of total light, respectively; (2) under favorable condition, increasing spike photosynthetic activity, rather than enlarging spike size or awn size, can enhance canopy photosynthesis; (3) covariation in tissue respiratory rate and photosynthetic rate may be a major factor responsible for less than expected increase in daily A_cnet; and (4) in general, erect leaves, lower spike position, shorter plant height, and proper plant densities can benefit daily A_cnet. Overall, the model, together with the facilities for quantifying plant architecture and tissue gas exchange, provides an integrated platform to study canopy photosynthesis and support rational design of photosynthetically efficient wheat crops.

Improving Soil Resource Uptake by Plants Through Capitalizing on Synergies Between Root Architecture and Anatomy and Root-Associated Microorganisms

Root architectural and anatomical phenotypes are highly diverse. Specific root phenotypes can be associated with better plant growth under low nutrient and water availability. Therefore, root ideotypes have been proposed as breeding targets for more stress-resilient and resource-efficient crops. For example, root phenotypes that correspond to the Topsoil Foraging ideotype are associated with better plant growth under suboptimal phosphorus availability, and root phenotypes that correspond to the Steep, Cheap and Deep ideotype are linked to better performance under suboptimal availability of nitrogen and water. We propose that natural variation in root phenotypes translates into a diversity of different niches for microbial associations in the rhizosphere, rhizoplane and root cortex, and that microbial traits could have synergistic effects with the beneficial effect of specific root phenotypes. Oxygen and water content, carbon rhizodeposition, nutrient availability, and root surface area are all factors that are modified by root anatomy and architecture and determine the structure and function of the associated microbial communities. Recent research results indicate that root characteristics that may modify microbial communities associated with maize include aerenchyma, rooting angle, root hairs, and lateral root branching density. Therefore, the selection of root phenotypes linked to better plant growth under specific edaphic conditions should be accompanied by investigating and selecting microbial partners better adapted to each set of conditions created by the corresponding root phenotype. Microbial traits such as nitrogen transformation, phosphorus solubilization, and water retention could have synergistic effects when correctly matched with promising plant root ideotypes for improved nutrient and water capture. We propose that elucidation of the interactive effects of root phenotypes and microbial functions on plant nutrient and water uptake offers new opportunities to increase crop yields and agroecosystem sustainability.

Integrated root phenotypes for improved rice performance under low nitrogen availability

Greater nitrogen efficiency would substantially reduce the economic, energy and environmental costs of rice production. We hypothesized that synergistic balancing of the costs and benefits for soil exploration among root architectural phenes is beneficial under suboptimal nitrogen availability. An enhanced implementation of the functional-structural model OpenSimRoot for rice integrated with the ORYZA_v3 crop model was used to evaluate the utility of combinations of root architectural phenes, namely nodal root angle, the proportion of smaller diameter nodal roots, nodal root number; and L-type and S-type lateral branching densities, for plant growth under low nitrogen. Multiple integrated root phenotypes were identified with greater shoot biomass under low nitrogen than the reference cultivar IR64. The superiority of these phenotypes was due to synergism among root phenes rather than the expected additive effects of phene states. Representative optimal phenotypes were predicted to have up to 80% greater grain yield with low N supply in the rainfed dry direct-seeded agroecosystem over future weather conditions, compared to IR64. These phenotypes merit consideration as root ideotypes for breeding rice cultivars with improved yield under rainfed dry direct-seeded conditions with limited nitrogen availability. The importance of phene synergism for the performance of integrated phenotypes has implications for crop breeding.

A hybrid kinetic and constraint-based model of leaf metabolism allows predictions of metabolic fluxes in different environments

While flux balance analysis (FBA) provides a framework for predicting steady-state leaf metabolic network fluxes, it does not readily capture the response to environmental variables without being coupled to other modelling formulations. To address this, we coupled an FBA model of 903 reactions of soybean (Glycine max) leaf metabolism with e-photosynthesis, a dynamic model that captures the kinetics of 126 reactions of photosynthesis and associated chloroplast carbon metabolism. Successful coupling was achieved in an iterative formulation in which fluxes from e-photosynthesis were used to constrain the FBA model and then, in turn, fluxes computed from the FBA model used to update parameters in e-photosynthesis. This process was repeated until common fluxes in the two models converged. Coupling did not hamper the ability of the kinetic module to accurately predict the carbon assimilation rate, photosystem II electron flux, and starch accumulation of field-grown soybean at two CO2 concentrations. The coupled model also allowed accurate predictions of additional parameters such as nocturnal respiration, as well as analysis of the effect of light intensity and elevated CO2 on leaf metabolism. Predictions included an unexpected decrease in the rate of export of sucrose from the leaf at high light, due to altered starch-sucrose partitioning, and altered daytime flux modes in the tricarboxylic acid cycle at elevated CO2 . Mitochondrial fluxes were notably different between growing and mature leaves, with greater anaplerotic, tricarboxylic acid cycle and mitochondrial ATP synthase fluxes predicted in the former, primarily to provide carbon skeletons and energy for protein synthesis.

Challenges for a Massive Implementation of Phenomics in Plant Breeding Programs

Due to climate change and expected food shortage in the coming decades, not only will it be necessary to develop cultivars with greater tolerance to environmental stress, but it is also imperative to reduce breeding cycle time. In addition to yield evaluation, plant breeders resort to many sensory assessments and some others of intermediate complexity. However, to develop cultivars better adapted to current/future constraints, it is necessary to incorporate a new set of traits, such as morphophysiological and physicochemical attributes, information relevant to the successful selection of genotypes or parents. Unfortunately, because of the large number of genotypes to be screened, measurements with conventional equipment are unfeasible, especially under field conditions. High-throughput plant phenotyping (HTPP) facilitates collecting a significant amount of data quickly; however, it is necessary to transform all this information (e.g., plant reflectance) into helpful descriptors to the breeder. To the extent that a holistic characterization of the plant (phenomics) is performed in challenging environments, it will be possible to select the best genotypes (forward phenomics) objectively but also understand why the said individual differs from the rest (reverse phenomics). Unfortunately, several elements had prevented phenomics from developing as desired. Consequently, a new set of prediction/validation methodologies, seasonal ambient information, and the fusion of data matrices (e.g., genotypic and phenotypic information) need to be incorporated into the modeling. In this sense, for the massive implementation of phenomics in plant breeding, it will be essential to count an interdisciplinary team that responds to the urgent need to release material with greater capacity to tolerate environmental stress. Therefore, breeding programs should (i) be more efficient (e.g., early discarding of unsuitable material), (ii) have shorter breeding cycles (fewer crosses to achieve the desired cultivar), and (iii) be more productive, increasing the probability of success at the end of the breeding process (percentage of cultivars released to the number of initial crosses).

Optimizing Crop Water Use for Drought and Climate Change Adaptation Requires a Multi-Scale Approach

Selection criteria that co-optimize water use efficiency and yield are needed to promote plant productivity in increasingly challenging and variable drought scenarios, particularly dryland cereals in the semi-arid tropics. Optimizing water use efficiency and yield fundamentally involves transpiration dynamics, where restriction of maximum transpiration rate helps to avoid early crop failure, while maximizing grain filling. Transpiration restriction can be regulated by multiple mechanisms and involves cross-organ coordination. This coordination involves complex feedbacks and feedforwards over time scales ranging from minutes to weeks, and from spatial scales ranging from cell membrane to crop canopy. Aquaporins have direct effect but various compensation and coordination pathways involve phenology, relative root and shoot growth, shoot architecture, root length distribution profile, as well as other architectural and anatomical aspects of plant form and function. We propose gravimetric phenotyping as an integrative, cross-scale solution to understand the dynamic, interwoven, and context-dependent coordination of transpiration regulation. The most fruitful breeding strategy is likely to be that which maintains focus on the phene of interest, namely, daily and season level transpiration dynamics. This direct selection approach is more precise than yield-based selection but sufficiently integrative to capture attenuating and complementary factors.

Root traits as drivers of plant and ecosystem functioning: current understanding, pitfalls and future research needs

The effects of plants on the biosphere, atmosphere and geosphere are key determinants of terrestrial ecosystem functioning. However, despite substantial progress made regarding plant belowground components, we are still only beginning to explore the complex relationships between root traits and functions. Drawing on the literature in plant physiology, ecophysiology, ecology, agronomy and soil science, we reviewed 24 aspects of plant and ecosystem functioning and their relationships with a number of root system traits, including aspects of architecture, physiology, morphology, anatomy, chemistry, biomechanics and biotic interactions. Based on this assessment, we critically evaluated the current strengths and gaps in our knowledge, and identify future research challenges in the field of root ecology. Most importantly, we found that belowground traits with the broadest importance in plant and ecosystem functioning are not those most commonly measured. Also, the estimation of trait relative importance for functioning requires us to consider a more comprehensive range of functionally relevant traits from a diverse range of species, across environments and over time series. We also advocate that establishing causal hierarchical links among root traits will provide a hypothesis-based framework to identify the most parsimonious sets of traits with the strongest links on functions, and to link genotypes to plant and ecosystem functioning.

X-ray computed tomography for 3D plant imaging

X-ray computed tomography (CT) is a valuable tool for 3D imaging of plant tissues and organs. Applications include the study of plant development and organ morphogenesis, as well as modeling of transport processes in plants. Some challenges remain, however, including attaining higher contrast for easier quantification, increasing the resolution for imaging subcellular features, and decreasing image acquisition and processing time for high-throughput phenotyping. In addition, phase contrast, multispectral, dark-field, soft X-ray, and time-resolved imaging are emerging. At the same time, a large amount of 3D image data are becoming available, posing challenges for data management. We review recent advances in the area of X-ray CT for plant imaging, and describe opportunities for using such images for studying transport processes in plants.

Shaping the Organ: A Biologist Guide to Quantitative Models of Plant Morphogenesis

Organ morphogenesis is the process of shape acquisition initiated with a small reservoir of undifferentiated cells. In plants, morphogenesis is a complex endeavor that comprises a large number of interacting elements, including mechanical stimuli, biochemical signaling, and genetic prerequisites. Because of the large body of data being produced by modern laboratories, solving this complexity requires the application of computational techniques and analyses. In the last two decades, computational models combined with wet-lab experiments have advanced our understanding of plant organ morphogenesis. Here, we provide a comprehensive review of the most important achievements in the field of computational plant morphodynamics. We present a brief history from the earliest attempts to describe plant forms using algorithmic pattern generation to the evolution of quantitative cell-based models fueled by increasing computational power. We then provide an overview of the most common types of "digital plant" paradigms, and demonstrate how models benefit from diverse techniques used to describe cell growth mechanics. Finally, we highlight the development of computational frameworks designed to resolve organ shape complexity through integration of mechanical, biochemical, and genetic cues into a quantitative standardized and user-friendly environment.

3D phytomer-based geometric modelling method for plants-the case of maize

Geometric plant modelling is crucial in in silico plants. Existing geometric modelling methods have focused on the topological structure and basic organ profiles, simplifying the morphological features. However, the models cannot effectively differentiate cultivars, limiting FSPM application in crop breeding and management. This study proposes a 3D phytomer-based geometric modelling method with maize (Zea Mays) as the representative plant. Specifically, conversion methods between skeleton and mesh models of 3D phytomer are specified. This study describes the geometric modelling of maize shoots and populations by assembling 3D phytomers. Results show that the method can quickly and efficiently construct 3D models of maize plants and populations, with the ability to show morphological, structural and functional differences among four representative cultivars. The method takes into account both the geometric modelling efficiency and 3D detail features to achieve automatic operation of geometric modelling through the standardized description of 3D phytomers. Therefore, this study provides a theoretical and technical basis for the research and application of in silico plants.

Single cell gene regulatory networks in plants: Opportunities for enhancing climate change stress resilience

Global warming poses major challenges for plant survival and agricultural productivity. Thus, efforts to enhance stress resilience in plants are key strategies for protecting food security. Gene regulatory networks (GRNs) are a critical mechanism conferring stress resilience. Until recently, predicting GRNs of the individual cells that make up plants and other multicellular organisms was impeded by aggregate population scale measurements of transcriptome and other genome-scale features. With the advancement of high-throughput single cell RNA-seq and other single cell assays, learning GRNs for individual cells is now possible, in principle. In this article, we report on recent advances in experimental and analytical methodologies for single cell sequencing assays especially as they have been applied to the study of plants. We highlight recent advances and ongoing challenges for scGRN prediction, and finally, we highlight the opportunity to use scGRN discovery for studying and ultimately enhancing abiotic stress resilience in plants.

Environment-coupled models of leaf metabolism

The plant leaf is the main site of photosynthesis. This process converts light energy and inorganic nutrients into chemical energy and organic building blocks for the biosynthesis and maintenance of cellular components and to support the growth of the rest of the plant. The leaf is also the site of gas–water exchange and due to its large surface, it is particularly vulnerable to pathogen attacks. Therefore, the leaf's performance and metabolic modes are inherently determined by its interaction with the environment. Mathematical models of plant metabolism have been successfully applied to study various aspects of photosynthesis, carbon and nitrogen assimilation and metabolism, aided suggesting metabolic intervention strategies for optimized leaf performance, and gave us insights into evolutionary drivers of plant metabolism in various environments. With the increasing pressure to improve agricultural performance in current and future climates, these models have become important tools to improve our understanding of plant–environment interactions and to propel plant breeders efforts. This overview article reviews applications of large-scale metabolic models of leaf metabolism to study plant–environment interactions by means of flux-balance analysis. The presented studies are organized in two ways — by the way the environment interactions are modelled — via external constraints or data-integration and by the studied environmental interactions — abiotic or biotic.

Overcoming the Challenges to Enhancing Experimental Plant Biology With Computational Modeling

The study of complex biological systems necessitates computational modeling approaches that are currently underutilized in plant biology. Many plant biologists have trouble identifying or adopting modeling methods to their research, particularly mechanistic mathematical modeling. Here we address challenges that limit the use of computational modeling methods, particularly mechanistic mathematical modeling. We divide computational modeling techniques into either pattern models (e.g., bioinformatics, machine learning, or morphology) or mechanistic mathematical models (e.g., biochemical reactions, biophysics, or population models), which both contribute to plant biology research at different scales to answer different research questions. We present arguments and recommendations for the increased adoption of modeling by plant biologists interested in incorporating more modeling into their research programs. As some researchers find math and quantitative methods to be an obstacle to modeling, we provide suggestions for easy-to-use tools for non-specialists and for collaboration with specialists. This may especially be the case for mechanistic mathematical modeling, and we spend some extra time discussing this. Through a more thorough appreciation and awareness of the power of different kinds of modeling in plant biology, we hope to facilitate interdisciplinary, transformative research.

Cherry Tomato Production in Intelligent Greenhouses-Sensors and AI for Control of Climate, Irrigation, Crop Yield, and Quality

Greenhouses and indoor farming systems play an important role in providing fresh and nutritious food for the growing global population. Farms are becoming larger and greenhouse growers need to make complex decisions to maximize production and minimize resource use while meeting market requirements. However, highly skilled labor is increasingly lacking in the greenhouse sector. Moreover, extreme events such as the COVID-19 pandemic, can make farms temporarily less accessible. This highlights the need for more autonomous and remote-control strategies for greenhouse production. This paper describes and analyzes the results of the second “Autonomous Greenhouse Challenge”. In this challenge, an experiment was conducted in six high-tech greenhouse compartments during a period of six months of cherry tomato growing. The primary goal of the greenhouse operation was to maximize net profit, by controlling the greenhouse climate and crop with AI techniques. Five international teams with backgrounds in AI and horticulture were challenged in a competition to operate their own compartment remotely. They developed intelligent algorithms and use sensor data to determine climate setpoints and crop management strategy. All AI supported teams outperformed a human-operated greenhouse that served as reference. From the results obtained by the teams and from the analysis of the different climate-crop strategies, it was possible to detect challenges and opportunities for the future implementation of remote-control systems in greenhouse production.

Voxel carving-based 3D reconstruction of sorghum identifies genetic determinants of light interception efficiency

Changes in canopy architecture traits have been shown to contribute to yield increases. Optimizing both light interception and light interception efficiency of agricultural crop canopies will be essential to meeting the growing food needs. Canopy architecture is inherently three-dimensional (3D), but many approaches to measuring canopy architecture component traits treat the canopy as a two-dimensional (2D) structure to make large scale measurement, selective breeding, and gene identification logistically feasible. We develop a high throughput voxel carving strategy to reconstruct 3D representations of sorghum from a small number of RGB photos. Our approach builds on the voxel carving algorithm to allow for fully automatic reconstruction of hundreds of plants. It was employed to generate 3D reconstructions of individual plants within a sorghum association population at the late vegetative stage of development. Light interception parameters estimated from these reconstructions enabled the identification of known and previously unreported loci controlling light interception efficiency in sorghum. The approach is generalizable and scalable, and it enables 3D reconstructions from existing plant high throughput phenotyping datasets. We also propose a set of best practices to increase 3D reconstructions' accuracy.

Simulating the effects of water limitation on plant biomass using a 3D functional-structural plant model of shoot and root driven by soil hydraulics

BACKGROUND AND AIMS

Improved modelling of carbon assimilation and plant growth to low soil moisture requires evaluation of underlying mechanisms in the soil, roots, and shoots. The feedback between plants and their local environment throughout the whole spectrum soil-root-shoot-environment is crucial to accurately describe and evaluate the impact of environmental changes on plant development. This study presents a 3D functional structural plant model, in which shoot and root growth are driven by radiative transfer, photosynthesis, and soil hydrodynamics through different parameterisation schemes relating soil water deficit and carbon assimilation. The new coupled model is used to evaluate the impact of soil moisture availability on plant productivity for two different groups of flowering plants under different spatial configurations.

METHODS

In order to address different aspects of plant development due to limited soil water availability, a 3D FSP model including root, shoot, and soil was constructed by linking three different well-stablished models of airborne plant, root architecture, and reactive transport in the soil. Different parameterisation schemes were used in order to integrate photosynthetic rate with root water uptake within the coupled model. The behaviour of the model was assessed on how the growth of two different types of plants, i.e. monocot and dicot, is impacted by soil water deficit under different competitive conditions: isolated (no competition), intra, and interspecific competition.

KEY RESULTS

The model proved to be capable of simulating carbon assimilation and plant development under different growing settings including isolated monocots and dicots, intra, and interspecific competition. The model predicted that (1) soil water availability has a larger impact on photosynthesis than on carbon allocation; (2) soil water deficit has an impact on root and shoot biomass production by up to 90 % for monocots and 50 % for dicots; and (3) the improved dicot biomass production in interspecific competition was highly related to root depth and plant transpiration.

CONCLUSIONS

An integrated model of 3D shoot architecture and biomass development with a 3D root system representation, including light limitation and water uptake considering soil hydraulics, was presented. Plant-plant competition and regulation on stomatal conductance to drought were able to be predicted by the model. In the cases evaluated here, water limitation impacted plant growth almost 10 times more than the light environment.

Two decades of functional-structural plant modelling: now addressing fundamental questions in systems biology and predictive ecology

BACKGROUND

Functional-structural plant models (FSPMs) explore and integrate relationships between a plant's structure and processes that underlie its growth and development. In the last 20 years, scientists interested in functional-structural plant modelling have expanded greatly the range of topics covered and now handle dynamical models of growth and development occurring from the microscopic scale, and involving cell division in plant meristems, to the macroscopic scales of whole plants and plant communities.

SCOPE

The FSPM approach occupies a central position in plant science; it is at the crossroads of fundamental questions in systems biology and predictive ecology. This special issue of Annals of Botany features selected papers on critical areas covered by FSPMs and examples of comprehensive models that are used to solve theoretical and applied questions, ranging from developmental biology to plant phenotyping and management of plants for agronomic purposes. Altogether, they offer an opportunity to assess the progress, gaps and bottlenecks along the research path originally foreseen for FSPMs two decades ago. This review also allows discussion of current challenges of FSPMs regarding (1) integration of multidisciplinary knowledge, (2) methods for handling complex models, (3) standards to achieve interoperability and greater genericity and (4) understanding of plant functioning across scales.

CONCLUSIONS

This approach has demonstrated considerable progress, but has yet to reach its full potential in terms of integration and heuristic knowledge production. The research agenda of functional-structural plant modellers in the coming years should place a greater emphasis on explaining robust emergent patterns, and on the causes of possible deviation from it. Modelling such patterns could indeed fuel both generic integration across scales and transdisciplinary transfer. In particular, it could be beneficial to emergent fields of research such as model-assisted phenotyping and predictive ecology in managed ecosystems.

Photosynthesis: basics, history and modelling

BACKGROUND

With limited agricultural land and increasing human population, it is essential to enhance overall photosynthesis and thus productivity. Oxygenic photosynthesis begins with light absorption, followed by excitation energy transfer to the reaction centres, primary photochemistry, electron and proton transport, NADPH and ATP synthesis, and then CO2 fixation (Calvin-Benson cycle, as well as Hatch-Slack cycle). Here we cover some of the discoveries related to this process, such as the existence of two light reactions and two photosystems connected by an electron transport 'chain' (the Z-scheme), chemiosmotic hypothesis for ATP synthesis, water oxidation clock for oxygen evolution, steps for carbon fixation, and finally the diverse mechanisms of regulatory processes, such as 'state transitions' and 'non-photochemical quenching' of the excited state of chlorophyll a.

SCOPE

In this review, we emphasize that mathematical modelling is a highly valuable tool in understanding and making predictions regarding photosynthesis. Different mathematical models have been used to examine current theories on diverse photosynthetic processes; these have been validated through simulation(s) of available experimental data, such as chlorophyll a fluorescence induction, measured with fluorometers using continuous (or modulated) exciting light, and absorbance changes at 820 nm (ΔA820) related to redox changes in P700, the reaction centre of photosystem I.

CONCLUSIONS

We highlight here the important role of modelling in deciphering and untangling complex photosynthesis processes taking place simultaneously, as well as in predicting possible ways to obtain higher biomass and productivity in plants, algae and cyanobacteria.

Multiscale computational models can guide experimentation and targeted measurements for crop improvement

Computational models of plants have identified gaps in our understanding of biological systems, and have revealed ways to optimize cellular processes or organ-level architecture to increase productivity. Thus, computational models are learning tools that help direct experimentation and measurements. Models are simplifications of complex systems, and often simulate specific processes at single scales (e.g. temporal, spatial, organizational, etc.). Consequently, single-scale models are unable to capture the critical cross-scale interactions that result in emergent properties of the system. In this perspective article, we contend that to accurately predict how a plant will respond in an untested environment, it is necessary to integrate mathematical models across biological scales. Computationally mimicking the flow of biological information from the genome to the phenome is an important step in discovering new experimental strategies to improve crops. A key challenge is to connect models across biological, temporal and computational (e.g. CPU versus GPU) scales, and then to visualize and interpret integrated model outputs. We address this challenge by describing the efforts of the international Crops in silico consortium.

Flux balance analysis of metabolism during growth by osmotic cell expansion and its application to tomato fruits

Cell expansion is a significant contributor to organ growth and is driven by the accumulation of osmolytes to increase cell turgor pressure. Metabolic modelling has the potential to provide insights into the processes that underpin osmolyte synthesis and transport, but the main computational approach for predicting metabolic network fluxes, flux balance analysis, often uses biomass composition as the main output constraint and ignores potential changes in cell volume. Here we present growth-by-osmotic-expansion flux balance analysis (GrOE-FBA), a framework that accounts for both the metabolic and ionic contributions to the osmotica that drive cell expansion, as well as the synthesis of protein, cell wall and cell membrane components required for cell enlargement. Using GrOE-FBA, the metabolic fluxes in dividing and expanding cells were analysed, and the energetic costs for metabolite biosynthesis and accumulation in the two scenarios were found to be surprisingly similar. The expansion phase of tomato fruit growth was also modelled using a multiphase single-optimization GrOE-FBA model and this approach gave accurate predictions of the major metabolite levels throughout fruit development, as well as revealing a role for transitory starch accumulation in ensuring optimal fruit development.

Towards a multiscale crop modelling framework for climate change adaptation assessment

Predicting the consequences of manipulating genotype (G) and agronomic management (M) on agricultural ecosystem performances under future environmental (E) conditions remains a challenge. Crop modelling has the potential to enable society to assess the efficacy of G × M technologies to mitigate and adapt crop production systems to climate change. Despite recent achievements, dedicated research to develop and improve modelling capabilities from gene to global scales is needed to provide guidance on designing G × M adaptation strategies with full consideration of their impacts on both crop productivity and ecosystem sustainability under varying climatic conditions. Opportunities to advance the multiscale crop modelling framework include representing crop genetic traits, interfacing crop models with large-scale models, improving the representation of physiological responses to climate change and management practices, closing data gaps and harnessing multisource data to improve model predictability and enable identification of emergent relationships. A fundamental challenge in multiscale prediction is the balance between process details required to assess the intervention and predictability of the system at the scales feasible to measure the impact. An advanced multiscale crop modelling framework will enable a gene-to-farm design of resilient and sustainable crop production systems under a changing climate at regional-to-global scales.

Using energy-efficient synthetic biochemical pathways to bypass photorespiration

Current crop yields will not be enough to sustain today's diets for a growing global population. As plant photosynthetic efficiency has not reached its theoretical maximum, optimizing photosynthesis is a promising strategy to enhance plant productivity. The low productivity of C3 plants is caused in part by the substantial energetic investments necessary to maintain a high flux through the photorespiratory pathway. Accordingly, lowering the energetic costs of photorespiration to enhance the productivity of C3 crops has been a goal of synthetic plant biology for decades. The use of synthetic bypasses to photorespiration in different plants showed an improvement of photosynthetic performance and growth under laboratory and field conditions, even though in silico predictions suggest that the tested synthetic pathways should confer a minimal or even negative energetic advantage over the wild type photorespiratory pathway. Current strategies increasingly utilize theoretical modeling and new molecular techniques to develop synthetic biochemical pathways that bypass photorespiration, representing a highly promising approach to enhance future plant productivity.

Crop Improvement from Phenotyping Roots: Highlights Reveal Expanding Opportunities

Root systems determine the water and nutrients for photosynthesis and harvested products, underpinning agricultural productivity. We highlight 11 programs that integrated root traits into germplasm for breeding, relying on phenotyping. Progress was successful but slow. Today's phenotyping technologies will speed up root trait improvement. They combine multiple new alleles in germplasm for target environments, in parallel. Roots and shoots are detected simultaneously and nondestructively, seed to seed measures are automated, and field and laboratory technologies are increasingly linked. Available simulation models can aid all phenotyping decisions. This century will see a shift from single root traits to rhizosphere selections that can be managed dynamically on farms and a shift to phenotype-based improvement to accommodate the dynamic complexity of whole crop systems.

Root Plasticity in the Pursuit of Water

One of the greatest challenges of terrestrial vegetation is to acquire water through soil-grown roots. Owing to the scarcity of high-quality water in the soil and the environment's spatial heterogeneity and temporal variability, ranging from extreme flooding to drought, roots have evolutionarily acquired tremendous plasticity regarding their geometric arrangement of individual roots and their three-dimensional organization within the soil. Water deficiency has also become an increasing threat to agriculture and dryland ecosystems due to climate change. As a result, roots have become important targets for genetic selection and modification in an effort to improve crop resilience under water-limiting conditions. This review addresses root plasticity from different angles: Their structures and geometry in response to the environment, potential genetic control of root traits suitable for water-limiting conditions, and contemporary and future studies of the principles underlying root plasticity post-Darwin's 'root-brain' hypothesis. Our increasing knowledge of different disciplines of plant sciences and agriculture should contribute to a sustainable management of natural and agricultural ecosystems for the future of mankind.

Practical steps to digital organism models, from laboratory model species to 'Crops in silico

A recent initiative named 'Crops in silico' proposes that multi-scale models 'have the potential to fill in missing mechanistic details and generate new hypotheses to prioritize directed engineering efforts' in plant science, particularly directed to crop species. To that end, the group called for 'a paradigm shift in plant modelling, from largely isolated efforts to a connected community'. 'Wet' (experimental) research has been especially productive in plant science, since the adoption of Arabidopsis thaliana as a laboratory model species allowed the emergence of an Arabidopsis research community. Parts of this community invested in 'dry' (theoretical) research, under the rubric of Systems Biology. Our past research combined concepts from Systems Biology and crop modelling. Here we outline the approaches that seem most relevant to connected, 'digital organism' initiatives. We illustrate the scale of experimental research required, by collecting the kinetic parameter values that are required for a quantitative, dynamic model of a gene regulatory network. By comparison with the Systems Biology Markup Language (SBML) community, we note computational resources and community structures that will help to realize the potential for plant Systems Biology to connect with a broader crop science community.

Functional phenomics: an emerging field integrating high-throughput phenotyping, physiology, and bioinformatics

The emergence of functional phenomics signifies the rebirth of physiology as a 21st century science through the use of advanced sensing technologies and big data analytics. Functional phenomics seeks to fill the significant knowledge gaps that still exist in the relationship of plant phenotype to function. Here, a general approach for the theory and practice of functional phenomics is outlined. The functional phenomics pipeline is proposed as a general method for conceptualizing, measuring, and validating utility of plant phenes, or elemental units of phenotype. The functional phenomics pipeline begins with ideotype development. Second, a phenotyping platform is developed to maximize the throughput of phene measurements. Target phenes and indicators of plant function, or performance, are measured in a mapping population. Forward genetics allows genetic mapping, while functional phenomics links phenes to plant performance. Based on these data, genotypes with contrasting phenotypes can be selected for smaller yet more intensive experiments to understand phene-environment interactions in depth. Simulation modeling is used to further understand the phenotypes, and all stages of the pipeline feed back to ideotype and phenotyping platform development. In total, functional phenomics represents an evolution of pre-existing disciplines, but the goals and unique methodologies constitute a novel research program.

Systems models, phenomics and genomics: three pillars for developing high-yielding photosynthetically efficient crops

Recent years witnessed a stagnation in yield enhancement in major staple crops, which leads plant biologists and breeders to focus on an urgent challenge to dramatically increase crop yield to meet the growing food demand. Systems models have started to show their capacity in guiding crops improvement for greater biomass and grain yield production. Here we argue that systems models, phenomics and genomics combined are three pillars for the future breeding for high-yielding photosynthetically efficient crops (HYPEC). Briefly, systems models can be used to guide identification of breeding targets for a particular cultivar and define optimal physiological and architectural parameters for a particular crop to achieve high yield under defined environments. Phenomics can support collection of architectural, physiological, biochemical and molecular parameters in a high-throughput manner, which can be used to support both model validation and model parameterization. Genomic techniques can be used to accelerate crop breeding by enabling more efficient mapping between genotypic and phenotypic variation, and guide genome engineering or editing for model-designed traits. In this paper, we elaborate on these roles and how they can work synergistically to support future HYPEC breeding.

Optimizing photorespiration for improved crop productivity

In C3 plants, photorespiration is an energy-expensive process, including the oxygenation of ribulose-1,5-bisphosphate (RuBP) by ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) and the ensuing multi-organellar photorespiratory pathway required to recycle the toxic byproducts and recapture a portion of the fixed carbon. Photorespiration significantly impacts crop productivity through reducing yields in C3 crops by as much as 50% under severe conditions. Thus, reducing the flux through, or improving the efficiency of photorespiration has the potential of large improvements in C3 crop productivity. Here, we review an array of approaches intended to engineer photorespiration in a range of plant systems with the goal of increasing crop productivity. Approaches include optimizing flux through the native photorespiratory pathway, installing non-native alternative photorespiratory pathways, and lowering or even eliminating Rubisco-catalyzed oxygenation of RuBP to reduce substrate entrance into the photorespiratory cycle. Some proposed designs have been successful at the proof of concept level. A plant systems-engineering approach, based on new opportunities available from synthetic biology to implement in silico designs, holds promise for further progress toward delivering more productive crops to farmer's fields.

Putting primary metabolism into perspective to obtain better fruits

Background

One of the key goals of fruit biology is to understand the factors that influence fruit growth and quality, ultimately with a view to manipulating them for improvement of fruit traits.

Scope

Primary metabolism, which is not only essential for growth but is also a major component of fruit quality, is an obvious target for improvement. However, metabolism is a moving target that undergoes marked changes throughout fruit growth and ripening.

Conclusions

Agricultural practice and breeding have successfully improved fruit metabolic traits, but both face the complexity of the interplay between development, metabolism and the environment. Thus, more fundamental knowledge is needed to identify further strategies for the manipulation of fruit metabolism. Nearly two decades of post-genomics approaches involving transcriptomics, proteomics and/or metabolomics have generated a lot of information about the behaviour of fruit metabolic networks. Today, the emergence of modelling tools is providing the opportunity to turn this information into a mechanistic understanding of fruits, and ultimately to design better fruits. Since high-quality data are a key requirement in modelling, a range of must-have parameters and variables is proposed.

Rightsizing root phenotypes for drought resistance

I propose that reduced root development would be advantageous for drought resistance in high-input agroecosystems. Selection regimes for crop ancestors and landraces include multiple stresses, intense competition, and variable resource distribution, which favored prolific root production, developmental plasticity in response to resource availability, and maintenance of unspecialized root tissues. High-input agroecosystems have removed many of these constraints to root function. Therefore, root phenotypes that focus on water capture at the expense of ancestral adaptations would be better suited to high-input agroecosystems. Parsimonious architectural phenotypes include fewer axial roots, reduced density of lateral roots, reduced growth responsiveness to local resource availability, and greater loss of roots that do not contribute to water capture. Parsimonious anatomical phenotypes include a reduced number of cortical cell files, greater loss of cortical parenchyma to aerenchyma and senescence, and larger cortical cell size. Parsimonious root phenotypes may be less useful in low-input agroecosystems, which are characterized by multiple challenges and trade-offs for root function in addition to water capture. Analysis of the fitness landscape of root phenotypes is a complex challenge that will be aided by the development of robust functional-structural models capable of simulating the dynamics of root-soil interactions.

Mathematical models of lignin biosynthesis

BACKGROUND

Lignin is a natural polymer that is interwoven with cellulose and hemicellulose within plant cell walls. Due to this molecular arrangement, lignin is a major contributor to the recalcitrance of plant materials with respect to the extraction of sugars and their fermentation into ethanol, butanol, and other potential bioenergy crops. The lignin biosynthetic pathway is similar, but not identical in different plant species. It is in each case comprised of a moderate number of enzymatic steps, but its responses to manipulations, such as gene knock-downs, are complicated by the fact that several of the key enzymes are involved in several reaction steps. This feature poses a challenge to bioenergy production, as it renders it difficult to select the most promising combinations of genetic manipulations for the optimization of lignin composition and amount.

RESULTS

Here, we present several computational models than can aid in the analysis of data characterizing lignin biosynthesis. While minimizing technical details, we focus on the questions of what types of data are particularly useful for modeling and what genuine benefits the biofuel researcher may gain from the resulting models. We demonstrate our analysis with mathematical models for black cottonwood (Populus trichocarpa), alfalfa (Medicago truncatula), switchgrass (Panicum virgatum) and the grass Brachypodium distachyon.

CONCLUSIONS

Despite commonality in pathway structure, different plant species show different regulatory features and distinct spatial and topological characteristics. The putative lignin biosynthes pathway is not able to explain the plant specific laboratory data, and the necessity of plant specific modeling should be heeded.

Allosteric Control of a Plant Receptor Kinase through S-Glutathionylation

Growing evidence supports the importance of protein S-glutathionylation as a regulatory post-translational modification with functional consequences for proteins. Discoveries of redox-state-dependent protein kinase S-glutathionylation have fueled discussion of redox-sensitive signaling. Following previously published experimental evidence for S-glutathionylation induced deactivation of the Arabidopsis thaliana kinase BRASSINOSTEROID INSENSITIVE 1-ASSOCIATED RECEPTOR-LIKE KINASE 1 (BAK1), we investigated the consequences of S-glutathionylation on the equilibrium conformational ensemble of BAK1 using all-atom molecular dynamics simulations. We found that glutathionylation of C408 allosterically destabilizes the active-like state of BAK1 and stabilizes an inactive conformation known to recur in protein kinases. Glutathionylation of C408 also has structural consequences throughout the BAK1 kinase domain, whereas glutathionylation of C353 in the N-lobe and C374 near the ATP-binding site have few notable effects on BAK1 compared with the unmodified protein. Our results suggest an allosteric mechanism for inhibition of BAK1 by C408 S-glutathionylation, and more generally, support the notion of protein kinase S-glutathionylation as a means of redox signaling in plant cells.