Rhytidome- and cork-type barks of holm oak, cork oak and their hybrids highlight processes leading to cork formation.

BACKGROUND: The periderm is basic for land plants due to its protective role during radial growth, which is achieved by the polymers deposited in the cell walls. In most trees, like holm oak, the first periderm is frequently replaced by subsequent internal periderms yielding a heterogeneous outer bark made of a mixture of periderms and phloem tissues, known as rhytidome. Exceptionally, cork oak forms a persistent or long-lived periderm which results in a homogeneous outer bark of thick phellem cell layers known as cork. Cork oak and holm oak distribution ranges overlap to a great extent, and they often share stands, where they can hybridize and produce offspring showing a rhytidome-type bark. RESULTS: Here we use the outer bark of cork oak, holm oak, and their natural hybrids to analyse the chemical composition, the anatomy and the transcriptome, and further understand the mechanisms underlying periderm development. We also include a unique natural hybrid individual corresponding to a backcross with cork oak that, interestingly, shows a cork-type bark. The inclusion of hybrid samples showing rhytidome-type and cork-type barks is valuable to approach cork and rhytidome development, allowing an accurate identification of candidate genes and processes. The present study underscores that abiotic stress and cell death are enhanced in rhytidome-type barks whereas lipid metabolism and cell cycle are enriched in cork-type barks. Development-related DEGs showing the highest expression, highlight cell division, cell expansion, and cell differentiation as key processes leading to cork or rhytidome-type barks. CONCLUSION: Transcriptome results, in agreement with anatomical and chemical analyses, show that rhytidome and cork-type barks are active in periderm development, and suberin and lignin deposition. Development and cell wall-related DEGs suggest that cell division and expansion are upregulated in cork-type barks whereas cell differentiation is enhanced in rhytidome-type barks.

Bioprospection of the bacterial ß-myrcene-biotransforming trait in the rhizosphere.

The biocatalysis of ß-myrcene into value-added compounds, with enhanced organoleptic/therapeutic properties, may be performed by resorting to specialized enzymatic machinery of ß-myrcene-biotransforming bacteria. Few ß-myrcene-biotransforming bacteria have been studied, limiting the diversity of genetic modules/catabolic pathways available for biotechnological research. In our model Pseudomonas sp. strain M1, the ß-myrcene catabolic core-code was identified in a 28-kb genomic island (GI). The lack of close homologs of this ß-myrcene-associated genetic code prompted a bioprospection of cork oak and eucalyptus rhizospheres, from 4 geographic locations in Portugal, to evaluate the environmental diversity and dissemination of the ß-myrcene-biotransforming genetic trait (Myr+). Soil microbiomes were enriched in ß-myrcene-supplemented cultures, from which ß-myrcene-biotransforming bacteria were isolated, belonging to Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, and Sphingobacteriia classes. From a panel of representative Myr+ isolates that included 7 bacterial genera, the production of ß-myrcene derivatives previously reported in strain M1 was detected in Pseudomonas spp., Cupriavidus sp., Sphingobacterium sp., and Variovorax sp. A comparative genomics analysis against the genome of strain M1 found the M1-GI code in 11 new Pseudomonas genomes. Full nucleotide conservation of the ß-myrcene core-code was observed throughout a 76-kb locus in strain M1 and all 11 Pseudomonas spp., resembling the structure of an integrative and conjugative element (ICE), despite being isolated from different niches. Furthermore, the characterization of isolates not harboring the Myr+-related 76-kb locus suggested that they may biotransform ß-myrcene via alternative catabolic loci, being thereby a novel source of enzymes and biomolecule catalogue for biotechnological exploitation. KEY POINTS: • The isolation of 150 Myr+ bacteria hints the ubiquity of such trait in the rhizosphere. • The Myr+ trait is spread across different bacterial taxonomic classes. • The core-code for the Myr+ trait was detected in a novel ICE, only found in Pseudomonas spp.

The influence of bioclimate on soil microbial communities of cork oak.

BACKGROUND: Soil microbiomes are important to maintain soil processes in forests and confer protection to plants against abiotic and biotic stresses. These microbiomes can be affected by environmental changes. In this work, soil microbial communities from different cork oak Portuguese forests under different edaphoclimatic conditions were described by using a metabarcoding strategy targeting ITS2 and 16S barcodes. RESULTS: A total of 11,974 fungal and 12,010 bacterial amplicon sequence variants (ASVs) were obtained, revealing rich and diverse microbial communities associated with different cork oak forests. Bioclimate was described as the major factor influencing variability in these communities (or bioclimates/cork oak forest for fungal community), followed by boron and granulometry. Also, pH explained variation of fungal communities, while C:N ratio contributed to bacterial variation. Fungal and bacterial biomarker genera for specific bioclimates were described. Their co-occurrence network revealed the existence of a complex and delicate balance among microbial communities. CONCLUSIONS: The findings revealed that bacterial communities are more likely to be affected by different edaphoclimatic conditions than fungal communities, also predicting a higher impact of climate change on bacterial communities. The integration of cork oak fungal and bacterial microbiota under different bioclimates could be further explored to provide information about useful interactions for increasing cork oak forest sustainability in a world subject to climate changes.

Nonlinear plant-plant interactions modulate impact of extreme drought and recovery on a Mediterranean ecosystem.

Interaction effects of different stressors, such as extreme drought and plant invasion, can have detrimental effects on ecosystem functioning and recovery after drought. With ongoing climate change and increasing plant invasion, there is an urgent need to predict the short- and long-term interaction impacts of these stressors on ecosystems. We established a combined precipitation exclusion and shrub invasion (Cistus ladanifer) experiment in a Mediterranean cork oak (Quercus suber) ecosystem with four treatments: (1) Q. suber control; (2) Q. suber with rain exclusion; (3) Q. suber invaded by shrubs; and (4) Q. suber with rain exclusion and shrub invasion. As key parameter, we continuously measured ecosystem water fluxes. In an average precipitation year, the interaction effects of both stressors were neutral. However, the combination of imposed drought and shrub invasion led to amplifying interaction effects during an extreme drought by strongly reducing tree transpiration. Contrarily, the imposed drought reduced the competitiveness of the shrubs in the following recovery period, which buffered the negative effects of shrub invasion on Q. suber. Our results demonstrate the highly dynamic and nonlinear effects of interacting stressors on ecosystems and urges for further investigations on biotic interactions in a context of climate change pressures.

Small molecule inhibitors of mammalian GSK-3ß promote in vitro plant cell reprogramming and somatic embryogenesis in crop and forest species.

Plant in vitro regeneration systems, such as somatic embryogenesis, are essential in breeding; they permit propagation of elite genotypes, production of doubled-haploids, and regeneration of whole plants from gene editing or transformation events. However, in many crop and forest species, somatic embryogenesis is highly inefficient. We report a new strategy to improve in vitro embryogenesis using synthetic small molecule inhibitors of mammalian glycogen synthase kinase 3ß (GSK-3ß), never used in plants. These inhibitors increased in vitro embryo production in three different systems and species, microspore embryogenesis of Brassica napus and Hordeum vulgare, and somatic embryogenesis of Quercus suber. TDZD-8, a representative compound of the molecules tested, inhibited GSK-3 activity in microspore cultures, and increased expression of embryogenesis genes FUS3, LEC2, and AGL15. Plant GSK-3 kinase BIN2 is a master regulator of brassinosteroid (BR) signalling. During microspore embryogenesis, BR biosynthesis and signalling genes CPD, GSK-3-BIN2, BES1, and BZR1 were up-regulated and the BAS1 catabolic gene was repressed, indicating activation of the BR pathway. TDZD-8 increased expression of BR signalling elements, mimicking BR effects. The findings support that the small molecule inhibitors promoted somatic embryogenesis by activating the BR pathway, opening up the way for new strategies using GSK-3ß inhibitors that could be extended to other species.

Bacteria could help ectomycorrhizae establishment under climate variations.

Rhizosphere microbiome is one of the main sources of plant protection against drought. Beneficial symbiotic microorganisms, such as ectomycorrhizal fungi (ECMF) and mycorrhiza helper bacteria (MHB), interact with each other for increasing or maintaining host plant fitness. This mutual support benefits all three partners and comprises a natural system for drought acclimation in plants. Cork oak (Quercus suber L.) tolerance to drought scenarios is widely known, but adaptation to climate changes has been a challenge for forest sustainability protection. In this work, ECMF and MHB communities from cork oak forests were cross-linked and correlated with climates. Cenococcum, Russula and Tuber were the most abundant ECMF capable of interacting with MHB (ECMF~MHB) genera in cork oak stands, while Bacillus, Burkholderia and Streptomyces were the most conspicuous MHB. Integrating all microbial data, two consortia Lactarius/Bacillaceae and Russula/Burkholderaceae have singled out but revealed a negative interaction with each other. Russula/Burkholderaceae might have an important role for cork oak forest sustainability in arid environments, which will be complemented by the lower drought adaptation of competitive Lactarius/Bacillaceae. These microbial consortia could play an essential role on cork oak forest resilience to upcoming climatic changes.

Genome-Wide Identification of Epigenetic Regulators in Quercus suber L.

Modifications of DNA and histones, including methylation and acetylation, are critical for the epigenetic regulation of gene expression during plant development, particularly during environmental adaptation processes. However, information on the enzymes catalyzing all these modifications in trees, such as Quercus suber L., is still not available. In this study, eight DNA methyltransferases (DNA Mtases) and three DNA demethylases (DDMEs) were identified in Q. suber. Histone modifiers involved in methylation (35), demethylation (26), acetylation (8), and deacetylation (22) were also identified in Q. suber. In silico analysis showed that some Q. suber DNA Mtases, DDMEs and histone modifiers have the typical domains found in the plant model Arabidopsis, which might suggest a conserved functional role. Additional phylogenetic analyses of the DNA and histone modifier proteins were performed using several plant species homologs, enabling the classification of the Q. suber proteins. A link between the expression levels of each gene in different Q. suber tissues (buds, flowers, acorns, embryos, cork, and roots) with the functions already known for their closest homologs in other species was also established. Therefore, the data generated here will be important for future studies exploring the role of epigenetic regulators in this economically important species.

A comparative transcriptomic approach to understanding the formation of cork.

KEY MESSAGE: The transcriptome comparison of two oak species reveals possible candidates accounting for the exceptionally thick and pure cork oak phellem, such as those involved in secondary metabolism and phellogen activity. Cork oak, Quercus suber, differs from other Mediterranean oaks such as holm oak (Quercus ilex) by the thickness and organization of the external bark. While holm oak outer bark contains sequential periderms interspersed with dead secondary phloem (rhytidome), the cork oak outer bark only contains thick layers of phellem (cork rings) that accumulate until reaching a thickness that allows industrial uses. Here we compare the cork oak outer bark transcriptome with that of holm oak. Both transcriptomes present similitudes in their complexity, but whereas cork oak external bark is enriched with upregulated genes related to suberin, which is the main polymer responsible for the protective function of periderm, the upregulated categories of holm oak are enriched in abiotic stress and chromatin assembly. Concomitantly with the upregulation of suberin-related genes, there is also induction of regulatory and meristematic genes, whose predicted activities agree with the increased number of phellem layers found in the cork oak sample. Further transcript profiling among different cork oak tissues and conditions suggests that cork and wood share many regulatory mechanisms, probably reflecting similar ontogeny. Moreover, the analysis of transcripts accumulation during the cork growth season showed that most regulatory genes are upregulated early in the season when the cork cambium becomes active. Altogether our work provides the first transcriptome comparison between cork oak and holm oak outer bark, which unveils new regulatory candidate genes of phellem development.

ChIP-Seq reveals that QsMYB1 directly targets genes involved in lignin and suberin biosynthesis pathways in cork oak (Quercus suber).

BACKGROUND: Gene activity is largely controlled by transcriptional regulation through the action of transcription factors and other regulators. QsMYB1 is a member of the R2R3-MYB transcription factor family related to secondary growth, and in particular, with the cork development process. In order to identify the putative gene targets of QsMYB1 across the cork oak genome we developed a ChIP-Seq strategy. RESULTS: Results provide direct evidence that QsMY1B targets genes encoding for enzymes involved in the lignin and suberin pathways as well as gene encoding for ABCG transporters and LTPs implicated in the transport of monomeric suberin units across the cellular membrane. These results highlight the role of QsMYB1 as a regulator of lignin and suberin biosynthesis, transport and assembly. CONCLUSION: To our knowledge, this work constitutes the first ChIP-Seq experiment performed in cork oak, a non-model plant species with a long-life cycle, and these results will contribute to deepen the knowledge about the molecular mechanisms of cork formation and differentiation.

Ectomycorrhizal fungal diversity and community structure associated with cork oak in different landscapes.

Cork oak (Quercus suber L.) forests play an important ecological and economic role. Ectomycorrhizal fungi (ECMF) are key components for the sustainability and functioning of these ecosystems. The community structure and composition of ECMF associated with Q. suber in different landscapes of distinct Mediterranean bioclimate regions have not previously been compared. In this work, soil samples from cork oak forests residing in different bioclimates (arid, semi-arid, sub-humid, and humid) were collected and surveyed for ectomycorrhizal (ECM) root tips. A global analysis performed on 3565 ECM root tips revealed that the ECMF community is highly enriched in Russula, Tomentella, and Cenoccocum, which correspond to the ECMF genera that mainly contribute to community differences. The ECMF communities from the rainiest and the driest cork oak forests were distinct, with soils from the rainiest climates being more heterogeneous than those from the driest climates. The analyses of several abiotic factors on the ECMF communities revealed that bioclimate, precipitation, soil texture, and forest management strongly influenced ECMF structure. Shifts in ECMF with different hyphal exploration types were also detected among forests, with precipitation, forest system, and soil texture being the main drivers controlling their composition. Understanding the effects of environmental factors on the structuring of ECM communities could be the first step for promoting the sustainability of this threatened ecosystem.

Ectomycorrhizal inoculation with Pisolithus tinctorius reduces stress induced by drought in cork oak.

We investigated whether the performance of cork oak under drought could be improved by colonization with the ectomycorrhizal fungus Pisolithus tinctorius. Results show that inoculation alone had a positive effect on plant height, shoot biomass, shoot basal diameter, and root growth. Under drought, root growth of mycorrhizal plants was significantly increased showing that inoculation was effective in increasing tolerance to drought. In accordance, mycorrhizal plants subjected to drought showed less symptoms of stress when compared to non-mycorrhizal plants, such as lower concentration of soluble sugars and starch, increased ability to maintain fatty acid content and composition, and increased unsaturation level of membrane lipids. After testing some of the mechanisms suggested to contribute to the enhanced tolerance of mycorrhizal plants to drought, we could not find any by which Pisolithus tinctorius could benefit cork oak, at least under the drought conditions imposed in our experiment. Inoculation did not increase photosynthesis under drought, suggesting no effect in sustaining stomatal opening at low soil water content. Similarly, plant water status was not affected by inoculation suggesting that P. tinctorius does not contribute to an increased plant water uptake during drought. Inoculation did increase nitrogen concentration in plants but it was independent of the water status. Furthermore, no significant mycorrhizal effect on drought-induced ROS production or osmotic adjustment was detected, suggesting that these factors are not important for the improved drought tolerance triggered by P. tinctorius.

Characterization of the cork formation and production transcriptome in Quercus cerris × suber hybrids.

Cork oak is the main cork-producing species worldwide, and plays a significant economic, ecological and social role in the Mediterranean countries, in particular in Portugal and Spain. The ability to produce cork is limited to a few species, hence it must involve specific regulation mechanisms that are unique to these species. However, to date, these mechanisms remain largely understudied, especially with approaches involving the use of high-throughput sequencing technology. In this study, the transcriptome of cork-producing and non-cork-producing Quercus cerris × suber hybrids was analyzed in order to elucidate the differences between the two groups of trees displaying contrasting phenotypes for cork production. The results revealed the presence of a significant number of genes exclusively associated with cork production, in the trees that developed cork. Moreover, several gene ontology subcategories, such as cell wall biogenesis, lipid metabolic processes, metal ion binding and apoplast/cell wall, were only detected in the trees with cork production. These results indicate the existence, at the transcriptome level, of mechanisms that seem to be unique and necessary for cork production, which is an advancement in our knowledge regarding the genetic regulation behind cork formation and production.

Landscape dynamics in Mediterranean oak forests under global change: understanding the role of anthropogenic and environmental drivers across forest types.

The Mediterranean region is projected to be extremely vulnerable to global change, which will affect the distribution of typical forest types such as native oak forests. However, our understanding of Mediterranean oak forest responses to future conditions is still very limited by the lack of knowledge on oak forest dynamics and species-specific responses to multiple drivers. We compared the long-term (1966-2006) forest persistence and land cover change among evergreen (cork oak and holm oak) and deciduous oak forests and evaluated the importance of anthropogenic and environmental drivers on observed changes for Portugal. We used National Forest Inventories to quantify the changes in oak forests and explored the drivers of change using multinomial logistic regression analysis and an information theoretical approach. We found distinct trends among oak forest types, reflecting the differences in oak economic value, protection status and management schemes: cork oak forests were the most persistent (62%), changing mostly to pines and eucalypt; holm oak forests were less persistent (53.2%), changing mostly to agriculture; and deciduous oak forests were the least persistent (45.7%), changing mostly to shrublands. Drivers of change had distinct importance across oak forest types, but drivers from anthropogenic origin (wildfires, population density, and land accessibility) were always among the most important. Climatic extremes were also important predictors of oak forest changes, namely extreme temperatures for evergreen oak forests and deficit of precipitation for deciduous oak forests. Our results indicate that under increasing human pressure and forecasted climate change, evergreen oak forests will continue declining and deciduous oak forests will be replaced by forests dominated by more xeric species. In the long run, multiple disturbances may change competitive dominance from oak forests to pyrophytic shrublands. A better understanding of forest dynamics and the inclusion of anthropogenic drivers on models of vegetation change will improve predicting the future of Mediterranean oak forests.

Oak protein profile alterations upon root colonization by an ectomycorrhizal fungus.

An increased knowledge on the real impacts of ectomycorrhizal symbiosis in forest species is needed to optimize forest sustainable productivity and thus to improve forest services and their capacity to act as carbon sinks. In this study, we investigated the response of an oak species to ectomycorrhizae formation using a proteomics approach complemented by biochemical analysis of carbohydrate levels. Comparative proteome analysis between mycorrhizal and nonmycorrhizal cork oak plants revealed no differences at the foliar level. However, the protein profile of 34 unique oak proteins was altered in the roots. Consistent with the results of the biochemical analysis, the proteome analysis of the mycorrhizal roots suggests a decreasing utilization of sucrose for the metabolic activity of mycorrhizal roots which is consistent with an increased allocation of carbohydrates from the plant to the fungus in order to sustain the symbiosis. In addition, a promotion of protein unfolding mechanisms, attenuation of defense reactions, increased nutrient mobilization from the plant-fungus interface (N and P), as well as cytoskeleton rearrangements and induction of plant cell wall loosening for fungal root accommodation in colonized roots are also suggested by the results. The suggested improvement in root capacity to take up nutrients accompanied by an increase of root biomass without apparent changes in aboveground biomass strongly re-enforces the potential of mycorrhizal inoculation to improve cork oak forest resistance capacity to cope with coming climate change.

Arabinogalactan proteins and pectin distribution during female gametogenesis in Quercus suber L.

BACKGROUND AND AIMS: Quercus suber L. (cork oak) is one of the most important monoecious tree species in semi-arid regions of Southern Europe, with a high ecological value and economic potential. However, as a result of its long reproductive cycle, complex reproductive biology and recalcitrant seeds, conventional breeding is demanding. In its complex reproductive biology, little is known about the most important changes that occur during female gametogenesis. Arabinogalactan proteins (AGPs) and pectins are the main components of plant cell walls and have been reported to perform common functions in cell differentiation and organogenesis of reproductive plant structures. AGPs have been shown to serve as important molecules in several steps of the reproductive process in plants, working as signalling molecules, associated with the sporophyte-gametophyte transition, and pectins have been implicated in pollen-pistil interactions before double fertilization. In this study, the distribution of AGP and pectin epitopes was assessed during female gametogenesis. METHODS: Immunofluorescence labelling of female flower cells was performed with a set of monoclonal antibodies (mAbs) directed to the carbohydrate moiety of AGPs (JIM8 and JIM13) and pectic homogalacturonans (HGs) (mAbs JIM5 and JIM7). KEY RESULTS: The selective labelling obtained with AGP and pectin mAbs JIM8, JIM13, JIM5 and JIM7 during Q. suber female gametogenesis shows that AGPs and pectic HG can work as markers for mapping gametophytic cell differentiation in this species. Pectic HG showed different distribution patterns, depending on their levels of methyl esterification. Methyl-esterified HGs showed a uniform distribution in the overall female flower cells before fertilization and a more specific pattern after fertilization. A low methyl-ester pectin distribution pattern during the different developmental stages appears to be related to the pathway that pollen tubes follow to reach the embryo sac. AGPs showed a more sparse distribution in early stages of development, but specific labelling is shown in the synergids and their filiform apparatus. CONCLUSIONS: The labelling obtained with anti-AGP and anti-pectin mAbs in Q. suber female flower cells showed a dynamic distribution of AGPs and pectic HGs, which may render these molecules useful molecular markers during female gametogenesis. Changes occurring during development will be determined in order to help describe cork oak ovule structural properties before and after fertilization, providing new insight to better understand Q. suber female gametogenesis.

Can behavioural differences in Platypus cylindrus (Coleoptera: Platypodinae) from Portugal and Tunisia be explained by genetic and morphological traits?

Platypus cylindrus is an important wood borer of cork oak trees (Quercus suber) in the Mediterranean region, namely Portugal, Morocco and Algeria where its presence has drastically increased in the past few decades. On the contrary, the insect is not a relevant pest in Tunisia. The aim of this work is to analyze morphological and genetic differences among Tunisian and Portuguese populations in order to understand their role in the diverse population dynamics (e.g., aggressiveness) of the insect. The information could be used as a novel tool to implement protective measures. Insects were collected from cork oak stands in Tunisia (Ain Beya, Babouch and Mzara) and Portugal (Chamusca and Crato). Morphological traits of female and male mycangial pits were determined, using scanning electron microscopy but no significant differences were found. Genetic differences were analyzed using nuclear (internal simple sequence repeat polymerase chain reaction) and mitochondrial (cytochrome oxidase I (COI)) molecular markers. The results showed a very low level of intraspecific polymorphism and genetic diversity. The alignment of COI sequences showed high percentage of identical sites (99%) indicating a very low variation in nucleotide composition. Other variables related with the ecology of the insect and its associated fungi must be studied for a better understanding of the differences in the insect population's dynamic in Mediterranean countries.

Temporal dynamics of spectral bioindicators evidence biological and ecological differences among functional types in a cork oak open woodland.

The application of spectral vegetation indices for the purpose of vegetation monitoring and modeling increased largely in recent years. Nonetheless, the interpretation of biophysical properties of vegetation through their spectral signature is still a challenging task. This is particularly true in Mediterranean oak forest characterized by a high spatial and temporal heterogeneity. In this study, the temporal dynamics of vegetation indices expected to be related with green biomass and photosynthetic efficiency were compared for the canopy of trees, the herbaceous layer, and two shrub species: cistus (Cistus salviifolius) and ulex (Ulex airensis). coexisting in a cork oak woodland. All indices were calculated from in situ measurements with a FieldSpec3 spectroradiometer (ASD Inc., Boulder, USA). Large differences emerged in the temporal trends and in the correlation between climate and vegetation indices. The relationship between spectral indices and temperature, radiation, and vapor pressure deficit for cork oak was opposite to that observed for the herbaceous layer and cistus. No correlation was observed between rainfall and vegetation indices in cork oak and ulex, but in the herbaceous layer and in the cistus, significant correlations were found. The analysis of spectral vegetation indices with fraction of absorbed PAR (fPAR) and quantum yield of chlorophyll fluorescence (ΔF/Fm') evidenced strongest relationships with the indices Normalized Difference Water Index (NDWI) and Photochemical Reflectance Index (PRI)512, respectively. Our results, while confirms the ability of spectral vegetation indices to represent temporal dynamics of biophysical properties of vegetation, evidence the importance to consider ecosystem composition for a correct ecological interpretation of results when the spatial resolution of observations includes different plant functional types.

Evaluation of the presence of arabinogalactan proteins and pectins during Quercus suber male gametogenesis.

BACKGROUND AND AIMS: Quercus suber (cork oak) is a dominant tree of the Fagaceae in forests of the south-west Iberian Peninsula. It is monoecious with a long progamic phase that provides a comprehensive system for comparative studies in development and sexual reproduction. In this study the distribution of arabinogalactan protein (AGPs) and pectin epitopes in anthers of Q. suber was assessed to map these hydroxyproline-rich glycoproteins and the galacturonate-rich acidic polysaccharides during pollen development. Methods Immunolocalization in male flowers was performed with a set of monoclonal antibodies directed against the carbohydrate moiety that recognizes AGPs and pectins. To identify AGP genes involved in cork oak male flower development, a search was conducted for annotated AGP genes in the available transcriptome data of the Cork Oak EST Consortium database (www.corkoakdb.org). KEY RESULTS: Ubiquitous labelling in all cell types was obtained with anti-homogalacturan antibodies for methyl-esterified pectins. In contrast, the antibody that labelled non-methyl-esterified homogalacturans had a preferential presence in microsporocyte cells walls at the beginning of pollen development. Intense labelling was obtained with anti-AGP antibodies both in the tapetum and in the intine wall near the pollen apertures and later in the generative cell wall and vegetative cell. Evaluation of the putative AGPs highly expressed in the male gametophyte was achieved by quantitative RT-PCR analysis in male and female cork oak flowers. CONCLUSIONS: Four putative AGP genes were identified that are preferentially expressed in the male flower compared with the female flower. The putative Arabidopsis thaliana orthologues of these genes are associated with preferential expression in pollen, suggesting that the AGPs probably play a significant role in cork oak reproduction.

Climatic drivers of cork growth depend on site aridity.

Cork is one of the main non-timber forest products in the world. Most of its production is concentrated in the Iberian Peninsula, a climate change hotspot. Climate warming may lead to increased aridification and reduce cork production in that region. However, we still lack assessments of climate-cork relationships across ample geographical and climatic gradients explicitly considering site aridity. We quantified cork growth by measuring cork ring width and related it to climate variables and a drought index using dendrochronology. Four cork oak (Quercus suber) forests located from north eastern Spain to south western Morocco (31.5-41.5° N) and subjected to different aridity levels were sampled. Warm conditions in spring to early summer, when cork is formed, reduced cork width, whereas high precipitation in winter and spring enhanced it. The response of cork to increased water availability in summer peaked (r = 0.89, p = 0.00002) in the most arid and continental site considering 14-month long droughts. A severe drought caused a disproportionate loss of cork production in this site, where for every five-fold decrease in the drought index, the cork-width index declined by a factor of thirteen. Therefore, site aridity determines the responses of cork growth to the soil water availability resulting from accumulated precipitation during winter and spring previous to cork growth and until summer. In general, this cumulative water balance, which is very dependent on temperature and evapotranspiration rate, is critical for cork production, especially in continental, dry sites. The precipitation during the hydrological year can be used as a proxy of cork production in similar sites. Assessments of climate-cork relationships in the western Mediterranean basin could be used as analogues to forecast the impacts of aridification on future cork production.