A Strategy for Studying Epigenetic Diversity in Natural Populations: Proof of Concept in Poplar and Oak.

These last 20 years, several techniques have been developed for quantifying DNA methylation, the most studied epigenetic marks in eukaryotes, including the gold standard method, whole-genome bisulphite sequencing (WGBS). WGBS quantifies genome-wide DNA methylation but has several inconveniences rendering it less suitable for population-scale epigenetic studies. The high cost of deep sequencing and the large amounts of data generated prompted us to seek an alternative approach. Restricting studies to parts of the genome would be a satisfactory alternative had there not been a major limitation: the need to select upstream targets corresponding to differentially methylated regions (DMRs) as targets. Given the need to study large numbers of samples, we propose a strategy for investigating DNA methylation variation in natural populations, considering the structural complexity of the genomes with their size and their content in unique as coding regions versus repeated regions as transposable elements. We first identified regions of highly variable DNA methylation in a representative subset of genotypes representative of the biological diversity in the population by WGBS. We then analysed the variations of DNA methylation in these targeted regions at the population level by Sequencing Capture Bisulphite (SeqCapBis). The entire strategy was then validated by applying it to another species. Our strategy was developed as a proof of concept on natural populations of two forest species: Populus nigra and Quercus petraea.

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.

Single-nucleotide polymorphism discovery and validation in high-density SNP array for genetic analysis in European white oaks.

An Illumina Infinium SNP genotyping array was constructed for European white oaks. Six individuals of Quercus petraea and Q. robur were considered for SNP discovery using both previously obtained Sanger sequences across 676 gene regions (1371 in vitro SNPs) and Roche 454 technology sequences from 5112 contigs (6542 putative in silico SNPs). The 7913 SNPs were genotyped across the six parental individuals, full-sib progenies (one within each species and two interspecific crosses between Q. petraea and Q. robur) and three natural populations from south-western France that included two additional interfertile white oak species (Q. pubescens and Q. pyrenaica). The genotyping success rate in mapping populations was 80.4% overall and 72.4% for polymorphic SNPs. In natural populations, these figures were lower (54.8% and 51.9%, respectively). Illumina genotype clusters with compression (shift of clusters on the normalized x-axis) were detected in ~25% of the successfully genotyped SNPs and may be due to the presence of paralogues. Compressed clusters were significantly more frequent for SNPs showing a priori incorrect Illumina genotypes, suggesting that they should be considered with caution or discarded. Altogether, these results show a high experimental error rate for the Infinium array (between 15% and 20% of SNPs potentially unreliable and 10% when excluding all compressed clusters), and recommendations are proposed when applying this type of high-throughput technique. Finally, results on diversity levels and shared polymorphisms across targeted white oaks and more distant species of the Quercus genus are discussed, and perspectives for future comparative studies are proposed.

A linkage disequilibrium perspective on the genetic mosaic of speciation in two hybridizing Mediterranean white oaks.

We analyzed the genetic mosaic of speciation in two hybridizing Mediterranean white oaks from the Iberian Peninsula (Quercus faginea Lamb. and Quercus pyrenaica Willd.). The two species show ecological divergence in flowering phenology, leaf morphology and composition, and in their basic or acidic soil preferences. Ninety expressed sequence tag-simple sequence repeats (EST-SSRs) and eight nuclear SSRs were genotyped in 96 trees from each species. Genotyping was designed in two steps. First, we used 69 markers evenly distributed over the 12 linkage groups (LGs) of the oak linkage map to confirm the species genetic identity of the sampled genotypes, and searched for differentiation outliers. Then, we genotyped 29 additional markers from the chromosome bins containing the outliers and repeated the multilocus scans. We found one or two additional outliers within four saturated bins, thus confirming that outliers are organized into clusters. Linkage disequilibrium (LD) was extensive; even for loosely linked and for independent markers. Consequently, score tests for association between two-marker haplotypes and the 'species trait' showed a broad genomic divergence, although substantial variation across the genome and within LGs was also observed. We discuss the influence of several confounding effects on neutrality tests and review the evolutionary processes leading to extensive LD. Finally, we examine how LD analyses within regions that contain outlier clusters and quantitative trait loci can help to identify regions of divergence and/or genomic hitchhiking in the light of predictions from ecological speciation theory.

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.

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.

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.

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.

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.

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.

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.

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.

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.