In plants, a complex mixture of solutes and macromolecules is transported by the phloem. Here, we examined how solutes and macromolecules are separated when they exit the phloem during the unloading process. We used a combination of approaches (non-invasive imaging, 3D-electron microscopy, and mathematical modelling) to show that phloem unloading of solutes in Arabidopsis roots occurs through plasmodesmata by a combination of mass flow and diffusion (convective phloem unloading). During unloading, solutes and proteins are diverted into the phloem-pole pericycle, a tissue connected to the protophloem by a unique class of ‘funnel plasmodesmata’. While solutes are unloaded without restriction, large proteins are released through funnel plasmodesmata in discrete pulses, a phenomenon we refer to as ‘batch unloading’. Unlike solutes, these proteins remain restricted to the phloem-pole pericycle. Our data demonstrate a major role for the phloem-pole pericycle in regulating phloem unloading in roots.https://doi.org/10.7554/eLife.24125.001
A mechanism called photosynthesis allows plants to use energy from sunlight to make sugars from carbon dioxide gas and water. These sugars can then be used as fuel, or as building blocks for wood and other plant structures. Every part of the plant requires sugars, but most photosynthesis happens in the leaves and stems, so the sugars need to be able to move around the plant to wherever they are needed.
Phloem tubes form a network that transports sugar, proteins and other molecules around the plant within a fluid known as sap. Because this network is so extensive, it is very difficult to study, which has left researchers with major questions about how it works. For example, it is not clear how the sugar and other molecules leave the phloem when they reach their destination.
Ross-Elliot et al. used a combination of microscopy and mathematical modeling to investigate how sugars and other molecules leave the phloem in the roots of a plant called Arabidopsis thaliana. The experiments show that these molecules move directly into cells within a neighboring tissue called the phloem-pole pericycle via pores known as funnel plasmodesmata.
Ross-Elliot et al. incorporated the experimental data into a mathematical model of phloem unloading. This model suggests that sugars and other small molecules move freely through the funnel plasmodesmata, but large proteins pass through these pores in pulses. Future challenges include finding out exactly how plants control phloem unloading and to investigate whether it is possible to modify the delivery of specific molecules to particular parts of the plant.https://doi.org/10.7554/eLife.24125.002
In plants, the products of photosynthesis in green tissues are delivered by the phloem to distant organs where they are utilized in growth or storage (Turgeon and Wolf, 2009; De Schepper et al., 2013). Modern agriculture aims to maximize the amounts of carbon-based products allocated to storage organs such as grains and tubers, structures that act as carbon ‘sinks’ for the assimilates delivered by the phloem. The process by which solutes exit the phloem is termed phloem unloading, and is a central target for regulating the flux of carbon into sink tissues (Ham and Lucas, 2014). In addition to assimilates, phloem sap also contains numerous proteins and RNAs (Kehr, 2006; Atkins et al., 2011; Batailler et al., 2012; Turnbull and Lopez-Cobollo, 2013). Unloading must therefore combine the seemingly antagonist functions of high selectivity with large permeability in order to allocate assimilates to growth zones while controlling the movement of macromolecules. The mechanisms that enable the precise and coordinated unloading of phloem-mobile compounds remain unknown.
Solute flow occurs through sieve elements (SEs), elongated cells connected to each other by perforated end walls called sieve plates. During differentiation SEs lose most of their cellular components, including their nucleus (Oparka and Turgeon, 1999; Furuta et al., 2014). Ontogenetically related, metabolically active companion cells (CCs) support the adjacent enucleate SEs throughout their lifespan (Lucas et al., 1996; Pritchard, 1996; van Bel and Knoblauch, 2000; Lalonde et al., 2001; Otero and Helariutta, 2017). The pressure flow model of phloem transport, originally proposed by Münch (1930), envisages that an osmotically generated pressure differential drives the bulk flow through the SEs that connect photosynthetic tissues (sources) with those in which carbon consumption occurs (sinks). The loading of solutes in source tissues results in a high osmotic potential within SEs. The reduction of turgor pressure in sink organs due to carbon consumption leads to a pressure gradient that provides the energy to overcome viscous resistance within the SEs, resulting in a passive phloem flow from source to sink (Münch, 1930). Phloem sap collected from excised aphid stylets contains a complex mixture of macromolecules and low-molecular weight solutes (Atkins et al., 2011). In a recent study, Paultre et al. (2016) showed that many proteins, including those with targeting sequences, can move across a graft union and be unloaded near the root tip. Significantly, these proteins entered a post-phloem domain beyond which their movement was restricted. This observation begs the question as to how the phloem is able to discriminate between macromolecules and solutes during unloading. Besides the central function in resource allocation, it is now well established that the phloem also serves as network for transmission of chemical (Kramer and Boyer, 1995a; Mullendore et al., 2015) and electrical (Hedrich et al., 2016) signals.
Phloem unloading in actively growing tissues such as the shoot or root apex occurs through the protophloem, a transient tissue that connects the conducting phloem with the receiver cells in sink tissues (Oparka et al., 1994). Initial investigations of phloem unloading in the root tip of Arabidopsis (Oparka et al., 1994, 1995; Wright and Oparka, 1996) provided evidence that unloading occurs through plasmodesmata (PD), the specialized pores that connect plant cells. Due to cell division and growth in the apical region of the root the demand for assimilates is high. Protophloem sieve elements (PSEs) become mature in such regions to accommodate this demand while the neighboring cells are still differentiating (Furuta et al., 2014). Because PSEs lose their nucleus, they cannot divide to keep pace with the growth of the neighboring cells. Therefore, they become increasingly elongated until they become inactive in transport and eventually obliterated (Erwin and Evert, 1967; Eleftheriou and Tsekos, 1982; Furuta et al., 2014). In the elongation zone of the root, solutes are transferred laterally from the metaphloem SEs (MSEs) to the PSEs, allowing phloem continuity between source and sink tissues (Stadler et al., 2005; Winter et al., 1992).
A currently favored hypothesis of phloem unloading is the ‘high-pressure manifold model’ proposed by Fisher (2000), recently evaluated by Patrick (2013). A central element of this model is that a low pressure gradient occurs along the flow path, with a steep drop in pressure between the PSEs and surrounding cells in the phloem unloading zone. In this scenario, allocation of carbon is controlled by the lateral hydraulic conductance in the unloading zone. Recent phloem turgor measurements in morning glory, however, did not support this model, as the bulk of the pressure is consumed by friction within the SEs, and only small pressure gradients are available for unloading (Knoblauch et al., 2016). The Fisher (2000) model also does not explain how both small solutes and macromolecules can leave the phloem simultaneously. Paultre et al., 2016 suggested recently that the removal of mobile proteins into a post-phloem domain may be necessary to prevent the terminal PSEs from becoming occluded, an event that would lead to dissipation of the turgor gradient between source and sink. Unfortunately, many of the factors that determine phloem unloading are not well studied. This is because the phloem in most sinks is difficult to access as it is embedded in an opaque layer of tissue. Phloem transport ceases immediately when the source is detached from the sink. Therefore, the function of the phloem can only be studied in situ, requiring new approaches for dissecting the factors that regulate phloem unloading.
The many unknowns surrounding phloem unloading in plants prompted us to conduct a detailed structure/function study of the terminal PSEs in the phloem unloading zone which, in Arabidopsis, is amenable to non-invasive imaging (Oparka et al., 1994; Knoblauch et al., 2015). We combined a detailed ultrastructural analysis of the cellular interfaces in the phloem unloading zone with the kinetics of phloem unloading of fluorescent solutes and macromolecules obtained by real-time imaging of growing roots. This analysis allowed us to derive new quantitative data on the factors that regulate phloem unloading in Arabidopsis. We report on the presence of a unique class of ‘funnel plasmodesmata’ that are involved specifically in the unloading of molecules into the phloem-pole pericycle. We show, by mathematical modelling, that phloem unloading of small solutes from PSEs is convective, i.e. it occurs continuously by a combination of mass flow and diffusion. In contrast, we find that macromolecules are unloaded in discrete pulses, a phenomenon we refer to as ‘batch unloading’. These macromolecules are diverted specifically into the two phloem-pole pericycle cells that abut each PSE where they are filtered out from the unloaded solutes.
In Arabidopsis roots, phloem unloading occurs exclusively from the protophloem, a short-lived tissue that is functional in the zone of root elongation (Oparka et al., 1994). Published images of phloem unloading in Arabidopsis give the impression that the unloading zone is a relatively broad region (Oparka et al., 1994; Wright and Oparka, 1996; Knoblauch et al., 2015). However, such images represent a ‘snapshot’ of unloading, taken at a defined time point following the application of fluorescent solutes to the leaf. Carboxyfluorescein diacetate (CFDA) is the most widely utilized phloem-mobile probe (Oparka et al., 1994; Wright and Oparka, 1996; Knoblauch and van Bel, 1998). This probe is non-fluorescent when applied to source leaves but is subsequently cleaved by endogenous esterases to produce fluorescent, membrane-impermeant, carboxyfluorescein (CF; Knoblauch et al., 2015). The dye travels with the phloem translocation stream to sink tissues, where it can be visualized (Figure 1A). When studied in real time, the dye characteristically shows preferential movement outwards into the cortex relative to the stele (Figure 1A; see also Oparka et al., 1994). Unloading of the dye in the root tip indicates a symplastic (plasmodesmata-mediated) pathway (Oparka et al., 1994; Wright and Oparka, 1996; Knoblauch et al., 2015). The cells in the unloading zone sequester the dye rapidly into their vacuoles where it becomes trapped. This results in an increasing fluorescence within cells over time (Wright and Oparka, 1996). Due to root growth, the cells initially involved in unloading move out of the phloem-unloading zone basipetally, but maintain their fluorescence due to the presence of dye in their vacuoles. Therefore, the apparent unloading zone broadens over time as root growth progresses, obscuring the current site of unloading. We therefore sought to define precisely the dimensions of the true phloem-unloading zone.
Esculin is a naturally fluorescent, glucosylated coumarin derivative recently introduced as phloem-mobile probe (Knoblauch et al., 2015). Unlike CFDA, it is loaded into the phloem by the sucrose transporter, SUC2, in the CCs of source tissues (Gora et al., 2012; Knoblauch et al., 2015). Sequestration of this probe is minimal and occurs only when high concentrations of the probe are applied (Figure 1B). This feature allows for extended acquisition of time-lapse movies. Unlike CF, esculin can also be detected clearly in lines expressing GFP (Knoblauch et al., 2015).
To define the developmental stage at which differentiating PSEs become integrated into the unloading zone, we used transgenic Arabidopsis lines expressing GFP targeted to the sieve-element ER (HDEL-GFP) under control of the MtSEO2 promoter, a SE-specific promoter (Froelich et al., 2011; Knoblauch and Peters, 2010); Figure 1B–E). This transgenic line clearly demarcates PSEs in the early stages of differentiation and is an excellent marker of the nuclear membrane (Figure 1C–E). We loaded source leaves of this line with esculin and acquired time-lapse movies of unloading in the terminal PSEs. Esculin did not enter differentiating PSEs that were symplastically isolated from the rest of the protophloem file. These cells still had a fully intact nuclear membrane (Figure 1C, Video 1). Once the nuclear membrane started to degenerate the sieve-plate pores opened rapidly, allowing esculin to enter the cell, which then became an integral member of the protophloem-unloading zone (Figure 1D and E). For convenience of reference, we refer to this protophloem cell as ‘PSE zero’ (yellow dashed arrow in Figure 1D). Remnants of the nuclear membrane remained for some time in PSEs that had been newly incorporated into the phloem-unloading zone (Figure 1C–E). These observations can be seen in Video 1. Our data suggest that nucleate, differentiating PSEs are isolated from the translocation stream. However, the degeneration of the nucleus and the opening of the sieve-plate pores are closely related events that lead to incorporation of PSE zero into the phloem-unloading zone. At this point the cell becomes competent to unload solutes (Figure 1C–E; Video 1).
In order to acquire functional data on the dimensions of the unloading zone, we conducted flow velocity measurements in individual PSE files using fluorescence recovery after photobleaching (FRAP; Froelich et al., 2011). Phloem flow velocities in the terminal region of the root (basipetal to the unloading zone) are in the range of 25 µm/s (Froelich et al., 2011). In this study, we photobleached CF after it had arrived in the phloem unloading zone of the root. In a tube of constant diameter and impermeable walls, the flow velocity is constant. When the walls become leaky, the flow velocity decreases because of loss of fluid. In the case of the protophloem, the tube becomes leaky when unloading occurs. To define the size of the unloading zone, we measured flow velocity along the phloem files and found that deceleration started at about 300–400 µm behind PSE zero (this dimension varied slightly, even within the two sieve tube files in the same root tip; Figure 2). Our data revealed that the phloem of the root is subdivided into distinct structural and functional domains. Translocation into the main root occurs through the metaphloem. In the elongation zone of the root, the metaphloem overlaps with the mature protophloem file, at which point solutes are transferred laterally from the metaphloem to the protophloem (‘transfer zone’; Figure 2). Movement then occurs through the protophloem towards the root tip (‘protophloem translocation zone’; Figure 2) and subsequently into the terminal functional PSEs (‘protophloem unloading zone’), at which point solutes are distributed laterally into the root tip. Apical to this zone lies the ‘protophloem differentiation zone’ that, as described above, plays no role in phloem unloading (Figure 2).
Arabidopsis roots grow at speeds of about 100–150 µm per hour (Beemster and Baskin, 1998). Consequently, after the sieve-plate pores open, a single PSE in the phloem unloading zone is active in unloading for only about 2–4 hr. Initially, cellular remnants such as the nuclear membrane, ribosomes and tonoplast are degraded and removed from the PSE. Subsequently, the PSEs become physically stretched in the elongation zone and finally move into an area that is active in translocation but inactive in unloading (‘protophloem translocation zone’; Figure 2). The PSEs that progress basipetally from the phloem unloading zone into the protophloem translocation zone must undergo a rapid cellular transformation, particularly in the PD on their lateral walls. However, the structural alterations that control the cessation of unloading are not known.
Callose is deposited in the neck region of plasmodesmata and restricts the size exclusion limit of the pore in response to various stimuli (Luna et al., 2011; Nakashima et al., 2003; Radford et al., 1998). To investigate if callose was involved in occluding PD in the translocation zone, we stained roots of intact seedlings with Sirofluor (Evans et al., 1984; Vatén et al., 2011). This fluorochrome has a strong affinity for callose and demarcates sieve plates and PD in developing cell walls (Stone et al., 1984; Vatén et al., 2011); Figure 3). The two protophloem files could be visualized easily with this stain (Figure 3A and B) and individual lateral PD were identified in the phloem unloading zone (Figure 3C and D). At the junction of the phloem unloading zone with the phloem translocation zone, callose deposits on the PSE wall increased in number (Figure 3E) until significant parts of the PSE walls were covered in callose (Figure 3F). This deposition of callose along the lateral walls of PSEs in the protophloem translocation zone correlated with the reduction in flow we observed in dye-loading experiments and may provide a structural basis for the functional isolation of this zone.
As well as sugars and amino acids, a large number of macromolecules including nucleic acids, and proteins occur in various amounts in the translocation stream (Fukumorita and Chino, 1982; Kehr, 2006; Turnbull and Lopez-Cobollo, 2013). Many of these may access the phloem sap via companion cells (CCs; Paultre et al., 2016). In addition, when new SEs develop in source regions of the plant, remnants of the nucleus, ribosomes, and tonoplast pass into the translocation stream. In both these scenarios, macromolecules must be removed from the terminal PSEs involved in unloading or they would impede the process of phloem unloading (Paultre et al., 2016). While small solutes and GFP (Stokes radius 2.82 nm; Terry et al., 1995) are unloaded from the PSEs into surrounding cells (Imlau et al., 1999), larger proteins unloaded from the PSE enter a unique post-phloem domain (Stadler et al., 2005; Paultre et al., 2016) and cannot cross the interface between the pericycle and endodermis (Paultre et al., 2016). We therefore sought to identify this domain as well as the routes taken by solutes and macromolecules during the unloading process.
In Arabidopsis, the root protophloem file is surrounded by five distinct cell files (Figure 4A and B). One immature metaphloem sieve element (MSE), located towards the center of the root, is in contact with the PSE. Two CCs, one to the left and one to the right, share a common cell wall area with the PSE and MSE, respectively. The complex is capped by two phloem pole pericycle (PPP) cell files (Figure 4A and B) that share cell walls with both the PSE and CCs. Small probes of the dimensions of sugars and amino acids (CFDA, esculin), and small proteins (GFP) move from cell to cell throughout the entire root tip without noticeable barriers (Figure 1; Oparka et al., 1994; Stadler et al., 2005; Wright and Oparka, 1996). However, translocated proteins of CC origin, are restricted to cells immediately adjacent to the mature protophloem file (Stadler et al., 2005; Paultre et al., 2016).
To monitor the fate of different phloem cargos, we compared the phloem-unloading pathway of small fluorescent probes (Knoblauch et al., 2015) with fluorescent fusion proteins (Froelich et al., 2011; Stadler et al., 2005), covering a molecular mass range of 340 Da to more than 100 kDa. Experiments were conducted in situ on living roots. As shown before, small probes were unloaded laterally towards the cortex and subsequently distributed amongst all root cells (Figure 1). Large probes of CC origin, however, were restricted to a specific domain (Figures 4 and 5). We wished to determine if proteins synthesized within SEs, rather than CCs, would show a similar pattern of unloading. Accordingly, we generated a transgenic Arabidopsis line expressing the Sieve Element Occlusion Related (SEOR) protein (Froelich et al., 2011; Pélissier et al., 2008) fused to YFP (SEOR-YFP). This protein is generated exclusively in young sieve elements and remains a structural component after maturation of the cell and its integration into the translocating sieve tube file. We crossed this line with a second reporter line in which HDEL-GFP was targeted to the ER (Figure 1). Both these fluorescent fusions were generated under the sieve-element specific promoters pSEOR; pMtSEO2. The promoters are active in young sieve elements only, but the gene products remain as structural components in mature sieve tubes (Pélissier et al., 2008; Froelich et al., 2011). GFP demarcated clearly the ER of the PSEs (Figures 1 and 4). Unlike HDEL-GFP, the signal from the SEOR-YFP fusion (112 kDa), which is expressed in the cytoplasm of immature SEs, was not restricted to SEs but entered two specific cell files significantly larger than the adjacent PSEs, and did not travel outward beyond this domain (Figure 4C). We used resin embedding of the phloem-unloading zone (Bell et al., 2013) to localize the YFP signal in transverse sections of the root. The YFP signal was confined to the two large PPP cells that abut the PSE files (Figure 4D and E). Time-lapse movies of the phloem unloading zone supported this observation (Figure 4F and G; Video 2). Initially, SEOR-YFP was detected only in the immature PSEs as small aggregates. Subsequently, groups of new PPP cells became fluorescent adjacent to the PSEs, presumably when PD connections opened up between the two cell types (Figure 4F and G). Using a fluorescent marker for CC nuclei (Zhang et al., 2008) we were unable to detect phloem unloading of SEOR-YFP into CCs, only into the PPP (Figure 4H,I). This result was unexpected as most models of phloem unloading assume that the exit of macromolecules occurs through CCs (Patrick, 1997; Thorne, 1985).
It was possible that SEOR-YFP was delivered directly to the PPP from neighboring immature PSEs in the root, rather than by long-distance transport of the protein from the shoot. To test if the latter occurred, we grafted SEOR-YFP scions onto non-transgenic rootstocks. 10 days after grafting, the roots of the grafted plants were indistinguishable from those of the native SEOR-YFP line (Figure 4I,J), demonstrating that SEOR-YFP protein, synthesized in SEs of the shoot, is unloaded from PSEs into the PPP.
To monitor the unloading of large proteins in real time, we used the transgenic lines pAtSUC2-GFP (27 kDa), pAtSUC2-ubiquitin-gfp (36 kDa), pAtSUC2-aequorin-gfp (48 kDa),, and pAtSEOR-AtSEOR-yfp (112 kDa; described above) and conducted fluorescence recovery after photobleaching (FRAP) experiments. Time-lapse movies revealed that unloading of large probes (GFP 27 kDa and above) did not occur at a constant rate, as observed with small solutes. Rather, these proteins were suddenly released into the adjoining PPP cells in distinct pulses (Figure 5, Video 3). We refer to this phenomenon as ‘batch unloading’, as the proteins are delivered in discrete pulses into individual cells. These cells then became highly fluorescent relative to their neighbors. Batch unloading did not occur at the same time into all cells, but independently into individual cells (Figure 5; Video 3). It appears that specific domains exist in the root tip that confer a size-dependent filtration of the macromolecules that arrive in the unloading zone. Small proteins such as GFP were batch unloaded and subsequently moved freely throughout the root (Figure 5A–F). However, large probes such as aequorin-GFP and SEOR-YFP became trapped inside the PPP after batch unloading (Figure 5G–J; Video 4).
To further test the role of the PPP in phloem unloading we inhibited the lateral movement of solutes into this cell layer using an inducible callose synthase system (Vatén et al., 2011). In this approach a modified, PD-specific, callose synthase (icals3m) is induced under estradiol treatment, inhibiting cell-cell movement into the induced cell layer of the root (Vatén et al., 2011). In Arabidopsis there are twelve CalS proteins, with diverse roles in callose production (Chen and Kim, 2009). During the course of this study, we found that CalS8 is expressed specifically in the PPP (Figure 6A,B), allowing us to generate an estradiol-inducible line expressing icals3m under the CalS8 promoter (pCALS8::icals3m). This line was generated specifically to block the connection between PSEs and the PPP. After callose induction we observed a significant arrest in primary root growth relative to roots that were transferred to a non-inducing medium (Figure 6—figure supplement 1). Next we labelled the cotyledons of induced and non-induced pCALS8::icals3m roots with CF at 8 hr and 24 hr after callose induction. Control roots showed the characteristic unloading pattern of wild-type plants (Figure 6C, n = 22) while induced plants showed a severe restriction in phloem unloading. At 8 hr induction, CF was restricted to PSE files with minimal lateral unloading into the root (Figure 6D, n = 20). At 24 hr induction, this effect was even stronger. CF did not enter the unloading zone and remained completely restricted to PSEs distal to the unloading zone (Figure 6E, n = 9). To determine the potential role played by CCs in phloem unloading, we used the sister of apple (sAPL) promoter to drive the production of the mutant CALS3 protein (psAPL::icals3m). We chose this promoter as other CC-specific promoters (e.g. SUC2) were expressed only weakly within the phloem-unloading zone (data not shown). Within the phloem unloading zone the psAPL promoter was expressed strongly in CCs and MSE, but not in the PPP (Figure 6F,G). As above, CF transport in these seedlings was monitored at 8 hr and 24 hr after estradiol induction. At both time points CF was unloaded without restriction, and induced and non-induced seedlings showed an identical pattern of unloading (Figure 6H–J, n = 3, n = 3, n = 2, respectively). Furthermore, we found that callose induction in the psAPL::icals3m seedlings did not significantly affect root growth relative to the uninduced control plants (Figure 6—figure supplement 1). To confirm that callose synthase was overexpressed in the predicted cell layers, roots were fixed, embedded and stained for callose after the CF transport assays. Control roots showed the characteristic callose labelling associated with PSE files (Figure 6K; see also Figure 3). Confocal imaging of pCALS8::icals3m seedlings confirmed that callose was deposited strongly in the PPP relative to other cells surrounding the PSE files (Figure 6L–O), while psAPL::icals3m roots showed strong callose labelling in the CC files, as predicted (Figure 6—figure supplement 6P–S). Collectively, these data provide strong evidence that phloem unloading occurs predominantly via the PPP, not CCs, and suggest that the PD between PSE and PPP are central players in the unloading process.
In order to understand the structural basis for the differential filtration of molecules in the unloading zone, we investigated the ultrastructure of PD at all cell interfaces between the protophloem file and its five interconnecting cells (See Figure 4A and B). We found different PD types in the unloading zone (Figure 7A). PD between the PSE and MP were simple, or occasionally branched (Figure 7B). Between PSEs and CCs, PD displayed the typical single pore on the PSE side and multiple branches towards the CC (Figure 7C), as shown often for connections between SEs and CCs in other tissues (Esau and Thorsch, 1985; Oparka and Turgeon, 1999). The structure of most of the PD between PSE and PPP, however, was quite different and showed an architecture not previously described. In addition to a few simple PD (10%; n = 20), funnel-shaped PD with apertures of up to 300 nm diameter were found at the PSE entrance, tapering towards the PPP entrance (Figure 7D–H).
In order to extract the number of PD available for unloading, we conducted serial block-face scanning electron microscopy followed by 3D reconstruction of the interface between PSEs and the adjacent cells (Denk and Horstmann, 2004; Furuta et al., 2014; Video 5). Using this method, we collected 2100 serial transverse sections that spanned six complete SEs (Figure 8A), including the junction between PSE zero and the adjacent immature PSE (Figure 8B). During reconstruction, we attributed the exact position of each PD to the different wall interfaces shared by the PSE, color-coded to ease identification (Video 5). Using 3D reconstruction, we could image PD on both the outer (Figure 8C) and inner (Figure 8D) walls of the PSE files. A movie showing a ‘fly through’ of the interior of the PSE file is shown in Video 5. PD relative distribution was 45.3%, 40.8%, and 13.9% between PSE-PPPs, PSE-CCs, and PSE-MSE, respectively, when corrected for the total interface between the cell types. Thus, PD leading from the PSE to PPP and CCs are equally abundant while those to the MSE are significantly lower. Next, we acquired data for the total numbers of PD along the entire phloem-unloading zone. We found that the total number averaged 527 ± 58 (n = 4) per protophloem unloading domain (since Arabidopsis roots are diarch, twice the number is available for unloading within the entire root tip). Based on the percentage distributions of PD obtained using serial sectioning, we calculated that in the phloem-unloading zone of the root there are approximately 240 plasmodesmata at the PSE-PPP interface, 215 plasmodesmata at the PSE-CC interface, and 73 plasmodesmata at the PSE-MP interface potentially available for unloading.
Current models of phloem unloading, such as the high-pressure-manifold model (Fisher, 2000; Patrick, 2013) assume high pressure differentials between PSEs and the unloading zone in order to drive fluids and solutes through the PD into the neighboring cells. Recently, direct sieve tube pressure, viscosity, flow velocity, and tube geometry measurements in morning glory showed, however, that the majority of energy provided by the pressure differential between source and sink is consumed to drive flow, and that a high pressure differential in the root system is unlikely to exist, calling the current model into question (Knoblauch et al., 2016). In principle, there are two possibilities of how the phloem sap might escape the PSEs. Bulk flow could move the entire volume of solvent and solutes through the PD into neighboring cells. Alternatively, the solutes could diffuse through PD while the solvent (water) could be removed from PSEs via membrane leakage, potentially via aquaporins (Doering-Saad et al., 2002). Finally, a combination of diffusion and bulk flow could account for the observed transport.
In order to evaluate the feasibility of the different routes, we gathered the necessary parameters to model flow in the phloem-unloading zone (Table 1). We measured the average phloem flow velocity in the protophloem translocation zone to be u = 22.6 ± 5.1 µm/s (n = 11) and the average sieve tube diameter to be: d = 3.6 ± 0.44 µm (n = 11) resulting in a volume flow rate of 230 femtoliter/s (230 µm3/s) to be unloaded from a single PSE file. Literature values on phloem sucrose concentration average about 500 mM (Hall and Baker, 1972; Kallarackal et al., 1989); Mendoza-Cózatl et al., 2008; Winter et al., 1992), hence the total amount of sugar to be unloaded is approximately . The solutes must leave the phloem through the PD, while some of the solvent (water) may be removed via membrane leakage. To facilitate transport, an average of 240 PD are available at the PSE/PPP interface, while the water may leak anywhere along the 350 µm-length of the unloading zone. For simple PD we assumed a cylinder of constant dimensions with a desmotubule of 15 nm diameter and a cytoplasmic sleeve width of 2.8 nm, corresponding to the hydrodynamic radius of GFP. Based on our observations that GFP is batch unloaded, but diffuses relatively freely throughout the root, this molecule appears to be slightly larger than the size exclusion limit imposed by the neck region of funnel PD. For funnel PD we assumed the same dimensions on the PPP side but an average funnel opening of 150 nm (Figure 7) towards the PSE. The location of the desmotubule in funnel PDs impacts the resistance which led us to assume two extremes which are discussed in detail in the appendix. The cell wall thickness (length of the simple and funnel PDs) was taken to be 500 nm.
To elucidate the role of funnel PD in unloading, we first considered the conditions necessary to facilitate unloading by simple PD. We found that the required pressure differential to drive unloading solely by bulk flow through simple PD would be 8.14 MPa (see Appendix for detailed calculations). This pressure differential has neither been measured in SEs, nor can it be considered as feasible. However, assuming funnel PD instead of simple PD, the required pressure would be as low as 0.05–0.2 MPa. Such a relatively low pressure differential would be facile to maintain between plant cells. Because the viscous resistance in wide pores is greatly reduced (Appendix), funnel PD are much more efficient in unloading compared to simple PD.
To estimate the solute concentration difference that would be required to account for unloading by pure diffusion, we assumed a length of the unloading zone of 350 µm. For the geometrical and physiological values listed above, a concentration difference of 276 mM would be required for diffusive unloading through funnel PDs. In this case, however, a second route for unloading of the solvent (water) would be necessary. The permeability of the plasma membrane varies significantly between cell types and values of 10−14 to 10−12 m/s/Pa have been reported (Kramer and Boyer, 1995a). Considering the parameters outlined above, a pressure of 0.059 to 5.9 MPa would be needed to remove the solvent from the PSE. Thus, at a high membrane permeability removal of the solvent over the membrane would be feasible.
Our analysis of unloading kinetics leads us to conclude that the distinctive funnel shape of the SE-to-PPP PD is crucial to enabling efficient unloading. In this case, convective phloem unloading, i.e. a combination of diffusion and bulk flow, is feasible, and that relatively moderate pressure and concentration differentials are necessary to drive transport. In our calculations (Appendix) we have assumed that PD connecting PSEs with CCs and MSEs play little or no direct role in the unloading process, consistent with our results using the icals3m system (Figure 6). In this scenario, the unique architecture of funnel PD could still accommodate the calculated unloading rates (Appendix)
Based on our structural and physiological data bulk flow is likely to be the dominant mechanism of unloading. Bulk flow would require a pressure differential of about 0.05–0.2 MPa, equal to an osmotic potential difference of about 20–80 mM. Solute unloading by diffusion alone, on the other hand, would require a concentration difference of 276 mM. While bulk flow predominates, the above mechanisms are not exclusive, and both may contribute to unloading. Bulk flow requires a pressure differential generated by osmosis which would lead to a concentration difference between PSEs and PPPs. This would induce diffusion, in parallel with bulk flow, even if its contribution is lower. The term for the combination of diffusion and bulk flow is convection, which leads us to propose a new model of Convective Phloem Unloading.
We have identified different domains within the root tip that provide a size dependent filtration of molecules. Cross sections of Arabidopsis roots (Figure 9A) show the typical diarch structure. The pentagonal architecture of the cell complex surrounding the PSEs and the different PD connecting the cell types is depicted in Figure 9B and C. Our results indicate that the PSE-PPP interface is the principal route for all solutes to be unloaded. Diffusion will, however, lead to a relatively quick redistribution of the solutes through simple PD connecting the cells within the post-phloem unloading zone, along a solute concentration gradient (Figure 9D,E). This gradient will also direct solutes to areas of highest demand, while higher consumption will lead to steeper gradients.
Unloading of small proteins such as GFP (27 kDa) occurs by batch unloading through funnel PD into PPPs. Once unloaded, GFP enters the post unloading zone and is evenly distributed throughout the root tip (Stadler et al., 2005). Large proteins of the size of aequorin-GFP (48 kDa) or SEOR-YFP (112 kDa), however, remain trapped within PPPs and cannot enter the post-phloem pathway to the meristem. Clearly, the different types of PD between PSEs and neighboring cells have major functional impacts on molecular flow in the root tip.
Many proteins have been found in extracted phloem sap, most showing no obvious function in long-distance signaling (Turnbull and Lopez-Cobollo, 2013; Batailler et al., 2012; Paultre, et al., 2016). Thus, the question arises as to how many of these molecules are components of the phloem sap by default rather than by design (Paultre, et al., 2016). To address this point we need to consider where such macromolecules originate. Many systemic macromolecules are thought to arise in CCs and pass into SEs through the pore-PD that connect SEs with CCs (Fisher et al., 1992; Lucas et al., 1996; Sjolund, 1997; Oparka and Turgeon, 1999). There appears to be selectivity to this movement, large proteins (above 70 kDa) are restricted from entering the translocation stream while those below this cutoff are not (Paultre, et al., 2016). Interestingly, in our current study SEOR-YFP (112 kDa), which is translated in immature SEs rather than CCs, was able to unload into the PPP, suggesting that the size exclusion limit that regulates macromolecular exchange between PPP and SE in sinks may be larger than between SEs and CCs in source phloem tissues. Another potential origin of macromolecules in phloem sap is young developing sieve elements. As shown in Figure 1, PSEs start to conduct when their nucleus degenerates. At this point the cytoplasmic content in the PSEs is still very dense. After integration into the mature PSE file, there is no other route left for degraded structures but phloem unloading, most likely into the PPP. Unlike the xylem pole pericycle (XPP), which is involved in lateral-root formation (Parizot et al., 2012) and apoplastic xylem loading (Takano et al., 2002), a specific role for the PPP has not been proposed, although its transcriptome shares similarity to the underlying protophloem (Parizot et al., 2012). We have now identified the ‘post-phloem domain’ previously described by Stadler et al. (2005) and Paultre, et al., 2016 as the PPP. It appears that while the XPP may be involved in apoplastic xylem loading, the PPP is intimately involved in symplastic phloem unloading, and may function as a repository for degraded PSE contents as well as systemic macromolecules. Although PD connect all cell types surrounding the PSEs, our data suggest that exit via the PPP is the major route of unloading in Arabidopsis roots. In mature source tissues, as well as in secondary sieve elements produced by the cambium, macromolecules will arrive in the phloem-unloading zone. This was demonstrated when we grafted SEOR-YFP scions onto wild-type rootstocks (Figure 4I,J). Without continuous macromolecule unloading, PSEs would fill up rapidly and phloem unloading would be impaired (Paultre, et al., 2016). PPPs are the major recipients of systemic macromolecules and it might be expected that they are specialized for the degradation of proteins and RNAs, an interesting area for future research. PPPs may also protect the downstream cells from receiving unwanted, potentially interfering mRNA molecules.
How does the root distinguish between true signaling molecules and those moving by default in the phloem? One possibility is that unloaded macromolecules destined to traffic further than the PPP may have specific sequences that interact with the PD that connect the PPP to the endodermis and beyond. Phloem-mobile viruses are able to move beyond the PPP suggesting that they can break the ‘barrier’ normally imposed by the PPP (Valentine et al., 2004). A clear future challenge will be to identify systemic macromolecules that can move beyond the PPP, and to analyze these macromolecules for specific motifs that, like transcription factors (Xu et al., 2011), allow them to interact with and modify PD.
Arabidopsis thaliana were grown in microROCs (Advanced Science Tools LLC. Pullman, WA) as described elsewhere (Froelich et al., 2011). Soil was kept saturated by placing the microROCs in 1–2 inches of water daily for 5 min creating a soil water saturation of >80%. Chamber conditions were 16 hr photoperiods with 100–200 μE m−2 s−1 at 18°C to 22°C.
Experiments were conducted on a Leica SP8 confocal microscope equipped with a supercontinuum laser (470–630 nm) and a pulsed 405 nm diode laser. 5(6)-Carboxyfluorescein diacetate (CFDA) stock solution, 5 mg/ml in acetone, was diluted 1:10 (v/v) in ddH2O and loaded into the phloem through application to cut cotyledons at 7–10 days post-germination as described earlier (Knoblauch et al., 2015). Sieve elements in the most apical root protophloem loaded with CFDA were photo-bleached with 480 nm, 488 nm, and 496 nm lasers concurrently at maximum power and 6.5x zoom with a 20x lens. Bleaching started apically and moved toward the hypocotyl until the entirety of CFDA in the sieve tube in the field of view was bleached. Subsequent recording of the refilling occurred immediately following photo-bleaching at 200 Hz and 0.75x digital zoom by excitation with the 488 nm laser line at 15% of continuous power. Emission was collected between 505 and 545 nm.
Seeds of Arabidopsis wt and the transgenic line pAtSEOR-AtSEOR-yfp were sterilized with 8% bleach and 1% TWEEN-20. After 5x wash in distilled water seeds were plated on MS with 1.2% agar and 0.2% sucrose, pH 5.7. After 48 hr stratification in the dark at 4°C seeds were oriented vertically at 23°C with 18 hr photoperiod. After 5–7 days, seedlings were grafted following the hypocotyl-grafting procedure of Turnbull et al. (2002) consisting of a transverse cut and butt alignment with silicon collars. The seedlings were cut transversely in the upper region of the hypocotyl with ultrafine microknives (Interfocus, n°10315–12). Scions were grafted onto wild type stocks using a short silicon collar for support on MS agar plates. The grafts were left to grow under LD with the plates still oriented vertically until new lateral roots of the stocks were fully established (~10 days). The grafts were imaged between 5 and 10 dag.
Ten days old seedlings of transgenic Arabidopsis thaliana, expressing AtSEO::SEO-YFP, were embedded in LR White following the method described by Bell et al. (2011). Tissue samples were fixed overnight at 4°C in a solution of 4% formaldehyde, 1% glutaraldehyde, 50 mM 1,4-piperazinediethanesulofonic acid (PIPES) and 1 mM CaCl2. The samples were then washed in buffer (50 mM PIPES, 1 mM CaCl2) three times for 10 min before dehydration in a graded ethanol series (50%, 70%, 2 × 90%). The tissue samples were then infiltrated at 4°C in medium grade LR White at 1:1, 1:2, 1:3, 1:9 ratios of 90% ethanol:resin for 45 min each before two 60 min changes in 100% LR White. The final embedding step was done at ambient temperature. The samples were then polymerised in gelatin capsules (TAAB) at 50°C for 24 hr.
Esculin and Sirofluor were excited with a pulsed 405 nm diode laser. Emission was collected with a hybrid detector between 420 and 480 nm or 420 and 600 nm respectively. GFP and esculin were imaged with sequential scan for fluorescence emission separation. For GFP, excitation was 488 nm and emission detection with a hybrid detector at 495–535 nm. Esculin excitation was 405 nm and emission collection with a hybrid detector at 420–470 nm. GFP and YFP were sequentially scanned for separation of emission fluorescence. GFP excitation was 484 nm at 50% continuous power. Emission collected with a hybrid detector at 489–505 nm. YFP excitation was 514 nm at 15% continuous power. Emission collected at 519–564 nm with a hybrid detector.
Samples were chemically fixed in the microwave with 2% glutaraldehye, 2% paraformaldehyde and 2% DMSO in 0.1 M cacodylate buffer. Microwave fixation time was 2 min, followed by a 2 min break and a final 2 min microwave period at maximum temperature of 25°C. Samples were then washed 5x for 10 min in 0.1 M cacodylate buffer. Postfixation took place in 1% osmium tetroxide (0.1 M cacodylate) overnight at 4°C. Samples were rinsed 3x for 10 min. in 0.1 M cacodylate buffer. Dehydration was carried out with methanol in the microwave for one minute followed by 5 min at room temperature for each step in 5% increments from 5% to 100% methanol. 100% methanol was replaced twice and the samples were then transferred to 100% propylene oxide with three exchanges. Samples were embedded in Spurr resin at the following propylene oxide:Spurr ratios: 2:1 for 2 hr, 1:1 for 2 hr, 1:2 for 2 hr, and 100% Spurr overnight 3X. The resin was cured at 65°C for 24 hr. Samples were sectioned and stained with 1% uranyl acetate for 6 min followed by a ddH2O rinse and stained in Reynold’s lead for 6 min. Following Reynold’s lead samples were quickly rinsed in 0.1 N NaOH and ddH2O. Samples were imaged with an FEI T20 at 200kv.
Material for SBFSM was prepared as described by Furuta et al. (2014). Wild-type Columbia roots were fixed with 2.5% glutaraldehyde, 2% formaldehyde in 0.1 M Na-Cacodylate buffer (pH 7.4) supplemented with 2 mM CaCl2 for 2–3 hr at room temperature and embedded in Durcupan ACM resin (Fluka, Sigma-Aldrich). After standard dehydration steps, samples were embedded in silicone holders filled with 100% Durcupan and infiltrated for at least 2 hr before polymerization at 60°C. The roots were trimmed to the desired starting point from the tip using an EM Ultracut UC6i ultramicrotome (Leica Mikrosysteme GmbH). Images were acquired with a FEG-SEM Quanta 250 (FEI, Hillsboro, OR), using a backscattered electron detector (Gatan Inc.). The block faces were sectioned at 40 nm increments. The images were initially processed and segmented using Microscopy Image Browser, a self-developed program written under Matlab environment (Mathworks Inc.) and available at http://mib.helsinki.fi/. A SBFSEM dataset of 2100 images was subsequently segmented and visualized in 3D using Amira 6.0 software (FEI Corp). The outer structure cell walls were discriminated in a hybrid method of interpolated interactive masking and grey-level segmentation. The individual plasmodesmata were identified and segmented sequentially by a manual interactive method and color-coded into groups that related to the particular cell wall interface along the protophloem files.
The promoters of CALS8 (At3g14570) and sAPL (At3g12730) were amplified using the primers CALS8-AscI-F AAGGCGCGCCCGGCAACATGAAATACGGGA and CALS8-XhoI-R ACAGCTCGAGGTTTTGGGAGAAAATCAATCAGAA, and SAPL-AscI-F AAGGCGCGCCAGCTAATAAGAAAGGGAGATCTCTG and SAPL-XhoI-R ACAGCTCGAGTTAACTAACAAAGTACTAAATGCCGA, respectively) and cloned into P4P1RpGEMt containing the estrogen receptor XVE (Zuo et al. 2000). Using the MultiSite Gateway system (Invitrogen), the inducible promoters were combined with the icals3m construct and the nopaline synthase terminator in destination vector pBm43GW. Arabidopsis Col-0 plants were dipped with the different constructs and positive transformants were selected using Basta (pCALS8) and Hygromycin (psAPL). Lines with single insertions were selected in T2 and homozygous plants were obtained in T3.
pCALS8::icals3m and psAPL::icals3m seeds were stratified for 2 days at 4°C in the dark and were sown on media containing 0.5 MS and 1.2% plant agar, pH 5.8. At 4 days post sowing, seedlings were induced in plates containing the same media plus 5 µM Beta estradiol or the corresponding amount of DMSO (mock). Plates were scanned at times 0, 3, 8 and 24 hr after the induction. Roots were measured using ImageJ and seedlings were subsequently used for CFDA loading. 1 mM CFDA in a solution of 72% acetone in water was mixed with 0.5% Adigor solution (Syngenta) to 10:1 proportions. 0.2 µL of this solution was applied to the adaxial side of the leaves, separating the leaves from the media to avoid undesired dye spreading. Approximately 10 roots were used for each treatment. 20 min after loading, roots were mounted in propidium iodide and imaged by confocal laser scanning microscopy (Zeiss LSM700).
Seedlings used for phloem unloading analysis subsequently underwent callose detection experiments through Sirofluor staining and callose immunolocalization allowing three independent experiments on the same biological material. Callose accumulation on the CFDA positive roots used for unloading experiments was confirmed first by Sirofluor staining. Whole seedlings were incubated 2 hr at room temperature in a 25 µg/mL (42.66 µM) solution of Sirofluor fluorochrome (Biosupplies, AU) in a 50 mM K2HPO4 buffer under vacuum (60 mPa). Roots were then dissected and mounted in a 1:1 solution of AF1 antifadent (Citifluor) and 67 mM K2HPO4 buffer, and subsequently imaged by confocal laser scanning microscopy (Zeiss LSM700). Following Sirofluor staining, roots were fixed for subsequent immunolocalization on cross sections in a fixative solution containing 4% formaldehyde (freshly prepared from paraformaldehyde powder, Sigma) and 0.5% glutaraldehyde (Sigma) in a 0.1 M phosphate buffer pH7. Fixation, dehydration and resin infiltration steps were done in a micro-wave using a PELCO BioWave Pro (Ted Pella, Redding, CA). Fixation was achieved at 150 W, under vacuum (20 Hg) and for 6 min (2’. Followed by a 2’break and a final 2’ microwave). Samples were left in the fixative overnight at 4°C and then washed 3 times 3 min (20 Hg, MW 150 W 1’ on – 1’ break – 1’ on). Roots were then aligned and embedded in 1% low-melting Agarose (Calbiochem) and processed through increasing dehydration steps (25%, 50%, 70%, 90%, 96%, 3 × 100% Ethanol, vacuum 20 Hg, MW 150 W 1’ on followed by 5’ break). Resin infiltration (LR White medium grade, Agar scientific) was achieved through increasing resin concentrations: 33% resin in ethanol 100%, 66% resin in ethanol 100%, and 3 times 100% resin (20 Hg, MW 200 W 2’ on – 2’ break – 2’ on). Polymerization was conducted overnight at 60°C. Semi-thin sections (0.5 µm) were taken with a Leica EM UC7 ultramicrotome. Callose immunolocalization on semi-thin sections was micro-wave assisted and performed as follows: blocking step (BSA 2% in PBS, 1 mL per slide; primary antibody (anti (1→3)-β-glucan (Biosupplies, AU), 1/200 in BSA 2% in PBS, 500 µL per slide, MW 170 W 2’ on – 2’ break – 2’ on); three washes in BSA 2% in PBS; secondary antibody (Alexa Fluor 488 anti-mouse IGG ThermoFisher Scientific, A-11017, 1/400 in BSA 2% in PBS (MW 170 W, 1’), 500 µL per slide; three washes in BSA 2% in PBS. Slides were finally mounted in a 1:1 solution of AF1 antifadent (Citifluor) with PBS, containing calcofluor as a cell wall counterstain, and imaged by confocal laser scanning microscopy (Zeiss LSM700).
The process of phloem unloading is mediated by plasmodesmata (PD) pores that link phloem sieve elements (SE) and neigboring phloem pole pericycle cells (PPP). Transport of solutes out of the phloem occurs by a combination of advection (bulk flow) and diffusion through PD (Appendix 1—figure 1a). The contribution from each mechanisms depends on the PD properties as well as the SE ability to maintain pressure and concentration gradients relative to the PPP.
In this supplement, we consider different unloading scenarios and estimate the pressure differentials required to drive phloem unloading by bulk flow in each case: First, we consider unloading by bulk flow through simple PD. Next, the pressure differential required to drive unloading through funnel PD is evaluated. Finally, we discuss the role of diffusive transport through through both simple and funnel PD, and aquaporin-mediated unloading of water. From these results, we deduce the most probable mode of phloem unloading able to account for our experimental observations.
Phloem unloading from SE to PPP is facilitated by two types of PD (Appendix 1—figure 1b). Simple PDs are annular pores of inner radius , slit width , and length . Funnel PDs taper from the SE towards the PPP from a wide approximately circular opening of radius . The parameters used in the following calculations are given in Appendix 1—table 1.
To ascertain the feasibility of different unloading scenarios we first compute the total amount of sugar unloaded from a single phloem file. Using measured values of the flow speed and SE radius yields a flow rate of
The sugar current, and the amount to be unloaded, is thus
where is the sap sugar concentration.
Assuming that all liquid must leave the SE via simple PD, the flow rate per PD is , while the maximum flow speed through each PD active in unloading is . Here, is the total number of PD and is the conductive area of one PD with inner radius and outer radius . The characteristic PD Reynolds number is , with viscosity and density . In the most general case, flow resistance through the PD has two contributions at low Reynolds numbers: internal viscous resistance to flow in the gap between the plasma membrane and the desmotubule - and - external viscous resistance associated with effects near the pore entrance. This leads to an expression for the pressure drop with two terms
Both pressure drops depend on details of the PD geometry, the nature of the desmotubule, and on the shape of the entrance region (Hasimoto, 1958; Blake, 1978; Kramer and Boyer, 1995a; Jensen et al., 2012). However, order-of-magnitude estimates can be obtained by assuming laminar Stokes flow inside the PD, and that the entrance region is free from physical obstructions.
Consider the case of unloading by bulk flow through simple PD. The pressure drops in Equation (3) can be estimated by assuming flow throught a straight slit of width and length . In this case, the two pressures are (Hasimoto, 1958; Bruus, 2008).
This pressure difference is relatively large, and most likely practically unfeasible, although we emphasize that this is an order-of-magnitude estimate. It is, however, probable that any objects (e.g. large proteins), which are present in the PD, will lead to an increase in .
We proceed to consider unloading through funnel PD, assuming are of this type. Since funnel PD posses a relatively wide opening facing the SE, we expect that the internal resistance of these is smaller than that computed for the straight narrow pores (Equation 4). There are two possible scenarios depending on the location of the desmotubule. First, we neglect the influence of the desmotubule, and consider flow in an open tapering cylinder. In the lubrication approximation, the hydraulic resistance per unit length of a cylinder of radius is . The pressure drop across a cylinder tapering linearly in radius from to is thus approximately
Second, we take into account the presence of the desmotubule and consider flow in the annular space between two concentric cylinders of inner radius and outer radius . The hydraulic resistance per unit length is (Pozrikidis, 2011)
The pressure drop for and is
The external entrance resistance is given by the expression in Equation (4)
In summary, a total pressure drop of
is required to drive unloading through funnel PD depending on the detailed location of the desmotubule. A cell-to-cell pressure differential of this magnitude does not appear unreasonable. Finally, we note that the pressure required for bulk unloading via funnel PD is at least 40 times smaller than that found for simple PD (Equation 4).
We proceed by considering the contribution to unloading from diffusion. A sugar molecule faces three obstacles when diffusing from SE to PPP (Bret-Harte and Silk, 1994). First, it must traverse the SE cell cytoplasm, i.e. cover a distance of a few microns. Second, it has to pass through the narrow constriction of the plasmodesmata. Finally, it must diffuse throughout the PPP cytoplasm.
We begin by considering the situation in the absence of PDs, i.e. assuming that the entire cell wall area is perfectly conductive. In that case, the unloading rate is
where is the SE-to-PPP sugar concentration difference and is the sum of the SE and PPP diameter and is the surface area of the SE. We assume here that , where is the length of the unloading zone. The parameters is the diffusion coefficient for sucrose. Equating the diffusive current in Equation (10) to the total incoming sap (Equation (2)) gives a required concentration gradient of
To ascertain the importance of the PDs in diffusive unloading, we consider opposite of the previous situation, i.e. the case where the sole resistance to diffusion is provided by the PDs. The diffusion rate through a short isolated pore of radius on a flat surface is governed by Stefans diameter law (Berg, 1993)
We proceed to consider the influence of the pore length. The diffusion rate through pores of radius and length is given by
Finally, we investigate the combined effects of cytoplasm and PD resistance to diffusion. We follow Berg (Berg, 1993) and write the total diffusive current as
where is the combined diffusion resistance. The three terms in the denominator are (the diffusion resistance of the PD opening) (the diffusion resistance of the PD interior) and (the diffusion resistance of the cell cytoplasm), corresponding to Equations (10) and (12). The total unloading rate is thus
where is the PD density. This leads to an expression for of the form
For funnel PD, the cell and pore opening resistances are similar to those compute above
Along the pore length, we take into account the varying cross-sectional area
This leads to a current of
Finally, this leads to an expression for the concentration difference required to drive diffusive unloading
If unloading occurs exclusively by diffusion through PDs, the water contained in the phloem sap must leave the SE by traversing the cell membrane. To estimate the pressure differential required to facilitate this, we assume that the volume flow rate move uniformly across the membrane of area (surface area of SEs). Taking as the permeability of the membrane leads to
It follows that the pressure is given by
The magnitude of the permeability varies considerably between cell types, with values in the range from to reported by Kramer and Boyer (see Table 3.1 of Kramer and Boyer ). This leads to pressures in the range
Soluble and filamentous proteins in Arabidopsis sieve elementsPlant, Cell & Environment 35:1258–1273.https://doi.org/10.1111/j.1365-3040.2012.02487.x
Correlative imaging of fluorescent proteins in resin-embedded plant materialPlant Physiology 161:1595–1603.https://doi.org/10.1104/pp.112.212365
Random Walks in BiologyPrinceton University Press.
Theoretical MicrofluidicsOxford: Oxford University Press.
Phloem transport: a review of mechanisms and controlsJournal of Experimental Botany 64:4839–4850.https://doi.org/10.1093/jxb/ert302
Use of aphid stylectomy and RT-PCR for the detection of transporter mRNAs in sieve elementsJournal of Experimental Botany 53:631–637.https://doi.org/10.1093/jexbot/53.369.631
Aspects of sieve element ontogeny and structure in smilax rotundifoliaBotanical Gazette 128:138–144.https://doi.org/10.1086/336389
Turnover of soluble proteins in the wheat sieve tubePlant Physiology 100:1433–1441.https://doi.org/10.1104/pp.100.3.1433
Biochemistry and Molecular Biology of Plants730–785, Long distance transport, Biochemistry and Molecular Biology of Plants.
Sugar, amino acid and inorganic contents in rice phloem sapPlant and Cell Physiology 23:273–283.https://doi.org/10.1093/oxfordjournals.pcp.a076347
The angiosperm phloem sieve tube system: a role in mediating traits important to modern agricultureJournal of Experimental Botany 65:1799–1816.https://doi.org/10.1093/jxb/ert417
On the flow of a viscous fluid past a thin screen at small reynolds numbersJournal of the Physical Society of Japan 13:633–639.https://doi.org/10.1143/JPSJ.13.633
Electrical wiring and Long-Distance plant communicationTrends in Plant Science 21:376–387.https://doi.org/10.1016/j.tplants.2016.01.016
Phloem sap proteins: their identities and potential roles in the interaction between plants and phloem-feeding insectsJournal of Experimental Botany 57:767–774.https://doi.org/10.1093/jxb/erj087
Münch, morphology, microfluidics - our structural problem with the phloemPlant, Cell & Environment 33:1439–1452.https://doi.org/10.1111/j.1365-3040.2010.02177.x
Multispectral phloem-mobile probes: properties and applicationsPlant Physiology 167:1211–1220.https://doi.org/10.1104/pp.114.255414
Water Relations of Plants and SoilsSan Diego: Academic Press.
eLSCompanion cells, eLS, John Wiley & Sons, Ltd, 10.1038/npg.els.0002087.
Callose deposition: a multifaceted plant defense responseMolecular Plant-Microbe Interactions 24:183–193.https://doi.org/10.1094/MPMI-07-10-0149
Identification of high levels of phytochelatins, glutathione and cadmium in the phloem sap of Brassica napus. A role for thiol-peptides in the long-distance transport of cadmium and the effect of cadmium on iron translocationThe Plant Journal 54:249–259.https://doi.org/10.1111/j.1365-313X.2008.03410.x
Die Stoffbewegung in Der PflanzeJena, Germany: Fischer.
Symplastic communication between primary and developing lateral roots of Arabidopsis thalianaJournal of Experimental Botany 46:187–197.https://doi.org/10.1093/jxb/46.2.187
In silico analyses of pericycle cell populations reinforce their relation with associated vasculature in ArabidopsisPhilosophical Transactions of the Royal Society B: Biological Sciences 367:1479–1488.https://doi.org/10.1098/rstb.2011.0227
PHLOEM UNLOADING: sieve element unloading and Post-Sieve element transportAnnual Review of Plant Physiology and Plant Molecular Biology 48:191–222.https://doi.org/10.1146/annurev.arplant.48.1.191
Introduction to Theoretical and Computational Fluid DynamicsOxford University Press.
Aphid stylectomy reveals an osmotic step between sieve tube and cortical cells in barley rootsJournal of Experimental Botany 47:1519–1524.https://doi.org/10.1093/jxb/47.10.1519
GFP tagging of sieve element occlusion (SEO) proteins results in green fluorescent forisomesPlant and Cell Physiology 49:1699–1710.https://doi.org/10.1093/pcp/pcn141
Molecular characterisation of recombinant green fluorescent protein by fluorescence correlation microscopyBiochemical and Biophysical Research Communications 217:21–27.https://doi.org/10.1006/bbrc.1995.2740
Phloem unloading of C and N assimilates in developing seedsAnnual Review of Plant Physiology 36:317–343.https://doi.org/10.1146/annurev.pp.36.060185.001533
Phloem transport: cellular pathways and molecular traffickingAnnual Review of Plant Biology 60:207–221.https://doi.org/10.1146/annurev.arplant.043008.092045
Sieve element and companion cell: the story of the comatose patient and the hyperactive nurseFunctional Plant Biology 27:477–487.https://doi.org/10.1071/PP99172
Callose biosynthesis regulates symplastic trafficking during root developmentDevelopmental Cell 21:1144–1155.https://doi.org/10.1016/j.devcel.2011.10.006
The fluorescent probe HPTS as a phloem-mobile, symplastic tracer: an evaluation using confocal laser scanning microscopyJournal of Experimental Botany 47:439–445.https://doi.org/10.1093/jxb/47.3.439
Christian S HardtkeReviewing Editor; University of Lausanne, Switzerland
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Phloem unloading in Arabidopsis roots is convective and regulated by the phloem-pole pericycle" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by Christian Hardtke as the Reviewing and Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Ulrich Z Hammes (Reviewer #1); David E Salt (Reviewer #2).
The reviewers have discussed the reviews with one another and the Senior Editor has drafted this decision to help you prepare a revised submission.
As you can see from the reviews below, all three reviewers appreciate your study and believe that your careful mapping of the root phloem unloading process and the identification of a new, specialized plasmodesmata type in this process are an exciting advance in the field. They also have a number of suggestions for further improvement of the manuscript, which we would like you to address one by one to the best of your capacity in a revised version.
The manuscript "Phloem unloading in Arabidopsis roots is convective and regulated by the phloem-pole pericycle" by Ross-Elliot and colleagues presents a breakthrough in our current understanding of phloem unloading- an agronomically important process. This manuscript addresses (and answers) a major question that arose from a previous manuscript of this (now extended) consortium which was published in eLife last year: how is phloem unloading achieved if the pressure gradient between the phloem and the surrounding tissue in sinks is low and the high-pressure manifold model that was believed to be working in roots does not work. Additionally, the authors found that phloem unloading occurs in a specialized cell type which is different from the tissue that was – so far – thought to be the major player in this context. An interesting consequence from this paper is also that it shows that possibly/likely many of the macromolecules, a very "hot" topic that is currently discussed controversially, found in the phloem represent "garbage" rather than signaling molecules. Last not least: I spent a while going through the literature: This is the first study that provides data on the anatomy and (sub)cellular organization of the unloading zone which led to the discovery of a new type of plasmodesmata (PD).
Taken together, I believe this paper contains many exciting findings that can be published pretty much without major changes.
In this manuscript the authors provide a wealth of elegant experimental data to help define the process of phloem unloading in roots of Arabidopsis thaliana.
I found the conclusions developed by the authors compelling and well supported by the data.
By combining confocal imaging of both GFP-tagged protophloem sieve elements (PSEs) and the fluorescent phloem mobile dye Esculin (which does not accumulate in vacuoles after unloading which is an advantage over carboxyfluorescein) the authors convincingly define the first transport competent PSE in a given file of PSE cells in the root. The authors use fluorescence recovery after photobleaching (FRAP) to measure phloem velocity and observe a reduction in this velocity in the phloem unloading zone. Using the first PSE identified as PSE-zero as a marker for the end of the unloading zone and flow rate reduction as a precise measure of phloem unloading the authors are able to precisely define the extent of the phloem unloading domain. A significant achievement. Strong evidence is presented using fluorescently tagged proteins and confocal imaging that proteins are unloaded in 'batches' and not as a steady stream. Furthermore, using cell-type specific expression of callose synthase the authors establish that unloading from the PSEs occurs into the phloem pole pericycle cells (PPP) and not into the companion cells (CC) as previously proposed. To further investigate the involvement of plasmodesmata (PD) in unloading from the PSEs into the PPP an in depth investigation using serial block-face scanning EM with 3D reconstruction is undertaken. This produces an incredible reconstruction of a file of PSEs and detailed structural information about PD that reveals a new type of 'funnel' PD that connect the PSEs and PPP. This is an exciting discovery. Modelling of phloem unloading is then used to predict that these new funnel PD provide a much more efficient pathway for phloem unloading, which presents a nice potential rational for the existence of such PD.
Overall I found this manuscript to be a fascinating read that describes a study that is rich with significant new insights, observations and testable hypothesis.
In this manuscript, by combination of non-invasive imaging, 3D electron microscopy and mathematical modeling, the authors show that phloem unloading in Arabidopsis root is driven by bulk flow and diffusion ("convection") through plasmodesmata (PD). Imaging supports the idea that phloem pole pericycle (PPP) is a predominant pathway of phloem unloading, rather than companion cells (CCs), supported by the artificial PD closure by callose deposition. Additionally, they also show that most PD interconnecting the protophloem sieve elements (PSEs) and neighboring PPP are shaped like a funnel, which could enable unloading with lower pressure differential between PSE and PPP. Furthermore, they indicate that, during unloading process, large macromolecules are trapped by PPP and restricted in further movement, whereas small solutes are unloaded through PPP and more freely move throughout the root tip, suggesting that PPP may regulate unloading process depending on molecular size, and act as a repository for unloaded macromolecules. Thus, in this study, authors reveal that PPP and its funnel PD act as central players in phloem unloading, which will adds to the regulatory mechanism of phloem unloading and PD function depending on their diverse structures.
This study is mostly descriptive, though the discovery of the phloem unloading route and funnel PDs are important advances. However the mechanism of unloading regulation (e.g. molecular size-dependent regulation, funnel PD formation) is still largely unknown.
I am a little concerned that the data on protein movement is only from 3 proteins, one of which, GFP, is freely mobile, and another, SEOR, is involved in making aggregates, and the third is a foreign protein (aequeorin). There are plenty of native plant fusions available, how would they behave?
The batch unloading is not clear to me, and only shown for GFP in video 3.
Results section, I don't understand: Once unloading occurred, the flow velocity decreased significantly due to loss of liquid volume during unloading.
Figure 6: Callose staining (I assume in green) appears weaker in the psAPL line (and not all the way around the cells) compared to the pCalS8 line- could this be a reason why the latter blocks transport and not the former? – Indeed, callose could be accumulated at CCs' cell wall by psAPL:iCalS3m, but is the function of PD between SE and CC really defective? To check PD function, authors could examine whether CC-expressed GFP probe can move to surrounding cells.
The references are often very old and recent refs on protein and mRNA movement are not mentioned- the discussion should be more balanced on this and avoid strong opinions on things not investigated in the manuscript, eg: "PPPs may also protect the downstream cells from receiving unwanted,potentially interfering mRNA molecules"https://doi.org/10.7554/eLife.24125.021
I am a little concerned that the data on protein movement is only from 3 proteins, one of which, GFP, is freely mobile, and another, SEOR, is involved in making aggregates, and the third is a foreign protein (aequeorin). There are plenty of native plant fusions available, how would they behave?
In the current study, we used three fusion proteins to cover a wide Mr size range. In the recent study by Paultre et al., 2016, it was shown that a host of native fusion proteins leave the terminal protophloem sieve elements and enter a unique ‘post-phloem domain’, identified in the current manuscript as the phloem-pole pericycle. These authors also used fusion proteins with organelle- targeting sequences and so our observations of protein exit from the protophloem are not just restricted to the three fusion proteins reported here.
The batch unloading is not clear to me, and only shown for GFP in Video 3.
As explained in the text (subsection “Batch Unloading of Proteins”), we refer to batch unloading when unloading is not continuous, but occurs spontaneously (in batches) into individual cells. This is the case for all probes of the size of GFP and larger. We provide examples in the form of a figure (Figure 5) and two videos (Video 3 and 4) which represent the extremes we tested (GFP, 27 kDa and SEOR-YFP 112 kDa). We have numerous videos with other probes (Aequorin GFP, Ubiquitin-GFP etc.), but we do not see a benefit in overloading the figure or adding more videos that basically show the same phenomenon. However, we have mentioned the additional fusion proteins that we tested in the text.
Results section, I don't understand: Once unloading occurred, the flow velocity decreased significantly due to loss of liquid volume during unloading.
We have rephrased the sentence to make it more easily understood.
Figure 6: Callose staining (I assume in green) appears weaker in the psAPL line (and not all the way around the cells) compared to the pCalS8 line- could this be a reason why the latter blocks transport and not the former?
Fluorescence immunolabelling is suitable for cellular localization of but not usually for quantification as this depends on various factors such as antigen exposure.
Indeed, callose could be accumulated at CCs' cell wall by psAPL:iCalS3m, but is the function of PD between SE and CC really defective? To check PD function, authors could examine whether CC-expressed GFP probe can move to surrounding cells.
We agree that it would be helpful to have an additional control. We used the suggested line for our batch unloading studies (pAtSUC2-GFP), but unfortunately it does not provide the results. GFP is constantly supplied to the sieve tube from the companion cells in the metaphloem and in the secondary phloem. Therefore we cannot discriminate between GFP that potentially leaks from protophloem CCs from GFP that is delivered to the root tip from metaphloem CCs. In order to make this experiment as the reviewers suggest, we should engineer a line which is expressing GFP in a CC-specific manner in the protophloem of the unloading zone only. Currently we don’t have such line and establishing such line will go beyond the two-month revision time frame.
The references are often very old and recent refs on protein and mRNA movement are not mentioned- the discussion should be more balanced on this and avoid strong opinions on things not investigated in the manuscript, eg: "PPPs may also protect the downstream cells from receiving unwanted,potentially interfering mRNA molecules"
We believe that our references are adequate and that up-to-date references on protein trafficking are indeed included. The suggestion that the PPP may filter unwanted mRNA molecules is simply a hypothesis that could be tested by others in future work.https://doi.org/10.7554/eLife.24125.022
- Michael Knoblauch
- Karl J Oparka
- Kaare H Jensen
- Kaare H Jensen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
KJO acknowledges the financial support of the BBSRC. We thank Ilya Belevich for preparing material for SBFSEM and Kirsten Knox and Andrea Paterlini for advice on tracer experiments. We thank Pawel Roszak for providing T2 seeds of psAPL::icals3m lines. We thank the Francheschi Microscopy and Imaging center for technical support. This work was supported by National Science Foundation grant IOS-1146500 (MK).
- Christian S Hardtke, University of Lausanne, Switzerland
© 2017, Ross-Elliott et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Downloads (link to download the article as PDF)
Download citations (links to download the citations from this article in formats compatible with various reference manager tools)
Open citations (links to open the citations from this article in various online reference manager services)
Quantitative analysis of plant and animal morphogenesis requires accurate segmentation of individual cells in volumetric images of growing organs. In the last years, deep learning has provided robust automated algorithms that approach human performance, with applications to bio-image analysis now starting to emerge. Here, we present PlantSeg, a pipeline for volumetric segmentation of plant tissues into cells. PlantSeg employs a convolutional neural network to predict cell boundaries and graph partitioning to segment cells based on the neural network predictions. PlantSeg was trained on 1xed and live plant organs imaged with confocal and light sheet microscopes. PlantSeg delivers accurate results and generalizes well across different tissues, scales, acquisition settings even on non plant samples. We present results of PlantSeg applications in diverse developmental contexts. PlantSeg is free and open-source, with both a command line and a user-friendly graphical interface (https://github.com/hci-unihd/plant-seg).
Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity.