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

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.

The Aquilegia genome reveals a hybrid origin of core eudicots

Background

Whole-genome duplications (WGDs) have dominated the evolutionary history of plants. One consequence of WGD is a dramatic restructuring of the genome as it undergoes diploidization, a process under which deletions and rearrangements of various sizes scramble the genetic material, leading to a repacking of the genome and eventual return to diploidy. Here, we investigate the history of WGD in the columbine genus Aquilegia, a basal eudicot, and use it to illuminate the origins of the core eudicots.

Results

Within-genome synteny confirms that columbines are ancient tetraploids, and comparison with the grape genome reveals that this tetraploidy appears to be shared with the core eudicots. Thus, the ancient gamma hexaploidy found in all core eudicots must have involved a two-step process: first, tetraploidy in the ancestry of all eudicots, then hexaploidy in the ancestry of core eudicots. Furthermore, the precise pattern of synteny sharing suggests that the latter involved allopolyploidization and that core eudicots thus have a hybrid origin.

Conclusions

Novel analyses of synteny sharing together with the well-preserved structure of the columbine genome reveal that the gamma hexaploidy at the root of core eudicots is likely a result of hybridization between a tetraploid and a diploid species.

Quantitative Trait Locus Analysis of Leaf Morphology Indicates Conserved Shape Loci in Grapevine

Leaf shape in plants plays important roles in water use, canopy structure, and physiological tolerances to abiotic stresses; all important traits for the future development and sustainability of grapevine cultivation. Historically, researchers have used ampelography, the study of leaf shape in grapevines, to differentiate Vitis species and cultivars based on finite leaf attributes. However, ampelographic measurements have limitations and new methods for quantifying shape are now available. We paired an analysis of finite trait attributes with a 17-point landmark survey and generalized Procrustes analysis (GPA) to reconstruct grapevine leaves digitally from five interspecific hybrid mapping families. Using the reconstructed leaves, we performed three types of quantitative trait loci (QTL) analyses to determine the genetic architecture that defines leaf shape. In the first analysis, we compared several important ampelographic measurements as finite trait QTL. In the second and third analyses, we identified significant shape variation via principal components analysis (PCA) and using a multivariate least squares interval mapping (MLSIM) approach. In total, we identified 271 significant QTL across the three measures of leaf shape and identified specific QTL hotspots in the grape genome which appear to drive major aspects of grapevine leaf shape.

Tissue-specific mouse mRNA isoform networks

Alternative Splicing produces multiple mRNA isoforms of genes which have important diverse roles such as regulation of gene expression, human heritable diseases, and response to environmental stresses. However, little has been done to assign functions at the mRNA isoform level. Functional networks, where the interactions are quantified by their probability of being involved in the same biological process are typically generated at the gene level. We use a diverse array of tissue-specific RNA-seq datasets and sequence information to train random forest models that predict the functional networks. Since there is no mRNA isoform-level gold standard, we use single isoform genes co-annotated to Gene Ontology biological process annotations, Kyoto Encyclopedia of Genes and Genomes pathways, BioCyc pathways and protein-protein interactions as functionally related (positive pair). To generate the non-functional pairs (negative pair), we use the Gene Ontology annotations tagged with "NOT" qualifier. We describe 17 Tissue-spEcific mrNa iSoform functIOnal Networks (TENSION) following a leave-one-tissue-out strategy in addition to an organism level reference functional network for mouse. We validate our predictions by comparing its performance with previous methods, randomized positive and negative class labels, updated Gene Ontology annotations, and by literature evidence. We demonstrate the ability of our networks to reveal tissue-specific functional differences of the isoforms of the same genes. All scripts and data from TENSION are available at: https://doi.org/10.25380/iastate.c.4275191 .

Post-transcriptional regulation through alternative splicing after infection with Flavobacterium columnare in channel catfish (Ictalurus punctatus)

Columnaris disease has long been recognized as a serious problem worldwide which affects both wild and cultured freshwater fish including the commercially important channel catfish (Ictalurus punctatus). The fundamental molecular mechanisms of the host immune response to the causative agent Flavobacterium columnare remain unclear, though gene expression analysis after the bacterial infection has been conducted. Alternative splicing, a post-transcriptional regulation process to modulate gene expression and increase the proteomic diversity, has not yet been studied in channel catfish following infection with F. columnare. In this study, genomic information and RNA-Seq datasets of channel catfish were used to characterize the changes of alternative splicing after the infection. Alternative splicing was shown to be induced by F. columnare infection, with 8.0% increase in alternative splicing event at early infection stage. Intriguingly, genes involved in RNA binding and RNA splicing themselves were significantly enriched in differentially alternatively spliced (DAS) gene sets after infection. This finding was consistent with our previous study in channel catfish following infection with Edwardsiella ictaluri. It was suggested to be a universal mechanism that genes involved in RNA binding and splicing were regulated to undergo differential alternative splicing after stresses in channel catfish. Moreover, many immune genes were observed to be differentially alternatively spliced after infection. Further studies need to be performed to get a deeper view of molecular regulation on alternative splicing after stresses, setting a foundation for developing catfish broodstocks with enhanced disease resistance.

Plant and fungus transcriptomic data from grapevine berries undergoing artificially-induced noble rot caused by Botrytis cinerea

Noble rot is a latent infection of grape berries caused by the necrotrophic fungus Botrytis cinerea, which develops under specific climatic conditions. The infected berries undergo biochemical and metabolic changes, associated with a rapid withering, which altogether offer interesting organoleptic features to sweet white wines. In this paper, we provide RNAseq datasets (raw and normalized counts as well as differentially expressed genes lists) of the transcriptome profiles of both grapevine berries (Vitis vinifera cv. Garganega) and B. cinerea during the establishment of noble rot, artificially induced in controlled conditions. The sequencing data are available in the NCBI GEO database under accession number GSE116741. These data were exploited in a comprehensive meta-analysis of gene expression during noble rot infection, gray mold and post-harvest withering. This highlighted an important common transcriptional reprogramming in different botrytized grape berry varieties and led to the identification of key genes specifically modulated during noble rot infection, which are described in the article entitled “Specific molecular interactions between Vitis vinifera and Botrytis cinerea are required for noble rot development in grape berries” Lovato et al., 2019.

Extensive changes in gene expression and alternative splicing due to homoeologous exchange in rice segmental allopolyploids

We report rampant homoeologous exchanges in progenies of a newly synthesized rice segmental allotetraploid and demonstrate their consequences to changes of gene expression and alternative splicing. Allopolyploidization is recurrent across the tree of angiosperms and known as a driving evolutionary force in both plants and animals. A salient feature of allopolyploidization is the induction of homoeologous exchange (HE) events between the constituent subgenomes, which may in turn cause changes in gene expression, transcript alternative splicing, and phenotypic novelty. However, this issue has been poorly studied, largely because lack of a system in which the exact parentage donating the subgenomes is known and the HE events are occurring in real time. Here, we employed whole-genome re-sequencing and RNA-seq-based transcriptome profiling in four randomly chosen progeny individuals (at the 10th-selfed generation) of segmental allotetraploids that were constructed by colchicine-mediated whole-genome doubling of F1 hybrids between the two subspecies (japonica and indica) of Asian cultivated Oryza sativa. We show that rampant HE events occurred in these tetraploid individuals, which converted most of the otherwise heterozygous genomic regions into a homogenized state of one parental subgenome. We demonstrate that genes within these homogenized genomic regions in the tetraploids showed high frequencies of altered expression and enhanced alternative splicing relative to their counterparts in the corresponding diploid parents in the embryo tissue. Intriguingly, limited overlaps between the differentially expressed genes and the differential alternative spliced genes were identified, which were partitioned to distinctly enriched gene ontology terms. Together, our results indicate that HE is a major mechanism to rapidly generate novelty in gene expression and transcriptome diversity, which may facilitate phenotypic innovation in nascent allopolyploids and relevant to allopolyploid crop breeding.

A QTL associated with leaf trichome traits has a major influence on the abundance of the predatory mite Typhlodromus pyri in a hybrid grapevine population

The abundance of predatory phytoseiid mites, Typhlodromus pyri, important biological control agents of spider mite pests in numerous crops, is positively influenced by the density of leaf trichomes and tuft-form domatia in vein axils. Identification of the genetic regions controlling both trophic levels could facilitate the improvement of predatory mite habitat in breeding programs. The abundance of T. pyri and non-glandular trichomes was measured in a segregating F₁ family derived from the cross of the complex Vitis hybrid, 'Horizon', with Illinois 547-1 (V. rupestris B38 × V. cinerea B9), finding positive correlation among traits. High density genetic maps were used to localize one major quantitative trait locus (QTL) on chromosome 1 of Illinois 547-1 associated with both predatory mite abundance and leaf trichomes. This QTL explained 23% of the variation in phytoseiid abundance and similar amounts of variance in domatia rating (21%), domatia size (16%), leaf bristle density (37% in veins and 33% in blades), and leaf hair density (20% in veins and 15% in blades). Another nine QTL distributed among chromosomes 1, 2, 5, 8, and 15 were associated solely with trichome density, and explained 7-17% of the phenotypic variation. Combined, our results provide evidence of the genetic architecture of non-glandular trichomes in Vitis, with a major locus influencing trichome densities, domatia size and predatory mite abundance. This information is relevant for breeding grapevines with a more favorable habitat for biological control agents.

Genome Wide Identification, Molecular Characterization, and Gene Expression Analyses of Grapevine NHX Antiporters Suggest Their Involvement in Growth, Ripening, Seed Dormancy, and Stress Response

Plant NHX antiporters are critical for cellular pH, Na⁺, and K⁺ homeostasis and salt tolerance. Even though their genomic and functional studies have been conducted in many species, the grapevine NHX family has not been described yet. Our work highlights the presence of six VvNHX genes whose phylogenetic analysis revealed their classification in two distinct groups: group I vacuolar (VvNHX1-5) and group II endosomal (VvNHX6). Several cis-acting regulatory elements related to tissue-specific expression, transcription factor binding, abiotic/biotic stresses response, and light regulation elements were identified in their promoter. Expression profile analyses of VvNHX genes showed variable transcription within organs and tissues with diverse patterns according to biochemical, environmental, and biotic treatments. All VvNHXs are involved in berry growth, except VvNHX5 that seems to be rather implicated in seed maturation. VvNHX4 would be more involved in floral development, while VvNHX2 and 3 display redundant roles. QPCR expression analyses of VvNHX1 showed its induction by NaCl and KNO₃ treatments, whereas VvNHX6 was induced by ABA application and strongly repressed by PEG treatment. VvNHX1 plays a crucial role in a bunch of grape developmental steps and adaptation responses through mechanisms of phyto-hormonal signaling. Overall, VvNHX family members could be valuable candidate genes for grapevine improvement.

Characterization of genes and alleles involved in the control of flowering time in grapevine

Grapevine (Vitis vinifera) is one of the most important perennial crop plants in worldwide. Understanding of developmental processes like flowering, which impact quality and quantity of yield in this species is therefore of high interest. This gets even more important when considering some of the expected consequences of climate change. Earlier bud burst and flowering, for example, may result in yield loss due to spring frost. Berry ripening under higher temperatures will impact wine quality. Knowledge of interactions between a genotype or allele combination and the environment can be used for the breeding of genotypes that are better adapted to new climatic conditions. To this end, we have generated a list of more than 500 candidate genes that may play a role in the timing of flowering. The grapevine genome was exploited for flowering time control gene homologs on the basis of functional data from model organisms like A. thaliana. In a previous study, a mapping population derived from early flowering GF.GA-47-42 and late flowering 'Villard Blanc' was analyzed for flowering time QTLs. In a second step we have now established a workflow combining amplicon sequencing and bioinformatics to follow alleles of selected candidate genes in the F1 individuals and the parental genotypes. Allele combinations of these genes in individuals of the mapping population were correlated with early or late flowering phenotypes. Specific allele combinations of flowering time candidate genes within and outside of the QTL regions for flowering time on chromosome 1, 4, 14, 17, and 18 were found to be associated with an early flowering phenotype. In addition, expression of many of the flowering candidate genes was analyzed over consecutive stages of bud and inflorescence development indicating functional roles of these genes in the flowering control network.

Transcriptome profiles of three Muscat table grape cultivars to dissect the mechanism of terpene biosynthesis

Vitis vinifera is widely grown worldwide for making wine and for use as table grapes. Of the existing cultivars, some have a floral and fruity flavour, referred to as a Muscat flavour. It is well-documented that this flavour originates from a series of terpene compounds, but the mechanism of terpene content differences among the Muscat-type cultivars remains unclear. Transcript and terpene metabolite profiles were integrated to elucidate the molecular mechanism of this phenomenon. In this research, three genotypes with different aromatic strengths were investigated by RNA sequencing. A total of 27 fruit samples from three biological replicates were sequenced on Illumina HiSeq2000 at three stages, corresponding to the veraison; berries had intermediate Brix value and were harvest-ripe. After quality assessment and data clearance, a total of 254.18 Gb of data with more than 97% Q20 bases were obtained, approximately 9.41 Gb data were generated per sample. These results will provide a valuable dataset for the discovery of the mechanism of terpene biosynthesis.

Transcriptomics of the grape berry shrivel ripening disorder

The lower expression at veraison of several ripening master regulators "switch genes" can play a central role in the induction of the berry shrivel ripening physiological disorder in grapevine. Berry shrivel (BS) is a ripening physiological disorder affecting grape berry with visible symptoms appearing after veraison. Berry shrivel leads to shrinking berries with a reduced weight and a lower content of sugars and anthocyanins. In this study, for the first time a transcriptomic analysis coupled with selected metabolites quantification was undertaken to understand the metabolic modifications induced by the disorder. Different stages of berry development were considered including pre- and symptomatic berries. No metabolic alterations in the berry transcriptome and in the metabolite content was observed in pre-symptomatic and pre-veraison samples. Interestingly, at veraison, with still not visible symptoms appearing on the berry, a subset of genes, called switch genes previously suggested as master regulators of the ripening onset in grape berries, were strongly lower expressed in BS. Later during the ripening phase and with visible symptoms of the disorder, more than 3000 genes were differentially expressed. The genes up-regulated were related to hormone biosynthesis, response to stress and the phenylpropanoid pathway, while the genes down-regulated during ripening belonged mainly to the flavonoid pathway, and the sugar metabolism. In agreement, BS berries showed lower content of sugars and anthocyanins from the onset of veraison onward, while the amount of acids was not significantly affected. In conclusion, these results highlight a pivotal role of the switch genes in grapevine ripening, as well as their possible contribution to induce the ripening disorder berry shrivel, although it remains unclear whether this is part of the cause or consequences of the BS disorder.

Characterization and Cloning of Grape Circular RNAs Identified the Cold Resistance-Related Vv-circATS1

Circular RNAs (circRNAs) are widely distributed and play essential roles in a series of developmental processes, although none have been identified or characterized in grapevines (Vitis vinifera). In this study, we characterized the function of grape circRNA and uncovered thousands of putative back-splicing sites by global transcriptome analysis. Our results indicated that several reported circRNA prediction algorithms should be used simultaneously to obtain comprehensive and reliable circRNA predictions in plants. Furthermore, the length of introns flanking grape circRNAs was closely related to exon circularization. Although the longer introns flanking grape circRNAs appeared to circularize more efficiently, a 20- to 50-nt region seemed large enough to drive grape circRNA biogenesis. In addition, the endogenous introns flanking circularized exon(s) in conjunction with reverse complementary sequences could support the accurate and efficient circularization of various exons in grape, which constitutes a new tool for exploring the functional consequences caused by circRNA expression. Finally, we identified 475 differentially expressed circRNAs in grape leaves under cold stress. Overexpression of Vv-circATS1, a circRNA derived from glycerol-3-P acyltransferase, improved cold tolerance in Arabidopsis (Arabidopsis thaliana), while the linear RNA derived from the same sequence cannot. These results indicate the functional difference between circRNA and linear RNA, and provide new insight into plant abiotic stress resistance.

Many rice genes are differentially spliced between roots and shoots but cytokinin has minimal effect on splicing

Alternatively spliced genes produce multiple spliced isoforms, called transcript variants. In differential alternative splicing, transcript variant abundance differs across sample types. Differential alternative splicing is common in animal systems and influences cellular development in many processes, but its extent and significance is not as well known in plants. To investigate differential alternative splicing in plants, we examined RNA-Seq data from rice seedlings. The data included three biological replicates per sample type, approximately 30 million sequence alignments per replicate, and four sample types: roots and shoots treated with exogenous cytokinin delivered hydroponically or a mock treatment. Cytokinin treatment triggered expression changes in thousands of genes but had negligible effect on splicing patterns. However, many genes were differentially spliced between mock-treated roots and shoots, indicating that our methods were sufficiently sensitive to detect differential splicing between data sets. Quantitative fragment analysis of reverse transcriptase-PCR products made from newly prepared rice samples confirmed 9 of 10 differential splicing events between rice roots and shoots. Differential alternative splicing typically changed the relative abundance of splice variants that co-occurred in a data set. Analysis of a similar (but less deeply sequenced) RNA-Seq data set from Arabidopsis showed the same pattern. In both the Arabidopsis and rice RNA-Seq data sets, most genes annotated as alternatively spliced had small minor variant frequencies. Of splicing choices with abundant support for minor forms, most alternative splicing events were located within the protein-coding sequence and maintained the annotated reading frame. A tool for visualizing protein annotations in the context of genomic sequence (ProtAnnot) together with a genome browser (Integrated Genome Browser) were used to visualize and assess effects of differential splicing on gene function. In general, differentially spliced regions coincided with conserved protein domains, indicating that differential alternative splicing is likely to affect protein function between root and shoot tissue in rice.

Comparative transcriptomics analysis uncovers alternative splicing events and molecular markers in cabbage (Brassica oleracea L.)

Alternative splicing (AS) events were identified and verified in cabbage by comparative transcriptome analysis. The corresponding markers were developed and the germplasm resources were identified. Alternative splicing (AS) is a central regulatory mechanism that greatly contributes to plant gene expression and transcriptome diversity. A large body of evidence has shown that AS complexity is relevant for plant development, evolution, complexity, and adaptation. Both insertion/deletion (InDel) and single nucleotide polymorphism (SNP) are typically co-dominant inheritance markers and have abundant polymorphisms. These have been widely used for marker-assisted selection, genetic mapping, and germplasm identification in plants. However, little is known about the molecular mechanisms underlying AS events and the development of markers including SNP and InDel from the cabbage transcriptome. In this study, three cabbage transcriptome datasets were collected and aligned to the cabbage reference genome to analyze AS events and marker development. 31,524 AS events were identified from three cabbage genotypes, accounting for 20.8% of the total cabbage genes. Alternative 3' splice site donor (A3SS) was the most frequent type of the four main AS events in cabbage. 70,475 InDels and 706,269 SNPs were identified with average frequencies of 1 InDel/6.9 kb and 1 SNP/0.7 kb, respectively. 71,942 potential SSRs were identified in 53,129 assembled unigenes with a density of 1 SSR/6.8 kb. The ratio of SNPs with synonymous/non-synonymous mutations was 1:0.65. 142 InDels and 36 SNPs were randomly selected and validated via Sanger sequencing and polymorphism was found among 66.2% of the InDels and 78.6% of the SNPs. Furthermore, 35 informative InDel markers were successfully used for genetic diversity analysis on 36 cabbage accessions. These results facilitate understanding of the molecular regulation mechanism underlying AS events in cabbage. They also provide molecular marker resource data for genetic mapping construction and germplasm identification, and facilitate the genetic improvement of cabbage via breeding.

Alternative splicing coupled to nonsense-mediated mRNA decay contributes to the high-altitude adaptation of maca (Lepidium meyenii)

Alpine plants remain the least studied plant communities in terrestrial ecosystems. However, how they adapt to high-altitude environments is far from clear. Here, we used RNA-seq to investigate a typical alpine plant maca (Lepidium meyenii) to understand its high-altitude adaptation at transcriptional and post-transcriptional level. At transcriptional level, we found that maca root significantly up-regulated plant immunity genes in day-time comparing to night-time, and up-regulated abiotic (cold/osmotic) stress response genes in Nov and Dec comparing to Oct. In addition, 17 positively selected genes were identified, which could be involved in mitochondrion. At post-transcriptional level, we found that maca had species-specific characterized alternative splicing (AS) profile which could be influenced by stress environments. For example, the alternative 3' splice site events (A3SS, 39.62%) were predominate AS events in maca, rather than intron retention (IR, 23.17%). Interestingly, besides serine/arginine-rich (SR) proteins and long non-coding RNAs (lncRNAs), a lot of components in nonsense-mediated mRNA decay (NMD) were identified under differential alternative splicing (DAS), supporting AS coupled to NMD as essential mechanisms for maca's stress responses and high-altitude adaptation. Taken together, we first attempted to unveil maca's high-altitude adaptation mechanisms based on transcriptome and post-transcriptome evidence. Our data provided valuable insights to understand the high-altitude adaptation of alpine plants.

Spaceflight-induced alternative splicing during seedling development in Arabidopsis thaliana

Plants grown in spaceflight experience novel environmental signals, including those associated with microgravity and ionizing radiation. Spaceflight triggers a response involving transcriptional re-programming and altered cell morphology, though many aspects of this response remain uncharacterized. We analyzed the spaceflight-induced transcriptome with a focus on genes that undergo alternative splicing to examine differential splicing associated with spaceflight—an unstudied characteristic of the molecular response to spaceflight exposure. RNA sequence data obtained during the APEX03 spaceflight experiment that was collected from two Arabidopsis thaliana ecotypes at two seedling stages grown onboard the International Space Station, or as ground controls at Kennedy Space Center, were re-examined to detect alternative splicing differences induced by spaceflight. Presence/absence variation analysis was used to identify putative expression-level differences in alternatively spliced isoforms between spaceflight and ground controls and was followed by analysis of significant differential alternative splicing. This study provides the first evidence of a role for alternative splicing in the molecular processes of physiological adaptation to the spaceflight environment.

Genome-wide identification and comparative analysis of alternative splicing across four legume species

MAIN CONCLUSION

Alternative splicing EVENTS were genome-wide identified for four legume species, and nitrogen fixation-related gene families and evolutionary analysis was also performed. Alternative splicing (AS) is a key regulatory mechanism that contributes to transcriptome and proteome diversity. Investigation of the genome-wide conserved AS events across different species will help with the understanding of the evolution of the functional diversity in legumes, allowing for genetic improvement. Genome-wide identification and characterization of AS were performed using the publically available mRNA, EST, and RNA-Seq data for four important legume species. A total of 15,165 AS genes in Glycine max, 6077 in Cicer arietinum, 7240 in Medicago truncatula, and 7358 in Lotus japonicus were identified. Intron retention (IntronR) was the dominant AS type among the identified events, with IntronR occurring from 53.76% in M. truncatula to 43.91% in C. arietinum. We identified 1159 AS genes that were conserved among four species. Furthermore, nine nitrogen fixation-related gene families with 237 genes were identified, and 80 of them were AS, accounting for the 43.48% in G. max and 27.78% in C. arietinum. An evolutionary analysis showed that these AS genes tended to be located adjacent to each other in the evolutionary tree and are unbalanced in the distribution in the sub-family. This study provides a foundation for future studies on transcription complexity, evolution, and the role of AS on plant functional regulation.

Flooding Responses on Grapevine: A Physiological, Transcriptional, and Metabolic Perspective

Studies on model plants have shown that temporary soil flooding exposes roots to a significant hypoxic stress resulting in metabolic re-programming, accumulation of toxic metabolites and hormonal imbalance. To date, physiological and transcriptional responses to flooding in grapevine are poorly characterized. To fill this gap, we aimed to gain insights into the transcriptional and metabolic changes induced by flooding on grapevine roots (K5BB rootstocks), on which cv Sauvignon blanc (Vitis vinifera L.) plants were grafted. A preliminary experiment under hydroponic conditions enabled the identification of transiently and steadily regulated hypoxia-responsive marker genes and drafting a model for response to oxygen deprivation in grapevine roots. Afterward, over two consecutive vegetative seasons, flooding was imposed to potted vines during the late dormancy period, to mimick the most frequent waterlogging events occurring in the field. Untargeted transcriptomic and metabolic profiling approaches were applied to investigate early responses of grapevine roots during exposure to hypoxia and subsequent recovery after stress removal. The initial hypoxic response was marked by a significant increase of the hypoxia-inducible metabolites ethanol, GABA, succinic acid and alanine which remained high also 1 week after recovery from flooding with the exception of ethanol that leveled off. Transcriptomic data supported the metabolic changes by indicating a substantial rearrangement of primary metabolic pathways through enhancement of the glycolytic and fermentative enzymes and of a subset of enzymes involved in the TCA cycle. GO and KEGG pathway analyses of differentially expressed genes showed a general down-regulation of brassinosteroid, auxin and gibberellin biosynthesis in waterlogged plants, suggesting a general inhibition of root growth and lateral expansion. During recovery, transcriptional activation of gibberellin biosynthetic genes and down-regulation of the metabolic ones may support a role for gibberellins in signaling grapevine rootstocks waterlogging metabolic and hormonal changes to the above ground plant. The significant internode elongation measured upon budbreak during recovery in plants that had experienced flooding supported this hypothesis. Overall integration of these data enabled us to draft a first comprehensive view of the molecular and metabolic pathways involved in grapevine's root responses highlighting a deep metabolic and transcriptomic reprogramming during and after exposure to waterlogging.

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.

The Rpv3-3 Haplotype and Stilbenoid Induction Mediate Downy Mildew Resistance in a Grapevine Interspecific Population

The development of new resistant varieties to the oomycete Plasmopara viticola (Berk.& Curt) is a promising way to combat downy mildew (DM), one of the major diseases threatening the cultivated grapevine (Vitis vinifera L.). Taking advantage of a segregating population derived from "Merzling" (a mid-resistant hybrid) and "Teroldego" (a susceptible landrace), 136 F1 individuals were characterized by combining genetic, phenotypic, and gene expression data to elucidate the genetic basis of DM resistance and polyphenol biosynthesis upon P. viticola infection. An improved consensus linkage map was obtained by scoring 192 microsatellite markers. The progeny were screened for DM resistance and production of 42 polyphenols. QTL mapping showed that DM resistance is associated with the herein named Rpv3-3 specific haplotype and it identified 46 novel metabolic QTLs linked to 30 phenolics-related parameters. A list of the 95 most relevant candidate genes was generated by specifically exploring the stilbenoid-associated QTLs. Expression analysis of 11 genes in Rpv3-3 +/- genotypes displaying disparity in DM resistance level and stilbenoid accumulation revealed significant new candidates for the genetic control of stilbenoid biosynthesis and oligomerization. These overall findings emphasized that DM resistance is likely mediated by the major Rpv3-3 haplotype and stilbenoid induction.

Metabolomic and transcriptomic changes underlying cold and anaerobic stresses after storage of table grapes

The currently accepted paradigm is that fruits and vegetables should be consumed fresh and that their quality deteriorates during storage; however, there are indications that some metabolic properties can, in fact, be improved. We examined the effects of low temperature and high-CO₂ conditions on table grapes, Vitis vinifera L. cv. 'Superior Seedless'. Berries were sampled at harvest (T0) and after low-temperature storage for 6 weeks under either normal atmosphere conditions (TC) or under an O₂ level of 5 kPa and elevated CO₂ levels of 5, 10 or 15 kPa (T5, T10, T15). Accumulation of 10 stilbenes, including E-ε-viniferin, E-miyabenol C and piceatannol, significantly increased under TC treatment as compared to T0 or T15. Sensory analysis demonstrated that elevated CO₂ elicited dose-dependent off-flavor accumulation. These changes were accompanied by an accumulation of 12 volatile metabolites, e.g., ethyl acetate and diacetyl, that imparted disagreeable flavors to fresh fruit. Transcriptome analysis revealed enrichment of genes involved in pyruvate metabolism and the phenylpropanoid pathway. One of the transcription factors induced at low temperature but not under high CO₂ was VvMYB14, which regulates stilbene biosynthesis. Our findings reveal the potential to alter the levels of targeted metabolites in stored produce through understanding the effects of postharvest treatments.

Swift metabolite changes and leaf shedding are milestones in the acclimation process of grapevine under prolonged water stress

BACKGROUND

Grape leaves provide the biochemical substrates for berry development. Thus, understanding the regulation of grapevine leaf metabolism can aid in discerning processes fundamental to fruit development and berry quality. Here, the temporal alterations in leaf metabolism in Merlot grapevine grown under sufficient irrigation and water deficit were monitored from veraison until harvest.

RESULTS

The vines mediated water stress gradually and involving multiple strategies: osmotic adjustment, transcript-metabolite alteration and leaf shedding. Initially stomatal conductance and leaf water potential showed a steep decrease together with the induction of stress related metabolism, e.g. up-regulation of proline and GABA metabolism and stress related sugars, and the down-regulation of developmental processes. Later, progressive soil drying was associated with an incremental contribution of Ca²⁺ and sucrose to the osmotic adjustment concomitant with the initiation of leaf shedding. Last, towards harvest under progressive stress conditions following leaf shedding, incremental changes in leaf water potential were measured, while the magnitude of perturbation in leaf metabolism lessened.

CONCLUSIONS

The data present evidence that over time grapevine acclimation to water stress diversifies in temporal responses encompassing the alteration of central metabolism and gene expression, osmotic adjustments and reduction in leaf area. Together these processes mitigate leaf water stress and aid in maintaining the berry-ripening program.

Expanding Alternative Splicing Identification by Integrating Multiple Sources of Transcription Data in Tomato

Tomato (Solanum lycopersicum) is an important vegetable and fruit crop. Its genome was completely sequenced and there are also a large amount of available expressed sequence tags (ESTs) and short reads generated by RNA sequencing (RNA-seq) technologies. Mapping transcripts including mRNA sequences, ESTs, and RNA-seq reads to the genome allows identifying pre-mRNA alternative splicing (AS), a post-transcriptional process generating two or more RNA isoforms from one pre-mRNA transcript. We comprehensively analyzed the AS landscape in tomato by integrating genome mapping information of all available mRNA and ESTs with mapping information of RNA-seq reads which were collected from 27 published projects. A total of 369,911 AS events were identified from 34,419 genomic loci involving 161,913 transcripts. Within the basic AS events, intron retention is the prevalent type (18.9%), followed by alternative acceptor site (12.9%) and alternative donor site (7.3%), with exon skipping as the least type (6.0%). Complex AS types having two or more basic event accounted for 54.9% of total AS events. Within 35,768 annotated protein-coding gene models, 23,233 gene models were found having pre-mRNAs generating AS isoform transcripts. Thus the estimated AS rate was 65.0% in tomato. The list of identified AS genes with their corresponding transcript isoforms serves as a catalog for further detailed examination of gene functions in tomato biology. The post-transcriptional information is also expected to be useful in improving the predicted gene models in tomato. The sequence and annotation information can be accessed at plant alternative splicing database (http://proteomics.ysu.edu/altsplice).

Increased Alternative Splicing as a Host Response to Edwardsiella ictaluri Infection in Catfish

Alternative splicing is the process of generating multiple transcripts from a single pre-mRNA used by eukaryotes to regulate gene expression and increase proteomic complexity. Although alternative splicing profiles have been well studied in mammalian species, they have not been well studied in aquatic species, especially after biotic stresses. In the present study, genomic information and RNA-Seq datasets were utilized to characterize alternative splicing profiles and their induced changes after bacterial infection with Edwardsiella ictaluri in channel catfish (Ictalurus punctatus). A total of 27,476 alternative splicing events, derived from 9694 genes, were identified in channel catfish. Exon skipping was the most abundant while mutually exclusive exon was the least abundant type of alternative splicing. Alternative splicing was greatly induced by E. ictaluri infection with 21.9% increase in alternative splicing events. Interestingly, genes involved in RNA binding and RNA splicing themselves were significantly enriched in differentially alternatively spliced genes after infection. Sequence analyses of splice variants of a representative alternatively spliced gene, splicing factor srsf2, revealed that certain spliced transcripts may undergo nonsense-mediated decay (NMD), suggesting functional significance of the induced alternative splicing. Although statistical analysis was not possible with such large datasets, results from quantitative real-time PCR from representative differential alternative splicing events provided general validation of the bacterial infection-induced alternative splicing. This is the first comprehensive study of alternative splicing and its changes in response to bacterial infection in fish species, providing insights into the molecular mechanisms of host responses to biotic stresses.

The interplay between membrane lipids and phospholipase A family members in grapevine resistance against Plasmopara viticola

Grapevine downy mildew, caused by the biotrophic oomycete Plasmopara viticola, is one of the most important diseases in modern viticulture. The search for sustainable disease control measure is of extreme importance, thus becoming imperative to fully characterize the mechanisms leading to an incompatible interaction. We have previously shown that lipid signalling events play an important role in grapevine's response to this pathogen, namely through changes in linolenic acid content, lipid peroxidation and jasmonic acid synthesis. Here, we have characterized the modulation of lipid metabolism in leaves from two V. vinifera cultivars (resistant and susceptible to P. viticola) in the first hours after pathogen inoculation. Prior to pathogen inoculation both genotypes present an inherently different fatty acid composition that is highly modulated in the resistant genotype after pathogen challenge. Such changes involve modulation of phospholipase A activity suggesting that the source of lipids mobilized upon pathogen infection are the chloroplast membranes. This work thus provides original evidence on the involvement of lipid signalling and phospholipases in grapevine immune responses to pathogen infection. The results are discussed considering the implications on the plant's physiological status and the use of discriminating lipid/fatty acids pattern in future selection procedures of cultivars.

Concealed ester formation and amino acid metabolism to volatile compounds in table grape (Vitis vinifera L.) berries

Volatile esters contribute to the aroma and flavor of many fruits but are normally absent in grape berries (Vitis vinifera L.). To examine the biosynthetic potential of grape berries to form volatile esters, berry sections were incubated with exogenous L-Phe, L-Leu or L-Met. In general, amino-acid incubation caused the accumulation of the respective aldehydes and alcohols. Moreover, L-Leu incubation resulted in the accumulation of 3-methylbutyl acetate and L-Phe incubation resulted in the accumulation 2-phenylethyl acetate in 'Muscat Hamburg' but not in the other grape accessions. Exogenous L-Met administration did not result in volatile esters accumulation but the accumulation of sulfur volatile compounds such as methional and dimethyl disulfide was prominent. Berry-derived cell-free extracts displayed differential alcohol acetyltransferase activities and supported the formation of 3-methylbutyl acetate and benzyl acetate. 2-Phenylethyl acetate was produced only in 'Muscat Hamburg' cell-free extracts. VvAAT2, a newly characterized gene, was preferentially expressed in 'Muscat Hamburg' berries and functionally expressed in E. coli. VvAAT2 possesses alcohol acetyltransferase activity utilizing benzyl alcohol, 2-phenylethanol, hexanol or 3-methylbutanol as substrates. Our study demonstrates that grape berries have a concealed potential to accumulate volatile esters and this process is limited by substrate availability.

Unique resistance traits against downy mildew from the center of origin of grapevine (Vitis vinifera)

The Eurasian grapevine (Vitis vinifera), an Old World species now cultivated worldwide for high-quality wine production, is extremely susceptible to the agent of downy mildew, Plasmopara viticola. The cultivation of resistant V. vinifera varieties would be a sustainable way to reduce the damage caused by the pathogen and the impact of disease management, which involves the economic, health and environmental costs of frequent fungicide application. We report the finding of unique downy mildew resistance traits in a winemaking cultivar from the domestication center of V. vinifera, and characterize the expression of a range of genes associated with the resistance mechanism. Based on comparative experimental inoculations, confocal microscopy and transcriptomics analyses, our study shows that V. vinifera cv. Mgaloblishvili, native to Georgia (South Caucasus), exhibits unique resistance traits against P. viticola. Its defense response, leading to a limitation of P. viticola growth and sporulation, is determined by the overexpression of genes related to pathogen recognition, the ethylene signaling pathway, synthesis of antimicrobial compounds and enzymes, and the development of structural barriers. The unique resistant traits found in Mgaloblishvili highlight the presence of a rare defense system in V. vinifera against P. viticola which promises fresh opportunities for grapevine genetic improvement.

Identification, characterization, and expression analysis of calmodulin and calmodulin-like genes in grapevine (Vitis vinifera) reveal likely roles in stress responses

Calcium (Ca²⁺) is an ubiquitous key second messenger in plants, where it modulates many developmental and adaptive processes in response to various stimuli. Several proteins containing Ca²⁺ binding domain have been identified in plants, including calmodulin (CaM) and calmodulin-like (CML) proteins, which play critical roles in translating Ca²⁺ signals into proper cellular responses. In this work, a genome-wide analysis conducted in Vitis vinifera identified three CaM- and 62 CML-encoding genes. We assigned gene family nomenclature, analyzed gene structure, chromosomal location and gene duplication, as well as protein motif organization. The phylogenetic clustering revealed a total of eight subgroups, including one unique clade of VviCaMs distinct from VviCMLs. VviCaMs were found to contain four EF-hand motifs whereas VviCML proteins have one to five. Most of grapevine CML genes were intronless, while VviCaMs were intron rich. All the genes were well spread among the 19 grapevine chromosomes and displayed a high level of duplication. The expression profiling of VviCaM/VviCML genes revealed a broad expression pattern across all grape organs and tissues at various developmental stages, and a significant modulation in biotic stress-related responses. Our results highlight the complexity of CaM/CML protein family also in grapevine, supporting the versatile role of its different members in modulating cellular responses to various stimuli, in particular to biotic stresses. This work lays the foundation for further functional and structural studies on specific grapevine CaMs/CMLs in order to better understand the role of Ca²⁺-binding proteins in grapevine and to explore their potential for further biotechnological applications.

Transcriptome analysis of alternative splicing in peanut (Arachis hypogaea L.)

BACKGROUND

Alternative splicing (AS) represents a mechanism widely used by eukaryotes for the post-transcriptional regulation of genes. The detailed exploration of AS in peanut has not been documented.

RESULTS

The strand-specific RNA-Seq technique was exploited to characterize the distribution of AS in the four samples of peanut (FH1-seed1, FH1-seed2, FH1-root and FH1-leaf). AS was detected as affecting around 37.2% of the full set of multi-exon genes. Some of these genes experienced AS throughout the plant, while in the case of others, the effect was organ-specific. Overall, AS was more frequent in the seed than in either the root or leaf. The predominant form of AS was intron retention, and AS in transcription start site and transcription terminal site were commonly identified in all the four samples. It is interesting that in genes affected by AS, the majority experienced only a single type of event. Not all of the in silico predicted transcripts appeared to be translated, implying that these are either degraded or sequestered away from the translation machinery. With respect to genes involved in fatty acid metabolism, about 61.6% were shown to experience AS.

CONCLUSION

Our report contributes significantly in AS analysis of peanut genes in general, and these results have not been mentioned before. The specific functions of different AS forms need further investigation.

The Major Origin of Seedless Grapes Is Associated with a Missense Mutation in the MADS-Box Gene VviAGL11

Seedlessness is greatly prized by consumers of fresh grapes. While stenospermocarpic seed abortion determined by the SEED DEVELOPMENT INHIBITOR (SDI) locus is the usual source of seedlessness in commercial grapevine (Vitis vinifera) cultivars, the underlying sdi mutation remains unknown. Here, we undertook an integrative approach to identify the causal mutation. Quantitative genetics and fine-mapping in two 'Crimson Seedless'-derived F1 mapping populations confirmed the major effect of the SDI locus and delimited the sdi mutation to a 323-kb region on chromosome 18. RNA-sequencing comparing seed traces of seedless and seeds of seeded F1 individuals identified processes triggered during sdi-determined seed abortion, including the activation of salicylic acid-dependent autoimmunity. The RNA-sequencing data set was investigated for candidate genes, and while no evidence for causal cis-acting regulatory mutations was detected, deleterious nucleotide changes in coding sequences of the seedless haplotype were predicted in two genes within the sdi fine-mapping interval. Targeted resequencing of the two genes in a collection of 124 grapevine cultivars showed that only the point variation causing the arginine-197-to-leucine substitution in the seed morphogenesis regulator gene AGAMOUS-LIKE11 (VviAGL11) was fully linked with stenospermocarpy. The concurrent postzygotic variation identified for this missense polymorphism and seedlessness phenotype in seeded somatic variants of the original stenospermocarpic cultivar supports a causal effect. We postulate that seed abortion caused by this amino acid substitution in VviAGL11 is the major cause of seedlessness in cultivated grapevine. This information can be exploited to boost seedless grape breeding.

Genomic signatures of different adaptations to environmental stimuli between wild and cultivated Vitis vinifera L

The application of population genetic methods in combination with gene mapping strategies can help to identify genes and mutations selected during the evolution from wild plants to crops and to explore the considerable genetic variation still maintained in natural populations. We genotyped a grapevine germplasm collection of 44 wild (Vitis vinifera subsp. sylvestris) and 48 cultivated (V. vinifera subsp. sativa) accessions at 54 K single-nucleotide polymorphisms (SNPs) to perform a whole-genome comparison of the main population genetic statistics. The analysis of Wright Fixation Index (FST) along the whole genome allowed us to identify several putative "signatures of selection" spanning over two thousand SNPs significantly differentiated between sativa and sylvestris. Many of these genomic regions included genes involved in the adaptation to environmental changes. An overall reduction of nucleotide diversity was observed across the whole genome within sylvestris, supporting a small effective population size of the wild grapevine. Tajima's D resulted positive in both wild and cultivated subgroups, which may indicate an ongoing balancing selection. Association mapping for six domestication-related traits was performed in combination with population genetics, providing further evidence of different perception and response to environmental stresses between sativa and sylvestris.

Annotation, classification, genomic organization and expression of the Vitis vinifera CYPome

Cytochromes P450 are enzymes that participate in a wide range of functions in plants, from hormonal signaling and biosynthesis of structural polymers, to defense or communication with other organisms. They represent one of the largest gene/protein families in the plant kingdom. The manual annotation of cytochrome P450 genes in the genome of Vitis vinifera PN40024 revealed 579 P450 sequences, including 279 complete genes. Most of the P450 sequences in grapevine genome are organized in physical clusters, resulting from tandem or segmental duplications. Although most of these clusters are small (2 to 35, median = 3), some P450 families, such as CYP76 and CYP82, underwent multiple duplications and form large clusters of homologous sequences. Analysis of gene expression revealed highly specific expression patterns, which are often the same within the genes in large physical clusters. Some of these genes are induced upon biotic stress, which points to their role in plant defense, whereas others are specifically activated during grape berry ripening and might be responsible for the production of berry-specific metabolites, such as aroma compounds. Our work provides an exhaustive and robust annotation including clear identification, structural organization, evolutionary dynamics and expression patterns for the grapevine cytochrome P450 families, paving the way to efficient functional characterization of genes involved in grapevine defense pathways and aroma biosynthesis.

In Vino Veritas: An Invitation for Ambitious, Collaborative Proteogenomics Campaigns on Plant and Animal Models

Vitis vinifera has been an emblematic plant for humans since the Neolithic period. Human civilization has been shaped by its domestication as both its medicinal and nutritional values were exploited. It is now cultivated on all habitable continents, and more than 5000 varieties have been developed. A global passion for the art of wine fuels innovation and a profound desire for knowledge on this plant. The genome sequence of a homozygotic cultivar and several RNA-seq datasets on other varieties have been released paving the way to gaining further insight into its biology and tailoring improvements to varieties. However, its genome annotation remains unpolished. In this issue of Proteomics, Chapman and Bellgard (Proteomics 2017, 17, 1700197) discuss how proteogenomics can help improve genome annotation. By mining shotgun proteomics data, they defined new protein-coding genes, refined gene structures, and corrected numerous mRNA splicing events. This stimulating study shows how large international consortia could work together to improve plant and animal genome annotation on a large scale. To achieve this aim, time should be invested to generate comprehensive, high-quality experimental datasets for a wide range of well-defined lineages and exploit them with pipelines capable of handling giant datasets.

Characterization of the ABA Receptor VlPYL1 That Regulates Anthocyanin Accumulation in Grape Berry Skin

ABA plays a crucial role in controlling several ripening-associated processes in grape berries. The soluble proteins named as PYR (pyrabactin resistant)/PYL (PYR-like)/RCAR (regulatory component of ABA receptor) family have been characterized as ABA receptors. Here, the function of a grape PYL1 encoding gene involved in the response to ABA was verified through heterologous expression. The expression level of VlPYL1 was highest in grape leaf and fruit tissues of the cultivar Kyoho, and the expression of VlPYL1 was increased during fruit development and showed a reduction in ripe berries. Over-expression of VlPYL1 enhances ABA sensitivity in Arabidopsis. Using the transient overexpression technique, the VlPYL1 gene was over-expressed in grape berries. Up-regulation of the VlPYL1 gene not only promoted anthocyanin accumulation but also induced a set of ABA-responsive gene transcripts, including ABF2 and BG3. Although tobacco rattle virus (TRV)-induced gene silencing (VIGS) was not successfully applied in the "Kyoho" grape, the application of the transient overexpression technique in grape fruit could be used as a novel tool for studying grape fruit development.

Single and multiple phenotype QTL analyses of downy mildew resistance in interspecific grapevines

KEY MESSAGE

Downy mildew resistance across days post-inoculation, experiments, and years in two interspecific grapevine F1 families was investigated using linear mixed models and Bayesian networks, and five new QTL were identified. Breeding grapevines for downy mildew disease resistance has traditionally relied on qualitative gene resistance, which can be overcome by pathogen evolution. Analyzing two interspecific F1 families, both having ancestry derived from Vitis vinifera and wild North American Vitis species, across 2 years and multiple experiments, we found multiple loci associated with downy mildew sporulation and hypersensitive response in both families using a single phenotype model. The loci explained between 7 and 17% of the variance for either phenotype, suggesting a complex genetic architecture for these traits in the two families studied. For two loci, we used RNA-Seq to detect differentially transcribed genes and found that the candidate genes at these loci were likely not NBS-LRR genes. Additionally, using a multiple phenotype Bayesian network analysis, we found effects between the leaf trichome density, hypersensitive response, and sporulation phenotypes. Moderate-high heritabilities were found for all three phenotypes, suggesting that selection for downy mildew resistance is an achievable goal by breeding for either physical- or non-physical-based resistance mechanisms, with the combination of the two possibly providing durable resistance.

Combinatorial Regulation of Stilbene Synthase Genes by WRKY and MYB Transcription Factors in Grapevine (Vitis vinifera L.)

Stilbene synthase (STS) is the key enzyme leading to the biosynthesis of resveratrol. Recently we reported two R2R3-MYB transcription factor (TF) genes that regulate the stilbene biosynthetic pathway in grapevine: VviMYB14 and VviMYB15. These genes are strongly co-expressed with STS genes under a range of stress and developmental conditions, in agreement with the specific activation of STS promoters by these TFs. Genome-wide gene co-expression analysis using two separate transcriptome compendia based on microarray and RNA sequencing data revealed that WRKY TFs were the top TF family correlated with STS genes. On the basis of correlation frequency, four WRKY genes, namely VviWRKY03, VviWRKY24, VviWRKY43 and VviWRKY53, were further shortlisted and functionally validated. Expression analyses under both unstressed and stressed conditions, together with promoter-luciferase reporter assays, suggested different hierarchies for these TFs in the regulation of the stilbene biosynthetic pathway. In particular, VviWRKY24 seems to act as a singular effector in the activation of the VviSTS29 promoter, while VviWRKY03 acts through a combinatorial effect with VviMYB14, suggesting that these two regulators may interact at the protein level as previously reported in other species.

Identification and characterization of evolutionarily conserved alternative splicing events in a mangrove genus Sonneratia

Alternative splicing (AS), which produces multiple mRNA transcripts from a single gene, plays crucial roles in plant growth, development and environmental stress responses. Functional significances of conserved AS events among congeneric species have not been well characterized. In this study, we performed transcriptome sequencing to characterize AS events in four common species of Sonneratia, a mangrove genus excellently adaptive to intertidal zones. 7,248 to 12,623 AS events were identified in approximately 25% to 35% expressed genes in the roots of the four species. The frequency of AS events in Sonneratia was associated with genomic features, including gene expression level and intron/exon number and length. Among the four species, 1,355 evolutionarily conserved AS (ECAS) events were identified from 1,170 genes. Compared with non-ECAS events, ECAS events are of shorter length and less possibility to introduce premature stop codons (PTCs) and frameshifts. Functional annotations of the genes containing ECAS events showed that four of the 26 enriched Gene Ontology (GO) terms are involved in proton transport, signal transduction and carbon metabolism, and 60 genes from another three GO terms are implicated in responses to osmotic, oxidative and heat stresses, which may contribute to the adaptation of Sonneratia species to harsh intertidal environments.

Capturing the Alternative Cleavage and Polyadenylation Sites of 14 NAC Genes in Populus Using a Combination of 3'-RACE and High-Throughput Sequencing

Detection of complex splice sites (SSs) and polyadenylation sites (PASs) of eukaryotic genes is essential for the elucidation of gene regulatory mechanisms. Transcriptome-wide studies using high-throughput sequencing (HTS) have revealed prevalent alternative splicing (AS) and alternative polyadenylation (APA) in plants. However, small-scale and high-depth HTS aimed at detecting genes or gene families are very few and limited. We explored a convenient and flexible method for profiling SSs and PASs, which combines rapid amplification of 3′-cDNA ends (3′-RACE) and HTS. Fourteen NAC (NAM, ATAF1/2, CUC2) transcription factor genes of Populus trichocarpa were analyzed by 3′-RACE-seq. Based on experimental reproducibility, boundary sequence analysis and reverse transcription PCR (RT-PCR) verification, only canonical SSs were considered to be authentic. Based on stringent criteria, candidate PASs without any internal priming features were chosen as authentic PASs and assumed to be PAS-rich markers. Thirty-four novel canonical SSs, six intronic/internal exons and thirty 3′-UTR PAS-rich markers were revealed by 3′-RACE-seq. Using 3′-RACE and real-time PCR, we confirmed that three APA transcripts ending in/around PAS-rich markers were differentially regulated in response to plant hormones. Our results indicate that 3′-RACE-seq is a robust and cost-effective method to discover SSs and label active regions subjected to APA for genes or gene families. The method is suitable for small-scale AS and APA research in the initial stage.

Intron retention in viruses and cellular genes: Detention, border controls and passports

Intron retention (IR), where one or more introns remain in the RNA after splicing, was long thought to be rare in mammalian cells, albeit common in plants and some viruses. Largely due to the development of better methods for RNA analysis, it has now been recognized that IR is much more common than previously thought and that this mechanism is likely to play an important role in mammalian gene regulation. To date, most publications and reviews about IR have described the resulting mRNAs as "dead end" products, with no direct consequence for the proteome. However, there are also many reports of mRNAs with retained introns giving rise to alternative protein isoforms. Although this was originally revealed in viral systems, there are now numerous examples of bona fide cellular proteins that are translated from mRNAs with retained introns. These new isoforms have sometimes been shown to have important regulatory functions. In this review, we highlight recent developments in this area and the research on viruses that led the way to the realization of the many ways in which mRNAs with retained introns can be regulated. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing RNA Export and Localization > Nuclear Export/Import RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Genome-wide sequencing of longan (Dimocarpus longan Lour.) provides insights into molecular basis of its polyphenol-rich characteristics

Longan (Dimocarpus longan Lour.), an important subtropical fruit in the family Sapindaceae, is grown in more than 10 countries. Longan is an edible drupe fruit and a source of traditional medicine with polyphenol-rich traits. Tree size, alternate bearing, and witches' broom disease still pose serious problems. To gain insights into the genomic basis of longan traits, a draft genome sequence was assembled. The draft genome (about 471.88 Mb) of a Chinese longan cultivar, "Honghezi," was estimated to contain 31 007 genes and 261.88 Mb of repetitive sequences. No recent whole-genome-wide duplication event was detected in the genome. Whole-genome resequencing and analysis of 13 cultivated D. longan accessions revealed the extent of genetic diversity. Comparative transcriptome studies combined with genome-wide analysis revealed polyphenol-rich and pathogen resistance characteristics. Genes involved in secondary metabolism, especially those from significantly expanded (DHS, SDH, F3΄H, ANR, and UFGT) and contracted (PAL, CHS, and F3΄5΄H) gene families with tissue-specific expression, may be important contributors to the high accumulation levels of polyphenolic compounds observed in longan fruit. The high number of genes encoding nucleotide-binding site leucine-rich repeat (NBS-LRR) and leucine-rich repeat receptor-like kinase proteins, as well as the recent expansion and contraction of the NBS-LRR family, suggested a genomic basis for resistance to insects, fungus, and bacteria in this fruit tree. These data provide insights into the evolution and diversity of the longan genome. The comparative genomic and transcriptome analyses provided information about longan-specific traits, particularly genes involved in its polyphenol-rich and pathogen resistance characteristics.

Global profiling of alternative splicing landscape responsive to drought, heat and their combination in wheat (Triticum aestivum L.)

Plant can acquire tolerance to environmental stresses via transcriptome reprogramming at transcriptional and alternative splicing (AS) levels. However, how AS coordinates with transcriptional regulation to contribute to abiotic stresses responses is still ambiguous. In this study, we performed genome-wide analyses of AS responses to drought stress (DS), heat stress (HS) and their combination (HD) in wheat seedlings, and further compared them with transcriptional responses. In total, we found 200, 3576 and 4056 genes exhibiting significant AS pattern changes in response to DS, HS and HD, respectively, and combined drought and heat stress can induce specific AS compared with individual one. In addition, wheat homeologous genes exhibited differential AS responses under stress conditions that more AS events occurred on B subgenome than on A and D genomes. Comparison of genes regulated at AS and transcriptional levels showed that only 12% of DS-induced AS genes were subjected to transcriptional regulation, whereas the proportion increased to ~40% under HS and HD. Functional enrichment analysis revealed that abiotic stress-responsive pathways tended to be highly overrepresented among these overlapped genes under HS and HD. Thus, we proposed that transcriptional regulation may play a major role in response to DS, which coordinates with AS regulation to contribute to HS and HD tolerance in wheat.

MicroRNA1917 targets CTR4 splice variants to regulate ethylene responses in tomato

Ethylene perception is regulated by receptors, and the downstream protein CONSTITUTIVE TRIPLE RESPONSE1 is a key suppressor of ethylene signalling. The non-conserved tomato (Solanum lycopersicum) microRNA1917 (Sly-miR1917) mediates degradation of SlCTR4 splice variants (SlCTR4sv) but the molecular details of this pathway remain unknown. Sly-miR1917 and the targeted SlCTR4sv are ubiquitously expressed in all tomato organs. Overexpression of Sly-miR1917 enhances ethylene responses, including the triple response in etiolated seedlings, in the absence of ethylene, as well as epinastic petiole growth, accelerated pedicel abscission, and fruit ripening. Enhanced ethylene signalling in Sly-miR1917-overexpressing plants (1917-OE) is accompanied by up-regulation of ethylene biosynthesis and signalling genes, and increased ethylene emission. These phenotypes were recovered by repressing the positive ethylene regulator EIN2. Moreover, the Sly-miR1917-targeted SlCTR4 splice variant SlCTR4sv3, expressed specifically in the abscission zone, exhibited the opposite expression pattern to Sly-miR1917. Complementation of the Arabidopsis thaliana ctr-1 mutant and yeast two-hybrid and bimolecular fluorescence complementation assays suggested that SlCTR4sv3 functions in ethylene signalling. Co-expression of Sly-miR1917 and SlCTR4sv3 in Nicotiana benthamiana further suggested that Sly-miR1917 cleaves SlCTR4sv3 in vivo. Database homology searching revealed a Solanum tuberosum CTR-like splice variant containing a Sly-miR1917 binding sequence, and a homologue of mature Sly-miR1917 in potato, indicating a conserved function for miR1917 and the regulatory miRNA-mediated ethylene network in solanaceous species.

Comparative epigenomics reveals evolution of duplicated genes in potato and tomato

The evolution of duplicated genes after polyploidization has been the subject of many evolutionary biology studies. Potato (Solanum tuberosum) and tomato (Solanum lycopersicum) are the first two sequenced genomes of asterids, and share a common polyploidization event. However, the epigenetic role of DNA methylation on the evolution of duplicated genes derived from polyploidization is not fully understood. Here, we explore the role of the DNA methylation in the evolution of duplicated genes in potato and tomato. The overall levels of DNA methylation are different, although patterns of DNA methylation are similar in potato and tomato. Different types of duplicated genes can display different methylation patterns in potato and tomato. In addition, we found that differences in the methylation levels between duplicated genes were associated with gene expression divergence. In particular, for the majority of duplicated gene pairs, one copy is always hyper- or hypo-methylated compared with the other copy across different tomato fruit ripening stages, and these genes are enriched for specific function related to transcription factor (TF) activity. Furthermore, transcription of hundreds of duplicated TFs was shown to be regulated by DNA methylation during fruit ripening stages in tomato, some of which are well-known fruit ripening TFs. Taken together, our results support the notion that DNA methylation may facilitate divergent evolution of duplicated genes and play roles in important biological processes such as tomato fruit ripening.

RNA-seq analysis reveals alternative splicing under salt stress in cotton, Gossypium davidsonii

BACKGROUND

Numerous studies have focused on the regulation of gene expression in response to salt stress at the transcriptional level; however, little is known about this process at the post-transcriptional level.

RESULTS

Using a diploid D genome wild salinity-tolerant cotton species, Gossypium davidsonii, we analyzed alternative splicing (AS) of genes related to salt stress by comparing high-throughput transcriptomes from salt-treated and well-watered roots and leaves. A total of 14,172 AS events were identified involving 6798 genes, of which intron retention (35.73%) was the most frequent, being detected in 3492 genes. Under salt stress, 1287 and 1228 differential alternative splicing (DAS) events were identified in roots and leaves, respectively. These DAS genes were associated with specific functional pathways, such as "responses to stress", "metabolic process" and "RNA splicing", implying that AS represents an important pathway of gene regulation in response to salt stress. Several salt response genes, such as pyrroline-5-carboxylate synthase (P5CS), K⁺ channel outward (KCO1), plasma membrane intrinsic protein (PIP) and WRKY33 which were involved in osmotic balance, ion homeostasis, water transportation and transcriptional regulation, respectively, were identified with differential alternative splicing under salt stress. Moreover, we revealed that 13 genes encoding Ser/Arg-rich (SR) proteins related to AS regulation were differentially alternatively spliced under salt stress.

CONCLUSION

This study first provide a comprehensive view of AS in G. davidsonii, and highlight novel insights into the potential roles of AS in plant responses to salt stress.

The comparison of alternative splicing among the multiple tissues in cucumber

BACKGROUND

Alternative splicing (AS) is an important post-transcriptional process. It has been suggested that most AS events are subject to tissue-specific regulation. However, the global dynamics of AS in different tissues are poorly explored.

RESULTS

To analyse global changes in AS in multiple tissues, we identified the AS events and constructed a comprehensive catalogue of AS events within each tissue based on the genome-wide RNA-seq reads from ten tissues in cucumber. First, we found that 58% of the multi-exon genes underwent AS. We further obtained 565 genes with significantly more AS events compared with random genes. These genes were found significant enrichment in biological processes related to the regulation of actin filament length. Second, significantly different AS event profiles among ten tissues were found. The tissues with the same origin of development are more likely to have a relatively similar AS profile. Moreover, 7370 genes showed tissue-specific AS events and were highly enriched in biological processes related to the positive regulation of cellular component organization. Root-specificity AS genes were related to the cellular response to DNA damage stimulus. Third, the genes with different intron retention (IR) patterns among the ten tissues showed significant difference in GC percentages of the retained intron, and the number of exons and FPKM of the major transcripts.

CONCLUSIONS

Our study provided a comprehensive view of AS in multiple tissues. We revealed novel insights into the patterns of AS in multiple tissues and the tissue-specific AS in cucumber.

Abiotic Stresses Modulate Landscape of Poplar Transcriptome via Alternative Splicing, Differential Intron Retention, and Isoform Ratio Switching

Abiotic stresses affect plant physiology, development, growth, and alter pre-mRNA splicing. Western poplar is a model woody tree and a potential bioenergy feedstock. To investigate the extent of stress-regulated alternative splicing (AS), we conducted an in-depth survey of leaf, root, and stem xylem transcriptomes under drought, salt, or temperature stress. Analysis of approximately one billion of genome-aligned RNA-Seq reads from tissue- or stress-specific libraries revealed over fifteen millions of novel splice junctions. Transcript models supported by both RNA-Seq and single molecule isoform sequencing (Iso-Seq) data revealed a broad array of novel stress- and/or tissue-specific isoforms. Analysis of Iso-Seq data also resulted in the discovery of 15,087 novel transcribed regions of which 164 show AS. Our findings demonstrate that abiotic stresses profoundly perturb transcript isoform profiles and trigger widespread intron retention (IR) events. Stress treatments often increased or decreased retention of specific introns - a phenomenon described here as differential intron retention (DIR). Many differentially retained introns were regulated in a stress- and/or tissue-specific manner. A subset of transcripts harboring super stress-responsive DIR events showed persisting fluctuations in the degree of IR across all treatments and tissue types. To investigate coordinated dynamics of intron-containing transcripts in the study we quantified absolute copy number of isoforms of two conserved transcription factors (TFs) using Droplet Digital PCR. This case study suggests that stress treatments can be associated with coordinated switches in relative ratios between fully spliced and intron-retaining isoforms and may play a role in adjusting transcriptome to abiotic stresses.

Comprehensive Workflow for the Genome-wide Identification and Expression Meta-analysis of the ATL E3 Ubiquitin Ligase Gene Family in Grapevine

Classification and nomenclature of genes in a family can significantly contribute to the description of the diversity of encoded proteins and to the prediction of family functions based on several features, such as the presence of sequence motifs or of particular sites for post-translational modification and the expression profile of family members in different conditions. This work describes a detailed protocol for gene family characterization. Here, the procedure is applied to the characterization of the Arabidopsis Tóxicos in Levadura (ATL) E3 ubiquitin ligase family in grapevine. The methods include the genome-wide identification of family members, the characterization of gene localization, structure, and duplication, the analysis of conserved protein motifs, the prediction of protein localization and phosphorylation sites as well as gene expression profiling across the family in different datasets. Such procedure, which could be extended to further analyses depending on experimental purposes, could be applied to any gene family in any plant species for which genomic data are available, and it provides valuable information to identify interesting candidates for functional studies, giving insights into the molecular mechanisms of plant adaptation to their environment.