Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 31 March 2015
Sec. Functional Plant Ecology

Patterns in hydraulic architecture from roots to branches in six tropical tree species from cacao agroforestry and their relation to wood density and stem growth

  • 1Plant Ecology and Ecosystems Research, Albrecht von Haller Institute for Plant Sciences, University of Göttingen, Göttingen, Germany
  • 2Faculty of Agriculture, Tadulaku University, Palu, Indonesia

For decades it has been assumed that the largest vessels are generally found in roots and that vessel size and corresponding sapwood area-specific hydraulic conductivity are acropetally decreasing toward the distal twigs. However, recent studies from the perhumid tropics revealed a hump-shaped vessel size distribution. Worldwide tropical perhumid forests are extensively replaced by agroforestry systems often using introduced species of various biogeographical and climatic origins. Nonetheless, it is unknown so far what kind of hydraulic architectural patterns are developed in those agroforestry tree species and which impact this exerts regarding important tree functional traits, such as stem growth, hydraulic efficiency and wood density (WD). We investigated wood anatomical and hydraulic properties of the root, stem and branch wood in Theobroma cacao and five common shade tree species in agroforestry systems on Sulawesi (Indonesia); three of these were strictly perhumid tree species, and the other three tree species are tolerating seasonal drought. The overall goal of our study was to relate these properties to stem growth and other tree functional traits such as foliar nitrogen content and sapwood to leaf area ratio. Our results confirmed a hump-shaped vessel size distribution in nearly all species. Drought-adapted species showed divergent patterns of hydraulic conductivity, vessel density, and relative vessel lumen area between root, stem and branch wood compared to wet forest species. Confirming findings from natural old-growth forests in the same region, WD showed no relationship to specific conductivity. Overall, aboveground growth performance was better predicted by specific hydraulic conductivity than by foliar traits and WD. Our study results suggest that future research on conceptual trade-offs of tree hydraulic architecture should consider biogeographical patterns underlining the importance of anatomical adaptation mechanisms to environment.

Introduction

The water transport pattern in trees is mainly determined by the plant hydraulic architecture, i.e., the spatial distribution of various xylem properties from roots to branches of a tree individual (McCulloh et al., 2010). The hydraulic efficiency of different compartments along the root-to-leaf flow path can be described by the sapwood area-specific hydraulic conductivity (KS), which is directly related to the hydraulic resistance of a given position (Tyree and Ewers, 1991; McElrone et al., 2004). According to the Hagen–Poiseuille law, even a small increase in mean vessel diameter causes an exponential increase of KS. This anatomical pattern represents the most economical way for a woody plant to enhance the path-length hydraulic conductivity. Independently of the efficiency of the hydraulic system water is transported in a metastable state below its vapor pressure in vascular plants, which makes them vulnerable to the formation of gas embolism. This can impair the transport of water from the soil to the leaves. Particularly wide vessels are not only most conductive but also most prone to the risk of hydraulic failure in form of xylem embolisms (Maherali et al., 2006; Awad et al., 2010; Cai et al., 2010; Hajek et al., 2014), resulting in a trade-off between hydraulic efficiency and cavitation resistance. As one of the basic organizing principles of tree hydraulic architecture it has been postulated that the mean vessel diameter in the xylem tissue generally decreases acropetally from roots to branches (‘vessel tapering’: Baas, 1982; Tyree and Zimmermann, 2002; Anfodillo et al., 2013). This principle has stimulated several conceptual models on plant hydraulic architecture during the past 15 years. They state that whole-plant hydraulic conductance dependent on distance to ground in support of the model by West et al. (1999) and Murray’s law (McCulloh et al., 2003). Consistent with these predictions it has indeed commonly been observed that the largest vessels along the water flow path are found in roots of trees from temperate or Mediterranean environments (Martinez-Vilalta et al., 2002; Pratt et al., 2007; Domec et al., 2009). However, recent studies in tropical forests in South America (Machado et al., 2007; Fortunel et al., 2013) and Indonesia (Schuldt et al., 2013) have produced contradictory results regarding the paradigm of continuous vessel tapering. Schuldt et al. (2013) supposed that mechanisms reducing cavitation risk may not have been evolved in these moist or perhumid environments where drought stress is normally not apparent.

Forested perhumid regions particularly in the tropics are underrepresented in studies so far and are moreover converted rapidly. Worldwide approximately 27.2 million ha of humid tropical forests have been cleared between 2000 and 2005 (Hansen et al., 2008) mainly for agricultural land use (Achard et al., 2002; FAO and JRC, 2012). In South-East Asia, a common driver of deforestation is the conversion of natural forests into cacao (Theobroma cacao) agroforestry systems. Cacao is native to tropical South America (Motamayor et al., 2008) and represents one of the commercially most important perennial cash crops worldwide. Traditionally cacao trees are planted under selectively thinned primary or older secondary forest in Indonesia, but nowadays cultivation is shifting to non-shaded monocultures or agroforests with introduced fast-growing legume tree species such as Gliricidia sepium to increase short-term income (Rice and Greenberg, 2000). Shade trees in cacao plantations enhance functional biodiversity, carbon sequestration, soil fertility and drought resistance and provide microclimatic benefits such as increased humidity and buffering temperature extremes (Schroth and Harvey, 2007; Tscharntke et al., 2011).

Considering the ecological relevance of the anatomical hydraulic properties described above, it is important to note that systematic studies on the ecological wood anatomy and hydraulic architecture of cacao and shade tree species are lacking so far. This is all the more unsatisfying since tropical agroforestry crop and shade tree species often originate from different biomes and possess distinct drought adaptations, but it is unknown so far if this implies differences in the hydraulic strategy of those crop and shade tree species. It is therefore unknown whether cacao and shade tree species in the agroforestry systems with different biogeographical origin have developed similar hydraulic properties as the tree species of the natural forest replaced by those.

A high aboveground biomass (AGB) production (including high crop yield) has been related to several plant functional traits like high stem hydraulic efficiency, high foliar nitrogen content, or low stem wood density (WD; Brodribb et al., 2002; Tyree, 2003; Zhang and Cao, 2009; Hoeber et al., 2014). Thereby low WD implying lower hydraulic safety is found to be associated with fast tree growth (Enquist et al., 1999; King et al., 2005; Poorter et al., 2010), while species with dense wood are considered to be more resistant to xylem cavitation due to the commonly assumed relation between WD and conduit size and thus xylem wall thickness and resistance to cell wall implosion under negative pressure (Jacobsen et al., 2005). Consequently, species with dense wood should show higher hydraulic safety at the cost of lower productivity (Meinzer et al., 2003; Bucci et al., 2004). Nevertheless, several studies, particularly from tropical environments, found WD decoupled of hydraulic efficiency traits and growth performance (Zhang and Cao, 2009; Russo et al., 2010; Fan et al., 2012; Schuldt et al., 2013). It would therefore be interesting to assess whether hydraulic properties and WD are related to the aboveground performance of crop and shade tree species in cacao agroforests.

In this study, we examined the inter-relationship between sapwood area-specific hydraulic conductivity of the root, stem, and branch xylem tissue with wood anatomical traits along the water flow path across six common cacao agroforestry tree species with different biogeographical origins from either seasonally dry or perhumid tropical environments growing in cacao agroforests in Central Sulawesi (Indonesia). We moreover wanted to relate aboveground growth performance to hydraulic efficiency, stem WD, foliar nitrogen content and foliar δ13C of these species. We hypothesized (i) that – in contrast to temperate tree species – the largest vessels along the water flow path are found in the stem xylem and not in the roots, (ii) that stem xylem hydraulic properties are unrelated to stem WD, and (iii) that aboveground productivity across species is positively related to vessel size and hydraulic conductivity.

Materials and Methods

Study site, Species, and Sampling

The study was carried out in a cacao agroforestry located in the Kulawi Valley, Bolabapu District, Central Sulawesi, Indonesia (S 01°55.9 E 120°02.2, elevation 571 m above see level) in May 2012. The climate of the study region is perhumid without a distinct dry season. Mean annual temperature recorded for the study area by Moser et al. (2010) is 25.5°C and mean annual precipitation is 2092 mm between 2002 and 2006. For the study, a cacao agroforestry plot with multi-species shade tree layer was selected from a larger number of preselected cacao agroforestry plots of a different investigation that were found to be representative in terms of management, aboveground structure and topographical patterns for this region. Caution was taken during the selection process that the plot was far enough above the groundwater table to guarantee that the trees had no direct access to this water source. All trees in the agroforest were planted simultaneously around 25 years ago.

Theobroma cacao L. (Malvaceae) originating from rainforests of lowland northern South America and five common shade tree species were studied: Leucaena leucocephala (Lam.) de Wit and Gliricidia sepium (Jacq.) Steud. (both Fabaceae), which are introduced species from seasonal dry forest areas of Central America. The three other species represent native origins: Gnetum gnemon L. (Gnetaceae), the short-term drought-tolerating Erythrina subumbrans (Hassk.) Merill (Fabaceae) and the strictly perhumid species Durio zibethinus. Murr. (Malvaceae). In the following we have grouped the species according to their drought tolerance as perhumid (T. cacao, D. zibethinus, G. gnemon) and seasonal (G. sepium, L. leucocephala, E. subumbrans). All species have diffuse-porous wood with G. gnemon being a gymnosperm bearing vessels structurally similar to angiosperms (Carlquist, 1994; Fisher and Ewers, 1995). We chose six tree replicates of each species with a diameter and height representative for the whole agroforestry (Table 2). For each tree three sun-exposed upper-crown branches and three topsoil root segments (diameter 6–14 mm; length 25–35 cm) were collected as well as one stem core of 5 cm length per tree taken with an increment corer (Haglöf, Långsele, Sweden) at 130 cm stem height. To ensure species identity the roots were traced back to the tree stem. In order to avoid microbial growth in the extracted tree organs, samples were stored in polyethylene tubes filled with water containing a sodium–silver chloride complex (Micropur Katadyn, Wallisellen, Switzerland). The samples were kept cool at 4°C and the conductivity measurements took place not more than 7 days after collection.

TABLE 1
www.frontiersin.org

TABLE 1. List of major variables with definition and units employed.

TABLE 2
www.frontiersin.org

TABLE 2. Tree height (H), diameter at breast height (DBH), wood density (WD), aboveground biomass (AGB), and basal area increment (BAI) of the six tree species in cocoa agroforests.

WD, AGB, and Productivity

Wood density, defined as oven-dry weight over wet volume, was measured for each stem core. The fresh volume of each sample was determined by Archimedes’ principle. Samples were then oven dried for 48 h at 105°C and dry mass recorded.

Aboveground biomass of the trees was calculated using the allometric equation of Chave et al. (2005) for tropical wet stands as: AGB = exp [-2.187 + 0.916 × ln (WD × DBH2 × H)], where AGB is the estimated aboveground biomass (kg), DBH the trunk diameter at 130 cm height (cm), H the total tree height (m), and WD the stem wood density (g cm-3). Since we obtained proper data on tree height only at the beginning of the study, we used stem basal area increment (BAI, cm2 yr-1) determined over a period of 12 months using dendrometer tapes (UMS GmbH, München, Germany) as indicator for aboveground productivity. However, it has been shown that AGB increment and BAI are very closely related in tropical trees (Hoeber et al., 2014). For T. cacao, G. sepium and D. zibethinus six tree replicates were monitored, whereas data from just three G. gnemon, two L. leucocephala and one of E. subumbrans were available for BAI.

Leaf Morphological and Chemical Properties

From each branch segment harvested for the hydraulic and anatomical measurements, all distal leaves were stripped off and oven-dried at 70°C for 48h to determine leaf dry weight. Specific leaf area (SLA, cm2 g-1) values were determined using data from nine additional branches per species where leaf surface areas were measured with the WinFOLIA software (Régent Instruments, Quebec, QC, Canada). Total leaf area per branch segment (AL, m2) was calculated by dividing dry weight through species-specific SLA values. Subsequently, leaf samples were grounded and analyzed for their foliar concentrations of C and N and for their foliar signatures of δ13C in the leaf bulk tissue with a Delta plus isotope mass spectrometer (Finnigan MAT, Bremen, Germany), a Conflo III interface (Thermo Electron Coorperation, Bremen, Germany) and a NA2500 elemental analyzer (CE-Instruments, Rodano, Milano, Italy) using standard δ notion: δ = (Rsample/Rstandard - 1) × 1000 (‰) in the laboratory for stable isotope measurements (KOSI) at the University of Göttingen.

Empirical Conductivity Measurements

Hydraulic conductivity of one to three root and branch segments per tree was empirically measured using the method described by Sperry et al. (1988). In total, 44 root and 39 branch segments were analyzed (mean root segment length ± SE: 291 ± 7.0 mm and diameter: 7.87 ± 0.25 mm; mean branch segment length: 308 ± 4.3 mm, and diameter: 9.12 ± 0.29 mm). All segments were recut under water with a razor blade, small lateral roots, and branches cut-off and sealed with quick-drying superglue (Loctite 431, Henkel, Düsseldorf, Germany) and activator (Loctite 7452 Aktivator, Henkel, Düsseldorf, Germany) that function on wet materials. Afterward, segments were attached under water to the tubing system of the conductivity apparatus, where the pressure difference of 6 kPa was generated by a 60 cm high water column. De-ionized water with a sodium–silver chloride complex (16 μg L-1 Ag, 8 mg L-1 NaCl, Micropur katadyn, Wallisellen, Switzerland) was used as measuring solution in order to avoid microbial growth in the tubing system, a common problem in tropical environments. While comparing our data with conductivities determined by other solutions, it has to be considered that different perfusion solutions can affect hydraulic conductivity (Espino and Schenk, 2011). The solution was passed through a 0.2 μm membrane filter (Maxi Capsule, Pall Corporation, USA) and each sample measured three times in row and flushed with the measuring solution for 5 min at 120 kPa in between each measurement to remove potential emboli. The hydraulic conductivity (Khemp, kg m s-1MPa-1) was calculated as Kh = (ΔV/Δt) × (l/ΔP), where l is the length of the segment (m), ΔP the pressure difference applied to the segment (MPa), ΔV the amount of water flowing out of the segment (kg), and Δt the time interval of measurement (s).

Segments of the branches and roots used for conductivity measurements were planed with a sliding microtome (G.S.L.1, WSL, Birmensdorf, Switzerland) to obtain high-quality top view images with a stereo-microscope (SteREOV20, Carl Zeiss MicroImaging GmbH, Göttingen, Gemany) and total cross-sectional (Across, mm2) and xylem cross-sectional area (Axylem, mm2) analyzed with ImageJ (v1.44p). Subsequently, for each species a regression analysis between Across and Axylem was carried out (Table A1). Empirical sapwood area-specific hydraulic conductivity (KSemp, kg m-1 MPa-1 s-1) was calculated by dividing Khemp by the calculated mean xylem cross-sectional area without pith and bark by applying the species-specific regression coefficients, and empirical leaf area-specific hydraulic conductivity (KLemp, kg m-1 MPa-1 s-1) by dividing Kh by the total supported leaf area (AL).

Vascular Anatomy

For the cross-sectional xylem anatomical analysis, 3 cm of the basipetal end of each root or branch segment used for empirical conductivity measurements was stained with safranin (1% in 50% ethanol, Merck, Darmstadt, Germany) and 10–20 μm semi-thin disks cut with a sliding microtome (G.S.L.1, WSL, Birmensdorf, Switzerland). For stem wood anatomy the outermost 4 cm of the increment core were used. Photographs of the cross-sectional cuts were taken with a stereo-microscope with an automatic stage equipped with a digital camera (SteREOV20, Carl Zeiss MicroImaging GmbH, Göttingen, Gemany) at 100× magnification. Per sample, 32 up to 107 single images were stitched together to obtain the whole cross-sectional area. Image processing was done with Adobe Photoshop CS6 (version 13.0.1, Adobe Systems Incorporated, USA) and ImageJ1 (version 1.47) using the particle analysis-function for estimating vessel density (VD, n mm-1), the idealized vessels diameter (d) from major (a) and minor (b) vessel radii using the equation given by White (1991) as d = [(32 × (a × b)3)/(a2 + b2)]1/4, and cumulative vessels lumen area (Alumen, m2). Single vessel diameters (d) were used to calculate the hydraulically weighted vessel diameter (dh) according to Sperry et al. (1994) as dh = ∑d4/∑d5. For these measurements all vessels of a cross section were analyzed, yielding 110 to 3,600 measured vessel per species and organ. The theoretical hydraulic conductivity (Khtheo) of a segment was calculated based on Hagen–Poiseuille’s law as Khtheo = ((π × ∑r4)/8η)×ρ, where r is the vessel radius, η the viscosity (1.002 × 10-3 Pa s) and ρ the density of water (998.2 kg m-3), both at 20°C. Theoretical sapwood area-specific hydraulic conductivity (KStheo, kg m-1 MPa-1 s-1) was obtained from Khtheo by dividing through the microscopically determined xylem cross-sectional area without bark and pit, and theoretical leaf area-specific hydraulic conductivity (KLtheo, kg m-1 MPa-1 s-1) by division of Khtheo by the total supported leaf area (AL).

Statistical Analyses

A principal-component analysis (PCA) was done to evaluate how aboveground growth performance, wood anatomical and leaf traits are associated among each other using the package CANOCO, version 4.5 (Biometris, Wageningen, the Netherlands). The matrix species factors were lumen area (Alumen), VD, hydraulically weighted vessel diameter (dh), stem basal increment (BAI), empirical (KSemp) as well as theoretical hydraulic conductivity (KStheo). All other statistical calculations were done with the R software package, version 3.1.0 (R Development Core Team, 2014). Pearson correlations were calculated for all pairwise combinations of wood anatomical properties, WD, and hydraulic traits. In case of non-linear relationships where the data are presented on a log-linear scale, the data were log10 transformed to achieve normal distribution before further statistical analyses were conducted. Comparisons of hydraulic and leaf traits among organs were conducted using mixed linear models (lme, package: ‘nlme’ and lm package: ‘stats’) with species as random factor to account for pseudoreplication. Predicted random effects and residuals of the models were checked for normal distribution and homoscedasticity using diagnosis plots and dependent variables were log-transformed and/or variance functions (varIdent or varExp) were used (Pinheiro and Bates, 2000) when necessary. Subsequently, multiple comparison tests between group means were tested post hoc with Tukey HSD tests (glht package: ‘multcomp’). In case of heteroscedasticity an adjusted statistical framework for simultaneous inference and robust covariance estimators (Herberich et al., 2010) was used to account for different variances between groups. To test the best predictor for aboveground growth performance we applied stepwise backward model selection (step.AIC, package: ‘MASS’) to identify the most parsimonious model, defined as the model with the lowest AIC (Akaike information criterion) score (Burnham and Anderson, 2002) including KStheo, WD, Nleaf, δ13C as well as species affiliation (whether it is perhumid or seasonal) and their interactions as explanatory variables.

Results

Tree Size and Aboveground Growth Performance

The variability in mean AGB between the studied species was high, ranging between 16.9 kg in T. cacao and 888.2 kg in L. leucocephala reflecting marked differences in height and diameter between the pruned T. cacao and G. sepium, and the other four shade tree species (Table 2). L. leucocephala was on average more than two times higher and larger compared to T. cacao; the other four species ranged between these two extremes (Table 2) even though all trees were planted at the same time. Stem WD varied by a factor of two across the six tree species with E. subumbrans showing the lowest WD and G. sepium and L. leucocephala showing the highest WD values. AGB was found to be a very good predictor for the annual BAI amongst all species (Figure 1). BAI numbers were thus very different across the six tree species and ranged from 6.5 and 10.2 cm2 yr-1 in T. cacao and E. subumbrans, respectively, to 68 and 88 cm2 yr-1 in D. zibethinus and L. leucocephala, respectively (Table 2).

FIGURE 1
www.frontiersin.org

FIGURE 1. Relationship between stem basal area increment (BAI) of cacao and four shade tree species and aboveground biomass (AGB). Each symbol represents mean values for each tree species (○ Th_ca; ∇ Du_zi; ⧫ Gl_se; □ Le_le; ● Gn_gn). Error bars indicate 1 SE.

Changes in Hydraulic Conductivity Along the Water Flow Path

The empirically determined sapwood area-specific hydraulic conductivity (KSemp) of root segments differed by a factor of 100 across species with E. subumbrans showing the highest values by far, whereas the smallest values were observed in roots of T. cacao (Figure 2; Table A2). The differences in KSemp of branches across species were much less pronounced (2.3 to 7.4 kg m-1 MPa-1 s-1) with G. gnemon showing the highest and T. cacao the lowest numbers. Overall, root segments always showed higher hydraulic conductivities than branches (‘lme’; p < 0.001). Furthermore, KSemp values (in both root and branch segments) were always smaller than the theoretically calculated hydraulic conductivity (KStheo) as derived from vessel diameters by Hagen–Poiseuille’s law. Even though branch and root segments around 30 cm lengths were used, probably open-cut vessels could not be avoided particularly for root segments of E. subumbrans. However, mean KSemp values reached 9–45% of respective KStheo values indicating that open-cut vessels were negligible for most species, except for G. gnemon where 50–81% of respective KStheo values were measured (Figure 2). Empirically measured and calculated specific conductivity in root segments showed a positive linear relationship (‘lme’; p < 0.001), but not for branch segments (‘lme’; p = 0.71).

FIGURE 2
www.frontiersin.org

FIGURE 2. Hydraulic characteristics – (A) empirical sapwood area-specific hydraulic conductivity (KSemp), (B) theoretically calculated sapwood area-specific hydraulic conductivity (KStheo), (C) vessel diameter (d), (D) hydraulically weighted vessel diameter (dh), (E) vessel density (VD), and (F) relative lumen area (Alumen) – of six cacao agroforestry species (Th_ca; Du_zi; Gn_gn; Gl_se; Le_le; Er_su) among root (white bars), stem (gray bars) and branch xylem (black bars). Values are means ± SE.

Anatomical Differences across Species in Root, Stem, and Branch Wood Properties

We found considerable variation in wood anatomical and derived hydraulic traits along the flow path from root, to stem and branch wood for all six species. Exemplary pictures for this variation from three of the species are given in Figure 3. In four of the six species average vessel diameter (d) was significantly largest in the stem and not in the root wood; in the remaining two species d was comparable between root and stem wood (Figure 2). Along the flow path smallest vessels were always observed in the branch wood of all species with the exception of D. zibethinus (Figure 2). The same pattern was observed for the hydraulically weighted vessel diameter (dh) for branch wood, while the differences in dh between root and stem wood were only significant in T. cacao and D. zibethinus. In general, several wood anatomical and derived hydraulic traits allowed a grouping between the three perhumid tree species originating from strictly wet tropical environments, and the three seasonal tree species reported to tolerate moderate droughts. As mentioned above, dh was not significantly higher in stem than in root wood for the three seasonal tree species, and VD was comparable between root and branch wood and did not differ significantly. On the other hand, highest vessel densities were observed in the branch wood of all perhumid tree species, although differences were only significant in two of the three species. However, when comparing the two groups (perhumid vs. seasonal) significant differences were found (‘lme’; p < 0.001). In general, VD varied considerably between the organs and species as well and was found to decrease in the order branch – root – stem across all six species (Figure 2). Variation in VD numbers was lowest (factor < 2) in the root xylem and highest (factor > 10) in the stem xylem. VD decreased exponentially with increasing vessels diameter; we therefore concentrate on changes in d along the flow path in the following (Figure 4).

FIGURE 3
www.frontiersin.org

FIGURE 3. Cross-sections of different tree parts along the flow path: branch (left row), roots (middle row), and stems (right row) for three common tree species from cocoa agroforests in Sulawesi, Indonesia. Erythrina subumbrans (upper line), Theobroma cacao (middle line), and Gliricidia sepium (lower line). The scale bars are presented in the figures and black bars represent 1000 μm.

FIGURE 4
www.frontiersin.org

FIGURE 4. Mean vessel diameter in relation to vessel density in tree organs (roots, stems, and branches) along the flow path for the six tree species.

Relative vessel lumen area (Alumen), i.e., the ratio of lumen to sapwood area, was lowest in the stem wood in all species except of G. gnemon (Figure 2). Root and branch wood showed therefore higher Alumen values that tended to show similar dimensions with the exception of the seasonal tree species that showed significantly higher Alumen values in the root segments.

When concentrating on hydraulic properties we found a similar pattern in KStheo, where hydraulic conductivity was highest in roots of seasonal tree species (‘lme’; p < 0.001), while there is no overall significant difference between the root and stem wood in the perhumid species (p > 0.05).

Leaf Morphological, Chemical, and Hydrological Properties

Specific leaf area of sun-exposed leaves was higher by roughly a factor of two in E. subumbrans, G. sepium, and L. leucocephala than in D. zibethinus, T. cacao, and G. gnemon (Table 3). The sapwood to leaf area ratio (‘Huber value,’ HV) of the sun-lit branch samples was lowest in G. sepium and T. cacao, and highest in E. subumbrans. Empirical leaf area-specific hydraulic conductivity (KLemp) in these branch samples showed a large variation across species ranging from 1.2 10-4 kg m-1 MPa-1 s-1 in T. cacao over 4.3–4.8 10-4 kg m-1 MPa-1 s-1 in D. zibethinus, G. sepium, E. subumbrans, and L. leucocephala to 14.2 10-4 kg m-1 MPa-1 s-1 in G. gnemon. The variation in theoretical leaf area-specific hydraulic conductivity (KLtheo) derived from the wood anatomical properties was less pronounced. Lowest values were found in T. cacao, followed by G. sepium, while the other four tree species had ca. 2–5 times higher KLtheo values (Table 3).

TABLE 3
www.frontiersin.org

TABLE 3. Leaf morphological, hydraulic, and chemical properties of the six investigated tree species.

Mass-specific foliar nitrogen concentration (Nleaf) was lowest in T. cacao, medium high in D. zibethinus and G. gnemon, and highest in the three seasonal species G. sepium, L. leucocephala, and E. subumbrans (Table 3). Variation in leaf carbon isotopic composition was rather small. The two species E. subumbrans and L. leucocephala revealed a ca. 1.0–1.8 higher δ13C value than the four other species that did not show significant differences in this variable.

Interrelationships between Vascular Properties, Tree Stem Growth, and Hydraulic Conductivity

A PCA on the inter-relationships between the investigated traits explained a large proportion of the total variance of the data set along the first four axes (Table 4). The first axis was strongly positively associated with all wood anatomical traits (including HV and δ13C), but negatively with WD. Axis 1 was furthermore positively related to DBH. Axis 2 was strongly associated with stem and branch lumen area as well as with the leaf traits (KLemp and Nleaf). In contrast to KSemp in root segments, branch KSemp was associated with axis 2 and therefore showed an inter-relationship with KLemp. BAI showed an only moderate association with the first axis and thus was only weakly correlated with the majority of wood anatomically and tree structural variables. BAI was correlated best with the third axis that was only associated with the variables AGB, DBH, and WD (positively), as well as root KStheo and KSemp (negatively). A Pearson’s coefficient of correlation analysis, however, revealed a strong relationship between BAI and KStheo on a species level for root, stem, and branch wood tissue (Figures 5A–C).

TABLE 4
www.frontiersin.org

TABLE 4. Results of a Principal Components Analysis (PCA) on the response of six agroforestry tree species with respect to stem BAI, anatomical properties of the coarse root, stem and branch wood as well as hydraulic and leaf traits.

FIGURE 5
www.frontiersin.org

FIGURE 5. Relationship between stem basal area increment (BAI) of cacao and four shade tree species and theoretically calculated cross sectional sapwood area-specific hydraulic conductivity (KStheo) in the root (A), stem (B), and branch wood (C). Each symbol represents mean values for each tree species (○ Th_ca; ∇ Du_zi; ⧫ Gl_se; □ Le_le; ● Gn_gn). Error bars indicate 1 SE.

A systematic correlation analysis of pairs of traits confirmed that most of the inter-relationships identified in the PCA on a species level were also valid on a tree individual level (Table 5). BAI was strongly interrelated with AGB as well as DBH and showed moreover a significant correlation with stem dh and stem KStheo. None of these three variables were correlated with WD that generally showed only few and moreover relatively low correlations with other variables in the Pearson correlation analysis. Mixed effect models incorporating the pseudoreplication due to species confirmed that WD does not explain neither KStheo, dh nor BAI in our data well (‘lme’; p > 0.05).

TABLE 5
www.frontiersin.org

TABLE 5. Pearson’s coefficients of correlation between pairs of traits.

Contrary to significant relationships of stem wood KStheo with AGB, BAI and DBH, no relationship of KSemp between any of these traits could be found, except for KSemp in root segments that were related to DBH.

As expected, all species and organs showed a positive relationship between KStheo and dh (Table 5). Foliar nitrogen content (Nleaf) as well as the carbon isotope signature (δ13C) did not show any relation with neither leaf area-specific hydraulic conductivity (KLtheo), KSemp nor HV, but a strong significant correlation within each other. Nleaf was unrelated to BAI among species (p > 0.1, r = 0.05) also when excluding the three seasonal species.

Stepwise model selection confirmed that KStheo is the best predictor for AGB together with WD and neither Nleaf, δ13C nor KSemp were explaining the variability in our data significantly.

Discussion

Patterns in Xylem Anatomy among Species in Stem, Root, and Branch Wood

Most studies on hydraulic anatomical properties in trees describe vessel sizes to be largest in roots and basipetally taper to the branches (Tyree and Zimmermann, 2002; McElrone et al., 2004; Goncalves et al., 2007; Domec et al., 2009; Lintunen and Kalliokoski, 2010). The first notification of this paradigm even goes back to observations by Nehemiah Grew in the 17th century (Baas, 1982). Generally, variation in conduit diameter is a compromise between hydraulic efficiency, safety, and the maximization of conductivity per growth investment due to conduit tapering (Sperry et al., 2008, 2012). As concluded by Tyree and Zimmermann (2002) the reason for conduit tapering is the control of water distribution, and more importantly to counter increases in flow resistance and gravimetrically forces with tree height to minimize the increasing risk of xylem dysfunction with path length (West et al., 1999; Anfodillo et al., 2006; Sperry et al., 2012). This is additionally mirrored in the meanwhile commonly observed relation between conduit size and vulnerability to cavitation (Wheeler et al., 2005; Maherali et al., 2006; Cai et al., 2010; Domec et al., 2010; Sterck et al., 2012). It is argued that in roots water stress will not be as great as in stems, since water potentials drop in going from root to stem to leaf (Tyree and Zimmermann, 2002). As long as soil water is still extractable, roots would then be less exposed to drought-induced embolism and might therefore afford larger vessels. Furthermore, with a small need of structural support and embedded in a soil matrix, biomechanical stress is unlikely to affect roots (McElrone et al., 2004; Pratt et al., 2007). Additionally, most plants have developed a mechanism to restore vessel functionality by refilling embolized vessels through living rays and paratracheal parenchyma. The contribution of paratracheal parenchyma was recently shown for grapevine by Brodersen et al. (2010), while the molecular and physiological paths were investigated by Chitarra et al. (2014). As coarse roots are located close to the water source it seems beneficial to restrict hydraulic failure to areas within the complex hydraulic network that are easily refilled. Embolism reversal is thought to occur by active transport of sugars into empty conduits, which are generally accumulated in high amounts within a trees rooting system (Würth et al., 2005).

Our results are in contrast to the common assumption as we found the largest vessels along the flow path in the stem xylem and not in the roots. Supporting our first hypothesis, our results are in accordance with the observations of a hump-shaped vessel size distribution along the flow path found in tropical trees of Indonesia (Schuldt et al., 2013), supported by findings from South America where the largest vessels were observed in the stem wood (Machado et al., 2007) or comparing just root and branch wood both organs showed similar vessel sizes (Fortunel et al., 2013). Our confirmative finding could represent a response to permanent water availability and low evaporative demand in this humid region, where trees without severe drought limitation might have developed roots with large relative lumen area and less structural tissue that can achieve sufficiently high axial conductivities in these organs. Thereby they would compensate for the smaller vessel diameters in roots than in the trunk in accordance with the pipe model theory by Shinozaki et al. (1964). Originally, this theory attempted to explain plant architecture in a quantitative way by proposing that photosynthetic organs should be supported by an adequate structure of non-photosynthetic organs in order to guarantee functionality (Chiba, 1998). Nevertheless, one has to keep in mind that the pipe model theory is not a hydraulic model, but should be viewed as a biomass allometry model with no particular implications concerning either hydraulics or biomechanics as proposed by McCulloh and Sperry (2005).

Machado et al. (2007) argued that shallow roots of moist tropical forest species, which is a common rooting pattern in tropical moist forests (Leuschner et al., 2006; Hertel et al., 2009), are subjected to variations in water availability and the narrower vessels in the root wood are a safety trade-off against cavitation. However, the vessel sizes found in coarse roots in the present study as well as in other tropical studies (Machado et al., 2007; Fortunel et al., 2013; Schuldt et al., 2013) are at least 30% larger compared to, e.g., temperate forest tree species (Köcher et al., 2012), and therefore might not directly be rated as an increased safety against cavitation compared to the stem or branch wood. It seems rational to assume that not the coarse root xylem, despite large vessel sizes, represents the most drought-sensitive organ, but rather that of fine roots with smaller diameter. In contrast to coarse roots, which are primarily responsible for axial water transport, fine roots represent the highest resistance for water transport within the rooting system due to radial water flow either along the apoplastic or cellular pathway (Steudle, 2000). As the most distal organs they are scarified in response to drought in order to avoid serious harm to coarse and large roots like it has been observed in various temperate and boreal forests (Gaul et al., 2008; Chenlemuge et al., 2013; Hertel et al., 2013). Fine roots might thereby act as a sort of ‘hydraulic fuse,’ which evolved from Zimmermann’s segmentation hypothesis (Tyree and Zimmermann, 2002) in analogy to the leaf petiole (Zufferey et al., 2011). At the root level this ‘hydraulic segmentation’ might additionally protect the below-ground system preventing the reverse water flow from main to lateral roots and back to the dryer soil as discussed for grapevine by Lovisolo et al. (2008). Woody plants would accordingly restrict hydraulic failure to redundant organs that are readily replaced (Sperry et al., 2002) although it has been argued that the term ‘hydraulic fuse’ should be reconsidered since roots are not necessarily an ‘expandable’ organ (Gonzalez-Benecke et al., 2010). The construction costs of fine roots and lignified small-diameter roots in term of carbon and nutrients may not be much smaller than for twigs and leaves, and the loss of roots is directly related to the loss of absorption capacity for nutrients and water. However, comparable data on fine root mortality and percentage loss of conductivity in coarse roots that would empirically support the idea that fine root are sacrificed in order to protect the hydraulic system are to our knowledge not available so far.

Concurrently with decreasing vessel size, conduit frequency is generally reported to increase from the roots to the branches (e.g., Lintunen and Kalliokoski, 2010). This commonly observed trade-off between VD and vessel diameter (Preston et al., 2006; Sperry et al., 2008; Zanne et al., 2010) could not be confirmed in our study where the stem wood showed by far the lowest VD compared to roots and branches. Since flow in capillary systems is proportional to the fourth power of vessel radius according to Hagen–Poiseuille law, variations in vessel diameter has a much greater effect on KStheo than variations in VD. However, due to the occurrence of a higher relative vessel lumen area and a few large vessels in roots we have observed that specific conductivity in the three seasonal species, i.e., G. sepium, E. subumbrans, and L. leucocephala, was higher in roots than in stems, even though the largest vessels were observed in the stem wood. At least for E. subumbrans and L. leucocephala this pattern could additionally be explained by the highest relative vessel lumen area found in roots, i.e., less space was occupied by tracheids or fibers compared to the stem or branch wood. Furthermore, disproportionately high empirical conductance measured in E. subumbrans might be attributed to the presence of open-cut vessels, which are highly conductive as water does not have to pass pit membranes, which generally account for more than 50% of the total hydraulic resistance (Choat et al., 2008). However, since tree hydraulic traits have been associated with general habitat preferences of various species (Sperry, 2000; Maherali et al., 2004), this finding could be due to the biogeographic background and could represent genetically determined adaptations to different water availability in the natural habitat of the species. While T. cacao, D. zibethinus and G. gnemon are known to be strictly wet tropical forest trees rather sensitive to drought and low air humidity (Brown, 1997; Carr and Lockwood, 2011), the other three species are reported to be fast-growing drought-resistant trees (Mrema et al., 1997; Fagbola et al., 2001). While drought resistance is recorded for some Erythrina species (da Silva et al., 2010; Manoharan et al., 2010), not many data are available on E. subumbrans, which is a species native to Indonesia. In a habitat where water stress is generally modest or absent such as the humid climate in Sulawesi, cavitation-avoiding mechanisms might be less beneficial than hydraulic efficiency and largest vessels can occur in stem xylem, thereby reducing the hydraulic resistance along the flow path. However, it remains speculative why the largest vessels along the flow path are observed in the root xylem only in biomes that frequently experience either drought- or frost stress. The size of a vessel is thought to be caused by the concentration of the plant hormone indole-3-acetic acid (IAA), an endogenous auxin, at the time of cell differentiation (Aloni, 1987; Lovisolo et al., 2002), which is also related to the cambial age and related cambial activity as seen by the radial increase in vessel size at the stem base of a tree (Spicer and Gartner, 2001). It would thus be of interest to extent the results of the present study to a quantification of IAA concentration in both coarse root and stem cambium in tropical and temperate trees; the latter should show higher concentrations in the root xylem independently of cambial age in agreement with the common paradigm that largest vessels are found in the rooting system.

Relationships between Vascular Properties, Tree Stem Growth, and Hydraulic Conductivity

Wood density is an easy to measure functional wood property that has been linked to various ecological and other functional traits. In species showing a relatively large fraction of vessels close to the hydraulically weighted mean vessel diameter (dh), KStheo should correlate negatively with WD (Bucci et al., 2004; Meinzer et al., 2008; Gonzalez-Benecke et al., 2010). Similar to observations on tropical forest trees from perhumid tropical environments (Poorter et al., 2010; Schuldt et al., 2013) we expected WD to be unrelated to wood anatomical and hydraulic properties. Even though we found a correlation between WD and dh on tree level, this relationship could not be confirmed on species level accounting for species pseudoreplication in mixed effect models. Also we found no significant relationship of WD to basal stem area increment, contradicting former results on a close relation between WD and growth for tropical trees (King et al., 2006; Poorter et al., 2010; Hietz et al., 2013). Several other studies report WD to be partially decoupled from hydraulic conductivity due to variation of frequency and size of fibers in angiosperms (Preston et al., 2006; Martinez-Cabrera et al., 2009; Zanne et al., 2010). Results on the relationship between WD and vascular properties as well as tree growth are thus partly conflicting; while some studies confirm that WD varies inversely with vessel size (Preston et al., 2006; Jacobsen et al., 2007; Thomas et al., 2007; McCulloh et al., 2011; Gleason et al., 2012), others did not support this finding (Martinez-Cabrera et al., 2009; Poorter et al., 2010; Russo et al., 2010; Zanne et al., 2010; Fan et al., 2012). These contradicting results are indicating that the relation between WD and growth or vessel traits is not necessarily interrelated and should be viewed separately. We further suspect that the relation between wood properties and tree hydraulics may depend as well on biogeographical origin and drought-adaptation strategy of the species investigated since convergent environmental factors such as water availability are known to lead to adaptations in functional wood anatomical properties (Swenson and Enquist, 2007; Gleason et al., 2013; Richardson et al., 2013).

We found wood anatomical and derived hydraulic properties to be a much better predictor for tree stem growth performance than WD as KStheo of the all tree organs studied were strongly positively correlated with stem BAI on a species level. This is in accordance with a growing body of studies showing strong links between growth rate and wood anatomical traits (Zhang and Cao, 2009; Poorter et al., 2010; Russo et al., 2010; Fan et al., 2012). In contrast, neither empirically measured branch and root KSemp, nor foliar δ13C or foliar nitrogen content were good predictors for aboveground growth performance.

We expected to find close correlations between functional leaf traits assumed to be associated with high aboveground productivity, i.e., high foliar N content and more negative foliar δ13C, and stem increment in our samples. However, no such correlation was found. This is most likely explained by the fact that our sampled species contained several N-fixing legume species, our relatively low species number as well as due to the fact that our study was conducted in a perhumid region were drought stress is not to be expected.

Conclusion

Our study results suggest that even though vessel traits, growth performance, and WD relations follow distinct conceptually determined trade-offs, some of these long-established paradigms might not be uniformly applicable to tree species from all biogeographic regions presumably due to their varying drought adaptation strategies. In moist tropical environments we could not confirm the paradigm of continuous conduit tapering from roots to branches although some traits (VD, relative vessel lumen area, and theoretical sapwood area-specific conductivity) enabled a clear separation between the three strictly wet tropical species and the three seasonal tree species. We therefore expect patterns in vessel traits along the flow path from roots to branches to be dependent on the long-term precipitation regime at the biogeographic origin of the investigated tree species. Furthermore and contrary to common knowledge, the investigated tree species did neither show a relationship between aboveground growth performance and WD nor foliar nitrogen content, nor between WD and vessel size. Instead, we found growth rate to be closely linked with wood anatomical and derived hydraulic traits. Future research should thus include a systematic approach to different biogeographic regions and cover a wider range of ecosystem types particularly underrepresented biomes.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This study was financed by the German Research Foundation (DFG) in the framework of the collaborative German – Indonesian research project CRC990 (‘Ecological and Socioeconomic Functions of Tropical Lowland Rainforest Transformation Systems (Sumatra, Indonesia),’ subproject B04 and ELUC (‘Environmental and land-use change in Sulawesi, Indonesia’), subproject ‘Shade trees in cacao agroforestry systems: influence on roots and net primary production.’ We acknowledge support by the Open Access Publication Funds of the Göttingen University. We further thank village heads and local plot owners for granting us access to and use of their properties as well as we are very thankful to our counterparts and local assistance for support of our study, and the three reviewers who provided helpful suggestions for improving the manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fpls.2015.00191/abstract

Footnote

  1. ^ http://rsb.info.nih.gov/ij

References

Achard, F., Eva, H. D., Stibig, H. J., Mayaux, P., Gallego, J., Richards, T.,et al. (2002). Determination of deforestation rates of the world’s humid tropical forests. Science 297, 999–1002. doi: 10.1126/science.1070656

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Aloni, R. (1987). Differentiation of Vascular Tissues. Annu. Rev. Plant Physiol. Plant Mol. Biol. 38, 179–204. doi: 10.1146/annurev.arplant.38.1.179

CrossRef Full Text | Google Scholar

Anfodillo, T., Carraro, V., Carrer, M., Fior, C., and Rossi, S. (2006). Convergent tapering of xylem conduits in different woody species. New Phytol. 169, 279–290. doi: 10.1111/j.1469-8137.2005.01587

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Anfodillo, T., Petit, G., and Crivellaro, A. (2013). Axial conduit widening in woody species: a still neglected anatomical pattern. IAWA J. 34, 352–364. doi: 10.1163/22941932-00000030

CrossRef Full Text | Google Scholar

Awad, H., Barigah, T., Badel, E., Cochard, H., and Herbette, S. (2010). Poplar vulnerability to xylem cavitation acclimates to drier soil conditions. Physiol. Plant. 139, 280–288. doi: 10.1111/j.1399-3054.2010.01367

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Baas, P. (1982). Systematic, phylogenetic and ecological wood anatomy-history and perspectives. Acta Bot. Neerl. 31, 501–501.

Google Scholar

Brodersen, C. R., McElrone, A. J., Choat, B., Matthews, M. A., and Shackel, K. A. (2010). The dynamics of embolism repair in xylem: in vivo visualizations using high-resolution computed tomography. Plant Physiol. 154, 1088–1095. doi: 10.1104/pp.110.162396

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Brodribb, T. J., Holbrook, N. M., and Gutierrez, M. V. (2002). Hydraulic and photosynthetic co-ordination in seasonally dry tropical forest trees. Plant Cell Environ. 25, 1435–1444. doi: 10.1046/j.1365-3040.2002.00919

CrossRef Full Text | Google Scholar

Brown, S. (1997). Estimating Biomass and Biomass Change of Tropical Forests: A Primer. Rome: FAO.

Google Scholar

Bucci, S. J., Goldstein, G., Meinzer, F. C., Scholz, F. G., Franco, A. C., and Bustamante, M. (2004). Functional convergence in hydraulic architecture and water relations of tropical savanna trees: from leaf to whole plant. Tree Physiol. 24, 891–899. doi: 10.1093/treephys/24.8.891

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Burnham, K. P., and Anderson, D. R. (2002). Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. New York, NY: Springer.

Google Scholar

Cai, J., Hacke, U., Zhang, S., and Tyree, M. T. (2010). What happens when stems are embolized in a centrifuge? Testing the cavitron theory. Physiol. Plant. 140, 311–320. doi: 10.1111/j.1399-3054.2010.01402

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Carlquist, S. (1994). Wood and bark anatomy of Gnetum gnemon L. Bot. J. Linn. Soc. 116, 203–221. doi: 10.1006/bojl.1994.1060

CrossRef Full Text | Google Scholar

Carr, M. K. V., and Lockwood, G. (2011). The water relations and irrigation requirements of cocoa (Theobroma cacao L.). Rev. Exp. Agric. 47, 653–676. doi: 10.1017/s0014479711000421

CrossRef Full Text | Google Scholar

Chave, J., Andalo, C., Brown, S., Cairns, M. A., Chambers, J. Q., Eamus, D.,et al. (2005). Tree allometry and improved estimation of carbon stocks and balance in tropical forests. Oecologia 145, 87–99. doi: 10.1007/s00442-005-0100

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Chenlemuge, T., Hertel, D., Dulamsuren, C., Khishigjargal, M., Leuschner, C., and Hauck, M. (2013). Extremely low fine root biomass in Larix sibirica forests at the southern drought limit of the boreal forest. Flora 208, 488–496. doi: 10.1016/j.flora.2013.08.002

CrossRef Full Text | Google Scholar

Chiba, Y. (1998). Architectural analysis of relationship between biomass and basal area based on pipe model theory. Ecol. Modell. 108, 219–225. doi: 10.1016/s0304-3800(98)00030-1

CrossRef Full Text | Google Scholar

Chitarra, W., Balestrini, R., Vitali, M., Pagliarani, C., Perrone, I., Schubert, A.,et al. (2014). Gene expression in vessel-associated cells upon xylem embolism repair in Vitis vinifera L. Planta 239, 887–899. doi: 10.1007/s00425-013-2017-7

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Choat, B., Cobb, A. R., and Jansen, S. (2008). Structure and function of bordered pits: new discoveries and impacts on whole-plant hydraulic function. New Phytol. 177, 608–625. doi: 10.1111/j.1469-8137.2007.02317.x

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

da Silva, E. C., Silva, M. F. A., Nogueira, R. J. M. C., and Albuquerque, M. B. (2010). Growth evaluation and water relations of Erythrina velutina seedlings in response to drought stress. Braz. J. Plant Physiol. 22, 225–233. doi: 10.1590/S1677-04202010000400002

CrossRef Full Text | Google Scholar

Domec, J.-C., Schaefer, K., Oren, R., Kim, H. S., and Mccarthy, H. R. (2010). Variable conductivity and embolism in roots and branches of four contrasting tree species and their impacts on whole-plant hydraulic performance under future atmospheric CO2 concentration. Tree Physiol. 30, 1001–1015. doi: 10.1093/treephys/tpq054

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Domec, J.-C., Warren, J. M., Meinzer, F. C., and Lachenbruch, B. (2009). Safety factors for xylem failure by implosion and air-seeding within roots, trunks and branches of young and old conifer trees. IAWA J. 30, 101–120. doi: 10.1163/22941932-90000207

CrossRef Full Text | Google Scholar

Enquist, B. J., West, G. B., Charnov, E. L., and Brown, J. H. (1999). Allometric scaling of production and life-history variation in vascular plants. Nature 401, 907–911. doi: 10.1038/44819

CrossRef Full Text | Google Scholar

Espino, S., and Schenk, H. J. (2011). Mind the bubbles: achieving stable measurements of maximum hydraulic conductivity through woody plant samples. J. Exp. Bot. 62, 1119–1132. doi: 10.1093/jxb/erq338

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Fagbola, O., Osonubi, O., Mulongoy, K., and Odunfa, S. A. (2001). Effects of drought stress and arbuscular mycorrhiza on the growth of Gliricidia sepium (Jacq). Walp, and Leucaena leucocephala (Lam.) de Wit. in simulated eroded soil conditions. Mycorrhiza 11, 215–223. doi: 10.1007/s005720100114

CrossRef Full Text | Google Scholar

Fan, Z.-X., Zhang, S.-B., Hao, G.-Y., Slik, J. W. F., and Cao, K.-F. (2012). Hydraulic conductivity traits predict growth rates and adult stature of 40 Asian tropical tree species better than wood density. J. Ecol. 100, 732–741. doi: 10.1111/j.1365-2745.2011.01939

CrossRef Full Text | Google Scholar

FAO, and JRC. (2012). “Global forest land-use change 1990–2005,” in FAO Forestry Paper No. 169, eds E. J. Lindquist, R. D’annunzio, A. Gerrand, K. Macdicken, F. Achard, R. Beuchle et al. (Rome: Food and Agriculture Organization of the United Nations and European Commission Joint Research Centre).

Google Scholar

Fisher, J. B., and Ewers, F. W. (1995). Vessel dimensions in liana and tree species of vessel dimensions Gnetum (Gnetales). Am. J. Bot. 82, 1350–1357. doi: 10.2307/2445861

CrossRef Full Text

Fortunel, C., Ruelle, J., Beauchêne, J., Fine, P. V. A., and Baraloto, C. (2013). Wood specific gravity and anatomy of branches and roots in 113 Amazonian rainforest tree species across environmental gradients. New Phytol. 202, 79–94. doi: 10.1111/nph.12632

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Gaul, D., Hertel, D., Borken, W., Matzner, E., and Leuschner, C. (2008). Effects of experimental drought on the fine root system of mature Norway spruce. For. Ecol. Manage. 256, 1151–1159. doi: 10.1016/j.foreco.2008.06.016

CrossRef Full Text | Google Scholar

Gleason, S. M., Butler, D. W., and Waryszak, P. (2013). Shifts in leaf and stem hydraulic traits across aridity gradients in Eastern Australia. Int. J. Plant Sci. 174, 1292–1301. doi: 10.1086/673239

CrossRef Full Text | Google Scholar

Gleason, S. M., Butler, D. W., Zieminska, K., Waryszak, P., and Westoby, M. (2012). Stem xylem conductivity is key to plant water balance across Australian angiosperm species. Funct. Ecol. 26, 343–352. doi: 10.1111/j.1365-2435.2012.01962

CrossRef Full Text | Google Scholar

Goncalves, B., Correia, C. M., Silva, A. P., Bacelar, E. A., Santos, A., Ferreira, H.,et al. (2007). Variation in xylem structure and function in roots and stems of scion-rootstock combinations of sweet cherry tree (Prunus avium L.). Trees 21, 121–130. doi: 10.1007/s00468-006-0102-2

CrossRef Full Text | Google Scholar

Gonzalez-Benecke, C. A., Martin, T. A., and Peter, G. F. (2010). Hydraulic architecture and tracheid allometry in mature Pinus palustris and Pinus elliottii trees. Tree Physiol. 30, 361–375. doi: 10.1093/treephys/tpp129

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Hajek, P., Leuschner, C., Hertel, D., Delzon, S., and Schuldt, B. (2014). Trade-offs between xylem hydraulic properties, wood anatomy and yield in Populus. Tree Physiol. 34, 744–756. doi: 10.1093/treephys/tpu048

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Hansen, M. C., Stehman, S. V., Potapov, P. V., Loveland, T. R., Townshend, J. R. G., Defries, R. S.,et al. (2008). Humid tropical forest clearing from 2000 to 2005 quantified by using multitemporal and multiresolution remotely sensed data. Proc. Natl. Acad. Sci. U.S.A. 105, 9439–9444. doi: 10.1073/pnas.0804042105

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Herberich, E., Sikorski, J., and Hothorn, T. (2010). A robust procedure for comparing multiple means under heteroscedasticity in unbalanced designs. PLoS ONE 5:e9788. doi: 10.1371/journal.pone.0009788

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Hertel, D., Moser, G., Culmsee, H., Erasmi, S., Horna, V., Schuldt, B.,et al. (2009). Below-and above-ground biomass and net primary production in a paleotropical natural forest (Sulawesi, Indonesia) as compared to neotropical forests. For. Ecol. Manage. 258, 1904–1912. doi: 10.1016/j.foreco.2009.07.019

CrossRef Full Text | Google Scholar

Hertel, D., Strecker, T., Müller-Haubold, H., and Leuschner, C. (2013). Fine root biomass and dynamics in beech forests across a precipitation gradient–is optimal resource partitioning theory applicable to water-limited mature trees? J. Ecol. 101, 1183–1200. doi: 10.1111/1365-2745.12124

CrossRef Full Text | Google Scholar

Hietz, P., Valencia, R., and Joseph Wright, S. (2013). Strong radial variation in wood density follows a uniform pattern in two neotropical rain forests. Funct. Ecol. 27, 684–692. doi: 10.1111/1365-2435.12085

CrossRef Full Text | Google Scholar

Hoeber, S., Leuschner, C., Köhler, L., Arias-Aguilar, D., and Schuldt, B. (2014). The importance of hydraulic conductivity and wood density to growth performance in eight tree species from a tropical semi-dry climate. For. Ecol. Manage. 330, 126–136. doi: 10.1016/j.foreco.2014.06.039

CrossRef Full Text | Google Scholar

Jacobsen, A. L., Ewers, F. W., Pratt, R. B., Paddock, W. A., and Davis, S. D. (2005). Do xylem fibers affect vessel cavitation resistance? Plant Physiol. 139, 546–556. doi: 10.1104/pp.104.058404

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Jacobsen, A. L., Pratt, R. B., Ewers, F. W., and Davis, S. D. (2007). Cavitation resistance among 26 chaparral species of southern California. Ecol. Monogr. 77, 99–115. doi: 10.1890/05–1879

CrossRef Full Text | Google Scholar

King, D. A., Davies, S. J., Supardi, M. N. N., and Tan, S. (2005). Tree growth is related to light interception and wood density in two mixed dipterocarp forests of Malaysia. Funct. Ecol. 19, 445–453. doi: 10.1111/j.1365-2435.2005.00982

CrossRef Full Text | Google Scholar

King, D. A., Davies, S. J., Tan, S., and Noor, N. S. M. (2006). The role of wood density and stem support costs in the growth and mortality of tropical trees. J. Ecol. 94, 670–680. doi: 10.1111/j.1365-2745.2006.01112

CrossRef Full Text | Google Scholar

Köcher, P., Horna, V., Beckmeyer, I., and Leuschner, C. (2012). Hydraulic properties and embolism in small-diameter roots of five temperate broad-leaved tree species with contrasting drought tolerance. Ann. For. Sci. 69, 693–703. doi: 10.1007/s13595-012-0189-0

CrossRef Full Text | Google Scholar

Leuschner, C., Wiens, M., Harteveld, M., Hertel, D., and Tjitrosemito, S. (2006). Patterns of fine root mass and distribution along a disturbance gradient in a tropical montane forest, central Sulawesi (Indonesia). Plant Soil 283, 163–174. doi: 10.1007/s11104-006-6638-5

CrossRef Full Text | Google Scholar

Lintunen, A., and Kalliokoski, T. (2010). The effect of tree architecture on conduit diameter and frequency from small distal roots to branch tips in Betula pendula, Picea abies and Pinus sylvestris. Tree Physiol. 30, 1433–1447. doi: 10.1093/treephys/tpq085

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Lovisolo, C., Perrone, I., Hartung, W., and Schubert, A. (2008). An abscisic acid-related reduced transpiration promotes gradual embolism repair when grapevines are rehydrated after drought. New Phytol. 180, 642–651. doi: 10.1111/j.1469-8137.2008.02592.x

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Lovisolo, C., Schubert, A., and Sorce, C. (2002). Are xylem radial development and hydraulic conductivity in downwardly-growing grapevine shoots influenced by perturbed auxin metabolism? New Phytol. 156, 65–74. doi: 10.1046/j.1469-8137.2002.00492.x

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Machado, S. R., Rodella, R. A., Angyalossy, V., and Marcati, C. R. (2007). Structural variations in root and stem wood of Styrax (Styracaceae) from Brazilian forest and cerrado. IAWA J. 28, 173–188. doi: 10.1163/22941932-90001632

CrossRef Full Text | Google Scholar

Maherali, H., Moura, C. F., Caldeira, M. C., Willson, C. J., and Jackson, R. B. (2006). Functional coordination between leaf gas exchange and vulnerability to xylem cavitation in temperate forest trees. Plant Cell Environ. 29, 571–583. doi: 10.1111/j.1365-3040.2005.01433

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Maherali, H., Pockman, W. T., and Jackson, R. B. (2004). Adaptive variation in the vulnerability of woody plants to xylem cavitation. Ecology 85, 2184–2199. doi: 10.1890/02-0538

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Manoharan, P. T., Shanmugaiah, V., Balasubramanian, N., Gomathinayagam, S., Sharma, M. P., and Muthuchelian, K. (2010). Influence of AM fungi on the growth and physiological status of Erythrina variegata Linn. grown under different water stress conditions. Eur. J. Soil Biol. 46, 151–156. doi: 10.1016/j.ejsobi.2010.01.001

CrossRef Full Text | Google Scholar

Martinez-Cabrera, H. I., Jones, C. S., Espino, S., and Schenk, H. J. (2009). Wood anatomy and wood density in shrubs: responses to varying aridity along transcontinental transects. Am. J. Bot. 96, 1388–1398. doi: 10.3732/ajb.0800237

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Martinez-Vilalta, J., Prat, E., Oliveras, I., and Pinol, J. (2002). Xylem hydraulic properties of roots and stems of nine Mediterranean woody species. Oecologia 133, 19–29. doi: 10.1007/s00442-002-1009-2

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

McCulloh, K. A., Meinzer, F. C., Sperry, J. S., Lachenbruch, B., Voelker, S. L., Woodruff, D. R.,et al. (2011). Comparative hydraulic architecture of tropical tree species representing a range of successional stages and wood density. Oecologia 167, 27–37. doi: 10.1007/s00442-011-1973-5

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

McCulloh, K. A., and Sperry, J. S. (2005). The evaluation of Murray’s law in Psilotum nudum (Psilotaceae), an analogue of ancestral vascular plants. Am. J. Bot. 92, 985–989. doi: 10.3732/ajb.92.6.985

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

McCulloh, K. A., Sperry, J. S., and Adler, F. R. (2003). Water transport in plants obeys Murray’s law. Nature 421, 939–942. doi: 10.1038/nature01444

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

McCulloh, K. A., Sperry, J. S., Lachenbruch, B., Meinzer, F. C., Reich, P. B., and Voelker, S. L. (2010). Moving water well: comparing hydraulic efficiency in twigs and trunks of coniferous, ring-porous, and diffuse-porous saplings from temperate and tropical forests. New Phytol. 186, 439–450. doi: 10.1111/j.1469-8137.2010.03181

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

McElrone, A. J., Pockman, W. T., Martinez-Vilalta, J., and Jackson, R. B. (2004). Variation in xylem structure and function in stems and roots of trees to 20 m depth. New Phytol. 163, 507–517. doi: 10.1111/j.1469-8137.2004.01127

CrossRef Full Text | Google Scholar

Meinzer, F. C., Campanello, P. I., Domec, J.-C., Gatti, M. G., Goldstein, G., Villalobos-Vega, R.,et al. (2008). Constraints on physiological function associated with branch architecture and wood density in tropical forest trees. Tree Physiol. 28, 1609–1617. doi: 10.1093/treephys/28.11.1609

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Meinzer, F. C., James, S. A., Goldstein, G., and Woodruff, D. (2003). Whole-tree water transport scales with sapwood capacitance in tropical forest canopy trees. Plant Cell Environ. 26, 1147–1155. doi: 10.1046/j.1365-3040.2003.01039

CrossRef Full Text | Google Scholar

Moser, G., Leuschner, C., Hertel, D., Hoelscher, D., Koehler, M., Leitner, D.,et al. (2010). Response of cocoa trees (Theobroma cacao) to a 13-month desiccation period in Sulawesi, Indonesia. Agroforest. Syst. 79, 171–187. doi: 10.1007/s10457-010-9303–9301

CrossRef Full Text | Google Scholar

Motamayor, J. C., Lachenaud, P., Mota, J., Loor, R., Kuhn, D. N., Brown, J. S.,et al. (2008). Geographic and genetic population differentiation of the Amazonian chocolate tree (Theobroma cacao L). PLoS ONE 3:e3311. doi: 10.1371/journal.pone.0003311

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Mrema, A. F., Granhall, U., and Sennerbyforsse, L. (1997). Plant growth, leaf water potential, nitrogenase activity and nodule anatomy in Leucaena leucocephala as affected by water stress and nitrogen availability. Trees 12, 42–48. doi: 10.1007/pl00009695

CrossRef Full Text | Google Scholar

Pinheiro, J. C., and Bates, D. M. (2000). Mixed-Effects Models in S and S-PLUS. New York, NY: Springer. doi: 10.1007/978-1-4419-0318-1

CrossRef Full Text | Google Scholar

Poorter, L., Mcdonald, I., Alarcon, A., Fichtler, E., Licona, J.-C., Pena-Claros, M.,et al. (2010). The importance of wood traits and hydraulic conductance for the performance and life history strategies of 42 rainforest tree species. New Phytol. 185, 481–492. doi: 10.1111/j.1469-8137.2009.03092

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Pratt, R. B., Jacobsen, A. L., Ewers, F. W., and Davis, S. D. (2007). Relationships among xylem transport, biomechanics and storage in stems and roots of nine Rhamnaceae species of the California chaparral. New Phytol. 174, 787–798. doi: 10.1111/j.1469-8137.2007.02061

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Preston, K. A., Cornwell, W. K., and Denoyer, J. L. (2006). Wood density and vessel traits as distinct correlates of ecological strategy in 51 California coast range angiosperms. New Phytol. 170, 807–818. doi: 10.1111/j.1469-8137.2006.01712

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

R Development Core Team. (2014). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. Available at: http://www.R-project.org

Google Scholar

Rice, R. A., and Greenberg, R. (2000). Cacao cultivation and the conservation of biological diversity. Ambio 29, 167–173.

PubMed Abstract | Full Text | Google Scholar

Richardson, S. J., Allen, R. B., Buxton, R. P., Easdale, T. A., Hurst, J. M., Morse, C. W.,et al. (2013). Intraspecific relationships among wood density, leaf structural traits and environment in four co-occurring species of Nothofagus in New Zealand. PLoS ONE 8:e58878. doi: 10.1371/journal.pone.0058878

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Russo, S. E., Jenkins, K. L., Wiser, S. K., Uriarte, M., Duncan, R. P., and Coomes, D. A. (2010). Interspecific relationships among growth, mortality and xylem traits of woody species from New Zealand. Funct. Ecol. 24, 253–262. doi: 10.1111/j.1365-2435.2009.01670

CrossRef Full Text | Google Scholar

Schroth, G., and Harvey, C. A. (2007). Biodiversity conservation in cocoa production landscapes: an overview. Biodivers. Conserv. 16, 2237–2244. doi: 10.1007/s10531-007-9195-1

CrossRef Full Text | Google Scholar

Schuldt, B., Leuschner, C., Brock, N., and Horna, V. (2013). Changes in wood density, wood anatomy and hydraulic properties of the xylem along the root-to-shoot flow path in tropical rainforest trees. Tree Physiol. 33, 161–174. doi: 10.1093/treephys/tps122

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Shinozaki, K., Yoda, K., Hozumi, K., and Kira, T. (1964). A quantitative analysis of plant form-the pipe model theory II. Further evidence of the theory and its application in forest ecology. Jpn. J. Ecol. 14, 133–139.

Google Scholar

Sperry, J. S. (2000). Hydraulic constraints on plant gas exchange. Agric. For. Meteorol. 104, 13–23. doi: 10.1016/s0168-1923(00)00144-1

CrossRef Full Text | Google Scholar

Sperry, J. S., Donnelly, J. R., and Tyree, M. T. (1988). A method for measuring hydraulic conductivity and embolism in xylem. Plant Cell Environ. 11, 35–40. doi: 10.1111/j.1365-3040.1988.tb01774

CrossRef Full Text | Google Scholar

Sperry, J. S., Hacke, U. G., Oren, R., and Comstock, J. P. (2002). Water deficits and hydraulic limits to leaf water supply. Plant Cell Environ. 25, 251–263. doi: 10.1046/j.0016-8025.2001.00799.x

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Sperry, J. S., Meinzer, F. C., and Mcculloh, K. A. (2008). Safety and efficiency conflicts in hydraulic architecture: scaling from tissues to trees. Plant Cell Environ. 31, 632–645. doi: 10.1111/j.1365-3040.2007.01765

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Sperry, J. S., Nichols, K. L., Sullivan, J. E. M., and Eastlack, S. E. (1994). Xylem embolism in ring-porous, diffuse-porous, and coniferous trees of northern Utah and interior Alaska. Ecology 75, 1736–1752. doi: 10.2307/1939633

CrossRef Full Text | Google Scholar

Sperry, J. S., Smith, D. D., Savage, V. M., Enquist, B. J., McCulloh, K. A., Reich, P. B.,et al. (2012). A species-level model for metabolic scaling in trees I. Exploring boundaries to scaling space within and across species. Funct. Ecol. 26, 1054–1065. doi: 10.1111/j.1365-2435.2012.02022

CrossRef Full Text | Google Scholar

Spicer, R., and Gartner, B. L. (2001). The effects of cambial age and position within the stem on specific conductivity in Douglas-fir (Pseudotsuga menziesii) sapwood. Trees 15, 222–229. doi: 10.1007/s004680100093

CrossRef Full Text | Google Scholar

Sterck, F. J., Martinez-Vilalta, J., Mencuccini, M., Cochard, H., Gerrits, P., Zweifel, R.,et al. (2012). Understanding trait interactions and their impacts on growth in Scots pine branches across Europe. Funct. Ecol. 26, 541–549. doi: 10.1111/j.1365-2435.2012.01963

CrossRef Full Text | Google Scholar

Steudle, E. (2000). Water uptake by roots: effects of water deficit. J. Exp. Bot. 51, 1531–1542. doi: 10.1093/jexbot/51.350.1531

CrossRef Full Text | Google Scholar

Swenson, N. G., and Enquist, B. J. (2007). Ecological and evolutionary determinants of a key plant functional trait: wood density and its community-wide variation across latitude and elevation. Am. J. Bot. 94, 451–459. doi: 10.3732/ajb.94.3.451

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Thomas, D. S., Montagu, K. D., and Conroy, J. P. (2007). Temperature effects on wood anatomy, wood density, photosynthesis and biomass partitioning of Eucalyptus grandis seedlings. Tree Physiol. 27, 251–260. doi: 10.1093/treephys/27.2.251

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Tscharntke, T., Clough, Y., Bhagwat, S. A., Buchori, D., Faust, H., Hertel, D.,et al. (2011). Multifunctional shade-tree management in tropical agroforestry landscapes-a review. J. Appl. Ecol. 48, 619–629. doi: 10.1111/j.1365-2664.2010.01939

CrossRef Full Text | Google Scholar

Tyree, M. T. (2003). Hydraulic limits on tree performance: transpiration, carbon gain and growth of trees. Trees 17, 95–100. doi: 10.1007/s00468-002-0227

CrossRef Full Text | Google Scholar

Tyree, M. T., and Ewers, F. W. (1991). The hydraulic architecture of trees and other woody plants. New Phytol. 119, 345–360. doi: 10.1111/j.1469-8137.1991.tb00035

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Tyree, M. T., and Zimmermann, M. H. (2002). Xylem Structure and The Ascent of Sap. Berlin: Springer. doi: 10.1007/978-3-662-04931-0

CrossRef Full Text | Google Scholar

West, G. B., Brown, J. H., and Enquist, B. J. (1999). A general model for the structure and allometry of plant vascular systems. Nature 400, 664–667. doi: 10.1038/23251

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Wheeler, J. K., Sperry, J. S., Hacke, U. G., and Hoang, N. (2005). Inter-vessel pitting and cavitation in woody Rosaceae and other vesselled plants: a basis for a safety versus efficiency trade-off in xylem transport. Plant Cell Environ. 28, 800–812. doi: 10.1111/j.1365-3040.2005.01330

CrossRef Full Text | Google Scholar

White, F. M. (1991). Viscous Fluid Flow, 2nd Edn. New York: McGraw-Hill, Inc.

Google Scholar

Würth, M. K. R., Peláez-Riedl, S., Wright, S. J., and Körner, C. (2005). Non-structural carbohydrate pools in a tropical forest. Oecologia 143, 11–24. doi: 10.1007/s00442-004-1773-2

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Zanne, A. E., Westoby, M., Falster, D. S., Ackerly, D. D., Loarie, S. R., Arnold, S. E. J.,et al. (2010). Angiosperm wood structure: global patterns in vessel anatomy and their relation to wood density and potential conductivity. Am. J. Bot. 97, 207–215. doi: 10.3732/ajb.0900178

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Zhang, J.-L., and Cao, K.-F. (2009). Stem hydraulics mediates leaf water status, carbon gain, nutrient use efficiencies and plant growth rates across dipterocarp species. Funct. Ecol. 23, 658–667. doi: 10.1111/j.1365-2435.2009.01552

CrossRef Full Text | Google Scholar

Zufferey, V., Cochard, H., Ameglio, T., Spring, J.-L., and Viret, O. (2011). Diurnal cycles of embolism formation and repair in petioles of grapevine (Vitis vinifera cv. Chasselas). J. Exp. Bot. 62, 3885–3894. doi: 10.1093/jxb/err081

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Keywords: shade tree, hydraulic conductivity, wood density, aboveground productivity, foliar nitrogen, perhumid climate, vessel diameter

Citation: Kotowska MM, Hertel D, Abou Rajab Y, Barus H and Schuldt B (2015) Patterns in hydraulic architecture from roots to branches in six tropical tree species from cacao agroforestry and their relation to wood density and stem growth. Front. Plant Sci. 6:191. doi: 10.3389/fpls.2015.00191

Received: 24 November 2014; Accepted: 09 March 2015;
Published online: 31 March 2015

Edited by:

Boris Rewald, University of Natural Resources and Life Sciences, Austria

Reviewed by:

Heidi Renninger, Mississippi State University, USA
Claudio Lovisolo, University of Turin, Italy

Copyright © 2015 Kotowska, Hertel, Rajab, Barus and Schuldt. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Martyna M. Kotowska, Plant Ecology and Ecosystems Research, Albrecht von Haller Institute for Plant Sciences, University of Göttingen, Untere Karspüle 2, 37073 Göttingen, Germany mkotows@gwdg.de

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.