A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

Responsive Alternative Splicing Events of Opisthopappus Species against Salt Stress

Salt stress profoundly affects plant growth, prompting intricate molecular responses, such as alternative splicing (AS), for environmental adaptation. However, the response of AS events to salt stress in Opisthopappus (Opisthopappus taihangensis and Opisthopappus longilobus) remains unclear, which is a Taihang Mountain cliff-dwelling species. Using RNA-seq data, differentially expressed genes (DEGs) were identified under time and concentration gradients of salt stress. Two types of AS, skipped exon (SE) and mutually exclusive exons (MXE), were found. Differentially alternative splicing (DAS) genes in both species were significantly enriched in "protein phosphorylation", "starch and sucrose metabolism", and "plant hormone signal transduction" pathways. Meanwhile, distinct GO terms and KEGG pathways of DAS occurred between two species. Only a small subset of DAS genes overlapped with DEGs under salt stress. Although both species likely adopted protein phosphorylation to enhance salt stress tolerance, they exhibited distinct responses. The results indicated that the salt stress mechanisms of both Opisthopappus species exhibited similarities and differences in response to salt stress, which suggested that adaptive divergence might have occurred between them. This study initially provides a comprehensive description of salt responsive AS events in Opisthopappus and conveys some insights into the molecular mechanisms behind species tolerance on the Taihang Mountains.

Boosting grapevine breeding for climate-smart viticulture: from genetic resources to predictive genomics

The multifaceted nature of climate change is increasing the urgency to select resilient grapevine varieties, or generate new, fitter cultivars, to withstand a multitude of new challenging conditions. The attainment of this goal is hindered by the limiting pace of traditional breeding approaches, which require decades to result in new selections. On the other hand, marker-assisted breeding has proved useful when it comes to traits governed by one or few genes with great effects on the phenotype, but its efficacy is still restricted for complex traits controlled by many loci. On these premises, innovative strategies are emerging which could help guide selection, taking advantage of the genetic diversity within the Vitis genus in its entirety. Multiple germplasm collections are also available as a source of genetic material for the introgression of alleles of interest via adapted and pioneering transformation protocols, which present themselves as promising tools for future applications on a notably recalcitrant species such as grapevine. Genome editing intersects both these strategies, not only by being an alternative to obtain focused changes in a relatively rapid way, but also by supporting a fine-tuning of new genotypes developed with other methods. A review on the state of the art concerning the available genetic resources and the possibilities of use of innovative techniques in aid of selection is presented here to support the production of climate-smart grapevine genotypes.

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

Deciphering the role of alternative splicing as modulators of defense response in the MYMIV- Vigna mungo pathosystem

Alternative splicing (AS) is a crucial regulatory mechanism that impacts transcriptome and proteome complexity under stressful situations. Although its role in abiotic stresses is somewhat understood, our understanding of the mechanistic regulation of pre-mRNA splicing in plant-pathogen interaction is meagre. To comprehend this unexplored immune reprogramming mechanism, transcriptome profiles of Mungbean Yellow Mosaic India Virus (MYMIV)-resistant and susceptible Vigna mungo genotypes were analysed for AS genes that may underlie the resistance mechanism. Results revealed a repertoire of AS-isoforms accumulated during pathogenic infestation, with intron retention being the most common AS mechanism. Identification of 688 differential alternatively spliced (DAS) genes in the resistant host elucidates its robust antiviral response, whereas 322 DAS genes were identified in the susceptible host. Enrichment analyses confirmed DAS transcripts pertaining to stress, signalling, and immune system pathways have undergone maximal perturbations. Additionally, a strong regulation of the splicing factors has been observed both at transcriptional and post-transcriptional levels. qPCR validation of candidate DAS transcripts with induced expression upon MYMIV-infection demonstrated a competent immune response in the resistant background. The AS-impacted genes resulted either in partial/complete loss of functional domains or altered sensitivity to miRNA-mediated gene silencing. A complex regulatory module, miR7517-ATAF2, has been identified in an aberrantly spliced ATAF2 isoform that exposes an intronic miR7517 binding site, thereby suppressing the negative regulator to enhance defense reaction. The present study establishes AS as a non-canonical immune reprogramming mechanism that operates in parallel, thereby offering an alternative strategy for developing yellow mosaic-resistant V. mungo cultivars. This article is protected by copyright. All rights reserved.

Characterization of Simple Sequence Repeat (SSR) Markers Mined in Whole Grape Genomes

SSR (simple sequence repeat) DNA markers are widely used for genotype DNA identification, QTL mapping, and analyzing genetic biodiversity. However, SSRs in grapes are still in their early stages, with a few primer pairs accessible. With the whole-genome sequencing (WGS) of several grape varieties, characterization of grape SSR changed to be necessary not only to genomics but to also help SSR development and utility. Based on this, we identified the whole-genome SSR of nine grape cultivars (‘PN40024’, ‘Cabernet Sauvignon’, ‘Carménère’, ‘Chardonnay’, ‘Merlot’, ‘Riesling’, ‘Zinfandel’, ‘Shine Muscat’, and ‘Muscat Hamburg’) with whole-genome sequences released publicly and found that there are great differences in the distribution of SSR loci in different varieties. According to the difference in genome size, the number of SSRs ranged from 267,385 (Cabernet Sauvignon) to 627,429 (Carménère), the density of the SSR locus in the genome of nine cultivars was generally 1 per Kb. SSR motif distribution characteristic analysis of these grape cultivars showed that the distribution patterns among grape cultivars were conservative, mainly enriched in A/T. However, there are some differences in motif types (especially tetranucleotides, pentanucleotides, and hexanucleotides), quantity, total length, and average length in different varieties, which might be related to the size of the assembled genome or the specificity of variety domestication. The distribution characteristics of SSRs were revealed by whole-genome analysis of simple repeats of grape varieties. In this study, 32 pairs of primers with lower polymorphism have been screened, which provided an important research foundation for the development of molecular markers of grape variety identification and the construction of linkage maps of important agronomic traits for crop improvement.

Dissecting the effect of ethylene in the transcriptional regulation of chilling treatment in grapevine leaves

Ethylene (ETH) plays important roles in various development programs and stress responses in plants. In grapevines, ETH increased dramatically under chilling stress and is known to positively regulate cold tolerance. However, the role of ETH in transcriptional regulation during chilling stress of grapevine leaves is still not clear. To address this gap, targeted hormone profiling and transcriptomic analysis were performed on leaves of Vitis amurensis under chilling stress with and without aminoethoxyvinylglycine (AVG, a inhibitor of ETH synthesis) treatment. APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) and WRKY transcription factors (TF) were only the two highly enriched TF families that were consistently up-regulated during chilling stress but inhibited by AVG. The comparison of leaf transcriptomes between chilling treatment and chilling with AVG allowed the identification of potential ETH-regulated genes. Potential genes that are positively regulated by ETH are enriched in solute transport, protein biosynthesis, phytohormone action, antioxidant and carbohydrate metabolism. Conversely, genes related to the synthesis and signaling of ETH, indole-3-acetic acid (IAA), abscisic acid (ABA) were up-regulated by chilling treatment but inhibited by AVG. The contents of ETH, ABA and IAA also paralleled with the transcriptome data, which suggests that the response of ABA and IAA during chilling stress may regulate by ETH signaling, and together may belong to an integrated network of hormonal signaling pathways underpinning chilling stress response in grapevine leaves. Together, these findings provide new clues for further studying the complex regulatory mechanism of ETH under low-temperature stress in plants more generally and new opportunities for breeding cold-resilient grapevines.

The Combination of Both Heat and Water Stresses May Worsen Botryosphaeria Dieback Symptoms in Grapevine

(1) Background: Grapevine trunk diseases (GTDs) have become a global threat to vineyards worldwide. These diseases share three main common features. First, they are caused by multiple pathogenic micro-organisms. Second, these pathogens often maintain a long latent phase, which makes any research in pathology and symptomatology challenging. Third, a consensus is raising to pinpoint combined abiotic stresses as a key factor contributing to disease symptom expression. (2) Methods: We analyzed the impact of combined abiotic stresses in grapevine cuttings artificially infected by two fungi involved in Botryosphaeria dieback (one of the major GTDs), Neofusicoccum parvum and Diplodia seriata. Fungal-infected and control plants were subjected to single or combined abiotic stresses (heat stress, drought stress or both). Disease intensity was monitored thanks to the measurement of necrosis area size. (3) Results and conclusions: Overall, our results suggest that combined stresses might have a stronger impact on disease intensity upon infection by the less virulent pathogen Diplodia seriata. This conclusion is discussed through the impact on plant physiology using metabolomic and transcriptomic analyses of leaves sampled for the different conditions.

Global Landscape of Alternative Splicing in Maize Response to Low Temperature

Maize (Zea mays L.) is an important food crop planted across the world, and low-temperature stress can affect maize germination. Alternative splicing (AS) is widely present in plants under abiotic stress; however, the response of AS to low-temperature stress in maize remains unclear. In this study, a genome-wide analysis of AS during maize response to low temperatures was performed. AS events were distributed on each chromosome, approximately 2.05-2.09 AS events per gene. Seven genes only had AS in low-temperature-resistant inbred lines. A total of 278 KEGGs and 46 GOs were enriched based on overlapping AS genes, which were associated with hormone and oxidoreductase activity. The mutant was used to verify the function of AS gene ZmWRKY48, and the RGR, RSL, RRL, and RRSA of the mutant decreased by 15.16%-19.87% compared with the normal line. These results contribute to subsequent analysis of the regulatory mechanism of maize in response to low-temperature stress.

Grape Berry Responses to Sequential Flooding and Heatwave Events: A Physiological, Transcriptional, and Metabolic Overview

Grapevine cultivation, such as the whole horticulture, is currently challenged by several factors, among which the extreme weather events occurring under the climate change scenario are the most relevant. Within this context, the present study aims at characterizing at the berry level the physiological response of Vitis vinifera cv. Sauvignon Blanc to sequential stresses simulated under a semi-controlled environment: flooding at bud-break followed by multiple summer stress (drought plus heatwave) occurring at pre-vèraison. Transcriptomic and metabolomic assessments were performed through RNASeq and NMR, respectively. A comprehensive hormone profiling was also carried out. Results pointed out a different response to the heatwave in the two situations. Flooding caused a developmental advance, determining a different physiological background in the berry, thus affecting its response to the summer stress at both transcriptional levels, with the upregulation of genes involved in oxidative stress responses, and metabolic level, with the increase in osmoprotectants, such as proline and other amino acids. In conclusion, sequential stress, including a flooding event at bud-break followed by a summer heatwave, may impact phenological development and berry ripening, with possible consequences on berry and wine quality. A berry physiological model is presented that may support the development of sustainable vineyard management solutions to improve the water use efficiency and adaptation capacity of actual viticultural systems to future scenarios.

Identification and analysis of HD-Zip genes involved in the leaf development of Liriodendron chinense using multidimensional analysis

Homeodomain-leucine zipper (HD-Zip) proteins are plant-specific transcription factors that play important roles in different biological processes, especially leaf development. However, no studies to date have identified the HD-Zip genes in Liriodendron chinense nor characterized their functions. We identified the HD-Zip genes in L. chinense by analysing the phylogeny, chromosome location, structure, conserved motif, cis-regulatory elements, synteny, post-transcriptional regulation and expression patterns of these genes during leaf development. A total of 36 LcHD-Zip genes were identified and divided into four subfamilies (HD-Zip I to IV). Synteny analysis revealed that segmental duplication was the main force driving the expansion of LcHD-Zip genes. These 36 LcHD-Zip genes exhibited 11 different expression patterns. Pattern 1, 2, 3, 4, 6, 7, 8 and 9 genes may play important roles in leaf development, such as leaf initiation, leaf polarity establishment, leaf shape development, phytohormone-mediated leaf growth and leaf epidermal structure formation. Four HD-Zip III genes were targeted by microRNAs (miRNAs), and the miR165/166a-HD-Zip regulatory module formed regulated leaf initiation and leaf polarity establishment. Overall, LcHD-Zip genes play key roles in leaf development of L. chinense. This work provides a foundation for the functional verification of HD-Zip genes identified in this study.

Genome-Wide Analysis of Multidrug and Toxic Compound Extruction Transporters in Grape

Grape (Vitis vinifera L.) is an important fruit crop in the world. It is used as a table grape and is also used for raisin and wine production. Grape berries accumulate secondary metabolites, such as anthocyanins, tannins, and resveratrol, which are known as functional compounds for human health. Multidrug and toxic compound extrusion transporter (MATEs) transport secondary metabolites. MATEs also transport other solutes, including organic acids, and toxic xenobiotics, depending on cation gradient and play various roles in plants. MATE comprises 300-500 amino acid residues and possesses a MATE domain and 8-12 transmembrane domains. In the present study, 59 MATE genes were identified in the grape genome, and phylogenetic analysis revealed the presence of four groups of grape MATEs (Group 1-4). Their information, such as gene structures, protein motifs, predicted subcellular localizations, and gene IDs of four genome annotations, that is, CRIBI v1, CRIBI v2, Genoscope, and Vcost v3, were annotated. The transport substrates and physiological functions of grape MATEs were estimated based on their homology with the analyzed MATEs in other plant species. Group 1 may transport toxic compounds and alkaloids, Group 2 may transport polyphenolic compounds, Group 3 may transport organic acids, and Group 4 may transport plant hormones related to signal transduction. In addition to the known anthocyanin transporters, VvMATE37 and VvMATE39, a novel anthocyanin transporter, VvMATE38 in Group 2, was suggested as a key transporter for anthocyanin accumulation in grape berry skin. VvMATE46, VvMATE47, and VvMATE49 in Group 3 may contribute to Al³⁺ detoxification and Fe²⁺/Fe³⁺ translocation via organic acid transport. This study provides helpful and fundamental information for grape MATE studies and resolves the confusion of gene IDs in different genome annotations.

The grapevine aspartic protease gene family: characterization and expression modulation in response to Plasmopara viticola

Grapevine aspartic proteases gene family is characterized and five VviAPs appear to be involved in grapevine defense against downy mildew. Grapevine (Vitis vinifera L.) is one of the most important crops worldwide. However, it is highly susceptible to the downy mildew disease caused by Plasmopara viticola (Berk. & Curt.) Berl. & De Toni. To minimize the use of fungicides used to control P. viticola, it is essential to gain a deeper comprehension on this pathosystem and proteases have gained particular interest in the past decade. Proteases were shown to actively participate in plant-pathogen interactions, not only in the processes that lead to plant cell death, stress responses and protein processing/degradation but also as components of the recognition and signalling pathways. The aim of this study was to identify and characterize the aspartic proteases (APs) involvement in grapevine defense against P. viticola. A genome-wide search and bioinformatics characterization of the V. vinifera AP gene family was conducted and a total of 81 APs proteins, coded by 65 genes, were found. VviAPs proteins can be divided into three categories, similar to those previously described for other plants. Twelve APs coding genes were selected, and expression analysis was conducted at several time-points after inoculation in both compatible and incompatible interactions. Five grapevine APs may be involved in grapevine tolerance against P. viticola. Our findings provide an overall understanding of the VviAPs gene family and establish better groundwork to further describe the roles of VviAPs in defense against P. viticola.

The bud dormancy disconnect: latent buds of grapevine are dormant during summer despite a high metabolic rate

Grapevine (Vitis vinifera L.) displays wide plasticity to climate; however, the physiology of dormancy along a seasonal continuum is poorly understood. Here we investigated the apparent disconnect between dormancy and the underlying respiratory physiology and transcriptome of grapevine buds, from bud set in summer to bud burst in spring. The establishment of dormancy in summer was pronounced and reproducible; however, this was coupled with little or no change in physiology, indicated by respiration, hydration, and tissue oxygen tension. The release of dormancy was biphasic; the depth of dormancy declined substantially by mid-autumn, while the subsequent decline towards spring was moderate. Observed changes in physiology failed to explain the first phase of dormancy decline, in particular. Transcriptome data contrasting development from summer through to spring also indicated that dormancy was poorly reflected by metabolic quiescence during summer and autumn. Gene Ontology and enrichment data revealed the prevailing influence of abscisic acid (ABA)-related gene expression during the transition from summer to autumn, and promoter motif analysis suggested that photoperiod may play an important role in regulating ABA functions during the establishment of dormancy. Transcriptomic data from later transitions reinforced the importance of oxidation and hypoxia as physiological cues to regulate the maintenance of quiescence and resumption of growth. Collectively these data reveal a novel disconnect between growth and metabolic quiescence in grapevine following bud set, which requires further experimentation to explain the phenology and dormancy relationships.

RNA-seq Analysis Reveals Alternative Splicing Under Heat Stress in Rainbow Trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss) is one of the most economically important cold-water farmed species in the world, and transcriptomic studies in response to heat stress have been conducted and will be studied in depth. Alternative splicing (AS), a post-transcriptional regulatory process that regulates gene expression and increases proteomic diversity, is still poorly understood in rainbow trout under heat stress. In the present study, 18,623 alternative splicing events were identified from 9936 genes using RNA transcriptome sequencing technology (RNA-Seq) and genomic information. A total of 2731 differential alternative splicing (DAS) events were found among 2179 differentially expressed genes (DEGs). Gene ontology analysis revealed that the DEGs were mainly enriched in cellular metabolic process, cell part, and organic cyclic compound binding under heat stress. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis displayed that the DEGs were enriched for 39 pathways, and some key pathways, such as lysine degradation, are involved in the regulation of heat stress in liver tissues of rainbow trout. The results were validated by qRT-PCR, confirming reliability of our bioinformatics analysis.

Extent of wild-to-crop interspecific introgression in grapevine (Vitis vinifera) as a consequence of resistance breeding and implications for the crop species definition

Over the past two centuries, introgression through repeated backcrossing has introduced disease resistance from wild grape species into the domesticated lineage Vitis vinifera subsp. sativa. Introgression lines are being cultivated over increasing vineyard surface areas, as their wines now rival in quality those obtained from preexisting varieties. There is, however, a lot of debate about whether and how wine laws defining commercial product categories, which are based on the classification of V. vinifera and interspecific hybrid grapes, should be revised to accommodate novel varieties that do not fit either category. Here, we developed a method of multilocus genotype analysis using short-read resequencing to identify haplotypic blocks of wild ancestry in introgression lines and quantify the physical length of chromosome segments free-of-introgression or with monoallelic and biallelic introgression. We used this genomic data to characterize species, hybrids and introgression lines and show that newly released resistant varieties contain 76.5-94.8% of V. vinifera DNA. We found that varietal wine ratings are not always commensurate with the percentage of V. vinifera ancestry and linkage drag of wild alleles around known resistance genes persists over at least 7.1-11.5 Mb, slowing down the recovery of the recurrent parental genome. This method also allowed us to identify the donor species of known resistance haplotypes, define the ancestry of wild genetic background in introgression lines with complex pedigrees, validate the ancestry of the historic varieties Concord and Norton, and unravel sample curation errors in public databases.

Systematic analysis of the Serine/Arginine-Rich Protein Splicing Factors (SRs) and focus on salt tolerance of PtSC27 in Populus trichocarpa

Serine/Arginine-Rich Protein Splicing Factors (SRs) are indispensable splicing factors, which play significant roles in spliceosome assembly, splicing regulation and regulation of plant stress. However, a comprehensive analysis and function research of SRs in the woody plant is still lacking. In this report, we conducted the identification and comprehensive analysis of the 71 SRs in poplar and three other dicots, including basic characterization, phylogenetic, conserved motifs, gene duplication, promoter and splice isoform of these genes. Based on the publicly available transcriptome data, expression pattern of SRs in poplar under low temperature, high temperature, drought and salt stress were further analyzed. Subsequently, a key candidate gene PtSC27 that responded to salt stress was screened. More importantly, overexpression of PtSC27 increased plant survival rate under salt stress, and enhanced salt tolerance by regulating malondialdehyde (MDA) content, peroxidase (POD) and catalase (CAT) enzyme activities in transgenic plants. Meanwhile, overexpression of PtSC27 made transgenic plants insensitive to exogenous ABA and improved the expression of some ABA signal-related genes under salt stress. Overall, our studies lay a foundation for understanding the structure and function of SRs in the poplar and provide useful gene resources for breeding through genetic engineering.

Multi-dimensional leaf phenotypes reflect root system genotype in grafted grapevine over the growing season

BACKGROUND

Modern biological approaches generate volumes of multi-dimensional data, offering unprecedented opportunities to address biological questions previously beyond reach owing to small or subtle effects. A fundamental question in plant biology is the extent to which below-ground activity in the root system influences above-ground phenotypes expressed in the shoot system. Grafting, an ancient horticultural practice that fuses the root system of one individual (the rootstock) with the shoot system of a second, genetically distinct individual (the scion), is a powerful experimental system to understand below-ground effects on above-ground phenotypes. Previous studies on grafted grapevines have detected rootstock influence on scion phenotypes including physiology and berry chemistry. However, the extent of the rootstock's influence on leaves, the photosynthetic engines of the vine, and how those effects change over the course of a growing season, are still largely unknown.

RESULTS

Here, we investigate associations between rootstock genotype and shoot system phenotypes using 5 multi-dimensional leaf phenotyping modalities measured in a common grafted scion: ionomics, metabolomics, transcriptomics, morphometrics, and physiology. Rootstock influence is ubiquitous but subtle across modalities, with the strongest signature of rootstock observed in the leaf ionome. Moreover, we find that the extent of rootstock influence on scion phenotypes and patterns of phenomic covariation are highly dynamic across the season.

CONCLUSIONS

These findings substantially expand previously identified patterns to demonstrate that rootstock influence on scion phenotypes is complex and dynamic and underscore that broad understanding necessitates volumes of multi-dimensional data previously unmet.

The genomes of 204 Vitis vinifera accessions reveal the origin of European wine grapes

In order to elucidate the still controversial processes that originated European wine grapes from its wild progenitor, here we analyse 204 genomes of Vitis vinifera and show that all analyses support a single domestication event that occurred in Western Asia and was followed by numerous and pervasive introgressions from European wild populations. This admixture generated the so-called international wine grapes that have diffused from Alpine countries worldwide. Across Europe, marked differences in genomic diversity are observed in local varieties that are traditionally cultivated in different wine producing countries, with Italy and France showing the largest diversity. Three genomic regions of reduced genetic diversity are observed, presumably as a consequence of artificial selection. In the lowest diversity region, two candidate genes that gained berry–specific expression in domesticated varieties may contribute to the change in berry size and morphology that makes the fruit attractive for human consumption and adapted for winemaking.

Reports on the origin of European wine grapes are controversial. Here, the authors perform population genetics analyses on a large set of representative wine-making varieties and reveal a single domestication event at the origin of the entire germplasm followed by repeated introgression from wild populations.

Grape Transcriptome Response to Powdery Mildew Infection: Comparative Transcriptome Profiling of Chinese Wild Grapes Provides Insights Into Powdery Mildew Resistance

Erysiphe necator, the fungal pathogen of grape powdery mildew disease, poses a great threat to the grape market and the wine industry. To better understand the molecular basis of grape responses to E. necator, we performed comparative transcriptome profiling on two Chinese wild grape accessions with varying degrees of resistance to E. necator. At 6-, 24-, and 96-h postinoculation of E. necator, 2,856, 2,678, and 1,542 differentially expressed genes (DEGs) were identified in the susceptible accession Vitis pseudoreticulata 'Hunan-1', and at those same time points, 1,921, 2,498, and 3,249 DEGs, respectively, were identified in the resistant accession V. quinquangularis 'Shang-24'. 'Hunan-1' had a substantially larger fraction of down-regulated genes than 'Shang-24' at every infection stage. Analysis of DEGs revealed that up-regulated genes were mostly associated with defense response and disease resistance-related metabolite biosynthesis, and such signaling genes were significantly suppressed in 'Hunan-1'. Interestingly, fatty acid biosynthesis- and elongation-related genes were suppressed by the fungus in the 'Shang-24' accession but somehow induced in the 'Hunan-1' accession, consistent with the concept that E. necator is likely to be a fatty acid auxotroph that requires lipids from the host. Moreover, genes involved in biosynthesis and signaling of phytohormones, such as jasmonic acid and cytokinin, as well as genes encoding protein kinases and nucleotide-binding domain leucine-rich repeat proteins, differentially responded to E. necator in the two wild grapes. The variation of gene regulation associated with nutrient uptake by the fungus and with signaling transduction and pathogen recognition suggests a multilayered regulatory network that works in concert to assist in the establishment of fungal pathogen infections.

Alternative splicing during fruit development among fleshy fruits

BACKGROUND

Alternative splicing (AS) is an important mechanism of posttranscriptional modification and dynamically regulates multiple physiological processes in plants, including fruit ripening. However, little is known about alternative splicing during fruit development in fleshy fruits.

RESULTS

We studied the alternative splicing at the immature and ripe stages during fruit development in cucumber, melon, papaya and peach. We found that 14.96-17.48% of multiexon genes exhibited alternative splicing. Intron retention was not always the most frequent event, indicating that the alternative splicing pattern during different developmental process differs. Alternative splicing was significantly more prevalent at the ripe stage than at the immature stage in cucumber and melon, while the opposite trend was shown in papaya and peach, implying that developmental stages adopt different alternative splicing strategies for their specific functions. Some genes involved in fruit ripening underwent stage-specific alternative splicing, indicating that alternative splicing regulates fruits ripening. Conserved alternative splicing events did not appear to be stage-specific. Clustering fruit developmental stages across the four species based on alternative splicing profiles resulted in species-specific clustering, suggesting that diversification of alternative splicing contributes to lineage-specific evolution in fleshy fruits.

CONCLUSIONS

We obtained high quality transcriptomes and alternative splicing events during fruit development across the four species. Dynamics and nonconserved alternative splicing were discovered. The candidate stage-specific AS genes involved in fruit ripening will provide valuable insight into the roles of alternative splicing during the developmental processes of fleshy fruits.

Alternative splicing regulation appears to play a crucial role in grape berry development and is also potentially involved in adaptation responses to the environment

BACKGROUND

Alternative splicing (AS) produces transcript variants playing potential roles in proteome diversification and gene expression regulation. AS modulation is thus essential to respond to developmental and environmental stimuli. In grapevine, a better understanding of berry development is crucial for implementing breeding and viticultural strategies allowing adaptation to climate changes. Although profound changes in gene transcription have been shown to occur in the course of berry ripening, no detailed study on splicing modifications during this period has been published so far. We report here on the regulation of gene AS in developing berries of two grapevine (Vitis vinifera L.) varieties, Gewurztraminer (Gw) and Riesling (Ri), showing distinctive phenotypic characteristics. Using the software rMATS, the transcriptomes of berries at four developmental steps, from the green stage to mid-ripening, were analysed in pairwise comparisons between stages and varieties.

RESULTS

A total of 305 differential AS (DAS) events, affecting 258 genes, were identified. Interestingly, 22% of these AS events had not been reported before. Among the 80 genes that underwent the most significant variations during ripening, 22 showed a similar splicing profile in Gw and Ri, which suggests their involvement in berry development. Conversely, 23 genes were subjected to splicing regulation in only one variety. In addition, the ratios of alternative isoforms were different in Gw and Ri for 35 other genes, without any change during ripening. This last result indicates substantial AS differences between the two varieties. Remarkably, 8 AS events were specific to one variety, due to the lack of a splice site in the other variety. Furthermore, the transcription rates of the genes affected by stage-dependent splicing regulation were mostly unchanged, identifying AS modulation as an independent way of shaping the transcriptome.

CONCLUSIONS

The analysis of AS profiles in grapevine varieties with contrasting phenotypes revealed some similarity in the regulation of several genes with developmental functions, suggesting their involvement in berry ripening. Additionally, many splicing differences were discovered between the two varieties, that could be linked to phenotypic specificities and distinct adaptive capacities. Together, these findings open perspectives for a better understanding of berry development and for the selection of grapevine genotypes adapted to climate change.

Severe Stunting Symptoms upon Nepovirus Infection Are Reminiscent of a Chronic Hypersensitive-like Response in a Perennial Woody Fruit Crop

Virus infection of plants can result in various degrees of detrimental impacts and disparate symptom types and severities. Although great strides have been made in our understanding of the virus-host interactions in herbaceous model plants, the mechanisms underlying symptom development are poorly understood in perennial fruit crops. Grapevine fanleaf virus (GFLV) causes variable symptoms in most vineyards worldwide. To better understand GFLV-grapevine interactions in relation to symptom development, field and greenhouse trials were conducted with a grapevine genotype that exhibits distinct symptoms in response to a severe and a mild strain of GFLV. After validation of the infection status of the experimental vines by high-throughput sequencing, the transcriptomic and metabolomic profiles in plants infected with the two viral strains were tested and compared by RNA-Seq and LC-MS, respectively, in the differentiating grapevine genotype. In vines infected with the severe GFLV strain, 1023 genes, among which some are implicated in the regulation of the hypersensitive-type response, were specifically deregulated, and a higher accumulation of resveratrol and phytohormones was observed. Interestingly, some experimental vines restricted the virus to the rootstock and remained symptomless. Our results suggest that GFLV induces a strain- and cultivar-specific defense reaction similar to a hypersensitive reaction. This type of defense leads to a severe stunting phenotype in some grapevines, whereas others are resistant. This work is the first evidence of a hypersensitive-like reaction in grapevine during virus infection.

Elevated transcription of transposable elements is accompanied by het-siRNA-driven de novo DNA methylation in grapevine embryogenic callus

Background

Somatic variation is a valuable source of trait diversity in clonally propagated crops. In grapevine, which has been clonally propagated worldwide for centuries, important phenotypes such as white berry colour are the result of genetic changes caused by transposable elements. Additionally, epiallele formation may play a role in determining geo-specific (‘terroir’) differences in grapes and thus ultimately in wine. This genomic plasticity might be co-opted for crop improvement via somatic embryogenesis, but that depends on a species-specific understanding of the epigenetic regulation of transposable element (TE) expression and silencing in these cultures. For this reason, we used whole-genome bisulphite sequencing, mRNA sequencing and small RNA sequencing to study the epigenetic status and expression of TEs in embryogenic callus, in comparison with leaf tissue.

Results

We found that compared with leaf tissue, grapevine embryogenic callus cultures accumulate relatively high genome-wide CHH methylation, particularly across heterochromatic regions. This de novo methylation is associated with an abundance of transcripts from highly replicated TE families, as well as corresponding 24 nt heterochromatic siRNAs. Methylation in the TE-specific CHG context was relatively low over TEs located within genes, and the expression of TE loci within genes was highly correlated with the expression of those genes.

Conclusions

This multi-‘omics analysis of grapevine embryogenic callus in comparison with leaf tissues reveals a high level of genome-wide transcription of TEs accompanied by RNA-dependent DNA methylation of these sequences in trans. This provides insight into the genomic conditions underlying somaclonal variation and epiallele formation in plants regenerated from embryogenic cultures, which is an important consideration when using these tissues for plant propagation and genetic improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07973-9.

A dynamic intron retention program regulates the expression of several hundred genes during pollen meiosis

Intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during early stages of pollen development. To improve our understanding of the gene regulatory mechanisms that drive developmental processes, we performed a genome-wide study of alternative splicing and isoform switching during five key stages of pollen development in field mustard, Brassica rapa. Surprisingly, for several hundred genes (12.3% of the genes analysed), isoform switching results in stage-specific expression of intron-retaining transcripts at the meiotic stage of pollen development. In such cases, we report temporally regulated switching between expression of a canonical, translatable isoform and an intron-retaining transcript that is predicted to produce a truncated and presumably inactive protein. The results suggest a new pervasive mechanism underlying modulation of protein levels in a plant developmental program. The effect is not based on gene expression induction but on the type of transcript produced. We conclude that intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during meiosis, especially genes related to ribosome biogenesis, mRNA transport and nuclear envelope architecture. We also propose that stage-specific expression of a non-functional isoform of Brassica rapa BrSDG8, a non-redundant member of histone methyltransferase gene family, linked to alternative splicing regulation, may contribute to the intron retention observed.

Fatty Acid Desaturases: Uncovering Their Involvement in Grapevine Defence against Downy Mildew

Grapevine downy mildew, caused by the biotrophic oomycete Plasmopara viticola, is one of the most severe and devastating diseases in viticulture. Unravelling the grapevine defence mechanisms is crucial to develop sustainable disease control measures. Here we provide new insights concerning fatty acid's (FA) desaturation, a fundamental process in lipid remodelling and signalling. Previously, we have provided evidence that lipid signalling is essential in the establishment of the incompatible interaction between grapevine and Plasmopara viticola. In the first hours after pathogen challenge, jasmonic acid (JA) accumulation, activation of its biosynthetic pathway and an accumulation of its precursor, the polyunsaturated α-linolenic acid (C18:3), were observed in the leaves of the tolerant genotype, Regent. This work was aimed at a better comprehension of the desaturation processes occurring after inoculation. We characterised, for the first time in Vitis vinifera, the gene family of the FA desaturases and evaluated their involvement in Regent response to Plasmopara viticola. Upon pathogen challenge, an up-regulation of the expression of plastidial FA desaturases genes was observed, resulting in a higher content of polyunsaturated fatty acids (PUFAs) of chloroplast lipids. This study highlights FA desaturases as key players in membrane remodelling and signalling in grapevine defence towards biotrophic pathogens.

Two sides of the same story in grapevine-pathogen interactions

Proteases are an integral part of plant defence systems, and their role in plant-pathogen interactions is unequivocal. Emerging evidence suggests that different protease families contribute to the establishment not only of hypersensitive response, priming, and signalling, but also of recognition events through complex proteolytic cascades. Moreover, they play a crucial role in pathogen/microbe-associated molecular pattern (PAMP/MAMP)-triggered immunity as well as in effector-triggered immunity. However, despite important advances in our understanding of the role of proteases in plant defence, the contribution of proteases to pathogen defence in grapevine remains poorly understood. In this review, we summarize current knowledge of the main grapevine pathosystems and explore the role of serine, cysteine, and aspartic proteases from both the host and pathogen point of views.

Molecular Tools for Adapting Viticulture to Climate Change

Adaptation of viticulture to climate change includes exploration of new geographical areas, new training systems, new management practices, or new varieties, both for rootstocks and scions. Molecular tools can be defined as molecular approaches used to study DNAs, RNAs, and proteins in all living organisms. We present here the current knowledge about molecular tools and their potential usefulness in three aspects of grapevine adaptation to the ongoing climate change. (i) Molecular tools for understanding grapevine response to environmental stresses. A fine description of the regulation of gene expression is a powerful tool to understand the physiological mechanisms set up by the grapevine to respond to abiotic stress such as high temperatures or drought. The current knowledge on gene expression is continuously evolving with increasing evidence of the role of alternative splicing, small RNAs, long non-coding RNAs, DNA methylation, or chromatin activity. (ii) Genetics and genomics of grapevine stress tolerance. The description of the grapevine genome is more and more precise. The genetic variations among genotypes are now revealed with new technologies with the sequencing of very long DNA molecules. High throughput technologies for DNA sequencing also allow now the genetic characterization at the same time of hundreds of genotypes for thousands of points in the genome, which provides unprecedented datasets for genotype-phenotype associations studies. We review the current knowledge on the genetic determinism of traits for the adaptation to climate change. We focus on quantitative trait loci and molecular markers available for developmental stages, tolerance to water stress/water use efficiency, sugar content, acidity, and secondary metabolism of the berries. (iii) Controlling the genome and its expression to allow breeding of better-adapted genotypes. High-density DNA genotyping can be used to select genotypes with specific interesting alleles but genomic selection is also a powerful method able to take into account the genetic information along the whole genome to predict a phenotype. Modern technologies are also able to generate mutations that are possibly interesting for generating new phenotypes but the most promising one is the direct editing of the genome at a precise location.

Functional conservation of the grapevine candidate gene INNER NO OUTER for ovule development and seed formation

Seedlessness represents a highly appreciated trait in table grapes. Based on an interesting case of seedless fruit production described in the crop species Annona squamosa, we focused on the Vitis vinifera INNER NO OUTER (INO) gene as a candidate. This gene encodes a transcription factor belonging to the YABBY family involved in the determination of abaxial identity in several organs. In Arabidopsis thaliana, this gene was shown to be essential for the formation and asymmetric growth of the ovule outer integument and its mutation leads to a phenotypic defect of ovules and failure in seed formation. In this study, we identified in silico the V. vinifera orthologue and investigated its phylogenetic relationship to INO genes from other species and its expression in different organs in seeded and seedless varieties. Applying cross-species complementation, we have tested its functionality in the Arabidopsis ino-1 mutant. We show that the V. vinifera INO successfully rescues the ovule outer integument growth and seeds set and also partially complements the outer integument asymmetric growth in the Arabidopsis mutant, differently from orthologues from other species. These data demonstrate that VviINO retains similar activity and protein targets in grapevine as in Arabidopsis. Potential implications for grapevine breeding are discussed.

Combined Metabolite and Transcriptome Profiling Reveals the Norisoprenoid Responses in Grape Berries to Abscisic Acid and Synthetic Auxin

The abscisic acid (ABA) increase and auxin decline are both indicators of ripening initiation in grape berry, and norisoprenoid accumulation also starts at around the onset of ripening. However, the relationship between ABA, auxin, and norisoprenoids remains largely unknown, especially at the transcriptome level. To investigate the transcriptional and posttranscriptional regulation of the ABA and synthetic auxin 1-naphthaleneacetic acid (NAA) on norisoprenoid production, we performed time-series GC-MS and RNA-seq analyses on Vitis vinifera L. cv. Cabernet Sauvignon grape berries from pre-veraison to ripening. Higher levels of free norisoprenoids were found in ABA-treated mature berries in two consecutive seasons, and both free and total norisoprenoids were significantly increased by NAA in one season. The expression pattern of known norisoprenoid-associated genes in all samples and the up-regulation of specific alternative splicing isoforms of VviDXS and VviCRTISO in NAA-treated berries were predicted to contribute to the norisoprenoid accumulation in ABA and NAA-treated berries. Combined weighted gene co-expression network analysis (WGCNA) and DNA affinity purification sequencing (DAP-seq) analysis suggested that VviGATA26, and the previously identified switch genes of myb RADIALIS (VIT_207s0005g02730) and MAD-box (VIT_213s0158g00100) could be potential regulators of norisoprenoid accumulation. The positive effects of ABA on free norisoprenoids and NAA on total norisoprenoid accumulation were revealed in the commercially ripening berries. Since the endogenous ABA and auxin are sensitive to environmental factors, this finding provides new insights to develop viticultural practices for managing norisoprenoids in vineyards in response to changing climates.

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.

The Grapevine E3 Ubiquitin Ligase VriATL156 Confers Resistance against the Downy Mildew Pathogen Plasmopara viticola

Downy mildew, caused by Plasmopara viticola, is one of the most severe diseases of grapevine (Vitis vinifera L.). Genetic resistance is an effective and sustainable control strategy, but major resistance genes (encoding receptors for specific pathogen effectors) introgressed from wild Vitis species, although effective, may be non-durable because the pathogen can evolve to avoid specific recognition. Previous transcriptomic studies in the resistant species Vitis riparia highlighted the activation of signal transduction components during infection. The transfer of such components to V. vinifera might confer less specific and therefore more durable resistance. Here, we describe the generation of transgenic V. vinifera lines constitutively expressing the V. riparia E3 ubiquitin ligase gene VriATL156. Phenotypic and molecular analysis revealed that the transgenic plants were less susceptible to P. viticola than vector-only controls, confirming the role of this E3 ubiquitin ligase in the innate immune response. Two independent transgenic lines were selected for detailed analysis of the resistance phenotype by RNA-Seq and microscopy, revealing the profound reprogramming of transcription to achieve resistance that operates from the earliest stages of pathogen infection. The introduction of VriATL156 into elite grapevine cultivars could therefore provide an effective and sustainable control measure against downy mildew.

Candidate genes and SNPs associated with stomatal conductance under drought stress in Vitis

BACKGROUND

Understanding the complexity of the vine plant's response to water deficit represents a major challenge for sustainable winegrowing. Regulation of water use requires a coordinated action between scions and rootstocks on which cultivars are generally grafted to cope with phylloxera infestations. In this regard, a genome-wide association study (GWAS) approach was applied on an 'ad hoc' association mapping panel including different Vitis species, in order to dissect the genetic basis of transpiration-related traits and to identify genomic regions of grape rootstocks associated with drought tolerance mechanisms. The panel was genotyped with the GrapeReSeq Illumina 20 K SNP array and SSR markers, and infrared thermography was applied to estimate stomatal conductance values during progressive water deficit.

RESULTS

In the association panel the level of genetic diversity was substantially lower for SNPs loci (0.32) than for SSR (0.87). GWAS detected 24 significant marker-trait associations along the various stages of drought-stress experiment and 13 candidate genes with a feasible role in drought response were identified. Gene expression analysis proved that three of these genes (VIT_13s0019g03040, VIT_17s0000g08960, VIT_18s0001g15390) were actually induced by drought stress. Genetic variation of VIT_17s0000g08960 coding for a raffinose synthase was further investigated by resequencing the gene of 85 individuals since a SNP located in the region (chr17_10,497,222_C_T) was significantly associated with stomatal conductance.

CONCLUSIONS

Our results represent a step forward towards the dissection of genetic basis that modulate the response to water deprivation in grape rootstocks. The knowledge derived from this study may be useful to exploit genotypic and phenotypic diversity in practical applications and to assist further investigations.

RNA-Seq Time Series of Vitis vinifera Bud Development Reveals Correlation of Expression Patterns with the Local Temperature Profile

Plants display sophisticated mechanisms to tolerate challenging environmental conditions and need to manage their ontogenesis in parallel. Here, we set out to generate an RNA-Seq time series dataset throughout grapevine (Vitis vinifera) early bud development. The expression of the developmental regulator VviAP1 served as an indicator of the progression of development. We investigated the impact of changing temperatures on gene expression levels during the time series and detected a correlation between increased temperatures and a high expression level of genes encoding heat-shock proteins. The dataset also allowed the exemplary investigation of expression patterns of genes from three transcription factor (TF) gene families, namely MADS-box, WRKY, and R2R3-MYB genes. Inspection of the expression profiles from all three TF gene families indicated that a switch in the developmental program takes place in July which coincides with increased expression of the bud dormancy marker gene VviDRM1.

Portuguese wild grapevine genome re-sequencing (Vitis vinifera sylvestris)

The first genome of Vitis vinifera vinifera (PN40024), published in 2007, boosted grapevine related studies. While this reference genome is a suitable tool for the overall studies in the field, it lacks the ability to unveil changes accumulated during V. v. vinifera domestication. The subspecies V. v. sylvestris preserves wild characteristics, making it a good material to provide insights into V. v. vinifera domestication. The difference in the reproductive strategy between both subspecies is one of the characteristics that set them apart. While V. v. vinifera flowers are hermaphrodite, V. v. sylvestris is mostly dioecious. In this paper, we compare the re-sequencing of the genomes from a male and a female individual of the wild sylvestris, against the reference vinifera genome (PN40024). Variant analysis reveals a low number but with high impact modifications in coding regions, essentially non-synonymous single nucleotide polymorphisms and frame shifts caused by insertions and deletions. The sex-locus was manually inspected, and the results obtained are in line with the most recent works related with wild grapevine sex. In this paper we also describe for the first time RNA editing in transcripts of 14 genes in the sex-determining region, including VviYABBY and VviPLATZ.

Plastid transit peptides-where do they come from and where do they all belong? Multi-genome and pan-genomic assessment of chloroplast transit peptide evolution

Subcellular relocalization of proteins determines an organism's metabolic repertoire and thereby its survival in unique evolutionary niches. In plants, the plastid and its various morphotypes import a large and varied number of nuclear-encoded proteins to orchestrate vital biochemical reactions in a spatiotemporal context. Recent comparative genomics analysis and high-throughput shotgun proteomics data indicate that there are a large number of plastid-targeted proteins that are either semi-conserved or non-conserved across different lineages. This implies that homologs are differentially targeted across different species, which is feasible only if proteins have gained or lost plastid targeting peptides during evolution. In this study, a broad, multi-genome analysis of 15 phylogenetically diverse genera and in-depth analyses of pangenomes from Arabidopsis and Brachypodium were performed to address the question of how proteins acquire or lose plastid targeting peptides. The analysis revealed that random insertions or deletions were the dominant mechanism by which novel transit peptides are gained by proteins. While gene duplication was not a strict requirement for the acquisition of novel subcellular targeting, 40% of novel plastid-targeted genes were found to be most closely related to a sequence within the same genome, and of these, 30.5% resulted from alternative transcription or translation initiation sites. Interestingly, analysis of the distribution of amino acids in the transit peptides of known and predicted chloroplast-targeted proteins revealed monocot and eudicot-specific preferences in residue distribution.

Molecular memory of Flavescence dorée phytoplasma in recovering grapevines

Flavescence dorée (FD) is a destructive phytoplasma disease of European grapevines. Spontaneous and cultivar-dependent recovery (REC) may occur in the field in FD-infected vines starting the year following the first symptoms. However, the biological underpinnings of this process are still largely unexplored. In this study, transcriptome sequencing (RNAseq), whole-genome bisulphite sequencing (WGBS) and metabolite analysis were combined to dissect molecular and metabolic changes associated to FD and REC in leaf veins collected in the field from healthy (H), FD and REC plants of the highly susceptible Vitis vinifera 'Barbera'. Genes involved in flavonoid biosynthesis, carbohydrate metabolism and stress responses were overexpressed in FD conditions, whereas transcripts linked to hormone and stilbene metabolisms were upregulated in REC vines. Accumulation patterns of abscisic acid and stilbenoid compounds analysed in the same samples confirmed the RNAseq data. In recovery conditions, we also observed the persistence of some FD-induced expression changes concerning inhibition of photosynthetic processes and stress responses. Several differentially expressed genes tied to those pathways also underwent post-transcriptional regulation by microRNAs, as outlined by merging our transcriptomic data set with a previously conducted smallRNAseq analysis. Investigations by WGBS analysis also revealed different DNA methylation marks between REC and H leaves, occurring within the promoters of genes tied to photosynthesis and secondary metabolism. The results allowed us to advance the existence of a "molecular memory" of FDp infection, involving alterations in the DNA methylation status of REC plants potentially related to transcriptional reprogramming events, in turn triggering changes in hormonal and secondary metabolite profiles.

Exploring transcriptional switches from pairwise, temporal and population RNA-Seq data using deepTS

Transcriptional switch (TS) is a widely observed phenomenon caused by changes in the relative expression of transcripts from the same gene, in spatial, temporal or other dimensions. TS has been associated with human diseases, plant development and stress responses. Its investigation is often hampered by a lack of suitable tools allowing comprehensive and flexible TS analysis for high-throughput RNA sequencing (RNA-Seq) data. Here, we present deepTS, a user-friendly web-based implementation that enables a fully interactive, multifunctional identification, visualization and analysis of TS events for large-scale RNA-Seq datasets from pairwise, temporal and population experiments. deepTS offers rich functionality to streamline RNA-Seq-based TS analysis for both model and non-model organisms and for those with or without reference transcriptome. The presented case studies highlight the capabilities of deepTS and demonstrate its potential for the transcriptome-wide TS analysis of pairwise, temporal and population RNA-Seq data. We believe deepTS will help research groups, regardless of their informatics expertise, perform accessible, reproducible and collaborative TS analyses of large-scale RNA-Seq data.

Two-omics data revealed commonalities and differences between Rpv12- and Rpv3-mediated resistance in grapevine

Plasmopara viticola is the causal agent of grapevine downy mildew (DM). DM resistant varieties deploy effector-triggered immunity (ETI) to inhibit pathogen growth, which is activated by major resistance loci, the most common of which are Rpv3 and Rpv12. We previously showed that a quick metabolome response lies behind the ETI conferred by Rpv3 TIR-NB-LRR genes. Here we used a grape variety operating Rpv12-mediated ETI, which is conferred by an independent locus containing CC-NB-LRR genes, to investigate the defence response using GC/MS, UPLC, UHPLC and RNA-Seq analyses. Eighty-eight metabolites showed significantly different concentration and 432 genes showed differential expression between inoculated resistant leaves and controls. Most metabolite changes in sugars, fatty acids and phenols were similar in timing and direction to those observed in Rpv3-mediated ETI but some of them were stronger or more persistent. Activators, elicitors and signal transducers for the formation of reactive oxygen species were early observed in samples undergoing Rpv12-mediated ETI and were paralleled and followed by the upregulation of genes belonging to ontology categories associated with salicylic acid signalling, signal transduction, WRKY transcription factors and synthesis of PR-1, PR-2, PR-5 pathogenesis-related proteins.

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.

Network aggregation improves gene function prediction of grapevine gene co-expression networks

Aggregation across multiple networks highlights robust co-expression interactions and improves the functional connectivity of grapevine gene co-expression networks. In recent years, the rapid accumulation of transcriptome datasets from diverse experimental conditions has enabled the widespread use of gene co-expression network (GCN) analysis in plants. In grapevine, GCN analysis has shown great promise for gene function prediction, however, measurable progress is currently lacking. Using accumulated microarray datasets from the grapevine whole-genome array (33 experiments, 1359 samples), we explored how meta-analysis through aggregation influences the functional connectivity (performance) of derived networks using guilt-by-association neighbor voting. Two annotation schemes, i.e. MapMan BIN and Pfam, at two sparsity thresholds, i.e. top 100 (stringent) and 300 (relaxed) ranked genes were evaluated. We observed that aggregating across multiple networks improves performance dramatically, with the aggregate outperforming the majority of functional terms across individual networks. Network sparsity and size (i.e. the number of samples and aggregates) were key factors influencing performance while the choice of annotation scheme had little. Systematic comparison with various state-of-the-art microarray and RNA-seq networks was also performed, however, none outperformed the aggregate microarray network despite having good predictive performance. Repeating these series of tests using a functional enrichment-based performance metric also showed remarkably consistent findings with guilt-by-association neighbor voting. To demonstrate its functionality, we explore the function and transcriptional regulation of grapevine EXPANSIN genes. We envisage that network aggregation will offer new and unique opportunities for gene function prediction in future grapevine functional genomics studies. To this end, we make the aggregate networks and associated metadata publicly available at VTC-Agg (https://sites.google.com/view/vtc-agg).

Drought stress modulates cuticular wax composition of the grape berry

Drought events are a major challenge for many horticultural crops, including grapes, which are often cultivated in dry and warm climates. It is not understood how the cuticle contributes to the grape berry response to water deficit (WD); furthermore, the cuticular waxes and the related biosynthetic pathways are poorly characterized in this fruit. In this study, we identified candidate wax-related genes from the grapevine genome by phylogenetic and transcriptomic analyses. Developmental and stress response expression patterns of these candidates were characterized across pre-existing RNA sequencing data sets and confirmed a high responsiveness of the pathway to environmental stresses. We then characterized the developmental and WD-induced changes in berry cuticular wax composition, and quantified differences in berry transpiration. Cuticular aliphatic wax content was modulated during development and an increase was observed under WD, with wax esters being strongly up-regulated. These compositional changes were related to up-regulated candidate genes of the aliphatic wax biosynthetic pathway, including CER10, CER2, CER3, CER1, CER4, and WSD1. The effect of WD on berry transpiration was not significant. This study indicates that changes in cuticular wax amount and composition are part of the metabolic response of the grape berry to WD, but these changes do not reduce berry transpiration.

Genome-Scale Characterization of Predicted Plastid-Targeted Proteomes in Higher Plants

Plastids are morphologically and functionally diverse organelles that are dependent on nuclear-encoded, plastid-targeted proteins for all biochemical and regulatory functions. However, how plastid proteomes vary temporally, spatially, and taxonomically has been historically difficult to analyze at a genome-wide scale using experimental methods. A bioinformatics workflow was developed and evaluated using a combination of fast and user-friendly subcellular prediction programs to maximize performance and accuracy for chloroplast transit peptides and demonstrate this technique on the predicted proteomes of 15 sequenced plant genomes. Gene family grouping was then performed in parallel using modified approaches of reciprocal best BLAST hits (RBH) and UCLUST. A total of 628 protein families were found to have conserved plastid targeting across angiosperm species using RBH, and 828 using UCLUST. However, thousands of clusters were also detected where only one species had predicted plastid targeting, most notably in Panicum virgatum which had 1,458 proteins with species-unique targeting. An average of 45% overlap was found in plastid-targeted protein-coding gene families compared with Arabidopsis, but an additional 20% of proteins matched against the full Arabidopsis proteome, indicating a unique evolution of plastid targeting. Neofunctionalization through subcellular relocalization is known to impart novel biological functions but has not been described before on a genome-wide scale for the plastid proteome. Further work to correlate these predicted novel plastid-targeted proteins to transcript abundance and high-throughput proteomics will uncover unique aspects of plastid biology and shed light on how the plastid proteome has evolved to influence plastid morphology and biochemistry.

Discovering the bacteriome of Vitis vinifera cv. Pinot Noir in a conventionally managed vineyard

The structure of the bacteriome associated with grapevine roots can affect plant development, health and grape quality. We previously investigated the bacterial biodiversity of the Vitis vinifera cv. Pinot Noir rhizosphere in a vineyard subjected to integrated pest management. The aim of this work is to characterize the bacteriome of V. vinifera cv. Pinot Noir in a conventionally managed vineyard using a metabarcoding approach. Comparisons between the microbial community structure in bulk soil and rhizosphere (variable space) were performed and shifts of bacteriome according to two sampling times (variable time) were characterized. Bacterial biodiversity was higher at the second than at the first sampling and did not differ according to the variable space. Actinobacteria was the dominant class, with Gaiella as the most represented genus in all the samples. Among Proteobacteria, the most represented classes were Alpha, Beta and Gamma-Proteobacteria, with higher abundance at the second than at the first sampling time. Bradyrhizobium was the most frequent genus among Alpha-Proteobacteria, while Burkholderia was the predominant Beta-Proteobacteria. Among Firmicutes, the frequency of Staphylococcus was higher than 60% in bulk soil and rhizosphere. Finally, the sampling time can be considered as one of the drivers responsible for the bacteriome variations assessed.

Grape-RNA: A Database for the Collection, Evaluation, Treatment, and Data Sharing of Grape RNA-Seq Datasets

Since its inception, RNA sequencing (RNA-seq) has become the most effective way to study gene expression. After more than a decade of development, numerous RNA-seq datasets have been created, and the full utilization of these datasets has emerged as a major issue. In this study, we built a comprehensive database named Grape-RNA, which is focused on the collection, evaluation, treatment, and data sharing of grape RNA-seq datasets. This database contains 1529 RNA-seq samples, 112 microRNA samples from the public platform, and 485 RNA-seq in-house datasets sequenced by our lab. We classified these data into 25 conditions and provide the sample information, cleaned raw data, expression level, assembled unigenes, useful tools, and other relevant information to the users. Thus, this study provides data and tools that should be beneficial for researchers by allowing them to easily use the RNA-seq. The provided information can greatly contribute to grape breeding and genomic and biological research. This study may improve the usage of RNA-seq.

Novel Aspects on The Interaction Between Grapevine and Plasmopara viticola: Dual-RNA-Seq Analysis Highlights Gene Expression Dynamics in The Pathogen and The Plant During The Battle For Infection

Mgaloblishvili, a Vitis vinifera cultivar, exhibits unique resistance traits against Plasmopara viticola, the downy mildew agent. This offers the unique opportunity of exploring the molecular responses in compatible and incompatible plant-pathogen interaction. In this study, whole transcriptomes of Mgaloblishvili, Pinot noir (a V. vinifera susceptible cultivar), and Bianca (a resistant hybrid) leaves, inoculated and non-inoculated with the pathogen, were used to identify P. viticola effector-encoding genes and plant susceptibility/resistance genes. Multiple effector-encoding genes were identified in P. viticola transcriptome, with remarkable expression differences in relation to the inoculated grapevine cultivar. Intriguingly, five apoplastic effectors specifically associated with resistance in V. vinifera. Gene coexpression network analysis identified specific modules and metabolic changes occurring during infection in the three grapevine cultivars. Analysis of these data allowed, for the first time, the detection in V. vinifera of a putative P. viticola susceptibility gene, encoding a LOB domain-containing protein. Finally, the de novo assembly of Mgaloblishvili, Pinot noir, and Bianca transcriptomes and their comparison highlighted novel candidate genes that might be at the basis of the resistant phenotype. These results open the way to functional analysis studies and to new perspectives in molecular breeding of grapevine for resistance to P. viticola.

Alternative Splicing Enhances the Transcriptome Complexity of Liriodendron chinense

Alternative splicing (AS) plays pivotal roles in regulating plant growth and development, flowering, biological rhythms, signal transduction, and stress responses. However, no studies on AS have been performed in Liriodendron chinense, a deciduous tree species that has high economic and ecological value. In this study, we used multiple tools and algorithms to analyze transcriptome data derived from seven tissues via hybrid sequencing. Although only 17.56% (8,503/48,408) of genes in L. chinense were alternatively spliced, these AS genes occurred in 37,844 AS events. Among these events, intron retention was the most frequent AS event, producing 1,656 PTC-containing and 3,310 non-PTC-containing transcripts. Moreover, 183 long noncoding RNAs (lncRNAs) also underwent AS events. Furthermore, weighted gene coexpression network analysis (WGCNA) revealed that there were great differences in the activities of transcription and post-transcriptional regulation between pistils and leaves, and AS had an impact on many physiological and biochemical processes in L. chinense, such as photosynthesis, sphingolipid metabolism, fatty acid biosynthesis and metabolism. Moreover, our analysis showed that the features of genes may affect AS, as AS genes and non-AS genes had differences in the exon/intron length, transcript length, and number of exons/introns. In addition, the structure of AS genes may impact the frequencies and types of AS because AS genes with more exons or introns tended to exhibit more AS events, and shorter introns tended to be retained, whereas shorter exons tended to be skipped. Furthermore, eight AS genes were verified, and the results were consistent with our analysis. Overall, this study reveals that AS and gene interaction are mutual-on one hand, AS can affect gene expression and translation, while on the other hand, the structural characteristics of the gene can also affect AS. This work is the first to comprehensively report on AS in L. chinense, and it can provide a reference for further research on AS in L. chinense.