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.

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 (CGs) (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 times of bud burst observed in common garden experiments (provenance tests). The CGs 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 (pi(tot)=6.15 x 10(-3); pi(silent)=11.2 x 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 CGs 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.

Contrasting relations between 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.

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.

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 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.

Inheritance and molecular variations of PCR-SSCP fragments in pedunculate oak (Quercus robur L.).

Single-strand conformaiton polymorphism (SSCP) profiles of six PCR-amplified fragments (250-800 bp) were analyzed in three full-sib families of pedunculate oak (Quercus robur L.) and their parents. Among the six fragments, four were polymorphic and one exhibited complex patterns that were not changed by varying the SSCP conditions. The number of bands for the analyzed fragments varied between two and four among individuals regardless of fragment size. As shown by segregation data, the variation in the number of bands between trees could only be attributed to the allelic composition (homozygotes vs heterozygotes): a genotype that exhibited two bands was presumptively homozygous, wheras a genotype exhibiting three or four bands was heterozygous. Mendelian proportions were observed in all crosses for each polymorphic fragment. In one cross, we could clearly identify a null allele due to a possible mutation at a primer site. Single-base mutations and short insertion-deletions were shown to be the molecular causes of the SSCP polymorphism observed between different alleles. The use of SSCP as a technique to identify co-dominant markers of PCR fragments (up to 800 bp) is recommended for gene diversity studies or for gene mapping.