Haplotype-resolved powdery mildew resistance loci reveal the impact of heterozygous structural variation on NLR genes in Muscadinia rotundifolia.

Muscadinia rotundifolia cv. Trayshed is a valuable source of resistance to grape powdery mildew. It carries 2 powdery mildew resistance-associated genetic loci, Run1.2 on chromosome 12 and Run2.2 on chromosome 18. The purpose of this study was to identify candidate resistance genes associated with each haplotype of the 2 loci. Both haplotypes of each resistance-associated locus were identified, phased, and reconstructed. Haplotype phasing allowed the identification of several structural variation events between haplotypes of both loci. Combined with a manual refinement of the gene models, we found that the heterozygous structural variants affected the gene content, with some resulting in duplicated or hemizygous nucleotide-binding leucine-rich repeat genes. Heterozygous structural variations were also found to impact the domain composition of some nucleotide-binding leucine-rich repeat proteins. By comparing the nucleotide-binding leucine-rich repeat proteins at Run1.2 and Run2.2 loci, we discovered that the 2 loci include different numbers and classes of nucleotide-binding leucine-rich repeat genes. To identify powdery mildew resistance-associated genes, we performed a gene expression profiling of the nucleotide-binding leucine-rich repeat genes at Run1.2b and Run2.2 loci with or without powdery mildew present. Several nucleotide-binding leucine-rich repeat genes were constitutively expressed, suggesting a role in powdery mildew resistance. These first complete, haplotype-resolved resistance-associated loci and the candidate nucleotide-binding leucine-rich repeat genes identified by this study are new resources that can aid the development of powdery mildew-resistant grape cultivars.

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.

Rootstock influences the effect of grapevine leafroll-associated viruses on berry development and metabolism via abscisic acid signalling.

Grapevine leafroll-associated virus (GLRaV) infections are accompanied by symptoms influenced by host genotype, rootstock, environment, and which individual or combination of GLRaVs is present. Using a dedicated experimental vineyard, we studied the responses to GLRaVs in ripening berries from Cabernet Franc grapevines grafted to different rootstocks and with zero, one, or pairs of leafroll infection(s). RNA sequencing data were mapped to a high-quality Cabernet Franc genome reference assembled to carry out this study and integrated with hormone and metabolite abundance data. This study characterized conserved and condition-dependent responses to GLRaV infection(s). Common responses to GLRaVs were reproduced in two consecutive years and occurred in plants grafted to different rootstocks in more than one infection condition. Though different infections were inconsistently distinguishable from one another, the effects of infections in plants grafted to different rootstocks were distinct at each developmental stage. Conserved responses included the modulation of genes related to pathogen detection, abscisic acid (ABA) signalling, phenylpropanoid biosynthesis, and cytoskeleton remodelling. ABA, ABA glucose ester, ABA and hormone signalling-related gene expression, and the expression of genes in several transcription factor families differentiated the effects of GLRaVs in berries from Cabernet Franc grapevines grafted to different rootstocks. These results support that ABA participates in the shared responses to GLRaV infection and differentiates the responses observed in grapevines grafted to different rootstocks.

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.

The genetic basis of sex determination in grapes.

It remains a major challenge to identify the genes and mutations that lead to plant sexual differentiation. Here, we study the structure and evolution of the sex-determining region (SDR) in Vitis species. We report an improved, chromosome-scale Cabernet Sauvignon genome sequence and the phased assembly of nine wild and cultivated grape genomes. By resolving twenty Vitis SDR haplotypes, we compare male, female, and hermaphrodite haplotype structures and identify sex-linked regions. Coupled with gene expression data, we identify a candidate male-sterility mutation in the VviINP1 gene and potential female-sterility function associated with the transcription factor VviYABBY3. Our data suggest that dioecy has been lost during domestication through a rare recombination event between male and female haplotypes. This work significantly advances the understanding of the genetic basis of sex determination in Vitis and provides the information necessary to rapidly identify sex types in grape breeding programs.

The genomic diversification of grapevine clones.

BACKGROUND: Vegetatively propagated clones accumulate somatic mutations. The purpose of this study was to better appreciate clone diversity and involved defining the nature of somatic mutations throughout the genome. Fifteen Zinfandel winegrape clone genomes were sequenced and compared to one another using a highly contiguous genome reference produced from one of the clones, Zinfandel 03. RESULTS: Though most heterozygous variants were shared, somatic mutations accumulated in individual and subsets of clones. Overall, heterozygous mutations were most frequent in intergenic space and more frequent in introns than exons. A significantly larger percentage of CpG, CHG, and CHH sites in repetitive intergenic space experienced transition mutations than in genic and non-repetitive intergenic spaces, likely because of higher levels of methylation in the region and because methylated cytosines often spontaneously deaminate. Of the minority of mutations that occurred in exons, larger proportions of these were putatively deleterious when they occurred in relatively few clones. CONCLUSIONS: These data support three major conclusions. First, repetitive intergenic space is a major driver of clone genome diversification. Second, clones accumulate putatively deleterious mutations. Third, the data suggest selection against deleterious variants in coding regions or some mechanism by which mutations are less frequent in coding than noncoding regions of the genome.

Iso-Seq Allows Genome-Independent Transcriptome Profiling of Grape Berry Development.

Transcriptomics has been widely applied to study grape berry development. With few exceptions, transcriptomic studies in grape are performed using the available genome sequence, PN40024, as reference. However, differences in gene content among grape accessions, which contribute to phenotypic differences among cultivars, suggest that a single reference genome does not represent the species' entire gene space. Though whole genome assembly and annotation can reveal the relatively unique or "private" gene space of any particular cultivar, transcriptome reconstruction is a more rapid, less costly, and less computationally intensive strategy to accomplish the same goal. In this study, we used single molecule-real time sequencing (SMRT) to sequence full-length cDNA (Iso-Seq) and reconstruct the transcriptome of Cabernet Sauvignon berries during berry ripening. In addition, short reads from ripening berries were used to error-correct low-expression isoforms and to profile isoform expression. By comparing the annotated gene space of Cabernet Sauvignon to other grape cultivars, we demonstrate that the transcriptome reference built with Iso-Seq data represents most of the expressed genes in the grape berries and includes 1,501 cultivar-specific genes. Iso-Seq produced transcriptome profiles similar to those obtained after mapping on a complete genome reference. Together, these results justify the application of Iso-Seq to identify cultivar-specific genes and build a comprehensive reference for transcriptional profiling that circumvents the necessity of a genome reference with its associated costs and computational weight.

Metabolite Profiling Reveals Developmental Inequalities in Pinot Noir Berry Tissues Late in Ripening.

Uneven ripening in Vitis vinifera is increasingly recognized as a phenomenon of interest, with substantial implications for fruit and wine composition and quality. This study sought to determine whether variation late in ripening (∼Modified Eichhorn-Lorenz stage 39) was associated with developmental differences that were observable as fruits within a cluster initiated ripening (véraison). Four developmentally distinct ripening classes of berries were tagged at cluster véraison, sampled at three times late in ripening, and subjected to untargeted HPLC-MS to measure variation in amino acids, sugars, organic acids, and phenolic metabolites in skin, pulp, and seed tissues separately. Variability was described using predominantly two strategies. In the first, multivariate analysis (Orthogonal Projections to Latent Structures-Discriminant Analysis, OPLS-DA) was used to determine whether fruits were still distinguishable per their developmental position at véraison and to identify which metabolites accounted for these distinctions. The same technique was used to assess changes in each tissue over time. In a second strategy and for each annotated metabolite, the variance across the ripening classes at each time point was measured to show whether intra-cluster variance (ICV) was growing, shrinking, or constant over the period observed. Indeed, berries could be segregated by OPLS-DA late in ripening based on their developmental position at véraison, though the four ripening classes were aggregated into two larger ripening groups. Further, not all tissues were dynamic over the period examined. Although pulp tissues could be segregated by time sampled, this was not true for seed and only moderately so for skin. Ripening group differences in seed and skin, rather than the time fruit was sampled, were better able to define berries. Metabolites also experienced significant reductions in ICV between single pairs of time points, but never across the entire experiment. Metabolites often exhibited a combination of ICV expansion, contraction and persistence. Finally, we observed significant differences in the abundance of some metabolites between ripening classes that suggest the berries that initiated ripening first remained developmentally ahead of the lagging fruit even late in the ripening phase. This presents a challenge to producers who would seek to harvest at uniformity or at a predefined level of variation.

The contribution of flowering time and seed content to uneven ripening initiation among fruits within Vitis vinifera L. cv. Pinot noir clusters.

MAIN CONCLUSION: Ripening initiation-associated hormonal changes and sugar accumulation for individual fruits differed by seed content and did not depend heavily on flowering time or duration from anthesis to clusters' onset of ripening. For Vitis vinifera, the ripening initiation of individual fruits in a cluster occurs unevenly. This developmental period is called véraison. Why individual fruits initiate ripening at different times is not well studied, though differences in seed content and unequal developmental durations that arise from asynchronous flowering within a cluster have been proposed. This study examined how much both variables contribute to individual fruits' ripening progress by mid-véraison, when half of berries in a cluster have initiated ripening, and whether either or both factors affect the timing of characteristic, ripening-initiation associated changes in abscisic acid and auxin before, at, and after véraison. Overall, developmental duration and flowering time did not sufficiently explain how far berries had progressed into the ripening stage because fruits did not require a fixed amount of time to initiate ripening. Fruits from early and late flowers but of similar seed content were able to initiate ripening at the same time despite differences in chronological age. This suggests either an early developmental enhancement occurred for late-initiated fruits or that flowering time is an inappropriate "day zero". Ultimately, only seed content was linked to the timing and magnitude of ripening-related hormone changes, supporting that seeds have a comparatively larger influence than flowering time on the ripening initiation of individual berries. More specifically, if the fraction of berry weight occupied by seed was high, then the initiation of ripening for that berry and its associated hormone changes were delayed relative to berries with less seed weight versus total berry weight.