- Research article
- Open Access
Costs and benefits of reticulate leaf venation
BMC Plant Biology volume 14, Article number: 234 (2014)
Recent theoretical and empirical work has identified redundancy as one of the benefits of the reticulate form in the evolution of leaf vein networks. However, we know little about the costs of redundancy or how those costs depend on vein network geometry or topology. Here, we examined both costs and benefits to redundancy in 339 individual reticulate leaf networks comprising over 3.5 million vein segments. We compared levels of costs and benefits within reticulate networks to those within analogous networks without loops known as Maximum Spanning Trees (MSTs).
We show that network robustness to varying degrees of simulated damage is positively correlated with structural indices of redundancy. We further show that leaf vein networks are topologically, geometrically and functionally more redundant than are MSTs. However, increased redundancy comes with minor costs in terms of increases in material allocation or decreases in conductance. We also show that full networks do not markedly decrease the distance to non-vein tissue in comparison to MSTs.
These results suggest the evolutionary transition to the reticulate type of networks found in modern Angiosperm flora involved a relatively minor increase in material and conductance costs with significant benefits in terms of network redundancy.
Hierarchical trees are considered to be the predominant type of physical distribution network in biology . Examples include the ramifying networks found in plants or mammalian cardiovascular or bronchial networks ,. However, not all biological networks are strictly hierarchical, and many networks exhibit both a hierarchical structure and loops that ostensibly allow for redundancy in the face of disturbance or perturbations, where redundancy is defined simply as the existence of multiple flow paths. This is perhaps most evident in the reticulate networks of the leaves of higher plants, notably most angiosperm lineages (Figure 1) −, but reticulate structures are also found in animal lineages such as mammalian capillary beds or some Gorgonian corals.
It has been suggested that the reticulate patterns found in higher leaves allow them to maintain supply of water and nutrients to photosynthetically active chloroplasts even when flow through some channels is lost −, as might be observed due to mechanical damage or herbivory. For example, recent work has shown that two broad classes of venation types, palmate and pinnate leaves, responded differently to network severing treatments . Palmate leaves suffer little loss in leaf hydraulic conductance, stomatal conductance, and photosynthetic rate, when compared to pinnate leaves, indicating that having multiple primary channels enables robustness to herbivory, embolism or other disturbance. Similarly, simulation work has shown that smaller leaves, with a higher vein length per area (VLA) of minor veins, are less vulnerable to embolism that larger leaves . Recent theoretical work has demonstrated that a combination of damage and fluctuating load favors the formation of loops in optimal transport networks ,, and loop formation in a fixed hierarchical tree necessarily increases VLA. The high VLA found in many reticulate angiosperm lineages are associated with high photosynthetic and transpiration rates and are thought to have facilitated the diversification and dominance of angiosperms ,. While certain lineages have retained networks with low, or no reticulation, such as some fern or gymnosperm clades, the overwhelming majority of broad leafed angiosperms have evolved reticulation, with multiple independent origins, suggesting strong selective pressure for this network arrangement ,,. The continued dominance of some communities and ecosystems by taxa such as ferns or gymnosperms, with little or no reticulation, suggests that these network strategies remain viable.
Thus while both empirical data and theoretical results support the idea that network redundancy and associated high VLA are advantageous for leaves, we know very little about the relative costs to leaves of having redundant venation. The costs of redundancy may have multiple origins, e.g., regulatory constraints, hydraulic costs or energetic demands, which need not be mutually exclusive. For example, the regulation of hormones and other factors that give rise to a reticulate vs. a strictly hierarchical network may be more prone to error or may be energetically more costly ,. Similarly, reticulate leaves may have higher resistance than hierarchical networks under certain flow regimes . Finally, the cost of redundancy may be energetic in that: (i) if non-photosynthetic veins (vessel bundles) displace photosynthetic tissue, reticulate networks may suffer from decreased total photosynthesis per unit area; (ii) the amount of materials necessary in the development of a reticulate network may exceed that of an analogous hierarchical network. For example, with estimates of 6.5 and 11.8 mmol glucose per g of cellulose and lignin respectively, xylem tissue has higher carbon costs than surrounding lamina ,. In addition, it has been demonstrated that the primary veins in leaves have lower nitrogen and carbon concentrations and higher density than surrounding lamina .
While there are hydraulic or energetic costs involved in the creation and maintenance of redundant networks, there are also clearly benefits to redundancy, otherwise this network form would not be so prevalent in leaves. The primary benefit to redundancy is the existence of multiple flow paths that maintain flow of water and nutrients to mesophyll tissue under moderate levels of disturbance . There may be additional benefits of redundancy, such as a higher VLA, or a reduction in the distance from veins to stomates and/or chloroplasts, as we discuss.
In this manuscript, we propose a combined empirical and computational approach to quantify the costs of minor vein redundancy in terms of hydraulic and material properties, and benefits in terms of robustness to disturbance and proximity of lamina tissue to veins. The first step in our approach is to extract the spatial structure of individual leaf networks utilizing a recently developed image segmentation and leaf network extraction software (LEAF GUI, www.leafgui.org) ,. Using LEAF GUI, we quantify the vein dimensions and connectivity in 339 leaves from 324 species in 72 angiosperm families, representing semi-automated measurement of 3,934,626 individual vein segments. We measure the level of redundancy of venation networks using established metrics of loopiness , meshedness , and VLA . Next, for each leaf network, we find the maximum spanning tree (MST), that is the network structure that most closely resembles the original leaf network, but that is strictly hierarchical (see Methods). The determination of the MST is equivalent to pruning veins computationally in such a way that the resulting network is both strictly hierarchical and has functional properties (such as estimates of hydraulic conductance) or material properties (such as total volume or vein length) that preserve network hierarchy, and are as close to the original network as possible (Figure 1).
We utilize the MST for inferring the cost of loops in networks based on the premise that bulk flow constraints necessitate a hierarchical tree  and that chloroplasts and/or stomates within leaves cannot be further than a minimum distance from the nearest vein/node ,. For example, open venation systems retain the characteristics of hierarchical trees, and thus any evolutionary transition to reticulate networks likely involved minor vein connections. We then evaluate the evidence that reticulate leaf networks can be pruned to become MSTs with minimal loss of (theoretical) conductance or other material properties. Lastly, to simulate damage to leaves we computationally prune both the reticulate nets and MSTs for each leaf which demonstrates that the overall loss of connectivity is expected to be greater in MSTs. To quantify this cost and evaluate it relative to the other measures we consider we introduce a new metric, "robustness" (see Methods), which represents the capacity of a reticulate network to remain connected to its source/sink node (point of petiole attachment) when subject to random removal of veins, relative to its MST counterpart. We close with a discussion the implications of such results for understanding the evolution and ecology of the reticulate leaf vein form.
We evaluated four metrics related to the reticulate structure and redundancy of leaves: VLA, loopiness, meshedness and robustness. VLA (mm−1) varied from a minimum value of 1.37 to a maximum of 10.76, with a mean of 3.17 (Figure 2). Loopiness (# of areoles/mm2) varied from a minimum value of 0.19 to a maximum of 5.20, with a mean of 1.40 (Figure 3). Meshedness which ranges from 0 (a "tree" structured graph without any loops) to 1 (a network which has a maximum number of loops. Meshedness varied from a minimum value of 0.06 (i.e. more like a tree) to a maximum value of 0.26 (more like a maximally connected planar graph) with a mean of 0.14 standard deviation of 0.04. Robustness, estimated based on the difference between reticulate nets and MSTs weighted for conductance, varied from a minimum of 0.026 to a maximum of 0.150 thus leaves varied in their robustness to disturbance (Figures 2 and 3, Additional file 1: Figure S1–S339). Hence, whereas VLA provides a measure of vein investment per unit area, loopiness is a better indicator of features like distance from non-photosynthetic tissue; meshedness provides a strong indicator of the shape of the network and its tendency to be redundant, and robustness provides a measure of a network's ability to remain connected under perturbations that damage the network. Note, our VLA values are lower on average, but overlap those previously reported ,. This is due to methodological differences, primarily the fact that our images were not magnified and of lower resolution (see discussion in; ,,).
As seen in Figure 2, all four network measures are correlated with one another. Robustness increases with increasing VLA, loopiness and meshedness (Figure 2), with meshedness being the best predictor of robustness (Additional file 2: Table S2). Thus as leaf networks increase their VLA and become more like planar graphs, and less like trees, their ability to buffer perturbations increases.
We evaluated the total cost of redundancy by comparing the total length, width, surface area, volume and conductance on a per-leaf basis (as estimated using the weighted graph extracted via LEAF GUI ) to the same property of the MST. First, as VLA, loopiness, meshedness or robustness increases, the fractional cost of redundancy for length, width and volume measures also increases (Figure 3).
We also find modest increases in network length (6.3%), width (12%), surface area (5.6%), volume (3.3%) or conductance (0.51%) for each reticulate network compared to the corresponding MST, i.e. that minimized for length, width, surface area, volume and conductance, respectively (Figure 4a). In other words, redundancy involves minimal investment in additional transport structures.
The mean distance from non-vein tissue to the venation network were statistically indistinguishable (p > 0.05 in all cases) whether evaluated using the reticulate network or any of its MST counterparts (Figure 4b). Further, we found that the maximum distance from a non-vein component to the network increased by a mean value of 2.73% with most leaves unchanged. Hence, the distances from vein to non-vein regions in MSTs and reticulate networks are not significantly different from one another.
Leaf vein networks display tremendous variety in their form ,, and vein network traits have been shown to be correlated with whole leaf conductance ,, photosynthetic rates , species diversification rates ,, and have been utilized as a proxy for climatic changes . Reticulate veined leaves first appear in the paleo-botanical record in the early Carboniferous  as simple cross linkages between semi-parallel veins. The subsequent divergence in reticulate form is vast and well documented with numerous morphological classes having been identified based primarily on the concept of vein order and arrangement ,. The physiological and theoretical consequences of this transition have only recently been investigated in the context of redundancy −.
We have examined the difference between reticulate nets and MSTs with respect to some of the costs and benefits of minor vein redundancy. However, given the demonstrated links between VLA and photosynthetic rates ,, it may be that there are benefits to the reticulate form over and above simple minor redundancy. As seen in Figure 4a, the mean network increased in length 6.3%. Given a constant area, this corresponds to a mean increase of 6.3% in VLA, which based on previously published empirical relationships , would lead to an increase in photosynthetic rate, all else being equal. Of course photosynthesis and gas exchange are dependent on numerous physiological traits in addition to VLA and a full understanding of the influence of VLA on leaf physiological rates is an important target for future research.
While we observe no statistical difference between MST and reticulate nets in the distance to non-vein areas, MSTs are indeed marginally further away (Figure 4a), thus the decrease in distance over which diffusion is the dominant flow regime that is enabled by reticulate nets may ultimately contribute to increased photosynthetic rate. Thus the reticulate form likely has multiple benefits that have led to the increase in its prominence, with natural selection likely acting not only on increased robustness, but also perhaps photosynthetic rate. Of course it is possible to increase VLA without becoming reticulate, and it is likely that some lineages have taken this course. Moreover, increasing the number of freely ending veinlets in leaves, as is found in many plant families, will also increase VLA without increasing loopiness ,, which may explain why our values for meshedness are not closer to those expected for fully planar networks.
We also find that the cost of redundancy increases with VLA, loopiness, meshedness and robustness (Figure 3). These costs are relatively minor and approach 5% by volume for the loopiest/densest leaves. Vessel bundles in leaves have a higher costs per unit mass due to the fraction of lignin and cellulose in their tissues ,. Thus, given the aforementioned relationship between VLA and measures of leaf performance such as photosynthetic rate, our results suggest that selection on high photosynthetic rates may have the added cost of an increased mass investment in vein structure, over and above that of a strictly hierarchical tree.
Similarly, the costs of redundancy for theoretical conductance are quite low, with a mean of 0.51% of the total. This highlights the fact that theoretical conductance scales approximately with the fourth power of vein radius ,. Large veins contribute much more to total conductance than small veins. Thus an increase in redundancy, by increasing the number and length of the minor most veins, does little to change the overall conductance or resistance (which is proportional to 1/conductance) of the leaf network. Detailed measurements of resistance in both xylem and mesophyll pathways indicate that the partitioning of resistance in leaves has both vein and mesophyll components which vary substantially between species, but are thought to be roughly equivalent on average .
Our estimates of conductance are based on the assumption that xylem conductance follows a Hagen-Poiseuille type scaling with conductance is proportional to radius to the 4th power and length to the 1st power, and further that there exists a consistent proportionality between vessel bundle dimensions and the dimensions of the xylem vessels they contain. While the use of the Hagen-Poiseuille relationship is well established in studies of plant hydraulics ,, due to the difficulty associated with sectioning and imaging small leaf veins, it is not currently known if a constant proportionality exists between vein diameters and xylem diameters. The Laplace-Young law states that for a conduit to resist transmural forces due to capillary tension, its thickness should be directly proportional to its internal radius, suggestive of such a proportionality , in veins that do not provide any additional biomechanical support to the leaf, which is likely true for the minor most veins we consider here. Recent work on tree branches has shown that the ratio of non-conducting to conducting area, remains approximately constant across branches of varying size due to an inverse relationship between xylem size and number, a so called "packing rule" for xylem . It is not known if this relationship holds in leaves, and an understanding of the relationship between vein diameter, the number and size of the xylem contained within veins, and their effect on vein conductance, is an important target for future research.
Overall, our results suggest that the transition from strictly hierarchical trees like those found in early ferns and gymnosperms to the reticulate networks found in subsequent tracheophyte lineages is unlikely to have resulted in a substantial cost either in terms of network resistance, linear dimensions, volume, and presumably mass . Moreover, assuming leaf area is fixed, redundancy has with it the added benefit of increased VLA which is known to increase photosynthetic rates. We suggest that the benefit of increased robustness in the face of disturbance and VLA increase outweighed the rather minor costs of redundancy in terms of material investment. Subsequent analyses will help to reveal how forms of redundancy differ between lineages or habitats, particularly those in which herbivory, high evaporative demand or other disturbances are prevalent.
Our dataset contains observations for leaves representing 339 individual cleared leaf images from 74 eudicot families from the National Cleared Leaf Collection housed at the Museum of Natural History, Smithsonian Institution (Additional file 2: Table S1). We went through the entire collection, selecting images for our analyses based on three criteria: 1) leaves were mostly intact, i.e. free from major tears or other damage; 2) image resolution was sufficient to resolve most of the highest order veins, and; 3) the contrast between leaf veins, areoles and background were significant enough for the LEAF GUI network extraction algorithms to resolve their structure. LEAF GUI is a recently developed software package designed specifically for the analysis of leaf vein images. Extensive descriptions of the underlying algorithms can be found in  and on (www.leafgui.org). Our analysis is based on network connectivity at the whole leaf level. We have used what is, to our knowledge, the most extensive, publically available source of entire images of leaf networks (note the entire collection is available at www.clearedleavesdb.org). We hope that future work on magnified images of entire leaf vein networks (which currently do not exist in large quantities) can confirm these results (see Discussion in ).
The LEAF GUI software returns a characterization of the leaf as a weighted graph comprised of nodes and edges, where an edge is defined as a vein segment, and nodes are defined as the intersection of two or more vein segments. In addition, LEAF GUI extracts metric and positional information for each vein segment, such that each vein segment has an associated vector of weights including length, width, surface area, volume or theoretical conductance, which is proportional to vein segment diameter to the fourth power assuming a proportionality between vein with and conduit width, and also assuming a constant conduit density and viscosity .
The degree of redundancy within each leaf was first estimated using two metrics: loopiness, and meshedness. Briefly, loopiness is defined simply as the number of areoles per unit area . Meshedness is meant to describe whether a network has a tree-like structure (meshedness = 0) or is a complete planar graph (meshedness = 1) regardless of how dense the veins are packed. It is therefore a purely topological index. The definition of meshedness is M = (m-n + 1)/(2n-5) where m is the number of edges (vein segments) and n is the number of nodes ,.
To compare the actual network to its non-reticulate counterpart, we utilized standard optimization routines to determine the MST within the extracted weighted graph for each leaf. To find the MST we employed Prim's algorithm  on the largest connected component and selected the node closest to the point of petiole attachment as the root node. The MST found for each leaf network differed depending on the measure being maximized, i.e. length, width, area, volume or theoretical conductance. The MST is a strictly hierarchical network (i.e., with no loops) which connects all vertices while maximizing some objective function (i.e., the sum of vein segment lengths, widths, surface area, volume or conductance) (Figure 1). Thus, the MST is that which connects all of these nodes without forming loops, thereby preserving vein hierarchy and ensuring supply to chloroplasts without being redundant.
To simulate the effects of network damage that might result from herbivory, embolism, mechanical damage, etc. we introduce an additional measure we term "robustness". To determine robustness, we iteratively pruned each vein network graph removing from one vein segment up to N vein , where N vein is simply the total number of vein segments removed. We performed this iterative operation for each vein network, for both the reticulate nets and their MST counter parts which were based on maximizing hydraulic conductance (Figure, 1, Additional file 1: Figure S1–S339). We then determined the fraction of the total number of nodes that were still connected to the petiole following each pruning iteration, and repeated this basic algorithm 1000 times. Plots of the mean fraction of nodes connected to the root node (the node closest to the point of petiole attachment) vs. the fraction of vein segments removed demonstrate that for a given fraction of vein segments removed, the mean reticulate network has a greater fraction of nodes connected compared to the mean MSTs (Figure 1, Additional file 1: Figure S1–S339). Moreover, the difference between the two curves (i.e.the larger curve minus the smaller curve) represents an additional measure of network redundancy we define as robustness, (Figures 2 and 3). Robustness is given by calculating the difference between the integral for each curve. Rather than trying to fit functions to each curve and then integrating those functions, we estimated the area under each curve through the use of Riemann sums, which is simply the sum of bin width (for example, 1 divided by the number of veins in the MST) times bin height (for example, the number veins remaining connected to the MST divided by the total number of veins in the MST) for all the curves in Additional file 1: Figure S1–S339.
LaBarbera M: Principles of design of fluid transport systems in biology. Science. 1990, 249: 992-1000. 10.1126/science.2396104.
Ball P: Branches: Nature's Patterns: A Tapestry In Three Parts. In Oxford: Oxford University Press; 2009.
Leopold LB: Trees and streams: the efficiency of branching patterns. J Theor Biol. 1971, 31: 339-354. 10.1016/0022-5193(71)90192-5.
Ellis B, Daly DC, Hickey LJ, Johnson KR, Mitchel JD, Wilf P, Wing SL: Manual of Leaf Architecture. morphological description and categorization of dicotyledonous and net-veined monocotyledonous angiosperms by Leaf Architecture. Cornell University Press, Washington D.C; 2009.
Brodribb TJ, Feild TS, Sack L: Viewing leaf structure and evolution from a hydraulic perspective. Funct Plant Biol. 2010, 37 (6): 488-498. 10.1071/FP10010.
Kull U, Herbig A: Leaf venation of angiosperms - form and evolution. Naturwissenschaften. 1995, 82 (10): 441-451. 10.1007/BF01131595.
Price CA, Enquist BJ: Scaling mass and morphology in leaves: an extension of the WBE model. Ecology. 2007, 88: 1132-1141. 10.1890/06-1158.
Couder Y, Pauchard L, Allain C, Adda-Bedia M, Douady S: The leaf venation as formed in a tensorial field. Eur Phys J B. 2002, 28 (2): 135-138. 10.1140/epjb/e2002-00211-1.
Nardini A, Tyree MT, Salleo S: Xylem cavitation in the leaf of Prunus laurocerasus and its impact on leaf hydraulics. Plant Physiol. 2001, 125 (4): 1700-1709. 10.1104/pp.125.4.1700.
Katifori E, Szollosi GJ, Magnasco MO: Damage and fluctuations induce loops in optimal transport networks. Phys Rev Lett. 2010, 104 (4): 4-10.1103/PhysRevLett.104.048704.
Sack L, Dietrich EM, Streeter CM, Sanchez-Gomez D, Holbrook NM: Leaf palmate venation and vascular redundancy confer tolerance of hydraulic disruption. Proc Natl Acad Sci U S A. 2008, 105 (5): 1567-1572. 10.1073/pnas.0709333105.
Corson F: Fluctuations and redundancy in optimal transport networks. Phys Rev Lett. 2010, 104 (4): 4-10.1103/PhysRevLett.104.048703.
Scoffoni C, Rawls M, McKown A, Cochard H, Sack L: Decline of leaf hydraulic conductance with dehydration: relationship to leaf size and venation architecture. Plant Physiol. 2011, 156 (2): 832-843. 10.1104/pp.111.173856.
Boyce CK, Brodribb TJ, Feild TS, Zwieniecki MA: Angiosperm leaf vein evolution was physiologically and environmentally transformative. Proc R Soc B Biol Sci. 2009, 276 (1663): 1771-1776. 10.1098/rspb.2008.1919.
Brodribb TJ, Feild TS: Leaf hydraulic evolution led a surge in leaf photosynthetic capacity during early angiosperm diversification. Ecol Lett. 2010, 13 (2): 175-183. 10.1111/j.1461-0248.2009.01410.x.
Roth-Nebelsick A, Uhl D, Mosbrugger V, Kerp H: Evolution and function of leaf venation architecture. Ann Bot. 2001, 87: 553-556. 10.1006/anbo.2001.1391.
Berleth T, Scarpella E, Prusinkiewicz P: Towards the systems biology of auxin-transport-mediated patterning. Trends Plant Sci. 2007, 12 (4): 151-159. 10.1016/j.tplants.2007.03.005.
Scarpella E, Marcos D, Friml J, Berleth T: Control of leaf vascular patterning by polar auxin transport. Genes Dev. 2006, 20 (8): 1015-1027. 10.1101/gad.1402406.
McKown AD, Cochard H, Sack L: Decoding leaf hydraulics with a spatially explicit model: principles of venation architecture and implications for its evolution. Am Nat. 2010, 175 (4): 447-460. 10.1086/650721.
Lambers H, Poorter H: Inherent variation in growth-rate between higher-plants - a search for physiological causes and ecological consequences. Adv Ecol Res. 1992, 23: 187-261. 10.1016/S0065-2504(08)60148-8.
Niinemets U, Portsmuth A, Tena D, Tobias M, Matesanz S, Valladares F: Do we underestimate the importance of leaf size in plant economics? Disproportional scaling of support costs within the spectrum of leaf physiognomy. Ann Bot. 2007, 100 (2): 283-303. 10.1093/aob/mcm107.
Price CA, Symonova O, Mileyko Y, Hilley T, Weitz JS: Leaf extraction and analysis framework graphical user interface: segmenting and analyzing the structure of leaf veins and areoles. Plant Physiol. 2010, 155: 236-245. 10.1104/pp.110.162834.
Price CA: LEAF GUI: Analyzing The Geometry Of Veins And Areoles Using Image Segmentation Algorithms. High Throughput Phenotyping In Plants. Edited by: Normanly J. Humana, New York; 2012.
Blonder B, Violle C, Bentley LP, Enquist BJ: Venation networks and the origin of the leaf economics spectrum. Ecol Lett. 2010, 14 (2): 91-100. 10.1111/j.1461-0248.2010.01554.x.
Buhl J, Gautrais J, Reeves N, Sole RV, Valverde S, Kuntz P, Theraulaz G: Topological patterns in street networks of self-organized urban settlements. Eur Phys J B. 2006, 49 (4): 513-522. 10.1140/epjb/e2006-00085-1.
Uhl D, Mosbrugger V: Leaf venation density as a climate and environmental proxy: a critical review and new data. Palaeogeogr Palaeoclimatol Palaeoecol. 1999, 149 (1–4): 15-26. 10.1016/S0031-0182(98)00189-8.
Raven JA: The significance of the distance from photosynthesizing cells to vascular tissue in extant and early vascular plants. Bot J Scotl. 1994, 47 (1): 65-81. 10.1080/03746609408684818.
Price CA, Wing SL, Weitz JS: Scaling and structure of dicotyledonous leaf venation networks. Ecol Lett. 2011, 15 (2): 87-95. 10.1111/j.1461-0248.2011.01712.x.
Sack L, Frole K: Leaf structural diversity is related to hydraulic capacity in tropical rain forest trees. Ecology. 2006, 87 (2): 483-491. 10.1890/05-0710.
Brodribb TJ, Feild TS, Jordan GJ: Leaf maximum photosynthetic rate and venation are linked by hydraulics. Plant Physiol. 2007, 144 (4): 1890-1898. 10.1104/pp.107.101352.
Price CA, Munro PR, Weitz JS: Estimates of leaf vein density are scale dependent. Plant Physiol In Press.
Sack L, Scoffoni C, McKown AD, Frole K, Rawls M, Havran JC, Tran H, Tran T: Developmentally based scaling of leaf venation architecture explains global ecological patterns. Nat Commun. 2012, 3: 1-10. 10.1038/ncomms1835.
Dengler N, Kang J: Vascular patterning and leaf shape. Curr Opin Plant Biol. 2001, 4 (1): 50-56. 10.1016/S1369-5266(00)00135-7.
Cochard H, Nardini A, Coll L: Hydraulic architecture of leaf blades: where is the main resistance?. Plant Cell Environ. 2004, 27 (10): 1257-1267. 10.1111/j.1365-3040.2004.01233.x.
Sack L, Holbrook NM: Leaf Hydraulics. Annu Rev Plant Biol. 2006, 57: 361-381. 10.1146/annurev.arplant.56.032604.144141.
Hickey LJ: Classification of the architecture of dicotyledonous leaves. Am J Bot. 1973, 60 (1): 17-33. 10.2307/2441319.
McKown AD, Dengler NG: Shifts in leaf vein density through accelerated vein formation in C-4 Flaveria (Asteraceae). Ann Bot. 2009, 104 (6): 1085-1098. 10.1093/aob/mcp210.
Lewis AM, Boose ER: Estimating volume flow-rates through xylem conduits. Am J Bot. 1995, 82 (9): 1112-1116. 10.2307/2446063.
Tyree MT, Ewers FW: The hydraulic architecture of trees and other woody plants. New Phytol. 1991, 119: 345-360. 10.1111/j.1469-8137.1991.tb00035.x.
Sack L, Scoffoni C: Leaf venation: structure, function, development, evolution, ecology and applications in the past, present and future. New Phytol. 2013, 198 (4): 983-1000. 10.1111/nph.12253.
Sherman TF: On connecting large vessels to small - the meaning of Murray Law. J Gen Physiol. 1981, 78 (4): 431-453. 10.1085/jgp.78.4.431.
Price CA, Knox S-JC, Brodribb TJ: The influence of branch order on optimal leaf vein geometries: Murray's law and area preserving branching. PLoS One, In Press.
Savage VM, Bentley LP, Enquist BJ, Sperry JS, Smith DD, Reich PB, von Allmen EI: Hydraulic trade-offs and space filling enable better predictions of vascular structure and function in plants. Proc Natl Acad Sci U S A. 2010, 107 (52): 22722-22727. 10.1073/pnas.1012194108.
Tyree MT, Ewers FW: Hydraulic Architecture Of Woody Tropical Plants. Tropical Forest Plant Ecophysiology. Edited by: Mulkey SS, Chazdon RL, Smith AP. 1996, Chapman and Hall, New York, 217-243. 10.1007/978-1-4613-1163-8_8.
Barthelemy M: Spatial Networks. Phys Rep. 2010, 499 (2011): 1-101.
Prim RC: Shortest connection networks and some generalizations. Bell Syst Tech J. 1957, 36 (6): 1389-1401. 10.1002/j.1538-7305.1957.tb01515.x.
CAP acknowledges the support of a Discovery Early Career Research Award (DECRA) from the Australian Research Council. JSW holds a Career Award at the Scientific Interface from the Burroughs Wellcome Fund. Scott Wing provided valuable assistance with the Smithsonian leaf image collection.
We declare that there are no competing interests with respect to this manuscript.
CAP and JSW conceived the study, performed the analysis and wrote, read, and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1-S339.: Plots of the leaf level mean fraction of the network disconnected from the petiole vs. the fraction of the vein segments removed (see Methods) for both full reticulate networks (blue symbols) and MSTs (red symbols) for all 339 leaves. (PDF 9 MB)
About this article
Cite this article
Price, C.A., Weitz, J.S. Costs and benefits of reticulate leaf venation. BMC Plant Biol 14, 234 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-014-0234-2
- Leaf veins
- Reticulate veins
- Network robustness