A mathematical model of biofilm growth and spread within plant xylem: Case study of Xylella fastidiosa in olive trees.

Xylem-limited bacterial pathogens cause some of the most destructive plant diseases. Though imposed measures to control these pathogens are generally ineffective, even among susceptible taxa, some hosts can limit bacterial loads and symptom expression. Mechanisms by which this resistance is achieved are poorly understood. In particular, it is still unknown how differences in vascular structure may influence biofilm growth and spread within a host. To address this, we developed a novel theoretical framework to describe biofilm behaviour within xylem vessels, adopting a polymer-based modelling approach. We then parameterised the model to investigate the relevance of xylem vessel diameters on Xylella fastidiosa resistance among olive cultivars. The functionality of all vessels was severely reduced under infection, with hydraulic flow reductions of 2-3 orders of magnitude. However, results suggest wider vessels act as biofilm incubators; allowing biofilms to develop over a long time while still transporting them through the vasculature. By contrast, thinner vessels become blocked much earlier, limiting biofilm spread. Using experimental data on vessel diameter distributions, we were able to determine that a mechanism of resistance in the olive cultivar Leccino is a relatively low abundance of the widest vessels, limiting X. fastidiosa spread.

A 3D image-based modelling approach for understanding spatiotemporal processes in phosphorus fertiliser dissolution, soil buffering and uptake by plant roots.

Phosphorus (P) is a key yield-limiting nutrient for crops, but the main source of P fertiliser is finite. Therefore, efficient fertilisation is crucial. Optimal P application requires understanding of the dynamic processes affecting P availability to plants, including fertiliser dissolution rate and soil buffer power. However, standard soil testing methods sample at fixed time points, preventing a mechanistic understanding of P uptake variability. We used image-based modelling to investigate the effects of fertiliser dissolution rate and soil buffer power on P uptake by wheat roots imaged using X-ray CT. We modelled uptake based on 1-day, 1-week, and 14-week dissolution of a fixed quantity of total P for two common soil buffer powers. We found rapid fertiliser dissolution increased short-term root uptake, but total uptake from 1-week matched 1-day dissolution. We quantified the large effects root system architecture had on P uptake, finding that there were trade-offs between total P uptake and uptake per unit root length, representing a carbon investment/phosphorus uptake balance. These results provide a starting point for predictive modelling of uptake from different P fertilisers in different soils. With the addition of further X-ray CT image datasets and a wider range of conditions, our simulation approach could be developed further for rapid trialling of fertiliser-soil combinations to inform field-scale trials or management.

Impact of root hairs on microscale soil physical properties in the field.

Aims: Recent laboratory studies revealed that root hairs may alter soil physical behaviour, influencing soil porosity and water retention on the small scale. However, the results are not consistent, and it is not known if structural changes at the small-scale have impacts at larger scales. Therefore, we evaluated the potential effects of root hairs on soil hydro-mechanical properties in the field using rhizosphere-scale physical measurements. Methods: Changes in soil water retention properties as well as mechanical and hydraulic characteristics were monitored in both silt loam and sandy loam soils. Measurements were taken from plant establishment to harvesting in field trials, comparing three barley genotypes representing distinct phenotypic categories in relation to root hair length. Soil hardness and elasticity were measured using a 3-mm-diameter spherical indenter, while water sorptivity and repellency were measured using a miniaturized infiltrometer with a 0.4-mm tip radius. Results: Over the growing season, plants induced changes in the soil water retention properties, with the plant available water increasing by 21%. Both soil hardness (P = 0.031) and elasticity (P = 0.048) decreased significantly in the presence of root hairs in silt loam soil, by 50% and 36%, respectively. Root hairs also led to significantly smaller water repellency (P = 0.007) in sandy loam soil vegetated with the hairy genotype (-49%) compared to the hairless mutant. Conclusions: Breeding of cash crops for improved soil conditions could be achieved by selecting root phenotypes that ameliorate soil physical properties and therefore contribute to increased soil health. Supplementary Information: The online version contains supplementary material available at 10.1007/s11104-022-05530-1.

Transition metal-catalysed carbene- and nitrene transfer to carbon monoxide and isocyanides.

Transition metal-catalysed carbene- and nitrene transfer to the C1-building blocks carbon monoxide and isocyanides provides heteroallenes (i.e. ketenes, isocyanates, ketenimines and carbodiimides). These are versatile and reactive compounds allowing in situ transformation towards numerous functional groups and organic compounds, including heterocycles. Both one-pot and tandem processes have been developed providing valuable synthetic methods for the organic chemistry toolbox. This review discusses all known transition metal-catalysed carbene- and nitrene transfer reactions towards carbon monoxide and isocyanides and in situ transformation of the heteroallenes hereby obtained, with a special focus on the general mechanistic considerations.

Quantifying citrate-enhanced phosphate root uptake using microdialysis.

AIMS: Organic acid exudation by plant roots is thought to promote phosphate (P) solubilisation and bioavailability in soils with poorly available nutrients. Here we describe a new combined experimental (microdialysis) and modelling approach to quantify citrate-enhanced P desorption and its importance for root P uptake. METHODS: To mimic the rhizosphere, microdialysis probes were placed in soil and perfused with citrate solutions (0.1, 1.0 and 10 mM) and the amount of P recovered from soil used to quantify rhizosphere P availability. Parameters in a mathematical model describing probe P uptake, citrate exudation, P movement and citrate-enhanced desorption were fit to the experimental data. These parameters were used in a model of a root which exuded citrate and absorbed P. The importance of soil citrate-P mobilisation for root P uptake was then quantified using this model. RESULTS: A plant needs to exude citrate at a rate of 0.73 µmol cm-1 of root h-1 to see a significant increase in P absorption. Microdialysis probes with citrate in the perfusate were shown to absorb similar quantities of P to an exuding root. CONCLUSION: A single root exuding citrate at a typical rate (4.3 × 10-5 µmol m-1 of root h-1) did not contribute significantly to P uptake. Microdialysis probes show promise for measuring rhizosphere processes when calibration experiments and mathematical modelling are used to decouple microdialysis and rhizosphere mechanisms.

Precipitation-optimised targeting of nitrogen fertilisers in a model maize cropping system.

Typically, half of the nitrogen (N) fertiliser applied to agricultural fields is lost to the wider environment. This inefficiency is driven by soil processes such as denitrification, volatilisation, surface run-off and leaching. Rainfall plays an important role in regulating these processes, ultimately governing when and where N fertiliser moves in soil and its susceptibility to gaseous loss. The interaction between rainfall, plant N uptake and N losses, however, remains poorly understood. In this study we use numerical modelling to predict the optimal N fertilisation strategy with respect to rainfall patterns and offer mechanistic explanations to the resultant differences in optimal times of fertiliser application. We developed a modelling framework that describes water and N transport in soil over a growing season and assesses nitrogen use efficiency (NUE) of split fertilisations within the context of different rainfall patterns. We used ninety rainfall patterns to determine their impact on optimal N fertilisation times. We considered the effects of root growth, root N uptake, microbial transformation of N and the effect of soil water saturation and flow on N movement in the soil profile. On average, we show that weather-optimised fertilisation strategies could improve crop N uptake by 20% compared to the mean uptake. In drier years, weather-optimising N applications improved the efficiency of crop N recovery by 35%. Further analysis shows that maximum plant N uptake is greatest under drier conditions due to reduced leaching, but it is harder to find the maximum due to low N mobility. The model could capture contrasting trends in NUE seen in previous arable cropping field trials. Furthermore, the model predicted that the variability in NUE seen in the field could be associated with precipitation-driven differences in N leaching and mobility. In conclusion, our results show that NUE in cropping systems could be significantly enhanced by synchronising fertiliser timings with both crop N demand and local weather patterns.

Significance of root hairs for plant performance under contrasting field conditions and water deficit.

BACKGROUND AND AIMS: Previous laboratory studies have suggested selection for root hair traits in future crop breeding to improve resource use efficiency and stress tolerance. However, data on the interplay between root hairs and open-field systems, under contrasting soils and climate conditions, are limited. As such, this study aims to experimentally elucidate some of the impacts that root hairs have on plant performance on a field scale. METHODS: A field experiment was set up in Scotland for two consecutive years, under contrasting climate conditions and different soil textures (i.e. clay loam vs. sandy loam). Five barley (Hordeum vulgare) genotypes exhibiting variation in root hair length and density were used in the study. Root hair length, density and rhizosheath weight were measured at several growth stages, as well as shoot biomass, plant water status, shoot phosphorus (P) accumulation and grain yield. KEY RESULTS: Measurements of root hair density, length and its correlation with rhizosheath weight highlighted trait robustness in the field under variable environmental conditions, although significant variations were found between soil textures as the growing season progressed. Root hairs did not confer a notable advantage to barley under optimal conditions, but under soil water deficit root hairs enhanced plant water status and stress tolerance resulting in a less negative leaf water potential and lower leaf abscisic acid concentration, while promoting shoot P accumulation. Furthermore, the presence of root hairs did not decrease yield under optimal conditions, while root hairs enhanced yield stability under drought. CONCLUSIONS: Selecting for beneficial root hair traits can enhance yield stability without diminishing yield potential, overcoming the breeder's dilemma of trying to simultaneously enhance both productivity and resilience. Therefore, the maintenance or enhancement of root hairs can represent a key trait for breeding the next generation of crops for improved drought tolerance in relation to climate change.

Mechanisms of root reinforcement in soils: an experimental methodology using four-dimensional X-ray computed tomography and digital volume correlation.

Vegetation on railway or highway slopes can improve slope stability through the generation of soil pore water suctions by plant transpiration and mechanical soil reinforcement by the roots. To incorporate the enhanced shearing resistance and stiffness of root-reinforced soils in stability calculations, it is necessary to understand and quantify its effectiveness. This requires integrated and sophisticated experimental and multi-scale modelling approaches to develop an understanding of the processes at different length scales, from individual root-soil interaction through to full soil-profile or slope scale. One of the challenges with multi-scale models is ensuring that they sufficiently closely represent real behaviour. This requires calibration against detailed high-quality and data-rich experiments. This study presents a novel experimental methodology, which combines in situ direct shear loading of a willow root-reinforced soil with X-ray computed tomography to capture the three-dimensional chronology of soil and root deformation within the shear zone. Digital volume correlation (DVC) analysis was applied to the computed tomography dataset to obtain full-field three-dimensional displacement and strain information. This paper demonstrates the feasibility and discusses the challenges associated with DVC experiments on root-reinforced soils.

Image-based quantification of soil microbial dead zones induced by nitrogen fertilization.

Microbial communities in agricultural soils underpin many ecosystem services including the maintenance of soil structure, food production, water purification and carbon storage. However, the impact of fertilization on the health of microbial communities is not well understood. This study investigates the spatial and temporal dynamics of nitrogen (N) transport away from a fertilizer granule with pore scale resolution. Specifically, we examined how soil structure and moisture content influence fertilizer derived N movement through the soil pore network and the subsequent impact of on soil microbial communities. We develop a mathematical model to describe N transport and reactions in soil at the pore-scale. Using X-ray Computed Tomography scans, we reconstructed a microscale description of a soil-pore geometry as a computational mesh. Solving two-phase water/air model produced pore-scale water distributions at 15, 30 and 70% water-filled pore volume. The N-speciation model considered ammonium (NH4+), nitrate (NO3-) and dissolved organic N (DON), and included N immobilization, ammonification and nitrification processes, as well as diffusion in soil solution. We simulated the dissolution of a fertilizer pellet and a pore scale N cycle at three different water saturations. To aid interpretation of the model results, microbial activity at a range of N concentrations was measured. The model showed that the diffusion and concentration of N in water films is critically dependent upon soil moisture and N species. We predict that the maximum NH4+ and NO3- concentrations in soil solution around the pellet under dry conditions are in the order of 1 × 103 and 1 × 104 mol m-3 respectively, and under wet conditions 2 × 102 and 1 × 103 mol m-3, respectively. Supporting experimental evidence suggests that these concentrations would be sufficient to reduce microbial activity in the short-term in the zone immediately around the fertilizer pellet (ranging from 0.9 to 3.8 mm), causing a major loss of soil biological functioning. This model demonstrates the importance of pore-scale processes in regulating N movement and their interactions with the soil microbiome.

Significance of root hairs at the field scale - modelling root water and phosphorus uptake under different field conditions.

BACKGROUND AND AIMS: Root hairs play a significant role in phosphorus (P) extraction at the pore scale. However, their importance at the field scale remains poorly understood. METHODS: This study uses a continuum model to explore the impact of root hairs on the large-scale uptake of P, comparing root hair influence under different agricultural scenarios. High vs low and constant vs decaying P concentrations down the soil profile are considered, along with early vs late precipitation scenarios. RESULTS: Simulation results suggest root hairs accounted for 50% of total P uptake by plants. Furthermore, a delayed initiation time of precipitation potentially limits the P uptake rate by over 50% depending on the growth period. Despite the large differences in the uptake rate, changes in the soil P concentration in the domain due to root solute uptake remains marginal when considering a single growth season. However, over the duration of 6 years, simulation results showed that noticeable differences arise over time. CONCLUSION: Root hairs are critical to P capture, with uptake efficiency potentially enhanced by coordinating irrigation with P application during earlier growth stages of crops.

Surface tension, rheology and hydrophobicity of rhizodeposits and seed mucilage influence soil water retention and hysteresis.

AIMS: Rhizodeposits collected from hydroponic solutions with roots of maize and barley, and seed mucilage washed from chia, were added to soil to measure their impact on water retention and hysteresis in a sandy loam soil at a range of concentrations. We test the hypothesis that the effect of plant exudates and mucilages on hydraulic properties of soils depends on their physicochemical characteristics and origin. METHODS: Surface tension and viscosity of the exudate solutions were measured using the Du Noüy ring method and a cone-plate rheometer, respectively. The contact angle of water on exudate treated soil was measured with the sessile drop method. Water retention and hysteresis were measured by equilibrating soil samples, treated with exudates and mucilages at 0.46 and 4.6 mg g-1 concentration, on dialysis tubing filled with polyethylene glycol (PEG) solution of known osmotic potential. RESULTS: Surface tension decreased and viscosity increased with increasing concentration of the exudates and mucilage in solutions. Change in surface tension and viscosity was greatest for chia seed exudate and least for barley root exudate. Contact angle increased with increasing maize root and chia seed exudate concentration in soil, but not barley root. Chia seed mucilage and maize root rhizodeposits enhanced soil water retention and increased hysteresis index, whereas barley root rhizodeposits decreased soil water retention and the hysteresis effect. The impact of exudates and mucilages on soil water retention almost ceased when approaching wilting point at -1500 kPa matric potential. CONCLUSIONS: Barley rhizodeposits behaved as surfactants, drying the rhizosphere at smaller suctions. Chia seed mucilage and maize root rhizodeposits behaved as hydrogels that hold more water in the rhizosphere, but with slower rewetting and greater hysteresis.

The effect of root exudates on rhizosphere water dynamics.

Most water and nutrients essential for plant growth travel across a thin zone of soil at the interface between roots and soil, termed the rhizosphere. Chemicals exuded by plant roots can alter the fluid properties, such as viscosity, of the water phase, potentially with impacts on plant productivity and stress tolerance. In this paper, we study the effects of plant exudates on the macroscale properties of water movement in soil. Our starting point is a microscale description of two fluid flow and exudate diffusion in a periodic geometry composed from a regular repetition of a unit cell. Using multiscale homogenization theory, we derive a coupled set of equations that describe the movement of air and water, and the diffusion of plant exudates on the macroscale. These equations are parametrized by a set of cell problems that capture the flow behaviour. The mathematical steps are validated by comparing the resulting homogenized equations to the original pore scale equations, and we show that the difference between the two models is ≲7% for eight cells. The resulting equations provide a computationally efficient method to study plant-soil interactions. This will increase our ability to predict how contrasting root exudation patterns may influence crop uptake of water and nutrients.

Determination of macro-scale soil properties from pore scale structures: image-based modelling of poroelastic structures.

We show how a combination of X-ray computed tomography (X-CT) and image-based modelling can be used to calculate the effect of moisture content and compaction on the macroscopic structural properties of soil. Our method is based on the equations derived in Daly & Roose (2018 Proc. R. Soc. A474, 20170141. (doi:10.1098/rspa.2017.0141)), which we have extended so they can be directly applied to the segmented images obtained from X-CT. We assume that the soils are composed of air-filled pore space, solid mineral grains and a mixed phase composed of both clay particles and water. We considered three different initial soil treatments, composed of two different compaction levels and two different moisture contents. We found that the effective properties of the soils were unaffected by compaction over the range tested in this paper. However, changing the moisture content significantly altered the hydraulic and mechanical properties of the soils. A key strength of this method is that it enables the optimization or even design of soils composed from different constituents, with specific mechanical and hydraulic properties.

Using high resolution X-ray computed tomography to create an image based model of a lymph node.

Lymph nodes are an important part of the immune system. They filter the lymphatic fluid as it is transported from the tissues before being returned to the blood stream. The fluid flow through the nodes influences the behaviour of the immune cells that gather within the nodes and the structure of the node itself. Measuring the fluid flow in lymph nodes experimentally is challenging due to their small size and fragility. In this paper, we present high resolution X-ray computed tomography images of a murine lymph node. The impact of the resulting visualized structures on fluid transport are investigated using an image based model. The high contrast between different structures within the lymph node provided by phase contrast X-ray computed tomography reconstruction results in images that, when related to the permeability of the lymph node tissue, suggest an increased fluid velocity through the interstitial channels in the lymph node tissue. Fluid taking a direct path from the afferent to the efferent lymphatic vessel, through the centre of the node, moved faster than the fluid that flowed around the periphery of the lymph node. This is a possible mechanism for particles being moved into the cortex.

Determination of macro-scale soil properties from pore-scale structures: model derivation.

In this paper, we use homogenization to derive a set of macro-scale poro-elastic equations for soils composed of rigid solid particles, air-filled pore space and a poro-elastic mixed phase. We consider the derivation in the limit of large deformation and show that by solving representative problems on the micro-scale we can parametrize the macro-scale equations. To validate the homogenization procedure, we compare the predictions of the homogenized equations with those of the full equations for a range of different geometries and material properties. We show that the results differ by [Formula: see text] for all cases considered. The success of the homogenization scheme means that it can be used to determine the macro-scale poro-elastic properties of soils from the underlying structure. Hence, it will prove a valuable tool in both characterization and optimization.

Fluid flow in porous media using image-based modelling to parametrize Richards' equation.

The parameters in Richards' equation are usually calculated from experimentally measured values of the soil-water characteristic curve and saturated hydraulic conductivity. The complex pore structures that often occur in porous media complicate such parametrization due to hysteresis between wetting and drying and the effects of tortuosity. Rather than estimate the parameters in Richards' equation from these indirect measurements, image-based modelling is used to investigate the relationship between the pore structure and the parameters. A three-dimensional, X-ray computed tomography image stack of a soil sample with voxel resolution of 6 µm has been used to create a computational mesh. The Cahn-Hilliard-Stokes equations for two-fluid flow, in this case water and air, were applied to this mesh and solved using the finite-element method in COMSOL Multiphysics. The upscaled parameters in Richards' equation are then obtained via homogenization. The effect on the soil-water retention curve due to three different contact angles, 0°, 20° and 60°, was also investigated. The results show that the pore structure affects the properties of the flow on the large scale, and different contact angles can change the parameters for Richards' equation.

Plant exudates may stabilize or weaken soil depending on species, origin and time.

We hypothesized that plant exudates could either gel or disperse soil depending on their chemical characteristics. Barley (Hordeum vulgare L. cv. Optic) and maize (Zea mays L. cv. Freya) root exudates were collected using an aerated hydroponic method and compared with chia (Salvia hispanica L.) seed exudate, a commonly used root exudate analogue. Sandy loam soil was passed through a 500-µm mesh and treated with each exudate at a concentration of 4.6 mg exudate g-1 dry soil. Two sets of soil samples were prepared. One set of treated soil samples was maintained at 4°C to suppress microbial processes. To characterize the effect of decomposition, the second set of samples was incubated at 16°C for 2 weeks at -30 kPa matric potential. Gas chromatography-mass spectrometry (GC-MS) analysis of the exudates showed that barley had the largest organic acid content and chia the largest content of sugars (polysaccharide-derived or free), and maize was in between barley and chia. Yield stress of amended soil samples was measured by an oscillatory strain sweep test with a cone plate rheometer. When microbial decomposition was suppressed at 4°C, yield stress increased 20-fold for chia seed exudate and twofold for maize root exudate compared with the control, whereas for barley root exudate decreased to half. The yield stress after 2 weeks of incubation compared with soil with suppressed microbial decomposition increased by 85% for barley root exudate, but for chia and maize it decreased by 87 and 54%, respectively. Barley root exudation might therefore disperse soil and this could facilitate nutrient release. The maize root and chia seed exudates gelled soil, which could create a more stable soil structure around roots or seeds. HIGHLIGHTS: Rheological measurements quantified physical behaviour of plant exudates and effect on soil stabilization.Barley root exudates dispersed soil, which could release nutrients and carbon.Maize root and chia seed exudates had a stabilizing effect on soil.Physical engineering of soil in contact with plant roots depends on the nature and origin of exudates.

Measurement of micro-scale soil deformation around roots using four-dimensional synchrotron tomography and image correlation.

This study applied time lapse (four-dimensional) synchrotron X-ray computed tomography to observe micro-scale interactions between plant roots and soil. Functionally contrasting maize root tips were repeatedly imaged during ingress into soil columns of varying water content and compaction. This yielded sequences of three-dimensional densiometric data, representing time-resolved geometric soil and root configurations at the micronmetre scale. These data were used as inputs for two full-field kinematic quantification methods, which enabled the analysis of three-dimensional soil deformation around elongating roots. Discrete object tracking was used to track rigid mineral grains, while continuum digital volume correlation was used to track grey-level patterns within local sub-volumes. These techniques both allowed full-field soil displacements to be quantified at an intra-rhizosphere spatial sampling scale of less than 300 µm. Significant differences in deformation mechanisms were identified around different phenotypes under different soil conditions. A uniquely strong contrast was observed between intact and de-capped roots grown in dry, compacted soil. This provides evidence that functional traits of the root cap significantly reduce the amount of soil disturbance per unit of root elongation, with this effect being particularly significant in drier soil.

Investigation of microvascular morphological measures for skeletal muscle tissue oxygenation by image-based modelling in three dimensions.

The supply of oxygen in sufficient quantity is vital for the correct functioning of all organs in the human body, especially for skeletal muscle during exercise. Traditionally, microvascular oxygen supply capability is assessed by the analysis of morphological measures on transverse cross-sections of muscle, e.g. capillary density or capillary-to-fibre ratio. In this work, we investigate the relationship between microvascular structure and muscle tissue oxygenation in mice. Phase contrast imaging was performed using synchrotron radiation computed tomography (SR CT) to visualize red blood cells (RBCs) within the microvasculature in mouse soleus muscle. Image-based mathematical modelling of the oxygen diffusion from the RBCs into the muscle tissue was subsequently performed, as well as a morphometric analysis of the microvasculature. The mean tissue oxygenation was then compared with the morphological measures of the microvasculature. RBC volume fraction and spacing (mean distance of any point in tissue to the closest RBC) emerged as the best predictors for muscle tissue oxygenation, followed by length density (summed RBC length over muscle volume). The two-dimensional measures of capillary density and capillary-to-fibre ratio ranked last. We, therefore, conclude that, in order to assess the states of health of muscle tissue, it is advisable to rely on three-dimensional morphological measures rather than on the traditional two-dimensional measures.

Phase contrast synchrotron radiation computed tomography of muscle spindles in the mouse soleus muscle.

Muscle spindles are skeletal muscle sensory organs involved in the sensation of position and movement of the body. We have explored the capability of phase contrast computed tomography to visualise muscle spindles in murine skeletal muscle. In particular, we have validated the visualisation of nerve fibres through phase contrast computed tomography using light microscopy on stained histological sections. We further present the first three-dimensional visualisation of muscle spindles in mouse soleus skeletal muscle in conjunction with the neurovascular bundle associated with it.