Disentangling the E ﬀ ects of Genotype and Environment on Growth and Wood Features of Balfourodendron riedelianum Trees by Common Garden Experiments in Brazil

: Intraspeciﬁc studies with populations replicated in di ﬀ erent sites allow the e ﬀ ects of genotype and environment on wood features and plant growth to be distinguished. Based on climate change predictions, this distinction is important for establishing future patterns in the distribution of tree species. We quantiﬁed the e ﬀ ects of genotype and environment on wood features and growth of 30-year-old Balfourodendron riedelianum trees. We used three provenances planted in two common garden experiments with di ﬀ erence in precipitation and temperature. We applied linear models to estimate the variability in wood and growth features and transfer functions to evaluate the responses of these features to temperature, precipitation, and the standardized precipitation evapotranspiration index (SPEI). Our results showed that genotype had an e ﬀ ect on vessels and rays, where narrower vessels with thinner walls and larger intervessel pits, and shorter, narrower and more numerous rays were observed in provenances from drier sites. We also observed the e ﬀ ect of the environment on wood features and growth. Trees growing in the wetter site were taller and larger, and they had wider vessels with thicker walls and lower ray density. Transfer functions indicated that an increase in temperature results in larger vessels with thicker walls, taller and denser rays, shorter and narrower ﬁbers with thinner walls, and lower wood density. From a functional perspective, these features make trees growing in warmer environments more prone to drought-induced embolisms and more vulnerable to mechanical damage and pathogen attacks. Tree growth varied with precipitation and SPEI, being negatively a ﬀ ected in the drier site. Overall, we demonstrated that both genotype and environment a ﬀ ect wood features, while tree growth is mainly inﬂuenced by the environment. Plastic responses in hydraulic characteristics could represent important functional traits to mitigate the consequences of ongoing climate change on the growth and survival of the species within its natural range.


Introduction
The effect of environmental variability on the wood features in long-lived plants has been studied for almost a century. These studies allowed establishment of the main ecological trends in wood anatomical features of woody plant populations grown in sites with different water availability, showing the effects of local conditions (e.g., precipitation, temperature, and soil physical properties) on their anatomical structure. However, besides the environmental conditions of the growing site, the wood features still remain under strong genetic control [1], making it a challenge to determine which anatomical features are determined by genotype, induced by environmental variability, or resulting from an interaction of these two factors [2].
Common garden experiments with different seed sources replicated in environments with variable water availability has proven to be a useful tool in disentangling the genetic from environmental effects on wood features in tree species with wide distributions [3,4]. For example, common garden experiments performed by Fisher et al. [5] and Schreiber et al. [6] showed that the place of origin of plants is an important selection factor for vessel diameter variability. In general, populations coming from places prone to cold or drought tend to have narrower vessels as a result of high selective pressures under water stress environments [1,7]. Instead, the responses of hydraulic traits to changes in environmental plant growth conditions varied among populations. In some of them, vessel diameters were sensitive to increased water availability, producing slightly wider vessels and consequently increasing water transport efficiency, with positive effects on plant growth and survival. However, in some populations, the genotype may buffer the effects of environmental variability on wood formation, demonstrating that the interaction between genotype and environment in determining their anatomical features is more complex than expected. The effects of the genotype and the environment are also observed in plant growth. In general, trees from drier sites are generally smaller than trees from wetter sites, according to a well-known positive effect of precipitation on growth [8,9].
Another important factor that should be considered in relation to the responses of wood to environmental variability is the involvement of cells in simultaneous multiple functions [10,11]. The numerous possibilities of arrangement, as well as the diverse interactions between the different cell types (water-conducting cells, fibers, and axial parenchyma and rays), have allowed long-lived plants to develop different strategies that ensure their growth and survival in periods of water stress. For example, in drier environments, to minimize hydraulic failures and avoid vessel collapse due to higher water flow stresses, narrower vessels are formed, with thicker cell walls [12,13], denser [12,13], more clustered [13,14], with smaller intervessel pits [15]. The narrowing of vessel diameter, associated with the thickening of cell walls, may confer greater safety to water transport, increasing the mechanical resistance of vessel walls to implosion [16]. Higher vessel density increases the possibility of water flow pathways in the case of embolism [17] at the same time as it increases the water conductive area [18]. Vessels in larger groupings provide greater connectivity between them, offering alternative pathways to water flow in cases of embolisms, which could keep the water column functional without substantial reduction in water transport [18].
Adaptation of the wood experiencing changes in water availability also involves adjustments between vessels, parenchyma cells (axial parenchyma and rays) and fibers. In a global approach, Morris et al. [19] showed that there is a positive correlation between vessel diameter and the amount of axial parenchyma in warmer sites. In these conditions, the interaction between parenchymal cells and vessels may act as an ancillary mechanism for greater water transport safety during periods of water stress, increasing the possibility of water refilling in embolized vessels [20] and water storage capacity in the stem [21]. Another possibility of wood construction to minimize the effects of water stress is through the interaction between vessels and fibers. Wood composed of a fiber matrix with thick walls in dry environments could provide greater mechanical support to the vessels and thus reduce the risks of vessel implosion in case of higher negative pressure [10,16]. The formation of a dense matrix of thick-walled fibers can also lead to an increase in the mechanical strength of the stem with a positive effect on wood density, while decreasing water conduction and storage capacity in the stem, reducing plant growth [22]. In addition, a dense fiber network could compensate for weaker areas of wood composed of large or denser vessels [16]. Such evidence shows that wood anatomical variations are related to the integrated changes between its different cell types to adjust to environmental variability [23], which could minimize, through different pathways, the effects of environmental variation on hydraulic traits. Thus, it is clear that detailed analyses of wood and its responses to drought events are essential for a detailed assessment of plants potential to withstand ongoing climate change. However, so far, little is known about the intraspecific variation and interaction between the cellular types of wood, as well as the effect of genotype and environment on these cells, in tree species with wide geographic distribution.
Under the current scenarios of changes in the seasonality of rainfall and warming [24], knowing which wood features are determined by genotype and which are more plastic will be essential to indicate the impacts of climate change on xylem formation. Furthermore, knowing the plastic responses of the hydraulic system and its effects on plant growth will contribute to improve the strategies of assisted migration of tree species from and towards forest sites worldwide, especially in neotropical forests. We therefore investigated the intraspecific variation in wood features and tree growth, estimated by stem height and diameter, in 30-year-old trees from three populations of Balfourodendron riedelianum (Engl.) Engl. (Rutaceae) grown in two common garden experiments as provenance tests. B. riedelianum is a long-lived and slow growing deciduous species reaching up to 32 m in height, with stems reaching 80 cm in diameter. Native of semideciduous seasonal forests, this species has a wide geographical distribution in Neotropical seasonally dry forests of Brazil, Argentina, and Paraguay [25]. In this study, we evaluated the effect of climate at the origin sites on wood features and tree growth through transfer functions. The following hypotheses were tested: (i) trees originating from drier sites exhibit lower growth and produce wood with higher hydraulic safety; (ii) trees planted in trial sites with greater water availability exhibit higher stem growth and wood with large vessels regardless of their origin.

Study Species, Common Gardens, and Sampling
This study was carried out in two common gardens with 30-year-old B. riedelianum trees located in Luís Antônio (LA; 21 •  Both sites are at an altitude of 500 m and have a dry season from May to September (rainfall ≤ 50 mm), as well as a hot and rainy season from November to March. The annual precipitation, mean annual temperature, summer mean temperature and winter mean temperature are described in Table 1. The soils differ between sites: dark red latosol with clay texture in LA [26] and sandy yellow latosol in PE [27].
The common gardens originated from seeds collected in three natural populations of B. riedelianum, Alvorada do Sul (AS), in the State of Paraná, Bauru (BA) and Gália (GA), both in the State of São Paulo ( Table 1). The seeds in AS were collected in the semideciduous seasonal forest. The seeds in BA were collected at the ecotone between the semideciduous seasonal forest and Cerrado. The seeds in GA were collected in the semideciduous seasonal forest. Both common gardens were established in 1985 using seedlings of B. riedelianum from open-pollinated seeds collected from 19 maternal trees per populations (AS, BA, and GA), and involving provenance and progeny trials [30]. Provenance refers to the place of origin of populations and progeny denotes individuals that share maternal genetic material (half-siblings). The common gardens were arranged as randomized complete blocks designed with six replications per provenance. Each block was divided into three plots, with one provenance per plot, and progeny allocated in sub-plots. In each row of the sub-plots, five seedlings of the same progeny were planted at a spacing of 3 m × 3 m, with two external border rows of the same species. We measured stem height and diameter at breast height, and collected wood samples from the main stem at 1.30 m above the ground from 36 trees for each common garden (12 trees per provenance). From each provenance, we selected three random progenies (AS, BA, and GA) and four random trees per progeny within the sub-plots. We collected trees of the same progeny in both trial sites.

Wood Anatomical Procedures
We collected one 10 cm disc from the main stem at 1.30 m above the ground per sampled tree. From each disc, we removed wood samples with approximately 3 cm 3 containing the sapwood nearest to the bark for wood-anatomical, scanning electron microscopy, and wood density analysis.
For anatomical analysis, we fixed wood samples in FAA 50 (formaldehyde, acetic acid, ethanol 50% 1:1:18) for 48 h and then transferred them to 70% ethanol [31] for storage. We performed histological sections using a sliding microtome on the radial longitudinal, tangential longitudinal, and cross directions. The thickness of the histological sections varied from 12 to 15 µm. Sections were stained with 1% aqueous of safranin [32] and 1% aqueous of astra blue [33] and mounted on permanent slides with Entellan ® synthetic resin. For maceration, we followed the Franklin method [34] modified by Kraus and Arduin [35]. We put small portions of the wood samples in a solution of glacial acetic acid and hydrogen peroxide (1:1 v:v) at 60 • C for 24 h, washed in water, stained with aqueous 1% safranin, and mounted semi-permanent slides in a glycerin and water solution (1:1 v:v). We analyzed the wood anatomy slides under a light microscope equipped with an AxioCAM MRC and Axiovision software. We used the recommendations of the International Association of Wood Anatomists Committee [36] for wood anatomical description.

Wood Data Collection
We measured wood anatomical features, including vessel density, vessel wall thickness, vessel element length; height and width of rays and numbers of rays per linear mm; and length, diameter, and lumen thickness of fibers. We performed 30 measurements for each wood anatomical features per sampled tree. We used images with 300 dpi resolution. To determine the vessel density, we counted the number of vessels in 30 areas of 1 mm 2 in transverse section per tree sampled. We calculated the fiber wall thickness (Fwt) by the following formula: where (Fd) = total fiber diameter, and (d) = fiber lumen diameter. We estimated vessel, fiber and parenchyma (rays + axial parenchyma) fractions (%) in the wood. We used 15 areas of 1 mm 2 in transverse sections images per tree. To assess the fraction (%) of area occupied by the parenchyma, we calculated by subtracting the area of fibers and vessels from total measured wood area. We performed all wood anatomical linear and area measurements using ImageJ 1.6.0 (https://imagej.nih.gov/ij/download.html; National Institutes of Health, Bethesda, MD, USA).
Since vessels are not exactly circular as seen in the transverse section, we calculated vessel diameter (Dv) of each vessel based on the vessel lumen area [37]. We calculated the vessel diameter as: where A is the vessel lumen area. The vessel lumen area was measured from 50 vessels per tree. We calculated the potential hydraulic conductivity (Kp) following the Hagen-Poiseuille law [22].
where Kp is the potential hydraulic conductivity (in Kg m −1 MPa −1 s −1 ), ρw is the density of water at 20 • C (998.2 kg m −3 at 20 • C), n is the viscosity of water at 20 • C (1.002 10 −3 Pa s at 20 • C), VD is the vessel density (n • vessels mm −2 ), and Dh is the hydraulically weighted vessel diameter (in m). The hydraulically weighted vessel diameter was calculated using the following formula: where d is the vessel diameter (in µm), and N is the number of vessels. The Dh represents the mean diameter that all vessels in stem would have in order to correspond to the same total conductivity for the same numbers of vessels [37,38].
In addition, we calculated the vessel grouping index as the total number of vessels divided by the total number of vessel groupings [37,39]. We selected 25 vessel groups and counted the total number of vessels in these groups and then divided the total number of vessel groups by 25. Both solitary vessels and grouped vessels were counted as vessel groups.

Scanning Electron Microscopy
We measured the intervessel pit area and intervessel pit aperture area on images taken with a scanning electron microscope. We selected three trees of each population (AS, BA, and GA) in each common garden (LA and PE). Wood samples fixed in FAA 50 with 2 mm lengths were split in a tangential plane, dehydrated in ethanol series from 50 to 100%, and then air dried for 24 h. The split wood samples were fixed to stubs with an electron-conductive carbon sticker, gold-coated with a sputter coater (Baltec SCD 050) for 3 min. We observed the intervessel pits with a scanning electron microscope at an accelerating voltage of 15 kV. The intervessel pit area and intervessel pit aperture were measured from 50 vessels per tree using ImageJ 1.6.0 (https://imagej.nih.gov/ij/download.html; National Institutes of Health, Bethesda, MD, USA).

Wood Density
We used the water displacement method to assess wood density (in g cm −3 ), corresponding to the wood dry mass over wood fresh volume [40]. Fresh volume of wood was estimated by immersion in water. We considered the weight of the displaced water equal to the volume of the sample (1 g = 1 cm 3 ). The wood samples were then dried at 105 • C for 72 h and weighed to measure dry mass.

Statistics and Transfer Functions
Multivariate analysis of variance (MANOVA) was performed based on the measurements collected on the 72 study trees, involving tree height, and diameter and wood features. The three provenance origins and two common gardens, as well as their interaction, were included as independent factors. Post-hoc tests based on contrasts were used to verify differences within each factor. We performed univariate analyses using mixed models by testing for differences between provenance origin, common gardens, and progenies using tree as a random factor. The factor progeny was nested in the provenance origin.
We extracted monthly temperatures and precipitations at the provenance origins and trial sites for a period covering 30 years (1970-2000) from the climatic datasets in WorldClim 2 [41]. Standardized precipitation evapotranspiration index (SPEI) [42] was calculated using SPEI package [43] in R software, version 1.2.5033 (R Foundation for Statistical Computing, Vienna, Austria) [44]. We performed Principal Component Analysis (PCA) to describe the climatic distance between sites based on the 19 bioclimatic parameters proposed by O'Donnell and Ignizio [45]. We identified the bioclimatic parameters with the stronger impact on the principal components based on Pearson's correlation coefficients. Response to climate shift was calculated using linear regressions as a function of the climate transfer difference, i.e., difference in temperature, precipitation, and SPEI between the provenance origin and trial sites. Unless otherwise specified, the statistics were performed in SAS 9.4 (SAS Institute Inc., Cary, NC, USA).  Stem height and diameter differed only between trial sites, where trees in LA trial site were taller and larger than those in PE trial site (Figure 2). In relation to trial sites, differences were observed in the following wood features: potential hydraulic conductivity, vessel diameter, vessel wall thickness, vessel element length, vessel grouping index, intervessel pit area, intervessel pit aperture area, fiber lumen diameter, fiber wall thickness, ray density, and parenchyma fraction. In general, trees in LA trial site had sparser and wider vessels (Figure 3), thicker walls, longer vessel elements, larger intervessel pits, fibers with thicker wall and larger lumen, lower ray density, lower parenchyma fraction, and higher potential hydraulic conductivity than in PE trial site ( Figure 4; Table 2; Supplementary Table S1). Wood density, vessel density, fiber length and diameter, ray height and width, vessel fraction, and fiber fraction did not differ significantly between the trial sites. Differences between provenances were observed for potential hydraulic conductivity, vessel diameter, vessel wall thickness, vessel density, vessel grouping index, intervessel pit area, intervessel pit aperture area, ray height, ray width, ray density, vessel fraction, fiber fraction, and parenchyma fraction.    In general, the trees from GA had larger vessels ( Figure 3) with thicker walls, less dense and clustered, and with larger intervessel pits, taller and wider and lower density rays, smaller vessel and parenchyma fraction, and larger fiber fraction, plus higher potential hydraulic conductivity than the trees from AS and BA (Figure 4; Table 2; Supplementary Table S1).

Variation in Wood Features and Growth
The analyses showed that vessel diameter, vessel wall thickness, vessel element length, and ray width differed between progenies. There was an interaction trial site × provenance in the potential hydraulic conductivity, vessel elements length, intervessel pit area, fiber lumen diameter, and fiber wall thickness, as well as a trial site × progeny interaction in vessel diameter and vessel wall thickness ( Table 2).

Climate Parameters at the Provenance Origins and Trial Sites
Bioclimatic parameters of provenance origins and trial sites showed a clear cluster including GA and BA provenances ( Figure 5). AS provenance is located at the periphery of the plot, indicating its distinct position and different climatic conditions in relation to the two other provenances. The trial sites are located in different positions in the plot. PE trial site is included in the same cluster as BA and GA provenances, indicating very similar climatic parameters. The first principal component of the PCA (PC1) explained 80.28% of the variation and was mainly connected to thermal conditions; more specifically it was positively correlated with mean temperature of the driest quarter (bio9) and mean temperature of the coldest quarter (bio 11). The second principal component of the PCA (PC2) explained 13.52% of the variation, was connected to precipitation, and, more specifically, it was correlated positively with annual precipitation (bio 12).

Transfer Functions
Temperature, precipitation, and SPEI of the trial sites had an effect on stem height and diameter at breast height and on most wood features in B. riedelianum trees (Table 3). Differences in stem height and diameter at breast height were related to variation in precipitation and SPEI with the increase in stem height and diameter when trees were transferred to a wetter location. However, stem height and diameter decrease when the provenances were transferred to the trial site with more drought (Figure 6).  The temperature of the trial sites had greater influence on wood features. Trees transferred to a warmer trial site had lower wood density, larger vessel diameter, thicker vessel wall, lower vessel grouping index, larger fiber diameter, lower fiber wall thickness, and higher height and ray density than trees transferred to a cooler trial site (Table 3; Figure 6). Precipitation and SPEI also had an effect on vessel diameter, potential hydraulic conductivity and length of vessel elements. The trees transferred to a wetter trial site and with less SPEI had a larger vessel diameter, higher potential hydraulic conductivity, longer vessel elements, lower ray density, and lower parenchyma fraction than trees transferred to a drier trial site. The fiber lumen diameter was affected only by the SPEI, with trees transferred to a drier trial site having a larger fiber lumen diameter. Moreover, there was an interaction between precipitation and temperature in the potential hydraulic conductivity, ray density, fiber lumen diameter, and fiber wall thickness ( Figure 6).

Discussion
This study investigated the intraspecific variation in wood features and tree growth of B. riedelianum trees originating from three populations and grown in two common garden experiments in São Paulo State, Brazil. The results show that tree growth, determined by stem height and diameter, was mostly influenced by the environment at the trial sites, specifically precipitation and SPEI, being promoted by wet conditions. Wood features vary among populations and are strongly influenced by the temperature of trial sites.

Genetic Difference between Populations
We raised the hypothesis that intraspecific variation in the wood features of B. riedelianum trees would be related to the climatic of populations origin site. We expected that populations from drier and warmer sites, with lower precipitation, would have wood characterized by higher water transport safety than trees from cooler sites with higher precipitation. Indeed, trees from Bauru, the site with lower annual precipitation and higher temperature, had narrower vessels and thinner walls, more clustered and denser, and occupying a larger fraction of wood than trees from other populations. Together, these anatomical features provide greater safety for water transport under higher negative tensions in drier environments [16,38]. However, unlike our expectation, trees from Alvorada do Sul, the site with higher annual precipitation than other provenances showed wood with water safety similar to that in trees from Bauru. We observed wood anatomical features with greater efficiency in water transport in trees from Gália, the site with intermediate precipitation between provenances.
We found a lower potential hydraulic conductivity, narrow vessels, with thin walls, in higher density and more clustered, with smaller intervessel pits and narrow openings, rays in higher density and higher parenchyma fraction (axial parenchyma + rays), features that together indicate greater water transport safety in trees from Alvorada do Sul. This region is characterized by high rainfall, fertile soils, low winter temperatures, and being prone to frost [28]. In general, narrow vessels are advantageous in sites with low temperatures to avoid freeze-induced cavitation, as well as to minimize the impacts of drought-induced cavitation [46]. According to Fisher et al. [5] and Schreiber et al. [6], the selection of narrow vessels may be an adaptation in populations that experience freezing periods that minimize the effects caused by frost induced cavitation, thus maintaining the xylem functional for the next growing season. Our findings, therefore, suggest that the formation of wood with high hydraulic safety in trees from Alvorada do Sul could reflect the low temperatures experienced at the origin site. However, experimental studies are necessary to confirm this statement.

Phenotypic Plasticity of Tree Growth and Wood Formation
In general, our trees responded to the increase in water availability through the formation of wood with a higher water transport efficiency. Regardless of origin, trees grown in a wetter environment (Luís Antônio) had wider vessels, with thicker walls, less clustered, with larger intervessel pits and wider openings, longer vessel elements, fibers with narrower lumen and thicker walls, less frequent rays and smaller fraction of parenchyma (axial parenchyma + rays). The increase in vessel diameter also had a positive effect on potential hydraulic conductivity. Our results showed that the vessels are the cells with the highest degree of plasticity in B. riedelianum, with variations in diameter, wall thickness, grouping, and size of intervessel pits. Vessels are responsible for long-distance water conduction in plants, performing a central role in hydraulic efficiency and safety [16]. Variations in vessel dimensions could represent a key trait for the adaptation or acclimation to changes in environmental conditions, such as increasing of temperatures and drought events, thus ensuring growth and survival of this species within its natural distribution range.
The plastic responses observed in wood of mature B. riedelianum trees in the wetter trial site reflect greater efficiency in conducting water and mechanical support on vessel walls. In this trial site, the formation of larger vessels provides low resistance to water flow and increases water transport ability [16,18]. The thickening of vessel walls makes larger vessels less prone to mechanical damage, such as micro-fractures and implosion of their walls [16]. Micro-fractures in the vessel wall may trigger the heterogeneous nucleation of gas particles between the cell walls or intercellular spaces, initiating the embolism process [38,47]. Large, thick-walled vessels also have larger intervessel pits with wider openings. These latter are potentially more efficient in conducting water but decrease the mechanical reinforcement of vessel cell walls [11,48]. In addition, intervessel pits with large openings could decrease support to the pit membrane during stretching [49]. Thus, the formation of fibers with narrower lumen and thicker walls found in B. riedelianum trees grown in the wetter trial site could compensate for potential cell wall weakening. Therefore, wide thicker-walled vessels immersed in a denser fiber matrix may confer greater mechanical strength to support larger intervessel pits and with wider openings in the wood.
We found plasticity in the size of vessel elements. In vessels, the length is determined by the size of initial fusiforms of the vascular cambium [50], and the diameter is linked to cell expansion. Variation in water availability to vascular cambium has an effect on the elongation of the initial fusiforms [51] and cellular expansion [52]. Therefore, the higher water availability in Luís Antônio may have favored the formation of longer and wider vessels. Similarly, the lower precipitation in Pederneiras may have been a limiting factor to elongation of the initial fusiform of the vascular cambium and cell expansion, thus resulting in the formation of shorter vessel elements and narrower vessels.
When testing the responses of the wood of B. reidelianum to variations in temperature, precipitation and SPEI among trial sites using transfer functions, we found that temperature had a greater effect on the anatomical features of vessels, fibers and rays. In general, the wood showed a greater capacity to water transport and vulnerability to mechanical damage, with larger vessels, less clustered and with thicker walls and longer fibers, wider and thinner walls in warmer environments. In the context of ongoing global changes, plants will experience warmer temperatures, often associated with less precipitation. The formation of wider vessels makes the wood more prone to drought-induced embolism, leading to hydraulic failure and even tree death [53]. Our results also indicate that taller rays and higher density occur in warmer environments. From the functional perspective, taller rays and higher density associated with larger vessels would be an important requirement to maintain hydraulic conductance in the stem, as well as hydraulic recovery through vessel refilling mechanisms [20,54].
Our results also indicate that wood density responds to an increase in temperature, with less dense wood occurring in warmer environments. The lower density may result from the combination of wider elements, vessels, and fibers, with thinner walls. Less dense wood leads to a decrease in the mechanical strength of the stem [55], making it more fragile and prone to breakage caused by animals or falling debris [56]. The decrease in wood density makes plants less resistant to damage caused by herbivores and pathogen attack, which can decrease their growth, life-span, and biomass allocation [57]. Therefore, our results suggest that an increase in temperature may stimulate B. riedelianum to form wood more prone to embolism, with lower mechanical resistance and higher susceptibility to pathogen attack.
The growth in stem height and diameter was sensitive to changes in the environmental conditions of the trial sites. Regardless of provenance, all trees responded similarly in both common garden experiments, with increased growth in stem height and diameter in the wetter trial site. In addition, the result of the transfer function on stem height and diameter growth demonstrates that tree growth benefits from an increase in precipitation and responds negatively to higher SPEI. Considering the climate predictions indicating 25% increases in precipitation across the natural range of B. riedelianum, a region that mainly covers southern and south-eastern Brazil [24], we may expect a positive impact on tree growth in a future scenario. However, the possible benefits of an increase in precipitation should be carefully examined. Firstly, possible changes in rainfall intensity, duration, and distribution throughout the year should be taken into account. Secondly, south and south-eastern Brazil is already exploited for crop production, and potential new areas are awaiting future agricultural expansion and development in the next decades, which will make the natural range of the studied species more prone to deforestation. At the same time, an important warming is predicted in South America [24], making it difficult to predict the growth responses of trees under complex climate change scenarios. Under these conditions, plastic responses in the growth of B. riedelianum could be essential to guarantee persistence of the species in the drier sites.

Conclusions
In this study, we investigated the intraspecific variation in wood features and tree growth of three populations of Balfourodendron riedelianum growing in two common garden experiments and evaluated the effect of climate at the origin sites through transfer functions. We observed that phenotypic variability in wood features results from a clear genetic differentiation between populations. We found that both genotype and environment have an effect on wood features, while variability in tree growth is mostly determined by the environment. Based on the wide variability in wood features and tree growth, B. riedelianum could be able to adjust growth and hydraulic traits to withstand the changes in temperature and precipitation predicted by the climatic models in the next decades, thus ensuring the survival of this species under warmer scenarios.
Supplementary Materials: The following are available online at http://www.mdpi.com/1999-4907/11/9/905/s1, Table S1: Height of stem, diameter at breast height of the stem, wood density, potential hydraulic conductivity and wood quantitative anatomical features of 30-year-old Balfourodendron riedelianum trees from three provenance (AS, BA and GA) grown in Luís Antônio and Pederneiras common garden experiments. AS = Alvorada do Sul; BA = Bauru; GA = Gália. The values are shown as mean ± standard deviation (minimum and maximum values).