Whole genome comparison between table and wine grapes reveals a comprehensive catalog of structural variants
BMC Plant Biology volume 14, Article number: 7 (2014)
Grapevine (Vitis vinifera L.) is the most important Mediterranean fruit crop, used to produce both wine and spirits as well as table grape and raisins. Wine and table grape cultivars represent two divergent germplasm pools with different origins and domestication history, as well as differential characteristics for berry size, cluster architecture and berry chemical profile, among others. ‘Sultanina’ plays a pivotal role in modern table grape breeding providing the main source of seedlessness. This cultivar is also one of the most planted for fresh consumption and raisins production. Given its importance, we sequenced it and implemented a novel strategy for the de novo assembly of its highly heterozygous genome.
Our approach produced a draft genome of 466 Mb, recovering 82% of the genes present in the grapevine reference genome; in addition, we identified 240 novel genes. A large number of structural variants and SNPs were identified. Among them, 45 (21 SNPs and 24 INDELs) were experimentally confirmed in ‘Sultanina’ and six SNPs in other 23 table grape varieties. Transposable elements corresponded to ca. 80% of the repetitive sequences involved in structural variants and more than 2,000 genes were affected in their structure by these variants. Some of these genes are likely involved in embryo development, suggesting that they may contribute to seedlessness, a key trait for table grapes.
This work produced the first structural variants and SNPs catalog for grapevine, constituting a novel and very powerful tool for genomic studies in this key fruit crop, particularly useful to support marker assisted breeding in table grapes.
Grapevine (Vitis vinifera L.) is one of the main fruit crops of the Mediterranean climate regions in the world ([1–3]). It is also one of the oldest agricultural crops, starting its cultivation during the Neolithic (6,000-5,000 BC) in the Transcaucasian region  from its wild progenitor V. vinifera L. subsp. sylvestris ([2, 5, 6]). It exhibits a high genetic diversity and heterozygosity level, as in most woody species. Genotypes of this species have been traditionally classified in two main groups: winemaking and table grape varieties. These last are mainly for fresh consumption but also some are intended for raisins production. Most of these cultivars are the result of centuries of selection and vegetative reproduction, representing a diverse knowledge, preferences and cultivation systems in different regions around the world . Indeed, the evaluation of the genetic diversity of collections of representative genotypes using SSR markers have revealed differences in the allelic composition between table grape and wine varieties, indicating that there is a substantial intra-specific differentiation . It is quite likely that table and wine grapes became what they are today through a divergent selection based on human preferences. On this regard, several traits such as thick pericarp, small berries with a larger number of seeds and high tannins and phenolic content have been selected for wine varieties, whereas thinner pericarps, seedlessness and larger rachis aiming to maximize the berry size, are the traits that have been selected for table grapes. There is increasing evidence that genetic diversity relies mainly in genomic structural variants such as SNPs, short sequence insertions and deletions (INDELs), inter- and intra-chromosomal translocations and inversions ([8–10]). Therefore, it is likely that differences observed between table and wine varieties are due to structural variants. Recently, a V. vinifera reference genome was assembled based on the sequencing of a nearly homozygous genotype (PN40024) . In addition, the genome sequence of the wine cultivar ‘Pinot noir’, a highly heterozygous genotype, was released . Furthermore, several genome sequencing initiatives in grapevine are in progress, most of them focused on the identification of polymorphisms related to traits of interest for wine production . However, as up to now no genomic sequence from a typical table grape variety has been released, it is not yet possible to establish at a genomic level how different are the two main groups of grapevine genotypes. This is a key aspect not only to increase the knowledge of the genome of the species but also for helping the breeding programs. Genetic variations and their associated genetic diversity are critical issues for obtaining new grape varieties. This is a labor-intensive task, where the use of marker-assisted selection (MAS) should expedite the selection process. The identification of markers to be used for MAS can be greatly improved when the structural variants present in the genome of the parents are known. Even though there is a reference genome from a wine variety, the genetic diversity observed in this species does not allow taking full advantage of this genomic tool.
‘Sultanina’ is one of the most important table grape varieties playing a pivotal role in modern breeding, mainly because of providing the seedlessness (stenospermocarpy) phenotype.
Genetic evidence indicates that Vitis vinifera is a highly heterozygous species, and the assembly of a heterozygous genome represents a bioinformatics challenge ([12, 14, 15]). In this work we sequenced and implemented a strategy for the de novo assembly of the highly heterozygous genome of ‘Sultanina’. Our results show that there are a number of structural variants with respect to the grapevine reference genome, including genome fragment translocations, INDELs and transposable elements relocalization. Moreover, a significant number of SNPs were detected and novel genes not present in the reference genome were also identified. Experimental validation of structural variants and SNPs predicted showed a high rate of success. This new assembled genome will allow us to get a better understanding of the genetics of the table grape group of cultivars, boosting its breeding based on a deeper understanding of the genomes used in the crossing blocks. We propose this assembly and its structural variants catalog as a genomics tool for this key fruit crop.
De novo assembly of ‘Sultanina’
We sequenced the diploid genome of ‘Sultanina’, a pivotal table grape genotype. The main challenge of the de novo assembly relied on its high heterozygosity ([1, 12, 14, 15]). The 25-mer analysis confirmed the highly heterozygous nature of this genome (Figure 1). To address this issue we used a novel approach called HAPLOIDIFY implemented in the ALLPATHS-LG  assembler. This is a decision process, based on statistics of the assembly, which chooses only one haplotype during the assembly. General features of the assembly are summarized in Table 1. We got a genome size of 466 Mb which is in agreement with the estimated size of the grapevine genomes ([11, 12]). The analysis of contig coverage (Additional file 1: Figure S1) showed that longest contigs are enriched in homozygosity . Indeed, 79% (326 Mb) of the total contig length was classified as homozygous and the remaining 21% (86.2 Mb) as heterozygous (see Methods for the detailed description of our classification strategy). In addition, using 95% of identity we were able to recover 82% (Table 1) of the genes present in the reference genome PN40024 and the anchoring of our scaffolds to the 19 chromosomes in this reference are well distributed, indicating again that our assembly is highly homozygous. Using 90% of identity, the recovery rate races to 86%. Thus, this de novo assembled genome offers a draft for the search of unique genomic features present in ‘Sultanina’. Interestingly, it allow us to seek for differences at the nucleotide resolution between table and wine grapes.
Novel genes found in ‘Sultanina’
From the whole genome comparison with the reference genome PN40024 we identified 240 novel genes in ‘Sultanina’ genome (Additional file 2: Table S1) that have EST support on the public NCBI EST database of Vitis sp. From them, 130 corresponded to transposon related genes and 88 to hypothetical genes. From the remaining 22 genes the most represented biological function was associated to disease resistance/defense response (13 genes). Other classes of novel genes that are represented in ‘Sultanina’ genome are related to proteolysis, embryo development, carbon-nitrogen bonds, methyltransferase and anthocyanin synthesis.
Structural variants (SVs) and SNPs catalog
We used both de novo assembly and reads mapping methods for the detection of structural variants (SVs) in the range of 1 bp to 50 kb between ‘Sultanina’ and the reference genome PN40024. We considered as SVs to INDELs, inversions and inter-intra chromosomal rearrangements; and SNPs were considered independently. We identified 310,855 insertions from 1 to 46,200 bp, 312,148 deletions from 1 to 9,993 bp and 5,871 complex SVs, defined as inversions or inter-intra chromosomal rearrangements from 10 to 41,402 bp. Also, 1,193,566 high quality SNPs were identified. Transposable elements are by far the most common genetic elements causing genomic variations in plants . In our study, we found that Gypsy-like and Copia retrotransposon elements are the most commonly found polymorphisms (Table 2), confirming previous findings in grapevine using a reduced part of the genome . We examined the whole genome distribution of SVs and SNPs. We identified homozygous and heterozygous SNPs and INDELs (Additional file 3: Table S2) which are distributed throughout the chromosomes (Figure 2). We found that around 70% of INDELs and SNPs are located in intergenic and intronic regions (Additional file 4: Table S3). Short SVs (below 50 bp) are the most abundant (Additional file 5: Figure S2). The higher frequency for those found in CDSs corresponds to SVs with lengths that are multiple of three nucleotides (Additional file 5: Figure S3) which is consistent with what has been described in other organisms . A significant number of genes exhibit homozygous INDELs, suggesting that the function of proteins encoded by these genes may be altered (Table 3). The whole genome distribution of polymorphisms revealed the existence of islands of homozygosity and heterozygosity. To further explore this phenomenon, the reference genome was divided into 4,256 disjoint intervals of length 100 kb and we counted the amount of heterozygous and homozygous variants on each interval (see Methods). We found 237 loci where both alleles were the same but diverged from the reference genome (highly homozygous variation) and 641 loci where only one allele diverged (highly heterozygous variation). The other loci could not be discriminated (Figure 2). Interestingly, among the loci that diverged between both genomes we found genes related to embryo development and it has been proposed that genomic regions with significantly high homozygosity have been related to domestication processes . About 3,700 genes showed a positive selection (based on dN/dS > 1). Among them, 540 genes had more than 10 SNPs (Additional file 6: Table S4). This suggests that around 2% of the genes present in the ‘Sultanina’ genome are undergoing a rapid divergence in protein coding regions. From these 540 genes, 410 presented a GO term associated. A GO enrichment analysis under Biological Process gave 59 categories that were statistically significantly overrepresented (Additional file 7: Table S5). Interestingly, genes related to response to stimulus, as well as anatomical and reproductive structure developments were within this group (Additional file 8: Figure S4).
Experimental confirmation of SNPs and INDELs predicted in ‘Sultanina’
Twenty seven INDELs predicted in the ‘Sultanina’ genome were selected for validation. Primer pairs amplifying fragments among 103 to 413 bp were selected and the amplicons were analyzed using capillary electrophoresis-laser-induced fluorescence (CE-LIF) assay (Additional file 9: Table S7). Twenty four INDELs, 11 deletions and 13 insertions, were confirmed (Additional file 10: Table S8). An example of these is shown in Figure 3. Interestingly, 22 out of the 24 confirmed INDELs fit the predicted homo- or heterozygous haplotype in ‘Sultanina’. A group of 23 heterozygous and homozygous SNPs predicted in the ‘Sultanina’ genome were selected to be confirmed by sequencing and qPCR-HRM (Additional file 9: Table S7). The group included 12 transitions and 11 transversions, with SNP-calling quality values distributed in the interval from 90.2 to 999. Twenty one of them (about 90%) were confirmed (Additional file 10: Table S8). Furthermore, robust and confident melting and HRM curves were optimized for six of such SNPs, and a group of 23 table grape varieties (Additional file 11: Table S9) were used to confirm the transferability of them. The average polymorphism information content value (PIC) for these six SNPs was 0.38, ranging from 0.12 to 0.5, suggesting their feasibility and transferability. As an example, the result for TSSNP820904 is shown in Figure 4. Three genotypic classes for this SNP were observed.
Seedlessness is a desirable trait in table grapes. A QTL for seedlessness has been mapped to chromosome 18 ([21–24]) and a polymorphic form of VvAGL11 (AGAMOUS-like 11) has been found to explain a high percentage of seedlessness variance in ‘Sultanina’ . Our SVs analysis confirmed a 15 bp heterozygous insertion in the 5’UTR of VvAGL11 (GSVIVT01025945001) gene in ‘Sultanina’ genome (Figure 5). This insertion is not present in this locus in the reference genome, which derives from a genotype that produces seeded fruits. In order to look for additional genes that may contribute to seedlessness in ‘Sultanina’, we searched for orthologous genes whose mutations in Arabidopsis lead to an embryo defective phenotype . Four hundred ninety six putative orthologous genes were identified in the ‘Sultanina’ genome. Forty two of these genes contained either INDELs in promoter and coding regions or no synonymous and frame shift SNPs in the coding region. Thirty of these genes were located in homozygous regions; therefore, we put more attention to these genes since they can be more tightly linked to seedlessness (Additional file 12: Table S6). Thirteen of these genes were also located in previously mapped QTLs for seedlessness in a progeny derived from the cross between ‘Ruby Seedless’ and ‘Sultanina’ . Therefore, it is likely that these genes affected by SVs or SNPs may be considered as main positional candidate genes responsible for seed development. Every SV or SNP present in each one of the 42 genes was confirmed by comparing the different reads used in the assembly of the respective contig.
‘Sultanina’ is an ancient seedless cultivar of unprecise geographical origin in old Persia. After it was brought to France and then popularized in America under the name ‘Thompson Seedless’, this cultivar has become key in the modern table grape breeding, being present in the pedigree of numerous modern varieties. It is also the main source of seedlessness used in breeding programs ([27, 28]), a prime trait for fresh consumption. Also, a number of somatic mutations exhibiting variations in berry size and seeds number have been described as derived from this genotype. However, no further studies have been done related to its phenotypic characteristics and its genetic constitution. Today, there is an increasing effort to establish the relationship between phenotypes and the genomic information of a species. In the case of the grapevine, the availability of a reference genome based on a wine-derived genotype (line PN40024) has not been as effective for table grape genetic studies as it would have been expected. This is probably due to the genetic divergence between wine and table grapes , phenotypically represented by traits such as the presence of seeds and their relationship with berry size , or the different content in phenolic compounds such as flavanols, flavonols and hydroxy-benzoic acids . In this work, we obtained the first draft of the highly heterozygous ‘Sultanina’ genome based entirely on NGS technologies and de novo assembly. The assembly of highly heterozygous genomes exhibits unique and difficult challenges. Moreover, there are few algorithmic ideas able to handle this kind of complexity ([12, 14, 15]). In plants, the most frequent strategy to build reference genomes has been based on the selection of highly homozygous individuals, what in most woody species is a very long process and seldomly addressed, not available for table grapes. Here we used ALLPATHS-LG assembler to tackle the heterozygotic nature of ‘Sultanina’. Our strategy led to a draft genome sharing similar metrics (size of the genome, number of contigs and scaffolds, as well as gene content) with the previously assembled genome of the heterozygous ‘Pinot noir’ , which was obtained through Sanger and 454 sequencing technologies. After a whole genome comparison between the ‘Sultanina’ genome and the grapevine reference genome PN40024, we succeeded to provide the first comprehensive catalog of SVs and SNPs between both genotypes, at the nucleotide level. This catalog contains about 1,800,000 variants including SNPs, INDELs, translocations and inversions. The SNP rate is in agreement with previous reports on this species . The classification of variants into homozygous and heterozygous revealed enriched islands of each kind distributed throughout the chromosomes. The experimental confirmation proved that about 90% of our SVs and SNPs predictions were true, showing the precision of the catalog. Indeed, our experimental validation of SVs can be considered as the first evidence suggesting the feasibility and transferability of SNP reported in ‘Sultanina’ catalog as useful tools for genetic studies in table grapes. We also found a set of rapidly evolving genes (540 genes with dN/dS ratios larger than one and 10 or more SNPs each) and 240 novel genes. Interestingly, GO terms related to pathogen resistance and quality traits were over-represented in rapidly evolving genes. This is likely due to a combination of natural selection by pressure of pathogens and artificial pressure due to the domestication process with selection of agronomically important traits (quality trait genes such as those related to cell wall metabolism and anatomical and reproductive structure development categories). A similar phenomenon has been observed in species such as rice, sorghum and maize when genomes of different varieties or landraces are compared ([10, 30–32]).
SVs and SNPs are a source of genetic variability; since, they are important in generating new genes or allelic variants that may be selected by natural or artificial means, if they confer an advantage to the fruit crop. The search for genes responsible for traits of interest has been tackled by seeking QTLs. However, the reduced size of the mapping populations commonly used in woody fruit crops renders too wide confidence intervals, corresponding to genomic regions of various cM harboring tens to hundreds of candidate genes per QTL . The availability of the ‘Sultanina’ genome would help to improve the saturation of the genomic region where a QTL has been identified, in a simpler and better way than it has been done until now based on the reference genome. This should reduce substantially the list of candidate genes to focus in subsequent analyses. In addition, the availability of a catalog of structural variants and SNPs can help in the identification of candidates genes related to traits of interest. In this work, we confirmed the INDEL previously described in the regulatory region of the VvAGL11 gene. This gene has been proposed as the main responsible for seedlessness , and this INDEL has been converted into an effective selection marker for seedlessness . Interestingly, the sequencing of ‘Sultanina’ highlighted other SVs and SNPs affecting the structure of genes related to embryo development, some of them located in other QTLs that explain the residual seedlessness phenotypic variance .
This ‘Sultanina’ draft genome should improve the efficiency of molecular assisted breeding in table grape and the search for genes associated to different traits could be better approached. In addition, the proposed SVs and SNPs catalog will become a powerful tool to improve and expedite processes such as synteny-based comparisons, mutations detection, transgenes localization, among other genetic studies and breeding-related applications in table grapes.
We produced a draft of the ‘Sultanina’ genome of size 466 Mb. Eighty-two percent of the genes present in the reference genome were recovered and 240 novel genes were identified. A large number of SVs and SNPs were found. Forty-five (21 SNPs and 24 INDELs) were experimentally confirmed in ‘Sultanina’ and among them six SNPs in other 23 table grape varieties. Two thousand genes were affected by these variants. The ‘Sultanina’ genome should improve the efficiency of molecular assisted breeding in table grape and the search for genes associated to different traits could be better approached. In addition, the proposed SVs and SNPs catalog will become a powerful tool to improve and expedite processes such as synteny-based comparisons, mutations detection, transgenes localization, among other genetic studies and breeding-related applications in table grapes.
The homozygous grapevine reference genome PN40024, mRNA and protein sequences were downloaded from the GENOSCOPE database . The heterozygous grape assembly, mRNA and protein sequences were downloaded from the IASMA database . Repeats libraries were downloaded from RepBase .
The sequenced vine was originally collected from a vineyard located in the vicinity of Santiago, Chile. It was confirmed as a true-to-type ‘Sultanina’ by using a standard set of microsatellite markers . It was planted in a pot in 2011 and has been maintained at INIA La Platina Experimental Station since then. The vine is clean of the most common grapevine viruses as tested by standard RT-PCR. A Total of 1,572 million reads were generated using Illumina sequencing. Three libraries were sequenced at different insert sizes (180 bp, 600 bp and 2000-3000 bp) using the Genome Analyzer II and HISeq 2000 platforms (Macrogen Inc. Seoul, Korea). The total sequencing represents a raw coverage of 327X, using an estimated genome size of 480 Mb for a highly heterozygous grape genome .
Before genome assembly, Illumina reads were corrected using Quake  with the following parameters: minimun length of reads 70 bp and minimum quality 20; 20% of the reads were thus eliminated. The genome assembly was performed by ALLPATHS-LG assembler  with a raw total coverage of 200X for overlap (180 bp) and jumping (600 bp, 2000-3000 bp) libraries. Since the genome is highly heterozygous , the HAPLOIDIFY variable was set. This setting examine mismatches in the graph of the assembly that result from single nucleotide variations (even those that are very close), selects one branch and discards the other following statistical criteria. Then, it replaces the reads from the discarded branches with the chosen ones, haploidifying the data set. A total of 47,863,057 reads were changed using this strategy. Then the assembly proceeded as described in . At the end of the assembly single nucleotide variations were reintroduced (by mapping back reads) and a mix of both haplotypes was obtained.
Identification of structural variants (SVs)
For the detection of short SVs (<50 bp) we aligned to the reference genome PN40024 the pair-end reads using BWA  with default parameters. Then we called the INDELs using the Dindel  program. To detect long (>50 bp) SVs (insertions, deletions and invertions) we applied a process similar to the one described previously in the literature using assembly methods . The assembled scaffold was pre-aligned to the reference genome using Nucmer  with mum option enabled. It counts matches that are unique in both the reference and the query. The matches were filtered with delta-filter allowing only one to one alignments, a minimum identity of 94% and a minimum alignment length of 1,000 bp. The scaffolds and best aligned regions were extracted and aligned using LASTZ  with ambiguous ‘N’ treatment, gap free extension tolerance up to 50 kb and high scoring segment pairs chaining options enabled. Scaffolds with no match in the pre-alignment were aligned to the whole reference genome with the same options. Finally, SV break points were extracted using all aligned regions between the assembly and the reference genome. To predict inter and intra chromosomal re-arrangements we used BWA alignment and BreakDancer  program with -t, -d and -g options enabled, allowing read tracking for each candidate SV.
In silico validation of SVs
Short SVs were validated using reads supported by Dindel program . Dindel uses a Bayesian approach to call short indels and genotypes by realigning reads to the candidate haplotype, avoiding homopolymer errors. Also, Dindel is optimized for Illumina sequencing technology. In order to validate long SVs, we implemented an approach similar to SoapSV . Our pipeline input is a modified version of the SoapSV output file that was produced after the alignment of the scaffolds against the PN40024 genome. This file contains all break points (coordinates) for each SV in our assembly and in the PN40024 genome. We splitted this process into four steps. First we removed all the SVs overlapping gap regions. Secondly, we divided the output file into two sets, insertions and deletions. Thirdly, we validated the SVs. We computed the coverage continuity in 500 bp up and down flanking regions of the SVs and inside the SVs using SAMtools  depth command. We considered valid a deletion in ‘Sultanina’ genome if the coverage dropped below a half or less in the reference genome and maintained constant ratio in the assembly. We considered valid an insertion, if a region contained half or less depth coverage in the reference genome and coverage maintained constant ratio in flanking regions of the break points of SVs in the assembly. For inter and intra chromosomal re-arrangements predicted by BreakDancer , we mapped the reads supporting the re-arrangements to our whole genome assembly. When at least three pair-end clones were mapped to the expected insert size, the inter or intra chromosomal re-arrangements were validated. The inversion predictions based on whole genome alignments were validated when they overlapped with an inversion prediction called by BreakDancer supported by at least three clones. Finally, the INDELs effect was predicted using SNPeff .
For high-quality SNPs, we excluded reads that were repeated (those that had more than one position in the genome) according to Bowtie  results. We initially called the SNPs using the mpileup function of SAMtools  with default parameters. Then, the candidate SNPs were filtered by VCFtools  using a window of 10 bp, a minimum depth of eight and a minimum quality of 40. Finally, the SNP effect was predicted by SNPeff  program.
SNPs and short indels were classified into homozygous or heterozygous by probabilistic methods implemented in SAMtools  and Dindel  programs. To define whether long and complex SVs where homozygous or heterozygous, we first classified the assembled contig into homozygous or heterozygous using the contig coverage . In order to do that, we took a total of 100X of reads and aligned them to the assembled contigs by BWA . By using intervals of different lengths, we could classify the homozygous and heterozygous contigs (Additional file 1: Figure S1). Contigs having coverage over 50X were considered homozygous, whereas those with coverage below 50X were considered heterozygous. Thus, the SVs genotype (homozygous or heterozygous) was defined based on the location of the variant within a given contig. To explore the island phenomenon (Figure 2), we performed a total of 4,256 Fisher exact tests with p < 0.01, corrected with FDR and fold change of ±2, using the rate between the amounts of homozygous against heterozygous variants in each window. Using these tests, we were able to infer the window genotype.
A total of 2,581 scaffolds of ‘Sultanina’ (total size of 6.5 Mb) could not be aligned to the reference genome and were used as input to search for putative novel genes. Identification of putative novel genes was performed using AUGUSTUS  with complete gene option enabled. A total of 1,113 candidate genes were found. Using MEGABLAST  we mapped the predicted mRNAs to the public EST Vitis sp. database downloaded from NCBI using as filters minimal-score equal to 100 and with a minimal-identity of 90%. It produced 327 genes with evidence in grape ESTs. Then, we eliminated all of those genes having a MEGABLAST match, using the same parameters, with the transcripts of the reference genome PN40024 . This process yielded 240 novel genes.
The functional annotation of the novel genes was done using the Non-Redundant database (pvalue 1e - 10) and Interpro .
Mapping the reference mRNAs of PN40024 reference genome into the ‘Sultanina’ genome
Using GMAP  we placed the public reference mRNAs from the reference PN40024 using parameter min-coverage 70% and min-identity 95%. We were able to place 82% of the reference genes in our genome assembly. With a less strict parameter of 90% of identity we mapped 86% of the reference transcripts in the ‘Sultanina’ assembly.
Repeat elements within SVs
In order to infer the most polymorphic elements causing variations in grape, we masked the reference and the assembled ‘Sultanina’ genomes using the last version of grape repeats from RepBase . We counted the amount of each kind of elements present within the long SVs (SVs >50 bps). If the deletion was in the reference, we counted the repeat elements present in the assembled ‘Sultanina’ genome. If an insertion appeared in the reference, we counted repeated elements in the reference (Additional file 13: Figure S5). The process was similar when the deletion or insertion was in the assembly.
The AgriGO tool  was used to detect enriched gene ontology terms through Singular Enrichment Analysis (SEA) coupled with available background data of the Arabidopsis TAIR 10 genome project . GO Term association was done by taking the grapevine genes and searching by blastp, the best homologue present in Arabidopsis . P-values for enrichment terms were calculated using hypergeometric distribution and statistical testing method with the Jekutieli multi-test adjustment method.
Total DNA was extracted from young leaves as has been described . The purified DNA was dissolved in TE buffer 1X and RNA was removed by incubating the sample with DNAase-free RNAase A. DNA concentration was measured using Qubit 2.0 digital fluorometer quantitation (Life Technologies). Samples with concentrations above 40 ng/uL were considered for the experiments.
A group of 27 INDELs identified in the ‘Sultanina’ genome were selected, including homozygous and heterozygous SVs. A total of 32 primers were designed, 27 of which were used to amplify them considering their performance (Additional file 9: Table S7). For capillary electrophoresis, PCR reactions were performed in a total volume of 10 uL, including the specific primers (0.1 uM) for each INDEL, approximately 40 ng of template DNA, 5 uL of GoTaq® Green Master Mix 2X (Promega) and 2 uL of ultrapure water (Applichem). PCR was performed at 95°C for 4 min, followed by 35 cycles of 95°C for 30 s, 57°C for 40 s, 72°C for 40 s; and a final cycle of 72°C for 5 min in a Thermo Electron’s Px2 Thermal Cycler (Thermo Electron Corp.).
INDELs confirmation by capillary electrophoresis-laser-induced fluorescence (CE-LIF) assay
An aliquot of 2 uL of PCR product was mixed with 22 uL of dsDNA Reagent Kit 35-500 bp buffer of Advanced Analytical, following conditions recommended, on Fragment Analyzer™ Automated CE System, using a 12-Capillary array cartridge (50 um [ID], 55 cm [EFF], 80 cm [TOT]), from Advanced Analytical. A pre-run was performed using 8.0 kV for 30 s, sample injection of 7.5 kV for 10 s, and a separation of 8.0 kV for 80 min. The analysis was conducted using the PROSize software and the results obtained were manually examinated.
Real time PCR and qPCR-HRM
A group of 23 SNPs predicted in the ‘Sultanina’ genome were selected, including homozygous and heterozygous SNPs. Specific primers were designed and used to amplify each SNP (Additional file 9: Table S7). Real time PCR reactions contained 5 uL of EvaGreen® Master Mix Dye 2X, 0.2 uM each primer and 0.5 ng of template DNA in a total reaction volume of 10 uL. The reactions were performed on a 72-Well Rotorgene-Q (Qiagen). Cycling conditions were 95°C for 2 min, and 50 cycles of 95°C for 5 s, 58°C for 10 s, and 72°C for 5 s. Following steps were 72°C for 2 min, 95°C for 5 s, 50°C for 30 s. The annealing temperature was optimized for each primer. Selected annealing temperature for primer TSSNP1037434 was 60°C and for TSSNP820904 and TSSNP820907 was 62°C; for all the other primers was 58°C. HRM was carried out from 65°C to 90°C, with 0.1°C increments each 2 s. Hold pre-melting at 65°C for 30 s and a final step at 65°C for 5 min were used. Raw HRM curves were recorded and normalized using the Rotorgene Q Series Software 2.0.2. HRM curve for each individual was visually scored. The data from low quality amplification were removed from HRM analysis. In particular, runs with CT value over 30 were considered not suitable for the analysis. Genotype assignations were done manually by examining normalized and derivatizes melt plots. Also, qPCR-HRM amplicons were quantified using Qubit 2.0 digital fluorometer quantitation (Life Technologies), and samples with concentrations above 20 ng/uL were sequenced. Alignments between reference genome sequence and SNPs amplified fragments were made using Sequencher software, in order to confirm these SNPs.
Single nucleotide polymorphism(s)
- SSR marker:
Simple sequence repeat marker
Insertion(s) and deletion(s)
Coding DNA sequence
Capillary electrophoresis-laser-induced fluorescence
Quantitative PCR high-resolution melting
Polymorphism information content value
Quantitative trait loci.
Myles S, Chia J-M, Hurwitz B, Simon C, Zhong GY, Buckler E, Ware D: Rapid genomic characterization of the genus Vitis. PLoS ONE. 2010, 5: e8219-10.1371/journal.pone.0008219.
Myles S, Boyko AR, Owens CL, Brown PJ, Grassi F, Aradhya MK, Prins B, Reynolds A, Chia J-M, Ware D, et al: Genetic structure and domestication history of the grape. Proc Natl Acad Sci USA. 2011, 108: 3530-3535. 10.1073/pnas.1009363108.
Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M: Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-seq. Plant Physiol. 2010, 152: 1787-1795. 10.1104/pp.109.149716.
McGovern PE: Ancient wine: the search for the origins of viticulture. Princeton, NJ: Princeton University Press; 2003.
Arroyo-García R, Ruiz-García L, Bolling L, Ocete R, López MA, Arnold C, Ergul A, Soylemezoglu G, Uzun HI, Cabello F, et al: Multiple origins of cultivated grapevine (Vitis vinifera L. ssp. sativa) based on chloroplast DNA polymorphisms. Mol Ecol. 2006, 15: 3707-3714. 10.1111/j.1365-294X.2006.03049.x.
Pelsy F: Molecular and cellular mechanisms of diversity within grapevine varieties. Heredity. 2010, 104: 331-340. 10.1038/hdy.2009.161.
Aradhya MK, Dangl GS, Prins BH, Boursiquot JM, Walker MA, Meredith CP, Simon CJ: Genetic structure and differentiation in cultivated grape Vitis vinifera L. Genet Res. 2003, 81: 179-192. 10.1017/S0016672303006177.
Huang X, Han B: A crop of maize variants. Nature Genetics. 2012, 44: 734-735. 10.1038/ng.2326.
Huang X, Kurata N, Wei X, Wang ZX, Wang A, Zhao Q, Zhao Y, Liu K, Lu H, Li W, et al: A map of rice genome variation reveals the origin of cultivated rice. Nature. 2012, 490: 497-501. 10.1038/nature11532.
Zheng L-Y, Guo X-S, He B, Sun L-J, Peng Y, Dong S-S, Liu T-F, Jiang S, Ramachandran S, Liu C-M, et al: Genome-wide patterns of genetic variation in sweet and grain sorghum (Sorghum bicolor). Genome Bio. 2011, 12: R114-10.1186/gb-2011-12-11-r114.
Jaillon O, Aury J-M, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.
Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro A, Pruss D, Pindo M, Fitzgerald LM, Vezzulli S, Reid J, et al: A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PLoS ONE. 2007, 2: e1326-10.1371/journal.pone.0001326.
International variant data. http://www.vitaceae.org.
Price JC, Udall JA, Bodily PM, Ward JA, Schatz MC, Page JT, Jensen JD, Snell QO, Clement MJ: De novo identification of "heterotigs" towards accurate and in-phase assembly of complex plant genomes. 2012, Las Vegas: Proceedings of BIOCOMP12.
Velasco R, Zharkikh A, Affourtit J, Dhingra A, Cestaro A, Kalyanaraman A, Fontana P, Bhatnagar SK, Troggio M, Pruss D, et al: The genome of the domesticated apple (Malus domestica Borkh.). Nat Genet. 2010, 42: 833-839. 10.1038/ng.654.
Gnerre S, Maccallum I, Przybylski D, Ribeiro FJ, Burton JN, Walker BJ, Sharpe T, Hall G, Shea TP, Sykes S, et al: High-quality draft assemblies of mammalian genomes from massively parallel sequence data. Proc Natl Acad Sci USA. 2011, 108: 1513-1508. 10.1073/pnas.1017351108.
Lisch D: How important are transposons for plant evolution?. Nat Rev Genet. 2012, 14: 49-61.
Carrier G, Cunff LL, Dereeper A, Legrand D, Sabot F, Bouchez O, Audeguin L, Boursiquot J-M, This P, Schonbach C, et al: Transposable elements are a major cause of somatic polymorphism in Vitis vinifera L. PLoS ONE. 2012, 7: e32973-10.1371/journal.pone.0032973.
Li Y, Zheng H, Luo R, Wu H, Zhu H, Li R, Cao H, Wu B, Huang S, Shao H, et al: Structural variation in two human genomes mapped at single-nucleotide resolution by whole genome de novo assembly. Nat Biotechnol. 2011, 29: 723-730. 10.1038/nbt.1904.
Doebley JF, Gaut BS, Smith BD: The molecular genetics of crop domestication. Cell. 2006, 127: 1309-1321. 10.1016/j.cell.2006.12.006.
Cabezas JA, Cervera MT, Ruiz-García L, Carreño J, Martínez-Zapater JM: A genetic analysis of seed and berry weight in grapevine. Genome. 2006, 49: 1572-1585. 10.1139/g06-122.
Doligez A, Bouquet A, Danglot Y, Lahogue F, Riaz S, Meredith P, Edwards J, This P: Genetic mapping of grapevine (Vitis vinifera L.) applied to the detection of QTLs for seedlessness and berry weight. Theor Appl Genet. 2002, 105: 780-795. 10.1007/s00122-002-0951-z.
Fanizza G, Lamaj F, Costantini L, Chaabane R, Grando MS: QTL analysis for fruit yield components in table grapes (Vitis vinifera). Theor Appl Genet. 2005, 111: 658-664. 10.1007/s00122-005-2016-6.
Mejía N, Soto B, Guerrero M, Casanueva X, Houel C, Miccono MA, Ramos R, Cunff LL, Boursiquot J-M, Hinrichsen P, et al: Molecular, genetic and transcriptional evidence for a role of VvAGL11 in stenospermocarpic seedlessness in grapevine. BMC Plant Biol. 2011, 11: 57-10.1186/1471-2229-11-57.
SeedGennes Project. http://www.seedgenes.org.
Mejía N, Gebauer M, Muñoz L, Hewstone N, Muñoz C, Hinrichsen P: Identification of QTLs for seedlessness, berry size, and ripening date in a seedless x seedless table grape progeny. Am J Enol Vitic. 2007, 58: 499-507.
Ibáñez J, Vargas AM, Palancar M, Borrego J, de Andrés MT: Genetic relationships among table-grape varieties. Am J Enol Vitic. 2009, 60: 35-42.
Vargas AM, Teresa de Andrés M, Borrego J, Ibáñez J: Pedigrees of fifty table-grape cultivars. Am J Enol Vitic. 2009, 60: 525-532.
Liang Z, Owens CL, Zhong G-Y, Cheng L: Polyphenolic profiles detected in the ripe berries of Vitis vinifera germplasm. Food Chem. 2011, 129: 940-950. 10.1016/j.foodchem.2011.05.050.
Hurwitz BL, Kudrna D, Yu Y, Sebastian A, Zuccolo A, Jackson SA, Ware D, Wing RA, Stein L: Rice structural variation: a comparative analysis of structural variation between rice and three of its closest relatives in the genus Oryza. Plant J. 2010, 63: 990-1003. 10.1111/j.1365-313X.2010.04293.x.
Li S, Wang S, Deng Q, Zheng A, Zhu J, Liu H, Wang L, Gao F, Zou T, Huang B, et al: Identification of genome-wide variations among three elite restorer lines for hybrid-rice. PLoS ONE. 2012, 7: e30952-10.1371/journal.pone.0030952.
Springer NM, Ying K, Fu Y, Ji T, Yeh C-T, Jia Y, Wu W, Richmond T, Kitzman J, Rosenbaum H, et al: Maize inbreds exhibit high levels of copy number variation (CNV) and presence/absence variation (PAV) in genome content. PLoS Genet. 2009, 5: e1000734-10.1371/journal.pgen.1000734.
Duchêne E, Butterlin G, Dumas V, Merdinoglu D: Towards the adaptation of grapevine varieties to climate change: QTLs and candidate genes for developmental stages. Theor Appl Genet. 2012, 124: 623-635. 10.1007/s00122-011-1734-1.
Karaagac E, Vargas AM, TeresaDe Andres M, Carreno I, Ibáñez J, Carreño J, Martínez-Zapater JM, Cabezas JA: Marker assisted selection for seedlessness in table grape breeding. Tree Genet Genomes. 2012, 8: 1003-1015. 10.1007/s11295-012-0480-0.
Genoscope database. http://www.genoscope.cns.fr
IASMA database. http://genomics.research.iasma.it
This P, Jung A, Boccacci P, Borrego J, Botta R, Costantini L, Crespan M, Dangl G, Eisenheld C, Ferreira-Monteiro F, et al: Development of a standard set of microsatellite reference alleles for identification of grape cultivars. Theor Appl Genet. 2004, 109: 1448-1458. 10.1007/s00122-004-1760-3.
Kelley DR, Schatz MC, Salzberg SL: Quake: quality-aware detection and correction of sequencing errors. Genome Bio. 2010, 11: R116-10.1186/gb-2010-11-11-r116.
Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.
Albers CA, Lunter G, MacArthur DG, McVean G, Ouwehand WH, Durbin R: Dindel: accurate indel calls from short-read data. Genome Res. 2011, 21: 961-973. 10.1101/gr.112326.110.
Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Bio. 2004, 5: R12-10.1186/gb-2004-5-2-r12.
Harris R: Improved pairwise alignment of genomic DNA. PhD thesisPennsylvania State University; 2007.
Chen K, Wallis JW, McLellan MD, Larson DE, Kalicki JM, Pohl CS, McGrath SD, Wendl MC, Zhang Q, Locke DP, et al: BreakDancer: an algorithm for high-resolution mapping of genomic structural variation. Nat Methods. 2009, 6: 677-681. 10.1038/nmeth.1363.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Subgroup GPDP, et al: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.
Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM: A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012, 6: 80-92. 10.4161/fly.19695.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Bio. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, et al: The variant call format and VCFtools. Bioinformatics. 2011, 27: 2156-2158. 10.1093/bioinformatics/btr330.
Stanke M, Diekhans M, Baertsch R, Haussler D: Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008, 24: 637-644. 10.1093/bioinformatics/btn013.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J. Mol. Biol. 1990, 3: 403-410.
Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, López R: InterProScan: protein domains identifier. Nucleic Acids Res. 2005, 33: W116-W120. 10.1093/nar/gki442.
Wu TD, Watanabe CK: GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005, 21: 1859-1875. 10.1093/bioinformatics/bti310.
AgriGO tool. http://bioinfo.cau.edu.cn/agriGO.
Arabidopsis TAIR 10 genome. http://arabidopsis.org/index.jsp.
Lodhi MA, Ye GN, Weeden NF, Reisch BI: A simple and efficient method for DNA extraction from grapevine cultivars and Vitis species. Plant Mol Biol Rep. 1994, 12: 6-13. 10.1007/BF02668658.
Botstein D, White RL, Skolnick M, Davis RW: Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet. 1980, 32: 314-331.
Chen H, He H, Zou Y, Chen W, Yu R, Liu X, Yang Y, Gao YM, Xu JL, Fan LM, Li Y, Li ZK, Deng XW: Development and application of a set of breeder-friendly SNP markers for genetic analyses and molecular breeding of rice (Oryza sativa L.). Theor Appl Genet. 2011, 123: 869-879. 10.1007/s00122-011-1633-5.
This project was supported by grants: Fondap 1509007, Basal programs PFB-03 and PFB-16, Genoma-Chile Fondef G07I-1002, CIRIC-INRIA Chile (line Natural Ressources) and Millennium Nucleus ICM-P10-062-F. We acknowledge the National Laboratory for High Performance Computing at the Center for Mathematical Modeling (PIA ECM-02- CONICYT).
The assembled genome and all of the associated variant analyses are freely available at http://vitisdb.cmm.uchile.cl/publicationmaterial. Reads can be downloaded from NCBI using STUDY accession SRP026420.
The authors have no conflicts of interest.
ADG, PH, AM and AO conceived the study. PH, AM and AO supervised the project. ADG, AMA, CME, PV, DT and CM participated in data analysis. CME designed and performed the experimental validation. ADG, AMA, PH, AM, CME and AO wrote the manuscript. All authors were involved in discussion of the manuscript. All authors read and approved the final manuscript.
Patricio Hinrichsen, Ariel Orellana and Alejandro Maass contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1: Histograms of contig coverage at 100X. The contig coverage is defined as the average depth at each position in a given contig. We depict histograms for different ranges of contig length (CL = Contig Length). Contigs with an average coverage out of the interval [20,120] are excluded. The largest contigs are mostly homozygous while smaller contigs are predominantly heterozygous. (PDF 54 KB)
Additional file 5: Figure S2: Length distribution of structural variants (SVs). Frequency of homozygous and heterozygous SVs in ‘Sultanina’ genome according to their length, and Figure S3 Structural variants in CDS. Frequency of homozygous and heterozygous SVs in coding sequences of ‘Sultanina’ genome according to their length. (PDF 164 KB)
Additional file 6: Table S4: SNPs in ‘Sultanina’ genes that present dN/dS ratio (nonsynonymous-to-synonymous substitutions) higher than 1 and their respective best homologue in Arabidopsis thaliana. (PDF 90 KB)
Additional file 8: Figure S4: GO enrichment on rapidly evolving genes under Biological Process Category. Only categories with FDR < 0.05 were considered as over represented. The analysis was done using the online Agrigo tool and the GO Slim plant category. The boxes contain the GO number, the category description, the p-value between parenthesis, the number of genes in each category out of the 410 that presented a GO term associated, the number of genes in each category out of 28,352 Arabidopsis genes. The arrows indicate the relationship among the GO categories. Black solid arrows mean that a GO category is also included in the other one, red solid arrows mean that one GO category positively regulates the other, green solid arrows mean that the GO category negatively regulates the other, black dashed arrows indicate that there are two significant nodes related to the GO category, black dotted arrows indicate that only one significant node is related to the GO category. (PDF 405 KB)
Additional file 9: Table S7: Primer sequences for the 50 selected structural variations (INDELs and SNPs) experimentally confirmed. (PDF 55 KB)
Additional file 10: Table S8: Experimental validation of 45 SVs (24 INDELs and 21 SNPs) identified in the ‘Sultanina’ genome. (PDF 39 KB)
Additional file 11: Table S9: Selected group of table grape varieties plus one used for wine production (‘Tokay’), representing different genetic backgrounds. (PDF 36 KB)
Additional file 12: Table S6: ‘Sultanina’ orthologous genes of Arabidopsis thaliana embryo development related genes containing SVs in promoter and coding regions. (PDF 53 KB)
Additional file 13: Figure S5: Identification of transposable elements within INDELs. We masked the repeat elements in the reference and the ‘Sultanina’ genomes using RepBase. Then, for each INDEL of length over than 50 bp we counted the total size in bp of the repeated elements contained within it. (PDF 43 KB)
About this article
Cite this article
Di Genova, A., Almeida, A.M., Muñoz-Espinoza, C. et al. Whole genome comparison between table and wine grapes reveals a comprehensive catalog of structural variants. BMC Plant Biol 14, 7 (2014). https://doi.org/10.1186/1471-2229-14-7
- ‘Sultanina’ genome
- Vitis vinifera L
- Structural variants