Host genotype interacts with aerial spore communities and influences the needle mycobiome of Norway spruce.

The factors shaping the composition of the tree mycobiome are still under investigation. We tested the effects of host genotype, site, host phenotypic traits, and air fungal spore communities on the assembly of the fungi inhabiting Norway spruce needles. We used Norway spruce clones and spore traps within the collection sites and characterized both needle and air mycobiome communities by high-throughput sequencing of the ITS2 region. The composition of the needle mycobiome differed between Norway spruce clones, and clones with high genetic similarity had a more similar mycobiome. The needle mycobiome also varied across sites and was associated with the composition of the local air mycobiome and climate. Phenotypic traits such as diameter at breast height or crown health influenced the needle mycobiome to a lesser extent than host genotype and air mycobiome. Altogether, our results suggest that the needle mycobiome is mainly driven by the host genotype in combination with the composition of the local air spore communities. Our work highlights the role of host intraspecific variation in shaping the mycobiome of trees and provides new insights on the ecological processes structuring fungal communities inhabiting woody plants.

Teasing apart the joint effect of demography and natural selection in the birth of a contact zone.

Vast population movements induced by recurrent climatic cycles have shaped the genetic structure of plant species. During glacial periods species were confined to low-latitude refugia from which they recolonized higher latitudes as the climate improved. This multipronged recolonization led to many lineages that later met and formed large contact zones. We utilize genomic data from 5000 Picea abies trees to test for the presence of natural selection during recolonization and establishment of a contact zone in Scandinavia. Scandinavian P. abies is today made up of a southern genetic cluster originating from the Baltics, and a northern one originating from Northern Russia. The contact zone delineating them closely matches the limit between two major climatic regions. We show that natural selection contributed to its establishment and maintenance. First, an isolation-with-migration model with genome-wide linked selection fits the data better than a purely neutral one. Second, many loci show signatures of selection or are associated with environmental variables. These loci, regrouped in clusters on chromosomes, are often related to phenology. Altogether, our results illustrate how climatic cycles, recolonization and selection can establish strong local adaptation along contact zones and affect the genetic architecture of adaptive traits.

Killing two enemies with one stone? Genomics of resistance to two sympatric pathogens in Norway spruce.

Trees must cope with the attack of multiple pathogens, often simultaneously during their long lifespan. Ironically, the genetic and molecular mechanisms controlling this process are poorly understood. The objective of this study was to compare the genetic component of resistance in Norway spruce to Heterobasidion annosum s.s. and its sympatric congener Heterobasidion parviporum. Heterobasidion root- and stem-rot is a major disease of Norway spruce caused by members of the Heterobasidion annosum species complex. Resistance to both pathogens was measured using artificial inoculations in half-sib families of Norway spruce trees originating from central to northern Europe. The genetic component of resistance was analysed using 63,760 genome-wide exome-capture sequenced SNPs and multitrait genome-wide associations. No correlation was found for resistance to the two pathogens; however, associations were found between genomic variants and resistance traits with synergic or antagonist pleiotropic effects to both pathogens. Additionally, a latitudinal cline in resistance in the bark to H. annosum s.s. was found; trees from southern latitudes, with a later bud-set and thicker stem diameter, allowed longer lesions, but this was not the case for H. parviporum. In summary, this study detects genomic variants with pleiotropic effects which explain multiple disease resistance from a genic level and could be useful for selection of resistant trees to both pathogens. Furthermore, it highlights the need for additional research to understand the evolution of resistance traits to multiple pathogens in trees.

Genome-wide association study identified novel candidate loci affecting wood formation in Norway spruce.

Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.

Effect of number of annual rings and tree ages on genomic predictive ability for solid wood properties of Norway spruce.

BACKGROUND: Genomic selection (GS) or genomic prediction is considered as a promising approach to accelerate tree breeding and increase genetic gain by shortening breeding cycle, but the efforts to develop routines for operational breeding are so far limited. We investigated the predictive ability (PA) of GS based on 484 progeny trees from 62 half-sib families in Norway spruce (Picea abies (L.) Karst.) for wood density, modulus of elasticity (MOE) and microfibril angle (MFA) measured with SilviScan, as well as for measurements on standing trees by Pilodyn and Hitman instruments. RESULTS: GS predictive abilities were comparable with those based on pedigree-based prediction. Marker-based PAs were generally 25-30% higher for traits density, MFA and MOE measured with SilviScan than for their respective standing tree-based method which measured with Pilodyn and Hitman. Prediction accuracy (PC) of the standing tree-based methods were similar or even higher than increment core-based method. 78-95% of the maximal PAs of density, MFA and MOE obtained from coring to the pith at high age were reached by using data possible to obtain by drilling 3-5 rings towards the pith at tree age 10-12. CONCLUSIONS: This study indicates standing tree-based measurements is a cost-effective alternative method for GS. PA of GS methods were comparable with those pedigree-based prediction. The highest PAs were reached with at least 80-90% of the dataset used as training set. Selection for trait density could be conducted at an earlier age than for MFA and MOE. Operational breeding can also be optimized by training the model at an earlier age or using 3 to 5 outermost rings at tree age 10 to 12 years, thereby shortening the cycle and reducing the impact on the tree.

Genotypic variation in Norway spruce correlates to fungal communities in vegetative buds.

The taxonomically diverse phyllosphere fungi inhabit leaves of plants. Thus, apart from the fungi's dispersal capacities and environmental factors, the assembly of the phyllosphere community associated with a given host plant depends on factors encoded by the host's genome. The host genetic factors and their influence on the assembly of phyllosphere communities under natural conditions are poorly understood, especially in trees. Recent work indicates that Norway spruce (Picea abies) vegetative buds harbour active fungal communities, but these are hitherto largely uncharacterized. This study combines internal transcribed spacer sequencing of the fungal communities associated with dormant vegetative buds with a genome-wide association study (GWAS) in 478 unrelated Norway spruce trees. The aim was to detect host loci associated with variation in the fungal communities across the population, and to identify loci correlating with the presence of specific, latent, pathogens. The fungal communities were dominated by known Norway spruce phyllosphere endophytes and pathogens. We identified six quantitative trait loci (QTLs) associated with the relative abundance of the dominating taxa (i.e., top 1% most abundant taxa). Three additional QTLs associated with colonization by the spruce needle cast pathogen Lirula macrospora or the cherry spruce rust (Thekopsora areolata) in asymptomatic tissues were detected. The identification of the nine QTLs shows that the genetic variation in Norway spruce influences the fungal community in dormant buds and that mechanisms underlying the assembly of the communities and the colonization of latent pathogens in trees may be uncovered by combining molecular identification of fungi with GWAS.

Association genetics identifies a specifically regulated Norway spruce laccase gene, PaLAC5, linked to Heterobasidion parviporum resistance.

It is important to improve the understanding of the interactions between the trees and pathogens and integrate this knowledge about disease resistance into tree breeding programs. The conifer Norway spruce (Picea abies) is an important species for the forest industry in Europe. Its major pathogen is Heterobasidion parviporum, causing stem and root rot. In this study, we identified 11 Norway spruce QTLs (Quantitative trait loci) that correlate with variation in resistance to H. parviporum in a population of 466 trees by association genetics. Individual QTLs explained between 2.1 and 5.2% of the phenotypic variance. The expression of candidate genes associated with the QTLs was analysed in silico and in response to H. parviporum hypothesizing that (a) candidate genes linked to control of fungal sapwood growth are more commonly expressed in sapwood, and; (b) candidate genes associated with induced defences are respond to H. parviporum inoculation. The Norway spruce laccase PaLAC5 associated with control of lesion length development is likely to be involved in the induced defences. Expression analyses showed that PaLAC5 responds specifically and strongly in close proximity to the H. parviporum inoculation. Thus, PaLAC5 may be associated with the lignosuberized boundary zone formation in bark adjacent to the inoculation site.

Rapid testing of red blood cell parameters in primary care patients using HemoScreen™ point of care instrument.

BACKGROUND: Patients with anemia are frequently encountered in primary care. Once anemia is detected, it is essential to define the type and identify the underlying cause prior to initiation of treatment. In most cases, the cause can be determined using information from the patient history, physical exam, and complete blood counts (CBC). Point of care testing of blood cell counts would speed up the work up of anemia patients. The aim of the present study was to evaluate if the HemoScreen™ instrument (PixCell Medical, Yokneam Ilit, Israel) could be used for primary care samples. It is a POCT instrument that utilizes single sample cuvettes and image analysis of full blood count including RBC, Hemoglobin, MCV, MCH, platelets, WBC, and WBC 5-part differential. METHODS: We compared the HemoScreen™ and the Sysmex XN instrument results of 100 primary care patient samples focusing on the total white blood cells, red blood cell parameters RBC, Hemoglobin, MCH, MCV and platelets. RESULTS: Deming correlations between the HemoScreen™ and the Sysmex XN instruments for the CBC were WBCHemoScreen™ = 1.016* WBCSysmex + 0.34; r = 0.981, RBCHemoScreen™ = 0.988* RBCSysmex + 0.015; r = 0.974, HemoglobinHemoScreen™ = 1.081* HemoglobinSysmex - 11.25; r = 0.964, MCHHemoScreen™ = 0.978* MCHSysmex + 0.78; r = 0.939, MCVHemoScreen™ = 0.963* MCVSysmex + 8.68; r = 0.946, PlateletsHemoScreen™ = 0.964* PlateletsSysmex + 25.7; r = 0.953. CONCLUSION: The HemoScreen™ instrument could provide rapid and accurate test results for evaluation of the red blood cell parameters in primary care. This new technology is interesting as it allows the analysis red blood cell parameters also at small primary care centers.

Accuracy of genomic selection for growth and wood quality traits in two control-pollinated progeny trials using exome capture as the genotyping platform in Norway spruce.

BACKGROUND: Genomic selection (GS) can increase genetic gain by reducing the length of breeding cycle in forest trees. Here we genotyped 1370 control-pollinated progeny trees from 128 full-sib families in Norway spruce (Picea abies (L.) Karst.), using exome capture as genotyping platform. We used 116,765 high-quality SNPs to develop genomic prediction models for tree height and wood quality traits. We assessed the impact of different genomic prediction methods, genotype-by-environment interaction (G × E), genetic composition, size of the training and validation set, relatedness, and number of SNPs on accuracy and predictive ability (PA) of GS. RESULTS: Using G matrix slightly altered heritability estimates relative to pedigree-based method. GS accuracies were about 11-14% lower than those based on pedigree-based selection. The efficiency of GS per year varied from 1.71 to 1.78, compared to that of the pedigree-based model if breeding cycle length was halved using GS. Height GS accuracy decreased to more than 30% while using one site as training for GS prediction and using this model to predict the second site, indicating that G × E for tree height should be accommodated in model fitting. Using a half-sib family structure instead of full-sib structure led to a significant reduction in GS accuracy and PA. The full-sib family structure needed only 750 markers to reach similar accuracy and PA, as compared to 100,000 markers required for the half-sib family, indicating that maintaining the high relatedness in the model improves accuracy and PA. Using 4000-8000 markers in full-sib family structure was sufficient to obtain GS model accuracy and PA for tree height and wood quality traits, almost equivalent to that obtained with all markers. CONCLUSIONS: The study indicates that GS would be efficient in reducing generation time of breeding cycle in conifer tree breeding program that requires long-term progeny testing. The sufficient number of trees within-family (16 for growth and 12 for wood quality traits) and number of SNPs (8000) are required for GS with full-sib family relationship. GS methods had little impact on GS efficiency for growth and wood quality traits. GS model should incorporate G × E effect when a strong G × E is detected.

Evaluation of the Alere D-dimer test for point of care testing.

The primary care regularly sees patients that have symptoms that could be due to thromboembolic diseases. It would be valuable to be able to rule out deep venous thrombosis or pulmonary embolism using Wells score and a negative D-dimer testing already at the primary care unit. This requires a validated D-dimer assay suitable for primary care use. We compared D-dimer results obtained with the new point of care analyzer Alere Triage(®) and the central hospital laboratory STA-R Evolution analyzer from the same patient samples (n = 102). We also calculated the total coefficient of variation (CV) for the Alere method. The two methods showed a good linear correlation (R(2) = 0.977) and a slope of 0.975. CV for the Alere D-dimer method was well below 10%. The study shows that the Alere D-dimer assay and the central laboratory standard assay show similar results. We suggest that the Alere D-dimer assay could be used in primary care in combination with Wells score to reduce referrals to the emergency unit.

Neuropeptide Y and measures of stress in a longitudinal study of women with the fibromyalgia syndrome.

OBJECTIVES: Neuropeptide Y is associated with stress in animal and human laboratory studies. However, data from clinical studies are scarce and no clinical longitudinal studies have been published. The aim of this clinical study was to assess the possible association between changes in the levels of pain, depression, and stress measures, on the one hand, and plasma neuropeptide Y levels, on the other. METHODS: Forty-four women with the fibromyalgia syndrome were exposed to a Cognitive Behavioral Therapy intervention. Levels of the plasma neuropeptide Y as well as pain, depression, and stress measures were obtained at the start and at the end of the intervention, and after a further six month follow-up. Based on these data, a before-and-after analysis was performed. RESULTS: Almost all measures of pain, depression, and stress improved during the study; specifically, variables measuring life control (coping), depression, and stress-related time urgency improved significantly. Moreover, during the same time period, the mean plasma neuropeptide Y level was reduced from 93.2 ± 38.8 fmol/mL before the Cognitive Behavioral Therapy to 75.6 ± 42.9 fmol/mL (p<0.001) at the end of the study. CONCLUSIONS: After exposure to a Cognitive Behavioral Therapy intervention, levels of most of the pain, depression, and stress measures improved, half of them significantly, as did the levels of neuropeptide Y. This circumstance indicates a possible functional relationship between pain-depression-stress and neuropeptide Y.

Divergent selection predating the Last Glacial Maximum mainly acted on macro-phenotypes in Norway spruce.

The current distribution and population structure of many species were, to a large extent, shaped by cycles of isolation in glacial refugia and subsequent population expansions. Isolation in and postglacial expansion through heterogeneous environments led to either neutral or adaptive divergence. Norway spruce is no exception, and its current distribution is the consequence of a constant interplay between evolutionary and demographic processes. We investigated population differentiation and adaptation of Norway spruce for juvenile growth, diameter of the stem, wood density, and tracheid traits at breast height. Data from 4461 phenotyped and genotyped Norway spruce from 396 half-sib families in two progeny tests were used to test for divergent selection in the framework of Q ST vs. F ST. We show that the macroscopic resultant trait (stem diameter), unlike its microscopic components (tracheid dimensions) and juvenile growth, was under divergent selection that predated the Last Glacial Maximum. Altogether, the current variation in these phenotypic traits in Norway spruce is better explained by local adaptation to ancestral environments than to current ones, where populations were partly preadapted, mainly through growth-related traits.

Inflammation and Interferon Signatures in Peripheral B-Lymphocytes and Sera of Individuals With Fibromyalgia.

Fibromyalgia (FM) is an idiopathic chronic disease characterized by widespread musculoskeletal pain, hyperalgesia and allodynia, often accompanied by fatigue, cognitive dysfunction and other symptoms. Autoimmunity and neuroinflammatory mechanisms have been suggested to play important roles in the pathophysiology of FM supported by recently identified interferon signatures in affected individuals. However, the contribution of different components in the immune system, such as the B-lymphocytes, in the progression to FM are yet unknown. Furthermore, there is a great need for biomarkers that may improve diagnostics of FM. Herein, we investigated the gene expression profile in peripheral B-cells, as well as a panel of inflammatory serum proteins, in 30 FM patients and 23 healthy matched control individuals. RNA sequence analysis revealed 60 differentially expressed genes when comparing the two groups. The group of FM patients showed increased expression of twenty-five interferon-regulated genes, such as S100A8 and S100A9, VCAM, CD163, SERPINA1, ANXA1, and an increased interferon score. Furthermore, FM was associated with elevated levels of 19 inflammatory serum proteins, such as IL8, AXIN1, SIRT2 and STAMBP, that correlated with the FM severity score. Together, the results shows that FM is associated with an interferon signature in B-cells and increased levels of a set of inflammatory serum proteins. Our findings bring further support for immune activation in the pathogenesis of FM and highlight candidate biomarkers for diagnosis and intervention in the management of FM.

Genetic architecture behind developmental and seasonal control of tree growth and wood properties in Norway spruce.

Genetic control of tree growth and wood formation varies depending on the age of the tree and the time of the year. Single-locus, multi-locus, and multi-trait genome-wide association studies (GWAS) were conducted on 34 growth and wood property traits in 1,303 Norway spruce individuals using exome capture to cover ~130K single-nucleotide polymorphisms (SNPs). GWAS identified associations to the different wood traits in a total of 85 gene models, and several of these were validated in a progenitor population. A multi-locus GWAS model identified more SNPs associated with the studied traits than single-locus or multivariate models. Changes in tree age and annual season influenced the genetic architecture of growth and wood properties in unique ways, manifested by non-overlapping SNP loci. In addition to completely novel candidate genes, SNPs were located in genes previously associated with wood formation, such as cellulose synthases and a NAC transcription factor, but that have not been earlier linked to seasonal or age-dependent regulation of wood properties. Interestingly, SNPs associated with the width of the year rings were identified in homologs of Arabidopsis thaliana BARELY ANY MERISTEM 1 and rice BIG GRAIN 1, which have been previously shown to control cell division and biomass production. The results provide tools for future Norway spruce breeding and functional studies.

Hydraulic efficiency compromises compression strength perpendicular to the grain in Norway spruce trunkwood.

The aim of this study was to investigate bending stiffness and compression strength perpendicular to the grain of Norway spruce (Picea abies (L.) Karst.) trunkwood with different anatomical and hydraulic properties. Hydraulically less safe mature sapwood had bigger hydraulic lumen diameters and higher specific hydraulic conductivities than hydraulically safer juvenile wood. Bending stiffness (MOE) was higher, whereas radial compression strength lower in mature than in juvenile wood. A density-based tradeoff between MOE and hydraulic efficiency was apparent in mature wood only. Across cambial age, bending stiffness did not compromise hydraulic efficiency due to variation in latewood percent and because of the structural demands of the tree top (e.g. high flexibility). Radial compression strength compromised, however, hydraulic efficiency because it was extremely dependent on the characteristics of the "weakest" wood part, the highly conductive earlywood. An increase in conduit wall reinforcement of earlywood tracheids would be too costly for the tree. Increasing radial compression strength by modification of microfibril angles or ray cell number could result in a decrease of MOE, which would negatively affect the trunk's capability to support the crown. We propose that radial compression strength could be an easily assessable and highly predictive parameter for the resistance against implosion or vulnerability to cavitation across conifer species, which should be topic of further studies.

Leveraging breeding programs and genomic data in Norway spruce (Picea abies L. Karst) for GWAS analysis.

BACKGROUND: Genome-wide association studies (GWAS) identify loci underlying the variation of complex traits. One of the main limitations of GWAS is the availability of reliable phenotypic data, particularly for long-lived tree species. Although an extensive amount of phenotypic data already exists in breeding programs, accounting for its high heterogeneity is a great challenge. We combine spatial and factor-analytics analyses to standardize the heterogeneous data from 120 field experiments of 483,424 progenies of Norway spruce to implement the largest reported GWAS for trees using 134 605 SNPs from exome sequencing of 5056 parental trees. RESULTS: We identify 55 novel quantitative trait loci (QTLs) that are associated with phenotypic variation. The largest number of QTLs is associated with the budburst stage, followed by diameter at breast height, wood quality, and frost damage. Two QTLs with the largest effect have a pleiotropic effect for budburst stage, frost damage, and diameter and are associated with MAP3K genes. Genotype data called from exome capture, recently developed SNP array and gene expression data indirectly support this discovery. CONCLUSION: Several important QTLs associated with growth and frost damage have been verified in several southern and northern progeny plantations, indicating that these loci can be used in QTL-assisted genomic selection. Our study also demonstrates that existing heterogeneous phenotypic data from breeding programs, collected over several decades, is an important source for GWAS and that such integration into GWAS should be a major area of inquiry in the future.

Combining transcriptomics and genetic linkage based information to identify candidate genes associated with Heterobasidion-resistance in Norway spruce.

The Heterobasidion annosum s.l species complex comprises the most damaging forest pathogens to Norway spruce. We revisited previously identified Quantitative Trait Loci (QTLs) related to Heterobasidion-resistance in Norway spruce to identify candidate genes associated with these QTLs. We identified 329 candidate genes associated with the resistance QTLs using a gene-based composite map for Pinaceae. To evaluate the transcriptional responses of these candidate genes to H. parviporum, we inoculated Norway spruce plants and sequenced the transcriptome of the interaction at 3 and 7 days post inoculation. Out of 298 expressed candidate genes 124 were differentially expressed between inoculation and wounding control treatment. Interestingly, PaNAC04 and two of its paralogs in the subgroup III-3 of the NAC family transcription factors were found to be associated with one of the QTLs and was also highly induced in response to H. parviporum. These genes are possibly involved in the regulation of biosynthesis of flavonoid compounds. Furthermore, several of the differentially expressed candidate genes were associated with the phenylpropanoid pathway including a phenylalanine ammonia-lyase, a cinnamoyl-CoA reductase, a caffeoyl-CoA O-methyltransferase and a PgMYB11-like transcription factor gene. Combining transcriptome and genetic linkage analyses can help identifying candidate genes for functional studies and molecular breeding in non-model species.

Primary adrenal lymphoma as a cause of adrenal insufficiency, a report of two cases.

SUMMARY: Primary adrenal lymphoma (PAL) is a rare cause of adrenal insufficiency. More than 90% is of B-cell origin. The condition is bilateral in up to 75% of cases, with adrenal insufficiency in two of three patients. We report two cases of adrenal insufficiency presenting at the age of 70 and 79 years, respectively. Both patients had negative 21-hydroxylase antibodies with bilateral adrenal lesions on CT. Biopsy showed B-cell lymphoma. One of the patients experienced intermittent disease regression on replacement dosage of glucocorticoids. LEARNING POINTS: Primary adrenal lymphoma (PAL) is a rare cause of adrenal insufficiency. Bilateral adrenal masses of unknown origin or in individuals with suspected extra-adrenal malignancy should be biopsied quickly when pheochromocytoma is excluded biochemically. Steroid treatment before biopsy may affect diagnosis. Adrenal insufficiency with negative 21-hydroxylase antibodies should be evaluated radiologically.

Genetic control of tracheid properties in Norway spruce wood.

Through the use of genome-wide association studies (GWAS) mapping it is possible to establish the genetic basis of phenotypic trait variation. Our GWAS study presents the first such effort in Norway spruce (Picea abies (L). Karst.) for the traits related to wood tracheid characteristics. The study employed an exome capture genotyping approach that generated 178 101 Single Nucleotide Polymorphisms (SNPs) from 40 018 probes within a population of 517 Norway spruce mother trees. We applied a least absolute shrinkage and selection operator (LASSO) based association mapping method using a functional multi-locus mapping approach, with a stability selection probability method as the hypothesis testing approach to determine significant Quantitative Trait Loci (QTLs). The analysis has provided 30 significant associations, the majority of which show specific expression in wood-forming tissues or high ubiquitous expression, potentially controlling tracheids dimensions, their cell wall thickness and microfibril angle. Among the most promising candidates based on our results and prior information for other species are: Picea abies BIG GRAIN 2 (PabBG2) with a predicted function in auxin transport and sensitivity, and MA_373300g0010 encoding a protein similar to wall-associated receptor kinases, which were both associated with cell wall thickness. The results demonstrate feasibility of GWAS to identify novel candidate genes controlling industrially-relevant tracheid traits in Norway spruce.

Shrinkage processes in standard-size Norway spruce wood specimens with different vulnerability to cavitation.

The aim of this study was to observe the radial shrinkage of Norway spruce [Picea abies (L. Karst.)] trunkwood specimens with different hydraulic vulnerability to cavitation from the fully saturated state until the overall shrinkage reaches a stable value, and to relate wood shrinkage and recovery from shrinkage to cavitations of the water column inside the tracheids. Radial shrinkage processes in standard-size sapwood specimens (6 mm x 6 mm x 100 mm; radial, tangential and longitudinal) obtained at different positions within the trunk, representing different ages of the cambium, were compared. Cavitation events were assessed by acoustic emission (AE) testing, hydraulic vulnerability by the AE feature analysis and shrinkage was calculated from the changes in contact pressure between the 150 kHz AE transducer and the wood specimen. Two shrinkage processes were observed in both juvenile (annual rings 1 and 2) and mature wood (annual rings 17-19), the first one termed tension shrinkage and the second one cell wall shrinkage process, which started when most of the tracheids reached relative water contents below fiber saturation. Maximum tension shrinkage coincided with high-energy AEs, and the periods of shrinkage recovery could be traced to tension release due to cavitation. Juvenile wood, which was less sensitive to cavitation, had lower earlywood tracheid diameters and was less prone to deformation due to tensile strain than mature wood, showed a lower cell wall shrinkage, and thus total shrinkage. Earlywood lumen diameters and maximum tension shrinkage were strongly positively related to each other, meaning that bigger tracheids are more prone to deformation at the same water tension than the smaller tracheids.