Catherine Bodénès

Comparison of microsatellites and amplified fragment length polymorphism markers for parentage analysis

This study compares the properties of dominant markers, such as amplified fragment length polymorphisms (AFLPs), with those of codominant multiallelic markers, such as microsatellites, in reconstructing parentage. These two types of markers were used to search for both parents of an individual without prior knowledge of their relationships, by calculating likelihood ratios based on genotypic data, including mistyping. Experimental data on 89 oak trees genotyped for six microsatellite markers and 159 polymorphic AFLP loci were used as a starting point for simulations and tests. Both sets of markers produced high exclusion probabilities, and among dominant markers those with dominant allele frequencies in the range 0.1–0.4 were more informative. Such codominant and dominant markers can be used to construct powerful statistical tests to decide whether a genotyped individual (or two individuals) can be considered as the true parent (or parent pair). Gene flow from outside the study stand (GFO), inferred from parentage analysis with microsatellites, overestimated the true GFO, whereas with AFLPs it was underestimated. As expected, dominant markers are less efficient than codominant markers for achieving this, but can still be used with good confidence, especially when loci are deliberately selected according to their allele frequencies.

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

Abstract.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® (low-acid peach) × cv. Fantasia (normally-acid nectarine) F2 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® × Fantasia population, gene mapping was performed on the almond [Prunus amygdalus (P. dulcis)] cv. Texas × peach cv. Earlygold F2 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® × Fantasia and Texas × Earlygold maps and data from a third year of trait assessment on the Ferjalou Jalousia® × 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.

A genetic linkage map of Quercus robur L. (pedunculate oak) based on RAPD, SCAR, microsatellite, minisatellite, isozyme and 5S rDNA markers

Abstract A genetic map of Pedunculate oak (Quercus robur) was constructed based on one 5S rDNA, 271 RAPD, ten SCAR, 18 microsatellite, one minisatellite, and six isozyme markers. A total of 94 individuals from a full-sib family was genotyped. Two maps, including 307 markers, were constructed according to the “two-way pseudo-testcross” mapping strategy. Testcross markers segregating in the 1 : 1 ratio were first used to establish separate maternal (893.2 cM, 12 linkage groups) and paternal (921.7 cM, 12 linkage groups) maps. Both maps provided 85–90% genome coverage. Homologies between the male and female linkage groups were then identified based on 74 intercross markers segregating in the 3 : 1, 1 : 2 : 1 and 1 : 1 : 1 : 1 ratios (RAPDs, SCARs, SSRs, 5S rDNA and isozymes) in the hybrid progeny. In each map, approximately 18% of the studied markers showed segregation distortion. More than 60% of the skewed markers were due to an excess of heterozygote genotypes. This map will be used for: (1) studying the molecular organisation of genomic regions involved in inter- and intraspecific differentiation in oaks and (2) identification of QTLs for adaptive traits.

A fast and cost-effective approach to develop and map EST-SSR markers: oak as a case study

BackgroundExpressed Sequence Tags (ESTs) are a source of simple sequence repeats (SSRs) that can be used to develop molecular markers for genetic studies. The availability of ESTs for Quercus robur and Quercus petraea provided a unique opportunity to develop microsatellite markers to accelerate research aimed at studying adaptation of these long-lived species to their environment. As a first step toward the construction of a SSR-based linkage map of oak for quantitative trait locus (QTL) mapping, we describe the mining and survey of EST-SSRs as well as a fast and cost-effective approach (bin mapping) to assign these markers to an approximate map position. We also compared the level of polymorphism between genomic and EST-derived SSRs and address the transferability of EST-SSRs in Castanea sativa (chestnut).ResultsA catalogue of 103,000 Sanger ESTs was assembled into 28,024 unigenes from which 18.6% presented one or more SSR motifs. More than 42% of these SSRs corresponded to trinucleotides. Primer pairs were designed for 748 putative unigenes. Overall 37.7% (283) were found to amplify a single polymorphic locus in a reference full-sib pedigree of Quercus robur. The usefulness of these loci for establishing a genetic map was assessed using a bin mapping approach. Bin maps were constructed for the male and female parental tree for which framework linkage maps based on AFLP markers were available. The bin set consisting of 14 highly informative offspring selected based on the number and position of crossover sites. The female and male maps comprised 44 and 37 bins, with an average bin length of 16.5 cM and 20.99 cM, respectively. A total of 256 EST-SSRs were assigned to bins and their map position was further validated by linkage mapping. EST-SSRs were found to be less polymorphic than genomic SSRs, but their transferability rate to chestnut, a phylogenetically related species to oak, was higher.ConclusionWe have generated a bin map for oak comprising 256 EST-SSRs. This resource constitutes a first step toward the establishment of a gene-based map for this genus that will facilitate the dissection of QTLs affecting complex traits of ecological importance.

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.

Contrasting relationships between the diversity of candidate genes and variation of bud burst in natural and segregating populations of European oaks

Nucleotide diversity was assessed within nine candidate genes (in total 4.6 kb) for the time of bud burst in nine sessile oak (Quercus petraea) populations distributed in central and northern Europe. The sampled populations were selected on the basis of their contrasting time of bud burst observed in common garden experiments (provenance tests). The candidate genes were selected according to their expression profiles during the transition from quiescent to developing buds and/or their functional role in model plants. The overall nucleotide diversity was large (πtot=6.15 × 10−3; πsilent=11.2 × 10−3), but population differentiation was not larger than for microsatellites. No outlier single-nucleotide polymorphism (SNP), departing from neutral expectation, was found among the total of 125 SNPs. These results contrasted markedly with the significant associations that were observed between the candidate genes and bud burst in segregating populations. Quantitative trait loci (QTLs) for bud burst were identified for 13 year*site seasonal observations in a cloned mapping pedigree. Nineteen QTLs were detected, and QTLs located on linkage groups 2, 5 and 9 contributed repeatedly to more than 12% of the phenotypic variation of the trait. Eight genes were polymorphic in the two parents of the pedigree and could be mapped on the existing genetic map. Five of them located within the confidence intervals of QTLs for bud burst. Interestingly, four of them located within the three QTLs exhibiting the largest contributions to bud burst.

A genetic map of maritime pine based on AFLP, RAPD and protein markers

Abstract TheAFLP (amplified fragment length polymorphism) technique was adapted to carry out genetic analysis in maritime pine, a species characterized by a large genome size (24 pg/C). A genetic linkage map was constructed for one F1 individual based on 239 AFLP and 127 RAPD (randomly amplified polymorphic DNA) markers. Markers were scored on megagametophytes (1n) from 200 germinated F2 seedlings. Polymorphism rate, labour time and cost of both AFLP and RAPD techniques were compared. The AFLP technique was found to be twice as fast and three-times less costly per marker than the RAPD technique. Thirteen linkage groups were identified with a LOD score ≥6 covering 1873 cM, which provided 93.4% of genome coverage. Proteins were extracted from needles (2n) of the F2 progeny and revealed by 2-DE (two-dimensional electrophoresis). Thirty one segregating proteins were mapped using a QTL detection strategy based on the quantification of protein accumulation. Two framework maps of the same F1 individual are now available. The first map (Plomion et al. 1996) uses RAPD markers and the second map, presented in this study, uses mostly AFLP markers. Although the total genetic length of both maps was almost identical, differences among homologous groups were observed.

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.

Comparative mapping in the Fagaceae and beyond with EST-SSRs

BackgroundGenetic markers and linkage mapping are basic prerequisites for comparative genetic analyses, QTL detection and map-based cloning. A large number of mapping populations have been developed for oak, but few gene-based markers are available for constructing integrated genetic linkage maps and comparing gene order and QTL location across related species.ResultsWe developed a set of 573 expressed sequence tag-derived simple sequence repeats (EST-SSRs) and located 397 markers (EST-SSRs and genomic SSRs) on the 12 oak chromosomes (2n = 2x = 24) on the basis of Mendelian segregation patterns in 5 full-sib mapping pedigrees of two species: Quercus robur (pedunculate oak) and Quercus petraea (sessile oak). Consensus maps for the two species were constructed and aligned. They showed a high degree of macrosynteny between these two sympatric European oaks. We assessed the transferability of EST-SSRs to other Fagaceae genera and a subset of these markers was mapped in Castanea sativa, the European chestnut. Reasonably high levels of macrosynteny were observed between oak and chestnut. We also obtained diversity statistics for a subset of EST-SSRs, to support further population genetic analyses with gene-based markers. Finally, based on the orthologous relationships between the oak, Arabidopsis, grape, poplar, Medicago, and soybean genomes and the paralogous relationships between the 12 oak chromosomes, we propose an evolutionary scenario of the 12 oak chromosomes from the eudicot ancestral karyotype.ConclusionsThis study provides map locations for a large set of EST-SSRs in two oak species of recognized biological importance in natural ecosystems. This first step toward the construction of a gene-based linkage map will facilitate the assignment of future genome scaffolds to pseudo-chromosomes. This study also provides an indication of the potential utility of new gene-based markers for population genetics and comparative mapping within and beyond the Fagaceae.

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.