De novo phased assembly of the Vitis riparia grape genome

Grapevine is one of the most important fruit species in the world. In order to better understand genetic basis of traits variation and facilitate the breeding of new genotypes, we sequenced, assembled, and annotated the genome of the American native Vitis riparia, one of the main species used worldwide for rootstock and scion breeding. A total of 164 Gb raw DNA reads were obtained from Vitis riparia resulting in a 225X depth of coverage. We generated a genome assembly of the V. riparia grape de novo using the PacBio long-reads that was phased with the 10x Genomics Chromium linked-reads. At the chromosome level, a 500 Mb genome was generated with a scaffold N50 size of 1 Mb. More than 34% of the whole genome were identified as repeat sequences, and 37,207 protein-coding genes were predicted. This genome assembly sets the stage for comparative genomic analysis of the diversification and adaptation of grapevine and will provide a solid resource for further genetic analysis and breeding of this economically important species.


Background & Summary
Since few decades and the development of sequencing technologies, the number of species whose genome has been totally sequenced has increased exponentially. There is a large variability for the quality of all the sequences assemblies. In 2017, 72 plant reference quality genome assemblies were reported in NCBI 1 . For plant breeding, the availability of a contiguous genome sequence provides a tool to better identify genes underlying traits and how they may be regulated by various environmental parameters in different genetic backgrounds. At the simplest, it allows for association of genetic markers for selection and introgression of traits across germplasm to enable the development of novel products for consumers 2,3 .
As an important crop, Vitis vinifera was one of the first higher plant species whose genome was sequenced by a French-Italian consortium 4 . The consortium decided to sequence a near homozygous V. vinifera cultivar related to Pinot Noir (PN40024) in order to facilitate the sequence assembly by limiting sequence variability. To date, this genome still stands as the reference for the grapevine community, but grapevine intra species and interspecies diversity makes using a single reference genome inadequate for studying the function of other genotypes. In order to address the variations in a cultivated V. vinifera variety, the Pinot Noir genome was sequenced using Sanger sequencing providing a high quality draft of the genome with about 10X coverage 5 . Next Generation Sequencing reads are too short to resolve abundant repeats in particular in plants genome, leading to incomplete or ambiguous assemblies 6 . Few attempts to produce high quality grapevine genomes were undertaken in grapevine and produced valuable data to study the genetic variations of V. vinifera cv. Tannat 7 and cv. Thompson seedless 8 through comparison with the reference genomes.
The last few years have seen rapid innovations in sequencing technologies and improvement in assembly algorithms that enabled the creation of highly contiguous genomes. The development of third generation sequencing technologies that deliver long reads from single molecules and carry the necessary information to phase haplotypes over several kilobases have greatly improved the feasibility of de novo assemblies [9][10][11] . Sequences of V. vinifera cv. Cabernet Sauvignon were first released 12 using PacBio sequencing and FALCON, and FALCON-Unzip pipeline 12 . This generated a 591 Mbp haplotype genome from a set of 718 primary contigs, and a set of correlated 2,037 haplotigs spanning 367 Mbp. The total p-contig size was larger than the estimated genome size of V. vinifera (~500 Mbp) suggesting that in some cases FALCON-Unzip underestimated the alternative haplotype sequences because of high heterozygosity between homologous regions, which is common in grapevine 13,14 . Later, the PacBio assembly and annotation of V. vinifera cv Chardonnay variety provided after curation of artefactual www.nature.com/scientificdata www.nature.com/scientificdata/ contig assignment, 854 p-tigs and 1883 h-tigs, totaling 490 Mb and 378 Mb 15 . More recently, another version of the Chardonnay genome was proposed with a different level of curation at 605 Mb 16 .
An evaluation of genetic diversity based on a panel of 783 V. vinifera varieties using 10 K SNPs revealed a high level of diversity (He = 0.32) and confirmed the close pedigree relationship within the cultivated grapevine due to the wide use of the most interesting parents during domestication and early selection by humans 17 . Considering that grape cultivation currently faces severe pathogen pressures and climate change, we assume that the exploitation of the natural genetic diversity may ensure the long-term sustainability of the grape and wine industries 18 . Grapes belong to the genus Vitis, which includes over 60 inter-fertile species. The most common grape cultivars derive their entire ancestry from the species V. vinifera, but wild relatives have also been exploited to create hybrid cultivars, often with increased disease resistances 19 .
To date, no wild Vitis genomes have been released so far and the only whole genome sequences for grape are from V. vinifera varieties and yet there is a clear need for genetic resources 20 . Here, we report the first de novo assembly and genome annotation of the North American native grape V. riparia. Using the latest sequencing technologies, we show that 10x Genomics Chromium data can be combined with long read PacBio sequencing to effectively determine genome phasing. The phased haplotypes of V. riparia genome will greatly contribute to give more insight into the functional consequences of genetic variants.

Methods
Sample collection, library construction and whole genome sequencing. The Vitis riparia Gloire de Montpellier (RGM) selection was obtained in 1880 by L. Vialla and R. Michel from North American collections and is the only commercially available pure V. riparia stock. RGM clone #1030 and the European native Vitis vinifera Cabernet sauvignon (CS) clone #15 were grown at INRA, Bordeaux (France). A F1 segregating population of 114 individuals named CSxRGM1995-1 was derived from the cross between CS and RGM 21 . This population was genotyped using the GBS approach 22 to create a high resolution genetic map to assist in anchoring and orienting the assembled V. riparia genome scaffolds.
Total DNA was isolated and extracted using QIAGEN Genomic-tips 100/G kit (Cat No./ID: 10243) following the tissue protocol extraction. Briefly, 1 g of young leaf material was ground in liquid nitrogen with mortar and pestle. After 3 h of lysis and one centrifugation step, the DNA was immobilized on the column. After several washing steps, DNA is eluted from the column, then desalted and concentrated by alcohol precipitation. The pellet is resuspended in TE buffer.
Three PacBio libraries with a 20-kb insert size were also constructed and sequenced on RSII platforms (97.71 Gb data; ~118-fold covering), following the standard PacBio protocol of Sequencing Kit 1.2.1 (Pacific Biosciences, USA). Four 10x Chromium Genomics libraries were constructed using the Chromium TM Genome Solution (10X Genomics, USA), and 2 × 150 bp sequenced on Illumina HiSeq3000, producing ~350 million paired-end linked-reads (~ 107-fold covering). Finally, 2 libraries for 2 × 100 bp sequencing were built with different insert sizes: 500 bp for paired-end (PE) and 6 kb for mate-pair (MP), based on the standard Illumina protocol and sequenced on the Illumina HiSeq2500. The raw reads were trimmed before being used for subsequent genome assembly. For Illumina HiSeq sequencing, the adaptor sequences, the reads containing more than 10% ambiguous nucleotides, as well as the reads containing more than 20% low-quality nucleotides (quality score less than 5), were all removed. After data cleaning and data preprocessing, we obtained a total of 164 Gb of clean data (52 Gb PacBio data, 59 Gb 10X Genomics, 33 Gb PE reads and 20 Gb MP reads,), representing 331X coverage of the V. riparia genome (Table 1).
Genome size and heterozygosity estimation. Lodhi and Reisch 23 estimated the genome size in grape to be approximately 475 Mb based on measurements using flow cytometry for 19 species including wild Vitis species, V. vinifera and V. labrusca cultivars. The measurements showed intraspecific variation in genome size between different varieties of Vitis vinifera ranging from 1C = 415 to 511 Mb, and between different North America Vitis species ranging from 1C = 411 to 541 Mb, with V. riparia around 470 Mb. Genome sequencing of different V. vinifera varieties gave values in the same range or greater depending on the methods of sequencing and assembly. In order to verify these values, we estimated genome size of V. riparia by the k-mer method 24,25 using data from pair-end and mate-pair Illumina sequencing. By analyzing the 21-mers depth distribution, a total of ~50 billion k-mers were estimated with a peak frequency of 100, corresponding to a genome size of 494 Mb and the estimated repeat sequencing ratio was 33.74%. In this study, V. riparia heterozygosity was estimated to be 0.46% (mean distance 1 SNP each 217 bp between heterozygous SNPs) from 10x Chromium Genomics data processing.
De novo Genome assembly and scaffolding of the Vitis riparia genome. We employed a hybrid de novo whole-genome assembly strategy, combining both short linked-reads and PacBio long reads data. www.nature.com/scientificdata www.nature.com/scientificdata/ Genome assembly was first performed on full PacBio cleaned reads using FALCON v0.3.0 26 . Error correction and pre-assembly were carried out with the FALCON/FALCON Unzip pipeline after evaluating the outcomes of using different parameters in FALCON during the pre-assembly process. Based on the contig N50 results, a length_cutoff of 5 kb and a length_cutoff_pr of 8 kb for the assembly step were ultimately chosen. The draft assembly was polished using Quiver 27 , which mapped the PacBio reads to the assembled genome with the BLASER pipeline 28 . Haplotypes were separated during assembly using FALCON-Unzip and the preliminary genome assembly was approximately 530 Mb (1,964 primary-contigs) and 317 Mb (3,344 haplotigs). A summary of the assembly statistics can be found in Table 1. Assembly was then processed with Purge Haplotigs 13 to investigate the proper assignment of contigs, followed by 2 rounds of polishing to correct residual SNP and INDELs errors with Pilon v1.22 software 29 using high-coverage (~106X) Illumina paired-end and mate pair data.
The 10x Chromium Genomics linked-reads were used to produce a separate V. riparia assembly using the Supernova assembler option-style = pseudohap2 and created two parallel pseudohaplotypes 30 . The mean input DNA molecule length reported by the Supernova assembler was 45 kb and the assembled genome size was 424 Mb with a N50 scaffold of 711 kb.
Subsequently, the PacBio assembly was scaffolded with the 10x Chromium Genomics one using the hybrid assembler LINKS 31 with 7 iterations, producing 870 scaffolds spanning 500 Mb with N50 = 964 kb and L50 = 255 (Table 2). Finally, genome phasing was reconstituted using Long Ranger analysis pipeline that processes Chromium sequencing output to align reads and call and phase SNPs, indels, and structural variants on the basis of molecular barcodes information.
Genotyping by Sequencing and genetic mapping. Two 96-plex GBS libraries (Keygene N.V. owns patents and patent applications protecting its Sequence Based Genotyping technologies) were constructed for the two parents (two replicates for each) and the 114 F1 plants of the cross CS × RGM. Raw reads were checked with FastQC 32 , demultiplexed with a custom script and cleaned with CutAdapt 33 . Cleaned reads were then mapped to the V. riparia RGM scaffolds previously obtained, the V. vinifera Cabernet Sauvignon contigs 12 and V. vinifera PN40024 genome assemblies 4 for SNP calling. Aligned on these genomes were performed using BWA 34 , SAMtools 35 and Picard tools 36 and SNP genotypes were detected with GATK 37 using the hardfilter parameters 38 . In the variant call format (VCF) output file only sites with less than 20% missing data and a minimum allele frequency (MAF) ≥ 0.2 were retained. The SNP set was parsed into two data sets based on a pseudo-test cross mapping strategy 39 using major_minor and get_pseudo_test_cross scripts from Hetmapps 40 . The segregation ratios of markers in the population were examined by Chi-square analysis. Markers with segregation ratios that differed from expected 1:1 at P < 0.05 were classified as segregation distortion markers and discarded. The RGM and CS sets contain 1591 and 2359 SNPs respectively. Linkage groups (LGs) were determined using software JoinMap ® 4.1 41,42 and Rqtl 43 .
LG were formed with a logarithm of odds (LOD) threshold of 6 and a maximum recombination frequency of 0.45. The 19 LGs that corresponded to the 19 chromosomes of grapevine were reconstructed and leaded to a total genetic map length of 2,268 cM and 2,514 cM for RGM and CS respectively.
Pseudo-molecule construction. The PacBio/10x Chromium Genomics hybrid scaffolding was organized into pseudo-molecules using GBS markers information from the CS × RGM genetic map. Scaffolds were anchored and oriented SNP using AllMaps 44 with the unequal weights2 parameters for a single run for the entire genome. Final pseudo-molecules were named according to Vitis vinifera PN40024 reference genome using SNP identification through SNP calling on this reference. Since PN40024 genome is the only one available who has been scaffolded into pseudo-molecules, collinearity with V. riparia was evaluated using D-GENIES 45 and showed extremely high conservation along the 19 chromosomes of the species (Fig. 1)  Genome annotation and gene prediction. Consistent with observations that long reads sequencing technologies are a better solution for resolving repeat sequences, we found that known repetitive elements accounted for 170 Mb (33.94%) of the genome in V. riparia. This is a lower proportion among grape genomes when comparing published values to date. However, when comparisons are performed with the same analysis workflow and tools 47,48 , the percentages obtained between the two genotypes were in the same range (Online-only www.nature.com/scientificdata www.nature.com/scientificdata/ Table 1). Similar to other grape genomes, long terminal repeat (LTR) elements constituted the highest proportion of all repeated elements in V. riparia, (21.44%) with Copia and Gypsy families accounting for 8.33% and 12.66% respectively. The Long Interspersed Nuclear Elements (LINEs) and Miniature Inverted-repeat Transposable Elements (MITEs) represented 3.61% and 6.02% of the whole genome respectively.
After repeat masking, the genome was ab initio annotated using MAKER-P pipeline 49,50 , SNAP 51 and Augustus 52 gene finder with 3 rounds of Maker and an Augustus prediction. Structural annotation was then followed with an Interproscan functional annotation and putative gene function assignation using BLAST on UniProtKB. MAKER-P quality metrics with a threshold of AED < 0.5 were chosen to retain the set of predicted genes. We finally generated a gene set of 37,207 protein-coding genes ( To facilitate genomic investigations for the community, a JBrowse Genome Browser 53 was set up for V. riparia pseudo-molecules and is available from https://www6.bordeaux-aquitaine.inra.fr/egfv/.

Data Records
The V. riparia genome project was deposited at NCBI under BioProject number PRJNA512170 and BioSample SAMN10662253. The DNA sequencing data from Illumina, PacBio and 10x Genomics have been deposited in the Sequence Read Archive (SRA) database under accession SRP174866 54 from SRX5189632 to SRX5189680. This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession SJAQ00000000 55 . The versions described in this paper are version SJAQ01000000 and genome ID 13150 (https:// www.ncbi.nlm.nih.gov/genome/13150). Genetic mapping data and structural and functional annotation file of the Vitis riparia assembly are available on figshare 56 .

technical Validation
To evaluate the accuracy and completeness of the V. riparia assembly, genome features were compared to those of V. vinifera (Table 2). We found that both contig and scaffold N50 lengths of Vitis riparia reached considerable continuity. The Guanine-Cytosine content (GC = 34.32%) was similar to those of V. vinifera Chardonnay (34.43%). www.nature.com/scientificdata www.nature.com/scientificdata/ To further assess the accuracy of the V. riparia genome assembly, the NGS-based short reads from whole-genome sequencing data were also aligned against the genome assembly using BWA mem 57 . We found that 98.4% of the reads were reliably aligned to the genome assembly, and 95.8% of the reads were properly aligned to the genome with their mates. Paired-end reads data were not used during the contig assembly, thus the high alignment ratio demonstrated the high quality of contig assembly.