Haplotyping the Vitis collinear core genome with rhAmpSeq improves marker transferability in a diverse genus

Wine metabolome and sensory analyses demonstrate the oenological potential of novel grapevine genotypes for sustainable viticulture in warm climates

BACKGROUND

Genetic breeding is essential to develop grapevine genotypes adapted to warm climates and resistant to pathogens. Traditionally cultivated Vitis vinifera is susceptible to biotic and abiotic stresses. Winemakers and consumers, however, perceive wines from non-vinifera or hybrid cultivars as inferior. In this study, sensory analyses and comprehensive metabolite profiling by targeted and untargeted approaches were used to investigate the oenological potential of wines from grapes of genotypes developed throughout four breeding cycles to improve climate adaptation, sugar contents and berry color.

RESULTS

Novel genotypes had higher yields and the wines exhibited increased contents of polyphenols, including anthocyanins. Volatile monoterpenes in the wines decreased throughout breeding cycles in the absence of selective pressure. Polyphenol contents were higher in intermediate wines, with hydroxytyrosol contents reaching up to three times reported values. Mouthfeel attributes astringency, leafy taste, flavor and body, and persistency showed significant correlation with untargeted features. Supervised model-based analyses of the metabolome effectively discriminate wines from distinct genetic origins.

CONCLUSION

Taken together, the results demonstrate the potential of novel grapevine genotypes to a more sustainable viticulture and quality wine production in warm climates. Comprehensive metabolite profiling of the wines reveals that genotype clustering is dependent on the chemical class and that traits not submitted to selective pressure are also altered by breeding. Supervised multivariate models were effective to predict the genetic origin of the wines based on the metabolic profile, indicating the potential of the technique to identify biomarkers for wines from sustainable genotypes. © 2024 Society of Chemical Industry.

A Sustainable Strategy for Marker-Assisted Selection (MAS) Applied in Grapevine (Vitis spp.) Breeding for Resistance to Downy (Plasmopara Viticola) and Powdery (Erysiphe Necator) Mildews

Plant breeders utilize marker-assisted selection (MAS) to identify favorable or unfavorable alleles in seedlings early. In this task, they need methods that provide maximum information with minimal input of time and economic resources. Grape breeding aimed at producing cultivars resistant to pathogens employs several resistance loci (Rpv, Ren, and Run) that are ideal for implementing MAS. In this work, a sustainable MAS protocol was developed based on non-purified DNA (crude), multiplex PCR of SSR markers, and capillary electrophoresis, and its application on grapevine seedlings to follow some main resistance loci was described. The optimized protocol was utilized on 8440 samples and showed high efficiency, reasonable throughput (2-3.2 min sample), easy handling, flexibility, and tolerable costs (reduced by at least 3.5 times compared to a standard protocol). The Rpv, Ren, and Run allelic data analysis did not show limitations to loci combination and pyramiding, but segregation distortions were frequent and displayed both low (undesired) and high rates of inheritance. The protocol and results presented are useful tools for grape breeders and beyond, and they can address sustainable changes in MAS. Several progenies generated have valuable pyramided resistance and will be the subject of new studies and implementation in the breeding program.

The wild side of grape genomics

With broad genetic diversity and as a source of key agronomic traits, wild grape species (Vitis spp.) are crucial to enhance viticulture's climatic resilience and sustainability. This review discusses how recent breakthroughs in the genome assembly and analysis of wild grape species have led to discoveries on grape evolution, from wild species' adaptation to environmental stress to grape domestication. We detail how diploid chromosome-scale genomes from wild Vitis spp. have enabled the identification of candidate disease-resistance and flower sex determination genes and the creation of the first Vitis graph-based pangenome. Finally, we explore how wild grape genomics can impact grape research and viticulture, including aspects such as data sharing, the development of functional genomics tools, and the acceleration of genetic improvement.

Biography of Vitis genomics: recent advances and prospective

The grape genome is the basis for grape studies and breeding, and is also important for grape industries. In the last two decades, more than 44 grape genomes have been sequenced. Based on these genomes, researchers have made substantial progress in understanding the mechanism of biotic and abiotic resistance, berry quality formation, and breeding strategies. In addition, this work has provided essential data for future pangenome analyses. Apart from de novo assembled genomes, more than six whole-genome sequencing projects have provided datasets comprising almost 5000 accessions. Based on these datasets, researchers have explored the domestication and origins of the grape and clarified the gene flow that occurred during its dispersed history. Moreover, genome-wide association studies and other methods have been used to identify more than 900 genes related to resistance, quality, and developmental phases of grape. These findings have benefited grape studies and provide some basis for smart genomic selection breeding. Moreover, the grape genome has played a great role in grape studies and the grape industry, and the importance of genomics will increase sharply in the future.

Metabolite profiling reveals the influence of grapevine genetic distance on the chemical signature of juices

BACKGROUND

Yield, disease tolerance, and climate adaptation are important traits in grapevine genetic breeding programs. Selection for these characteristics causes unpredictable changes in primary and specialized metabolism, affecting the physicochemical properties and chemical composition of the berries and their processed products, juice, and wine. In this study, we investigated the influence of the genetic distance between grapevine genotypes on the chemical signatures of the juices, by integrating comprehensive metabolic profiling to genetic analyses.

RESULTS

The studied grapevine cultivars exhibited low genetic diversity. Breeding for agronomic traits promoted higher contents of soluble sugars, total phenolics, and anthocyanins in the juices. Untargeted juice metabolomics identified a total of 147 metabolites, consisting of 30 volatiles, 21 phenolics, and 96 ultra-high-performance liquid chromatography-mass spectrometry (UHPLC-MS) features. Juices from grapes of the most recent cultivars exhibited increased levels of trans-resveratrol, catechin, and luteolin. The blend of volatiles from juices of later cultivars was also more complex, consisting of 29 distinct metabolites in 'BRS Magna'. Grapes from 'BRS Carmem', an intermediate cultivar, gave the most divergent UHPLC-MS juice profile.

CONCLUSION

Contents of soluble solids, total phenolics, and anthocyanins in grape juices were increased by controlled crosses and hybrid selection. Integrative analyses demonstrated that the juices' metabolic profiles accurately represent the cultivars' genetic distances. Juices from 'BRS Violeta' and 'BRS Magna' show relevant positive association with health-related phenolics and a distinct set of odor volatiles, although these characteristics were specifically sought by breeding. © 2023 Society of Chemical Industry.

A super-pangenome of the North American wild grape species

BACKGROUND

Capturing the genetic diversity of wild relatives is crucial for improving crops because wild species are valuable sources of agronomic traits that are essential to enhance the sustainability and adaptability of domesticated cultivars. Genetic diversity across a genus can be captured in super-pangenomes, which provide a framework for interpreting genomic variations.

RESULTS

Here we report the sequencing, assembly, and annotation of nine wild North American grape genomes, which are phased and scaffolded at chromosome scale. We generate a reference-unbiased super-pangenome using pairwise whole-genome alignment methods, revealing the extent of the genomic diversity among wild grape species from sequence to gene level. The pangenome graph captures genomic variation between haplotypes within a species and across the different species, and it accurately assesses the similarity of hybrids to their parents. The species selected to build the pangenome are a great representation of the genus, as illustrated by capturing known allelic variants in the sex-determining region and for Pierce's disease resistance loci. Using pangenome-wide association analysis, we demonstrate the utility of the super-pangenome by effectively mapping short reads from genus-wide samples and identifying loci associated with salt tolerance in natural populations of grapes.

CONCLUSIONS

This study highlights how a reference-unbiased super-pangenome can reveal the genetic basis of adaptive traits from wild relatives and accelerate crop breeding research.

A multitiered haplotype strategy to enhance phased assembly and fine mapping of a disease resistance locus

Fine mapping of quantitative trait loci (QTL) to dissect the genetic basis of traits of interest is essential to modern breeding practice. Here, we employed a multitiered haplotypic marker system to increase fine mapping accuracy by constructing a chromosome-level, haplotype-resolved parental genome, accurate detection of recombination sites, and allele-specific characterization of the transcriptome. In the first tier of this system, we applied the preexisting panel of 2,000 rhAmpSeq core genome markers that is transferable across the entire Vitis genus and provides a genomic resolution of 200 kb to 1 Mb. The second tier consisted of high-density haplotypic markers generated from Illumina skim sequencing data for samples enriched for relevant recombinations, increasing the potential resolution to hundreds of base pairs. We used this approach to dissect a novel Resistance to Plasmopara viticola-33 (RPV33) locus conferring resistance to grapevine downy mildew, narrowing the candidate region to only 0.46 Mb. In the third tier, we used allele-specific RNA-seq analysis to identify a cluster of 3 putative disease resistance RPP13-like protein 2 genes located tandemly in a nonsyntenic insertion as candidates for the disease resistance trait. In addition, combining the rhAmpSeq core genome haplotype markers and skim sequencing-derived high-density haplotype markers enabled chromosomal-level scaffolding and phasing of the grape Vitis × doaniana 'PI 588149' assembly, initially built solely from Pacific Biosciences (PacBio) high-fidelity (HiFi) reads, leading to the correction of 16 large-scale phasing errors. Our mapping strategy integrates high-density, phased genetic information with individual reference genomes to pinpoint the genetic basis of QTLs and will likely be widely adopted in highly heterozygous species.

Basic leucine zipper gene VvbZIP61 is expressed at a quantitative trait locus for high monoterpene content in grape berries

The widely appreciated muscat flavor of grapes and wine is mainly attributable to the monoterpenes that accumulate in ripe grape berries. To identify quantitative trait loci (QTL) for grape berry monoterpene content, an F1 mapping population was constructed by a cross between two grapevine genotypes, one with neutral aroma berries (cv. 'Beifeng') and the other with a pronounced muscat aroma (elite Vitis vinifera line '3-34'). A high-density genetic linkage map spanning 1563.7 cM was constructed using 3332 SNP markers that were assigned to 19 linkage groups. Monoterpenes were extracted from the berry of the F1 progeny, then identified and quantified by gas chromatography-mass spectrometry. Twelve stable QTLs associated with the amounts of 11 monoterpenes in berries were thus identified. In parallel, the levels of RNA in berries from 34 diverse cultivars were estimated by RNA sequencing and compared to the monoterpene content of the berries. The expression of five genes mapping to stable QTLs correlated well with the monoterpene content of berries. These genes, including the basic leucine zipper VvbZIP61 gene on chromosome 12, are therefore considered as potentially being involved in monoterpene metabolism. Overexpression of VvbZIP61 in Vitis amurensis callus through Agrobacterium-mediated transformation significantly increased the accumulation of several monoterpenes in the callus, including nerol, linalool, geranial, geraniol, β-myrcene, and D-limonene. It is hypothesized that VvbZIP61 expression acts to increase muscat flavor in grapes. These results advance our understanding of the genetic control of monoterpene biosynthesis in grapes and provide important information for the marker-assisted selection of aroma compounds in grape breeding.

Vitis labrusca genome assembly reveals diversification between wild and cultivated grapevine genomes

Wild grapevines are important genetic resources in breeding programs to confer adaptive fitness traits and unique fruit characteristics, but the genetics underlying these traits, and their evolutionary origins, are largely unknown. To determine the factors that contributed to grapevine genome diversification, we performed comprehensive intragenomic and intergenomic analyses with three cultivated European (including the PN40024 reference genome) and two wild North American grapevine genomes, including our newly released Vitis labrusca genome. We found the heterozygosity of the cultivated grapevine genomes was twice as high as the wild grapevine genomes studied. Approximately 30% of V. labrusca and 48% of V. vinifera Chardonnay genes were heterozygous or hemizygous and a considerable number of collinear genes between Chardonnay and V. labrusca had different gene zygosity. Our study revealed evidence that supports gene gain-loss events in parental genomes resulted in the inheritance of hemizygous genes in the Chardonnay genome. Thousands of segmental duplications supplied source material for genome-specific genes, further driving diversification of the genomes studied. We found an enrichment of recently duplicated, adaptive genes in similar functional pathways, but differential retention of environment-specific adaptive genes within each genome. For example, large expansions of NLR genes were discovered in the two wild grapevine genomes studied. Our findings support variation in transposable elements contributed to unique traits in grapevines. Our work revealed gene zygosity, segmental duplications, gene gain-and-loss variations, and transposable element polymorphisms can be key driving forces for grapevine genome diversification.

Multigenic resistance to Xylella fastidiosa in wild grapes (Vitis sps.) and its implications within a changing climate

Xylella fastidiosa is a bacterium that infects crops like grapevines, coffee, almonds, citrus and olives. There is little understanding of the genes that contribute to plant resistance, the genomic architecture of resistance, and the potential role of climate in shaping resistance, in part because major crops like grapevines (Vitis vinifera) are not resistant to the bacterium. Here we study a wild grapevine species, V. arizonica, that segregates for resistance. Using genome-wide association, we identify candidate resistance genes. Resistance-associated kmers are shared with a sister species of V. arizonica but not with more distant species, suggesting that resistance evolved more than once. Finally, resistance is climate dependent, because individuals from low ( < 10 °C) temperature locations in the wettest quarter were typically susceptible to infection, likely reflecting a lack of pathogen pressure in colder climates. In fact, climate is as effective a predictor of resistance phenotypes as some genetic markers. We extend our climate observations to additional crops, predicting that increased pathogen pressure is more likely for grapevines and almonds than some other susceptible crops.

Discovery and genome-guided mapping of REN12 from Vitis amurensis, conferring strong, rapid resistance to grapevine powdery mildew

Powdery mildew resistance genes restrict infection attempts at different stages of pathogenesis. Here, a strong and rapid powdery mildew resistance phenotype was discovered from Vitis amurensis 'PI 588631' that rapidly stopped over 97% of Erysiphe necator conidia, before or immediately after emergence of a secondary hypha from appressoria. This resistance was effective across multiple years of vineyard evaluation on leaves, stems, rachises, and fruit and against a diverse array of E. necator laboratory isolates. Using core genome rhAmpSeq markers, resistance mapped to a single dominant locus (here named REN12) on chromosome 13 near 22.8-27.0 Mb, irrespective of tissue type, explaining up to 86.9% of the phenotypic variation observed on leaves. Shotgun sequencing of recombinant vines using skim-seq technology enabled the locus to be further resolved to a 780 kb region, from 25.15 to 25.93 Mb. RNASeq analysis indicated the allele-specific expression of four resistance genes (NLRs) from the resistant parent. REN12 is one of the strongest powdery mildew resistance loci in grapevine yet documented, and the rhAmpSeq sequences presented here can be directly used for marker-assisted selection or converted to other genotyping platforms. While no virulent isolates were identified among the genetically diverse isolates and wild populations of E. necator tested here, NLR loci like REN12 are often race-specific. Thus, stacking of multiple resistance genes and minimal use of fungicides should enhance the durability of resistance and could enable a 90% reduction in fungicides in low-rainfall climates where few other pathogens attack the foliage or fruit.

Genetic analysis of grapevine root system architecture and loci associated gene networks

Own-rooted grapevines and grapevine rootstocks are vegetatively propagated from cuttings and have an adventitious root system. Unraveling the genetic underpinnings of the adventitious root system architecture (RSA) is important for improving own-rooted and grafted grapevine sustainability for a changing climate. Grapevine RSA genetic analysis was conducted in an Vitis sp. 'VRS-F2' population. Nine root morphology, three total root system morphology, and two biomass traits that contribute to root anchorage and water and nutrient uptake were phenotyped. Quantitative trait loci (QTL) analysis was performed using a high density integrated GBS and rhAmpSeq genetic map. Thirty-one QTL were detected for eleven of the RSA traits (surface area, root volume, total root length, fresh weight, number of tips, forks or links, longest root and average root diameter, link length, and link surface area) revealing many small effects. Several QTL were colocated on chromosomes 1, 9, 13, 18, and 19. QTL with identical peak positions on chromosomes 1 or 13 were enriched for AP2-EREBP, AS2, C2C2-CO, HMG, and MYB transcription factors, and QTL on chromosomes 9 or 13 were enriched for the ALFIN-LIKE transcription factor and regulation of autophagy pathways. QTL modeling for individual root traits identified eight models explaining 13.2 to 31.8% of the phenotypic variation. 'Seyval blanc' was the grandparent contributing to the allele models that included a greater surface area, total root length, and branching (number of forks and links) traits promoting a greater root density. In contrast, V. riparia 'Manitoba 37' contributed the allele for greater average branch length (link length) and diameter, promoting a less dense elongated root system with thicker roots. LATERAL ORGAN BOUNDARY DOMAIN (LBD or AS2/LOB) and the PROTODERMAL FACTOR (PFD2 and ANL2) were identified as important candidate genes in the enriched pathways underlying the hotspots for grapevine adventitious RSA. The combined QTL hotspot and trait modeling identified transcription factors, cell cycle and circadian rhythm genes with a known role in root cell and epidermal layer differentiation, lateral root development and cortex thickness. These genes are candidates for tailoring grapevine root system texture, density and length in breeding programs.

Characterization of fatty acid compositions in longissimus thoracis muscle and identification of candidate gene and SNPs related to polyunsaturated fatty acid in Hu sheep

Fatty acid (FA) composition contributes greatly to the quality and nutritional value of lamb meat. In the present study, FA was measured in longissimus thoracis (LT) muscles of 1,085 Hu sheep using gas chromatography. Comparative transcriptomic analysis was conducted in LT muscles to identify differentially expressed genes (DEGs) between six individuals with high polyunsaturated fatty acids (H-PUFA, 15.27% ± 0.42%) and six with low PUFA (L-PUFA, 5.22% ± 0.25%). Subsequently, the single nucleotide polymorphisms (SNPs) in a candidate gene PLIN2 were correlated with FA traits. The results showed a total of 29 FA compositions and 8 FA groups were identified, with the highest content of monounsaturated fatty acids (MUFA, 46.54%, mainly C18:1n9c), followed by saturated fatty acids (SFA, 44.32%, mainly C16:0), and PUFA (8.72%, mainly C18:2n6c), and significant correlations were observed among the most of FA traits. Transcriptomic analyses identified 110 upregulated and 302 downregulated DEGs between H-PUFA and L-PUFA groups. The functional enrichment analysis revealed three significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and 17 gene ontology (GO) terms, in which regulation of lipolysis in adipocytes, the AMPK signaling pathway, and the PPAR signaling pathway may play important roles in FA metabolism and biosynthesis. In addition, weighted gene co-expression network analysis (WGCNA) identified 37 module genes associated with PUFA-related traits. In general, PLIN1, LIPE, FABP4, LEP, ACACA, ADIPOQ, SCD, PCK2, FASN, PLIN2, LPL, FABP3, THRSP, and ACADVL may have a great impact on PUFA metabolism and lipid deposition. Four SNPs within PLIN2 were significantly associated with FA. Of those, SNP1 (g.287 G>A) was significantly associated with C18:1n9c and MUFA, and SNP4 (g.7807 T>C) was significantly correlated with PUFA (C18:3n3). In addition, the combined genotype of SNP1 (g.287 G>A), SNP3 (g.7664 T>C), and SNP4 (g.7807 T>C) were significantly correlated with C16:1, C17:0, C18:1C6, PUFA (C18:3n3, C22:6n3), and n-6/n-3 PUFA. These results contribute to the knowledge of the biological mechanisms and genetic markers involved in the composition of FA in Hu sheep.

A high-density integrated map for grapevine based on three mapping populations genotyped by the Vitis18K SNP chip

We present a high-density integrated map for grapevine, allowing refinement and improved understanding of the grapevine genome, while demonstrating the applicability of the Vitis18K SNP chip for linkage mapping. The improvement of grapevine through biotechnology requires identification of the molecular bases of target traits by studying marker-trait associations. The Vitis18K SNP chip provides a useful genotyping tool for genome-wide marker analysis. Most linkage maps are based on single mapping populations, but an integrated map can increase marker density and show order conservation. Here we present an integrated map based on three mapping populations. The parents consist of the well-known wine cultivars 'Cabernet Sauvignon', 'Corvina' and 'Rhine Riesling', the lesser-known wine variety 'Deckrot', and a table grape selection, G1-7720. Three high-density population maps with an average inter-locus gap ranging from 0.74 to 0.99 cM were developed. These maps show high correlations (0.9965-0.9971) with the reference assembly, containing only 93 markers with large order discrepancies compared to expected physical positions, of which a third is consistent across multiple populations. Moreover, the genetic data aid the further refinement of the grapevine genome assembly, by anchoring 104 yet unanchored scaffolds. From these population maps, an integrated map was constructed which includes 6697 molecular markers and reduces the inter-locus gap distance to 0.60 cM, resulting in the densest integrated map for grapevine thus far. A small number of discrepancies, mainly of short distance, involve 88 markers that remain conflictual across maps. The integrated map shows similar collinearity to the reference assembly (0.9974) as the single maps. This high-density map increases our understanding of the grapevine genome and provides a useful tool for its further characterization and the dissection of complex traits.

HiFi chromosome-scale diploid assemblies of the grape rootstocks 110R, Kober 5BB, and 101-14 Mgt

Cultivated grapevines are commonly grafted on closely related species to cope with specific biotic and abiotic stress conditions. The three North American Vitis species V. riparia, V. rupestris, and V. berlandieri, are the main species used for breeding grape rootstocks. Here, we report the diploid chromosome-scale assembly of three widely used rootstocks derived from these species: Richter 110 (110R), Kober 5BB, and 101–14 Millardet et de Grasset (Mgt). Draft genomes of the three hybrids were assembled using PacBio HiFi sequences at an average coverage of 53.1 X-fold. Using the tool suite HaploSync, we reconstructed the two sets of nineteen chromosome-scale pseudomolecules for each genome with an average haploid genome size of 494.5 Mbp. Residual haplotype switches were resolved using shared-haplotype information. These three reference genomes represent a valuable resource for studying the genetic basis of grape adaption to biotic and abiotic stresses, and designing trait-associated markers for rootstock breeding programs.

Measurement(s)	Genome Assembly Sequence
Technology Type(s)	PacBio Sequel System
Sample Characteristic - Organism	Vitis cinerea var. helleri x Vitis rupestris • Vitis riparia x Vitis rupestris • Vitis cinerea var. helleri x Vitis riparia
Sample Characteristic - Location	State of California

High throughput saliency-based quantification of grape powdery mildew at the microscopic level for disease resistance breeding

Imaging-based high throughput phenotyping (HTP) systems have demonstrated promising solutions to enhance genetic understanding of grapevine powdery mildew (PM) resistance and have accelerated PM-resistant cultivar breeding. The accuracy and throughput of extracting phenotypic traits from images are still the bottleneck of modern HTP systems, especially at the microscopic level. The goal of this study was to develop a saliency-based processing pipeline for the quantification of PM infection in microscopic images and comprehensively evaluate its performance for genetic analyses. An input image was segregated into subimages that were classified as infected or healthy by a pretrained CNN classifier. Saliency maps from the classification were generated post-hoc and used for the quantification of PM infection in the input image at the pixel level without the use of mask annotations. A total of seven phenotypic traits were extracted from images collected for a biparental population. Experimental results showed that optimal combinations of convolutional neural network and saliency methods achieved strong measurement correlations (r = 0.74 to 0.75) with human assessments at the image patch level, and the traits calculated by the saliency-based processing pipeline were highly correlated (r = 0.87 to 0.88) with reference PM infection ratings at the leaf image level. The high quantification accuracy of the saliency-based pipeline led to the increased explanation of phenotypic variance and reliable identification of quantitative trait loci. Therefore, the saliency-based processing pipeline can be used as an effective and efficient analysis tool for PM disease research and breeding programs in the future, especially agricultural and life science studies requiring microscopic image analysis.

A multi-locus genome-wide association study reveals the genetics underlying muscadine antioxidant in berry skin

Muscadine berries display enhanced nutraceutical value due to the accumulation of distinctive phytochemical constituents with great potential antioxidant activity. Such nutritional and health merits are not only restricted to muscadine, but muscadine berries accumulate higher amounts of bioactive polyphenolics compared with other grape species. For the genetic study of the antioxidant trait in muscadine, a multi-locus genome-wide association study (GWAS) with 350 muscadine genotypes and 1,283 RNase H2 enzyme-dependent amplicon sequencing (rhAmpSeq) markers was performed. Phenotyping was conducted with several antioxidant-related traits, including total phenolic content (TPC), total flavonoid content (TFC), 1,1-diphenyl-2-picrylhydrazyl (DPPH) free radical scavenging activity, and FRAP antioxidant assay in muscadine berry skin. The correlation coefficient analysis revealed that the TPC, and DPPH/FRAP activities were significantly correlated. Through the GWAS analysis, 12 QTNs were identified from the four traits, of which six were pleiotropic QTNs. Two pleiotropic QTNs, chr2_14464718 and chr4_16491374, were commonly identified from the TPC and DPPH/FRAP activities. Co-located genes with the two pleiotropic QTNs were isolated, and two candidate genes were identified with transcriptome analysis. UDP-glycosyltransferase and 4-hydroxy-4-methyl-2-oxoglutarate aldolase were the candidate genes that are positively and negatively correlated to the quantitative property of traits, respectively. These results are the first genetic evidence of the quantitative property of antioxidants in muscadine and provide genetic resources for breeding antioxidant-rich cultivars for both Muscadinia and Euvitis species.

Assembly of complete diploid-phased chromosomes from draft genome sequences

De novo genome assembly is essential for genomic research. High-quality genomes assembled into phased pseudomolecules are challenging to produce and often contain assembly errors because of repeats, heterozygosity, or the chosen assembly strategy. Although algorithms that produce partially phased assemblies exist, haploid draft assemblies that may lack biological information remain favored because they are easier to generate and use. We developed HaploSync, a suite of tools that produces fully phased, chromosome-scale diploid genome assemblies, and performs extensive quality control to limit assembly artifacts. HaploSync scaffolds sequences from a draft diploid assembly into phased pseudomolecules guided by a genetic map and/or the genome of a closely related species. HaploSync generates a report that visualizes the relationships between current and legacy sequences, for both haplotypes, and displays their gene and marker content. This quality control helps the user identify misassemblies and guides Haplosync’s correction of scaffolding errors. Finally, HaploSync fills assembly gaps with unplaced sequences and resolves collapsed homozygous regions. In a series of plant, fungal, and animal kingdom case studies, we demonstrate that HaploSync efficiently increases the assembly contiguity of phased chromosomes, improves completeness by filling gaps, corrects scaffolding, and correctly phases highly heterozygous, complex regions.

Berry Anthocyanin, Acid, and Volatile Trait Analyses in a Grapevine-Interspecific F2 Population Using an Integrated GBS and rhAmpSeq Genetic Map

Increased map density and transferability of markers are essential for the genetic analysis of fruit quality and stress tolerance in interspecific grapevine populations. We used 1449 GBS and 2000 rhAmpSeq markers to develop a dense map for an interspecific F₂ population (VRS-F₂) that was derived by selfing a single F₁ from a Vitis riparia x ‘Seyval blanc’ cross. The resultant map contained 2519 markers spanning 1131.3 cM and was highly collinear with the Vitis vinifera ‘PN40024’ genome. Quantitative trait loci (QTL) for berry skin color and flower type were used to validate the map. Four rhAmpSeq transferable markers were identified that can be used in pairs (one pistillate and one hermaphroditic) to predict pistillate and hermaphrodite flower type with ≥99.7% accuracy. Total and individual anthocyanin diglucoside QTL mapped to chromosome 9 near a 5-O-GLUCOSYLTRANSFERASE candidate gene. Malic acid QTL were observed on chromosome 1 and 6 with two MALATE DEHYRDROGENASE CYTOPLASMIC 1 and ALUMINUM-ACTIVATED MALATE TRANSPORTER 2-LIKE (ALMT) candidate genes, respectively. Modeling malic acid identified a potential QTL on chromosome 8 with peak position in proximity of another ALMT. A first-ever reported QTL for the grassy smelling volatile (E)-2-hexenal was found on chromosome 2 with a PHOSPHOLIPID HYDROPEROXIDE GLUTATHIONE PEROXIDASE candidate gene near peak markers.

Stable QTL for malate levels in ripe fruit and their transferability across Vitis species

Malate is a major contributor to the sourness of grape berries (Vitis spp.) and their products, such as wine. Excessive malate at maturity, commonly observed in wild Vitis grapes, is detrimental to grape and wine quality and complicates the introgression of valuable disease resistance and cold hardy genes through breeding. This study investigated an interspecific Vitis family that exhibited strong and stable variation in malate at ripeness for five years and tested the separate contribution of accumulation, degradation, and dilution to malate concentration in ripe fruit in the last year of study. Genotyping was performed using transferable rhAmpSeq haplotype markers, based on the Vitis collinear core genome. Three significant QTL for ripe fruit malate on chromosomes 1, 7, and 17, accounted for over two-fold and 6.9 g/L differences, and explained 40.6% of the phenotypic variation. QTL on chromosomes 7 and 17 were stable in all and in three out of five years, respectively. Variation in pre-veraison malate was the major contributor to variation in ripe fruit malate (39%), and based on two and five years of data, respectively, their associated QTL overlapped on chromosome 7, indicating a common genetic basis. However, use of transferable markers on a closely related Vitis family did not yield a common QTL across families. This suggests that diverse physiological mechanisms regulate the levels of this key metabolite in the Vitis genus, a conclusion supported by a review of over a dozen publications from the past decade, showing malate-associated genetic loci on all 19 chromosomes.

Candidate resistance genes to foliar phylloxera identified at Rdv3 of hybrid grape

The foliage of the native grape species Vitis riparia and certain cold-hardy hybrid grapes are particularly susceptible to the insect pest phylloxera, Daktulosphaira vitifoliae Fitch. A previous study using a cold-hardy hybrid grape biparental F1 population (N~125) detected the first quantitative trait locus (QTL) for foliar resistance on chromosome 14, designated as resistance to Daktulosphaira vitifoliae 3 (Rdv3). This locus spans a ~7-Mbp (10-20 cM) region and is too wide for effective marker-assisted selection or identification of candidate genes. Therefore, we fine mapped the QTL using a larger F1 population, GE1783 (N~1023), and genome-wide rhAmpSeq haplotype markers. Through three selective phenotyping experiments replicated in the greenhouse, we screened 184 potential recombinants of GE1783 using a 0 to 7 severity rating scale among other phylloxera severity traits. A 500-kb fine mapped region at 4.8 Mbp on chromosome 14 was identified. The tightly linked rhAmpSeq marker 14_4805213 and flanking markers can be used for future marker-assisted breeding. This region contains 36 candidate genes with predicted functions in disease resistance (R genes and Bonzai genes) and gall formation (bifunctional 3-dehydroquinate dehydratase/shikimate dehydrogenase). Disease resistance genes suggest a traditional R-gene-mediated resistance mechanism often accompanied by a hypersensitive response, which has been widely studied in the plant pathology field. A novel resistance mechanism, non-responsiveness to phylloxera gall formation is proposed as a function of the bifunctional dehydratase gene, which plays a role in gallic acid biosynthesis and is important in gall formation. This study has implications for improvement of foliar phylloxera resistance in cold-hardy hybrid germplasm and is a starting place to understand the mechanism of resistance in crops to gall-forming insects.

Chromosome-level genome sequence assembly and genome-wide association study of Muscadinia rotundifolia reveal the genetics of 12 berry-related traits

Vitis has two subgenera: Euvitis, which includes commercially important Vitis vinifera and interspecific hybrid cultivars, and Muscadinia. Of note, the market for Muscadinia grapes remains small, and only Muscadinia rotundifolia is cultivated as a commercial crop. To establish a basis for the study of Muscadinia species, we generated chromosome-level whole-genome sequences of Muscadinia rotundifolia cv. Noble. A total of 393.8 Mb of sequences were assembled from 20 haploid chromosomes, and 26 394 coding genes were identified from the sequences. Comparative analysis with the genome sequence of V. vinifera revealed a smaller size of the M. rotundifolia genome but highly conserved gene synteny. A genome-wide association study of 12 Muscadinia berry-related traits was performed among 356 individuals from breeding populations of M. rotundifolia. For the transferability of markers between Euvitis and Muscadinia, we used 2000 core genome rhAmpSeq markers developed to allow marker transferability across Euvitis species. A total of 1599 (80%) rhAmpSeq markers returned data in Muscadinia. From the GWAS analyses, we identified a total of 52 quantitative trait nucleotides (QTNs) associated with the 12 berry-related traits. The transferable markers enabled the direct comparison of the QTNs with previously reported results. The whole-genome sequences along with the GWAS results provide a new basis for the extensive study of Muscadinia species.

The Genetic Basis of Anthocyanin Acylation in North American Grapes (Vitis spp.)

Hydroxycinnamylated anthocyanins (or simply 'acylated anthocyanins') increase color stability in grape products, such as wine. Several genes that are relevant for anthocyanin acylation in grapes have been previously described; however, control of the degree of acylation in grapes is complicated by the lack of genetic markers quantitatively associated with this trait. To characterize the genetic basis of anthocyanin acylation in grapevine, we analyzed the acylation ratio in two closely related biparental families, Vitis rupestris B38 × 'Horizon' and 'Horizon' × Illinois 547-1, for 2 and 3 years, respectively. The acylation ratio followed a bimodal and skewed distribution in both families, with repeatability estimates larger than 0.84. Quantitative trait locus (QTL) mapping with amplicon-based markers (rhAmpSeq) identified a strong QTL from 'Horizon' on chromosome 3, near 15.85 Mb in both families and across years, explaining up to 85.2% of the phenotypic variance. Multiple candidate genes were identified in the 14.85-17.95 Mb interval, in particular, three copies of a gene encoding an acetyl-CoA-benzylalcohol acetyltransferase-like protein within the two most strongly associated markers. Additional population-specific QTLs were found in chromosomes 9, 10, 15, and 16; however, no candidate genes were described. The rhAmpSeq markers reported here, which were previously shown to be highly transferable among the Vitis genus, could be immediately implemented in current grapevine breeding efforts to control the degree of anthocyanin acylation and improve the quality of grapes and their products.

Genomic mating in outbred species: predicting cross usefulness with additive and total genetic covariance matrices

Diverse crops are both outbred and clonally propagated. Breeders typically use truncation selection of parents and invest significant time, land, and money evaluating the progeny of crosses to find exceptional genotypes. We developed and tested genomic mate selection criteria suitable for organisms of arbitrary homozygosity level where the full-sibling progeny are of direct interest as future parents and/or cultivars. We extended cross variance and covariance variance prediction to include dominance effects and predicted the multivariate selection index genetic variance of crosses based on haplotypes of proposed parents, marker effects, and recombination frequencies. We combined the predicted mean and variance into usefulness criteria for parent and variety development. We present an empirical study of cassava (Manihot esculenta), a staple tropical root crop. We assessed the potential to predict the multivariate genetic distribution (means, variances, and trait covariances) of 462 cassava families in terms of additive and total value using cross-validation. Most variance (89%) and covariance (70%) prediction accuracy estimates were greater than zero. The usefulness of crosses was accurately predicted with good correspondence between the predicted and the actual mean performance of family members breeders selected for advancement as new parents and candidate varieties. We also used a directional dominance model to quantify significant inbreeding depression for most traits. We predicted 47,083 possible crosses of 306 parents and contrasted them to those previously tested to show how mate selection can reveal the new potential within the germplasm. We enable breeders to consider the potential of crosses to produce future parents (progeny with top breeding values) and varieties (progeny with top own performance).

Introgression among North American wild grapes (Vitis) fuels biotic and abiotic adaptation

BACKGROUND

Introgressive hybridization can reassort genetic variants into beneficial combinations, permitting adaptation to new ecological niches. To evaluate evolutionary patterns and dynamics that contribute to introgression, we investigate six wild Vitis species that are native to the Southwestern United States and useful for breeding grapevine (V. vinifera) rootstocks.

RESULTS

By creating a reference genome assembly from one wild species, V. arizonica, and by resequencing 130 accessions, we focus on identifying putatively introgressed regions (pIRs) between species. We find six species pairs with signals of introgression between them, comprising up to ~ 8% of the extant genome for some pairs. The pIRs tend to be gene poor, located in regions of high recombination and enriched for genes implicated in disease resistance functions. To assess potential pIR function, we explore SNP associations to bioclimatic variables and to bacterial levels after infection with the causative agent of Pierce's disease (Xylella fastidiosa). pIRs are enriched for SNPs associated with both climate and bacterial levels, suggesting that introgression is driven by adaptation to biotic and abiotic stressors.

CONCLUSIONS

Altogether, this study yields insights into the genomic extent of introgression, potential pressures that shape adaptive introgression, and the evolutionary history of economically important wild relatives of a critical crop.

Diploid chromosome-scale assembly of the Muscadinia rotundifolia genome supports chromosome fusion and disease resistance gene expansion during Vitis and Muscadinia divergence

Muscadinia rotundifolia, the muscadine grape, has been cultivated for centuries in the southeastern United States. M. rotundifolia is resistant to many of the pathogens that detrimentally affect Vitis vinifera, the grape species commonly used for winemaking. For this reason, M. rotundifolia is a valuable genetic resource for breeding. Single-molecule real-time reads were combined with optical maps to reconstruct the two haplotypes of each of the 20 M. rotundifolia cv. Trayshed chromosomes. The completeness and accuracy of the assembly were confirmed using a high-density linkage map. Protein-coding genes were annotated using an integrated and comprehensive approach. This included using full-length cDNA sequencing (Iso-Seq) to improve gene structure and hypothetical spliced variant predictions. Our data strongly support that Muscadinia chromosomes 7 and 20 are fused in Vitis and pinpoint the location of the fusion in Cabernet Sauvignon and PN40024 chromosome 7. Disease-related gene numbers in Trayshed and Cabernet Sauvignon were similar, but their clustering locations were different. A dramatic expansion of the Toll/Interleukin-1 Receptor-like Nucleotide-Binding Site Leucine-Rich Repeat (TIR-NBS-LRR) class was detected on Trayshed chromosome 12 at the Resistance to Uncinula necator 1 (RUN1)/Resistance to Plasmopara viticola 1 (RPV1) locus, which confers strong dominant resistance to powdery and downy mildews. A genome browser, annotation, and Blast tool for Trayshed are available at www.grapegenomics.com.

Multiple independent recombinations led to hermaphroditism in grapevine

We studied the grape sex-determining region (SDR) in 12 Vitis genomes and demonstrated its conservation across 556 genotypes including 193 accessions from 47 world-wide wild grapevine species and 363 accessions of cultivated grapevine. Although the grape SDR is recombination free in all wild species, we found two distinct hermaphrodite (H) haplotypes (H1 and H2) among the cultivated grapevines, both chimeras of male (M) and female (f) haplotypes. The two independent recombinations carry different genetic signatures which long predate the domestication of grapevine, suggesting independent evolutions of this trait in wild European grapevine gene pools prior to human domestication.

Hermaphroditic (perfect) flowers were a key trait in grapevine domestication, enabling a drastic increase in yields due to the efficiency of self-pollination in the domesticated grapevine (Vitis vinifera L. ssp. vinifera). In contrast, all extant wild Vitis species are dioecious, each plant having only male or female flowers. In this study, we identified the male (M) and female (f) haplotypes of the sex-determining region (SDR) in the wild grapevine species V. cinerea and confirmed the boundaries of the SDR. We also demonstrated that the SDR and its boundaries are precisely conserved across the Vitis genus using shotgun resequencing data of 556 wild and domesticated accessions from North America, East Asia, and Europe. A high linkage disequilibrium was found at the SDR in all wild grape species, while different recombination signatures were observed along the hermaphrodite (H) haplotype of 363 cultivated accessions, revealing two distinct H haplotypes, named H1 and H2. To further examine the H2 haplotype, we sequenced the genome of two grapevine cultivars, 'Riesling' and 'Chardonnay'. By reconstructing the first two H2 haplotypes, we estimated the divergence time between H1 and H2 haplotypes at ∼6 million years ago, which predates the domestication of grapevine (∼8,000 y ago). Our findings emphasize the important role of recombination suppression in maintaining dioecy in wild grape species and lend additional support to the hypothesis that at least two independent recombination events led to the reversion to hermaphroditism in grapevine.

Genome-Wide Identification of Sequence Variations and SSR Marker Development in the Munake Grape Cultivar

The Munake grape cultivar produces uniquely flavored and high-quality fruits. Despite the numerous beneficial agronomic traits of Munake, there are few genetic resources available for this cultivar. To address this knowledge gap, the entire genome was sequenced using whole-genome sequencing approaches and compared with a Vitis vinifera L. reference genome. This study describes more than 3 million single nucleotide polymorphism (SNP), 300,000 insertion and deletion (InDel), 14,000 structural variation (SV), and 80,000 simple sequence repeat (SSR) markers (one SSR per 4.23 kb), as well as their primers. Among the SSRs, 44 SSR primer pairs were randomly selected and validated by polymerase chain reaction (PCR), allowing discrimination between the different Munake cultivar genotypes. The genetic data provided allow a deeper understanding of Munake cultivar genomic sequence and contribute to better knowledge of the genetic basis behind its key agronomic traits.

Fine Mapping of Leaf Trichome Density Revealed a 747-kb Region on Chromosome 1 in Cold-Hardy Hybrid Wine Grape Populations

Segregation for leaf trichome density was observed in a cold-hardy hybrid grape population GE1025 (N = ∼125, MN1264 × MN1246) that was previously used to detect a quantitative trait locus (QTL) underlying foliar phylloxera resistance on chromosome 14. Our hypothesis was that high trichome density was associated with resistance to phylloxera. Existing literature found trichome density QTL on chromosomes 1 and 15 using a hybrid grape population of "Horizon" × Illinois 547-1 and suggested a few candidate genes. To validate the reported QTL and our hypothesis, interval mapping was conducted in GE1025 with previous genotyping-by-sequencing (GBS) single nucleotide polymorphism (SNP) genotype data and phenotypic scores collected using a 0-6 trichome density scale at several leaf positions. Evaluations were done on replicated forced dormant cuttings in 2 years and on field-grown leaves in 1 year. There was no strong relationship between trichome density and phylloxera resistance except for a Pearson's correlation (r) of about -0.2 between a few trichome density traits and phylloxera severity traits at 2 and 3 weeks after infestation. Two genetic regions were repeatedly detected for multiple trichome density traits: from 10 to 20.7 Mbp (∼10 Mbp) on chromosome 1 for ribbon and simple density traits and from 2.4 to 8.9 Mbp on chromosome 10 for ribbon density traits, explaining 12.1-48.2 and 12.6-27.5% of phenotypic variation, respectively. To fine map, we genotyped a larger population, GE1783 (N = ∼1,023, MN1264 × MN1246), with conserved rhAmpSeq haplotype markers across multiple Vitis species and phenotyped 233 selected potential recombinants. Evaluations were conducted on field-grown leaves in a single year. The QTL for ribbon trichome density on adaxial vein and adaxial leaf and simple density on abaxial vein was fine mapped to 12.63-13.38 Mbp (747 kb) on chromosome 1. We found variations of MN1264 and MN1246 at candidate genes NAC transcription factor 29, EF-hand protein, and MYB140 in this region and three other surrounding candidate genes proposed previously. Even though no strong relationship between foliar phylloxera resistance and trichome density was found, this study validated and fine mapped a major QTL for trichome density using a cold-hardy hybrid grape population and shed light on a few candidate genes that have implications for different breeding programs.

Identification of SNPs associated with magnesium and sodium uptake and the effect of their accumulation on micro and macro nutrient levels in Vitis vinifera

Macro and micro nutrient accumulation affects all stages of plant growth and development. When nutrient deficiencies or excesses occur, normal plant growth is altered resulting in symptoms such as leaf chlorosis, plant stunting or death. In grapes, few genomic regions associated with nutrient accumulation or deficiencies have been identified. Our study evaluated micro and macro nutrient concentrations in Vitis vinifera L. to identify associated SNPs using an association approach with genotype by sequencing data. Nutrient concentrations and foliar symptoms (leaf chlorosis and stunting) were compared among 249 F₁Vitis vinifera individuals in 2015 and 2016. Foliar symptoms were consistent (≥90%) between years and correlated with changes in nutrient concentrations of magnesium (r = 0.65 and r = 0.38 in 2015 and 2016, respectively), aluminum (r = 0.24 and r = 0.49), iron (r = 0.21 and r = 0.49), and sodium (r = 0.32 and r = 0.21). Single nucleotide polymorphisms associated with symptoms, sodium, and magnesium were detected on each chromosome with the exception of 5, 7 and 17 depending on the trait and genome used for analyses explaining up to 40% of the observed variation. Symptoms and magnesium concentration were primarily associated with SNPs on chromosome 3, while SNPs associated with increased sodium content were primarily found on chromosomes 11 and 18. Mean concentrations for each nutrient varied between years in the population between symptomatic and asymptomatic plants, but relative relationships were mostly consistent. These data suggest a complex relationship among foliar symptoms and micro and macro nutrients accumulating in grapevines.

Discovery of the REN11 Locus From Vitis aestivalis for Stable Resistance to Grapevine Powdery Mildew in a Family Segregating for Several Unstable and Tissue-Specific Quantitative Resistance Loci

Race-specific resistance loci, whether having qualitative or quantitative effects, present plant-breeding challenges for phenotypic selection and deciding which loci to select or stack with other resistance loci for improved durability. Previously, resistance to grapevine powdery mildew (GPM, caused by Erysiphe necator) was predicted to be conferred by at least three race-specific loci in the mapping family B37-28 × C56-11 segregating for GPM resistance from Vitis aestivalis. In this study, 9 years of vineyard GPM disease severity ratings plus a greenhouse and laboratory assays were genetically mapped, using a rhAmpSeq core genome marker platform with 2,000 local haplotype markers. A new qualitative resistance locus, named REN11, on the chromosome (Chr) 15 was found to be effective in nearly all (11 of 12) vineyard environments on leaves, rachis, berries, and most of the time (7 of 12) stems. REN11 was independently validated in a pseudo-testcross with the grandparent source of resistance, "Tamiami." Five other loci significantly predicted GPM severity on leaves in only one or two environments, which could indicate race-specific resistance or their roles in different timepoints in epidemic progress. Loci on Chr 8 and 9 reproducibly predicted disease severity on stems but not on other tissues and had additive effects with REN11 on the stems. The rhAmpSeq local haplotype sequences published in this study for REN11 and Chr 8 and 9 stem quantitative trait locus (QTL) can be used directly for marker-assisted selection or converted to SNP assays. In screening for REN11 in a diversity panel of 20,651 vines representing the diversity of Vitis, this rhAmpSeq haplotype had a false positive rate of 0.034% or less. The effects of the other foliar resistance loci detected in this study seem too unstable for genetic improvement regardless of quantitative effect size, whether due to race specificity or other environmental variables.

Transcriptome-Derived Amplicon Sequencing Markers Elucidate the U.S. Podosphaera macularis Population Structure Across Feral and Commercial Plantings of Humulus lupulus

Obligately biotrophic plant pathogens pose challenges in population genetic studies due to their genomic complexities and elaborate culturing requirements with limited biomass. Hop powdery mildew (Podosphaera macularis) is an obligately biotrophic ascomycete that threatens sustainable hop production. P. macularis populations of the Pacific Northwest (PNW) United States differ from those of the Midwest and Northeastern United States, lacking one of two mating types needed for sexual recombination and harboring two strains that are differentially aggressive on the cultivar Cascade and able to overcome the Humulus lupulus R-gene R6 (V6), respectively. To develop a high-throughput marker platform for tracking the flow of genotypes across the United States and internationally, we used an existing transcriptome of diverse P. macularis isolates to design a multiplex of 54 amplicon sequencing markers, validated across a panel of 391 U.S. samples and 123 international samples. The results suggest that P. macularis from U.S. commercial hop yards form one population closely related to P. macularis of the United Kingdom, while P. macularis from U.S. feral hop locations grouped with P. macularis of Eastern Europe. Included in this multiplex was a marker that successfully tracked V6-virulence in 65 of 66 samples with a confirmed V6-phenotype. A new qPCR assay for high-throughput genotyping of P. macularis mating type generated the highest resolution distribution map of P. macularis mating type to date. Together, these genotyping strategies enable the high-throughput and inexpensive tracking of pathogen spread among geographical regions from single-colony samples and provide a roadmap to develop markers for other obligate biotrophs.

Population Genomics of Filamentous Plant Pathogens-A Brief Overview of Research Questions, Approaches, and Pitfalls

With ever-decreasing sequencing costs, research on the population biology of plant pathogens is transitioning from population genetics-using dozens of genetic markers or polymorphism data of several genes-to population genomics-using several hundred to tens of thousands of markers or whole-genome sequence data. The field of population genomics is characterized by rapid theoretical and methodological advances and by numerous steps and pitfalls in its technical and analytical workflow. In this article, we aim to provide a brief overview of topics relevant to the study of population genomics of filamentous plant pathogens and direct readers to more extensive reviews for in-depth understanding. We briefly discuss different types of population genomics-inspired research questions and give insights into the sampling strategies that can be used to answer such questions. We then consider different sequencing strategies, the various options available for data processing, and some of the currently available tools for population genomic data analysis. We conclude by highlighting some of the hurdles along the population genomic workflow, providing cautionary warnings relative to assumptions and technical challenges, and presenting our own future perspectives of the field of population genomics for filamentous plant pathogens.

Draft genome of the Native American cold hardy grapevine Vitis riparia Michx. 'Manitoba 37'

Vitis riparia, a critically important Native American grapevine species, is used globally in rootstock and scion breeding and contributed to the recovery of the French wine industry during the mid-19th century phylloxera epidemic. This species has abiotic and biotic stress tolerance and the largest natural geographic distribution of the North American grapevine species. Here we report an Illumina short-read 369X coverage, draft de novo heterozygous genome sequence of V. riparia Michx. 'Manitoba 37' with the size of ~495 Mb for 69,616 scaffolds and a N50 length of 518,740 bp. Using RNAseq data, 40,019 coding sequences were predicted and annotated. Benchmarking with Universal Single-Copy Orthologs (BUSCO) analysis of predicted gene models found 96% of the complete BUSCOs in this assembly. The assembly continuity and completeness were further validated using V. riparia ESTs, BACs, and three de novo transcriptome assemblies of three different V. riparia genotypes resulting in >98% of respective sequences/transcripts mapping with this assembly. Alignment of the V. riparia assembly and predicted CDS with the latest V. vinifera 'PN40024' CDS and genome assembly showed 99% CDS alignment and a high degree of synteny. An analysis of plant transcription factors indicates a high degree of homology with the V. vinifera transcription factors. QTL mapping to V. riparia 'Manitoba 37' and V. vinifera PN40024 has identified genetic relationships to phenotypic variation between species. This assembly provides reference sequences, gene models for marker development and understanding V. riparia's genetic contributions in grape breeding and research.