Exp Neurobiol 2020; 29(4): 285-299
Published online August 31, 2020
© The Korean Society for Brain and Neural Sciences
Jungryul Ahn1, Huu Lam Phan2, Seongkwang Cha1, Kyo-in Koo2, Yongseok Yoo3* and Yong Sook Goo1*
1Department of Physiology, Chungbuk National University School of Medicine, Cheongju 28644,
2Department of Biomedical Engineering, University of Ulsan, Ulsan 44610,
3Department of Electronics Engineering, Incheon National University, Incheon 22012, Korea
Correspondence to: *To whom correspondence should be addressed.
Yongseok Yoo, TEL: 82-32-835-8453, FAX: 82-32-835-0774
Yong Sook Goo, TEL: 82-43-261-2870, FAX: 82-43-272-1603
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License
(http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and
reproduction in any medium, provided the original work is properly cited.
Neurons communicate with other neurons in response to environmental changes. Their goal is to transmit information to their targets reliably. A burst, which consists of multiple spikes within a short time interval, plays an essential role in enhancing the reliability of information transmission through synapses. In the visual system, retinal ganglion cells (RGCs), the output neurons of the retina, show bursting activity and transmit retinal information to the lateral geniculate neuron of the thalamus. In this study, to extend our interest to the population level, the burstings of multiple RGCs were simultaneously recorded using a multi-channel recording system. As the first step in network analysis, we focused on investigating the pairwise burst correlation between two RGCs. Furthermore, to assess if the population bursting is preserved across species, we compared the synchronized bursting of RGCs between marmoset monkey (callithrix jacchus), one species of the new world monkeys and mouse (C57BL/6J strain). First, monkey RGCs showed a larger number of spikes within a burst, while the inter-spike interval, burst duration, and inter-burst interval were smaller compared with mouse RGCs. Monkey RGCs showed a strong burst synchronization between RGCs, whereas mouse RGCs showed no correlated burst firing. Monkey RGC pairs showed significantly higher burst synchrony and mutual information than mouse RGC pairs did. Comprehensively, through this study, we emphasize that two species have a different bursting activity of RGCs and different burst synchronization suggesting two species have distinctive retinal processing.
Neurons communicate with other neurons in response to changes in the internal and external environment. Primarily, spiking neurons transmit information to other neurons in the temporal pattern of action potentials called spike train due to its all-or-none nature [1, 2]. Their goal is to transmit information to other neurons faithfully. Accordingly, it is very crucial that this spike train successfully carries a large amount of information without any loss to their targets.
Regarding the strategy of transferring information, using only one spike could be very vulnerable as it is sensitive to noise contamination and requires very accurate timing. If a spike is lost or mistimed during synaptic transmission, it can lead to a fatal loss of information. Alternatively, multiple spike trains can carry even much information without any loss, assuring the robust transmission [3, 4].
A burst, which consists of multiple spikes within a short time interval, plays an essential role in improving the reliability of information transmission through synapses. A single spike of a presynaptic neuron may not evoke the excitatory postsynaptic potential (EPSP), while a burst improves the probability of EPSP in a postsynaptic neuron [5-7]. This synaptic transmission through the burst has been widely observed in various neural systems, including Purkinje cells [8, 9], thalamocortical neurons [10, 11], and pyramidal neurons in hippocampus [12, 13]. In these neural systems, synaptic transmission through the burst can induce synaptic plasticity. For instance, in rat hippocampal pyramidal cells, the burst induced EPSP and long-term potentiation (LTP), whereas the single spike did not induce LTP . This result suggests that burst firing may be associated with the modulation of the neural plasticity.
Unfortunately, our knowledge that bursts increase information transmission is limited to a single synapse. However, since most computations in the nervous system are performed by a network of synapses rather than a single synapse, so studying the burst activity at the network level will improve our understanding of the contribution of bursts to information transmission.
In the visual system, retinal ganglion cells (RGCs), the output neurons of the retina, show bursting activity in response to changes in the internal circuit or external visual environment [14-16]. Mostly, this burst information of the RGCs reaches to lateral geniculate neuron (LGN) of the thalamus, which is known as the relay neuron between RGC and the primary visual cortex in the visual pathway. The bursts of RGCs augment the transmission of retinal signals to LGN and induce its synapse plasticity [17-19]. These findings emphasize the functional role of RGC bursts for information transmission in retinogeniculate communication. However, most studies are still limited to a single unit of RGC.
To extend our interest to the network level, in this study, we recorded the population bursting of RGCs using a multi-channel recording system, and the temporal correlation of bursts among RGCs was investigated through correlation analysis. As the first step in network analysis, we focused on the pairwise burst correlation between two RGCs. Furthermore, to find out if the phenomena of population bursting are preserved across species, we compared the synchronized bursting activity of RGCs between marmoset monkey and mouse. However, marmoset monkey is one species of the new world monkeys, and it is not the most representative model for primates like macaque monkeys.
This study emphasizes that two species have a different bursting activity of RGCs, and different burst synchronization. In particular, mouse RGCs hardly show synchronized bursting activity, but monkey RGCs share their bursting activity with nearby RGCs, suggesting two species have distinctive retinal processing.
Male mice at postnatal day 56 and above (>P56) (C57BL/6J strain; The Jackson Lab., ME, USA) were used for this study (n=6). Institutional animal care committee of Chungbuk National University approved animal use protocols (approval number: CBNUA-1172-18-02). The mice were anesthetized with intramuscular injection of 30 mg/kg tiletamine-zolazepam hydroxide (Zoletil 50; Virbac, Sao Paulo, Brazil), 10 mg/kg zylazine hydrochloride (Rumpun; Bayer Korea, Seoul, South Korea), and 5000 IU heparine sodium (Heparin; JW Pharmaceutical Corp., Seoul, South Korea). Monkey eyeballs were obtained from marmoset monkeys (callithrix jacchus, n=4) used by the Laboratory Animal Center of the OSONG KBIO HEALTH, following institutional guidelines for the care and use of animals. In the process of transportation of monkey eyeballs after enucleation, CO2-independent medium, well known as transport buffer, was used to avoid tissue damage caused by pH changes in the hypoxic buffer (CO2-Independent Medium; Thermo Fisher Scientific, MA, USA). Four retinas at postnatal year two and higher were used. The detailed procedures for the preparation of the ex-vivo retinal patch were described in previous studies [20, 21]. Briefly, the eyeball was enucleated, and then the retina was isolated from the sclera and the retinal pigment epithelium (RPE) and cut into a 2×2 mm2 patch. The retinal patch was prepared under illumination of 4.3 nW/cm2 in an artificial cerebrospinal fluid (ACSF) solution (124 mM NaCl, 10 mM Glucose, 1.15 mM KH2PO4, 25 mM NaHCO3, 1.15 mM MgSO4, 2.5 mM CaCl2, and 5 mM KCl) bubbled with 95% O2 and 5% CO2 to maintain a pH of 7.3~7.4 and a temperature of 32°C. The isolated retina was mounted the retinal ganglion cell (RGC) layer down on a planar multi-electrode array (MEA), and it was continuously perfused with oxygenated solution (flow rate: 1~3 ml/min) during the experiment.
The data acquisition system (MEA60 system; Multichannel Systems GmbH, Reutlingen, Germany) included a planar 64-channel perforated MEA (60pMEA200/30iR), an amplifier (MEA1060), temperature control units (TC01), data acquisition hardware (Mc_Card) and software (Mc_Rack). The MEA contained 64 circular electrodes in an 8×8 grid layout with electrode diameters of 30 µm and inter-electrode distances of 200 µm (Fig. 1A). The electrodes were coated with porous titanium nitride (TiN) and were embedded in a perforated polyimide foil, which allows sufficient oxygen and nutrient supply to the retina . Multi-electrode recordings of the retinal activity were obtained from 59 electrodes except for one reference electrode and four inactive electrodes, with a bandwidth ranging from 1 to 3,000 Hz at a gain of 1,200. The data sampling rate was 25 kHz per electrode. From the raw waveform of retinal recording, RGC spikes were isolated using a 100 Hz high-pass filter. The threshold for spike detection was set to four times the standard deviation of the background noise. The recording data were processed with spike sorting software (Offline SorterTM; Plexon Inc., TX, USA) to separate multiunit activities containing different spike waveforms into individual cell units using principal component analysis . Data analysis was processed utilizing a commercial analyzing software (NeuroExplorerTM; Nex Technologies, CO, USA) and a custom-made MatlabTM (MathWorks, MA, USA) code.
For burst analysis, we identified the burst quantitatively using the Poisson surprise algorithm [24, 25]. First, the mean of inter-spike interval (ISI) of the whole spike train,
Starting from the initial burst of three spikes, the following spike is included in the burst if the probability of the null hypothesis (independent spikes) decreases by adding the spike of question. Specifically, the null hypothesis is that
When the probability of kth consecutive spikes (
Finally, any bursts with
The burst-correlation of two simultaneously recorded burst trains between RGCs can be quantified using the cross-correlogram, which displays the amount of synchronized burst firing between two cells as a function of the difference of the median time of each burst. First, we binned burst trains from RGC pairs to generate a binary burst count for each RGC as a function of time. Next, to quantify the strength of a correlation between two RGCs, we used the burst synchrony index (BSI). It is defined as the ratio between the probability of synchronized burst firing (i.e., two cells fire together during a time bin) and the probability expected from statistically independent firing . The BSI was calculated using the following equation:
where 1 denotes the occurrence of a burst. The BSI value of two independent burst trains is zero, and a greater BSI value indicates higher synchrnony.
The time interval for calculating probabilities in Eq. 2 was chosen to be 100 ms, considering two components (the deviation of median times for each burst and IBI values). For burst analysis, average IBI values were around 200 ms from both two species. Thus, the time interval should be shorter than 200 ms. The BSI was calculated and plotted for the inter-electrode distances ranging from the nearest distance of 200 µm to a maximum of 1,400 µm by step of 200 µm.
We measured mutual information (MI) of bursting of two RGCs in bits, which quantifies the amount of information obtained about the burst of an RGC through observation of the burst of the other RGC. In general, the MI of two variables
For the control experiment, MI analysis was further performed assuming two simulated conditions of independent firing and perfectly synchronized firing between RGCs. First, to mimic the independent firing state, a surrogate burst train was created by shuffling the burst train of the reference cell (A). The shuffled train (A') was obtained by randomly rearranging the original burst train of the reference cell, destroying any relationship between the origin and the surrogate trains. Then, the MI between the reference cell and the surrogate cell (A↔A') was calculated. Second, for perfectly synchronized firing, a duplicate train of the reference cell was created and MI analysis was performed (A↔A). Each of the shuffled condition and the perfectly synchronized condition represents the baseline and the upper bound of the mutual information between RGCs on pairwise condition.
We applied light stimulation to the retina to functionally classify RGCs and found the receptive fields (RFs) of RGCs. Briefly, we classified RGCs into two types, considering the spatiotemporal pattern of RGC RFs, where ON RGCs are considered as cells displaying a response when the flash is ‘on’, while OFF RGCs show a response when the flash is ‘off’. Visual stimuli were generated by custom-made software, written in MatlabTM with Psychtoolbox [27, 28], and presented using a digital light processing (DLP) projector (ep7122; Hewlett-Packard, CA, USA) with a maximum resolution of 1,400×1,050 pixels and a refresh rate of 60 Hz. The visual stimulus was projected onto the photoreceptor layer with a size of 1.7×1.7 mm2 through the water dipping 20× objective lens (UMPLFLN20XW; Olympus Corp., Tokyo, Japan), with a beamsplitter (32-506; Edmund Optics Inc., NJ, USA) inserted into the light path for monitoring. Light stimuli were attenuated using a neutral density filter (NE220B; Thorlabs Inc., NJ, USA). For characterization of RFs of RGCs, we applied a spatiotemporal white noise stimulation, temporally updated at a frame rate of 30 Hz, and spatially arranged in a checkerboard layout with a pixel width of 67 µm (Fig. 1C). Each pixel was independently modulated according to the Bernoulli distribution with equal probabilities for +1 (white) and -1 (dark) (100% contrast and mean intensity of 3.8 µW/cm2). Light intensity corresponding to photopic vision was measured with a calibrated radiometer (ILT-5000; International Light Technologies, MA, USA). This white noise stimulation was applied for 15 min.
For an estimation of RGC RFs, the spatial maps of spike-triggered average (STA) were manipulated (Fig. 1C). Our previous study showed spatiotemporal patterns of mouse RFs through STA analysis . Briefly, for STA, spikes from each RGC were counted to produce a train of spike counts binned at 33.3 ms (30 Hz). The STA was obtained by averaging the stimuli preceding a spike with the time window of 400 ms.
All results shown here, including burst parameters, BSI, and MI were represented by averaging for analyzed neurons. Their error bars were indicated with ±standard error of the mean (SEM). In Fig. 3, the BSI was calculated for the inter-electrode distances. Their curves were plotted with the exponential fit. Significant difference between two species in a certain distance was verified with paired T-test (⁎: p<0.05, ⁎⁎: p<0.01, ⁎⁎⁎: p<0.001). To determine the statistical difference of BSIs among inter-electrode distances, ANOVA analysis (p<0.05) for each mouse and monkey retina was performed with post hoc tests of Turkey, Duncan, and Scheffe.
For burst analysis, we used 601 RGCs from 6 mouse retinas and 436 RGCs from 4 marmoset monkey retinas. First, we investigated the bursting patterns of RGCs from two species. Fig. 2A showed each representative burst from mouse and monkey retina. Compared with the mouse, the bursts of monkey RGCs showed remarkable bursting properties such as more spikes in a burst (Fig. 2A). Therefore, between mouse and monkey retinas, we compared burst parameters such as spike counts in a burst, inter-spike interval (ISI) in a burst, burst duration, burst frequency, and inter-burst interval (IBI) (Fig. 2B). For spike counts in a burst and burst frequency, monkey RGCs showed significantly more spike counts (4.21±0.01) within a burst and burst frequency (4.46±0.05 Hz) compared with those of mouse RGCs (3.67±0.01) and (2.42±0.02 Hz), respectively. In monkey RGCs, averaged inter-spike interval in a burst (35.56±0.29 ms) and burst duration (114.20±1.56 ms) are significantly shorter than in mouse RGCs (ISI: 51.67±0.26 ms, burst duration: 138.40±1.36 ms). For IBI, monkey RGCs showed shorter IBIs (222.97±1.06 ms) compared with mouse RGCs (272.57±3.33 ms). Overall, the bursts of monkey RGCs have more spikes over a short period than those of mouse RGCs, suggesting that compared with mouse RGC's bursts, monkey bursts could deliver retinal activity to the higher-order LGN neuron more reliably.
We investigated the synchronized bursting between RGC pairs from two species. Fig. 3A showed each representative RGC bursts in monkey and mouse retinas. Monkey RGCs showed a strong burst correlation between two RGCs compared with mouse RGCs. In the raster plot of two RGCs, the bursts of monkey RGCs were relatively well-synchronized in time. This synchronized pattern was well represented as a significant single peak through burst cross-correlogram, in which a single peak appeared within the time ranges of 100 ms from zero interval. On the other hand, mouse RGCs showed no correlated bursting between two RGCs, even with the nearest spacing of 200 µm. No significant peak in burst cross-correlogram suggests two RGCs have little correlated bursting.
To assess the strength of the burst correlation between RGC pairs, BSI was calculated for each RGC pair recorded on separate channels of MEA. For BSI calculation, we included only one RGC for each channel of MEA to exclude possible effects due to different RGC density or RF diameter across species. We used 193 RGCs (out of total 601 RGCs) from 6 mouse retinas and 156 RGCs (out of 436 RGCs) from 4 monkey retinas. More details are written in discussion section. Fig. 3B showed the BSI as a function of inter-electrode distance. A significant difference in BSI between two species appeared only for RGC pairs with 200 µm spacing (***: p<0.001). The BSI of monkey RGCs decreased with the increment of the distance between RGC pairs, while mouse RGCs showed a near-zero BSIs regardless of distances between RGCs. Noticeably, the BSI of monkey RGC pairs with 200 µm spacing was found to deviate from the fitting curve extrapolated only at the other distances.
To find out the statistical difference of BSIs among inter-electrode distances, we performed ANOVA analysis for mouse and monkey retina, respectively. In monkey retina, two statistically different groups (group 1: 200 µm, group 2: 400~1,400 µm, p<0.05) were found, while the mouse retina showed no significant differences along with the inter-electrode distances (p>0.05).
How much information a pair of RGCs share through pairwise synchronized bursting? To answer this question, we performed MI analysis on synchronized bursting for a pair of RGCs. MI analysis was performed for RGC pairs with 200 µm spacing, since a significant difference in BSI between two species appeared only for RGCs with 200 µm apart (Fig. 3B).
As described in the method section for mutual information, we made shuffled train and perfect sync train. Each of the shuffled- and the perfect conditions represents the baseline and the upper bound of the mutual information between RGCs on pairwise condition. The probability distribution of four possible firing patterns between two RGCs in shuffled, pairwise, and perfect sync train was observed (Fig. 4A). Across six histograms, the probability of (0,0) firing pattern was always dominant. Also, each structure of the shuffled histogram and perfectly synchronized histogram is similar between monkey and mouse. The pairwise histogram of monkey RGCs showed an increase in probability for the (1,1) pattern and a decrease in probability for the (0,1), (1,0) patterns compared with the pairwise histogram of mouse RGCs. The pairwise histogram of mouse RGCs for the (1,1) pattern appears to be very similar to the shuffled histogram, which is the firing pattern distribution expected from independent firing. Based on these probability histograms, the following MIs were calculated respectively: For monkey,
For statistical analysis of MIs, 6726 pairs of monkey RGCs and 8482 pairs of mouse RGCs with 200 µm spacing were used (Fig. 4B). Each group of the shuffled and perfect sync did not show any differences between two species (
We showed a comparison of spontaneous bursting activity across species (Figs. 2, 3). However, a comparison of bursting activity for different physiological cell types was not shown. For this purpose, we performed a functional classification of RGCs using light stimulus. RGCs were classified into two types, ON and OFF RGCs (Detailed information can be found in the Method section, light stimulation and Spike-triggered average and receptive field estimation).
Out of six mouse retinas, 153 ON- and 174 OFF-RGCs were harvested. Fig. 5A shows the mosaics of receptive fields (RFs) of RGCs in a representative mouse retinal patch. Exemplary spatial and temporal RF profiles of RGCs are shown on the bottom side of the RF mosaic figures. Each spatial map of RFs was obtained from a specific recording channel (electrode number of MEA for ON cell in the mouse retina: 54, OFF cell in the mouse retina: 45). The ON RGCs shows a cluster of pixels with white color around the selected electrode, which has a positive contrast temporally. On the other hand, the OFF RGCs shows a negative contrast with black pixels. Each RF size of mouse ON and OFF RGCs was 259.0±63.9 µm and 276.6±67.0 µm, respectively. Additionally, there was no significant difference in RF size between mouse ON and OFF RGCs (p>0.05).
Based on the classified ON and OFF RGCs, especially in the mouse retina, we investigated whether there is a difference in spontaneous bursting activity for each cell type. Fig. 5B and 5C show no significant difference between ON and OFF RGCs at all values, including bursting parameters and burst synchronization.
It was impossible to compare the bursting activity between ON and OFF RGCs in monkeys since we only harvested three ON RGCs without any OFF RGCs out of monkey retinas. The difficulty of harvesting the light-sensitive RGCs in monkey retina will be mentioned in the discussion section.
In the nervous system, the coding scheme relates to how information is reliably transmitted through relays to target neurons. In the case of the visual system, it is associated with transmission through neural synaptic relays between RGCs and LGNs. Bursts are well known for their role in improving the reliability of information transfer (i.e., synaptic transmission) [3, 4]. While a single spike may not pass through a synapse, multiple spikes in a burst will cross the synapse more reliably, improving the possibility of eliciting postsynaptic spikes. Besides, some studies reported that the correlation strength between neurons could be enhanced through the bursts rather than individual spikes [29, 30]. The bursting activity can efficiently enhance the transfer of information either by storing information in the burst or by strengthening the connection with other neurons.
Our results showed that mouse RGCs are not preferential to bursting activity compared with marmoset monkey RGCs (Fig. 2). Given the low spike density, the importance of burst activity in the mouse retina is ambiguous. On the other hand, the burst statistics of monkey RGCs highlight that monkey RGCs exhibit higher spike counts in a burst. These findings indicate monkey RGCs transmit spike information much densely through burst carrier into the higher-order visual system than mouse RGCs do. Comprehensively, our results emphasize RGCs of mouse and monkey retinas have different bursting patterns, suggesting two species have distinctive retinal processing.
The main reasons for different bursting patterns between two species can be explained mainly by two factors, such as (1) Intrinsic RGC bursting created by intrinsic mechanism without any synaptic input , and (2) Input-driven bursting by presynaptic neurons of RGCs such as bipolar cells .
First, for intrinsic bursting, there are many discovered mechanisms of intrinsic bursting, including voltage-[33-35] and Ca2+- gated ion channels. Typically, T-type voltage-gated calcium channel (T-type VGCC; transient opening calcium channels called Cav3) is thought to be involved in the bursting process. T-type VGCC allows the influx of calcium ions. Thus, T-type VGCC induces more frequent depolarization of membrane potential and causes high-frequency firing . Due to its transient activation properties, T-type VGCC plays a role in the rhythmic bursts. In RGCs, T-type VGCC is known to be involved in the maintenance of intracellular Ca2+ concentration and the generation of RGC bursting . T-type VGCC in RGCs has been identified in various species including the salamander , rat [31, 40, 41], turtle , and cat .
As another candidate channel for the intrinsic bursting process, Ca2+-activated K+ channel (KCa) has been reported to be essential for the modulation of bursting activity [44-46]. In particular, KCa is involved in the control of bursts by affecting the hyperpolarization stage of the spike [47-50]. Activation of this channel leads to a decrease of sustained burst firing by increasing the spike threshold time. Blocking of KCa with charybdotoxin results in burst firing by increasing the depolarization of the membrane potentials. Similarly, in RGCs, the KCa channel contributes to the modulation of burst firing in RGCs , and blocking the KCa channel increases the burst firing rate of RGCs . These observations indicate that the modulation of KCa is deeply involved in the spontaneous bursting of RGCs.
Comprehensively, the adjustment between excitation through T-type VGCCs and inhibition through KCa leads to the intrinsic bursting of RGCs. The previous report in intact mouse retina showed that T-type VGCC amplitudes decreased from embryonic through postnatal stages . Unfortunately, however, for the best knowledge of authors, there has been no study for the difference of expression levels of these channels according to the developmental stage on monkey RGC. However, we believe that the different channel density of these burst-associated channels in RGCs may provide a reasonable guess for understanding the distinct bursting activity between two species.
Second, for input-driven bursting by presynaptic neurons of RGC, a difference between rodent retina and primate retina has been reported recently . Using immunohistochemistry in macaque and human retinas, Puthussery et al. showed that voltage-gated sodium channel (NaV1.1) known to be essential for the spike generation, was observed in an axon initial segment (AIS) region of bipolar cells. In other vertebrates, the NaV1.1 was not seen in the AIS of bipolar cells. In rat and mouse retinas, NaV1.1 was found in AISs of the amacrine cell but not in bipolar cell axons. Thus, the high-density expression of NaV1.1 in bipolar cell axons seems to be a unique characteristic of the primate retina. Here, NaV1.1 is not a burst-associated ion channel, so it might not be directly related to the bursting activity of RGCs. However, it can be expected that different presynaptic inputs generated by channel specialization in bipolar cells will result in the different bursting activity of RGCs between two species.
In this study, we showed comparisons of synchronized burst firing of RGCs in mouse and monkey retinas. For burst firing, we found that mouse RGCs showed little correlated firing, while monkey RGCs showed a strong correlation of bursts between RGC pairs.
Burst synchrony index (BSI), the strength of correlated firing between monkey RGCs decreased with the increment of the distance between RGC pairs. There are significant differences in BSI between mouse and monkey RGCs with a distance of 200 µm. In mouse retina, there is no significant difference in the synchronized bursting between RGC pairs irrespective of inter-electrode distances, even for the nearest spacing of 200 µm.
Since we used adult mice (>P56), our finding may not be generalized that mouse RGCs do not show the synchronized bursts throughout their lives. In fact, during the early development of visual circuits before eye-opening, mouse RGCs show strongly correlated bursts called retinal waves in the retina [53, 54]. Their correlated activity works for the refinement of immature retinal circuits and the retinal projection between the retina and higher-order visual system, i.e., LGN and visual cortex. Burst correlation between RGCs appeared most strongly in early postnatal days (P9), became weakened during retinal development, and eventually reached to little synchronization at P42 close to adult age (P56) . In conclusion, at infant age, the burst correlation contributes to the pruning of neural circuits, whereas it disappears at adult age. Thus, here in our study, as shown in Fig. 3, mice over P56 showed no correlated firing between a pair of RGCs and appeared to be independent bursts of each RGC. Unfortunately, for comparison, there have been no monkey studies related to burst correlation according to retinal development in our best knowledge. This issue is to be studied further in the future.
The reason why two species have different burst synchrony can be considered in two different aspects related to RGC: (1) Difference in RGC density, and (2) Different RF size of RGC. First, for RGC density, previous studies showed that the RGC density in primate fovea is higher than that of mouse RGCs [55-57]. The higher the RGC density, the shorter the average distance between RGCs. Second, for RF size, the average RF diameter of monkey RGCs is about 125 µm in range of 50 to 200 µm [58, 59]. In contrast, the average RF diameter of mouse RGCs is about 250 µm ranging from 100 to 400 µm [60-62], which mostly agrees with our result of Fig 5A (RF diameter of mouse ON and OFF RGCs: 259.0±63.9 µm and 276.6±67.0 µm). Thus, monkey RGCs showed smaller RF diameter than mouse RGCs on average. In our study, since we always compare the burst synchrony between two RGCs recorded on a separate channel of MEA, at least 200 µm apart, the RF overlap may happen in the mouse, not the monkey retina. If there is RF overlap, it might increase the burst synchrony in the mouse, not the monkey retina. But we found higher burst synchrony in monkey retina. Therefore, we believe that the difference in burst synchrony is not due to different RGC density or RF diameter between the mouse and the monkey retina.
The other possibility for the inter-species difference in burst synchrony could be found in the different expression levels of gap junction channels. For spontaneous bursting, synchronized bursts are thought to be mainly generated through direct connections via a gap-junction channel called connexin-36 between RGCs . The difference in the expression level of the gap junction channel between monkey and mouse would be one good candidate to explain the difference in synchronized bursts. But in our best knowledge, there has been no report for direct comparison of the expression level of the gap junction channel between mouse and monkey retina.
We found that while mouse RGCs have little information of 2.0×10-5, monkey RGCs share around 2.7×10-3 of information with one neighboring RGC through pairwise synchronized bursting. One can question the physiological meaning of such a burst correlation. At first glance, burst correlation may not seem to be efficient for transmission due to the critical problem of information redundancy. What are the benefits of burst correlation between neurons at the expense of redundancy?
Regarding population bursting, there can be two main benefits: (1) Compensation. Timely synchronized bursts would improve the reliability of information through shared information. During synaptic transmission, if one synapse fails to transmit information, the other synapse could serve as a complement. Thus, the success of information transmission would be guaranteed. (2) Amplification. The synchronized bursts would amplify the signal input by synaptic integration of multiple RGCs. It helps the EPSP in a postsynaptic of LGN to exceed the threshold potential, making it easier for LGN to fire output spike or burst. Besides, based on our knowledge that the RGC burst induces single synapse plasticity of LGN [17-19], we could assume that the synchronized population bursting of multiple synapses will result in more dynamic changes in synaptic plasticity than by a single synapse.
As a result, mouse RGCs are likely to transmit information independently, suggesting that mouse RGCs select a coding scheme to minimize the redundancy of information. On the other hand, monkey RGCs seem to transmit information through mutual dependence. This finding suggests that the higher the evolutionary stage of the species, the more burst synchronization may be required for stable transmission and dynamic synaptic changes.
The next question is that how many RGCs are synchronized, that is, what is the optimal structure of the primate retinal network to satisfy both reliable information transmission and low redundancy. Greater synchrony between RGCs would enhance reliable information transmission to the LGN. However, if the synchrony between RGCs is too high, the amount of transmitted information would be lower. In future works, we would like to investigate more complex networks (i.e., burst correlation among three or more neurons). By this, we could answer how many neurons are needed to maximize mutual information, which will improve our understanding of the functional role of the population bursting for information transmission.
Harvesting light-sensitive RGCs in monkey retina was problematic due to some constraints. It usually takes at least about 40 minutes to transport monkey eyeballs after enucleation from OSONG KBIO HEALTH to our lab. During the retinal preparation, we could often observe the retinal pigment epithelium (RPE) layer, supporting essential nutrients to the retina, was detached severely from the retina. RPE detachment from the retina can induce severe damage to the nearest photoreceptors (PRs) rather than the relatively distant RGCs. In some experiments of RPE detachment, after removal of RPE, the visual pathway from PR to RGC was functionally inactive while RGC's spontaneous firing was well maintained and robust [64-66]. To improve the tissue quality during the time window from enucleation to recording, we used a CO2-independent medium well known as transport buffer to avoid tissue damage caused by pH changes of the buffer solution. However, despite these efforts, there was substantial tissue damage. A fundamental modification of the experimental procedure is required, such as shorter transfer time, to harvest light-responsive RGCs reliably.
Here, we investigated the bursting patterns and burst synchronization of RGC bursts without light stimulus across species. As a next step, we would like to examine the bursting activity of light-stimulated RGCs across species. Many studies address that the burst can code different stimulus features and expand the additional coding space using the number of spikes and temporal patterns (i.e., ISI) within a burst [67-70]. Recently, some studies showed that thalamocortical relay (TCR) neurons or pyramidal neurons use a multiplexed neural code, which distinguishes multiple input streams through burst parameters, including the burst duration and ISI within a burst [71, 72]. Accordingly, burst parameters, including the number of spikes and ISIs in a burst, and burst duration, etc., could be employed to discriminate the multiple stimulus features about what and when of the stimuli.
In this context, with a visual stimulus, the bursting activity may differ from spontaneous bursting. For the additional coding space, some burst parameters, including spike count and burst duration, can increase further compared with spontaneous bursting. Besides, the increase of burst duration may induce the increase of IBI for full recovery during the quiescence period before the neuron can burst again.
It will be essential to investigate the change in burst synchronization with a visual stimulus. For spontaneous bursting, synchronized bursts are mainly generated through direct connections by a gap-junction channel called connexin-36 between RGCs . With a visual stimulus, the light activates a series of retinal circuits, including photoreceptors and bipolar cells. Then, light-driven presynaptic inputs induce RGC firing. Network-mediated input through a chemical synapse by a bipolar cell relay to RGCs will further increase burst synchronization. This network-mediated input to RGCs may augment the burst correlation in mouse retina as well as in monkey retina. It would be interesting to observe if the species difference in the burst synchrony index will be preserved even with visual stimuli.Need for
Our study showed an
Our future work is to investigate how correlated bursting activity contributes to the changes in LGN. To do this, we must record RGC and LGN signals simultaneously. Some studies showed that
We thank the Laboratory Animal Center of the OSONG KBIO HEALTH for providing us the marmoset monkey eyeballs. This research was supported in part by the Bio & Medical Technology Development Program (NRF-2017M3A9E2056460 & 2017M3A9E2062707), the Basic Science Research Program (NRF-2018R1A2B6003917) through the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (MSIP), and an Incheon National University Research Grant in 2017. This work was financially supported by the research year of Chungbuk National University in 2018, in part.