Research Article | Open Access

Sam Sultan, Joseph Snider, Adam Conn, Mao Li, Christopher N. Topp, Saket Navlakha, "A Statistical Growth Property of Plant Root Architectures", *Plant Phenomics*, vol. 2020, Article ID 2073723, 11 pages, 2020. https://doi.org/10.34133/2020/2073723

# A Statistical Growth Property of Plant Root Architectures

#### Abstract

Numerous types of biological branching networks, with varying shapes and sizes, are used to acquire and distribute resources. Here, we show that plant root and shoot architectures share a fundamental design property. We studied the spatial density function of plant architectures, which specifies the probability of finding a branch at each location in the 3-dimensional volume occupied by the plant. We analyzed 1645 root architectures from four species and discovered that the spatial density functions of all architectures are population-similar. This means that despite their apparent visual diversity, all of the roots studied share the same basic shape, aside from stretching and compression along orthogonal directions. Moreover, the spatial density of all architectures can be described as variations on a single underlying function: a Gaussian density truncated at a boundary of roughly three standard deviations. Thus, the root density of any architecture requires only four parameters to specify: the total mass of the architecture and the standard deviations of the Gaussian in the three growth directions. Plant shoot architectures also follow this design form, suggesting that two basic plant transport systems may use similar growth strategies.

#### 1. Introduction

Understanding how biological branching networks grow and distribute branches in space has long fascinated mathematical biologists [1]. These networks serve as the basis for acquiring and distributing resources, for example, neural arbors that process information [2–4], vascular networks that circulate blood [5, 6], and plant architectures that transport sugars and nutrients [7–9]. Due to basic similarities in their design and function [10, 11], it is natural to ask whether there are quantitative features shared by several of these networks [8, 12].

Many structural features of plant root architectures have been studied, including root length and root depth [13], the distribution of root hairs [14–16], the size and number of lateral branches [17–20], and biomass allocation to fine roots [21]. More global properties of root shapes have also been analyzed [22], including allometric and fractal scaling [23–25], root foraging precision [26, 27], and topological morphology via persistent homology methods [28]. These properties can affect numerous biological functions performed by root architectures, including anchorage, carbon sequestration, and search for water and nutrients in the soil [29–33]. Identifying the molecular mechanisms that drive the formation of different shapes can aid in uncovering genotype-to-phenotype relationships [34] and in breeding specific traits of interest in crops [35–39].

We focus here on branch density [40, 41] and the spatial distribution of roots [42–44], which have numerous functional implications. Water and nutrients in the soil are ephemeral, distributed heterogeneously within a complex chemical and physical environment and intensely competed for by other organisms. Thus, when and where a plant grows its roots have strong implications for its success in resource capture and therefore its viability. In this foraging process, a plant’s limited carbon resources are continually partitioned between the growth of existing roots (elongation) and the generation of new roots (branching), which is conditioned by genetic mechanisms and both internal and external signals. The outcome of these iterative processes at any point in time is reflected by the global properties of the root system architecture, which is a complex 3D shape intimately tied to root system function.

Spatial distributions of root length, biomass, and density reflect the topological and geometric properties of root system architecture [45] and have been used extensively as quantitative summaries to compare the genetic, environmental, and genetic x environmental differences in how roots forage. Significant understanding of root system architecture has developed around mathematical models that fit observations of root systems to spatial density functions, beginning with the simple linear model of Gerwitz and Page [46]. However, these models do not fit all data well; in part, because observations from the field have almost exclusively consisted of sparse, lower dimensional representations of much more complex 3D root structures. Classic and still widely-used methods include analysis of root fragments from slabs taken from trench or monolith excavations, or from soil cores or minirhizotron data [47]. Efforts to better approximate the true root structure from these limited samples have used angular distributions [48] or 3D structural simulations [49] to produce better quantitative relationships between the observable arrangement of roots as they intersect the 2D plane of a slab, to the actual density distribution in 3D space. Current architectural (reviewed by Dunbabin et al. [50]) and density continuum models (reviewed by Dupuy et al. [51]) have been greatly aided by continued improvements in root phenotyping and computer simulation, including the addition of soil chemical and physical parameters that shape outcomes of root structure and resource foraging functions [52, 53]. Furthermore, describing root systems in terms of spatial densities allow a direct comparison to shoot systems.

##### 1.1. Our Contributions

Here, we use tools from probability theory to study the spatial density function of root architectures, which describes the probability of finding a branch at each point in the 3D volumetric space occupied by the root system. Our goal is to uncover statistical properties of this function, to study how this function varies across 4 species comprising 10 genotypes, and to determine if there are features of this function shared by all architectures studied. Specifically, we analyzed 1645 3D root architectures from four species (rice, corn, tomato, and *Arabidopsis*), including several genotypes per species. Using statistical moments to characterize the spatial density function of each architecture, we found that all functions are population-similar, i.e., all 1645 architectures analyzed had the same underlying shape and thus could be superimposed on top of each other, modulo compression and stretching along each direction. We then derived an analytical form of the spatial density function and discovered that all root architectures can be closely approximated by a 3D Gaussian function truncated at approximately 3.3 standard deviations. This means that, for each architecture, root density is highest at the center of mass and then decays like a Gaussian function until the boundary of the root system is reached at 3.3 standard deviations from the center. The truncation corresponds to the boundary of the plant, outside of which there is zero density, since plants have a finite spatial domain. This result also implies that the root density of any architecture requires only four parameters to specify: the total mass of the architecture and the standard deviations of the Gaussian in the three growth directions. We also find species-specific variation in this property, which suggests functional specialization across genotypes despite following the same general design template. Finally, prior work found that a Gaussian function truncated at 2.0 standard deviations closely approximated the spatial density function of hundreds of 3D shoot architectures [54]. To our knowledge, this highlights one of the first shared structural features of plant architectures above and below ground.

#### 2. Results

##### 2.1. Dataset of 3-Dimensional Root System Architectures

Phenotyping root system architectures in their natural habitat remain technologically challenging due to soil and particulates that obfuscate visibility. To overcome these physical barriers, several phenotyping methods have been proposed, each of which employs a different trade-off in terms of accuracy, coverage, throughput, and cost [18, 34]. These methods are based, for example, on ground-penetrating radar [55–57], electric resistivity tomography [56, 58, 59], minirhizotrons [60–64], shovelomics [65, 66], soil coring [67], X-ray tomography [68–71], magnetic resonance imaging [42, 72, 73], and 3D laser scanning [74, 75].

Here, we used a controlled environment, gel-based optical imaging platform that allows the 3D architectures of freely-growing root systems to be measured noninvasively in large numbers [35, 75–78] (Supplementary Methods (available here)). This method outputs a point-cloud representation of an architecture (Figure 1(a)). Overall, we collected 1645 architectures from four species, with multiple genotypes per species.
(1)Rice [35]. We collected two genotypes of *Oryza sativa* (Bala and Azucena), as well as a third group representing the F6 recombinant inbred lines of the parental cross Azucena x Bala with 171 different families. The recombinant inbred lines provide additional genetic and phenotypic variation within the rice population. In total, there were 36, 29, and 421 individual Azucena, Bala, and Azucena x Bala plants, respectively. Most individual plants were imaged roughly three times on days 12, 14, and 16, respectively, totaling 1358 root architectures.(2)Corn. We collected four genotypes of maize (*Zea mays* L. ssp. *mays*): B73, Mo17, IHP, and ILP. In total, we analyzed 106 architectures, each captured on growth day 6.(3)*Arabidopsis thaliana* [79]. We collected two genotypes: Col-0 and tob1-1 (containing a mutation in an auxin biosynthesis pathway gene). In total, we analyzed 133 architectures representing 31 individual plants, each imaged on approximately days 9, 11, 13, 15, 17, 21, and 25.(4)Tomato. We collected two genotypes of *Solanum lycopersicum*: wildtype cultivar “Moneymaker” and an RNAi line in a nitrogen metabolism gene [80]. In total, we analyzed 48 architectures representing 19 individual plants, each imaged on approximately days 6, 8, 10, 12, and 14.

**(a) Root system architectures**

**(b) Gaussian spatial density function**

The plant species chosen represent two major forms of root systems—tap-rooted (tomato, Arabidopsis) and fibrous-rooted (corn, rice)—with fundamentally different developmental patterning. Further, three of the species are important agriculturally (tomato, rice, corn), and one is a traditional model system (*Arabidopsis*). Thus, these species are phylogenetically and developmentally diverse and relevant economically and scientifically. Moreover, the dataset contains multiple genotypes per species and multiple time-points per plant, and it contains architectures ranging over two orders of magnitude in size, from 2,011 to 317,933 cloud points. Thus, this dataset represents a strong benchmark for testing whether a statistical property is shared by a broad class of root architectures.

##### 2.2. Describing Root Architectures Using Statistical Moments

To characterize how roots are distributed in space, we studied the root system’s spatial density function (Figure 1(b)). This function describes the density of points in the 3D volumetric territory occupied by the roots. The methods used to study the spatial density function are briefly summarized below; full technical details can be found in the Supplement and in prior published reports [54, 81].

Finding properties of the spatial density function that can be compared across many diverse architectures is a formidable challenge due to noise and variability in architecture sizes. For example, grid-counting techniques, which define density as the total sum of points in each voxel, can be highly sensitive to noise and other shifts in form (e.g., rotations and translations). Further, using a small size for each voxel will generate density functions that are sparse and difficult to compare, whereas large-sized voxels may not provide sufficient spatial resolution to make detailed comparisons.

To overcome these challenges, we define the spatial density function by its statistical moments [81]. Knowing all of the moments of a distribution is equivalent to knowing the function that generates the distribution, and the accuracy of the description increases with more moments calculated [82]. Lower order moments capture coarse features of the architecture—e.g., the 0th moment is the total mass of the root system, the 1st moment divided by the total mass is the center of mass, and the 2nd moment is the rotational inertia—whereas higher order moments capture more intricate details, including the shape of individual branches. Further, using moments to study the spatial density function avoids reliance on template matching or other assumptions about the number of functional forms for architectures [83, 84].

To calculate moments, we start with a 3D point cloud of the architecture (Figure 1(a)). A point cloud is a set of points in a Cartesian coordinate system, where each point is the 3D position of the point on the surface or interior of the object (root system). Following standard probability theory, the product moment is defined as:
where denotes the center of mass of the root system in the -direction (i.e., the mean of the -coordinates of each point), and is the volume of a voxel in mm^{3}. The resolution of the camera system and the resulting density of points varied from root to root depending on the experimental setup. We calculate moments for even values of (to ensure all moments are positive) and thus leave out the absolute value sign. Following the definition, the moment is equal to the total mass —— which is used to normalize for size.

##### 2.3. Are Root System Architectures Population-Similar?

We first describe the goal of our analysis intuitively, and then, we formalize these notions mathematically.

Our first goal is to use architecture moments to test if all spatial density functions are population-similar. This asks: do all density functions have the same shape, modulo stretching and compression of the architecture along one or more directions? Whether two root systems share the same density function can be difficult to assess visually. For example, two root systems may have the same shape, but one may be stretched along the (up-down) direction to forage deeper in search for water, whereas another root system may develop lateral roots that stretch along the direction. Similarly, two architectures scaled to be the same size could superimpose exactly, but this may be difficult to tell when comparing a large versus a small root system.

Population similarity is different than the commonly studied property of self-similarity [1, 8, 84]. “Self-similarity” is often associated with having the property that a single structure looks similar at all magnification scales (e.g., it is fractal). Population-similarity considers not a single structure, but rather a population of structures, and asks whether all of them can be viewed as variations, via stretching or compression, of a single form [81]. In other words, different root spatial density functions are population-similar if they can be “transformed” into one another by expanding or contracting root density along orthogonal spatial dimensions [54].

The concept of equivalence across stretching translates mathematically in terms of a density function that measures the probability that an object of size has mass at 3D position . The density function depends explicitly on spatial location and the overall scale, , of the architecture. Without loss of generality, we conventionally define spatial location for each architecture with respect to its center of mass. However, we do not know the length scale, , of each architecture because, unless we are generating the architectures ourselves, we only measure the consequences of changing , not itself.

To rewrite everything in terms of measurable quantities, consider the mathematical definition of equivalence across stretching:

This formulation expresses all possible architectures in terms of a single reference architecture, , that is stretched (divided) by a factor of and scaled (multiplied) by a power law, . The key advantage of this formulation is the separation between how the architectures change size and shape, reflected in the power law and terms, respectively. As detailed above, the density function itself is experimentally inaccessible, but we can measure its moments from the point cloud data using Eqn. (1). Moments of the theoretical density function may also be calculated as where the information about population similarity is encoded in the first term, and the form of the reference architecture is encoded in the single parameter, . See Supplementary Information for further derivations of the population similarity test and power-law scaling.

Eqn. (3) forms the basis to determine the degree to which root spatial density functions are population similar. This involves two steps: (1)The first step is to plot, versus for various values of and for each architecture. If, for each value of , the two have a linear relationship (with a different slope for each ), then the architectures share the same function. The term represents a typical measure of size that is computable using only the moments themselves, and no other quantity (Supplementary Methods). This measure is also proportional to another common measure of size, the convex hull volume, i.e., the smallest convex polytope that encloses all the cloud points (Figure 2(a)).(2)The second step is to plot the slope of the lines generated in the first step versus , the moment order. The difference between the slope of this line and 1 denotes the degree to which the architectures are all population similar [81].

**(a)**

**(b)**

**(c)**

**(d)**

Even if roots are population-similar, they may be generated by different population-similar functions, such as a uniform, an exponential, or a Gaussian spatial density. Each species or genotype may also belong to its own functional class; e.g., tomato architectures may have a uniform spatial density, whereas corn architectures may have a Gaussian density. This test can determine how many different classes of population-similar functions are required to describe the root architectures because architectures from one class will fall on one line, and architectures for another class will fall on a different line, for each value of . Thus, this test quantifies the degree to which architectures are population-similar and the number of functional classes of architectures.

For the first test, we plotted versus for for all 1645 architectures (Figure 2(b)). Strikingly, we found that for each moment order , all architectures lied on the same line, suggesting that all architectures share the same functional form. As a reminder, the function describes the spatial density of roots, which is a probability distribution that specifies, at each point in 3D space, the probability of finding a branch at that point. Deviation from the line increases with higher moment orders because higher moment orders have a much larger range of values (from approximately to ), and because they capture increasingly fine branching structures. Nonetheless, each line closely approximated the data ( for each of the 11 least-squares regression lines).

For the second test, we plotted the slopes of the lines calculated in the first test versus , the moment order (Figure 2(c)). The slope of this line was , meaning that the root architectures deviate by only from being truly population-similar. We tested the robustness of this result in three ways. First, we repeated the second test using both even and odd moments between 0 and 20; this yielded a nearly identical slope of , suggesting that enough moments are being considered. Second, we computed the slope for all species together, while only sampling 100 random Azucena x Bala plants, instead of using all 1175 of them; we found a similar slope of , suggesting that the strong representation of Azucena x Bala plants in our dataset does not unduly bias our results. Third, we computed the slope for each species with at least 100 architectures independently and found some species-specific differences in scaling, as expected, but qualitatively similar results: *Arabidopsis* (), rice (), and corn ().

Notably, we grew tomato plants in “nonoptimal” growth conditions and found significant deviation from population similarity (Figure S1), suggesting that population similarity is not inevitable under every growth regime and may indeed be an actively generated property.

##### 2.4. What Is the Functional Form of Root Spatial Densities?

The tests above indicate that all 1645 architectures are population-similar and share a single density function. Next, we seek to find a function with few parameters that gives us a sufficient statistical description of the root spatial density function. For all population-similar functions, the plot of moment order versus slope (Figure 2(c)) will have a slope of 1 [81]. Thus, the only other parameter of the line (the intercept) must provide information about the actual form of the function.

To determine the form of the function, we compared the intercepts from Figure 2(c) to the intercepts of two simple density functions: a 3D uniform function and a 3D Gaussian function truncated at a spherical boundary. A uniform function implies that root density is equal everywhere within the volume occupied by the root system. A Gaussian function implies that root density is highest at the center of mass and then decays as you move outward, with the rate of decay controlled by the standard deviations in each orthogonal direction. The truncation (boundary) of both functions corresponds to the edge of the volume that the root system occupies, outside of which there is zero density.

Strikingly, we found that all the root architectures can be described as a 3D Gaussian truncated at approximately 3.3 standard deviations from the center of mass (Figure 2(d)). We compared the goodness of fit of this function by comparing the intercepts of the root architectures to the intercepts of the uniform function and to a Gaussian function with a higher and lower standard deviation—both of which are poorer fits (Figure 2(d)). Indeed, a difference of only one standard deviation (2.3 or 4.3) significantly departs from the data.

In summary, we derived a simple, statistical description of the spatial density function of all 1645 architectures. A truncated Gaussian density function (Figure 1(b)) implies that only 4 parameters are required to describe statistically how roots are distributed in space: the center of mass and the standard deviations in each of the three spatial growth directions.

##### 2.5. Is There Variation in This Property across Species or Genotypes?

Our results suggest a statistical similarity across many diverse architectures, but it is also to be expected that there exists some variation in spatial densities across species or genotypes. Differences in the apparent visual similarities of root architectures can have profoundly different impacts on a plant’s growth.

To explore this, we focused on scatter around the regression lines for each moment order to test if this scatter represented noise or systematic differences in architectures from different species or genotypes (Supplementary Methods). We studied scatter around the line, corresponding to the total mass of the root system () versus its volume (). Recall that this measure of volume is derived entirely from the moments and highly correlates with the convex hull volume (Figure 2(a)). Mass and volume showed a linear relationship on a log-log plot, indicating that as total mass increases, the volumetric space occupied increases according to power law. Thus, roots below the regression line have a smaller total mass for the same volumetric space, and roots above the regression line have a larger total mass for the same volume. Scatter in any of the 11 moment lines can be studied in this manner, but we chose mass and volume because they are relatively easy to interpret.

We found systematic differences in the mass-to-volume relationship across species and across some genotypes within the same species. For example, Figures 3(a) and 3(b) show that 100.0% of tomato WT/Moneymaker plants lay below the regression line, indicating that they had a smaller mass than other root systems, despite occupying the same volume. Figures 3(c) and 3(d) shows that 88.46% of the root systems from the corn B73 and Mo17 genotypes lie above the regression line, whereas the other two corn genotypes (IHP and ILP) lie roughly equally above and below the line. We found similar differences across genotypes of rice, with one genotype (family 182 of Azucena x Bala) laying entirely above the regression line and another genotype (family 16 of Azucena x Bala) laying almost entirely below the regression line (Figures 3(e) and 3(f)).

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

#### 3. Discussion

We studied the spatial density functions of 1645 root systems and found that all architectures were population-similar. Population similarity means that the root density function of all architectures can be transformed into one another by stretching or compressing along orthogonal spatial directions. We also found that a single function, a 3D Gaussian truncated at roughly 3.3 standard deviations, closely describes the spatial densities of the architectures studied. This, as opposed to having different functions for different genotypes or species, was not an expected result. There are clearly other, more complicated functions that may also fit this data, but our goal here was to find a simple function with relatively few parameters that provided a robust fit. Evolutionarily, developing the regulatory circuits to generate and fine-tune a single functional form appears to be parsimonious, and uncovering the molecular mechanisms driving this form is a natural next question.

The density property revealed in this work was identified by a large-scale analysis of entire, nonsimulated, 3D root system architectures across multiple species. Therefore, the growth property incorporates actual heterogeneities in root growth, since it was biologically produced. We pooled data from different species to test if this property may be shared by a diverse set of architectures and thus a useful root phenotyping trait; a per-species analysis, where each species was analyzed separately, also showed similar trends. While all root systems we examined fit under the Gaussian function truncated at approximately 3 s.d., the clustering of species and/or genotypes in groups that deviated from this function is of great interest and highlights potential differences in underlying genetic mechanisms that pattern growth. We consider this a new trait, or phenotype, from which to evaluate root system architecture comprehensively and to constrain architecture modeling schemes.

A previous analysis of 557 shoot architectures—comprising three species (tomato, tobacco, sorghum), each grown in several growth conditions (ambient light, high-light, high-heat, drought, shade), and across roughly 20 developmental time-points per plant—found that the spatial density of shoot architecture branches could also be described by a Gaussian function truncated at approximately 2.0 standard deviations [54], instead of 3.3 standard deviations found here. The truncation amount indicates how many standard deviations away from the center of mass there exists nonzero root density. The higher truncation means that root architectures are more spread-out (relative to their mass) than shoot architectures. There is substantial interest in understanding the growth strategies that root systems use to forage (e.g., [24]), and our results provide a new constraint on this process and perhaps a property common with shoot foraging.

Finally, there are many biological branching structures where Gaussian spatial densities are not observed (e.g., blood vessel networks [6], sand dune morphologies [85, 86], and dendrites of retinal ganglion cells [87]), indicating that these properties should not be considered the null hypothesis. Indeed, given that along with providing support, the main drive for a plant to construct roots is to gather nutrients from the soil, the most obvious design would be space filling (i.e., a uniform spatial density). When an engineer designs anchors for a building, they extend supports equally in a compact space. When a landscape designer spreads irrigation, they do so evenly over the space. Why then have evolutionary processes selected Gaussian distributions in these species, with density concentrated near the center and fewer roots spaced out over an extended region? One possibility is that Gaussian functions balance maximal reach of supporting structures with efficiently gathering resources, e.g., as observed in river systems [88]. Alternatively, the Gaussian densities may simply reflect the competition between roots for space since overlapping Gaussians generate a maximally “flat” combined distribution, as observed in retinal ganglion cells [89]. The well-known central limit theorem that generates Gaussians for random, independent processes does not apply here because the roots are self-avoiding and cannot overlap each other. In such a case, the expected distribution for a random process is strongly non-Gaussian [90]. While the details of the growth process remain to be discovered, our results show that the roots expend resources to array themselves in a specific structure.

Prior work also did not find a similarity between plant-fungus interaction networks below ground and plant-animal interaction networks above ground [91]. Of course, we do not expect every architecture to abide by this rule (e.g., vine-like root systems that do not branch), and it remains to be seen how these properties are affected by changing environmental conditions and across the thousands of other species of plants in the natural world. Nonetheless, this property appears to be relevant to at least 4 species—including two monocots and two dicots, comprising 10 genotypes—and thus may represent a broader principle of plant architecture design.

#### Data Availability

Point cloud data for all 1645 3D root architectures are available online at: http://plant3d.navlakhalab.net.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Authors’ Contributions

J.S., C.N.T., and S.N. conceived the study. All authors performed the analysis and wrote the manuscript. Sam Sultan and Joseph Snider are co-first authors.

#### Acknowledgments

The authors thank Zhengbin Liu for contributing unpublished data of tomato and corn root architectures. S.S. thanks Dr. and Mrs. Larry Greenfield for their generous lab support. S.N. was supported by the Pew Charitable Trusts, the National Science Foundation under award CAREER DBI-1846554, and the NIDCD of the National Institutes of Health under award number 1R01DC017695.

#### Supplementary Materials

Supplementary Methods. Mathematical definition of population similarity, Mathematical definition of power-law scaling, Example of deviation from population similarity, Analysis of skeleton architectures, Analysis of rotated architectures.* (Supplementary Materials)*

#### References

- B. Mandelbrot and M. Novak,
*Thinking in Patterns: Fractals and Related Phenomena in Nature*, World Scientic, 2004. - H. Cuntz, F. Forstner, A. Borst, and M. Hausser, “One rule to grow them all: a general theory of neuronal branching and its practical application,”
*PLoS Computational Biology*, vol. 6, no. 8, p. e1000877, 2010. View at: Publisher Site | Google Scholar - J. M. L. Budd, K. Kovács, A. S. Ferecskó, P. Buzás, U. T. Eysel, and Z. F. Kisvárday, “Neocortical axon arbors Trade-off material and conduction delay conservation,”
*PLoS Computational Biology*, vol. 6, no. 3, article e1000711, 2010. View at: Publisher Site | Google Scholar - A. Chandrasekhar and S. Navlakha, “Neural arbors are Pareto optimal,”
*Proceedings of the Biological Sciences*, vol. 286, no. 1902, article 20182727, 2019. View at: Publisher Site | Google Scholar - P. Blinder, A. Y. Shih, C. Rafie, and D. Kleinfeld, “Topological basis for the robust distribution of blood to rodent neocortex,”
*Proceedings od the National Academy of Sciences of the United States of America*, vol. 107, no. 28, pp. 12670–12675, 2010. View at: Publisher Site | Google Scholar - P. Blinder, P. S. Tsai, J. P. Kaufhold, P. M. Knutsen, H. Suhl, and D. Kleinfeld, “The cortical angiome: an interconnected vascular network with noncolumnar patterns of blood flow,”
*Nature Neuroscience*, vol. 16, no. 7, pp. 889–897, 2013. View at: Publisher Site | Google Scholar - P. Prusinkiewicz and A. Lindenmayer,
*The Algorithmic Beauty of Plants*, Springer-Verlag New York, Inc., New York, NY, USA, 1996. - G. B. West, J. H. Brown, and B. J. Enquist, “A general model for the structure and allometry of plant vascular systems,”
*Nature*, vol. 400, no. 6745, pp. 664–667, 1999. View at: Publisher Site | Google Scholar - K. D. Farnsworth and K. J. Niklas, “Theories of optimization, form and function in branching architecture in plants,”
*Functional Ecology*, vol. 9, no. 3, pp. 355–363, 1995. View at: Publisher Site | Google Scholar - A. Conn, U. V. Pedmale, J. Chory, and S. Navlakha, “High-resolution laser scanning reveals plant architectures that reflect universal network design principles,”
*Cell Systems*, vol. 5, no. 1, pp. 53–62.e3, 2017. View at: Publisher Site | Google Scholar - K. J. Niklas, “Evolutionary walks through a land plant morphospace,”
*Journal of Experimental Botany*, vol. 50, no. 330, pp. 39–52, 1999. View at: Publisher Site | Google Scholar - J. H. Brown, V. K. Gupta, B. L. Li, B. T. Milne, C. Restrepo, and G. B. West, “The fractal nature of nature: power laws, ecological complexity and biodiversity,”
*Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences*, vol. 357, no. 1421, pp. 619–626, 2002. View at: Publisher Site | Google Scholar - J. Lynch, “Root architecture and plant productivity,”
*Plant Physiology*, vol. 109, no. 1, pp. 7–13, 1995. View at: Publisher Site | Google Scholar - Y. J. Zhang, J. P. Lynch, and K. M. Brown, “Ethylene and phosphorus availability have interacting yet distinct effects on root hair development,”
*Journal of Experimental Botany*, vol. 54, no. 391, pp. 2351–2361, 2003. View at: Publisher Site | Google Scholar - Z. Ma, D. G. Bielenberg, K. M. Brown, and J. P. Lynch, “Regulation of root hair density by phosphorus availability inArabidopsis thaliana,”
*Plant, Cell & Environment*, vol. 24, no. 4, pp. 459–467, 2001. View at: Publisher Site | Google Scholar - N. Koebernick, K. R. Daly, S. D. Keyes et al., “High-resolution synchrotron imaging shows that root hairs influence rhizosphere soil structure formation,”
*The New Phytologist*, vol. 216, no. 1, pp. 124–135, 2017. View at: Publisher Site | Google Scholar - B. Péret, B. de Rybel, I. Casimiro et al., “Arabidopsis lateral root development: an emerging story,”
*Trends in Plant Science*, vol. 14, no. 7, pp. 399–408, 2009. View at: Publisher Site | Google Scholar - J. Zhu, P. A. Ingram, P. N. Benfey, and T. Elich, “From lab to field, new approaches to phenotyping root system architecture,”
*Current Opinion in Plant Biology*, vol. 14, no. 3, pp. 310–317, 2011. View at: Publisher Site | Google Scholar - J. P. Lynch and K. M. Brown, “Topsoil foraging - an architectural adaption of plants to low phosphorous availability,”
*Plant and Soil*, vol. 237, no. 2, pp. 225–237, 2001. View at: Publisher Site | Google Scholar - J. Zhu, S. M. Kaeppler, and J. P. Lynch, “Topsoil foraging and phosphorus acquisition efficiency in maize (Zea mays),”
*Functional Plant Biology*, vol. 32, no. 8, pp. 749–762, 2005. View at: Publisher Site | Google Scholar - H. Poorter, K. J. Niklas, P. B. Reich, J. Oleksyn, P. Poot, and L. Mommer, “Biomass allocation to leaves, stems and roots: meta-analyses of interspecific variation and environmental control,”
*The New Phytologist*, vol. 193, no. 1, pp. 30–50, 2012. View at: Publisher Site | Google Scholar - A. Bucksch, A. Atta-Boateng, A. F. Azihou et al., “Morphological plant modeling: unleashing geometric and topological potential within the plant sciences,”
*Frontiers in Plant Science*, vol. 8, p. 900, 2017. View at: Publisher Site | Google Scholar - G. M. Berntson, “Topological scaling and plant root system architecture: developmental and functional hierarchies,”
*New Phytologist*, vol. 135, no. 4, pp. 621–634, 1997. View at: Publisher Site | Google Scholar - H. Ozier-Lafontaine, F. Lecompte, and J. F. Sillon, “Fractal analysis of the root architecture of Gliricidia sepium for the spatial prediction of root branching, size and mass: model development and evaluation in agroforestry,”
*Plant and Soil*, vol. 209, no. 2, pp. 167–179, 1999. View at: Publisher Site | Google Scholar - M. van Noordwijk, L. Y. Spek, and P. de Willigen, “Proximal root diameter as predictor of total root size for fractal branching models,”
*Plant and Soil*, vol. 164, no. 1, pp. 107–117, 1994. View at: Publisher Site | Google Scholar - H. De Kroon, E. J. W. Visser, H. Huber, L. Mommer, and M. J. Hutchings, “A modular concept of plant foraging behaviour: the interplay between local responses and systemic control,”
*Plant, Cell & Environment*, vol. 32, no. 6, pp. 704–712, 2009. View at: Publisher Site | Google Scholar - S. W. Kembel and J. F. Cahill Jr., “Plant phenotypic plasticity belowground: a phylogenetic perspective on root foraging trade-offs,”
*The American Naturalist*, vol. 166, no. 2, pp. 216–230, 2005. View at: Publisher Site | Google Scholar - M. Li, K. Duncan, C. N. Topp, and D. H. Chitwood, “Persistent homology and the branching topologies of plants,”
*American Journal of Botany*, vol. 104, no. 3, pp. 349–353, 2017. View at: Publisher Site | Google Scholar - A. H. Fitter, “An architectural approach to the comparative ecology of plant root systems,”
*New Phytologist*, vol. 106, pp. 61–77, 1987. View at: Publisher Site | Google Scholar - D. Robinson, A. Hodge, B. S. Griffiths, and A. H. Fitter, “Plant root proliferation in nitrogen-rich patches confers competitive advantage,”
*Proceedings of the Royal Society of London. Series B: Biological Sciences*, vol. 266, no. 1418, pp. 431–435, 1999. View at: Publisher Site | Google Scholar - J. Lopez-Bucio, A. Cruz-Ramirez, and L. Herrera-Estrella, “The role of nutrient availability in regulating root architecture,”
*Current Opinion in Plant Biology*, vol. 6, no. 3, pp. 280–287, 2003. View at: Publisher Site | Google Scholar - D. B. Kell, “Breeding crop plants with deep roots: their role in sustainable carbon, nutrient and water sequestration,”
*Annals of Botany*, vol. 108, no. 3, pp. 407–418, 2011. View at: Publisher Site | Google Scholar - M. Balduzzi, B. M. Binder, A. Bucksch et al., “Reshaping plant biology: qualitative and quantitative descriptors for plant morphology,”
*Frontiers in Plant Science*, vol. 8, p. 117, 2017. View at: Publisher Site | Google Scholar - A. Paez-Garcia, C. Motes, W. R. Scheible, R. Chen, E. Blancaflor, and M. Monteros, “Root traits and phenotyping strategies for plant improvement,”
*Plants*, vol. 4, no. 2, pp. 334–355, 2015. View at: Publisher Site | Google Scholar - C. N. Topp, A. S. Iyer-Pascuzzi, J. T. Anderson et al., “3D phenotyping and quantitative trait locus mapping identify core regions of the rice genome controlling root architecture,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 110, no. 18, pp. E1695–E1704, 2013. View at: Publisher Site | Google Scholar - H. Tian, I. De Smet, and Z. Ding, “Shaping a root system: regulating lateral versus primary root growth,”
*Trends in Plant Science*, vol. 19, no. 7, pp. 426–431, 2014. View at: Publisher Site | Google Scholar - B. Hufnagel, S. M. Sousa, L. Assis et al., “Duplicate and conquer: multiple homologs of PHOSPHORUS-STARVATION TOLERANCE1 enhance phosphorus acquisition and sorghum performance on low-phosphorus soils,”
*Plant Physiology*, vol. 166, no. 2, pp. 659–677, 2014. View at: Publisher Site | Google Scholar - G. Wachsman, E. E. Sparks, and P. N. Benfey, “Genes and networks regulating root anatomy and architecture,”
*The New Phytologist*, vol. 208, no. 1, pp. 26–38, 2015. View at: Publisher Site | Google Scholar - S. H. Su, N. M. Gibbs, A. L. Jancewicz, and P. H. Masson, “Molecular mechanisms of root gravitropism,”
*Current Biology*, vol. 27, no. 17, pp. R964–R972, 2017. View at: Publisher Site | Google Scholar - S. V. Hatzig, S. Schiessl, A. Stahl, and R. J. Snowdon, “Characterizing root response phenotypes by neural network analysis,”
*Journal of Experimental Botany*, vol. 66, no. 18, pp. 5617–5624, 2015. View at: Publisher Site | Google Scholar - A. D. Bird and H. Cuntz, “Dissecting sholl analysis into its functional components,”
*Cell Reports*, vol. 27, no. 10, pp. 3081–3096.e5, 2019. View at: Publisher Site | Google Scholar - D. Dusschoten, R. Metzner, J. Kochs et al., “Quantitative 3D analysis of plant roots growing in soil using magnetic resonance imaging,”
*Plant Physiology*, vol. 170, no. 3, pp. 1176–1188, 2016. View at: Publisher Site | Google Scholar - R. J. Flavel, C. N. Guppy, S. M. R. Rabbi, and I. M. Young, “An image processing and analysis tool for identifying and analysing complex plant root systems in 3D soil using non-destructive analysis: Root1,”
*PLoS One*, vol. 12, no. 5, article e0176433, 2017. View at: Publisher Site | Google Scholar - F. Danjon, D. H. Barker, M. Drexhage, and A. Stokes, “Using three-dimensional plant root architecture in models of shallow-slope stability,”
*Annals of Botany*, vol. 101, no. 8, pp. 1281–1293, 2007. View at: Publisher Site | Google Scholar - P. Prusinkiewicz and P. B. de Reuille, “Constraints of space in plant development,”
*Journal of Experimental Botany*, vol. 61, no. 8, pp. 2117–2129, 2010. View at: Publisher Site | Google Scholar - A. Gerwitz and E. R. Page, “An empirical mathematical model to describe plant root systems,”
*The Journal of Applied Ecology*, vol. 11, no. 2, pp. 773–781, 1974. View at: Publisher Site | Google Scholar - W. Böhm, “Methods of studying root systems,” in
*Ecological Studies*, vol. 33, Springer-Verlag, 1979. View at: Google Scholar - A. R. G. Lang and F. M. Melhuish, “Lengths and diameters of plant roots in non-random populations by analysis of plane surfaces,”
*Biometrics*, vol. 26, no. 3, pp. 421–431, 1970. View at: Publisher Site | Google Scholar - P. Grabarnik, L. Pages, and A. G. Bengough, “Geometrical properties of simulated maize root systems: consequences for length density and intersection density,”
*Plant and Soil*, vol. 200, no. 2, pp. 157–167, 1998. View at: Publisher Site | Google Scholar - V. M. Dunbabin, J. A. Postma, A. Schnepf et al., “Modelling root–soil interactions using three–dimensional models of root growth, architecture and function,”
*Plant and Soil*, vol. 372, no. 1-2, pp. 93–124, 2013. View at: Publisher Site | Google Scholar - L. Dupuy, P. J. Gregory, and A. G. Bengough, “Root growth models: towards a new generation of continuous approaches,”
*Journal of Experimental Botany*, vol. 61, no. 8, pp. 2131–2143, 2010. View at: Publisher Site | Google Scholar - S. Schlüter, S. R. G. A. Blaser, M. Weber, V. Schmidt, and D. Vetterlein, “Quantification of root growth patterns from the soil perspective via root distance models,”
*Frontiers in Plant Science*, vol. 9, article 1084, 2018. View at: Publisher Site | Google Scholar - T. Roose, S. D. Keyes, K. R. Daly et al., “Challenges in imaging and predictive modeling of rhizosphere processes,”
*Plant and Soil*, vol. 407, no. 1-2, pp. 9–38, 2016. View at: Publisher Site | Google Scholar - A. Conn, U. V. Pedmale, J. Chory, C. F. Stevens, and S. Navlakha, “A statistical description of plant shoot architecture,”
*Current Biology*, vol. 27, no. 14, pp. 2078–2088.e3, 2017. View at: Publisher Site | Google Scholar - D. B. Stover, F. P. Day, J. R. Butnor, and B. G. Drake, “Effect of elevated CO2 on coarse-root biomass in Florida scrub detected by ground-penetrating radar,”
*Ecology*, vol. 88, no. 5, pp. 1328–1334, 2007. View at: Publisher Site | Google Scholar - T. Zenone, G. Morelli, M. Teobaldelli et al., “Preliminary use of ground-penetrating radar and electrical resistivity tomography to study tree roots in pine forests and poplar plantations,”
*Functional Plant Biology*, vol. 35, no. 10, pp. 1047–1058, 2008. View at: Publisher Site | Google Scholar - J. R. Butnor, J. A. Doolittle, K. H. Johnsen, L. Samuelson, T. Stokes, and L. Kress, “Utility of ground-penetrating radar as a root biomass survey tool in forest systems,”
*Soil Science Society of America Journal*, vol. 67, no. 5, pp. 1607–1615, 2003. View at: Publisher Site | Google Scholar - M. Amato, B. Basso, G. Celano, G. Bitella, G. Morelli, and R. Rossi, “In situ detection of tree root distribution and biomass by multi-electrode resistivity imaging,”
*Tree Physiology*, vol. 28, no. 10, pp. 1441–1448, 2008. View at: Publisher Site | Google Scholar - M. Amato, G. Bitella, R. Rossi, J. A. Gómez, S. Lovelli, and J. J. F. Gomes, “Multi-electrode 3D resistivity imaging of alfalfa root zone,”
*European Journal of Agronomy*, vol. 31, no. 4, pp. 213–222, 2009. View at: Publisher Site | Google Scholar - G. H. Bates, “A device for the observation of root growth in the soil,”
*Nature*, vol. 139, no. 3527, pp. 966-967, 1937. View at: Publisher Site | Google Scholar - J. L. Maeght, B. Rewald, and A. Pierret, “How to study deep roots-and why it matters,”
*Frontiers in Plant Science*, vol. 4, p. 229, 2013. View at: Publisher Site | Google Scholar - C. M. Iversen, M. T. Murphy, M. F. Allen et al., “Advancing the use of minirhizotrons in wetlands,”
*Plant and Soil*, vol. 352, no. 1-2, pp. 23–39, 2012. View at: Publisher Site | Google Scholar - R. L. Hendrick and K. S. Pregitzer, “Spatial variation in tree root distribution and growth associated with minirhizotrons,”
*Plant and Soil*, vol. 143, pp. 283–288, 1992. View at: Publisher Site | Google Scholar - R. Rellán-Álvarez, G. Lobet, H. Lindner et al., “GLO-Roots: an imaging platform enabling multidimensional characterization of soil-grown root systems,”
*eLife*, vol. 4, 2015. View at: Publisher Site | Google Scholar - A. Bucksch, J. Burridge, L. M. York et al., “Image-based high-throughput field phenotyping of crop roots,”
*Plant Physiology*, vol. 166, no. 2, pp. 470–486, 2014. View at: Publisher Site | Google Scholar - S. Trachsel, S. M. Kaeppler, K. M. Brown, and J. P. Lynch, “Shovelomics: high throughput phenotyping of maize (Zea Mays L.) root architecture in the field,”
*Plant and Soil*, vol. 341, no. 1-2, pp. 75–87, 2011. View at: Publisher Site | Google Scholar - A. P. Wasson, G. J. Rebetzke, J. A. Kirkegaard, J. Christopher, R. A. Richards, and M. Watt, “Soil coring at multiple field environments can directly quantify variation in deep root traits to select wheat genotypes for breeding,”
*Journal of Experimental Botany*, vol. 65, no. 21, pp. 6231–6249, 2014. View at: Publisher Site | Google Scholar - S. Mairhofer, S. Zappala, S. Tracy et al., “Recovering complete plant root system architectures from soil via X-ray
*μ*-Computed tomography,”*Plant Methods*, vol. 9, no. 1, p. 8, 2013. View at: Publisher Site | Google Scholar - A. Kaestner, M. Schneebeli, and F. Graf, “Visualizing three-dimensional root networks using computed tomography,”
*Geoderma*, vol. 136, no. 1-2, pp. 459–469, 2006. View at: Publisher Site | Google Scholar - S. J. Mooney, T. P. Pridmore, J. Helliwell, and M. J. Bennett, “Developing X-ray computed tomography to non-invasively image 3D root systems architecture in soil,”
*Plant and Soil*, vol. 352, no. 1-2, pp. 1–22, 2012. View at: Publisher Site | Google Scholar - S. R. Tracy, J. A. Roberts, C. R. Black, A. McNeill, R. Davidson, and S. J. Mooney, “The X-factor: visualizing undisturbed root architecture in soils using X-ray computed tomography,”
*Journal of Experimental Botany*, vol. 61, no. 2, pp. 311–313, 2010. View at: Publisher Site | Google Scholar - L. Borisjuk, H. Rolletschek, and T. Neuberger, “Surveying the plant’s world by magnetic resonance imaging,”
*The Plant Journal*, vol. 70, no. 1, pp. 129–146, 2012. View at: Publisher Site | Google Scholar - H. Van As, “Intact plant MRI for the study of cell water relations, membrane permeability, cell-to-cell and long distance water transport,”
*Journal of Experimental Botany*, vol. 58, no. 4, pp. 743–756, 2007. View at: Publisher Site | Google Scholar - H. Gartner, B. Wagner, I. Heinrich, and C. Denier, “3D-laser scanning: a new method to analyze coarse tree root systems,”
*Forest Snow and Landscape Research*, vol. 82, no. 1, pp. 95–106, 2009. View at: Google Scholar - S. Fang, X. Yan, and H. Liao, “3D reconstruction and dynamic modeling of root architecture in situ and its application to crop phosphorus research,”
*The Plant Journal*, vol. 60, no. 6, pp. 1096–1108, 2009. View at: Publisher Site | Google Scholar - R. T. Clark, R. B. MacCurdy, J. K. Jung et al., “Three-dimensional root phenotyping with a novel imaging and software platform,”
*Plant Physiology*, vol. 156, no. 2, pp. 455–465, 2011. View at: Publisher Site | Google Scholar - S. Fang, R. T. Clark, Y. Zheng et al., “Genotypic recognition and spatial responses by rice roots,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 110, no. 7, pp. 2670–2675, 2013. View at: Publisher Site | Google Scholar - P. R. Zurek, C. N. Topp, and P. N. Benfey, “Quantitative trait locus mapping reveals regions of the maize genome controlling root system architecture,”
*Plant Physiology*, vol. 167, no. 4, pp. 1487–1496, 2015. View at: Publisher Site | Google Scholar - M. Michniewicz, C. H. Ho, T. A. Enders et al., “TRANSPORTER OF IBA1 links auxin and cytokinin to Influence root architecture,”
*Developmental Cell*, vol. 50, no. 5, pp. 599–609.e4, 2019. View at: Publisher Site | Google Scholar - W. L. Araújo, T. Tohge, S. Osorio et al., “Antisense inhibition of the 2-oxoglutarate dehydrogenase complex in tomato demonstrates its importance for plant respiration and during leaf senescence and fruit maturation,”
*Plant Cell*, vol. 24, no. 6, pp. 2328–2351, 2012. View at: Publisher Site | Google Scholar - J. Snider, A. Pillai, and C. F. Stevens, “A universal property of axonal and dendritic arbors,”
*Neuron*, vol. 66, no. 1, pp. 45–56, 2010. View at: Publisher Site | Google Scholar - E. Jaynes and G. Bretthorst,
*Probability Theory: The Logic of Science*, Cambridge University Press, 2012. - P. Ferraro, C. Godin, and P. Prusinkiewicz, “Toward a quantification of self-similarity in plants,”
*Fractals*, vol. 13, no. 2, pp. 91–109, 2011. View at: Publisher Site | Google Scholar - C. Godin and P. Ferraro, “Quantifying the degree of self-nestedness of trees: Application to the structural analysis of plants,”
*IEEE/ACM Transactions on Computational Biology and Bioinformatics*, vol. 7, no. 4, pp. 688–703, 2010. View at: Publisher Site | Google Scholar - D. Hunt and V. M. Savage, “Asymmetries arising from the space-filling nature of vascular networks,”
*Physical Review E*, vol. 93, no. 6, article 062305, 2016. View at: Publisher Site | Google Scholar - J. D. Pelletier, “Deviations from self-similarity in barchan form and flux: the case of the Salton Sea dunes, California,”
*Journal of Geophysical Research - Earth Surface*, vol. 118, no. 4, pp. 2406–2420, 2013. View at: Publisher Site | Google Scholar - J. Panico and P. Sterling, “Retinal neurons and vessels are not fractal but space-filling,”
*The Journal of Comparative Neurology*, vol. 361, no. 3, pp. 479–490, 1995. View at: Publisher Site | Google Scholar - B. Milne and V. Gupta, “Horton ratios link self-similarity with maximum entropy of eco-geomorphological properties in stream networks,”
*Entropy*, vol. 19, no. 6, p. 249, 2017. View at: Publisher Site | Google Scholar - Y. S. Liu, C. F. Stevens, and T. O. Sharpee, “Predictable irregularities in retinal receptive fields,”
*Proceedings. National Academy of Sciences. United States of America*, vol. 106, no. 38, pp. 16499–16504, 2009. View at: Publisher Site | Google Scholar - P.-G. De Gennes and P.-G. Gennes,
*Scaling concepts in polymer physics*, Cornell university press, 1979. - H. Toju, P. R. Guimarães Jr., J. M. Olesen, and J. N. Thompson, “Below-ground plant-fungus network topology is not congruent with above-ground plant-animal network topology,”
*Science Advances*, vol. 1, no. 9, article e1500291, 2015. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2020 Sam Sultan et al. Exclusive Licensee Nanjing Agricultural University. Distributed under a Creative Commons Attribution License (CC BY 4.0).