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.

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.

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.