Loïc Le Cunff

Quantitative Genetic Bases of Anthocyanin Variation in Grape ( Vitis vinifera L. ssp. sativa ) Berry: A Quantitative Trait Locus to Quantitative Trait Nucleotide Integrated Study

The combination of QTL mapping studies of synthetic lines and association mapping studies of natural diversity represents an opportunity to throw light on the genetically based variation of quantitative traits. With the positional information provided through quantitative trait locus (QTL) mapping, which often leads to wide intervals encompassing numerous genes, it is now feasible to directly target candidate genes that are likely to be responsible for the observed variation in completely sequenced genomes and to test their effects through association genetics. This approach was performed in grape, a newly sequenced genome, to decipher the genetic architecture of anthocyanin content. Grapes may be either white or colored, ranging from the lightest pink to the darkest purple tones according to the amount of anthocyanin accumulated in the berry skin, which is a crucial trait for both wine quality and human nutrition. Although the determinism of the white phenotype has been fully identified, the genetic bases of the quantitative variation of anthocyanin content in berry skin remain unclear. A single QTL responsible for up to 62% of the variation in the anthocyanin content was mapped on a Syrah × Grenache F1 pseudo-testcross. Among the 68 unigenes identified in the grape genome within the QTL interval, a cluster of four Myb-type genes was selected on the basis of physiological evidence (VvMybA1, VvMybA2, VvMybA3, and VvMybA4). From a core collection of natural resources (141 individuals), 32 polymorphisms revealed significant association, and extended linkage disequilibrium was observed. Using a multivariate regression method, we demonstrated that five polymorphisms in VvMybA genes except VvMybA4 (one retrotransposon, three single nucleotide polymorphisms and one 2-bp insertion/deletion) accounted for 84% of the observed variation. All these polymorphisms led to either structural changes in the MYB proteins or differences in the VvMybAs promoters. We concluded that the continuous variation in anthocyanin content in grape was explained mainly by a single gene cluster of three VvMybA genes. The use of natural diversity helped to reduce one QTL to a set of five quantitative trait nucleotides and gave a clear picture of how isogenes combined their effects to shape grape color. Such analysis also illustrates how isogenes combine their effect to shape a complex quantitative trait and enables the definition of markers directly targeted for upcoming breeding programs.

Construction of nested genetic core collections to optimize the exploitation of natural diversity in Vitis vinifera L. subsp. sativa

BackgroundThe first high quality draft of the grape genome sequence has just been published. This is a critical step in accessing all the genes of this species and increases the chances of exploiting the natural genetic diversity through association genetics. However, our basic knowledge of the extent of allelic variation within the species is still not sufficient. Towards this goal, we constructed nested genetic core collections (G-cores) to capture the simple sequence repeat (SSR) diversity of the grape cultivated compartment (Vitis vinifera L. subsp. sativa) from the worlds largest germplasm collection (Domaine de Vassal, INRA Hérault, France), containing 2262 unique genotypes.ResultsSub-samples of 12, 24, 48 and 92 varieties of V. vinifera L. were selected based on their genotypes for 20 SSR markers using the M-strategy. They represent respectively 58%, 73%, 83% and 100% of total SSR diversity. The capture of allelic diversity was analyzed by sequencing three genes scattered throughout the genome on 233 individuals: 41 single nucleotide polymorphisms (SNPs) were identified using the G-92 core (one SNP for every 49 nucleotides) while only 25 were observed using a larger sample of 141 individuals selected on the basis of 50 morphological traits, thus demonstrating the reliability of the approach.ConclusionThe G-12 and G-24 core-collections displayed respectively 78% and 88% of the SNPs respectively, and are therefore of great interest for SNP discovery studies. Furthermore, the nested genetic core collections satisfactorily reflected the geographic and the genetic diversity of grape, which are also of great interest for the study of gene evolution in this species.

A candidate gene association study on muscat flavor in grapevine (Vitis vinifera L.)

BackgroundThe sweet, floral flavor typical of Muscat varieties (Muscats), due to high levels of monoterpenoids (geraniol, linalool and nerol), is highly distinct and has been greatly appreciated both in table grapes and in wine since ancient times. Muscat flavor determination in grape (Vitis vinifera L.) has up to now been studied by evaluating monoterpenoid levels through QTL analysis. These studies have revealed co-localization of 1-deoxy-D-xylulose 5-phosphate synthase (VvDXS) with the major QTL positioned on chromosome 5.ResultsWe resequenced VvDXS in an ad hoc association population of 148 grape varieties, which included muscat-flavored, aromatic and neutral accessions as well as muscat-like aromatic mutants and non-aromatic offsprings of Muscats. Gene nucleotide diversity and intragenic linkage disequilibrium (LD) were evaluated. Structured association analysis revealed three SNPs in moderate LD to be significantly associated with muscat-flavored varieties. We identified a putative causal SNP responsible for a predicted non-neutral substitution and we discuss its possible implications for flavor metabolism. Network analysis revealed a major star-shaped cluster of reconstructed haplotypes unique to muscat-flavored varieties. Moreover, muscat-like aromatic mutants displayed unique non-synonymous mutations near the mutated site of Muscat genotypes.ConclusionsThis study is a crucial step forward in understanding the genetic regulation of muscat flavor in grapevine and it also sheds light on the domestication history of Muscats. VvDXS appears to be a possible human-selected locus in grapevine domestication and post-domestication. The putative causal SNP identified in Muscat varieties as well as the unique mutations identifying the muscat-like aromatic mutants under study may be immediately applied in marker-assisted breeding programs aimed at enhancing fragrance and aroma complexity respectively in table grape and wine cultivars.

Molecular, genetic and transcriptional evidence for a role of VvAGL11 in stenospermocarpic seedlessness in grapevine

BackgroundStenospermocarpy is a mechanism through which certain genotypes of Vitis vinifera L. such as Sultanina produce berries with seeds reduced in size. Stenospermocarpy has not yet been characterized at the molecular level.ResultsGenetic and physical maps were integrated with the public genomic sequence of Vitis vinifera L. to improve QTL analysis for seedlessness and berry size in experimental progeny derived from a cross of two seedless genotypes. Major QTLs co-positioning for both traits on chromosome 18 defined a 92-kb confidence interval. Functional information from model species including Vitis suggested that VvAGL11, included in this confidence interval, might be the main positional candidate gene responsible for seed and berry development.Characterization of VvAGL11 at the sequence level in the experimental progeny identified several SNPs and INDELs in both regulatory and coding regions. In association analyses performed over three seasons, these SNPs and INDELs explained up to 78% and 44% of the phenotypic variation in seed and berry weight, respectively. Moreover, genetic experiments indicated that the regulatory region has a larger effect on the phenotype than the coding region. Transcriptional analysis lent additional support to the putative role of VvAGL11s regulatory region, as its expression is abolished in seedless genotypes at key stages of seed development. These results transform VvAGL11 into a functional candidate gene for further analyses based on genetic transformation.For breeding purposes, intragenic markers were tested individually for marker assisted selection, and the best markers were those closest to the transcription start site.ConclusionWe propose that VvAGL11 is the major functional candidate gene for seedlessness, and we provide experimental evidence suggesting that the seedless phenotype might be caused by variations in its promoter region. Current knowledge of the function of its orthologous genes, its expression profile in Vitis varieties and the strong association between its sequence variation and the degree of seedlessness together indicate that the D-lineage MADS-box gene VvAGL11 corresponds to the Seed Development Inhibitor locus described earlier as a major locus for seedlessness. These results provide new hypotheses for further investigations of the molecular mechanisms involved in seed and berry development.

SNiPlay: a web-based tool for detection, management and analysis of SNPs. Application to grapevine diversity projects

BackgroundHigh-throughput re-sequencing, new genotyping technologies and the availability of reference genomes allow the extensive characterization of Single Nucleotide Polymorphisms (SNPs) and insertion/deletion events (indels) in many plant species. The rapidly increasing amount of re-sequencing and genotyping data generated by large-scale genetic diversity projects requires the development of integrated bioinformatics tools able to efficiently manage, analyze, and combine these genetic data with genome structure and external data.ResultsIn this context, we developed SNiPlay, a flexible, user-friendly and integrative web-based tool dedicated to polymorphism discovery and analysis. It integrates:1) a pipeline, freely accessible through the internet, combining existing softwares with new tools to detect SNPs and to compute different types of statistical indices and graphical layouts for SNP data. From standard sequence alignments, genotyping data or Sanger sequencing traces given as input, SNiPlay detects SNPs and indels events and outputs submission files for the design of Illuminas SNP chips. Subsequently, it sends sequences and genotyping data into a series of modules in charge of various processes: physical mapping to a reference genome, annotation (genomic position, intron/exon location, synonymous/non-synonymous substitutions), SNP frequency determination in user-defined groups, haplotype reconstruction and network, linkage disequilibrium evaluation, and diversity analysis (Pi, Wattersons Theta, Tajimas D).Furthermore, the pipeline allows the use of external data (such as phenotype, geographic origin, taxa, stratification) to define groups and compare statistical indices.2) a database storing polymorphisms, genotyping data and grapevine sequences released by public and private projects. It allows the user to retrieve SNPs using various filters (such as genomic position, missing data, polymorphism type, allele frequency), to compare SNP patterns between populations, and to export genotyping data or sequences in various formats.ConclusionsOur experiments on grapevine genetic projects showed that SNiPlay allows geneticists to rapidly obtain advanced results in several key research areas of plant genetic diversity. Both the management and treatment of large amounts of SNP data are rendered considerably easier for end-users through automation and integration. Current developments are taking into account new advances in high-throughput technologies.SNiPlay is available at: http://sniplay.cirad.fr/.

Dissecting genetic architecture of grape proanthocyanidin composition through quantitative trait locus mapping

BackgroundProanthocyanidins (PAs), or condensed tannins, are flavonoid polymers, widespread throughout the plant kingdom, which provide protection against herbivores while conferring organoleptic and nutritive values to plant-derived foods, such as wine. However, the genetic basis of qualitative and quantitative PA composition variation is still poorly understood. To elucidate the genetic architecture of the complex grape PA composition, we first carried out quantitative trait locus (QTL) analysis on a 191-individual pseudo-F1 progeny. Three categories of PA variables were assessed: total content, percentages of constitutive subunits and composite ratio variables. For nine functional candidate genes, among which eight co-located with QTLs, we performed association analyses using a diversity panel of 141 grapevine cultivars in order to identify causal SNPs.ResultsMultiple QTL analysis revealed a total of 103 and 43 QTLs, respectively for seed and skin PA variables. Loci were mainly of additive effect while some loci were primarily of dominant effect. Results also showed a large involvement of pairwise epistatic interactions in shaping PA composition. QTLs for PA variables in skin and seeds differed in number, position, involvement of epistatic interaction and allelic effect, thus revealing different genetic determinisms for grape PA composition in seeds and skin. Association results were consistent with QTL analyses in most cases: four out of nine tested candidate genes (VvLAR1, VvMYBPA2, VvCHI1, VvMYBPA1) showed at least one significant association with PA variables, especially VvLAR1 revealed as of great interest for further functional investigation. Some SNP-phenotype associations were observed only in the diversity panel.ConclusionsThis study presents the first QTL analysis on grape berry PA composition with a comparison between skin and seeds, together with an association study. Our results suggest a complex genetic control for PA traits and different genetic architectures for grape PA composition between berry skin and seeds. This work also uncovers novel genomic regions for further investigation in order to increase our knowledge of the genetic basis of PA composition.

An efficient and rapid protocol for plant nuclear DNA preparation suitable for next generation sequencing methods.

PREMISE OF THE STUDY In this study, we developed a nuclear DNA extraction protocol for Next Generation Sequencers (NGS). METHODS AND RESULTS We applied this extraction method to grapevines and coffee trees, which are known to contain many secondary metabolites. The nuclear DNA obtained was sequenced by the 454/GS-FLX method. We obtained excellent results, with less than 4% cytoplasmic DNA, in a similar way to a BAC (Bacterial Artificial Chromosome)-building protocol. We also compared our protocol with a classic DNA extraction using specific cytoplasmic DNA amplification. Results showed a lower cytoplasmic DNA contamination with the new protocol. CONCLUSIONS The method presented here is fast and economical. The DNA obtained is of high quality, with a low level of cytoplasmic DNA contamination, and very efficient for the construction of sequencing libraries.

Selection of candidate genes for grape proanthocyanidin pathway by an integrative approach

Proanthocyanidins (PA) play a major role in plant protection against biotic and abiotic stresses. Moreover these molecules are known to be beneficial for human health and are responsible for astringency of foods and beverages such as wine and thus have a great impact on the final quality of the product. Genes playing a role in the PA pathway are only partially known. The amount of available transcriptomic and genetic data to select candidate genes without a priori knowledge from orthologous function increases every day. However, the methods used so far generate so many candidate genes that it is impossible to validate all of them. In this study, we used an integrative strategy based on different screening methods to select a reduced list of candidate genes. We have crossed results from different screening methods including QTL mapping and three transcriptomic studies to select 20 candidate genes, located in QTL intervals and fulfilling at least two transcriptomic screenings. This list includes three glucosyltransferases, already suspected to have a role in the PA biosynthetic pathway. Among the 17 remaining genes, we selected three genes to perform further analysis by association genetic studies. For each of these genes, we found a polymorphism linked to PA variation. The three genes (VvMybC2-L1, VvGAT-like and VvCob-like), not previously known to play a role in PA synthesis, are promising candidates for further molecular physiology studies.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

VvGAI1 polymorphisms associate with variation for berry traits in grapevine

Gibberellins have important effects over different developmental processes in plants and the exogenous application of GA3 has been extensively used in grapevine to improve certain quality traits (to control berry size, set and weight, to decrease cluster compactness and to induce the absence of seeds in grapes). In this study, we performed association analyses using a core collection of 127 table grape accessions in order to identify causal SNPs in the gene VvGAI1, a negative regulator of gibberellins response. Seventeen quantitative descriptors, related with fertility, bunch and berry traits were used for the association analysis with VvGAI1. Fifteen polymorphisms and 15 haplotypes were identified, and some of them associated significantly with important traits for table grape: berry texture, juice yield, and bunch weight. The genetic structure found in the collection of varieties was consistent with their origins, and was included in the association test. In addition, nucleotide and haplotypic diversity, linkage disequilibrium and the possible existence of selection over the gene were also evaluated to support the association test.