Grégoire Le Provost

Range‐wide phylogeography and gene zones in Pinus pinaster Ait. revealed by chloroplast microsatellite markers

Some 1339 trees from 48 Pinus pinaster stands were characterized by five chloroplast microsatellites, detecting a total of 103 distinct haplotypes. Frequencies for the 16 most abundant haplotypes (pk > 0.01) were spatially interpolated over a lattice made by 430 grid points. Fitting of spatially interpolated values on raw haplotype frequencies at the same geographical location was tested by regression analysis. A range‐wide ‘diversity map’ based on interpolated haplotype frequencies allowed the identification of one ‘hotspot’ of diversity in central and southeastern Spain, and two areas of low haplotypic diversity located in the western Iberian peninsula and Morocco. Principal component analysis (PCA) carried out on haplotypes frequency surfaces allowed the construction of a colour‐based ‘synthetic’ map of the first three PC components, enabling the detection of the main range‐scale genetic trends and the identification of three main ‘gene pools’ for the species: (i) a ‘southeastern’ gene pool, including southeastern France, Italy, Corsica, Sardinia, Pantelleria and northern Africa; (ii) an ‘Atlantic’ gene pool, including all the western areas of the Iberian peninsula; and (iii) a ‘central’ gene pool, located in southeastern Spain. Multivariate and amova analyses carried out on interpolated grid point frequency values revealed the existence of eight major clusters (‘gene zones’), whose genetic relationships were related with the history of the species. In addition, demographic models showed more ancient expansions in the eastern and southern ranges of maritime pine probably associated to early postglacial recolonization. The delineation of the gene zones provides a baseline for designing conservation areas in this key Mediterranean pine.

Bioinformatic analysis of ESTs collected by Sanger and pyrosequencing methods for a keystone forest tree species: oak

BackgroundThe Fagaceae family comprises about 1,000 woody species worldwide. About half belong to the Quercus family. These oaks are often a source of raw material for biomass wood and fiber. Pedunculate and sessile oaks, are among the most important deciduous forest tree species in Europe. Despite their ecological and economical importance, very few genomic resources have yet been generated for these species. Here, we describe the development of an EST catalogue that will support ecosystem genomics studies, where geneticists, ecophysiologists, molecular biologists and ecologists join their efforts for understanding, monitoring and predicting functional genetic diversity.ResultsWe generated 145,827 sequence reads from 20 cDNA libraries using the Sanger method. Unexploitable chromatograms and quality checking lead us to eliminate 19,941 sequences. Finally a total of 125,925 ESTs were retained from 111,361 cDNA clones. Pyrosequencing was also conducted for 14 libraries, generating 1,948,579 reads, from which 370,566 sequences (19.0%) were eliminated, resulting in 1,578,192 sequences. Following clustering and assembly using TGICL pipeline, 1,704,117 EST sequences collapsed into 69,154 tentative contigs and 153,517 singletons, providing 222,671 non-redundant sequences (including alternative transcripts). We also assembled the sequences using MIRA and PartiGene software and compared the three unigene sets. Gene ontology annotation was then assigned to 29,303 unigene elements. Blast search against the SWISS-PROT database revealed putative homologs for 32,810 (14.7%) unigene elements, but more extensive search with Pfam, Refseq_protein, Refseq_RNA and eight gene indices revealed homology for 67.4% of them. The EST catalogue was examined for putative homologs of candidate genes involved in bud phenology, cuticle formation, phenylpropanoids biosynthesis and cell wall formation. Our results suggest a good coverage of genes involved in these traits. Comparative orthologous sequences (COS) with other plant gene models were identified and allow to unravel the oak paleo-history. Simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were searched, resulting in 52,834 SSRs and 36,411 SNPs. All of these are available through the Oak Contig Browser http://genotoul-contigbrowser.toulouse.inra.fr:9092/Quercus_robur/index.html.ConclusionsThis genomic resource provides a unique tool to discover genes of interest, study the oak transcriptome, and develop new markers to investigate functional diversity in natural populations.

A micromethod for high throughput RNA extraction in forest trees.

A large quantity of high quality RNA is often required in the analysis of gene expression. However, RNA extraction from samples taken from woody plants is generally complex, and represents the main limitation to study gene expression, particularly in refractory species like conifers. Standard RNA extraction protocols are available but they are highly time consuming, and not adapted to large scale extraction. Here we present a high-throughput RNA extraction protocol. This protocol was adapted to a micro-scale by modifying the classical cetyltrimethylammonium (CTAB) protocol developed for pine: (i) quantity of material used (100-200 mg of sample), (ii) disruption of samples in microtube using a mechanical tissue disrupter, and (iii) the use of SSTE buffer. One hundred samples of woody plant tissues/organs can be easily treated in two working days. An average of 15 \ig of high quality RNA per sample was obtained. The RNA extracted is suitable for applications such as real time reverse transcription polymerase chain reaction, cDNA library construction or synthesis of complex targets for microarray analysis.

Development and implementation of a highly- multiplexed SNP array for genetic mapping in maritime pine and comparative mapping with loblolly pine

BackgroundSingle nucleotide polymorphisms (SNPs) are the most abundant source of genetic variation among individuals of a species. New genotyping technologies allow examining hundreds to thousands of SNPs in a single reaction for a wide range of applications such as genetic diversity analysis, linkage mapping, fine QTL mapping, association studies, marker-assisted or genome-wide selection. In this paper, we evaluated the potential of highly-multiplexed SNP genotyping for genetic mapping in maritime pine (Pinus pinaster Ait.), the main conifer used for commercial plantation in southwestern Europe.ResultsWe designed a custom GoldenGate assay for 1,536 SNPs detected through the resequencing of gene fragments (707 in vitro SNPs/Indels) and from Sanger-derived Expressed Sequenced Tags assembled into a unigene set (829 in silico SNPs/Indels). Offspring from three-generation outbred (G2) and inbred (F2) pedigrees were genotyped. The success rate of the assay was 63.6% and 74.8% for in silico and in vitro SNPs, respectively. A genotyping error rate of 0.4% was further estimated from segregating data of SNPs belonging to the same gene. Overall, 394 SNPs were available for mapping. A total of 287 SNPs were integrated with previously mapped markers in the G2 parental maps, while 179 SNPs were localized on the map generated from the analysis of the F2 progeny. Based on 98 markers segregating in both pedigrees, we were able to generate a consensus map comprising 357 SNPs from 292 different loci. Finally, the analysis of sequence homology between mapped markers and their orthologs in a Pinus taeda linkage map, made it possible to align the 12 linkage groups of both species.ConclusionsOur results show that the GoldenGate assay can be used successfully for high-throughput SNP genotyping in maritime pine, a conifer species that has a genome seven times the size of the human genome. This SNP-array will be extended thanks to recent sequencing effort using new generation sequencing technologies and will include SNPs from comparative orthologous sequences that were identified in the present study, providing a wider collection of anchor points for comparative genomics among the conifers.

Transcriptional profiling of bud dormancy induction and release in oak by next-generation sequencing

BackgroundIn temperate regions, the time lag between vegetative bud burst and bud set determines the duration of the growing season of trees (i.e. the duration of wood biomass production). Dormancy, the period during which the plant is not growing, allows trees to avoid cold injury resulting from exposure to low temperatures. An understanding of the molecular machinery controlling the shift between these two phenological states is of key importance in the context of climatic change. The objective of this study was to identify genes upregulated during endo- and ecodormancy, the two main stages of bud dormancy. Sessile oak is a widely distributed European white oak species. A forcing test on young trees was first carried out to identify the period most likely to correspond to these two stages. Total RNA was then extracted from apical buds displaying endo- and ecodormancy. This RNA was used for the generation of cDNA libraries, and in-depth transcriptome characterization was performed with 454 FLX pyrosequencing technology.ResultsPyrosequencing produced a total of 495,915 reads. The data were cleaned, duplicated reads removed, and sequences were mapped onto the oak UniGene data. Digital gene expression analysis was performed, with both R statistics and the R-Bioconductor packages (edgeR and DESeq), on 6,471 contigs with read numbers ≥ 5 within any contigs. The number of sequences displaying significant differences in expression level (read abundance) between endo- and ecodormancy conditions ranged from 75 to 161, depending on the algorithm used. 13 genes displaying significant differences between conditions were selected for further analysis, and 11 of these genes, including those for glutathione-S-transferase (GST) and dehydrin xero2 (XERO2) were validated by quantitative PCR.ConclusionsThe identification and functional annotation of differentially expressed genes involved in the “response to abscisic acid”, “response to cold stress” and “response to oxidative stress” categories constitutes a major step towards characterization of the molecular network underlying vegetative bud dormancy, an important life history trait of long-lived organisms.

Decoding the oak genome: public release of sequence data, assembly, annotation and publication strategies

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole‐genome shotgun (WGS) approach, without the use of costly and time‐consuming methods, such as fosmid or BAC clone‐based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS‐FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired‐end and mate‐pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired‐end reads and contaminants detected, resulting in a total of 17 910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.

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

BackgroundThe availability of a large expressed sequence tags (EST) resource and recent advances in high-throughput genotyping technology have made it possible to develop highly multiplexed SNP arrays for multi-objective genetic applications, including the construction of meiotic maps. Such approaches are particularly useful in species with a large genome size, precluding the use of whole-genome shotgun assembly with current technologies.ResultsIn this study, a 12 k-SNP genotyping array was developed for maritime pine from an extensive EST resource assembled into a unigene set. The offspring of three-generation outbred and inbred mapping pedigrees were then genotyped. The inbred pedigree consisted of a classical F2 population resulting from the selfing of a single inter-provenance (Landes x Corsica) hybrid tree, whereas the outbred pedigree (G2) resulted from a controlled cross of two intra-provenance (Landes x Landes) hybrid trees. This resulted in the generation of three linkage maps based on SNP markers: one from the parental genotype of the F2 population (1,131 markers in 1,708 centimorgan (cM)), and one for each parent of the G2 population (1,015 and 1,110 markers in 1,447 and 1,425 cM for the female and male parents, respectively). A comparison of segregation patterns in the progeny obtained from the two types of mating (inbreeding and outbreeding) led to the identification of a chromosomal region carrying an embryo viability locus with a semi-lethal allele. Following selfing and segregation, zygote mortality resulted in a deficit of Corsican homozygous genotypes in the F2 population. This dataset was also used to study the extent and distribution of meiotic recombination along the length of the chromosomes and the effect of sex and/or genetic background on recombination. The genetic background of trees in which meiotic recombination occurred was found to have a significant effect on the frequency of recombination. Furthermore, only a small proportion of the recombination hot- and cold-spots were common to all three genotypes, suggesting that the spatial pattern of recombination was genetically variable.ConclusionThis study led to the development of classical genomic tools for this ecologically and economically important species. It also identified a chromosomal region bearing a semi-lethal recessive allele and demonstrated the genetic variability of recombination rate over the genome.

Molecular and phenotypic profiling from the base to the crown in maritime pine wood‐forming tissue

Environmental, developmental and genetic factors affect variation in wood properties at the chemical, anatomical and physical levels. Here, the phenotypic variation observed along the tree stem was explored and the hypothesis tested that this variation could be the result of the differential expression of genes/proteins during wood formation. Differentiating xylem samples of maritime pine (Pinus pinaster) were collected from the top (crown wood, CW) to the bottom (base wood, BW) of adult trees. These samples were characterized by Fourier transform infrared spectroscopy (FTIR) and analytical pyrolysis. Two main groups of samples, corresponding to CW and BW, could be distinguished from cell wall chemical composition. A genomic approach, combining large-scale production of expressed sequence tags (ESTs), gene expression profiling and quantitative proteomics analysis, allowed identification of 262 unigenes (out of 3512) and 231 proteins (out of 1372 spots) that were differentially expressed along the stem. A good relationship was found between functional categories from transcriptomic and proteomic data. A good fit between the molecular mechanisms involved in CW-BW formation and these two types of wood phenotypic differences was also observed. This work provides a list of candidate genes for wood properties that will be tested in forward genetics.

Analysis of BAC end sequences in oak, a keystone forest tree species, providing insight into the composition of its genome

BackgroundOne of the key goals of oak genomics research is to identify genes of adaptive significance. This information may help to improve the conservation of adaptive genetic variation and the management of forests to increase their health and productivity. Deep-coverage large-insert genomic libraries are a crucial tool for attaining this objective. We report herein the construction of a BAC library for Quercus robur, its characterization and an analysis of BAC end sequences.ResultsThe Eco RI library generated consisted of 92,160 clones, 7% of which had no insert. Levels of chloroplast and mitochondrial contamination were below 3% and 1%, respectively. Mean clone insert size was estimated at 135 kb. The library represents 12 haploid genome equivalents and, the likelihood of finding a particular oak sequence of interest is greater than 99%. Genome coverage was confirmed by PCR screening of the library with 60 unique genetic loci sampled from the genetic linkage map. In total, about 20,000 high-quality BAC end sequences (BESs) were generated by sequencing 15,000 clones. Roughly 5.88% of the combined BAC end sequence length corresponded to known retroelements while ab initio repeat detection methods identified 41 additional repeats. Collectively, characterized and novel repeats account for roughly 8.94% of the genome. Further analysis of the BESs revealed 1,823 putative genes suggesting at least 29,340 genes in the oak genome. BESs were aligned with the genome sequences of Arabidopsis thaliana, Vitis vinifera and Populus trichocarpa. One putative collinear microsyntenic region encoding an alcohol acyl transferase protein was observed between oak and chromosome 2 of V. vinifera.ConclusionsThis BAC library provides a new resource for genomic studies, including SSR marker development, physical mapping, comparative genomics and genome sequencing. BES analysis provided insight into the structure of the oak genome. These sequences will be used in the assembly of a future genome sequence for oak.

Role of waterlogging-responsive genes in shaping interspecific differentiation between two sympatric oak species

Pedunculate (Quercus robur L.) and sessile oak (Quercus petreae Matt. Liebl.) are closely related species with a widely sympatric distribution in Europe. These two oak species are also known to display different ecological features, particularly related to their adaptation to soil waterlogging. Pedunculate oak grows in humid areas and can withstand high moisture content of the soil, whereas sessile oak requires drier soil with better drainage. The main goal of this study was to explore the role of gene expression contributing to differences in terms of waterlogging tolerance between these two species. We implemented a series of experiments aimed at evaluating whether differentially expressed genes between species are associated with their ecological preferences and underlie adaptive genetic divergence. Rooted cuttings of both species were grown in hydroponic conditions and subjected to gradual root hypoxia. White roots were sampled after 6, 12, 24 and 48 h. Real-time polymerase chain reaction (qPCR) was first used to monitor the expression of 10 known waterlogging-responsive genes, to identify discriminating sampling time points along the kinetics of hypoxia. Secondly, four subtractive suppressive hybridization libraries (sessile vs. pedunculate, pedunculate vs. sessile for early and late responses) were generated to isolate differentially expressed genes between species. A total of 2160 high-quality expressed sequence tags were obtained and annotated, and a subset of 45 genes were selected for qPCR analysis in a second independent factorial experimental design applying two stress durations per two species. Significant differences of gene expression between pedunculate and sessile oaks were detected, suggesting species-specific molecular strategies to respond to hypoxia. This study revealed that the ability of pedunculate oak to maintain glycolysis and fermentation under hypoxic conditions may help explain its tolerance to waterlogging.