Next Article in Journal
Vernonia britteniana Root Phytochemical Studies, In Vitro Cercaricidal Activity on the Larval Stage of Schistosoma mansoni and Antioxidant Activities
Previous Article in Journal
Cnidium monnieri (L.) Cusson Flower as a Supplementary Food Promoting the Development and Reproduction of Ladybeetles Harmonia axyridis (Pallas) (Coleoptera: Coccinellidae)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Depth Topological Analysis of Arabidopsis Mid-SUN Proteins and Their Interaction with the Membrane-Bound Transcription Factor MaMYB

1
Department of Biological and Molecular Sciences, Faculty of Health and Life Sciences, Oxford Brookes University, Oxford OX3 0BP, UK
2
CNRS, Inserm, GReD Clermont-Ferrand, Université Clermont Auvergne, 63001 Clermont-Ferrand, France
*
Author to whom correspondence should be addressed.
Plants 2023, 12(9), 1787; https://doi.org/10.3390/plants12091787
Submission received: 16 March 2023 / Revised: 20 April 2023 / Accepted: 24 April 2023 / Published: 27 April 2023
(This article belongs to the Section Plant Cell Biology)

Abstract

:
Mid-SUN proteins are a neglected family of conserved type III membrane proteins of ancient origin with representatives in plants, animals, and fungi. Previous higher plant studies have associated them with functions at the nuclear envelope and the endoplasmic reticulum (ER). In this study, high-resolution confocal light microscopy is used to explore the localisation of SUN3 and SUN4 in the perinuclear region, to explore topology, and to study the role of mid-SUNs on endoplasmic reticulum morphology. The role of SUN3 in the ER is reinforced by the identification of a protein interaction between SUN3 and the ER membrane-bound transcription factor maMYB. The results highlight the importance of mid-SUNs as functional components of the ER and outer nuclear membrane.

1. Introduction

Sad1/Unc84 (SUN) domain proteins are conserved across all eukaryotes and can be classified into two subgroups, Cter-SUN proteins, with a C-terminal SUN domain, and mid-SUNs, in which the SUN domain occupies a central location in the protein [1,2]. Cter-SUNs are inner nuclear membrane (INM)-localised and interact with outer nuclear membrane (ONM)-localised plant Klarsicht/Anc1/Syne homology (KASH) domain proteins to form Linker of Nucleoskeleton and Cytokseleton (LINC) complexes with multiple functions [3,4,5,6,7,8,9]. On the other hand, mid-SUN proteins are type III transmembrane proteins and information on their functions is limited to a small number of systems [2,10,11,12,13,14,15,16,17,18,19]. In yeast and metazoans, they are predominantly ER-located with potential functions in integral membrane protein folding in the ER [12]. In maize, three mid-SUN proteins, ZmSUN3-5, were identified, two (ZmSUN3 & 4) being expressed ubiquitously while ZmSUN5 is most abundant in pollen [2]. This expression pattern was also found to be true for the three Arabidopsis mid-SUNs (SUN3, SUN4, and SUN5) [13] and for GbSUN5 in cotton [20]. The Arabidopsis mid-SUNs have been shown to interact with LINC components–Cter-SUNs and KASH proteins–and to play a role in nuclear morphology and shape [13]. Arabidopsis SUN3 and SUN4 were found to localise in both the NE and ER, which was also shown for chickpea mid-SUN protein CaSUN1 [14]. CaSUN1 was identified in a screen investigating potential components of the dehydration response pathway, which involves the ER-localised unfolded protein response pathway [14]. Recently, a triple sun3/sun4/sun5 mutant was shown to have impaired pollen tubes and plant growth with leucine-rich receptor kinases normally destined for the plasma membrane missorted to acidic compartments [18]. Mid-SUNs were identified as part of a POD1–SUN–CRT3 chaperone complex guarding the ER sorting of the LRR receptor. Xue et al. [17] found that loss of SUN3/4/5 caused tubule-to-sheet transition and overexpression that led to ER deformation, suggesting that the coiled-coil domain contributes to the maintenance of ER morphology.
Taken together, the previous work described above leaves a number of areas of uncertainty about the location, topology, and function of mid-SUN proteins in plants. A NE localisation has been suggested, with evidence of interaction with NE components, while there is mounting evidence for roles in the endoplasmic reticulum, including as part of a chaperone complex involved in the ER sorting of a plasma membrane receptor and in the determination of the morphology of the ER. Here, we explore in detail the location of mid-SUN proteins at the outer nuclear membrane and ER. We explore the role of mid-SUNs in determining the structure and properties of the ER and provide novel insights into the topology of mid-SUN protein domains. We finally provide evidence demonstrating that mid-SUN proteins interact with an ER-bound transcription factor, maMYB.

2. Results

2.1. Mid-SUN Proteins Are Enriched in the ONM and ER

With previous reports of ER and NE localisation of plant, yeast, and metazoan mid-SUN proteins, we firstly aimed to better understand the subcellular localisation of Arabidopsis SUN3 and SUN4. We therefore compared their ER/NE ratio using fluorescence intensity (FI) of ER and NE markers in transiently expressing Nicotiana benthamiana (N. benthamiana) leaf epidermal cells by establishing the intensity of fluorescence in regions of interest (ROI) in the NE and cortical ER. To achieve this, AtSUN2 and RETICULON1 (RTN1) were used as NE and ER markers, respectively. SUN2-YFP was predominantly localised to the NE (NE/ER ratio = 3.3 ± 0.28) while the ER marker GFP-RTN1 was predominantly ER-localised (1.2 ± 0.09) as expected. RFP-SUN3 (1.4 ± 0.09) and SUN4-RFP (1.2 ± 0.07) were also found to be significantly enriched in the ER with values very similar to RTN1 (Figure 1A). To dissect the subcellular localisation of the two mid-SUN proteins more precisely, we resolved their membrane enrichment by determining fluorescence intensity (FI) on line profiles. Line profiles as illustrated in Figure 1 and Supplementary Figure S1 were defined using calnexin-mCherry (CXN-mCherry) as a marker of the ER [21]. In this case, SUN2 fluorescence was significantly nearer the centre of the nucleus (0.067 ± 0.012 µm) in good agreement with its predominant INM localization, while SUN3 peak fluorescence was external to that of CXN-mCherry (GFP-SUN3: −0.036 ± 0.007 µm and SUN3-GFP: −0.028 ± 0.007 µm) with SUN4 colocalised with it (GFP-SUN4: 0 ± 0.009 µm) suggesting strongly that the mid-SUNs are predominantly localised at the ONM/ER (Figure 1B).

2.2. Quantifying the Effect of SUN3 and 4 OverExpression on ER Architecture

SUN3 and 4 have been implicated in maintaining ER morphology [18]. To further explore the potential function of SUN3 and 4 as ER morphogens, we compared ER morphology of control cells expressing GFP-HDEL alone (Figure 2A) with N- and C-terminally tagged SUN3 and SUN4 (Figure 2B–E) using the ER structure quantification software AnalyzER [22]. An initial analysis by MANOVA revealed significant changes in ER structure on overexpression of both SUN proteins (Pillias trace, p-value = 9.31 × 10−12 and Roy’s largest root, p-value = 1.90 × 10−22). Plotting the first two discriminate canonical variables (Figure 2F) and a dendrogram plot of group means (Figure 2G) indicated that SUN3 and SUN4 overexpression results in different modifications to ER structure. The changes induced by overexpression were significant but not large in magnitude and many of the results were within the 95% confidence limit of the control group (Figure 2F). Which terminal was tagged appeared to have limited impact on the effect of SUN3 and 4 on ER structure. ANOVA comparisons, with an applied Bonferroni correction for multiple comparisons, revealed changes in tubule length, cisternae area and shape, and polygonal region size and shape (Figure 2H). Interestingly, SUN3 and 4 appeared to have largely opposing impacts on ER structure, with SUN3 overexpression significantly increasing mean cisternae area (by at least 85.3%) whilst SUN4 overexpression reduces the mean cisternal area by at least 21.15%, though this change is not statistically significant (Supplementary Table S1).
SUN4 overexpression also caused GFP-HDEL to localise preferentially to small pockets spread throughout the ER. This phenotype is similar to what has been observed previously in overexpression of membrane constricting members of the reticulon protein family [23]. Using high-resolution confocal microscopy imaging, we were able to observe that whilst CXN-mCherry, which normally localises to ER and ONM membranes, has an even distribution of lumenal contents along its tubules (visualised by GFP-HDEL). This was not the case for RFP-SUN4, SUN4-RFP, and RFP-RTN1 (Figure 3A–D). Tracing the intensity of the luminal marker GFP-HDEL along tubules also containing RFP-SUN4 and SUN4-RFP showed inconsistent luminal contents, with patches of more or less GFP-HDEL visible. However, when compared to RFP-RTN1 overexpression, there were no long stretches of tubule completely lacking GFP-HDEL. This suggests that membrane-bound SUN4 was able to induce some ER lumen constriction resulting in GFP-HDEL patches—but not to the same extent and potential mechanism as RFP-RTN1 (Figure 3E). In addition, peaks in CXN-mCherry and GFP-HDEL show some co-occurrence along the length of tubules. The same is true of SUN4 and GFP-HDEL, whilst RFP-RTN1 and GFP-HDEL have peaks in intensity that do not co-occur, as GFP-HDEL is forced into areas with minimal RFP-RTN1 expression (Figure 3F). This suggests that areas of high SUN4 in the ER membrane do not correlate well with areas of constriction, as they do on RFP-RTN1 overexpression.

2.3. Using roGFP to Determine the Membrane Topology of SUN3 and 4

Mid-SUN proteins contain three transmembrane domains, coiled-coil regions, and a SUN domain with the CC and SUN domains involved in protein–protein interaction [13]. The membrane topology of both SUN3 and SUN4 proteins was examined using the redox sensitive fluorescent tag roGFP2 (hereafter referred to as roGFP; [24]). roGFP has two excitation peaks at 405 nm and 488 nm [24], and the efficiency of each excitation peak is dependent on the redox state of the surrounding environment. In the more oxidised environment of the ER lumen, cysteine residues in roGFP form disulphide bonds, improving excitation at 405 nm [25,26,27]. In the more reduced environment of the cytosol, excitation is most efficient at 488 nm [26,27]. Similar to previous studies, cytosolic roGFP transiently expressed in N. benthamiana showed almost no excitation at 405 nm and significant excitation at 488 nm (Figure 4A), a ratio of 0.1 ± 0.009 (Figure 4K). When roGFP was localised to the lumen of the ER, using the ER retention signal HDEL, approximately equivalent excitation at both 405 nm and 488 nm (Figure 4B), a ratio of 1.2 ± 0.03 (Figure 4K), was observed. The roGFP-SUN3 (Figure 4C) had a similar excitation profile to cytosolic roGFP, ratio 0.2 ± 0.03, while SUN3-roGFP (Figure 4D), had a similar excitation profile to ER lumenal roGFP (ratio 8.9 ± 0.08). SUN4 followed the same patterns, with roGFP-SUN4 showing similar excitation to cytosolic roGFP (ratio 0.1 ± 0.008; Figure 4E), and SUN4-roGFP showing excitation similar to the ER lumen (ratio 1.1 ± 0.08; Figure 4F). These results indicate that the N-terminal tails of SUN3 and SUN4 are located in the cytosol, the SUN domain and coiled-coil domain are present in the ER lumen, followed by two transmembrane domains, leaving the C-terminal tail of SUN3 and SUN4 located in the ER lumen (Figure 4G).
This topology analysis implies that the SUN domain of SUN3 and SUN4 is not accessible to cytosolic and soluble nuclear proteins, whereas it can interact with proteins located in the ER lumen and nuclear periplasm via its SUN and coiled-coil domain and more limited in the cytoplasm via its N-terminus and a cytosolic loop.

2.4. Arabidopsis Mid-SUN3 Interacts with the ER-Localised Membrane Transcription Factor maMYB

To gain further insights into the possible function of SUN3 at the ER, a Membrane Yeast Two-Hybrid (MYTH) screen was performed to search for SUN3 interactors (Supplementary Table S2). A transcription factor localized at the ER, maMYB [28], was detected. maMYB contains two transmembrane domains at the N-terminus and two SANT domains at the C-terminus [28]. maMYB had the strongest and most reproducible interaction with SUN3. Weaker interactions were detected between maMYB and SUN1, as well as with SUN2 and SUN5, but no interaction was recorded with SUN4 (Figure 5A) suggesting a specific role for the interaction. Interactions were then confirmed in planta using apFRET. Transiently expressed fluorescent protein fusions of SUN proteins and maMYB were colocalised in the ER and NE. Protein–protein interactions were determined by calculating the FRET efficiency [EF] between each pair of proteins measured in ROI selected at the ER and the NE. Protein interactions were detected between CFP–maMYB or maMYB–CFP and YFP–SUN3 while no interaction was detected between maMYB and SUN1, SUN2, or SUN4 in planta (Supplementary Table S3). Despite co-localisation at the NE and ER, only SUN3 interacted with maMYB in planta whilst SUN4 did not, in agreement with the MYTH results. Deletion analysis performed in MYTH on maMYB and SUN3 suggested the interaction requires the N-terminal domains of both SUN3 and maMYB (Figure 5B,C).
To confirm the relevance of the interaction between maMYB and SUN3, co-occurrence of the two proteins was searched across a set of species chosen to represent the green lineage (Supplementary Table S4). According to the Orthologous Matrix (OMA) and EggNOG databases [29,30], SUN3 and maMYB have a large number of orthologues across plant species, although the mid-SUN domain was already present in unicellular green algae while maMYB appears later in mosses (Supplementary Figure S2). Interestingly, the transmembrane domains of maMYB and SUN3 as predicted by TMHMM-2.0 [31] are well conserved in all orthologues. Altogether, the localization and interaction between maMYB and SUN3 and the conservation of transmembrane domains during evolution suggest that the two proteins may have an important function in the ER and potentially the outer nuclear membrane.

2.5. Disruption of SUN3 and maMYB Affects Root Hair Elongation

To understand whether SUN3 and maMYB could act in the same pathway, we studied the genetic interaction between SUN3 and maMYB. To this end, a mamyb-1 mutant was generated using CRISPR/Cas9 [32] by targeting the 5′ maMYB coding sequence. Transcripts from the mamyb-1 mutant allele are produced but given the stop codon occurring after 26 amino acids, no maMYB protein was revealed by Western Blot suggesting that it is a true null mutation (Supplementary Figure S3). Homozygous mamyb-1 plants were crossed with the sun3-1 mutant [13] to obtain the double mutant, sun3-1 mamyb-1. Both single and double mutants did not show any obvious growth defects. Disruption of maMYB was previously shown to affect root hair length but not root hair initiation using RNAi lines [28]. As previously described, we did not observe significant differences in root hair numbers among the tested genotypes. Therefore, root hair elongation was investigated using a similar experimental procedure as described by Slabaugh et al. [28] (Figure 6A). In this assay, each mutant line was normalised with the wild-type Col0 line (mutant/Col0 ratio set as 1) and root hair length expressed as a ratio between a given mutant and Col0 plants grown simultaneously. Root hair length was reduced in the sun3-1 (64.46 ± 0.26%), mamyb-1 (49.38 ± 0.25%), and in sun3-1 mamyb-1 (53.11 ± 0.21%) lines while overexpression of maMYB-1-YFP under a 35S promoter had no significant impact on root hair length (89.44 ± 0.34%) (Figure 6B). These results suggested that SUN3 and maMYB are acting in the same pathway to promote root hair elongation as no additive effect was observed in double mutants versus single mutants.

3. Discussion

Knowledge of the localisation and function of mid-SUN proteins is expanding rapidly, with growing evidence from numerous studies, especially of higher plant systems. Here, we showed that mid-SUN proteins are functional constituents of the ER and ONM with luminal SUN and coiled-coil domains, as well as functions in ER shaping and specialised transcription factor binding.
Using fluorescence intensity of ER and NE markers, we showed the enrichment of SUN3 and SUN4 in the ONM and in the ER, but not in the INM (Figure 1). As the ONM is continuous with the ER through junctional regions, this is not unexpected and other ER markers (such as Calnexin) are found in both membranes. However, mid-SUNs appear to have specific functionality in both membranes as they previously have been shown to interact with (ONM)-localised KASH proteins and Cter-SUNs [13]. In yeast, the sole mid-SUN protein, SLP1, is part of a complex with the YERP65 protein recruiting the Cter-SUN protein Mps3 to the NE [12]. Whether plant mid-SUN proteins have a similar chaperoning function in transporting NE components to their destined membrane is currently unknown but an enticing concept to investigate.
While both are localised in the ER and ONM, SUN3, but not SUN4, was previously shown to have different mobile behaviours in the ER and NE indicating different protein interactions in these membranes [13]. This mirrors our findings that SUN3 but not SUN4 is able to interact with maMYB (Figure 5) providing an example of such differing interactions.
The transcription factor maMYB is a member of the R2R3-MYB subgroup of the homeodomain-like MYB superfamily. This gene family is a large, plant-specific subfamily of the MYB gene family in eukaryotes that all have DNA binding functionality in the MYB domain [33]. Arabidopsis maMYB is a membrane-bound transcription factor involved in plant development and responses to environmental stress, including cell division [34], hormone signalling [35], and root hair elongation [28]. Here, we found that maMyb and SUN3 interact both in vitro and in planta in the ER (Figure 5 and Supplementary Table S3). Evolutionary analysis of this SUN3–maMYB complex shows the presence of SUN3 through the green lineage while maMYB emerges in the land plants (Supplementary Figure S2), reinforcing the idea of a fundamental function in eukaryotic plant cells for mid-SUN proteins and a later specialised function for maMYB. This specialisation is further suggested by the fact maMYB was only shown to interact with SUN3 and not SUN4 (Figure 5 and Supplementary Table S3). To investigate whether there is a functional as well as evolutionary link between the two, we examined root hair length. It is known that silencing of maMYB using the RNAi approach results in a significant reduction of root hair length in Arabidopsis but similar root hair density compared to wild-type [28]. We were able to replicate this by showing that mamyb-1 had a severely reduced root hair length (Figure 6). We also observed a similar significant reduction in root hair length in the sun3-1 mutant compared to Col0 (Figure 6). The simultaneous inactivation of SUN3 and maMYB in the double mutant sun3-1 mamyb-1 resulted in slightly shorter root hairs compared to the single sun3-1 mutant indicating that SUN3 and MaMYB are both involved in root hair elongation. Reduced root hair length was previously observed in a triple sun3 sun4 sun5 Arabidopsis mutant, in which the phenotype could be rescued by expression of SUN4-GFP [18]. This redundancy was not observed in our study with the shorter root hair length phenotype resulting from the absence of SUN3 only.
Protein interactions are mediated by domains and motifs and so to understand better how SUN3 and SUN4 are functional in the ER and ONM, it is essential to look at their topology. The roGFP assays in this work (Figure 4) as well as previous BiFC analysis [18] suggest that mid-SUN proteins are oriented with the SUN domain, coiled-coil domain, and C-terminus in the ER and NE lumen, and with the N-terminus in the cytosol (Figure 4) with a potential for the protein to bridge opposite ER membranes. This topology analysis implies that the SUN domain of SUN3 and SUN4 is not accessible to cytosolic and soluble nuclear proteins, whereas it can interact with proteins located in the ER and NE lumen. This lends support to previous observations showing mid-SUNs to interact with KASH proteins where the luminal KASH and SUN proteins mediate this interaction [13]. In addition, it was shown that the luminal coiled-coil domain of SUN2 is required for SUN2–SUN3 interactions. It is therefore likely that the SUN2–SUN3 interactions also occur in the NE lumen mediated by their respective coiled coils.
While the bulk of the mid-SUN protein is localised in the ER/NE lumen, protein interactions are also possible on the cytoplasmic side and via the TMDs in the membrane. Indeed, protein domain deletion analysis using MYTH suggested the interaction requires the N-terminal domains of both SUN3 and maMYB (Figure 5). As the N-terminal domains of both interactors used in the study are composed of a cytosolic N-terminal tail and transmembrane domains, the interaction can take place either in the cytosol and/or in the ER membrane. The MYB domain, located toward the cytosol, does not interact with SUN3, and remains free for a proteolytic activation and transport to the nucleus.
Their functional presence in the ER indicates that mid-SUN proteins might have an effect on ER structure. Indeed, previously Xue et al. [17] showed that loss of SUN3, SUN4, and SUN5 leads to increased cisternal area, and that overexpression of SUN4 results in the formation of enlarged isolated bubbles of lumenal contents. To explore their role as potential ER morphogens, we used AnalyzER to quantify effects on ER structure. Our results demonstrated that SUN3 and SUN4 have different, opposing effects on ER structure, with SUN3 inducing an increase in cisternal area and SUN4 resulting in an abundance of tubules with a partial constriction phenotype (Figure 2). Climp63, a mammalian lumenal spacer, has the effect of inducing cisternae proliferation on overexpression and acts to stabilise across larger cisternae [36]. This is similar to the impact of SUN3 on ER structure and of potential interest as no plant Climp63 has been reported so far. However, mammalian Climp63 localises preferentially to ER cisternae, as tubule formation by membrane bending is reliant on the reticulon proteins that shape the ER through positive membrane curvature.
The partial constriction phenotype of SUN4 is particularly interesting. It was previously shown that a loss of sun3 sun4 sun5 could not be compensated for by RTN overexpression [18] suggesting that the function of the mid-SUN proteins in membrane shaping is different to known ER morphogens. Our results further support this by showing that SUN4 overexpression, while appearing to cause a mild constriction phenotype, does not act in a similar manner to RTN1 (Figure 3). The magnitude of the modifications to ER structure also suggest that although SUN3 and 4 act to shape the ER, their impact when compared to known plant ER morphogens such as the Lunapark proteins [37], the reticulon protein family [23], and ROOT HAIR DEFECTIVE 3 (RHD3) [38] is relatively low. However, this more subtle effect on membrane morphology paired with their functional presence at the NE may indicate that the mid-SUNs could have potential functions in NE morphology and dynamics.
Taken together, this work suggests important roles for mid-SUN proteins in the ER and ONM. Studies to date, including this, reveal interactions with Cter-SUNs and KASH proteins (potentially at the ONM [13]), a subtle but distinctive role in ER morphology, interactions with a membrane transcription factor, the unfolded protein response pathway (Jaiswal et al., 2014) in pollen tube and root hair growth, and a chaperone complex guarding the ER sorting of LRR receptor [17]. Given their primitive evolutionary origins, it is a diverse set of functions and it will be interesting to unravel the molecular processes underlying them, especially for the roles of plant mid-SUN proteins.

4. Materials and Methods

4.1. Plant Material and Growth Conditions

Arabidopsis thaliana Col0 plants were grown from seeds surface sterilized using 0.01% SDS and 70% ethanol for 5 min, and washed twice with pure ethanol, then dried on Whatman paper. Dried seeds were transferred to half strength Murashige and Skoog media including vitamins, and supplemented with 1% sucrose, 1% Phytagel, and adjusted at pH 5.7, then stratified for 3 days at 4 °C in the dark. Arabidopsis seedlings were grown at 23 °C, 16 h photoperiod under illumination of 100 μE. N. benthamiana plants were grown on soil in a greenhouse in the same conditions.

4.2. Accessions Numbers, T-DNA Mutants and 35S::MaMYB-vYFP Line

Accession numbers for this study include SUN1 (AT5G04990), SUN2 (AT3G10730), SUN3 (AT1G22882), SUN4 (AT1G71360), SUN5 (AT4G24950), and maMYB (AT5G45420). The sun3 insertion T-DNA line Flag_374A03 (sun3-1) from the Salk collection was obtained from the European Arabidopsis Stock Centre. As Flag_374A03 line was in the Ws (Wassilewskija) background [13], three successive generations of outcross with wild-type Col0 plants were used to generate a near isogenic Flag_374A03 Col0 line. The T-DNA insertion was genotyped by PCR. A line overexpressing maMYB fused upstream to venus YFP under the 35S cauliflower promoter was provided by Prof F Brandizzi (maMYB-vYFP) [28].

4.3. Molecular Cloning and CRISPR/Cas9-Mediated Genome Editing

Arabidopsis SUN3 and SUN4 cDNA were cloned into the N-terminal roGFP2 (hereafter referred to as roGFP) vector pCM01 and C-terminal roGFP vector pSS01 created by Brach et al. (2009) using Gateway cloning technology (Invitrogen, Waltham, MA, USA).
For fluorescent fusions SUN3 and SUN4 cDNA was cloned into the vectors pB7RWG2 (C-terminal mRFP), pB7WGR2 (N-terminal mRFP), pB7FWG2 (C-terminal eGFP), pB7WGF2 (N-terminal eGFP) [39].
The mamyb-1 mutant was generated using the CRISPR/Cas9 system [32]. A 20-nucleotide guide targeting the 5′ of the mamyb coding sequence (5′-GCAGAAGATGAACTCCGTTC-3′) was selected using Crispr-Plant (http://omap.org/crispr2/) [40] and cloned into pDe-Cas9 by Gateway LR clonase (Thermo Fisher Scientific, Waltham, MA, USA). The new binary vector was transformed into Agrobacterium tumefaciens (GV3101) and then into A. thaliana Col0 by the floral-dip method [41]. The selected knockdown mamyb-1 mutant line results from an 8 nucleotides deletion leading to a premature stop codon and a truncated protein of 26 amino acids (Figure 1). This CRISPR-Cas9 generated mutant was crossed with the Flag_374A03 Col0 line to generate the sun3-1 mamyb-1 double mutant. The CRISPr/Cas9 mutants are genotyped by PCR amplification and sequencing.

4.4. Subcellular Localisation of SUN3 and 4 Analysed via Proximity Distance

The subcellular localisation of SUN3 and SUN4 was measured using radial line profiles and confocal microscopy. Briefly, z-stacks of nuclei were collected and the medial point of each nucleus was selected for analysis, to limit the impact of membrane curvature on the analysis. Using the Radial Profile Extended plugin on ImageJ, ROIs of 30° for each nucleus were drawn and the fluorescence intensity (FI) was measured along each radial line profile. FI values of all the lines were averaged to produce an average FI plotted as distances (expressed in µm) along the radius. The mean distance for the peak intensities of SUN2-CFP, GFP-SUN3, SUN3-GFP, and GFP-SUN4 was then subtracted with the peak intensity of the ONM/ER marker calnexin-(CXN) mCherry to give proximity distances (in µm). Negative distance referred to markers enriched in ER while positive values are enriched in the nucleus. Significant differences between group means were determined first by a Kruskal–Wallis H test followed by Dunn’s test with Bonferroni correction as a post-hoc test, χ2(6) = 36.078, p = 2.662 × 10−6.

4.5. Quantifying Changes in ER Structure

N. benthamiana was infiltrated as described by [42]. Briefly, transformed Agrobacteria (GV3101) were grown overnight at room temperature. Cultures were pelleted by centrifugation at 2200× g for 3 min at room temperature, resuspended in infiltration buffer (50 mM MES, 2 mM Na3PO4.12H2O, 5mg mL−1 glucose, and 0.1 mM acetosyringone) and diluted to the appropriate optical density (OD600 0.15 for SUN3 and 4 roGFP constructs; OD600 0.1 for roGFP-HDEL and roGFP cytosolic constructs). All constructs were infiltrated alongside the silencing inhibitor P19 (OD600 = 0.1) (Lakatos et al., 2004) to improve expression. After infiltration by injection, plants were transferred to an incubator and maintained at 23 °C for 72 h before imaging.
ER structure analysis was performed on images collected using confocal microscopy on N. benthamiana transiently expressing fluorescently tagged proteins of interest. Analysis of ER structure was performed using the AnalyzER software, as described in [22]. In brief, confocal images were collected and imported into AnalyzER. Imaging conditions, including pixel size and channel order, were input into the software, together with estimated minimum and maximum tubule width. Images were upsampled to increase pixel density and to prevent pixelation error later in the analysis. A region of background outside the cell was measured, with the intensity of this region subtracted from the overall images. Image intensity was then enhanced by deploying a guided filter, to increase the intensity of the ER network. The area of analysis was masked out using thresholding, and within this region cisternae were identified using an opening function, followed by active contouring, to map the outline of cisternae onto the underlying image intensity. Any cisternae with an area of less than 0.3 µm2 were reclassified as puncta as defined by [43]. The tubular network was then enhanced using an edge preserving phase congruency filter [44]. The width of each tubule was calculated using integrated intensity granulometry. Multiple structure metrics were then calculated for each specific subregion of the ER, including: tubule length, the length of each tubule between junction points and cisternae in the ER; cisternae area, the area of each cisternae; polygonal region area, the area of each polygonal region within the ER, where the tubules have been thinned to a single pixel wide skeleton; elongation, the ratio of the major and minor axes; circularity, the ratio of the radius measured for the area of the region to the radius measured from the perimeter of the region and roughness, the ratio of the regions perimeter squared, to the area.

4.6. Tubule Morphology Quantification

Tubule morphology was calculated by measuring the diameter of tubules along their length, with the position determined by the tubular skeleton. The intensity of pixels along the tubular skeleton was measured and peaks in intensity identified. Troughs were identified using the same method, using an inverted intensity profile. Once peaks and troughs had been identified, the width of the tubule at that point was measured and the distance between peaks calculated. Using this data, it is possible to calculate the following measurements: peak/trough width, the estimated width of each peak/trough; peak/trough intensity, the mean intensity at the peak/trough; peak/trough covariance, the mean covariance of intensities of both channels, measured at the peak/trough intensity of each tubule. Positive values are an indicator of similar behaviour, whilst negative values suggest that the tubule intensities move in opposite directions.

4.7. Determination of SUN3 and SUN4 Topology Using roGFP

Samples of the leaf epidermis were imaged using a Zeiss PlanApo x100/1.46 NA oil immersion objective on a Zeiss 880 confocal microscope with Airyscan detector. Samples were excited using 405 and 488 nm lasers with line switching. The detection wavelength was set to between 500–588 nm. For improved signal to noise ratio, images were collected using 4× averaging. The mean intensity of an image was quantified by measuring the intensity in an ROI on the ER (or cytoplasm in the case of the cytoplasmic control) using FIJI. A region of background intensity was also measured in the area surrounding the ER, or outside the cell when using cytoplasmic controls. The background intensity was subtracted from the measured fluorophore intensity in both channels. Then, the intensity of the 405 nm image was divided by the intensity of the 488 nm image.

4.8. Membrane Yeast Two-Hybrid (MYTH)

Protein interaction of membrane proteins SUN3 and maMYB were performed using Split-Ubiquitin-based Membrane Yeast Two-Hybrid (MYTH) system as previously described [45]. Yeast media used were a standard Yeast Nitrogen Base (YNB) supplemented with amino acids and bases as required. Interactions were tested on Permissive Medium (PM: YNB without Leu and Trp) and Test Medium (TM: YNB without Leu, Trp, Ade, and His) using the NMY51 strain (MATa, his3Δ200, trp1-901, leu2-3,112, ade2, LYS2::(lexAop)4-HIS3, ura3::(lexAop)8-lacZ, ade2::(lexAop)8-ADE2, GAL4) and the following prey plasmid pPR3N (, TRP1, AmpR), and bait plasmid pBT3N (CEN, LEU2, KanR). Full-length SUN constructs were as in [13]. SUN and MaMYB deletion constructs used for the MYTH were generated by ‘gap-repair’ homologous recombination in yeast [46] from cDNAs amplified using chimeric primer pairs containing 5′ ends with 35 bp of homology to the linearized yeast target plasmid. After digestion by SfI1, bait and prey plasmids and PCR-amplified cDNAs were cotransformed in yeast with a vector/insert of 1/3. All constructs were checked and sequenced. Screen for new SUN3 interacting proteins was performed using a cDNA library (DUALSYSTEM Biotech) cloned into the prey vector pDSL-Nx (, TRP1, AmpR) as described in [13].

4.9. Conservation of MaMYB and mid-SUN3 and Phylogenetic Reconstruction

maMYB and SUN3 protein orthologues were selected in 16 species from green algae, basal land plants and angiosperms: Arabidopsis thaliana, Arabidopsis lyrata, Brassica oleracea, Populus trichocarpa, Theobroma cacao, Prunus persica, Glycine max, Oryza sativa sub sp indica, Musa acuminata, Zea mays, Amborella trichopoda, Picea glauca, Selaginella moelledorffii, Physcomitrella patens, Chlamydomonas reinhardtii, Ostreococcus lucimarinus. Orthologues were collected using OMA and EggNOG databases [29,30]. Graphic representation was performed using ITOL [47]. Signal peptide were predicted by SignalP-5.0 [48] and Phobius [49], transmembrane domain by TMHMM-2.0 [32]. The accession numbers of all orthologues are provided in Supplementary Table S4.

4.10. Root Hair Measurements

Arabidopsis thaliana wild-type (Col0), 3 mutant lines (sun3-1, mamyb-1, sun3-1 mamyb-1) and a 35S::maMYB-YFP line were used for the root hair analysis. The 11-day-old seedlings, grown upright, were imaged using an Olympus XC30 digital camera mounted on an Olympus SZX12 stereomicroscope. Digital images were analysed with ImageJ software. Root hair measurements were taken according to Slabaugh et al. [28]. Root hairs were measured within the differentiation zone approximately 5 mm from the root tip and for each plant the five longest root hairs in focus were measured in five independent replicates.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12091787/s1, Figure S1: Co-localisation of CXN-mCherry with AtSUN2-CFP, AtSUN3-GFP, GFP-AtSUN3 and GFP-AtSUN4; Figure S2: SUN3 is present all across the green lineage when maMYB is only conserved in land plants; Figure S3: Characterisation of maMYB-1 mutant generated by CRISPR/Cas9; Table S1: Percentage change in ER structure variables; Table S2: Screen for SUN3 interactors using Membrane Yeast Two-Hybrid (MYTH) system; Table S3: FRET efficiency (EF) indicates interactions, or lack thereof. Positive protein–protein interactions are displayed in green; Table S4: Protein orthologues of SUN3 and maMYB in 16 plant species.

Author Contributions

A.B.-S., C.T., E.V., D.E.E., V.K. and K.G. supervised the study. E.V., M.V. and S.M. performed MYTH experiments. E.V., M.V. and S.T. established and characterized the mutant lines. A.B.-S. and C.C. performed the root hair experiments; C.T. and C.C. the bioinformatic analysis. B.A. carried out FRET, ER/NE ratio, and subcellular localization studies. C.P. carried out AnalyzER and topology analysis. A.B.-S., C.T., D.E.E., V.K. and K.G. contributed to the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Centre National de la Recherche Scientifique (CNRS), Institut National de la Santé et de la Recherche Médicale (INSERM), Université Clermont Auvergne (UCA), 16-IDEX-0001 CAP20-25 CIR 1 and the COST-Action INDEPTH (CA16212) to A.B.-S., C.T., C.C., E.V., S.M. and S.T. B.A. was supported by a Groome PhD studentship from Oxford Brookes University.

Data Availability Statement

All authors agree with MDPI Research Data Policies.

Acknowledgments

The authors thank Federica Brandizzi for providing the 35S-maMYB-vYFP transgenic lines that were used as controls in these experiments. C.T. and A.B.-S. thank the CLIC microscopy platform for their help in collecting root hair images.

Conflicts of Interest

The author declare no conflict of interest.

References

  1. Graumann, K.; Runions, J.; Evans, D.E. Characterization of SUN-domain Proteins at the Higher Plant Nuclear Envelope. Plant J. 2010, 61, 134–144. [Google Scholar] [CrossRef] [PubMed]
  2. Murphy, S.P.; Simmons, C.R.; Bass, H.W. Structure and Expression of the Maize (Zea mays L.) SUN-Domain Protein Gene Family: Evidence for the Existence of Two Divergent Classes of SUN Proteins in Plants. BMC Plant Biol. 2010, 10, 269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Oda, Y.; Fukuda, H. Dynamics of Arabidopsis SUN Proteins during Mitosis and Their Involvement in Nuclear Shaping. Plant J. 2011, 66, 629–641. [Google Scholar] [CrossRef] [PubMed]
  4. Zhou, Z.; Du, X.; Cai, Z.; Song, X.; Zhang, H.; Mizuno, T.; Suzuki, E.; Yee, M.R.; Berezov, A.; Murali, R.; et al. Structure of Sad1-UNC84 Homology (SUN) Domain Defines Features of Molecular Bridge in Nuclear Envelope. J. Biol. Chem. 2012, 287, 5317–5326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Tamura, K.; Iwabuchi, K.; Fukao, Y.; Kondo, M.; Okamoto, K.; Ueda, H.; Nishimura, M.; Hara-Nishimura, I. Myosin XI-i Links the Nuclear Membrane to the Cytoskeleton to Control Nuclear Movement and Shape in Arabidopsis. Curr. Biol. 2013, 23, 1776–1781. [Google Scholar] [CrossRef] [Green Version]
  6. Zhou, X.; Groves, N.R.; Meier, I. Plant Nuclear Shape Is Independently Determined by the SUN-WIP-WIT2-Myosin XI-i Complex and CRWN1. Nucleus 2015, 6, 144–153. [Google Scholar] [CrossRef] [Green Version]
  7. Murphy, S.P.; Gumber, H.K.; Mao, Y.; Bass, H.W. A Dynamic Meiotic SUN Belt Includes the Zygotene-Stage Telomere Bouquet and Is Disrupted in Chromosome Segregation Mutants of Maize (Zea mays L.). Front. Plant Sci. 2014, 5, 314. [Google Scholar] [CrossRef] [Green Version]
  8. Varas, J.; Graumann, K.; Osman, K.; Pradillo, M.; Evans, D.E.; Santos, J.L.; Armstrong, S.J. Absence of SUN1 and SUN2 Proteins in Arabidopsis Thaliana Leads to a Delay in Meiotic Progression and Defects in Synapsis and Recombination. Plant J. 2015, 81, 329–346. [Google Scholar] [CrossRef] [Green Version]
  9. Pradillo, M.; Evans, D.; Graumann, K. The Nuclear Envelope in Higher Plant Mitosis and Meiosis. Nucleus 2019, 10, 55–66. [Google Scholar] [CrossRef] [Green Version]
  10. Sohaskey, M.L.; Jiang, Y.; Zhao, J.J.; Mohr, A.; Roemer, F.; Harland, R.M. Osteopotentia Regulates Osteoblast Maturation, Bone Formation, and Skeletal Integrity in Mice. J. Cell Biol. 2010, 189, 511–525. [Google Scholar] [CrossRef] [Green Version]
  11. Shimada, N.; Inouye, K.; Sawai, S.; Kawata, T. SunB, a Novel Sad1 and UNC-84 Domain-Containing Protein Required for Development of Dictyostelium Discoideum. Dev. Growth Differ. 2010, 52, 577–590. [Google Scholar] [CrossRef] [PubMed]
  12. Friederichs, J.M.; Gardner, J.M.; Smoyer, C.J.; Whetstine, C.R.; Gogol, M.; Slaughter, B.D.; Jaspersen, S.L. Genetic Analysis of Mps3 SUN Domain Mutants in Saccharomyces Cerevisiae Reveals an Interaction with the SUN-Like Protein Slp1. G3 2012, 2, 1703–1718. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Graumann, K.; Vanrobays, E.; Tutois, S.; Probst, A.V.; Evans, D.E.; Tatout, C. Characterization of Two Distinct Subfamilies of SUN-Domain Proteins in Arabidopsis and Their Interactions with the Novel KASH-Domain Protein AtTIK. J. Exp. Bot. 2014, 65, 6499–6512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Jaiswal, D.K.; Mishra, P.; Subba, P.; Rathi, D.; Chakraborty, S.; Chakraborty, N. Membrane-Associated Proteomics of Chickpea Identifies Sad1/UNC-84 Protein (CaSUN1), a Novel Component of Dehydration Signaling. Sci. Rep. 2014, 4, 4177. [Google Scholar] [CrossRef] [Green Version]
  15. Poulet, A.; Probst, A.V.; Graumann, K.; Tatout, C.; Evans, D. Exploring the Evolution of the Proteins of the Plant Nuclear Envelope. Nucleus 2017, 8, 46–59. [Google Scholar] [CrossRef] [Green Version]
  16. Vasnier, C.; de Muyt, A.; Zhang, L.; Tessé, S.; Kleckner, N.E.; Zickler, D.; Espagne, E. Absence of SUN-Domain Protein Slp1 Blocks Karyogamy and Switches Meiotic Recombination and Synapsis from Homologs to Sister Chromatids. Proc. Natl. Acad. Sci. USA 2014, 111, E4015–E4023. [Google Scholar] [CrossRef] [Green Version]
  17. Xue, Y.; Meng, J.-G.; Jia, P.-F.; Zhang, Z.-R.; Li, H.-J.; Yang, W.-C. POD1-SUN-CRT3 Chaperone Complex Guards the ER Sorting of LRR Receptor Kinases in Arabidopsis. Nat. Commun. 2022, 13, 2703. [Google Scholar] [CrossRef]
  18. Xue, Y.; Jia, P.-F.; Li, H.-J. SUN3/4/5 Proteins Regulate Endoplasmic Reticulum Tubule Formation and Luminal Spacing in Arabidopsis. J. Genet. Genom. 2022; Online ahead of print. [Google Scholar] [CrossRef]
  19. Zhang, S.; Xu, C.; Larrimore, K.E.; Ng, D.T.W. Slp1-Emp65: A Guardian Factor That Protects Folding Polypeptides from Promiscuous Degradation. Cell 2017, 171, 346–357.e12. [Google Scholar] [CrossRef] [Green Version]
  20. Yuan, L.; Pan, J.; Zhu, S.; Li, Y.; Yao, J.; Li, Q.; Fang, S.; Liu, C.; Wang, X.; Li, B.; et al. Evolution and Functional Divergence of SUN Genes in Plants. Front. Plant Sci. 2021, 12, 646622. [Google Scholar] [CrossRef]
  21. Liu, D.Y.T.; Smith, P.M.C.; Barton, D.A.; Day, D.A.; Overall, R.L. Characterisation of Arabidopsis Calnexin 1 and Calnexin 2 in the Endoplasmic Reticulum and at Plasmodesmata. Protoplasma 2017, 254, 125–136. [Google Scholar] [CrossRef] [PubMed]
  22. Pain, C.; Kriechbaumer, V.; Kittelmann, M.; Hawes, C.; Fricker, M. Quantitative Analysis of Plant ER Architecture and Dynamics. Nat. Commun. 2019, 10, 984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sparkes, I.; Tolley, N.; Aller, I.; Svozil, J.; Osterrieder, A.; Botchway, S.; Mueller, C.; Frigerio, L.; Hawes, C. Five Arabidopsis Reticulon Isoforms Share Endoplasmic Reticulum Location, Topology, and Membrane-Shaping Properties. Plant Cell 2010, 22, 1333–1343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Dooley, C.T.; Dore, T.M.; Hanson, G.T.; Jackson, W.C.; Remington, S.J.; Tsien, R.Y. Imaging Dynamic Redox Changes in Mammalian Cells with Green Fluorescent Protein Indicators. J. Biol. Chem. 2004, 279, 22284–22293. [Google Scholar] [CrossRef] [Green Version]
  25. Hanson, G.T.; Aggeler, R.; Oglesbee, D.; Cannon, M.; Capaldi, R.A.; Tsien, R.Y.; Remington, S.J. Investigating Mitochondrial Redox Potential with Redox-Sensitive Green Fluorescent Protein Indicators. J. Biol. Chem. 2004, 279, 13044–13053. [Google Scholar] [CrossRef] [Green Version]
  26. Meyer, A.J.; Brach, T.; Marty, L.; Kreye, S.; Rouhier, N.; Jacquot, J.-P.; Hell, R. Redox-Sensitive GFP in Arabidopsis Thaliana Is a Quantitative Biosensor for the Redox Potential of the Cellular Glutathione Redox Buffer. Plant J. 2007, 52, 973–986. [Google Scholar] [CrossRef] [Green Version]
  27. Schwarzländer, M.; Fricker, M.D.; Müller, C.; Marty, L.; Brach, T.; Novak, J.; Sweetlove, L.J.; Hell, R.; Meyer, A.J. Confocal Imaging of Glutathione Redox Potential in Living Plant Cells. J. Microsc. 2008, 231, 299–316. [Google Scholar] [CrossRef]
  28. Slabaugh, E.; Held, M.; Brandizzi, F. Control of Root Hair Development in Arabidopsis Thaliana by an Endoplasmic Reticulum Anchored Member of the R2R3-MYB Transcription Factor Family. Plant J. 2011, 67, 395–405. [Google Scholar] [CrossRef]
  29. Huerta-Cepas, J.; Szklarczyk, D.; Heller, D.; Hernández-Plaza, A.; Forslund, S.K.; Cook, H.; Mende, D.R.; Letunic, I.; Rattei, T.; Jensen, L.J.; et al. EggNOG 5.0: A Hierarchical, Functionally and Phylogenetically Annotated Orthology Resource Based on 5090 Organisms and 2502 Viruses. Nucleic Acids Res. 2019, 47, D309–D314. [Google Scholar] [CrossRef]
  30. Altenhoff, A.M.; Train, C.-M.; Gilbert, K.J.; Mediratta, I.; Mendes de Farias, T.; Moi, D.; Nevers, Y.; Radoykova, H.-S.; Rossier, V.; Warwick Vesztrocy, A.; et al. OMA Orthology in 2021: Website Overhaul, Conserved Isoforms, Ancestral Gene Order and More. Nucleic Acids Res. 2021, 49, D373–D379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Möller, S.; Croning, M.D.R.; Apweiler, R. Evaluation of Methods for the Prediction of Membrane Spanning Regions. Bioinformatics 2001, 17, 646–653. [Google Scholar] [CrossRef] [PubMed]
  32. Fauser, F.; Schiml, S.; Puchta, H. Both CRISPR/Cas-Based Nucleases and Nickases Can Be Used Efficiently for Genome Engineering in Arabidopsis Thaliana. Plant J. 2014, 79, 348–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Dubos, C.; Stracke, R.; Grotewold, E.; Weisshaar, B.; Martin, C.; Lepiniec, L. MYB Transcription Factors in Arabidopsis. Trends Plant Sci. 2010, 15, 573–581. [Google Scholar] [CrossRef]
  34. Kim, Y.-S.; Kim, S.-G.; Park, J.-E.; Park, H.-Y.; Lim, M.-H.; Chua, N.-H.; Park, C.-M. A Membrane-Bound NAC Transcription Factor Regulates Cell Division in Arabidopsis. Plant Cell 2006, 18, 3132–3144. [Google Scholar] [CrossRef]
  35. Park, J.; Kim, Y.-S.; Kim, S.-G.; Jung, J.-H.; Woo, J.-C.; Park, C.-M. Integration of Auxin and Salt Signals by the NAC Transcription Factor NTM2 during Seed Germination in Arabidopsis. Plant Physiol. 2011, 156, 537–549. [Google Scholar] [CrossRef] [Green Version]
  36. Shibata, Y.; Shemesh, T.; Prinz, W.A.; Palazzo, A.F.; Kozlov, M.M.; Rapoport, T.A. Mechanisms Determining the Morphology of the Peripheral ER. Cell 2010, 143, 774–788. [Google Scholar] [CrossRef] [Green Version]
  37. Kriechbaumer, V.; Breeze, E.; Pain, C.; Tolmie, F.; Frigerio, L.; Hawes, C. Arabidopsis Lunapark Proteins Are Involved in ER Cisternae Formation. New Phytol. 2018, 219, 990–1004. [Google Scholar] [CrossRef] [Green Version]
  38. Stefano, G.; Brandizzi, F. Unique and Conserved Features of the Plant ER-Shaping GTPase RHD3. Cell. Logist. 2014, 4, e28217. [Google Scholar] [CrossRef] [Green Version]
  39. Karimi, M.; Depicker, A.; Hilson, P. Recombinational Cloning with Plant Gateway Vectors. Plant Physiol. 2007, 145, 1144–1154. [Google Scholar] [CrossRef]
  40. Minkenberg, B.; Zhang, J.; Xie, K.; Yang, Y. CRISPR-PLANT v2: An Online Resource for Highly Specific Guide RNA Spacers Based on Improved off-Target Analysis. Plant Biotechnol. J. 2019, 17, 5–8. [Google Scholar] [CrossRef] [Green Version]
  41. Clough, S.J.; Bent, A.F. Floral Dip: A Simplified Method for Agrobacterium-Mediated Transformation of Arabidopsis Thaliana. Plant J. 1998, 16, 735–743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Sparkes, I.A.; Runions, J.; Kearns, A.; Hawes, C. Rapid, Transient Expression of Fluorescent Fusion Proteins in Tobacco Plants and Generation of Stably Transformed Plants. Nat. Protoc. 2006, 1, 2019–2025. [Google Scholar] [CrossRef]
  43. Sparkes, I.A.; Frigerio, L.; Tolley, N.; Hawes, C. The Plant Endoplasmic Reticulum: A Cell-Wide Web. Biochem. J. 2009, 423, 145–155. [Google Scholar] [CrossRef] [PubMed]
  44. Kovesi, P. Phase Congruency: A Low-Level Image Invariant. Psychol. Res. 2000, 64, 136–148. [Google Scholar] [CrossRef] [PubMed]
  45. Voisin, M.; Vanrobays, E.; Tatout, C. Investigation of Nuclear Periphery Protein Interactions in Plants Using the Membrane Yeast Two-Hybrid (MbY2H) System. Methods Mol. Biol. 2018, 1840, 221–235. [Google Scholar] [CrossRef]
  46. Oldenburg, K.R.; Vo, K.T.; Michaelis, S.; Paddon, C. Recombination-Mediated PCR-Directed Plasmid Construction in Vivo in Yeast. Nucleic Acids Res. 1997, 25, 451–452. [Google Scholar] [CrossRef] [Green Version]
  47. Letunic, I.; Bork, P. Interactive Tree of Life (ITOL) v5: An Online Tool for Phylogenetic Tree Display and Annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
  48. Nielsen, H.; Tsirigos, K.D.; Brunak, S.; von Heijne, G. A Brief History of Protein Sorting Prediction. Protein J. 2019, 38, 200–216. [Google Scholar] [CrossRef] [Green Version]
  49. Käll, L.; Krogh, A.; Sonnhammer, E.L.L. A Combined Transmembrane Topology and Signal Peptide Prediction Method. J. Mol. Biol. 2004, 338, 1027–1036. [Google Scholar] [CrossRef]
Figure 1. NE/ER fluorescence intensity (FI) ratios of ER and NE proteins. (A) NE/ER FI as measured for NE marker SUN2-YFP, ER marker GFP-RTN1, as well as for RFP-SUN3 and SUN4-RFP, when transiently coexpressed with p19 in N. benthamiana leaf epidermal cells. There were 30 cells analysed for each construct tested. There were statistically significant differences between group means as determined by Kruskal–Wallis H test, χ2(3) = 48.95, p = 1.377 × 10−10. (B) Proximity distances (µm) of fluorescent fusion proteins were computed with respect to the ER/ONM marker protein CXN-mCherry used as a reference (distance = 0 µm). SUN2-CFP (0.067 ± 0.012 µm), GFP-SUN3 (−0.036 ± 0.007 µm), SUN3-GFP (−0.028 ± 0.007 µm), and GFP-SUN4: (0 ± 0.009 µm) were measured by collecting 20 nuclei from each construct in three biological replicates.
Figure 1. NE/ER fluorescence intensity (FI) ratios of ER and NE proteins. (A) NE/ER FI as measured for NE marker SUN2-YFP, ER marker GFP-RTN1, as well as for RFP-SUN3 and SUN4-RFP, when transiently coexpressed with p19 in N. benthamiana leaf epidermal cells. There were 30 cells analysed for each construct tested. There were statistically significant differences between group means as determined by Kruskal–Wallis H test, χ2(3) = 48.95, p = 1.377 × 10−10. (B) Proximity distances (µm) of fluorescent fusion proteins were computed with respect to the ER/ONM marker protein CXN-mCherry used as a reference (distance = 0 µm). SUN2-CFP (0.067 ± 0.012 µm), GFP-SUN3 (−0.036 ± 0.007 µm), SUN3-GFP (−0.028 ± 0.007 µm), and GFP-SUN4: (0 ± 0.009 µm) were measured by collecting 20 nuclei from each construct in three biological replicates.
Plants 12 01787 g001
Figure 2. Overexpression of SUN3 and SUN4 affect ER structure. (AE) Transient overexpression of GFP-HDEL (A) shown in green and (B) RFP-SUN3, (C) SUN3-RFP, (D) RFP-SUN4, and (E) SUN4-RFP shown in magenta. (F) Pairwise scatterplot of the first two canonical variables from the MANOVA analysis grouped by expressed proteins, the dotted line represents the 95% confidence limit. (G) Dendrogram of class means using MANOVA. (H) ANOVA results of planned comparisons of ER structure variables after ER structure was modified by SUN3 and SUN4 overexpression, compared to control ER structure. The F statistic and p-values (with applied Bonferroni correction) are shown. Images collected from five biological repeats, control, n = 35; RFP-SUN3, n = 13; SUN3-RFP, n = 5; RFP-SUN4, n = 26; SUN4-RFP, n = 27. Scale bars 5 µm.
Figure 2. Overexpression of SUN3 and SUN4 affect ER structure. (AE) Transient overexpression of GFP-HDEL (A) shown in green and (B) RFP-SUN3, (C) SUN3-RFP, (D) RFP-SUN4, and (E) SUN4-RFP shown in magenta. (F) Pairwise scatterplot of the first two canonical variables from the MANOVA analysis grouped by expressed proteins, the dotted line represents the 95% confidence limit. (G) Dendrogram of class means using MANOVA. (H) ANOVA results of planned comparisons of ER structure variables after ER structure was modified by SUN3 and SUN4 overexpression, compared to control ER structure. The F statistic and p-values (with applied Bonferroni correction) are shown. Images collected from five biological repeats, control, n = 35; RFP-SUN3, n = 13; SUN3-RFP, n = 5; RFP-SUN4, n = 26; SUN4-RFP, n = 27. Scale bars 5 µm.
Plants 12 01787 g002
Figure 3. Comparison of ER constriction induced by SUN4 and RTN1 overexpression. High-resolution representative images of transmembrane proteins (magenta) and the ER luminal marker GFP-HDEL (green), alongside intensity traces along a single tubule, highlighted around the edge in cyan and along the centre of the tubule with a white dashed line. Images and line profiles shown for (A) CXN-mCherry, (B) RFP-SUN4, (C) SUN4-RFP, and (D) RFP-RTN1. Comparison of the mean trough intensity of the GFP-HDEL channel is shown (E), and the covariance of intensity peaks along all tubules in each image (F). Comparisons between groups using the Tukey HSD shown by asterisks, with *** representing p-value ≤ 0.001, ** representing p-value = 0.01–0.001, and * representing p-value = 0.05–0.01. Nonsignificant differences between the groups are denoted with ns. Results taken from three biological repeats with CXN-mCherry, n = 9; RFP-SUN4, n = 12; SUN4-RFP, n = 10; and RFP-RTN1, n = 9 cells analysed. Scale bars 2 µm.
Figure 3. Comparison of ER constriction induced by SUN4 and RTN1 overexpression. High-resolution representative images of transmembrane proteins (magenta) and the ER luminal marker GFP-HDEL (green), alongside intensity traces along a single tubule, highlighted around the edge in cyan and along the centre of the tubule with a white dashed line. Images and line profiles shown for (A) CXN-mCherry, (B) RFP-SUN4, (C) SUN4-RFP, and (D) RFP-RTN1. Comparison of the mean trough intensity of the GFP-HDEL channel is shown (E), and the covariance of intensity peaks along all tubules in each image (F). Comparisons between groups using the Tukey HSD shown by asterisks, with *** representing p-value ≤ 0.001, ** representing p-value = 0.01–0.001, and * representing p-value = 0.05–0.01. Nonsignificant differences between the groups are denoted with ns. Results taken from three biological repeats with CXN-mCherry, n = 9; RFP-SUN4, n = 12; SUN4-RFP, n = 10; and RFP-RTN1, n = 9 cells analysed. Scale bars 2 µm.
Plants 12 01787 g003
Figure 4. Characterising SUN3 and SUN4 membrane topology using a redox sensitive fluorescent protein (roGFP). Representative images of N. benthamiana epidermal cells transient expressing the redox sensitive fluorophore, roGFP, (A) localised to the cytosol, (B) the ER/NE lumen and attached to the N- and C-terminal of (C,D) SUN3 and (E,F) SUN4. All images collected in the emission range of 500–550 nm, with excitation by either a 405 nm laser (cyan) or a 488 nm laser (green). (G–J) Schematic representation of the membrane topology of SUN3 and SUN4. (K) A plot of the ratio of roGFP intensity after activation by 405 nm and 488 nm. Images collected across three biological repeats, (cytosolic roGFP, n = 10; roGFP-HDEL, n = 8; roGFP-SUN3, n = 12; SUN3-roGFP, n = 8; roGFP-SUN4, n = 9; SUN4-roGFP, n = 9). Scale bars 5 µm.
Figure 4. Characterising SUN3 and SUN4 membrane topology using a redox sensitive fluorescent protein (roGFP). Representative images of N. benthamiana epidermal cells transient expressing the redox sensitive fluorophore, roGFP, (A) localised to the cytosol, (B) the ER/NE lumen and attached to the N- and C-terminal of (C,D) SUN3 and (E,F) SUN4. All images collected in the emission range of 500–550 nm, with excitation by either a 405 nm laser (cyan) or a 488 nm laser (green). (G–J) Schematic representation of the membrane topology of SUN3 and SUN4. (K) A plot of the ratio of roGFP intensity after activation by 405 nm and 488 nm. Images collected across three biological repeats, (cytosolic roGFP, n = 10; roGFP-HDEL, n = 8; roGFP-SUN3, n = 12; SUN3-roGFP, n = 8; roGFP-SUN4, n = 9; SUN4-roGFP, n = 9). Scale bars 5 µm.
Plants 12 01787 g004
Figure 5. Interaction between maMYB and SUNs proteins in MYTH. All MYTH assays have been done on permissive and test media (PM and TM, respectively). Results from drop tests and schematic representation of protein structures drawn at scale are shown. (A) MYTH assays between full-length maMYB as bait and the SUN family proteins or maMYB as preys. Empty vector is used as a negative control. No interaction was observed with SUN4. (B) MYTH assays between full-length SUN3 as bait and deletions of maMYB as preys. Empty prey vector is used as a negative control. Weaker interaction was observed with maMYB-deltaNter. (C) MYTH assays between full-length maMYB as bait and deletions of SUN3 as preys. Empty prey vector is used as a negative control.
Figure 5. Interaction between maMYB and SUNs proteins in MYTH. All MYTH assays have been done on permissive and test media (PM and TM, respectively). Results from drop tests and schematic representation of protein structures drawn at scale are shown. (A) MYTH assays between full-length maMYB as bait and the SUN family proteins or maMYB as preys. Empty vector is used as a negative control. No interaction was observed with SUN4. (B) MYTH assays between full-length SUN3 as bait and deletions of maMYB as preys. Empty prey vector is used as a negative control. Weaker interaction was observed with maMYB-deltaNter. (C) MYTH assays between full-length maMYB as bait and deletions of SUN3 as preys. Empty prey vector is used as a negative control.
Plants 12 01787 g005
Figure 6. Mutant lines for sun3-1 and mamyb-1 show shorter root hair. (A) Representative images from root hairs for each genotype. Digital images were analysed with imageJ. At least 25 roots from seedlings were investigated and the five longest root hairs were measured at approximately 5 mm to the root tip, within the differentiation zone, for each plant as in Slabaugh et al. (2007). (B) Root hairs lengths are expressed as mutant/Col0 ratio (Col0 is set at 1, red line) for 35S-maMYB (n = 170), sun3-1 (n = 95), mamyb-1 (n = 40), and the corresponding double mutant line sun3-1 mamyb-1 (n = 85). Data are representative of five experimental replicates. Scale bar: 0.5 mm. p < 0.01 (***).
Figure 6. Mutant lines for sun3-1 and mamyb-1 show shorter root hair. (A) Representative images from root hairs for each genotype. Digital images were analysed with imageJ. At least 25 roots from seedlings were investigated and the five longest root hairs were measured at approximately 5 mm to the root tip, within the differentiation zone, for each plant as in Slabaugh et al. (2007). (B) Root hairs lengths are expressed as mutant/Col0 ratio (Col0 is set at 1, red line) for 35S-maMYB (n = 170), sun3-1 (n = 95), mamyb-1 (n = 40), and the corresponding double mutant line sun3-1 mamyb-1 (n = 85). Data are representative of five experimental replicates. Scale bar: 0.5 mm. p < 0.01 (***).
Plants 12 01787 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Andov, B.; Boulaflous-Stevens, A.; Pain, C.; Mermet, S.; Voisin, M.; Charrondiere, C.; Vanrobays, E.; Tutois, S.; Evans, D.E.; Kriechbaumer, V.; et al. In Depth Topological Analysis of Arabidopsis Mid-SUN Proteins and Their Interaction with the Membrane-Bound Transcription Factor MaMYB. Plants 2023, 12, 1787. https://doi.org/10.3390/plants12091787

AMA Style

Andov B, Boulaflous-Stevens A, Pain C, Mermet S, Voisin M, Charrondiere C, Vanrobays E, Tutois S, Evans DE, Kriechbaumer V, et al. In Depth Topological Analysis of Arabidopsis Mid-SUN Proteins and Their Interaction with the Membrane-Bound Transcription Factor MaMYB. Plants. 2023; 12(9):1787. https://doi.org/10.3390/plants12091787

Chicago/Turabian Style

Andov, Bisa, Aurelia Boulaflous-Stevens, Charlotte Pain, Sarah Mermet, Maxime Voisin, Camille Charrondiere, Emmanuel Vanrobays, Sylvie Tutois, David E. Evans, Verena Kriechbaumer, and et al. 2023. "In Depth Topological Analysis of Arabidopsis Mid-SUN Proteins and Their Interaction with the Membrane-Bound Transcription Factor MaMYB" Plants 12, no. 9: 1787. https://doi.org/10.3390/plants12091787

APA Style

Andov, B., Boulaflous-Stevens, A., Pain, C., Mermet, S., Voisin, M., Charrondiere, C., Vanrobays, E., Tutois, S., Evans, D. E., Kriechbaumer, V., Tatout, C., & Graumann, K. (2023). In Depth Topological Analysis of Arabidopsis Mid-SUN Proteins and Their Interaction with the Membrane-Bound Transcription Factor MaMYB. Plants, 12(9), 1787. https://doi.org/10.3390/plants12091787

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop