Comparative transcriptomics of Central Asian Vitis vinifera accessions reveals distinct defense strategies against powdery mildew

The grapevine LysM receptor-like kinase VvLYK5-1 recognizes chitin oligomers through its association with VvLYK1-1

The establishment of defense reactions to protect plants against pathogens requires the recognition of invasion patterns (IPs), mainly detected by plasma membrane-bound pattern recognition receptors (PRRs). Some IPs, also termed elicitors, are used in several biocontrol products that are gradually being developed to reduce the use of chemicals in agriculture. Chitin, the major component of fungal cell walls, as well as its deacetylated derivative, chitosan, are two elicitors known to activate plant defense responses. However, recognition of chitooligosaccharides (COS) in Vitis vinifera is still poorly understood, hampering the improvement and generalization of protection tools for this important crop. In contrast, COS perception in the model plant Arabidopsis thaliana is well described and mainly relies on a tripartite complex formed by the cell surface lysin motif receptor-like kinases (LysM-RLKs) AtLYK1/CERK1, AtLYK4 and AtLYK5, the latter having the strongest affinity for COS. In grapevine, COS perception has for the moment only been demonstrated to rely on two PRRs VvLYK1-1 and VvLYK1-2. Here, we investigated additional players by overexpressing in Arabidopsis the two putative AtLYK5 orthologs from grapevine, VvLYK5-1 and VvLYK5-2. Expression of VvLYK5-1 in the atlyk4/5 double mutant background restored COS sensitivity, such as chitin-induced MAPK activation, defense gene expression, callose deposition and conferred non-host resistance to grapevine downy mildew (Erysiphe necator). Protein-protein interaction studies conducted in planta revealed a chitin oligomer-triggered interaction between VvLYK5-1 and VvLYK1-1. Interestingly, our results also indicate that VvLYK5-1 mediates the perception of chitin but not chitosan oligomers showing a part of its specificity.

The grape powdery mildew resistance loci Ren2, Ren3, Ren4D, Ren4U, Run1, Run1.2b, Run2.1, and Run2.2 activate different transcriptional responses to Erysiphe necator

Multiple grape powdery mildew (PM) genetic resistance (R) loci have been found in wild grape species. Little is known about the defense responses associated with each R locus. In this study, we compare the defense mechanisms associated with PM resistance in interspecific crosses segregating for a single R locus from Muscadinia rotundifolia (Run1, Run1.2b, Run2.1, Run2.2), Vitis cinerea (Ren2), V. romanetii (Ren4D and Ren4U), and the interspecific hybrid Villard blanc (Ren3). By combining optical microscopy, visual scoring, and biomass estimation, we show that the eight R loci confer resistance by limiting infection at different stages. We assessed the defense mechanisms triggered in response to PM at 1 and 5 days post-inoculation (dpi) via RNA sequencing. To account for the genetic differences between species, we developed for each accession a diploid synthetic reference transcriptome by incorporating into the PN40024 reference homozygous and heterozygous sequence variants and de novo assembled transcripts. Most of the R loci exhibited a higher number of differentially expressed genes (DEGs) associated with PM resistance at 1 dpi compared to 5 dpi, suggesting that PM resistance is mostly associated with an early transcriptional reprogramming. Comparison of the PM resistance-associated DEGs showed a limited overlap between pairs of R loci, and nearly half of the DEGs were specific to a single R locus. The largest overlap of PM resistance-associated DEGs was found between Ren3 ⁺, Ren4D ⁺, and Ren4U ⁺ genotypes at 1 dpi, and between Ren4U ⁺ and Run1 ⁺ accessions at 5 dpi. The Ren3 ⁺, Ren4D ⁺, and Ren4U ⁺ were also found to have the highest number of R locus-specific DEGs in response to PM. Both shared and R locus-specific DEGs included genes from different defense-related categories, indicating that the presence of E. necator triggered distinct transcriptional responses in the eight R loci.

Transcriptomic and metabolomic integration as a resource in grapevine to study fruit metabolite quality traits

Transcriptomics and metabolomics are methodologies being increasingly chosen to perform molecular studies in grapevine (Vitis vinifera L.), focusing either on plant and fruit development or on interaction with abiotic or biotic factors. Currently, the integration of these approaches has become of utmost relevance when studying key plant physiological and metabolic processes. The results from these analyses can undoubtedly be incorporated in breeding programs whereby genes associated with better fruit quality (e.g., those enhancing the accumulation of health-promoting compounds) or with stress resistance (e.g., those regulating beneficial responses to environmental transition) can be used as selection markers in crop improvement programs. Despite the vast amount of data being generated, integrative transcriptome/metabolome meta-analyses (i.e., the joint analysis of several studies) have not yet been fully accomplished in this species, mainly due to particular specificities of metabolomic studies, such as differences in data acquisition (i.e., different compounds being investigated), unappropriated and unstandardized metadata, or simply no deposition of data in public repositories. These meta-analyses require a high computational capacity for data mining a priori, but they also need appropriate tools to explore and visualize the integrated results. This perspective article explores the universe of omics studies conducted in V. vinifera, focusing on fruit-transcriptome and metabolome analyses as leading approaches to understand berry physiology, secondary metabolism, and quality. Moreover, we show how omics data can be integrated in a simple format and offered to the research community as a web resource, giving the chance to inspect potential gene-to-gene and gene-to-metabolite relationships that can later be tested in hypothesis-driven research. In the frame of the activities promoted by the COST Action CA17111 INTEGRAPE, we present the first grapevine transcriptomic and metabolomic integrated database (TransMetaDb) developed within the Vitis Visualization (VitViz) platform (https://tomsbiolab.com/vitviz). This tool also enables the user to conduct and explore meta-analyses utilizing different experiments, therefore hopefully motivating the community to generate Findable, Accessible, Interoperable and Reusable (F.A.I.R.) data to be included in the future.

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.

Identification of powdery mildew resistance in wild grapevine (Vitis vinifera subsp. sylvestris Gmel Hegi) from Croatia and Bosnia and Herzegovina

Wild grapevine (Vitis vinifera subsp. sylvestris) is widely recognized as an important source of resistance or tolerance genes for diseases and environmental stresses. Recent studies revealed partial resistance to powdery mildew (Erysiphe necator, PM) in V. sylvestris from Central Asia. Here, we report resistance to PM of V. sylvestris collected from different regions of Croatia and in seedling populations established from in situ V. sylvestris accessions. Ninety-one in situ individuals and 67 V. sylvestris seedlings were evaluated for PM resistance according to OIV 455 descriptor. Three SSR markers (SC47-18, SC8-071-0014, and UDV-124) linked to PM resistance locus Ren1 were used to decipher allelic structure. Nine seedlings showed resistance in in vivo evaluations while leaf disk assays revealed three PM-resistant accessions. One V. vinifera cultivar used as a control for PM evaluations also showed high phenotypic resistance. Based on the presence of one or two resistance alleles that are linked to the Ren1 locus, 32 resistant seedlings and 41 resistant in situ genotypes were identified in the investigated set. Eight seedlings showed consistent phenotypic PM resistance, of which seven carried one or two alleles at the tested markers. This study provides the first evidence of PM resistance present within the eastern Adriatic V. sylvestris germplasm.

Resistance to Powdery Mildew in Qingke Involves the Accumulation of Aromatic Phenolamides Through Jasmonate-Mediated Activation of Defense-Related Genes

Powdery mildew (PM) leads to severe yield reduction in qingke (Hordeum vulgare L. var. nudum). Although studies have focused on identifying PM-related resistance genes, mechanistic insights into the metabolic regulation networks of resistance against PM have rarely been explored in qingke. Here, we integrated transcriptomic, proteomic and metabolomic data using PM-susceptible (G72) and PM-resistant (K69) accessions to systemically explore the mechanisms of PM resistance. The integrated results show that a rapidly transduction of jasmonic acid (JA) and (+)-7-iso-jasmonoyl-L-isoleucine (JA-Ile), and importantly, a inducing accumulation of aromatic PAs conferred qingke-specific resistance for PM stress. Functional analysis revealed that the four BAHD N-acyltransferase genes were responsible for the synthesis of aliphatic and aromatic PAs. The expression of the four genes are induced by methyl jasmonate (MeJA) and PM treatment. Co-expression network analysis shows that a histone lysine demethylase, JMJ705 gene, also induced by MeJA and PM treatment, had highly correlation with PAs biosynthesis. Chromatin immunoprecipitation (ChIP)-seq assays revealed that the level of trimethylated histone H3 lysine 27 (H3K27me3) of the four genes in MeJA and PM-treated plants was significantly reduced. Overall, our results suggest that a novel strategy for jasmonic acid signal-mediated demethylation controlling the accumulation of aromatic PAs to enhance plant immune resistance through removal of H3K27me3 and activating defense-related gene expression.

PRGdb 4.0: an updated database dedicated to genes involved in plant disease resistance process

The Plant Resistance Genes database (PRGdb; http://prgdb.org/prgdb4/) has been greatly expanded, keeping pace with the increasing amount of available knowledge and data (sequenced proteomes, cloned genes, public analysis data, etc.). The easy-to-use style of the database website has been maintained, while an updated prediction tool, more data and a new section have been added. This new section will contain plant resistance transcriptomic experiments, providing additional easy-to-access experimental information. DRAGO3, the tool for automatic annotation and prediction of plant resistance genes behind PRGdb, has been improved in both accuracy and sensitivity, leading to more reliable predictions. PRGdb offers 199 reference resistance genes and 586.652 putative resistance genes from 182 sequenced proteomes. Compared to the previous release, PRGdb 4.0 has increased the number of reference resistance genes from 153 to 199, the number of putative resistance genes from 177K from 76 proteomes to 586K from 182 sequenced proteomes. A new section has been created that collects plant-pathogen transcriptomic data for five species of agricultural interest. Thereby, with these improvements and data expansions, PRGdb 4.0 aims to serve as a reference to the plant scientific community and breeders worldwide, helping to further study plant resistance mechanisms that contribute to fighting pathogens.

Construction of a high-density genetic map and detection of a major QTL of resistance to powdery mildew (Erysiphe necator Sch.) in Caucasian grapes (Vitis vinifera L.)

BACKGROUND

Vitis vinifera L. is the most cultivated grapevine species worldwide. Erysiphe necator Sch., the causal agent of grape powdery mildew, is one of the main pathogens affecting viticulture. V. vinifera has little or no genetic resistances against E. necator and the grape industry is highly dependent on agrochemicals. Some Caucasian V. vinifera accessions have been reported to be resistant to E. necator and to have no genetic relationships to known sources of resistance to powdery mildew. The main purpose of this work was the study and mapping of the resistance to E. necator in the Caucasian grapes 'Shavtsitska' and 'Tskhvedianis tetra'.

RESULTS

The Caucasian varieties 'Shavtsitska' and 'Tskhvedianis tetra' showed a strong partial resistance to E. necator which segregated in two cross populations: the resistant genotypes delayed and limited the pathogen mycelium growth, sporulation intensity and number of conidia generated. A total of 184 seedlings of 'Shavtsitska' x 'Glera' population were genotyped through the Genotyping by Sequencing (GBS) technology and two high-density linkage maps were developed for the cross parents. The QTL analysis revealed a major resistance locus, explaining up to 80.15% of the phenotypic variance, on 'Shavtsitska' linkage group 13, which was associated with a reduced pathogen infection as well as an enhanced plant necrotic response. The genotyping of 105 Caucasian accessions with SSR markers flanking the QTL revealed that the resistant haplotype of 'Shavtsitska' was shared by 'Tskhvedianis tetra' and a total of 25 Caucasian grape varieties, suggesting a widespread presence of this resistance in the surveyed germplasm. The uncovered QTL was mapped in the region where the Ren1 locus of resistance to E. necator, identified in the V. vinifera 'Kishmish vatkana' and related grapes of Central Asia, is located. The genetic analysis conducted revealed that the Caucasian grapes in this study exhibit a resistant haplotype different from that of Central Asian grape accessions.

CONCLUSIONS

The QTL isolated in 'Shavtsitska' and present in the Caucasian V. vinifera varieties could be a new candidate gene of resistance to E. necator to use in breeding programmes. It co-localizes with the Ren1 locus but shows a different haplotype from that of grapevines of Central Asia. We therefore consider that the Caucasian resistance locus, named Ren1.2, contains a member of a cluster of R-genes, of which the region is rich, and to be linked with, or possibly allelic, to Ren1.

Transcriptional, hormonal, and metabolic changes in susceptible grape berries under powdery mildew infection

Grapevine (Vitis vinifera) berries are extremely sensitive to infection by the biotrophic pathogen Erysiphe necator, causing powdery mildew disease with deleterious effects on grape and wine quality. The combined analysis of the transcriptome and metabolome associated with this common fungal infection has not been previously carried out in any fruit. In order to identify the molecular, hormonal, and metabolic mechanisms associated with infection, healthy and naturally infected V. vinifera cv. Carignan berries were collected at two developmental stages: late green (EL33) and early véraison (EL35). RNA sequencing combined with GC-electron impact ionization time-of-flight MS, GC-electron impact ionization/quadrupole MS, and LC-tandem MS analyses revealed that powdery mildew-susceptible grape berries were able to activate defensive mechanisms with the involvement of salicylic acid and jasmonates and to accumulate defense-associated metabolites (e.g. phenylpropanoids, fatty acids). The defensive strategies also indicated organ-specific responses, namely the activation of fatty acid biosynthesis. However, defense responses were not enough to restrict fungal growth. The fungal metabolic program during infection involves secretion of effectors related to effector-triggered susceptibility, carbohydrate-active enzymes and activation of sugar, fatty acid, and nitrogen uptake, and could be under epigenetic regulation. This study also identified potential metabolic biomarkers such as gallic, eicosanoic, and docosanoic acids and resveratrol, which can be used to monitor early stages of infection.

Succession of Fungal Communities at Different Developmental Stages of Cabernet Sauvignon Grapes From an Organic Vineyard in Xinjiang

Fungi present on grape surface considerably impact grape growth and quality. However, information of the fungal community structures and dynamics on the worldwide cash crop, the Cabernet Sauvignon grape, from the budding to ripening stages remains limited. Here, we investigated the succession of fungal communities on Cabernet Sauvignon grapes from an organic vineyard in Xinjiang, China at different developmental stages via high-throughput sequencing combined with multivariate data analysis. In total, 439 fungal amplicon sequence variants (ASVs) from six phyla were identified. The fungal communities differed over the budding to the berry stages. Moreover, Aspergillus, Malassezia, Metschnikowia, and Udeniomyces were predominant during the unripe stage, whereas Erysiphe, Cryptococcus, Vishniacozyma, and Cladosporium were dominant in the ripe stages. Notably, Vishniacozyma was the most abundant genus, conserved in all development stages. Moreover, network analysis resulted in 171 edges—96 negative and 75 positive. Moreover, fungal genera such as Vishniacozyma, Sporobolomyces, Aspergillus, Alternaria, Erysiphe, Toxicodendron, and Metschnikowia were present in the hubs serving as the main connecting nodes. Extensive mutualistic interactions potentially occur among the fungi on the grape surface. In conclusion, the current study expounded the characteristics of the Cabernet Sauvignon grape fungal community during the plant growth process, and the results provided essential insights into the potential impacts of fungal communities on grape growth and health.

Mining Grapevine Downy Mildew Susceptibility Genes: A Resource for Genomics-Based Breeding and Tailored Gene Editing

Several pathogens continuously threaten viticulture worldwide. Until now, the investigation on resistance loci has been the main trend to understand the interaction between grapevine and the mildew causal agents. Dominantly inherited gene-based resistance has shown to be race-specific in some cases, to confer partial immunity, and to be potentially overcome within a few years since its introgression. Recently, on the footprint of research conducted in Arabidopsis, putative genes associated with downy mildew susceptibility have been discovered also in the grapevine genome. In this work, we deep-sequenced four putative susceptibility genes-namely VvDMR6.1, VvDMR6.2, VvDLO1, VvDLO2-in 190 genetically diverse grapevine genotypes to discover new sources of broad-spectrum and recessively inherited resistance. Identified Single Nucleotide Polymorphisms were screened in a bottleneck analysis from the genetic sequence to their impact on protein structure. Fifty-five genotypes showed at least one impacting mutation in one or more of the scouted genes. Haplotypes were inferred for each gene and two of them at the VvDMR6.2 gene were found significantly more represented in downy mildew resistant genotypes. The current results provide a resource for grapevine and plant genetics and could corroborate genomic-assisted breeding programs as well as tailored gene editing approaches for resistance to biotic stresses.

Population structure of Erysiphe necator on domesticated and wild vines in the Middle East raises questions on the origin of the grapevine powdery mildew pathogen

Plant pathogens usually originate and diversify in geographical regions where hosts and pathogens co-evolve. Erysiphe necator, the causal agent of grape powdery mildew, is a destructive pathogen of grapevines worldwide. Although Eastern US is considered the centre of origin and diversity of E. necator, previous reports on resistant native wild and domesticated Asian grapevines suggest Asia as another possible origin of the pathogen. By using multi-locus sequencing, microsatellites and a novel application of amplicon sequencing (AmpSeq), we show that the population of E. necator in Israel is composed of three genetic groups: Groups A and B that are common worldwide, and a new group IL, which is genetically differentiated from any known group in Europe and Eastern US. Group IL showed distinguished ecological characteristics: it was dominant on wild and traditional vines (95%); its abundance increased along the season; and was more aggressive than A and B isolates on both wild and domesticated vines. The low genetic diversity within group IL suggests that it has invaded Israel from another origin. Therefore, we suggest that the Israeli E. necator population was founded by at least two invasions, of which one could be from a non-East American source, possibly from Asian origin.

Identification of SNPs and InDels associated with berry size in table grapes integrating genetic and transcriptomic approaches

BACKGROUND

Berry size is considered as one of the main selection criteria in table grapes breeding programs, due to the consumer preferences. However, berry size is a complex quantitive trait under polygenic control, and its genetic determination of berry weight is not yet fully understood. The aim of this work was to perform marker discovery using a transcriptomic approach, in order to identify and characterize SNP and InDel markers associated with berry size in table grapes. We used an integrative analysis based on RNA-Seq, SNP/InDel search and validation on table grape segregants and varieties with different genetic backgrounds.

RESULTS

Thirty SNPs and eight InDels were identified using a transcriptomic approach (RNA-Seq). These markers were selected from SNP/InDel found among segregants from a Ruby x Sultanina population with contrasting phenotypes for berry size. The set of 38 SNP and InDel markers was distributed in eight chromosomes. Genotype-phenotype association analyses were performed using a set of 13 RxS segregants and 41 table grapes varieties with different genetic backgrounds during three seasons. The results showed several degrees of association of these markers with berry size (10.2 to 30.7%) as other berry-related traits such as length and width. The co-localization of SNP and /or InDel markers and previously reported QTLs and candidate genes associated with berry size were analysed.

CONCLUSIONS

We identified a set of informative and transferable SNP and InDel markers associated with berry size. Our results suggest the suitability of SNPs and InDels as candidate markers for berry weight in seedless table grape breeding. The identification of genomic regions associated with berry weight in chromosomes 8, 15 and 17 was achieved with supporting evidence derived from a transcriptome experiment focused on SNP/InDel search, as well as from a QTL-linkage mapping approach. New regions possibly associated with berry weight in chromosomes 3, 6, 9 and 14 were identified.

VviERF6Ls: an expanded clade in Vitis responds transcriptionally to abiotic and biotic stresses and berry development

Background

VviERF6Ls are an uncharacterized gene clade in Vitis with only distant Arabidopsis orthologs. Preliminary data indicated these transcription factors may play a role in berry development and extreme abiotic stress responses. To better understand this highly duplicated, conserved clade, additional members of the clade were identified in four Vitis genotypes. A meta-data analysis was performed on publicly available microarray and RNA-Seq data (confirmed and expanded with RT-qPCR), and Vitis VviERF6L1 overexpression lines were established and characterized with phenotyping and RNA-Seq.

Results

A total of 18 PN40024 VviERF6Ls were identified; additional VviERF6Ls were identified in Cabernet Sauvignon, Chardonnay, and Carménère. The amino acid sequences of VviERF6Ls were found to be highly conserved. VviERF6L transcripts were detected in numerous plant organs and were differentially expressed in response to numerous abiotic stresses including water deficit, salinity, and cold as well as biotic stresses such as red blotch virus, N. parvum, and E. necator. VviERF6Ls were differentially expressed across stages of berry development, peaking in the pre-veraison/veraison stage and retaining conserved expression patterns across different vineyards, years, and Vitis cultivars. Co-expression network analysis identified a scarecrow-like transcription factor and a calmodulin-like gene with highly similar expression profiles to the VviERF6L clade. Overexpression of VviERF6L1 in a Seyval Blanc background did not result in detectable morphological phenotypes. Genes differentially expressed in response to VviERF6L1 overexpression were associated with abiotic and biotic stress responses.

Conclusions

VviERF6Ls represent a large and distinct clade of ERF transcription factors in grapevine. The high conservation of protein sequence between these 18 transcription factors may indicate these genes originate from a duplication event in Vitis. Despite high sequence similarity and similar expression patterns, VviERF6Ls demonstrate unique levels of expression supported by similar but heterogeneous promoter sequences. VviERF6L gene expression differed between Vitis species, cultivars and organs including roots, leaves and berries. These genes respond to berry development and abiotic and biotic stresses. VviERF6L1 overexpression in Vitis vinifera results in differential expression of genes related to phytohormone and immune system signaling. Further investigation of this interesting gene family is warranted.

Genetic mapping and survey of powdery mildew resistance in the wild Central Asian ancestor of cultivated grapevines in Central Asia

Cultivated grapevines (Vitis vinifera) lack resistance to powdery mildew (PM) with few exceptions. Resistance to this pathogen within V. vinifera has been reported in earlier studies and identified as the Ren1 locus in two Central Asian table grape accessions. Other PM-resistant cultivated varieties and accessions of the wild ancestor V. vinifera subsp. sylvestris were soon identified raising questions regarding the origin of the resistance. In this study, F1 breeding populations were developed with a PM susceptible V. vinifera subsp. vinifera breeding line and a PM-resistant subsp. sylvestris accession. Genotyping was carried out with five Ren1 locus linked SSR markers. A PM resistance locus explaining up to 96% of the phenotypic variation was identified in the same genomic position, where the Ren1 locus was previously reported. New SSR marker alleles linked with the resistance locus were identified. We report results of PM resistance in multiple accessions of subsp. sylvestris collected as seed lots or cuttings from five countries in the Caucasus and Central Asia. A total of 20 females from 11 seed lots and 19 males from nine seed lots collected from Georgia, Armenia, and Azerbaijan were resistant to PM. Three male and one female plant collected as cuttings from Afghanistan and Iran were also resistant to PM. Allelic analysis of markers linked with the Ren1 locus in conjunction with disease evaluation data found a high diversity of allelic haplotypes, which are only possible via recombination events occurring over a long time period. Sequence analysis of two alleles of the SSR marker that cosegregates with the resistance found SNPs that were present in the wild progenitor and in cultivated forms. Variable levels of PM resistance among the tested accessions were also observed. These lines of evidence suggest that the powdery mildew fungus may have been present in Asia for a longer time than currently thought, giving the wild progenitor V. vinifera subsp. sylvestris time to coevolve with and develop resistance to this pathogen.

Advanced Modeling for the Identification of Different Pathogen Tolerant Vines to Reduce Fungicides and Energy Consumption

The aim of this study is the application of advanced modeling techniques to identify powdery mildew tolerant cultivars and reduce fungicides and energy consumption. The energy savings resulting from the increased efficiency of the use of fungicides is an innovative aspect investigated within the project AGROENER researching on energy efficiency. In this preliminary study, investigations through phenotyping methods could represent a potential solution, especially if they are used in combination with tools and algorithms able to extract and convert a large amount of data. Twelve different grapevine cultivars were tested. The construction of an artificial model, characterized by absolute optima of response to a pathogen (i.e., low values of disease incidence and severity and first day of the pathogen appearance), allowed us to cover the potential variability of a real dataset. To identify the cultivars that tolerate powdery mildew the most, two Soft Independent Modeling of Class Analogy (SIMCA) models were built. The modeling efficiencies, indicated by sensitivity value, were equal to 100%. These statistical multivariate classifications identified some of these tolerant cultivars, as the best responding to the pathogen.

Inferring RPW8-NLRs's evolution patterns in seed plants: case study in Vitis vinifera

Genomic and transcriptomic studies in plants and, more in deep, in grapevine reveal that the disease-resistance RNL gene family is highly variable. RNLs (RPW8-NLRs) are a phylogenetically distinct class of nucleotide oligomerization domain (NOD)-like receptors (NLRs) identified in plants. Two RNLs, namely, the NRG1 (N Requirement Gene 1) and the ADR1 (Activated Disease Resistance 1), have been characterized; however, little is known about the RNL evolutionary history in higher plants. To trace the diversification of RNL gene subfamily, we scanned the NLR proteins of 73 plant genomes belonging to 29 taxa, revealing a noticeable diversification across species and within the same genus or botanic family together with a conspicuous expansion in important crop species. To explore the RNL variability in Vitis vinifera and gain information with respect to their structure, evolutionary diversification of five grape genomes ('Aglianico', 'Falanghina', 'Sultanina', 'Tannat', and 'Nebbiolo') has been compared to the reference genome ('Pinot Noir'). The number of RNLs ranged from 6 ('Sultanina') to 14 ('Nebbiolo'), in contrast to the 10 'Pinot Noir' RNLs. The phylogenetic study on grapevine RNLs revealed that all collapsed into NRG1-clade, rather than four. To investigate more in depth the means of intraspecific variability of grape RNL copies, a transcriptomic profiling in response to powdery mildew (PM) infection was carried out through qRT-PCRs and public databases interrogation. The RNL expression variability identified in transcriptome data sets supports the hypothesis of a functional expansion/contraction in grapevine varieties. Although no direct correlations between grapevine PM-resistance and RNL expression was identified, our work can provide good candidates for functional studies able to elucidate the putative "helper" role of RNLs in grape immune signalling.

Emergent Ascomycetes in Viticulture: An Interdisciplinary Overview

The reduction of pesticide usage is a current imperative and the implementation of sustainable viticulture is an urgent necessity. A potential solution, which is being increasingly adopted, is offered by the use of grapevine cultivars resistant to its main pathogenic threats. This, however, has contributed to changes in defense strategies resulting in the occurrence of secondary diseases, which were previously controlled. Concomitantly, the ongoing climate crisis is contributing to destabilizing the increasingly dynamic viticultural context. In this review, we explore the available knowledge on three Ascomycetes which are considered emergent and causal agents of powdery mildew, black rot and anthracnose. We also aim to provide a survey on methods for phenotyping disease symptoms in fields, greenhouse and lab conditions, and for disease control underlying the insurgence of pathogen resistance to fungicide. Thus, we discuss fungal genetic variability, highlighting the usage and development of molecular markers and barcoding, coupled with genome sequencing. Moreover, we extensively report on the current knowledge available on grapevine-ascomycete interactions, as well as the mechanisms developed by the host to counteract the attack. Indeed, to better understand these resistance mechanisms, it is relevant to identify pathogen effectors which are involved in the infection process and how grapevine resistance genes function and impact the downstream cascade. Dealing with such a wealth of information on both pathogens and the host, the horizon is now represented by multidisciplinary approaches, combining traditional and innovative methods of cultivation. This will support the translation from theory to practice, in an attempt to understand biology very deeply and manage the spread of these Ascomycetes.

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.

Genome-wide characterization revealed role of NBS-LRR genes during powdery mildew infection in Vitis vinifera

NBS-LRR comprises a large class of disease resistance (R) proteins that play a widespread role in plant protection against pathogens. In grapevine, powdery mildew cause significant losses in its productivity and efforts are being directed towards finding of resistance loci or genes imparting resistance/tolerance against such fungal diseases. In the present study, we performed genome-wide analysis of NBS-LRR genes during PM infection in grapevine. We identified 18, 23, 12, 16, 10, 10, 9, 20 and 14 differentially expressed NBS-LRR genes in response to PM infection in seven partially PM-resistant (DVIT3351.27, Husseine, Karadzhandal, Khalchili, Late vavilov, O34-16, Sochal) and 2 PM-susceptible (Carignan and Thompson seedless) V. vinifera accessions. Further, the identified sequences were characterized based on chromosomal locations, physicochemical properties, gene structure and motif analysis, and functional annotation by Gene Ontology (GO) mapping. The NBS-LRR genes responsive to powdery mildew could potentially be exploited to improve resistance in grapes.

High-throughput sequencing of microbial community diversity in soil, grapes, leaves, grape juice and wine of grapevine from China

In this study Illumina MiSeq was performed to investigate microbial diversity in soil, leaves, grape, grape juice and wine. A total of 1,043,102 fungal Internal Transcribed Spacer (ITS) reads and 2,422,188 high quality bacterial 16S rDNA sequences were used for taxonomic classification, revealed five fungal and eight bacterial phyla. At the genus level, the dominant fungi were Ascomycota, Sordariales, Tetracladium and Geomyces in soil, Aureobasidium and Pleosporaceae in grapes leaves, Aureobasidium in grape and grape juice. The dominant bacteria were Kaistobacter, Arthrobacter, Skermanella and Sphingomonas in soil, Pseudomonas, Acinetobacter and Kaistobacter in grape and grapes leaves, and Oenococcus in grape juice and wine. Principal coordinate analysis showed structural separation between the composition of fungi and bacteria in all samples. This is the first study to understand microbiome population in soil, grape, grapes leaves, grape juice and wine in Xinjiang through High-throughput Sequencing and identify microorganisms like Saccharomyces cerevisiae and Oenococcus spp. that may contribute to the quality and flavor of wine.

Co-expression network analysis and cis-regulatory element enrichment determine putative functions and regulatory mechanisms of grapevine ATL E3 ubiquitin ligases

Arabidopsis thaliana Toxicos en Levadura (ATL) proteins are a subclass of the RING-H2 zinc finger binding E3 ubiquitin ligases. The grapevine (Vitis vinifera) ATL family was recently characterized, revealing 96 members that are likely to be involved in several physiological processes through protein ubiquitination. However, the final targets and biological functions of most ATL E3 ligases are still unknown. We analyzed the co-expression networks among grapevine ATL genes across a set of transcriptomic data related to defense and abiotic stress, combined with a condition-independent dataset. This revealed strong correlations between ATL proteins and diverse signal transduction components and transcriptional regulators, in particular those involved in immunity. An enrichment analysis of cis-regulatory elements in ATL gene promoters and related co-expressed genes highlighted the importance of hormones in the regulation of ATL gene expression. Our work identified several ATL proteins as candidates for further studies aiming to decipher specific grapevine resistance mechanisms activated in response to pathogens.

Neofusicoccum parvum Colonization of the Grapevine Woody Stem Triggers Asynchronous Host Responses at the Site of Infection and in the Leaves

Grapevine trunk diseases cause important economic losses in vineyards worldwide. Neofusicoccum parvum, one of the most aggressive causal agents of the trunk disease Botryosphaeria dieback, colonizes cells and tissues of the grapevine wood, leading to the formation of an internal canker. Symptoms then extend to distal shoots, with wilting of leaves and bud mortality. Our aim was to characterize the transcriptional dynamics of grapevine genes in the woody stem and in the leaves during Neofusicoccum parvum colonization. Genome-wide transcriptional profiling at seven distinct time points (0, 3, and 24 hours; 2, 6, 8, and 12 weeks) showed that both stems and leaves undergo extensive transcriptomic reprogramming in response to infection of the stem. While most intense transcriptional responses were detected in the stems at 24 hours, strong responses were not detected in the leaves until the next sampling point at 2 weeks post-inoculation. Network co-expression analysis identified modules of co-expressed genes common to both organs and showed most of these genes were asynchronously modulated. The temporal shift between stem vs. leaf responses affected transcriptional modulation of genes involved in both signal perception and transduction, as well as downstream biological processes, including oxidative stress, cell wall rearrangement and cell death. Promoter analysis of the genes asynchronously modulated in stem and leaves during N. parvum colonization suggests that the temporal shift of transcriptional reprogramming between the two organs might be due to asynchronous co-regulation by common transcriptional regulators. Topology analysis of stem and leaf co-expression networks pointed to specific transcription factor-encoding genes, including WRKY and MYB, which may be associated with the observed transcriptional responses in the two organs.

Elucidation of the molecular responses of a cucumber segment substitution line carrying Pm5.1 and its recurrent parent triggered by powdery mildew by comparative transcriptome profiling

BACKGROUND

Powdery mildew (PM) is one of the most severe fungal diseases of cucurbits, but the molecular mechanisms underlying PM resistance in cucumber remain elusive. In this study, we developed a PM resistant segment substitution line SSL508-28 that carried a segment on chromosome five representing the Pm5.1 locus from PM resistant donor Jin5-508 using marker-assisted backcrossing of an elite PM susceptible cucumber inbred line D8.

RESULTS

Whole-genome resequencing of SSL508-28, Jin5-508 and D8 was performed to identify the exact boundaries of the breakpoints for this introgression because of the low density of available single sequence repeat markers. This led to the identification of a ~6.8 Mb substituted segment predicted to contain 856 genes. RNA-seq was used to study gene expression differences in PM treated (plants harvested 48 h after inoculation) and untreated (control) SSL508-28 and D8 lines. Exactly 1,248 and 1,325 differentially expressed genes (DEGs) were identified in SSL508-28 and D8, respectively. Of those, 88 DEGs were located in the ~6.8 Mb segment interval. Based on expression data and annotation, we identified 8 potential candidate genes that may participate in PM resistance afforded by Pm5.1, including two tandemly arrayed genes encoding receptor protein kinases, two transcription factors, two genes encoding remorin proteins, one gene encoding a P-type ATPase and one gene encoding a 70 kDa heat shock protein. The transcriptome data also revealed a complex regulatory network for Pm5.1-mediated PM resistance that may involve multiple signal regulators and transducers, cell wall modifications and the salicylic acid signaling pathway.

CONCLUSION

These findings shed light on the cucumber PM defense mechanisms mediated by Pm5.1 and provided valuable information for the fine mapping of Pm5.1 and breeding of cucumber with enhanced resistance to PM.

Genome-wide characterisation and expression profile of the grapevine ATL ubiquitin ligase family reveal biotic and abiotic stress-responsive and development-related members

The Arabidopsis Tóxicos en Levadura (ATL) protein family is a class of E3 ubiquitin ligases with a characteristic RING-H2 Zn-finger structure that mediates diverse physiological processes and stress responses in plants. We carried out a genome-wide survey of grapevine (Vitis vinifera L.) ATL genes and retrieved 96 sequences containing the canonical ATL RING-H2 domain. We analysed their genomic organisation, gene structure and evolution, protein domains and phylogenetic relationships. Clustering revealed several clades, as already reported in Arabidopsis thaliana and rice (Oryza sativa), with an expanded subgroup of grapevine-specific genes. Most of the grapevine ATL genes lacked introns and were scattered among the 19 chromosomes, with a high level of duplication retention. Expression profiling revealed that some ATL genes are expressed specifically during early or late development and may participate in the juvenile to mature plant transition, whereas others may play a role in pathogen and/or abiotic stress responses, making them key candidates for further functional analysis. Our data offer the first genome-wide overview and annotation of the grapevine ATL family, and provide a basis for investigating the roles of specific family members in grapevine physiology and stress responses, as well as potential biotechnological applications.

Expression of a Grapevine NAC Transcription Factor Gene Is Induced in Response to Powdery Mildew Colonization in Salicylic Acid-Independent Manner

Tissue colonization by grape powdery mildew (PM) pathogen Erysiphe necator (Schw.) Burr triggers a major remodeling of the transcriptome in the susceptible grapevine Vitis vinifera L. While changes in the expression of many genes bear the signature of salicylic acid (SA) mediated regulation, the breadth of PM-induced changes suggests the involvement of additional regulatory networks. To explore PM-associated gene regulation mediated by other SA-independent systems, we designed a microarray experiment to distinguish between transcriptome changes induced by E. necator colonization and those triggered by elevated SA levels. We found that the majority of genes responded to both SA and PM, but certain genes were responsive to PM infection alone. Among them, we identified genes of stilbene synthases, PR-10 proteins, and several transcription factors. The microarray results demonstrated that the regulation of these genes is either independent of SA, or dependent, but SA alone is insufficient to bring about their regulation. We inserted the promoter-reporter fusion of a PM-responsive transcription factor gene into a wild-type and two SA-signaling deficient Arabidopsis lines and challenged the resulting transgenic plants with an Arabidopsis-adapted PM pathogen. Our results provide experimental evidence that this grape gene promoter is activated by the pathogen in a SA-independent manner.

Identification of two novel powdery mildew resistance loci, Ren6 and Ren7, from the wild Chinese grape species Vitis piasezkii

BACKGROUND

Grapevine powdery mildew Erysiphe necator is a major fungal disease in all grape growing countries worldwide. Breeding for resistance to this disease is crucial to avoid extensive fungicide applications that are costly, labor intensive and may have detrimental effects on the environment. In the past decade, Chinese Vitis species have attracted attention from grape breeders because of their strong resistance to powdery mildew and their lack of negative fruit quality attributes that are often present in resistant North American species. In this study, we investigated powdery mildew resistance in multiple accessions of the Chinese species Vitis piasezkii that were collected during the 1980 Sino-American botanical expedition to the western Hubei province of China.

RESULTS

A framework genetic map was developed using simple sequence repeat markers in 277 seedlings of an F1 mapping population arising from a cross of the powdery mildew susceptible Vitis vinifera selection F2-35 and a resistant accession of V. piasezkii DVIT2027. Quantitative trait locus analyses identified two major powdery mildew resistance loci on chromosome 9 (Ren6) and chromosome 19 (Ren7) explaining 74.8 % of the cumulative phenotypic variation. The quantitative trait locus analysis for each locus, in the absence of the other, explained 95.4 % phenotypic variation for Ren6, while Ren7 accounted for 71.9 % of the phenotypic variation. Screening of an additional 259 seedlings of the F1 population and 910 seedlings from four pseudo-backcross populations with SSR markers defined regions of 22 kb and 330 kb for Ren6 and Ren7 in the V. vinifera PN40024 (12X) genome sequence, respectively. Both R loci operate post-penetration through the induction of programmed cell death, but vary significantly in the speed of response and degree of resistance; Ren6 confers complete resistance whereas Ren7 confers partial resistance to the disease with reduced colony size. A comparison of the kinetics of induction of powdery mildew resistance mediated by Ren6, Ren7 and the Run1 locus from Muscadinia rotundifolia, indicated that the speed and strength of resistance conferred by Ren6 is greater than that of Run1 which, in turn, is superior to that conferred by Ren7.

CONCLUSIONS

This is the first report of mapping powdery mildew resistance in the Chinese species V. piasezkii. Two distinct powdery mildew R loci designated Ren6 and Ren7 were found in multiple accessions of this Chinese grape species. Their location on different chromosomes to previously reported powdery mildew resistance R loci offers the potential for grape breeders to combine these R genes with existing powdery mildew R loci to produce grape germplasm with more durable resistance against this rapidly evolving fungal pathogen.

cDNA-AFLP analysis reveals differential gene expression in incompatible interaction between infected non-heading Chinese cabbage and Hyaloperonospora parasitica

Non-heading Chinese cabbage (Brassica rapa ssp. chinensis) is one of the main green leafy vegetables in the world, especially in China, with significant economic value. Hyaloperonospora parasitica is a fungal pathogen responsible for causing downy mildew disease in Chinese cabbage, which greatly affects its production. The objective of this study was to identify transcriptionally regulated genes during incompatible interactions between non-heading Chinese cabbage and H. parasitica using complementary DNA-amplified fragment length polymorphism (cDNA-AFLP). We obtained 129 reliable differential transcript-derived fragments (TDFs) in a resistant line 'Suzhou Qing'. Among them, 121 upregulated TDFs displayed an expression peak at 24-48 h post inoculation (h.p.i.). Fifteen genes were further selected for validation of cDNA-AFLP expression patterns using quantitative reverse transcription PCR. Results confirmed the altered expression patterns of 13 genes (86.7%) revealed by the cDNA-AFLP. We identified four TDFs related to fungal resistance among the 15 TDFs. Furthermore, comparative analysis of four TDFs between resistant line 'Suzhou Qing' and susceptible line 'Aijiao Huang' showed that transcript levels of TDF14 (BcLIK1_A01) peaked at 48 h.p.i. and 25.1-fold increased in the resistant line compared with the susceptible line. Similarly, transcript levels of the other three genes, TDF42 (BcCAT3_A07), TDF75 (BcAAE3_A06) and TDF88 (BcAMT2_A05) peaked at 24, 48 and 24 h.p.i. with 25.1-, 100- and 15.8-fold increases, respectively. The results suggested that the resistance genes tended to transcribe at higher levels in the resistance line than in the susceptible line, which may provide resistance against pathogen infections. The present study might facilitate elucidating the molecular basis of the infection process and identifying candidate genes for resistance improvement of susceptible cultivars.

Transcriptome profiling of grapevine seedless segregants during berry development reveals candidate genes associated with berry weight

BACKGROUND

Berry size is considered as one of the main selection criteria in table grape breeding programs. However, this is a quantitative and polygenic trait, and its genetic determination is still poorly understood. Considering its economic importance, it is relevant to determine its genetic architecture and elucidate the mechanisms involved in its expression. To approach this issue, an RNA-Seq experiment based on Illumina platform was performed (14 libraries), including seedless segregants with contrasting phenotypes for berry weight at fruit setting (FST) and 6-8 mm berries (B68) phenological stages.

RESULTS

A group of 526 differentially expressed (DE) genes were identified, by comparing seedless segregants with contrasting phenotypes for berry weight: 101 genes from the FST stage and 463 from the B68 stage. Also, we integrated differential expression, principal components analysis (PCA), correlations and network co-expression analyses to characterize the transcriptome profiling observed in segregants with contrasting phenotypes for berry weight. After this, 68 DE genes were selected as candidate genes, and seven candidate genes were validated by real time-PCR, confirming their expression profiles.

CONCLUSIONS

We have carried out the first transcriptome analysis focused on table grape seedless segregants with contrasting phenotypes for berry weight. Our findings contributed to the understanding of the mechanisms involved in berry weight determination. Also, this comparative transcriptome profiling revealed candidate genes for berry weight which could be evaluated as selection tools in table grape breeding programs.

Grapevine Pathogenic Microorganisms: Understanding Infection Strategies and Host Response Scenarios

Grapevine (Vitis vinifera L.) is one of the most important fruit crop worldwide. Commercial cultivars are greatly affected by a large number of pathogenic microorganisms that cause diseases during pre- and/or post-harvest periods, affecting production, processing and export, along with fruit quality. Among the potential threats, we can find bacteria, fungi, oomycete, or viruses with different life cycles, infection mechanisms and evasion strategies. While plant-pathogen interactions are cycles of resistance and susceptibility, resistance traits from natural resources are selected and may be used for breeding purposes and for a sustainable agriculture. In this context, here we summarize some of the most important diseases affecting V. vinifera together with their causal agents. The aim of this work is to bring a comprehensive review of the infection strategies deployed by significant types of pathogens while understanding the host response in both resistance and susceptibility scenarios. New approaches being used to uncover grapevine status during biotic stresses and scientific-based procedures needed to control plant diseases and crop protection are also addressed.

Transcriptomic Crosstalk between Fungal Invasive Pathogens and Their Host Cells: Opportunities and Challenges for Next-Generation Sequencing Methods

Fungal invasive infections are an increasing health problem. The intrinsic complexity of pathogenic fungi and the unmet clinical need for new and more effective treatments requires a detailed knowledge of the infection process. During infection, fungal pathogens are able to trigger a specific transcriptional program in their host cells. The detailed knowledge of this transcriptional program will allow for a better understanding of the infection process and consequently will help in the future design of more efficient therapeutic strategies. Simultaneous transcriptomic studies of pathogen and host by high-throughput sequencing (dual RNA-seq) is an unbiased protocol to understand the intricate regulatory networks underlying the infectious process. This protocol is starting to be applied to the study of the interactions between fungal pathogens and their hosts. To date, our knowledge of the molecular basis of infection for fungal pathogens is still very limited, and the putative role of regulatory players such as non-coding RNAs or epigenetic factors remains elusive. The wider application of high-throughput transcriptomics in the near future will help to understand the fungal mechanisms for colonization and survival, as well as to characterize the molecular responses of the host cell against a fungal infection.