# Stabilization of membrane topologies by proteinaceous remorin scaffolds

Chao Su, Marta Rodriguez-Franco, Beatrice Lace, Nils Nebel, Casandra Hernandez-Reyes, Pengbo Liang, Eija Schulze, Evgeny V Mymrikov, Nikolas M Gross, Julian Knerr, et al.

► To cite this version:

Chao Su, Marta Rodriguez-Franco, Beatrice Lace, Nils Nebel, Casandra Hernandez-Reyes, et al. Stabilization of membrane topologies by proteinaceous remorin scaffolds. Nature Communications, 2023, 14 (1), pp.323. 10.1038/s41467-023-35976-5. hal-03955426

HAL Id: hal-03955426

https://ut3-toulouseinp.hal.science/hal-03955426

Submitted on 25 Jan 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Stabilization of membrane topologies by proteinaceous remorin scaffolds

Received: 14 March 2022

Accepted: 10 January 2023

Published online: 19 January 2023

<figure><img src="image_3.png" /></figure>

Check for updates

Chao Su<fnref n="1" />, Marta Rodriguez-Franco Ⓒ<fnref n="1" />, Beatrice Lace<fnref n="1" />, Nils Nebel Ⓒ<fnref n="1" />,

Casandra Hernandez-Reyes<fnref n="1" />,², Pengbo Liang<fnref n="1" />, Eija Schulze<fnref n="1" />,

Evgeny V. Mymrikov Ⓒ²,³, Nikolas M. Gross Ⓒ²,³,⁴, Julian Knerr Ⓒ⁵,

Hong Wang Ⓒ²,⁵, Lina Siukstaite<fnref n="1" />,⁶, Jean Keller Ⓒ⁷, Cyril Libourel Ⓒ⁷,

Alexandra A. M. Fischer Ⓒ²,⁴,⁶,⁸, Katharina E. Gabor⁹, Eric Mark¹⁰, Claudia Popp¹¹,

Carola Hunte Ⓒ²,³,⁶, Wilfried Weber Ⓒ²,⁶,⁸, Petra Wendler Ⓒ¹⁰, Thomas Stanislas⁹,

Pierre-Marc Delaux⁷, Oliver Einsle¹², Robert Grosse Ⓒ²,⁵, Winfried Römer Ⓒ<fnref n="1" />,²,⁶ &

Thomas Ott Ⓒ<fnref n="1" />,²

<figure><img src="image_4.png" /></figure>

In plants, the topological organization of membranes has mainly been attributed to the cell wall and the cytoskeleton. Additionally, few proteins, such as plant-specific remorins have been shown to function as protein and lipid organizers. Root nodule symbiosis requires continuous membrane re-arrangements, with bacteria being finally released from infection threads into membrane-confined symbiosomes. We found that mutations in the symbiosis-specific SYMREM1 gene result in highly disorganized perimicrobial membranes. AlphaFold modelling and biochemical analyses reveal that SYMREM1 oligomerizes into antiparallel dimers and may form a higher-order membrane scaffolding structure. This was experimentally confirmed when expressing this and other remorins in wall-less protoplasts is sufficient where they significantly alter and stabilize *de novo* membrane topologies ranging from membrane blebs to long membrane tubes with a central actin filament. Reciprocally, mechanically induced membrane indentations were equally stabilized by SYMREM1. Taken together we describe a plant-specific mechanism that allows the stabilization of large-scale membrane conformations independent of the cell wall.

Intracellular colonization of host cells represents a characteristic feature of mutualistic associations like the root nodule symbiosis (RNS) and the arbuscular mycorrhiza symbiosis (AMS) that occur between a host plant and soil-borne bacteria, such as rhizobia, or

Glomeromycotean fungi, respectively. Molecularly, RNS-related infection is tightly controlled and requires the perception of strain-specific rhizobial lipochitooligosaccharides (Nod factors) by host LysM-type receptor-like kinases<fnref n="1" />⁻³ and of rhizobial exopolysaccharides⁴. Nod

<footnote n="1">¹Faculty of Biology, University of Freiburg, 79104 Freiburg, Germany. ²CIBSS - Centre of Integrative Biological Signalling Studies, University of Freiburg, 79104 Freiburg, Germany. ³Institute for Biochemistry and Molecular Biology, ZBMZ, Faculty of Medicine, University of Freiburg, 79104 Freiburg, Germany. ⁴Spemann Graduate School of Biology and Medicine (SGBM), University of Freiburg, 79104 Freiburg, Germany. ⁵Institute of Pharmacology, Medical Faculty, University of Freiburg, 79104 Freiburg, Germany. ⁶BIOSS - Centre for Biological Signalling Studies, University of Freiburg, 79104 Freiburg, Germany. ⁷Laboratoire de Recherche en Sciences Végétales (LRSV), Université de Toulouse, CNRS, UPS, INP Toulouse, Castanet Tolosan, France. ⁸Division of Synthetic Biology, Faculty of Biology, University of Freiburg, 79104 Freiburg, Germany. ⁹Center for Plant Molecular Biology (ZMBP), University of Tübingen, 72076 Tübingen, Germany. ¹⁰Institute of Biochemistry and Biology, Department of Biochemistry, University of Potsdam, 14476 Potsdam-Golm, Germany. ¹¹Ludwig-Maximilians-University (LMU) Munich, Institute of Genetics, 82152 Martinsried, Germany. ¹²Institute of Biochemistry, Faculty of Chemistry, University of Freiburg, 79104 Freiburg, Germany. e-mail: Thomas.Ott@biologie.uni-freiburg.de</footnote>
factor perception by the host and heteromeric receptor complex assembly trigger a symbiotic signaling cascade that results in peri-nuclear calcium spiking⁵. This calcium signature is, in turn, decoded by the calcium-calmodulin-dependent kinase CCaMK/DMI3 and the transcriptional activator CYCLOPS/IPD3⁶,⁷. Upon phosphorylation of CYCLOPS/IPD3 by CCaMK/DMI3, RNS-specific gene expression is triggered by the activation of specific transcription factors such as NODULE INCEPTION (NIN)⁸,⁹. Consequently, mutations in these genes result in the inability of the host to maintain intracellular infections and bacterial release.

Intruding rhizobia mostly infect legume hosts such as *Medicago truncatula* and *Lotus japonicus* via young, growing root hairs that swell and later curl around surface attached rhizobia to entrap them just below the root hair tip. This process is driven by cellular repolarization of the actin and microtubule cytoskeleton¹⁰⁻¹³ and involves, among others, the actin polymerizing formin protein SYFO1¹⁴. Root hair curling results in an entrapment of the symbiont in a so-called infection chamber¹⁵. A re-polarization of the cell towards the infection chamber leads to the targeted secretion of proteins and membrane constituents, which enables the formation of a membrane-surrounded tunnel called the 'infection thread' (IT) that emerges from the infection chamber, transcellularly progresses though the root cortex and finally branches inside the nodule primordium¹⁶,¹⁷. Within these primordia and later in indeterminate nodules of *M. truncatula*, rhizobia are continuously released from bulges of nodular ITs (infection droplets) into these infected cells. Upon release, rhizobia differentiate into nitrogen-fixing bacteroids that remain encapsulated by the host-derived symbiosome/peribacteroid membrane¹⁸,¹⁹. It should be noted that pre-infection membrane invaginations and constitutive IT stabilization during RNS are most likely host-driven processes. In addition, initial membrane invaginations, young IT segments around the growing tip, infection droplets as precursors of bacterial release sites, and symbiosome membranes encapsulating the nitrogen-fixing symbiont are devoid of a rigid cell wall that could provide structural support to these sites²⁰,²¹. Although this is crucial to the success of the infection, our understanding of symbiotic membrane stabilization in the absence of a primary cell wall is rather sparse.

Such roles could alternatively be maintained by oligomeric scaffold proteins like clathrins as described during endocytosis²²,²³ or Bar-domain proteins like amphisin or BIN1 in human cells²⁴. During RNS, the scaffold SYMREM1, a member of the plant-specific remorin family²⁵, recruits and stabilizes the symbiotic receptor LYK3 in membrane nanodomains²⁶. This function might be explained by remorin-induced alterations in membrane fluidity²⁷,²⁸, higher order protein oligomerization²⁹⁻³¹ or maintenance of membrane-associated and phase-separated condensates³². Receptor scaffolding additionally relies on the Flotillin protein FLOT4³³. While FLOT4 is already expressed prior to rhizobial infection, SYMREM1 expression is exclusively triggered in the presence of the symbiont²⁶,³⁴,³⁵. Furthermore, SYMREM1 localization to membrane nanodomains, but not to the plasma membrane per se, depends on FLOT4²⁶. Like all other remorins, SYMREM1 is comprised of a conserved alpha-helical C-terminal region, while the intrinsically disordered N-terminal region (IDR) of SYMREM1 is highly variable in sequence and conformation³⁶. Membrane association of remorins is mediated by the Remorin C-terminal Anchor (RemCA) and, in most cases, assisted by palmitoylation³⁶,³⁷. Furthermore, it has been shown that remorins oligomerize at the plasma membrane in planta³⁸ and can form higher order filamentous structures or protein lattices in vitro³¹.

Since remorins can alter membrane fluidity²⁸ and *symrem1* mutants largely fail to release rhizobia into nodule cells³⁵, these proteins may have greater impact on membrane topology than currently envisioned. Furthermore, SYMREM1 accumulates predominantly on cell wall-devoid symbiotic membranes such as IT tips²⁶, nodular infection droplets that precede bacterial release and the symbiosome

membrane, which surrounds the released and nitrogen-fixing rhizobia inside the nodule³⁵. Therefore, we investigated putative roles of SYMREM1 in membrane dynamics and shape.

## Results

### SYMREM1 controls symbiotic membrane topologies

To visualize symbiotic membranes in great detail and to examine their precise morphology, we labeled phosphatidylserine (PS), a central phospholipid of biological membranes, by expressing a LactC2 biosensor³⁹. This allowed clear imaging of membrane patterns at ITs, infection droplets and symbiosome membranes (Supplementary Fig. 1a-f). As previously shown and recapitulated here, all three membrane sites were also targeted by SYMREM1 (Fig. 1a-f; Supplementary Fig. 1g-i)³⁵. While ITs and infection droplet structures were mostly filled with rhizobia, we frequently observed additional empty tube-like membrane structures on enlarged ITs associated with infection droplets in nodule cortex cells of wild-type (WT) plants (Fig. 1g). In line with this accumulation of SYMREM1 at ITs, on infection droplets and on symbiosome membranes, and in contrast to WT plants (Supplementary Fig. 2a, d), *symrem1* mutants are greatly impaired in releasing bacteria to cells of the inner nodule cortex and exhibited bulky ITs, as revealed by light and electron microscopy (Supplementary Fig. 2b-c, e-f)³⁵. To assess possible alterations in membrane topology, we also expressed the PS-biosensor in the *symrem1-1* mutant. In contrast to WT plants, ITs observed in the *symrem1* mutant did not display signs of membrane tubulation (Fig. 1h). While differentiated bacteroids with a tightly aligned symbiosome membrane were observed in WT nodules (Fig. 1f, k; Supplementary Fig. 2g, j), symbiosome membranes were either loosely aligned with bacterial shapes (Fig. 1i, l-m; Supplementary Fig. 2h-i, k-l) or released as empty membrane spheres (Fig. 1j, m, n). While loosely attached symbiosome membranes were observed in all *symrem1* mutant nodules, only some nodule cells contained empty membrane spheres. These structures rather resembled symbiosome membranes than tonoplasts as they contained PS (Fig. 1j), a lipid that does not occur on vacuolar membranes⁴⁰.

So far, SYMREM1 has been described as a molecular scaffold that interacts with and stabilizes the entry receptor LYK3 at the plasma membrane²⁶,³⁵. While LYK3 has been localized to ITs²⁶,³³, no persistent accumulation on the infection droplet nor the symbiosome membrane has been observed⁴¹. This is in contrast to SYMREM1 that remained and even accumulated on these unwalled symbiotic membranes (Fig. 1d-f; Supplementary Fig. 1g-i) implying that SYMREM1 serves additional functions. Given the lack of stabilized tubular outgrowths from ITs close to infection droplets, our findings that symbiosome membranes appeared less supported and consequently loosely organized around the bacteroid and the occurrence of empty membrane spheres in *symrem1* mutants, we hypothesized that SYMREM1 functions in maintaining distinct topologies of unwalled membranes (summarized in Supplementary Fig. 3).

### SYMREM1 functions as a structural membrane scaffold

The hypothesis of SYMREM1 being a structural membrane scaffold was further supported by the fact that other members of the remorin family have been shown to form higher order proteinaceous lattices and/or filamentous structures³⁰,³¹,⁴²,⁴³ in vitro. We also confirmed such filamentous structures for purified, recombinant SYMREM1 using TEM analysis followed by 2D classification of manually picked filamentous segments. Here, we detected auto-assembled amorphous protein filaments that were partially branched or scrambled (Fig. 2a). Systematic inspection of 389 filament fragments revealed an average width between 84 and 125 Å for the filamentous particles with a few class averages showing helical features such as twists (Fig. 2b). Besides this, we mostly observed irregular filaments and protein bodies (Fig. 2a) that may represent filament seeds. Analytical gel-filtration of
<figure><img src="image_5.png" /></figure>

**Fig. 1 | SYMREM1 stabilizes confined membrane tubes during infection.**
**a–c** Illustrations indicate an infection thread (IT) with an infection droplet (ID; **a**), symbiosome (S) formation (**b**) and a symbiosome-filled infected cell (**c**); CW cell wall, V vacuole. **d–f** Expression of a mEosEM-SYMREM1 fusion protein (green) specifically labels IT membranes (**d**), accumulates at bacterial droplet structures (ID) (**e**), and symbiosome membranes (**f**). *S. meliloti* expressing a mCherry marker is shown in magenta. The white arrow (in **e**) indicates the bacterial release site; the black circle line (in **f**) indicates the central vacuole. Single focal plane images for maximum projections of z-stacks (**d–f**) can be found in Supplementary Fig. 1. Experiments were performed with three biological replicates with at least 10 nodules being imaged per replicate. **g–j** Membranes were visualized by expressing

the phosphatidylserine (PS) biosensor LactC2. Spatially confined membrane tubes (arrow) were found on wild-type (WT) IT containing rhizobia (magenta, **g**) but not on ITs in *symrem1* mutants (**h**). Symbiosome membranes are loosely associated with released rhizobia (arrows, **i**) or appear as empty spheres (arrowheads, **j**) in *symrem1* mutants. Data were collected based on three biological replicates. Images (**d–j**) were taken as z-stacks (internal distance is 0.5 µm) and shown as 3D projections generated by using Imaris. These patterns were confirmed by transmission electron microscopy for WT (**k**) and *symrem1* mutants (**l–n**). Scale bars indicate 5 µm (**d–j**) and 2 µm (**k–n**). IT infection thread, ID infection droplet. The sketches (**a–c**) were drawn with Inkscape.

fresh SYMREM1 protein extracts revealed an apparent molecular weight of 57 kDa that matches the calculated size for a dimeric remorin protein (Fig. 2c). Since the filaments were too amorphous for further structural assessment by cryo-EM, we conducted 3D modeling using AlphaFold⁴⁴.

In line with previous analyses, AlphaFold predicted a largely disordered N-terminal region encompassing residues 1–69, followed by an α-helical segment from residues 70–187 and a disordered C-terminus at residues 188–205. Interestingly, the search for multimers of SYMREM1 consistently yielded a head-to-tail dimer as a recurrent structural unit (Fig. 2d). Here, the slightly bent α-helical segments of two monomers interacted through extended hydrophobic patches, with the N-termini pointing outwards. Using AlphaFold²⁴⁵ predictions for higher-order oligomers we repeatedly obtained alignments in flexible sheets that can be extended to helical structures (Fig. 2e–g).

Our model further indicated a consistent interaction between two dimers that was not mediated by hydrophobic interactions, but rather by a few selected hydrogen bonds. Calculating electrostatic surface potentials on these structures revealed a strong enrichment of positive potentials at one site of the sheet (Fig. 2f). Interestingly, the propagation of these dimeric interactions resulted in the prediction of helical superstructures with flexible inner diameters (Fig. 2g). The disordered N-termini are extending to each side of the helix and are not visualized in the model.

Next, we asked whether the predicted sheet-like and potential super-helical organization of the SYMREM1 oligomer has an impact on membrane topologies along infection threads. For this, we ectopically expressed fluorescently labeled SYMREM1 in WT *M. truncatula* plants, a setup that was previously shown to elevate nodulation levels in *Lotus japonicus*⁴⁶. Here, we observed that the frequency of stabilized
<figure><img src="image_6.png" /><figcaption>**Fig. 2** | SYMREM1 forms oligomeric alpha helical assemblies. **a** Representative raw electron micrograph of purified, recombinant SYMREM1 stained with 2% uranyl acetate. Arrowheads indicate irregular protein bodies. Scale bar indicates 100 nm, experiments were performed twice with independently isolated recombinant SYMREM1 protein. **b** 2D class averages derived from multivariate statistical analysis of all 389 particle images. Each class contains on average ten images. Class averages that show twisted features are marked with a white asterisk. Scale bar indicates 100 Å. **c** Elution profile of recombinant His-SYMREM1. Molecular masses (in kDa) and positions of elution peaks for standard proteins are indicated with triangles on the top. Insert: SDS-PAGE after Coomassie staining of the purified His-SYMREM1 (labeled by an arrow); molecular masses of the pre-stained protein standards are indicated on the left in kDa. **d** AlphaFold predictions for homodimers of SYMREM1</figcaption></figure>

(left) and MtREM2.1 (right). One monomer is colored from blue at the N-terminus to red at the C-terminus, the other in white. The Extended helical regions of both remorins form highly similar, antiparallel dimers. **e** Prediction of higher-order oligomers using AlphaFold2. The remorin homodimers form flexible sheets that can be extended into helical structures. **f** Electrostatic surface potential maps for the two faces of the sheets formed by SYMREM1 (above) and MtREM2.1 (below), contoured from -5k_BT (red) to +5k_BT (blue). In both cases, the convex faces show a positive electrostatic potential, while that of the concave faces is negative. **g** A predicated helical super-structure based on the oligomerization of SYMREM1 dimers. It is important to note that the diameter of this higher order structural prediction is highly variable and depends on the alignment of the individual dimers.

membrane tubes remained unaltered (Fig. 3a) compared to those observed in LactC2-labeled WT cells (Fig. 1g), indicating a temporary nature of these structures at membrane interfaces that continuously release rhizobia. To test this hypothesis, we made use of the release-compromised *ipd3-1* mutant, which is defective in the transcriptional activator CYCLOPS that regulates, among other genes, the expression of endogenous *SYMREM1*⁴⁷. This reported transcriptional regulation also translated into reduced SYMREM1 protein levels in the *ipd3-1* mutant (Supplementary Fig. 4). In line with this and similar to WT plants, PS-labeling revealed membrane tube formation in 36% of all assessed ITs in the *ipd3-1* mutant (Fig. 3b). However, ectopic expression of fluorophore-tagged SYMREM1 in *ipd3-1* significantly increased IT-associated membrane tubulation to 75% with several tubes per IT being frequently observed (Fig. 3c). These data further support a function of SYMREM1 in tubulation of membranes that are not supported by a rigid cell wall.

To address the types of negatively charged lipids that can be bound by the positively charged surfaces of SYMREM1 oligomers, we purified recombinantly expressed and His-tagged SYMREM1 from *E. coli* and hybridized it to a lipid strip containing 15 different spotted lipids. Clear binding of SYMREM1 was observed for phosphatidylinositol (PI), PI(3)P, PI(3,4)P₂, PI(3,5)P₂, and PI(4,5)P₂ while PI monophosphates such as PI(4)P and PI(5)P were not bound by SYMREM1 (Supplementary Fig. 5a). To verify these results independently, we generated giant unilamellar vesicles (GUVs) comprised of PI(4)P, PI(3,5)P₂ or PI(4,5)P₂ and incubated these GUVs with recombinant SYMREM1. While no association was observed with GUVs containing PI(4)P (Supplementary Fig. 5b), SYMREM1 bound to GUVs with PI(3,5)P₂ and PI(4,5)P₂, which resulted in lipid accumulations at these sites

(Supplementary Fig. 5c, d). Although recombinant SYMREM1 associated with these GUVs it should be noted that the protein is not palmitoylated (S-acylated) when being isolated from *E. coli*. This lipidation has, however, been reported to greatly enhance tight membrane binding of SYMREM1³⁶. To test the effect of tight membrane association of SYMREM1, we generated GUVs comprised of Ni-NTA lipids (18:1 DGS-NTA(Ni)) that should tightly associate with His-tagged SYMREM1. While His-tagged GFP (control) did not alter the GUVs at all (Fig. 3d), SYMREM1 application resulted in dramatically altered topologies with many membrane blebs being formed (Fig. 3e) that frequently detached as membrane spheres (Supplementary Fig. 6a). We also occasionally observed stabilized membrane invaginations (Supplementary Fig. 6b). We then expanded this analysis and generated recombinant His-tagged SYMREM1 variants comprised of the sole N-terminal intrinsically disordered region (SYMREM1¹D*, residues 1-73) and the C-terminal region (SYMREM1Cterm; residues 74-205). While the association of SYMREM1¹D* did not change GUV topologies (Fig. 3f), incubation of GUVs with SYMREM1Cterm had similar effects as observed for full-length SYMREM1 with membrane blebs (Fig. 3g; Supplementary Fig. 6c) and membrane invaginations (Supplementary Fig. 6d) being formed.

To assess whether we can recapitulate these findings in an in vivo system, we expressed these SYMREM1 variants as YFP fusion proteins in cell wall-free *Nicotiana benthamiana* mesophyll protoplasts. Such cells adopt a spherical shape under mild hypotonic conditions as visualized by expressing a generic membrane marker probe (mCitrine-LTI6b; Supplementary Fig. 7a). However, when expressing full-length SYMREM1, we observed numerous tubular outgrowths developing on 64 out of 112 (57%) inspected protoplasts (Fig. 4a-c, Supplementary Fig. 7b). The formed tubular structures greatly varied in length with the
<figure><img src="image_7.png" /><figcaption>**Fig. 3** | SYMREM1 can induce membrane morphology changes. a–c Membrane tubes were found on nodular ITs in the *ipd3-1* mutant (arrow, b), while ectopic expression of SYMREM1 greatly increased these tubular outgrowths in the *ipd3-1* mutant (arrows, c) but not on WT (arrow, a). Membranes were visualized by expressing mEosEM-SYMREM1 (a, c) or the phosphatidylserine (PS) biosensor LactC2 (b), both shown in green. IT infection thread. *S. meliloti* expressing a mCherry marker is shown in magenta. Data were collected based on three biological replicates. All images were taken as z-stacks (internal distance is 0.5 µm) and</figcaption></figure>

are shown as 3D projection generated by using Imaris. **d–g** Application of His-tagged GFP (control, green) did not alter the morphology of Ni-NTA containing GUVs (Atto 647N-DOPE was used as membrane marker; magenta) (**d–d''**). By contrast, application of His-GFP tagged SYMREM1 (green, **e–e''**) and His-GFP tagged SYMREM1<sup>Cterm</sup> (green, **g–g''**) resulted in dramatically altered topologies with many membrane blebs being formed. This was not observed when applying the His-GFP-tagged SYMREM1<sup>IDR</sup> variant (green, **f–f''**). GUV experiments were performed twice and showed similar results. Scale bars indicate 5 µm (**a–c**) and 10 µm (**d–g**).
<figure><img src="image_8.png" /><figcaption>**Fig. 4** | SYMREM1 stabilizes membrane tubulation and curvature in a cell wall-independent manner. a–c *N. benthamiana* leaf protoplasts ectopically expressing YFP-SYMREM1 develop multiple membrane tubes as shown by confocal laser-scanning microscopy (a) and scanning electron microscopy (b) and close-up in (c). d–i Protoplasts expressing different SYMREM1 variants: only IDR domain (SYMREM1<sup>IDR</sup>, d), only C-terminal including RemCA (SYMREM1<sup>Cterm</sup>, e), full-length sequence lacking the RemCA domain (SYMREM1<sup>FLΔRemCA</sup>, f), RemCA domain only</figcaption></figure>

(SYMREM1<sup>RemCA</sup>, g, h), and the RemCA domain bearing a mutation of the palmi-
toylated Cys197 residue (SYMREM1<sup>RemCA C197A</sup>, i). Scale bars indicate 10 µm. All the
confocal images are shown as maximal projections. Magenta signals derive from
chlorophyll autofluorescence within chloroplasts. Numbers indicate frequencies of
all observations. Protoplast isolation experiments were at least performed in four
biological replicates, with similar tube formation frequencies found in each
replicate.

longest ones being more than 70 mm and occasionally branched
(Supplementary Fig. 7c). Scanning electron microscopy (SEM) (Fig. 4b,
c) revealed an average diameter of these protrusions of
0.18 ± 0.03 mm (Supplementary Fig. 7d). Expression of SYMREM1<sup>IDR</sup> in
protoplasts failed to induce and stabilize these tubes as verified by
confocal laser-scanning microscopy (Fig. 4d) and SEM (Supplementary
Fig. 7e, f). By contrast, expression of the helical C-terminal region
(SYMREM1<sup>Cterm</sup>) phenocopied patterns observed for the full-length
protein (Fig. 4e) even though the tubes were mostly shorter. Taken
together, these observations are consistent with those from GUVs.

Since SYMREM1 associates with the plasma membrane via the
amphipathic and often palmitoylated RemCA peptide (amino acids
171-205)³⁶, we tested the impact of these C-terminal 35 amino acids on
membrane topologies. Truncating RemCA from full-length SYMREM1
(SYMREM1<sup>FLΔRemCA</sup>, residues 1-170) resulted in a cytosolic protein that
was unable to induce and/or stabilize these membrane outgrowths
(Fig. 4f). Expression of the RemCA membrane anchor alone (SYM-
REM1<sup>RemCA</sup>) was sufficient to drive membrane topology changes. While
13 out of 121 protoplasts developed tubes on the protoplast surface
(Fig. 4g), 23 out of 121 protoplasts exhibited less confined and more
bulky membrane blebs (Fig. 4h). Mutating the palmitoylated Cys197
residue in RemCA (SYMREM1<sup>RemCA C197A</sup>) fully abolished membrane
association of the peptide and its impact on membrane topology

(Fig. 4i). These data show that the RemCA peptide alone is sufficient for
initiating membrane tubulation but cannot maintain and/or stabilize
long protrusions at high-frequency.

SYMREM1 indirectly associates with actin

As tubular membrane outgrowths such as pollen tubes and root hairs
are usually actin-dependent, we co-expressed full-length SYMREM1
with the actin marker LifeAct in protoplasts and observed that all
tubes contained central actin filaments (Fig. 5a). The importance of
this was additionally supported by the fact that application of cyto-
chalasin D, an actin depolymerizing agent, dramatically decreased the
number of tubular outgrowths (11 out of 118 protoplasts) in proto-
plasts expressing SYMREM1 (Fig. 5b). Instead, upon actin depoly-
merization we observed the same kind of membrane blebs (Fig. 5c) as
induced upon expression of the SYMREM1<sup>RemCA</sup> peptide (Fig. 4h). By
contrast, the application of the microtubule de-polymerizing agent
oryzalin did not alter membrane tube formation (Fig. 5d). Further-
more, expression of the symbiotic and actin-associated formin pro-
tein SYFO1<sup>14</sup> resulted in a tip-localized signal at these protrusions
(Fig. 5e) with SYFO1 being recruited to SYMREM1 foci over time
(Supplementary Fig. 8a). Interestingly, SYFO1 was also present in
membrane blebs induced by either expressing SYMREM1 and sub-
sequent cytochalasin D treatment (Fig. 5f) or when solely expressing
<figure><img src="image_9.png" /><figcaption>**Fig. 5** | The formation of membrane tubes is actin-dependent. a–a'' *N. benthamiana* leaf protoplasts ectopically expressing mCherry-SYMREM1 (magenta) developed multiple membrane tubes that comprised a central actin strand as labeled by LifeAct (green). b–d Protoplasts ectopically expressing YFP-SYMREM1 and treated with Cytochalasin D (b, c) and Oryzalin treatment (d). Magenta signals in (b–d) derived from chlorophyll autofluorescence within chloroplasts. Numbers indicate frequencies of observations and images are shown as maximal projections.</figcaption></figure>

e–f' *N. benthamiana* protoplasts ectopically expressing mCherry-SYMREM1 (magenta) developed multiple membrane tubes with a tip localized formin protein SYFO1 (green, white arrowheads, e–e'). SYFO1 remained to be recruited into membrane blebs after Cytochalasin D treatment of these protoplasts (green, white arrowheads, f–f'). Scale bars indicate 10 µm. Data were obtained from four independent biological replicates.

the RemCA peptide (Supplementary Fig. 8b) while no actin strands were found in these blebs (Supplementary Fig. 8c). Since cytochalasin D treatment of SYMREM1-expressing protoplasts severely inhibited membrane tubulation we concluded that SYMREM1 stabilizes rather than actively drives outgrowth of these membrane tubes in an actin-dependent manner.

To investigate whether SYMREM1 and actin directly interact, we performed co-sedimentation assays. To demonstrate functionality of the setup, we first incubated actin with the known actin-binding formin protein mDial-ct from human cells⁴⁸. In the control samples G-actin and mDial-ct were predominantly present in the pellet and the supernatant, respectively, while incubation of both proteins resulted in an enrichment of mDial-ct in the pellet, indicative of actin-binding (Supplementary Fig. 9a). This phenomenon was neither observed when mixing 10 or 30 µM SYMREM1 with actin since SYMREM1 remained in the soluble fraction (Supplementary Fig. 9b, c). These data

suggest that actin and SYMREM1, at least in vitro, do not interact directly with each other.

To test the association of SYMREM1 with actin under symbiotic conditions, we conducted co-labeling experiments on nodule sections, where we labeled actin with either Phalloidin or by using stable transgenic plants expressing an actin reporter (ABD-mCherry). SYMREM1 was labeled by a specific antibody³⁵ or by expressing a fluorescently tagged SYMREM1 protein. In these experiments, SYMREM1 was found along infection threads within the nodular fixation zone (Fig. 6a, Supplementary Fig. 10) with actin being ubiquitously present in all cells. To gain a more detailed view on this, we applied structured illumination super-resolution microscopy (SIM) to these samples. Actin and SYMREM1 tightly aligned but were slightly displaced to each other along nodular ITs (Fig. 6b–d). A similar situation was found in the nodular fixation zone, where both proteins showed the same association pattern (Fig. 6e–h; Supplementary Figs. 10, 11).
<figure><img src="image_10.png" /><figcaption>**Fig. 6 | SYMREM1 and F-actin colocalize at infection structures inside nodules.**
Sections of two weeks old nodules from a stable transgenic *Medicago truncatula* plant expressing an ABD-mCherry actin reporter (green). *a*, *e* SYMREM1 was immunolabelled (magenta) using a specific α-SYMREM1 antibody. Maximum intensity projection of an infection thread (IT) (*b*) and symbiosomes (*f*) immunolabelled for SYMREM1 (magenta), counterstained with Phalloidin (green) and imaged using super-resolution structured illumination microscopy (SIM). c–c′, g–g′′ Close-ups of areas indicated by white boxes in (*b*, *f*). *d*, *h* 2D linescans of the</figcaption></figure>

transect indicated by the white arrow in (c", g"). Images are maximum projections (*a*, *b*, *e*) or single plane SIM images (*f*, c–c", g–g"). Structures depicted in (f–g") derived from symbiosome membranes for which the raw images are provided in Supplementary Fig. 11. Scale bars indicate 5 µm (*a*, *b*, *e*, *f*) and 200 nm in (*c*, *g*). For confocal analysis (*a*, *e*), data were collected based on three biological replicates, with at least eight nodules being assessed for each replicate. For SIM analysis, two biological replicates were performed and three nodules were assessed for each replica.

## Membrane scaffolding is conserved among remorin proteins

To further test whether SYMREM1 is able to stabilize exogenously induced membrane deformations and functions as a stabilizing scaffold exclusively for positive curvatures, we assessed whether negative curvatures, as underlying IT-associated membrane tubes, were also maintained in the presence of this protein. Therefore, we isolated SYMREM1-expressing protoplasts and immediately indented them with a micro-capillary for 30 min. Here, 10/13 protoplasts expressing the LTI6b membrane marker (control) re-inflated directly after releasing the pressure (Supplementary Fig. 12a). This indicates that the cytoskeleton alone is not sufficient to dynamically maintain these short-term indentations as shown for protoplasts that had been confined in controlled geometries⁴⁰. By contrast, full re-inflation was observed in only 4/14 protoplasts expressing SYMREM1, while microcapillary-induced membrane deformations were maintained in the majority (10/14) of these protoplasts (Supplementary Fig. 12b). Such stabilization was also observed in most protoplasts expressing the truncated SYMREM1Cterm variant (Supplementary Fig. 12c) while the RemCA peptide alone (SYMREM1RemCA) was not sufficient to mediate this phenomenon (Supplementary Fig. 12d). These data suggest that topology scaffolding of SYMREM1 is highly flexible and does not rely on an invariable protein structure.

As a last set of experiments, we assessed whether membrane topology scaffolding is limited to SYMREM1 or a more wide-spread feature within the remorin protein family. In analogy to the *Medicago* SYMREM1 protein, expression of the orthologous gene from *Lotus*

*japonicus* (*LjSYMREM1*; Lj4g3v2928720) resulted in stabilized membrane tubes at high frequency (Supplementary Fig. 13a). Evolutionary, the presence of SYMREM1 proteins is mainly restricted to species that maintained the root nodule symbiosis and species that do not support intracellular infection of symbionts, either fungal or bacterial, through their epidermis have lost *SYMREM 1* (Supplementary Fig. 14). In addition, the only non-symbiotic clade that retained *SYMREM1*, the Caryophyllales, experienced bursts of diversifying positive selection indicative of neofunctionalization (Supplementary Data 1), suggesting that the purifying selection maintaining *SYMREM1* in plant genomes is linked to the maintenance of symbiotic abilities. Interestingly, a closely related protein (MedtruChr5g0397781; MtREM2.1) shows an expanded phylogenetic distribution, while both *SYMREM1* and *MtREM2.1* derive from a Papilionoideae duplication. To assess whether MtREM2.1 has similar impact on membrane topology as shown for SYMREM1, we expressed a MtREM2.1 fusion protein in mesophyll protoplasts. Indeed, the presence of this protein resulted in the stabilization of numerous membrane tubes (Supplementary Fig. 13b), and this function is supported by AlphaFold predictions that revealed a putative structure that is highly similar to the one predicted for SYMREM1 (Fig. 2d–f). To address whether MtREM2.1 localizes to highly curved membranes, we first surveyed its gene expression profile by using global transcriptome data, which revealed an induced expression during endosymbiotic interactions⁵⁰. We independently confirmed this by quantitative Real-Time PCR on cDNA generated from nodulated and mycorrhized roots (Supplementary Fig. 15a). To increase the
<figure><img src="image_11.png" /><figcaption>**Fig. 7** | Remorin-induced membrane tubulation in protoplasts. a–h *N. ben-thamiana* leaf protoplasts ectopically expressing different remorins from the non-symbiotic plant *Arabidopsis thaliana* covering all sub-groups of the Arabidopsis remorin family. AtREMI.3 (a, b), AtREMI.4 (c), AtREM2.1 (d), AtREM3.2 (e), AtREM4.2 (f), AtREM5.1 (g), AtREM6.1 (h). i–k Membrane tubulation activity was also assessed for a N-terminally truncated variant of AtREMI.3 (AtREMI.3<sup>Cterm</sup>, i), a chimeric protein comprised of the AtREMI.3 N-terminal and the SYMREM1</figcaption></figure>

C-terminal region (AtREMI.3<sup>Nterm</sup>::SYMREM1<sup>Cterm</sup>, j), a dual-legume SYMREM1 chimera (LjSYMREM1<sup>Nterm</sup>::SYMREM1<sup>Cterm</sup>, k) and SYMREM1 (l). Scale bars indicate 10 µm. All images are shown as maximal projections. The magenta signal indicates chlorophyll autofluorescence from chloroplasts. Numbers indicate the frequencies of the presented observations across at least three independent biological replicates.

spatio-temporal resolution of this analysis, we generated a transcriptional reporter with a 1 kb long fragment of the putative *MtREM2.1* promoter driving the expression of a β-glucuronidase gene (GUS). *MtREM2.1* promoter activity was detected in young and maintained in the apical region of mature nodules (Supplementary Fig. 15b–d). In addition, we found increased GUS activity in cells containing arbuscules when inoculating roots with the arbuscular mycorrhizal fungus *Rhizophagus irregularis* (Supplementary Fig. 15e, f). To assess the localization of the MtREM2.1 protein, we generated *Medicago truncatula* stable lines (WT R108 background) expressing an mCherry-MtREM2.1 fusion protein under the control of the endogenous promoter. In roots inoculated with *R. irregularis*, MtREM2.1 predominantly localized to the host-derived and highly curved membrane (periarbuscular membrane) surrounding intracellular hyphae, from their early penetration into inner cortical cells (Supplementary Fig. 15g) until the development of fully branched arbuscules (Supplementary Fig. 15h). Interestingly this arbuscular membrane is, similarly to infection droplets or symbiosome membranes, not supported by a rigid cell wall⁵.

While all three tested group 2 remorins (Medicago SYMREM1, LjSYMREM1, and MtREM2.1) stabilized membrane tubes on protoplasts, systematic expression of different members of the

Arabidopsis remorin family, confirmed this ability for AtREMI.3, AtREMI.4, and AtREM4.2 even though membrane tubulation was stabilized with much lower frequencies in all cases (Fig. 7a–c, f). In addition, and similar to the effect observed in the presence of RemCA (Fig. 4h) or cytochalasin D (Fig. 5c), expression of AtREMI.3 resulted in membrane bleb formation in 11% of the cases. No such effects were observed for a so far uncharacterized remorin (AtREM2.1; AT1G63295), which was mostly cytosolic (Fig. 7d), whereas expression of AtREM3.2, AtREM5.1, and AtREM6.1 led to focal accumulations and the formation of shorter membrane tubes (Fig. 7e, g, h). These results indicate that the ability of remorins to act as membrane topology scaffolds is conserved in land plants. Besides this, it was intriguing to note that expression of an N-terminally truncated variant of AtREMI.3 lacking its IDR significantly increased the ability to stabilize membrane tubes from 6% to 24% (Fig. 7i), suggesting an inhibitory role of the disordered N-terminal region. This hypothesis was supported by expressing a chimeric protein comprised of the AtREMI.3 N-terminal and the Medicago SYMREM1 C-terminal region (AtREMI.3<sup>N-term</sup>::SYMREM1<sup>C-term</sup>) that reduced membrane tubulation frequencies from 57% (Figs. 4a, 7l) to 17% (Fig. 7j), while expression of a LjSYMREM1/SYMREM1 chimera maintained high tubulation levels at 47% (Fig. 7k).
Taken together, these results show that the largely legume-
confined remorins SYMREM1 and MtREM2.1 represent evolutionary
optimizations to structurally support curvatures of cell wall-devoid
membranes throughout endosymbiotic infection processes.

Discussion

Plant membranes are usually supported by a rigid cell wall that coun-
teracts the intracellular turgor pressure and maintains, together with
cortical microtubules and actin, cellular shapes over their
lifespan⁴⁹,⁵²⁻⁵⁴. As all cell wall constituents are either secreted to the
apoplast or synthesized directly onto the extracellular leaflet of the
plasma membrane, these responses require time. Furthermore, sub-
sequent loosening of the rigidified cell wall to adopt temporal changes
in membrane topologies mostly account for protrusion-like shapes. In
contrast, larger scale and negatively curved membrane invaginations
as found during microbial infections of host cells may require other
types of scaffolding prior to de novo cell wall apposition. Such func-
tions can be maintained by scaffold proteins as exemplified during
clathrin-mediated endocytosis, where the adaptor complex AP2
together with clathrin light and heavy chains stabilizes endocytotic
vesicles²²,²³. In addition, multimeric and helical forming ESCRT-III
subunits predominantly contribute to membrane remodeling and
cytokinetic abscission⁵⁵,⁵⁶. In plants, further evidence has been pre-
sented that endocytosis of membrane nanodomain-localized proteins
might occur via a flotillin- and remorin-dependent but clathrin-
independent pathway⁵⁷⁻⁵⁹. However, compelling evidence for topolo-
gical scaffolding functions of these proteins has been missing, while
the impact of remorins on membrane fluidity is experimentally
supported²⁸.

Here, we demonstrate that plant remorin proteins can stabilize
membrane topologies independent of a cell wall with legume group 2
remorins, which are mostly found in plants that have maintained a
symbiotic association with (both) rhizobia and arbuscular mycorrhiza
fungi, showing the strongest effect on membrane tubulation (Fig. 4;
Supplementary Fig. 13). This makes sense since both endosymbioses
require large-scale membrane invaginations during intracellular colo-
nization and stabilization of the symbiotic membranes. Here, SYM-
REM1 is required for stabilizing the IT and the symbiosome membrane
(Fig. 1d–f; Supplementary Fig. 1g–i). By contrast, MtREM2.1 is poten-
tially involved in stabilizing the periарbuscular membrane enclosing
mycorrhizal hyphae (Supplementary Fig. 15g, h). Although cell wall
components are deposited within the older parts of ITS²⁰, the initial
curvature needs to be differently stabilized and additionally main-
tained around cell wall-free bacterial release sites. In contrast to
membrane bending proteins belonging to the Bin/Amphiphysin/RVS
(BAR) domain family that mediate either positive of negative
curvature⁶⁰ and that are, amongst others, involved in autophagosome
formation⁶¹ and cell plate formation⁶² in plants, SYMREM1 seems
capable of stabilizing both types of directions (Fig. 4; Supplementary
Figs. 6; 12). This feature is mainly mediated by the evolutionary con-
served C-terminal region of the protein, which is additionally required
for protein oligomerization into dimers via hydrogen bonding and
hydrophobic interactions (Fig. 2d). Membrane association is further
assisted by a palmitoylation of a C-terminal cysteine residue³⁶ (Fig. 4i).
We propose that the required forces are provided by a high membrane
binding energy and internal rigidity of the SYMREM1 scaffold that
exceed the energy required for membrane bending and membrane
rigidity, respectively, as suggested for other scaffold proteins⁶³. This
assumption is supported by our AlphaFold model predicting curved
dimers and multimers (Fig. 2d–g) and by the fact that accumulated
shallow insertions of membrane scaffolds, as also mediated by the
RemCA anchor²⁷ and palmitoylation for SYMREM1, have been gen-
erally shown to be effective to induce membrane curvatures⁶⁴. Fur-
thermore, the formation of higher order protein sheets or filaments
has been shown to occur in vitro³¹ (Fig. 2a, b). The impact of such

protein lattices on membrane curvature has also been demonstrated for the human N-Bar protein Endophilin that can adopt higher-order arrangements⁶⁵. Similar observations have also been reported for the F-Bar protein Imp2 from fission yeast, where a helical alignment of Imp2-subunits results in membrane tubulation of human Cos-7 cells⁶⁶. Whether SYMREM1 can form such helical superstructures as predicted (Fig. 2g) and shown for Imp2 and ESCRT-III proteins⁵⁵,⁵⁶ remains to be investigated. Potential hints have recently been provided for the group 1 remorin AtREM1.3 when being expressed in human Cos-7 cells where AtREM1.3 was found to helically wrap around actin strands⁶⁷. Evidence for such tight association of actin and SYMREM1 is also provided in this study, where we showed that membrane tube elongation in protoplasts is dependent on an intact actin cytoskeleton (Fig. 5) and likely driven by formin-mediated actin polymerization as also shown for filopodia growth and retraction⁶⁸. During innate immune responses in Arabidopsis the control of formin condensation is also mediated by remorins³⁸. Furthermore, SYMREM1 and actin tightly align alongside symbiotic membranes (Fig. 6; Supplementary Fig. 10). However, a direct interaction between SYMREM1 and filamentous actin could not be demonstrated in vitro (Supplementary Fig. 9), whereas such interactions have been shown for the group 6 remorin GSD1 from rice⁶⁹.

In summary, we hypothesize that SYMREM1 provides a scaffolding
structure that enable the stabilization of large-scale membrane
topologies as found during bacterial release at nodular infection
threads (Fig. 1g). We favor a model in which shallow insertions of oli-
gomeric and rigid remorin scaffolds provide the forces for membrane
bending and stabilization that may, however, be assisted but other, yet
unknown, proteins. In our hands, a membrane-bending model exclu-
sively based on protein crowding is less likely. While it was reported
that sole crowding of any membrane-associated protein when
exceeding a membrane coverage of more than 20% induces membrane
tubulation in GUVs²⁰, this was not confirmed in our study as neither
membrane-immobilized GFP nor SYMREM1<sup>IDR</sup> (Fig. 3d, f) or insertion of
the transmembrane protein LTI6b into protoplasts (Supplementary
Fig. 7a) induced these effects. In vivo, such tubulations and invagina-
tions, are sufficient to drive cytoskeleton polarization and conse-
quently targeted secretion of proteins to these sites. In *symrem1*
mutant nodules, the lack of such polarization may therefore prevent
the secretion of proteins required for bacterial release.

Taken together we unraveled a novel mechanism that allows plant cells to stabilize distinct membrane topologies in the absence of a cell wall by an interplay between oligomeric remorin scaffolds and the cytoskeleton. Consequently, remorins may serve roles additional to those of BAR-domain, ESCRT or other scaffold proteins<sup>71</sup>.

Methods

Plant growth and Rhizobia inoculation

Seeds of *M. truncatula* were surface sterilized by covering them with pure sulfuric acid (H₂SO₄) for 10 min, followed by 4–6 times washing with sterile water. Seeds were then covered with bleaching solution (12% NaOCl, 0.1% SDS) for no longer than 60 seconds and washed 4 to 6 times with sterile water. After surface sterilization, seeds were transferred to 1% agar plates and stratified at 4 °C for 3 days in darkness. Germination was allowed for up to 24 h at 24 °C in darkness. The seed coat was removed and seedlings were transferred to plates containing Fahräeus medium supplemented with 0.5 mM NH₄NO₃. One week later, they were transferred onto fresh Fahräeus medium without nitrate, but containing AVG (0.1 µM). Inoculations were performed after 4 days of growth on plates without nitrate. For inoculation of *M. truncatula* roots, a *S. meliloti* (Sm2011) liquid culture was centrifuged (3 min, 3000 rpm), washed once with liquid Fahräeus medium and resuspended in liquid Fahräeus medium to a final OD₆₀₀ = 0.03. Each root was covered with 1 ml of rhizobia suspension, which was removed after 6 minutes. Afterwards, the plants were placed in a controlled environment chamber at 24 °C with a 16/8 h light/dark photoperiod,
keeping the roots in the dark, for 3 weeks before harvesting the nodules.

## Hairy root transformation and rhizobial inoculation

M. truncatula hairy root transformation was performed as previously described<sup>72</sup>. Briefly, transgenic *Agrobacterium rhizogenes* (ARqua I), carrying the plasmid of interest, was grown in LB liquid culture for one day and 300 µl of the culture were spread on LB agar plates supplemented with the corresponding antibiotics for selection, and incubated on plates for two more days before transformation. M. truncatula seeds were prepared as mentioned above. After germination the seed coat was removed from the cotyledons of the seedlings under water, and the root meristem was cut off with a scalpel. Cut seedlings were dipped on the Agrobacterium plates and transferred onto solid Fahraeus medium (containing 0.5 mM NH<sub>4</sub>NO<sub>3</sub>). Transformed seedlings were incubated for three days at 22 °C in darkness, following 4 days at 22 °C in white light but keeping the roots in the dark. One week after transformation, seedlings were transferred onto new Fahraeus medium (0.5 mM NH<sub>4</sub>NO<sub>3</sub>) and grown for another 10 days at 24 °C in a controlled environment chamber with 16 h/8 h light/dark photoperiod. Afterwards, the roots were screened to examine the transformation efficiency using a stereomicroscope to detect the corresponding fluorescent signal. Untransformed roots were cut off and plants showing fluorescence roots were transferred to pots (2 plants per pot) having a mixture of equal volume of quartz-sand and vermiculite. All pots were individually watered with liquid Fahraeus medium (without nitrate) and tap water once a week. After 3–5 days, the pots were inoculated with S. meliloti (OD<sub>600</sub> = 0.003).

## Evolutionary analysis

Homologs of MtSYMREM1 (Medtr8g097320.2 or MtrunA17Chr8g0386521) were retrieved using the tBLASTn v2.11.0+<sup>73</sup> against a database of 189 species covering all Viridiplantae lineages with default parameters and e-value threshold of 1e-10. Coding sequences of putative homologs were aligned using MAFFT v7.471<sup>74</sup> with default parameters. Positions with more than 80% of gaps were removed from the subsequent alignment using trimAl v1.4rev15<sup>75</sup> and cleaned alignment subjected to phylogenetic analysis using Maximum Likelihood approach. Prior tree reconstruction, best-fitting evolutionary model was tested using ModelFinder<sup>76</sup> according to the Bayesian Information Criteria. Maximum likelihood analysis was conducted using IQ-TREE2 v2.0.3<sup>77</sup>. Branches support were estimated using 10,000 replicates of both sh-aLRT<sup>78</sup> and UltraFast Bootstraps<sup>79</sup>. Analysis of the main tree revealed that MtSYMREM1 clade derives from Eudicots duplication and then, a subtree corresponding to the Eudicots clade of MtSYMREM1 has been reconstructed based on protein sequences aligned using MUSCLE v3.8.1551<sup>80</sup> and phylogenetic procedure described above. Trees were visualized using the iTOL v6 platform<sup>81</sup>.

To look for specific selective pressure acting on the Caryophyllales clade, which is composed only of species not forming infection threads, we conducted branch and branch site analysis. We estimated relaxation (*K* < 1) or intensification (*K* > 1) parameters and also positive selection acting on Caryphyllales (Supplementary Data 1), we used the RELAX and aBSREL programs implemented in the HYHPY software<sup>82,83</sup>. These methods calculate different synonymous and nonsynonymous substitution rates (ω = dN/dS) using the phylogenetic tree topology for both foreground and background branches. Translated CDS of MtSYMREM1 homologs were aligned using DECIPHER<sup>84</sup>. In total, 191 sequences and 162 codons were analyzed (Supplementary Data 2).

## Nodule sections

For analysing the subcellular localization of SYMREM1 and the LactC2 biosensor in WT and *ipd3* (*sym1-TE7*) mutant plants, the corresponding

constructs were used to transform Medicago plants by hairy root transformation as described above. Nodules were harvested 2 weeks after being inoculated with S. meliloti (mCherry) in open pots and directly embedded in 7% low melting agarose. Semi-thin (70 µm) longitudinal sections were obtained using a vibratome microtome (VT1000S, Leica) and the sections were analyzed using a confocal microscope (Leica TCS SP8).

## Construct design

Constructs used for SYMREM1 localization assays, protoplast analyses, and yeast transformation are described in ref. ³⁶. The PS reporter (2x LactC2 domain), the plasma membrane marker (LTI6b), the SYMREM1 coding sequence and all the chimeric sequences were synthesized by Life Technologies, then cloned to expression vector by Golden Gate cloning<sup>85</sup>. For SYMREM1 expression and purification, the SYMREM1 coding sequence of *Medicago truncatula* was recombined into the Gateway (GW) compatible pDEST17 vector via LR-reaction. A pET303-EGFP-10xHis empty vector (provided by Dr. Nicole Gensch, University of Freiburg) was used for expression and purification of recombinant SYMREM1 (details see below). For AtREM2.1, the coding sequence template was provided by Julien Gronnier (University of Tübingen), then amplified and cloned into the expression vector by Golden Gate cloning<sup>86</sup>. Other AtREM expression vectors were built in a previous study<sup>36</sup>. All the used constructs and all primers used are listed in Supplementary Data 3 and Supplementary Data 4, respectively. The sequence data from this article can be found in phytozome (https://phytozome.jgi.doe.gov/) with gene IDs: SYMREM1 (Medtr8g097320), SYFO1 (Medtr5g036540). All gene IDs are listed in Supplementary Data 3.

For promoter activation studies, a sequence 1 kb upstream of the start codon of the *Medicago truncatula REM2.1* gene was amplified via PCR from genomic DNA (A17) and cloned into a pENTR/D-TOPO vector. The *ProMtREM2.1*:GUS construct was created by a LR-reaction of pKGWFS7-vector and pENTR/D-TOPO: *ProREM2.1*.

## Histochemical promoter analysis (GUS-staining), WGA staining, and microscopy

The activation patterns of the *MtREM2.1* promoter were analyzed via β-Glucuronidase (GUS) activity. Transgenic roots were stained in GUS-staining solution (0.1 M NaPO<sub>4</sub>; 1 mM EDTA; 1 mM K<sub>3</sub>Fe (CN)<sub>6</sub>; 1 mM K<sub>4</sub>Fe (CN)<sub>6</sub>; 1% Triton-X 100; 1 mM X-Gluc) at 37 °C for 4 h in the dark. For fluorescent visualization of fungal structures, colonized roots were fixed in 50% ethanol for at least 12 h and afterwards cleared for 2 days at room temperature in 10% KOH. After a washing step with distilled water, roots were incubated in 0.1 M HCl for 1 h at RT. Prior to the final staining, roots were washed with distilled water and rinsed once with 1x PBS (phosphate buffered saline; pH7.4). Roots were placed in a 1x PBS-WGA-AlexaFluor594 staining solution (0.2 µg/mL WGA-AlexaFluor594; Thermo Fisher Scientific, Germany) for at least 6 h at 4 °C in dark.

## Expression analysis

Total RNA extraction was performed according to the Spectrum Plant Total RNA Kit (Sigma-Aldrich, Germany) manual. Root material was grinded in liquid nitrogen and 100 mg per root sample was used for extraction. Extracted RNA was treated with DNase I, Amp Grade (Invitrogen, Germany). The absence of genomic DNA was verified via PCR. Synthesis of cDNA was performed with 700 ng of RNA in a total reaction volume of 20 µl using the Superscript III kit (Invitrogen, Germany). For qRT-PCR analysis, a Fast SYBR Green Master Mix (Applied Biosystems, Germany) was used in a 10 µl reaction volume. A CFX96TM Real-Time system (Bio-Rad, Germany) was used for PCR reactions and detection. Expression was normalized to *Ubiquitin*. At least three biological replicates were analyzed in technical duplicates per treatment.
## Transformation of *Nicotiana benthamiana* leaves and protoplast isolation

Transgenic *Agrobacterium tumefaciens* carrying the plasmids of interest were grown in LB liquid culture overnight at 28 °C with the appropriate antibiotics. The culture was centrifuged (4000 rpm, 2 min) and the pellet was resuspended in Agromix (10 mM MgCl₂; 10 mM MES/KOH pH 5.6; 150 µM Acetosyringone) to an OD₆₀₀ of 0.3. Bacteria were mixed with the silencing suppressor p19 before being incubated for 2 h at 25 °C in darkness and then infiltrated at the lower site of *N. benthamiana* leaves. Two days after infiltration, the transformed leaves were harvested to isolate protoplasts according as described earlier⁸⁷ with a small modification (PNT solution contained 300 mg/l of CaCl₂). All experiments using protoplasts were done at least three times independently.

## Microcapillary assay

Isolated protoplasts were embedded in 0.5% agarose on a cover of a Petri Dish. The injection set-up consisted of an inverted microscope (Zeiss Axiovert 135 TV) with a motor driven micromanipulator (LANG GmbH & Co. KG, Type: STM3) mounted at the right side of the stage. Femtotips injection needles (Eppendorf) were adapted by removing the sharp-pointed tip of the needle by hand, until obtaining a needle that could not penetrate the protoplast plasma-membrane.

## Confocal laser-scanning microscopy

Protoplast images were obtained using a Leica TCS SP8 confocal microscope equipped with a HC PL APO 20x/0.75 IMM CORR CS2 objective (Leica Microsystems, Mannheim, Germany), with the exception of images in Supplementary Fig. 7b and Fig. 5a that were taken with a Zeiss LSM880 Airyscan using a 63x/1.4 oil immersion lens. GFP was excited with a White Light Laser (WLL) at 488 nm and the emission detected at 500–550 nm. YFP was excited with a 514 nm laser line and detected at 520–555 nm. mCherry fluorescence was excited at 561 nm and emission was detected between 575–630 nm. Samples, co-expressing two fluorophores were imaged in sequential mode between frames. All images were taken as z-stacks (internal distance is 0.5 µm) and further analyses and projections were performed with either ImageJ/(Fiji) software⁸⁸ or Imaris.

For GUV experiments: Images were acquired by confocal fluorescence microscopy (Nikon Eclipse Ti-E inverted microscope using a Nikon A1R confocal laser scanning system with laser lines: 405 nm, 488 nm, 561 nm, 640 nm; 60x/1.49 oil immersion objective; Nikon Instruments, Inc.) and analyzed using the corresponding NIS Elements software (NIS Elements Confocal 5.20, Nikon Instruments, Inc.) and ImageJ (Fiji win64, Open source).

## Super resolution microscopy

Super resolution images were generated using an ELYRA7, structured-illumination-microscope with 3D Lattice SIM (Zeiss), equipped with a 63x/1.4 Oil DIC objective, a Pecon incubation chamber ensuring stable temperature and processed using the Zen black software package (Zeiss). Samples were stained as indicated and imaged as z-stack with 0,091 µm interval. Each stack was processed as single plane or maximum intensity projection (MIP) using Imaris software 64×9.0. Brightness and contrast were applied equally to all images prior to quantification. For colocalization analysis, we performed 2D line scan analysis on single plane and MIP images as indicated using Zen black software and Imaris.

## Protein expression and purification

*E. coli* BL21(DE3) cells were transformed with the plasmid pDEST17 encoding His-SYMREM1 (lipid strips, electron microscopy and actin sedimentation assay) and His-GFP-SYMREM1 (GUVs assembled with negatively charged lipids) proteins and the plasmid pET303-EGFP-10xHis encoding SYMREM1-GFP-His, SYREM1<sup>IDR</sup>-GFP-His and

SYMREM1<sup>C-term</sup>-GFP-His (GUVs assembled with Ni-NTA lipids). A single colony of transformed *E. coli* was transferred into LB medium and grown overnight for obtaining a pre-culture. Then, 40 ml of pre-culture was inoculated in 2 l of LB media and further grown at 37 °C. Protein expression was induced by 1 mM IPTG at an OD₆₀₀ of 0.5–0.6. Afterwards, cells were incubated overnight (about 20 h) at 25 °C. Cells were harvested by centrifugation at 6000 × g for 15 min. The cell pellet was resuspended in 100–150 ml Lysis buffer (20 mM HEPES, 500 mM NaCl, 20 mM imidazole, 10% glycerol, 1 mM EDTA, 1 mM Pefabloc, pH 7.2) and cells were passed through Constant Cell Disrupter (Constant Systems Limited). Cell debris was removed by centrifugation at 120,000 × g for 45 min. The cleared cell lysate was loaded onto IMAC columns (5 ml HisTrap FF) pre-equilibrated with Loading buffer A (20 mM HEPES, 500 mM NaCl, 20 mM imidazole, pH 7.2) and washed with 10 column volumes (CV) of Loading buffer A. Proteins were eluted with a linear gradient of imidazole from 20 to 450 mM in 15 CV. The eluted fractions were pooled and concentrated by spin filtration to 5 ml. Precipitated proteins were removed by an additional centrifugation for 10 min at 10,000 × g before loading onto gel-filtration column (HiLoad 16/60 Superdex 200 pg) equilibrated with PBS. Eluted fractions after gel-filtration were analyzed with SDS-PAGE. Those fractions containing pure His-SYMREM1 were pooled and concentrated by spin filtration to the working concentration. For Immunoblot Analysis: 2 weeks old nodules were harvested. Extracted proteins were separated on a 12% (w/v) SDS-PAGE gel and transferred overnight at 30 V to a PVDF (polyvinylidene difluoride) membrane. Membranes were then blocked with 5% (w/v) milk for 1 hour at room temperature before being hybridized with the SYREM1 peptide antibodies³⁵ at a dilution of 1:500 for 1.5 h at room temperature. Membranes were then washed with TBST three times for 10 min before incubating with a second antibody (anti-rabbit (Sigma), 1:2000 dilution, 1 h at room temperature). Prior to image acquisition, the membranes were washed again with TBST three time for 10 min.

mDial-ct proteins were purified from the *Escherichia coli* strain BL21 CP using the protocol for C-terminal constructs (mDial-ct) as previously described⁸⁹.

## Preparation of giant unilamellar vesicles (GUVs)

1,2-dioleoyl-*sn*-glycero-3-phosphocholine (DOPC), cholesterol and 1,2-dioleoyl-*sn*-glycero-3-[(*N*-5-amino-1-carboxypentyl) iminodiacetic acid] succinyl] [18:1 DGS-NTA(Ni); abbreviated Ni-NTA] were obtained from Avanti Polar Lipids (Alabaster, AL, USA). Atto 647 N 1,2-dioleoyl-*sn*-glycero-3-phosphoethanolamine (Atto 647N-DOPE) was purchased from Sigma-Aldrich Chemie GmbH (Darmstadt, Germany). The negatively charged lipids, brain PI(4,5)P₂ (L-α-phosphatidylinositol-4,5-bisphosphate) from porcine, 18:1 PI(3,5)P₂ (1,2-dioleoyl-*sn*-glycero-3-phospho-(1'-myo-inositol-3',5'-bisphosphate)), 18:1 PI(4)P (1,2-dioleoyl-*sn*-glycero-3-phospho-(1'-myo-inositol-4'-phosphate)) and Soy PI (L-α-phosphatidylinositol)) were purchased from Sigma-Aldrich.

As previously described⁹⁰ GUVs were formed by a classical electroformation protocol. In brief, solutions of lipids (at a concentration of 0.5 mg/mL), composed of DOPC, cholesterol, Atto 647N-DOPE (as membrane marker), and either the NiNTA or desired negatively charged lipid of choice (at 59.7:30:0.3:10 mol%, respectively), were prepared in chloroform and spread on indium tin oxide (ITO-) covered glass slides. For negatively charged lipids, lipid solutions were preheated to 50 °C before spreading on ITO slides.

To remove residual solvent, the slides were incubated under vacuum for at least one hour, for the negatively charged lipids, 50 °C temperature was applied. A chamber was assembled with two slides, filled with 318 mOsm L⁻¹ sucrose solution as formation buffer, and an AC electrical field with a voltage of 1 V mm⁻¹ was applied (to the chamber) for 2 h at room temperature for Ni-NTA solutions, while for negatively charged lipids this was done at 50 °C temperature, to ensure a homogeneous distribution of lipids⁹¹. GUVs
were observed in hand-made chambers⁹⁰ using 318 mOsm L⁻¹ PBS as imaging buffer.

### Actin co-sedimentation assay

Purified SYMREM1 (10 µM or 30 µM) or mDial-ct (10 µM) proteins were incubated with 10 µM G-actin (Cytoskeleton Inc.) for 1 h at RT in an actin-polymerization buffer (F-buffer: 5 mM Tris-HCl pH 7.5, 100 mM KCl, 1 mM MgCl₂, 0.2 mM CaCl₂, 0.2 mM EGTA, 0.2 mM ATP and 0.5 mM DTT). After ultracentrifugation at 80,000 × g for 30 min the supernatant was removed, as the non-pelleted fraction. The remaining pellet was then washed one time with F-buffer and subsequently resuspended in an equal volume compared to the supernatant. Laemml buffer was added to all lysates, boiled at 96 °C for 10 min and conducted to SDS-page analysis.

### PIP lipid strips

The recombinant His-SYMREM1 protein was used for lipid binding assays as previously described⁹².

### Phalloidin staining and SYMREM1 immunolocalization

Nodules (10 dpi) of *M. truncatula* were harvested in fixative solution (4% paraformaldehyde, 0.1% Triton X-100 and 10% DMSO in ASB solution (100 mM PIPES pH = 6.8, 1 mM MgCl₂, 1 mM CaCl₂)). Samples were then twice incubated under vacuum for 15 min and kept at room temperature for 4–6 h. Nodules were then embedded in 7% low melting agarose, sectioned into semi-thin (70 µm) longitudinal sections using a vibratome (VT1000S, Leica) and subsequently kept in a fixative solution for an additional hour at room temperature. After fixation, nodules sections were washed with ASB twice and kept in blocking solution (ASB supplied with 3% bovine serum albumin (BSA), 0.2% cold fish gelatin) for 30 min (room temperature). After blocking, nodule sections were incubated with ASB buffer additionally supplied with the SYMREM1 antibody (1:250) overnight. After incubation, sections were washed with ASB 3–5 times for 5 min before incubating with goat anti-rabbit Alexa488 (Molecular Probes) secondary antibodies (1:500 in ASB) and Phalloidin (Alexa Fluor™ 568) for at least 2 h at room temperature. Samples were then finally washed 3–5 times with ASB prior to imaging.

### Electron microscopy

Transmission electron microscopy (TEM) was performed on nodules harvested at (3 wpi). Nodules were cut longitudinally in half, immediately fixed in MTSB buffer⁹³ containing 2.5% glutaraldehyde and 4% p-Formaldehyde under vacuum for 15 min (twice), and stored at 4 °C in fixative solution until used further. After washing five times for 10 min each with buffer, nodules were post-fixed with 1% OsO₄ in H₂O at 4 °C for 2 h and again washed five times (10 min each) with H₂O at room temperature. The tissue was in block stained with 1% Uranyl Acetate for 1 hour in darkness, washed three times (10 min each) in H₂O, and dehydrated in EtOH/H₂O graded series (30%, 50%, 70%, 80%, 90%, 95% 15 min each). Final dehydration was achieved by incubating the samples twice in absolute EtOH (30 min each) followed by incubation in dehydrated acetone twice (30 min each). Embedding of the samples was performed by gradually infiltrating them with Epoxy resin (Agar 100) mixed with acetone at 1:3, 1:1 and 3:1 ratio for 12 h each, and finally in pure Epoxy resin for 48 h with resin changes every 12 h. Polymerization was carried out at 60 °C for 48 h. Ultrathin sections of -70 nm were obtained with a Reichert-Jung ultra-microtome and collected in TEM slot grids. Images were acquired with a Philips CM 10 transmission electron microscope coupled to a Gatan BioScan 792 CCD camera at 80 kV acceleration voltage.

Scanning Electron Microscopy was carried out on freshly isolated protoplasts and on longitudinal vibratome sections (70 µm) of nodules collected after 3 wpi. The material was immediately fixed, dehydrated in graded EtOH series until 100% EtOH, and critical point dried in

absolute EtOH-CO₂. Dried material was mounted on carbon tabs and coated with 5 nm platinum. Imaging of samples was performed using a Hitachi S-4800 microscope at 5 kV acceleration voltage.

Negative staining of purified SYMREM1 protein was performed by applying 5 µl protein solution to glow-discharged 400 Cu mesh carbon grids for 10 min, blotting and negatively staining using 2% (w/v) uranyl acetate. Images were recorded under low-dose conditions on a Talos F200C transmission electron microscope operated at 200 kV and equipped with a Ceta 16 M camera. Micrographs were taken at a nominal magnification of 73,000x. A total of 389 segments were manually selected using RELION-3.1.0⁹⁴. The defocus and astigmatism of the images were determined with CTFFIND4.1⁹⁵ and numerical phase-flipping was done to correct for effects of the contrast transfer function using RELION-3.1.0. Image processing was done using IMAGIC-5⁹⁶. Particle images were band pass filtered between 400 and 10 Å, normalized and centered by iteratively aligning them to a vertically oriented class average. Class averages containing 5–10 images were obtained by four rounds of classification based on multivariate statistical analysis, followed by multi-reference alignment using homogenous classes as new references.

### Protein 3D prediction

In-silico structure predictions were made with DeepMind AlphaFold 2.1.1 in multimer mode⁴⁴,⁴⁵. Electrostatic surface potentials were generated with ABPS⁹⁷, figures were created with PyMOL (Schrödinger LLC).

### Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

### Data availability

All data are available in the main text or the supplementary materials.

Source Data are provided Source data are provided with this paper.

### References

1. Amor, B. B. et al. The NFP locus of Medicago truncatula controls an early step of Nod factor signal transduction upstream of a rapid calcium flux and root hair deformation. *Plant J.* **34**, 495–506 (2003).

2. Limpens, E. et al. LysM domain receptor kinases regulating rhizobial Nod factor-induced infection. *Science* **302**, 630–633 (2003).

3. Smit, P. et al. Medicago LYK3, an entry receptor in rhizobial nodulation factor signaling. *Plant Physiol.* **145**, 183–191 (2007).

4. Kawaharada, Y. et al. Receptor-mediated exopolysaccharide perception controls bacterial infection. *Nature* **523**, 308–312 (2015).

5. Wais, R. J. et al. Genetic analysis of calcium spiking responses in nodulation mutants of Medicago truncatula. *Proc. Natl Acad. Sci. USA* **97**, 13407–13412 (2000).

6. Levy, J. et al. A putative Ca2+ and calmodulin-dependent protein kinase required for bacterial and fungal symbioses. *Science* **303**, 1361–1364 (2004).

7. Messinese, E. et al. A novel nuclear protein interacts with the symbiotic DMI3 calcium- and calmodulin-dependent protein kinase of Medicago truncatula. *Mol. Plant Microbe Interact.* **20**, 912–921 (2007).

8. Singh, S., Katzer, K., Lambert, J., Cerri, M. & Parniske, M. CYCLOPS, a DNA-binding transcriptional activator, orchestrates symbiotic root nodule development. *Cell Host Microbe* **15**, 139–152 (2014).

9. Yano, K. et al. CYCLOPS, a mediator of symbiotic intracellular accommodation. *Proc. Natl Acad. Sci. USA* **105**, 20540–20545 (2008).

10. Timmers, A. C., Auriac, M. C. & Truchet, G. Refined analysis of early symbiotic steps of the Rhizobium-Medicago interaction in relationship with microtubular cytoskeleton rearrangements. *Development* **126**, 3617–3628 (1999).
11. Yokota, K. et al. Rearrangement of actin cytoskeleton mediates invasion of Lotus japonicus roots by Mesorhizobium loti. *Plant Cell* **21**, 267–284 (2009).

12. Zepeda, I. et al. Visualization of highly dynamic F-actin plus ends in growing phaseolus vulgaris root hair cells and their responses to Rhizobium etli nod factors. *Plant Cell Physiol.* **55**, 580–592 (2014).

13. Vassileva, V. N., Kouchi, H. & Ridge, R. W. Microtubule dynamics in living root hairs: transient slowing by lipochitin oligosaccharide nodulation signals. *Plant Cell* **17**, 1777–1787 (2005).

14. Liang, P. et al. Formin-mediated bridging of cell wall, plasma membrane, and cytoskeleton in symbiotic infections of Medicago truncatula. *Curr. Biol.* **31**, 2712–2719.e2715 (2021).

15. Fournier, J. et al. Remodeling of the infection chamber before infection thread formation reveals a two-step mechanism for rhizobial entry into the host legume root hair. *Plant Physiol.* **167**, 1233–1242 (2015).

16. Gage, D. J. Infection and invasion of roots by symbiotic, nitrogen-fixing rhizobia during nodulation of temperate legumes. *Microbiol. Mol. Biol. Rev.* **68**, 280–300 (2004).

17. Monahan-Giovanelli, H., Pinedo, C. A. & Gage, D. J. Architecture of infection thread networks in developing root nodules induced by the symbiotic bacterium Sinorhizobium meliloti on Medicago truncatula. *Plant Physiol.* **140**, 661–670 (2006).

18. Mellor, R. B. Bacteroids in the rhizobium-legume symbiosis inhabit a plant internal lytic compartment - implications for other microbial endosymbioses. *J. Exp. Bot.* **40**, 831–839 (1989).

19. Verma, D. P. S. & Hong, Z. L. Biogenesis of the peribacteroid membrane in root nodules. *Trends Microbiol.* **4**, 364–368 (1996).

20. Brewin, N. J. Plant cell wall remodelling in the rhizobium-legume symbiosis. *Crit. Rev. Plant Sci.* **23**, 293–316 (2004).

21. Robledo, M. et al. Rhizobium cellulase CelC2 is essential for primary symbiotic infection of legume host roots. *Proc. Natl Acad. Sci. USA* **105**, 7064–7069 (2008).

22. Kovtun, O., Dickson, V. K., Kelly, B. T., Owen, D. J. & Briggs, J. A. G. Architecture of the AP2/clathrin coat on the membranes of clathrin-coated vesicles. *Sci. Adv.* **6**, eaba8381 (2020).

23. Robinson, M. S. Forty years of clathrin-coated vesicles. *Traffic* **16**, 1210–1238 (2015).

24. Adam, J., Basnet, N. & Mizuno, N. Structural insights into the cooperative remodeling of membranes by amphiphysin/BIN1. *Sci. Rep.* **5**, 15452 (2015).

25. Raffaele, S., Mongrand, S., Gamas, P., Niebel, A. & Ott, T. Genome-wide annotation of remorins, a plant-specific protein family: evolutionary and functional perspectives. *Plant Physiol.* **145**, 593–600 (2007).

26. Liang, P. et al. Symbiotic root infections in Medicago truncatula require remorin-mediated receptor stabilization in membrane nanodomains. *Proc. Natl Acad. Sci. USA* **115**, 5289–5294 (2018).

27. Gronnier, J. et al. Structural basis for plant plasma membrane protein dynamics and organization into functional nanodomains. *eLife* **6**, e26404 (2017).

28. Huang, D. et al. Salicylic acid-mediated plasmodesmal closure via Remorin-dependent lipid organization. *Proc. Natl Acad. Sci. USA* **116**, 21274–21284 (2019).

29. Raffaele, S. et al. Remorin, a Solanaceae protein resident in membrane rafts and plasmodesmata, impairs potato virus X movement. *Plant Cell* **21**, 1541–1555 (2009).

30. Marin, M., Thallmair, V. & Ott, T. The intrinsically disordered N-terminal region of AtREM1.3 remorin protein mediates protein-protein interactions. *J. Biol. Chem.* **287**, 39982–39991 (2012).

31. Martinez, D. et al. Coiled-coil oligomerization controls localization of the plasma membrane REMORINS. *J. Struct. Biol.* **206**, 12–19 (2019).

32. Jaillais, Y. & Ott, T. The nanoscale organization of the plasma membrane and its importance in signaling: a proteolipid perspective. *Plant Physiol.* **182**, 1682–1696 (2020).

33. Haney, C. H. et al. Symbiotic Rhizobia bacteria trigger a change in localization and dynamics of the Medicago truncatula receptor kinase LYK3. *Plant Cell* **23**, 2774–2787 (2011).

34. Haney, C. H. & Long, S. R. Plant flotillins are required for infection by nitrogen-fixing bacteria. *Proc. Natl Acad. Sci. USA* **107**, 478–483 (2010).

35. Lefebvre, B. et al. A remorin protein interacts with symbiotic receptors and regulates bacterial infection. *Proc. Natl Acad. Sci. USA* **107**, 2343–2348 (2010).

36. Konrad, S. S. et al. S-acylation anchors remorin proteins to the plasma membrane but does not primarily determine their localization in membrane microdomains. *N. Phytol.* **203**, 758–769 (2014).

37. Perraki, A. et al. Plasma membrane localization of solanum tuberosum remorin from group 1, homolog 3 is mediated by conformational changes in a novel C-terminal anchor and required for the restriction of potato virus X movement. *Plant Physiol.* **160**, 624–637 (2012).

38. Ma, Z. et al. Membrane nanodomains modulate formin condensation for actin remodeling in Arabidopsis innate immune responses. *Plant Cell* **34**, 374–394 (2022).

39. Platre, M. P. et al. A combinatorial lipid code shapes the electrostatic landscape of plant endomembranes. *Dev. Cell* **45**, 465–480.e411 (2018).

40. Noack, L. C. & Jaillais, Y. Functions of anionic lipids in plants. *Annu. Rev. Plant Biol.* **71**, 71–102 (2020).

41. Moling, S. et al. Nod factor receptors form heteromeric complexes and are essential for intracellular infection in medicago nodules. *Plant Cell* **26**, 4188–4199 (2014).

42. Bariola, P. A. et al. Remorins form a novel family of coiled coil-forming oligomeric and filamentous proteins associated with apical, vascular and embryonic tissues in plants. *Plant Mol. Biol.* **55**, 579–594 (2004).

43. Yue, J., Li, C., Liu, Y. & Yu, J. A remorin gene SiREM6, the target gene of SiARDP, from foxtail millet (Setaria italica) promotes high salt tolerance in transgenic Arabidopsis. *PLoS ONE* **9**, e100772 (2014).

44. Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. *Nature* **596**, 583–589 (2021).

45. Evans, R. et al. Protein complex prediction with AlphaFold-Multimer. Preprint at bioRxiv https://doi.org/10.1101/2021.10.04.463034 (2022).

46. Tóth, K. et al. Functional domain analysis of the remorin protein LjSYMREM1 in Lotus japonicus. *PLoS ONE* **7**, e30817 (2012).

47. Ovchinnikova, E. et al. IPD3 controls the formation of nitrogen-fixing symbiosomes in pea and Medicago Spp. *Mol. Plant Microbe Interact.* **24**, 1333–1344 (2011).

48. Lee, C. Y. et al. Regulation of actin catch-slip bonds with a RhoA-formin module. *Sci. Rep.* **6**, 35058 (2016).

49. Durand-Smet, P., Spelman, T. A., Meyerowitz, E. M. & Jonsson, H. Cytoskeletal organization in isolated plant cells under geometry control. *Proc. Natl Acad. Sci. USA* **117**, 17399–17408 (2020).

50. Carrere, S., Verdier, J. & Gamas, P. MtExpress, a comprehensive and curated RNAseq-based gene expression atlas for the model legume Medicago truncatula. *Plant Cell Physiol.* **62**, 1494–1500 (2021).

51. Gutjahr, C. & Parniske, M. Cell and developmental biology of arbuscular mycorrhiza symbiosis. *Annu. Rev. Cell Dev. Biol.* **29**, 593–617 (2013).

52. Colin, L. et al. Cortical tension overrides geometrical cues to orient microtubules in confined protoplasts. *Proc. Natl Acad. Sci. USA* **117**, 32731–32738 (2020).

53. Cosgrove, D. J. Plant cell wall extensibility: connecting plant cell growth with cell wall structure, mechanics, and the action of wall-modifying enzymes. *J. Exp. Bot.* **67**, 463–476 (2016).
54. Szymanski, D. & Staiger, C. J. The actin cytoskeleton: functional arrays for cytoplasmic organization and cell shape control. *Plant Physiol.* **176**, 106–118 (2018).

55. McCullough, J. et al. Structure and membrane remodeling activity of ESCRT-III helical polymers. *Science* **350**, 1548–1551 (2015).

56. Moser von Filseck, J. et al. Anisotropic ESCRT-III architecture governs helical membrane tube formation. *Nat. Commun.* **11**, 1516 (2020).

57. Fan, L., Li, R., Pan, J., Ding, Z. & Lin, J. Endocytosis and its regulation in plants. *Trends Plant Sci.* **20**, 388–397 (2015).

58. Glebov, O. O., Bright, N. A. & Nichols, B. J. Flotillin-1 defines a clathrin-independent endocytic pathway in mammalian cells. *Nat. Cell Biol.* **8**, 46–54 (2006).

59. Li, R. et al. A membrane microdomain-associated protein, Arabidopsis Flot1, is involved in a clathrin-independent endocytic pathway and is required for seedling development. *Plant Cell* **24**, 2105–2122 (2012).

60. Simunovic, M., Evergren, E., Callan-Jones, A. & Bassereau, P. Curving cells inside and out: roles of BAR domain proteins in membrane shaping and its cellular implications. *Annu. Rev. Cell Dev. Biol.* **35**, 111–129 (2019).

61. Zhuang, X. et al. A BAR-domain protein SH3P2, which binds to phosphatidylinositol 3-phosphate and ATG8, regulates autophagosome formation in Arabidopsis. *Plant Cell* **25**, 4596–4615 (2013).

62. Ahn, G. et al. SH3 domain-containing protein 2 plays a crucial role at the step of membrane tubulation during cell plate formation. *Plant Cell* **29**, 1388–1405 (2017).

63. Kozlov, M. M. et al. Mechanisms shaping cell membranes. *Curr. Op. Cell Biol.* **29**, 53–60 (2014).

64. Campelo, F., McMahon, H. T. & Kozlov, M. M. The hydrophobic insertion mechanism of membrane curvature generation by proteins. *Biophys. J.* **95**, 2325–2339 (2008).

65. Mim, C. et al. Structural basis of membrane bending by the N-BAR protein endophilin. *Cell* **149**, 137–145 (2012).

66. McDonald, N. A. et al. The tubulation activity of a fission yeast F-BAR protein is dispensable for its function in cytokinesis. *Cell Rep.* **14**, 534–546 (2016).

67. Wei, Z. et al. Plasmodesmata-like intercellular connections by plant remorin in animal cells. Preprint at *bioRxiv* https://doi.org/10.1101/791137 (2020).

68. Bornschlogl, T. How filopodia pull: what we know about the mechanics and dynamics of filopodia. *Cytoskeleton* **70**, 590–603 (2013).

69. Gui, J., Zheng, S., Shen, J. & Li, L. Grain setting defect1 (GSD1) function in rice depends on S-acylation and interacts with actin 1 (OsACT1) at its C-terminal. *Front. Plant Sci.* **6**, 804 (2015).

70. Stachowiak, J. C. et al. Membrane bending by protein-protein crowding. *Nat. Cell Biol.* **14**, 944–949 (2012).

71. Jarsch, I. K., Daste, F. & Gallop, J. L. Membrane curvature in cell biology: an integration of molecular mechanisms. *J. Cell Biol.* **214**, 375–387 (2016).

72. Boisson-Dernier, A. et al. Agrobacterium rhizogenes-transformed roots of Medicago truncatula for the study of nitrogen-fixing and endomycorrhizal symbiotic associations. *Mol. Plant Microbe Interact.* **14**, 695–700 (2001).

73. Camacho, C. et al. BLAST+: architecture and applications. *BMC Bioinforma.* **10**, 421 (2009).

74. Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. *Mol. Biol. Evol.* **30**, 772–780 (2013).

75. Capella-Gutierrez, S., Silla-Martinez, J. M. & Gabaldon, T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. *Bioinformatics* **25**, 1972–1973 (2009).

76. Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., von Haeseler, A. & Jermiin, L. S. ModelFinder: fast model selection for accurate phylogenetic estimates. *Nat. Methods* **14**, 587 (2017).

77. Minh, B. Q. et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era (vol 37, pg 1530, 2020). *Mol. Biol. Evol.* **37**, 2461–2461 (2020).

78. Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. *Syst. Biol.* **59**, 307–321 (2010).

79. Hoang, D. T., Chernomor, O., von Haeseler, A., Minh, B. Q. & Vinh, L. S. UFBoot2: improving the ultrafast bootstrap approximation. *Mol. Biol. Evol.* **35**, 518–522 (2018).

80. Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. *Nucl. Acids Res.* **32**, 1792–1797 (2004).

81. Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. *Nucl. Acids Res.* **49**, W293–W296 (2021).

82. Wertheim, J. O., Murrell, B., Smith, M. D., Pond, S. L. K. & Scheffler, K. RELAX: detecting relaxed selection in a phylogenetic framework. *Mol. Biol. Evol.* **32**, 820–832 (2015).

83. Smith, M. D. et al. Less is more: an adaptive branch-site random effects model for efficient detection of episodic diversifying selection. *Mol. Biol. Evol.* **32**, 1342–1353 (2015).

84. Wright, E. S. Using DECIPHER v2.0 to Analyze Big Biological Sequence Data in R. *R. J.* **8**, 352–359 (2016).

85. Weber, E., Engler, C., Gruetzner, R., Werner, S. & Marillonnet, S. A modular cloning system for standardized assembly of multigene constructs. *PLoS ONE* **6**, e16765 (2011).

86. Binder, A. et al. A modular plasmid assembly kit for multigene expression, gene silencing and silencing rescue in plants. *PLoS ONE* **9**, e88218 (2014).

87. Schnorf, M. et al. An improved approach for transformation of plant cells by microinjection: molecular and genetic analysis. *Transgenic Res.* **1**, 23–30 (1991).

88. Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. *Nat. Methods* **9**, 676–682 (2012).

89. Li, F. & Higgs, H. N. Dissecting requirements for auto-inhibition of actin nucleation by the formin, mDia1. *J. Biol. Chem.* **280**, 6986–6992 (2005).

90. Madl, J., Villringer, S. & Römer, W. In *Chemical and Synthetic Approaches in Membrane Biology* (ed. Shukla, A. K.) 17–36 (Humana Press, 2016).

91. Roux, A. et al. Membrane curvature controls dynamin polymerization. *Proc. Natl Acad. Sci. USA* **107**, 4141–4146 (2010).

92. Su, C. et al. The Medicago truncatula DREPP protein triggers microtubule fragmentation in membrane nanodomains during symbiotic infections. *Plant Cell* **32**, 1689–1702 (2020).

93. Miart, F. et al. Spatio-temporal analysis of cellulose synthesis during cell plate formation in Arabidopsis. *Plant J.* **77**, 71–84 (2014).

94. Zivanov, J. et al. New tools for automated high-resolution cryo-EM structure determination in RELION-3. *eLife* **7**, e42166 (2018).

95. Rohou, A. & Grigorieff, N. CTFFIND4: fast and accurate defocus estimation from electron micrographs. *J. Struct. Biol.* **192**, 216–221 (2015).

96. vanHeel, M., Harauz, G., Orlova, E. V., Schmidt, R. & Schatz, M. A new generation of the IMAGIC image processing system. *J. Struct. Biol.* **116**, 17–24 (1996).

97. Jurrus, E. et al. Improvements to the APBS biomolecular solvation software suite. *Protein Sci.* **27**, 112–128 (2018).

## Acknowledgements

We would like to thank Carmen Schubert and Rosula Hinnenberg for their excellent technical help and the entire Ott lab team for fruitful
discussions and providing their individual expertise throughout the course of the project. We also thank the staff of the Life Imaging Center (LIC) in the Hilde Mangold House (HMH) of the Albert-Ludwigs-University of Freiburg for help with their confocal microscopy resources, and the excellent support in image recording. Special thanks also to Norbert Roos and Jens Wohlmann at the Electron Microscopy Facility, Department of Biosciences, University of Oslo, Norway, for helping with the SEM sample preparation and imaging and to Falk Tauber (University of Freiburg, Cluster of Excellence livMatS) for designing and printing a tentative indentation support. Many thanks also to Nicole Gensch (BIOSS, University of Freiburg) for providing the pET303-EGFP-10xHis vector, Maria Harrison (Boyce Thompson Institute, Ithaca, USA) for providing the Medicago ABD-mCherry line and Julien Gronnier (University of Tübingen) for providing the AtREM2.1 CDS and critically reading the manuscript. J.K., C.L., and P.M.D. belong to the LRSV, which is part of the TULIP LABEX (ANR-10-LABX-41). The microscopes are operated by the Microscopy and Image Analysis Platform (MIAP) and the Life Imaging Center (LIC), Freiburg. Engineering Nitrogen Symbiosis for Africa (ENSA) project currently supported through a grant to the University of Cambridge by the Bill & Melinda Gates Foundation (OPP1172165) and UK government's Department for International Development (DFID) (T.O., P.M.D.). Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) 431626755 (T.O.), 442219341 (P.W.). DFG under Germany's Excellence Strategy grant CIBSS - EXC-2189 - Project ID 39093984 (T.O./C.H./W.W./R.G./W.R.). China Scholarship Council (CSC) grants 201708080016 and 201506350004 (C.S./P.L.). DFG project number 414136422 (CLSM; T.O.), DFG project number 426849454 (TEM; T.O.) and DFG project number 406260942 (cryo-TEM; P.W.). European Union's Horizon 2020 Research and Innovation Program under the Marie Skłodowska-Curie grant agreement synBIOcarb (No. 814029; W.R.). Ministry for Science, Research and Arts of the State of Baden-Württemberg (Az: 33-7532.20; W.R.)

## Author contributions

Conceptualization, C.S., W.R. and T.O.; Investigation, C.S., M.R.-F., B.L., N.N., C.H.-R., P.L., E.S., E.V.M., N.M.G., J. Kn, H.W., L.S., J. Ke, C.L., A.A.M.F., K.E.G., E.M., C.P., P.W., T.S., P.M.D., O.E. and T.O.; Writing—Original Draft, C.S., M.R.-F., and T.O.; Writing—Review & Editing, C.S., M.R.-F., B.L., N.N., C.H.-R., P.L., E.S., E.V.M., N.M.G., J. Kn, H.W., L.S., J. Ke, C.L., A.A.M.F., K.E.G., E.M., C.P., C.H., W.W., P.W., T.S., P.M.D., O.E., R.G., W.R. and T.O.; Funding Acquisition, C.H., W.W., P.W., T.S., P.M.D., O.E., R.G., W.R. and T.O.; Supervision, C.H., W.W., P.W., T.S., P.M.D., O.E., R.G., W.R. and T.O.

## Funding

Open Access funding enabled and organized by Projekt DEAL.

## Competing interests

A patent application has been submitted (U.S. Provisional Application Serial No.: 63/245662) on behalf of C.S., L.S., W.R. and T.O. concerning the use of SYMREM1 and other proteins for regulating symbiotic infections by altering membrane topologies. All other authors declare no further competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-35976-5.

**Correspondence** and requests for materials should be addressed to Thomas Ott.

**Peer review information** *Nature Communications* thanks Patricia Bas-sereau, Dong Wang and Xiaohong Zhuang for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdic-tional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023