Comparison of Three Assembly Strategies for a Heterozygous Seedless Grapevine Genome Assembly

BackgroundDe novo heterozygous assembly is an ongoing challenge requiring improved assembly approaches. In this study, three strategies were used to develop de novo Vitis vinifera ‘Sultanina’ genome assemblies for comparison with the inbred V. vinifera (PN40024 12X.v2) reference genome and a published Sultanina ALLPATHS-LG assembly (AP). The strategies were: 1) a default PLATANUS assembly (PLAT_d) for direct comparison with AP assembly, 2) an iterative merging strategy using METASSEMBLER to combine PLAT_d and AP assemblies (MERGE) and 3) PLATANUS parameter modifications plus GapCloser (PLAT*_GC).ResultsThe three new assemblies were greater in size than the AP assembly. PLAT*_GC had the greatest number of scaffolds aligning with a minimum of 95% identity and ≥1000 bp alignment length to V. vinifera (PN40024 12X.v2) reference genome. SNP analysis also identified additional high quality SNPs. A greater number of sequence reads mapped back with zero-mismatch to the PLAT_d, MERGE, and PLAT*_GC (\u3e94%) than was found in the AP assembly (87%) indicating a greater fidelity to the original sequence data in the new assemblies than in AP assembly. A de novo gene prediction conducted using seedless RNA-seq data predicted \u3e 30,000 coding sequences for the three new de novo assemblies, with the greatest number (30,544) in PLAT*_GC and only 26,515 for the AP assembly. Transcription factor analysis indicated good family coverage, but some genes found in the VCOST.v3 annotation were not identified in any of the de novo assemblies, particularly some from the MYB and ERF families.ConclusionsThe PLAT_d and PLAT*_GC had a greater number of synteny blocks with the V. vinifera (PN40024 12X.v2) reference genome than AP or MERGE. PLAT*_GC provided the most contiguous assembly with only 1.2% scaffold N, in contrast to AP (10.7% N), PLAT_d (6.6% N) and Merge (6.4% N). A PLAT*_GC pseudo-chromosome assembly with chromosome alignment to the reference genome V. vinifera, (PN40024 12X.v2) provides new information for use in seedless grape genetic mapping studies. An annotated de novo gene prediction for the PLAT*_GC assembly, aligned with VitisNet pathways provides new seedless grapevine specific transcriptomic resource that has excellent fidelity with the seedless short read sequence data

Transcriptomic Response is More Sensitive to Water Deficit in Shoots than Roots of \u3cem\u3eVitis riparia\u3c/em\u3e (Michx.)

Background: Drought is an important constraint on grapevine sustainability. Vitis riparia, widely used in rootstock and scion breeding, has been studied in isolated leaf drying response studies; however, it is essential to identify key root and shoot water deficit signaling traits in intact plants. This information will aid improved scion and rootstock selection and management practices in grapevine. RNAseq data were generated from V. riparia roots and shoots under water deficit and well-watered conditions to determine root signaling and shoot responses to water deficit. Results: Shoot elongation, photosynthetic rate, and stomatal conductance were significantly reduced in water deficit (WD) treated than in well-watered grapevines. RNAseq analysis indicated greater transcriptional differences in shoots than in roots under WD, with 6925 and 1395 genes differentially expressed, respectively (q-value \u3c 0.05). There were 50 and 25 VitisNet pathways significantly enriched in WD relative to well-watered treatments in grapevine shoots and roots, respectively. The ABA biosynthesis genes beta-carotene hydroxylase, zeaxanthin epoxidase, and 9-cis-epoxycarotenoid dioxygenases were up-regulated in WD root and WD shoot. A positive enrichment of ABA biosynthesis genes and signaling pathways in WD grapevine roots indicated enhanced root signaling to the shoot. An increased frequency of differentially expressed reactive oxygen species scavenging (ROS) genes were found in the WD shoot. Analyses of hormone signaling genes indicated a strong ABA, auxin, and ethylene network and an ABA, cytokinin, and circadian rhythm network in both WD shoot and WD root. Conclusions: This work supports previous findings in detached leaf studies suggesting ABA-responsive binding factor 2 (ABF2) is a central regulator in ABA signaling in the WD shoot. Likewise, ABF2 may have a key role in V. riparia WD shoot and WD root. A role for ABF3 was indicated only in WD root. WD shoot and WD root hormone expression analysis identified strong ABA, auxin, ethylene, cytokinin, and circadian rhythm signaling networks. These results present the first ABA, cytokinin, and circadian rhythm signaling network in roots under water deficit. These networks point to organ specific regulators that should be explored to further define the communication network from soil to shoo

VitisNet: ‘‘Omics’’ Integration through Grapevine Molecular Networks

Background: Genomic data release for the grapevine has increased exponentially in the last five years. The Vitis vinifera genome has been sequenced and Vitis EST, transcriptomic, proteomic, and metabolomic tools and data sets continue to be developed. The next critical challenge is to provide biological meaning to this tremendous amount of data by annotating genes and integrating them within their biological context. We have developed and validated a system of Grapevine Molecular Networks (VitisNet). Methodology/Principal Findings: The sequences from the Vitis vinifera (cv. Pinot Noir PN40024) genome sequencing project and ESTs from the Vitis genus have been paired and the 39,424 resulting unique sequences have been manually annotated. Among these, 13,145 genes have been assigned to 219 networks. The pathway sets include 88 ‘‘Metabolic’’, 15 ‘‘Genetic Information Processing’’, 12 ‘‘Environmental Information Processing’’, 3 ‘‘Cellular Processes’’, 21 ‘‘Transport’’, and 80 ‘‘Transcription Factors’’. The quantitative data is loaded onto molecular networks, allowing the simultaneous visualization of changes in the transcriptome, proteome, and metabolome for a given experiment. Conclusions/Significance: VitisNet uses manually annotated networks in SBML or XML format, enabling the integration of large datasets, streamlining biological functional processing, and improving the understanding of dynamic processes in systems biology experiments. VitisNet is grounded in the Vitis vinifera genome (currently at 8x coverage) and can be readily updated with subsequent updates of the genome or biochemical discoveries. The molecular network files can be dynamically searched by pathway name or individual genes, proteins, or metabolites through the MetNet Pathway database and web-portal at http://metnet3.vrac.iastate.edu