Identification and Characterization of the EXO70 Gene Family in Polyploid Wheat and Related Species

Characterization of sucrose nonfermenting-1-related protein kinase 2 (SnRK2) gene family in Haynaldia villosa demonstrated SnRK2.9-V enhances drought and salt stress tolerance of common wheat

BACKGROUND

The sucrose nonfermenting-1-related protein kinase 2 (SnRK2) plays a crucial role in responses to diverse biotic/abiotic stresses. Currently, there are reports on these genes in Haynaldia villosa, a diploid wild relative of wheat.

RESULTS

To understand the evolution of SnRK2-V family genes and their roles in various stress conditions, we performed genome-wide identification of the SnRK2-V gene family in H. villosa. Ten SnRK2-V genes were identified and characterized for their structures, functions and spatial expressions. Analysis of gene exon/intron structure further revealed the presence of evolutionary paths and replication events of SnRK2-V gene family in the H. villosa. In addition, the features of gene structure, the chromosomal location, subcellular localization of the gene family were investigated and the phylogenetic relationship were determined using computational approaches. Analysis of cis-regulatory elements of SnRK2-V gene members revealed their close correlation with different phytohormone signals. The expression profiling revealed that ten SnRK2-V genes expressed at least one tissue (leave, stem, root, or grain), or in response to at least one of the biotic (stripe rust or powdery mildew) or abiotic (drought or salt) stresses. Moreover, SnRK2.9-V was up-regulated in H. villosa under the drought and salt stress and overexpressing of SnRK2.9-V in wheat enhanced drought and salt tolerances via enhancing the genes expression of antioxidant enzymes, revealing a potential value of SnRK2.9-V in wheat improvement for salt tolerance.

CONCLUSION

Our present study provides a basic genome-wide overview of SnRK2-V genes in H. villosa and demonstrates the potential use of SnRK2.9-V in enhancing the drought and salt tolerances in common wheat.

Heterologous expression of rice OsEXO70FX1 confers tolerance to cadmium in Arabidopsis thaliana and fission yeast

Cadmium (Cd) is an environmental toxicant that accumulates in grains, which greatly increases the risk of human exposure to Cd via food chain. The exocytosis of Cd is one of the essential detoxification mechanisms in plants. OsEXO70s, which facilitate the fusion of secretory vesicles and target membranes, has undergone significant expansion in rice. Here, we uncovered 40 OsEXO70 genes characterized by genome-wide profiling and focused on the potential functions of OsEXO70s, especially OsEXO70FX1, in Cd stress. Overexpression of OsEXO70FX1 enhanced both diamide and Cd tolerances in Schizosaccharomyces pombe (S. pombe), and in Arabidopsis resulted in 11% more seedlings survival rate and about 70% longer primary roots under Cd treatment compared with WT (empty vector). Meanwhile, Cd treatment upregulated the expression levels of some exocyst subunits in overexpression lines. Trichomes isolated from overexpression lines were observed to accumulate more Cd. Also, reactive oxygen species (ROS) induced by Cd stress reflected less sensitivity of OsEXO70FX1 overexpression lines to Cd stress, which was evidenced in the Cd determination assay. These results provide the fundament to future research on rice EXO70 family and suggest that it may have evolved a specialized role in response to Cd stress.

The functional and structural characterization of Xanthomonas campestris pv. campestris core effector XopP revealed a new kinase activity

Exo70B1 is a protein subunit of the exocyst complex with a crucial role in a variety of cell mechanisms, including immune responses against pathogens. The calcium-dependent kinase 5 (CPK5) of Arabidopsis thaliana (hereafter Arabidopsis), phosphorylates AtExo70B1 upon functional disruption. We previously reported that, the Xanthomonas campestris pv. campestris effector XopP compromises AtExo70B1, while bypassing the host's hypersensitive response, in a way that is still unclear. Herein we designed an experimental approach, which includes biophysical, biochemical, and molecular assays and is based on structural and functional predictions, utilizing AplhaFold and DALI online servers, respectively, in order to characterize the in vivo XccXopP function. The interaction between AtExo70B1 and XccXopP was found very stable in high temperatures, while AtExo70B1 appeared to be phosphorylated at XccXopP-expressing transgenic Arabidopsis. XccXopP revealed similarities with known mammalian kinases and phosphorylated AtExo70B1 at Ser107, Ser111, Ser248, Thr309, and Thr364. Moreover, XccXopP protected AtExo70B1 from AtCPK5 phosphorylation. Together these findings show that XccXopP is an effector, which not only functions as a novel serine/threonine kinase upon its host target AtExo70B1 but also protects the latter from the innate AtCPK5 phosphorylation, in order to bypass the host's immune responses. Data are available via ProteomeXchange with the identifier PXD041405.

Genome-Wide Analysis of Exocyst Complex Subunit Exo70 Gene Family in Cucumber

Cucumber (Cucumis sativus L.) is an important vegetable worldwide, but its yield is affected by a wide range of pathogens and pests. As the major subunit of the exocyst complex, the roles of Exo70 members have been shown in Arabidopsis and rice, but their function are unknown in cucumber. Here, we identified 18 CsExo70 members in cucumber, which were divided into three groups (Exo70.1-Exo70.3) and nine subgroups (Exo70A-Exo70I) based on the phylogenetic tree. Subsequently, systematical analyses were performed, including collinearity, gene structure, cis-acting elements, conserved motifs, expression patterns, and subcellular localization. Our results showed that CsExo70 genes were generally expressed in all tissues, and CsExo70C1 and CsExo70C2 were highly expressed in the stamen. Moreover, the expression levels of most CsExo70 genes were induced by Pseudomonas syringae pv. lachrymans (Psl) and Fusarium oxysporum f. sp. cucumerinum Owen (Foc), especially CsExo70E2 and CsExo70H3. In addition, these CsExo70s displayed similar location patterns with discrete and punctate signals in the cytoplasm. Together, our results indicate that CsExo70 members may be involved in plant development and resistance, and provide a reference for future in-depth studies of Exo70 genes in cucumber.

Genome-Wide Association Analysis of Senescence-Related Traits in Maize

Senescence is a programmed process that involves the destruction of the photosynthesis apparatus and the relocation of nutrients to the grain. Identifying senescence-associated genes is essential to adapting varieties for the duration of the cultivation cycle. A genome-wide association study (GWAS) was performed using 400 inbred maize lines with 156,164 SNPs to study the genetic architecture of senescence-related traits and their relationship with agronomic traits. We estimated the timing of senescence to be 45 days after anthesis in the whole plant and specifically in the husks. A list of genes identified in a previous RNAseq experiment as involved in senescence (core senescence genes) was used to propose candidate genes in the vicinity of the significant SNPs. Forty-six QTLs of moderate to high effect were found for senescence traits, including specific QTLs for husk senescence. The allele that delayed senescence primarily increased grain yield and moisture. Seven and one significant SNPs were found in the coding and promoter regions of eight core senescence genes, respectively. These genes could be potential candidates for generating a new variation by genome editing for functional analysis and breeding purposes, particularly Zm00001d014796, which could be responsible for a QTL of senescence found in multiple studies.

A Combination of a Genome-Wide Association Study and a Transcriptome Analysis Reveals circRNAs as New Regulators Involved in the Response to Salt Stress in Maize

Salinization seriously threatens the normal growth of maize, especially at the seedling stage. Recent studies have demonstrated that circular RNAs (circRNAs) play vital roles in the regulation of plant stress resistance. Here, we performed a genome-wide association study (GWAS) on the survival rate of 300 maize accessions under a salt stress treatment. A total of 5 trait-associated SNPs and 86 candidate genes were obtained by the GWAS. We performed RNA sequencing for 28 transcriptome libraries derived from 2 maize lines with contrasting salt tolerance under normal and salt treatment conditions. A total of 1217 highly expressed circRNAs were identified, of which 371 were responsive to a salt treatment. Using PCR and Sanger sequencing, we verified the reliability of these differentially expressed circRNAs. An integration of the GWAS and RNA-Seq analyses uncovered two differentially expressed hub genes (Zm00001eb013650 and Zm00001eb198930), which were regulated by four circRNAs. Based on these results, we constructed a regulation model of circRNA/miRNA/mRNA that mediated salt stress tolerance in maize. By conducting hub gene-based association analyses, we detected a favorable haplotype in Zm00001eb198930, which was responsible for high salt tolerance. These results help to clarify the regulatory relationship between circRNAs and their target genes as well as to develop salt-tolerant lines for maize breeding.

Genome-wide association study for grain yield and component traits in bread wheat (Triticum aestivum L.)

Genomic regions governing days to heading (DH), grain filling duration (GFD), grain number per spike (GNPS), grain weight per spike (GWPS), plant height (PH), and grain yield (GY) were investigated in a set of 280 diverse bread wheat genotypes. The genome-wide association studies (GWAS) panel was genotyped using a 35K Axiom Array and phenotyped in five environments. The GWAS analysis showed a total of 27 Bonferroni-corrected marker-trait associations (MTAs) on 15 chromosomes representing all three wheat subgenomes. The GFD showed the highest MTAs (8), followed by GWPS (7), GY (4), GNPS (3), PH (3), and DH (2). Furthermore, 20 MTAs were identified with more than 10% phenotypic variation. A total of five stable MTAs (AX-95024590, AX-94425015, AX-95210025 AX-94539354, and AX-94978133) were identified in more than one environment and associated with the expression of DH, GFD, GNPS, and GY. Similarly, two novel pleiotropic genomic regions with associated MTAs i.e. AX-94978133 (4D) and AX-94539354 (6A) harboring co-localized QTLs governing two or more traits were also identified. In silico analysis revealed that the SNPs were located on important putative candidate genes such as F-box-like domain superfamily, Lateral organ boundaries, LOB, Thioredoxin-like superfamily Glutathione S-transferase, RNA-binding domain superfamily, UDP-glycosyltransferase family, Serine/threonine-protein kinase, Expansin, Patatin, Exocyst complex component Exo70, DUF1618 domain, Protein kinase domain involved in the regulation of grain size, grain number, growth and development, grain filling duration, and abiotic stress tolerance. The identified novel MTAs will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection (MAS).

Expression Analysis of TCP Transcription Factor Family in Autopolyploids of Chrysanthemum nankingense

Autopolyploids often exhibit plant characteristics different from their diploid ancestors and are frequently associated with altered genes expression controlling growth and development. TCP is a unique transcription factor family in plants that is closely related to plant growth and development. Based on transcriptome sequencing of Chrysanthemum nankingense, 23 full-length TCP genes were cloned. The expression of CnTCP9 was most variable in tetraploids, at least threefold greater than diploids. Due to the lack of a C. nankingense transgenic system, we overexpressed CnTCP9 in Arabidopsis thaliana (Col-0) and Chrysanthemum morifolium. Overexpression of CnTCP9 caused enlargement of leaves in A. thaliana and petals in C. morifolium, and the expression of genes downstream of the GA pathway in C. morifolium were increased. Our results suggest that autopolyploidization of C. nankingense led to differential expression of TCP family genes, thereby affecting plant characteristics by the GA pathway. This study improves the understanding of enlarged plant size after autopolyploidization.

Functional analysis of ARF1 from Cymbidium goeringii in IAA response during leaf development

Background

Cymbidium is an economically important genus of flowering orchids cultivated in China because of showing graceful leaf shapes and elegant flower coloration. However, the deterioration of the ecological environment and the difficulty of conservation management have become increasing challenges for maintaining its germplasm resources. ARFs are critical transcription factors in the auxin signaling pathway and have been found to play pivotal roles in leaf growth and development in previous studies. However, their functions and mechanisms in Cymbidium goeringii remain to be clarified.

Methods

The sequence of the CgARF1 gene was analyzed by bioinformatics. The expression of this gene in different tissues and under IAA treatment was detected by quantitative real-time PCR analysis. The CgARF1 gene was overexpressed in wild-type Arabidopsis and Nicotiana benthamiana via the Agrobacterium infection method. Acetone-ethanol solvent extraction was applied for the determination of chlorophyll contents, and the contents of endogenous hormones were determined using the enzyme-linked immunosorbent assay technique.

Results

CgARF1 cloned from C. goeringii 'Songmei' was 2,049 bp in length and encoded 682 amino acids containing three typical domains: a B3 DNA binding domain, an Aux_resp domain and an AUX/IXX family domain. The expression pattern of CgARF1 in different tissues of C. goeringii showed that its expression was highest in the leaves and changed greatly under IAA treatment. Subcellular localization studies showed that the protein was mainly localized in the cell nucleus. CgARF1-overexpressing lines exhibited leaf senescence and a decreased chlorophyll content. Under IAA treatment, CgARF1 regulates the rooting length, rooting number and rooting rate from detached leaves. The levels of endogenous hormones in transgenic leaves were also significantly changed.

Conclusion

These results indicated that CgARF1 overexpression is responsive to IAA treatment during leaf development. This study provides a foundation for future research on the function of the ARF gene family in C. goeringii.

Genome-Wide Identification of GDSL-Type Esterase/Lipase Gene Family in Dasypyrum villosum L. Reveals That DvGELP53 Is Related to BSMV Infection

GDSL-type esterase/lipase proteins (GELPs) characterized by a conserved GDSL motif at their N-terminus belong to the lipid hydrolysis enzyme superfamily. In plants, GELPs play an important role in plant growth, development and stress response. The studies of the identification and characterization of the GELP gene family in Triticeae have not been reported. In this study, 193 DvGELPs were identified in Dasypyrum villosum and classified into 11 groups (clade A–K) by means of phylogenetic analysis. Most DvGELPs contain only one GDSL domain, only four DvGELPs contain other domains besides the GDSL domain. Gene structure analysis indicated 35.2% DvGELP genes have four introns and five exons. In the promoter regions of the identified DvGELPs, we detected 4502 putative cis-elements, which were associated with plant hormones, plant growth, environmental stress and light responsiveness. Expression profiling revealed 36, 44 and 17 DvGELPs were highly expressed in the spike, the root and the grain, respectively. Further investigation of a root-specific expressing GELP, DvGELP53, indicated it was induced by a variety of biotic and abiotic stresses. The knockdown of DvGELP53 inhibited long-distance movement of BSMV in the tissue of D. villosum. This research provides a genome-wide glimpse of the D. villosum GELP genes and hints at the participation of DvGELP53 in the interaction between virus and plants.

Identification and Comprehensive Structural and Functional Analyses of the EXO70 Gene Family in Cotton

The EXO70 gene is a vital component of the exocytosis complex and participates in biological processes ranging from plant cell division to polar growth. There are many EXO70 genes in plants and their functions are extensive, but little is known about the EXO70 gene family in cotton. Here, we analyzed four cotton sequence databases, identified 165 EXO70 genes, and divided them into eight subgroups (EXO70A–EXO70H) based on their phylogenetic relationships. EXO70A had the most exons (≥11), whereas the other seven each had only one or two exons. Hence, EXO70A may have many important functions. The 84 EXO70 genes in Asian and upland cotton were expressed in the roots, stems, leaves, flowers, fibers, and/or ovules. Full-length GhEXO70A1-A cDNA was homologously cloned from upland cotton (Gossypium hirsutum, G. hirsutum). Subcellular analysis revealed that GhEXO70A1-A protein was localized to the plasma membrane. A yeast two-hybrid assay revealed that GhEXO70A1-A interacted with GhEXO84A, GhEXO84B, and GhEXO84C. GhEXO70A1-A silencing significantly altered over 4000 genes and changed several signaling pathways related to metabolism. Thus, the EXO70 gene plays critical roles in the physiological functions of cotton.

Genome-wide identification and characterization of the Lateral Organ Boundaries Domain (LBD) gene family in polyploid wheat and related species

BACKGROUND

Wheat (Triticum aestivum) originated from three different diploid ancestral grass species and experienced two rounds of polyploidization. Exploring how certain wheat gene subfamilies have expanded during the evolutionary process is of great importance. The Lateral Organ Boundaries Domain (LBD) gene family encodes plant-specific transcription factors that share a highly conserved LOB domain and are prime candidates for this, as they are involved in plant growth, development, secondary metabolism and stress in various species.

METHODS

Using a genome-wide analysis of high-quality polyploid wheat and related species genome sequences, a total of 228 LBD members from five Triticeae species were identified, and phylogenetic relationship analysis of LBD members classified them into two main classes (classes I and II) and seven subgroups (classes I a-e, II a and II b).

RESULTS

The gene structure and motif composition analyses revealed that genes that had a closer phylogenetic relationship in the same subgroup also had a similar gene structure. Macrocollinearity and microcollinearity analyses of Triticeae species suggested that some LBD genes from wheat produced gene pairs across subgenomes of chromosomes 4A and 5A and that the complex evolutionary history of TaLBD4B-9 homologs was a combined result of chromosome translocation, polyploidization, gene loss and duplication events. Public RNA-seq data were used to analyze the expression patterns of wheat LBD genes in various tissues, different developmental stages and following abiotic and biotic stresses. Furthermore, qRT-PCR results suggested that some TaLBDs in class II responded to powdery mildew, regulated reproductive growth and were involved in embryo sac development in common wheat.

The Exocyst Complex Subunit EXO70E1-V From Haynaldia villosa Interacts With Wheat Powdery Mildew Resistance Gene CMPG1-V

EXO70 belongs to the exocyst complex subunit that plays a critical role in regulating plant cell polarity establishment and defense response. A previous study proved that the E3 ligase CMPG1-V from Haynaldia villosa, a diploid wheat relative, positively regulates the resistance to wheat powdery mildew (Pm), caused by fungus Blumeria graminis f.sp tritici (Bgt). In this study, a member of EXO70 superfamily named EXO70E1-V was isolated from H. villosa, and EXO70E1-V interacted with CMPG1-V were shown by yeast two-hybrid (Y2H), pull-down assay, bimolecular fluorescence complementation (BiFC) assay, and luciferase complementation imaging (LCI) assay. It is localized in various subcellular organs, i.e., plasma membrane (PM) and endoplasmic reticulum. Co-expression of EXO70E1-V and CMPG1-V showed dot-like structure fluorescence signals that were mainly in PM and nucleus. Expression of EXO70E1-V was relatively higher in leaf and was significantly induced by Bgt infection and exogenous application of hormones such as salicylic acid. Transient or stable overexpression of EXO70E1-V could not enhance/decrease the Pm resistance level, suggesting overexpression of EXO70E1-V alone has no impact on Pm resistance in wheat.

Genome-wide identification and expression analysis of the EXO70 gene family in grape (Vitis vinifera L)

EXO70 is the pivotal protein subunit of exocyst, which has a very crucial role in enhancing the shielding effect of the cell wall, resisting abiotic and hormonal stresses. This experiment aims to identify family members of the EXO70 gene family in grape and predict the characteristics of this gene family, so as to lay the foundation of further exploring the mechanism of resisting abiotic and hormone stresses of VvEXO70s. Therefore, the Vitis vinifera ‘Red Globe’ tube plantlet were used as materials. Bioinformatics was used to inquire VvEXO70 genes family members, gene structure, system evolution, cis-acting elements, subcellular and chromosomal localization, collinearity, selective pressure, codon bias and tissue expression. All of VvEXO70s had the conserved pfam03081 domain which maybe necessary for interacting with other proteins. Microarray analysis suggested that most genes expressed to varying degrees in tendrils, leaves, seeds, buds, roots and stems. Quantitative Real-Time PCR (qRT-PCR) showed that the expression levels of all genes with 5 mM salicylic acid (SA), 0.1 mM methy jasmonate (MeJA), 20% PEG6000 and 4 °C for 24 h were higher than for 12 h. With 20% PEG6000 treatment about 24 h, the relative expression of VvEXO70-02 was significantly up-regulated and 361 times higher than CK. All genes’ relative expression was higher at 12 h than that at 24 h after treatment with 7 mM hydrogen peroxide (H₂O₂) and 0.1 mM ethylene (ETH). In conclusion, the expression levels of 14 VvEXO70 genes are distinguishing under these treatments, which play an important role in the regulation of anti-stress signals in grape. All of these test results provide a reference for the future research on the potential function analysis and plant breeding of VvEXO70 genes.

Identification of QTLs/Defense Genes Effective at Seedling Stage Against Prevailing Races of Wheat Stripe Rust in India

Resistance in modern wheat cultivars for stripe rust is not long lasting due to the narrow genetic base and periodical evolution of new pathogenic races. Though nearly 83 Yr genes conferring resistance to stripe rust have been cataloged so far, few of them have been mapped and utilized in breeding programs. Characterization of wheat germplasm for novel sources of resistance and their incorporation into elite cultivars is required to achieve durable resistance and thus to minimize the yield losses. Here, a genome-wide association study (GWAS) was performed on a set of 391 germplasm lines with the aim to identify quantitative trait loci (QTL) using 35K Axiom® array. Phenotypic evaluation disease severity against four stripe rust pathotypes, i.e., 46S119, 110S119, 238S119, and 47S103 (T) at the seedling stage in a greenhouse providing optimal conditions was carried out consecutively for 2 years (2018 and 2019 winter season). We identified, a total of 17 promising QTl which passed FDR criteria. Moreover these 17 QTL identified in the current study were mapped at different genomic locations i.e. 1B, 2A, 2B, 2D, 3A, 3B, 3D, 4B, 5B and 6B. These 17 QTLs identified in the present study might play a key role in marker-assisted breeding for developing stripe rust resistant wheat cultivars.

Comprehensive analysis of polygalacturonase genes offers new insights into their origin and functional evolution in land plants

Polygalacturonase (PG) is a hydrolase that participates in pectin degradation, pod shattering and fruit softening. Here, we identified 2786 PG genes across 54 plants, which could be divided into three groups. Evolutionary analysis suggested that PG family originated from the charophyte green algae, and Subgroups A2-A4 evolved from the Subgroup A1 after the tracheophyte-angiosperm split. Whole-genome duplication was the major force leading to PG gene expansion. Interestingly, the PG genes continuously expanded in eudicots, whereas it contracted in monocots after the eudicot-monocot split. PG genes in Group A are expressed at high levels in floral organs, whereas genes in Groups B and C are expressed at high levels in various tissues. Moreover, three BnaPG15 members were found for their potential possibility in pod shattering in Brassica napus. Our results provide new insight into the evolutionary history of PG family, and their potentially functional role in plants.

OsExo70B1 Positively Regulates Disease Resistance to Magnaporthe oryzae in Rice

The exocyst, an evolutionarily conserved octameric protein complex, mediates tethering of vesicles to the plasma membrane in the early stage of exocytosis. Arabidopsis Exo70, a subunit of the exocyst complex, has been found to be involved in plant immunity. Here, we characterize the function of OsExo70B1 in rice. OsExo70B1 mainly expresses in leaf and shoot and its expression is induced by pathogen-associated molecular patterns (PAMPs) and rice blast fungus Magnaporthe oryzae (M. oryzae). Knocking out OsExo70B1 results in significantly decreased resistance and defense responses to M. oryzae compared to the wild type, including more disease lesions and enhanced fungal growth, downregulated expression of pathogenesis-related (PR) genes, and decreased reactive oxygen species accumulation. In contrast, the exo70B1 mutant does not show any defects in growth and development. Furthermore, OsExo70B1 can interact with the receptor-like kinase OsCERK1, an essential component for chitin reception in rice. Taken together, our data demonstrate that OsExo70B1 functions as an important regulator in rice immunity.

Characterization of the Heavy-Metal-Associated Isoprenylated Plant Protein (HIPP) Gene Family from Triticeae Species

Heavy-metal-associated (HMA) isoprenylated plant proteins (HIPPs) only exist in vascular plants. They play important roles in responses to biotic/abiotic stresses, heavy-metal homeostasis, and detoxification. However, research on the distribution, diversification, and function of HIPPs in Triticeae species is limited. In this study, a total of 278 HIPPs were identified from a database from five Triticeae species, and 13 were cloned from Haynaldia villosa. These genes were classified into five groups by phylogenetic analysis. Most HIPPs had one HMA domain, while 51 from Clade I had two, and all HIPPs had good collinear relationships between species or subgenomes. In silico expression profiling revealed that 44 of the 114 wheat HIPPs were dominantly expressed in roots, 43 were upregulated under biotic stresses, and 29 were upregulated upon drought or heat treatment. Subcellular localization analysis of the cloned HIPPs from H. villosa showed that they were expressed on the plasma membrane. HIPP1-V was upregulated in H. villosa after Cd treatment, and transgenic wheat plants overexpressing HIPP1-V showed enhanced Cd tolerance, as shown by the recovery of seed-germination and root-growth inhibition by supplementary Cd. This research provides a genome-wide overview of the Triticeae HIPP genes and proved that HIPP1-V positively regulates Cd tolerance in common wheat.

Three subfamilies of exocyst EXO70 family subunits in land plants: early divergence and ongoing functional specialization

Localized delivery of plasma membrane and cell wall components is an essential process in all plant cells. The vesicle-tethering complex, the exocyst, an ancient eukaryotic hetero-octameric protein cellular module, assists in targeted delivery of exocytosis vesicles to specific plasma membrane domains. Analyses of Arabidopsis and later other land plant genomes led to the surprising prediction of multiple putative EXO70 exocyst subunit paralogues. All land plant EXO70 exocyst subunits (including those of Bryophytes) form three distinct subfamilies-EXO70.1, EXO70.2, and EXO70.3. Interestingly, while the basal well-conserved EXO70.1 subfamily consists of multiexon genes, the remaining two subfamilies contain mostly single exon genes. Published analyses as well as public transcriptomic and proteomic data clearly indicate that most cell types in plants express and also use several different EXO70 isoforms. Here we sum up recent advances in the characterization of the members of the family of plant EXO70 exocyst subunits and present evidence that members of the EXO70.2 subfamily are often recruited to non-canonical functions in plant membrane trafficking pathways. Engagement of the most evolutionarily dynamic EXO70.2 subfamily of EXO70s in biotic interactions and defence correlates well with massive proliferation and conservation of new protein variants in this subfamily.