How dense can you be? New automatic measures of vein density in angiosperm leaves

Premise

Because of the trade-off between water loss and carbon dioxide assimilation, the conductivity of the transpiration path in a leaf is an important limit on photosynthesis. Closely packed veins correspond to short paths and high assimilation rates while widely spaced veins are associated with higher resistance to flow and lower maximum photosynthetic rates. Vein length per area (VLA) has become the standard metric for comparing leaves with different vein densities; its measurement typically utilizes digital image processing with varying amounts of human input.

Methods and Results

Here, we propose three new ways of measuring vein density using image analysis that improve on currently available procedures: (1) areole area distributions, (2) a sizing transform, and (3) a distance map. Each alternative has distinct practical, statistical, and biological limitations and advantages. In particular, we advocate the log-transformed modal distance map of a vein mask as an estimator to replace VLA as a standard metric for vein density.

Conclusions

These methods, for which open-source code appropriate for high-throughput automation is provided, improve on VLA by producing determinate measures of vein density as distributions rather than point estimates. Combined with advances in image quality and computational efficiency, these methods should help clarify the physiological and evolutionary significance of vein density.

CO2 starvation experiments provide support for the carbon-limited hypothesis on the evolution of CAM-like behaviour in Isoëtes

BACKGROUND AND AIMS

Crassulacean acid metabolism (CAM) is an adaptation to increase water use efficiency in dry environments. Similar biochemical patterns occur in the aquatic lycophyte genus Isoëtes. It has long been assumed and accepted that CAM-like behaviour in these aquatic plants is an adaptation to low daytime carbon levels in aquatic ecosystems, but this has never been directly tested.

METHODS

To test this hypothesis, populations of Isoëtes engelmannii and I. tuckermanii were grown in climate-controlled chambers and starved of atmospheric CO2 during the day while pH was measured for 24 h.

KEY RESULTS

We demonstrate that terrestrial plants exposed to low atmospheric CO2 display diel acidity cycles similar to those in both xerophytic CAM plants and submerged Isoëtes.

CONCLUSIONS

Daytime CO2 starvation induces CAM-like nocturnal acid accumulation in terrestrial Isoëtes, substantiating the hypothesis that carbon starvation is a selective pressure for this physiological behaviour.

Global plant trait relationships extend to the climatic extremes of the tundra biome

The majority of variation in six traits critical to the growth, survival and reproduction of plant species is thought to be organised along just two dimensions, corresponding to strategies of plant size and resource acquisition. However, it is unknown whether global plant trait relationships extend to climatic extremes, and if these interspecific relationships are confounded by trait variation within species. We test whether trait relationships extend to the cold extremes of life on Earth using the largest database of tundra plant traits yet compiled. We show that tundra plants demonstrate remarkably similar resource economic traits, but not size traits, compared to global distributions, and exhibit the same two dimensions of trait variation. Three quarters of trait variation occurs among species, mirroring global estimates of interspecific trait variation. Plant trait relationships are thus generalizable to the edge of global trait-space, informing prediction of plant community change in a warming world.

TRY plant trait database - enhanced coverage and open access

Plant traits-the morphological, anatomical, physiological, biochemical and phenological characteristics of plants-determine how plants respond to environmental factors, affect other trophic levels, and influence ecosystem properties and their benefits and detriments to people. Plant trait data thus represent the basis for a vast area of research spanning from evolutionary biology, community and functional ecology, to biodiversity conservation, ecosystem and landscape management, restoration, biogeography and earth system modelling. Since its foundation in 2007, the TRY database of plant traits has grown continuously. It now provides unprecedented data coverage under an open access data policy and is the main plant trait database used by the research community worldwide. Increasingly, the TRY database also supports new frontiers of trait-based plant research, including the identification of data gaps and the subsequent mobilization or measurement of new data. To support this development, in this article we evaluate the extent of the trait data compiled in TRY and analyse emerging patterns of data coverage and representativeness. Best species coverage is achieved for categorical traits-almost complete coverage for 'plant growth form'. However, most traits relevant for ecology and vegetation modelling are characterized by continuous intraspecific variation and trait-environmental relationships. These traits have to be measured on individual plants in their respective environment. Despite unprecedented data coverage, we observe a humbling lack of completeness and representativeness of these continuous traits in many aspects. We, therefore, conclude that reducing data gaps and biases in the TRY database remains a key challenge and requires a coordinated approach to data mobilization and trait measurements. This can only be achieved in collaboration with other initiatives.

Reading the leaves: A comparison of leaf rank and automated areole measurement for quantifying aspects of leaf venation

The reticulate venation that is characteristic of a dicot leaf has excited interest from systematists for more than a century, and from physiological and developmental botanists for decades. The tools of digital image acquisition and computer image analysis, however, are only now approaching the sophistication needed to quantify aspects of the venation network found in real leaves quickly, easily, accurately, and reliably enough to produce biologically meaningful data. In this paper, we examine 120 leaves distributed across vascular plants (representing 118 genera and 80 families) using two approaches: a semiquantitative scoring system called "leaf ranking," devised by the late Leo Hickey, and an automated image-analysis protocol. In the process of comparing these approaches, we review some methodological issues that arise in trying to quantify a vein network, and discuss the strengths and weaknesses of automatic data collection and human pattern recognition. We conclude that subjective leaf rank provides a relatively consistent, semiquantitative measure of areole size among other variables; that modal areole size is generally consistent across large sections of a leaf lamina; and that both approaches-semiquantitative, subjective scoring; and fully quantitative, automated measurement-have appropriate places in the study of leaf venation.

Reading the leaves: A comparison of leaf rank and automated areole measurement for quantifying aspects of leaf venation

The reticulate venation that is characteristic of a dicot leaf has excited interest from systematists for more than a century, and from physiological and developmental botanists for decades. The tools of digital image acquisition and computer image analysis, however, are only now approaching the sophistication needed to quantify aspects of the venation network found in real leaves quickly, easily, accurately, and reliably enough to produce biologically meaningful data. In this paper, we examine 120 leaves distributed across vascular plants (representing 118 genera and 80 families) using two approaches: a semiquantitative scoring system called "leaf ranking," devised by the late Leo Hickey, and an automated image-analysis protocol. In the process of comparing these approaches, we review some methodological issues that arise in trying to quantify a vein network, and discuss the strengths and weaknesses of automatic data collection and human pattern recognition. We conclude that subjective leaf rank provides a relatively consistent, semiquantitative measure of areole size among other variables; that modal areole size is generally consistent across large sections of a leaf lamina; and that both approaches-semiquantitative, subjective scoring; and fully quantitative, automated measurement-have appropriate places in the study of leaf venation.

The function of the aerenchyma in arborescent lycopsids: evidence of an unfamiliar metabolic strategy

Most species of the modern family Isoëtaceae (Quillworts) some other modern hydrophytes, use a metabolic pathway for carbon fixation that involves uptake of sedimentary carbon and enrichment of CO(2) in internal gas spaces as a carbon-concentrating mechanism. This metabolism, which is related to 'aquatic CAM', is characterized by morphological, physiological and biochemical adaptations for decreasing photorespirative loss, aerating roots and maintaining high growth rates in anoxic, oligotrophic, stressed environments. Some of the closest relatives of the Isoëtaceae were the 'arborescent lycopsids', which were among the dominant taxa in the coal swamps found in lowland ecosystems during the Carboniferous and Permian periods (approx. 300 Ma). Morphological, ecological and geochemical evidence supports the hypothesis that the arborescent lycopsids had an unusual metabolism similar to that of modern Isoëtaceae and processed a biogeochemically significant proportion of organically fixed carbon over a period of about 100 million years in the late Palaeozoic. The temporal coincidence between the dominance of plants with this metabolism and an anomalous global atmosphere (high O(2); low CO(2)) supports the idea that biosphere feedbacks are important in regulating global climatic homeostasis. The potential influence of this metabolism on the global carbon cycle and its specific adaptive function suggest that it should perhaps be considered a fourth major photosynthetic pathway.

Transient response of quasi-isotropic fiber composite laminates to internal line sources

This article considers the propagation of elastic waves in an eight-ply quasi-isotropic laminate arising from line sources of dislocation located at each of the seven interfaces in turn. The laminate is composed of identical layers of a fiber composite material which is modeled as a homogeneous transversely isotropic elastic continuum with the axis of transverse isotropy along the fiber direction. The line source sets up a straight crested wave traveling along the laminate in the direction normal to the load line and the elastodynamic equations within each layer are solved by taking the Laplace transform with respect to time and the Fourier transform with respect to the spatial coordinate in the direction of propagation. The resulting system of six first-order differential equations in each layer is solved to obtain the transforms of the displacement and stress components throughout the laminate. The time history of any displacement or stress component at any location may then be recovered by numerical inversion of the double transform. The graphs presented show the time history of the normal displacement of the top surface of the laminate at distances of 1 and 20 plate thicknesses from the plane of action of the sources. These graphs are for four different orientations of the line of action of the sources, namely, at angles 0, 30, 60, and 90 degrees to the fiber direction in the surface layer. The numerical inversion involves a summation over different modes of Rayleigh-Lamb waves in the laminate and results are also presented showing the contributions to the overall response from some of the individual modes.

A computational method for wave propagation from a point load in an anisotropic material

One approach which is employed to solve dynamic point load problems in plates and laminates is to take integral transforms to reduce the governing equations to a system of ordinary differential equations with respect to the depth variable. The solution of this system leads to expressions for the transforms of the displacement and stress components at any level in the plate and the transient response at any location may then be recovered by inversion of the multiple transforms. The formal transform inversion involves a double infinite integral but by making a change of variable this may be replaced by an infinite integral associated with a line source and a finite integral with respect to the orientation of the line. A first attempt at applying this approach to obtain the point load response of quasi-isotropic fibre composite laminate led to a non-causal predicted signal. This paper deals with an investigation of this proposed method applied to the simpler model problem of wave propagation in a two-dimensional anisotropic medium. Results are obtained for two different time histories of point loads, namely: a delta function; and a single period of a sine function. In the case of the delta function source a comparison is made with the analytic solution and the errors arising from the numerical approach are discussed. Graphs are also presented showing the non-causal contributions to the overall response which arise at individual angles of orientation of the line source.

Elastic waves in fibre composite laminates

This paper considers the propagation of elastic waves in an eight-ply quasi-isotropic laminate arising from line sources of dislocation located at each of the seven interfaces in turn. The line source sets up a straight crested wave travelling along the laminate in a direction normal to the load line, and the elastodynamic equations within each layer are solved by taking the Laplace transform with respect to time and the Fourier transform with respect to the spatial coordinate in the direction of propagation. The resulting system of six first-order differential equations in each layer is solved to obtain the transforms of the displacement and stress components throughout the laminate. The time history of any displacement or stress component at any location may then be recovered by numerical inversion of the double transform. Examples are shown of the time history of the normal displacement of the top surface of the laminate at a distance of 20 plate thicknesses from the plane of action of the sources. The numerical inversion involves a summation over different modes of Rayleigh-Lamb waves in the laminate and contributions to the overall response from some of the individual modes are displayed.