Charting the native architecture of Chlamydomonas thylakoid membranes with single-molecule precision
Abstract
Thylakoid membranes scaffold an assortment of large protein complexes that work together to harness the energy of light. It has been a longstanding challenge to visualize how the intricate thylakoid network organizes these protein complexes to finely tune the photosynthetic reactions. Previously, we used in situ cryo-electron tomography to reveal the native architecture of thylakoid membranes (Engel et al., 2015). Here, we leverage technical advances to resolve the individual protein complexes within these membranes. Combined with a new method to visualize membrane surface topology, we map the molecular landscapes of thylakoid membranes inside green algae cells. Our tomograms provide insights into the molecular forces that drive thylakoid stacking and reveal that photosystems I and II are strictly segregated at the borders between appressed and non-appressed membrane domains. This new approach to charting thylakoid topology lays the foundation for dissecting photosynthetic regulation at the level of single protein complexes within the cell.
Introduction
Membranes orchestrate cellular life. In addition to compartmentalizing the cell into organelles, membranes can organize their embedded proteins into specialized domains, concentrating molecular partners together to drive biological processes (Malinsky et al., 2013; Case et al., 2019). Due to the labile and often transient nature of these membrane domains, it remains a challenge to study how individual protein complexes are organized within them.
The question of membrane domain organization is especially pertinent to thylakoids, sheet-like membrane-bound compartments that produce oxygen while converting light energy into biochemical energy, thereby sustaining most of the life on Earth. These light-dependent photosynthetic reactions are driven by the concerted actions of four large protein complexes within the thylakoid membrane. Photosystem II (PSII), cytochrome b6f (cytb6f), and photosystem I (PSI) form an electron transport chain that produces NADPH and pumps protons into the thylakoid lumen. ATP synthase then uses this electrochemical gradient across the membrane to generate ATP. In most photosynthetic eukaryotes, the thylakoid membranes are subdivided into appressed regions that face adjacent membranes within thylakoid stacks (called grana in higher plants) and non-appressed regions that freely face the stroma (Austin and Staehelin, 2011; Nevo et al., 2012; Pribil et al., 2014; Engel et al., 2015). In both plants and algae, PSII and PSI appear to be segregated to the appressed and non-appressed membranes, respectively (Goodenough and Staehelin, 1971; Wollman et al., 1980; Pribil et al., 2014; Flori et al., 2017). This lateral heterogeneity is believed to coordinate photosynthesis by concentrating different reactions within specialized membrane domains, while the redistribution of light-harvesting complex II (LHCII) antennas between these domains may enable adaptation to changing environmental conditions (Minagawa and Tokutsu, 2015; Nawrocki et al., 2016).
Much of what is known about the molecular organization of thylakoids comes from freeze-fracture electron microscopy (and the related deep etch technique), which provides views of membrane-embedded protein complexes within the cell (Heuser, 2011). Combined with biochemical fractionation (Andersson and Anderson, 1980; Albertsson, 2001), these membrane panoramas helped describe the lateral heterogeneity of thylakoids (Goodenough and Staehelin, 1971; Staehelin, 1976; Olive et al., 1979; Wollman et al., 1980; Staehelin and Arntzen, 1983). However, the platinum replicas produced by this technique have limited resolution and only provide access to random fracture planes through the membranes. More recently, atomic force microscopy (AFM) has been used to map the macromolecular organization of thylakoids (Sznee et al., 2011; Johnson et al., 2014; Wood et al., 2018). AFM can very precisely measure topology and thus, can distinguish between each of the major photosynthetic complexes within a hydrated membrane. However, the membranes must first be removed from the cell, and only one membrane surface can be visualized at a time. As an alternative, cryo-electron tomography (cryo-ET) has been shown to resolve PSII complexes within hydrated thylakoids (Daum et al., 2010; Kouřil et al., 2011; Levitan et al., 2019). Although multiple overlapping membranes can be imaged with this technique, the thylakoids in these studies were isolated from the chloroplast in order to produce sufficiently thin samples.
Dissecting the interrelationship between membrane domains and thylakoid architecture requires a molecular view of intact thylakoid networks within native cells. In pursuit of this goal, we combined cryo-focused ion beam milling (Marko et al., 2007; Schaffer et al., 2017) with cryo-ET (Asano et al., 2016) to image thylakoid membranes within vitreously-frozen Chlamydomonas reinhardtii cells. Several years ago, we demonstrated how this approach can capture the undisturbed membrane architecture of this green alga, revealing an elaborate system of stacked thylakoids (Engel et al., 2015). However, due to the lower resolution of the CCD cameras used at that time, this investigation was limited to a description of membrane architecture, without visualizing the protein complexes embedded within these membranes. Here, we leverage advances in direct detector cameras and the contrast-enhancing Volta phase plate (Danev et al., 2014) to resolve the thylakoid-embedded complexes, enabling us to describe the molecular organization of thylakoids in situ, within their native cellular context.
Results and discussion
Our tomograms show numerous densities corresponding to ribosomes and photosynthetic complexes decorating the appressed and non-appressed surfaces of the thylakoid network (Figure 1A–B, Figure 1—figure supplement 1, Video 1). Close visual inspection revealed the unmistakable shapes of ATP synthase and PSII protruding into the stroma and lumen, respectively (Figure 1C–D). In order to map these decorating densities into the native thylakoid architecture, we developed a visualization approach called a ‘membranogram’, where densities from the tomogram are projected onto the surface of a segmented membrane (Figure 1E). The result is a topological view of the membrane surface that resembles AFM data. However, unlike AFM, membranograms can display the topology of both sides of each membrane within the cellular volume. By dynamically growing and shrinking the segmentations, the membranograms allow us to track how densities change as they extend from the membrane surface and compare these densities to mapped in molecular models of different complexes (Video 2). This enabled the manual assignment of thylakoid-associated complexes by considering whether the densities projected into the stroma or thylakoid lumen (Figure 2A) and then comparing the densities to known structures of each complex (Figure 2—figure supplement 1). We used manually picked positions to generate subtomogram averages of PSII and ATP synthase, structurally confirming that membranograms can correctly identify these complexes (Figure 2—figure supplement 2).
We used membranograms to examine the molecular organization within appressed and non-appressed domains of the thylakoid network (Figure 2). The complexes in appressed membranes (M2 and M3 in Figure 2B–D) could be reliably assigned. The stromal surfaces of these membranes showed no large densities, whereas two clearly distinguished classes of densities were seen on luminal surfaces: large PSII dimers and smaller cytb6f dimers that extend ~4 nm and ~3 nm into the lumen, respectively. Non-appressed membranes (M1 in Figure 2B–D) were more complicated to analyze. We first assigned thylakoid-associated ribosomes and ATP synthases with high certainty by growing the segmentation ~10 nm away from the membrane. Assignment of PSI was more difficult: we marked single densities that protruded ~3 nm into the stroma and were not directly under a ribosome or ATP synthase. On the luminal side, almost no large densities were observed that could correspond to PSII. Although we did see frequent small dimers that resembled cytb6f, they were more difficult to identify due to increased background signal on the luminal surfaces of non-appressed membranes. We generated membrane models from the assigned particles (Figure 2D) to analyze how the complexes are arranged within the plane of the membrane. In total, we quantified 84 membrane regions from four tomograms (Table 1) and found clear evidence of lateral heterogeneity: PSII was almost exclusively found in the appressed regions, whereas PSI, ATP synthase, and ribosomes were restricted to the non-appressed regions. Cytb6f was observed with almost equal abundance in both regions, with the caveat that we had lower confidence in our assignment of this complex in non-appressed membranes. Poly-ribosome chains were clearly resolved decorating the non-appressed thylakoids (Figures 1A and 2C–D, Figure 2—figure supplement 3). In contrast, PSI, PSII, cytb6f, and ATP synthase were distributed relatively evenly along their respective membrane regions (nearest-neighbor distances plotted in Figure 2—figure supplement 4).
The concentrations we measured for different photosynthetic complexes were similar to previous biochemical estimates (Allred and Staehelin, 1986; Vallon et al., 1991; Albertsson, 2001) as well as AFM measurements (Johnson et al., 2014) and counts of ‘membrane-embedded particles’ from freeze-fracture (Wollman et al., 1980; Staehelin, 2003), indicating that direct comparisons can be made between cryo-ET data and these earlier studies. In particular, the concentration of PSII dimers that we observed in appressed regions (1122 particles/μm2) was in agreement with previous freeze-fracture measurements from Chlamydomonas (1250 particles/μm2) (Wollman et al., 1980). We observed very few PSII complexes in the non-appressed membranes (Figure 2—figure supplement 1B), but it remains to be tested whether different conditions, such as damaging high light intensities, may increase the abundance of PSII outside the appressed domains. Our observation that cytb6f is distributed relatively evenly between appressed and non-appressed membranes is also consistent with previous reports, although this ratio varies widely in the literature and may be modulated by physiological conditions (Vallon et al., 1991; Staehelin, 2003; Kirchhoff et al., 2017). The total number of densities that we counted on non-appressed membranes (4907 particles/μm2) was similar to the number of particles previously observed by freeze-fracture in Chlamydomonas (4500 particles/μm2) (Wollman et al., 1980). However, we are not aware of any previous direct measurement of PSI complexes to which we could compare. Spectroscopic estimates for the stoichiometry of PSII/PSI monomers in thylakoids vary widely from 0.8 to 2.1 (Danielsson et al., 2004; Fan et al., 2007); by cryo-ET, we measured a PSII/PSI monomer ratio of 2.1, but we note that PSI was not identified with high confidence (Table 1). The concentration of ATP synthase that we measured on non-appressed membranes (1652 particles/μm2) was very similar to previous cryo-ET analysis of isolated spinach thylakoids (1770 particles/μm2) (Daum et al., 2010) but more than twice the concentration initially proposed from freeze-fracture (710 particles/μm2) (Miller and Staehelin, 1976; Staehelin, 2003).
At the boundary between the appressed and non-appressed membranes of higher plants (i.e., the grana and stroma lamellae), it is widely believed that there is a specialized domain called the ‘grana margin’ where PSII and PSI intermix (Anderson, 1989; Wollenberger et al., 1995; Kitmitto et al., 1997; Albertsson, 2001). Although green algae do not have traditional grana, we were curious if such margin regions could be observed in Chlamydomonas. Using membranograms, we visualized the molecular organization of native thylakoids that transitioned between appressed and non-appressed regions (Figure 3, Figure 3—figure supplement 1). Strikingly, we observed a sharp boundary between PSII and PSI that exactly matched the division in thylakoid architecture. The photosystems did not intermix, and we saw no clear difference in particle abundance near the boundary. We therefore find no evidence for ‘grana margins’ in Chlamydomonas cells grown under these conditions (mixotrophic growth in moderate light, see Materials and methods).
What drives the strict lateral heterogeneity that we observe between appressed and non-appressed domains? PSI is presumably excluded from appressed membranes because its ~3 nm stromal density is too bulky to fit into the ~3 nm space between stacked thylakoids (Daum et al., 2010; Kirchhoff et al., 2011; Engel et al., 2015). Conversely, PSII and its associated LHCII antennas may induce thylakoid stacking, a causal relationship that would precisely limit PSII to appressed membranes. Several studies have observed semi-crystalline arrays of C2S2-type (Boekema et al., 2000; Daum et al., 2010) or C2S2M2-type (Kouřil et al., 2012) PSII-LHCII supercomplexes in thylakoids isolated from higher plants, and it has been proposed that the overlap of LHCII or PSII between membranes mediates thylakoid stacking (McDonnel and Staehelin, 1980; Boekema et al., 2000; Standfuss et al., 2005; Daum et al., 2010; Albanese et al., 2017; Albanese et al., 2020). Although we observed randomly oriented PSII complexes instead of ordered arrays, we nonetheless looked for evidence of supercomplex interactions across native thylakoid stacks (Figure 4). We first created membranogram overlays of adjacent membranes spanning either the thylakoid lumen or stromal gap (Figure 4B). Then we generated membrane models by using the positions and rotational orientations of PSII luminal densities seen in the membranograms to place structures of C2S2M2L2-type PSII-LHCII supercomplexes (Burton-Smith et al., 2019; Shen et al., 2019; Sheng et al., 2019), the largest supercomplexes that have been isolated from Chlamydomonas (Figure 4C–D). Note that because LHCII barely protrudes from the membrane surface (McDonnel and Staehelin, 1980; Standfuss et al., 2005; Johnson et al., 2014) and thus is not well resolved in membranograms, we relied solely on the orientations of the PSII core complexes to place the supercomplex models. While the majority of C2S2M2L2-type supercomplexes fit within the plane of the membrane, we observed a ~3% in-plane overlap between the models (Figure 4C), indicating that some PSII may form smaller supercomplexes under the moderate light conditions that we examined (~90 µmol photons m−2s−1). It should be noted that the previously characterized C2S2M2L2-type supercomplexes were isolated from cells grown under lower light (20–50 µmol photons m−2s−1), which should favor larger supercomplex assemblies. Nevertheless, we observed that there is ample space within the appressed regions of Chlamydomonas thylakoids to accommodate large PSII-LHCII supercomplexes. Mapping in C2S2M2-type supercomplexes, a slightly smaller arrangement that has been purified from higher plants (Su et al., 2017; van Bezouwen et al., 2017), resulted in almost no in-plane overlap between the models. Supercomplex models of different sizes occupied 47.3 ± 6.0% (C2S2M2L2), 40.7 ± 5.2% (C2S2M2), and 29.1 ± 3.7% (C2S2) of the membrane surface area (Figure 4E), with cytb6f occupying an additional 5.8 ± 1.6%. Thylakoids are ~70% protein (Kirchhoff et al., 2002), suggesting that other complexes such as extra LHCII antennas may occupy up to 20% of the surface area. This spacing should also allow room for rapid diffusion of plastoquinone between PSII and cytb6f within appressed membranes.
The membranogram overlays showed very little overlap of PSII between adjacent membranes. In particular, PSII complexes were almost completely interdigitated across the thylakoid lumen, with only ~4% overlap of EM density (Figure 4B and F, Figure 4—figure supplement 1). PSII interdigitation has previously been predicted from spinach thylakoids with a contracted ~4.5 nm lumen, which sterically could only permit interdigitation (Daum et al., 2010). However, the thylakoids in our light-adapted Chlamydomonas cells have a ~9 nm lumen, which is enough space to permit face-to-face interactions of PSII across the lumen, yet the PSII complexes remain interdigitated. It has been suggested that interdigitation could facilitate diffusion of the soluble electron carrier plastocyanin through the lumen, while also enabling the lumen to contract and block diffusion in response to changing light conditions (Kirchhoff et al., 2011). To investigate interactions across the stromal gap, we measured the intermembrane overlap from the membrane models of C2S2M2L2-type supercomplexes and compared it to simulated data where the same number of supercomplexes were randomly positioned (Figure 4D and F, Figure 4—figure supplement 2). We found that the overlap of PSII-to-PSII and LHCII-to-LHCII were both statistically indistinguishable from random. It remains to be tested whether randomly distributed interactions covering ~8% of the PSII surface area and ~13% of the S2M2L2 LHCII surface area (disregarding the ~3% in-plane overlap) can mediate the tight adhesion between stacked thylakoid membranes.
In situ cryo-ET has enabled us to visualize the native architecture of thylakoids with molecular resolution, revealing several striking findings that raise questions about how thylakoids partition their proteins and form appressed stacks. We observe sharp boundaries between PSII and PSI that are tightly coupled to membrane architecture, with no evidence of intermixing at transitions between appressed and non-appressed membrane domains. This argues against the presence of ‘grana margins’. However, more definitive conclusions will require investigation of the true grana stacks found in higher plants. Semi-crystalline arrays of PSII have been reported in thylakoids isolated from spinach, pea, and diatoms (Staehelin, 1976; Boekema et al., 2000; Daum et al., 2010; Sznee et al., 2011; Goral et al., 2012; Kouřil et al., 2012; Levitan et al., 2019). However, randomly arranged PSII has been observed even more frequently (Goodenough and Staehelin, 1971; Staehelin, 1976; Olive et al., 1979; Wollman et al., 1980; Daum et al., 2010; Kouřil et al., 2011; Goral et al., 2012; Johnson et al., 2014; Wood et al., 2018; Levitan et al., 2019), and thus the physiological relevance of PSII arrays remains debated (Tsvetkova et al., 1995; Ruban and Johnson, 2015; Tietz et al., 2015; Charuvi et al., 2016). In our first look at PSII organization within native membranes, we found that PSII complexes are randomly oriented within the appressed thylakoids of Chlamydomonas, with ample spacing to accommodate C2S2M2-type and C2S2M2L2-type supercomplexes. This observation supports the physiological relevance of these larger supercomplexes, which have been primarily characterized in vitro (Caffarri et al., 2009; Su et al., 2017; van Bezouwen et al., 2017; Burton-Smith et al., 2019; Shen et al., 2019; Sheng et al., 2019). However, it remains to be proven whether semi-crystalline PSII arrays form inside the native cellular environment or whether they are a consequence of isolating thylakoid membranes. Unlike earlier techniques, our combination of in situ cryo-ET with membranograms has enabled us to visualize the organization of PSII across multiple adjacent stacked thylakoids within native cells. We determined that PSII supercomplexes randomly overlap between appressed membranes, challenging the idea that thylakoid stacking is mediated by specific static interactions between LHCII or PSII proteins across the stromal gap. Beyond these observations, our study establishes the methodological foundation to explore the conservation of thylakoid molecular architecture with other ecologically important phototrophs and dissect the mechanisms that adapt this architecture to changing environmental conditions.
Materials and methods
Cell culture
Request a detailed protocolWe used the Chlamydomonas reinhardtii strains mat3-4 (CC-3994) (Umen, 2001) and wild-type CC-125, provided by the Chlamydomonas Resource Center, University of Minnesota, MN. The mat3-4 strain has smaller cells that vitrify better by plunge freezing. Comparative 77K fluorescence measurements showed that mat3-4 and wild-type cells had similar spectra under the conditions that the cells were grown and frozen onto EM grids (Figure 1—figure supplement 2). Thus, our description of the arrangement of photosynthetic complexes within mat3-4 thylakoids is likely comparable to wild-type thylakoids.
For cryo-ET, mat3-4 cells were grown to mid-log phase (1000–2000 cells/μL) in Tris-acetate-phosphate (TAP) medium under constant light conditions (~90 µmol photons m−2s−1) and bubbling with normal air. For 77K measurements, both mat3-4 and wild-type strains were grown under these same conditions, and measurements were made both immediately and after allowing cultures to sit in low light for 30 min.
77K measurements
Request a detailed protocolFresh colonies of CC-125 and CC-3944 were picked from agar plates and suspended in sterile air-bubbled flasks containing 50 mL TAP medium. The cultures were grown at room temperature under 90 µmol photons m−2s−1 light until reaching 1,000 cells/μL. Samples containing 15 mL culture suspension were transferred into duplicate 15 mL Falcon tubes. From one of these tubes, five 1 mL aliquots were transferred to Eppendorf tubes and frozen immediately in liquid nitrogen. The second Falcon tube was kept under dim light of ~10 μmol photons m−2s−1 for 30 min, with gentle mixing by inversion, and five 1 mL samples were then frozen as above. Fluorescence emission spectra were acquired as previously described (Kanazawa et al., 2014) under liquid nitrogen, using a custom-made apparatus. Excitation was from a 440 nm diode laser, the output of which was sent through one branch of a bifurcated optical fiber (diameter of 1 mm) to the sample. The emission light was collected through the second fiber and measured with a spectrometer (Ocean Optics HR200 + ER). The average 77K profile for each condition was composed of three independent biological replicates, each with five technical replicates (15 measurements total per condition).
Vitrification and Cryo-FIB milling
Request a detailed protocolUsing a Vitrobot Mark 4 (FEI, Thermo Fisher Scientific), 4 µL of mat3-4 cell culture was blotted onto R2/1 carbon-coated 200-mesh copper EM grids (Quantifoil Micro Tools) and plunge frozen in a liquid ethane/propane mixture. Grids were stored in liquid nitrogen until used for FIB milling. Cryo-FIB milling was performed following a previously reported procedure (Schaffer et al., 2015; Schaffer et al., 2017), using a Quanta dual-beam FIB/SEM instrument (FEI, Thermo Fisher Scientific), equipped with a Quorum PP3010 preparation chamber. Grids were clipped into Autogrid support rings modified with a cut-out on one side (FEI, Thermo Fisher Scientific). In the preparation chamber, the frozen grids were sputtered with a fine metallic platinum layer to make the sample conductive. Once loaded onto the dual-beam microscope’s cryo-stage, the grids were coated with a thicker layer of organometallic platinum using a gas injection system to protect the sample surface. After coating, the cells were milled with a gallium ion beam to produce ~100 nm-thick lamellas (Figure 1 = 60–90 nm, Figure 1—figure supplement 1A = 100–140 nm, Figure 1—figure supplement 1C = 60–80 nm, Figure 1—figure supplement 1E = 110–140 nm). During transfer out of the microscope, the finished lamellas were coated with another fine layer of metallic platinum to prevent detrimental charging effects during Volta phase plate cryo-ET imaging, as previously described (Mahamid et al., 2016; Schaffer et al., 2017).
Cryo-ET
Request a detailed protocolAfter FIB milling, grids were transferred into a 300 kV Titan Krios microscope (FEI, Thermo Fisher Scientific), equipped with a Volta phase plate (Danev et al., 2014), a post-column energy filter (Quantum, Gatan), and a direct detector camera (K2 summit, Gatan). Prior to tilt-series acquisition, the phase plate was conditioned to about 0.5π phase shift (Danev et al., 2017). Using SerialEM software (Mastronarde, 2005), bidirectional tilt-series (separated at 0°) were acquired with 2° steps between −60° and +60°. Individual tilts were recorded in movie mode with 12 frames per second, at an object pixel size of 3.42 Å and a target defocus of −0.5 µm. The total accumulated dose for the tilt-series was kept below ~100 e-/Å2. Each tomogram was acquired from a separate cell and thus is both a biological and technical replicate.
Tomogram reconstruction
Request a detailed protocolFrame alignment was performed with K2Align (https://github.com/dtegunov/k2align). Using IMOD software (Kremer et al., 1996), tilt-series were aligned with patch tracking, and bin4 reconstructions (13.68 Å pixel size) were created by weighted back projection. Of the 13 Volta phase plate tomograms acquired, four tomograms (Figure 1, Figure 1—figure supplement 1) were selected for analysis of photosynthetic complexes based on good IMOD tilt-series alignment scores and visual confirmation of well-resolved complexes at the thylakoid membranes.
Membrane segmentation
Request a detailed protocolSegmentation of chloroplast membranes was performed in Amira software (FEI, Thermo Fisher Scientific), aided by automated membrane detection from the TomoSegMemTV package (Martinez-Sanchez et al., 2014). Bin4 tomograms were processed with TomoSegMemTV to generate correlation volumes with high pixel intensity corresponding to membrane positions. The original bin4 tomograms and correlation volumes were imported into Amira, and the correlation volumes were segmented by 3D threshold-based selection, producing one-voxel-wide segmentations at the centers of the membranes. Using the 3D lasso selection tool, these segmentations were subdivided into appressed and non-appressed membrane regions (Figure 1B, Figure 1—figure supplement 1B,D and F). To generate membranograms visualizing complexes directly protruding from the membrane surface (PSII, PSI, cytb6f), the segmentations were grown by two voxels in all directions to produce a five-voxel-wide segmentation with a surface that matched the surface of the membrane. To generate membranograms visualizing complexes ~10 nm above the membrane surface (ATP synthase, membrane-bound ribosomes), segmentations were grown by 10 voxels in all directions. To produce smooth 3D surfaces, the segmented voxels were transformed into a polygonal mesh with the ‘generate surface’ command, decimated to 10% triangle density with the ‘remesh surface’ command, and smoothed with the ‘smooth surface’ command (50 iterations, 0.4 lambda). These surfaces were exported in the OBJ 3D model format.
Membranograms and particle picking
Request a detailed protocolBin4 tomograms and corresponding membrane segmentations (OBJ models) were loaded into Membranorama software (https://github.com/dtegunov/membranorama; copy archived at https://github.com/elifesciences-publications/membranorama). This software projects tomographic density onto the surface of a 3D membrane segmentation to create a membranogram. Segmentations that had been grown by two voxels in Amira had surfaces that intersected densities protruding directly from the membrane surface (PSII, PSI, cytb6f). Segmentations that had been grown by 10 voxels in Amira had surfaces that intersected densities that were ~10 nm above the membrane surface (ATP synthase, membrane-bound ribosomes). Importantly, the Membranorama software can dynamically grow and shrink 3D membrane segmentations in real time, enabling the user to interactively track how densities appear at different distances from the membrane surface (see Video 2). The software also allows to-scale 3D models of each molecular complex (PDB: 6IJJ for PSI, 6KAD for PSII, 1Q90 for cytb6f, 6FKF for ATP synthase, 5MMM for ribosome) (Stroebel et al., 2003; Bieri et al., 2017; Hahn et al., 2018; Sheng et al., 2019; Su et al., 2019) to be mapped onto the membrane and compared to the tomographic densities. This is accomplished interactively by clicking the surface of the 3D membrane segmentation and using the mouse wheel to rotate each particle in the plane of the membrane. Using these features, we manually assigned membrane-associated densities to different classes of macromolecular complexes based on their positions relative to the membrane and their characteristic structural features (Figure 2—figure supplement 1). For the luminal side of the membrane, we exclusively used segmentations that had been grown by two voxels. PSII was assigned to large dimeric densities projecting ~4 nm from the membrane surface, and cytb6f was assigned to small dimeric densities projecting ~3 nm from the surface. For the stromal side of the membrane, we started with segmentations that had been grown by 10 voxels, assigning large round densities with ~25 nm diameters to ribosomes and smaller round densities with ~10 nm diameters to the F1 subunit of ATP synthase. The stator of ATP synthase was often observed as a small density adjacent to the larger F1 density (see Figure 2—figure supplement 1B). After assigning the positions of these two complexes, we next loaded the corresponding segmentation that had been grown by two voxels and assigned PSI to small round densities projecting ~3 nm from the surface that were not positioned directly under a ribosome or ATP synthase. This order of particle picking prevented misassignment of PSI to the stalk of ATP synthase or the translocon structures that attach ribosomes to thylakoid membranes. The clarity of membrane-associated densities varied between different membranes within the same tomogram, likely due to effects of the tomographic missing wedge on different membrane curvatures and orientations, as well as local differences in the quality of tilt-series alignment. Therefore, only larger complexes (ribosomes, ATP synthase, PSII) were assigned for membranes with lower-clarity densities. Of the 51 non-appressed membranes quantified in this study, all complexes were assigned in 28 membranes, all complexes except for cytb6f were assigned in two membranes, only ATP synthase and ribosomes were assigned in four membranes, and only ribosomes were assigned in 17 membranes. All complexes were assigned in the 33 appressed membranes quantified in this study. Within appressed membranes, PSII and cytb6f were assigned with high confidence. Within non-appressed membranes, ribosomes and ATP synthase were assigned with high confidence, whereas PSI and cytb6f were assigned with lower confidence (denoted by an asterisk in Table 1).
Subtomogram averaging
Request a detailed protocolWe used subtomogram averaging as a structural confirmation of the membranogram-picked positions for PSII and ATP synthase (Figure 2—figure supplement 2). Manually assigned positions and orientations were exported from Membranorama and used as starting parameters for real space subtomogram alignment in PyTom software (Hrabe et al., 2012). No classification was performed, and all membranogram-picked subvolumes were included in the averages (396 PSII particles, 639 ATP synthase particles).
Analysis of protein complex organization
Nearest-neighbor distances within the plane of a membrane
Request a detailed protocolTo measure nearest-neighbor distances between molecular complexes within the thylakoids (Figure 2—figure supplement 4), segmentations of essentially flat membrane regions were exported from Amira (FEI, Thermo Fisher Scientific) as MRC volumes. Coordinates of the particles (PSII, cytb6f, ATP synthase, ribosomes) assigned on each membrane region were exported from Membranorama software. Each membrane region was projected with its corresponding particles onto a flat plane to generate a 2D surface. Nearest-neighbor distances between the particles were then measured using Matlab scripts calculating the shortest path between objects. Clustered poly-ribosomes exclude large regions of the stromal surface (see Figure 2C–D, Figure 2—figure supplement 3), which can cause misleading nearest-neighbor measurements for ATP synthase. To avoid this, the membranes were cropped to exclude regions containing poly-ribosome clusters before measuring ATP synthase distances.
Overlap between adjacent membranes using EM densities
Request a detailed protocolTo calculate overlap between EM densities from two adjacent membranes (Figure 4F), 11 membrane pairs separated by the stromal gap and six membrane pairs separated by the thylakoid lumen were segmented in Amira and imported into Membranorama software. Using the tools in Membranorama, surfaces from each membrane were selected and overlaid along a vector orthogonal to both membrane surfaces, ensuring a geometrically accurate superposition of the membrane densities. Images of the overlaid membranograms were then analyzed in Fiji software (Schindelin et al., 2012) as described in Figure 4—figure supplement 1. Thresholding and cropping the membranograms were necessary to decrease noise and avoid edge effects, respectively.
Overlap between adjacent membranes using membrane models
Request a detailed protocolLHCII light-harvesting antennas are poorly visualized by cryo-ET because they are almost entirely embedded within the thylakoid membrane. Nevertheless, we incorporated hypothetical PSII-associated LHCII complexes into our analysis by using the PSII core positions and orientations visualized in membranograms to generate membrane models containing C2S2M2L2-type PSII-LHCII supercomplexes. First, 3D segmentations of appressed membrane regions were exported from Amira (FEI, Thermo Fisher Scientific). In the Membranorama software, we manually aligned a correctly scaled 3D model of the C2S2M2L2-type PSII-LHCII supercomplex (PDB: 6KAD) (Sheng et al., 2019) with the EM densities observed for each PSII core particle. Coordinates and orientations of all the particles were exported from Membranorama. Using Matlab scripts, the PSII-LHCII supercomplex structure was filtered to 25 Å resolution and then mapped at the assigned positions and orientations into the 3D volume of the membrane segmentation, generating a 3D model of a one-voxel-thick membrane with embedded PSII-LHCII supercomplexes. To compare these experimentally determined models with simulated models containing randomly distributed supercomplexes, the same number of PSII-LHCII supercomplexes were placed one-by-one at random positions in the membrane (randomly selected voxels of the one-voxel-thick 3D membrane segmentation) and each time rotated by a random in-plane angle before placing the next supercomplex particle. Whenever a newly placed particle overlapped with a preexisting particle, rotation of the new particle was first attempted to avoid overlap, and if this failed, the particle was moved to a new random position (see Figure 4—figure supplement 2 for a 2D schematic representation of this 3D procedure). 100 random models were generated for each appressed membrane region.
To measure the relative membrane area covered by the placed supercomplex structures (Figure 4E), the one-voxel-thick membrane segmentation was masked where it was intersected by the 3D structure volumes (PSII core alone, or PSII-LHCII supercomplexes of increasing size: C2S2-type, C2S2M2-type, and C2S2M2L2-type). This masked ‘occupied area’ was then divided by the total number of voxels in the membrane segmentation to yield a percentage. To calculate the overlap of ‘occupied area’ between adjacent membranes (Figure 4F), the C2S2M2L2-type supercomplex structure was subdivided into separate PSII (C2) and LHCII (S2M2L2) regions using Chimera software (Goddard et al., 2007). Next, volumes of the complexes in one membrane were extended along a vector orthogonal to both membrane surfaces until they also intersected the adjacent membrane. These extended complex densities were then used to mask ‘occupied area’ on both membranes. Similar to the overlap calculation for EM density (detailed in Figure 4—figure supplement 1), the percent of overlap was calculated as the number of membrane voxels that were masked by complexes from both membranes (equivalent to white in Figure 4—figure supplement 1) divided by the sum of all masked voxels on the membrane (equivalent to white + green + magenta in Figure 4—figure supplement 1). ‘White’ overlap voxels masked by complexes from both membranes were only counted once in this calculation. The schematic models in Figure 4C–D provide a 2D representation of this 3D analysis.
Data availability
The four tomograms analyzed in this study have been deposited in the Electron Microscopy Data Bank (EMD-10780-10783). The membranorama software is freely available at: https://github.com/dtegunov/membranorama (copy archived at https://github.com/elifesciences-publications/membranorama).
-
Electron Microscopy Data BankID EMD-10780. In Situ Cryo-Electron Tomogram of the Chlamydomonas Chloroplast.
-
Electron Microscopy Data BankID EMD-10781. In Situ Cryo-Electron Tomogram of the Chlamydomonas Chloroplast.
-
Electron Microscopy Data BankID EMD-10782. In Situ Cryo-Electron Tomogram of the Chlamydomonas Chloroplast.
-
Electron Microscopy Data BankID EMD-10783. In Situ Cryo-Electron Tomogram of the Chlamydomonas Chloroplast.
References
-
A quantitative model of the domain structure of the photosynthetic membraneTrends in Plant Science 6:349–354.https://doi.org/10.1016/s1360-1385(01)02021-0
-
Spatial organization of the cytochrome b6-f complex within chloroplast thylakoid membranesBiochimica Et Biophysica Acta (BBA) - Bioenergetics 849:94–103.https://doi.org/10.1016/0005-2728(86)90100-3
-
The grana margins of plant thylakoid membranesPhysiologia Plantarum 76:243–248.https://doi.org/10.1111/j.1399-3054.1989.tb05640.x
-
Lateral heterogeneity in the distribution of chlorophyll-protein complexes of the thylakoid membranes of spinach chloroplastsBiochimica et Biophysica Acta (BBA) - Bioenergetics 593:427–440.https://doi.org/10.1016/0005-2728(80)90078-X
-
In situ Cryo-Electron tomography: a Post-Reductionist approach to structural biologyJournal of Molecular Biology 428:332–343.https://doi.org/10.1016/j.jmb.2015.09.030
-
Arrangement of photosystem II supercomplexes in crystalline macrodomains within the thylakoid membrane of green plant chloroplastsJournal of Molecular Biology 301:1123–1133.https://doi.org/10.1006/jmbi.2000.4037
-
Structural determination of the large photosystem II-light-harvesting complex II supercomplex of Chlamydomonas reinhardtii using nonionic amphipolJournal of Biological Chemistry 294:15003–15013.https://doi.org/10.1074/jbc.RA119.009341
-
Functional architecture of higher plant photosystem II supercomplexesThe EMBO Journal 28:3052–3063.https://doi.org/10.1038/emboj.2009.232
-
Regulation of transmembrane signaling by phase separationAnnual Review of Biophysics 48:465–494.https://doi.org/10.1146/annurev-biophys-052118-115534
-
Quantification of photosystem I and II in different parts of the thylakoid membrane from spinachBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1608:53–61.https://doi.org/10.1016/j.bbabio.2003.10.005
-
The stoichiometry of the two photosystems in higher plants revisitedBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1767:1064–1072.https://doi.org/10.1016/j.bbabio.2007.06.001
-
Plastid thylakoid architecture optimizes photosynthesis in diatomsNature Communications 8:15885.https://doi.org/10.1038/ncomms15885
-
Visualizing density maps with UCSF ChimeraJournal of Structural Biology 157:281–287.https://doi.org/10.1016/j.jsb.2006.06.010
-
The origins and evolution of freeze-etch electron microscopyMicroscopy 60 Suppl 1:S3–S29.https://doi.org/10.1093/jmicro/dfr044
-
PyTom: a python-based toolbox for localization of macromolecules in cryo-electron tomograms and subtomogram analysisJournal of Structural Biology 178:177–188.https://doi.org/10.1016/j.jsb.2011.12.003
-
The site of regulation of light capture in Symbiodinium: does the peridinin–chlorophyll a–protein detach to regulate light capture?Biochimica Et Biophysica Acta (BBA) - Bioenergetics 1837:1227–1234.https://doi.org/10.1016/j.bbabio.2014.03.019
-
Sublocalization of cytochrome b6f Complexes in Photosynthetic MembranesTrends in Plant Science 22:574–582.https://doi.org/10.1016/j.tplants.2017.04.004
-
Two-dimensional crystals of photosystem I in higher plant grana marginsJournal of Biological Chemistry 272:19497–19501.https://doi.org/10.1074/jbc.272.31.19497
-
Fine structure of granal thylakoid membrane organization using cryo electron tomographyBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1807:368–374.https://doi.org/10.1016/j.bbabio.2010.11.007
-
Supramolecular organization of photosystem II in green plantsBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1817:2–12.https://doi.org/10.1016/j.bbabio.2011.05.024
-
Computer visualization of three-dimensional image data using IMODJournal of Structural Biology 116:71–76.https://doi.org/10.1006/jsbi.1996.0013
-
Membrane microdomains, rafts, and detergent-resistant membranes in plants and fungiAnnual Review of Plant Biology 64:501–529.https://doi.org/10.1146/annurev-arplant-050312-120103
-
Robust membrane detection based on tensor voting for electron tomographyJournal of Structural Biology 186:49–61.https://doi.org/10.1016/j.jsb.2014.02.015
-
Automated electron microscope tomography using robust prediction of specimen movementsJournal of Structural Biology 152:36–51.https://doi.org/10.1016/j.jsb.2005.07.007
-
Analysis of the thylakoid outer surface. Coupling factor is limited to unstacked membrane regionsThe Journal of Cell Biology 68:30–47.https://doi.org/10.1083/jcb.68.1.30
-
Dynamic regulation of photosynthesis in Chlamydomonas reinhardtiiThe Plant Journal : For Cell and Molecular Biology 82:413–428.https://doi.org/10.1111/tpj.12805
-
Structure and dynamics of thylakoids in land plantsJournal of Experimental Botany 65:1955–1972.https://doi.org/10.1093/jxb/eru090
-
Optimized cryo-focused ion beam sample preparation aimed at in situ structural studies of membrane proteinsJournal of Structural Biology 197:73–82.https://doi.org/10.1016/j.jsb.2016.07.010
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Reversible particle movements associated with unstacking and restacking of chloroplast membranes in vitroThe Journal of Cell Biology 71:136–158.https://doi.org/10.1083/jcb.71.1.136
-
Regulation of chloroplast membrane function: protein phosphorylation changes the spatial organization of membrane componentsThe Journal of Cell Biology 97:1327–1337.https://doi.org/10.1083/jcb.97.5.1327
-
Jumping mode atomic force microscopy on grana membranes from spinachJournal of Biological Chemistry 286:39164–39171.https://doi.org/10.1074/jbc.M111.284844
-
Functional implications of photosystem II crystal formation in photosynthetic membranesJournal of Biological Chemistry 290:14091–14106.https://doi.org/10.1074/jbc.M114.619841
-
Factors influencing PS II particle array formation in Arabidopsis thaliana chloroplasts and the relationship of such arrays to the thermostability of PS IIBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1228:201–210.https://doi.org/10.1016/0005-2728(94)00175-5
-
Control of cell division by a retinoblastoma protein homolog in ChlamydomonasGenes & Development 15:1652–1661.https://doi.org/10.1101/gad.892101
-
Further characterization of the chloroplast grana margins: the non-detergent preparation of granal photosystem I cannot reduce ferredoxin in the absence of NADP+ reductionBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1230:10–22.https://doi.org/10.1016/0005-2728(95)00027-G
-
Organization of the photosystem II centers and their associated antennae in the thylakoid membranes: a comparative ultrastructural, biochemical, and biophysical study of Chlamydomonas wild type and mutants lacking in photosystem II reaction centersThe Journal of Cell Biology 87:728–735.https://doi.org/10.1083/jcb.87.3.728
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (FOR 2092, EN 1194/1-1)
- Benjamin D Engel
Basic Energy Sciences (DE-FG02-91ER20021)
- Atsuko Kanazawa
Max-Planck-Gesellschaft
- Wolfgang Baumeister
Helmholtz-Gemeinschaft
- Benjamin D Engel
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Radostin Danev, Günter Pfeifer, and Matthias Pöge for technical and computational assistance, Stefan Pfeffer and Roberta Croce for helpful discussions, and Karin Engel for critically reading the manuscript. This work was supported by a grant from the Deutsche Forschungsgemeinschaft to BDE (EN 1194/1–1 as part of FOR 2092). AK was supported by the U.S. Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES) under Award number DE-FG02-91ER20021. Additional funding was provided by the Helmholtz Zentrum München and the Max Planck Society.
Copyright
© 2020, Wietrzynski et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 7,634
- views
-
- 911
- downloads
-
- 101
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Microbiology and Infectious Disease
- Plant Biology
Programmed cell death occurring during plant development (dPCD) is a fundamental process integral for plant growth and reproduction. Here, we investigate the connection between developmentally controlled PCD and fungal accommodation in Arabidopsis thaliana roots, focusing on the root cap-specific transcription factor ANAC033/SOMBRERO (SMB) and the senescence-associated nuclease BFN1. Mutations of both dPCD regulators increase colonization by the beneficial fungus Serendipita indica, primarily in the differentiation zone. smb-3 mutants additionally exhibit hypercolonization around the meristematic zone and a delay of S. indica-induced root-growth promotion. This demonstrates that root cap dPCD and rapid post-mortem clearance of cellular corpses represent a physical defense mechanism restricting microbial invasion of the root. Additionally, reporter lines and transcriptional analysis revealed that BFN1 expression is downregulated during S. indica colonization in mature root epidermal cells, suggesting a transcriptional control mechanism that facilitates the accommodation of beneficial microbes in the roots.
-
- Cell Biology
- Plant Biology
Plants distribute many nutrients to chloroplasts during leaf development and maturation. When leaves senesce or experience sugar starvation, the autophagy machinery degrades chloroplast proteins to facilitate efficient nutrient reuse. Here, we report on the intracellular dynamics of an autophagy pathway responsible for piecemeal degradation of chloroplast components. Through live-cell monitoring of chloroplast morphology, we observed the formation of chloroplast budding structures in sugar-starved leaves. These buds were then released and incorporated into the vacuolar lumen as an autophagic cargo termed a Rubisco-containing body. The budding structures did not accumulate in mutants of core autophagy machinery, suggesting that autophagosome creation is required for forming chloroplast buds. Simultaneous tracking of chloroplast morphology and autophagosome development revealed that the isolation membranes of autophagosomes interact closely with part of the chloroplast surface before forming chloroplast buds. Chloroplasts then protrude at the site associated with the isolation membranes, which divide synchronously with autophagosome maturation. This autophagy-related division does not require DYNAMIN-RELATED PROTEIN 5B, which constitutes the division ring for chloroplast proliferation in growing leaves. An unidentified division machinery may thus fragment chloroplasts for degradation in coordination with the development of the chloroplast-associated isolation membrane.