CC BY
44
Functional Connectivity of Nei Hippocampal Culture Grown on Mici
DOI: 10.17691/stm2017.9.2.07 Received June 16, 2016
A.A. Gladkov, Junior Researcher, Neuroengineering Department, Center of Tran Researcher, Molecular and Cellular Technology Department, Central Research Labo V.N. Kolpakov, Laboratory Technician, Neuroengineering Department, Center of Trans Y.I. Pigareva, Junior Researcher, Neuroengineering Department, Center of Translational V.B. Kazantsev, DSc, Professor, Vice-Rector for Research and Innovations1; I.V. Mukhina, DSc, Professor, Head of the Central Research Laboratory2; Head of the Departme of Normal Physiology named after N.Y. Belenkov2; Professor, Neurotechnology Department, Instituti of Biology and Biomedicine1; Head of the Center of Translational Technology1 A.S. Pimashkin, PhD, Senior Teacher, Neurotechnology Department1
1Lobachevsky State University of Nizhni Novgorod, 23 Prospekt Gagarina, Nizhny Novgorod, 603950, Russian Federation;
2Nizhny Novgorod State Medical Academy, 10/1 Minin and Pozharsky Square, Nizhny Novgorod, 603005, Russian Federation
Dissociated neural cultures are used as convenient experimental models to study basic mechanisms of brain signal processing, memory, learning and synaptic plasticity in neuronal networks. Evaluation of short-term and long-term memory in hippocampal cultures requires simultaneous multisite recording of signals and bioelectrical stimulation.
The aim of the investigation was to reveal the characteristic features of neuronal network plasticity in the model of primary hippocampal cell cultures, to study the stability of spontaneous and stimulus-induced changes in spiking patterns of cultures.
Materials and Methods. Long-term changes in functional connectivity characteristics of spikes propagating in the neural network of hippocampal cultures grown on microelectrode arrays were evaluated. It was investigated how low-frequency stimulation (0.1-0.5 Hz) lying in the frequency band of spontaneous bursting activity altered functional connections in the network at different time scales.
Results. Low-frequency stimulation of hippocampal cultures grown on microelectrode arrays was found to induce reconfiguration of the network connectivity. This effect could be preserved during tens of minutes. On the time scale of hours, stimulation-induced connectivity pattern disappeared due to spontaneous changes.
Key words: microelectrode array; electrical stimulation of neurons in vitro; hippocampal culture; functional connectivity of neural networks; synaptic plasticity.
Dissociated neural cultures grown on microelectrode arrays are widely used as convenient models to study the fundamental mechanisms of brain signal processing, memory, learning and synaptic plasticity in neuronal networks [1-4]. Low-frequency stimulation of neuronal cultures is known to evoke network bursts whose patterns depend on the location of the stimulating electrode [5, 6]. The induced reverberating bursts play a key role in changing functional connectivity of neural networks [7-11]. Cultures of dissociated neural cells are known to display dominant spontaneous activity patterns that enhance the existing connectivity [12]. Low-frequency stimulation by one electrode induces generation of network responses organized as reverberating activity bursts different from the dominant spontaneous activity pattern. Such stimulation alters signal propagation pathways and induces connectivity changes [12]. Several research groups have shown that
low-frequency stimulation may induce network plasticity [13-15].
However, spontaneous electrical activity in cultured neural networks demonstrates a high level of variability in the functioning of mature cultures for more than 30 days in vitro [16]. This variability can be associated with spontaneous synaptic plasticity changes due to the presence of cyclic connections in the neural network. The spontaneous changes of functional connectivity in the network may interfere with stimulus-induced changes. To reveal stimulus-induced synaptic connectivity, it is important to measure the dynamics of spontaneous changes in the network.
The aim of the investigation was to evaluate long-term changes in functional characteristics of spikes propagating in the neural network of hippocampal cultures grown on microelectrode array.
In this study we performed long-term experiments
For contacts: Alexander S. Pimashkin, e-mail: [email protected]
Connectivity of Neural Network in Dissociated Hippocampal Culture Grown on Microelectrode Array СТМ j 2017 — vol. 9, No.2 61
recording spontaneous activity of the neural network. Low-frequency voltage stimulation was applied through 12 electrodes to induce bursting responses in certain locations of the network. We hypothesized that changes in functional connectivity of the network were induced by synaptic plasticity. Functional connectivity was estimated by cross-correlation analysis of spiking patterns of spontaneous network activity recorded through different microelectrodes (by analogy with [17]). Connectivity changes induced by electrical stimulation were found on the time scale of 20 min. Besides, low-frequency stimulation changed burst initiation pattern (a propagation pathway of the signal inducing burst initiation) and spiking rate of spontaneous network activity.
Materials and Methods
Cell culture. Hippocampal cells were dissociated from embryonic mice (E18) and plated on 60-electrode microelectrode arrays (MEAs) (Multi Channel Systems MCS, GmbH, Germany) at final density of approximately 9,000-10,000 cells/mm2 [18]. One half of the medium was replaced by new one every 2 days. The cells were cultured under constant conditions at 35.5°C, 5% CO2 and 95% air humidity in cell culture incubator (MCO-18AIC, Sanyo, Japan). Experiments were performed when the cultures were 3-5 weeks in vitro.
Recording and electrical stimulation of cultured cells. Signals were recorded from microelectrode array with 59 planar electrodes of 30 |jm diameter and 200 |jm inter-electrode distance and 1 reference electrode. Spikes were detected using threshold estimated for each channel separately based on the noise statistics. The amplitudes of detected spikes were in the range of 10-80 |jV. Signal analysis and statistic calculations
were performed using custom-made software developed in Matlab [5]. Voltage pulse stimulation (±600 mV, 260 |js per phase, the first phase positive) was applied. Experimental protocol consisted of 3 h of spontaneous activity recordings before and after stimulation. The 3-hour recordings were divided into 60- and 20-minute intervals for analysis. Low-frequency voltage stimulation protocol included 30 pulses per each of 12-15 electrodes with 5-seconds inter-stimulus interval (Figure 1).
Connectivity analysis. Cross-correlation of activity observed in each electrode pair (i, j) was calculated as a percentage of spikes from i-th electrode recorded on j-th electrode with various time delays. We used the term 'connection delay' to denote the time delay corresponding to the maximum of cross-correlation function which was regarded as connection strength. If connection delay was longer than 3 ms (Figure 2), the pair of electrodes was treated as synaptically connected [17]. Connectivity map characterized by connection strength and connection delay was constructed for each 20- or 60-minutes raster. To measure connectivity changes over time, each recording was compared with the first recording and two parameters were calculated: connectivity coefficient (ADC%) calculated as the percentage of appeared and disappeared connections and stable connection strength change (SCSC) coefficient determined as the mean difference between connection strengths detected in both recordings. SCSC was measured in percentage accounting for the difference between the maxima of cross-correlation function characterizing the number of transmitted spikes.
Bursting activity comparison. Mean spiking rates per burst were compared for each electrode and each pair of recordings (rasters), respectively. The differential was determined as the difference between the mean number of spikes per burst for each pair of recordings (rasters), for each recording electrode separately (Figure 3 (a)). The difference in spiking rates per burst for each pair of recordings (rasters) for each electrode was quantified using the parameter denoted as overlap in burst spiking rate distributions. To calculate this parameter, we used the methods of clustering two sets of burst spiking rates. K-means clustering method was applied to all bursts for each pair of recordings (rasters A and B). As a result, the bursts were divided into two data sets (A' and B'). The sets of bursts with known raster origins (sets A and B) we compared to the set of bursts after cluster analysis (A' and B'). In this case, classification error was represented as the number of the patterns identified incorrectly as to the recording source (Na=A'HB and Nb=B'HA), where Na
Control (pre-stimulation) C1 C2 C3
Post-stimulation PS1 PS2 PS3
20 min 20 min 20 min 20 min 20 min 20 min
Spontaneous activity
recording, 1 h
LFS
Spontaneous activity
recording, 1 h
C1
C2
C3
PS1
PS2
PS3
1 h 1 h 1 h 1 h 1 h 1 h
Spontaneous activity recording, 3 h LFS Spontaneous activity recording, 3 h
AlSiA«^
Tib*?::
Figure 1. Experimental protocol of electrical stimulation and spontaneous activity recording; LFS — low-frequency stimulation
//////////////////////^^^^^
62 CTM j 2017 — vol. 9, No.2 A.A. Gladkov, V.N. Kolpakov, YI. Pigareva, V.B. Kazantsev, I.V. Mukhina, A.S. Pimashkin
o o a
I
s
s
o
0
3
1 i
S?
1
I
o &
s I
o %
a.
m .
" I
I
I
o
g
o
isj
«
w
CA oT
° £ a- o
CD o 03 c
& ÏÏ §8
o1
3 -r
03
CD
CD
CD
o.
CD 03
Kl =1
i < iÂ
5 % °2 03 es
03
s <
s
03
CA
o o
3
ÇA
ca'
I—I-
CD
CD O
■ o
3'i
CD CD __o
S" < a-
CD W o'cd
o
CD
5; 03. o c » CD
C/) ^
"i 3
^ CD
3 03
ca 3
s. "2.
CD
OCË
CD —1
03 T3
CA
—! . 03
cd CA 1—1-
03 3'
3 03
q. 1—1-
cd
03 q.
3!
cd -i 5'
ç/j o-
1—h c
cd 1—1-
CA 0
3
'tî ><
cq' c t. 1—1-
3"
cd 3'
w 03
^_^ CA
3
03.
^_^
0 cq
—■ o
c
H t3
3" 0
cd —h
1—1-
3 3"
c cd
3 cd
o- cd
cd o
-i
o —h o q.
cd
cd CA
cd o 3"
r-l- cd
0 q. r 0
cd
CA tt
CD CD 03
CL 3
cn o
I—I-
I—I-
o 5-
o o c
03 Q.
,2 ®
CQ û CD
CA =;■
O- CD C
3 §
= a.
S <8
s ®
= §
çd S
iv) 3 œ 0
0 t3 3
3 0 3 3
1—1-
00 00 03 3 0'
(—t- 3" cd 0
3' c/) o ^ cd
c 03 1—1-03 3. CA cd o r-l-—1
0 3 o- 0
3 cq c q.
cd
w —h c/3
1—1-
3" 03 3! o 3 q. 3' cq 03 03 ><
cd 03 03 CA
-i >< o
CA iv) 1—1-<' 03
3' 0 .nr 3
D
c/) C/)
o si
s. I
CD S-
□ ÛJ 3
i- w S
o
Q.
O
3 o"
03 ® 3 o
s
(—t-
2" «fl CD ci.
^ I
c/) E.
q3
O =!:
V.s
CD -h § £
Q. O
lit
T3 T3
o o
03
3
T3 03.
o-g
ig £ 3
3 «Q
C/)
m ^
->■ fi)
co ¡s-
'—- o
— ' s" i—t- w
S* (D ^ 3- -3
VI (D (D 9
C 3
ï* (n S02
w ® 3 w
"o №
» 3 -=?
S 3 03
® çd
o- c/5
cd
5' ®
cd
03 O
c
=!: o
< q.
03 cd
o'
3 ca cd
3
3 03
3 l'as 3
x cd
3 £ c/) 3
■g
03
a
cq £
g CD s cd cd ç/)
o!
cd
-?
3 S
cd cd
03 3^ -rj
tn 0) ä
S" 3 TT
çd =;■
o
^ tq
s
o 3-fit 3
03 (q
03
I—I-
cd
o 03
cd
O CA o
S g-
® s
03 cd
rt.
3"
3 %
It
c/) t3
cl 03
n3 cd
0 o
1 -
il
cd (fl
03 cd c5' cd
cd o o ><
«A CA
03 e.
3 03 (q
cd o
03
0 CA
1—1- 1—1-
<■ 3'
1—h >< c
S 03 1—1-
cd 0'
3
cd
q. H
<' 3"
q! cd
cd cd
q. -i
__ cd
3 o
1—1- o o -i
q.
o
3 3'
cd 1 cq
c/)
O 03 3" 3 03 q.
3 3 ?
s
sr M
o^ cfl
3 5'
03 o"
âR
v>
03
cd t3
03
s
cd ? s 03
o q. c,
œ
^ 3-
3 w 3
o 2, 3
3. ®
tq -1
v>
çd. 03
cd o
s
3 Si
03. œ 3
3= -g â ¡■S
I—I- ^
cd c/)
o
O" t/) cd 0:
» g- 3
«e- T < <-
cd
cd 03
cd <.
I—I-
3"
a- _.
! tq
03
1—1-
cd 3 §
m
03 S 3 cq
m s
■o
25: 3
ca 7
< A cd
CA —■
3.
o CD
03 3.
o
2 -Q
■ tq
Si a
cd c
o £ CA
03 o
<
cd CA
cq'
03
I—I-
cd q.
cd ' -q
CD 03 03 3 < O (O 2 3" CD CD
Si
i m
q. 1—1- -s —■
cd cd 03 (—t- ><
>< cd
—h o "o cd
H cd -s 03
3" o
cd cd —^ 3"
03 cd
o 3 c/)
3" q. c/)
t3 —h cd* 0 3 en o
03 -î
cd cd
03
0 —h o cd o 3" o —h
CA
CD
o o
CD 03
3 CD CD Q-q.
03
03 O (0 ^ 03
< o
CD "O < ^ ^
CD 3 O M (û 3 03
' Si
CD CD
«a -.
03
^ CA
CA CD ~ CA
CD
O CD
=î 03
^ 1-1-
CD CD
o ^ 1—1- 3" c/) 1—1-
cd cd 03 1—1-
—1 CÄ
03 q. 1—1-
t3 cd 0'
cq 03.
03 -i cd 1—1-
3 cd cd
cd CA 1—1-
-q o
c —h c/j
03_ c
0 q. 0 3"
cd -i 3
cd 3 cd 03
03 o CA
cd n.
< cd
03
c cd CA H 3" c -i cd t3 -i
c/3 cd
o J c/)
—h cd
3 1—1-
cd cd CA
q.
ca'
o —h
CD
£1 o -1
O CD
sO O
o^ o
O (D
g œ
3 a
T3
q) o q. (/) (D
(D q.
q) q. T3
0
ä
1
(J)
f—K
3'
c
a
o
?
o -h
o ° o
=J (D
rt- Q)
° §-
0 -1
1 s, ® s
in «P
o|
(D
t] 8
Sä I S
8 o
§ N3
8' o
o
o —h
T3
s
CL en
en' "a
1—1-
CT c CD en
(—1-
0' S 0
0 0
—h Q.
0 < CD Q.
CD
i
Q) CD
T3
S
5'
O" CD
c 3
c/5 c-1- 0 < CD
c/> Q.
"a
—h
O
CQ 3
—,
Q) c-1- Q)
CD a
Q)
CL
en' en
1—1- en'
CT
c (—1-
0 11
en
C/) 0
"a
CD c-1-
c
CD s
Q) en
O
S
0
0 <
Q. Q)
C
CQ CD
q. q.
tn iÎ
(D O
-Q (D d
(d 5
3
o"
(D
o "g
a- ^
!S ®
ä %
-, (D
° 8 o" =h
S Ä
8 8 =j
ô' =r
q) o
3 -n
0 ? cq' c
3 -t -1 (d
tt (i)
d o H
r-f 3"
0
r-t- >< 5' â
c
0
3
0
0
0
-h
O
¥
-h
o
c
(d
3
0
<
tn
r+
3'
c
d
r-f
o'
3
0
3
r-f
3"
0
in
■5
3'
cq
-t
d
r-f
0
0
-h
CQ
o m
O 03
3 3
=t o.
03 <
a I
rt- cd'
o q.
CD CO CA
? 5'
CD -, 3. 03
03 CA CD
— cn 3 >
a- os
c 3 i q.
cfl
s
cfl ^ T3 5
- 03
CA —
«Si
w m w w
3- CD
CD ^
U CD
rt- CD ^
O ^ CD
03 CD — -s
3. o
= 03 en C tQ ^ ^ CD rt-
O" —1 CD ? a o-o ^
S. « ^
-h o 1 -h 3
? CA CD E C 03
CA o 3 & 3 <*
f
zä«
03 CD CD + 3 3.
z w (g
o" 10
" I—I-
o
CD C/)
03 CD =î
w œ
2 0 3
CD p
CA E
3 CD 3 = ~ o" 3 > ÇD
a^rr
C CD C/)
w Êj 5T
(—t- 3
«a o
^ o 3
§ E" 3
3 a
3 s ®
CD —■ (—t-
t- 3
>ÎP CD
Number of electrodes (%)
(D
en o
Number of electrodes (%)
O J^ 00 N3 O)
(D Q)
8 n>
o. "a
(D
Q)
Q) <
So»
3 c ID </l W CD
O Ö:
(/) "O ^ "Ç O CD
^ 2. — =: CD =f CD O (/)
O CD CD -1
ÇD. Cn C® -n
O ^ H- CD
O ""
<" — 0' "Ç c± o =J TT — CD
~ en
SS -Ä —
x t3 -5
Q) -:
-s
ô" 8 -
2-2-1
3D® ?
"15 8«
-H ^ 3
I - S"
O _. ?
o i O
3
CD
(/)
Q)
„ (-1-O CD
3
< J?
en 'S en
"a. o 2
s: o
3 CD 3
en x a.
S-1 o < "D
o S
?"
§ 3
o o
Si
B-. * u
§ s
3 o
o t» 5T ?
Electrodes
Burst profile, spikes per bin
Spikes transferred (%)
œ
cû_ CD O
CD.
O" c
c/5
c-1-
T3
o
œ"
en
O o
î? (D
r+ _^
S' 0
ho o
ai o
«-O œ
en tt
MÄ»
"O C/) CD
en
>
a
£
z n P3
a
PO
P3 cc
P3 >
P3 n ffi
s
stimulation. Some electrodes were found to show the increase of spiking rate in the bursts during 2 h after stimulation.
Burst activation changes induced by low-frequency stimulation. The dominant pathway of signal propagation was determined based on delays of burst activation for different parts of the network. Burst activation middle was analyzed for each channel before and after stimulation. Changes in burst activation induced by low-frequency stimulation were more significant than spontaneous changes on a small group of electrodes (network sites) (Figure 4 (a), (c)).
Connectivity reorganization in the network.
Connectivity changes induced by stimulation were studied by analyzing the sequence of three-hour activity: two one-hour raster plots recorded before stimulation and one-hour plot recorded after stimulation. We analyzed connectivity coefficient (ADC%) (the number of appeared/disappeared connections) during each hour comparing each raster plot to the connectivity at the first hour of recording (C1: pre-stimulation (control), PS1: post-stimulation). Statistical significance of the found difference was estimated using Mann-Whitney U test, p<0.01, n=5 (Figure 5 (a)). Median ADC% values
30
<D
■o
O
O ^
CO s_
■g 10
C1-C2 C1-C3 C1-PS1 C1-PS2
-20 -10 0 10 Inter-channel "activation middle" change (ms) a
—C1-C2 —C1-C3 —C1-PS1 C1-PS2
-20 -10 Inter-channel "activation middle" change (ms) b
Figure 4. The influence of low-frequency stimulation on network burst activation:
(a) distribution of difference in delays of burst activation middle; (b) enlarged left part of the distribution in Figure 4 (a). C1-C2, C1-C3 — comparison of two control recordings; C1-PS1, C1-PS2 — comparison of pre-stimulation (control) and post-stimulation recordings
60 min time scale
Coefficient of connectivity changes (%)
Coefficient of common connection strength changes (%)
C2-C3 C2-PS1 C2-C3 C2-PS1
a b
20 min time scale
Coefficient of connectivity changes (%)
Coefficient of common connection strength changes (%)
p=0.0003
y
1
0 C2-C3 C2-PS1
c
p= 0.0043
_1_
C2-C3 C2-PS1 d
Figure 5. Network connectivity changes at 60-and 20-minute time intervals after low-frequency stimulation:
(a), (c) connectivity coefficient and (b), (d) stable connection strength coefficient before and after stimulation at 60-minute time interval (a), (b) and 20-minute time interval (c), (d). Red line depicts median values, blue color depicts the boundaries of 25 and 75 percentiles; C2-C3 — comparison of two control recordings; C2-PS1 — comparison of pre-stimulation (control) and post-stimulation recordings
////////////////////^^^^
64 CTM i 2017 — V0l. 9, N0.2 A.A. Gladkov, V.N. Kolpakov, Y.I. Pigareva, V.B. Kazantsev, I.V. Mukhina, A.S. Pimashkin
for pre- and post-stimulation data were 31.0 and 54.4, respectively (i.e. they changed by 1.75 times). Differences (Mann-Whitney U test, p<0.05, n=5) were also found when SCSC. The medians were 0.006 and 0.0125, their relative difference was 2.08. However, at the time scales of 20 min the differences in SCSC coefficient were more significant (Man-Whitney U test, p<0.001). The medians of ADC% coefficient were 0.5 in control and 33.1 when comparing pre- and post-stimulation activity (Figure 5 (c)). The relative change was 66.2 in contrast to 1.75 for 60-minute recordings. SCSC coefficient medians equaled 0.0001 in control and 0.0095 comparing pre- and post-stimulation activity. The relative change was 95.0 in contrast to 2.08 for 60-minute recordings.
The total number of connections at 60-minute time interval was found to remain unchanged after the stimulation (Figure 6). Significant changes in SCSC coefficient and change in connectivity coefficient ADC% clearly demonstrated reorganization of functional connectivity in the network after low-frequency stimulation.
Discussion. In this study, dissociated hippocampal neurons were cultured on microelectrode arrays. Spontaneous bursting activity was recorded for several hours (up to 8 h), then low-frequency electrical stimulation was applied. Two stimulation protocols were applied: low- and high-frequency stimulation. Low-frequency stimulation was performed through 12 randomly selected electrodes at 0.1-0.5 Hz, high-frequency stimulation — through 2 channels and 8 electrodes (4 electrodes per one stimulation channel). There were applied 150 series, each consisting of 20 stimuli, the time interval between stimuli was 100 ms, between series — 6 s, the delay between the stimuli of the first and the second channels equaled 10 ms. The recorded activity was divided into 60- and 20-minute intervals and followed by spike-transfer function analysis to estimate connectivity.
Functional connectivity was estimated by cross-correlation analysis of spiking patterns of spontaneous network activity recorded through different microelectrodes (by analogy with [17]). The paper by le Feber [12] demonstrates that the integral indicator of neuronal network connectivity reorganization (Euclidean space between connectivity strength matrices for all pairs of recording electrodes before and after stimulation) changes after electrical stimulation.
We have developed the indicators of neuronal network activity alteration: connectivity coefficient, stable connection strength change coefficient, mean spiking rate change, overlap in burst spiking rate distributions, the difference in delays of burst activation middle, which provide the possibility to evaluate functional connectivity of neuronal network. Changes in connectivity strength and connectivity reorganization (number of appeared/
300
250
Connectivity dynamics
-•-Experiment 1 -•-Experiment 2 -•-Experiment 3 -•-Experiment 4 -•-Experiment 5 -•- Experiment 6 -•-Experiment 7
— Mean
— Stimulation
Time (h)
Figure 6. Total number of correlations before and after low-frequency stimulation in dynamics
disappeared connections) indicate significant functional changes in culture network due to low-frequency stimulation. It should be noted that, at the same time, the total number of connections remained unchanged.
Using the introduced indicators, it has been established for the first time that neuronal network activity alteration observed on the time scale of 20 min is much more profound than that on the time scale of 1 h. This suggests that spontaneous network connectivity changes may inhibit the changes induced by stimulation.
The introduced indicators of neuronal network activity alteration provide efficient and simple tools for evaluation of functional state of neural networks in vitro when lengthy electrophysiological recordings on microelectrode arrays are made. These indicators can be applied in studies on testing pharmacological drugs intended to improve or rehabilitate the functional state of neuronal tissue.
Conclusion. Electrical stimulation by short two-phase square-wave pulses (±800 mV, 260 ps per phase, positive first) of low frequency (at 3-second inter-stimulus interval during 20 min) causes functional changes in the network of durably cultured dissociated brain cells, particularly, the changes in connectivity architecture of network activity with stable total quantity of cross-correlation connections. The developed indicators, connectivity coefficient and stable connection strength change coefficient appear to be reliable indicators of activity-dependent plasticity in neuronal network. Changes in the parameters of mean spiking rate, overlap in burst spiking rate distributions, difference in delays of burst activation middle after the stimulation are observed only within a few network sites (recording electrodes). Neuronal network activity alteration induced by low-frequency stimulation can be preserved for tens of minutes. At a long time interval (1 h), spontaneous connectivity changes inhibit stimulation-induced pattern.
Study Funding. This study was supported by the Russian Science Foundation (No.14-19-01381).
Connectivity of Neural Network in Dissociated Hippocampal Culture Grown on Microelectrode Array СТМ j 2017 — vol. 9, No.2 65
Conflicts of Interest. The authors have no conflicts of interest to disclose.
References
1. Chiappalone M., Massobrio P., Martinoia S. Network plasticity in cortical assemblies. Eur J Neurosci 2008; 28(1): 221-237, https://doi.org/10.1111/j.1460-9568.2008.06259.x.
2. Bakkum D.J., Chao Z.C., Potter S.M. Long-term activity-dependent plasticity of action potential propagation delay and amplitude in cortical networks. PLoS One 2008; 3(5): e2088, https://doi.org/10.1371/journal.pone.0002088.
3. Eytan D., Brenner N., Marom S. Selective adaptation in networks of cortical neurons. J Neurosci 2003; 23(28): 93499356.
4. Marom S., Shahaf G. Development, learning and memory in large random networks of cortical neurons: lessons beyond anatomy. Q Rev Biophys 2002; 35(01): 63-87, https:// doi.org/10.1017/s0033583501003742.
5. Pimashkin A., Kastalskiy I., Simonov A., Koryagina E., Mukhina I., Kazantsev V. Spiking signatures of spontaneous activity bursts in hippocampal cultures. Front Comput Neurosci 2011; 5: 46, https://doi.org/10.3389/fncom.2011.00046.
6. Shahaf G., Eytan D., Gal A., Kermany E., Lyakhov V., Zrenner C., Marom S. Order-based representation in random networks of cortical neurons. PLoS Comput Biol 2008; 4(11): e1000228, https://doi.org/10.1371/journal.pcbi.1000228.
7. Maccione A., Garofalo M., Nieus T., Tedesco M., Berdondini L., Martinoia S. Multiscale functional connectivity estimation on low-density neuronal cultures recorded by high-density CMOS Micro Electrode Arrays. J Neurosci Methods 2012; 207(2): 161-171, https://doi.org/10.1016/j. jneumeth.2012.04.002.
8. Makarov V.A., Panetsos F., de Feo O. A method for determining neural connectivity and inferring the underlying network dynamics using extracellular spike recordings. J Neurosci Methods 2005; 144(2): 265-279, https://doi. org/10.1016/j.jneumeth.2004.11.013.
9. Ullo S., Nieus T.R., Sona D., Maccione A., Berdondini L., Murino V. Functional connectivity estimation over large networks at cellular resolution based on electrophysiological recordings and structural prior. Front
Neuroanat 2014; 8: 137, https://doi.org/10.3389/fnana.2014. 00137.
10. Wagenaar D.A. Controlling bursting in cortical cultures with closed-loop multi-electrode stimulation. J Neurosci 2005; 25(3): 680-688, https://doi.org/110.1523/ jneurosci.4209-04.2005.
11. Poli D., Pastore V.P., Massobrio P. Functional connectivity in in vitro neuronal assemblies. Front Neural Circuits 2015; 9: 57, https://doi.org/10.3389/fncir.2015.00057.
12. le Feber J., Witteveen T., van Veenendaal T.M., Dijkstra J. Repeated stimulation of cultured networks of rat cortical neurons induces parallel memory traces. Learn Mem 2015; 22(12): 594-604, https://doi.org/10.1101/lm.039362.115.
13. Bologna L.L., Nieus T., Tedesco M., Chiappalone M., Benfenati F., Martinoia S. Low-frequency stimulation enhances burst activity in cortical cultures during development. Neuroscience 2010; 165(3): 692-704, https://doi.org/10.1016/j. neuroscience.2009.11.018.
14. Maeda E., Kuroda Y., Robinson H.P.C., Kawana A. Modification of parallel activity elicited by propagating bursts in developing networks of rat cortical neurones. Eur J Neurosci 1998; 10(2): 488-496, https://doi.org/10.1046/j.1460-9568.1998.00062.x.
15. Vajda I., van Pelt J., Wolters P., Chiappalone M., Martinoia S., van Someren E., van Ooyen A. Low-frequency stimulation induces stable transitions in stereotypical activity in cortical networks. Biophys J 2008; 94(12): 5028-5039, https:// doi.org/10.1529/biophysj.107.112730.
16. Wagenaar D.A., Pine J., Potter S.M. An extremely rich repertoire of bursting patterns during the development of cortical cultures. BMC Neurosci 2006; 7(1): 11, https://doi. org/10.1186/1471-2202-7-11.
17. le Feber J., Rutten W.L., Stegenga J., Wolters P.S., Ramakers G.J., van Pelt J. Conditional firing probabilities in cultured neuronal networks: a stable underlying structure in widely varying spontaneous activity patterns. J Neural Eng 2007; 4(2): 54-67, https://doi.org/10.1088/1741-2560/4/2/006.
18. Pimashkin A., Gladkov A., Agrba E., Mukhina I., Kazantsev V. Selectivity of stimulus induced responses in cultured hippocampal networks on microelectrode arrays. Cogn Neurodyn 2016; 10(4): 287-299, https://doi.org/10.1007/ s11571-016-9380-6.
//////////////////////^^^^
66 СТМ i 2017 — v0l. 9, N0.2 A.A. Gladkov, V.N. Kolpakov, Y.I. Pigareva, V.B. Kazantsev, I.V. Mukhina, A.S. Pimashkin
