Ecological and social factors influence interspecific pathogens occurrence among bees.

The interspecific transmission of pathogens can occur frequently in the environment. Among wild bees, the main spillover cases are caused by pathogens associated with Apis mellifera, whose colonies can act as reservoirs. Due to the limited availability of data in Italy, it is challenging to accurately assess the impact and implications of this phenomenon on the wild bee populations. In this study, a total of 3372 bees were sampled from 11 Italian regions within the BeeNet project, evaluating the prevalence and the abundance of the major honey bee pathogens (DWV, BQCV, ABPV, CBPV, KBV, Nosema ceranae, Ascosphaera apis, Crithidia mellificae, Lotmaria passim, Crithidia bombi). The 68.4% of samples were positive for at least one pathogen. DWV, BQCV, N. ceranae and CBPV showed the highest prevalence and abundance values, confirming them as the most prevalent pathogens spread in the environment. For these pathogens, Andrena, Bombus, Eucera and Seladonia showed the highest mean prevalence and abundance values. Generally, time trends showed a prevalence and abundance decrease from April to July. In order to predict the risk of infection among wild bees, statistical models were developed. A low influence of apiary density on pathogen occurrence was observed, while meteorological conditions and agricultural management showed a greater impact on pathogen persistence in the environment. Social and biological traits of wild bees also contributed to defining a higher risk of infection for bivoltine, communal, mining and oligolectic bees. Out of all the samples tested, 40.5% were co-infected with two or more pathogens. In some cases, individuals were simultaneously infected with up to five different pathogens. It is essential to increase knowledge about the transmission of pathogens among wild bees to understand dynamics, impact and effects on pollinator populations. Implementing concrete plans for the conservation of wild bee species is important to ensure the health of wild and human-managed bees within a One-Health perspective.

Anisotropic fluid flows in black phosphorus nanochannels.

With the development of advanced micro/nanoscale technologies, two-dimensional materials have emerged from laboratories and have been applied in practice. To investigate the mechanisms of solid-liquid interactions in potential applications, molecular dynamics simulations are employed to study the flow behavior of n-dodecane (C12) molecules confined in black phosphorus (BP) nanochannels. Under the same external conditions, a significant difference in the velocity profiles of fluid molecules is observed when flowing along the armchair and zigzag directions of the BP walls. The average velocity of C12 molecules flowing along the zigzag direction is 9-fold higher than that along the armchair direction. The friction factor at the interface between C12 molecules and BP nanochannels and the orientations of C12 molecules near the BP walls are analyzed to explain the differences in velocity profiles under various flow directions, external driving forces, and nanochannel widths. The result shows that most C12 molecules are oriented parallel to the flow direction along the zigzag direction, leading to a relatively smaller friction factor hence a higher average velocity. In contrast, along the armchair direction, most C12 molecules are oriented perpendicular to the flow direction, leading to a relatively larger friction factor and thus a lower average velocity. This work provides important insights into understanding the anisotropic liquid flows in nanochannels.

A Data-Driven Approach to Refine Predictions of Differentiated Thyroid Cancer Outcomes: A Prospective Multicenter Study.

CONTEXT: The risk stratification of patients with differentiated thyroid cancer (DTC) is crucial in clinical decision making. The most widely accepted method to assess risk of recurrent/persistent disease is described in the 2015 American Thyroid Association (ATA) guidelines. However, recent research has focused on the inclusion of novel features or questioned the relevance of currently included features. OBJECTIVE: To develop a comprehensive data-driven model to predict persistent/recurrent disease that can capture all available features and determine the weight of predictors. METHODS: In a prospective cohort study, using the Italian Thyroid Cancer Observatory (ITCO) database (NCT04031339), we selected consecutive cases with DTC and at least early follow-up data (n = 4773; median follow-up 26 months; interquartile range, 12-46 months) at 40 Italian clinical centers. A decision tree was built to assign a risk index to each patient. The model allowed us to investigate the impact of different variables in risk prediction. RESULTS: By ATA risk estimation, 2492 patients (52.2%) were classified as low, 1873 (39.2%) as intermediate, and 408 as high risk. The decision tree model outperformed the ATA risk stratification system: the sensitivity of high-risk classification for structural disease increased from 37% to 49%, and the negative predictive value for low-risk patients increased by 3%. Feature importance was estimated. Several variables not included in the ATA system significantly impacted the prediction of disease persistence/recurrence: age, body mass index, tumor size, sex, family history of thyroid cancer, surgical approach, presurgical cytology, and circumstances of the diagnosis. CONCLUSION: Current risk stratification systems may be complemented by the inclusion of other variables in order to improve the prediction of treatment response. A complete dataset allows for more precise patient clustering.

Online revenue maximization for server pricing.

Efficient and truthful mechanisms to price resources on servers/machines have been the subject of much work in recent years due to the importance of the cloud market. This paper considers revenue maximization in the online stochastic setting with non-preemptive jobs and a unit capacity server. One agent/job arrives at every time step, with parameters drawn from the underlying distribution. We design a posted-price mechanism which can be efficiently computed and is revenue-optimal in expectation and in retrospect, up to additive error. The prices are posted prior to learning the agent's type, and the computed pricing scheme is deterministic, depending only on the length of the allotted time interval and on the earliest time the server is available. We also prove that the proposed pricing strategy is robust to imprecise knowledge of the job distribution and that a distribution learned from polynomially many samples is sufficient to obtain a near-optimal truthful pricing strategy.

Assessment with clinical data of a coupled bio-hemodynamics numerical model to predict leukocyte adhesion in coronary arteries.

Numerical simulations of coupled hemodynamics and leukocyte transport and adhesion inside coronary arteries have been performed. Realistic artery geometries have been obtained for a set of four patients from intravascular ultrasound and angiography images. The numerical model computes unsteady three-dimensional blood hemodynamics and leukocyte concentration in the blood. Wall-shear stress dependent leukocyte adhesion is also computed through agent-based modeling rules, fully coupled to the hemodynamics and leukocyte transport. Numerical results have a good correlation with clinical data. Regions where high adhesion is predicted by the simulations coincide to a good approximation with artery segments presenting plaque increase, as documented by clinical data from baseline and six-month follow-up exam of the same artery. In addition, it is observed that the artery geometry and, in particular, the tortuosity of the centerline are a primary factor in determining the spatial distribution of wall-shear stress, and of the resulting leukocyte adhesion patterns. Although further work is required to overcome the limitations of the present model and ultimately quantify plaque growth in the simulations, these results are encouraging towards establishing a predictive methodology for atherosclerosis progress.

Individual reproductive success in Norway spruce natural populations depends on growth rate, age and sensitivity to temperature.

Quantifying the individual reproductive success and understanding its determinants is a central issue in evolutionary research for the major consequences that the transmission of genetic variation from parents to offspring has on the adaptive potential of populations. Here, we propose to distil the myriad of information embedded in tree-ring time series into a set of tree-ring-based phenotypic traits to be investigated as potential drivers of reproductive success in forest trees. By using a cross-disciplinary approach that combines parentage analysis and a thorough dendrophenotypic characterisation of putative parents, we assessed sex-specific relationships between such dendrophenotypic traits (i.e., age, growth rate and parameters describing sensitivity to climate and to extreme climatic events) and reproductive success in Norway spruce. We applied a full probability method for reconstructing parent-offspring relationships between 604 seedlings and 518 adult trees sampled within five populations from southern and central Europe. We found that individual female and male reproductive success was positively associated with tree growth rate and age. Female reproductive success was also positively influenced by the correlation between growth and the mean temperature of the previous vegetative season. Overall, our results showed that Norway spruce individuals with the highest fitness are those who are able to keep high-growth rates despite potential growth limitations caused by reproductive costs and climatic limiting conditions. Identifying such functional links between the individual ecophysiological behaviour and its evolutionary gain would increase our understanding on how natural selection shapes the genetic composition of forest tree populations over time.

New Tools for Conservation Biological Control: Testing Ant-Attracting Artificial Nectaries to Employ Ants as Plant Defenders.

Knowledge of the role of ants in many agroecosystems is relatively scarce, and in temperate regions the possibility to exploit ants as biocontrol agents for crop protection is still largely unexplored. Drawing inspiration from mutualistic ant-plant relationships mediated by extrafloral nectaries (EFNs), we tested the use of artificial nectaries (ANs) in order to increase ant activity on pear trees and to evaluate the effects on the arthropods, plant health and fruit production. While EFNs secrete a complex solution mainly composed of sugars and amino acids, ANs were filled with water and sucrose only. The results suggest that ANs can be used as manipulative instruments to increase ant activity over long periods of time. High ant activity was significantly linked to lower incidence of the pathogen fungus Venturia pyrina (pear scab) on pear leaves, and of the presence of Cydia pomonella (codling moth) caterpillars on pear fruit production. These results further encourage exploring underrated possibilities in the development of new tools for conservation biological control (CBC).

Designing Cost-Sharing Methods for Bayesian Games.

We study the design of cost-sharing protocols for two fundamental resource allocation problems, the Set Cover and the Steiner Tree Problem, under environments of incomplete information (Bayesian model). Our objective is to design protocols where the worst-case Bayesian Nash equilibria have low cost, i.e. the Bayesian Price of Anarchy (PoA) is minimized. Although budget balance is a very natural requirement, it puts considerable restrictions on the design space, resulting in high PoA. We propose an alternative, relaxed requirement called budget balance in the equilibrium (BBiE). We show an interesting connection between algorithms for Oblivious Stochastic optimization problems and cost-sharing design with low PoA. We exploit this connection for both problems and we enforce approximate solutions of the stochastic problem, as Bayesian Nash equilibria, with the same guarantees on the PoA. More interestingly, we show how to obtain the same bounds on the PoA, by using anonymous posted prices which are desirable because they are easy to implement and, as we show, induce dominant strategies for the players.

Disentangling the effects of spatial proximity and genetic similarity on individual growth performances in Norway spruce natural populations.

Cambial growth is a phenotypic trait influenced by various physiological processes, numerous biotic and abiotic drivers, as well as by the genetic background. By archiving the outcome of such complex interplay, tree-rings are an exceptional resource for addressing individual long-term growth responses to changing environments and climate. Disentangling the effects of the different drivers of tree growth, however, remains challenging because of the lack of multidisciplinary data. Here, we combine individual dendrochronological, genetic and spatial data to assess the relative importance of genetic similarity and spatial proximity on Norway spruce (Picea abies (L.) Karst.) growth performances. We intensively sampled five plots from two populations in southern and central Europe, characterizing a total of 482 trees. A two-step analytical framework was developed. First, the effects of climate and tree age on tree-ring width (TRW) were estimated for each individual using a random slope linear mixed-effects model. Individual parameters were then tested against genetic and spatial variables by Mantel tests, partial redundancy analyses and variance partitioning. Our modelling approach successfully captured a large fraction of variance in TRW (conditional R2 values up to 0.94) which was largely embedded in inter-individual differences. All statistical approaches consistently showed that genetic similarity was not related to variation in the individual parameters describing growth responses. In contrast, up to 29% of the variance of individual parameters was accounted by spatial variables, revealing that microenvironmental features are more relevant than genetic similarity in determining similar growth patterns. Our study highlights both the advantages of modelling dendrochronological data at the individual level and the relevance of microenvironmental variation on individual growth patterns. These two aspects should be carefully considered in future multidisciplinary studies on growth dynamics in natural populations.

Dependence of leukocyte capture on instantaneous pulsatile flow.

Atherosclerosis, an artery disease, is currently the leading cause of death in the United States in both men and women. The first step in the development of atherosclerosis involves leukocyte adhesion to the arterial endothelium. It is broadly accepted that blood flow, more specifically wall shear stress (WSS), plays an important role in leukocyte capture and subsequent development of an atherosclerotic plaque. What is less known is how instantaneous WSS, which can vary by up to 5 Pa over one cardiac cycle, influences leukocyte capture. In this paper we use direct numerical simulations (DNS), performed using an in-house code, to illustrate that leukocyte capture is different whether as a function of instantaneous or time-averaged blood flow. Specifically, a stenotic plaque is modeled using a computational fluid dynamics (CFD) solver through fully three-dimensional Navier-Stokes equations and the immersed boundary method. Pulsatile triphasic inflow is used to simulate the cardiac cycle. The CFD is coupled with an agent-based leukocyte capture model to assess the impact of instantaneous hemodynamics on stenosis growth. The computed wall shear stress agrees well with the results obtained with a commercial software, as well as with theoretical results in the healthy region of the artery. The analysis emphasizes the importance of the instantaneous flow conditions in evaluating the leukocyte rate of capture. That is, the capture rate computed from mean flow field is generally underpredicted compared to the actual rate of capture. Thus, in order to obtain a reliable estimate, the flow unsteadiness during a cardiac cycle should be taken into account.

Nitrogen deposition outweighs climatic variability in driving annual growth rate of canopy beech trees: Evidence from long-term growth reconstruction across a geographic gradient.

In this study, we investigated the role of climatic variability and atmospheric nitrogen deposition in driving long-term tree growth in canopy beech trees along a geographic gradient in the montane belt of the Italian peninsula, from the Alps to the southern Apennines. We sampled dominant trees at different developmental stages (from young to mature tree cohorts, with tree ages spanning from 35 to 160 years) and used stem analysis to infer historic reconstruction of tree volume and dominant height. Annual growth volume (GV ) and height (GH ) variability were related to annual variability in model simulated atmospheric nitrogen deposition and site-specific climatic variables, (i.e. mean annual temperature, total annual precipitation, mean growing period temperature, total growing period precipitation, and standard precipitation evapotranspiration index) and atmospheric CO2 concentration, including tree cambial age among growth predictors. Generalized additive models (GAM), linear mixed-effects models (LMM), and Bayesian regression models (BRM) were independently employed to assess explanatory variables. The main results from our study were as follows: (i) tree age was the main explanatory variable for long-term growth variability; (ii) GAM, LMM, and BRM results consistently indicated climatic variables and CO2 effects on GV and GH were weak, therefore evidence of recent climatic variability influence on beech annual growth rates was limited in the montane belt of the Italian peninsula; (iii) instead, significant positive nitrogen deposition (Ndep ) effects were repeatedly observed in GV and GH ; the positive effects of Ndep on canopy height growth rates, which tended to level off at Ndep values greater than approximately 1.0 g m-2  y-1 , were interpreted as positive impacts on forest stand above-ground net productivity at the selected study sites.

Long-term temporal relationships between environmental conditions and xylem functional traits: a meta-analysis across a range of woody species along climatic and nitrogen deposition gradients.

The objectives of this study were to provide a quantitative description of the long-term effects of environmental variability on xylem functional traits, in order to better assess xylem capacity to change in response to climate change. Twenty-six sites throughout the world, primarily in Europe, were chosen where results from long-term measurements of anatomical traits were previously published. Published data on long-term xylem anatomy (conduit size and density) and ring width variability were compiled across a range of tree species, which was subsequently related to variability in temperature, precipitation and nitrogen deposition rates across the study sites using generalized additive models and Bayesian methods. We found some appreciable relationships between xylem traits (conduit area Ac and conduit density Dc) and environmental variables; whereas combined trait indices (lumen fraction: Ac × Dc and vessel composition: Ac/Dc) were found to be rather constant across a wide range of environmental conditions and to be decoupled from tree growth rates. Overall, results suggested xylem traits coordinated towards a homeostasis in xylem function, which appeared to act across a wide range of environmental conditions. Results showed also nitrogen deposition was associated with xylem traits and vessel composition: increased nitrogen availability due to nitrogen deposition might facilitate construction of a xylem structure efficient for water transport, and concurrently provide capacity to withstand the risks of drought-induced embolism.

Functional adjustments of xylem anatomy to climatic variability: insights from long-term Ilex aquifolium tree-ring series.

The present study assessed the effects of climatic conditions on radial growth and functional anatomical traits, including ring width, vessel size, vessel frequency and derived variables, i.e., potential hydraulic conductivity and xylem vulnerability to cavitation in Ilex aquifolium L. trees using long-term tree-ring time series obtained at two climatically contrasting sites, one mesic site in Switzerland (CH) and one drought-prone site in Italy (ITA). Relationships were explored by examining different xylem traits, and point pattern analysis was applied to investigate vessel clustering. We also used generalized additive models and bootstrap correlation functions to describe temperature and precipitation effects. Results indicated modified radial growth and xylem anatomy in trees over the last century; in particular, vessel frequency increased markedly at both sites in recent years, and all xylem traits examined, with the exception of xylem cavitation vulnerability, were higher at the CH mesic compared with the ITA drought site. A significant vessel clustering was observed at the ITA site, which could contribute to an enhanced tolerance to drought-induced embolism. Flat and negative relationships between vessel size and ring width were observed, suggesting carbon was not allocated to radial growth under conditions which favored stem water conduction. Finally, in most cases results indicated that climatic conditions influenced functional anatomical traits more substantially than tree radial growth, suggesting a crucial role of functional xylem anatomy in plant acclimation to future climatic conditions.

A global assessment of forest surface albedo and its relationships with climate and atmospheric nitrogen deposition.

We present a global assessment of the relationships between the short-wave surface albedo of forests, derived from the MODIS satellite instrument product at 0.5° spatial resolution, with simulated atmospheric nitrogen deposition rates (Ndep ), and climatic variables (mean annual temperature Tm and total annual precipitation P), compiled at the same spatial resolution. The analysis was performed on the following five forest plant functional types (PFTs): evergreen needle-leaf forests (ENF); evergreen broad-leaf forests (EBF); deciduous needle-leaf forests (DNF); deciduous broad-leaf forests (DBF); and mixed-forests (MF). Generalized additive models (GAMs) were applied in the exploratory analysis to assess the functional nature of short-wave surface albedo relations to environmental variables. The analysis showed evident correlations of albedo with environmental predictors when data were pooled across PFTs: Tm and Ndep displayed a positive relationship with forest albedo, while a negative relationship was detected with P. These correlations are primarily due to surface albedo differences between conifer and broad-leaf species, and different species geographical distributions. However, the analysis performed within individual PFTs, strengthened by attempts to select 'pure' pixels in terms of species composition, showed significant correlations with annual precipitation and nitrogen deposition, pointing toward the potential effect of environmental variables on forest surface albedo at the ecosystem level. Overall, our global assessment emphasizes the importance of elucidating the ecological mechanisms that link environmental conditions and forest canopy properties for an improved parameterization of surface albedo in climate models.

Genetic Drift Linked to Heterogeneous Landscape and Ecological Specialization Drives Diversification in the Alpine Endemic Columbine Aquilegia thalictrifolia.

The European Alpine system is an extensive mountain range, whose heterogeneous landscape together with Quaternary climatic oscillations significantly affected organismal diversity and distribution in Europe. The model genus Aquilegia represents a textbook example of a rapid and recent radiation through the Northern hemisphere, with the majority of the European taxa occuring in the Alpine system. However, the processes governing genetic differentiation of the genus in this complex geographic area are still widely unexplored. In this work, we used 9 microsatellite loci to study the genetic structure and diversity of 11 populations of Aquilegia thalictrifolia Schott & Kotschy, an alpine taxon characterized by a marked ecological specificity. We found that, despite the endemic and fragmented distribution, A. thalictrifolia has overall high levels of heterozygosity, which is consistent to the substantial inbreeding depression that characterizes the genus. Strong spatial genetic structuring of populations suggests a historical prevalence of genetic drift over gene flow, with natural barriers and ecological niche hindering migration. An analytical comparison of fixation and population differentiation indexes allowed us to infer hypotheses of the postglacial history and more recent demographic events that have influenced the genetics of the species. Overall, our results indicate allopatry as a major force of differentiation in the European scenario, likely to underlie the development of taxonomic boundaries in a broader geographic context. This adds to previous notions on the primary evolutionary forces shaping the Aquilegia radiation in Europe.

Within-population genetic structure in beech (Fagus sylvatica L.) stands characterized by different disturbance histories: does forest management simplify population substructure?

The fine-scale assessment of both spatially and non-spatially distributed genetic variation is crucial to preserve forest genetic resources through appropriate forest management. Cryptic within-population genetic structure may be more common than previously thought in forest tree populations, which has strong implications for the potential of forests to adapt to environmental change. The present study was aimed at comparing within-population genetic structure in European beech (Fagus sylvatica L.) plots experiencing different disturbance levels. Five plot pairs made up by disturbed and undisturbed plots having the same biogeographic history were sampled throughout Europe. Overall, 1298 individuals were analyzed using four highly polymorphic nuclear microsatellite markers (SSRs). Bayesian clustering within plots identified 3 to 11 genetic clusters (within-plot θ ST ranged from 0.025 to 0.124). The proportion of within-population genetic variation due to genetic substructuring (F CluPlot = 0.067) was higher than the differentiation among the 10 plots (F PlotTot = 0.045). Focusing on the comparison between managed and unmanaged plots, disturbance mostly explains differences in the complexity of within-population genetic structure, determining a reduction of the number of genetic clusters present in a standardized area. Our results show that: i) genetic substructuring needs to be investigated when studying the within-population genetic structure in forest tree populations, and ii) indices describing subtle characteristics of the within-population genetic structure are good candidates for providing early signals of the consequences of forest management, and of disturbance events in general.

Effect of habitat fragmentation on the genetic diversity and structure of peripheral populations of beech in Central Italy.

Fragmentation can affect the demographic and genetic structure of populations near the boundary of their biogeographic range. Higher genetic differentiation among populations coupled with lower level of within-population variability is expected as a consequence of reduced population size and isolation. The effects of these 2 factors have been rarely disentangled. Given their high gene flow, anemophilous forest trees should be more affected, in terms of loss of genetic diversity, by small population size rather than geographic isolation alone. We studied the impact of distance from the main range (a measure of isolation) and reduced population size on the within-population and among population components of genetic variability. We assayed 11 isozyme loci in a total of 856 individuals in 27 marginal populations of European beech (Fagus sylvatica L.) in Central Italy. Populations were divided into 3 groups with an increasing level of fragmentation. In the most fragmented group, the within-population genetic variability was slightly smaller and the among population differentiation significantly larger than in the other 2 groups. Isolation-by-distance was lost when only pairs of populations involving at least one from the most fragmented group were considered and maintained in the other groups. These results support the role of random genetic drift having a larger impact on the most fragmented group, whereas gene flow seems to balance genetic drift in the 2 less fragmented ones. Given that average distance from the main range is not different between the intermediate and the most fragmented group, but average population size is smaller, we can conclude that gene flow is effective, even at relatively long distances, in balancing the effect of fragmentation if population size is not too small.

Gene flow and simulation of transgene dispersal from hybrid poplar plantations.

Gene flow is a primary determinant of potential ecological impacts of transgenic trees. However, gene flow is a complex process that must be assessed in the context of realistic genetic, management, and environmental conditions. We measured gene flow from hybrid poplar plantations using morphological and genetic markers, and developed a spatially explicit landscape model to simulate pollination, dispersal, establishment, and mortality in the context of historical and projected disturbance and land-use regimes. Most pollination and seed establishment occurred within 450 m of the source, with a very long tail. Modeled transgene flow was highly context-dependent, strongly influenced by the competitive effects of transgenes, transgenic fertility, plantation rotation length, disturbance regime, and spatial and temporal variation in selection. The use of linked infertility genes even if imperfect, substantially reduced transgene flow in a wide range of modeled scenarios. The significance of seed and vegetative dispersal was highly dependent on plantation size. Our empirical and modeling studies suggest that transgene spread can be spatially extensive. However, the amount of spread is highly dependent on ecological and management context, and can be greatly limited or prevented by management or mitigation genes such as those that cause sexual infertility.

Assessing the effects of nitrogen deposition and climate on carbon isotope discrimination and intrinsic water-use efficiency of angiosperm and conifer trees under rising CO2 conditions.

The objective of this study is to globally assess the effects of atmospheric nitrogen deposition and climate, associated with rising levels of atmospheric CO2 , on the variability of carbon isotope discrimination (Δ(13) C), and intrinsic water-use efficiency (iWUE) of angiosperm and conifer tree species. Eighty-nine long-term isotope tree-ring chronologies, representing 23 conifer and 13 angiosperm species for 53 sites worldwide, were extracted from the literature, and used to obtain long-term time series of Δ(13) C and iWUE. Δ(13) C and iWUE were related to the increasing concentration of atmospheric CO2 over the industrial period (1850-2000) and to the variation of simulated atmospheric nitrogen deposition and climatic variables over the period 1950-2000. We applied generalized additive models and linear mixed-effects models to predict the effects of climatic variables and nitrogen deposition on Δ(13) C and iWUE. Results showed a declining Δ(13) C trend in the angiosperm and conifer species over the industrial period and a 16.1% increase of iWUE between 1850 and 2000, with no evidence that the increased rate was reduced at higher ambient CO2 values. The temporal variation in Δ(13) C supported the hypothesis of an active plant mechanism that maintains a constant ratio between intercellular and ambient CO2 concentrations. We defined linear mixed-effects models that were effective to describe the variation of Δ(13) C and iWUE as a function of a set of environmental predictors, alternatively including annual rate (Nrate ) and long-term cumulative (Ncum ) nitrogen deposition. No single climatic or atmospheric variable had a clearly predominant effect, however, Δ(13) C and iWUE showed complex dependent interactions between different covariates. A significant association of Nrate with iWUE and Δ(13) C was observed in conifers and in the angiosperms, and Ncum was the only independent term with a significant positive association with iWUE, although a multi-factorial control was evident in conifers.

Variability of stomatal conductance in a small and isolated population of silver fir (Abies alba Mill.).

We analyzed the response to drought of 420 individuals from eight half-sib families from a small and isolated population of silver fir (Abies alba Mill.): 105 of them were kept in well-watered conditions as control while the remaining 315 were exposed to drought for 27 days. A model describing stomatal behavior derived from Monteith and developed in beech by Leonardi et al. was fitted to experimental transpiration data obtained simply from the difference between two daily pot weighings. The estimated parameters were maximum stomatal conductance, maximum transpiration in well-watered conditions and sensitivity to soil water deficit. The model worked well: convergence for all but four individuals and concordance between experimental and fitted data were good (R(2)=0.86). Inter-individual variability for all three estimated parameters was high and two of them (maximum stomatal conductance and sensitivity to soil water deficit) were significantly different among families, suggesting genetic control. Our results validate the simplified method used to evaluate individual stomatal parameters. We also show that in the small and isolated population of our study substantial adaptive variability remains, a crucial prerequisite to endure environmental conditions determined by climatic change foreseen for the next decades.