- 1Department of Neurotechnology, Institute of Biology and Biomedicine, National Research Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod, Russia
- 2Molecular and Cell Technologies Group, Central Scientific Research Laboratory, Privolzhsky Research Medical University, Nizhny Novgorod, Russia
Alzheimer’s disease (AD) is a widespread chronic neurodegenerative pathology characterized by synaptic dysfunction, partial neuronal death, cognitive decline and memory impairments. The major hallmarks of AD are extracellular senile amyloid plaques formed by various types of amyloid proteins (Aβ) and the formation and accumulation of intracellular neurofibrillary tangles. However, there is a lack of relevant experimental models for studying changes in neural network activity, the features of intercellular signaling or the effects of drugs on the functional activity of nervous cells during AD development. In this work, we examined two experimental models of amyloidopathy using primary hippocampal cultures. The first model involves the embryonic brains of 5xFAD mice; the second uses chronic application of amyloid beta 1-42 (Aβ1-42). The model based on primary hippocampal cells obtained from 5xFAD mice demonstrated changes in spontaneous network calcium activity characterized by a decrease in the number of cells exhibiting Ca2+ activity, a decrease in the number of Ca2+ oscillations and an increase in the duration of Ca2+ events from day 21 of culture development in vitro. Chronic application of Aβ1-42 resulted in the rapid establishment of significant neurodegenerative changes in primary hippocampal cultures, leading to marked impairments in neural network calcium activity and increased cell death. Using this model and multielectrode arrays, we studied the influence of amyloidopathy on spontaneous bioelectrical neural network activity in primary hippocampal cultures. It was shown that chronic Aβ application decreased the number of network bursts and spikes in a burst. The spatial structure of neural networks was also disturbed that characterized by reduction in both the number of key network elements (hubs) and connections between network elements. Moreover, application of brain-derived neurotrophic factor (BDNF) recombinant protein and BDNF hyperexpression by an adeno-associated virus vector partially prevented these amyloidopathy-induced neurodegenerative phenomena. BDNF maintained cell viability and spontaneous bioelectrical and calcium network activity in primary hippocampal cultures.
Introduction
Alzheimer’s disease (AD) studies are becoming more relevant each year due to the increase in the life expectancy of the population and the accumulation of information regarding AD polyetiology (Dubois et al., 2014; Hampel et al., 2015). The features of AD pathological processes and the development of new strategies to prevent neurodegeneration are actively pursued worldwide (Hadar and Gurwitz, 2018; Cao et al., 2018; Wang N. et al., 2018). Nevertheless, there is no clear therapeutic solution for highly accelerated neurodegeneration, even if it is diagnosed at an early stage.
Investigations of AD processes have raised questions about the possibility of using endogenous regulatory molecules, such as neurotrophic factors, to correct neurodegeneration at different stages of pathology development. The content of neurotrophic factors decreases with neurodegeneration, and this process correlates with AD stages (Peng et al., 2005; Wang and Holsinger, 2018). Brain-derived neurotrophic factor (BDNF) is a potent biological agent that maintains cell viability and functional neuron activity in various pathological states, including severe genetically determined neurodegenerative diseases (Criscuolo et al., 2015; de Pins et al., 2019; Choi et al., 2018). Viral constructs carrying the BDNF gene are a promising therapeutic strategy to restore BDNF levels in the brain. Several studies have indicated the efficacy of viral vectors carrying neurotrophic factor genes in the treatment of Parkinson’s disease (Lim et al., 2010; Cheng et al., 2018; Tereshchenko et al., 2014). The establishment of approaches to use viral constructs carrying the BDNF gene in AD has been carrying out for the past 10 years (Nagahara et al., 2009; Nagahara et al., 2013; Jiao et al., 2016).
AD is characterized by significant variability in the age of disease manifestation and the rate of disease progression. The major hallmarks of AD are the pathological accumulation of amyloid beta (Aβ) protein in the form of extracellular plaques in brain parenchyma and capillaries and the abnormal phosphorylation of tau protein, which forms neurofibrillary tangles (Bourdenx et al., 2017). Aggregation of Aβ and phosphorylated tau occurs gradually; monomers are aggregated into oligomers in neurons and then collected into fibrils, leading to the formation of amyloid plaques and neurofibrillary tangles (Walker et al., 2013; Cline et al., 2018).
The elaboration of experimental models of AD is a key for better understanding AD pathogenesis and assessing the potential of new therapeutic approaches for effective neurodegenerative process correction (Drummond and Wisniewski, 2017). In vivo models are currently the most frequently used experimental models of AD and are mostly based on transgenic mice that overexpress human genes associated with the familial form of AD (such as the familial AD (FAD) lines), resulting in the formation of amyloid plaques (Mucke et al., 2000; Webster et al., 2014).
The adequacy of any biological or mathematical model depends on specific tasks and possible approaches to its solution. Investigation of neural networks as the minimal functional unit of the nervous system responsible for the processes of reconsolidation and storage of information is considered one of the principal aspects of studies on the neurodegeneration processes. A neural network is not only a functionally connected complex of neurons but also a single functional ensemble capable of responding in a consolidated manner to changes caused by both external and internal stimuli (Yuste, 2015; Mishchenko et al., 2019). A single neural network in the native brain is extremely difficult to study and cannot be examined in a comprehensive manner at this time. Primary hippocampal cultures are considered an adequate biological model that allows the study of individual cellular and network reactions under stress and the effects of neuroprotectants in a chronic experiment with the possibility of multiple measurements of neural network activity (Johnstone et al., 2010; Vedunova et al., 2015; Hasan and Berdichevsky, 2016). To investigate the changes in functional neural network activity caused by AD development, we used a protocol for creating primary neuronal cultures obtained from 5xFAD murine embryos. Notably, however, models using transgenic animals mostly simulate familial forms of AD, which account for only 5% of all cases of this pathology (Bilkei-Gorzo, 2014; Karch et al., 2014; Kim et al., 2014). These models often lack a complex of pathological traits exhibited by patients with AD. These transgenic mice are characterized by amyloid plaques, compromised synaptic transmission and memory impairment, but these symptoms are not always accompanied by neuronal loss and, most importantly, neurofibrillary tangle formation (Oakley et al., 2006). The poor correlation between preclinical research of new therapeutic drugs and clinical trials is probably associated with this issue (Banik et al., 2015; Cummings et al., 2018). Therefore, the development of relevant experimental models will provide a more complete view of pathogenic processes in AD. Such opportunities may be possible in an in vitro amyloidosis model based on synthetic amyloid peptide application (Stancu et al., 2014; Villalobos Acosta et al., 2018; Mango et al., 2019).
Our present study is devoted to adapting an in vitro amyloidopathy model that allows investigation of the functional activity of neural networks. Using our model, we also studied the influence of BDNF on cell viability and the reorganization of neural networks in AD development.
Materials and Methods
Ethics Statement
All experimental protocols used in this study were approved by the Bioethics Committee of Lobachevsky University and carried out in accordance with Act708n (23.08.2010) of the Russian Federation National Ministry of Public Health, which states the rules of laboratory practice for the care and use of laboratory animals, and the Council Directive 2010/63 EU of the European Parliament (September 22, 2010) on the protection of animals used for scientific purposes. C57BL/6J mice were killed by cervical vertebra dislocation, and their embryos were then surgically removed and sacrificed by decapitation.
Primary Neuronal Cultures
Primary hippocampal cells were obtained from murine embryos (day 18 of gestation). A detailed protocol for culture preparation is described in Vedunova et al., 2015. Hippocampi were surgically isolated. Cell dissociation was achieved through mechanical dissection followed by incubation for 20 min in 0.25% trypsin-EDTA solution (Gibco, 25200056, United States). The obtained cell suspension was centrifuged at 1000 rpm for 3 min. Then, the cell pellet was resuspended in NeurobasalTM medium (Gibco, 21103049, United States) supplemented with 2% B27 (Gibco, 175040446, United States), 0.5 mM L-glutamine (Gibco, 25030024, United States) and 5% fetal bovine serum (FBS) (PanEco, K055, Russia). To perform a viability assessment, immunocytochemical analysis and registration of functional calcium activity, we placed cells on coverslips (18x18 mm) pretreated with polyethyleneimine solution (1 mg/mL) (Sigma-Aldrich, Germany). For electrophysiological experiments, cells were cultured on multielectrode arrays (MEAs; MEA60, Multichannel, Germany). The initial density of cells was 9000 cells/mm2. Half of the medium containing 0.4% FBS was replaced every third day. Cell viability was maintained under constant conditions of 35.5°C, 5% CO2 and a humidified atmosphere in a CO2 incubator (Sheldon Manufacturing, United States).
5xFAD Embryo Genotyping
To obtain primary hippocampal cultures from 5xFAD mice, we performed genotyping using polymerase chain reaction (PCR) for wild-type and mutant embryos. During genotyping procedures, cell viability was maintained by placing embryonic hippocampal tissue in warm Neurobasal medium on a thermoshaker (750 rpm, 37°C).
The 5xFAD mouse line expresses both mutant human APP695, which harbors the Swedish mutation (K670N, M671L), the Florida mutation (I716V), and the London mutation (V717I), and human PSEN1, which harbors two FAD mutations (M146L and L286V). Both transgenes are expressed under the control of the mouse Thy1 promoter to induce overexpression in the brain. Primers for PCR are shown in Table 1. PCR was performed with Taq polymerase in a C-1000 Thermal Cycler (Bio-Rad) with the following steps:
Lid: 110° C
Volume: 20 μl
(1). Denaturation: 94°C, 3:00;
(2). Denaturation: 94°C, 0: 30;
(3). Annealing: 64°C, 1:00;
(4). Elongation: 72°C, 1:00;
(5). Beginning at step 2, repeat 35 times;
(6). 72°C, 2:00;
(7). Storage: 4°C, ∞.
PCR was carried out in 0.2-ml disposable tubes with optically transparent caps. The reaction mixture for each gene was prepared in accordance with the following protocol (Table 2):
In a test tube, 19.5 μl of the prepared PSEN1 reaction mixture, 0.5 μl of a DNA template sample, 19 μl of the prepared amyloid precursor protein (APP) reaction mixture and 1 μl of a DNA template sample were added.
To discriminate the genotype of individual mice, we used 2% agarose gel electrophoresis.
The expected results for the transgenes were a 610-bp fragment for PSEN1 and 377-bp fragment for APP, with a 325-bp fragment for PSEN1 and 324-bp fragment for APP as the internal positive controls.
Preparation of Amyloid β Treated With Hexafluoroisopropanol (HFIP)
To prepare 1 mM Aβ, we added HFIP solution (220 μl) directly to the lyophilized β-amyloid peptide powder (Aβ42) (InnovaGen, Sweden) and incubated the mixture at room temperature for 30 min. The obtained solution was transferred to microcentrifuge tubes and left in a hood overnight until a transparent film formed and HFIP evaporated. Next, dimethylsulfoxide (DMSO) was added to the tubes, which were mixed on a vortex for approximately 30 s and then centrifuged at 1000 rpm for 1 min. The Aβ-DMSO solution (5 mM) was placed in a sonicator for 10 min. To prepare a fibrillar β-amyloid peptide, we added 10 mM HCl (98 μL) to the Aβ-DMSO solution, mixed the solution for 15 s and then incubated it at 37°C for 24 h.
β-Amyloidopathy Model Based on Synthetic Aβ42
We used two protocols for in vitro β-amyloidopathy modeling. The first model involved one application of Aβ42 to the culture medium at a final concentration of 3.5 μM on day 10 of primary hippocampal culture development in vitro (DIV 10).
The second protocol involved chronic application of the same Aβ concentration. The obtained fibrillar amyloid-β was added to the culture medium every 48 h (i.e., after each change of culture medium) at a final concentration of 3.5 μM from DIV 10 to DIV 28 (Figure 1).
Figure 1. Scheme of the chronic amyloid-β application. The fibrillar amyloid-β was added to the culture medium every 48 h (i.e., after each change of culture medium) at a final concentration of 3.5 μM from DIV 10 to DIV 28.
AAV-Syn-BDNF-EGFP Virus Vector
To obtain a viral construct encoding the BDNF gene, we used the following plasmids: AAV-Syn-EGFP and helper plasmids pDP5, DJvector and pHelper. The same plasmids with cDNA EGFP were used to produce a control viral construct - AAV-Syn-EGFP.
The bacterial pUC19 plasmid served as the basis for the AAV-Syn-EGFP plasmid. This plasmid carries the sequences of the human synapsin (hSyn) promoter, woodchuck hepatitis posttranscriptional regulatory element (WPRE) enhancer, and SV40 polyA signal sequence flanked by inverted terminal repeats (ITRs) from adeno-associated serotype 2 virus (AAV2). The developed AAV-Syn-BDNF-EGFP included the following sequences: (1) the hSyn promoter, allowing expression of the gene of interest only in neuronal cells; (2) the regulatory WPRE enhancer, which markedly strengthens hSyn function; (3) a multilinker for open reading frame (ORF) cloning of the embedded gene; (4) the EGFP gene; (5) the SV40 polyA signal sequence flanked by ITRs from AAV2; (6) a gene cassette encoding ampicillin resistance (AmpR promoter and AmpR gene) for positive selection of colonies carrying this plasmid; and (7) a sequence corresponding to the nucleotide sequence encoding the functional BDNF protein.
The designed primers mBDNF-EcoRI-fw (5′-ATTGAATT CATGGGCCACATGCTGTCC-3′) and mBDNF-BamHI-rv (5′-AATGGATCCAATCTTCCCCTTTTAATGGTCAGTG-3′) were used. Detailed procedures for the generation and isolation of the viral construct are described in Mitroshina et al. (2018).
Primary hippocampal cultures were infected with AAV-Syn-BDNF-EGFP or vehicle (AAV-Syn-EGFP) on DIV 7. To infect the cultures, we mixed 3 μL of the viral sample with 50 μL of fresh culture medium. The medium was temporarily removed from the culture dishes, and the working solution of the viral vector was directly added to the cells. The cultures were incubated at 35.5°C with 5% CO2 for 20 min; then, the culture medium was returned to the cultures.
Cell Viability Analysis
To identify dead cell nuclei and the total number of cell nuclei, we stained primary hippocampal cultures with propidium iodide (Sigma, Germany) and bisbenzimide (Thermo Fisher, United States) according to Vedunova et al., 2015. Propidium iodide and bisbenzimide at concentrations of 5 μg/mL and 1 μg/mL, respectively, were added to the culture medium 30 min before viability measurements. Visualization of stained cells was carried out on a Leica DMIL HC inverted fluorescence microscope (Leica, Germany). We estimated the ratio of the number of propidium iodide-positive cells to the number of bisbenzimide-positive cells.
Immunocytochemical Analysis
The presence of Aβ in dissociated hippocampal cultures was detected by using primary chicken antibodies to Aβ (1:1000, Abcam, ab2539, United Kingdom) and secondary antibodies conjugated to a goat anti-chicken fluorescent marker (1: 100, Alexa Fluor 555, Invitrogen, 1719602, United States). Primary guinea pig antibodies to βIII-tubulin (1:1000, Synaptic systems, 302304, Germany) and secondary antibodies conjugated to a goat anti-guinea pig fluorescent label (1:100, Alexa Fluor 647, Invitrogen, 1711474, United States) were used as neuronal markers. The cultures were fixed in 4% paraformaldehyde in PBS for 20 min at room temperature; 0.2% Triton X-100/PBS was used for cell permeabilization. Immunocytochemically stained cultures were imaged using a Zeiss 510 NLO fluorescent confocal microscope (Carl Zeiss, Germany). The obtained images were analyzed using a custom ImageJ plugin. We conducted a comparative assessment of observation fields with equal densities of cells and imaged by the same laser power and photodetector settings. The average fluorescence intensity in the yellow channel, corresponding to the presence of Aβ42 in the observation field, was estimated.
Ca2+ Imaging
For imaging studies of functional Ca2+ activity in primary hippocampal cultures, we used a Zeiss 510 NLO fluorescent confocal microscope (Carl Zeiss, Germany) with a W Plan-Apochromat 20 × /1.0 objective. This method allows visualization of the functional neural network architecture at the cellular level. Oregon Green 488 BAPTA-1 AM (OGB-1) (0.4 μM, Thermo Fisher, United States), which was used as a calcium sensor, was dissolved in DMSO (Sigma, Germany) with 4% Pluronic F-127 (Thermo Fisher, United States) and then added to the culture medium for 40 min at 37°C and 5% CO2. OGB-1 was excited at 488 nm and recorded in the range of 500–530 nm. Time series of 512 × 512 pixel images of 420 × 420-μm fields of view were recorded at 2 Hz. A confocal pinhole of 1 airy unit was used to obtain an axial optical slice resolution of 1.6 μm. Detection and further analysis of Ca2+ oscillations was performed in the Astroscanner program. A more detailed description of the image analysis is provided in our previous articles (Vedunova et al., 2013; Zakharov et al., 2013). The following parameters of spontaneous Ca2+ activity were taken into account: the percentage of functional active cells and the duration (s) and frequency (the amount of Ca2+ events/min) of Ca2+ oscillations.
Electrophysiological Methods and Cross-Correlation Analysis
Spontaneous bioelectrical activity of neural networks in primary hippocampal cultures under chronic Aβ application was measured on days 14, 21, and 28 of cultivation. Extracellular action potentials were detected by MEAs (MEA60) and the USB-MEA-120 system (Multichannel system, Germany). The MEAs consisted of 59 planar TIN electrodes 8 × 8 grid) with a diameter of 30 μm and spaced 200 μm apart. Electrophysiological data were recorded simultaneously from 59 channels at a sampling rate of 20 kHz/channel. All signaling and statistical analyses were performed using custom-made software (MATLAB®6.0, United States).
Small network bursts were detected by calculating the total spiking rate (TSR), which considered the total number of spikes from all electrodes within 50-ms time bins. The criterion of a small network burst was the rapid appearance of a large number of spikes over four electrodes within a small (50-ms) time bin (Pimashkin et al., 2011; Vedunova et al., 2013). A more detailed description of the method for spikes and small burst detection is provided in our previous article (Mishchenko et al., 2019).
The following parameters of spontaneous bioelectrical activity of neural networks were analyzed: the number of small network bursts and the number of spikes per burst.
For cross-correlation analysis, the dataset obtained from electrophysiological recordings is presented as a raster plot. The network graph method was then used to detect the neuronal groups.
To assess the degree of synchronization between all pairs of cells, considering axonal delays, we calculated the proportion of transmitted spikes. The number of delayed synchronous spikes was normalized by the number of spikes received by the postsynaptic neuron nj. The cross-correlation matrix was calculated using the following formula:
Next, we selected the largest 5% of Cij coefficients and defined a set of indices, i.e., hubs of cells with a maximum number of functionally active connections. In addition, for each hub “i,” we calculated the number of connections to index i within the array Cij.
Next, the graph was constructed. The vertex size was proportional to the number of significant connections, and the edge of the graph corresponded to the functional connections of spikes transferred from one neuron to another at individual time points for each pair of axonal delays, i.e., τ ± δ/2 (Shishkina et al., 2018).
Real-Time PCR
Quantitative real-time PCR was used to analyze the levels of TrkB-FL receptor (TrkB gene) expression. Total RNA was isolated from primary hippocampal cell cultures on DIV 21 under chronic Aβ application using an ExtractRNA kit (eUROGEN, Russia). Then, cDNA was synthesized by Moloney murine leukemia virus (MMLV) reverse transcriptase (eUROGEN, Russia) and a random primer.
Quantitative real-time PCR was performed with qPCRmix-HS SYBR (eUROGEN, Russia) and an Applied Biosystems 7500 RT-PCR thermal cycler. The following primers were used:
TrkB-fw1, 5′-TTTCCGCCACCTTGACTTGTCT-3′;
TrkB-rv1, 5′-GTCGGGGCTGGATTTAGTCTCC-3′;
Oaz1_fw, 5′- AAGGACAGTTTTGCAGCTCTCC -3′; and
Oaz1_rv, 5′- TCTGTCCTCACGGTTCTTGGG-3′.
Data processing was carried out using the ΔΔCt method and a reference sample in which the target gene level was taken as a unit. Normalization was performed relative to the reference gene (Oaz1).
Statistical Analysis
All quantified data are presented as the mean ± standard error of the mean (SEM). Statistical analyses were performed using two-way ANOVA implemented in Sigma Plot 11.0 software (Systat Software, Inc.). The Student–Newman–Keuls (SNK) test was used as a post hoc test following ANOVA. Differences between groups were considered significant if the corresponding p-value was less than 0.05.
Results
Features of the Morphology and Spontaneous Calcium Activity of Neuron-Glial Networks in Primary Hippocampal Cultures Obtained From 5xFAD Murine Embryos
First, we adapted an in vitro amyloidopathy model, allowing studies of changes in neural network activity. To obtain a valid model, we conducted a single and chronic application of Aβ1-42 to primary hippocampal cultures obtained from C57BL/6 murine embryos and investigated the features of long-term cultivation of primary hippocampal cultures obtained from 5xFAD murine embryos.
Comparative morphological assessment did not reveal significant changes between primary hippocampal cultures obtained from wild-type and 5xFAD murine embryos over 28 DIV (Figure 2, Supplementary Figure S1). There was also no decrease in cell viability in either experimental group (Figure 2).
Figure 2. (A) Representative light field images of primary hippocampal cultures obtained from wild-type and 5xFAD mice. Scale bar – 100 μm. Comparative morphological assessment did not reveal significant changes between primary hippocampal cultures obtained from wild-type and 5xFAD murine embryos over 28 DIV. (B) Analysis of cell viability in primary hippocampal cultures obtained from wild-type and 5xFAD murine embryos. There was no decrease in cell viability in either experimental group.
Immunocytochemical analysis of Aβ accumulation in primary hippocampal cells obtained from C57Bl/6 and 5xFAD mice revealed that both types of cultures underwent endogenous Aβ synthesis, which was detected intraneuronally. The level of Aβ expression did not change significantly throughout the entire observation period (DIV 14, 17, 21, and 28) (Figure 3). Thus, an increase in the production of endogenous amyloid in the 5xFAD murine brain may occur in later stages. Numerous studies have shown that amyloid plaques are formed in the 5xFAD murine brain after eight months of age, although the first cognitive impairments are observed beginning at four months of age (Radde et al., 2006; Bilkei-Gorzo, 2014; Liu et al., 2017).
Figure 3. Representative confocal images of immunocytochemical staining of primary hippocampal cultures on DIV 28. (A) Primary hippocampal cultures obtained from C57BL/6 murine embryos. (B) Primary hippocampal cultures obtained from 5xFAD murine embryos. Red: fluorescence channel of βIII tubulin (neuronal microtubule marker); yellow: fluorescence channel of Aβ. Scale bar – 20 μm. The level of Aβ expression did not change significantly.
No morphological changes in primary 5xFAD murine hippocampal cultures were shown, and Aβ accumulations were not observed. However, the functional Ca2+ activity of primary 5xFAD murine hippocampal cultures was significantly altered in comparison with that of neuronal cultures obtained from wild-type murine embryos (Figure 4). Spontaneous network Ca2+ activity was detected beginning at DIV 10, consistent with previous data on primary hippocampal culture development in vitro (Shirokova et al., 2013). However, the percentage of cells that exhibited Ca2+ activity in the cultures obtained from 5xFAD murine embryos was significantly lower than that in the control cultures (DIV 21: control, 78.2 ± 10.02%; 5xFAD, 44.3 ± 7.02%; DIV 28: control, 65 ± 4.6%; 5xFAD, 35.5 ± 7.91%). The frequency of Ca2+ oscillations in the 5xFAD group of cultures was also significantly lower than the control values (DIV 21: control, 1.94 ± 0.14 oscillations/min (osc/min); 5xFAD, 1.16 ± 0.22 osc/min; DIV 28: control, 2.64 ± 0.46 osc/min; 5xFAD, 1.44 ± 0.36 osc/min). Additionally, the duration of Ca2+ oscillations in 5xFAD primary cultures was 1.55 (DIV 21) and 1.75 (DIV 28) times higher than those in the control cultures. The identified alterations in spontaneous calcium activity of 5xFAD primary cultures in the absence of pronounced morphological changes suggest that functional changes in nervous cells occur much earlier than visible neurodegenerative changes.
Figure 4. Main parameters of spontaneous calcium activity in primary hippocampal cultures obtained from wild-type and 5xFAD murine embryos (E18) during development in vitro. (A) Proportion of cells exhibiting Ca2+ activity; (B) number of Ca2+ oscillations per min; (C) duration of Ca2+ oscillations. ∗vs. “Control,” p < 0.05, ANOVA. The functional Ca2+ activity of 5xFAD primary neuronal cultures was significantly altered in comparison with that of wild-type murine primary cultures. The percentage of cells that exhibited Ca2+ activity and the frequency of Ca2+ oscillations in the cultures obtained from 5xFAD murine embryos were significantly lower than that in control cultures. Additionally, the duration of Ca2+ oscillations in 5xFAD primary cultures was 1.55 (DIV 21) and 1.75 (DIV 28) times higher than those in control cultures.
Features of the Morphology and Spontaneous Calcium Activity of Neuron-Glial Networks in Primary Hippocampal Cultures in the β-Amyloidopathy Model in vitro
Next, we used an in vitro β-amyloidopathy model established by adding synthetic Aβ1-42 to the culture medium. Two amyloidopathy modeling protocols were adapted, with the first involving one application of Aβ at a final concentration of 3.5 μM on DIV 10, and the second (chronic application) involving the introduction of Aβ to the culture medium at a final concentration of 3.5 μM every 48 h (i.e., after each change of culture medium) from DIV 10 to DIV 28. Notably, DIV 10 is a time of culture development characterized by a large number of chemical synapses and the formation of spontaneous bioelectrical and calcium network activity (Shirokova et al., 2013).
To verify the effectiveness of the experimental protocols used, we performed a cell viability assessment (Figure 5A) and immunohistochemical analysis of Aβ aggregate formation in primary hippocampal cultures (Figure 6). One application of Aβ did not affect primary culture viability. By contrast, chronic application of Aβ led to a significant decrease in the number of viable cells. On DIV 21, the percentage of living cells in the chronic Aβ group was 83.57 ± 6.95% (sham: 96.74 ± 0.74%) and continued to decrease at DIV 28 (sham: 97.47 ± 0.62%, chronic Aβ: 65.87 ± 8.39%).
Figure 5. (A) Analysis of cell viability in primary hippocampal cultures under single and chronic exogenous application of Aβ. * vs. “Sham,” p < 0.05, ANOVA. One application of Aβ did not affect primary culture viability, whereas its chronic application significant decreased the number of viable cells. (B–D) Main parameters of spontaneous calcium activity in primary hippocampal cultures under single and chronic Aβ application (B) Proportion of cells exhibiting Ca2+ activity; (C) number of Ca2+ oscillations per min; (D) duration of Ca2+ oscillations. * vs. “Sham,” p < 0.05, ANOVA. One Aβ application did not cause a pronounced effect on Ca2+ activity in neuron-glial networks. Severe Ca2+ activity suppression under chronic Aβ application was observed from DIV 21.
Figure 6. Representative confocal images of immunocytochemical staining of primary hippocampal cultures on DIV 28. (A) Sham culture, (B) primary culture with single Aβ application, (C) primary culture with chronic Aβ application. Red: fluorescence channel of βIII tubulin (neuronal microtubule marker); yellow: fluorescence channel of Aβ. Scale bar – 20 μm. A significant amount of Aβ associated with cells as protein globules was detected only in cultures with chronic Aβ application.
On DIV 28, immunocytochemical analysis revealed a significant amount of Aβ associated with cells as protein globules only in cultures with chronic Aβ application (Figure 6).
Moreover, the effect of Aβ application using the two protocols on functional calcium activity of neuron-glial networks in primary hippocampal cultures was studied. In addition to the lack of influence on cell viability, one Aβ application did not cause a pronounced effect on Ca2+ activity in neuron-glial networks. No significant changes in the number of cells exhibiting Ca2+ activity or the frequency of Ca2+ oscillations throughout the observation period were revealed (Figures 5B–D). The detected decrease in the duration of Ca2+ oscillations on DIV 14 (sham: 8.03 ± 0.73 s, single Aβ: 5.84 ± 0.85) was normalized to sham values by DIV 21 (sham: 8.26 ± 0.81 s, single Aβ: 7.82 ± 0.61).
Severe Ca2+ activity suppression under chronic Aβ application was observed from DIV 21, consistent with the cell viability analysis. Compared to the sham group, the chronic Aβ group had a significantly lower number of cells that exhibited Ca2+ activity on DIV 21 and DIV 28 (DIV 21: sham, 80.06 ± 3.40%; chronic Aβ, 29.01 ± 16.38%; DIV 28: sham, 68.70 ± 9.53%; chronic Aβ, 17.18 ± 9.48%). On DIV 28, the chronic Aβ group exhibited a significantly lower Ca2+ oscillation frequency (sham: 2.49 ± 0.34 osc/min, chronic Aβ: 0.35 ± 0.13 osc/min) and a significantly higher Ca2+ oscillation duration (sham: 6.2 ± 0.23 s, chronic Aβ: 13.24 ± 2.37 s) than the sham group.
In further studies, we used the AD model based on chronic Aβ application because it caused the most pronounced neurodegenerative changes.
Effects of Chronic BDNF Application on Cell Viability, Bioelectrical and Calcium Activity of Primary Hippocampal Cultures in the Chronic β-Amyloidopathy Model
To investigate the possible neuroprotective effect of chronic BDNF application in our amyloidopathy model, we analyzed primary hippocampal cell viability. On DIV 21, when Aβ application caused a significant decrease in culture viability, both recombinant BDNF and AAV-Syn-BDNF preserved the number of viable cells (sham: 87.63 ± 3.72%, Aβ: 72.52 ± 4.18%, Aβ+BDNF: 77.11 ± 6.22%, Aβ+AAV-Syn-BDNF: 82.34 ± 5.60%). The use of control virus vector AAV-Syn-EGFP did not preserve the viability of neuronal cultures. The viability parameter did not differ from that in the Aβ group and was significantly lower than the values in the sham group (Aβ+AAV-Syn-EGFP: 69.41 ± 5.87%). On DIV 28, the percentage of viable cells in the Aβ and Aβ+AAV-Syn-EGFP groups was decreased to 61.34 ± 3.38% and 58.72 ± 5.31% respectively, whereas in the Aβ+BDNF and Aβ+AAV-Syn-BDNF groups, this parameter was significantly higher at 76.96 ± 2.20% and 72.01 ± 2.95%, respectively (Figure 7A).
Figure 7. (A) Analysis of cell viability in primary hippocampal cell cultures under exogenous administration of Aβ and BDNF. * vs. “Sham,” # vs. “Aβ,” p < 0.05, ANOVA. On DIV 21 and 28 when Aβ and AAV-Syn-EGFP application caused a significant decrease in culture viability, both recombinant BDNF and AAV-Syn-BDNF preserved the number of viable cells. (B–D) Main parameters of spontaneous calcium activity in primary hippocampal cell cultures under chronic exogenous application of Aβ and BDNF during development in vitro. (B) Proportion of cells exhibiting Ca2+ activity; (C) number of Ca2+ oscillations per min; (D) duration of Ca2+ oscillations. ∗vs. “Sham,” # vs. “Aβ,” p < 0.05, ANOVA. On DIV 28, the activity of cells in BDNF-treated cultures was significantly higher compared to those in cultures exposed to chronic Aβ. The frequency of Ca2+ oscillations in the Aβ+BDNF and Aβ+AAV-Syn-BDNF groups was significantly higher than that in the Aβ group. Moreover, the duration of Ca2+ oscillations in the BDNF-treated group did not differ from that in the sham group on DIV 28. The parameters of spontaneous calcium activity in the “Aβ+AAV-Syn-EGFP” group did not differ from the values in the “Aβ” group.
Analysis of the main parameters of spontaneous Ca2+ activity identified significant changes in the functional state of neuronal cells beginning at DIV 21. The number of cells that exhibited Ca2+ activity in the Aβ+BDNF and Aβ+AAV-Syn-BDNF groups was significantly higher than that in the Aβ and Aβ+AAV-Syn-EGFP groups (DIV 21: sham: 88.27 ± 5.56%, Aβ: 37.74 ± 3.27%, Aβ+AAV-Syn-EGFP: 39.21 ± 5.97%; Aβ+AAV-Syn-BDNF: 77.83 ± 4.89%, Aβ+BDNF: 69.31 ± 4.35%). On DIV 28, the neuroprotective effect remained; the activity of cells in primary cultures treated with BDNF was significantly higher than that of those in AAV-Syn-EGFP-treated cultures and in cultures exposed to chronic Aβ (Figures 7B–D).
A decreased frequency of Ca2+ oscillations and an increased duration of Ca2+ events caused by chronic Aβ application were observed on DIV 21. Additionally, the frequency of Ca2+ oscillations in the Aβ+BDNF and Aβ+AAV-Syn-BDNF groups was significantly higher than that in the Aβ group group (DIV 21: Aβ, 1.14 ± 0.11 osc/min; Aβ+AAV-Syn-BDNF, 1.76 ± 0.17 osc/min; Aβ+BDNF, 2.38 ± 0.16 osc/min). The frequency of Ca2+ oscillations in the Aβ+AAV-Syn-EGFP group did not differ from that in the Aβ group (1.08 ± 0.27 osc/min). Moreover, the duration of Ca2+ oscillations in the BDNF-treated group did not differ from that in the sham group on DIV 28 (Figures 7C,D).
Thus, an AAV construct containing the BDNF gene, but not an AAV-Syn-EGFP virus vector, prevented impairments in functional neural network Ca2+ activity in our amyloidopathy model in vitro. The effect of BDNF hyperexpression was comparable to the effects of chronic recombinant protein application.
Electrophysiological data analysis revealed that Aβ disrupted neural network formation on DIV 14, and this effect was characterized by a significant decrease in the number of spikes in a small network burst (Table 3).
Table 3. Main parameters of spontaneous bioelectrical activity in primary hippocampal cell cultures under chronic Aβ application.
Further suppression of spontaneous bioelectrical activity in the cultures with chronic Aβ application was observed at a later stage. Compared to the sham group, the Aβ group exhibited fewer small network bursts and spikes in a burst on DIV 21 and DIV 28 (Table 3). BDNF partially negated the decreased spontaneous bioelectrical activity of primary hippocampal cultures. On DIV 21, the number of small network bursts in the Aβ+BDNF group was significantly higher than that in the Aβ group. Moreover, the number of spikes in a network burst in the Aβ+BDNF and Aβ+AAV-Syn-BDNF groups did not differ from that in the sham group (Table 3).
According to the classical concept, a network burst is an event comprising no fewer than four spikes simultaneously recorded from different electrodes in a 50-ms interval (Wagenaar et al., 2006; Pimashkin et al., 2011; Vedunova et al., 2013). However, for a more detailed analysis of the complex structure of the neural network, we selected events that simultaneously captured the prevailing part of functionally active cells. In this regard, all neural network bursts were conditionally divided into small (from 4 to 100 spikes in 50 ms) and large (101 or more spikes in 50 ms) groups (Mishchenko et al., 2019).
Our studies revealed that Aβ almost completely inhibited large burst formation in primary hippocampal cultures at early stages (Table 3C, Figure 8). Large network bursts were first measured in the sham group at DIV 14, whereas large network events in the Aβ group were completely absent on both DIV 14 and DIV 21 and did not begin to form until DIV 28. However, the formation of large network bursts was in cultures treated with recombinant BDNF and AAV-Syn-BDNF on DIV 21.
Figure 8. Number of spikes/50 ms and representative raster diagrams of spontaneous bioelectrical activity in primary hippocampal cultures under chronic Aβ application (DIV 21). The use of Aβ completely inhibited large network burst formation in cultures until DIV 28. In contrast, in cultures treated with recombinant BDNF and AAV-syn-BDNF the formation of complex network bursts were observed on DIV 21.
For a better understanding of neural network structure, we performed a cross-correlation analysis and correlation graph reconstruction to study the structure of functional interconnections in a network with a definition of its key activity elements–hubs (Shishkina et al., 2018). On DIV 14, a neuronal network in an intact primary hippocampal culture has a complex architecture, which includes several active hubs that form at least 10 connections between nearby electrodes (Figure 9). AD modeling by chronic application of synthetic Aβ led to significant simplification of the internal neural network structure on DIV 14. At this time, the number of connections between the network elements was 5 times lower than that in the sham group (Table 4). Simplification of the internal functional structure of neural networks in the Aβ group continued during culture development (Figure 9). Active centers and elements forming at least 10 connections were nearly absent on DIV 21. By DIV 28, neural network activity was almost completely abolished.
Figure 9. Internal functional structure of neural networks in the primary hippocampal cultures on DIV 21. Graphical representation of the correlated connections among neurons in the network. The electrode number is presented as “Cell X.” The number of connections on the electrode is indicated in square brackets. The vertex size is proportional to the number of significant connections. (A) Sham; (B) Aβ; (C) BNDF; (D) Aβ+AAV-syn-BDNF. A neuronal network in an intact primary hippocampal culture has a complex architecture, which includes several active hubs. Chronic application of synthetic Aβ simplified the internal neural network structure. BDNF preserved the complexity of neural network architecture and maintained the active centers in the network throughout the entire observation period.
Table 4. Average number of connections in a hub in primary hippocampal neural networks on different days of development in vitro.
Brain-derived neurotrophic factor preserved the complexity of neural network architecture. The active centers in the network were maintained throughout the entire observation period. The number of connections in the hubs in the BDNF-treated groups did not differ from the values in the sham group (DIV 21: sham, 16.5 ± 4.3; Aβ, 1.1 ± 1.3; Aβ+BDNF, 2.25 ± 4.1; Aβ+AAV-Syn-BDNF, 11.4 ± 3.2).
Next, we assumed that a molecular mechanism of BDNF neuroprotective action could be associated with activation of TrkB receptors, having a high affinity to BDNF. For this purpose, using RT-PCR, we performed quantitative evaluation of TrkB-FL receptors. We showed no significant alterations in TrkB-FL mRNA levels in primary hippocampal cultures on DIV21 under chronic Aβ1-42 application. The obtained data are presented in Supplementary Figure S2.
Discussion
Currently, animal models of AD are commonly used to study various changes associated with neurodegenerative processes. However, the ongoing in vivo studies do not allow for in-depth analysis of functional neural network activity. There have been a few experimental studies on neuronal cultures obtained from 5xFAD transgenic mice (ex. Park et al., 2016), but they did not investigate the features of functional neuron-glial network activity during AD development. Our data on functional calcium activity inhibition during the cultivation of primary hippocampal cells obtained from 5xFAD murine embryos are of interest since the proposed model may serve as a way to assess functional changes and identify mechanisms of neural network reorganization at the earliest stages of a familial form of AD. However, our model is not suitable for studying the effect of amyloidopathy.
In this study, we used experimental model of amyloidopathy to investigate the functional activity of neural networks in vitro. Amyloid protein oligomers are currently considered the most neurotoxic in AD; their accumulation is much better correlated with the severity of cognitive symptoms than the presence of plaques or neurofibrillary tangles (DiChiara et al., 2017; Bourdenx et al., 2017; Cline et al., 2018). Amyloid plaques mainly consist of aggregated Aβ, the most common forms of which are Aβ1-40 and Aβ1-42 (Scarano et al., 2016; Cline et al., 2018). Aβ accumulation initially occurs intraneuronally, mainly in synapses (Pignataro and Middei, 2017; Forner et al., 2017); plaques and tangles form in the brain parenchyma at later stages. Therefore, in our amyloidopathy model, we used Aβ1-42. We demonstrated neurodegenerative effects in this model, characterized by an increase in cell death, and immunocytochemically confirmed the formation of extracellular amyloid aggregates.
A pronounced suppressive action of chronic amyloidopathy on spontaneous functional Ca2+ activity in primary hippocampal cultures was also shown. The negative effect of Aβ was demonstrated by decreases in the number of cells exhibiting Ca2+ activity and frequency of Ca2+ oscillations. Studies on the features of Ca2+ network activity in AD are currently limited. Previous in vitro experiments have shown that Aβ oligomers can increase cytosolic Ca2+ levels, eventually leading to mitochondrial Ca2+ overload and partial neuronal death (Kuchibhotla et al., 2008; Villalobos et al., 2012; Tu et al., 2014). Interestingly, this effect of Aβ42 oligomer application was not revealed at the early stages of primary neuronal culture development (2-10 DIV), but Ca2+ activity was sharply increased beginning at DIV 14, at which point different age-related signs were observed in neurons. In this regard, the cytosolic and mitochondrial Ca2+ responses induced by single Aβ42 oligomer administration increased dramatically over time (adult cultures DIV 21-28) (Calvo-Rodríguez et al., 2016; Núñez et al., 2018). Our data indicated that chronic Aβ42 application increased the suppression of functional Ca2+ activity in primary hippocampal cultures and may correlate with in vivo studies of Ca2+ activity, indicating neural dysfunction in the brains of aging mice with modeled AD. These pathological changes result in the silencing of some neurons, whereas other neurons in areas enriched with Aβ plaques are hyperactive (Busche et al., 2012; Busche, 2018).
Currently, many research groups are interested in calcium activity changes during AD. However, there is still very limited data on the molecular mechanisms underlying these changes. One possible neurodegenerative mechanism in AD is through the hyperstimulation of glutamate receptors, mainly NMDA, resulting in an excessive increase in Ca2+ levels, causing excitotoxicity and further neuronal death (Kabir et al., 2019). Moreover, because numerous studies have demonstrated that Aβ protein oligomer stimulates calcium influx via the NMDA receptor in the pathogenesis of AD, NMDA receptor antagonists are considered potential therapeutic substances for presymptomatic AD (Kodis et al., 2018; Müller et al., 2018). Increased IP3R function is another expected cause of changes in intracellular transmission of Ca2+ signals (Mak et al., 2015). Overall, the aspect of calcium activity changes are extremely interesting and merit further detailed investigation.
We also studied the neuroprotective effects of the neurotrophic factor BDNF. We were especially interested in examining these effects in a model that permitted the evaluation of rapid neurodegeneration. Therefore, we used the model based on chronic Aβ application.
Cognitive dysfunction in AD is associated with impairments in neurotrophic factors [such as BDNF, glial cell line-derived neurotrophic factor (GDNF), and nerve growth factor (NGF)] levels in the blood (Budni et al., 2015). Yasutake et al., 2006, showed a decrease in BDNF levels in blood in the late stages of AD. Additionally, Aβ has been shown to directly inhibit BDNF proteolysis from pro-BDNF (Zheng et al., 2010) and reduce the retrograde axonal transport of the BDNF/TrkB system through a mechanism including ubiquitin carboxy-terminal hydrolase L1 (Poon et al., 2013). Thus, a decrease in BDNF levels is likely one of the important molecular mechanisms of AD pathogenesis. In this regard, we considered BDNF a potential neuroprotectant in amyloidosis-induced neurodegenerative processes.
As neurodegeneration has a chronic irreversible pattern in AD, we studied the effects of constant BDNF hyperexpression induced by an AVV vector carrying the BDNF gene. Chronic application of recombinant BDNF protein (1 ng/ml) served as a positive control. Similar viral constructs carrying neurotrophic factor genes are considered promising agents that are being actively studied as gene therapies for neurodegenerative diseases, such as Parkinson’s disease (Cheng et al., 2018; Tereshchenko et al., 2014). However, there is a lack of data on the effectiveness of these approaches in AD. For instance, using a transgenic model of amyloidosis in vivo, the team led by Prof. M.H. Tuszynski has shown that viral delivery of BDNF gene after the beginning of pathological processes development restrains the loss of synapses, partially normalizes the aberrant expression of the Aβ gene, improves synaptic transmission and restores learning and memory (Nagahara et al., 2009). Jiao SS. et al. demonstrated the influence of AAV-BDNF construct application in vivo (Jiao et al., 2016). The recovery of the BDNF level attenuated behavioral impairments and prevented neuronal loss but did not affect the level of tau hyperphosphorylation in P301L mouse brains. These data indicate that delivery of the BDNF gene is a promising method for neurodegeneration therapy. However, there have been no studies of neural network activity following treatment with genetically engineered constructs in AD.
Our in vitro studies using a β-amyloidopathy model showed a pronounced neuroprotective effect of BDNF. Notably, we reported cell viability preservation as well as maintenance of the functional integrity of neural networks, characterized by the normalization of spontaneous Ca2+ activity.
Electrophysiological data analysis revealed significant impairments in primary hippocampal neural network activity under chronic Aβ application. Our previous studies have shown several stages of functional activity formation in neuronal cultures. The first single small network bursts are detected on DIV 7; their numbers and the number of spikes in a burst gradually increase by DIV 14. Stable spontaneous bioelectrical activity in primary hippocampal cultures is measured beginning at DIV 14 (Shirokova et al., 2013; Mishchenko et al., 2019).
Our present study revealed that Aβ decreased the number of small network bursts and spikes in a burst on DIV 14 and resulted in a lack of large network events up to DIV 28. A significant simplification in the functional neural network architecture was also observed. These changes are potentially associated with neurodegenerative processes that occur in primary neuronal cultures, with a reduction in cellular outgrowth and synaptic loss, as well as the death of some functionally significant network elements. Changes in spontaneous bioelectrical neural network activity in AD have been very poorly investigated. The available data were mainly obtained by the patch-clamp method and describe the features of synaptic transmission in individual neurons and synaptic endings (Lasala et al., 2019; Mondragón-Rodríguez et al., 2018), and studies using multichannel electroencephalography only indirectly describe neural network activity in AD models (Ahnaou et al., 2017).
Therefore, the electrophysiological data presented in the manuscript characterizing the network activity and functional architecture of neural networks in the AD model are of great interest and may be considered a fundamental basis for the further study of neural network activity, including the application of different stimulation protocols.
A recently published study based on MEA recordings showed that one application of Aβ oligomers initially led to activation followed by further (after 12-24 h) suppression and desynchronization in neural network activity, ultimately resulting in network destruction (Gao et al., 2019). In our work, we studied more prolonged effects caused by amyloid on neural network activity and identified a significant decrease. Our results are consistent with previous data from in vivo experiments and studies performed on brain slices obtained from different AD animal models, indicating that the early preclinical stages of AD are characterized by hyperactivation of neuronal activity, whereas the later stages are characterized by its suppression (de Haan et al., 2017; Heggland et al., 2019).
In could be assumed that activation of TrkB receptors is a backbone pathway for BDNF neuroprotective action, and alterations in receptors functions are one of the pathogenetic aspects of AD. A number of studies have shown a decrease in TkrB-FL expression in both in vitro and in vivo models of AD (Kemppainen et al., 2012; Lei et al., 2018). For instance, rat intracerebroventricular injection of Aβ25-35 solution (5 μL) was shown to result in a significant decrease in TrkB levels (Wang K. et al., 2018). Immunohistochemical studies of postmortem tissues showed a significantly lower staining density of TrkB receptor in the hippocampus of patients with AD than in that of healthy controls. In addition, TrkB staining was inversely correlated with both Amylo-Glo and pTau staining in the same region. These observations strongly confirm that changes in the BDNF-TrkB system are involved in the pathology of AD (Bharani et al., 2019). However, we showed no significant alterations in TrkB-FL mRNA levels in primary hippocampal cultures 2 weeks after the beginning of Aβ application.
In summary, in this work, we examined two experimental models of amyloidopathy using primary hippocampal cultures, one employing chronic application of Aβ1-42 and the other using the embryonic brains of 5xFAD mice. Chronic application of Aβ1-42 resulted in the rapid establishment of significant neurodegenerative changes in primary hippocampal cultures, leading to marked impairments in neural network calcium activity and increased cell death. We studied the influence of amyloidopathy on spontaneous bioelectrical neural network activity in primary hippocampal cultures. Chronic Aβ application led to a decrease in the number of network bursts and spikes in a burst and disturbed the spatial structure of neural networks, reducing the number of key network elements (hubs) and the number of connections between network elements. The model based on primary hippocampal cells obtained from 5xFAD mice demonstrated changes in spontaneous network calcium activity characterized by a decrease in the number of cells exhibiting Ca2+ activity, a decrease in the frequency of Ca2+ oscillations and an increase in the duration of Ca2+ events beginning at DIV 21.
Moreover, application of BDNF recombinant protein and BDNF hyperexpression by an AAV vector partially prevented these amyloidopathy-induced neurodegenerative phenomena. BDNF maintained cell viability and spontaneous bioelectrical and calcium network activity in primary hippocampal cultures. The internal functional structure of neural networks, including the number of hubs and connections between active elements in the network, was also partially preserved.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
Ethics Statement
The animal study was reviewed and approved by the Bioethics Committee of Lobachevsky University.
Author Contributions
EM, TM, and RY carried out the experiments on primary cultures with support from VK. EM and MV wrote the manuscript with input from all other authors. AB and MG designed and developed the virus. EM, TM, and MV ensured the financing of the project. MV supervised the project, conceptualized the original idea, and was in charge of the overall direction of the study. EE performed Real-Time PCR analysis. All authors read and approved the final manuscript.
Funding
This work was supported by the state projects “Provision of Scientific Research” (Nos. 6.6379.2017/8.9, 17.3335.2017/4.6, 6.6659.2017/6.7, and 0729-2020-0061) and partially by RFBR (Project Nos. 18-015-00391 and 18-315-20003) and Grant of the President of the Russian Federation (MK-1485.2019.4).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank Alexey S. Pimashkin for technical support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2020.00582/full#supplementary-material
FIGURE S1 | Representative light field images of primary hippocampal cultures obtained from wild-type and 5xFAD murine embryos (DIV 21), Scale bar – 20 μm. Comparative morphological assessment did not reveal significant changes between cultures obtained from wild-type and 5xFAD murine embryos over 28 DIV.
FIGURE S2 | Features of TrkB-FL gene expression on DIV 21 under chronic exogenous administration of Aβ and BDNF. Data are normalized to the reference gene (Oaz1). The data represent the mean values ± SEMs from three independent experiments. We showed no significant alterations in TrkB-FL mRNA levels in primary hippocampal cultures on DIV21 under chronic Aβ application.
References
Ahnaou, A., Moechars, D., Raeymaekers, L., Biermans, R., Manyakov, N. V., Bottelbergs, A., et al. (2017). Emergence of early alterations in network oscillations and functional connectivity in a tau seeding mouse model of Alzheimer’s disease pathology. Sci. Rep. 7:14189. doi: 10.1038/s41598-017-13839-6
Banik, A., Brown, R. E., Bamburg, J., Lahiri, D. K., Khurana, D., Friedland, R. P., et al. (2015). Translation of pre-clinical studies into successful clinical trials for Alzheimer’s disease: what are the roadblocks and how can they be overcome? J. Alzheimers Dis. 47, 815–843. doi: 10.3233/JAD-150136
Bharani, K. L., Ledreux, A., Gilmore, A., Carroll, S. L., and Granholm, A. C. (2019). Serum pro-BDNF levels correlate with phospho-tau staining in Alzheimer’s disease. Neurobiol. Aging 87, 49–59. doi: 10.1016/j.neurobiolaging.2019.11.010
Bilkei-Gorzo, A. (2014). Genetic mouse models of brain ageing and Alzheimer’s disease. Pharmacol. Ther. 142, 244–257. doi: 10.1016/j.pharmthera.2013.12.009
Bourdenx, M., Koulakiotis, N. S., Sanoudou, D., Bezard, E., Dehay, B., and Tsarbopoulos, A. (2017). Protein aggregation and neurodegeneration in prototypical neurodegenerative diseases: examples of amyloidopathies, tauopathies and synucleinopathies. ProgNeurobiol 155, 171–193. doi: 10.1016/j.pneurobio.2015.07.003
Budni, J., Bellettini-Santos, T., Mina, F., Garcez, M. L., and Zugno, A. I. (2015). The involvement of BDNF, NGF and GDNF in aging and Alzheimer’s disease. Aging Dis. 6, 331–341. doi: 10.14336/AD.2015.0825
Busche, M. A. (2018). In vivo two-photon calcium imaging of hippocampal neurons in Alzheimer mouse models. Methods Mol. Biol. 1750, 341–351. doi: 10.1007/978-1-4939-7704-8_23
Busche, M. A., Chen, X., Henning, H. A., Reichwald, J., Staufenbiel, M., Sakmann, B., et al. (2012). Critical role of soluble amyloid-β for early hippocampal hyperactivity in a mouse model of Alzheimer’s disease. Proc. Natl. Acad. Sci. U.S.A. 109, 8740–8745. doi: 10.1073/pnas.1206171109
Calvo-Rodríguez, M., García-Durillo, M., Villalobos, C., and Núñez, L. (2016). Aging enables Ca2+ overload and apoptosis induced by amyloid-β oligomers in rat hippocampal neurons: neuroprotection by non-steroidal anti-inflammatory drugs and R-flurbiprofen in aging neurons. J. Alzheimers Dis. 54, 207–221. doi: 10.3233/JAD-151189
Cao, J., Hou, J., Ping, J., and Cai, D. (2018). Advances in developing novel therapeutic strategies for Alzheimer’s disease. Mol. Neurodegener. 13:64. doi: 10.1186/s13024-018-0299-8
Cheng, S., Tereshchenko, J., Zimmer, V., Vachey, G., Pythoud, C., Rey, M., et al. (2018). Therapeutic efficacy of regulable GDNF expression for Huntington’s and Parkinson’s disease by a high-induction, background-free “GeneSwitch” vector. Exp. Neurol. 309, 79–90. doi: 10.1016/j.expneurol.2018.07.017
Choi, S. H., Bylykbashi, E., Chatila, Z. K., Lee, S. W., Pulli, B., Clemenson, G. D., et al. (2018). Combined adult neurogenesis and BDNF mimic exercise effects on cognition in an Alzheimer’s mouse model. Science 361:eaan8821. doi: 10.1126/science.aan8821
Cline, E. N., Bicca, M. A., Viola, K. L., and Klein, W. L. (2018). The amyloid-β oligomer hypothesis: beginning of the third decade. J. Alzheimers Dis. 64(Suppl. 1), S567–S610. doi: 10.3233/JAD-179941
Criscuolo, C., Fabiani, C., Bonadonna, C., Origlia, N., and Domenici, L. (2015). BDNF prevents amyloid-dependent impairment of LTP in the entorhinal cortex by attenuating p38 MAPK phosphorylation. Neurobiol. Aging 36, 1303–1309. doi: 10.1016/j.neurobiolaging.2014.11.016
Cummings, J., Lee, G., Ritter, A., and Zhong, K. (2018). Alzheimer’s disease drug development pipeline: 2018. Alzheimers Dement 4, 195–214. doi: 10.1016/j.trci.2018.03.009
de Pins, B., Cifuentes-Díaz, C., Farah, A. T., López-Molina, L., Montalban, E., Sancho-Balsells, A., et al. (2019). Conditional BDNF delivery from astrocytes rescues memory deficits, spine density, and synaptic properties in the 5xFAD mouse model of Alzheimer disease. J. Neurosci. 39, 2441–2458. doi: 10.1523/JNEUROSCI.2121-18.2019
de Haan, W., van Straaten, E. C. W., Gouw, A. A., and Stam, C. J. (2017). Altering neuronal excitability to preserve network connectivity in a computational model of Alzheimer’s disease. PLoS Comput. Biol. 13:e1005707. doi: 10.1371/journal.pcbi.1005707
DiChiara, T., DiNunno, N., Clark, J., Bu, R. L., Cline, E. N., Rollins, M. G., et al. (2017). Alzheimer’s toxic amyloid beta oligomers: unwelcome visitors to the Na/K ATPase alpha3 docking station. Yale J. Biol. Med. 90, 45–61.
Drummond, E., and Wisniewski, T. (2017). Alzheimer’s disease: experimental models and reality. Acta Neuropathol. 133, 155–175. doi: 10.1007/s00401-016-1662-x
Dubois, B., Feldman, H. H., Jacova, C., Hampel, H., Molinuevo, J. L., Blennow, K., et al. (2014). Advancing research diagnostic criteria for Alzheimer’s disease: the IWG-2 criteria. Lancet Neurol. 6, 614–629. doi: 10.1016/S1474-4422(14)70090-0
Forner, S., Baglietto-Vargas, D., Martini, A. C., Trujillo-Estrada, L., and LaFerla, F. M. (2017). Synaptic impairment in Alzheimer’s disease: a dysregulated symphony. Trends Neurosci. 40, 347–357. doi: 10.1016/j.tins.2017.04.002
Gao, F., Gao, K., He, C., Liu, M., Wan, H., and Wang, P. (2019). Multi-site dynamic recording for Aβ oligomers-induced Alzheimer’s disease in vitro based on neuronal network chip. Biosens. Bioelectron. 133, 183–191. doi: 10.1016/j.bios.2019.03.025
Hadar, A., and Gurwitz, D. (2018). Peripheral transcriptomic biomarkers for early detection of sporadic Alzheimer disease? Dialogues Clin. Neurosci. 4, 293–300. doi: 10.31887/dcns.2018.20.4/dgurwitz
Hampel, H., Schneider, L. S., Giacobini, E., Kivipelto, M., Sindi, S., Dubois, B., et al. (2015). Advances in the therapy of Alzheimer’s disease: targeting amyloid beta and tau and perspectives for the future. Expert Rev. Neurother. 1, 83–105. doi: 10.1586/14737175.2015.995637
Hasan, M. F., and Berdichevsky, Y. (2016). Neural circuits on a chip. Micromachines (Basel) 7:E157. doi: 10.3390/mi7090157
Heggland, I., Kvello, P., and Witter, M. P. (2019). Electrophysiological characterization of networks and single cells in the hippocampal region of a transgenic rat model of Alzheimer’s disease. eNeuro 6:ENEURO.0448-17.2019. doi: 10.1523/ENEURO.0448-17.2019
Jiao, S. S., Shen, L. L., Zhu, C., Bu, X. L., Liu, Y. H., Liu, C. H., et al. (2016). Brain-derived neurotrophic factor protects against tau-related neurodegeneration of Alzheimer’s disease. Transl. Psychiatry 6:e907. doi: 10.1038/tp.2016.186
Johnstone, A. F., Gross, G. W., Weiss, D. G., Schroeder, O. H., Gramowski, A., and Shafer, T. J. (2010). Microelectrode arrays: a physiologically based neurotoxicity testing platform for the 21st century. Neurotoxicology 4, 331–350. doi: 10.1016/j.neuro.2010.04.001
Kabir, M. T., Sufian, M. A., Uddin, M. S., Begum, M. M., Akhter, S., Islam, A., et al. (2019). NMDA receptor antagonists: repositioning of memantine as a multitargeting agent for Alzheimer’s therapy. Curr. Pharm. Des. 25, 3506–3518. doi: 10.2174/138161282566619101110244
Karch, C. M., Cruchaga, C., and Goate, A. M. (2014). Alzheimer’s disease genetics: from the bench to the clinic. Neuron 83, 11–26. doi: 10.1016/j.neuron.2014.05.041
Kemppainen, S., Rantamäki, T., Jerónimo-Santos, A., Lavasseur, G., Autio, H., Karpova, N., et al. (2012). Impaired TrkB receptor signaling contributes to memory impairment in APP/PS1 mice. Neurobiol. Aging 33:1122.e23-39. doi: 10.1016/j.neurobiolaging.2011.11.006
Kim, D. H., Yeo, S. H., Park, J. M., Choi, J. Y., Lee, T. H., Park, S. Y., et al. (2014). Genetic markers for diagnosis and pathogenesis of Alzheimer’s disease. Gene 545, 185–193. doi: 10.1016/j.gene.2014.05.031
Kodis, E. J., Choi, S., Swanson, E., Ferreira, G., and Bloom, G. S. (2018). N-methyl-D-aspartate receptor-mediated calcium influx connects amyloid-β oligomers to ectopic neuronal cell cycle reentry in Alzheimer’s disease. Alzheimers Dement. 14, 1302–1312. doi: 10.1016/j.jalz.2018.05.017
Kuchibhotla, K. V., Goldman, S. T., Lattarulo, C. R., Wu, H. Y., Hyman, B. T., and Bacskai, B. J. (2008). A plaques lead to aberrant regulation of calcium homeostasis in vivo resulting in structural and functional disruption of neuronal networks. Neuron 59, 214–225. doi: 10.1016/j.neuron.2008.06.008
Lasala, M., Fabiani, C., Corradi, J., Antollini, S., and Bouzat, C. (2019). Molecular modulation of human α7 nicotinic receptor by Amyloid-β peptides. Front. Cell. Neurosci. 13:37. doi: 10.3389/fncel.2019.00037
Lei, H., Zhang, Y., Huang, L., Xu, S., Li, J., Yang, L., et al. (2018). L-3-n-butylphthalide regulates proliferation, migration, and differentiation of neural stem cell in vitro and promotes neurogenesis in APP/PS1 Mouse model by regulating BDNF/TrkB/CREB/Akt pathway. Neurotox. Res. 34, 477–488. doi: 10.1007/s12640-018-9905-3
Lim, S., Airavaara, M., and Harvey, B. K. (2010). Viral vectors for neurotrophic factor delivery: a gene therapy approach for neurodegenerative diseases of the CNS. Pharmacol. Res. 61, 14–26. doi: 10.1016/j.phrs.2009.10.002
Liu, P., Reichl, J. H., Rao, E. R., McNellis, B. M., Huang, E. S., Hemmy, L. S., et al. (2017). Quantitative comparison of dense-core amyloid plaque accumulation in amyloid-β precursor protein transgenic mice. J. Alzheimers Dis. 56, 743–761. doi: 10.3233/JAD-161027
Mak, D. O., Cheung, K. H., Toglia, P., Foskett, J. K., and Ullah, G. (2015). Analyzing and quantifying the gain-of-function enhancement of IP3 receptor gating by familial Alzheimer’s disease-causing mutants in presenilins. PLoS Comput. Biol. 11:e1004529. doi: 10.1371/journal.pcbi.1004529
Mango, D., Saidi, A., Cisale, G. Y., Feligioni, M., Corbo, M., and Nisticò, R. (2019). Targeting synaptic plasticity in experimental models of Alzheimer’s disease. Front Pharmacol. 10:778. doi: 10.3389/fphar.2019.00778
Mishchenko, T. A., Mitroshina, E. V., Usenko, A. V., Voronova, N. V., Astrakhanova, T. A., Shirokova, O. M., et al. (2019). Features of neural network formation and their functions in primary hippocampal cultures in the context of chronic TrkB receptor system influence. Front. Physiol. 9:1925. doi: 10.3389/fphys.2018.01925
Mitroshina, E. V., Mishchenko, T. A., Usenko, A. V., Epifanova, E. A., Yarkov, R. S., Gavrish, M. S., et al. (2018). AAV-Syn-BDNF-EGFP virus construct exerts neuroprotective action on the hippocampal neural network during hypoxia in vitro. Int. J. Mol. Sci. 19:2295. doi: 10.3390/ijms19082295
Mondragón-Rodríguez, S., Gu, N., Manseau, F., and Williams, S. (2018). Alzheimer’s transgenic model is characterized by very early brain network alterations and β-CTF fragment accumulation: reversal by β-secretase inhibition. Front. Cell. Neurosci. 12:121. doi: 10.3389/fncel.2018.00121
Mucke, L., Masliah, E., Yu, G. Q., Mallory, M., Rockenstein, E. M., Tatsuno, G., et al. (2000). High-level neuronal expression of aβ 1-42 in wildtype human amyloid protein precursor transgenic mice: synaptotoxicity without plaque formation. J. Neurosci. 20, 4050–4058. doi: 10.1523/jneurosci.20-11-04050.2000
Müller, M. K., Jacobi, E., Sakimura, K., Malinow, R., and von Engelhardt, J. (2018). NMDA receptors mediate synaptic depression, but not spine loss in the dentate gyrus of adult amyloid Beta (Aβ) overexpressing mice. Acta Neuropathol. Commun. 6:110. doi: 10.1186/s40478-018-0611-4
Nagahara, A. H., Mateling, M., Kovacs, I., Wang, L., Eggert, S., Rockenstein, E., et al. (2013). Early BDNF treatment ameliorates cell loss in the entorhinal cortex of APP transgenic mice. J. Neurosci. 33, 15596–15602. doi: 10.1523/JNEUROSCI.5195-12.2013
Nagahara, A. H., Merrill, D. A., Coppola, G., Tsukada, S., Schroeder, B. E., Shaked, G. M., et al. (2009). Neuroprotective effects of brain-derived neurotrophic factor in rodent and primate models of Alzheimer’s disease. Nat. Med. 15, 331–337. doi: 10.1038/nm.1912
Núñez, L., Calvo-Rodríguez, M., Caballero, E., García-Durillo, M., and Villalobos, C. (2018). Neurotoxic Ca2+ signaling induced by amyloid-β oligomers in aged hippocampal neurons in vitro. Methods Mol. Biol. 1779, 341–354. doi: 10.1007/978-1-4939-7816-8_20
Oakley, H., Cole, S. L., Logan, S., Maus, E., Shao, P., Craft, J., et al. (2006). Intraneuronal beta-amyloid aggregates, neurodegeneration, and neuron loss in transgenic mice with five familial Alzheimer’s disease mutations: potential factors in amyloid plaque formation. J. Neurosci. 26, 10129–10140. doi: 10.1523/JNEUROSCI.1202-06.2006
Park, S. E., Lee, J., Chang, E. H., Kim, J. H., Sung, J. H., Na, D. L., et al. (2016). Activin A secreted by human mesenchymal stem cells induces neuronal development and neurite outgrowth in an in vitro model of Alzheimer’s disease: neurogenesis induced by MSCs via activin A. Arch. Pharm. Res. 39, 1171–1179. doi: 10.1007/s12272-016-0799-4
Peng, S., Wuu, J., Mufson, E. J., and Fahnestock, M. (2005). Precursor form of brain-derived neurotrophic factor and mature brain-derived neurotrophic factor are decreased in the pre-clinical stages of Alzheimer’s disease. J. Neurochem. 93, 1412–1421. doi: 10.1111/j.1471-4159.2005.03135.x
Pignataro, A., and Middei, S. (2017). Trans-synaptic spread of amyloid-β in alzheimer’s disease: paths to β-amyloidosis. Neural. Plast. 2017:5281829. doi: 10.1155/2017/5281829
Pimashkin, A., Kastalskiy, I., Simonov, A., Koryagina, E., Mukhina, I., and Kazantsev, V. (2011). Spiking signatures of spontaneous activity bursts in hippocampal cultures. Front. Comput. Neurosci. 5:46. doi: 10.3389/fncom.2011.00046
Poon, W. W., Carlos, A. J., Aguilar, B. L., Berchtold, N. C., Kawano, C. K., Zograbyan, V., et al. (2013). β-Amyloid (Aβ) oligomers impair brain-derived neurotrophic factor retrograde trafficking by down-regulating ubiquitin C-terminal hydrolase, UCH-L1. J. Biol. Chem. 288, 169–182. doi: 10.1074/jbc.M113.463711
Radde, R., Bolmont, T., Kaeser, S. A., Coomaraswamy, J., Lindau, D., Stoltze, L., et al. (2006). Aβ42-driven cerebral amyloidosis in transgenic mice reveals early and robust pathology. EMBO Rep. 7, 940–946. doi: 10.1038/sj.embor.7400784
Scarano, S., Lisi, S., Ravelet, C., Peyrin, E., and Minunni, M. (2016). Detecting Alzheimer’s disease biomarkers: from antibodies to new bio-mimetic receptors and their application to established and emerging bioanalytical platforms – a critical review. Anal. Chim. Acta 940, 21–37. doi: 10.1016/j.aca.2016.08.008
Shirokova, O. M., Frumkina, L. E., Vedunova, M. V., Mitroshina, E. V., Zakharov, Y. N., Khaspekov, L., et al. (2013). Morphofunctional patterns of neuronal network developing in dissociated hippocampal cell cultures. Sovrem. Tehnol. Med. 5, 6–13.
Shishkina, T. V., Mishchenko, T. A., Mitroshina, E. V., Shirokova, O. M., Pimashkin, A. S., Kastalskiy, I. A., et al. (2018). Glial cell line-derived neurotrophic factor (GDNF) counteracts hypoxic damage to hippocampal neural network function in vitro. Brain Res. 1678, 310–321. doi: 10.1016/j.brainres.2017.10.023
Stancu, I. C., Vasconcelos, B., Terwel, D., and Dewachter, I. (2014). Models of β-amyloid induced Tau-pathology: the long and “folded” road to understand the mechanism. Mol. Neurodegener. 9:51. doi: 10.1186/1750-1326-9-51
Tereshchenko, J., Maddalena, A., Bähr, M., and Kügler, S. (2014). Pharmacologically controlled, discontinuous GDNF gene therapy restores motor function in a rat model of Parkinson’s disease. Neurobiol. Dis. 65, 35–42. doi: 10.1016/j.nbd.2014.01.009
Tu, S., Okamoto, S., Lipton, S. A., and Xu, H. (2014). Oligomeric Aβ-induced synaptic dysfunction in Alzheimer’s disease. Mol. Neurodegener. 9:48. doi: 10.1186/1750-1326-9-48
Vedunova, M., Sakharnova, T., Mitroshina, E., Perminova, M., Pimashkin, A., Zakharov, Y., et al. (2013). Seizure-like activity in hyaluronidase-treated dissociated hippocampal cultures. Front. Cell. Neurosci. 7:149. doi: 10.3389/fncel.2013.00149
Vedunova, M. V., Mishchenko, T. A., Mitroshina, E. V., and Mukhina, I. V. (2015). TrkB-mediated neuroprotective and antihypoxic properties of Brain-derived neurotrophic factor. Oxid. Med. Cell. Longev. 2015:453901. doi: 10.1155/2015/453901
Villalobos, C., Caballero, E., Sanz-Blasco, S., and Nunez, L. (2012). Study of neurotoxic intracellular calcium signalling triggered by amyloids. Methods Mol. Biol. 849, 289–302. doi: 10.1007/978-1-61779-551-0_20
Villalobos Acosta, D. M. Á, Chimal Vega, B., Correa Basurto, J., Fragoso Morales, L. G., and Rosales Hernández, M. C. (2018). Recent advances by in silico and in vitro studies of amyloid-β 1-42 fibril depicted a s-shape conformation. Int. J. Mol. Sci. 19:2415. doi: 10.3390/ijms19082415
Wagenaar, D. A., Pine, J., and Potter, S. M. (2006). An extremely rich repertoire of bursting patterns during the development of cortical cultures. BMC Neurosci. 7:11. doi: 10.1186/1471-2202-7-11
Walker, L. C., Diamond, M. I., Duff, K. E., and Hyman, B. T. (2013). Mechanisms of protein seeding in neurodegenerative diseases. JAMA Neurol. 70, 304–310. doi: 10.1001/jamaneurol.2013.1453
Wang, K., Sun, W., Zhang, L., Guo, W., Xu, J., Liu, S., et al. (2018). Oleanolic acid ameliorates Aβ25-35 injection-induced memory deficit in Alzheimer’s disease model rats by maintaining synaptic plasticity. CNS Neurol. Disord. Drug Targets 17, 389–399. doi: 10.2174/1871527317666180525113109
Wang, N., Qiu, P., Cui, W., Zhang, B., Yan, X., and He, S. (2018). Recent advances in multi-target anti-Alzheimer disease compounds (2013 up to present). Curr. Med. Chem. 26, 5684–5710. doi: 10.2174/0929867326666181203124102
Wang, R., and Holsinger, R. (2018). Exercise-induced brain-derived neurotrophic factor expression: therapeutic implications for Alzheimer’s dementia. Ageing Res. Rev. 46, 109–121. doi: 10.1016/j.arr.2018.10.002
Webster, S. J., Bachstetter, A. D., Nelson, P. T., Schmitt, F. A., and Van Eldik, L. J. (2014). Using mice to model Alzheimer’s dementia: an overview of the clinical disease and the preclinical behavioral changes in 10 mouse models. Front. Genet. 5:88. doi: 10.3389/fgene.2014.00088
Yasutake, C., Kuroda, K., Yanagawa, T., Okamura, T., and Yoneda, H. (2006). Serum BDNF, TNF-alpha and IL-1beta levels in dementia patients: comparison between Alzheimer’s disease and vascular dementia. Eur. Arch. Psychiatry Clin. Neurosci. 256, 402–406. doi: 10.1007/s00406-006-0652-8
Yuste, R. (2015). From the neuron doctrine to neural networks. Nat. Rev. Neurosci. 16, 487–497. doi: 10.1038/nrn3962
Zakharov, Y. N., Mitroshina, E. V., Shirokova, O., and Mukhina, I. V. (2013). “Calcium transient imaging as tool for neuronal and glial network interaction study,” in Models, Algorithms, and Technologies for Network Analysis. Springer Proceedings in Mathematics & Statistics, Vol. 32, eds B. Goldengorin, V. Kalyagin, and P. Pardalos, (New York, NY: Springer), 225–232. doi: 10.1007/978-1-4614-5574-5_12
Keywords: neural networks, Alzheimer’s disease, β-amyloidopathy, brain-derived neurotrophic factor, microelectrode arrays, calcium imaging, neuroprotection
Citation: Mitroshina EV, Yarkov RS, Mishchenko TA, Krut’ VG, Gavrish MS, Epifanova EA, Babaev AA and Vedunova MV (2020) Brain-Derived Neurotrophic Factor (BDNF) Preserves the Functional Integrity of Neural Networks in the β-Amyloidopathy Model in vitro. Front. Cell Dev. Biol. 8:582. doi: 10.3389/fcell.2020.00582
Received: 02 November 2019; Accepted: 16 June 2020;
Published: 08 July 2020.
Edited by:
Adelaide Fernandes, University of Lisbon, PortugalReviewed by:
Sandra Henriques Vaz, University of Lisbon, PortugalLaurent Roybon, Lund University, Sweden
Copyright © 2020 Mitroshina, Yarkov, Mishchenko, Krut’, Gavrish, Epifanova, Babaev and Vedunova. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elena V. Mitroshina, aGVsZW5taXRyb3NoaW5hQGdtYWlsLmNvbQ==; Maria V. Vedunova, bXZlZHVub3ZhQHlhbmRleC5ydQ==