The Arabidopsis V-ATPase is localized to the TGN/EE via a seed plant-specific motif
Abstract
The V-ATPase is a versatile proton-pump found in a range of endomembrane compartments yet the mechanisms governing its differential targeting remain to be determined. In Arabidopsis, VHA-a1 targets the V-ATPase to the TGN/EE whereas VHA-a2 and VHA-a3 are localized to the tonoplast. We report here that the VHA-a1 targeting domain serves as both an ER-exit and as a TGN/EE-retention motif and is conserved among seed plants. In contrast, Marchantia encodes a single VHA-isoform that localizes to the TGN/EE and the tonoplast in Arabidopsis. Analysis of CRISPR/Cas9 generated null alleles revealed that VHA-a1 has an essential function for male gametophyte development but acts redundantly with the tonoplast isoforms during vegetative growth. We propose that in the absence of VHA-a1, VHA-a3 is partially re-routed to the TGN/EE. Our findings contribute to understanding the evolutionary origin of V-ATPase targeting and provide a striking example that differential localization does not preclude functional redundancy.
Introduction
Compartmentalization into distinct membrane-bound organelles that offer different chemical environments optimized for the biological processes that occur within them is a hallmark of eukaryotic cells. Vacuolar-type H+-ATPases (V-ATPases) are rotary engines that couple the energy released by ATP hydrolysis to the transport of protons across membranes of several intracellular compartments. They consist of two subcomplexes: The cytosolic V1 subcomplex responsible for ATP hydrolysis composed of eight different subunits (A, B, C, D, E, F, G, H) and the membrane-integral VO subcomplex consisting of 6 subunits (a, d, e, c, c’, c’) required for proton translocation (Sze et al., 2002). Although all eukaryotic V-ATPases are strikingly similar regarding their structure and biochemical activity, their biological reach has been greatly diversified by cell-type specific expression and differential subcellular localization (Cotter et al., 2015). Differential targeting of the V-ATPase is mediated by isoforms of subunit a, the largest of the V-ATPase subunits that consists of a C- terminal hydrophobic domain with eight transmembrane domains and a large N-terminal domain that is accessible for cytosolic interaction partners (Zhao et al., 2015; Vasanthakumar et al., 2019). Previous studies in yeast and plants have shown that the targeting information is contained in the N-terminal domain (Dettmer et al., 2006; Kawasaki-Nishi et al., 2001a), however the responsible targeting domain has so far only been addressed for the yeast isoforms Vhp1p and Stv1p. Whereas Vph1p targets the yeast V-ATPase to the vacuole, a motif containing an aromatic residue (WKY) in the N-terminus of Stv1p targets the V-ATPase to the Golgi/endosomal network (Finnigan et al., 2012). The localization of Stv1p at the Golgi is dependent on interaction of the WKY motif with phosphatidylinositol-4-phosphate (PI(4)P; Banerjee and Kane, 2017). Mammals possess four subunit a isoforms (a1, a2, a3, and a4) which steer the V-ATPase to different endomembranes or the plasma membrane (Forgac, 2007; Marshansky et al., 2014; Futai et al., 2019). In Arabidopsis VHA-a2 and VHA-a3 target the V-ATPase to the tonoplast (Dettmer et al., 2006) where the combined action of V-ATPase and V-PPase energizes secondary active transport and maintains the acidic environment required for the lytic function of plant vacuoles (Sze et al., 1999; Kriegel et al., 2015). In contrast, VHA-a1 targets the V-ATPase to the TGN/EE, a highly dynamic organelle that receives and sorts proteins from the endocytic, recycling and secretory pathways (Dettmer et al., 2006; Viotti et al., 2010). The combined action of the V-ATPase and proton-coupled antiporters including the NHX-type cation proton exchangers (Bassil et al., 2011a and Bassil et al., 2011b; Dragwidge et al., 2019) and the ClC (Chloride channel) Cl-/NO- proton antiporters (von der Fecht-Bartenbach et al., 2007) is responsible for generating and maintaining the acidic pH of the TGN/EE (Luo et al., 2015). Genetic and pharmacological inhibition of the V-ATPase interferes with endocytic and secretory trafficking (Dettmer et al., 2006; Viotti et al., 2010; Luo et al., 2015) and causes defects in both cell division and cell expansion (Dettmer et al., 2006; Brüx et al., 2008). Despite the TGN/EE being the central hub for protein sorting, it is still not clear how the identity of this compartment is specified and how resident proteins required for functionality are maintained while at the same time very similar proteins are sorted and leave the TGN/EE as cargo.
We have shown previously that after assembly involving dedicated ER-chaperones (Neubert et al., 2008), VHA-a1 and VHA-a3 containing Vo-subcomplexes leave the ER via different trafficking routes. Whereas VHA-a1 containing V-ATPases are exported in a coat protein complex II (COPII)-dependent manner (Viotti et al., 2013), VHA-a3 containing complexes have been shown to be delivered in a Golgi-independent manner to the tonoplast (Viotti et al., 2013). It remains to be addressed how V-ATPases containing the different VHA-a isoforms are sorted in the ER and how they are targeted to their final destinations in the cell. Here, we focus on VHA-a1 to understand how the V-ATPase is targeted to and retained in the TGN/EE. Targeting of TGN-resident proteins has been shown to be mediated by groups of acidic amino acids (acidic clusters), di-leucine and tyrosine-based motifs in mammals (Bos et al., 1993; Schäfer et al., 1995; Alconada et al., 1996; Xiang et al., 2000). In yeast, retrieval to the late Golgi (equivalent to the TGN) is dependent on the concerted action of aromatic-based amino acid (aa) motifs in the cytoplasmic tails of proteins and slow anterograde transport to the late endosome (Wilcox et al., 1992; Nothwehr et al., 1993; Cooper and Stevens, 1996; Cereghino et al., 1995; Bryant and Stevens, 1997). We have previously shown that the TGN/EE targeting information is contained within the first 228 aa of VHA-a1 (Dettmer et al., 2006). Here, by using chimeric proteins, 3D homology modeling, site-directed mutagenesis and live cell imaging, we identified a region that is required for ER-export as well as TGN/EE-retention. Although most plant genomes encode multiple VHA-a isoforms, the VHA-a1 targeting domain (a1-TD) is conserved only among seed plants. The liverwort Marchantia polymorpha with a single gene encoding subunit a (MpVHA-a) that is mostly present at the tonoplast in both Marchantia and Arabidopsis might represent the ancestral state in which acidification of the TGN/EE was achieved by transitory V-ATPase complexes. The identification of homozygous mutants lacking VHA-a1 indicates that the VHA-a isoforms, despite their differential localization, act redundantly. Based on the results presented here, we propose a model in which differential targeting of the VHA-a isoforms is based on competition for entry into COPII vesicles.
Results
Identification of a region in the VHA-a1 N-terminus (a1-TD) that is necessary and sufficient for TGN/EE localization
The tri-peptide motif that mediates Golgi-localization of Stv1p in yeast (Finnigan et al., 2012) is not conserved in VHA-a1 (Figure 1—figure supplement 1) implying that TGN/EE localization evolved independently. We thus used chimeric proteins consisting of increasing lengths of the cytosolic N-terminal domain of VHA-a1 (a1NT 37, 85, 131, 179 and 228 aa) fused to decreasing lengths of the C-terminal domain of VHA-a3, to further narrow down the region required for TGN/EE targeting. Constructs encoding the chimeric proteins fused to GFP were expressed in Arabidopsis under control of the UBIQUITIN10 promoter (UBQ10; Grefen et al., 2010). Whereas a1NT37a3-GFP, a1NT85a3-GFP (Figure 1—figure supplement 2) and a1NT131a3-GFP (Figure 1A) all localized to the tonoplast, a1NT179a3-GFP and a1NT228a3-GFP were detectable at the tonoplast but also in a punctate pattern reminiscent of the TGN/EE (Figure 1A). TGN/EE localization was confirmed by colocalization with the endocytic tracer FM4-64 (Dettmer et al., 2006) and treatment with the fungal drug Brefeldin A (BFA) that causes the aggregation of the TGN/EE to form BFA compartments (Nebenführ et al., 2002; Dettmer et al., 2006). After 3 hr of BFA treatment, the core of BFA compartments was labeled by a1NT179a3-GFP and a1NT228a3-GFP (Figure 1B). From these observations, we concluded that the region necessary for TGN/EE localization of VHA-a1 (a1-targeting domain; a1-TD) is located between residues L132 and E179.
To visualize the a1-TD in the three-dimensional (3D) structure, we used I-TASSER (Roy et al., 2010) to generate 3D models of the VHA-a1 and VHA-a3 N-termini with the atomic models of the yeast isoforms Stv1p and Vph1p (Vasanthakumar et al., 2019) serving as templates, respectively (Figure 1C, Supplementary file 1 and Figure 1—figure supplement 3). Superimposition of the obtained models for the VHA-a1 and VHA-a3 N-termini revealed that the a1-TD is one of the regions predicted to be different (Figure 1C). To visualize the orientation of the VHA-a1 N-terminus within the V-ATPase complex, we aligned our model for the VHA-a1 N-terminus model with the structure of the yeast V-ATPase (PDB6O7V; Vasanthakumar et al., 2019) and found that the a1-TD would be accessible for recognition in the fully assembled V-ATPase complex (Figure 1D).
To test if the VHA-a1-TD is sufficient to target VHA‐a3 to the TGN/EE, we introduced the 34 aa region between K140 and S174 into the N-terminus of VHA-a3. The resulting GFP fusion construct (UBQ10:VHA-a3-a1-TD-GFP) was transformed into wildtype plants and VHA-a3-a1-TD-GFP was detectable at the tonoplast and at the TGN/EE where it colocalized with VHA-a1-mRFP (Figure 1E). It remains to be determined why TGN/EE localization of VHA-a3-a1-TD is only partial, however it demonstrates that the a1-TD is sufficient to partially re-route VHA-a3 to the TGN/EE.
Mutation of the a1-TD in VHA-a1 leads to mislocalization to the tonoplast
Phylogenetic analysis of VHA-a sequences from mono- and dicot species showed that VHA-a1 and VHA-a3 represent two clearly separated clades (Figure 2—figure supplement 1). The a1-TD includes a number of amino acids conserved among all of the sequences in the VHA-a1 clade (Figure 2A) but not in the VHA-a3 clade that could either be part of a di-acidic ER-export motif or an acidic cluster involved in TGN/EE-retention. Next, we used site-directed mutagenesis to identify which amino acids are required for correct localization of VHA-a1. Conserved amino acids in the a1-TD were exchanged individually or in combinations of two or three to the corresponding amino acid of VHA-a3. To assess the effect that each mutation had on the localization of VHA-a1-GFP, we performed colocalization experiments with VHA-a1-mRFP. Unmutated VHA-a1-GFP coexpressed with VHA-a1-mRFP was used as a control. Nine to ten images were subjected to quantitative image analysis. The overall degree of colocaliaztion was measured using the standard Pearson’s correlation coefficient (Manders et al., 1992). In addition, to determine the fraction of VHA-a1-GFP in the VHA-a1-mRFP location, the Manders’ Colocalization Coefficients (Manders et al., 1993) were also calculated for each individual image.
Mutations with high Pearson’s correlation coefficients predominately localized to the TGN/EE (E161S and F134Y; Figure 2C,D and K). Three mutations led to a dual TGN/EE and tonoplast localization (L159T + E161S, L159T, E155 + E156 + I157 deletion (ΔEEI) and Figure 2E,F,G and H) and had intermediate Pearson’s and Mander’s coefficients (M2) (Figure 2K). Mutation E156Q + L159T + E161S (ELE) also led to a dual TGN/EE and tonoplast localization but displayed fewer TGN/EEs therefore had a lower Mander’s coefficient compared to the other mutations that had a dual TGN/EE and tonoplast localization. Two mutations led to VHA-a1 being localized to an unknown compartment (E156Q and E156Q + L159T; Figure 2I and J) and had the lowest Pearson’s and Mander’s coefficients (Figure 2K).
The a1-TD is required for COPII-mediated ER-export
As we have shown previously that VHA-a3 containing complexes are delivered in a Golgi-independent manner from the ER to the tonoplast (Viotti et al., 2013), we next asked if the tonoplast signal observed for the mutated VHA-a1 proteins was due to poor retention at the TGN/EE or partial entry into the provacuolar route caused by reduced ER-exit.
To address if the a1-TD contains an ER-export motif, we made use of a dominant negative mutation of Sar1 (Sar1-GTP) that has been shown to block COPII-mediated ER-export (daSilva et al., 2004). We expressed AtSar1b-GTP-CFP under a dexamethasone (DEX) inducible promoter (pUBQ10 >GR > AtSar1b-GTP-CFP; Moore et al., 1998) and used electron microscopy (EM) to confirm that COPII-mediated ER-export was indeed blocked. Upon induction of AtSar1b-GTP-CFP, we observed large clusters of vesicles in the periphery of Golgi stack remnants and ER cisternae appeared swollen (Figure 3—figure supplement 1). Furthermore, in Arabidopsis transgenic lines expressing the Golgi-marker Sialyl transferase (ST; Boevink et al., 1998) or the brassinosteroid receptor (BRI1; Geldner et al., 2007) as GFP fusion proteins, induction of AtSar1b-GTP-CFP caused their ER-retention and colocalization with AtSar1b-GTP-CFP in bright and often large punctae that might represent the clusters of uncoated COPII vesicles as observed by EM (Figure 3—figure supplement 1).
Similarly, the punctate pattern of VHA-a1-GFP disappeared after induction of AtSar1b-GTP-CFP and was replaced by a characteristic ER pattern with dense punctae that colocalized with AtSar1b-GTP-CFP (Figure 3A). In contrast, VHA-a3-GFP was not retained in the ER after the induction of AtSar1b-GTP-CFP (Figure 3B). To exclude that these differences in ER retention between VHA-a1-GFP and VHA-a3-GFP were due to differences in induction strength of AtSar1b-GTP-CFP, we also blocked ER exit in a line co-expressing VHA-a1-GFP and VHA-a3-mRFP. In cells in which VHA-a1-GFP was retained in the ER, VHA-a3-mRFP was not affected (Figure 3C). Conversely, VHA-a3-a1-TD-GFP is partially retained in the ER upon expression of AtSar1b-GTP-CFP (Figure 3D). Taken together, our data confirms that the ER-export of VHA-a1 is COPII-dependent and reveals that the a1-TD contains an ER-export signal.
When AtSar1b-GTP was used to block ER exit, we observed a significant increase in the GFP fluorescence intensity at the tonoplast for all the mutated VHA-a1 proteins (Figure 4A and B) indicating that blocking ER-export causes an increase in the trafficking to the tonoplast via a Golgi-independent route.
Based on our observation that VHA-a3 carrying a mutation (VHA-a3-R729N-GFP) that renders the complex inactive is retained in the ER when expressed in wildtype but found at the tonoplast in the vha-a2 vha-a3 mutant (Figure 5—figure supplement 1; Krebs et al., 2010), we next tested if tonoplast localization of proteins carrying mutations in the a1-TD increases in the absence of competing complexes containing VHA-a2 and VHA-a3. Indeed, when the mutated VHA-a1 proteins were expressed in the vha-a2 vha-a3 mutant, the ratio of TGN/EE-to-tonoplast fluorescence intensity increased as compared to the wildtype background (Figure 5A and B). Similarly, VHA-a1 mutants that showed only TGN localization (E161S) or a punctate pattern that is different from VHA-a1-RFP (E156Q and E156Q+ L159T) in the wildtype background displayed a tonoplast signal in the vha-a2 vha-a3 double mutant (Figure 5A and B). These results strongly support the idea that all mutations affect ER export and that the observed tonoplast signal is not a result of poor retention at the TGN/EE but due to a Golgi-independent pathway to the tonoplast.
The a1-TD is conserved in the plant kingdom and originates with the gymnosperms
To trace the evolutionary origin of the a1-TD in the plant kingdom, VHA-a sequences from selected species were subjected to phylogenetic analysis. Our analysis revealed that all VHA-a isoforms from seed plants including the gymnosperm Pinus taeda cluster into two distinct clades containing VHA-a1 and VHA-a3 respectively (Figure 6—figure supplement 1 and Figure 6A). Interestingly, the genomes of most basal plants encode multiple isoforms of VHA-a, however they do not fall into either of these two clades implying that duplication of VHA-a occurred independently. Closer analysis of the a1-TD sequence revealed that the hallmarks of the a1-TD (acidic residues flanking a critical leucine residue) are conserved throughout the angiosperms but are absent in bryophytes, lycophytes, ferns and hornworts. Interestingly, the gymnosperm sequence from P. taeda contains the a1-TD sequence with the exception of one flanking glutamic residue (Figure 6B). To test if the a1-TD is functionally conserved we generated chimeric proteins consisting of the N-terminal domain of VHA-a from Pinus taeda and Amborella trichopoda as representatives of the seed plants and Selaginella moellendorffii as a non-seed plant fused to the C-terminal domain of VHA-a1 (Figure 6C). Whereas the Pine and Amborella chimeric proteins colocalized with VHA-a1 at the TGN/EE, the Selaginella chimeric protein localized to the tonoplast. These results suggest that the a1-TD domain is conserved within the spermatophytes but not in other green plants.
Mislocalized VHA-a1 can replace the tonoplast V-ATPase
Mislocalization of VHA-a1 at the tonoplast provides the opportunity to address if the isoforms belonging to the VHA-a1 and VHA-a3 clades are functionally divergent. We used the ability to complement the dwarfed phenotype of the vha-a2 vha-a3 double mutant (Krebs et al., 2010) in standard long day conditions (Figure 7—figure supplement 1) as well as in short day conditions (Figure 7A) as a proxy. We confirmed that the proteins were expressed (Figure 7—figure supplement 2) and observed that complementation of the growth phenotype correlated with the intensity of the tonoplast signal with VHA-a1E161S-GFP conferring the lowest and VHA-a1ΔEEI-GFP the highest degree of rescue (Figure 7A and B). The ability of V-ATPases with incorporated mutated VHA-a1-GFP subunits to complement the vha-a2 vha-a3 double mutant rosette phenotype is an indication that these complexes localize to the tonoplast and are functional in all plant tissues and not only in root cells. Three VHA-a1-GFP mutation lines showing high (ΔEEI), intermediate (E156Q + L159T + E161S) and low degrees (L159T + E161S) of rescue were selected for further analysis. Cell sap pH measurements were performed as a proxy of vacuolar pH. Whereas all VHA-a1-GFP mutation lines had more acidic vacuoles than the vha-a2 vha-a3 mutant none of them reached wildtype vacuolar pH (Figure 7C). Taken together, our data indicates that although VHA-a1 containing complexes are functional when mistargeted to the tonoplast their enzymatic properties are different and they cannot fully replace the vacuolar isoforms.
The Marchantia V-ATPase is dual localized at the TGN/EE and tonoplast and is functional at the tonoplast in Arabidopsis
Most plants outside of the spermatophytes possess multiple VHA-a isoforms with the exception of the liverwort Marchantia polymorpha and the fern Salvinia cucullata (Supplementary file 1). Focusing on Marchantia, we addressed if the TGN/EE or the tonoplast V-ATPase represent the ancestral state or if dual localization of the V-ATPase can also be achieved in the absence of differentially localized isoforms. To examine the subcellular localization of the single, M. polymorpha VHA-a protein (MpVHA-a), mVenus-tagged MpVHA-a (MpVHA-a-mVenus) was co-expressed with the TGN marker, mRFP-MpSYP6A (Kanazawa et al., 2016) in M. polymorpha thallus cells. MpVHA-a-mVenus driven by the CaMV35S or MpEF1α promoter predominantly localized to the vacuolar membrane, but some additional punctate structures were detectable (Figure 8—figure supplement 1).
Next, we expressed a MpVHA-a-mVenus fusion construct (UBQ10:MpVHA-a-mVenus) in the Arabidopsis wildtype background. Whereas MpVHA-a-mVenus localization at the tonoplast of root cells was clearly visible (Figure 8—figure supplement 2), a clear punctate pattern was not observed. However, after BFA treatment, the core of BFA compartments was labeled with MpVHA-a-mVenus (Figure 8—figure supplement 2) indicating that MpVHA-a is also present at the TGN/EE. To test for functionality, MpVHA-a-mVenus was expressed in the vha-a2 vha-a3 mutant background. MpVHA-a-mVenus localized to the TGN/EE and tonoplast and was found in the core of FM4-64 labeled BFA compartments (Figure 8A). A growth assay conducted in short day conditions (SD; 22 °C and 10 hr light) revealed that V-ATPases that incorporate MpVHA-a can complement the vha-a2 vha-a3 double mutant to wildtype levels in terms of rosette size and cell sap pH (Figure 8B and C). We also confirmed that the MpVHA-a-mVenus proteins are expressed (Figure 8D).
VHA-a1 has a unique and essential function during pollen development which cannot be fulfilled by VHA-a2, VHA-a3 or MpVHA-a
Next, we wanted to test if MpVHA-a can also replace VHA-a1 at the TGN/EE. Null alleles of single copy-encoded VHA-subunits cause male gametophyte lethality (Dettmer et al., 2005) and it has thus not been possible to identify homozygous T-DNA mutants for VHA-a1. To avoid T-DNA related silencing problems in complementation experiments with heterozygous vha-a1/+ mutants, we used CRISPR/Cas9 under control of an egg cell-specific promoter to generate mutant alleles (Wang et al., 2015). Different regions of VHA-a1 were targeted using four different guide RNAs (gRNAs 1–4, Figure 9A). Wildtype plants as well as VHA-a1-GFP expressing transgenic plants were transformed and T1 plants were analyzed by sequencing of PCR products spanning the CRISPR site. Four alleles were selected for further studies: vha-a1-1 containing a 260 bp deletion which eliminates the start codon and vha-a1-2,−3 and −4 that each contain single base pair insertions leading to frameshifts and early stop codons (Figure 9A).
Surprisingly, we did not only identify heterozygous vha-a1/+ but also homozygous and bi-allelic vha-a1 individuals (Supplementary file 1). The latter two were indistinguishable from wildtype during vegetative growth (Figure 9B), but failed to produce seeds due to a defect in pollen development (Figure 9C). To confirm that the defect in pollen development is indeed caused by a lack of VHA-a1, we analyzed if it is rescued by VHA-a1-GFP. Although VHA-a1:VHA-a1-GFP (Dettmer et al., 2006) is also targeted by the chosen gRNAs and UBQ10:VHA-a1-GFP is targeted by gRNAs 2,3 and 4, the use of appropriate primer combinations allowed us to distinguish between mutations in the endogenous locus and the transgene (Figure 9—figure supplement 1). Plants that carried mutations corresponding to vha-a1-2 in the transgene did not express VHA-a1-GFP confirming that it is indeed a null allele (Figure 9—figure supplement 1). Importantly, all plants with no wildtype allele of either VHA-a1 or VHA-a1-GFP were defective in pollen development, whereas vha-a1 mutants expressing VHA-a1-GFP under the control of UBQ10 promoter showed normal pollen development (Figure 9C).
We next established stable lines expressing UBQ10:VHA-a1-GFP in the vha-a1-1 background and performed reciprocal crosses to determine if transmission via the female gametophyte is also affected. As expected, vha-a1-1 was only transmitted via the male gametophyte in the presence of VHA-a1-GFP. In contrast, transmission of vha-a1-1 via the female gametophyte did not require the presence of the transgene indicating that development of the female gametophyte is not affected (Supplementary file 1).
As homozygous vha-a1 mutants cannot be transformed, rescue experiments would require crosses with a pollen donor that carries a wildtype allele of VHA-a1 so that complementation could only be determined in the F2. For crosses we thus used vha–a1 mutants that still contained the CRISPR T-DNA (Cas9+) so that the incoming wildtype allele would be mutated and homozygous vha-a1 mutants could be obtained in the resulting F1. Using this strategy, we identified bi–allelic and homozygous vha-a1 mutants expressing UBQ10:MpVHA-a-mVenus and found that it was localized both at the TGN/EE and the tonoplast (Figure 9—figure supplement 2). However, MpVHA-a-mVenus did not rescue the pollen phenotype of vha-a1 (Figure 9—figure supplement 2) indicating that the TGN/EE-localized V-ATPase of seed plants has acquired a unique function during pollen development.
During vegetative growth V-ATPases containing VHA-a2 and VHA-a3 compensate for the lack of VHA-a1
Given the importance of acidification of the TGN/EE for endomembrane trafficking and the fact that we have previously reported that RNA-mediated knock-down of VHA-a1 causes reduced cell expansion (Brüx et al., 2008), the observation that the vha-a1 mutant displays no different phenotypes from wildtype during vegetative growth is unexpected. However, vha-a1 is hypersensitive to the V-ATPase inhibitor Concanamycin-A (ConcA). Analysis of the progeny of vha-a1 VHA-a1-GFP (Cas9+) mutants revealed that vha-a1 vha-a1-GFP etiolated seedlings displayed shorter hypocotyls and had deformed roots in contrast to vha-a1 mutants expressing VHA-a1-GFP. This result indicates that a target of ConcA is present at the TGN/EE in vha-a1 albeit in small amounts or with a higher sensitivity to ConcA. (Figure 10A and Figure 10—figure supplement 1). To test if VHA-a2 or VHA-a3-containing V-ATPases might compensate for the lack of VHA-a1, vha-a1-1 (Cas9+) was crossed with the vha-a2 vha-a3 double mutant. vha-a1 vha-a2/+ vha-a3/+ mutants obtained from this cross were significantly smaller than vha-a2/+ vha-a3/+ individuals indicating that tonoplast V-ATPases indeed compensate for the lack of VHA-a1 during vegetative growth (Figure 10B and Figure 10—figure supplement 2).
Subsequently, vha-a1 (Cas9+) was crossed with vha-a2 and vha-a3 single mutants to determine if both are able to compensate for the lack of VHA-a1. vha-a1 vha-a2/+ and vha-a1 vha-a3/+ plants both had reduced rosette sizes, however the latter showed a stronger reduction which is in accordance with VHA-a3 being expressed at higher levels than VHA-a2 (Figure 10—figure supplement 2). To analyze mutants with even lower numbers of VHA-a wildtype alleles vha-a1 (Cas9+) vha-a2/+ vha-a3/+ was crossed with the vha-a2 vha-a3 double mutant. Consistent with the results of the single mutant crosses, vha-a1/+ vha-a2/+ vha-a3 was smaller than vha-a1/+ vha-a2 vha-a3/+. The smallest plants obtained were vha-a1 vha-a2 vha-a3/+ (Figure 10C). Neither did we find vha-a1 vha-a2/+ vha-a3 nor the homozygous triple mutant vha-a1 vha-a2 vha-a3, suggesting that they are not viable.
Discussion
The information needed to deliver and keep the V-ATPase at the TGN/EE is contained in the a1-TD
The compartments of the eukaryotic endomembrane system are acidified to varying degrees by the activity of the vacuolar H+-ATPase (V-ATPase). Acidification enables specific biochemical reactions as well as secondary active transport and can thus be considered a central component of compartmentation. Throughout eukaryotes, isoforms of the membrane-integral VO-subunit a are used to achieve differential targeting of the V-ATPase and although we know very little about the underlying mechanisms, it is clear that differential targeting has arisen independently in different groups. Whereas the targeting information is contained in the N-terminus of the Golgi-localized yeast isoform Stv1p (Finnigan et al., 2012), it was shown to be provided by the C-terminal part for some of the 17 isoforms of Paramecium tetraurelia that localize the V-ATPase to at least seven different compartments (Wassmer et al., 2006). Moreover, the targeting motif identified in Stv1p (Finnigan et al., 2012) is neither conserved in plant nor mammalian subunit a isoforms implying that differential targeting mediated by subunit a evolved independently in the different eukaryotic lineages. Given the essential role that the V-ATPase plays in endocytic and secretory trafficking (Dettmer et al., 2006; Luo et al., 2015), we set out to identify the mechanism underlying its highly specific localization at the TGN/EE (Dettmer et al., 2006). Based on our previous studies, in which we had shown that the targeting information is located within the 228 amino acids of VHA-a1 (Dettmer et al., 2006), we first narrowed down this region to a 30 aa region (a1-TD) which we show here to be both necessary and sufficient for TGN/EE localization. Structural modeling revealed that the a1-TD that contains several highly conserved acidic amino acids is exposed and accessible for interaction partners. Acidic clusters have been shown to be involved in TGN localization (Schäfer et al., 1995; Alconada et al., 1996; Xiang et al., 2000) and we thus used site-directed mutagenesis and found that mutation of conserved aa in the a1-TD indeed leads to partial mislocalization of VHA-a1 to the tonoplast as expected if retention at the TGN/EE would be affected. However, the acidic residues in the a1-TD could also represent di-acidic ER-export motifs and as we have shown previously that VHA-a3 is targeted to the tonoplast in a Golgi-independent manner (Viotti et al., 2013), we needed to differentiate between the a1-TD serving as an ER exit motif and/or a TGN-retention motif. We used inducible expression of the dominant negative GTPase; Sar1BH74L to block ER exit and observed that the mutated VHA-a1 proteins were not retained in the ER and importantly that the signal at the tonoplast was increased. Assuming that the a1-TD contains one or several overlapping diacidic-ER-exit motifs and no TGN-retention signal, it should re-route VHA-a3 into COPII- and Golgi-dependent trafficking to the tonoplast. However, a dual localization was observed for VHA-a3-a1-TD indicating either that the a1-TD contains a TGN-retention motif or that VHA-a3 has a yet unknown feature that actively sorts it into the provacuolar route.
Which VHA-a isoform came first, TGN/EE or tonoplast and are they redundant?
Phylogenetic analysis revealed that the VHA-a isoforms of all angiosperms as well as the gymnosperm P. taeda fall into two distinct clades. The a1-TD is conserved within the VHA-a1 clade and exchange of the N-terminal domains within this clade resulted in TGN/EE localization whereas the N-terminal domain of Selaginella that lacks the a1-TD resulted in tonoplast localization indicating that the duplication leading to differential targeting occured in a common ancestor of all seed plants. With the exception of Marchantia and Salvina culculatta, all other plants outside of the spermatophytes possess multiple VHA-a isoforms, however the a1-TD is not conserved in these plants. Evidence for differential targeting of VHA-a isoforms is missing but it seems reasonable to speculate that differential localization arose independently. Based on this, Marchantia might represent the ancestral ground state with a single VHA-isoform predominantly localized at the tonoplast. The trafficking machinery of Marchantia has distinct features and it remains to be determined if acidification of the TGN/EE plays a similar role as in Arabidopsis. If so, Marchantia would indeed represent a unique model for a single VHA-a isoform being sufficient for V-ATPase activity in two locations. A comparable situation has been proposed for the ancestral form of Stv1p and Vph1p that also served both Golgi- and vacuole acidification based on slow anterograde post-Golgi trafficking (Finnigan et al., 2011). In support of this notion, it has recently been shown that V-ATPase complexes containing Vph1p contribute more to Golgi-acidification while en route to the vacuole than the Golgi-localized complexes containing Stv1p (Deschamps et al., 2020). The identity of the VHA-a isoform does not only determine the subcellular localization, it can influence a number of biochemical properties including assembly, coupling efficiency, protein abundance and reversible dissociation (Leng et al., 1998; Kawasaki-Nishi et al., 2001b, Kawasaki-Nishi et al., 2001a).
Given the central role of the V-ATPase in many cellular functions it is important to not only understand the sorting of this molecular machine but also to address the functional diversification of plant VHA-a isoforms. The fact that VHA-a1 subunits with mutations in the a1-TD which are mislocalized to the tonoplast can complement the dwarf phenotype of the vha-a2 vha-a3 mutant to varying degrees argues strongly that they have retained their basal proton-pumping activity and do not require the specific lipid-environment of the TGN/EE for functionality. However, VHA-a1 containing V-ATPases at the vacuole cannot acidify the vacuole to wildtype levels and it remains to be determined if this is caused eg. by a difference in coupling efficiency (Kawasaki-Nishi et al., 2001b) or pH-dependent feedback regulation (Rienmüller et al., 2012). In contrast, when the Marchantia VHA-a was expressed in vha-a2 vha-a3, the resulting transgenic lines were indistinguishable from wildtype. As null alleles of single copy-encoded VHA-subunits cause male gametophyte lethality (Dettmer et al., 2005), we assumed that we would have to transform heterozygous vha-a1/+ plants to analyse in their progeny if MpVHA-a and VHA-a3-a1-TD can functionally replace VHA-a1 at the TGN/EE. As we encountered T-DNA related silencing problems in earlier complementation experiments, we used CRISPR/Cas9 under control of an egg cell-specific promoter (Wang et al., 2015) to generate new vha-a1 mutant alleles. Although vha-a1 mutants were nearly indistinguishable from wildtype during vegetative growth they were found to be completely sterile due to a block in male gametophyte development. Given that we have previously reported that inducible knock-down of VHA-a1 via either RNAi or amiRNA causes strong cell expansion defects (Brüx et al., 2008), the fact that plants homozygous for several independent null alleles are viable could be due to genetic compensation (El-Brolosy and Stainier, 2017). Indeed, when we systematically reduced the number of wildtype alleles for VHA-a2 and VHA-a3 in the vha-a1 null background, a vegetative phenotype became manifest indicating that the two tonoplast-isoforms are able to complement for the lack of VHA-a1 during vegetative growth. It remains to be determined if VHA-a2 or VHA-a3 are up-regulated in the vha-a1 null background, however we propose a model in which compensation is explained by rerouting of the tonoplast isoforms to theTGN/EE.
VHA-a1 and VHA-a3 compete for entry into COPII vesicles
We assume that Marchantia reflects the ancestral state in which ER-export is mediated via COPII and acidification of the TGN/EE is either not required or is provided by V-ATPase complexes passing through en route to the tonoplast. Similar to yeast, the duplication of the ancestral VHA-a isoform provided the basis for a dedicated TGN/EE-isoform which in turn allowed the evolution of a Golgi-independent trafficking route from the ER to the tonoplast, which we have shown previously to exist in meristematic root cells (Viotti et al., 2013). Alternatively, a single VHA-isoform could be able to enter both trafficking routes and it will thus be of importance to determine the evolutionary origin of the provacuolar route.
We propose here that in Arabidopsis, cells that require large quantities of tonoplast membranes such as the cells of the root tip, the VHA-a isoforms compete for entry into COPII vesicles. Assembly of all VO-complexes takes place at the ER with the help of dedicated assembly factors and based on RNA expression and enzymatic activity the ratio of VHA-a1 to VHA-a3 is roughly 1:10 (Hanitzsch et al., 2007; Neubert et al., 2008). We assume that due to the presence of the a1-TD, the affinity of VHA-a1 for Sec24, the cargo receptor for COPII vesicles, is much higher than the affinity of VHA-a3. Small quantities of VHA-a3 can enter COPII vesicles but the majority is predominantly sorted into the Golgi-independent route (Figure 11).The fact that VHA-a3 is strongly trafficked to the TGN/EE by the insertion of the a1-TD provides further evidence for our competition model which also predicts that blocking of COPII-mediated export by DN-Sar1-GTP would lead to re-routing of VHA-a1 to the tonoplast when its affinity for COPII is reduced by point mutations in the a1-TD. The fact that tonoplast localization of a1-TD mutants is enhanced in the vha-a2 vha-a3 background provides further support for a competition between the isoforms for the two ER-exits. Similarly, in the absence of VHA-a1, VHA-a2/VHA-a3 would be able to enter the COPII-route more efficiently and would thus be able to compensate for the lack of V-ATPases at the TGN/EE providing an explanation for the lack of a vegetative phenotype of the vha-a1 mutant. We have shown previously that despite their seemingly strict spatial separation, the activities of the V-ATPase at the TGN/EE and the tonoplast are coordinated. In the vha-a2 vha-a3 mutant in which VHA-a1 is the only remaining target, treatment with the V-ATPase inhibitor ConcA causes an additional increase in vacuolar pH (Kriegel et al., 2015) although the localization of VHA-a1 is not affected. Why then, does VHA-a3 not enter the COPII-route and rescue the reduced activity at the TGN/EE in the RNA-mediated knock-down lines of VHA-a1? Based on our model, we propose that in the RNA-mediated knock-down lines, although levels of VHA-a1 are reduced, they are still sufficient to outcompete VHA-a3 for entry into COPII vesicles.
In conclusion, our study provides the first report on a targeting motif for subunit a from a multicellular eukaryotic V-ATPase. Interaction partners that recognize and bind to the a1-TD need to be identified. ER-exit motifs as well as provacuole directing motifs in VHA-a3 should also be uncovered to better understand its trafficking via COPII vesicles and provacuoles. Our results incite further investigations to determine additional functions of the endosomal V-ATPase during vegetative growth and to learn what new functions it acquired during evolution that are essential for pollen development.
Materials and methods
Plant materials and growth conditions
Request a detailed protocolArabidopsis thaliana, Columbia 0 (Col-0) ecotype was used in all experiments in this study. The vha-a2 vha-a3 double mutant was characterized by Krebs et al., 2010. VHA-a1p:VHA-a1-GFP, VHA-a3p:VHA-a3-GFP and VHA-a1p:VHA-a1-RFP lines were established by Dettmer et al., 2006. BRI1-GFP lines were previously established by Geldner et al., 2007.
Growth of Arabidopsis seedlings for confocal microscopy was performed on plates. The standard growth medium used contained 1/2 Murashige and Skoog (MS), 0.5% sucrose, 0.5% phyto agar, 10 mM MES and the pH was set to 5.8 using KOH. Agar and MS basal salt mixture were purchased from Duchefa. Seeds were surface sterilized with ethanol and stratified for 48 hr at 4°C. Plants were grown in long day conditions (LD; 16 hr light/8 hr dark) for 5 days.
For the rosette phenotype assays, seeds were stratified for 48 hr at 4°C and then placed on soil. Seedlings were transferred to individual pots at 7 days after germination (DAG). Plants were grown either in LD or short day conditions (SD; 10 hr light/14 hr dark) for the required time.
M. polymorpha accession Takaragaike-1 (Tak-1, male; Ishizaki et al., 2008) was used in this study. The growth condition and the transformation method were described previously (Kubota et al., 2013; Kanazawa et al., 2016). M. polymorpha expressing mRFP-MpSYP6A was generated in the previous work (Kanazawa et al., 2016).
Construct design and plant transformation
VHA-a1/VHA-a3 Chimeras
Request a detailed protocolFive chimeric proteins were made which consisted of increasing lengths of the VHA-a1 N-terminus (37aa, 85aa, 131aa, 179 aa and 228 aa) fused to decreasing lengths of the C-terminal domain of VHA-a3. The GreenGate cloning system was used (Lampropoulos et al., 2013). Unique overhangs for each chimera were designed to allow seamless fusion of the VHA-a1 and VHA-a3 cDNA sequences (Supplementary file 1). The primers used are listed in Supplementary file 1. All PCR products were blunt end cloned into the pJET1.2 vector (ThermoFisher Scientific) and verified by sequencing (Eurofins). Verified clones were digested using either Bglll or NotI and ClaI to release the fragments. To combine the VHA-a1 and VHA-a3 fragments, the GreenGate cloning system was applied using modules described in Supplementary file 1.
UBQ:VHA-a3-a1-TD-GFP
Request a detailed protocolThe targeting domain of VHA-a1 (a1-TD) was introduced into VHA-a3 by PCR techniques. The N and C-termini of VHA-a3 and the a1-TD were amplified using primers in Supplementary file 1. All PCR products were blunt end cloned into the pJET1.2 vector and verified by sequencing (Eurofins). Verified clones were digested using either Bglll or NotI and ClaI to release the fragments. The destination vector was made by combining the fragments with GreenGate modules described in Supplementary file 1.
Site-directed mutagenesis of VHA-a1
Request a detailed protocolThe pJET1.2 clone carrying VHA-a1 NT with 179 aa from the chimeras was used as a template for site-directed mutagenesis of VHA-a1. PCR mediated site-directed mutagenesis was performed using primers indicated in Supplementary file 1 and according to Wang and Malcolm, 1999. Five cycles of single-primer extension reactions were carried out before the standard QuikChange Site-Directed Mutagenesis Protocol (Stratagene, La Jolla, CA, USA). Arabidopsis VHA-a1 intron 10 was amplified from genomic DNA, a connecting VHA-a1 fragment (Y2) and a C-terminal fragment (VHA-a1-CT) were amplified from VHA-a1 cDNA using primers indicated in Supplementary file 1. All PCR fragments were subcloned into the pJet1.2 vector and were verified by sequencing (Eurofins). Verified clones were digested using Bglll to release the fragments. The final destination vectors were made with modules described in Supplementary file 1.
UBQ:VHA-a3R729N-GFP
Request a detailed protocolTo generate UBQ:VHA-a3R729N-GFP, the full length coding sequence of VHA-a3 without Stop codon, was amplified from Arabidopsis Col-0 wild type cDNA using primer pairs a3SacII-Fw and a3SalI-Rv (Supplementary file 1). The resulting 2482 bp fragment was subcloned into pJET2.1blunt to receive pJET1.2blunt-cVHA-a3-STOP. Next, pJET1.2blunt-cVHA-a3-STOP was used as a template to perform site-directed mutagenesis PCR using primer combination cVHA-a3-R729N-Fw and cVHA-a3-R729N-Rv (Supplementary file 1) to receive pJET1.2blunt-cVHA-a3R729N-STOP. The coding sequence of GFP5(S65T) and an interconnecting linker was amplified from VHA-a3-GFP (Dettmer et al., 2006) with primer combination GFP5-SalI-Fw2 and GFP5-PstI-Rv2 (Supplementary file 1). The resulting 764 bp fragment was subcloned into pJET1.2blunt to receive pJET1.2blunt-GFP(SalI/PstI). The coding sequences of VHA-a3, GFP5 and the mutation R729N were verified by sequencing (Eurofins). To obtain the binary vector UBQ10-VHA-a3-R729N-GFP, pUTKan (Krebs et al., 2012) and pJET1.2blunt-cVHA-a3R729N-STOP were digested with SacII and SalI. The resulting 10847 and 2473 bp fragments were ligated to receive UBQ10-VHA-a3-R729N. Finally, UBQ10-VHA-a3-R729N and pJET1.2blunt-GFP(SalI/PstI) were digested with SalI and PstI. The 13291 and 756 bp fragments were ligated to obtain UBQ10-VHA-a3-R729N-GFP.
Dex:Sar1BH74L-CFP
Request a detailed protocolThe Arabidopsis Sar1B cDNA sequence with the H74L mutation was synthesized by Eurofins. The synthesized fragment contained Eco31I at its 5' and 3' ends for subcloning. PCR was performed on the synthesized fragment using primers indicated in Supplementary file 1 and the PCR product was blunt end cloned into the pJET1.2 vector. The SarB1H74L sequence was verified by sequencing (Eurofins). Verified clones were digested using NotI and ClaI to release the SarB1H74L fragments. The dexamethasone-inducible construct was made by using the LhG4/pOp system combined with the ligand-binding domain of the rat glucocorticoid receptor (GR; Moore et al., 1998; Craft et al., 2005; Samalova et al., 2005). Two GreenGate reactions were performed to create two intermediate vectors that were later combined on one T-DNA. The first intermediate vector (pKSM002) contained the GR-LhG4 transcription factor expressed under the UBQ10 promoter and the second intermediate vector contained Sar1BH74L under the pOP6 promoter (pKSN009). The two intermediate vectors were combined on one final destination vector pGGZ003 (Supplementary file 1).
UBQ10:MpVHA-a-mVenus
Request a detailed protocolThe cDNA sequence of MpVHA-a was amplified by PCR using female Marchantia polymorpha (ecotype BoGa) thallus cDNA as a template. The sequence was amplified in two parts (MpVHA-aNT and MpVHA-aCT) separated at the exon 10-exon 11 junction using primers indicated in Supplementary file 1. All PCR products were blunt end cloned into the pJET1.2 vector and verified by sequencing (Eurofins). Verified clones were digested using Bglll to release the fragments. The final destination vector was made with the MpVHA-aNT, VHA-a1-intron10 and MpVHA-aCT fragments and modules described in Supplementary file 1.
CaMV35S:MpVHA-a-mVenus and MpEF1α:MpVHA-a-mVenus
Request a detailed protocolThe genome sequence of MpVHA-a containing all exons and introns was amplified by PCR using the Tak-1 genome as a template with primers (Supplementary file 1). The amplified products were subcloned into pENTR D-TOPO (Invitrogen), and cDNA for mVenus was inserted into the AscI site of the pENTR vectors using In-Fusion (Clontech) according to the manufacturer’s instructions. The resultant entry sequences were transferred to pMpGWB302 or pMpGWB303 (Ishizaki et al., 2015) using LR Clonase II (Invitrogen) according to the manufacturer’s instructions.
UBQ10:A.trichopoda-VHA-aNT-VHA-a1-mCherry,UBQ10:S.moellendorffii-VHA-aNT-VHA-a1-mCherry and UBQ10:P.taeda-VHA-aNT-VHA-a1-mVenus
Request a detailed protocolThe first 682 bp of the sequence identified as evm_27.model.AmTr_v1.0_scaffold00080.37 (A. trichopoda VHA-a), 649 bp of the sequence identified as 182335 (S.moellendorffii-VHA-a) and a 703 bp fragment starting from the 163rd bp of the sequence identified as 5A_I15_VO_L_1_T_29156/41278 (P.taeda-VHA-a) were synthesized (Eurofins). The synthesized fragments contained Eco31I at their 5' and 3' ends for subcloning. PCR was performed on the synthesized fragments using primers indicated in Supplementary file 1 and the PCR products were blunt end cloned into the pJET1.2 vector. A connecting VHA-a1 fragment (dodo) was also amplified from VHA-a1 cDNA using primers in Supplementary file 1. All PCR fragments were verified by sequencing (Eurofins). Verified clones were digested using either Bglll or NotI and ClaI to release the fragments. The destination vectors were made by combining the following components: A.trichopoda-VHA-aNT or P.taeda-VHA-aNT or S.moellendorffii-VHA-aNT with VHA-a1 dodo, VHA-a1-intron10, VHA-a1-CT and modules described in Supplementary file 1.
UBQ10::ST-GFP
Request a detailed protocolRat Sialyltransferase (ST) sequence was amplified from pre-existing p16:ST-pHusion plasmid (Luo et al., 2015) using GG-ST-C-fwd and GG-ST-C-rev primers (Supplementary file 1), sub-cloned into pGGC000 after Eco31I digest and sequenced to verify correct insert sequence. pGGC-ST was used in a GreenGate reaction to generate UBQ10:ST-GFP using modules described in Supplementary file 1.
UBQ10:VHA-a3-pmScarlet-I
Request a detailed protocolThe VHA-a3 cDNA sequence was amplified from cDNA with primers listed in Supplementary file 1. After digest with Eco31I, VHA-a3 was subcloned into pGGC000. pGGC-VHA-a3 was used in a GreenGate reaction to assemble the final destination vector as described in Supplementary file 1.
UBQ10:VHA-a1-GFP
Request a detailed protocolThe VHA-a1 cDNA sequence was split into two parts by PCR to generate the fragments VHA-a1 NT (Y5) and VHA-a1 CT separated at the exon 10-exon 11 junction using primers indicated in Supplementary file 1. These fragments were digested with Eco31I and then subcloned into pGGC000 with VHA-a1 intron 10. pGGC-VHA-a1-intron 10 was sequenced and used in a GreenGate reaction with modules described in Supplementary file 1 to generate the final destination vector.
Destination vector pGGZ004
Request a detailed protocolThe vector backbone pGGZ004 is a GreenGate compatible plant binary vector which, unlike previous GreenGate vector backbones (Lampropoulos et al., 2013), contains the replicase required for plasmid replication in A. tumefaciens and therefore does not require the pSOUP helper plasmid (Hellens et al., 2000). pGGZ004 is based on the plant binary vector pTKan (Krebs et al., 2012), which originated from pPZP212 (Hajdukiewicz et al., 1994). To generate pGGZ004, pTKan was cut at its multiple cloning site with KpnI followed by an incomplete digest with Eco31I, which opened the plasmid next to the T-DNA left border. The resulting 6971 bp fragment comprised the T-DNA border regions, the pBR322 bom site and the ColE1 and pVS1 plasmid origins for replication in E. coli and in Agrobacterium, respectively (Hajdukiewicz et al., 1994). GreenGate overhangs A and G were added to the opened vector backbone by sticky end ligation of annealed oligonucleotides Eco31I-w-AG-Fw and Eco31I-w-AG-Rv which contained KpnI and Eco31I compatible overhangs. Lastly, a site-directed mutagenesis using primer combination pGGZ004SDM1-Fw and pGGZ004SDM1-Rv was performed to remove an Eco31I site in the vector backbone (Supplementary file 1). The final pGGZ004 vector was fully sequenced before use for further cloning.
CRISPR/Cas9 constructs
View detailed protocolCRISPR target sites in VHA-a1 were selected using CHOPCHOP (https://chopchop.cbu.uib.no; Labun et al., 2016) and CCtop; (https://crispr.cos.uni-heidelberg.de, Stemmer et al., 2015). gRNAs with at least 4 bp difference to every other region in the Arabidopsis thaliana genome were selected as precaution against off–target mutations. gRNA sequences (Supplementary file 1) were inserted into the plasmid pHEE401E as described by the authors (Wang et al., 2015). Three CRISPR plasmids were cloned, one containing two gRNAs (gRNA1 and gRNA2), aiming at deleting the region between the two CRISPR sites, and two plasmids containing one gRNA each, gRNA3 and gRNA4.
Depending on the final destination vector used, all constructs were transformed into either of two Agrobacterium tumefaciens strains. The strain GV3101:pMP90 was used if the final destination vector was pGGZ004. Selection was done on 5 mg/ml rifampicin, 10 mg/ml gentamycin, and 100 mg/ml spectinomycin. The strain ASE1(pSOUP+) was used if pGGZ001/3 were used. Selection was done on 100 μg/ml spectinomycin, 5 μg/ml tetracycline (for pSOUP), 25 μg/ml chloramphenicol, and 50 μg/ml kanamycin. Arabidopsis plants were transformed using standard procedures. Transgenic plants were selected on MS medium containing appropriate antibiotics.
Analysis of mutations at CRISPR sites and genotyping of mutants
Request a detailed protocolGenomic DNA was extracted from rosette leaves and amplified by PCR using primers flanking the CRISPR sites (Supplementary file 1). PCR products were sequenced by Eurofins or analyzed by agarose gel electrophoresis. Presence/absence of CRISPR T-DNA was monitored by PCR with Cas9 specific primers (Supplementary file 1) followed by agarose gel electrophoresis. VHA-a2 and VHA–a3 wildtype and mutant alleles were identified by PCR using specific primers (Supplementary file 1) and agarose gel electrophoresis. Screening for vha-a1-1 was done by detecting the PCR product spanning the CRISPR sites (260 bp shorter for vha-a1-1 compared to wildtype) by agarose gel electrophoresis (Figure 9—figure supplement 1). Screening for vha-a1 mutants with other vha-a1 alleles was performed by observation of the pollen phenotype after having established that the defect in pollen development is caused by knockout of VHA-a1 (Supplementary file 1).
Pharmacological treatments and stains
Request a detailed protocolArabidopsis seedlings were incubated in liquid 1/2 MS medium with 0.5% sucrose, pH 5.8 (KOH), containing 50 µM BFA, 1 µM FM4-64, 60 µM DEX and 125 nM ConcA or the equivalent amount of DMSO in control samples for the required time at room temperature. Stock solutions were prepared in DMSO.
Hypocotyl length measurement
Request a detailed protocolFor the hypocotyl length measurements, seeds were sown on agar plates made with distilled water without nutrients containing 2.5 mM MES, pH 5.8 adjusted with KOH, 0.8% phytoagar and 125 nM ConcA (Santa Cruz Biotechnology) dissolved in DMSO. After stratification for 48 hr at 4°C, a 4 hr light stimulus was given at 22°C. With the plates wrapped in aluminium foil the seedlings were grown in the dark for 4 days, then lain on overhead transparency sheets in drops of water and scanned. Quantification of hypocotyl length from the scans was done using Fiji (Fiji is just image J; Schindelin et al., 2012).
Alexander staining for pollen viability
Request a detailed protocolAnthers from flowers at stage 13 (Smyth et al., 1990) were incubated in Alexander’s stain (Alexander, 1969) on objective slides covered with coverslips at room temperature for 20 hr. Bright field images were taken using a Zeiss Axio Imager M1 microscope. Adjustment of brightness and contrast of whole images was done with Fiji.
Confocal microscopy image acquisition
Request a detailed protocolArabidopsis root cells of 6 day-old seedlings were analyzed by confocal laser scanning microscopy (CLSM) using a Leica TCS SP5II microscope equipped with a Leica HCX PL APO lambda blue 63.0 × 1.20 UV water immersion objective. Pictures were taken above the division zone. CFP, GFP and mVenus were excited at 458 nm, 488 nm and 514 nm with a VIS-argon laser respectively. mRFP, mCherry, pmScarlet-I and FM4-64 were excited at 561 nm with a VIS-DPSS 561 laser diode. For image acquisition, the Leica Application Suite Advanced Fluorescence software was used. Image post processing operations (Gaussian blur between 0.6 and 0.7; adjustment of brightness and contrast of the whole image) were done with Fiji.
CLSM on Marchantia thalli was performed according to Kanazawa et al., 2016 with a minor modification. Briefly, the dorsal cells of 5-day-old thalli were observed using LSM780 confocal microscope (Carl Zeiss) equipped with an oil immersion lens (63×, N.A. = 1.4). For the spectral imaging, The samples were excited at 488 nm (Argon 488) and 561 nm (DPSS 561–10), and the emission was collected by 20 GaAsP detectors between 482 and 659 nm. Spectral unmixing and processing of the obtained images were conducted using ZEN2012 software (Carl Zeiss). For semi-superresolutional imaging, the samples were observed by LM780 confocal microscope equipped with Airyscan (Carl Zeiss). The samples were excited at 488 nm and 561 nm, and the emission was separated by BP420-480+BP495-550 (for mVenus), BP420-480+BP495-620 (for mRFP), and BP570-620+LP645 (for chlorophyll autofluorescence) filters, and detected by Airyscan. The semi-superresolutional processing was conducted using ZEN 2.3 SP1 software (Carl Zeiss).
Image analysis and quantifications
Colocalization analyses
Request a detailed protocolColocalization analysis was performed with the JACoP plugin on Image J (Bolte, 2006). Nine to ten images per condition were subjected to a mean blur of 2 pixels and a background subtraction with a rolling ball radius of 20 pixels. The Pearson correlation coefficient and the Mander’s correlation coefficients M1 (VHA-a1-mRFP) and M2 (mutated VHA-a1-GFP) were determined (Dunn et al., 2011). For Mander’s analysis a uniform threshold was set to all pictures to avoid background contribution. OriginPro was used to determine if the data were normally distributed (Shapiro-Wilk test) and to perform statistical tests. The data was found to be normally distributed and the student T-test was used to test for significance between means.
Tonoplast signal intensity
Request a detailed protocolAverage tonoplast intensities were measured in roots of Arabidopsis seedlings expressing DEX:Sar1BH74L-CFP and UBQ10:VHA-a1-GFP with mutations or VHA-a3:VHA-a3-GFP after 6 hr induction with 60 µM DEX (Sigma-Aldrich). Images were acquired sequentially. In the first sequential scan, GFP was excited at 488 nm and emission detected at 500–557 nm. In the second sequential scan CFP was exited at 458 nm and emission detected between 468–481 nm. Intensity measurements were done using plot profiles in ImageJ. The maximum intensity values along line profiles across the TGN/EE and tonoplast were recorded. OriginPro was used to perform statistical tests. The data was found not to be normally distributed and the Mann-Whitney test was used to test for significance between means.
Ratio of TGN/EE-to-tonoplast signal intensity
Request a detailed protocolTo determine the ratio of TGN/EE-to-tonoplast fluorescence intensity, images of root cells expressing UBQ10:VHA-a1-GFP with mutations in the wildtype and vha-a2 vha-a3 background were acquired.The maximum intensity values along line profiles across the TGN/EE and tonoplast were measured using ImageJ. The average TGN/EE and tonoplast intensities were calculated for each image and the TGN/EE-to-tonoplast fluorescence intensity ratio was calculated. The TGN/EE-to-tonoplast fluorescence intensity ratios for n ≥ 10 images for each mutation were averaged and statistical tests were performed with OriginPro. The student T-test was used to test for significance between means.
High-pressure freezing, freeze substitution, and electronmicroscopy
Request a detailed protocolSeven-day-old Arabidopsis wildtype seedlings expressing DEX:AtSar1b-GTP were induced for 6 hr with 60 μM DEX. Seedlings were then processed as previously described (Scheuring et al., 2011). Freeze substitution was performed in a Leica EM AFS2 freeze substitution unit in dry acetone supplemented with 0.3% uranyl acetate as previously described (Hillmer et al., 2012). Root tips were cut axially with a Leica Ultracut S microtome to obtain ultrathin sections. Sections were examined in a JEM1400 transmission electron microscope (JEOL) operating at 80 kV. Micrographs were recorded with a TemCam F416 digital camera (TVIPS, Gauting, Germany).
pH measurements
Request a detailed protocolCell sap pH measurements were performed as previously described (Krebs et al., 2010).
Tonoplast membrane preparation
View detailed protocolTonoplast membranes were prepared from rosette leaves of 6-week-old plants grown under short day conditions as previously described by Barkla et al., 1999; Leidi et al., 2010.
SDS-PAGE and immunoblotting
Request a detailed protocolMicrosomal membrane and tonoplast membrane proteins were analyzed by SDS-PAGE and subsequent immunoblotting. Upon gel electrophoresis, the proteins were transferred to a PVDF membrane (Bio-Rad). The primary antibodies against VHA-a1 (AS142822) and VHA-a3 (AS204369) were purchased from Agrisera and were used in a dilution of 1:1000 in 2% BSA-TBS-T. The primary antibody against VHA-B (Ward et al., 1992) and anti-GFP (Roth et al., 2018) were previously described. Antigen on the membrane was visualized with horseradish peroxidase-coupled anti-rabbit IgG (Promega) for VHA-a1 and VHA-a3, anti-mouse IgG (Sigma) for VHA-B and chemiluminescent substrate (Peqlab). Immunostained bands were analyzed using a cooled CCD camera system (Intas).
Phylogenetic analysis
Request a detailed protocolFor phylogenetic reconstruction in a first step the best molecular evolutionary model was determined by running the program PartitionFinder 2 (Lanfear et al., 2017). Phylogenetic reconstruction was then performed by running raxml-ng (Kozlov et al., 2019) setting the model to JTT+I+G and starting the analysis from 10 most parsimonious and 10 random trees. In order to estimate the reliability of the phylogenetic reconstruction 500 bootstrap replicates were run.
Multiple sequence alignments
Request a detailed protocolMultiple sequence alignments were performed using Clustal omega (Madeira et al., 2019). Aligned sequences were analyzed in Geneious 10.1.3.
Homology modeling
Request a detailed protocol3D models of the VHA-a1 and VHA-a3 N-termini were obtained through homology modeling with cryo-EM derived models of Stv1p-VO subcomplex (PDB6O7U) and Vph1p-VO subcomplex (PDB6O7T) as templates (Vasanthakumar et al., 2019). Homology modeling was performed according to Roy et al., 2010.
Accession numbers
Request a detailed protocolSequence data from this article can be found in the Arabidopsis Genome Initiative or GenBank/EMBL databases under the following accession numbers: VHA-a1, At2g28520; VHA-a3, At4g39080; Sar1B, AT1G56330.1. Data for Marchantia polymorpha VHA-a can be found on the Marchantia genome database with the following ID: Mp3g15140.1. Amborella trichopoda and Selaginella moellendorffii VHA-a sequences can be found on the Phytozome platform (Goodstein et al., 2012) with the following identifiers: A.trichopoda VHA-a; evm_27.TU.AmTr_v1.0_scaffold00080.37 and S.moellendorffii VHA-a; 182335. Pinus taeda sequence data can be found on the PineRefSeq project on the TreeGenes platform (Falk et al., 2019; Wegrzyn et al., 2008) with the following identifier: 5A_I15_VO_L_1_T_29156/41278.
Appendix 1
Appendix 2
Legend to Supplementary file 1
Table 1 | Statistical support for the predicted structural homology models obtained from I-TASSER |
Table 2 | Overview of the number of VHA-a isoforms in plants with VHA-a isoforms that do not contain the a1-TD |
Table 3 | Overview T1 analysis CRISPR/Cas9 VHA-a1 lines |
Table 4 | Segregation of vha-a1-1 |
Table 5 | Chimeric proteins overhangs and primers |
Table 6 | GreenGate modules and constructs |
Table 7 | CRISPR/Cas9 gRNA sequences and genotyping primers |
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data is provided for Figures 2,4B, 5B, 7B and C, 8C and D, 10A, Supplemental Figure 8B, Supplemental Figure 9 and Supplemental Figure 15A.
-
Dryad Digital RepositoryThe hornwort genome and early land plant evolution.https://doi.org/10.5061/dryad.msbcc2ftv
References
-
Differential staining of aborted and nonaborted pollenStain Technology 44:117–122.https://doi.org/10.3109/10520296909063335
-
Direct interaction of the golgi V-ATPase a-subunit isoform with PI(4)P drives localization of golgi V-ATPases in yeastMolecular Biology of the Cell 28:2518–2530.https://doi.org/10.1091/mbc.e17-05-0316
-
Stacks on tracks: the plant golgi apparatus traffics on an actin/ER networkThe Plant Journal : For Cell and Molecular Biology 15:441–447.https://doi.org/10.1046/j.1365-313x.1998.00208.x
-
JACoP v2.0: improving the user experience with co-localization studiesImageJ User Dev 224:213–232.
-
Recent insights into the structure, regulation, and function of the V-ATPasesTrends in Biochemical Sciences 40:611–622.https://doi.org/10.1016/j.tibs.2015.08.005
-
Essential role of the V-ATPase in male gametophyte developmentThe Plant Journal 41:117–124.https://doi.org/10.1111/j.1365-313X.2004.02282.x
-
A practical guide to evaluating colocalization in biological microscopyAmerican Journal of Physiology-Cell Physiology 300:C723–C742.https://doi.org/10.1152/ajpcell.00462.2010
-
The reconstructed ancestral subunit a functions as both V-ATPase isoforms Vph1p and Stv1p in Saccharomyces cerevisiaeMolecular Biology of the Cell 22:3176–3191.https://doi.org/10.1091/mbc.e11-03-0244
-
Vacuolar ATPases: rotary proton pumps in physiology and pathophysiologyNature Reviews Molecular Cell Biology 8:917–929.https://doi.org/10.1038/nrm2272
-
Vacuolar-type ATPase: a proton pump to lysosomal traffickingProceedings of the Japan Academy, Series B 95:261–277.https://doi.org/10.2183/pjab.95.018
-
Endosomal signaling of plant steroid receptor kinase BRI1Genes & Development 21:1598–1602.https://doi.org/10.1101/gad.1561307
-
Phytozome: a comparative platform for green plant genomicsNucleic Acids Research 40:D1178–D1186.https://doi.org/10.1093/nar/gkr944
-
The small, versatile pPZP family of Agrobacterium binary vectors for plant transformationPlant Molecular Biology 25:989–994.https://doi.org/10.1007/BF00014672
-
pGreen: a versatile and flexible binary ti vector for Agrobacterium-mediated plant transformationPlant Molecular Biology 42:819–832.https://doi.org/10.1023/a:1006496308160
-
SNARE molecules in Marchantia polymorpha: unique and conserved features of the membrane fusion machineryPlant and Cell Physiology 57:307–324.https://doi.org/10.1093/pcp/pcv076
-
The Amino-terminal Domain of the Vacuolar Proton-translocating ATPase a Subunit Controls Targeting and in Vivo Dissociation, and the Carboxyl-terminal Domain Affects Coupling of Proton Transport and ATP HydrolysisJournal of Biological Chemistry 276:47411–47420.https://doi.org/10.1074/jbc.M108310200
-
Yeast V-ATPase Complexes Containing Different Isoforms of the 100-kDa a-subunit Differ in Coupling Efficiency and in Vivo DissociationJournal of Biological Chemistry 276:17941–17948.https://doi.org/10.1074/jbc.M010790200
-
Efficient Agrobacterium-mediated transformation of the liverwort Marchantia polymorpha using regenerating thalliBioscience, Biotechnology, and Biochemistry 77:167–172.https://doi.org/10.1271/bbb.120700
-
CHOPCHOP v2: a web tool for the next generation of CRISPR genome engineeringNucleic Acids Research 44:W272–W276.https://doi.org/10.1093/nar/gkw398
-
PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analysesMolecular Biology and Evolution 34:772–773.https://doi.org/10.1093/molbev/msw260
-
Function of the COOH-terminal domain of Vph1p in activity and assembly of the yeast V-ATPaseJournal of Biological Chemistry 273:6717–6723.https://doi.org/10.1074/jbc.273.12.6717
-
The EMBL-EBI search and sequence analysis tools APIs in 2019Nucleic Acids Research 47:W636–W641.https://doi.org/10.1093/nar/gkz268
-
Dynamics of three-dimensional replication patterns during the S-phase, analysed by double labelling of DNA and confocal microscopyJournal of Cell Science 103:857–862.
-
Measurement of co-localization of objects in dual-colour confocal imagesJournal of Microscopy 169:375–382.https://doi.org/10.1111/j.1365-2818.1993.tb03313.x
-
Eukaryotic V-ATPase: novel structural findings and functional insightsBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1837:857–879.https://doi.org/10.1016/j.bbabio.2014.01.018
-
Brefeldin A: deciphering an enigmatic inhibitor of secretionPlant Physiology 130:1102–1108.https://doi.org/10.1104/pp.011569
-
Luminal and Cytosolic pH Feedback on Proton Pump Activity and ATP Affinity of V-type ATPase from ArabidopsisJournal of Biological Chemistry 287:8986–8993.https://doi.org/10.1074/jbc.M111.310367
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Early flower development in ArabidopsisThe Plant Cell 2:755–767.https://doi.org/10.1105/tpc.2.8.755
-
A simple nomenclature for a complex proton pump: vha genes encode the vacuolar H+-ATPaseTrends in Plant Science 7:157–161.https://doi.org/10.1016/S1360-1385(02)02240-9
-
Function of the anion transporter AtCLC-d in the trans-Golgi networkThe Plant Journal : For Cell and Molecular Biology 50:466–474.https://doi.org/10.1111/j.1365-313X.2007.03061.x
-
Seventeen a-subunit isoforms of Paramecium V-ATPase provide high specialization in localization and functionMolecular Biology of the Cell 17:917–930.https://doi.org/10.1091/mbc.e05-06-0511
-
TreeGenes: a forest tree genome databaseInternational Journal of Plant Genomics 2008:1–7.https://doi.org/10.1155/2008/412875
-
The PC6B cytoplasmic domain contains two acidic clusters that direct sorting to distinct trans-Golgi network/endosomal compartmentsMolecular Biology of the Cell 11:1257–1273.https://doi.org/10.1091/mbc.11.4.1257
Article and author information
Author details
Funding
DFG (SFB1101 TPA02)
- Karin Schumacher
Ministry of Education, Culture, Sports, Science and Technology (19H05675)
- Takashi Ueda
Ministry of Education, Culture, Sports, Science and Technology (18H02470)
- Takashi Ueda
Ministry of Education, Culture, Sports, Science and Technology (18K14738)
- Takehiko Kanazawa
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Prof. Dr. Sabine Zachgo (University of Osnabrück) for kindly providing us with Marchantia polymorpha (ecotype BoGa) cDNA. We also thank Prof. Dr. Stephan Rensing (University of Marburg) for providing us with Selaginella moellendorffii leaf tissue. We thank Dr. Rainer Waadt (COS, University of Heidelberg) for kindly providing GreenGate entry vectors. We thank Dr. Alyona Minina (COS, University of Heidelberg) for assisting with Alexander staining of pollen. We thank Dr. Takayuki Kohchi, Dr. Ryuichi Nishihama (Kyoto University), and Dr. Kimitsune Ishizaki (Kobe University) for vectors. We also thank Dr. Kazuo Ebine, Mayuko Yamamoto, Koji Hayashi (NIBB) and Fabian Fink (COS, University of Heidelberg) for supporting experiments. The support in M. polymorpha cultivation was provided by the Model Plant Research Facility of National Institute for Basic Biology. This work was supported by the Deutsche Forschungsgemeinschaft within SFB1101 to KS and by Grants-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (to TU, 19H05675, and 18H02470, and TK, 18K14738). Electron microscopy was performed by S Hillmer at the Electron Microscopy Core Facility of Heidelberg University with technical assistance of S Gold.
Copyright
© 2020, Lupanga et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,064
- views
-
- 493
- downloads
-
- 27
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Plant Biology
Plants distribute many nutrients to chloroplasts during leaf development and maturation. When leaves senesce or experience sugar starvation, the autophagy machinery degrades chloroplast proteins to facilitate efficient nutrient reuse. Here, we report on the intracellular dynamics of an autophagy pathway responsible for piecemeal degradation of chloroplast components. Through live-cell monitoring of chloroplast morphology, we observed the formation of chloroplast budding structures in sugar-starved leaves. These buds were then released and incorporated into the vacuolar lumen as an autophagic cargo termed a Rubisco-containing body. The budding structures did not accumulate in mutants of core autophagy machinery, suggesting that autophagosome creation is required for forming chloroplast buds. Simultaneous tracking of chloroplast morphology and autophagosome development revealed that the isolation membranes of autophagosomes interact closely with part of the chloroplast surface before forming chloroplast buds. Chloroplasts then protrude at the site associated with the isolation membranes, which divide synchronously with autophagosome maturation. This autophagy-related division does not require DYNAMIN-RELATED PROTEIN 5B, which constitutes the division ring for chloroplast proliferation in growing leaves. An unidentified division machinery may thus fragment chloroplasts for degradation in coordination with the development of the chloroplast-associated isolation membrane.
-
- Plant Biology
Photosynthetic organisms exhibit remarkable diversity in their light-harvesting complexes (LHCs). LHCs are associated with photosystem I (PSI), forming a PSI-LHCI supercomplex. The number of LHCI subunits, along with their protein sequences and pigment compositions, has been found to differ greatly among the PSI-LHCI structures. However, the mechanisms by which LHCIs recognize their specific binding sites within the PSI core remain unclear. In this study, we determined the cryo-electron microscopy structure of a PSI supercomplex incorporating fucoxanthin chlorophyll a/c-binding proteins (FCPs), designated as PSI-FCPI, isolated from the diatom Thalassiosira pseudonana CCMP1335. Structural analysis of PSI-FCPI revealed five FCPI subunits associated with a PSI monomer; these subunits were identified as RedCAP, Lhcr3, Lhcq10, Lhcf10, and Lhcq8. Through structural and sequence analyses, we identified specific protein–protein interactions at the interfaces between FCPI and PSI subunits, as well as among FCPI subunits themselves. Comparative structural analyses of PSI-FCPI supercomplexes, combined with phylogenetic analysis of FCPs from T. pseudonana and the diatom Chaetoceros gracilis, underscore the evolutionary conservation of protein motifs crucial for the selective binding of individual FCPI subunits. These findings provide significant insights into the molecular mechanisms underlying the assembly and selective binding of FCPIs in diatoms.