Differentiation and characterization of excitatory and inhibitory ...

11 downloads 0 Views 5MB Size Report
Jan 8, 2018 - Differentiation and characterization of excitatory and inhibitory synapses by cryo-electron tomography and correlative microscopy. Chang-Lu ...
This Accepted Manuscript has not been copyedited and formatted. The final version may differ from this version. A link to any extended data will be provided when the final version is posted online.

Research Articles: Cellular/Molecular

Differentiation and characterization of excitatory and inhibitory synapses by cryo-electron tomography and correlative microscopy 1

1

1

2

2

1

2

Chang-Lu Tao , Yun-Tao Liu , Rong Sun , Bin Zhang , Lei Qi , Sakar Shivakoti , Chong-Li Tian , Peijun 3

2

1,4,5

Zhang , Pak-Ming Lau , Z. Hong Zhou

and Guo-Qiang Bi

1,6

1

National Laboratory for Physical Sciences at the Microscale, and School of Life Sciences, University of Science and Technology of China, Hefei, Anhui 230026, China 2

CAS Key Laboratory of Brain Function and Disease, and School of Life Sciences, University of Science and Technology of China, Hefei, Anhui 230026, China 3

Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Roosevelt Drive Oxford OX37BN, UK 4

The California NanoSystems Institute, University of California, Los Angeles, CA 90095, USA

5

Department of Microbiology, Immunology and Molecular Genetics, University of California, Los Angeles, CA 90095, USA

6

CAS Center for Excellence in Brain Science and Intelligence Technology, and Innovation Center for Cell Signaling Network, University of Science and Technology of China, Hefei, Anhui 230026, China DOI: 10.1523/JNEUROSCI.1548-17.2017 Received: 2 June 2017 Revised: 17 December 2017 Accepted: 24 December 2017 Published: 8 January 2018

Author contributions: C.-L. Tao, Y.-T.L., P.Z., P.-M.L., Z.H.Z., and G.-Q.B. designed research; C.-L. Tao, Y.-T.L., R.S., B.Z., and L.Q. performed research; C.-L. Tao, Y.-T.L., R.S., S.S., C.-L. Tian, and G.-Q.B. analyzed data; C.-L. Tao, Y.-T.L., S.S., P.Z., P.-M.L., Z.H.Z., and G.-Q.B. wrote the paper. Conflict of Interest: The authors declare no competing financial interests. This work was supported in part by grants from the CAS (XDB02050000 to GQB), NSFC (30725017 and 91232722 to GQB; 31070935 to PML), MOST (2009CB941300 to GQB), and NIH (GM071940 to ZHZ). We acknowledge use of instruments at the Center for Integrative Imaging of Hefei National Laboratory for Physical Sciences at the Microscale and those at the Electron Imaging Center for Nanomachines of UCLA supported by NIH (S10RR23057 and S10OD018111) and NSF (DBI-133813). We thank Xiaokang Zhang and Peng Ge for technical advice on imaging and data processing, Jay He for help with design and construction of the cryoCLEM platform, Weidong Yao and Ann Marie Craig for sharing the PSD-95 and mCherrry-gephyrin plasmids, Xiaobing Chen for insightful discussions, Chunhong Qiu for help with illustrations, and the anonymous reviewers for their insightful comments and constructive suggestions to improve the paper. Correspondence: Guo-Qiang Bi, School of Life Science, USTC, 443 Huangshan Rd, Hefei, China, 230027. E-mail: [email protected] and Z. Hong Zhou, The California NanoSystems Institute, UCLA, 570 Westwood Plaza, CA 90095, USA. E-mail: [email protected] Cite as: J. Neurosci ; 10.1523/JNEUROSCI.1548-17.2017 Alerts: Sign up at www.jneurosci.org/cgi/alerts to receive customized email alerts when the fully formatted version of this article is published.

Accepted manuscripts are peer-reviewed but have not been through the copyediting, formatting, or proofreading process. Copyright © 2018 the authors

1

Differentiation and characterization of excitatory and inhibitory synapses by

2

cryo-electron tomography and correlative microscopy

3 4

Abbreviated title: Visualizing synapse types in 3D by cryoCLEM

5 6

Chang-Lu Tao1,*, Yun-Tao Liu1,*, Rong Sun1, Bin Zhang2, Lei Qi2, Sakar Shivakoti1,

7

Chong-Li Tian2, Peijun Zhang3, Pak-Ming Lau2, Z. Hong Zhou1,4,5,#, Guo-Qiang Bi1,6,#

8

1

9

University of Science and Technology of China, Hefei, Anhui 230026, China; 2CAS Key

National Laboratory for Physical Sciences at the Microscale, and School of Life Sciences,

10

Laboratory of Brain Function and Disease, and School of Life Sciences, University of

11

Science and Technology of China, Hefei, Anhui 230026, China; 3Division of Structural

12

Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Roosevelt

13

Drive Oxford OX37BN, UK; 4The California NanoSystems Institute, University of

14

California, Los Angeles, CA 90095, USA; 5Department of Microbiology, Immunology and

15

Molecular Genetics, University of California, Los Angeles, CA 90095, USA; 6CAS Center

16

for Excellence in Brain Science and Intelligence Technology, and Innovation Center for

17

Cell Signaling Network, University of Science and Technology of China, Hefei, Anhui

18

230026, China

19

*These authors contributed equally to this work.

20

#

21

Hefei, China, 230027. E-mail: [email protected] and Z. Hong Zhou, The California

22

NanoSystems Institute, UCLA, 570 Westwood Plaza, CA 90095, USA. E-mail:

23

[email protected].

24

Number of text pages: 51; Number of figures: 9; Number of movies: 5;

Correspondence: Guo-Qiang Bi, School of Life Science, USTC, 443 Huangshan Rd,

1

25

Number of words: Abstract (215); Significance Statement (118); Introduction (640);

26

Discussion (1207).

27 28

Conflict of Interest: The authors declare no competing financial interests.

29 30

Acknowledgements: This work was supported in part by grants from the CAS

31

(XDB02050000 to GQB), NSFC (30725017 and 91232722 to GQB; 31070935 to PML),

32

MOST (2009CB941300 to GQB), and NIH (GM071940 to ZHZ). We acknowledge use of

33

instruments at the Center for Integrative Imaging of Hefei National Laboratory for Physical

34

Sciences at the Microscale and those at the Electron Imaging Center for Nanomachines

35

of UCLA supported by NIH (S10RR23057 and S10OD018111) and NSF (DBI-133813).

36

We thank Xiaokang Zhang and Peng Ge for technical advice on imaging and data

37

processing, Jay He for help with design and construction of the cryoCLEM platform,

38

Weidong Yao and Ann Marie Craig for sharing the PSD-95 and mCherrry-gephyrin

39

plasmids, Xiaobing Chen for insightful discussions, Chunhong Qiu for help with

40

illustrations, and the anonymous reviewers for their insightful comments and constructive

41

suggestions to improve the paper.

42 43

2

44

Abstract

45

As key functional units in neural circuits, different types of neuronal synapses play distinct

46

roles in brain information processing, learning and memory. Synaptic abnormalities are

47

believed to underlie various neurological and psychiatric disorders. Here, by combining

48

cryo-electron tomography and cryo-correlative light and electron microscopy, we

49

distinguished intact excitatory and inhibitory synapses of cultured hippocampal neurons,

50

and visualized the in situ three-dimensional organization of synaptic organelles and

51

macromolecules in their native state. Quantitative analyses of over a hundred synaptic

52

tomograms reveal that excitatory synapses contain a mesh-like postsynaptic density

53

(PSD) with thickness ranging from 20-50 nm. In contrast, the PSD in inhibitory synapses

54

assumes a thin sheet-like structure ~12 nm from the postsynaptic membrane. On the

55

presynaptic side, spherical synaptic vesicles (SVs) of 25-60 nm diameter and

56

discus-shaped ellipsoidal SVs of various sizes coexist in both synaptic types, with more

57

ellipsoidal ones in inhibitory synapses. High-resolution tomograms obtained using a Volta

58

phase plate and electron filtering and counting reveal glutamate receptor-like and GABAA

59

receptor-like structures that interact with putative scaffolding and adhesion molecules,

60

reflecting details of receptor anchoring and PSD organization. These results provide an

61

updated view of the ultrastructure of excitatory and inhibitory synapses, and demonstrate

62

the potential of our approach to gain insight into the organizational principles of cellular

63

architecture underlying distinct synaptic functions.

64

3

65

Significance Statement

66

To understand functional properties of neuronal synapses, it is desirable to analyze their

67

structure at molecular resolution. We have developed an integrative approach combining

68

cryoET and correlative fluorescence microscopy to visualize 3D ultrastructural features of

69

intact excitatory and inhibitory synapses in their native state. Our approach shows that

70

inhibitory synapses contain uniform thin sheet-like PSDs, while excitatory synapses

71

contain previously known mesh-like PSDs. We discovered “discus-shaped” ellipsoidal

72

synaptic vesicles, and their distributions along with regular spherical vesicles in synaptic

73

types are characterized. High-resolution tomograms further allowed identification of

74

putative neurotransmitter receptors and their heterogeneous interaction with synaptic

75

scaffolding proteins. The specificity and resolution of our approach enables precise in situ

76

analysis of ultrastructural organization underlying distinct synaptic functions.

77

4

78

Introduction

79

Chemical synapses are basic functional units in neural circuits for information

80

transmission, processing and storage (Eccles, 1964; Sudhof and Malenka, 2008; Mayford

81

et al., 2012). From the enormous number of synapses in the brain, the plasticity of each

82

synapse, and the molecular and functional diversity across these synapses, arise the

83

brain’s remarkable computational power and cognitive capacity (Milner et al., 1998; Bi

84

and Poo, 2001). Glutamatergic and GABAergic synapses, the two main types of central

85

synapses, play opposite roles in excitation and inhibition. They have been shown by

86

biochemical and electrophysiological studies to contain different sets of molecular and

87

cellular components and to exhibit distinct functional properties and plasticity rules (Craig

88

and Boudin, 2001; Sudhof and Malenka, 2008; Vogels and Abbott, 2009;

89

Sassoe-Pognetto et al., 2011). How are these different components organized into the

90

intricate machinery to perform distinct synaptic functions? Electron microscopy (EM) has

91

been a primary probe to address this question by visualizing the ultrastructure of various

92

synapses (Siksou et al., 2007; Harris and Weinberg, 2012).

93

Classical EM uses chemical fixation, dehydration and plastic embedding, followed by

94

sectioning and heavy-metal staining to image brain tissues and cultured neurons. Such

95

meticulous processing has allowed electron beams to image various biological

96

specimens at high contrast. Indeed, classical EM observations have shaped much of our

97

current knowledge about synaptic ultrastructure (Sorra and Harris, 2000; Harris and

98

Weinberg, 2012). For example, prominent subcellular features such as the postsynaptic

99

density (PSD) and synaptic vesicles (SVs) are well documented, especially for excitatory

100

synapses (Gray, 1959; Colonnier, 1968; Schikorski and Stevens, 1997; Harris and 5

101

Weinberg, 2012). The 3D resolving capability of electron tomography (ET) has yielded

102

better views of synaptic ultrastructure (Harlow et al., 1998; Ress et al., 2004; Burette et

103

al., 2012). The improved structural preservation provided by high-pressure freezing with

104

freeze substitution (HPF-FS) (Tatsuoka and Reese, 1989) combined with ET (Rostaing et

105

al., 2006) has allowed studies of the organization and dynamics of SVs (Siksou et al.,

106

2007; Watanabe et al., 2013; Imig et al., 2014; Jung et al., 2016) and the 3D organization

107

of macromolecular complexes in individual synapses, such as the PSD-95/glutamate

108

receptor complex at the PSD (Chen et al., 2008; Chen et al., 2015).

109

It remains challenging to characterize the structure and organization of cellular and

110

molecular machinery of specific synaptic types at higher resolution (Hurbain and Sachse,

111

2011; Tao et al., 2012), because damage or deformation from sample preparation

112

procedures can complicate structural interpretation (Hurbain and Sachse, 2011).

113

Cryo-electron tomography (cryoET), which aims to overcome this limitation, has been

114

used to visualize distribution of SVs and other ultrastructural features in isolated

115

synaptosomes and cryo-sections of brain tissues (Fernandez-Busnadiego et al., 2010;

116

Shi et al., 2014; Wilhelm et al., 2014; Perez de Arce et al., 2015). However, cryoET alone

117

cannot unambiguously identify synapse types due to a lack of specific labeling such as

118

immunogold staining or photoconversion of diaminobenzidine for classical EM (Megias et

119

al., 2001; Schikorski and Stevens, 2001; Rostaing et al., 2006; Harris and Weinberg,

120

2012). One way to overcome this is to take advantage of the molecular specificity of

121

fluorescence labeling and sample preservation of cryoET in cryo-correlative light and

122

electron microscopy (cryoCLEM), as suggested earlier (Lucic et al., 2007), although the

123

ability of this approach to distinguish different synapse types is yet to be realized. In the

124

current study, we developed an efficient cryoCLEM platform to identify different types of 6

125

synapses in cultured hippocampal neurons, and to define pre- and postsynaptic

126

ultrastructural features of excitatory and inhibitory synapses in their native state. By

127

high-resolution cryoET with cutting-edge direct electron detection (Li et al., 2013), Volta

128

phase plate (VPP) (Danev et al., 2014; Fukuda et al., 2015), and electron energy filter

129

(Verbeeck et al., 2004) technologies, we could also visualize putative glutamate and

130

GABAA receptors and their organization at the postsynaptic membrane of excitatory and

131

inhibitory synapses.

132

Materials and Methods

133

The overall workflow of experimental procedures is illustrated in Figure 1A. Primary

134

neuronal cultures were grown on EM grids and then plunge-frozen for cryoET imaging

135

followed by 3D reconstruction. For some cultures transfected with constructs of

136

fluorescent protein-tagged synaptic proteins, cryo-fluorescence microscopy was

137

performed before cryoET for correlative imaging. All animal procedures were performed

138

following the guidelines of the Animal Experiments Committee at the University of

139

Science and Technology of China.

140

Primary culture of hippocampal neurons

141

Low-density cultures of dissociated embryonic rat hippocampal neurons were prepared

142

as previously described (Bi and Poo, 1998) with modifications. Quantifoil R2/2 gold EM

143

grids (200 mesh with holey carbon film of 2 μm hole size and 2 μm spacing) or Quantifoil

144

R2/2 gold NH2 finder grids (100 mesh with holey carbon film of 2 μm hole size and 2 μm

145

spacing) were plasma cleaned with H2 and O2 for 10 s using a plasma cleaning system

146

(Gatan), and sterilized with UV light for 30 min. These grids were then coated with 7

147

poly-L-lysine (Sigma) overnight, followed by washing with Hank’s balanced salt solution

148

and ddH2O for about 12 hours each. Hippocampi were removed from embryonic day-18

149

(E18) rats (without distinguishing sex differences) and were treated with trypsin for 15 min

150

at 37 °C, followed by washing and gentle trituration. The dissociated cells were plated on

151

the poly-L-lysine coated EM grids in 35-mm Petri dishes at a density of 40,000-60,000

152

cells/ml, and maintained in incubators at 37 °C in 5% CO2. The culture medium was

153

NeuroBasal (NB) (Invitrogen) supplemented with 5% heat-inactivated bovine calf serum

154

(PAA) plus 5% heat-inactivated fetal bovine serum (HyClone), 1X Glutamax (Invitrogen)

155

and 1X B27 (Invitrogen). Twenty-four hours after plating, half of the medium was replaced

156

by serum-free culture medium. Subsequently, one third of the culture medium was

157

replaced with fresh culture medium twice a week. For correlative microscopy, cultures

158

were co-infected with lentiviruses encoding PSD-95-EGFP and mCherry-gephyrin

159

constructs (see below) for 5-7 days in vitro before vitrification of the grid. Twelve hours

160

after the infection, half of the culture medium was replaced by fresh medium.

161

To prevent overgrowth of glial cells, the cultures were treated with cytosine

162

arabinoside (Sigma) at various stages. Cultures were used for cryoEM imaging at 14-18

163

days in vitro (DIV), when healthy, low-density cultures formed patches of monolayer

164

neuronal processes (Figure 1B1). We judge whether the culture is healthy based on

165

morphological criteria, e.g. smooth soma and dendrites with multiple branches viewed

166

under phase-contrast light microscopy, and at least one probable synapse each few

167

micrometers along the dendrites of transfected neurons, as viewed under fluorescence

168

microscopy. According to our experience, such criteria predict retention of functional

169

properties of synaptic transmission and plasticity evaluated with patch-clamp recording

170

and calcium imaging. 8

171

DNA constructs and lentivirus preparation

172

The PSD-95 cDNA was amplified from GW-PSD-95-EGFP plasmid (a generous gift from

173

Dr. Weidong Yao) and subcloned into pLenti-CaMKII-mKate2 vector to produce

174

pLenti-CaMKII-PSD-95-mKate2. The GGGS-EGFP cDNA was amplified from the

175

pEGFPN1 plasmid, and then subcloned into pLenti-CaMKII-PSD-95-mKate2 to produce

176

the pLenti-CaMKII-PSD-95-EGFP plasmid. The mCherry-gephyrin lentiviral construct

177

(Dobie and Craig, 2011) was a generous gift from Dr. Ann Marie Craig. Both

178

PSD-95-EGFP and mCherry-gephyrin lentiviral constructs were packaged into lentivirus

179

following a protocol from Dr. Karl Deisseroth’s lab (Zhang et al., 2010b).

180

Frozen-hydrated sample preparation

181

After being removed from the CO2 incubator, low-density neuronal cultures (DIV 14-18)

182

on EM grids were first placed in extracellular solution (ECS, containing 150 mM NaCl, 3

183

mM KCl, 3 mM CaCl2, 2 mM MgCl2, 10 mM HEPES and 5 mM glucose, pH 7.3), then

184

mounted on a Vitrobot IV (FEI). Protein A-coated colloidal gold beads (15-nm size, CMC)

185

were added to the grid (4 μl each, stock solution washed in ECS and diluted 10 times

186

after centrifugation) as fiducial markers. The grids were then plunged into liquid ethane

187

for rapid vitrification of the samples, which were then stored in liquid nitrogen until use.

188

CryoET imaging

189

CryoET data were collected with single-axis tilt using either a Tecnai F20 transmission

190

electron microscope (Tecnai F20, FEI) equipped with an Eagle 4K×4K multiport CCD

191

camera (FEI), or a Titan Krios (FEI) with a K2 Summit direct electron detector (K2 camera,

192

Gatan). The Tecnai F20 was operated at an acceleration voltage of 200 KV. Tilt series 9

193

were collected from -60° to +60° at of 2° intervals using FEI Xplore 3D software, with the

194

defocus value set at -12 to -18 μm, and the total electron dosage of ~100 e-/Å2. The final

195

pixel size was 0.755 nm. The Titan Krios were operated at an acceleration voltage of 300

196

KV, with or without Volta phase plate (VPP) and Gatan image filter (GIF). In either

197

configuration, images were collected by the K2 camera in counting mode. In the absence

198

of VPP and GIF, tilt series were acquired from -64° to +64° at 2° intervals using Leginon

199

(Suloway et al., 2005), with the defocus value maintained at -10 μm, and the total

200

accumulated dose of ~120 e-/Å2. The final pixel size was 0.765 nm. When VPP and GIF

201

were used, the energy filter slit was set at 20 eV, and VPP was conditioned by

202

pre-irradiation for 60 s to achieve an initial phase shift of ~0.3π (Fukuda et al., 2015). Tilt

203

series were acquired from -66° to +60° at an interval of 2° or 3° using SerialEM

204

(Mastronarde, 2005) with the defocus value maintained at -1 μm and the total

205

accumulated dose of ~150 e-/Å2. The final pixel size was 0.435 nm.

206

For this study, we examined 78 grids, of which 12 were used for data collection. The

207

rest were discarded because the grids were damaged during transfer or freezing, cultures

208

were too dense and/or too thick, or cultures appeared not healthy with few or no

209

synapses found. Usually 3-5 grid squares (each about 100×100 μm2) per grid were

210

selected for imaging. To obtain high quality cryoET images, it is critical to choose thin

211

culture areas with healthy yet relatively low-density dendrites (Figure 1B2). Generally,

212

areas thicker than ~500 nm were ignored. At this thickness, subcellular structures such

213

as mitochondria, microtubules and synaptic vesicles could not be distinguished in single

214

projection images.

215

Cryo-correlative light and electron microscopy imaging

10

216

The hardware of our cryo-light microscope system includes a custom-built cryo-chamber

217

with liquid nitrogen supply, a Gatan 626 EM cryo-holder, and an Olympus IX71 inverted

218

fluorescence microscope (Figure 1C). The inside channel of the cryo-chamber was

219

precooled to -190 °C by liquid nitrogen, and maintained below -180 °C, as monitored by a

220

thermoelectric sensor. Nitrogen gas flowed through the objective lens and light source

221

windows during the experiment to prevent frost accumulation. Then, an EM grid with

222

frozen-hydrated sample was loaded onto an EM cryo-holder, which was subsequently

223

inserted into the cryo-chamber.

224

For cryoCLEM imaging, fluorescence images were taken using a 40X air objective

225

lens (Olympus LUCPLFLN 40X, NA 0.6) and an ANDOR NEO sCMOS camera (Andor)

226

attached to the fluorescence microscope. For each field of view, three images were

227

collected, one in bright-field, another in the EGFP channel (Ex:470/40, DM: 495, Em:

228

525/50; Chroma, 49002) and the third in the mCherry channel (Ex: 562/40, DM: 593, Em:

229

641/75; Semrock, mCherry-B-000). Typically about 10 sets of images were sufficient to

230

cover all good areas (about 40 squares) on each grid, which took about 20 min to

231

complete; a "good area" was defined as a grid square (~100×100 μm2) of appropriate

232

sample thickness that displayed multiple dendritic branches, usually containing dozens of

233

PSD-95-EGFP puncta or multiple mCherry-gephyrin puncta under fluorescence

234

microscopy (Figure 1D1).

235

Immediately after the light microscopy (LM) imaging, the EM cryo-holder with grid

236

was directly transferred into a Tecnai F20 scope. Areas of the sample imaged in cryo-light

237

microscope were identified in the EM using the indexes of the finder grids (Figure 1D1,

238

D2). Low magnification (330X) EM images were collected and roughly aligned with

239

bright-field LM images using Midas program in IMOD package (RRID: SCR_003297) 11

240

(Kremer et al., 1996). After rough alignment, a set of holes on the carbon layer of the grid

241

were picked from both the low magnification EM images and their corresponding

242

fluorescence images using 3dmod in IMOD package. Transformation functions between

243

the EM and LM images was calculated based on the selected positions by minimizing the

244

mean squared error.

245

When the low magnification EM image and LM image were optimally aligned (Figure

246

1D3), ~15 holes on carbon (in one square) in the low magnification EM image were

247

selected, with their pixel positions recorded. The same holes were identified at 5,000X

248

magnification and the mechanical coordinates were recorded. Afterwards, the

249

transformation function from the pixel positions to EM mechanical coordinates was

250

determined using similar linear regression methods. With the transformation functions,

251

positions of selected fluorescent puncta (putative excitatory or inhibitory synapses) were

252

converted into corresponding EM mechanical coordinates, which were used to guide EM

253

image acquisition. Tilt series were collected on the area with selected fluorescent signals.

254

Finally, tomographic slices were fine-aligned and merged with the fluorescence images to

255

identify each synapse (Figure 1D4) using Midas and ImageJ (RRID: SCR_003070).

256

Python scripts (RRID: SCR_008394) to integrate the correlation procedures are available

257

to interested readers upon request.

258

3D Reconstruction and rendering

259

Tilt series were aligned and reconstructed using IMOD. Gold beads added to the sample

260

before plunge freezing were used as fiducial markers to align the tilt series.

261

Reconstruction was performed using a simultaneous iterative reconstruction technique

262

with 5 or 15 iterations. A fraction of the data collected were discarded during 12

263

reconstruction for technical reasons: for Tecnai F20, ~50% of the data were discarded

264

because of issues such as stage drift, beam blockade at high tilt angles, and occasional

265

auto-focus failure; for Titan Krios, 20 nm) PSD structures 21

469

right next to the plasma membrane were easily identifiable in 13 out of the 14 excitatory

470

synapses (Figure 3A) and spanned nearly the entire area of the uniform synaptic cleft

471

(Figure 3A1). The existence of such thick PSDs is consistent with the common belief that

472

excitatory synapses are “asymmetric”, with dense molecular scaffolds on the

473

postsynaptic side (Colonnier, 1968; Peters and Palay, 1996). In contrast, nearly all (7 out

474

of the 8) inhibitory synapses identified by cryoCLEM (Figure 3B) had distinct thin

475

sheet-like PSD (thin PSD for short), positioned in parallel and close to the postsynaptic

476

membrane (Figure 3B1). Thus, the previously-termed “symmetric” inhibitory synapses

477

are in fact asymmetric under cryoET.

478

Quantitative analyses of PSD structures in excitatory and inhibitory synapses.

479

We quantified the PSD profiles of the above 20 identified synapses together with 90

480

additional synapses with visible PSDs from the 101 synapses obtained with cryoET but

481

not CLEM, by plotting the mean pixel density as a function of its distance to the

482

postsynaptic membrane (Figure 4A, B). The curve contains two major peaks representing

483

the pre- and postsynaptic membranes. There are also two smaller peaks on the curve,

484

one between the pre- and postsynaptic membrane, representing an electron-dense band

485

within the synaptic cleft reported previously as an intermediate band (Gray, 1959), and

486

the other to the right of the postsynaptic membrane defined as “PSD peak,” presumably

487

indicating a postsynaptic proteinaceous layer. This quantitative approach allowed us to

488

identify d1, the peak position of the PSD, and d2, which provides a measure of the

489

thickness of the PSD (Figure 4A2, B2). Scatter plot of PSD thickness and PSD peak

490

positions of all synapses that contain visible PSDs shows two well-defined clusters, thick

491

and thin, which overlap with the two distinct populations formed by the CLEM-identified

492

excitatory and inhibitory synapses, respectively (Figure 4C). Thus, the distinct PSD 22

493

patterns detected by cryoET can be used as hallmarks to distinguish the two types of

494

synapses.

495

With the PSD pattern as an unequivocal criterion, we systematically characterized

496

pre- and postsynaptic features in all 110 synapses with visible PSDs analyzed above,

497

including 85 excitatory synapses and 25 inhibitory synapses. Their postsynaptic densities

498

exhibit distinct PSD peak positions (14.7±3.0 nm, n=85 for excitatory synapse; 9.1±1.1

499

nm, n=25 for inhibitory synapse) and thickness (32.7±7.7 nm, n=85 for excitatory synapse;

500

12.3±1.8 nm, n=25 for inhibitory synapse). Compared to the uniformly thin PSDs of

501

inhibitory synapses, the thick PSDs of excitatory synapses exhibit substantial variability

502

(Figure 4C, D). By averaging all density profiles for excitatory and inhibitory synapses,

503

respectively, we found that the two types of synapses have similar cleft width (~26 nm)

504

(Figure 4E), in contrast to the previous report of narrower cleft for inhibitory synapses

505

(Peters and Palay, 1996). Inside the cleft, distinctive band-like structures are visible in all

506

synapses (Figure 4A, B), as reported previously (Gray, 1959; Zuber et al., 2005). We

507

speculate that these structures are protein complexes involved in cell adhesion (Missler

508

et al., 2012). Intriguingly, the density profile around the presynaptic membrane peak is

509

asymmetric; the density values on the cytoplasmic side are slightly higher than that on the

510

cleft side, especially in excitatory synapses. This presumably reflects extra protein on the

511

cytoplasmic side of the presynaptic membrane, forming a weak version of the active zone

512

commonly observed in conventional EM (Phillips et al., 2001; Sudhof, 2012). Additional

513

studies with better molecular identifications are needed to gain more insight into these

514

structures.

515 516

In addition to synapses with visible PSD, we also observed 13 structures that met our criteria for synapses but exhibited no visible PSD (Figure 5A). They differ from the more 23

517

frequently observed nonsynaptic boutons similar to those reported previously (Bourne et al.,

518

2013), because they had distinct uniform synaptic cleft structures that were absent in the

519

latter. These “PSD-free synapses" could be either excitatory or inhibitory as evidenced

520

from the cryoCLEM data (Figure 5B, C). They might reflect a special transient stage, e.g.,

521

at an early phase of synaptogenesis or on the way towards elimination (Klemann and

522

Roubos, 2011). The variation in the existence and thickness of PSDs reflects potentially

523

diverse synaptic subtypes and associated signaling and structural mechanisms.

524

Quantitative analyses of synaptic vesicles in excitatory and inhibitory synapses

525

On the presynaptic side, we found that in both excitatory and inhibitory synapses, the

526

majority of SVs are perfectly spherical (Figure 6A-C), with an average diameter of about

527

40 nm (40.9±5.0 nm, 16476 vesicles in excitatory synapses; 41.3±6.6 nm, 4766 vesicles

528

in inhibitory synapses, Figure 6D). This differs from the classical EM observation that

529

inhibitory vesicles tend to be smaller compared to excitatory ones (Peters and Palay,

530

1996), but is consistent with findings using HPF-FS (Tatsuoka and Reese, 1989; Korogod

531

et al., 2015), and with measurements from synaptosomes using cryoEM

532

(Fernandez-Busnadiego et al., 2010). Intriguingly, we also observed some apparently

533

ellipsoidal vesicles in both types of synapses (Figure 6A, B). To quantify the shape of

534

different vesicles, we first performed two-dimensional (2D) analysis, based on the

535

maximal cross-section of each SV in the x-y plane, of 19,056 SVs in both excitatory and

536

inhibitory synapses. This analysis identified a relatively small population of ellipsoidal

537

vesicles with major/minor axis ratio significantly greater than 1 (Figure 6E). Interestingly,

538

although ellipsoidal vesicles were found in both excitatory and inhibitory synapses, more

539

were found in the latter (8.5% ± 9%, n=35 in excitatory synapses and 16.9% ± 11%, n=15

540

in inhibitory synapses; p=0.006, two-sample Kolmogorov-Smirnov test). Indeed, 5 out of 24

541

15 inhibitory synapses analyzed (in contrast to 3 out of 35 excitatory synapses) had 20%

542

or more vesicles that were ellipsoidal (Figure 6F).

543

A limitation of 2D analysis is that the true 3D shape of the vesicles could not be

544

determined. For example, if the 2D cross-sections of a vesicle is circular, its 3D shape

545

can be either spherical, “discus-shaped” (i.e., oblate spheroid), or “olive-shaped” (i.e.,

546

prolate spheroid) (Figure 6G). Therefore, we measured the three principal axes of each

547

vesicle using a 3D fitting program (see Materials and Methods) similar to a method

548

described previously (Kukulski et al., 2012). Measurement of axes with this program on

549

70 ellipsoidal vesicles and 70 of their neighboring spherical vesicles in 5 excitatory and 5

550

inhibitory synapses (in high-resolution tomograms obtained with VPP, electron filtering

551

and counting) revealed that the majority of the ellipsoidal vesicles were discus-shaped

552

rather than olive-shaped for both excitatory and inhibitory synapses (Figure 6H). The

553

ellipticity of the ellipsoidal vesicles exhibit large variability (Figure 6H), which may reflect

554

their different compositions and functional roles in synaptic transmission.

555

Visualization of putative receptors and scaffolding proteins in individual synapses

556

New tools that allow cryoET to achieve higher resolution, including direct electron

557

detection, Volta phase plate, and electron energy filter have facilitated characterization of

558

molecular complexes such as proteasome in intact cultured neurons (Asano et al., 2015).

559

Using these tools, we obtained high-quality tilt series (Figure 7A, B and Movie 1) and

560

tomograms with high contrast and high resolution, permitting visualization of features

561

such as the two leaflets of the membrane bilayer, microtubule protofilaments, and

562

putative proteasomes (Figure 7C-E). Two high-quality synaptic tomograms obtained

563

using Volta phase plate, and electron filtering and counting, were selected for further 25

564

study. One is a spine excitatory synapse with thick PSD structure. The other is a shaft

565

inhibitory synapse with thin sheet-like PSD.

566

In the excitatory synapse (Figure 8A), large features such as membranous

567

organelles and ribosomes, as well as actin and microtubule filaments with ultrastructural

568

details were readily identified and segmented (Figure 8B and Movie 2). Furthermore,

569

numerous particles and filamentous structures of various sizes and shapes were

570

visualized within and across pre- and postsynaptic compartments (Figure 8C). These

571

structures were presumably individual protein molecules and complexes. Of special

572

interest were particles near the postsynaptic membrane, some of which have shapes

573

similar to that known for glutamate receptors, including N-methyl-D-aspartate receptors

574

(NMDARs) and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors

575

(AMPARs), which constitute a major fraction of the postsynaptic membrane proteins

576

(Valtschanoff and Weinberg, 2001; Chen et al., 2008; Dani et al., 2010; Jacob and

577

Weinberg, 2015). We thus visually classified these particles as putative “glutamate

578

receptors”, and defined the remaining particles visible on the cleft side of the postsynaptic

579

membrane as “non-receptor” particles (Figure 8C). Plotting the length and width of all

580

these particles reveals that the visually identified putative glutamate receptors form a

581

cluster, although not well separated from the “non-receptor” particles (Figure 8D). The

582

average length (12.1±1.4 nm, n=81) and width (8.6±1.4 nm, n=81) of particles in this

583

cluster is similar to that of the extracellular domain of AMPARs (length: 12.0±0.2 nm,

584

width: 10.5±2.4 nm) and NMDARs (length: 10.5±0.2 nm, width: 10.3±1.4 nm) based on

585

their crystal structures [see detailed calculation of averaged dimensions in Materials and

586

Methods] (Figure 8D-D2). In total, this synapse contained 81 putative glutamate

587

receptors, intermingled with other membrane proteins to occupy the surface of the 26

588

postsynaptic membrane area (Figure 8E, F and Movie 3). This number agrees with

589

estimates of the total number of AMPARs and NMDARs in a glutamatergic synapse

590

based on quantitative immunoEM (Nusser et al., 1998b; Takumi et al., 1999), quantitative

591

mass spectrometry (Sheng and Hoogenraad, 2007; Lowenthal et al., 2015) and visual

592

identification with ET after HPF-FS (Chen et al., 2008; Chen et al., 2015).

593

Among the 81 receptor-like structures in the excitatory synapse, 16 displayed a clear

594

globular density (~10 nm in diameter) on the cytoplasmic side (Figure 8G1). Such

595

globular densities are unlikely to belong to the MAGUK-family proteins, which have

596

filamentous shapes (Nakagawa et al., 2004; Chen et al., 2008; Chen et al., 2011). We

597

thus suspect that these 16 structures are likely NMDARs, known to have much larger

598

cytoplasmic domains than AMPARs (Chen et al., 2008). Of the remaining 65 receptor-like

599

structures, 11 had relatively low image quality and thus prevented classification based on

600

their cytoplasmic structures, whereas 54 structures could be visually identified as

601

AMPAR-like structures based on the lack of large globular density. Among them, 44 were

602

found to each link to one or two filamentous structure that might represent PSD-95 or

603

similar MAGUK-family proteins (Figure 8G2, G3). These putative scaffolding structures

604

appeared to contact the cytoplasmic side of AMPAR-like structures (Figure 8G2, G3),

605

reminiscent of PSD-95 anchoring AMPAR through its interaction with stargazin, which

606

binds to the side of the AMPAR (Meyer et al., 2004; Nakagawa et al., 2006). We also

607

observed 10 AMPAR-like structures not associated with any PSD-95-like structures

608

(Figure 8G4). The majority of these PSD-95-like structures linking to AMPAR-like

609

structures were in near-perpendicular orientation with respect to the postsynaptic

610

membrane (Figure 8G). Together with about two hundred similar filaments connecting

611

directly to the membrane, they form a set of “vertical pillars” to shape an overall core 27

612

structure of the PSD (Movie 3), as also seen by ET of samples prepared using HPF-FS

613

(Chen et al., 2008).

614

A high-resolution tomogram of an inhibitory synapse also revealed rich ultrastructural

615

details (Figure 9A-C and Movie 4). On the cleft side of the postsynaptic membrane, many

616

particles (Figure 9C) were found with shapes similar to that of the type-A γ-aminobutyric

617

acid receptors (GABAAR), the primary inhibitory transmitter receptor in these

618

hippocampal neurons (Bi and Poo, 1998; Nusser et al., 1998a). With visual inspection,

619

we provisionally identified about 143 particles as GABAAR, and about 109 other particles

620

visible on the extracellular side of the postsynaptic membrane as “non-receptor” particles

621

that likely represent other synaptic proteins such as adhesion molecules. Plotting the

622

length and width of all these particles revealed that the visually-identified putative

623

GABAARs formed a cluster, and that the sizes of these putative GABAARs (length:

624

7.1±0.9 nm, width: 5.9±0.9 nm, n=143) were similar to that of the extracellular domain of

625

GABAAR based on its crystal structure (length: 6.2±0.1 nm, width: 6.4±0.1 nm) (Figure

626

9D-D2). Note that the averaged putative GABAAR, similar to the averaged putative

627

glutamate receptor, is sounded by a “halo” (Figure 8D1, 9D1), which could be partially

628

due to fringes arising from uncorrected contrast transfer function. However, for individual

629

particles, such effects appeared to be minimal and did not affect the measurements of

630

particle sizes. The majority of the putative “non-receptor” particles were uniformly skinny

631

but with variable lengths (Figure 9D). Within this inhibitory synapse, the 143 putative

632

GABAARs lying amidst 109 “non-receptor” membrane protein molecules covered the

633

entire ~0.1 μm2 postsynaptic membrane (Figure 9E, F and Movie 5). The number and

634

density of GABAAR-like particles are consistent with a previous estimate of 30-200

635

GABAARs per GABAergic synapse (1,250 receptors/μm2) (Nusser et al., 1997). 28

636

The majority of the GABAAR-like particles were associated with one or two

637

“hammer-shaped” structures on the cytoplasmic side, each with a dense “head” and a

638

thin “neck” (Figure 9G). The heads of these hammer-shaped structures were consistently

639

located at ~12 nm from the postsynaptic membrane, and the necks bridged the

640

transmembrane domain between the GABAAR-like structure and the dense head (Figure

641

9G1-3). We speculate that these hammer-shaped structures are protein complexes

642

containing gephyrin molecules, the major postsynaptic scaffolding component of

643

inhibitory synapse (Tretter et al., 2012). Some of the GABAAR-like structures had only

644

thin necks on their cytoplasmic side, lacking the head density (Figure 9G4), suggesting

645

that the neck might be the cytosolic domain of GABAAR. These putative receptor-linked

646

gephyrin-like structures, together with similar particles not linked to receptor-like particles

647

but also lying ~12 nm from the postsynaptic membrane, appeared to form a cross-linked

648

matrix, which could provide anchoring sites and structural support for GABAARs, as

649

previously proposed (Tyagarajan and Fritschy, 2014). Interestingly, we also found many

650

GABAAR-like structures also linked to densities on the cleft side (Figure 9G3, G4). These

651

densities might represent the cell adhesion molecule neurexin, previously reported to

652

bind directly to GABAAR (Zhang et al., 2010a).

653

Discussion

654

The complex and highly organized molecular machinery inside neuronal synapses

655

provides the structural basis for synaptic transmission and plasticity. In this study, we

656

have developed an approach of cryo-correlative microscopy to distinguish excitatory and

657

inhibitory synapses in intact neurons in culture and to visualize their 3D structures in their

658

native state. By quantifying ultrastructural features of more than a hundred hippocampal 29

659

synapses, we have characterized the ultrastructural features across excitatory and

660

inhibitory central synapses. Because the neurons we used were from embryonic brain

661

and cultured for only a couple of weeks, their synapses may not be as mature as those in

662

the adult brain. Nonetheless, such cultured neurons has been shown to exhibit basic

663

physiological properties of synaptic transmission and plasticity similar to those in more

664

intact preparations such as brain slices. Thus, it is likely that the basic ultrastructural

665

features we observed also reflect synaptic architecture in the brain, at least during its

666

early development.

667

Our approach allows for unequivocal differentiation of the ultrastructure of excitatory

668

and inhibitory synapses in hippocampal neurons. Our results show that excitatory

669

synapses have distinct thick PSDs, and inhibitory synapses have more uniformly thin

670

PSDs. These corroborate classic findings regarding the ultrastructure of excitatory and

671

inhibitory synapses (Colonnier, 1968; Gray, 1969), while providing an updated

672

description of different PSD types in their native state. Quantitative analysis reveals that

673

excitatory PSDs have a broad distribution of thickness. This is consistent with the idea

674

that the dynamically interacting PSD molecules may be in a mixed gel/liquid phase (Zeng

675

et al., 2016), and suggests the existence of multiple structural configurations, perhaps

676

reflecting different states of their activation (Dosemeci et al., 2001) and plasticity (Bi and

677

Poo, 1998; Montgomery and Madison, 2004). The existence of multiple functional and

678

plasticity states in excitatory synapses could be critical for optimal learning and memory

679

storage in neuronal circuits, as suggested by theoretical studies (Fusi et al., 2005; Fusi

680

and Abbott, 2007). In contrast, inhibitory PSDs are much more uniformly and regularly

681

organized, consistent with the meager evidence for structural plasticity in such synapses.

682

Intriguingly, there is a “gap” area with relatively lower electron density between the 30

683

postsynaptic membrane and the PSD peak in both excitatory and inhibitory synapses

684

(see the density profiles of Figure 4A2, B2). We suspect that much of this gap reflects the

685

relatively lower protein content here, compared to the PSD peak where more extensive

686

interactions between protein molecules may occur.

687

High-resolution cryoET has allowed visualization of the distinct molecular

688

organization underlying the different functional properties of excitatory and inhibitory

689

synapses. In the excitatory synapse, a set of PSD-95-like filamentous structures formed

690

“vertical pillars” immediately underneath the postsynaptic membrane to organize the thick

691

PSD meshwork. These filaments have not been observed in classic studies with

692

conventional EM (Gray, 1959; Colonnier, 1968; Peters and Palay, 1996), or later studies

693

by ET with brain slices prepared using HPF-FS (Rostaing et al., 2006; Siksou et al., 2007),

694

but is consistent with more recent observation using HPF-FS and ET in cultures (Chen et

695

al., 2008). The discrepancies could be due to differences in image resolution, or could

696

reflect differences between synapses in slice and culture. Consistent with previous

697

observations (Chen et al., 2008; Chen et al., 2015), we also found that some of the

698

vertical filaments were linked to putative NMDA and AMPA type glutamate receptors.

699

However, the cytoplasmic domains of the putative receptors we observed are much

700

smaller. We speculate that the previously-seen large cytoplasmic domains were due to

701

heavy metal staining.

702

In the inhibitory synapse, the thin sheet-like PSD is likely to represent putative

703

gephyrin complexes forming a layer of interacting “heads”, which connect to putative

704

GABAARs and the postsynaptic membrane through filamentous thin “necks”. This is

705

consistent with the hypothesis that gephyrin molecules form a single layer comprising a

706

hexagonal planar lattice that provides docking sites of GABAARs (Tretter et al., 2012; 31

707

Heine et al., 2013). Such interactions could provide a stable structural matrix to ensure

708

efficient synaptic transmission.

709

On the presynaptic side, fine details of synaptic vesicle organization were revealed,

710

including filamentous structures tethering vesicles to the presynaptic membrane, as

711

observed previously in slices using ET with HPF-FS (Siksou et al., 2007), as well as in

712

slices and isolated synaptosomes using cryoET (Fernandez-Busnadiego et al., 2010;

713

Fernandez-Busnadiego et al., 2013). Of considerable interest is the presence of

714

discus-shaped ellipsoidal vesicles in both excitatory and inhibitory synapses.

715

“Pleomorphic” vesicles have been observed since early EM studies of the synapse, and

716

were generally considered an indicator of inhibitory synapses (Uchizono, 1965). The

717

nature of such vesicles has been debated as more recent studies indicated that their

718

occurrence was associated with specific conditions of sample processing (Tatsuoka and

719

Reese, 1989; Peters and Palay, 1996; Korogod et al., 2015). Although synapses of

720

cultured hippocampal neurons may have different characteristics compared to mature

721

synapses in brain slices, our results suggest that ellipsoidal vesicles could exist in both

722

excitatory and inhibitory native synapses, but their existence may not be used as a

723

definitive criterion to classify synapse types.

724

A technological advantage of our approach was the ability to use cryoCLEM for

725

unambiguous identification of excitatory and inhibitory synapse types. This is potentially

726

extendable to broader applications especially in light of functional heterogeneity of

727

synapses in neuronal circuits (Dobrunz and Stevens, 1997; Bi and Poo, 2001; Craig and

728

Boudin, 2001; Vogels and Abbott, 2009; Letellier et al., 2016). Compared to immuno-EM

729

labeling, the use of fluorescent protein tagging ensured high label density, while avoiding

730

additional staining steps that can cause significant structural distortions and structural 32

731

artifacts. The platform we developed here uses the same EM cryo-holder for shuttling the

732

sample between the light and electron microscopes. Besides being more convenient for

733

reliable correlation between light microscopy and EM, this method also avoids repeated

734

grid transfers, thus protecting the sample from potential damage and contamination. This

735

together with our method for accurate correlation between light and electron microscopy

736

greatly improved the efficiency of our approach, and was key to the success of our

737

cryoCLEM experiments.

738

Another advantage of our approach is that direct plunge-freezing of neurons cultured

739

on EM grids at low density prevents unwanted disturbance to the synapses unavoidable

740

in the synaptosome preparations previously used for cryoET studies

741

(Fernandez-Busnadiego et al., 2010; Shi et al., 2014; Perez de Arce et al., 2015).

742

Therefore, the ultrastructure of intact synapses can be preserved near their physiological

743

state. However, plunge-freezing is limited to monolayer cultured neurons and synapses

744

that are no more than a few hundred nanometer thick. High-pressure freezing and

745

cryo-sectioning could provide an appropriate tool to extend cryoET to native circuits in

746

brain tissue (Zuber et al., 2005). By implementing the latest cryoET technologies

747

including VPP, electron energy filter and direct electron detection, which greatly improve

748

resolution and signal-to-noise ratio (Danev et al., 2014; Fukuda et al., 2015), we were

749

able to visualize synaptic ultrastructural features. Individual molecules in the synapses,

750

such as GABAA receptors that were previously not accessible, could be identified,

751

localized and counted, thus providing a straightforward way to study key synaptic proteins

752

in individual synapse. With further technical improvements along the lines outlined here,

753

future studies with 3D classification and sub-tomogram averaging could identify additional

33

754

synaptic proteins more confidently, and reveal finer structural details of protein

755

complexes in situ.

756

References

757 758 759

Asano S, Fukuda Y, Beck F, Aufderheide A, Forster F, Danev R, Baumeister W (2015) Proteasomes. A molecular census of 26S proteasomes in intact neurons. Science 347:439-442.

760 761

Bi G, Poo M (2001) Synaptic modification by correlated activity: Hebb's postulate revisited. Annu Rev Neurosci 24:139-166.

762 763

Bi GQ, Poo MM (1998) Synaptic modifications in cultured hippocampal neurons: dependence on spike timing, synaptic strength, and postsynaptic cell type. J Neurosci 18:10464-10472.

764 765

Bourne JN, Chirillo MA, Harris KM (2013) Presynaptic ultrastructural plasticity along CA3-->CA1 axons during long-term potentiation in mature hippocampus. J Comp Neurol 521:3898-3912.

766 767

Brandt F, Carlson LA, Hartl FU, Baumeister W, Grunewald K (2010) The three-dimensional organization of polyribosomes in intact human cells. Molecular cell 39:560-569.

768 769

Burette AC, Lesperance T, Crum J, Martone M, Volkmann N, Ellisman MH, Weinberg RJ (2012) Electron Tomographic Analysis of Synaptic Ultrastructure. J Comp Neurol 520:2697-2711.

770 771

Chen X, Winters C, Azzam R, Li X, Galbraith JA, Leapman RD, Reese TS (2008) Organization of the core structure of the postsynaptic density. Proc Natl Acad Sci U S A 105:4453-4458.

772 773 774

Chen X, Levy JM, Hou A, Winters C, Azzam R, Sousa AA, Leapman RD, Nicoll RA, Reese TS (2015) PSD-95 family MAGUKs are essential for anchoring AMPA and NMDA receptor complexes at the postsynaptic density. Proc Natl Acad Sci U S A 112:E6983-6992.

775 776 777

Chen X, Nelson CD, Li X, Winters CA, Azzam R, Sousa AA, Leapman RD, Gainer H, Sheng M, Reese TS (2011) PSD-95 is required to sustain the molecular organization of the postsynaptic density. J Neurosci 31:6329-6338.

778 779

Colonnier M (1968) Synaptic patterns on different cell types in the different laminae of the cat visual cortex. An electron microscope study. Brain Res 9:268-287.

780 781

Craig AM, Boudin H (2001) Molecular heterogeneity of central synapses: afferent and target regulation. Nature neuroscience 4:569-578. 34

782 783 784

Danev R, Buijsse B, Khoshouei M, Plitzko JM, Baumeister W (2014) Volta potential phase plate for in-focus phase contrast transmission electron microscopy. Proc Natl Acad Sci U S A 111:15635-15640.

785 786

Dani A, Huang B, Bergan J, Dulac C, Zhuang X (2010) Superresolution imaging of chemical synapses in the brain. Neuron 68:843-856.

787 788 789

Dobie FA, Craig AM (2011) Inhibitory synapse dynamics: coordinated presynaptic and postsynaptic mobility and the major contribution of recycled vesicles to new synapse formation. J Neurosci 31:10481-10493.

790 791

Dobrunz LE, Stevens CF (1997) Heterogeneity of release probability, facilitation, and depletion at central synapses. Neuron 18:995-1008.

792 793 794

Dosemeci A, Tao-Cheng JH, Vinade L, Winters CA, Pozzo-Miller L, Reese TS (2001) Glutamate-induced transient modification of the postsynaptic density. Proc Natl Acad Sci U S A 98:10428-10432.

795 796 797

Durr KL, Chen L, Stein RA, De Zorzi R, Folea IM, Walz T, McHaourab HS, Gouaux E (2014) Structure and dynamics of AMPA receptor GluA2 in resting, pre-open, and desensitized states. Cell 158:778-792.

798

Eccles JC (1964) The physiology of synapses. Berlin: Springers.

799 800 801

Fernandez-Busnadiego R, Zuber B, Maurer UE, Cyrklaff M, Baumeister W, Lucic V (2010) Quantitative analysis of the native presynaptic cytomatrix by cryoelectron tomography. The Journal of cell biology 188:145-156.

802 803 804

Fernandez-Busnadiego R, Asano S, Oprisoreanu AM, Sakata E, Doengi M, Kochovski Z, Zurner M, Stein V, Schoch S, Baumeister W, Lucic V (2013) Cryo-electron tomography reveals a critical role of RIM1alpha in synaptic vesicle tethering. The Journal of cell biology 201:725-740.

805 806

Fitzgibbon A, Pilu M, Fisher RB (1999) Direct least square fitting of ellipses. IEEE transactions on pattern analysis and machine intelligence 21:476-480.

807 808

Fukuda Y, Laugks U, Lucic V, Baumeister W, Danev R (2015) Electron cryotomography of vitrified cells with a Volta phase plate. J Struct Biol 190:143-154.

809 810

Fusi S, Abbott LF (2007) Limits on the memory storage capacity of bounded synapses. Nature neuroscience 10:485-493.

35

811 812

Fusi S, Drew PJ, Abbott LF (2005) Cascade models of synaptically stored memories. Neuron 45:599-611.

813 814

Gray EG (1959) Axo-somatic and axo-dendritic synapses of the cerebral cortex: an electron microscope study. Journal of anatomy 93:420-433.

815 816

Gray EG (1969) Electron microscopy of excitatory and inhibitory synapses: a brief review. Progress in brain research 31:141-155.

817 818

Harlow M, Ress D, Koster A, Marshall RM, Schwarz M, McMahan UJ (1998) Dissection of active zones at the neuromuscular junction by EM tomography. J Physiol-Paris 92:75-78.

819 820

Harris KM, Weinberg RJ (2012) Ultrastructure of synapses in the mammalian brain. Cold Spring Harbor perspectives in biology 4.

821 822

Heine M, Karpova A, Gundelfinger ED (2013) Counting gephyrins, one at a time: a nanoscale view on the inhibitory postsynapse. Neuron 79:213-216.

823 824 825

Hrabe T, Chen Y, Pfeffer S, Cuellar LK, Mangold AV, Forster F (2012) PyTom: a python-based toolbox for localization of macromolecules in cryo-electron tomograms and subtomogram analysis. J Struct Biol 178:177-188.

826 827 828

Hurbain I, Sachse M (2011) The future is cold: cryo-preparation methods for transmission electron microscopy of cells. Biology of the cell / under the auspices of the European Cell Biology Organization 103:405-420.

829 830 831

Imig C, Min SW, Krinner S, Arancillo M, Rosenmund C, Sudhof TC, Rhee J, Brose N, Cooper BH (2014) The morphological and molecular nature of synaptic vesicle priming at presynaptic active zones. Neuron 84:416-431.

832 833

Jacob AL, Weinberg RJ (2015) The organization of AMPA receptor subunits at the postsynaptic membrane. Hippocampus 25:798-812.

834 835

Jung JH, Szule JA, Marshall RM, McMahan UJ (2016) Variable priming of a docked synaptic vesicle. Proc Natl Acad Sci U S A 113:E1098-1107.

836 837

Klemann CJ, Roubos EW (2011) The gray area between synapse structure and function-Gray's synapse types I and II revisited. Synapse 65:1222-1230.

838 839

Korogod N, Petersen CC, Knott GW (2015) Ultrastructural analysis of adult mouse neocortex comparing aldehyde perfusion with cryo fixation. eLife 4.

36

840 841

Kremer JR, Mastronarde DN, McIntosh JR (1996) Computer visualization of three-dimensional image data using IMOD. J Struct Biol 116:71-76.

842 843

Kukulski W, Schorb M, Kaksonen M, Briggs JA (2012) Plasma membrane reshaping during endocytosis is revealed by time-resolved electron tomography. Cell 150:508-520.

844 845

Lee CH, Lu W, Michel JC, Goehring A, Du J, Song XQ, Gouaux E (2014) NMDA receptor structures reveal subunit arrangement and pore architecture. Nature 511:191-+.

846 847 848

Letellier M, Park YK, Chater TE, Chipman PH, Gautam SG, Oshima-Takago T, Goda Y (2016) Astrocytes regulate heterogeneity of presynaptic strengths in hippocampal networks. Proc Natl Acad Sci U S A 113:E2685-2694.

849 850 851

Li XM, Mooney P, Zheng S, Booth CR, Braunfeld MB, Gubbens S, Agard DA, Cheng YF (2013) Electron counting and beam-induced motion correction enable near-atomic-resolution single-particle cryo-EM. Nature methods 10:584.

852 853 854

Liu B, Xue Y, Zhao W, Chen Y, Fan C, Gu L, Zhang Y, Zhang X, Sun L, Huang X, Ding W, Sun F, Ji W, Xu T (2015) Three-dimensional super-resolution protein localization correlated with vitrified cellular context. Scientific reports 5:13017.

855 856 857

Lowenthal MS, Markey SP, Dosemeci A (2015) Quantitative mass spectrometry measurements reveal stoichiometry of principal postsynaptic density proteins. Journal of proteome research 14:2528-2538.

858 859 860

Lucic V, Kossel AH, Yang T, Bonhoeffer T, Baumeister W, Sartori A (2007) Multiscale imaging of neurons grown in culture: from light microscopy to cryo-electron tomography. J Struct Biol 160:146-156.

861 862

Mastronarde DN (2005) Automated electron microscope tomography using robust prediction of specimen movements. J Struct Biol 152:36-51.

863 864

Mayford M, Siegelbaum SA, Kandel ER (2012) Synapses and memory storage. Cold Spring Harbor perspectives in biology 4.

865 866

Megias M, Emri Z, Freund TF, Gulyas AI (2001) Total number and distribution of inhibitory and excitatory synapses on hippocampal CA1 pyramidal cells. Neuroscience 102:527-540.

867 868 869

Meyer G, Varoqueaux F, Neeb A, Oschlies M, Brose N (2004) The complexity of PDZ domain-mediated interactions at glutamatergic synapses: a case study on neuroligin. Neuropharmacology 47:724-733.

37

870 871

Miller PS, Aricescu AR (2014) Crystal structure of a human GABAA receptor. Nature 512:270-275.

872 873

Milner B, Squire LR, Kandel ER (1998) Cognitive neuroscience and the study of memory. Neuron 20:445-468.

874 875

Missler M, Sudhof TC, Biederer T (2012) Synaptic cell adhesion. Cold Spring Harbor perspectives in biology 4:a005694.

876 877

Montgomery JM, Madison DV (2004) Discrete synaptic states define a major mechanism of synapse plasticity. Trends Neurosci 27:744-750.

878 879

Nakagawa T, Cheng Y, Sheng M, Walz T (2006) Three-dimensional structure of an AMPA receptor without associated stargazin/TARP proteins. Biological chemistry 387:179-187.

880 881 882

Nakagawa T, Futai K, Lashuel HA, Lo I, Okamoto K, Walz T, Hayashi Y, Sheng M (2004) Quaternary structure, protein dynamics, and synaptic function of SAP97 controlled by L27 domain interactions. Neuron 44:453-467.

883 884

Nusser Z, Cull-Candy S, Farrant M (1997) Differences in synaptic GABA(A) receptor number underlie variation in GABA mini amplitude. Neuron 19:697-709.

885 886

Nusser Z, Hajos N, Somogyi P, Mody I (1998a) Increased number of synaptic GABA(A) receptors underlies potentiation at hippocampal inhibitory synapses. Nature 395:172-177.

887 888 889

Nusser Z, Lujan R, Laube G, Roberts JD, Molnar E, Somogyi P (1998b) Cell type and pathway dependence of synaptic AMPA receptor number and variability in the hippocampus. Neuron 21:545-559.

890 891 892

Perez de Arce K, Schrod N, Metzbower SW, Allgeyer E, Kong GK, Tang AH, Krupp AJ, Stein V, Liu X, Bewersdorf J, Blanpied TA, Lucic V, Biederer T (2015) Topographic Mapping of the Synaptic Cleft into Adhesive Nanodomains. Neuron 88:1165-1172.

893

Peters A, Palay SL (1996) The morphology of synapses. Journal of neurocytology 25:687-700.

894 895 896

Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE (2004) UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem 25:1605-1612.

897 898 899

Phillips GR, Huang JK, Wang Y, Tanaka H, Shapiro L, Zhang WD, Shan WS, Arndt K, Frank M, Gordon RE, Gawinowicz MA, Zhao YM, Colman DR (2001) The presynaptic particle web: Ultrastructure, composition, dissolution, and reconstitution. Neuron 32:63-77. 38

900 901

Ress DB, Harlow ML, Marshall RM, McMahan UJ (2004) Methods for generating high-resolution structural models from electron microscope tomography data. Structure 12:1763-1774.

902 903 904

Rostaing P, Real E, Siksou L, Lechaire JP, Boudier T, Boeckers TM, Gertler F, Gundelfinger ED, Triller A, Marty S (2006) Analysis of synaptic ultrastructure without fixative using high-pressure freezing and tomography. Eur J Neurosci 24:3463-3474.

905 906

Sassoe-Pognetto M, Frola E, Pregno G, Briatore F, Patrizi A (2011) Understanding the molecular diversity of GABAergic synapses. Frontiers in cellular neuroscience 5:4.

907 908

Schikorski T, Stevens CF (1997) Quantitative ultrastructural analysis of hippocampal excitatory synapses. J Neurosci 17:5858-5867.

909 910

Schikorski T, Stevens CF (2001) Morphological correlates of functionally defined synaptic vesicle populations. Nature neuroscience 4:391-395.

911 912

Schorb M, Briggs JA (2014) Correlated cryo-fluorescence and cryo-electron microscopy with high spatial precision and improved sensitivity. Ultramicroscopy 143:24-32.

913 914

Sheng M, Hoogenraad CC (2007) The postsynaptic architecture of excitatory synapses: a more quantitative view. Annu Rev Biochem 76:823-847.

915 916 917

Shi L, Du X, Zhou H, Tao C, Liu Y, Meng F, Wu G, Xiong Y, Xia C, Wang Y, Bi G, Zhou JN (2014) Cumulative effects of the ApoE genotype and gender on the synaptic proteome and oxidative stress in the mouse brain. Int J Neuropsychopharmacol 17:1863-1879.

918 919 920

Siksou L, Rostaing P, Lechaire JP, Boudier T, Ohtsuka T, Fejtova A, Kao HT, Greengard P, Gundelfinger ED, Triller A, Marty S (2007) Three-dimensional architecture of presynaptic terminal cytomatrix. J Neurosci 27:6868-6877.

921 922

Sorra KE, Harris KM (2000) Overview on the structure, composition, function, development, and plasticity of hippocampal dendritic spines. Hippocampus 10:501-511.

923

Sudhof TC (2012) The presynaptic active zone. Neuron 75:11-25.

924 925

Sudhof TC, Malenka RC (2008) Understanding Synapses: Past, Present, and Future. Neuron 60:469-476.

926 927

Suloway C, Pulokas J, Fellmann D, Cheng A, Guerra F, Quispe J, Stagg S, Potter CS, Carragher B (2005) Automated molecular microscopy: the new Leginon system. J Struct Biol 151:41-60.

928 929

Takumi Y, Ramirez-Leon V, Laake P, Rinvik E, Ottersen OP (1999) Different modes of expression of AMPA and NMDA receptors in hippocampal synapses. Nature neuroscience 2:618-624. 39

930 931

Tao CL, Xia CL, Chen XB, Zhou ZH, Bi GQ (2012) Ultrastructural analysis of neuronal synapses using state-of-the-art nano-imaging techniques. Neuroscience bulletin 28:321-332.

932 933

Tatsuoka H, Reese TS (1989) New structural features of synapses in the anteroventral cochlear nucleus prepared by direct freezing and freeze-substitution. J Comp Neurol 290:343-357.

934 935

Tretter V, Mukherjee J, Maric HM, Schindelin H, Sieghart W, Moss SJ (2012) Gephyrin, the enigmatic organizer at GABAergic synapses. Frontiers in cellular neuroscience 6:23.

936 937

Tyagarajan SK, Fritschy JM (2014) Gephyrin: a master regulator of neuronal function? Nat Rev Neurosci 15:141-156.

938 939

Uchizono K (1965) Characteristics of excitatory and inhibitory synapses in the central nervous system of the cat. Nature 207:642-643.

940 941

Valtschanoff JG, Weinberg RJ (2001) Laminar organization of the NMDA receptor complex within the postsynaptic density. J Neurosci 21:1211-1217.

942 943

Verbeeck J, Van Dyck D, Van Tendeloo G (2004) Energy-filtered transmission electron microscopy: an overview. Spectrochim Acta B 59:1529-1534.

944 945

Vogels TP, Abbott LF (2009) Gating multiple signals through detailed balance of excitation and inhibition in spiking networks. Nature neuroscience 12:483-491.

946 947 948

Watanabe S, Rost BR, Camacho-Perez M, Davis MW, Sohl-Kielczynski B, Rosenmund C, Jorgensen EM (2013) Ultrafast endocytosis at mouse hippocampal synapses. Nature 504:242-247.

949 950 951

Wilhelm BG, Mandad S, Truckenbrodt S, Krohnert K, Schafer C, Rammner B, Koo SJ, Classen GA, Krauss M, Haucke V, Urlaub H, Rizzoli SO (2014) Composition of isolated synaptic boutons reveals the amounts of vesicle trafficking proteins. Science 344:1023-1028.

952 953 954

Zeng M, Shang Y, Araki Y, Guo T, Huganir RL, Zhang M (2016) Phase Transition in Postsynaptic Densities Underlies Formation of Synaptic Complexes and Synaptic Plasticity. Cell 166:1163-1175 e1112.

955 956

Zhang C, Atasoy D, Arac D, Yang X, Fucillo MV, Robison AJ, Ko J, Brunger AT, Sudhof TC (2010a) Neurexins physically and functionally interact with GABA(A) receptors. Neuron 66:403-416.

957 958 959

Zhang F, Gradinaru V, Adamantidis AR, Durand R, Airan RD, de Lecea L, Deisseroth K (2010b) Optogenetic interrogation of neural circuits: technology for probing mammalian brain structures. Nature protocols 5:439-456. 40

960 961 962

Zuber B, Nikonenko I, Klauser P, Muller D, Dubochet J (2005) The mammalian central nervous synaptic cleft contains a high density of periodically organized complexes. P Natl Acad Sci USA 102:19192-19197.

963

964

Figure/Movie Legends

965

Figure 1. Imaging primary rat hippocampal neurons with cryoET / cryoCLEM. (A)

966

Illustration of the workflow of cryoET/cryoCLEM imaging of neurons grown on gold EM

967

grids. (B) Representative results from different stages of the workflow. (B1) Light

968

microscopy image of cultured neurons (red arrows indicating cell bodies). (B2) CryoEM

969

image of neuronal processes in one grid square. (B3) A single cryoEM projection image

970

of the boxed area in (B2) showing a synapse-like structure with a presynaptic bouton

971

(Bouton) containing a dense population of SVs (green circles), a postsynaptic spine

972

(Spine), and a relatively uniform cleft (yellow arrow). Inset shows a zoomed-in view of

973

the synaptic bouton area with a dense population of SVs. (B4) A tomographic slice

974

showing fine structure of the same synapse in (B3), which was identified as a spine

975

synapse by following through the tomogram in 3D, with mitochondrion (Mit),

976

microtubules (MT), and SVs (green circles) and superposed with segmented presynaptic

977

membrane (green) and postsynaptic membrane (red). Inset shows a zoomed-in view of

978

the synaptic cleft area with transcleft structures. (C) Schematics depicting main

979

components of cryo-fluorescence light microscope with an EM cryo-holder. (D) Pipeline 41

980

of imaging synapse with cryoCLEM. (D1) Merged cryo-fluorescence and cryo-bright-field

981

light images. (D2) Low magnification cryoEM image including the same grid square. (D3)

982

Merged images of boxed area in (D1) and (D2) after fine alignment. (D4) Tomographic

983

slice of the boxed area in (D3) superimposed with aligned fluorescence image showing

984

the structure of a synapse with a green fluorescent punctum.

985 986

Figure 2. Synapses of various sizes, shapes and ultrastructural details imaged

987

with cryoET. (A-C) Three tomographic slices showing structures of different synapses.

988

In the synapses, structures such as SVs and dense core vesicle in presynaptic boutons

989

(Bouton), microtubules (MT) in boutons and dendritic shaft (Shaft), mitochondria (Mit) in

990

presynaptic bouton and postsynaptic spine (Spine) are clearly visible. (A1-C1)

991

Zoomed-in views of corresponding boxed areas from (A-C) showing thick (dashed

992

parallel lines in (A1 and B1)) and thin (dashed parallel lines in (C1)) PSDs, as well as

993

SVs attached (cyan arrowheads in (A1)) or fused (pink arrowheads in (B1)) to the

994

presynaptic membrane. (D and E) Two synapses sharing the same presynaptic axon

995

(determined by following through their tomograms in 3D), both with thick PSDs (D1 and

996

D2) or both with thin PSDs (E1 and E2), respectively. (F and G) Two synapses sharing

997

the same postsynaptic spine, both with thick PSDs (F1 and F2), or one with thin PSD

998

(G1) and the other with thick PSD (G2). 

999 42

1000

Figure 3. Identification of excitatory and inhibitory synapses with cryoCLEM. (A

1001

and B) Tomographic slices of an excitatory (A) and inhibitory (B) synapse colocalized

1002

with PSD-95-EGFP and mCherry-gephyrin puncta, respectively. (A1 and B1) Zoomed-in

1003

views of the boxed area in (A) and (B) showing the synapse with thick and thin PSD,

1004

respectively. Red dashed lines indicate the range of the PSD.

1005 1006

Figure 4. Quantitative and statistical characterization of excitatory and inhibitory

1007

PSDs. (A and B) Tomographic slices of two synapses with thick and thin PSD

1008

respectively. (A1 and B1) Zoomed-in views of the marked areas in (A) and (B). (A2 and

1009

B2) Normalized density profiles of the two synapses in (A) and (B) respectively with

1010

cross-sectional mean density plotted against distance to postsynaptic membrane (see

1011

Materials and Methods). On the x-axis of this plot, 0 was set to be the position of the

1012

postsynaptic membrane, and positive values are on the postsynaptic side. The density

1013

profiles were normalized against the density values at distance ranging from 100 to 200

1014

nm such that the average density value in this range is zero and their standard deviation

1015

is unity. d1: PSD peak position, d2 is the sum of d1 and the length constant obtained

1016

from the exponential fit of the profile from d1 to the flattened background, to provide a

1017

measure of the thickness of the PSD (see Materials and Methods). (C) Scatter plot of

1018

PSD peak position and PSD thickness of all synapses show two well-defined clusters. (D)

1019

Histogram shows the PSD thickness distribution of all excitatory and inhibitory synapses 43

1020

respectively. (E) Averaged density curve of all excitatory synapses and all inhibitory

1021

synapses respectively.

1022 1023

Figure 5. Synapses without visible PSDs. (A) A 15 nm-thick tomographic slice of an

1024

unidentified synapse imaged by cryoET only. (B and C) A 15 nm-thick tomographic slice

1025

of an excitatory (B) and an inhibitory (C) synapse identified by cryoCLEM superposed

1026

with the fluorescence image of colocalized PSD-95-EGFP and mCherry-gephyrin,

1027

respectively. (A1-C1) Zoomed-in views of the boxed area in (A-C) respectively showing

1028

that no PSD structures are visible in these synapses.

1029 1030

Figure 6. Heterogeneity of synaptic vesicles in excitatory and inhibitory synapses.

1031

(A and B) Tomographic slices of an excitatory and an inhibitory synapse respectively.

1032

Insets are zoomed-in views showing thick and thin PSDs from (A) and (B) respectively.

1033

(C) Scatter plot showing the major and minor axes of SVs in the two synapses in (A and

1034

B) measured by 2D fitting. (D) Distribution of vesicle sizes in excitatory and inhibitory

1035

synapses (16,476 vesicles in 35 excitatory synapses and 4,766 vesicles in 15 inhibitory

1036

synapses). (E) Distribution of ellipticity of SVs (major to minor axis ratio) in excitatory and

1037

inhibitory synapses. A threshold (dashed line) was set at major/minor = 1.14, which is

1038

about twice the peak position (major/minor = 1.07) from perfect circle (major/minor = 1)

1039

to separate ellipsoidal from spherical vesicles. Coincidentally, this threshold is also close 44

1040

to the cross point of the two distribution curves. (F) Cumulative frequency of the fraction

1041

of ellipsoidal vesicles in excitatory and inhibitory synapses. (G) Tomographic slices of

1042

spherical (left column), discus-shaped (i.e., oblate spheroid, middle column) and

1043

olive-shaped (i.e., prolate spheroid, right column) SVs viewed in three orthogonal planes

1044

that were rotated so that each of the three principle axes (a, b, and c) of the vesicles can be

1045

measured horizontally in the corresponding plane. For the spherical vesicles, aҮbҮc; for

1046

discus-shaped ones, a