Pauline Garnier-Géré

Automated SNP detection in expressed sequence tags: statistical considerations and application to maritime pine sequences.

We developed an automated pipeline for the detection of single nucleotide polymorphisms (SNPs) in expressed sequence tag (EST) data sets, by combining three DNA sequence analysis programs:Phred, Phrap and PolyBayes. This application requires access to the individual electrophoregram traces. First, a reference set of 65 SNPs was obtained from the sequencing of 30 gametes in 13 maritime pine (Pinus pinaster Ait.) gene fragments (6671 bp), resulting in a frequency of 1 SNP every 102.6 bp. Second, parameters of the three programs were optimized in order to retrieve as many true SNPs, while keeping the rate of false positive as low as possible. Overall, the efficiency of detection of true SNPs was 83.1%. However, this rate varied largely as a function of the rare SNP allele frequency: down to 41% for rare SNP alleles (frequency ` 10%), up to 98% for allele frequencies above 10%. Third, the detection method was applied to the 18498 assembled maritime pine (Pinus pinaster Ait.) ESTs, allowing to identify a total of 1400 candidate SNPs, in contigs containing between 4 and 20 sequence reads. These genetic resources, described for the first time in a forest tree species, were made available at http://www.pierroton.inra/genetics/Pinesnps. We also derived an analytical expression for the SNP detection probability as a function of the SNP allele frequency, the number of haploid genomes used to generate the EST sequence database, and the sample size of the contigs considered for SNP detection. The frequency of the SNP allele was shown to be the main factor influencing the probability of SNP detection.

Influences of triploidy, parentage and genetic diversity on growth of the Pacific oyster Crassostrea gigas reared in contrasting natural environments.

An increasing number of hypotheses are being proposed to explain the faster growth potential of triploids in molluscs, including their partial sterility or their higher heterozygosity compared to diploids. Triploid advantage however, remains controversial for poorer sites, because of a potential trade‐off with survival. These questions were addressed in Crassostrea gigas by deploying meiosis II triploids and their diploid siblings from a single mass spawning of three males and seven females, in two contrasting locations for their trophic resources. One hundred and fifty individuals were sampled at each site after nine months, measured for weight and biochemical composition, and genotyped using three microsatellite and seven allozyme loci. Higher performance was observed at the fast‐growing site for all traits except shell weight, and triploids had greater weights and biochemical contents than diploids at harvest. Triploids also grew faster at the poorer site, and showed similar survival rates to diploids at both sites. Triploids had significantly higher average allozyme and microsatellite diversity. However, they performed better for a wide range of individual heterozygosity values, arguing for an advantage of the triploid state per se, that could be due to positive effects on growth of both sterility of triploids with subsequent resource re‐allocation and possible faster transcription with three copies of each gene. Despite evidence of very low or no inbreeding in the diploid sample, positive associations between individual allozyme diversity and growth were detected, which explained little but significant amounts of phenotypic variation. These associations were interpreted as direct effects of allozymes, either alone or including epistatic interactions with other loci. In addition, measures of individual distance (mean‐d2) specific to microsatellites, were negatively correlated with growth in diploids, indicating possible effects of outbreeding depression between more distant genomes of parents from distinct populations.

Standardizing for microsatellite length in comparisons of genetic diversity

Mutation rates at microsatellites tend to increase with the number of repeats of the motif, leading to higher levels of polymorphism at long microsatellites. To standardize levels of diversity when microsatellites differ in size, we investigate the relationship between tract length and variation and provide a formula to adjust allelic richness to a fixed mean number of repeats in the specific case of chloroplast microsatellites. A comparison between 39 loci from eight species of conifers (where chloroplast DNA is paternally inherited) and 64 loci from 12 species of angiosperms (where chloroplast DNA is generally predominantly maternally inherited) indicates that the greater allelic richness found in conifers remains significant after controlling for number of repeats. The approach stresses the advantage of reporting variation in number of repeats instead of relative fragment sizes.

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.

High-quality SNPs from genic regions highlight introgression patterns among European white oaks (Quercus petraea and Q. robur).

In the post-genomics era, non-model species like most Fagaceae still lack operational diversity resources for population genomics studies. Sequence data were produced from over 800 gene fragments covering ~530 kb across the genic partition of European oaks, in a discovery panel of 25 individuals from western and central Europe (11 Quercus petraea, 13 Q. robur, one Q. ilex as an outgroup). Regions targeted represented broad functional categories potentially involved in species ecological preferences, and a random set of genes. Using a high-quality dedicated pipeline, we provide a detailed characterization of these genic regions, which included over 14500 polymorphisms, with ~12500 SNPs −218 being triallelic-, over 1500 insertion-deletions, and ~200 novel di- and tri-nucleotide SSR loci. This catalog also provides various summary statistics within and among species, gene ontology information, and standard formats to assist loci choice for genotyping projects. The distribution of nucleotide diversity (θπ) and differentiation (FST) across genic regions are also described for the first time in those species, with a mean n θπ close to ~0.0049 in Q. petraea and to ~0.0045 in Q. robur across random regions, and a mean FST ~0.13 across SNPs. The magnitude of diversity across genes is within the range estimated for long-term perennial outcrossers, and can be considered relatively high in the plant kingdom, with an estimate across the genome of 41 to 51 million SNPs expected in both species. Individuals with typical species morphology were more easily assigned to their corresponding genetic cluster for Q. robur than for Q. petraea, revealing higher or more recent introgression in Q. petraea and a stronger species integration in Q. robur in this particular discovery panel. We also observed robust patterns of a slightly but significantly higher diversity in Q. petraea, across a random gene set and in the abiotic stress functional category, and a heterogeneous landscape of both diversity and differentiation. To explain these patterns, we discuss an alternative and non-exclusive hypothesis of stronger selective constraints in Q. robur, the most pioneering species in oak forest stand dynamics, additionally to the recognized and documented introgression history in both species despite their strong reproductive barriers. The quality of the data provided here and their representativity in terms of species genomic diversity make them useful for possible applications in medium-scale landscape and molecular ecology projects. Moreover, they can serve as reference resources for validation purposes in larger-scale resequencing projects. This type of project is preferentially recommended in oaks in contrast to SNP array development, given the large nucleotide variation and the low levels of linkage disequilibrium revealed.

Species-specific alleles at a β-tubulin gene show significant associations with leaf morphological variation within Quercus petraea and Q. robur populations

Quercus petraea and Q. robur are largely sympatric oak species in western and central Europe and known for their intensive genetic exchange which has made the discovery of species-diagnostic markers a huge challenge. Various natural white oak populations (Q. petraea/Q. robur including mixed stands) were investigated for their variability and differentiation patterns at a β-tubulin gene (qutub8) in a European-wide survey. This gene was chosen as a possible candidate among loci subjected to selection and maintaining integrity between species. Two frequent alleles depicted as indels within qutub8’s first intron showed remarkably high interspecific genetic differentiation, with Weir and Cockerham’s theta per allele values ranging from 0.17 to 0.30 for one allele and from 0.04 to 0.19 for the other allele in such mixed oak stands where the multi-allelic qutub8 locus showed significant interspecific FST. For three mixed stands, qutub8’s FST significantly departed from the expected neutral differentiation patterns (FST ranging from 0.063 to 0.080 for this multi-allelic marker) and thus could be influenced by selection. Significant associations were found between genotypic variation and leaf dimensions as well as leaf structure patterns, after having accounted for species and stand effects. Qutub8 represents a locus that exhibits significant species differentiation and is linked to morphological discriminant traits. Consequently, qutub8 likely contributes to species divergence within the European white oak complex.