Bacteriome dataset from the rhizosphere of trees in a Pinus pinaster and Pinus halepensis dominated forest subjected to drought conditions.

The Mediterranean basin is drastically affected by intense and frequent droughts, which jeopardize the diversity and survival of its forest, for example, Pinus pinaster forests. The dynamics of the bacterial communities inhabiting the rhizosphere of Pinus pinaster and other plants from a pine dominated forest under contrasting hydric conditions was monitored. The forest was located in Sierra de Oria (southeast Spain), and it was mainly composed by P. pinaster, P. halepensis, woody shrub species and herbaceous plants. 18 trees visually belonging to P. pinaster located along the perimeter and across the forest were selected for the analysis. All the trees were separated at least 50 m each other. Although all of them belonged to P. pinaster morphologically according to visual identification, the genotyping of the roots confirmed that they corresponded to P. pinaster, P. halepensis, and other plant species different from genus Pinus, although in the last case it was not possible to identify the plant species. At a distance less than 50 cm from the trunk, the litter and topsoil were removed, and the soil closely attached to non-suberified roots (rhizosphere soil) was collected (depth of 5-25 cm). Sampling was carried out in two seasons with contrasting temperature and rainfall patterns: on July 18, 2017 (summer) and April 24, 2018 (spring). After rhizosphere soil DNA and RNA extraction (and cDNA synthesis), a metabarcoding approach was followed by sequencing the V3-V4 hypervariable regions of the 16S rRNA gene and its derived transcripts by Illumina MiSeq platform. Sequencing reads were bioinformatically processed; specifically, they were filtered, trimmed, clustered into ASV (Amplicon Sequence Variants), and taxonomically identified. As a result, a total of 1,123,209 and 1,089,359 quality sequences were obtained from DNA and RNA-derived libraries, which resulted in 5,241 and 5,231 ASVs, respectively. Total communities (DNA) were mainly dominated by phyla Proteobacteria, Acidobacteria, Actinobacteria, Verrucomicrobia and Bacteroidetes in summer and spring, while potentially active populations (RNA libraries) were rich in Proteobacteria, Acidobacteria, Candidate division WPS-1, Actinobacteria and Verrucomicrobia both in summer and spring. On the other hand, DNA libraries were mainly dominated by genera Sphingomonas and acidobacterial groups Gp4 and Gp6, while potentially active bacteria (RNA) were rich in acidobacterial Gp3, Gp4, Gp6 and Phenylobacterium, although their relative abundance depended on the considered season. This dataset can provide valuable information about bacterial candidates which could be used as bioindicators of drought conditions. In addition to shifts in the bacterial relative abundance due to seasonal changes, the ratio RNA-based cDNA:DNA could be calculated as proxy of the potential activity of bacterial taxa. Moreover, these data could aid in developing bioformulations based on microorganisms which could be resistant to desiccation and involved in the drought resistance mechanisms of the host plant.

Correlating the above- and belowground genotype of Pinus pinaster trees and rhizosphere bacterial communities under drought conditions.

Increasing temperatures along with severe droughts are factors that may jeopardize the survival of the forests in the Mediterranean basin. In this region, Pinus pinaster is a common conifer species, that has been used as a model species in evolutionary studies due to its adaptive response to changing environments. Although its drought tolerance mechanisms are already known, knowledge about the dynamics of its root microbiota is still scarce. We aimed to decipher the structural (bacterial abundance), compositional, functional and associative changes of the P. pinaster rhizosphere bacterial communities in spring and summer, at DNA and RNA level (environmental DNA, live and dead cells, and those synthesizing proteins). A fundamental aspect of root microbiome-based approaches is to guarantee the correct origin of the samples. Thus, we assessed the genotype of host needles and roots from which rhizosphere samples were obtained. For more than 50% of the selected trees, genotype discrepancies were found and in three cases the plant species could not be determined. Rhizosphere bacterial communities were homogeneous with respect to diversity and structural levels regardless of the host genotype in both seasons. Nonetheless, significant changes were seen in the taxonomic profiles depending on the season. Seasonal changes were also evident in the bacterial co-occurrence patterns, both in DNA and RNA libraries. While spring communities switched to more complex networks, summer populations resulted in more compartmentalized networks, suggesting that these communities were facing a disturbance. These results may mirror the future status of bacterial communities in a context of climate change. A keystone hub was ascribed to the genus Phenylobacterium in the functional network calculated for summer. Overall, it is important to validate the origin and identity of plant samples in any plant-microbiota study so that more reliable ecological analyses are performed.

Comprehensive analysis of the isomiRome in the vegetative organs of the conifer Pinus pinaster under contrasting water availability.

An increasing number of microRNAs (miRNAs) and miRNA-related sequences produced during miRNA biogenesis, comprising the isomiRome, have been recently highlighted in different species as critical mediators of environmental stress responses. Conifers have some of the largest known genomes but an extensive characterization of the isomiRome from any conifer species has been lacking. We provide here a comprehensive overview of the Pinus pinaster isomiRome expressed in roots, stem and needles under well-watered and drought conditions. From the 13,441 unique small RNA sequences identified, 2,980 were annotated as canonical miRNAs or miRNA* and the remaining were classified as isomiRNA or miRNA-like sequences. A survey of their expression patterns highlighted roots as the most responsive organ under drought, where specific sequences of which a 24-nt novel miRNA stood out, were strongly down-regulated. Given the putative roles of the miRNA-targeted transcripts validated specifically in root tissues, some of the miRNAs, conserved and novel, are shortlisted as potential regulators of drought response. These results provide a valuable resource for comparative studies between gymnosperms and angiosperms. Furthermore, it evidences high transferability of the isomiRome between pine species being a useful basis for further molecular regulation and physiological studies, and especially those focused on adaptation to drought conditions.

Molecular study of drought response in the Mediterranean conifer Pinus pinaster Ait.: Differential transcriptomic profiling reveals constitutive water deficit-independent drought tolerance mechanisms.

Adaptation of long-living forest trees to respond to environmental changes is essential to secure their performance under adverse conditions. Water deficit is one of the most significant stress factors determining tree growth and survival. Maritime pine (Pinus pinaster Ait.), the main source of softwood in southwestern Europe, is subjected to recurrent drought periods which, according to climate change predictions for the years to come, will progressively increase in the Mediterranean region. The mechanisms regulating pine adaptive responses to environment are still largely unknown. The aim of this work was to go a step further in understanding the molecular mechanisms underlying maritime pine response to water stress and drought tolerance at the whole plant level. A global transcriptomic profiling of roots, stems, and needles was conducted to analyze the performance of siblings showing contrasted responses to water deficit from an ad hoc designed full-sib family. Although P. pinaster is considered a recalcitrant species for vegetative propagation in adult phase, the analysis was conducted using vegetatively propagated trees exposed to two treatments: well-watered and moderate water stress. The comparative analyses led us to identify organ-specific genes, constitutively expressed as well as differentially expressed when comparing control versus water stress conditions, in drought-sensitive and drought-tolerant genotypes. Different response strategies can point out, with tolerant individuals being pre-adapted for coping with drought by constitutively expressing stress-related genes that are detected only in latter stages on sensitive individuals subjected to drought.

Single-Copy Genes as Molecular Markers for Phylogenomic Studies in Seed Plants.

Phylogenetic relationships among seed plant taxa, especially within the gymnosperms, remain contested. In contrast to angiosperms, for which several genomic, transcriptomic and phylogenetic resources are available, there are few, if any, molecular markers that allow broad comparisons among gymnosperm species. With few gymnosperm genomes available, recently obtained transcriptomes in gymnosperms are a great addition to identifying single-copy gene families as molecular markers for phylogenomic analysis in seed plants. Taking advantage of an increasing number of available genomes and transcriptomes, we identified single-copy genes in a broad collection of seed plants and used these to infer phylogenetic relationships between major seed plant taxa. This study aims at extending the current phylogenetic toolkit for seed plants, assessing its ability for resolving seed plant phylogeny, and discussing potential factors affecting phylogenetic reconstruction. In total, we identified 3,072 single-copy genes in 31 gymnosperms and 2,156 single-copy genes in 34 angiosperms. All studied seed plants shared 1,469 single-copy genes, which are generally involved in functions like DNA metabolism, cell cycle, and photosynthesis. A selected set of 106 single-copy genes provided good resolution for the seed plant phylogeny except for gnetophytes. Although some of our analyses support a sister relationship between gnetophytes and other gymnosperms, phylogenetic trees from concatenated alignments without 3rd codon positions and amino acid alignments under the CAT + GTR model, support gnetophytes as a sister group to Pinaceae. Our phylogenomic analyses demonstrate that, in general, single-copy genes can uncover both recent and deep divergences of seed plant phylogeny.

Analysis of DNA Cytosine Methylation Patterns Using Methylation-Sensitive Amplification Polymorphism (MSAP).

Different molecular techniques have been developed to study either the global level of methylated cytosines or methylation at specific gene sequences. One of them is the methylation-sensitive amplified polymorphism technique (MSAP) which is a modification of amplified fragment length polymorphism (AFLP). It has been used to study methylation of anonymous CCGG sequences in different fungi, plants, and animal species. The main variation of this technique resides on the use of isoschizomers with different methylation sensitivity (such as HpaII and MspI) as a frequent-cutter restriction enzyme. For each sample, MSAP analysis is performed using both EcoRI/HpaII- and EcoRI/MspI-digested samples. A comparative analysis between EcoRI/HpaII and EcoRI/MspI fragment patterns allows the identification of two types of polymorphisms: (1) methylation-insensitive polymorphisms that show common EcoRI/HpaII and EcoRI/MspI patterns but are detected as polymorphic amplified fragments among samples and (2) methylation-sensitive polymorphisms which are associated with the amplified fragments that differ in their presence or absence or in their intensity between EcoRI/HpaII and EcoRI/MspI patterns. This chapter describes a detailed protocol of this technique and discusses the modifications that can be applied to adjust the technology to different species of interest.

Organ-specific metabolic responses to drought in Pinus pinaster Ait.

Drought is an important driver of plant survival, growth, and distribution. Water deficit affects different pathways of metabolism, depending on plant organ. While previous studies have mainly focused on the metabolic drought response of a single organ, analysis of metabolic differences between organs is essential to achieve an integrated understanding of the whole plant response. In this work, untargeted metabolic profiling was used to examine the response of roots, stems, adult and juvenile needles from Pinus pinaster Ait. full-sib individuals, subjected to a moderate and long lasting drought period. Cyclitols content showed a significant alteration, in response to drought in all organs examined, but other metabolites increased or decreased differentially depending on the analyzed organ. While a high number of flavonoids were only detected in aerial organs, an induction of the glutathione pathway was mainly detected in roots. This result may reflect different antioxidant mechanisms activated in aerial organs and roots. Metabolic changes were more remarkable in roots than in the other organs, highlighting its prominent role in the response to water stress. Significant changes in flavonoids and ascorbate metabolism were also observed between adult and juvenile needles, consistent with previously proven differential functional responses between the two developmental stages. Genetic polymorphisms in candidate genes coding for a Myb1 transcription factor and a malate dehydrogenase (EC 1.1.1.37) were associated with different concentration of phenylalanine, phenylpropanoids and malate, respectively. The results obtained will support further research on metabolites and genes potentially involved in functional mechanisms related to drought tolerance in trees.

Evidence of intense chromosomal shuffling during conifer evolution.

While recent advances have been gained on genome evolution in angiosperm lineages, virtually nothing is known about karyotype evolution in the other group of seed plants, the gymnosperms. Here we used high density gene-based linkage mapping to compare the karyotype structure of two families of conifers (the most abundant group of gymnosperms) separated around 290 million years ago: Pinaceae and Cupressaceae. We propose for the first time a model based on the fusion of 20 ancestral chromosomal blocks that may have shaped the modern karyotpes of Pinaceae (with n=12) and Cupressaceae (with n=11). The considerable difference in modern genome organization between these two lineages contrasts strongly with the remarkable level of synteny already reported within the Pinaceae. It also suggests a convergent evolutionary mechanism of chromosomal block shuffling that has shaped the genomes of the spermatophytes.

Massive sequencing of Ulmus minor's transcriptome provides new molecular tools for a genus under the constant threat of Dutch elm disease.

Elms, especially Ulmus minor and U. americana, are carrying out a hard battle against Dutch elm disease (DED). This vascular wilt disease, caused by Ophiostoma ulmi and O. novo-ulmi, appeared in the twentieth century and killed millions of elms across North America and Europe. Elm breeding and conservation programmes have identified a reduced number of DED tolerant genotypes. In this study, three U. minor genotypes with contrasted levels of tolerance to DED were exposed to several biotic and abiotic stresses in order to (i) obtain a de novo assembled transcriptome of U. minor using 454 pyrosequencing, (ii) perform a functional annotation of the assembled transcriptome, (iii) identify genes potentially involved in the molecular response to environmental stress, and (iv) develop gene-based markers to support breeding programmes. A total of 58,429 putative unigenes were identified after assembly and filtering of the transcriptome. 32,152 of these unigenes showed homology with proteins identified in the genome from the most common plant model species. Well-known family proteins and transcription factors involved in abiotic, biotic or both stresses were identified after functional annotation. A total of 30,693 polymorphisms were identified in 7,125 isotigs, a large number of them corresponding to single nucleotide polymorphisms (SNPs; 27,359). In a subset randomly selected for validation, 87% of the SNPs were confirmed. The material generated may be valuable for future Ulmus gene expression, population genomics and association genetics studies, especially taking into account the scarce molecular information available for this genus and the great impact that DED has on elm populations.

Nucleotide polymorphisms in a pine ortholog of the Arabidopsis degrading enzyme cellulase KORRIGAN are associated with early growth performance in Pinus pinaster.

We have carried out a candidate-gene-based association genetic study in Pinus pinaster Aiton and evaluated the predictive performance for genetic merit gain of the most significantly associated genes and single nucleotide polymorphisms (SNPs). We used a second generation 384-SNP array enriched with candidate genes for growth and wood properties to genotype mother trees collected in 20 natural populations covering most of the European distribution of the species. Phenotypic data for total height, polycyclism, root-collar diameter and biomass were obtained from a replicated provenance-progeny trial located in two sites with contrasting environments (Atlantic vs Mediterranean climate). General linear models identified strong associations between growth traits (total height and polycyclism) and four SNPs from the korrigan candidate gene, after multiple testing corrections using false discovery rate. The combined genomic breeding value predictions assessed for the four associated korrigan SNPs by ridge regression-best linear unbiased prediction (RR-BLUP) and cross-validation accounted for up to 8 and 15% of the phenotypic variance for height and polycyclic growth, respectively, and did not improve adding SNPs from other growth-related candidate genes. For root-collar diameter and total biomass, they accounted for 1.6 and 1.1% of the phenotypic variance, respectively, but increased to 15 and 4.1% when other SNPs from lp3.1, lp3.3 and cad were included in RR-BLUP models. These results point towards a desirable integration of candidate-gene studies as a means to pre-select relevant markers, and aid genomic selection in maritime pine breeding programs.

Epigenetic variability in the genetically uniform forest tree species Pinus pinea L.

There is an increasing interest in understanding the role of epigenetic variability in forest species and how it may contribute to their rapid adaptation to changing environments. In this study we have conducted a genome-wide analysis of cytosine methylation pattern in Pinus pinea, a species characterized by very low levels of genetic variation and a remarkable degree of phenotypic plasticity. DNA methylation profiles of different vegetatively propagated trees from representative natural Spanish populations of P. pinea were analyzed with the Methylation Sensitive Amplified Polymorphism (MSAP) technique. A high degree of cytosine methylation was detected (64.36% of all scored DNA fragments). Furthermore, high levels of epigenetic variation were observed among the studied individuals. This high epigenetic variation found in P. pinea contrasted with the lack of genetic variation based on Amplified Fragment Length Polymorphism (AFLP) data. In this manner, variable epigenetic markers clearly discriminate individuals and differentiates two well represented populations while the lack of genetic variation revealed with the AFLP markers fail to differentiate at both, individual or population levels. In addition, the use of different replicated trees allowed identifying common polymorphic methylation sensitive MSAP markers among replicates of a given propagated tree. This set of MSAPs allowed discrimination of the 70% of the analyzed trees.

Genetic control of functional traits related to photosynthesis and water use efficiency in Pinus pinaster Ait. drought response: integration of genome annotation, allele association and QTL detection for candidate gene identification.

BACKGROUND: Understanding molecular mechanisms that control photosynthesis and water use efficiency in response to drought is crucial for plant species from dry areas. This study aimed to identify QTL for these traits in a Mediterranean conifer and tested their stability under drought. RESULTS: High density linkage maps for Pinus pinaster were used in the detection of QTL for photosynthesis and water use efficiency at three water irrigation regimes. A total of 28 significant and 27 suggestive QTL were found. QTL detected for photochemical traits accounted for the higher percentage of phenotypic variance. Functional annotation of genes within the QTL suggested 58 candidate genes for the analyzed traits. Allele association analysis in selected candidate genes showed three SNPs located in a MYB transcription factor that were significantly associated with efficiency of energy capture by open PSII reaction centers and specific leaf area. CONCLUSIONS: The integration of QTL mapping of functional traits, genome annotation and allele association yielded several candidate genes involved with molecular control of photosynthesis and water use efficiency in response to drought in a conifer species. The results obtained highlight the importance of maintaining the integrity of the photochemical machinery in P. pinaster drought response.

De novo assembly of maritime pine transcriptome: implications for forest breeding and biotechnology.

Maritime pine (Pinus pinasterAit.) is a widely distributed conifer species in Southwestern Europe and one of the most advanced models for conifer research. In the current work, comprehensive characterization of the maritime pine transcriptome was performed using a combination of two different next-generation sequencing platforms, 454 and Illumina. De novo assembly of the transcriptome provided a catalogue of 26 020 unique transcripts in maritime pine trees and a collection of 9641 full-length cDNAs. Quality of the transcriptome assembly was validated by RT-PCR amplification of selected transcripts for structural and regulatory genes. Transcription factors and enzyme-encoding transcripts were annotated. Furthermore, the available sequencing data permitted the identification of polymorphisms and the establishment of robust single nucleotide polymorphism (SNP) and simple-sequence repeat (SSR) databases for genotyping applications and integration of translational genomics in maritime pine breeding programmes. All our data are freely available at SustainpineDB, the P. pinaster expressional database. Results reported here on the maritime pine transcriptome represent a valuable resource for future basic and applied studies on this ecological and economically important pine species.

Molecular response to water stress in two contrasting Mediterranean pines (Pinus pinaster and Pinus pinea).

Adaptation to water stress has determined the evolution and diversification of vascular plants. Water stress is forecasted to increase drastically in the next decades in certain regions, such as in the Mediterranean basin. Consequently, a proper knowledge of the response and adaptations to drought stress is essential for the correct management of plant genetic resources. However, most of the advances in the understanding of the molecular response to water stress have been attained in angiosperms, and are not always applicable to gymnosperms. In this work we analyse the transcriptional response of two emblematic Mediterranean pines, Pinus pinaster and Pinus pinea, which show noticeable differences in their performance under water stress. Using microarray analysis, up to 113 genes have been detected as significantly induced by drought in both species. Reliability of expression patterns has been confirmed by RT-PCR. While induced genes with similar profiles in both species can be considered as general candidate genes for the study of drought response in conifers, genes with diverging expression patterns can underpin the differences displayed by these species under water stress. Most promising candidate genes for drought stress response include genes related to carbohydrate metabolism, such as glycosyltransferases or galactosidases, sugar transporters, dehydrins and transcription factors. Additionally, differences in the molecular response to drought and polyethylene-glycol-induced water stress are also discussed.

Epigenetic regulation of adaptive responses of forest tree species to the environment.

Epigenetic variation is likely to contribute to the phenotypic plasticity and adaptative capacity of plant species, and may be especially important for long-lived organisms with complex life cycles, including forest trees. Diverse environmental stresses and hybridization/polyploidization events can create reversible heritable epigenetic marks that can be transmitted to subsequent generations as a form of molecular "memory". Epigenetic changes might also contribute to the ability of plants to colonize or persist in variable environments. In this review, we provide an overview of recent data on epigenetic mechanisms involved in developmental processes and responses to environmental cues in plant, with a focus on forest tree species. We consider the possible role of forest tree epigenetics as a new source of adaptive traits in plant breeding, biotechnology, and ecosystem conservation under rapid climate change.

DNA sequence variation of drought-response candidate genes in Austrocedrus chilensis

Background: Austrocedrus chilensis (D. Don) Pic. Ser. et Bizzarri commonly known as Patagonian cypress is a member of the Cupressaceae family, characterized by a high adaptive potential for growing in marginal areas and good timber quality. The species grows over a wide area and under a wide range of rainfall. This study assessed adaptive genetic variation at SNP level in candidate genes involved in response to drought stress. Results: A total of 18 single nucleotide polymorphisms (SNPs) were found among 1,428 bp. Average nucleotide diversity value (π = 0.00312) was similar to those previously reported in other Cupressaceae. The Fst average among genes and populations was 0.163 and the lowest differentiation was observed in continuous and humid populations. A number of neutrality tests were applied to find evidence of positive selection in our candidate gene set, but only AcAQP2 gene in Pedregoso and San Ramón populations revealed significant departures from neutrality with positive values suggesting balancing selection. Conclusions: In this study we report the levels of nucleotide diversity searched in some drought stress candidate genes in Austrocedrus chilensis and the selective factors that may be acting on this species.

Annotated genetic linkage maps of Pinus pinaster Ait. from a Central Spain population using microsatellite and gene based markers.

BACKGROUND: Pinus pinaster Ait. is a major resin producing species in Spain. Genetic linkage mapping can facilitate marker-assisted selection (MAS) through the identification of Quantitative Trait Loci and selection of allelic variants of interest in breeding populations. In this study, we report annotated genetic linkage maps for two individuals (C14 and C15) belonging to a breeding program aiming to increase resin production. We use different types of DNA markers, including last-generation molecular markers. RESULTS: We obtained 13 and 14 linkage groups for C14 and C15 maps, respectively. A total of 211 and 215 markers were positioned on each map and estimated genome length was between 1,870 and 2,166 cM respectively, which represents near 65% of genome coverage. Comparative mapping with previously developed genetic linkage maps for P. pinaster based on about 60 common markers enabled aligning linkage groups to this reference map. The comparison of our annotated linkage maps and linkage maps reporting QTL information revealed 11 annotated SNPs in candidate genes that co-localized with previously reported QTLs for wood properties and water use efficiency. CONCLUSIONS: This study provides genetic linkage maps from a Spanish population that shows high levels of genetic divergence with French populations from which segregating progenies have been previously mapped. These genetic maps will be of interest to construct a reliable consensus linkage map for the species. The importance of developing functional genetic linkage maps is highlighted, especially when working with breeding populations for its future application in MAS for traits of interest.

Towards decoding the conifer giga-genome.

Several new initiatives have been launched recently to sequence conifer genomes including pines, spruces and Douglas-fir. Owing to the very large genome sizes ranging from 18 to 35 gigabases, sequencing even a single conifer genome had been considered unattainable until the recent throughput increases and cost reductions afforded by next generation sequencers. The purpose of this review is to describe the context for these new initiatives. A knowledge foundation has been acquired in several conifers of commercial and ecological interest through large-scale cDNA analyses, construction of genetic maps and gene mapping studies aiming to link phenotype and genotype. Exploratory sequencing in pines and spruces have pointed out some of the unique properties of these giga-genomes and suggested strategies that may be needed to extract value from their sequencing. The hope is that recent and pending developments in sequencing technology will contribute to rapidly filling the knowledge vacuum surrounding their structure, contents and evolution. Researchers are also making plans to use comparative analyses that will help to turn the data into a valuable resource for enhancing and protecting the world's conifer forests.

Identification of water stress genes in Pinus pinaster Ait. by controlled progressive stress and suppression-subtractive hybridization.

Climate change is a major challenge particularly for forest tree species, which will have to face the severe alterations of environmental conditions with their current genetic pool. Thus, an understanding of their adaptive responses is of the utmost interest. In this work we have selected Pinus pinaster as a model species. This pine is one of the most important conifers (for which molecular tools and knowledge are far more scarce than for angiosperms) in the Mediterranean Basin, which is characterised in all foreseen scenarios as one of the regions most drastically affected by climate change, mainly because of increasing temperature and, particularly, by increasing drought. We have induced a controlled, increasing water stress by adding PEG to a hydroponic culture. We have generated a subtractive library, with the aim of identifying the genes induced by this stress and have searched for the most reliable expressional candidate genes, based on their overexpression during water stress, as revealed by microarray analysis and confirmed by RT-PCR. We have selected a set of 67 candidate genes belonging to different functional groups that will be useful molecular tools for further studies on drought stress responses, adaptation, and population genomics in conifers, as well as in breeding programs.

Functional and genetic characterization of gas exchange and intrinsic water use efficiency in a full-sib family of Pinus pinaster Ait. in response to drought.

Drought is an important environmental factor in Mediterranean ecosystems affecting seedling recruitment, productivity or susceptibility to fires and pathogens. Studying water use efficiency in these environments is crucial due to its adaptive value allowing trees to cope with low water availability. We studied the phenotypic variability and genetic control of intrinsic water use efficiency (WUE(i)) and related traits in a full-sib family of Pinus pinaster under drought imposition. We detected significant differences in WUE(i) between clones of the same family and moderate heritability estimates that indicate some degree of genetic control over this trait. Stomatal conductance to water vapor was the trait most affected by drought imposition and it showed the strongest influence in WUE(i). Stomatal conductance to water vapor and specific leaf area (SLA) were the traits with highest heritabilities and they showed a significant genetic correlation with WUE(i), suggesting that selection of needles with low SLA values will improve WUE(i) in this species by reducing water losses through stomatal control.