Search

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.

Xylella fastidiosais a bacterium that infects crops like grapevines, coffee, almonds, citrus and olives, causing economically devastating damage. There is, however, little understanding of the genes that contribute to resistance, the genomic architecture of resistance, and the potential role of climate in shaping resistance, in part because major crops like grapevines (V. vinifera) are not resistant to the bacterium. Here we studied a wild grapevine species,Vitis arizonica, that segregates for resistance toX. fastidiosa. Using genome-wide association, we identified candidate genes that mediate the host response toX. fastidiosainfection. We uncovered evidence that resistance requires genes from multiple genomic regions, based on data from breeding populations and from additionalVitisspecies. We also inferred that resistance evolved more than once in the wild, suggesting that wildVitisspecies may be a rich source for resistance alleles and mechanisms. Finally, resistance inV. arizonicawas 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 these climates. Surprisingly, climate was nearly as effective a predictor of resistance phenotypes as some genetic markers. This work underscores that pathogen pressure is likely to increase with climate, but it also provides genetic insight and tools for breeding and transforming resistant crops.

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.

Muscadinia rotundifolia cv. Trayshed is a valuable source of resistance to grape powdery mildew. It carries two powdery mildew resistance-associated genetic loci, Run1.2 on chromosome 12 and Run2.2 on chromosome 18. In this study, we identified, phased, and reconstructed the two haplotypes of each resistance-associated locus. Haplotype phasing allowed the identification of several structural variation events between haplotypes of both loci. Combined with 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 (NLR) genes. Structural variations also impacted the domain composition of some NLRs. These findings emphasize the need of generating haplotype-resolved sequences instead of using consensus sequences for identifying haplotype-specific candidate genes. Comparison of the NLRs in the Run1.2 and Run2.2 loci indicated that the two loci are composed of a different number and classes of NLR genes. We provide a list of candidate NLR genes from the Run1.2b and Run2.2 loci, whose expression suggests a role in powdery mildew resistance in Trayshed. These first complete and haplotype-resolved resistance-associated loci, and their candidate NLR genes, represent new resources to develop powdery mildew-resistant grape cultivars.

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.

Background ‘Nebbiolo’ is a grapevine cultivar typical of north-western Italy, appreciated for producing high-quality red wines. Grapevine cultivars are characterized by possessing highly heterozygous genomes, including a great incidence of genomic rearrangements larger than 50 bp, so called structural variations (SVs). Even though abundant, SVs are an under-explored source of genetic variation mainly due to methodological limitations at their detection. Results We employed a multiple platform approach to produce long-range genomic data for two different ‘Nebbiolo’ clones, namely: optical mapping, long-reads and linked-reads. We performed a haplotype-resolved de novo assembly for cultivar ‘Nebbiolo’ (clone CVT 71) and used an ab-initio strategy to annotate it. The annotated assembly enhanced our ability to detect SVs, enabling the study of genomic regions not present in the grapevines’ reference genome and accounting for their functional implications. We performed variant calling analyses at three different organizational levels: i) between haplotypes of clone CVT 71 (primary assembly vs haplotigs), ii) between ‘Nebbiolo’ and ‘Cabernet Sauvignon’ assemblies and iii) between clones CVT 71 and CVT 185, representing different ‘Nebbiolo’ biotypes. The cumulative size of non-redundant merged SVs indicated a total of 79.6 Mbp for the first comparison and 136.1 Mbp for the second one, while no SVs were detected for the third comparison. Interestingly, SVs differentiating cultivars and haplotypes affected similar numbers of coding genes. Conclusions Our results suggest that SVs accumulation rate and their functional implications in ‘Nebbiolo’ genome are highly-dependent on the organizational level under study. SVs are abundant when comparing ‘Nebbiolo’ to a different cultivar or the two haplotypes of the same individual, while they turned absent between the two analysed clones.

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., GLRaVs elicit changes in gene expression, ripening‐related metabolites, and hormone metabolism in berries; one set of responses is observed across infections and host genotypes while other responses differ based on rootstock. ​

Structural variants (SVs) are a largely unexplored feature of plant genomes. Little is known about the type and size of SVs, their distribution among individuals and, especially, their population dynamics. Understanding these dynamics is critical for understanding both the contributions of SVs to phenotypes and the likelihood of identifying them as causal genetic variants in genome-wide associations. Here, we identify SVs and study their evolutionary genomics in clonally propagated grapevine cultivars and their outcrossing wild progenitors. To catalogue SVs, we assembled the highly heterozygous Chardonnay genome, for which one in seven genes is hemizygous based on SVs. Using an integrative comparison between Chardonnay and Cabernet Sauvignon genomes by whole-genome, long-read and short-read alignment, we extended SV detection to population samples. We found that strong purifying selection acts against SVs but particularly against inversion and translocation events. SVs nonetheless accrue as recessive heterozygotes in clonally propagated lineages. They also define outlier regions of genomic divergence between wild and cultivated grapevines, suggesting roles in domestication. Outlier regions include the sex-determination region and the berry colour locus, where independent large, complex inversions have driven convergent phenotypic evolution. This study identified structural variants in grapevine populations, including wild progenitors, using a de novo assembly and comparative genomics approach, and examined their evolutionary genomics and roles in domestication and phenotypic evolution.

Significance 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.

Muscadinia rotundifolia, the muscadine grape, has been cultivated for centuries in the southeastern United States.M. rotundifoliais resistant to many of the pathogens that detrimentally affectVitis vinifera, the grape species commonly used for winemaking. For this reason,M. rotundifoliais 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 20M. rotundifoliacv. Trayshed chromosomes. The completeness and accuracy of the assembly were confirmed using a high-density linkage map ofM. rotundifolia.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 thatMuscadiniachromosomes 7 and 20 are fused inVitisand 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 theResistance 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 for Trayshed, its annotation, and an associated Blast tool are available atwww.grapegenomics.com.

The Botryosphaeriaceae is a fungal family that includes many destructive vascular pathogens of woody plants (e.g., Botryosphaeria dieback of grape, Panicle blight of pistachio). Species in the genera Botryosphaeria, Diplodia, Dothiorella, Lasiodiplodia, Neofusicoccum, and Neoscytalidium attack a range of horticultural crops, but they vary in virulence and their abilities to infect their hosts via different infection courts (flowers, green shoots, woody twigs). Isolates of seventeen species, originating from symptomatic apricot, grape, pistachio, and walnut were tested for pathogenicity on grapevine wood after 4 months of incubation in potted plants in the greenhouse. Results revealed significant variation in virulence in terms of the length of the internal wood lesions caused by these seventeen species. Phylogenomic comparisons of the seventeen species of wood-colonizing fungi revealed clade-specific expansion of gene families representing putative virulence factors involved in toxin production and mobilization, wood degradation, and nutrient uptake. Statistical analyses of the evolution of the size of gene families revealed expansions of secondary metabolism and transporter gene families in Lasiodiplodia and of secreted cell wall degrading enzymes (CAZymes) in Botryosphaeria and Neofusicoccum genomes. In contrast, Diplodia, Dothiorella, and Neoscytalidium generally showed a contraction in the number of members of these gene families. Overall, species with expansions of gene families, such as secreted CAZymes, secondary metabolism, and transporters, were the most virulent (i.e., were associated with the largest lesions), based on our pathogenicity tests and published reports. This study represents the first comparative phylogenomic investigation into the evolution of possible virulence factors from diverse, cosmopolitan members of the Botryosphaeriaceae.

The Botryosphaeriaceae is a fungal family that includes many destructive vascular pathogens of woody plants (e.g., Botryosphaeria dieback of grape, Panicle blight of pistachio). Species in the genera Botryosphaeria, Diplodia, Dothiorella, Lasiodiplodia, Neofusicoccum, and Neoscytalidium attack a range of horticultural crops, but they vary in virulence and in their abilities to infect their hosts via different infection courts (flowers, green shoots, woody twigs). Isolates of seventeen species, originating from symptomatic apricot, grape, pistachio, and walnut were tested for pathogenicity to grapevine wood after four months of incubation in potted plants in the greenhouse. Results revealed significant variation in virulence in terms of the length of the internal wood lesions caused by these seventeen species. Phylogenomic comparisons of the seventeen species of wood-colonizing fungi revealed clade-specific expansion of gene families representing putative virulence factors involved in toxin production and mobilization, wood degradation, and nutrient uptake. Statistical analyses of the evolution of the size of gene families revealed expansions of secondary metabolism and transporter gene families in Lasiodiplodia and of secreted cell wall degrading enzymes (CAZymes) in Botryosphaeria and Neofusicoccum genomes. In contrast, Diplodia, Dothiorella, and Neoscytalidium generally showed a contraction in the number of members of these gene families. Overall, species with expansions of gene families, such as secreted CAZymes, secondary metabolism, and transporters, were the most virulent (i.e., were associated with the largest lesions), based on our pathogenicity tests and published reports. This study represents the first comparative phylogenomic investigation into the evolution of possible virulence factors from diverse, cosmopolitan members of the Botryosphaeriaceae.

Structural Variants (SVs) are a widely unexplored source of genetic variation, both due to methodological limitations and because they are generally associated to deleterious effects. However, with the advent of long-range genomic platforms, it has become easier to directly detect SVs. In the same direction, clonally propagated crops provide a unique opportunity to study SVs, offering a suitable genomic environment for their accumulation in heterozygosis. In particular, it has been reported that SVs generate drastic levels of heterozygosity in grapevines. ‘Nebbiolo’ (Vitis vinifera L.) is a grapevine cultivar typical of north-western Italy, appreciated for its use in producing high-quality red wines. Here, we aimed to analyze the frequency of SVs in ‘Nebbiolo’, at three different organizational levels. For this purpose, we generated genomic data based on long-reads, linked-reads and optical mapping. We assembled a reference genome for this cultivar and compared two different clones, including V. vinifera reference genome (PN40024) in our comparisons. Our results indicate that SVs differentially occurring between ‘Nebbiolo’ clones might be rare, while SVs differentiating haplotypes of the same individual are as abundant as those that occur differentially between cultivars.

Phytophthora megakarya and P. palmivora are oomycete pathogens that cause black pod rot of cacao (Theobroma cacao), the most economically important disease on cacao globally. While P. palmivora is a cosmopolitan pathogen, P. megakarya, which is more aggressive on cacao than P. palmivora, has been reported only in West and Central Africa where it has been spreading and devastating cacao farms since the 1950s. In this study, we reconstructed the complete diploid genomes of multiple isolates of both species using single-molecule real-time sequencing. Thirty-one additional genotypes were sequenced to analyze inter- and intra-species genomic diversity. The P. megakarya genome is exceptionally large (222 Mbp) and nearly twice the size of P. palmivora (135 Mbp) and most known Phytophthora species (∼100 Mbp on average). Previous reports pointed toward a whole-genome duplication (WGD) in P. palmivora. In this study, we demonstrate that both species underwent independent and relatively recent WGD events. In P. megakarya we identified a unique combination of WGD and large-scale transposable element driven genome expansion, which places this genome in the upper range of Phytophthora genome sizes, as well as effector pools with 1,382 predicted RxLR effectors. Finally, this study provides evidence of adaptive evolution of effectors like RxLRs and Crinklers, and discusses the implications of effector expansion and diversification.

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., Grapevine is one of a few ancestrally dioecious crops that are reverted to hermaphroditism during domestication. Here, the authors identify candidate genes related to male- and female-sterility in grapes and describe the genetic process that led to hermaphroditism during domestication.

Sex determination in grapevine evolved through a complex succession of switches in sexual systems. Phased genomes built with single molecule real-time sequencing reads were assembled for eleven accessions of cultivated hermaphrodite grapevines and dioecious males and females, including the ancestor of domesticated grapevine and other related wild species. By comparing the phased sex haplotypes, we defined the sex locus of theVitisgenus and identified polymorphisms spanning regulatory and coding sequences that are in perfect association with each sex-type throughout the genus. These findings identified a novel male-fertility candidate gene,INP1, and significantly refined the model of sex determination inVitisand its evolution.

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.

‘Nebbiolo’ (Vitis vinifera) is among the most ancient and prestigious wine grape varieties characterised by a wide genetic variability exhibited by a high number of clones (vegetatively propagated lines of selected mother plants). However, limited information is available for this cultivar at the molecular and genomic levels. The whole-genomes of three ‘Nebbiolo’ clones (CVT 71, CVT 185 and CVT 423) were re-sequenced and a de novo transcriptome assembly was produced. Important remarks about the genetic peculiarities of ‘Nebbiolo’ and its intra-varietal variability useful for clonal identification were reported. In particular, several varietal transcripts identified for the first time in ‘Nebbiolo’ were disease resistance genes and single-nucleotide variants (SNVs) identified in ‘Nebbiolo’, but not in other cultivars, were associated with genes involved in the stress response. Ten newly discovered SNVs were successfully employed to identify some periclinal chimeras and to classify 98 ‘Nebbiolo’ clones in seven main genotypes, which resulted to be linked to the geographical origin of accessions. In addition, for the first time it was possible to discriminate some ‘Nebbiolo’ clones from the others.

Phytophthora megakarya (Pmeg) and P. palmivora (Ppal) are oomycete pathogens that cause black pod rot of cacao (Theobroma cacao), the most economically important disease on cacao globally. While Ppal is a cosmopolitan pathogen, Pmeg, which is more aggressive on cacao than Ppal, has been reported only in West and Central Africa where it has been spreading and devastating cacao farms since the 1950s. In this study, we reconstructed the complete diploid genomes of multiple isolates of both species using single-molecule sequencing. Thirty-one additional genotypes were sequenced to analyze inter- and intra-species genomic diversity. These resources make it possible to better understand the molecular basis of virulence differences in closely related and consequential pathogens and study their evolutionary history. The Pmeg genome is exceptionally large (222 Mbp) and nearly twice the size Ppal (135 Mbp) and most known Phytophthora species (∼100 Mbp on average). We show that the genomes of both species recently expanded by independent whole-genome duplications (WGD). WGD and the dramatic transposable element associated expansion of a few gene families led to the exceptionally large genome and transcriptome of Pmeg and the diversification of virulence-related genes including secreted RxLR effectors. Finally, this study provides evidence of adaptive evolution among well-known effectors and discusses the implications of effector expansion and diversification.

Vegetatively propagated clones accumulate somatic mutations. The purpose of this study was to better understand the consequences of clonal propagation 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.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 genic and non-repetitive intergenic spaces, likely because of higher levels of methylation in the region and the increased likelihood of methylated cytosines to 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.These data support three major conclusions. First, repetitive intergenic space is a major driver of clone genome diversification. Second, clonal propagation is associated with the accumulation of putatively deleterious mutations. Third, the data suggest selection against deleterious variants in coding regions such that mutations are less frequent in coding than noncoding regions of the genome.

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.

With approximately 450 species, spiny Solanum species constitute the largest monophyletic group in the Solanaceae family, but a high-quality genome assembly from this group is presently missing. We obtained a chromosome-anchored genome assembly of eggplant (Solanum melongena), containing 34,916 genes, confirming that the diploid gene number in the Solanaceae is around 35,000. Comparative genomic studies with tomato (S. lycopersicum), potato (S. tuberosum) and pepper (Capsicum annuum) highlighted the rapid evolution of miRNA:mRNA regulatory pairs and R-type defense genes in the Solanaceae, and provided a genomic basis for the lack of steroidal glycoalkaloid compounds in the Capsicum genus. Using parsimony methods, we reconstructed the putative chromosomal complements of the key founders of the main Solanaceae clades and the rearrangements that led to the karyotypes of extant species and their ancestors. From 10% to 15% of the genes present in the four genomes were syntenic paralogs (ohnologs) generated by the pre-γ, γ and T paleopolyploidy events, and were enriched in transcription factors. Our data suggest that the basic gene network controlling fruit ripening is conserved in different Solanaceae clades, and that climacteric fruit ripening involves a differential regulation of relatively few components of this network, including CNR and ethylene biosynthetic genes.

Though grape transcriptomics has expanded dramatically over the last ten years, few additional novel genomic resources were developed since the release of the PN40024 reference genome in 2007. This is partly because of the difficulty associated with assembling grape genomes. Despite a relatively small genome size of ~500 Mb and modest repeat content, high sequence and structural heterozygosity makes assembling grape genomes particularly challenging. Without assemblies representative of the genetic diversity within the cultivated germplasm, identifying cultivar-specific functions not represented in the PN40024 genome has remained elusive. Now, third-generation sequencing technologies and long-range scaffolding methods have made it possible to relatively inexpensively and rapidly generate highly contiguous and complete grape genomes. This chapter will describe the challenges associated with the isolation of high-quality nucleic acids suitable for long-read sequencing and provide an overview of the sequencing and assembling approaches that can be used to successfully reconstruct grape genomes.

Structural variants (SVs) affect plant phenotypes, but they are a largely unexplored feature of plant genomes. Little is known about the type and size of SVs, their distribution among individuals or their evolutionary dynamics. Here we identify SVs and study their evolutionary dynamics in clonally propagated grapevine cultivars and their outcrossing wild relatives. To catalog SVs, we assembled the highly heterozygous Chardonnay genome, for which one in seven genes is hemizygous. Using genomic inference as the standard, we extended SV detection to population samples. We found that negative selection acts against SVs, but particularly against inversion and translocation events. SVs nonetheless accrue as recessive heterozygotes in clonal lineages. They also define outlier regions of genomic divergence between wild and cultivated grapevines, suggesting roles in domestication. Outlier regions include the sex determination region and the berry color locus, where independent large, complex inversions drive convergent phenotypic evolution.

The role of plant tyrosyl-DNA phosphodiesterase 1α in genome stability is studied using a Medicago truncatula MtTdp1α-depleted line. Lack of MtTdp1α results in a 39% reduction of methylated cytosines as compared to control. RNA-Seq analyses revealed that 11 DNA transposons and 22 retrotransposons were differentially expressed in the Tdp1α-2a line. Among them all, DNA transposons (MuDR, hAT, DNA3-11_Mad) and seven retrotransposons (LTR (Long Terminal Repeat)/Gipsy, LTR/Copia, LTR and NonLTR/L1) were down-regulated, while the 15 retrotransposons were up-regulated. Results suggest that the occurrence of stress-responsive cis-elements as well as changes in the methylation pattern at the LTR promoters might be responsible for the enhanced retrotransposon transcription.

The Ascomycete fungusPhaeoacremonium minimumis one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence amongPm. minimumisolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference ofPm. minimumcomprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included thePm. minimumpan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

The first draft genome sequencing of the non-model fungal pathogen Pyrenochaeta lycopersici showed an expansion of gene families associated with heterokaryon incompatibility and lacking of mating-type genes, providing insights into the genetic basis of this "imperfect" fungus which lost the ability to produce the sexual stage. However, due to the Illumina short-read technology, the draft genome was too fragmented to allow a comprehensive characterization of the genome, especially of the repetitive sequence fraction. In this work, the sequencing of another P. lycopersici isolate using long-read Single Molecule Real-Time sequencing technology was performed with the aim of obtaining a gapless genome. Indeed, a gapless genome assembly of 62.7 Mb was obtained, with a fraction of repetitive sequences representing 30% of the total bases. The gene content of the two P. lycopersici isolates was very similar, and the large difference in genome size (about 8 Mb) might be attributable to the high fraction of repetitive sequences detected for the new sequenced isolate. The role of repetitive elements, including transposable elements, in modulating virulence effectors is well established in fungal plant pathogens. Moreover, transposable elements are of fundamental importance in creating and re-modelling genes, especially in imperfect fungi. Their abundance in P. lycopersici, together with the large expansion of heterokaryon incompatibility genes in both sequenced isolates, suggest the presence of possible mechanisms alternative to gene re-assorting mediated by sexual recombination. A quite large fraction (~9%) of repetitive elements in P. lycopersici, has no homology with known classes, strengthening this hypothesis. The availability of a gapless genome of P. lycopersici allowed the in-depth analysis of its genome content, by annotating functional genes and TEs. This goal will be an important resource for shedding light on the evolution of the reproductive and pathogenic behaviour of this soilborne pathogen and the onset of a possible speciation within this species.

Dinoflagellates are phytoplanktonic organisms found in both freshwater and marine habitats. They are often studied because related to harmful algal blooms but they are also known to produce bioactive compounds for the treatment of human pathologies. The aim of this study was to sequence the full transcriptome of the dinoflagellate Amphidinium carterae in both nitrogen-starved and -replete culturing conditions (1) to evaluate the response to nitrogen starvation at the transcriptional level, (2) to look for possible polyketide synthases (PKSs) in the studied clone (genes that may be involved in the synthesis of bioactive compounds), (3) if present, to evaluate if nutrient starvation can influence PKS expression, (4) to look for other possible enzymes of biotechnological interest and (5) to test strain cytotoxicity on human cell lines. Results showed an increase in nitrogen metabolism and stress response in nitrogen-starved cells and confirmed the presence of a type I β-ketosynthase. In addition, L-asparaginase (used for the treatment of Leukemia and for acrylamide reduction in food industries) and cellulase (useful for biofuel production and other industrial applications) have been identified for the first time in this species, giving new insights into possible biotechnological applications of dinoflagellates.

Background The environment has a profound influence on the organoleptic quality of tomato (Solanum lycopersicum) fruit, the extent of which depends on a well-regulated and dynamic interplay among genes, metabolites and sensorial attributes. We used a systems biology approach to elucidate the complex interacting mechanisms regulating the plasticity of sensorial traits. To investigate environmentally challenged transcriptomic and metabolomic remodeling and evaluate the organoleptic consequences of such variations we grown three tomato varieties, Heinz 1706, whose genome was sequenced as reference and two ?local? ones, San Marzano and Vesuviano in two different locations of Campania region (Italy). Results Responses to environment were more pronounced in the two ?local? genotypes, rather than in the Heinz 1706. The overall genetic composition of each genotype, acting in trans, modulated the specific response to environment. Duplicated genes and transcription factors, establishing different number of network connections by gaining or losing links, play a dominant role in shaping organoleptic profile. The fundamental role of cell wall metabolism in tuning all the quality attributes, including the sensorial perception, was also highlighted. Conclusions Although similar fruit-related quality processes are activated in the same environment, different tomato genotypes follow distinct transcriptomic, metabolomic and sensorial trajectories depending on their own genetic makeup. Electronic supplementary material The online version of this article (doi:10.1186/s12870-017-1008-4) contains supplementary material, which is available to authorized users.