Skip to main content

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

Abstract

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.

Background

Drought resistance is crucial for growth and survival of species living in water scarce environments [1]. Unraveling the molecular mechanisms that control functional traits, such as photosynthesis and water use efficiency in response to drought, is especially relevant in view of its implication in survival, growth and biomass production. However, carbon uptake in response to drought is a complex process with many mechanisms acting in coordination in final CO2 fixation [2]. From stomatal and mesophyll resistances to diffusion of CO2 to biochemical processes within chloroplast, complex mechanisms are involved in net carbon fixation [25]. The functional bases that control carbon uptake under water stress have been largely studied [6, 7], but less information is available about its genetic regulation.

Complex functional trait dissection can be achieved through two approaches: association studies and QTL (Quantitative Trait Loci) mapping [8]. The resolution power of association studies is higher than QTL mapping [9]. However, the rapid decay in linkage disequilibrium of conifers [10] makes the development of genome wide association studies in these species laborious and advocates in favor of candidate gene approaches [11].

In the past, identification of candidate genes underlying QTL was difficult due to the use of anonymous markers and limited sequences with functional information, and thus restricted the approach to model plant species [12]. Nowadays, gene-based markers are easily developed and much more functional information is available for a wide range of organisms [1318], allowing to integrate functional annotation with QTL studies [19]. Moreover, the development and application of high throughput genotyping technologies have allowed the construction of dense genetic maps [2027], http://dendrome.ucdavis.edu/cmap/]. The use of highly saturated genetic maps allows to narrow down the position of loci involved in the genetic control of the targeted trait and the combination of high density gene based maps with functional annotation allows to identify positional candidate genes for these QTL [19, 28]. Suggested candidate genes are suitable for association studies that can validate marker-trait associations [29]. Therefore, identification of positional candidate genes within QTL confidence intervals, some of them with known function in other species, could be considered as a preliminary step that contributes to the detection of genes underlying traits of interest [30]. Additionally, QTL mapping allows the evaluation of the genetic basis for potential adaptation in natural populations [31, 32] and to extend the understanding of relationships between different morpho-functional traits [33]. The identification of the main QTL involved in drought response could be a first step to develop marker assisted selection (MAS) strategies for these traits [11].

Consequently, the detection of QTL involved in photosynthesis and water use efficiency in the context of drought response is a first attempt to understand the genetic basis regulating the expression of these traits. QTL studies on functional drought response have been largely implemented for non-forest model species [3439]. Some of these QTL studies in crop species have recently identified genomic regions controlling photochemistry of carbon uptake [40, 41]. Breeding programs implemented in crops have reported yield improvement associated with increased photosynthesis [42]. However, fewer QTL analyses on functional drought response of forest tree species have been performed [11, 4346] and to our knowledge none of them has focused on the photochemical machinery.

QTL studies involve development of a segregant progeny for target traits, phenotypic and molecular characterization of the progeny and construction of genetic maps [47]. The power to resolve the location of a QTL is related to the size of the studied population and the mapping coverage [48]. Additionally, forest tree species are characterized by long generation times which hinder development of backcross or three-generation pedigrees by controlled crosses. On the other hand, replication of each genotype is needed for a reliable phenotypic evaluation [49], especially when working with physiological parameters that are extremely sensible to environmental conditions [41, 50, 51].

Mediterranean species are particularly threatened by drought [5254], especially in the context of climatic change predictions [55]. Pinus pinaster Ait. is an important conifer in Mediterranean region with a high ecological and socio-economical value [5658]. Although P. pinaster shows evidence of drought adaptation [59, 60], recurrent or severe drought periods can limit its growth [61, 62]. Understanding the molecular basis of drought tolerance is of high importance for a suitable management of the available genetic resources of P. pinaster in conservation, afforestation or breeding programs. QTL and association studies of drought tolerance traits have been developed in several tree species, such as P. taeda[23, 63, 64], Populus sp. [19, 65, 66] or Quercus robur[43, 67]. Several QTL and association studies in P. pinaster have analyzed the molecular basis of different processes related to growth or wood quality traits [6873], terpenes [74] and serotiny [75]. However, to date only association studies based on a few potential candidate genes [59, 60] and one QTL study have analyzed the molecular basis of drought tolerance in P. pinaster[44].

The main objective of this work was to unravel the genetic basis of different functional parameters related to carbon uptake and water use efficiency in response to drought for P. pinaster. For this purpose a QTL analysis using vegetatively propagated genotypes in order to improve the reliability of phenotypic estimates was designed. Several specific objectives were outlined: 1) construction of dense gene-based linkage maps with functional information; 2) identification of genomic regions underlying photosynthesis and water use efficiency in response to drought through QTL analysis; and 3) identification of a set of promising candidate genes in targeted genomic regions that may be involved in the genetic regulation of photosynthesis and water use efficiency in response to drought.

Methods

Plant material, experimental setup and phenotypic evaluation

Plant material, experimental setup and phenotypic evaluation are explained in detail in de Miguel et al. [76]. Briefly, 162 seedlings from a F1 full-sib family of P. pinaster obtained from a controlled cross between a male parent (Oria6) from Oria, a natural population from South-East Spain (37° 31 ’N 2° 21 ’W) and a female parent (Gal1056) from a breeding program established in Pontevedra, North-West Spain (42° 10 ’N 8° 30 ’W), were vegetatively replicated and established in an incomplete block design in a greenhouse at Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA). Phenotypic evaluations were conducted on the 103 clones for which at least three ramets were obtained. For phenotypic characterization three time-points of measurement were carried out starting in October 2009. During the 1st time-point, plants were watered close to full holding capacity. Then, watering was withdrawn and during the 2nd time-point of measurement plants were left 7 days without watering. The final third batch of measurements was carried out after plants have been 14 days without watering. Net photosynthetic rate (An, μmol CO2m-2 s-1), stomatal conductance to water vapour (gsw, molH2Om-2 s-1), intrinsic water use efficiency (WUEi, μmol CO2 molH2O-1), specific leaf area (SLA, m2Kg-1), maximum efficiency of photosystem II under light conditions (Fv’Fm’) and quantum yield (ΦPSII) were measured for all plants. Chlorophyll fluorescence parameters were measured following the procedure described in Cano et al. [6].

In the 1st time-point of measurement four adult needles were collected for each plant, dried and ground into a fine homogeneous powder. Carbon isotope composition was measured with a PDZ Europa ANCA-GSL elemental analyzer interfaced to a PDZ Europa 20–20 continuous flow isotope ratio mass spectrometer (Sercon Ltd., Cheshire, UK) at Stable Isotope Facility UC Davis, California, USA. The isotopic composition of 13C (‰) was expressed as [77]:

δ 13 C = R s - R b R b * 1000

Where Rs and Rb refer to the 13C/12C ratio in the sample and in the Pee Dee Belemnite standard, respectively.

Broad-sense heritability estimates and genetic correlations were calculated for the analyzed traits according to de Miguel et al. [76].

DNA extraction and marker genotyping

The mapping progeny was genotyped with nuclear microsatellites (single sequence repeats, nSSR), selective amplification of microsatellite polymorphic loci (SAMPL) and single nucleotide polymorphism (SNP) markers. Different DNA extraction methods were used in needles: a modified protocol from Dellaporta et al. [78] for nSSRs, SAMPLs and SNP array D (detailed below); the commercial kit Invisorb DNA plants HTS 96 kit (Invitek GmbH, Berlin, Germany) for SNP arrays A and C (detailed below) and the commercial kit DNeasy Plant mini kit (Qiagen, Düsseldorf, Germany) for SNP array B (detailed below).

A total of twenty nine primer pairs designed for amplification of nSSR loci in P. pinaster and P. taeda[79, 80] were tested for their segregations in the mapping population and both progenitors and six progeny individuals were genotyped. The whole mapping progeny was then genotyped only for polymorphic microsatellite loci. PCRs were performed in 10 μl containing 10 ng of DNA, 1x PCR reaction buffer (Invitrogen, Grand Island, NY, USA), 250 μM of each dNTP (Invitrogen, Grand Island, NY, USA), 0.25 U Taq polymerase (Invitrogen, Grand Island, NY, USA), 4 mM MgCl2 (Invitrogen, Grand Island, NY, USA) except for A6D12 where 2 mM MgCl2 was used, 0.2 μM of forward primer and 0.2 μM of reverse primer labeled on its 5’ end with IRD800. The PCR profile used was 94°C 4 min, 2 cycles of 94°C 45 s, 60°C 45 s, 72°C 45 s, 18 touchdown cycles of 94°C 45 s, 59.5°C 45 s (-0.5°C/cycle), 72°C 45 s, 20 cycles 94°C 30 s, 50°C 30 s, 72°C 45 s and final extension at 72°C 5 min. PCR reactions were carried out with a Perkin-Elmer GenAmp 9700 thermal cycler (Perkin Elmer Inc., Waltham, Massachusetts, USA). Amplified products were separated on denaturing gels containing 6% (w/v) acrylamide/bisacrylamide (19:3), 7 M urea and 1 x TBE and visualized in a 4300 DNA Analyzer (LI-COR Biosciences, Lincoln, NE, USA). Fragments were scored visually as codominant markers.

SAMPL genotyping was performed as in de Miguel et al. [81] with several modifications. Preamplification was carried out using EcoRI + A / MseI + C primer combination. In order to identify the most informative selective primer combinations (those with a higher number of informative polymorphic fragments) different primer combinations were tested using DNA from the progenitors and 6 offspring. A total of five CATA/EcoRI and three GATA/EcoRI primer combinations were used for the selective amplification. The whole mapping progeny was then genotyped for the eight selected SAMPL primer combinations. Primers CATA and GATA were IRDye 700 and IRDye 800 5’end labeled, respectively. Samples were loaded into denaturing gels containing 16% (w/v) Long Ranger® 50% (w/v) Gel Solution (Lonza, Basel, Switzerland), 7 M urea and 1 x TBE. Fragment detection was carried out on a 4300 DNA Analyzer (LI-COR Biosciences, Lincoln, NE, USA). Each gel was visually scored twice independently by two different people.

In this study, four SNP genotyping assays were used, three of which were Golden Gate assays (Illumina Inc., San Diego, CA, USA): SNP arrays A and C, which were two different 1,536 BeadArray™ experiments; and SNP array B, which was a 384 BeadXpress®. The SNP array D was a 12 K Infinium assay (Illumina Inc., San Diego, CA, USA). SNP arrays B and C were used to genotype the whole mapping progeny, whereas A and D could be used only on 83 and 70 progeny individuals, respectively. SNP array B was developed including many of the SNPs targeted in array A [82] and 14 additional SNPs from candidate genes for drought resistance [59] in order to complete the information for a set of genes of special interest (see de Miguel et al. [81] for further details). When the same SNP was successfully genotyped in both assays only the data of SNP array B was used because of the higher number of individuals genotyped in this assay. SNP array C was designed using a P. pinaster gene catalog obtained from 454 sequencing of cDNA libraries constructed with different tissues from 9 siblings of the mapping progeny submitted to different growing conditions (i.e. drought stress versus control plants; E Sáez-Laguna et al., unpublished). The genotyping of SNP array C was developed at CNIO, Madrid, Spain. Finally, SNP array D contained 10,593 SNPs identified from unigene set “PineContig_v2” of P. pinaster[20]. The four genotyping assays were carried out according to the manufacturer’s instructions (Illumina Inc., San Diego, CA, USA) and SNPs clusters revised manually with Illumina Genome Studio v.1.9.4 software with a GenCall score cutoff of 0.15. SNP clusters were modified manually to refine cluster positions when necessary. For the SNP array D (12 K Infinium) SNPs with Gen-Train values lower than 0.25 were discarded, with values between 0.25 and 0.5 were manually scored and with values higher than 0.5 were automatically scored.

Construction of dense linkage maps

For the construction of two genetic maps, one for each progenitor (Gal1056 and Oria6), the “two-way-pseudo-testcross” mapping strategy was applied [83]. The consensus map for the cross, combining markers informative for both parents, was also developed (GxO). Linkage analyses and map estimations were performed using the regression mapping algorithm implemented in the software JoinMap® v4.1 [84] with the CP population type and using a recombination fraction < 0.35 and a LOD > 3 as mapping parameters. Map distances were calculated using Kosambi mapping function [85]. For map building a goodness-of-fit jump threshold of 5 was established. JoinMap suggests three genetic maps with increasing number of markers (map1, map2 and map3). In map2, new markers were added because more pair wise data were available but statistical support is the same as in map1. In map3, the remaining loci were added by increasing the goodness-of fit jump threshold. In these cases map2 was kept for further analyses. Mean χ2 contribution to the goodness of fit and number of double recombinants were inspected in order to remove not reliably positioned markers from the estimated maps. When a pair of markers was considered identical based on the lack of recombination between them, only one of the markers was selected for mapping (see Additional file 1). Segregation ratios were tested using χ2 test (p ≤ 0.01) after Bonferroni correction. Framework maps for Gal1056, Oria6 and GxO were also built. For this purpose, only the most informative markers with very reliable positions and inter marker distance of circa 10 cM were kept. Total genome length was calculated as the sum of all mapped marker intervals. Estimated genome length was determined from the partial linkage data according to Hulbert et al. [86] modified by Chakravarti et al. [87] (Method 3). To estimate genome length using framework maps, a minimum LOD score of three was chosen. Observed map coverage was calculated as the ratio of total genome length to estimated genome length. To estimate the number of different mapped genes a BlastN was performed between gene sequences contained in the different SNP genotyping arrays. Sequences with a percentage of identity higher than 98% were considered the same gene. To test whether the mapped genes were evenly distributed between linkage groups χ2 tests (p < 0.05) were performed by comparing observed and estimated numbers of genes per linkage group (LG). The expected number of genes for each LG was obtained by multiplying the ratio size of LG to total genome length by the total number of mapped genes. Linkage maps were compared with previously developed P. pinaster maps [20, 81, 82] based on common SNPs and SSRs.

QTL mapping

In order to avoid errors in marker order that may have some impact on the precision and accuracy of QTL placement, QTL analyses were performed using the framework linkage maps. QTL detection was carried out using the regression algorithm implemented in the software MapQTL® v6.0 [88]. Interval mapping was applied followed by multiple QTL mapping (MQM) when more than one QTL was found for a trait. Analyses were performed using a mapping step size of one. The thresholds (95% and 99% confidence) for QTL significance were determined using a chromosome and genome wide permutation test with 10,000 iterations. Support intervals for the detected QTL were estimated based on the observed decrease of LOD value in one and two units. QTL identified with only the 95% significance at chromosome level were considered as suggestive of putative QTL. Each detected QTL received an identification name indicating the measured trait, the time-point of measurement, the linkage group (LG) and the map (“f” and “m” for female and male parents respectively and “i” for consensus map) where the QTL was detected.

Candidate genes search

Functional annotation for gene based markers of SNP arrays A, B and D were described by Chancerel et al. [20, 82]. Functional annotation for SNP array C was obtained with Blast2GO software [89]. For high-scoring segment pair (HSP) a restrictive E-20 e-value was chosen finding in a Blastx search against a set of 88,516 reference proteins from UniProt (http://www.uniprot.org/). In order to update annotation information for the sequences of the four SNP arrays that mapped to the QTL, a second round of annotation was performed using Blast2GO software with a restrictive E-25 e-value for HSP, and annotation was completed with InterPro (http://www.ebi.ac.uk/interpro/) and Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/) searches. For those genes within the significant QTL confidence intervals (±2 LOD), functional annotations were queried to identify functional relationship between the positional candidate genes and each analyzed trait. In order to cover all mapped genes, QTL confidence intervals in framework maps were extrapolated to maps with all the mapped markers.

Association between phenotypes and alleles at candidate loci was further studied by ANOVA using the traits as dependent variables and the SNP genotypes as factors. Thereof nineteen traits (seven different traits measured at three water irrigation regimes, except δ13C measured only at 1st time-point) and 73 SNPs located in 58 identified candidate genes were inspected. False discovery rate (FDR) was calculated using the package qvalue. Association analyses were carried out in R version 2.15.2 (R Development Core Team, 2012).

Results

Phenotypic evaluation

Descriptive statistics of all analyzed traits are shown in Table 1. Almost all traits showed a close to normal distribution with low levels of skewness and kurtosis. Although normal distribution is an assumption in interval mapping, this method and MQM are quiet robust against deviations from normality [88]. Water stress produced a decrease in mean values for almost all variables except for WUEi and SLA that showed higher and very similar mean values, respectively, for the three time-points of measurement. Coefficients of variation were progressively higher with the imposition of drought stress being gsw the trait that showed the higher coefficient of variation in the 3rd time-point of measurement.

Table 1 Descriptive statistics of measured traits in the F 1 full sib family Gal1056xOria6 (n = 103)

Phenotypic correlations between the studied traits are presented in Table 2. An was correlated with gsw and with chlorophyll fluorescence parameters (Fv’Fm’ and ΦPSII). The magnitude of the correlation coefficients was very similar for the 1st time-point of measurement. However, under drought stress An showed a higher correlation coefficient with chlorophyll fluorescence parameters than with gsw. Besides, An and chlorophyll fluorescence parameters showed a tight genetic correlation (see Additional file 2). For WUEi and δ13C, a significant phenotypic (Table 2) and broad sense genetic correlation (see Additional file 2) was found. Both traits had higher phenotypic correlation coefficients with gsw than with An. SLA was moderately correlated with An, gsw, WUEi and δ13C for the 1st and 2nd time-points of measurement (Table 2). Broad sense heritability estimates for the analyzed traits are presented in Additional file 3. All of them presented moderate to low values of heritability being the higher estimates for gsw, WUEi and δ13C.

Table 2 Pearson correlation coefficients and statistical significance for measured traits in the F 1 full sib family Gal1056xOria6 (n = 103)

Highly saturated linkage maps

For Gal1056, Oria6 and consensus map, 17, 16 and 13 linkage groups (LG) were obtained, respectively (Table 3). The three constructed genetic linkage maps had in total 2,107 markers representing 1,314 mapped genes (Table 3). Genes were evenly distributed between linkage groups (χ2 test p > 0.05 for the three linkage maps). Map coverage was 65–100% and average distance between two adjacent markers was smaller than 2 cM (Table 3). The vast majority of markers with distorted segregations were discarded because of insufficient linkage information to be mapped (Table 3). Out of the six distorted markers, five mapped in the first 10 and 20 cM of LG 5 in Oria6 and consensus maps, respectively (see Additional files 4 and 5).

Table 3 Mapping features of the two parental linkage maps (Gal1056 and Oria6) and consensus map for the cross (GxO)

Through comparisons between both parental maps, as well as with previously developed maps for P. pinaster[20, 81, 82] based on 654 common markers, 12 groups could be identified for the three maps, which is in agreement with the haploid number of chromosomes for the species. Common markers among the different genetic maps compared mapped always in the same homologous LG excepting three markers (see Additional file 6): contigs FN696780 and AL749831 that mapped in LG 9 and LG 4 in Chancerel et al. [20] and in LG 7 and LG 9 in this study, respectively (see Additional files 4 and 5); contig CT577280 that mapped in LG 7 and LG 4 in the two different maps obtained in Chancerel et al. [20] while in Gal1056 and the consensus map it was mapped in LG 4.

For the 82% and 86% of contigs with more than one mapped SNP, they mapped at less than 1 cM in Gal1056 and Oria6 respectively. There was a significant exception for contig BX249015 that had one SNP mapped in LG 5 (BX249015-204) in Gal1056, Oria6 and the consensus map and the other SNP mapped in LG 8 (BX249015-289) in Oria6 and the consensus map (see Additional files 4 and 5), whereas this contig was mapped in LG 5 in Chancerel et al. [20].

The consensus linkage map is available at Dendrome (http://dendrome.ucdavis.edu/cmap/).

QTL detection

Of the 55 detected QTL (Table 4, Figures 1 and 2), 28 were highly significant QTL, whereas the remaining 27 could be considered as suggestive or putative QTL. QTL were detected for all traits but the higher number of QTL were detected for Fv’Fm’ and ΦPSII (Table 4, Figures 1 and 2). The total phenotypic variance explained for a single QTL ranged from 4.6% (WUEi) to 20.9% (Fv’Fm’). The higher percentage of total phenotypic variance explained by all the QTL detected for a trait in a time-point of measurement was 44% (Fv’Fm’).

Table 4 Identified QTL in Gal1056, Oria6 and GxO maps
Figure 1
figure 1

Graphical representation of the QTL identified on the parental (Gal1056 and Oria6) and consensus (GxO) framework linkage maps (LGs 1 to 6). Only linkage groups for which QTL have been detected are presented. The one and two LOD confidence intervals are indicated by squares and lines, respectively. Colored QTL are the significant QTL (significant at 99% at chromosome level or 95% at genome wide level), with each color representing a different trait, and black QTL are the suggestive QTL (significant at 95% confidence at chromosome level). An = net photosynthetic rate (μmol CO2m-2 s-1); gsw = stomatal conductance to water vapour (molH2Om-2 s-1); WUEi = intrinsic water use efficiency (μmol CO2 molH2O-1); δ13C = isotopic composition of 13C (‰); SLA = specific leaf area (m2Kg-1); Fv’Fm’ = maximum efficiency of PSII under light conditions; ΦPSII = quantum yield. S1, S2 and S3 stand for 1st time-point of measurement (well watered plants), 2nd time-point of measurement (seven days without watering) and 3rd time-point of measurement (14 days without watering) respectively.

Figure 2
figure 2

Graphical representation of the QTL identified on the parental (Gal1056 and Oria6) and consensus (GxO) framework linkage maps (LGs 7 to 12). Only linkage groups for which QTL have been detected are presented. The one and two LOD confidence intervals are indicated by squares and lines, respectively. Colored QTL are the significant QTL (significant at 99% at chromosome level or 95% at genome wide level), with each color representing a different trait, and black QTL are the suggestive QTL (significant at 95% confidence at chromosome level). An = net photosynthetic rate (μmol CO2m-2 s-1); gsw = stomatal conductance to water vapour (molH2Om-2 s-1); WUEi = intrinsic water use efficiency (μmol CO2 molH2O-1); δ13C = isotopic composition of 13C (‰); SLA = specific leaf area (m2Kg-1); Fv’Fm’ = maximum efficiency of PSII under light conditions; ΦPSII = quantum yield. S1, S2 and S3 stand for 1st time-point of measurement (well watered plants), 2nd time-point of measurement (seven days without watering) and 3rd time-point of measurement (14 days without watering) respectively.

Consequently, four QTL hotspots could be identified in LG 5, LG 6, LG 7 and LG 12 (Table 4, Figures 1 and 2) due to the co-localization of QTL for different traits (Figures 1 and 2): SLA co-localized with gsw in LG 5 and LG 7; Fv’Fm’ co-localized with WUEi in LG 5 and LG 12 with δ13C in LG 6, with An in LG 6 and LG 12 and with ΦPSII in LG 7 and LG 12. QTL for SLA, Fv’Fm’ and ΦPSII were detected for the three time-points of measurement. Some of them co-localized in the same region, such as the identified for SLA in LG 5, LG 7 and LG 12 and for Fv’Fm’ in LG 6 and LG 7. Co-localization of QTL for the same traits at different levels of water stress highlights the stability of QTL with the imposition of drought stress. QTL for An could only be detected for the 2nd and 3rd time-point of measurement while QTL for gsw and WUEi were only detected in the 1st and 2nd time-point of measurement (Table 4).

Candidate gene identification

The 74% of the mapped sequences (991 out of 1,348) were annotated. Gene annotations and co-localization with the detected QTL lead to the identification of 58 positional candidate genes that could be involved in the expression of the targeted traits (see Additional file 7). Genes related with oxidative stress, ATPase family proteins or proteins of the light harvesting centers were found in the confidence intervals of QTL for net photosynthesis or chlorophyll fluorescence traits. Genes related with stomatal regulation, ABA signaling pathways or cell wall composition were found in QTL for gsw and WUEi. Genes expressed under drought conditions co-localized with QTL identified in the 2nd or 3rd time-point of measurements but not in the first one, which could be pointing out the induced drought functional response of P. pinaster. Other remarkable co-localizations were found for two QTL for SLA with an enzyme involved in auxin biosynthesis, or between a QTL for δ13C and a member of the aquaporin family (see Additional file 7 for a detailed list of candidate genes). ANOVA test developed for the 73 tested SNPs in candidate genes resulted in 43 significant associations with at least one of the analyzed traits (data not shown). After corrections using the false discovery rate estimated, only three SNPs of the gene MYB1 (m746, m747 and m751) remained significantly associated with Fv’Fm’ measured in the 1st time-point (well watered plants) and SLA measured in the 3rd time-point (14 days without watering; Figure 3). SNP m746 was located in an intron but m747 and m751 were located in exon regions. The base substitution in SNP m747 was a non-synonymous change between a threonine (when a cytosine is present) and an isoleucine (when a thymine is present) while in m751 was a synonymous change. SNPs m747 and m751 explained 14.4% and 12.6% of the phenotypic variance for Fv’Fm’ measured in the 1st time-point and SLA measured in the 3rd time-point respectively. SNP m746 explained 9% of the total phenotypic variance for Fv’Fm’ measured in the 1st time-point.

Figure 3
figure 3

Boxplots for SNPs in candidate gene MYB 1 significantly associated with traits. For each one of the three SNPs (gene name and SNP position between brackets) is shown the p-value of ANOVA and false discovery rate (q-value) estimated for Fv’Fm’ (maximum efficiency of PSII under light conditions) and SLA (specific leaf area). S1, S2 and S3 stand 1st, 2nd and 3rd time-point of measurements respectively.

Discussion

Highly saturated linkage maps

Combining different types of markers three highly dense linkage maps were constructed. They include more than 1,000 genes scattered throughout the genome of P. pinaster and distributed in 12 groups that match the chromosome number of the species. The aforementioned highly saturated maps, with less than 2 cM mean distance between markers, are in the range of recently published linkage maps for other conifer species [20, 22, 24, 27, 51]. Estimated map length was higher in the female than in the male parent. Differences in genome length between parental maps are usually found in conifer species [9093] and it may be a consequence of differences in the recombination rate between parental trees [94, 95].

The accuracy of the SNP genotyping assays previously proved [20, 82] has been confirmed in this study by genes with more than one SNP that mapped in almost all cases within a distance lower than 1 cM. The single exception of contig BX249015 could be attributed to the existence of two paralogous genes for this sequence placed in different LGs. Indeed, high levels of synteny and colinearity were observed between female and male parental maps. The fact that four out of the five distorted markers mapped in the same region suggests that segregation distortion could be due to pre or post-zygotic selection rather than to genotyping errors.

The construction of dense genetic maps for different conifers provides additional tools for studying conifer genomes organization and evolution at a finer scale [27]. In addition, high density linkage maps can be used to position scaffolds along linkage groups contributing to the assembly of a reference genome sequence [24, 96]. P. pinaster genome sequencing is currently in progress, and it should be noted that Oria6, the male progenitor of the mapping family, is the genotype from which the haploid line was selected and its DNA used as template [97].

Additionally, development of dense genetic maps from individuals belonging to two Spanish natural populations (from Northwest coast and Southeast mountains) that show high levels of genetic divergence with the Spanish (from the Castilian Plateau; [81]) and French populations (from Landes and Corsica; [20]), from which segregating progenies have been previously mapped, is important to explore the genetic organization and evolution of the species. Synteny and colinearity were highly conserved when compared with 654 common markers with previous studies [20, 81, 82]. Only three discrepancies were found that supposed just a 0.46% over all the common markers analyzed: Contig CT577280 was mapped in LG 7 and LG 4 in two of three obtained maps for P. pinaster in Chancerel et al. [20] and it was suggested the existence of two paralogous genes for this sequence. In this study the position of CT577280 in LG 4 was confirmed. Contigs FN696780 and AL749831 mapped in LG 9 and LG 4 in Chancerel et al. [20] and in LG 7 and LG 9 in this study respectively, which suggest also the existence of two paralogous genes for these sequences.

The high level of synteny and colinearity observed between the genetic maps developed for individuals that belong to populations with very different genetic backgrounds [98] points out the high reliability in the marker order obtained. Thus, it is possible the development of a composite genetic map for the species by integrating the genetic maps developed by de Miguel et al. [81], Chancerel et al. [20] and those obtained in this work, which is currently in progress. Parental maps are the most accurate regarding both, marker order and marker distances; since they have been constructed through separated information of the meiosis occurred in each progenitor. Accuracy is also related with the presence of genotyping errors, missing values and segregation distortion in the molecular marker data used for the construction of linkage maps [99]. In this study it was achieved by the thorough genotypic data integrity obtained by using highly stringent thresholds to consider SNPs for mapping. Also, the position of the SNPs genotyped in fewer individuals (SNP array D) was validated through the comparison with previously developed P. pinaster maps [20]. In addition, only a few distorted markers have been mapped and almost all located in a narrow region of a single linkage group (Oria6 and GxO LG 5), which points towards a probable biological origin. Even so, QTL analyses have been developed using framework linkage maps to minimize the problems that possible errors in marker order could cause.

QTL detection

One of the main goals of this work was to identify QTL for leaf functional traits related to photosynthesis and water use efficiency in response to drought. QTL analysis in forest tree species is challenging by its long generation times which hinder the development of classical mapping populations like backcross, F2 or recombinant inbred lines. In order to overcome this shortcoming alternative strategies are usually developed for QTL detection in trees, such as the two-way pseudo-testcross [83] used in this study. In this work, two parents from contrasting populations in their drought response were selected to maximize the variability of the F1 obtained progeny, at molecular and functional levels. Although some recent QTL studies in trees worked with larger progenies [51, 100, 101], the 162 obtained siblings in this study are in the range or higher than other QTL analysis in trees [31, 46, 69, 70, 102]. On the other side, gas exchange parameters are extremely sensible to variations in the environmental conditions. To cope with the problem of environmental noise in phenotypic evaluation, different strategies have been used for QTL analysis in the literature. For example the implementation of statistical and physiological models to adjust phenotypic values for microclimatic differences [41, 50] or the development of inbred line populations for annual species [103, 104]. In this study, four vegetative copies of each genotype were established in a completely randomized block design in a greenhouse. All these efforts allowed the identification of significant and suggestive QTL for important traits related to water use efficiency and carbon uptake in P. pinaster.

For all the analyzed traits several QTL were identified with moderate effects rather than a single or few QTL with large effect, as expected for complex functional traits. The only exception was δ13C for which only one QTL could be identified. For gas exchange parameters, the percentage of observed phenotypic variance explained when taking into account all significant QTL detected in a single trait in each time-point of measurements was in average 20%. In addition, some suggestive QTL for water use efficiency estimated from different approaches have been reported, although their existence should be further tested using larger population sizes or analyzing their stability in different genetic backgrounds [105]. For chlorophyll fluorescence parameters, a large number of QTL were identified. The significant QTL detected in each time-point of measurement for Fv’Fm’ explained together up to 44% of the observed phenotypic variance. The results achieved point towards a tight genetic control of photochemical traits, as previously described in crop species like maize [103, 106], barley [104] or soybean [40].

For δ13C only one QTL was detected in LG 6, in agreement with Brendel et al. [44] that found a QTL in the same region of LG 6 for P. pinaster. No co-localization of QTL for WUEi and δ13C was found in spite of the significant phenotypic and genetic correlation between both traits. In this study δ13C was measured only in the 1st time-point of measurement, thus its value probably reflected the water use efficiency in well watered conditions. As the genotypes analyzed have showed high genetic variability in water use efficiency in response to drought [76], it could be expected to observe higher variation in δ13C in needles developed under water limiting conditions maintained in a long-lasting water stress period [107109]. Higher variability on δ13C would enhance the detection of QTL for this trait and it might be possible to find other QTL as reported by Brendel et al. [44], who found four significant and four suggestive QTL. Differences in the number of detected QTL for δ13C between both studies could also be explained because Brendel et al. [44] measured δ13C in tree growth rings from 15 years-old trees while in this study δ13C was measured in needles of one year-old seedlings. Nevertheless, the co-localization of this QTL for δ13C between two genetically unrelated progenies from wide geographic origins (Landes x Landes versus Galicia x Oria) and growing under different environmental conditions supports its stability.

Interaction between QTL and environmental conditions was tested performing the QTL analyses using three different time-points of measurements corresponding to different water irrigation regimes. In general, most of the detected QTL were environment-specific, suggesting that genes are differentially activated during maritime pine drought response [110]. Nevertheless, several QTL for Fv’Fm’, ΦPSII and SLA were less sensitive to environmental conditions and maintained the same location with drought imposition, confirming the stability of these QTL across different levels of water stress endured by plants.

This way, four clusters of QTL were identified in LG 5, LG 6, LG 7 and LG 12. Clustering of QTL could be related with the pleiotropic effect of one or a few genes affecting different traits rather to the existence of rich gene regions, as genes were homogeneously distributed between LGs. Chancerel et al. [20] detected higher number of genes in LG 6 and LG 12 than in the other linkage groups, however the maps developed in this study could not confirm these results.

QTL for photosynthesis measured through gas exchange and chlorophyll fluorescence parameters co-localized in LG 6 and LG 12, accordingly with the high broad-sense genetic correlation found between both traits. However, additional and no co-localizing QTL were identified for these traits in other LGs, suggesting that CO2 fixation and electron transport were not entirely coupled, in agreement with Gu et al. [41]. Uncoupling of these two processes may be due to drought effects on stomatal conductance, biochemical alterations of carbon fixation enzymes, or photoinhibition affecting electron transport rate [7, 8].

Under drought stress An and gsw showed a lower level of phenotypic correlation while the correlation coefficient between An and Fv’Fm’ or ΦPSII increased with water stress, which suggests that under stomatal closure the differences that can be observed between genotypes in carbon fixation could be due to differences in electron transport through PSII rather than to differences in stomatal conductance, as previously observed in other species [106].

SLA showed a significant phenotypic and genetic correlation with WUEi and Fv’Fm’. The identification of relationships between two traits using phenotypic correlations may not distinguish whether the traits could be causally related or simply varying in association. However, the coincidence of QTL for two traits is strong evidence that they could be functionally related [36]. QTL co-localization of SLA with WUEi and Fv’Fm’ was found in LG 5 and LG 7 pointing towards a strong inter-relationship between SLA, WUEi and Fv’Fm’. The aforementioned co-localization could indicate that plants with lower SLA are more efficient in water use but had a lower efficiency of electron transport through photosystem II that could be explained because of the higher importance of gsw over An in determining WUEi in this species [76, 111113].

Most of the detected QTL were found only in one of the two progenitors. The parental trees were selected from two distant populations, showing high level of genetic differentiation, and with a different degree of drought tolerance. Oria6 came from the southeast of Spain governed by a Mediterranean climate with long, hard and frequent summer dry periods, while Gal1056 came from the northwest of Spain where Atlantic climate is present. Consequently, a higher degree of drought adaptation is expected in Oria6 than in Gal1056. Controlled crosses performed with so different parental trees in their response to drought are very useful to compare QTL identified in individuals with different genetic backgrounds.

Candidate genes within QTL

The identification of the gene or genes underlying a trait has been described as one of the greatest challenge for geneticists during this century [114]. The development of high density linkage maps using gene-based markers selected, in some cases, for their known implication in drought response allowed the identification of potential candidate genes for the quantitative multigenic traits analyzed in this study. Due to the lack of sequence annotation, a considerable number of mapped sequences showing high homology with cDNA sequences from other conifers could not be functionally inspected. Thus, some QTL with large effect had no obvious candidate genes but hold great promise to identify unknown genes underlying the corresponding processes in the future. For other QTL, positional candidate genes with known function in other species that were selected according to their functional similarity with genes involved in processes related with the studied trait were identified. A MIXTA-LIKE TRANSCRIPTION FACTOR (MYB) and a HISTONE CHAPERONE were found at 25 and 12 cM from the LOD peak of one of the four most clearly detected QTL, Fv’Fm’S2LG7i. MYB transcription factors are a wide group related with multiple physiological processes such as photosynthesis signaling [115]. The HISTONE CHAPERONE acts as a heat protection protein [116]. The increase of leaf temperature could be an important consequence under drought stress conditions due to reduced transpiration caused by stomatal closure. In this sense, the gene encoding the MYB transcription factor and the HISTONE CHAPERONE also co-localized with Fv’Fm’S3LG7i, both QTL measured under water stress. Another gene of the MYB family encoding the MYB 1 transcription factor, co-localized with several QTL for Fv’ Fm’ and ΦPSII measured in well watered conditions (Fv’Fm’S1LG7i, Fv’Fm’S1LG7_2f, ΦPSII S1LG7_2f) and SLA measured under water stress (SLAS3LG7i). In this sense, three SNPs positioned in MYB 1 gene resulted in a significant association with Fv’Fm’ measured in the 1st time-point of measurements (well watered plants) and SLA measured in the 3rd time-point of measurements (14 days without watering). Lepoittevin et al. [117] found that the gene MYB 1 showed complete linkage disequilibrium in P. pinaster over a distance of 1,304 bp. Together with their intron/exon location and base substitution types, this points towards association of SNPs m751 and m746 with target traits could be the consequence of genetic linkage with m747, that had higher chance to influence Fv’Fm’ and SLA. The expression of MYB 1 regulates genes of the phenylalanine pathway in white spruce [118] and maritime pine [119]. Increase of isoprenoid related compounds has been described to be related with photoprotection mechanisms triggering under abiotic stresses [120]. In this respect, some of the SNPs observed for MYB 1 at present could be related with enhancing maintenance of photochemistry function as higher Fv’Fm’ during drought. These associations should be further validated analyzing, i.e. specific nucleotide variants in a panel of unrelated genotypes [121].

Several genes related with oxidative stress co-localized with QTL for photosynthesis under water stress conditions inferred both by gas exchange or chlorophyll fluorescence. For example, 5 -ADENYLSULFATE REDUCTASE-LIKE 4-LIKE that was implicated in the cell redox homeostasis [122], co-localized with QTL ΦPSIIS1LG8_2m; PROLYL 4-HYDROXYLASE ALPHA SUBUNIT-LIKE PROTEIN that has oxidoreductase activity [123], co-localized with QTL AnS2LG9m; or CINNAMOYL- REDUCTASE 1-LIKE and PEROXIREDOXIN- CHLOROPLASTIC-LIKE that are enzymes from the flavonoid and phenylpropanoid biosynthesis pathways, respectively [124], were on the confidence interval of QTL Fv’Fm’S2LG1_2m. Overall, gene annotation seems to point out to an important role of maintenance photochemical integrity machinery in the drought response of P. pinaster.

Several genes that have been described to be related with regulation of stomatal aperture were found in the range of QTL for gsw and WUEi. For example, MALATE DEHYDROGENASE catalyzes the reaction which converts malate to oxalacetate and a reduction in malate before stomatal closure was observed [125, 126]. Also, PHOSPHOLIPASE C 3-LIKE is required for the control of stomatal aperture by ABA [127, 128]. Genes encoding these enzymes co-localized with QTL gswS2LG5f, WUEiS2LG5f, WUEiS2LG5i. AQUAPORIN NIP1-2-LIKE co-localized with δ13CS1LG6i, which was found interesting because of the importance of aquaporins in determining the leaf water status [129] and the proved stability of this QTL.

Conclusions

The in-depth analysis of genetic control of the CO2 fixation process in response to drought was possible after measuring different functional parameters using complementary techniques, such as gas exchange and chlorophyll fluorescence, that measure final carbon capacity uptake. The use of maritime pine replicated genotypes and a suitable experimental design have made possible to identify genetic control for functional and morphological leaf traits, measured under three water irrigation regimes as they are highly dependent on environmental conditions. Several genomic regions implicated in the genetic control of drought resistance traits have been identified. The identification of potential candidate genes leads this project a step beyond the simple detection of QTL. Nonetheless, further association studies with proposed candidate genes are needed in order to validate detected SNP marker-trait associations.

References

  1. Chaves M, Davies B: Drought effects and water use efficiency: improving crop production in dry environments. Funct Plant Biol. 2010, 37: 3-5.

    Google Scholar 

  2. Pinheiro C, Chaves MM: Photosynthesis and drought: can we make metabolic connections from available data?. J Exp Bot. 2011, 62: 869-882.

    CAS  PubMed  Google Scholar 

  3. Jones H: Partitioning stomatal and non stomatal limitations to photosynthesis. Plant Cell Environ. 1985, 8: 95-104.

    Google Scholar 

  4. Foyer CH, Neukermans J, Queval G, Noctor G, Harbinson J: Photosynthetic control of electron transport and the regulation of gene expression. J Exp Bot. 2012, 63: 1637-1661.

    CAS  PubMed  Google Scholar 

  5. Cano FJ, Sánchez-Gómez D, Rodríguez-Calcerrada J, Warren CR, Gil L, Aranda I: Effects of drought on mesophyll conductance and photosynthetic limitations at different tree canopy layers. Plant Cell Environ. 2013, 36: 1961-1980. doi:10.1111/pce.12103

    PubMed  Google Scholar 

  6. Flexas J, Medrano H: Drought inhibition of photosynthesis in C3 plants: stomatal and non-stomatal limitations revisited. Ann Botany. 2002, 89: 183-189.

    CAS  Google Scholar 

  7. Flexas J, Barbour MM, Brendel O, Cabrera HM, Carriquí M, Diaz-Espejo A, Douthe C, Dreyer E, Jp JP, Gago J: Mesophyll diffusion conductance to CO2: An unappreciated central player in photosynthesis. Plant Sci. 2012, 193–194: 70-84.

    PubMed  Google Scholar 

  8. González-Martínez SC, Krutovsky KV, Neale DB: Forest-tree population genomics and adaptive evolution. New Phytol. 2006, 170: 227-238.

    PubMed  Google Scholar 

  9. Neale DB, Kremer A: Forest tree genomics: growing resources and applications. Nat Rev Genet. 2011, 12: 111-122.

    CAS  PubMed  Google Scholar 

  10. Neale D, Savolainen O: Association genetics of complex traits in conifers. Trends Plant Sci. 2004, 9: 325-330.

    CAS  PubMed  Google Scholar 

  11. Thavamanikumar S, Southerton SG, Bossinger G, Thumma BR: Dissection of complex traits in forest trees: opportunities for marker-assisted selection. TGG. 2013, 9: 627-639.

    Google Scholar 

  12. Mochida K, Shinozaki K: Genomics and bioinformatics resources for crop improvement. Plant Cell Physiol. 2010, 51: 497-523.

    CAS  PubMed Central  PubMed  Google Scholar 

  13. Deschamps S, Campbell MA: Utilization of next-generation sequencing platforms in plant genomics and genetic variant discovery. Mol Breed. 2010, 25: 553-570.

    CAS  Google Scholar 

  14. Rigault P, Boyle B, Lepage P, Cooke JEK, Bousquet J, MacKay JJ: A white spruce gene catalog for conifer genome analyses. Plant Physiol. 2011, 157: 14-28.

    CAS  PubMed Central  PubMed  Google Scholar 

  15. Zhou Y, Gao F, Liu R, Feng J, Li H: De novo sequencing and analysis of root transcriptome using 454 pyrosequencing to discover putative genes associated with drought tolerance in Ammopiptanthus mongolicus. BMC Genomics. 2012, 13: 266-

    CAS  PubMed Central  PubMed  Google Scholar 

  16. Qiu Q, Ma T, Hu Q, Liu B, Wu Y, Zhou H, Wang Q, Wang J, Liu J: Genome-scale transcriptome analysis of the desert poplar. Populus euphratica. Tree Physiol. 2011, 31: 452-461.

    PubMed  Google Scholar 

  17. Li X, Wu HX, Dillon SK, Southerton SG: Generation and analysis of expressed sequence tags from six developing xylem libraries in Pinus radiata D Don. BMC Genomics. 2009, 10: 41-

    CAS  PubMed Central  PubMed  Google Scholar 

  18. Perdiguero P, Collada C, Barbero MC, Casado GG, Cervera MT, Soto A: Identification of water stress genes in Pinus pinaster Ait. by controlled progressive stress and suppression-subtractive hybridization. Plant Physiol Bioch. 2012, 50: 44-53.

    CAS  Google Scholar 

  19. Monclus R, Leplé JC, Bastien C, Bert PF, Villar M, Marron N, Brignolas F, Jorge V, Orléans F, Champenoux F: Integrating genome annotation and QTL position to identify candidate genes for productivity, architecture and water-use efficiency in Populus spp. BMC Plant Biol. 2012, 12: 173-

    CAS  PubMed Central  PubMed  Google Scholar 

  20. Chancerel E, Lamy J-B, Lesur I, Noirot C, Klopp C, Ehrenmann F, Boury C, Le Provost G, Label P, Lalanne C: High-density linkage mapping in a pine tree reveals a genomic region associated with inbreeding depression and provides clues to the extent and distribution of meiotic recombination. BMC Biol. 2013, 11: 50-

    CAS  PubMed Central  PubMed  Google Scholar 

  21. Echt C, Saha S, Krutovsky K, Wimalanathan K, Erpelding J, Liang C: An annotated genetic map of loblolly pine based on microsatellite and cDNA markers. BMC Genet. 2011, 12: 17-

    CAS  PubMed Central  PubMed  Google Scholar 

  22. Eckert AJ, Pande B, Ersoz ES, Wright MH, Rashbrook VK, Nicolet CM, Neale DB: High-throughput genotyping and mapping of single nucleotide polymorphisms in loblolly pine (Pinus taeda L.). TGG. 2009, 5: 225-234. 61

    Google Scholar 

  23. Eckert AJ, van Heerwaarden J, Wegrzyn JL, Nelson CD, Ross-Ibarra J, Gonzalez-Martinez SC, Neale D: Patterns of population structure and environmental associations to aridity across the range of loblolly pine (Pinus taeda L., Pinaceae). Genetics. 2010, 185: 969-982.

    CAS  PubMed Central  PubMed  Google Scholar 

  24. Martínez-García PJ, Stevens KA, Wegrzyn JL, Liechty J, Crepeau M, Langley CH, Neale DB: Combination of multipoint maximum likelihood (MML) and regression mapping algorithms to construct a high-density genetic linkage map for loblolly pine (Pinus taeda L.). TGG. 2013, 9: 1529-1535.

    Google Scholar 

  25. Moriguchi Y, Ujino-Ihara T, Futamura N, Saito M, Ueno S, Matsumoto A, Tani N, Taira H, Shinohara K, Tsumura Y: The construction of a high-density linkage map for identifying SNP markers that are tighly linked to a nuclear-recessive major gene for male sterility in Cryptomeria japonica D.Don. BMC Genomics. 2012, 19: 95-

    Google Scholar 

  26. Pavy N, Pelgas B, Beauseigle S, Blais S, Gagnon F, Gosselin I, Lamothe M, Isabel N, Bousquet J: Enhancing genetic mapping of complex genomes through the design of highly-multiplexed SNP arrays: application to the large and unsequenced genomes of white spruce and black spruce. BMC Genomics. 2008, 9: 21-

    PubMed Central  PubMed  Google Scholar 

  27. Pavy N, Pelgas B, Laroche J, Rigault P, Isabel N, Bousquet J: A spruce gene map infers ancient plant genome reshuffling and subsequent slow evolution in the gymnosperm lineage leading to extant conifers. BMC Biol. 2012, 10: 84-

    CAS  PubMed Central  PubMed  Google Scholar 

  28. Sehgal D, Rajaram V, Armstead IP, Vadez V, Yadav YP, Hash CT, Yadav RS: Integration of gene-based markers in a pearl millet genetic map for identification of candidate genes underlying drought tolerance quantitative trait loci. BMC Plant Biol. 2012, 12: 9-

    CAS  PubMed Central  PubMed  Google Scholar 

  29. Prunier J, Pelgas B, Gagnon F, Desponts M, Isabel N, Beaulieu J, Bousquet J: The genomic architecture and association genetics of adaptive characters using a candidate SNP approach in boreal black spruce. BMC Genomics. 2013, 14: 368-

    CAS  PubMed Central  PubMed  Google Scholar 

  30. Pflieger S, Lefebvre V, Causse M: The candidate gene approach in plant genetics: a review. Mol Breed. 2001, 7: 275-291.

    CAS  Google Scholar 

  31. Casasoli M, Pot D, Plomion C, Monteverdi MC, Barreneche T, Lauteri M, Villani F: Identification of QTLs affecting adaptive traits in Castanea sativa Mill Plant. Cell Environ. 2004, 27: 1088-1101.

    CAS  Google Scholar 

  32. Hurme P, Sillanpaa MJ, Arjas E, Repo T, Savolainen O: Genetic basis of climatic adaptation in scots pine by bayesian quantitative trait locus analysis. Genetics. 2000, 156: 1309-1322.

    CAS  PubMed Central  PubMed  Google Scholar 

  33. Lebreton C, Lazić-Jančić V, Steed A, Pekić S, Quarrie S: Identification of QTL for drought responses in maize and their use in testing causal relationships between traits. J Exp Bot. 1995, 46: 853-865.

    CAS  Google Scholar 

  34. Whan A, Robinson N, Lakshmanan P, Schmidt S, Aitken K: A quantitative genetics approach to nitrogen use efficiency in sugarcane. Funct Plant Biol. 2010, 37: 448-454.

    Google Scholar 

  35. Price A, Courtois B: Mapping QTLs associated with drought resistance in rice: Progress, problems and prospects. Plant Growth Regul. 1999, 29: 123-133.

    CAS  Google Scholar 

  36. Thumma BR, Naidu BP, Chandra A, Cameron DF, Bahnisch LM, Liu C: Identification of causal relationships among traits related to drought resistance in Stylosanthes scabra using QTL analysis. J Exp Bot. 2001, 52: 203-

    CAS  PubMed  Google Scholar 

  37. Juenger TE, Mckay JK, Hausmann N, Keurentjes JJB, Sen S, Stowe KA, Dawson TE, Simms EL, Richards JH: Identification and characterization of QTL underlying whole plant physiology in Arabidopsis thaliana: δ13C, stomatal conductance and transpiration efficiency. Plant Cell Environ. 2005, 28: 697-708.

    CAS  Google Scholar 

  38. Specht J, Chase K, Macrander M, Graef G, Chung J, Markwell J, Germann M, Orf J, Lark K: Soybean response to water: A QTL analysis of drought tolerance. Crop Sci. 2001, 41: 493-509.

    CAS  Google Scholar 

  39. Agbicodo E, Fatokun C, Muranaka S, Visser R: Breeding drought tolerant cowpea: constraints, accomplishments, and future prospects. Euphytica. 2009, 167: 353-370.

    Google Scholar 

  40. Yin Z, Meng F, Song H, He X, Xu X, Yu D: Mapping quantitative trait loci associated with chlorophyll a fluorescence parameters in soybean (Glycine max (L.) Merr.). Planta. 2010, 231: 875-885.

    CAS  PubMed  Google Scholar 

  41. Gu J, Yin X, Struik PC, Stomph TJ, Wang H: Using chromosome introgression lines to map quantitative trait loci for photosynthesis parameters in rice (Oryza sativa L.) leaves under drought and well-watered field conditions. J Exp Bot. 2012, 63: 455-469.

    CAS  PubMed Central  PubMed  Google Scholar 

  42. Fischer R, Edmeades GO: Breeding and cereal yield progress. Crop Sci. 2010, 50: 85-98.

    Google Scholar 

  43. Brendel O, Le Thiece D, Scotti-Saintagne C, Bodénès C, Kremer A, Guehl J-M: Quantitative trait loci controlling water use efficiency and related traits in Quercus robur L. TGG. 2008, 4: 263-278.

    Google Scholar 

  44. Brendel O, Pot D, Plomion C, Rozenberg P, Guehl JM: Genetic parameters and QTL analysis of δ13C and ring width in maritime pine. Plant Cell Environ. 2002, 25: 945-953.

    CAS  Google Scholar 

  45. Muchero W, Sewell MM, Ranjan P, Gunter LE, Tschaplinski TJ, Yin T, Tuskan GA: Genome anchored QTLs for biomass productivity in hybrid Populus grown under contrasting environments. PLoS One. 2013, 8: e54468-

    CAS  PubMed Central  PubMed  Google Scholar 

  46. Freeman JS, Potts BM, Downes GM, Pilbeam D, Thavamanikumar S, Vaillancourt R: Stability of quantitative trait loci for growth and wood properties across multiple pedigrees and environments in Eucalyptus globulus. New Phytol. 2013, 198: 1121-1134.

    CAS  PubMed  Google Scholar 

  47. Paterson AH, Lander ES, Hewitt JD, Peterson S, Lincoln S, Tanksley S: Resolution of quantitative traits into Mendelian factors by using a complete linkage map of restriction fragment lenght polymorphisms. Nature. 1988, 335: 721-726.

    CAS  PubMed  Google Scholar 

  48. Kruglyak L, Lander ES: High-resolution genetic mapping of complex traits. Am J Hum Genet. 1995, 56: 1212-

    CAS  PubMed Central  PubMed  Google Scholar 

  49. Paran I, Zamir D: Quantitative traits in plants: beyond the QTL. Trends Genet. 2003, 19: 303-306.

    CAS  PubMed  Google Scholar 

  50. Yin X, Kropff MJ, Stam P: The role of ecophysiological models in QTL analysis: the example of specific leaf area in barley. Heredity. 1999, 82: 415-421.

    PubMed  Google Scholar 

  51. Pelgas B, Bousquet J, Meirmans P, Ritland K, Isabel N: QTL mapping in white spruce: gene maps and genomic regions underlying adaptive traits across pedigrees, years and environments. BMC Genomics. 2011, 12: 145-

    PubMed Central  PubMed  Google Scholar 

  52. Granda E, Camarero JJ, Gimeno TE, Martínez-Fernández J, Valladares F: Intensity and timing of warming and drought differentially affect growth patterns of co-occurring Mediterranean tree species. Eur J Forest Res. 2013, 132: 469-480.

    Google Scholar 

  53. Körner C, Sarris D, Christodoulakis D: Long-term increase in climatic dryness in the East-Mediterranean as evidenced for the island of Samos. Regional Environ Change. 2005, 5: 27-36.

    Google Scholar 

  54. Sarris D, Christodoulakis D, KÖRNER C: Recent decline in precipitation and tree growth in the eastern Mediterranean. Glob Chang Biol. 2007, 13: 1187-1200.

    Google Scholar 

  55. Christensen JH, Hewitson B, Busuioc A, Chen A, Gao X, Held R, Jones R, Kolli RK, Kwon W, Laprise R: Regional climate projections. Contribution of Working group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Climate Change, 2007: The Physical Science Basis. Edited by: Solomon S, Qin D, Manning M, Marquis M, Averyt K, Tignor MMB, Miller HLR, Chen Z. 2007, Cambridge, UK and New York, USA: Cambridge University Press, 847-940.

    Google Scholar 

  56. Blanco E, Casado MA, Costa M, Escribano R, García-Anton M, Génova M, Gómez-Manzaneque A, Gómez-Manzaneque F, Moreno JC, Morla C, Regato P, Sainz H: Los bosques ibéricos. Una interpretación geobotánica. 2005, Barcelona: Planeta, 4

    Google Scholar 

  57. Tadesse W, Nanos N, Aunon F, Arrabal C, Garcia C, Gil L, Alia R, Pardos J: Genetic improvement of resin yield from maritime pine in Spain. Forest Chem Rev. 2001, 111: 11-

    Google Scholar 

  58. INF3: Inventario Nacional Forestal. Edited by http://www.magrama.gob.es/es/biodiversidad/temas/inventarios-nacionales/inventario-forestal-nacional/; 2007

  59. Grivet D, Sebastiani F, Alía R, Bataillon T, Torre S, Zabal-Aguirre M, Vendramin GG, González-Martínez SC: Molecular footprints of local adaptation in two Mediterranean conifers. Mol Biol Evol. 2011, 28: 101-

    CAS  PubMed  Google Scholar 

  60. Eveno E, Collada C, Guevara MA, Léger V, Soto A, Díaz L, Léger P, González-Martínez SC, Cervera MT, Plomion C, Garnier-Géré P: Contrasting patterns of selection at Pinus pinaster Ait. drought stress candidate genes as revealed by genetic differentiation analyses. Mol Biol Evol. 2008, 25: 417-437.

    CAS  PubMed  Google Scholar 

  61. Martínez-Vilalata J, Piñol J: Drought-induced mortality and hydraulic architecture in pine populations of the NE Iberian Peninsula. For Ecol Manag. 2002, 161: 247-256.

    Google Scholar 

  62. Sabaté S, Gracia CA, Sánchez A: Likely effects of climate change on growth of Quercus ilex, Pinus halepensis, Pinus pinaster, Pinus sylvestris and Fagus sylvatica forests in the Mediterranean region. For Ecol Manag. 2002, 162: 23-37.

    Google Scholar 

  63. González-Martínez SC, Huber D, Ersoz E, Davis J, Neale D: Association genetics in Pinus taeda L. II Carbon isotope discrimination. Heredity. 2008, 101: 19-26.

    PubMed  Google Scholar 

  64. González-Martínez SC, Ersoz E, Brown GR, Wheeler NC, Neale DB: DNA sequence variation and selection of tag single-nucleotide polymorphisms at candidate genes for drought-stress response in Pinus taeda L. Genetics. 2006, 172: 1915-

    PubMed Central  PubMed  Google Scholar 

  65. Street NR, Skogstrom O, Sjodin A, Tucker J, Rodriguez-Acosta M, Nilsson P, Jansson S, Taylor G: The genetics and genomics of the drought response in Populus. Plant J. 2006, 48: 321-341.

    CAS  PubMed  Google Scholar 

  66. Tschaplinski TJ, Tuskan GA, Sewell MM, Gebre GM, Todd DE, Pendley CD: Phenotypic variation and quantitative trait locus identification for osmotic potential in an interspecific hybrid inbred F2 poplar pedigree grown in contrasting environments. Tree Physiol. 2006, 26: 595-604.

    PubMed  Google Scholar 

  67. Parelle J, Zapater M, Scotti-Saintagne C, Kremer A, Jolivet Y, Dreyer E, Brendel O: Quantitative trait loci of tolerance to waterlogging in a European oak (Quercus robur L.): physiological relevance and temporal effect patterns. Plant Cell Environ. 2007, 30: 422-434.

    PubMed  Google Scholar 

  68. Gerber S, Lascoux M, Kremer A: Relation between protein markers and quantitative traits in maritime pine (Pinus pinaster Ait.). Silvae Genet. 1997, 46: 286-291.

    Google Scholar 

  69. Markussen T, Fladung M, Achere V, Favre JM, Faivre-Rampant P, Aragones A, DA Silva Pérez D, Havengt L, Ritter E: Identification of QTLs controlling growth, chemical and physical wood property traits in Pinus pinaster, Ait. Silvae Genet. 2003, 52: 8-15.

    Google Scholar 

  70. Pot D, Rodrigues J, Rozenberg P, Chantre G, Tibbits J, Cahalan C, Pichavant F, Plomion C: QTLs and candidate genes for wood properties in maritime pine (Pinus pinaster Ait.). TGG. 2006, 2: 10-24.

    Google Scholar 

  71. Plomion C, Durel CE, O'Malley DM: Genetic dissection of height in maritime pine seedlings raised under accelerated growth conditions. Theor Appl Genet. 1996, 93: 849-858.

    CAS  PubMed  Google Scholar 

  72. Chagné D, Brown GR, Lalanne C, Madur D, Pot D, Neale D, Plomion C: Comparative genome and QTL mapping between maritime and loblolly pines. Mol Breed. 2003, 12: 185-195.

    Google Scholar 

  73. Lepoittevin C, Harvengt L, Plomion C, Garnier-Géré P: Association mapping for growth, straightness and wood chemistry traits in the Pinus pinaster Aquitaine breeding population. TGG. 2012, 8: 113-126.

    Google Scholar 

  74. Plomion C, Yani A, Marpeau A: Genetic determinism of delta 3-carene in maritime pine using RAPD markers. Genome. 1996, 39: 1123-1127.

    CAS  PubMed  Google Scholar 

  75. Budde KB, Heuertz M, Hernández-Serrano A, Pausas JG, Vendramin GG, Verdú M, González-Martínez SC: In situ genetic association for serotiny, a fire-related trait, in Mediterranean maritime pine (Pinus pinaster). New Phytol. 2014, 201: 230-241. doi:10.1111/nph.12483

    CAS  PubMed  Google Scholar 

  76. de Miguel M, Sánchez-Gómez D, Cervera MT, Aranda I: 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. Tree Physiol. 2012, 32: 94-103.

    PubMed  Google Scholar 

  77. Sharp Z: Principles of stable isotope geochemistry. 2007, Upper Saddle River: NJ Pearson Education Inc.

    Google Scholar 

  78. Dellaporta SL, Wood J, Hicks JB: A plant DNA minipreparation: version II. Plant Mol Biol Report. 1983, 1: 19-21.

    CAS  Google Scholar 

  79. Guevara MA, Chagné D, Almeida H, Byrnes M, Collada C, Favre JM, Harvengt L, Jeandroz S, Orazio C, Plomion C, Ramboer A, Rocheta M, Sebastiani F, Soto A, Vendramin GG, Cervera MT: Isolation and characterization of nuclear microsatellite loci in Pinus pinaster Ait. Mol Ecol Notes. 2005, 5: 57-59.

    CAS  Google Scholar 

  80. Chagné D: Développement de marqueurs moléculaires chez le pin maritime (Pinus pinaster Ait.) et cartographie génétique comparée des conifères. PhD. 2004, Nancy I: Nancy, France

    Google Scholar 

  81. de Miguel M, de María N, Guevara MA, Díaz L, Sáez-Laguna E, Sánchez-Gómez D, Chancerel E, Aranda I, Collada C, Plomion C, Cabezas JA, Cervera MT: Annotated genetic linkage maps of Pinus pinaster Ait. from a Central Spain population using microsatellite and gene based markers. BMC Genomics. 2012, 13: 527-

    CAS  PubMed Central  PubMed  Google Scholar 

  82. Chancerel E, Lepoittevin C, Le Provost G, Lin Y-C, Jaramillo-Correa J, Eckert A, Wegrzyn J, Zelenika D, Boland A, Frigerio J-M, Chaumeil P, Garnier-Gere P, Boury C, Grivet D, Gonzalez-Martinez S, Rouze P, Van de Peer Y, Neale D, Cervera M, Kremer A, Plomion C: Development and implementation of a highly-multiplexed SNP array for genetic mapping in maritime pine and comparative mapping with loblolly pine. BMC Genomics. 2011, 12: 368-

    CAS  PubMed Central  PubMed  Google Scholar 

  83. Grattapaglia D, Sederoff R: Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-testcross mapping strategy and RAPD markers. Genetics. 1994, 137: 1121-1137.

    CAS  PubMed Central  PubMed  Google Scholar 

  84. Van Ooijen JW: Joinmap 4, software for the calculation of genetic maps in experimental populations. Edited by: Kiazma BV. 2006, Wageningen

    Google Scholar 

  85. Kosambi D: The estimation of map distances from recombination values. Annals Eugen. 1944, 12: 172-175.

    Google Scholar 

  86. Hulbert S, Ilott T, Legg E, Lincoln S, Lander E, Michelmore R: Genetic analysis of the fungus, Bremia lactucae, using restriction fragment length polymorphisms. Genetics. 1988, 120: 947-

    CAS  PubMed Central  PubMed  Google Scholar 

  87. Chakravarti A, Lasher LK, Reefer JE: A maximum likelihood method for estimating genome length using genetic linkage data. Genetics. 1991, 128: 175-182.

    CAS  PubMed Central  PubMed  Google Scholar 

  88. Van Ooijen JW: MapQTL 6.0, software for the mapping of quantitative trait loci in experimental populations of diploid species. Edited by: Kiazma BV. 2009, Wageningen

    Google Scholar 

  89. Conesa A, Götz S, García-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676.

    CAS  PubMed  Google Scholar 

  90. Kuramoto N, Kondo T, Fujisawa Y, Nakata R, Hayashi E, Goto Y: Detection of quantitative trait loci for wood strength in Cryptomeria japonica. Can J Forest Res-Revue Canadienne De Recherche Forestiere. 2000, 30: 1525-1533.

    CAS  Google Scholar 

  91. Nikaido A, Ujino T, Iwata H, Yoshimura K, Yoshimura H, Suyama Y, Murai M, Nagasaka K, Tsumura Y: AFLP and CAPS linkage maps of Cryptomeria japonica. Theor Appl Genet. 2000, 100: 825-831.

    CAS  Google Scholar 

  92. Scotti I, Burelli A, Cattonaro F, Chagné D, Fuller J, Hedley PE, Jansson G, Lalanne C, Madur D, Neale D, Plomion C, Powell W, Troggio M, Morgante M: Analysis of the distribution of marker classes in a genetic linkage map: a case study in Norway spruce (Picea abies Karst). TGG. 2005, 1: 93-102.

    Google Scholar 

  93. Komulainen P, Brown GR, Mikkonen M, Karhu A, Garcia-Gil MR, O'Malley D, Lee B, Neale DB, Savolainen O: Comparing EST-based genetic maps between Pinus sylvestris and Pinus taeda. Theor Appl Genet. 2003, 107: 667-678.

    CAS  PubMed  Google Scholar 

  94. Sewell MM, Sherman BK, Neale DB: A consensus map for loblolly pine (Pinus taeda L.). I. Construction and integration of individual linkage maps from two outbred three-generation pedigrees. Genetics. 1999, 151: 321-330.

    CAS  PubMed Central  PubMed  Google Scholar 

  95. Plomion C, O'Malley D: Recombination rate differences for pollen parents and seed parents in pine. Heredity. 1996, 77: 341-350.

    CAS  Google Scholar 

  96. Mackay J, Dean JFD, Plomion C, Peterson DG, Cánovas FM, Pavy N, Ingvarsson PK, Savolainen O, Guevara MÁ, Fluch S, Vinceti B, Abarca D, Díaz-Sala C, Cervera M-T: Towards decoding the conifer giga-genome. Plant Mol Biol. 2012, 80: 555-569.

    CAS  PubMed  Google Scholar 

  97. Arrillaga I, Guevara M, Muñoz-Bertomeu J, Lázaro-Gimeno D, Sáez-Laguna E, Díaz L, Torralba L, Mendoza-Poudereux I, Segura J, Cervera M: Selection of haploid cell lines from megagametophyte cultures of maritime pine as a DNA source for massive sequencing of the species. Plant Cell Tiss Org Cult. 2014, doi:10.1007/s11240-014-0470-z

    Google Scholar 

  98. Bucci G, González-Martínez SC, Le Provost G, Plomion C, Ribeiro MM, Sebastiani F, Alía R, Vendramin GG: Range-wide phylogeography and gene zones in Pinus pinaster Ait. revealed by chloroplast microsatellite markers. Molec Ecol. 2007, 16: 2137-2153.

    CAS  Google Scholar 

  99. Hackett CA, Broadfoot LB: Effects of genotyping errors, missing values and segregation distortion in molecular marker data on the construction of linkage maps. Heredity. 2003, 90: 33-38.

    CAS  PubMed  Google Scholar 

  100. Ukrainetz N, Ritland K, Mansfield S: Identification of quantitative trait loci for wood quality and growth across eight full-sib coastal Douglas-fir families. TGG. 2008, 4: 159-170.

    Google Scholar 

  101. Novaes E, Osorio L, Drost DR, Miles BL, Boaventura-Novaes CRD, Benedict C, Dervinis C, Yu Q, Sykes R, Davis M, Martin TA, Peter GF, Kirst M: Quantitative genetic analysis of biomass and wood chemistry of Populus under different nitrogen levels. New Phytol. 2009, 182: 878-890.

    CAS  PubMed  Google Scholar 

  102. Bartholomé J, Salmon F, Vigneron P, Bouvet J-M, Plomion C, Gion J-M: Plasticity of primary and secondary growth dynamics in Eucalyptus hybrids: a quantitative genetics and QTL mapping perspective. BMC Plant Biol. 2013, 13: 120-

    PubMed Central  PubMed  Google Scholar 

  103. Fracheboud Y, Jompuk C, Ribaut J, Stamp P, Leipner J: Genetic analysis of cold-tolerance of photosynthesis in maize. Plant Mol Biol. 2004, 56: 241-253.

    CAS  PubMed  Google Scholar 

  104. Guo P, Baum M, Varshney RK, Graner A, Grando S, Ceccarelli S: QTLs for chlorophyll and chlorophyll fluorescence parameters in barley under post-flowering drought. Euphytica. 2008, 163: 203-214.

    CAS  Google Scholar 

  105. Kearsey M, Farquhar A: QTL analysis in plants; where are we now?. Heredity. 1998, 80: 137-142.

    PubMed  Google Scholar 

  106. Fracheboud Y, Ribaut JM, Vargas M, Messmer R, Stamp P: Identification of quantitative trait loci for cold-tolerance of photosynthesis in maize (Zea mays L.). J Exp Bot. 2002, 53: 1967-1977.

    CAS  PubMed  Google Scholar 

  107. Aranda I, Alía R, Ortega U, Dantas AK, Majada J: Intra-specific variability in biomass partitioning and carbon isotopic discrimination under moderate drought stress in seedlings from four Pinus pinaster populations. TGG. 2010, 6: 169-178.

    Google Scholar 

  108. Correia I, Almeida MH, Aguiar A, Alía R, David TS, Pereira JS: Variations in growth, survival and carbon isotope composition (δ13C) among Pinus pinaster populations of different geographic origins. Tree Physiol. 2008, 28: 1545-1552.

    PubMed  Google Scholar 

  109. Corcuera L, Gil-Pelegrin E, Notivol E: Phenotypic plasticity in Pinus pinaster delta C13: environment modulates genetic variation. Ann For Sci. 2010, 67: 812-

    Google Scholar 

  110. Watkinson JI, Sioson AA, Vasquez-Robinet C, Shukla M, Kumar D, Ellis M, Heath LS, Ramakrishnan N, Chevone B, Watson LT: Photosynthetic acclimation is reflected in specific patterns of gene expression in drought-stressed loblolly pine. Plant Physiol. 2003, 133: 1702-1716.

    CAS  PubMed Central  PubMed  Google Scholar 

  111. Fernández M, Gil L, Pardos JA: Effects of water supply on gas exchange in Pinus pinaster Ait. provenances during their first growing season. Ann For Sci. 2000, 57: 9-16.

    Google Scholar 

  112. Fernandez M, Novillo C, Pardos JA: Effects of water and nutrient availability in Pinus pinaster Ait. Open pollinated families at an early age: Growth, gas exchange and water relations. New For. 2006, 31: 321-342.

    Google Scholar 

  113. Guehl J, Fort C, Ferhi A: Differential response of leaf conductance, carbon isotope discrimination and water-use efficiency to nitrogen deficiency in maritime pine and pedunculate oak plants. New Phytol. 1995, 131: 149-157.

    Google Scholar 

  114. Price AH: Believe it or not, QTLs are accurate!. Trends Plant Sci. 2006, 11: 213-216.

    CAS  PubMed  Google Scholar 

  115. Saibo NJM, Lourenço T, Oliveira MM: Transcription factors and regulation of photosynthetic and related metabolism under environmental stresses. Ann Botany. 2009, 103: 609-623.

    CAS  Google Scholar 

  116. Kotak S, Larkindale J, Lee U, von Koskull-Döring P, Vierling E, Scharf K-D: Complexity of the heat stress response in plants. Curr Opin Plant Biol. 2007, 10: 310-316.

    CAS  PubMed  Google Scholar 

  117. Lepoittevin C, Garnier-Gere P, Hubert F, Plomion C: Strong linkage disequilibrium and balanced selection in Pinus pinaster transcription factors putatively involved in wood formation. Oral presentation to IUFRO-CTIA Joint Conference “Adaptation, Breeding and Conservation in the Era of Forest Tree Genomics and Environmental Change”. 2008, IUFRO-CTIA: Québec City (Canada), 26-28.

    Google Scholar 

  118. Bomal C, Duval I, Giguère I, Fortin E, Caron S, Stewart D, Boyle B, Séguin A, MacKay JJ: Opposite action of R2R3-MYBs from different subgroups on key genes of the shikimate and monolignol pathways in spruce. J Exp Bot. 2014, 65: 495-508.

    CAS  PubMed Central  PubMed  Google Scholar 

  119. Craven-Bartle B, Pascual MB, Cánovas FM, Ávila C: A Myb transcription factor regulates genes of the phenylalanine pathway in maritime pine. Plant J. 2013, 74: 755-766.

    CAS  PubMed  Google Scholar 

  120. Peñuelas J, Munné-Bosc S: Isoprenoids: an evolutionary pool for photoprotection. Trends Plant Sci. 2005, 10: 166-169.

    PubMed  Google Scholar 

  121. Harjes CE, Rocheford TR, Bai L, Brutnell TP, Kandianis CB, Sowinski SG, Stapleton AE, Vallabhaneni R, Williams M, Wurtzel ET: Natural genetic variation in lycopene epsilon cyclase tapped for maize biofortification. Science. 2008, 319: 330-

    CAS  PubMed Central  PubMed  Google Scholar 

  122. Martin MN, Tarczynski MC, Shen B, Leustek T: The role of 5'-adenylylsulfate reductase in controlling sulfate reduction in plants. Photosynth Res. 2005, 86: 309-323.

    CAS  PubMed  Google Scholar 

  123. Vlad F, Spano T, Vlad D, Bou Daher F, Ouelhadj A, Kalaitzis P: Arabidopsis prolyl-hydroxylases are differentially expressed in response to hypoxia, anoxia and mechanical wounding. Physiol Plant. 2007, 130: 471-483.

    CAS  Google Scholar 

  124. Baier M, Noctor G, Foyer CH, Dietz K-J: Antisense suppression of 2-cysteine peroxiredoxin in Arabidopsis specifically enhances the activities and expression of enzymes associated with ascorbate metabolism but not glutathione metabolism. Plant Physiol. 2000, 124: 823-832.

    CAS  PubMed Central  PubMed  Google Scholar 

  125. Riccardi F, Gazeau P, de Vienne D, Zivy M: Protein changes in response to progressive water deficit in maize quantitative variation and polypeptide identification. Plant Physiol. 1998, 117: 1253-1263.

    CAS  PubMed Central  PubMed  Google Scholar 

  126. Cramer GR, Van Sluyter SC, Hopper DW, Pascovici D, Keighley T, Haynes PA: Proteomic analysis indicates massive changes in metabolism prior to the inhibition of growth and photosynthesis of grapevine (Vitis vinifera L.) in response to water deficit. BMC Plant Biol. 2013, 13: 49-

    CAS  PubMed Central  PubMed  Google Scholar 

  127. Chaves MM, Pereira JS, Maroco JP: Understanding plant responses to drought—from genes to the whole plant. Funct Plant Biol. 2003, 30: 239-264.

    CAS  Google Scholar 

  128. Hunt L, Mills LN, Pical C, Leckie CP, Aitken FL, Kopka J, Mueller-Roeber B, McAinsh MR, Hetherington AM, Gray JE: Phospholipase C is required for the control of stomatal aperture by ABA. Plant J. 2003, 34: 47-55.

    CAS  PubMed  Google Scholar 

  129. Prado K, Boursiac Y, Tournaire-Roux C, Monneuse J-M, Postaire O, Da Ines O, Schäffner AR, Hem S, Santoni V, Maurel C: Regulation of Arabidopsis Leaf Hydraulics Involves Light-Dependent Phosphorylation of Aquaporins in Veins. The Plant Cell Online. 2013, 25: 1029-1039.

    CAS  Google Scholar 

  130. Ajjawi I, Coku A, Froehlich JE, Yang Y, Osteryoung KW, Benning C, Last RL: A J-like protein influences fatty acid composition of chloroplast lipids in Arabidopsis. PLoS One. 2011, 6: e25368-

    CAS  PubMed Central  PubMed  Google Scholar 

  131. Andrews M, Huizinga DH, Crowell DN: The CaaX specificities of Arabidopsis protein prenyltransferases explain era1 and ggb phenotypes. BMC Plant Biol. 2010, 10: 118-

    PubMed Central  PubMed  Google Scholar 

  132. Banci L, Bertini I, Luchinat C, Turano P: Electron Transfer, Respiration, and Photosynthesis. Biological Inorganic Chemistry: Structure and Reactivity. Edited by: Bertini I, Gray H, Stiefel EI, Valentine JS. 2007, USA: University Science Books, 229-261.

    Google Scholar 

  133. Barkla BJ, Vera-Estrella R, Maldonado-Gama M, Pantoja O: Abscisic acid induction of vacuolar H + -ATPase activity in Mesembryanthemum crystallinum is developmentally regulated. Plant Physiol. 1999, 120: 811-820.

    CAS  PubMed Central  PubMed  Google Scholar 

  134. Bazakos C, Manioudaki ME, Therios I, Voyiatzis D, Kafetzopoulos D, Awada T, Kalaitzis P: Comparative Transcriptome Analysis of Two Olive Cultivars in Response to NaCl-Stress. PLoS One. 2012, 7: e42931-

    CAS  PubMed Central  PubMed  Google Scholar 

  135. Beligni MV, Mayfield SP: Arabidopsis thaliana mutants reveal a role for CSP41a and CSP41b, two ribosome-associated endonucleases, in chloroplast ribosomal RNA metabolism. Plant Mol Biol. 2008, 67: 389-401.

    CAS  PubMed  Google Scholar 

  136. Brandt U: Proton-translocation by membrane-bound NADH: ubiquinone-oxidoreductase (complex I) through redox-gated ligand conduction. Biochim Biophys Acta. 1997, 1318: 79-91.

    CAS  PubMed  Google Scholar 

  137. Bray EA: Genes commonly regulated by water-deficit stress in Arabidopsis thaliana. J Exp Bot. 2004, 55: 2331-2341.

    CAS  PubMed  Google Scholar 

  138. Çakır B, Olcay AC: Molecular cloning, phylogenetic analysis, and expression profiling of a grape CMP-sialic acid transporter-like gene induced by phytohormone and abiotic stress. Gen Genom. 2013, 35: 225-238.

    Google Scholar 

  139. Dietz KJ, Tavakoli N, Kluge C, Mimura T, Sharma SS, Harris GC, Chardonnens AN, Golldack D: Significance of the V-type ATPase for the adaptation to stressful growth conditions and its regulation on the molecular and biochemical level. J Exp Bot. 2001, 52: 1969-1980.

    CAS  PubMed  Google Scholar 

  140. Fujimoto SY, Ohta M, Usui A, Shinshi H, Ohme-Takagi M: Arabidopsis ethylene-responsive element binding factors act as transcriptional activators or repressors of GCC box-mediated gene expression. Plant Cell Online. 2000, 12: 393-404.

    CAS  Google Scholar 

  141. Golldack D, Li C, Mohan H, Probst N: Gibberellins and abscisic acid signal crosstalk: living and developing under unfavorable conditions. Plant Cell Rep. 2013, 32: 1007-1016.

    CAS  PubMed  Google Scholar 

  142. Gómez J, Sanchez-Martínez D, Stiefel V, Rigau J, Puigdomenech P, Pages M: A gene induced by the plant hormone abscisic acid in response to water stress encodes a glycine-rich protein. Nature. 1988, 334: 262-264.

    PubMed  Google Scholar 

  143. Goyer A, Collakova E, de la Garza RD, Quinlivan EP, Williamson J, Gregory JF, Shachar-Hill Y, Hanson AD: 5-Formyltetrahydrofolate is an inhibitory but well tolerated metabolite in Arabidopsis leaves. J Biol Chem. 2005, 280: 26137-26142.

    CAS  PubMed  Google Scholar 

  144. Gross EL: Plastocyanin: Structure, Location, Diffusion and Electron Transfer Mechanisms. Oxygenic Photosynthesis: The Light Reactions. Edited by: Ort DR, Yocum CF, Heichel IF. 2004, Netherlands: Springer, 413-429.

    Google Scholar 

  145. Guo L, Devaiah SP, Narasimhan R, Pan X, Zhang Y, Zhang W, Wang X: Cytosolic glyceraldehyde-3-phosphate dehydrogenases interact with phospholipase Dδ to transduce hydrogen peroxide signals in the Arabidopsis response to stress. Plant Cell Online. 2012, 24: 2200-2212.

    CAS  Google Scholar 

  146. Ham BK, Park JM, Lee SB, Kim MJ, Lee IJ, Kim KJ, Kwon CS, Paek KH: Tobacco Tsip1, a DnaJ-Type Zn Finger Protein, Is Recruited to and Potentiates Tsi1-Mediated Transcriptional Activation. Plant Cell. 2006, 18: 2005-2020.

    CAS  PubMed Central  PubMed  Google Scholar 

  147. Hardie DG: Plant protein serine/threonine kinases: Classification and Functions. Annu Rev Plant Physiol Plant Mol Biol. 1999, 50: 97-131.

    CAS  PubMed  Google Scholar 

  148. Hernández I, Van Breusegem F: Opinion on the possible role of flavonoids as energy escape valves: Novel tools for nature's Swiss army knife?. Plant Sci. 2010, 179: 297-301.

    Google Scholar 

  149. Ioannidis NE, Cruz JA, Kotzabasis K, Kramer DM: Evidence That Putrescine Modulates the Higher Plant Photosynthetic Proton Circuit. PLoS One. 2012, 7: e29864-

    CAS  PubMed Central  PubMed  Google Scholar 

  150. Janicka-Russak M: Plant plasma membrane H + -ATPase in adaptation of plants to abiotic stresses. Abiotic stress response in plants - physiological, biochemical and genetic perspectives. Edited by: Shanker A, Venkateswarlu B. 2011, Rijeka, Croatia: Intech, 197-218.

    Google Scholar 

  151. Jansson S, Gustafsson P: Evolutionary conservation of the chlorophyll a/b binding proteins cDNAs encoding Type I, II and III LHC I polypeptides from the gymnosperm Scots pine. Mol Genet Genom. 1991, 229: 67-76.

    CAS  Google Scholar 

  152. Kawasaki T, Koita H, Nakatsubo T, Hasegawa K, Wakabayashi K, Takahashi H, Umemura K, Umezawa T, Shimamoto K: Cinnamoyl-CoA reductase, a key enzyme in lignin biosynthesis, is an effector of small GTPase Rac in defense signaling in rice. Proc Nat Acad Sci U S A. 2006, 103: 230-235.

    CAS  Google Scholar 

  153. Kim J, Malladi A, van Iersel MW: Physiological and molecular responses to drought in Petunia: the importance of stress severity. J Exp Bot. 2012, 63: 6335-6345.

    CAS  PubMed Central  PubMed  Google Scholar 

  154. Kim JS, Jung HJ, Lee HJ, Kim KA, Goh C-H, Woo Y, Oh SH, Han YS, Kang H: Glycine-rich RNA-binding protein7 affects abiotic stress responses by regulating stomata opening and closing in Arabidopsis thaliana. Plant J. 2008, 55: 455-466.

    CAS  PubMed  Google Scholar 

  155. Kim S, Choi HI, Ryu HJ, Park JH, Kim MD, Kim SY: ARIA, an Arabidopsis arm repeat protein interacting with a transcriptional regulator of abscisic acid-responsive gene expression, is a novel abscisic acid signaling component. Plant Physiol. 2004, 136: 3639-3648.

    CAS  PubMed Central  PubMed  Google Scholar 

  156. Kiyosue T, Yamaguchi-Shinozaki K, Shinozaki K: Cloning of cDNAs for genes that are early-responsive to dehydration stress (ERDs) in Arabidopsis thaliana L.: identification of three ERDs as HSP cognate genes. Plant Mol Biol. 1994, 25: 791-798.

    CAS  PubMed  Google Scholar 

  157. Kushwaha HR, Singh AK, Sopory SK, Singla-Pareek SL, Pareek A: Genome wide expression analysis of CBS domain containing proteins in Arabidopsis thaliana (L.) Heynh and Oryza sativa L. reveals their developmental and stress regulation. BMC Genomics. 2009, 10: 200-

    PubMed Central  PubMed  Google Scholar 

  158. Li AD, Anderson LE: Expression and Characterization of Pea Chloroplastic Glyceraldehyde-3-Phosphate Dehydrogenase Composed of Only the B-Subunit. Plant Physiol. 1997, 115: 1201-1209.

    CAS  PubMed Central  PubMed  Google Scholar 

  159. Liu F, Guo J, Bai P, Duan Y, Wang X, Chen Y, Feng H, Huang L, Kang Z: Wheat TaRab7 GTPase is part of the signaling pathway in responses to stripe rust and abiotic stimuli. PLoS One. 2012, 7: e37146-

    CAS  PubMed Central  PubMed  Google Scholar 

  160. Miao Y, Lv D, Wang P, Wang XC, Chen J, Miao C, Songa CP: An Arabidopsis glutathione peroxidase functions as both a redox transducer and a scavenger in abscisic acid and drought stress responses. Plant Cell. 2006, 18: 2749-2766.

    CAS  PubMed Central  PubMed  Google Scholar 

  161. Mishra M, Das R, Pandey GK: Role of ethylene responsive factors (ERFs) in abiotic stress mediated signaling in plants. e-J BiolSci. 2009, 1: 133-146.

    Google Scholar 

  162. Munnik T, Irvine R, Musgrave A: Phospholipid signalling in plants. Biochim Biophys Acta-Lipids and Lipid Metabolism. 1998, 1389: 222-272.

    CAS  Google Scholar 

  163. Noir S, Bömer M, Takahashi N, Ishida T, Tsui TL, Balbi V, Shanahan H, Sugimoto K, Devoto A: Jasmonate controls leaf growth by repressing cell proliferation and the onset of endoreduplication while maintaining a potential stand-by mode. Plant Physiol. 2013, 161: 1930-1951.

    CAS  PubMed Central  PubMed  Google Scholar 

  164. País SM, Tellez-Iñón M, Capiati DA: Serine/threonine protein phosphatases type 2A and their roles in stress signaling. Plant Signal Behav. 2009, 4: 1013-1015.

    PubMed Central  PubMed  Google Scholar 

  165. Pedone KH, Der CJ: Small GTPase. Encyclopedia of Biological Chemistry. Edited by: Lennarz WJ, Lane MD. 2013, USA: Academic Press, 242-248.

    Google Scholar 

  166. Pichersky E, Jansson S: The Light-Harvesting Chlorophyll a/b-Binding Polypeptides and Their Genes in Angiosperm and Gymnosperm Species. Oxygenic Photosynthesis: The Light Reactions. Edited by: Ort D, Yocum C, Heichel I. 2004, Netherlands: Springer, 507-521.

    Google Scholar 

  167. Quilliam LA: Ras Family. Encyclopedia of Biological Chemistry. Edited by: Lennarz WJ, Lane MD. 2013, London, UK: Academic Press, 12-16.

    Google Scholar 

  168. Ribas-Carbo M, Taylor NL, Giles L, Busquets S, Finnegan PM, Day DA, Lambers H, Medrano H, Berry JA, Flexas J: Effects of water stress on respiration in soybean leaves. Plant Physiol. 2005, 139: 466-473.

    CAS  PubMed Central  PubMed  Google Scholar 

  169. Sakamoto H, Maruyama K, Sakuma Y, Meshi T, Iwabuchi M, Shinozaki K, Yamaguchi-Shinozaki K: Arabidopsis Cys2/His2-type zinc finger proteins function as transcription repressors under drought, cold, and high-salinity stress conditions. Plant Physiol. 2004, 136: 2734-2746.

    CAS  PubMed Central  PubMed  Google Scholar 

  170. Saxena M, Bisht R, Roy SD, Sopory SK, Bhalla-Sarin N: Cloning and characterization of a mitochondrial glyoxalase II from Brassica juncea that is upregulated by NaCl, Zn, and ABA. Biochem Biophys Res Commun. 2005, 336: 813-819.

    CAS  PubMed  Google Scholar 

  171. Seidel T: Structure and Regulation of Plant Vacuolar H + -ATPase. Prog Bot. 2009, 70: 93-126.

    CAS  Google Scholar 

  172. Shimazaki K, Doi M, Assmann SM, Kinoshita T: Light regulation of stomatal movement. Annu Rev Plant Biol. 2007, 58: 219-247.

    CAS  PubMed  Google Scholar 

  173. Sugano S, Kaminaka H, Rybka Z, Catala R, Salinas J, Matsui K, Ohme-Takagi M, Takatsuji H: Stress-responsive zinc finger gene ZPT2-3 plays a role in drought tolerance in petunia. Plant J. 2003, 36: 830-841.

    CAS  PubMed  Google Scholar 

  174. Thidholm E, Lindstrom V, Tissier C, Robinson C, Schroder W, Funk C: Novel approach reveals localization and assembly pathway of the PsbS and PsbW proteins into the photosystem II dimer. FEBS Lett. 2002, 513: 217-222.

    CAS  PubMed  Google Scholar 

  175. Torres GA, Gimenes MA, de Rosa VE, Quecini V: Identifying water stress-response mechanisms in citrus by in silico transcriptome analysis. Genet Mol Biol. 2007, 30: 888-905.

    CAS  Google Scholar 

  176. Tovar-Méndez A, Miernyk JA, Randall DD: Regulation of pyruvate dehydrogenase complex activity in plant cells. Eur J Biochem. 2003, 270: 1043-1049.

    PubMed  Google Scholar 

  177. Van Nocker S, Ludwig P: The WD-repeat protein superfamily in Arabidopsis: conservation and divergence in structure and function. BMC Genomics. 2003, 4: 50-

    PubMed Central  PubMed  Google Scholar 

  178. Wang B, Luttge U, Ratajczak R: Effects of salt treatment and osmotic stress on V-ATPase and V-PPase in leaves of the halophyte Suaeda salsa. J Exp Bot. 2001, 52: 2355-2365.

    CAS  PubMed  Google Scholar 

  179. Xu C, Jing R, Mao X, Jia X, Chang X: A wheat (Triticum aestivum) protein phosphatase 2A catalytic subunit gene provides enhanced drought tolerance in tobacco. Ann Bot. 2007, 99: 439-450.

    CAS  PubMed Central  PubMed  Google Scholar 

  180. Yao X, Xiong W, Ye T, Wu Y: Overexpression of the aspartic protease ASPG1 gene confers drought avoidance in Arabidopsis. J Exp Bot. 2012, 63: 2579-2593.

    CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgments

This work was supported by the Spanish projects MAPINSEQ (AGL2009-10496; Spanish Ministry of Science and Innovation), PinCoxSeq (AGL2012-35175; Ministry of Economy and Competitiveness) and the Plant-KBBE project SUSTAINPINE (PLE2009-0016). The research leading to these results has also received funding from the European Union’s Seventh Framework Programme (FP7/2007-2013) under grant agreement n° 289841 (ProCoGen). L Alté and S Ferrándiz are gratefully acknowledged for their assistance and Rose Daniels for the revision of the quality of written English.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to María-Teresa Cervera.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MdM: SSR genotyping, phenotypic evaluation, genetic maps, QTL analysis, candidate genes search and wrote the first draft of the manuscript. JAC: genetic maps, QTL analysis and candidate genes search. NdM: SAMPL genotyping, design of SNP array C and candidate genes search. DS-G: phenotypic evaluation. MAG: SSR genotyping, SAMPL genotyping, design of SNP array C and candidate genes search. MDV: SAMPL and SNP genotyping and candidate genes search. ES-L: DNA extraction and functional annotations. LD: SSR genotyping, SAMPL genotyping. JAM: phenotypic evaluation. MCB: DNA extraction, SAMPL genotyping. CC: design of SNP array C and candidate genes search. CD-S: design of SNP array C and candidate genes search. M-TC and IA: conceived and designed the experiments and collected funding. All authors have read and approved the final version of the manuscript.

Electronic supplementary material

12864_2013_6163_MOESM1_ESM.xlsx

Additional file 1: Identical markers based on recombination rate. Not positioned markers correspond to unlinked markers or markers which position could not be reliably estimated. (XLSX 31 KB)

Additional file 2: Broad sense genetic correlations (±standard error) between the analyzed traits.(DOCX 13 KB)

Additional file 3: Broad sense heritability (estimate ± standard error).(DOCX 13 KB)

Additional file 4: Parental linkage maps for Gal1056, Oria6 and consensus map for both progenitors (GxO).(PDF 281 KB)

12864_2013_6163_MOESM5_ESM.xlsx

Additional file 5: Mapped markers in parental linkage maps for Gal1056, Oria6 and consensus map for both progenitors (GxO).(XLSX 213 KB)

12864_2013_6163_MOESM6_ESM.pdf

Additional file 6: Marker order comparison with maps obtained by Chancerel et al. [20].(PDF 63 KB)

12864_2013_6163_MOESM7_ESM.xlsx

Additional file 7: Candidate genes within QTL [130180].(XLSX 30 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

de Miguel, M., Cabezas, JA., de María, N. et al. 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. BMC Genomics 15, 464 (2014). https://doi.org/10.1186/1471-2164-15-464

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-464

Keywords