In the mouse retina, three different types of photoreceptors provide input to 14 bipolar cell (BC) types. Classically, most BC types are thought to contact all cones within their dendritic field; ON-BCs would contact cones exclusively via so-called invaginating synapses, while OFF-BCs would form basal synapses. By mining publically available electron microscopy data, we discovered interesting violations of these rules of outer retinal connectivity: ON-BC type X contacted only ~20% of the cones in its dendritic field and made mostly atypical non-invaginating contacts. Types 5T, 5O and 8 also contacted fewer cones than expected. In addition, we found that rod BCs received input from cones, providing anatomical evidence that rod and cone pathways are interconnected in both directions. This suggests that the organization of the outer plexiform layer is more complex than classically thought.https://doi.org/10.7554/eLife.20041.001
Parallel visual processing already starts at the very first synapse of the visual system, where photoreceptors distribute the signal onto multiple types of bipolar cells. In the mouse retina, two types of cone photoreceptors differing in their spectral properties – short (S-) and medium wavelength-sensitive (M-) cones – and rod photoreceptors provide input to 14 types of bipolar cell (reviewed in Euler et al., 2014). The precise connectivity rules between photoreceptors and bipolar cell (BC) types determine which signals are available to downstream circuits. Therefore, the connectome of the outer retina is essential for a complete picture of visual processing in the retina.
For some mouse BC types, specific connectivity patterns have already been described: For example, based on electrical recordings and immunohistochemistry cone bipolar cell type 1 (CBC1) have been suggested to selectively contact M-cones, whereas CBC9 exclusively contacts S-cones (Haverkamp et al., 2005; Breuninger et al., 2011). The other BC types are thought to contact all M-cones within their dendritic field, but their connectivity to S-cones is unclear (Wässle et al., 2009). In addition, two fundamental cone-BC contact shapes have been described: invaginating contacts with the dendritic tips extending into the cone pedicle and flat (basal) contacts that touch the cone pedicle base, commonly associated with ON- and OFF-BCs, respectively (Dowling and Boycott, 1966; Kolb, 1970; Hopkins and Boycott, 1995).
Rod bipolar cells (RBCs) are thought to exclusively receive rod input and to feed this signal into the cone pathway via AII amacrine cells (reviewed by Bloomfield and Dacheux, 2001). However, physiological data indicate that RBCs may receive cone photoreceptor input as well (Pang et al., 2010). Also, types CBC3A, CBC3B and CBC4 have been reported to receive direct rod input (Mataruga et al., 2007; Haverkamp et al., 2008; Tsukamoto and Omi, 2014), suggesting that rod and cone pathways are much more interconnected than their names implicate.
Here, we analyzed an existing electron microscopy dataset (Helmstaedter et al., 2013) to quantify the connectivity between photoreceptors and bipolar cells in the mouse. We did not find evidence for additional M- or S-cone selective CBC types in addition to the reported CBC1 and 9. However, we found interesting violations of established rules of outer retinal connectivity: The newly discovered CBCX (Helmstaedter et al., 2013), likely an ON-CBC (Ichinose et al., 2014), had unexpectedly few and mostly atypical basal contacts to cones. CBC5T, CBC5O and CBC8 also contacted fewer cones than expected from their dendritic field. In addition, we provide anatomical evidence that rod and cone pathways are connected in both directions: Not only OFF-types CBC3A, CBC3B and CBC4 get direct input from rods but also RBCs from cones.
We used the serial block-face electron microscopy (SBEM) dataset e2006 published by Helmstaedter et al. (2013) to analyze the connectivity between photoreceptors and bipolar cells in the outer plexiform layer (OPL) of the mouse retina (Figure 1A). To this end, we reconstructed the volume of all cone axon terminals (cone pedicles; n = 163) in the dataset as well as the dendritic trees of all BCs (n = 451; Figure 1B, see Materials and methods).
To identify S- and M-cones, we used the fact that type nine cone bipolar cells selectively target S-cones (Figure 1C,D) (Mariani, 1984; Kouyama and Marshak, 1992; Haverkamp et al., 2005; Breuninger et al., 2011). We found 48 contacts of CBC9s and cones, involving 43 cones (Figure 1—figure supplement 1A). We visually assessed all contacts and found that 29 of these were in the periphery of the cone pedicle, where no synapses are expected (Figure 1—figure supplement 1B) (Dowling and Boycott, 1966; Chun et al., 1996). This left 14 potential S-cones with invaginating contacts from at least one CBC9. It has been shown that S-cones are contacted by the dendrites of all neighboring CBC9s and that these contacts occur mostly at the tip of dendritic branches (Haverkamp et al., 2005). Out of the 14 candidate cones, eight cones had only one CBC9 contact. Some of these cones were contacted by a CBC9 dendritic branch that continued past the contact site. Other cones – although contact by one CBC9 – were not contacted by passing dendrites from other CBC9s. The other six cones had at least two invaginating contacts from CBC9s. These originated from two different CBC9s or – in case they originated only from a single CBC9 – at least one of them was formed by a dendritic branch ending at the cone (Figure 1E). We labeled the eight cones that featured only a single CBC9 contact as M-cones (Figure 1—figure supplement 1C), and defined the remaining six candidate cones as S-cones (Figure 1D and Figure 1—figure supplement 1C and D, see Materials and methods). This corresponds to a fraction of 4.8% S-cones (6/124 cones within the dendritic field of at least one CBC9), matching the 3–5% reported in previous studies (Röhlich et al., 1994; Haverkamp et al., 2005).
An alternative scheme for identifying S-cones would have been to classify all cones with at least one invaginating contact from CBC9 as S-cones. This would have resulted in a total of 14 S-cones out of 124 cones (Figure 3—figure supplement 2A) or a fraction of 11.3%. Because this S-cone percentage is much larger than the 3–5% reported earlier (Haverkamp et al., 2005), we consider this scenario as very unlikely (p=0.0037, binomial test, null hypothesis: 5% S-cones, n = 124).
We next developed an automatic method to distinguish contacts likely corresponding to synaptic connections from false contacts. As the tissue in the dataset was stained to enhance cell-surface contrast to enable automatic reconstruction, it is not possible to distinguish between synaptic contacts based on explicit ultrastructural synaptic markers, such as vesicles, synaptic ribbons or postsynaptic densities (see also discussion in Helmstaedter et al., 2013). In contrast to the synaptic contacts in the inner plexiform layer studied by Helmstaedter et al. (Helmstaedter et al., 2013), the highly stereotypical morphology of synapses at photoreceptor axon terminals allowed us to classify the contacts (Haverkamp et al., 2000): The ribbon synapses of the cones are placed exclusively in the presynaptic area at the bottom of the cone pedicles. Here, ON-cone bipolar cells (ON-CBCs) make invaginating contacts, where the dendritic tips reach a few hundred nanometers into the presynaptic area of cone pedicles (Figure 2A) (Dowling and Boycott, 1966). In contrast, OFF-cone BCs (OFF-CBCs) make basal contacts in the same area (Figure 2B). These ‘true’ contacts have to be distinguished from contacts in the periphery or at the (out)sides of the cone pedicle as well as contacts between dendrites and cone telodendria, which can happen, for instance as dendrites pass by (Figure 2C).
In total, we found n = 20,944 contacts in n = 2620 pairs of cones and BCs. We trained a support vector machine (SVM) classifier to distinguish whether or not an individual BC obtains input from a cone (as opposed to classifying each individual contact site, see Materials and methods). To this end, we defined a set of seven features, including contact area, eccentricity and contact height, which allowed distinguishing between potential synaptic contacts and ‘false’ contacts (Figure 2D–F, Figure 2—figure supplement 1). For training of the classifier, we manually labeled a randomly selected set of contacts (n = 50 for OFF-CBCs, n = 108 for ON-CBCs and n = 67 for RBCs). Given the highly stereotypical anatomy of the photoreceptor-BC synapse, labeling performed by an experienced human observer is expected to be very accurate – we here consider therefore the human labels as ‘ground-truth’. We trained separate classifiers for ON-CBCs, OFF-CBCs and RBCs and found that they could reliably distinguish between true and false contacts, with a success rate of ~90% (leave-one-out cross-validation accuracy, Figure 2G–I). Deviations between the labels of the automatic classifier and the human labels did not vary systematically with BC type (see Materials and methods). Such deviations typically occurred when human labels were assigned based on more global structural features of a contact; such more contextual features were not included in the features used for automatic classification (examples of misclassified contacts are shown in Figure 2—figure supplement 2 and Videos 1–3).
We analyzed contacts between CBCs and S- and M-cones in the center of the EM stack where cones were covered by a complete set of all BC types. There was no difference in the number of CBCs contacted by S- and M-cones with 12.2 ± 1.5 CBCs (n = 5 cones, mean ± SEM) for S-cones and 12.2 ± 0.4 CBCs (n = 71 cones) for M-cones, respectively. Similarly, the total number of contact points per cone was almost identical for S- and M-cones with an average of 108 ± 24 per S- and 105 ± 5 per M-cone.
We first studied the convergent connectivity onto the different CBC types and studied how many cones provide input to a single BC of each type (Figure 3A and B). To this end, we classified type 5 CBCs, which had not been further subdivided by Helmstaedter et al. (Helmstaedter et al., 2013), into three types (Figure 3—figure supplement 1, see Materials and methods) in agreement with recent reports (Greene et al., 2016).
Most CBC types were contacted predominantly by M-cones, with an average of 2–6 cones contacting individual CBCs. One exception was the CBC9 that – by our definition of S-cones – received considerable S-cone input. We also detected a few contacts between CBC9s and M-cones; these are a consequence of our restrictive definition of S-cone and originate from those cones for which we found only single CBC9 contacts, such that they were classified as M-cone (see above, Figure 1; see also Figure 3—figure supplement 2).
We next evaluated the divergent connectivity from S- and M-cones to CBCs and studied how many individual BCs of each type were contacted by a single cone (Figure 3C). We found that each M-cone contacted on average a little less than one CBC1, while S-cones contacted almost no CBC1, consistent with previous reports (Breuninger et al., 2011). Conversely, we found that M-cones almost never contacted CBC9s (see above), but S-cones contacted on average two. Both cone types contacted all other CBC types (Figure 3D), with each cone making contact with at least one CBC2, 3B, 4, 5I, 6 and 7. In contrast, some cones did not contact ON-CBC types 5T, 5O, X and 8, such that they were contacted by considerably less than one cone on average.
In addition, we tested the hypothesis that CBCs other than type 1 and 9 unselectively contact all cones within their dendritic field (Wässle et al., 2009). To this end, we compared the number of contacted cones and the number of cones that are in reach of the BC dendrites (Figure 3E–G). OFF-CBCs (types 1–4) contacted on average 65–75% of the cones in their dendritic field, with very similar numbers across types (Figure 3G). In contrast, ON-CBCs showed greater diversity: The connectivity pattern of types 5I, 6 and 7 was similar to that observed in the OFF types (Figure 3G); these cells sampled from the majority of cones within their dendritic field (60–80%). CBC5T, 5O, X and 8, however, contacted less than half of the cones within their dendritic field (Figure 3G), with the lowest fraction contacted by CBCX (~20%). This result is not due to a systematic error in our contact classification: We manually checked volume-reconstructed dendritic trees of the respective types for completeness and frequently found dendrites passing underneath a cone with a distance of 1–3 µm without contacting it (Figure 3—figure supplement 3).
Finally, we studied the contact density along CBC dendrites (Figure 3H and I). To check for systematic variation independent of the absolute size of the CBC dendritic tree, we normalized the cone contact density for the dendritic field size of each CBC type (Figure 3I). Almost all CBC types received input at a very similar location relative to their soma, except for CBCX, which received the majority of inputs closer to the soma than all other types relative to its dendritic field size.
As a control, we also ran the connectivity analysis with the set of S-cones from our alternative, more liberal classification (Figure 3—figure supplement 2B,C). In this analysis, CBC9 was the only color specific BC type whereas all other BC types, including CBC1, contacted both S- and M-cones without preferences (Figure 3—figure supplement 2C). This contradicts the result of a previous analysis based on physiology, which implies that CBC1 does not receive S-cone input (Breuninger et al., 2011).
CBCX had an atypical connectivity pattern compared to other CBC types, so we decided to study its connections in more detail. This BC type has only recently been identified by (Helmstaedter et al., 2013; Shekhar et al., 2016). It has a compact dendritic tree but a relatively wide axonal terminal system that stratifies narrowly at approximately the same depth as CBC5O and 5I do. Interestingly, CBCX seems to sample the cone input very sparsely, with input from only two cones on average, and contacting only about 20% of the cones available in its dendritic field (Figure 3C,D and G). In fact, dendrites of CBCX oftentimes passed underneath cones or even stopped shortly before cone pedicles without making contacts at all (Figure 4A and B). It is unlikely that this resulted from incomplete skeletons for these BCs, as all skeletons were independently verified for this study and corrected where necessary (see Materials and methods).
We re-examined all detected contacts between CBCXs and cones and found that very few of those were ‘classical’ invaginating ON-CBC contacts (3 out of 19 contacts, n = 7 cells, Figure 4B–D). The vast majority were ‘tip’ contacts (16 out of 19 contacts, n = 7 cells), which were similar to basal contacts made by OFF-CBC dendrites (Figure 4B–D). The available data was not conclusive with regards to the question whether these tip contacts of CBCX are smaller than those of OFF-CBCs (median area: 0.05 µm2 for n = 22 CBCX contacts; 0.10 µm2 for n = 23 OFF-CBC contacts, but p=0.17, Wilcoxon ranksum test).
In contrast to the CBCX, the other ON-CBC types made mostly invaginating contacts (71 out of 81 contacts, n = 12 cells, two cells per BC type, Figure 4D), indicating a significant effect of cell type on contact type (GLM with Poisson output distribution, n = 38, interaction: p=3.6x10−7, see Materials and methods). We checked if CBCX receive rod input instead but did not observe any rod contacts (see below). Thus, the CBCX appears to be an ON-CBC with both very sparse and atypical cone contacts similar to those made by OFF-CBCs. Still, based on the axonal stratification depth (Helmstaedter et al., 2013) and recent electrophysiological and functional recordings (Ichinose et al., 2014; Franke et al., 2016) this BC type is most likely an ON-CBC, supported by its mGluR6 expression (Shekhar et al., 2016).
We next analyzed the connectivity between photoreceptors and rod bipolar cells (RBCs) to test the hypothesis that RBCs may contact cones directly (Pang et al., 2010). Three cells labeled as RBC in Helmstaedter et al. (Helmstaedter et al., 2013) were excluded from this analysis, as they did not contact any rods (Figure 5—figure supplement 1). We also found some rods not contacted by any RBC, which is likely due to incomplete tracing of the fine dendritic tips of some RBCs.
In fact, RBCs did not only contact rod spherules but also cone pedicles (Figure 5A,B). These contacts were typical ON-CBC contacts with invaginating dendritic tips into the cone pedicles (Figure 5B). To quantify the cone-to-RBC connectivity in more detail, we counted the number of cones contacted by an individual RBC. While the vast majority (75%) contacted at least one cone, only 25% of all RBCs (n = 141) did not contact any (Figure 5C). However, we did not find a preference of RBCs to connect S- or M-cones (Figure 5D). Conversely, 45% of cones contacted a single RBC, ~35% spread their signal to two to four RBCs, and only 20% of the cones did not make any contact with an RBC (Figure 5E). Our finding provides an anatomical basis to the physiologically postulated direct cone input into a subset of RBCs (Pang et al., 2010). Next, we evaluated whether RBCs contacting only rods or both cone(s) and rods represent two types of RBC, as hypothesized by Pang et al. (2010). However, the two groups of RBCs did not differ regarding the stratification depth of their axonal arbor (Figure 5—figure supplement 2A), number of rod contacts (Figure 5—figure supplement 2B) or potential connectivity to AII amacrine cells (Figure 5—figure supplement 2C), and did not form independent mosaics (Figure 5—figure supplement 2D). In addition, the dendritic field size (116 vs 131 µm², p=0.1, n = 139 RBCs) and the number of dendritic tips (46 vs. 45.5, p=0.8, n = 12 manually counted RBCs) did not differ significantly between the two groups. Therefore, the available anatomical data argue against two types of RBC.
Analogous to the analysis above, we skeletonized and volume rendered a complete set of over 2000 neighboring rod spherules (about 50% of the EM field, Figure 6A, Figure 6—figure supplement 1) and identified rod-to-bipolar cell connections. In addition to the well-described invaginating rod-to-RBC connections (Figure 6B), we also found basal contacts between OFF-CBCs and rods close to the invaginating RBC dendrites (Figure 6C), as described earlier (Hack et al., 1999; Mataruga et al., 2007; Haverkamp et al., 2008; Tsukamoto and Omi, 2014). We did not find any contacts between ON-CBCs and rods (in agreement with Tsukamoto and Omi, 2014; but see Tsukamoto et al., 2007).
A single RBC contacted about 35 rods (Figure 6D), which is slightly more than what was recently reported (~25 rods, Tsukamoto and Omi, 2013). A single rod contacted one or two RBCs, but very rarely no RBC or more than two (Figure 6E). In all cases with two invaginating dendrites, the dendrites belonged to two different RBCs (n = 30 rods). The rods without RBC contacts were mainly located at the border of the reconstructed volume, where we could not recover all RBCs. The number of rods contacting OFF-CBCs was much lower: Whereas CBC1 and CBC2 did not receive considerable rod input, CBC3A, CBC3B and CBC4 were contacted by 5–10 rods, with CBC3B receiving the strongest rod input (Figure 6D).
We analyzed an existing electron microscopy dataset (Helmstaedter et al., 2013) to quantify the connectivity between photoreceptors and bipolar cells. We found interesting violations of established principles of outer retinal connectivity: The newly discovered CBCX (Helmstaedter et al., 2013), likely an ON-CBC (Ichinose et al., 2014; Franke et al., 2016), had unexpectedly few and mostly atypical basal contacts to cones. While CBC types 5T, 5O and 8 also contacted fewer cones than expected from their dendritic field, they exhibited ‘standard’ invaginating synapses. In addition, we provide anatomical evidence that rod and cone pathways are interconnected, showing frequent cone-RBC contacts. The emerging picture of BC types with their input profiles are summarized in Figure 7.
Since the dataset we used was not labeled for synaptic structures, we used automatic classifiers based on structural criteria to identify putative synaptic contacts between BCs and photoreceptors. Due to the highly stereotypical anatomy of the photoreceptor-BC synapse, these criteria allow unambiguous identification of synaptic sites for trained humans (see also Results). For example, we used the proximity of the closest contact to the center of the cone pedicle region as a feature, where presynaptic ribbons have been reported at the ultrastructural level (Dowling and Boycott, 1966; Chun et al., 1996).
The overall accuracy of the classifiers evaluated with human annotated labels was high (~90%). Nevertheless, it is possible that a few contacts were misclassified. Manual quality control, however, revealed no systematic errors. Therefore it is unlikely that classification errors affected our main conclusions. Mismatches between human and classifier labels usually occurred when the human used the global context for the assessment of a contact, knowledge that is not easily transferred into an algorithm (Figure 2—figure supplement 2). Still, human error cannot be ruled out: For instance, contact points labeled by humans as non-synaptic may feature a gap junction and are therefore functional rather than random. For reference, all data including software for classifying and examining BC-cone contacts is available online.
Unfortunately, the retinal location of the EM stack used here is unknown (Helmstaedter et al., 2013); it may originate from the ventral retina, where M-cones co-express S-opsin (Röhlich et al., 1994; Baden et al., 2013) However, as ‘true’ S-cones were shown to be evenly distributed across the retina (Haverkamp et al., 2005), CBC9 connectivity can be used for identification of S-cones independent of location. Nevertheless, it cannot be excluded that opsin co-expression in M-cones in the ventral retina may influence the connectivity patterns between the M-cones and the remaining bipolar cell types.
We found that ON-CBCs 5T, 5O, X and 8 contact fewer cones than expected from the size of their dendritic field. We observed that many of their dendrites passed by the cone pedicles with a distance of 1–3 µm or even ended under a cone pedicle without contacting it (Figure 3—figure supplement 1). This is in agreement with a recent study reporting that CBC8 does not contact all cones within its dendritic field (Dunn and Wong, 2012), but in contrast to earlier studies that concluded that different diffuse BC types receive input from all cones within their dendritic field (Boycott and Wässle, 1991; Wässle et al., 2009). However, a crucial difference with the earlier studies and our study is spatial resolution: Conventional light microscopy can resolve depth with a resolution of several hundreds of nanometers, whereas the EM dataset we used has a resolution of 25 nm, allowing us to more accurately assess whether pre- and postsynaptic structures are in contact with each other.
Could diffusion-based synaptic signaling play a role in ON-CBCs with sparse contacts? ‘Diffusion-based’ synaptic contacts have been proposed for OFF-CBCs (DeVries et al., 2006) and between cones (Szmajda and Devries, 2011). However, although diffusion-based transmission may be present in the outer retina, there is no evidence so far that cone signals reach bipolar cells that neither make invaginating or basal contacts (i.e. with dendrites just passing by a cone pedicle).
As shown above, the CBCX makes the fewest contacts with cones of all ON CBCs. On average, they contacted only about two cones, representing only 20% of the cones within the area of their dendrites. This finding is in agreement with a recent single-cell RNA-seq study, which found that CBCXs feature lower expression levels for metabotropic glutamate receptor mGluR6 (grm6) – the hallmark of ON-BCs – compared to other ON-CBC types (Shekhar et al., 2016). The mouse CBCX is reminiscent of the giant CBC in macaque retina (Joo et al., 2011, Tsukamoto and Omi, 2016) with respect to several features: Also the latter has a very large and sparsely branched dendritic tree and a relatively large axonal arbor that stratifies in the middle of the IPL and contacts only about 50% of the cones in its dendritic field.
In contrast to all other ON-CBCs, we found that the vast majority of CBCX contacts were not invaginating but rather resembled basal OFF-CBC contacts. It is unclear if these tip contacts are indeed functional synaptic sites. This is not the first finding to challenge the traditional view that ON-CBCs form only invaginating and OFF-CBCs only basal synaptic contacts. In the primate fovea, diffuse ON-CBCs (DBs) form basal contacts with foveal cones since almost all invaginating sites are taken by midget bipolar cell dendrites (Calkins et al., 1996). This spatial limitation is less evident in mid-peripheral primate retina. At 3–4 mm eccentricity, diffuse ON-CBCs receive 10% (DB5) to 40% (DB4 and DB6) of their cone input through basal synapses (Hopkins and Boycott, 1996).
Interestingly, CBCX contacts in the IPL also appear to be distinct from those of other BC types: First, the majority of cells contacted by CBCX in the IPL are amacrine cells rather than ganglion cells (Helmstaedter et al., 2013). Second, they form sparse contacts relative to their axon terminal size with comparatively few cells. Thus, the CBCX seems to be an exception, an unusual BC type in many respects in addition to its sparse and atypical connectivity properties in the OPL, reminiscent of a recently described dendrite-less interneuron type that expresses BC-specific genes (Shekhar et al., 2016) and was named GluMI (glutamatergic monopolar interneuron) (Della Santina et al., 2016). It is conceivable that – similar to the GluMI – the CBCX is evolutionary on its way to retracting its dendrites. Alternatively, CBCX develops a tad later than other CBC types and most potential synaptic sites at cone ribbons are already occupied, such that the CBCX can only form few connections – reminiscent of the situation in primate (see above).
We found that cones connect to 75% of RBCs; in many cases, one cone contacted multiple RBCs. In turn, 35% of RBCs received converging input from several cones. This massive cone input via invaginating synapses to RBCs suggests a prominent use of the primary rod pathway (Bloomfield and Dacheux, 2001) during photopic conditions. Consistent with our findings, RBCs can be activated under photopic light conditions (Franke et al., 2016). However, since rods recover at high light levels (Blakemore and Rushton, 1965), the functional significance of cone input to RBC remains unclear. In principle, it is possible that the observed cone-RBC synapses are developmental ‘leftovers’ without physiological relevance, but as the cone-RBC contacts look like standard cone-CBC invaginating synapses, we think that it is more likely that they contribute to RBC activation under photopic conditions, especially in mid-range light intensities where cones are active but rods not yet recovering (Tikidji-Hamburyan et al., 2015).
If this was the case, cone activation of RBCs could indirectly inhibit OFF-CBCs via AII amacrine cells. This suggests that RBCs may contribute to crossover inhibition (Molnar and Werblin, 2007). On the other hand, it is unclear whether gap junctions between AIIs and ON-CBCs are in an open or closed state under light-adapted conditions (Bloomfield et al., 1997; Kuo et al., 2016). With open gap junctions, activating RBCs may boost the signal in ON-CBC axon terminals and therefore enhance contrast (in complement with the OFF-CBC inhibition).
Based on the physiological finding that only a subset of RBCs receive input from cones, Pang et al. (2010) suggested that there may be two distinct RBC types, with the rod-only one having axon terminals ending closer to the ganglion cell layer. Our data do not provide evidence for two RBC types based on the connectivity in the outer retina. This agrees well with recent findings from single-cell RNA-seq experiments, where all RBCs fell into a single genetic cluster with little heterogeneity (Shekhar et al., 2016).
We quantified the number of rods contacting the five OFF-CBC types. Whereas CBC1 and 2 received almost no rod input, we observed flat/basal contacts between rods and types CBC3A, 3B and 4, providing a quantitative confirmation of this finding (Mataruga et al., 2007; Haverkamp et al., 2008; Tsukamoto and Omi, 2014). CBC3A and four received input from ~5 rods in addition to the ~5 cones contacted by them. CBC3B sampled from the same number of cones but was contacted by about twice as many rods. Since these basal contacts between rods and OFF CBCs have been shown to express AMPA receptors (Hack et al., 1999), rods likely provide considerable input to OFF-CBCs, possibly representing a distinct scotopic OFF channel complementing the scotopic ON channel via RBCs. Interestingly, the morphologically similar CBC3A and 3B may obtain their (functional) differences not only from the expression of different ionotropic glutamate receptors (Puller et al., 2013) but also from their connectivity with rods.
Here, we performed a systematic quantitative analysis of the photoreceptor-to-bipolar cell synapse. We showed that there are exceptions to several established principles of outer retinal connectivity. In particular, we found several ON-BC types that contacted only a relatively small fraction of the cones in their dendritic field. We also find that rod and cone pathways already interact strongly in the outer plexiform layer. Whether these are general features of mammalian retinas or evolutionary specializations unique to the mouse remains to be seen.
We used the SBEM dataset e2006 published by (Helmstaedter et al., 2013) for our analysis (http://www.neuro.mpg.de/connectomics). The dataset has a voxel resolution of 16.5×16.5×25 nm with dimensions 114 µm × 80 µm × 132 µm. We performed volume segmentation of the outer plexiform layer (OPL) using the algorithms of (Helmstaedter et al., 2013). The preprocessing of the data consisted of three steps: (i) Segmentation of the image stack, (ii) merging of the segmented regions and (iii) collection of regions into cell volumes based on traced skeletons.
We modified the segmentation algorithm to prevent merging of two segments if the total volume was above a threshold (>50,000 voxels), as sometimes the volumes of two cone pedicles could not be separated with the original algorithm. Although this modification resulted in overall smaller segments, these were collected and correctly assigned to cells based on the skeletons in the last step of the preprocessing.
We identified 163 cone pedicles and created skeletons spanning their volume using the software KNOSSOS ([Helmstaedter et al., 2012], www.knossostool.org,). We typically traced the center of the cone pedicle coarsely and added the individual telodendria for detailed reconstruction. In addition, we traced 2177 rod spherules covering half of the dataset (Figure 6). For our analysis, we used all photoreceptors for which at least 50% of the volume had been reconstructed (resulting in 147 cones and 1799 rods). We used the BC skeletons published by Helmstaedter et al. (2013), with the following exceptions: We completed the dendritic trees of three XBCs (CBCXs), which were incompletely traced in the original dataset. In addition, we discarded three BCs originally classified as RBCs because they were lacking rod contacts as well as the large axonal boutons typical for RBCs (Supp. Figure 6A–C), and one BC classified as a CBC9 because its dendritic field was mostly outside of the data stack (Supp. Figure 6D).
Next, we used the algorithm by (Helmstaedter et al., 2013) to detect and calculate the position and area of 20,944 contact points between cone pedicles and BC dendrites and 7993 contact points between rod spherules and BC dendrites. To simplify the later visual inspection of contacts, we used the reconstructed cell volumes to generate colored overlays for the raw data to highlight the different cells in KNOSSOS.
We detected 169 contacts in 51 pairs of CBC9s and cones. Upon manual inspection, we found a total of 32 invaginating (potentially synaptic) contacts between 6 CBC9s and 14 cone pedicles.
Based on immunocytochemistry, it has been shown that S-cones are contacted by all CBC9 within reach and that CBC9 contacts to S-cones are mostly at the tips of the dendrites (Haverkamp et al., 2005). For all 14 contacted cones, we analyzed the number of invaginating CBC9 contacts, the number of contacting CBC9s, the fraction of CBC9 with dendrites close to the cone that make contact and whether the dendrites end at the cone or continue beyond it (Figure 1E). Based on these criteria, we classified 6 out of these 14 cones as S-cones (see also Figure 1—figure supplement 1). In addition to our main analysis, we present an alternative analysis that considers the case if all 14 cones were counted as S-cones (Figure 3—figure supplement 2).
CBC5s were classified initially based on their connectivity to ganglion cells and amacrine cells into types 5A and 5R, where 5R was a group containing multiple types (Helmstaedter et al., 2013). In addition, some CBC5s could not be classified due to a lack of axonal overlap with the reconstructed ganglion cells of the types used for classification. Considering the separate coverage factors for dendritic and axonal overlap of all CBC5s together (OPL: 3.14, IPL: 2.89), dividing them into three subtypes is conceivable considering the numbers for other CBC types (Table 1). This has already been suggested by (Greene et al., 2016), who divide CBC5s into three subtypes based on axonal density profiles (using a different EM dataset that includes only the inner retina).
We followed the classification approach suggested by Greene et al. (Greene et al., 2016): First, we calculated the densities of both ON- and OFF-starburst amacrine cells (SACs) dendrites along the optical axis. We fitted the peak of these profiles with a surface using bivariate B-splines of third order. Next, we corrected the density profiles of CBC5 axonal trees by mapping the SAC surfaces to parallel planes. We then applied principal component analysis (Figure 3—figure supplement 1A) to obtain a first clustering into three groups by fitting a Gaussian mixture model (GMM) (Bishop, 2006) with three components onto the first three principal components of the axon density profiles. The resulting density profiles of the three clusters matches those found by (Greene et al., 2016) (Figure 3—figure supplement 1B). As we noted a few violations of the postulated tiling of the retina by each type (Seung and Sümbül, 2014), we implemented a heuristic to shift cells to a different cluster or swap pairs of cells optimizing a cost function including both overlap in IPL and OPL as well as the GMM clustering (Figure 3—figure supplement 1):
with the parameter vector of cell , the mixture component cell is assigned to, the mean of the mixture component , the covariance matrix of the mixture component , the Kronecker delta, the area of the dendritic field/axonal tree of cell and the overlap of cell and in the OPL/IPL. The overlap of two cells is calculated as the intersection of the convex hull of the dendritic fields/axonal trees. Likely, our CBCX corresponds to CBC5D from (Shekhar et al., 2016) and CBC5T to CBC5C. Possibly, CBC5I corresponds to CBC5A and CBC5O to CBC5B (see discussion in Greene et al., 2016; Shekhar et al., 2016).
To distinguish potential synaptic contacts between photoreceptors and BCs from accidental contacts, we developed an automatic classification procedure exploiting the stereotypical anatomy of cone-BC synapses (triads, Dowling and Boycott, 1966). First, we grouped all contacts for a specific cone-BC pair, in the following referred to as a contact-set. We obtained a training data set by randomly selecting 10 contact-sets per CBC type and 50 RBC-cone contact-sets. We excluded CBCX from the training data because of their atypical contacts. To increase classifier performance we added 17 additional RBC-cone contact-sets manually classified as invaginating contacts as well as all 48 CBC9-cone contact-sets classified for the S-cone identification. For those contact-sets, we visually inspected each individual contact point in the raw data combined with volume segmentation overlay using KNOSSOS. Then we classified it either as a central basal contact (potentially synaptic) or peripheral contact (e.g. at the side of a cone or contact with telodendria, likely non-synaptic) for OFF-CBCs or as invaginating contact vs. peripheral contact for ON-CBCs and RBCs. Next, we extracted a set of seven parameters for each contact (see Figure 2—figure supplement 1):
Contact area: The total contact area aggregated over all contact points between a BC and a cone
Eccentricity: The distance between the cone center and the closest contact point in the plane perpendicular to the optical axis
Contact height: The distance of the contact point with minimal eccentricity from the bottom of the cone pedicle (measured along the optical axis, normalized by the height of the cone pedicle).
Distance to branch point: Minimal distance between a contact point and the closest branch point, measured along the dendrite
Distance to tip: Minimal distance between a contact point and the closest dendritic tip. A large distance occurs for example for a contact between a passing dendrite and a cone.
Smallest angle between the dendrite and the optical axis at a contact point
Number of contact points between cone and BC
Based on those parameters we trained a support vector machine classifier with radial basis functions (C-SVM) for each OFF-CBC, ON-CBC and RBC cone contact using the Python package scikit-learn. Optimal parameters were determined using leave-one-out cross validation (see Table 2 for scores and error rates). Typically, 0–2 errors for 10 labeled training samples occurred per BC type (three in one case, CBC3A; 7/48 for CBC9).
As the reconstructed rod spherules cover only half of the EM dataset, we restricted the analysis to bipolar cells with their soma position inside this area. To automatically classify the contacts to rods, we followed a similar scheme as for the cones. Again, we grouped the contacts for each pair of BC and rod spherule. As training data, we selected all putative contact sites with CBC1s (n = 5) and CBC2s (n = 32), 20 random contacts to CBC types 3A, 3B and four as well as 100 random contacts to RBCs. Again, we classified these contacts by visual inspection in KNOSSOS using the raw data with a colored segmentation overlay. In addition, we manually inspected all 132 contact points between rod spherules and ON-CBCs, but could not identify a single potential synaptic contact. We trained SVM classifiers for contacts between rods and RBCs/OFF-CBCs using the same parameters as for the contacts to cones. As synaptic contacts between OFF-CBCs and rod spherules are basal contacts situated close to the invaginating RBC contacts, we added the minimum distance to the next (synaptic) RBC contact as an additional classification parameter for OFF-CBCs. As a consequence, we restricted the analysis of OFF-CBC-to-rod contacts to those rods were RBC contacts could be identified (n = 1685). See Table 3 for scores and error rates from the leave-one-out cross-validation.
Error bars in all plots are 95% confidence intervals (CI) calculated as percentiles of the bootstrap distribution obtained via case resampling. In Figure 4D, we used a generalized linear mixed model with Poisson output distribution and fixed effects contact type and cell type and random effect cell identity (R package lme4). The model yielded a significant intercept (z = 8.72, p<0.0001), a significant main effect of cell type (z = 4.11, p=4x10−4), a significant main effect of contact type (z = 2.66, p=0.008) and a significant interaction cell x contact type (z = −5.09, p<3.6x10−7).
Jupyter notebooks and data for reproducing our analysis and main figures are available online at https://github.com/berenslab/pr_bc_connectivity.
Pattern Recognition and Machine LearningNew York: Springer.
Dark adaptation and increment threshold in a rod monochromatThe Journal of Physiology 181:612–628.https://doi.org/10.1113/jphysiol.1965.sp007786
Rod vision: pathways and processing in the mammalian retinaProgress in Retinal and Eye Research 20:351–384.https://doi.org/10.1016/S1350-9462(00)00031-8
Morphological classification of bipolar cells of the primate retinaEuropean Journal of Neuroscience 3:1069–1088.https://doi.org/10.1111/j.1460-9568.1991.tb00043.x
Chromatic bipolar cell pathways in the mouse retinaJournal of Neuroscience 31:6504–6517.https://doi.org/10.1523/JNEUROSCI.0616-11.2011
Organization of the Primate Retina: Electron MicroscopyProceedings of the Royal Society B: Biological Sciences 166:80–111.https://doi.org/10.1098/rspb.1966.0086
Retinal bipolar cells: elementary building blocks of visionNature Reviews Neuroscience 15:507–519.https://doi.org/10.1038/nrn3783
Type 4 OFF cone bipolar cells of the mouse retina express calsenilin and contact cones as well as rodsThe Journal of Comparative Neurology 507:1087–1101.https://doi.org/10.1002/cne.21612
The primordial, blue-cone color system of the mouse retinaJournal of Neuroscience 25:5438–5445.https://doi.org/10.1523/JNEUROSCI.1117-05.2005
High-accuracy neurite reconstruction for high-throughput neuroanatomyNature Neuroscience 14:1081–1088.https://doi.org/10.1038/nn.2868
Synapses between cones and diffuse bipolar cells of a primate retinaJournal of Neurocytology 24:680–694.https://doi.org/10.1007/BF01179818
Roles of ON cone bipolar cell subtypes in temporal coding in the mouse retinaJournal of Neuroscience 34:8761–8771.https://doi.org/10.1523/JNEUROSCI.3965-13.2014
Organization of the outer plexiform layer of the primate retina: electron microscopy of Golgi-impregnated cellsPhilosophical Transactions of the Royal Society B: Biological Sciences 258:261–283.https://doi.org/10.1098/rstb.1970.0036
Type 3a and type 3b OFF cone bipolar cells provide for the alternative rod pathway in the mouse retinaThe Journal of Comparative Neurology 502:1123–1137.https://doi.org/10.1002/cne.21367
Inhibitory feedback shapes bipolar cell responses in the rabbit retinaJournal of Neurophysiology 98:3423–3435.https://doi.org/10.1152/jn.00838.2007
Glutamate spillover between mammalian cone photoreceptorsJournal of Neuroscience 31:13431–13441.https://doi.org/10.1523/JNEUROSCI.2105-11.2011
Retinal output changes qualitatively with every change in ambient illuminanceNature Neuroscience 18:66–74.https://doi.org/10.1038/nn.3891
Functional allocation of synaptic contacts in microcircuits from rods via rod bipolar to AII amacrine cells in the mouse retinaJournal of Comparative Neurology 521:3541–3555.https://doi.org/10.1002/cne.23370
Cone contacts, mosaics, and territories of bipolar cells in the mouse retinaJournal of Neuroscience 29:106–117.https://doi.org/10.1523/JNEUROSCI.4442-08.2009
Fred RiekeReviewing Editor; Howard Hughes Medical Institute, University of Washington, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Connectivity map of bipolar cells and photoreceptors in the mouse retina" for consideration by eLife. Your article has been favorably evaluated by David Van Essen (Senior Editor) and three reviewers, one of whom, Fred Rieke (Reviewer #1), is a member of our Board of Reviewing Editors.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
We all agreed that the paper provides important and interesting data on connectivity in the outer retina, and that the results were generally well presented. We all also agreed on two key points that need to be strengthened in the paper:
1) Separation of false vs. true contacts. The EM dataset used in the paper does not permit unambiguous identification of synapses, and hence the authors rely on an automated procedure to separate true synaptic contacts from other membrane appositions between two cells. This procedure is verified by comparison with results from human annotation. More details of this procedure and its validity are needed. Some specific items that would help are: (1) examples of contacts misidentified by the algorithm and corrected by human observers; (2) a discussion from the beginning of the results of the impact of errors in identification for the conclusions of the paper; and, (3) some discussion of the fact that there is no unambiguous "ground truth" identification available (i.e. that human observers also likely make some errors) and how that could impact the results.
2) Discussion of novel contacts. A strength of the paper is quantification of unexpected synaptic contacts. We felt it would help readers interpret the importance of those results if they were put in the context of physiological signaling (both established and predicted), and in the context of development (e.g. perhaps these contacts play a role in development rather than signaling in the mature retina – see comments from reviewer 2 in particular).
These points and other points are described in more detail in the individual reviews below.
This paper used an existing electron microscopy data set to provide a comprehensive analysis of connectivity in the outer retina in mouse. The results provide direct evidence for several unexpected connectivity features that had been suggested previously but not proven definitively; further, the paper shows that a newly discovered cone bipolar type exhibits properties very distinct from the others. These results are quite interesting and are generally presented clearly. I have a few concerns about the analysis and some suggestions about presentation:
1) Identification of S cones.
S cones were identified based on the assumption that they alone were contacted by type 9 bipolar cells at invaginating contacts made at the end of the bipolar dendrites. These assumptions result in an approximately correct number of S cones given past estimates of their density. Nonetheless, the process seems almost certain to misidentify some cones. It is not entirely clear how important this is for the conclusions of the paper. The section in the Discussion presenting an alternative identification of S cones is quite helpful in this regard, but in my view, comes too late. I would like to see a more thorough analysis of the potential errors in identification and the impact of those errors on the conclusions of the paper. Doing that from the start should help build a reader's confidence. Related to this point, it would help to lay out the criteria very directly (e.g. I was not entirely clear from the second paragraph of the subsection “Identification of S- and M-cones” if contacts only at the end of the dendrites was a hard criterion).
2) Identification of synaptic contacts.
The original dataset does not permit clear identification of synapses, so the raw data the authors have to work with consists of points of close apposition of two cells. Synapses were identified automatically based on extracted features of each close apposition. These identifications were compared with those from manual identification to verify the accuracy of the automated process. At present, it is hard for a reader to evaluate and confirm the accuracy of identification of true contacts vs. false contacts, since this relies on a "ground truth" that itself has some uncertainty at least for the reader. I think it is important to construct as strong an argument as possible that manual inspection and classification indeed provides nearly unambiguous identification of true synapses. As for the point above, I would suggest doing this from the start of the paper. It would be very helpful to have a better sense for the impact of classification errors (which are not so rare from the analysis described in the last paragraph of the subsection “Classification of photoreceptor-BC contacts”) on the conclusions. Related to this point is the possibility of distant contacts or near contacts detecting diffusing transmitter. This is brought up late in the paper (subsection “Sparse contacts between some ON CBC types and cones”, last paragraph), but it would help to deal with the issue earlier (e.g. around subsection “Classification of photoreceptor-BC contacts”, first paragraph). Another point related to this is the uniqueness of identification of basal contacts – e.g. how certain can we be that rod-Off CB basal contacts are real synapses?
On the whole this is a fine analysis of photoreceptor – bipolar cell contacts in the mouse retina based on an existing medium resolution SBFSEM connectome. The key new conclusions are that many rod BCs actually contact a few cones (though I think the significance is developmental rather than physiological) and many OFF cone BCs contact a few or several rods based on subtype.
Despite the impoverished resolution inherent in SEM datasets (no ribbons, vesicles, gap junctions etc.), it is possible to use proximity to suggest possible connectivity. From my experience, there will be few systematic errors in this process, unlike the inner plexiform layer where proximity can rarely be used, if ever.
The findings emphasize some new bipolar cell classes and will be very useful in guiding higher resolution TEM connectomics in process in other labs.
On the whole it is a nicely illustrated paper.
In this elegant and well-executed study, Behrens and colleagues perform a detailed volume reconstruction of previously available and published serial EM data from the mouse retina. The original reconstruction by Helmstaedter and colleagues identified a new bipolar cell type (the CBCX) and focused on reconstructing the inner plexiform layer. The authors of the current study perform a reconstruction of the synaptic connections in the outer plexiform layer. This study offers significant novel findings that add to and also challenge the existing knowledge of outer retina connectivity. Most importantly, they show that there is a potential cross-over of input from cones to rod bipolar cells with the majority of RBCs (~75%) contacting cones. This further challenges the status quo of separate and isolated channels for parallel processing in the vertebrate retina. Furthermore, the authors show that the newly discovered CBCX makes surprisingly few contacts with the cones (~20%) available within the cell's dendritic reach. The work also challenges the assumption that there is equal divergence between one cone and all cone bipolar cell types.
An issue of concern, as the authors themselves point out, is that the tissue is devoid of typical synaptic markers like presynaptic ribbons, synaptic vesicles, and postsynaptic densities, which introduces a challenge when trying to decide whether contacts between pre and postsynaptic cells represent a functioning synapse. The authors' statements about connectivity should be qualified for this reason. The authors describe "true" and "false" synapses, but these judgements are unclear given that the dataset did not contain synaptic information that could be used to verify contacts. Furthermore, the authors did not illustrate "false" contacts, making the interpretation of the "true" contacts ambiguous. To this point, the subsection “Does a ‘contact’ represent a synaptic connection?” of the Discussion mentions contacts and synaptic connections and "errors" in classifying a contact. Can the authors show examples in a supplementary figure of what these human and algorithm errors look like? Such clarification would be useful for establishing common guidelines for contact classification in EM data of the retina.https://doi.org/10.7554/eLife.20041.024
- Thomas Euler
- Philipp Berens
- Philipp Berens
- Thomas Euler
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Helmstaedter et al. (2013) for making their data available. This work was funded by the DFG (EXC 307 and BE 5601/1–1) and the BMBF through the BCCN Tübingen (FKZ 01GQ1002) and the Bernstein Award to PB (FKZ 01GQ1601).
- Fred Rieke, Reviewing Editor, Howard Hughes Medical Institute, University of Washington, United States
© 2016, Behrens et al.