Sage Journals: Discover world-class research

Abstract

Communication systems within and between plant cells involve the transfer of ions and molecules between compartments, and are essential for development and responses to biotic and abiotic stresses. This in turn requires the regulated movement and fusion of membrane systems with their associated cargo. Recent advances in genomics has provided new resources with which to investigate the evolutionary relationships between membrane proteins across plant species. Members of the soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) are known to play important roles in vesicle trafficking across plant, animal and microbial species. Using recent public expression and transcriptomic data from 9 representative green plants, we investigated the evolution of the SNARE classes and linked protein changes to functional specialization (expression patterns). We identified an additional 3 putative SNARE genes in the model plant Arabidopsis. We found that all SNARE classes have expanded in number to a greater or lesser degree alongside the evolution of multicellularity, and that within-species expansions are also common. These gene expansions appear to be associated with the accumulation of amino acid changes and with sub-functionalization of SNARE family members to different tissues. These results provide an insight into SNARE protein evolution and functional specialization. The work provides a platform for hypothesis-building and future research into the precise functions of these proteins in plant development and responses to the environment.

Keywords

SNAREs vesicle trafficking membrane fusion phylogenetics Viridiplantae

Introduction

Communication between cells, and between compartments within cells, depends on membrane structure and function, and is fundamental to plant and animal growth and development and responses to environmental stresses. The regulated dynamics of membranes, and the associated proteins, is mediated by the membrane trafficking system, and is of great current interest to those interested in understanding the cell as an integrated system of signals and molecular responses. The membrane trafficking system in plant cells comprises the biosynthetic secretory pathway, the endocytic pathway and the vesicle transport pathway.¹ The organelles involved in these processes use small membrane-enclosed transport vesicles to exchange molecular information.

There are 4 essential steps in the vesicle trafficking system, namely budding, vesicle movement, tethering and fusion (Figure 1). During budding, coat proteins and dynamin-related GTPases in a donor compartment are used to form a vesicle and deform the local membrane until a vesicle is freed by scission. Cargo and vesicle (v)-SNAREs (soluble N-ethylmaleimide-sensitive factor attachment protein receptors) are incorporated into the budding vesicle by binding to coat subunits. In the movement phase, the freed vesicle moves towards the acceptor compartment by association with cytoskeletal motors. Molecular motors including kinesin and myosin have all been shown to be involved in this process.^2,3 Then tethering and docking factors work in conjunction with Rab GTPases and SNARE proteins to tether the vesicle to their acceptor membrane.⁴ In the final, fusion step, a tetrameric target (t)-SNARE complex is formed from a single v-SNARE molecule (members of the synaptobrevin or VAMP family of proteins) and a trimeric target membrane t-SNARE complex (members of the syntaxin and SNAP-25 families) that allows vesicles to identify their target compartment and complete membrane fusion and cargo delivery.^5-9

Figure 1.

The vesicle trafficking pathway.

Previous reports have suggested there are 61 SNARE genes represented in the Arabidopsis genome,^10-12 and the encoded proteins are located in different subcellular compartments and are essential for membrane fusion (Supplemental Figure 1; Supplemental Tables 1 and 2). SNARE proteins form a superfamily of small proteins of 100 to 300 amino acids that are characterized by a conserved protein motif called the SNARE domain, a membrane-anchor domain in the C terminus (TM) and a variable N-terminal domain.⁶ The SNARE domain contains a SNARE motif of 60 to 70 amino acids, and consists of a heptad repeat forming a coiled-coil structure, via hetero-oligomeric interactions.^6,8,13 SNAREs can be classified on either a functional or structural basis. Classification on the basis of function puts SNAREs into vesicle-associated (v-SNARE) and target membrane-associated (t-SNARE) groupings.¹⁴ An alternative type of classification relies on amino acid sequence, though this results in a similar grouping to the v/t classification. Here, SNAREs are classified as either Q-SNAREs or R-SNAREs according to whether they contain a conserved glutamine (Q) or arginine (R) residue in the SNARE domain.¹⁵ t-SNAREs broadly correspond to Q-SNAREs, and v-SNAREs to R-SNAREs. SNAREs can be further subdivided into Qa-, Qb-, Qc-, Qb+c- and R-SNAREs, each with a hydrophobic SNARE domain.^6,16 SNAREs include the vesicle-associated membrane protein (VAMP) known as synaptobrevin in animal systems, with a role in fusing synaptic vesicles to the plasma membrane.^17,18

Most R-SNAREs localize to trafficking vesicles, and are attached via the C-terminal transmembrane (TM) domain. In plants, R-SNAREs form 3 major subfamilies, the Sec22-, YKT6- and VAMP7-like R-SNAREs.^12,19 Vesicle resident R-SNAREs are typically designated VAMPs (vesicle-associated membrane proteins), and similar to other SNARE proteins have 3 functional domains, the N terminal domain, a SNARE motif and the transmembrane domain.²⁰ VAMPs can be further classified into 2 subgroups, that is, with either a short or a long N-terminal sequence: the shorter VAMPs are termed ‘brevins’ and longer VAMPs, ‘longins’.²¹ Based on functions in synaptic exocytosis in mammalian systems, the shorter R-SNAREs are also termed synaptobrevins. However, all plant R-SNAREs are longins.^22,23

In Arabidopsis, VAMP7-like R-SNAREs have been subjected to duplication into 2 gene families of 4 and 8 members respectively. The 4 VAMP71 group proteins are involved in endosomal trafficking;^10,24 while the 8 VAMP72 group proteins are likely specific to the green plant lineage and represents R-SNARE components involved in secretion.^12,19 VAMPs have roles as key proteins in abiotic stress tolerance (VAMP711, VAMP712),^25,26 in gravitropic responses,²⁷ in cell plate formation (VAMP721, VAMP722),^28-33 in cytokinesis,^29,34 in defence responses,^28,31,35 in the transport of phytohormones (eg, Abscisic acid (ABA) and auxin) and in root formation and hormone response in plants.^36,37 These functions are therefore dependent on correct membrane vesicle trafficking and its regulation. However, R-SNAREs may play a key role in determining the specificity in vesicle budding, and an important mechanism for SNARE localization is interaction with vesicle coats. For example, Springer & Schekman³⁸ have shown that R-SNAREs may be components of the COPII vesicles that are involved in ER-Golgi transport. Rein et al³⁹ found that R-SNAREs must be packaged into COPI vesicles during retrieval from the Golgi.

In this article, we analyze phylogenetic relationships between different classes of SNARE proteins across plant species across the plant kingdom, and summarize recent findings on the role of SNARE proteins in the model plant species Arabidopsis thaliana, to provide some insight into the evolution of these important proteins and to consider directions of research to understand better these proteins.

Materials and Methods

Phylogenetic analysis

To understand the phylogenetic relationships between different classes of SNAREs across plant species, we built a database of Arabidopsis SNARE-encoding genes and their homologues in 8 well annotated plant genomes (Brachypodium distachyon (grass), Chlamydomonas reinhardtii (alga), Glycine max (soybean), Oryza sativa (rice), Physcomitrella patens (moss), Populus trichocarpa (poplar), Sorghum bicolor (sorghum), Vitis vinifera (grapevine)). These species span the green plant kingdom (Viridiplantae), capturing various major taxonomic groupings in plants including: Embryophyta (all land plants excluding alga), Trachaeophyta (all vascular plants excluding alga and moss), Moncotyledoneae (rice, sorghum and grass), and Dicotelydoneae (Arabidopsis, soybean, poplar, grapevine). Each Arabidopsis gene on the TAIR database (https://www.arabidopsis.org) was searched to find its list of plant gene homologues that have previously identified by sequence-based methods through the PANTHER (Protein ANalysis THrough Evolutionary Relationships) curated database.⁴⁰ The list of protein accession codes were searched on the UniProt database (https://www.uniprot.org) and fasta formatted amino acid coding sequences were downloaded (listed in Suppl. Table 3). Phylogenetic analysis was performed in R v3.5.3 (https://www.r-project.org) using the R packages: ape v5.3,⁴¹ ggtree,⁴² msa v1.4.3,⁴³ phytools v0.6-99,⁴⁴ phangorn v2.5.5⁴⁵ and seqinr v3.6.1.⁴⁶ The protein lists for each gene class were filtered by gene name to remove duplicated sequences, protein fragments, and allelic variants of the same gene. The remaining amino acid sequences for each gene class were aligned using the Expresso structural alignment algorithm of T-Coffee alignment suite (http://tcoffee.crg.cat/apps/tcoffee/⁴⁷) that prioritizes alignment of protein structural motifs. The best amino acid evolutionary model for the aligned amino acid sequences was identified using the phangorn modelTest function, which was identified to be JTT+G+I in all cases. Maximum likelihood distances were calculated for the model using the phangorn dist.ml function and a best fit tree estimated using the phangorn optim.pml function with options appropriate for the JTT+G+I model. The resulting maximum likelihood trees were tested with 100 bootstrap iterations of the phangorn bootstrap.pml function. Trees were visualized using phytools and ggtree. Clades of related genes that were identified in each tree are highlighted with different colours. To understand the evolutionary relationships between the functional domains of different Arabidopsis SNAREs, the Arabidopsis amino acid sequences were analyzed to detect protein motifs with P < .05 using MEME software (http://meme-suite.org/tools/meme⁴⁸). The amino acid sequences of the genes from other species were searched for regions that were significantly similar to these protein domains using FIMO software (http://meme-suite.org/doc/fimo.html⁴⁹). Some short domains gave multiple hits for each gene so the significance threshold was adjusted from P < 1E × 04 to between P < 1E × 06 and P <1E × 08 that gave no more than 2 hits per domain per gene. The putative protein domains were identified when possible by manually checking the locations of these putative domains against the UniProt protein domain annotation of each Arabidopsis gene. Identified protein domain information was visualized alongside the phylogenetic trees using the gheatmap function of the R package ggtree to highlight the evolutionary gain and loss of motifs in each gene class.

The expression of SNARE protein-coding genes in Arabidopsis

The expression data for 64 SNARE protein-coding genes in 79 organs and developmental stages of Arabidopsis were downloaded from TraVA (Transcriptome Variation Analysis, http://travadb.org/; Supplemental Tables 5 and Table 6).⁵⁰ Expression levels of each gene across samples were scaled to have a mean of zero and standard deviation of one. The expression heatmap was generated using dChip software⁵¹ with the default settings of paired Pearson correlation distances as distance and centroid linkage method for building clusters.

Relationship between genetic distance and expression differences

Paired Pearson correlation distances were calculated from the TraVA gene expression data using the cor.dist function of the R package bioDist v1.58.0.⁵² Coding sequences of each class of SNARE genes were aligned using webPRANK,⁵³ a phylogeny-aware alignment algorithm that maintains the correct reading frame. Maximum likelihood trees were constructed for each class as described above but using the optimal GTR+G+I DNA model of evolution instead. Paired gene maximum likelihood genetic distances were extracted from this tree using the R ape function cophenetic.phylo. Sub-functionalization within each gene class was tested with Mantel tests to investigate relationships between genetic and expression distances using the R phytools multi.mantel function. For the between gene class sub-functionalization tests, mean, standard deviation, and sample size values were calculated for the genetic and expression distance matrices for each gene class. The relationship between mean gene class genetic and expression distance weighted by sample size was then tested with linear models in R. Neo-functionalization was tested by measuring paired gene non-synonymous to synonymous base change ratios (Dn/Ds) for each pair of genes within each webPRANK aligned gene class using the R seqinr kaks function assuming that amino acid changes contribute to changes in protein function. Gene pairs that gave valid Dn/Ds values greater than zero and less than 9 were considered. Mantel tests were used to investigate relationships between Dn/Ds and expression distances. Between gene class neo-functionalization tests were performed with linear models of gene class means weighted by sample size.

Results

Phylogenetic analysis of plant SNARE proteins

A search of PANTHER homologues of Arabidopsis SNARE genes led to the identification of 3 new putative Arabidopsis SNARE homologues; 2 Qc-SNARE homologues, AT1G16225 and AT1G16230; and an R-SNARE homologue, AT3G25013, and many other similar genes across the 8 other target species. The phylogenetic trees of each SNARE gene class are presented in Figure 2. Identified protein domains are summarized in the trees and in Supplemental Table 4. The main evolutionary patterns that the phylogenetic analyses reveal in each gene class are discussed in the following sections.

Figure 2.

Maximum likelihood trees of all 5 classes of SNARE genes in 9 different species representing all green plants.

For the Qa-SNAREs, we identified 8 distinct clades that are highly conserved across plants (Figure 2). Some clades can be traced back into more basal green plant representatives than others and clades also reveal a progressive increase in protein domain complexity. Each of the clades labelled SYP31-32, SYP41-43, and SYP81 could be traced back to the most ancient Viridiplantae grouping including alga. The basal SYP81 clade gains 1 recognisable tSNARE domain from the Trachaeophyta onwards, missing this domain in alga and moss. The basal SYP31-32 clade gains the same tSNARE domain earlier from the Embryophyta onwards, missing it only in the alga. The SYP41-43 clade has this original tSNARE domain throughout and gains a second tSNARE domain from the embryophyta onwards. Some later gene duplication has taken place in these clades with, for example, 2 gene copies in the moss for SYP31-32 and SYP41-43 and up to 4 gene copies of SYP32 present in soybean. The clade SYP21-24 is the next most ancient as it can be traced back to all Embryophyta. This clade is characterized by 4 gene copies in both moss and Arabidopsis and mostly contains the 2 tSNARE domains and a new coiled coil protein domain. These 4 basal clades show sporadic evidence for a transmembrane domain. The remaining clades, SYP131-132, PHYPA, SYP111-112 and SYP121-125, are more closely related to each other and are only found in Embryophyta. There is a diverse 7 member moss-specific PHYPA clade, suggesting the evolution of multiple specialist functions for this SNARE gene family in moss. These results are consistent with the recent phylogeny study of SYP1 subfamily of Qa-SNAREs by Slane et al⁵⁴ The other clades, SYP111-112, SYP121-125, SYP131-132 are only found in Trachaeophyta. These clades have diversified into many gene members with evidence of gene duplication within soybean with 8 and 5 gene copies in the SYP121-125 and SYP111-112 clades, respectively. These 4 clades mostly all contain the transmembrane domain, the 2 tSNARE domains and a 3 preceding coiled coil domains, the SYP-121-125 gaining an extra fourth coiled coil domain. Genes that lack several of these domains, such as B9I5W3_POPTR and K7KQY1_SOYBN have relatively long branch lengths and might be recent pseudogenes copies. Coiled coil domains enable protein interactions and interestingly, the Arabidopsis genes SYP121 and SYP122 function as negative regulators of the programmed cell death reaction, the salicylic acid, jasmonic acid and ethylene signalling pathways,⁵⁵ and the membrane traffic mediated by SYP121 and SYP122 are associated differentially with lots of cargo proteins.⁵⁶ The Arabidopsis genes SYP123, SYP124 and SYP125 that function to specify root hair and pollen tissues, tissues not found in algae or moss, are specific to Trachaeophyta.^57-60

The Qb-SNAREs cluster into 6 clades (Figure 2). All of the clades are present across all Viridiplantae, except GOS11 that is present in all Embryophyta, suggesting that the functions of each of the Qb-SNARE clades are ancient and essential. Gene diversification is evident in these clades. Multiple moss gene copies are present in all clades, possibly reflecting a need for more Qb-SNAREs with increasing multicellularity and adaptation to land. The largest clade is VTI11-14, with multiple copies found in Arabidopsis (4 gene copies plus an isoform), soybean (5 copies), poplar and moss (4 copies each), grapevine (3 copies), and grass (2 copies). The smallest group of Qb-SNAREs is the MEMB11-12 clade, which is nonetheless found in multiple copies in Arabidopsis (2 copies), soybean (2 copies) and moss (3 copies). This clade appears to be missing from grass and poplar, but this may be an artefact of database incompleteness. In terms of protein domains, the clades GOS11, GOS12 and SEC20, show no recognisable domains. The MEMB11-12 clade contains the tSNARE domain from the Qa class, and Trachaeophyta members also contain a shared coiled coil domain. Only 2 moss genes show a recognisable transmembrane domain. The VTI11-14 clade contains this tSNARE and coiled coil domain with an extra coiled coil domain in Embryophyta members of the clade. Some members have a transmembrane domain. The NPSN11-13 clade has an algal member with only the tSNARE domain, 5 moss members that show a mix of 2 tSNARE and 2 coiled coil domains that are distinct from the other clades in this class, and a remaining Trachaeophyta members all show these 4 domains and a transmembrane domain.

For the Qc-SNAREs, there are 6 clades, all except USE11-12 and SYP51-52 that can be traced back to all Viridiplantae (Figure 2). The USE11-12 clade traces to Embryophyta but the moss gene is not at the expected basal position. This clade contains just 1 recognised coiled coil domain in the Dicotyledoneae; Arabidopsis, grapevine, and soybean. The SFT11-12 clade contains 2 tSNARE domains starting from the Trachaeophyta. The BSA14A-B clade contains one of these tSNARE domains in all Viridiplantae and a second distinct tSNARE domain in Trachaeophyta members only. The SYP61 clade contains the same set of domains as the prior clade but not all domains are reliably identified in all clade members. This clade also sometimes shows a transmembrane domain. Interestingly, the SYP51-52 clade is found only in some Dicotyledoneae including soybean, grapevine, and Arabidopsis. There are 4 Arabidopsis gene copies in this clade, including 2 previously uncharacterized genes, At1g16225 and At1g16230. These genes are mainly expressed in root tissue. This clade contains the widespread tSNARE domain, the tSNARE domain of BSA14A-B and SYP61 and a third distinct tSNARE domain in most members. The largest clade is SYP71-73 with 2 copies each for alga, moss, and grass, and poplar, 3 copies in grapevine, and Arabidopsis and 5 copies in soybean. This clade contains a distinct coiled coil domain and the widespread tSNARE domain in all Viridiplantae and the SYP61-like tSNARE domain in all Embryophyta. An occasional transmembrane domain is identified in this clade.

The Qb+cSNAP SNARES can be divided into 3 clades, CHLRE SNAP29, SNAP29+SNAP33 but these clades are relatively poorly supported (Figure 2). The CHLRE clade contains highly divergent alga gene and a grapevine gene. The small SNAP30 clade contains Dicotyledoneae Arabidopsis, poplar and 2 soybean genes. The final larger SNAP29+33 group contains Embryophyta members but the divergent 2 moss genes are not in the expected basal position. As with all other gene classes examined, multiple gene copies were found in many species, most notably soybean with 6 copies. Most class members contained 3 tSNARE domains corresponding to Qb and Qc class tSNARE domains.

The R-SNAREs form 3 large clades distinguished by YKT61-62+SEC22, VAMP711-714, and VAMP721-728 encoding genes that are each conserved in all Viridiplantae (Figure 2). The smallest clade is YKT61-62+SEC22, which is represented at least once in all studied species. The base of this clade is defined by a pair of divergent alga genes, Arabidopsis SEC22, and a soybean gene C6TGR3. These genes except for soybean contain a widespread vSNARE domain. The remaining YKT61-62 subclade is found in all Viridiplantae, suggesting its essential cell trafficking role. This clade contains mostly the widespread vSNARE domain, a second vSNARE domain and up to 3 longin domains. The algal gene A8J6T1 contains just 1 of the 3 longin domains. VAMP711-714 is the next largest clade, again with representatives in all studied species, including probable gene duplicates in poplar (5 copies), soybean (4 copies), Arabidopsis (4 copies), and moss (3 copies). This clade contains the widespread vSNARE domain and 3 longin domains that are distinct from the previous clade. Only the algal gene A8J924 contains a single central longin domain. Dicotyledoneae members of the VAMP711-13 subclade mostly contain a transmembrane domain also. The largest clade is represented by VAMP721-728 with multiple members from all Viridiplantae. Gene expansion is likely to have occurred within Arabidopsis for this clade with VAMP-721, 722, 723, 725, and 726 all forming a monophyletic subclades. The basal alga and moss also show monophyletic clusters of duplicated genes (clusters of 4 genes each). It might be that the role of this clade of genes in SNARE complex formation has allowed for repeated pattern of redundant gene duplication within species. This clade includes a divergent Arabidopsis gene AT3G25013 that was designated as a new putative R-SNARE gene by PANTHER. This clade contains the same vSNARE and 3 longin domains as the VAMP711-14 as well as up to 2 transmembrane domains from Embryophyta onwards. Genes missing several of these domains have relatively long branch lengths and are probably pseudogenes.

The expression of SNARE protein-coding genes in Arabidopsis

The expression of 64 SNARE protein-coding genes in 79 organs and developmental stages of Arabidopsis is shown in Figure 3 (see Methods for details). SNARE-encoding genes show a wide range of expression levels and distinct regulation during Arabidopsis development. Arabidopsis SNARE genes are clustered according to their expression patterns in root, leaf, meristem, anther, flower, silique and seeds. Two SNARE-encoding genes (SYP112 and VTI13) are abundantly expressed in flower. Five SNARE transcripts (SYP121, VAMP722, VAMP723, SYP122 and SNAP33) are abundant in leaf and internode. Twenty SNARE genes (SYP123, SYP31, SYP43, VAMP711, GOS11, SYP124, SYP72, SYP131, YKT62, SYP81, USE12, VAMP725, VAMP726, SYP125, SNAP30, SEC20, SEC22, SYP61, NPSN11 and NPSN13) are highly expressed in the anther. Interestingly, the Arabidopsis genes SYP123, SYP124, SYP125 and SYP131 might be important for exocytosis during pollen tube growth.^54,57-60 SEC22 is essential for gametophyte development and maintenance of Golgi-stack integrity.⁶¹ Seven SNARE gene transcripts (VTI11, YKT61, SYP21, MEMB11, SYP52, SYP51 and VAMP727) are abundant in anther and seed. SNAP29 is abundant in the inflorescence meristem and seeds. Nine (SYP41, SYP73, SFT11, BS14a, USE11, VAMP712, VAMP724, SYP23, VAMP714) are strongly expressed in the silique and seeds, and 5 (GOS12, BS14b, SYP32, SYP71 and SYP132) are abundant in meristem and flower. VTI14 is abundant in inflorescence meristem. Four SNARE genes (SYP22, SFT12, SYP42 and MEMB12) are strongly expressed in anther, silique and seeds. Seven (VAMP721, VAMP728, NPSN12, VAMP713, SYP24, SYP111 and VTI12) are strongly expressed in meristems.

Figure 3.

Expression heat map of SNARE genes in different organs and developmental stages of Arabidopsis.

Relationship between genetic distance and expression differences

Large gene families like SNAREs usually arise through gene duplication.⁶² Redundant duplicated genes tend to have a limited evolutionary lifespan unless they acquire new useful functions. Gene paralogues can acquire functional diversity either through sub-functionalization when different tissues or developmental stages come to predominantly express only 1 of the available gene paralogues, or neo-functionalization, when relaxed evolution during the redundant phase allows gene copies to explore adaptive space and acquire novel functions.^63,64

We explored the first of these hypotheses of gene family diversification for SNARE genes by comparing genetic distances to expression differences within and between each gene class. If duplicated genes gradually accumulate expression differences under weak or absent divergent selection, then a significant positive relationship between genetic distance and expression distance would be expected within each gene class. More rapid accumulation of expression differences between duplicated genes would mask this null relationship. For comparisons between gene classes, we might expect that larger families would show more sub-functionalization compared to smaller families. Analysis of the relationships between genetic and expression distances suggested that no significant positive relationships were found for any gene class suggesting that the hypothesis of rapid sub-functionalization of gene duplicates cannot be ruled out within gene classes (Supplemental Figure 2). No relationship was found between genetic and expression distance across gene families failing to support the hypothesis of sub-functionalization at this level (Figure 4).

Figure 4.

Relationships between (A, left) genetic distance versus expression distances, and (B, right) Dn/Ds ratios versus expression distances across SNARE gene classes in Arabidopsis.

We explored the second hypothesis of gene diversification by neo-functionalization. According to the neo-functionalization hypothesis, divergent selection for amino acid changes between duplicated genes would weaken an otherwise longer-term positive relationship for amino acid changes to accumulate over longer periods of time. The results of Mantel tests showed that, except for Qa, there were no significant positive relationships between Dn/Ds and expression distances/genetic distances at either the within-or between gene class levels (Supplemental Figures 3–5). The largest gene classe, Qa, showed a positive linear relationship between genetic distance and Dn/Ds (r² = 0.61, β = 6.44, P = 0.005) supporting a null hypothesis of gradual functional change (Supplemental Figure 4). The other gene families did not show significant positive linear relationships suggesting that some gene pairs in these families might accumulate putatively functional amino acid differences at a faster rate than others. The need for coordinated interactions between different R class genes in the SNARE core complex might contribute to their relatively rapid amino acid change.

Discussion

During the last few years there have been advances in our understanding of the functions of SNARE proteins, especially in the membrane fusion and vesicle transport pathways. However, the functional redundancy between SNAREs makes it hard to define precisely their fusion specificity in relation to specific cellular or developmental events. For example, functional redundancy exists between VTI11 and VTI12,^65,66 SYP121 and SYP122,^67,68 and VAMP721 and VAMP722 R-SNAREs.²⁸ Homozygous double mutants of each show embryo lethal or severe growth defects, but single gene knockout mutations do not show any obvious phenotypes. SYP22 and SYP23,⁶⁹ VAMP711, VAMP712 and VAMP713⁷⁰ may either form distinct SNARE complexes or they could have overlapping functions for specific developmental roles. The redundancy between SNAREs makes it difficult to identify specific SNARE functions if we use traditional experimental approaches.⁷¹

Our phylogenetic analysis provides insights into patterns and drivers of SNARE family gene diversification across green plants. Each SNARE class could be broken down into subclades, many of which had genes representing all studied species, attesting to their ancient evolutionary origins. Some sets of homologous genes including: Qa: SYP21-24, SYP31-32, SYP41-43, SYP81; Qb: SEC20, MEMB11-12, VTI11-14, NPSN11-13; Qc: SYP61, SYP71-73, BS14a-b, SFT11-12; and R: SEC22 plus YKT61-62, VAMP711-714, VAMP721-728 are found in all species examined, indicating that the SNARE complex was inherited vertically throughout all Viridiplantae.¹² Many of these subclades were diverse, with relatively high numbers of gene duplications within some species. Species specific gene duplications frequently occurred in plants. The phylogenetic trees show that many SNARE subclades including: Qa SYP1; Qc USE11-12; and Qc SYP51-52, appeared during the evolution of multicellular land plants (Embryophyta), with a similar story for syntaxins in higher animals.¹² These more recent subclades contain an increased number of protein domains involved in protein-protein interactions and associated with cell membranes. Therefore 1 reason for the diversification of SNAREs might be related to the increased complexity of inter- and intracellular communication systems that are associated with multicellularity. For other SNARE classes, for example: Qa-SNAREs, Qb-NPSN11-13, and R-VAMP721-728, it is less clear why basal green plants such as moss have species-specific subclades with as many as 7 genes. Understanding these phylogenetic relationships will help address our understanding of the evolution of SNARE function. Gene duplications in Arabidopsis seem to be more associated with sub-functionalization to different tissue types than neo-functionalization in terms of the accumulation of amino acid changes, as no associations were found between genetic and expression distance.

Some information on the functions of SNAREs in different tissues is available. Depending on the specific intracellular trafficking process, different v-/t-SNARE complexes are correspondingly formed.¹⁸ This specificity of t-SNARE and v-SNARE complex formation ensures that vesicles are targeted to the correct compartment and induce membrane fusion.⁷² For example, SYP111 takes part in membrane fusion events forming the cell plate and the transport of secretory vesicles at the plasma membrane.^34,73-75 SYP122 plays a role in cell wall deposition and in tethering of donor and target membrane.⁷⁶ SYP81, SYP31, SYP32, GOS11, GOS12, MEMB11, SFT11, SFT12, BS14a, BS14b, SEC22, YKT61 and YKT62 mediate anterograde traffic between ER and Golgi and retrograde traffic within the Golgi apparatus.^5,61,77-82 SEC20, MEMB12, USE11, USE12, SYP71, SYP72 and SYP73 mediate retrograde traffic from Golgi to ER for protein recycling and balance maintenance.^12,29,83-86 SYP51 and SYP52 take part in direct membrane transport from ER to tonoplast and Golgi to mediate vesicle trafficking.^87,88 VAMP714 interacts with SYP121 and SYP22 to mediate vesicle transport from the Golgi apparatus to the vacuole.^89,90

These cell trafficking functions are important mediators of environmental responses. Thus, Q-SNAREs have roles as key proteins in gravitropic responses (VTI11, VTI14),^27,65 in cytokinesis (SYP31, SYP32, NPSN11, NPSN12, NPSN13),^4,10,91-93 in auxin homeostasis (VTI11, SYP41, SYP42, SYP43),^94-99 and in cell plate formation and pathogen resistance (SYP121, SYP131, SYP132, SNAP33).^{35,37,67,100-106}

SYP121 mediates vesicle fusion at the Arabidopsis plasma membrane, and binds the K⁺ channels voltage sensors to coordinate membrane trafficking with K⁺ uptake for growth.¹⁰⁷ Karnik et al.^30,108,109 and Zhang et al.¹¹⁰ found that SM protein SEC11 binds and selectively regulates secretory traffic mediated by N terminus of SYP121 and is important for assembling and recycling of the SNARE during membrane fusion. And Sec1/Munc18 proteins interact with some SNARE proteins to form the exocyst complex, which initiates membrane fusion.^111,112 Moreover, Waghmare et al.⁵⁶ found that the membrane traffic mediated by SYP121 and SYP122 is associated differentially with lots of cargo proteins.

Xue et al²⁶ found that VAMP711 regulates ABA-mediated inhibition of PM H+-ATPase activity and drought stress response by regulating stoma closure. Zhang et al^19,75 found that VAMP721 and VAMP722 interact with the same K+ channels and that this interaction suppresses channel activity, and VAMP721 assembles with SYP121 to coordinate K+ channel gating during SNARE assembly and vesicle fusion. Kim et al.¹¹³ found that CALRETICULIN 1 (CRT1) and CRT2 are critical components in the accumulation of VAMP721 and VAMP722 during ER stress responses. Zhang et al.¹¹⁴ found that SYP22 and VAMP727 mediate the BRI trafficking to PM.

These studies show that SNARE proteins likely function in trafficking at most membranes, but especially at the vacuole and plasma membrane. And VAMP727 forms a complex with SYP22, VTI11, SYP121 and SYP51 to drive membrane fusion and vesicle transport pathways between PVCs and vacuoles.^32,33,115 Löfke et al.⁹⁹ found that auxin can increase the amount of SNARE proteins in the vacuolar membrane. Shirakawa et al.¹¹⁶ found that SYP22 functions in the polarized localization of the auxin efflux carrier PIN1 in the Arabidopsis leaf. Xia et al¹¹⁷ found that SYP132 and associated endocytosis play significant roles in auxin-regulated H⁺-ATPase traffic and associated functions at the plasma membrane. Moreover, we have recently found that AtVAMP714 is also required for the exocytic localization of PIN1 and PIN2 proteins to the plasma membrane via the Golgi, and for polar auxin transport; and the related genes VAMP711, VAMP712 and VAMP713 are also inducible by auxin.¹¹⁸

Conclusion

The analysis presented in this paper shows that SNARE proteins likely play very important roles in plant development, with functions being identified in auxin transport and signalling. Phylogenetic studies show that SNARE diversification is associated with the evolution of multicellularity in plants, therefore providing new insights into the functions of SNARE proteins in vesicle transport pathways. This phylogenetic context can be combined with future molecular genetic findings for greater insights. Knowing where and how those membrane protein complexes originate and are destined in the cell, and their dynamics and relationship with membrane composition that determines for example transporter polarity in a tissue-specific manner,¹¹⁹ will improve our knowledge of the importance of SNARE complexes in vesicle trafficking pathways and in a range of developmental and stress response pathways in higher plants.

Supplemental Material

00f8b8c046636_EvoBioRevSupplTable5y_xyz466362fa485fb – Supplemental material for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins

Supplemental material, 00f8b8c046636_EvoBioRevSupplTable5y_xyz466362fa485fb for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins by Xiaoyan Gu, Adrian Brennan, Wenbin Wei, Guangqin Guo and Keith Lindsey in Evolutionary Bioinformatics

Supplemental Material

2873cc8b46636_EvoBioRevSupplTable4_xyz466366e8025bf – Supplemental material for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins

Supplemental material, 2873cc8b46636_EvoBioRevSupplTable4_xyz466366e8025bf for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins by Xiaoyan Gu, Adrian Brennan, Wenbin Wei, Guangqin Guo and Keith Lindsey in Evolutionary Bioinformatics

Supplemental Material

3185449246636_EvoBioRevSupplTable1_xyz466364bc83175 – Supplemental material for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins

Supplemental material, 3185449246636_EvoBioRevSupplTable1_xyz466364bc83175 for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins by Xiaoyan Gu, Adrian Brennan, Wenbin Wei, Guangqin Guo and Keith Lindsey in Evolutionary Bioinformatics

Supplemental Material

529d24f646636_EvoBioRevSupplTable6_xyz466362abe5173 – Supplemental material for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins

Supplemental material, 529d24f646636_EvoBioRevSupplTable6_xyz466362abe5173 for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins by Xiaoyan Gu, Adrian Brennan, Wenbin Wei, Guangqin Guo and Keith Lindsey in Evolutionary Bioinformatics

Supplemental Material

ef0a8cf646636_EvoBioRevSupplTable3y_xyz466366d196b73 – Supplemental material for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins

Supplemental material, ef0a8cf646636_EvoBioRevSupplTable3y_xyz466366d196b73 for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins by Xiaoyan Gu, Adrian Brennan, Wenbin Wei, Guangqin Guo and Keith Lindsey in Evolutionary Bioinformatics

Supplemental Material

f763b0b446636_EvoBioRevSupplTable2_xyz46636386aed71 – Supplemental material for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins

Supplemental material, f763b0b446636_EvoBioRevSupplTable2_xyz46636386aed71 for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins by Xiaoyan Gu, Adrian Brennan, Wenbin Wei, Guangqin Guo and Keith Lindsey in Evolutionary Bioinformatics

Supplemental Material

Supplementary_Figures_1-5_Gu_et_al. – Supplemental material for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins

Supplemental material, Supplementary_Figures_1-5_Gu_et_al. for Vesicle Transport in Plants: A Revised Phylogeny of SNARE Proteins by Xiaoyan Gu, Adrian Brennan, Wenbin Wei, Guangqin Guo and Keith Lindsey in Evolutionary Bioinformatics

Footnotes

Acknowledgements

Not applicable.

Funding:

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: XG thanks the China Scholarship Council for the funding to support her studies at Durham University. KL acknowledges The Biotechnology and Biological Sciences Research Council (BBS/B/0773X) and Durham University for funding.

Declaration of conflicting interests:

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Author contributions

Phylogenetic analysis was carried out by XG and AB, expression analysis by XG and WW. All authors contributed to the writing of the manuscript.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files.

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

ORCID iDs

Adrian Brennan

Wenbin Wei

Keith Lindsey

Supplemental material

Supplemental material for this article is available online.

References

Uemura

Ueda

Plant vacuolar trafficking driven by RAB and SNARE proteins. Curr Opin Plant Biol. 2014;22:116-121.

Matanis

Akhmanova

Wulf

, et al. Bicaudal-D regulates COPI-independent Golgi–ER transport by recruiting the dynein–dynactin motor complex. Nat Cell Biol. 2002;4:986-992.

Shorter

Beard

Seemann

Dirac-Svejstrup

Warren

Sequential tethering of Golgins and catalysis of SNARE-pin assembly by the vesicle-tethering protein p115. J Cell Biol. 2002;157:45-62.

Bassham

Brandizzi

Otegui

Sanderfoot

AA.

The secretory system of Arabidopsis. Arabidopsis Book. 2008;6:e0116.

Sanderfoot

Raikhel

. The secretory system of Arabidopsis. In: Somerville

Meyerowitz

, eds. The Arabidopsis Book. Rockville, MD: American Society of Plant Biologists; 2003:24.

Jahn

Scheller

RH.

SNAREs-engines for membrane fusion. Nat Rev Mol Cell Biol. 2006;7:631-643.

Cai

Reinisch

Ferro-Novick

Coats, tethers, Rabs, and SNAREs work together to mediate the intracellular destination of a transport vesicle. Dev Cell. 2007;12:671-682.

Bassham

Blatt

MR.

SNAREs: cogs and coordinators in signaling and development. Plant Physiol. 2008;147:1504-1515.

Fujimoto

Ueda

Conserved and plant-unique mechanisms regulating plant post-Golgi traffic. Front Plant Sci. 2012;3:197.

10.

Uemura

Ueda

Ohniwa

Nakano

Takeyasu

Sato

MH.

Systematic analysis of SNARE molecules in Arabidopsis: dissection of the post-Golgi network in plant cells. Cell Struct Func. 2004;29:49-65.

11.

Lipka

Kwon

Panstruga

SNARE-ware: the role of SNARE-domain proteins in plant biology. Annu Rev Cell Dev Biol. 2007;23:147-174.

12.

Sanderfoot

AA.

Increases in the number of SNARE genes parallels the rise of multicellularity among the green plants. Plant Physiol. 2007;144:6-17.

13.

Rehman

Sansebastiano

. SNARE proteins as signaling elements. In: Hakeem

Rehman

Tahir

, eds. Plant Signaling: Understanding the Molecular Crosstalk. New Delhi: Springer; 2014:39-49.

14.

Söllner

Whiteheart

Brunner

, et al. SNAP receptors implicated in vesicle targeting and fusion. Nature. 1993;362:318-324.

15.

Fasshauer

Sutton

Brunger

Jahn

Conserved structural features of the synaptic fusion complex: SNARE proteins reclassified as Q- and R-SNAREs. Proc Natl Acad Sci U S A. 1998;95:15781-15786.

16.

Bock

Matern

Peden

Scheller

RH.

A genomic perspective on membrane compartment organization. Nature. 2001;409:839-841.

17.

Walch-Solimena

Takei

Marek

, et al. Synaptotagmin: a membrane constituent of neuropeptide-containing large dense-core vesicles. J Neurosci. 1993;13:3895-3903.

18.

Bonifacino

Glick

BS.

The mechanisms of vesicle budding and fusion. Cell. 2004;116:153-166.

19.

Zhang

Karnik

Wang

Wallmeroth

Blatt

Grefen

The Arabidopsis R-SNARE VAMP721 interacts with KAT1 and KC1 K⁺ channels to moderate K⁺ current at the plasma membrane. Plant Cell. 2015;27:1697-1717.

20.

Zhang

Wang

Zheng

Requirement of R-SNAREs VAMP721 and VAMP722 for the gametophyte activity, embryogenesis and seedling root development in Arabidopsis. Plant Growth Regul. 2015;77:57-65.

21.

Rossi

Banfield

Vacca

, et al. Longins and their longin domains: regulated SNAREs and multifunctional SNARE regulators. Trends Biochem Sci. 2004;29:682-688.

22.

Martinez-Arca

Rudge

Vacca

, et al. A dual mechanism controlling the localization and function of exocytic v-SNAREs. Proc Natl Acad Sci U S A. 2003;100:9011-9016.

23.

Uemura

Sato

Takeyasu

The longin domain regulates subcellular targeting of VAMP7 in Arabidopsis thaliana. FEBS Lett. 2005;579:2842-2846.

24.

Hong

SNAREs and traffic. Biochim Biophys Acta. 2005;1744:120-144.

25.

Leshem

Golani

Kaye

Levine

Reduced expression of the v-SNAREs AtVAMP71/AtVAMP7C gene family in Arabidopsis reduces drought tolerance by suppression of abscisic acid-dependent stomatal closure. J Exp Bot. 2010;61:2615-2622.

26.

Xue

Yang

Wang

Guo

VAMP711 is required for abscisic acid-mediated inhibition of plasma membrane H⁺-ATPase activity. Plant Physiol. 2018;178:1332-1343.

27.

Yano

Sato

Saito

Sato

Morita

Tasaka

A SNARE complex containing SGR3/AtVAM3 and ZIG/VTI11 in gravity-sensing cells is important for Arabidopsis shoot gravitropism. Proc Natl Acad Sci U S A. 2003;100:8589-8594.

28.

Zhang

Liu

Hao

Jin

Lin

Arabidopsis R-SNARE proteins VAMP721 and VAMP722 are required for cell plate formation. PLoS One. 2011;6:e26129.

29.

Karnik

Grefen

Bayne

, et al. Arabidopsis Sec1/Munc18 protein SEC11 is a competitive and dynamic modulator of SNARE binding and SYP121-dependent vesicle traffic. Plant Cell. 2013;25:1368-1382.

30.

Yun

Kwaaitaal

Kato

, et al. Requirement of vesicle-associated membrane protein 721 and 722 for sustained growth during immune responses in Arabidopsis. Mol Cells. 2013;35:481-488.

31.

Zhang

Zhao

Zheng

Wang

JY.

The analysis of mutant phenotypes and tissue expression reveals a role of SNAREs VAMP721 and VAMP722 in seedling growth. Biol Plant. 2017;61:275-283.

32.

Uemura

Nakano

Takagi

, et al. A Golgi-released subpopulation of the trans-Golgi network mediates protein secretion in Arabidopsis. Plant Physiol. 2019;179:519-532.

33.

EI-Kasmi

Krause

Hiller

, et al. SNARE complexes of different composition jointly mediate membrane fusion in Arabidopsis cytokinesis. Mol Biol Cell. 2013;24:1593-1601.

34.

Collins

Thordal-Christensen

Lipka

, et al. SNARE-protein-mediated disease resistance at the plant cell wall. Nature. 2003;425:973-977.

35.

Kwon

Neu

Pajonk

, et al. Co-option of a default secretory pathway for plant immune responses. Nature. 2008;451:835-840.

36.

Dacks

Doolittle

WF.

Novel syntaxin gene sequences from Giardia, Trypanosoma and algae: implications for the ancient evolution of the eukaryotic endomembrane system. J Cell Sci. 2002;115:1635-1642.

37.

Enami

Ichikawa

Uemura

, et al. Differential expression control and polarized distribution of plasma membrane-resident SYP1 SNAREs in Arabidopsis thaliana. Plant Cell Physiol. 2009;50:280-289.

38.

Springer

Schekman

Nucleation of COPII vesicular coat complex by endoplasmic reticulum to Golgi vesicle SNAREs. Science. 1998;281:698-700.

39.

Rein

Andag

Duden

Schmitt

Spang

ARF-GAP–mediated interaction between the ER-Golgi v-SNAREs and the COPI coat. J Cell Biol. 2002;157:395-404.

40.

Muruganujan

Ebert

Huang

Thomas

PD.

PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 2019;47:D419-D42. doi:10.1093/nar/gky1038.

41.

Paradis

Claude

Strimmer

APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289-290.

42.

Smith

Zhu

Guan

Lam

TTY

. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8:28-36.

43.

Bodenhofer

Bonatesta

Horejš-Kainrath

Hochreiter

msa: an R package for multiple sequence alignment. Bioinformatics. 2015;31:3997-3999.

44.

Revell

LJ.

phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217-223.

45.

Schliep

KP.

phangorn: phylogenetic analysis in R. Bioinformatics. 2011;27:592-593.

46.

Charif

Lobry

. SeqinR 1.0-2: a contributed package to the R project for statistical computing devoted to biological sequences retrieval and analysis. In: Bastolla

Porto

Roman

Vendruscolo

, eds. Structural Approaches to Sequence Evolution: Molecules, Networks, Populations. Series Biological and Medical Physics, Biomedical Engineering. New York: Springer Verlag; 2007:207-232.

47.

Notredame

Higgins

Heringa

T-Coffee: a novel method for multiple sequence alignments. J Mol Biol. 2000;302:205-217.

48.

Bailey

Elkan

. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. In: Proceedings of the Second International Conference on Intelligent Systems for Molecular Biology. Menlo Park, CA: AAAI Press; 1994:28-36.

49.

Grant

Bailey

Noble

WS.

FIMO: scanning for occurrences of a given motif. Bioinformatics. 2011;27:1017-1018.

50.

Klepikova

Kasianov

Gerasimov

Logacheva

Penin

AA.

A high resolution map of the Arabidopsis thaliana developmental transcriptome based on RNA-seq profiling. Plant J. 2016;88:1058-1070.

51.

Wong

WH.

Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci U S A. 2001;98:31-36.

52.

Ding

Gentleman

Carey

bioDist: different distance measures. 2019; R package version 1.58.0.

53.

Löytynoja

Goldman

webPRANK: a phylogeny-aware multiple sequence aligner with interactive alignment browser. BMC Bioinformatics. 2010;11:579.

54.

Slane

Reichardt

El-Kasmi

Bayer

Jürgens

Evolutionarily diverse SYP1 Qa-SNAREs jointly sustain pollen tube growth in Arabidopsis. Plant J. 2017;92:375-385.

55.

Zhang

Feechan

Pedersen

, et al. A SNARE-protein has opposing functions in penetration resistance and defence signalling pathways. Plant J. 2007;49:302-312.

56.

Waghmare

Lileikyte

Karnik

Goodman

Blatt

Jones

AME

. SNAREs SYP121 and SYP122 mediate the secretion of distinct cargo subsets. Plant Physiol. 2018;178:1679-1688.

57.

Kato

Steger

AP.

A systems model of vesicle trafficking in Arabidopsis pollen tubes. Plant Physiol. 2010;152:590-601.

58.

Silva

Ul-Rehman

Rato

Di Sansebastiano

Malhó

Asymmetric localization of Arabidopsis SYP124 syntaxin at the pollen tube apical and sub-apical zones is involved in tip growth. BMC Plant Biol. 2010;10:179.

59.

Ul-Rehman

Silva

Malho

Localization of Arabidopsis SYP125 syntaxin in the plasma membrane sub-apical and distal zones of growing pollen tubes. Plant Signal Behav. 2011;6:665-670.

60.

Ichikawa

Hirano

Enami

, et al. Syntaxin of plant proteins SYP123 and SYP132 mediate root hair tip growth in Arabidopsis thaliana. Plant Cell Physiol. 2014;55:790-800.

61.

El-Kasmi

Pacher

Strompen

, et al. Arabidopsis SNARE protein SEC22 is essential for gametophyte development and maintenance of Golgi-stack integrity. Plant J. 2011;66:268-279.

62.

Voordeckers

Pougach

Verstrepen

KJ.

How do regulatory networks evolve and expand throughout evolution?

Curr Opin Biotechnol. 2015;34:180-188.

63.

Conant

Wolfe

KH.

Turning a hobby into a job: how duplicated genes find new functions. Nat Rev Genet. 2008;9:938-950.

64.

Viaene

Vekemans

Irish

, et al. Pistillata—duplications as a mode for floral diversification in (basal) asterids. Mol Biol Evol. 2009;26:2627-45.

65.

Kato

Morita

Fukaki

, et al. SGR2, a phospholipase-like protein, and ZIG/SGR4, a SNARE, are involved in the shoot gravitropism of Arabidopsis. Plant Cell. 2002;14:33-46.

66.

Surpin

Zheng

Morita

, et al. The VTI family of SNARE proteins is necessary for plant viability and mediates different protein transport pathways. Plant Cell. 2003;15:2885-2899.

67.

Assaad

Qiu

Youngs

, et al. The PEN1 syntaxin defines a novel cellular compartment upon fungal attack and is required for the timely assembly of papillae. Mol Biol Cell. 2004;15:5118-5129.

68.

Zhang

Lenk

Andersson

, et al. A lesion-mimic syntaxin double mutant in Arabidopsis reveals novel complexity of pathogen defense signaling. Mol Plant. 2008;1:510-527.

69.

Ohtomo

Ueda

Shimada

, et al. Identification of an allele of VAM3/SYP22 that confers a semi-dwarf phenotype in Arabidopsis thaliana. Plant Cell Physiol. 2005;46:1358-1365.

70.

Kim

Brandizzi

News and views into the SNARE complexity in Arabidopsis. Front Plant Sci. 2012;3:28.

71.

Di Sansebastiano

. Defining new SNARE functions: the i-SNARE. Front Plant Sci. 2013;4:99.

72.

Scales

Chen

Yoo

Patel

Doung

Scheller

RH.

SNAREs contribute to the specificity of membrane fusion. Neuron. 2000;26:457-464.

73.

Lukowitz

Mayer

Jürgens

Cytokinesis in the Arabidopsis embryo involves the syntaxin related KNOLLE gene product. Cell. 1996;84:61-71.

74.

Kim

O’Connell

Maekawa-Yoshikawa

Uemura

Neumann

Schulze-Lefert

The powdery mildew resistance protein RPW8.2 is carried on VAMP721/722 vesicles to the extrahaustorial membrane of haustorial complexes. Plant J. 2014;79:835-847.

75.

Zhang

Karnik

Waghmare

Donald

Blatt

MR.

VAMP721 conformations unmask an extended motif for K⁺ channel binding and gating control. Plant Physiol. 2017;173:536-551.

76.

Rehman

Stigliano

Lycett

, et al. Tomato Rab11a characterization evidenced a difference between SYP121-dependent and SYP122-dependent exocytosis. Plant Cell Physiol. 2008;49:751-766.

77.

McNew

Søgaard

Lampen

, et al. Ykt6p, a prenylated SNARE essential for endoplasmic reticulum-Golgi transport. J Biol Chem. 1997;272:17776-17783.

78.

Tsui

Banfield

DK.

Yeast Golgi SNARE interactions are promiscuous. J Cell Sci. 2000;113:145-152.

79.

Parlati

McNew

Fukuda

Miller

Söllner

Rothman

JE.

Topological restriction of SNARE dependent membrane fusion. Nature. 2000;407:194-198.

80.

Chatre

Brandizzi

Hocquellet

Hawes

Moreau

Sec22 and Memb11 are v-SNAREs of the anterograde endoplasmic reticulum-Golgi pathway in tobacco leaf epidermal cells. Plant Physiol. 2005;139:1244-1254.

81.

Bubeck

Scheuring

Hummel

, et al. The syntaxins SYP31 and SYP81 control ER–Golgi trafficking in the plant secretory pathway. Traffic. 2008;9:1629-1652.

82.

Tarte

Seok

Woo

, et al. Arabidopsis Qc-SNARE gene AtSFT12 is involved in salt and osmotic stress responses and Na⁺ accumulation in vacuoles. Plant Cell Rep. 2015;34:1127-1138.

83.

Burri

Varlamov

Doege

, et al. A SNARE required for retrograde transport to the endoplasmic reticulum. Proc Natl Acad Sci U S A. 2003;100:9873-9877.

84.

Dilcher

Veith

Chidambaram

Hartmann

Schmitt

von Mollard

GF.

Use1p is a yeast SNARE protein required for retrograde traffic to the ER. EMBO J. 2003;22:3664-3674.

85.

Suwastika

Uemura

Shiina

Sato

Takeyasu

SYP71, a plant-specific Qc-SNARE protein, reveals dual localization to the plasma membrane and the endoplasmic reticulum in Arabidopsis. Cell Struct Funct. 2008;33:185-192.

86.

Kanazawa

Era

Minamino

, et al. SNARE molecules in Marchantia polymorpha: unique and conserved features of the membrane fusion machinery. Plant Cell Physiol. 2016;57:307-324.

87.

Sanderfoot

Kovaleva

Bassham

Raikhel

NV.

Interactions between syntaxins identify at least five SNARE complexes within the Golgi/prevacuolar system of the Arabidopsis cell. Mol Biol Cell. 2001;12:3733-3743.

88.

Barozzi

Papadia

Stefano

, et al. Variation in membrane trafficking linked to SNARE AtSYP51 interaction with aquaporin NIP1;1. Front Plant Sci. 2019;9:1949.

89.

Szponarski

Sommerer

Boyer

Rossignol

Gibrat

Large-scale characterization of integral proteins from Arabidopsis vacuolar membrane by two-dimensional liquid chromatography. Proteomics. 2004;4:397-406.

90.

Cao

, et al. OsSYP121 accumulates at fungal penetration sites and mediates host resistance to rice blast. Plant Physiol. 2019;179:1330-1342. doi:10.1104/pp.18.01013.

91.

Heese

Gansel

Sticher

, et al. Functional characterization of the KNOLLE-interacting t-SNARE AtSNAP33 and its role in plant cytokinesis. J Cell Biol. 2001;155:239-249.

92.

Zheng

Bednarek

Sanderfoot

Alonso

Ecker

Raikhel

NV.

NPSN11 is a cell plate-associated SNARE protein that interacts with the syntaxin KNOLLE. Plant Physiol. 2002;129:530-539.

93.

Rancour

Dickey

Park

Bednarek

SY.

Characterization of AtCDC48. Evidence for multiple membrane fusion mechanisms at the plane of cell division in plants. Plant Physiol. 2002;130:1241-1253.

94.

Zheng

von-Mollard

Kovaleva

Stevens

Raikhel

NV.

The Plant vesicle-associated SNARE AtVTI1a likely mediates vesicle transport from the trans-Golgi network to the prevacuolar compartment. Mol Biol Cell. 1999;10:2251-2264.

95.

Bassham

Sanderfoot

Kovaleva

Zheng

Raikhel

NV.

AtVPS45 complex formation at the trans-Golgi network. Mol Biol Cell. 2000;11:2251-2265.

96.

Uemura

Kim

Saito

, et al. Qa-SNAREs localized to the trans-Golgi network regulate multiple transport pathways and extracellular disease resistance in plants. Proc Natl Acad Sci U S A. 2012;109:1784-1789.

97.

Ebine

Okatani

Uemura

, et al. A SNARE complex unique to seed plants is required for protein storage vacuole biogenesis and seed development of Arabidopsis thaliana. Plant Cell. 2008;20:3006-3021.

98.

Fujiwara

Uemura

Ebine

, et al. Interactomics of Qa-SNARE in Arabidopsis thaliana. Plant Cell Physiol. 2014;55:781-789.

99.

Löfke

Dünser

Scheuring

Vehn

JK.

Auxin regulates SNARE-dependent vacuolar morphology restricting cell size. ELife. 2015;4:e05868.

100.

Kargul

Gansel

Tyrrell

Sticher

Blatt

MR.

Protein-binding partners of the tobacco syntaxin NtSyr1. FEBS Lett. 2001;508:253-258.

101.

Sutter

Campanoni

Tyrrell

Blatt

MR.

Selective mobility and sensitivity to SNAREs is exhibited by the Arabidopsis KAT1 K1 channel at the plasma membrane. Plant Cell. 2006;18:935-954.

102.

Dhonukshe

Baluska

Schlicht

, et al. Endocytosis of cell surface material mediates cell plate formation during plant cytokinesis. Devel Cell. 2006;10:137-150.

103.

Kalde

Nuhse

Findlay

Peck

SC.

The syntaxin SYP132 contributes to plant resistance against bacteria and secretion of pathogenesis-related protein 1. Proc Natl Acad Sci U S A. 2007;104:11850-11855.

104.

Grefen

Blatt

MR.

SNAREs-molecular governors in signalling and development. Curr Opin Plant Biol. 2008;11:600-609.

105.

Reichardt

Stierhof

Mayer

, et al. Plant cytokinesis requires de novo secretory trafficking but not endocytosis. Curr Biol. 2007;17:2047-2053.

106.

Reichardt

Slane

El-Kasmi

, et al. Mechanisms of functional specificity among plasma-membrane syntaxins in Arabidopsis. Traffic. 2011;12:1269-1280.

107.

Grefen

Karnik

Larson

, et al. A vesicle-trafficking protein commandeers Kv channel voltage sensors for voltage-dependent secretion. Nat Plants. 2015;1:1-12.

108.

Karnik

Zhang

Waghmare

Aderhold

Grefen

Blatt

MR.

Binding of SEC11 indicates its role in SNARE recycling after vesicle fusion and identifies two pathways for vesicular traffic to the plasma membrane. Plant Cell. 2015;27:675-694. doi:10.1105/tpc.114.134429.

109.

Karnik

Waghmare

Zhang

, et al. Commandeering channel voltage sensors for secretion, cell turgor, and volume control. Trends Plant Sci. 2017;22:81-95.

110.

Zhang

Karnik

Alvim

Donald

Blatt

MR.

Dual sites for SEC11 on the SNARE SYP121 implicate a binding exchange during secretory traffic. Plant Physiol. 2019;180:228-239.

111.

Žárský

Cvrcková

Potocký

Hála

Exocytosis and cell polarity in plants - exocyst and recycling domains. New Phytol. 2009;183:255-272.

112.

Vukašinovic

Žárský

Tethering complexes in the Arabidopsis endomembrane system. Front Cell Dev Biol. 2016;4:1-13.

113.

Kim

Choi

Kwon

Yun

HS.

Endoplasmic reticulum stress-induced accumulation of VAMP721/722 requires CALRETICULIN 1 and CALRETICULIN 2 in Arabidopsis. J Integr Plant Biol. 2019;61:974-980.

114.

Zhang

Liu

Zhu

, et al. SYP22 and VAMP727 regulate BRI1 plasma membrane targeting to control plant growth in Arabidopsis. New Phytol. 2019;223:1059-1065.

115.

Takemotoa

Ebinea

Askanid

, et al. Distinct sets of tethering complexes, SNARE complexes, and Rab GTPases mediate membrane fusion at the vacuole in Arabidopsis. Proc Natl Acad Sci U S A. 2018;115:E2457-E2466.

116.

Shirakawa

Ueda

Shimada

Nishiyama

Hara-Nishimura

Vacuolar SNAREs function in the formation of the leaf vascular network by regulating auxin distribution. Plant Cell Physiol. 2009;50:1319-1328.

117.

Xia

Marquès-Bueno

Bruce

Karnik

Unusual roles of secretory SNARE SYP132 in plasma membrane H⁺-ATPase traffic and vegetative plant growth. Plant Physiol. 2019;180:837-858.

118.

Fonseka

Casson

, et al. PIN protein localization by an R-SNARE-dependent exocytic pathway. BIORXIV/2019/821272.

119.

Short

Leighton

Imriz

, et al. Epidermal expression of a sterol biosynthesis gene regulates root growth by a non-cell autonomous mechanism in Arabidopsis. Development. 2018;145:dev160572.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.33 MB

0.00 MB

0.12 MB

0.07 MB

0.10 MB

4.57 MB

0.24 MB

0.06 MB