High-density SNP assay development for genetic analysis in maritime pine (Pinus pinaster)

Maritime pine provides essential ecosystem services in the south-western Mediterranean basin, where it covers around 4 million ha. Its scattered distribution over a range of environmental conditions makes it an ideal forest tree species for studies of local adaptation and evolutionary responses to climatic change. Highly multiplexed single nucleotide polymorphism (SNP) genotyping arrays are increasingly used to study genetic variation in living organisms and for practical applications in plant and animal breeding and genetic resource conservation. We developed a 9k Illumina Infinium SNP array and genotyped maritime pine trees from (i) a three-generation inbred (F2) pedigree, (ii) the French breeding population and (iii) natural populations from Portugal and the French Atlantic coast. A large proportion of the exploitable SNPs (2052/8410, i.e. 24.4%) segregated in the mapping population and could be mapped, providing the densest ever gene-based linkage map for this species. Based on 5016 SNPs, natural and breeding populations from the French gene pool exhibited similar level of genetic diversity. Population genetics and structure analyses based on 3981 SNP markers common to the Portuguese and French gene pools revealed high levels of differentiation, leading to the identification of a set of highly differentiated SNPs that could be used for seed provenance certification. Finally, we discuss how the validated SNPs could facilitate the identification of ecologically and economically relevant genes in this species, improving our understanding of the demography and selective forces shaping its natural genetic diversity, and providing support for new breeding strategies.

Transcriptional analysis of differentially expressed genes in response to stem inclination in young seedlings of pine

The gravitropic response in trees is a widely studied phenomenon, however understanding of the molecular mechanism involved remains unclear. The purpose of this work was to identify differentially expressed genes in response to inclination using a comparative approach for two conifer species. Young seedlings were subjected to inclination and samples were collected at four different times points. First, suppression subtractive hybridisation (SSH) was used to identify differentially regulated genes in radiata pine (Pinus radiata D. Don). cDNA libraries were constructed from the upper and lower part of inclined stems in a time course experiment, ranging from 2.5 h to 1 month. From a total of 3092 sequences obtained, 2203 elements were assembled, displaying homology to a public database. A total of 942 unigene elements were identified using bioinformatic tools after redundancy analysis. Of these, 614 corresponded to known function genes and 328 to unknown function genes, including hypothetical proteins. Comparative analysis between radiata pine and maritime pine (Pinus pinaster Ait.) was performed to validate the differential expression of relevant candidate genes using qPCR. Selected genes were involved in several functional categories: hormone regulation, phenylpropanoid pathway and signal transduction. This comparative approach for the two conifer species helped determine the molecular gene pattern generated by inclination, providing a set of Pinus gene signatures that may be involved in the gravitropic stress response. These genes may also represent relevant candidate genes involved in the gravitropic response and potentially in wood formation.

Genetic determinism of δ3-carene in maritime pine using RAPD markers

An F2 progeny of maritime pine (Pinus pinaster Aïton) was used to investigate the mode of inheritance of δ3-carene using a quantitative and a qualitative approach. A previously reported genetic map constructed with random amplified polymorphic DNA (RAPD) markers made it possible to locate one major quantitative trait locus (QTL) accounting for most of the phenotypic variation of this trait on linkage group 5. In the qualitative approach, the "C" locus that controls the relative quantity of δ3-carene (C+ for the richness allele and C− for the poorness allele) was found to be strongly associated with RAPD markers in the same genomic region of linkage group 5. The colocation between the QTL and the "C" locus suggests that a major gene or closely linked loci affect the variation in δ3-carene expression. Key words : Pinus pinaster, terpenes, QTL, RAPD, linkage analysis.

Characterisation and natural variation of a dehydrin gene in Quercus petraea (Matt.) Liebl

For the first time in sessile oak [Quercus petraea (Matt.) Liebl.], the isolation and characterisation of a full-length dehydrin gene and its promoter region, as well as its allelic variation in natural populations, is reported. Dehydrins (Dhn) are stress-related genes important for the survival of perennial plants in a seasonal climate. A full-length dehydrin gene (Dhn3) was characterised at the nucleotide level and the protein structure was modelled. Additionally, the allelic variation was analysed in five natural populations of Quercus petraea (Matt.) Liebl. sampled along an altitudinal gradient in the French Pyrenees. The analysed sequences contain typical domains of the K(n) class of dehydrins in the coding region. Also, the 5'untranslated region (promoter) of the gene was amplified, which shows typical motifs essential for drought- and cold-responsive gene expression. Single nucleotide substitutions and indels (insertions/deletions) within the coding region determine large biochemical differences at the protein level. However, only low levels of genetic differentiation between populations from different altitudes were detectable.

Plasticity of maritime pine (Pinus pinaster) wood-forming tissues during a growing season

The seasonal effect is the most significant external source of variation affecting vascular cambial activity and the development of newly divided cells, and hence wood properties. Here, the effect of edapho-climatic conditions on the phenotypic and molecular plasticity of differentiating secondary xylem during a growing season was investigated. Wood-forming tissues of maritime pine (Pinus pinaster) were collected from the beginning to the end of the growing season in 2003. Data from examination of fibre morphology, Fourier-transform infrared spectroscopy (FTIR), analytical pyrolysis, and gas chromatography/mass spectrometry (GC/MS) were combined to characterize the samples. Strong variation was observed in response to changes in edapho-climatic conditions. A genomic approach was used to identify genes differentially expressed during this growing season. Out of 3512 studied genes, 19% showed a significant seasonal effect. These genes were clustered into five distinct groups, the largest two representing genes over-expressed in the early- or late-wood-forming tissues, respectively. The other three clusters were characterized by responses to specific edapho-climatic conditions. This work provides new insights into the plasticity of the molecular machinery involved in wood formation, and reveals candidate genes potentially responsible for the phenotypic differences found between early- and late-wood.

Distribution of genomic regions differentiating oak species assessed by QTL detection

Pedunculate oak and sessile oak are two sympatric interfertile species that exhibit leaf morphological differences. We aimed to detect quantitative trait loci (QTLs) of these traits in order to locate genomic regions involved in species differentiation. A total of 15 leaf morphological traits were assessed in a mixed forest stand composed of Quercus petraea and Q. robur and in a full-sib pedigree of Q. robur. The progeny of the full-sib family were vegetatively propagated in two successive experiments comprising 174 and 216 sibs, and assessments were made on two leaves collected on each of the 1080 and 1530 cuttings corresponding to the two experiments. Traits that exhibited strong species differences in the mixed stand tended also to have higher repeatability values in the mapping population, thus indicating higher genetic control. A genetic map was constructed for QTL detection. Composite interval mapping with the one QTL model was used for QTL detection. From one to three QTLs were detected for 13 traits. In-depth analysis of the QTLs, controlling the five morphological traits that exhibited the highest interspecific differences in the mixed stand, indicated that they were distributed on six linkage groups, with two clusters comprising QTLs of at least two discriminant traits. These results were reinforced when error 1 for QTL detection was set at 5% at the chromosome level, as up to nine clusters could be identified. In conclusion, traits involved in interspecific differentiation of oaks are under polygenic control and widespread in clusters across the genome.

Cross-species transferability and mapping of genomic and cDNA SSRs in pines

Two unigene datasets of Pinus taeda and Pinus pinaster were screened to detect di-, tri- and tetranucleotide repeated motifs using the SSRIT script. A total of 419 simple sequence repeats (SSRs) were identified, from which only 12.8% overlapped between the two sets. The position of the SSRs within their coding sequences were predicted using FrameD. Trinucleotides appeared to be the most abundant repeated motif (63 and 51% in P. taeda and P. pinaster, respectively) and tended to be found within translated regions (76% in both species), whereas dinucleotide repeats were preferentially found within the 5'- and 3'-untranslated regions (75 and 65%, respectively). Fifty-three primer pairs amplifying a single PCR fragment in the source species (mainly P. taeda), were tested for amplification in six other pine species. The amplification rate with other pine species was high and corresponded with the phylogenetic distance between species, varying from 64.6% in P. canariensis to 94.2% in P. radiata. Genomic SSRs were found to be less transferable; 58 of the 107 primer pairs (i.e. 54%) derived from P. radiata amplified a single fragment in P. pinaster. Nine cDNA-SSRs were located to their chromosomes in two P. pinaster linkage maps. The level of polymorphism of these cDNA-SSRs was compared to that of previously and newly developed genomic-SSRs. Overall, genomic SSRs tend to perform better in terms of heterozygosity and number of alleles. This study suggests that useful SSR markers can be developed from pine ESTs.

Comparative mapping between quercus and castanea using simple-sequence repeats (SSRs)

Simple sequence repeat (SSR) markers from Quercus and Castanea were used for comparative mapping between Quercus robur (L.) and Castanea sativa (Mill.). We tested the transferability of SSRs developed in Quercus to Castanea and vice-versa. In total, 47% (25) of the Quercus SSRs and 63% (19) of the Castanea SSRs showed a strong amplification product in the non-source species. From these 44 putative comparative anchor tags, 19 (15 from Quercus and 4 from Castanea) were integrated in two previously established genetic linkage maps for the two genera. SSR loci were sequenced to confirm the orthology of the markers. The combined information from both genetic mapping and sequence analysis were used to determine the homeology between seven linkage groups, aligned on the basis of pairs or triplets of common markers, while two additional groups were matched using a single microsatellite marker. Orthologous loci identified between Q. robur and C. sativa will be useful as anchor loci for comparative mapping studies within the Fagaceae family.

Candidate genes and QTLs for sugar and organic acid content in peach [ Prunus persica (L.) Batsch]

The identification of genes involved in variation of peach fruit quality would assist breeders in creating new cultivars with improved fruit quality. Major genes and quantitative trait loci (QTLs) for physical and chemical components of fruit quality have already been detected, based on the peach [ Prunus persica (L.) Batsch] cv. Ferjalou Jalousia((R)) (low-acid peach) x cv. Fantasia (normally-acid nectarine) F(2) intraspecific cross. Our aim was to associate these QTLs to structural genes using a candidate gene/QTL approach. Eighteen cDNAs encoding key proteins in soluble sugar and organic acid metabolic pathways as well as in cell expansion were isolated from peach fruit. A single-strand conformation polymorphism strategy based on specific cDNA-based primers was used to map the corresponding genes. Since no polymorphism could be detected in the Ferjalou Jalousia((R)) x Fantasia population, gene mapping was performed on the almond [ Prunus amygdalus ( P. dulcis)] cv. Texas x peach cv. Earlygold F(2) interspecific cross from which a saturated map was available. Twelve candidate genes were assigned to four linkage groups of the peach genome. In a second step, the previous QTL detection was enhanced by integrating anchor loci between the Ferjalou Jalousia((R)) x Fantasia and Texas x Earlygold maps and data from a third year of trait assessment on the Ferjalou Jalousia((R)) x Fantasia population. Comparative mapping allowed us to detect a candidate gene/QTL co-location. It involved a cDNA encoding a vacuolar H(+)-pyrophosphatase ( PRUpe;Vp2) that energises solute accumulation, and QTLs for sucrose and soluble solid content. This preliminary result may be the first step in the future development of marker-assisted selection for peach fruit sucrose and soluble solid content.

Comparison of genetic diversity estimates within and among populations of maritime pine using chloroplast simple-sequence repeat and amplified fragment length polymorphism data

We compared the genetic variation of Pinus pinaster populations using amplified fragment length polymorphism (AFLP) and chloroplast simple-sequence repeat (cpSSR) loci. Populations' levels of diversity within groups were found to be similar with AFLPs, but not with cpSSRs. The high interlocus variance associated with the AFLP loci could account for the lack of differences in the former. Although AFLPs revealed much lower genetic diversity than cpSSRs, the levels of among-population differentiation found with the two types of marker were similar, provided that loci showing fewer than four null-homozygotes, in any population, were pruned from the AFLP data. Moreover, the French and Portuguese populations were clearly differentiated from each other, with both markers. The Mantel test showed that the genetic distance matrix calculated using the AFLP data was correlated with the matrix derived from the cpSSRs. Because of the concordance found between markers we conclude that gene flow was indeed the predominant force shaping nuclear and chloroplastic genetic variation of the populations within regions, at the geographical scale studied.

Drought differentially affects expression of a PR-10 protein, in needles of maritime pine (Pinus pinaster Ait.) seedlings

A full length cDNA encoding a PR-10 protein was isolated from maritime pine drought-stressed seedlings. The predicted protein contained 150 amino acids, has a molecular mass of 16.7 kDa and an isoelectric point of 5.28. The transcript level of PR-10 displayed a transient accumulation in needles of drought-stressed plants, and was not detectable in root and stem tissues.

Genetic diversity within and among Pinus pinaster populations: comparison between AFLP and microsatellite markers

Twenty-three populations of Pinus pinaster (13 Aquitaine populations and 10 Corsican populations) were analysed at three microsatellite loci and 122 AFLP loci. The aims of the study were: (i) to compare levels of within-population and among-population diversity assessed with both kinds of markers; (ii) to compare Aquitaine and Corsican provenances of P. pinaster; and (iii) to know if both markers gave the same information for conservation purposes. Classical population genetics statistics were estimated and the ranking of populations obtained using different markers and/or parameters were compared by computing Spearman's rank correlations. Even though microsatellites showed a higher within-population diversity, they showed the same level of differentiation as AFLP markers. Moreover, both markers also showed a higher genetic diversity in the Aquitaine provenance and a higher differentiation among Corsican populations. AFLPs and microsatellites gave different population diversity rankings. Consequently, the results do not support the potential population identification within each provenance for conservation purposes.

Two genes encoding distinct cytosolic glutamine synthetases are closely linked in the pine genome

The major isoenzyme of glutamine synthetase found in leaves of angiosperms is the chloroplastic form. However, pine seedlings contain two cytosolic glutamine synthetases in green cotyledons: GS1a, the predominant isoform, and GS1b, a minor enzyme whose relative amount is increased following phosphinotricin treatment. We have cloned a GS1b cDNA, and comparison with the previously reported GS1a cDNA sequence indicated that they correspond to separate cytosolic GS genes encoding distinct protein products. Phylogenetic analysis showed that the newly reported sequence is closer to cytosolic angiosperm GS than to GS1a, suggesting therefore that GS1a could be a divergent gymnospermous GS1 gene. Gene mapping using a F2 family of maritime pine showed co-localization of both GS genes on group 2 of the genetic linkage map. This result supports the proposed origin of different members of the GS1 family by adjacent gene duplication. The implications for gymnosperm genome organization are discussed.

Compression wood-responsive proteins in developing xylem of maritime pine (Pinus pinaster ait.)

When a conifer shoot is displaced from its vertical position, compression wood (CW) is formed on the under side and can eventually return the shoot to its original position. Changes in cell wall structure and chemistry associated with CW are likely to result from differential gene/protein expression. Two-dimensional polyacrylamide gel electrophoresis of differentiating xylem proteins was combined with the physical characterization of wooden samples to identify and characterize CW-responsive proteins. Differentiating xylem was harvested from a 22-year-old crooked maritime pine (Pinus pinaster Ait.) tree. Protein extracted from different samples were revealed by high-resolution silver stained two-dimensional polyacrylamide gel electrophoresis and analyzed with a computer-assisted system for single spot quantification. Growth strain (GS) measurements allowed xylem samples to be classified quantitatively from normal wood to CW. Regression of lignin and cellulose content on GS showed that an increase in the percentage of lignin and a decrease of the percentage of cellulose corresponded to increasing GS values, i.e. CW. Of the 137 studied spots, 19% were significantly associated with GS effect. Up-regulated proteins included 1-aminocyclopropane-1-carboxylate oxidase (an ethylene forming enzyme), a putative transcription factor, two lignification genes (caffeic O-methyltransferase and caffeoyl CoA-O-methyltransferase), members of the S-adenosyl-L-methionine-synthase gene family, and enzymes involved in nitrogen and carbon assimilation (glutamine synthetase and fructokinase). A clustered correlation analysis was performed to study simultaneously protein expression along a gradient of gravistimulated stressed xylem tissue. Proteins were found to form "expression clusters" that could identify: (a) Gene product under similar control mechanisms, (b) partner proteins, or (c) functional groups corresponding to specialized pathways. The possibility of obtaining regulatory correlations and anticorrelations between proteins provide us with a new category of homology (regulatory homology) in tracing functional relationships.

Proteomics for genetic and physiological studies in plants

Proteomics is becoming a necessity in plant biology, as it is in medicine, zoology and microbiology, for deciphering the function and role of the genes that are or will be sequenced. In this review we focus on the various, mainly genetic, applications of the proteomic tools that have been developed in recent years: characterization of individuals or lines, estimation of genetic variability within and between populations, establishment of genetic distances that can be used in phylogenetic studies, characterization of mutants and localization of the genes encoding the revealed proteins. Improvements in specifically devoted software have permitted precise quantification of the variation in amounts of proteins, leading to the concept of "protein quantity loci" which, combined with the "quantitative trait loci" approach, results in testable hypotheses regarding the role of "candidate proteins" in the metabolism or phenotype under study. This new development is exemplified by the reaction of plants to drought, a trait of major agronomic interest. The accumulation of data regarding genomic and cDNA sequencing will be connected to the protein databases currently developed in plants.

Separation and characterization of needle and xylem maritime pine proteins

Two-dimensional gel electrophoresis (2-DE) and image analysis are currently used for proteome analysis in maritime pine (Pinus pinaster Ait.). This study presents a database of expressed proteins extracted from needles and xylem, two important tissues for growth and wood formation. Electrophoresis was carried out by isoelectric focusing (IEF) in the first dimension and sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) in the second. Silver staining made it possible to detect an average of 900 and 600 spots on 2-DE gels from needles and xylem, respectively. A total of 28 xylem and 35 needle proteins were characterized by internal peptide microsequencing. Out of these 63 proteins, 57 (90%) could be identified based on amino acid similarity with known proteins, of which 24 (42%) have already been described in conifers. Overall comparison of both tissues indicated that 29% and 36% of the spots were specific to xylem and needles, respectively, while the other spots were of identical molecular weight and isoelectric point. The homology of spot location in 2-DE patterns was further validated by sequence analysis of proteins present in both tissues. A proteomic database of maritime pine is accessible on the internet (http://www.pierroton.inra.fr/genetics/2D/).

Water-deficit-responsive proteins in maritime pine

Two-dimensional polyacrylamide gel electrophoresis (2D-PAGE) was used to identify drought-responsive proteins during progressive water deprivation of two-year old maritime pine seedlings. Stress was applied by withholding water during vegetative growth. Needles were sampled before, during and after the stress. Out of about 1000 spots that were quantified by computer analysis, 38 responded during stress. Some proteins were accumulated while others were suppressed. One to three internal microsequences were obtained for 11 proteins, 10 of which were identified on the basis of sequence homologies. These proteins are quite diverse and are involved in photosynthesis, cell elongation, antioxidant metabolism and lignification.

Genetic analysis using trans-dominant linked markers in an F2 family

Trans-dominant linked markers pairs (trans referring to the repulsion linkage phase) provide a model for inferring the F2 progeny genotype based upon both the conditional probabilities of F2 genotypes, given the F2 phenotype, and prior information on marker arrangement. Prior information of marker arrangement can be readily obtained from a linkage analysis performed on marker segregation data in a family resulting by crossing the F1 individual to a "tester" parent or else can be obtained directly from the gametes of the F1, or from recombinant inbred lines. We showed that a trans-dominant linked marker (TDLM) pair can be recoded as a "co-dominant megalocus" when the recombination fraction, r1, for apair of TDLMs is less than 0.05. We obtained a maximum-likelihood estimator (MLE) of the recombination frequency, r2, between a TDLM pair and a co-dominant marker in an F2 family using the EM algorithm. The MLE was biased. Mean bias increased as r1 and r2 increased, and decreased as sample size increased. The information content for r2 was compared to the information content of dominant and co-dominant markers segregating in an F2 family. It was almost identical with two co-dominant markers when r1≤0.01 and r2≥0.05. For larger values of r1, (0.05≤r1≤0.15) a TDLM pair provided 75%-66% of the information content of two co-dominant markers. Although dominant markers can be converted to co-dominant markers by a laborious process of cloning, sequencing, and PCR, TDLM pairs could easily substitute for co-dominant markers in order to detect quantitative trait loci (QTLs) and estimate gene action in an F2 family.

Genetic dissection of height in maritime pine seedlings raised under accelerated growth conditions

Random Amplified Polymorphic DNAs (RAPDs) were used to investigate quantitative trait loci (QTL) for traits related to height growth on 126 F2 seedlings of maritime pine (Pinuspinaster Ait). The haploid megagametophyte was used to determine the maternal genotype of each F2 individual. The seedlings were raised for 2 years in a greenhouse under accelerated growth conditions consisting of intense fertilization combined with continuous light treatments. Total height was measured at different developmental stages, and height growth components were measured after the second growth period. QTLs were identified for each trait. For total height, QTLs of different developmental stages were located on distinct linkage groups. However, rather than a complete temporal change in QTL expression, our results showed that maturation may induce a progressive shift of the genetic control of height growth. This may provide an explanation for a low juvenile-mature phenotypic correlation previously reported for height. Height growth components related to the initiation (controlled by the apical meristem) and elongation of shoot cycles (controlled by the subapical meristem) were mapped to different chromosomes, suggesting that the activity of these meristems is controlled by separate genetic mechanisms.

Genomic analysis in maritime pine (Pinus pinaster). Comparison of two RAPD maps using selfed and open-pollinated seeds of the same individual

Two genomic maps were constructed for one individual tree of maritime pine, Pinus pinaster Ait., using a common set of 263 RAPD markers (random amplified polymorphic DNA). The RAPD markers were chosen from a larger number of polymorphic RAPD fragments on the basis of repeatability and inheritance in a three-generation pedigree. The maps were constructed from two independent mapping samples of 62 megagametophytes (In) from a self cross and from an open-pollinated cross. The markers were grouped (LOD≥4; θ≤0.25) and assigned to 13 major and 5 minor linkage groups. Two framework maps were constructed using the ordering criterion of interval support≥3. Comparison of the two framework maps suggested that the locus order was incorrect for 2% of the framework markers. A bootstrap analysis showed that this error rate was representative for our data set. The results showed that framework maps constructed using RAPD markers were repeatable and that differences in locus order for maps of different genotypes or species could result from chance. The total map distance was 1380 cM, and the map provided coverage of approximately 90% of the genome.