Are high-throughput root phenotyping platforms suitable for informing root system architecture models with genotype-specific parameters? An evaluation based on the root model ArchiSimple and a small panel of wheat cultivars.

Given the difficulties in accessing plant roots in situ, high-throughput root phenotyping (HTRP) platforms under controlled conditions have been developed to meet the growing demand for characterizing root system architecture (RSA) for genetic analyses. However, a proper evaluation of their capacity to provide the same estimates for strictly identical root traits across platforms has never been achieved. In this study, we performed such an evaluation based on six major parameters of the RSA model ArchiSimple, using a diversity panel of 14 bread wheat cultivars in two HTRP platforms that had different growth media and non-destructive imaging systems together with a conventional set-up that had a solid growth medium and destructive sampling. Significant effects of the experimental set-up were found for all the parameters and no significant correlations across the diversity panel among the three set-ups could be detected. Differences in temperature, irradiance, and/or the medium in which the plants were growing might partly explain both the differences in the parameter values across the experiments as well as the genotype × set-up interactions. Furthermore, the values and the rankings across genotypes of only a subset of parameters were conserved between contrasting growth stages. As the parameters chosen for our analysis are root traits that have strong impacts on RSA and are close to parameters used in a majority of RSA models, our results highlight the need to carefully consider both developmental and environmental drivers in root phenomics studies.

Identification of quantitative trait loci controlling nitrogen use efficiency-related traits in rice at the seedling stage under salt condition by genome-wide association study.

Rice cultivation is facing both salt intrusion and overuse of nitrogen fertilizers. Hence, breeding new varieties aiming to improve nitrogen use efficiency (NUE), especially under salt conditions, is indispensable. We selected 2,391 rice accessions from the 3K Rice Genomes Project to evaluate the dry weight under two N concentrations [2.86 mM - standard N (SN), and 0.36 mM - low N (LN)] crossed with two NaCl concentrations [0 (0Na) and 60 mM (60Na)] at the seedling stage. Genome-wide association studies for shoot, root, and plant dry weight (DW) were carried out. A total of 55 QTLs - 32, 16, and 7 in the whole, indica, and japonica panel - associated with one of the tested traits were identified. Among these, 27 QTLs co-localized with previously identified QTLs for DW-related traits while the other 28 were newly detected; 24, 8, 11, and 4 QTLs were detected in SN-0Na, LN-0Na, SN-60Na, and LN-60Na, respectively, and the remaining 8 QTLs were for the relative plant DW between treatments. Three of the 11 QTLs in SN-60Na were close to the regions containing three QTLs detected in SN-0Na. Eleven candidate genes for eight important QTLs were identified. Only one of them was detected in both SN-0Na and SN-60Na, while 5, 0, 3, and 2 candidate genes were identified only once under SN-0Na, LN-0Na, SN-60Na, and LN-60Na, respectively. The identified QTLs and genes provide useful materials and genetic information for future functional characterization and genetic improvement of NUE in rice, especially under salt conditions.

Hydraulic flux-responsive hormone redistribution determines root branching.

Plant roots exhibit plasticity in their branching patterns to forage efficiently for heterogeneously distributed resources, such as soil water. The xerobranching response represses lateral root formation when roots lose contact with water. Here, we show that xerobranching is regulated by radial movement of the phloem-derived hormone abscisic acid, which disrupts intercellular communication between inner and outer cell layers through plasmodesmata. Closure of these intercellular pores disrupts the inward movement of the hormone signal auxin, blocking lateral root branching. Once root tips regain contact with moisture, the abscisic acid response rapidly attenuates. Our study reveals how roots adapt their branching pattern to heterogeneous soil water conditions by linking changes in hydraulic flux with dynamic hormone redistribution.

Physiological adaptive traits are a potential allele reservoir for maize genetic progress under challenging conditions.

Combined phenomic and genomic approaches are required to evaluate the margin of progress of breeding strategies. Here, we analyze 65 years of genetic progress in maize yield, which was similar (101 kg ha-1 year-1) across most frequent environmental scenarios in the European growing area. Yield gains were linked to physiologically simple traits (plant phenology and architecture) which indirectly affected reproductive development and light interception in all studied environments, marked by significant genomic signatures of selection. Conversely, studied physiological processes involved in stress adaptation remained phenotypically unchanged (e.g. stomatal conductance and growth sensitivity to drought) and showed no signatures of selection. By selecting for yield, breeders indirectly selected traits with stable effects on yield, but not physiological traits whose effects on yield can be positive or negative depending on environmental conditions. Because yield stability under climate change is desirable, novel breeding strategies may be needed for exploiting alleles governing physiological adaptive traits.

Investigating Soil-Root Interactions with the Numerical Model R-SWMS.

In this chapter, we present the Root and Soil Water Movement and Solute transport model R-SWMS, which can be used to simulate flow and transport in the soil-plant system. The equations describing water flow in soil-root systems are presented and numerical solutions are provided. An application of R-SWMS is then briefly discussed, in which we combine in vivo and in silico experiments in order to decrypt water flow in the soil-root domain. More precisely, light transmission imaging experiments were conducted to generate data that can serve as input for the R-SWMS model. These data include the root system architecture, the soil hydraulic properties and the environmental conditions (initial soil water content and boundary conditions, BC). Root hydraulic properties were not acquired experimentally, but set to theoretical values found in the literature. In order to validate the results obtained by the model, the simulated and experimental water content distributions were compared. The model was then used to estimate variables that were not experimentally accessible, such as the actual root water uptake distribution and xylem water potential.

Non-invasive hydrodynamic imaging in plant roots at cellular resolution.

A key impediment to studying water-related mechanisms in plants is the inability to non-invasively image water fluxes in cells at high temporal and spatial resolution. Here, we report that Raman microspectroscopy, complemented by hydrodynamic modelling, can achieve this goal - monitoring hydrodynamics within living root tissues at cell- and sub-second-scale resolutions. Raman imaging of water-transporting xylem vessels in Arabidopsis thaliana mutant roots reveals faster xylem water transport in endodermal diffusion barrier mutants. Furthermore, transverse line scans across the root suggest water transported via the root xylem does not re-enter outer root tissues nor the surrounding soil when en-route to shoot tissues if endodermal diffusion barriers are intact, thereby separating 'two water worlds'.

Modification of the Expression of the Aquaporin ZmPIP2;5 Affects Water Relations and Plant Growth.

The plasma membrane intrinsic protein PIP2;5 is the most highly expressed aquaporin in maize (Zea mays) roots. Here, we investigated how deregulation of PIP2;5 expression affects water relations and growth using maize overexpression (OE; B104 inbred) or knockout (KO; W22 inbred) lines. The hydraulic conductivity of the cortex cells of roots grown hydroponically was higher in PIP2;5 OE and lower in pip2;5 KO lines compared with the corresponding wild-type plants. While whole-root conductivity decreased in the KO lines compared to the wild type, no difference was observed in OE plants. This paradox was interpreted using the MECHA hydraulic model, which computes the radial flow of water within root sections. The model hints that the plasma membrane permeability of the cells is not radially uniform but that PIP2;5 may be saturated in cell layers with apoplastic barriers, i.e. the endodermis and exodermis, suggesting the presence of posttranslational mechanisms controlling the abundance of PIP in the plasma membrane in these cells. At the leaf level, where the PIP2;5 gene is weakly expressed in wild-type plants, the hydraulic conductance was higher in the PIP2;5 OE lines compared with the wild-type plants, whereas no difference was observed in the pip2;5 KO lines. The temporal trend of leaf elongation rate, used as a proxy for that of xylem water potential, was faster in PIP2;5 OE plants upon mild stress, but not in well-watered conditions, demonstrating that PIP2;5 may play a beneficial role in plant growth under specific conditions.

Lateral Roots: Random Diversity in Adversity.

Lateral roots are essential for soil foraging and uptake of minerals and water. They feature a large morphological diversity that results from divergent primordia or root growth and development patterns. Besides a structured diversity, resulting from the hierarchical and developmental organization of root systems, there exists a random diversity, occurring between roots of similar age, of the same hierarchical order, and exposed to uniform conditions. The physiological bases and functional consequences of this random diversity are largely ignored. Here we review the evidence for such random diversity throughout the plant kingdom, present innovative approaches based on statistical modeling to account for such diversity, and set the list of its potential benefits in front of a variable and unpredictable soil environment.

Connecting the dots between computational tools to analyse soil-root water relations.

In recent years, many computational tools, such as image analysis, data management, process-based simulation, and upscaling tools, have been developed to help quantify and understand water flow in the soil-root system, at multiple scales (tissue, organ, plant, and population). Several of these tools work together or at least are compatible. However, for the uninformed researcher, they might seem disconnected, forming an unclear and disorganized succession of tools. In this article, we show how different studies can be further developed by connecting them to analyse soil-root water relations in a comprehensive and structured network. This 'explicit network of soil-root computational tools' informs readers about existing tools and helps them understand how their data (past and future) might fit within the network. We also demonstrate the novel possibilities of scale-consistent parameterizations made possible by the network with a set of case studies from the literature. Finally, we discuss existing gaps in the network and how we can move forward to fill them.

Going with the Flow: Multiscale Insights into the Composite Nature of Water Transport in Roots.

As water often limits crop production, a more complete understanding of plant water capture and transport is necessary. Here, we developed MECHA, a mathematical model that computes the flow of water across the root at the scale of walls, membranes, and plasmodesmata of individual cells, and used it to test hypotheses related to root water transport in maize (Zea mays). The model uses detailed root anatomical descriptions and a minimal set of experimental cell properties, including the conductivity of plasma membranes, cell walls, and plasmodesmata, which yield quantitative and scale-consistent estimations of water pathways and root radial hydraulic conductivity (k r). MECHA revealed that the mainstream hydraulic theories derived independently at the cell and root segment scales are compatible only if osmotic potentials within the apoplastic domains are uniform. The results suggested that the convection-diffusion of apoplastic solutes explained most of the offset between estimated k r in pressure clamp and osmotic experiments, while the contribution of water-filled intercellular spaces was limited. Furthermore, sensitivity analyses quantified the relative impact of cortex and endodermis cell conductivity of plasma membranes on root k r and suggested that only the latter contributed substantially to k r due to the composite nature of water flow across roots. The explicit root hydraulic anatomy framework brings insights into contradictory interpretations of experiments from the literature and suggests experiments to efficiently address questions pertaining to root water relations. Its scale consistency opens avenues for cross-scale communication in the world of root hydraulics.

The Xerobranching Response Represses Lateral Root Formation When Roots Are Not in Contact with Water.

Efficient soil exploration by roots represents an important target for crop improvement and food security [1, 2]. Lateral root (LR) formation is a key trait for optimizing soil foraging for crucial resources such as water and nutrients. Here, we report an adaptive response termed xerobranching, exhibited by cereal roots, that represses branching when root tips are not in contact with wet soil. Non-invasive X-ray microCT imaging revealed that cereal roots rapidly repress LR formation as they enter an air space within a soil profile and are no longer in contact with water. Transcript profiling of cereal root tips revealed that transient water deficit triggers the abscisic acid (ABA) response pathway. In agreement with this, exogenous ABA treatment can mimic repression of LR formation under transient water deficit. Genetic analysis in Arabidopsis revealed that ABA repression of LR formation requires the PYR/PYL/RCAR-dependent signaling pathway. Our findings suggest that ABA acts as the key signal regulating xerobranching. We conclude that this new ABA-dependent adaptive mechanism allows roots to rapidly respond to changes in water availability in their local micro-environment and to use internal resources efficiently.

Advanced Backcross QTL Analysis of Fiber Strength and Fineness in a Cross between Gossypium hirsutum and G. mustelinum.

The molecular genetic basis of cotton fiber strength and fineness in crosses between Gossypium mustelinum and Gossypium hirsutum (Upland cotton) was dissected using 21 BC3F2 and 12 corresponding BC3F2:3 and BC3F2:4 families. The BC3F2 families were genotyped with simple sequence repeat markers from a G. hirsutum by G. mustelinum linkage map, and the three generations of BC3-derived families were phenotyped for fiber strength (STR) and fineness (Micronaire, MIC). A total of 42 quantitative trait loci (QTLs) were identified through one-way analysis of variance, including 15 QTLs for STR and 27 for MIC, with the percentage of variance explained by individual loci averaging 13.86 and 14.06%, respectively. Eighteen of the 42 QTLs were detected at least twice near the same markers in different generations/families or near linked markers in the same family, and 28 of the 42 QTLs were identified in both mixed model-based composite interval mapping and one-way variance analyses. Alleles from G. mustelinum increased STR for eight of 15 and reduced MIC for 15 of 27 QTLs. Significant among-family genotypic effects (P < 0.001) were detected in 13 and 10 loci for STR and MIC respectively, and five loci showed significant (P < 0.001) genotype × family interaction for MIC. These results support the hypothesis that fiber quality improvement for Upland cotton could be realized by introgressing G. mustelinum alleles although complexities due to the different effects of genetic background on introgressed chromatin might be faced. Building on prior work with G. barbadense, G. tomentosum, and G. darwinii, QTL mapping involving introgression of G. mustelinum alleles offers new allelic variation to Upland cotton germplasm.

Segregation distortion and genome-wide digenic interactions affect transmission of introgressed chromatin from wild cotton species.

KEY MESSAGE: This study reports transmission genetics of chromosomal segments into Gossypium hirsutum from its most distant euploid relative, Gossypium mustelinum . Mutilocus interactions and structural rearrangements affect introgression and segregation of donor chromatin. Wild allotetraploid relatives of cotton are a rich source of genetic diversity that can be used in genetic improvement, but linkage drag and non-Mendelian transmission genetics are prevalent in interspecific crosses. These problems necessitate knowledge of transmission patterns of chromatin from wild donor species in cultivated recipient species. From an interspecific cross, Gossypium hirsutum × Gossypium mustelinum, we studied G. mustelinum (the most distant tetraploid relative of Upland cotton) allele retention in 35 BC3F1 plants and segregation patterns in BC3F2 populations totaling 3202 individuals, using 216 DNA marker loci. The average retention of donor alleles across BC3F1 plants was higher than expected and the average frequency of G. mustelinum alleles in BC3F2 segregating families was less than expected. Despite surprisingly high retention of G. mustelinum alleles in BC3F1, 46 genomic regions showed no introgression. Regions on chromosomes 3 and 15 lacking introgression were closely associated with possible small inversions previously reported. Nonlinear two-locus interactions are abundant among loci with single-locus segregation distortion, and among loci originating from one of the two subgenomes. Comparison of the present results with those of prior studies indicates different permeability of Upland cotton for donor chromatin from different allotetraploid relatives. Different contributions of subgenomes to two-locus interactions suggest different fates of subgenomes in the evolution of allotetraploid cottons. Transmission genetics of G. hirsutum × G. mustelinum crosses reveals allelic interactions, constraints on fixation and selection of donor alleles, and challenges with retention of introgressed chromatin for crop improvement.

QTL analysis of cotton fiber length in advanced backcross populations derived from a cross between Gossypium hirsutum and G. mustelinum.

KEY MESSAGE: QTLs for fiber length mapped in three generations of advanced backcross populations derived from crossing Gossypium hirsutum and Gossypium mustelinum showed opportunities to improve elite cottons by introgression from wild relatives. The molecular basis of cotton fiber length in crosses between Gossypium hirsutum and Gossypium mustelinum was dissected using 21 BC3F2 and 12 corresponding BC3F2:3 and BC3F2:4 families. Sixty-five quantitative trait loci (QTLs) were detected by one-way analysis of variance. The QTL numbers detected for upper-half mean length (UHM), fiber uniformity index (UI), and short fiber content (SFC) were 19, 20, and 26 respectively. Twenty-three of the 65 QTLs could be detected at least twice near adjacent markers in the same family or near the same markers across different families/generations, and 32 QTLs were detected in both one-way variance analyses and mixed model-based composite interval mapping. G. mustelinum alleles increased UHM and UI and decreased SFC for five, one, and one QTLs, respectively. In addition to the main-effect QTLs, 17 epistatic QTLs were detected which helped to elucidate the genetic basis of cotton fiber length. Significant among-family genotypic effects were detected at 18, 16, and 16 loci for UHM, UI, and SFC, respectively. Six, two, and two loci showed genotype × family interaction for UHM, UI and SFC, respectively, illustrating complexities that might be faced in introgression of exotic germplasm into cultivated cotton. Co-location of many QTLs for UHM, UI, and SFC accounted for correlations among these traits, and selection of these QTLs may improve the three traits simultaneously. The simple sequence repeat (SSR) markers associated with G. mustelinum QTLs will assist breeders in transferring and maintaining valuable traits from this exotic source during cultivar development.

Gravimetric phenotyping of whole plant transpiration responses to atmospheric vapour pressure deficit identifies genotypic variation in water use efficiency.

There is increasing interest in rapidly identifying genotypes with improved water use efficiency, exemplified by the development of whole plant phenotyping platforms that automatically measure plant growth and water use. Transpirational responses to atmospheric vapour pressure deficit (VPD) and whole plant water use efficiency (WUE, defined as the accumulation of above ground biomass per unit of water used) were measured in 100 maize (Zea mays L.) genotypes. Using a glasshouse based phenotyping platform with naturally varying VPD (1.5-3.8kPa), a 2-fold variation in WUE was identified in well-watered plants. Regression analysis of transpiration versus VPD under these conditions, and subsequent whole plant gas exchange at imposed VPDs (0.8-3.4kPa) showed identical responses in specific genotypes. Genotype response of transpiration versus VPD fell into two categories: 1) a linear increase in transpiration rate with VPD with low (high WUE) or high (low WUE) transpiration rate at all VPDs, 2) a non-linear response with a pronounced change point at low VPD (high WUE) or high VPD (low WUE). In the latter group, high WUE genotypes required a significantly lower VPD before transpiration was restricted, and had a significantly lower rate of transpiration in response to VPD after this point, when compared to low WUE genotypes. Change point values were significantly positively correlated with stomatal sensitivity to VPD. A change point in stomatal response to VPD may explain why some genotypes show contradictory WUE rankings according to whether they are measured under glasshouse or field conditions. Furthermore, this novel use of a high throughput phenotyping platform successfully reproduced the gas exchange responses of individuals measured in whole plant chambers, accelerating the identification of plants with high WUE.

Determination by near infrared microscopy of the nitrogen and carbon content of tomato (Solanum lycopersicum L.) leaf powder.

Near infrared microscopy (NIRM) has been developed as a rapid technique to predict the chemical composition of foods, reduce analytical costs and time and ease sample preparation. In this study, NIRM has been evaluated as an alternative to classical chemical analysis to determine the nitrogen and carbon content of small samples of tomato (Solanum lycopersicum L.) leaf powder. Near infrared spectra were obtained by NIRM for independent leaf samples collected on 216 plants grown under six different levels of nitrogen. From these, 30 calibration and 30 validation samples covering the spectral range of the whole set were selected and their nitrogen and carbon contents were determined by a reference method. The calibration model obtained for nitrogen content proved to be excellent, with a coefficient of determination in calibration (R(2)c) higher than 0.9 and a ratio of performance to deviation (RPDc) higher than 3. Statistical indicators of prediction using the validation set were also very high (R(2)p values > 0.90). However, the calibration model obtained for carbon content was much less satisfactory (R(2)c < 0.50). NIRM appears as a promising and suitable tool for a rapid, non-destructive and reliable determination of nitrogen content of tiny samples of tomato leaf powder.

Genome-Wide Analysis of Yield in Europe: Allelic Effects Vary with Drought and Heat Scenarios.

Assessing the genetic variability of plant performance under heat and drought scenarios can contribute to reduce the negative effects of climate change. We propose here an approach that consisted of (1) clustering time courses of environmental variables simulated by a crop model in current (35 years × 55 sites) and future conditions into six scenarios of temperature and water deficit as experienced by maize (Zea mays L.) plants; (2) performing 29 field experiments in contrasting conditions across Europe with 244 maize hybrids; (3) assigning individual experiments to scenarios based on environmental conditions as measured in each field experiment; frequencies of temperature scenarios in our experiments corresponded to future heat scenarios (+5°C); (4) analyzing the genetic variation of plant performance for each environmental scenario. Forty-eight quantitative trait loci (QTLs) of yield were identified by association genetics using a multi-environment multi-locus model. Eight and twelve QTLs were associated to tolerances to heat and drought stresses because they were specific to hot and dry scenarios, respectively, with low or even negative allelic effects in favorable scenarios. Twenty-four QTLs improved yield in favorable conditions but showed nonsignificant effects under stress; they were therefore associated with higher sensitivity. Our approach showed a pattern of QTL effects expressed as functions of environmental variables and scenarios, allowing us to suggest hypotheses for mechanisms and candidate genes underlying each QTL. It can be used for assessing the performance of genotypes and the contribution of genomic regions under current and future stress situations and to accelerate breeding for drought-prone environments.

RBOH-mediated ROS production facilitates lateral root emergence in Arabidopsis.

Lateral root (LR) emergence represents a highly coordinated process in which the plant hormone auxin plays a central role. Reactive oxygen species (ROS) have been proposed to function as important signals during auxin-regulated LR formation; however, their mode of action is poorly understood. Here, we report that Arabidopsis roots exposed to ROS show increased LR numbers due to the activation of LR pre-branch sites and LR primordia (LRP). Strikingly, ROS treatment can also restore LR formation in pCASP1:shy2-2 and aux1 lax3 mutant lines in which auxin-mediated cell wall accommodation and remodeling in cells overlying the sites of LR formation is disrupted. Specifically, ROS are deposited in the apoplast of these cells during LR emergence, following a spatiotemporal pattern that overlaps the combined expression domains of extracellular ROS donors of the RESPIRATORY BURST OXIDASE HOMOLOGS (RBOH). We also show that disrupting (or enhancing) expression of RBOH in LRP and/or overlying root tissues decelerates (or accelerates) the development and emergence of LRs. We conclude that RBOH-mediated ROS production facilitates LR outgrowth by promoting cell wall remodeling of overlying parental tissues.

Analysis of root growth from a phenotyping data set using a density-based model.

Major research efforts are targeting the improved performance of root systems for more efficient use of water and nutrients by crops. However, characterizing root system architecture (RSA) is challenging, because roots are difficult objects to observe and analyse. A model-based analysis of RSA traits from phenotyping image data is presented. The model can successfully back-calculate growth parameters without the need to measure individual roots. The mathematical model uses partial differential equations to describe root system development. Methods based on kernel estimators were used to quantify root density distributions from experimental image data, and different optimization approaches to parameterize the model were tested. The model was tested on root images of a set of 89 Brassica rapa L. individuals of the same genotype grown for 14 d after sowing on blue filter paper. Optimized root growth parameters enabled the final (modelled) length of the main root axes to be matched within 1% of their mean values observed in experiments. Parameterized values for elongation rates were within ±4% of the values measured directly on images. Future work should investigate the time dependency of growth parameters using time-lapse image data. The approach is a potentially powerful quantitative technique for identifying crop genotypes with more efficient root systems, using (even incomplete) data from high-throughput phenotyping systems.

Loss of function of 1-FEH IIb has more impact on post-harvest inulin degradation in Cichorium intybus than copy number variation of its close paralog 1-FEH IIa.

Key Message: The loss of mini-exon 2 in the 1-FEH IIb glycosyl-hydrolase results in a putative non-functional allele. This loss of function has a strong impact on the susceptibility to post-harvest inulin depolymerization. Significant variation of copy number was identified in its close paralog 1-FEH IIa, but no quantitative effect of copy number on carbohydrates-related phenotypes was detected. Inulin polyfructan is the second most abundant storage carbohydrate in flowering plants. After harvest, it is depolymerized by fructan exohydrolases (FEHs) as an adaptive response to end-season cold temperatures. In chicory, the intensity of this depolymerization differs between cultivars but also between individuals within a cultivar. Regarding this phenotypic variability, we recently identified statistically significant associations between inulin degradation and genetic polymorphisms located in three FEHs. We present here new results of a systematic analysis of copy number variation (CNV) in five key members of the chicory (Cichorium intybus) GH32 multigenic family, including three FEH genes and the two inulin biosynthesis genes: 1-SST and 1-FFT. qPCR analysis identified a significant variability of relative copy number only in the 1-FEH IIa gene. However, this CNV had no quantitative effect. Instead, cloning of the full length gDNA of a close paralogous sequence (1-FEH IIb) identified a 1028 bp deletion in lines less susceptible to post-harvest inulin depolymerization. This region comprises a 9 bp mini-exon containing one of the three conserved residues of the active site. This results in a putative non-functional 1-FEH IIb allele and an observed lower inulin depolymerization. Extensive genotyping confirmed that the loss of mini-exon 2 in 1-FEH IIb and the previously identified 47 bp duplication located in the 3'UTR of 1-FEH IIa belong to a single haplotype, both being statistically associated with reduced susceptibility to post-harvest inulin depolymerization. Emergence of these haplotypes is discussed.