- Open Access
Serotonin 5-HT6 receptors affect cognition in a mouse model of Alzheimer’s disease by regulating cilia function
Alzheimer's Research & Therapyvolume 9, Article number: 76 (2017)
Serotonin receptor 5-HT6 is involved in cognition and Alzheimer’s disease (AD) development. However, the mechanism of 5-HT6 in AD pathology is not clear.
Since 5-HT6 is almost exclusively expressed in the primary cilia, using immunostaining we examined the number of cilia in the hippocampus of AD animal model APP/PS1 mice. By overexpressing and knocking down 5-HT6 in the primary cultured hippocampal neurons, we investigated the roles of 5-HT6 in alternating ciliary morphology. Furthermore, 5-HT6 antagonist was applied to confirm its roles in cognition using the Morris water maze test, Y maze, and fear conditioning.
In the present study, we found that the primary cilia were elongated in the hippocampus of APP/PS1 mice compared with WT mice. 5-HT6 regulated cilia length, influenced cilia and axon initial segment (AIS) morphology, and affected localization of ARL13B and AnkG. We also found that, by changing cilia morphology, the AIS was elongated, branched, and more proximal to the cell body in both WT and APP/PS1 mouse neurons. Alterations of cilia also decreased the axonal length in WT and APP/PS1 neurons. Furthermore, in the water maze test, Y maze, and fear conditioning test, 5-HT6 antagonist SB271046 recovered the cognitive impairment of APP/PS1 mice.
We suggest that 5-HT6 plays a critical role in AD development through regulating the morphology and function of neuronal primary cilia, which is possibly related to the AIS and axon alterations in AD development.
Alzheimer’s disease (AD) is the most common form of dementia in people aged 65 years and older. AD is associated with impairments in memory, cognition, language, behavior, and personality . The pathology of AD in the brain includes extracellular senile plaques composed of amyloid beta (Aβ), intracellular tau aggregates known as neurofibrillary tangles (NFT), and neuronal and synaptic loss . Many studies have shown that individuals with the 267C allele of 5-HT6 have increased risk of AD [2, 3]. Several studies have reported that 5-HT6 has great impact in cognition, especially in AD models . Up to now, there is no effective treatment for AD. During the past 30 years, many studies have indicated that 5-HT6 may be a potential target for cognitive improvement in AD . Some 5-HT6 antagonists are currently being tested in clinical trials . However, although there is a strong link between 5-HT6 and AD development and treatment, the mechanism of how 5-HT6 is associated with AD is still not well understood. Besides, several 5-HT6 antagonist drugs, such as Pfizer's SAM-760 and SAM-531, Abbott's A-964324, Avineuro Pharmaceuticals' AVN-322, Suven's SUVN-502, Synosia Therapeutics' SYN-120, idalorpidine, and others, showed no significant positive effects for AD patients . The clinical failures made the pathophysiology of 5-HT6 antagonist more uncertain and indicated the necessity to clarify the inner function of 5-HT6 in AD. In the present study, we aimed to determine the mechanism of 5-HT6 in AD pathology.
Previous studies have shown that 5-HT6 can influence dendritic growth and neuronal migration [8, 9], suggesting that 5-HT6 might influence neuritic and axonal development. Furthermore, it has been reported that AD is associated with axonal abnormalities, including axonal swelling and abnormal accumulation of axonal proteins, which are related to dystrophic neurites and amyloid plaque formation . Abundant age-dependent axonal spheroids and myelin ovoids have been reported in the spinal cord of APPSwe/London/PS1M146V mice, an established familial AD (FAD) transgenic mouse model . However, the causative association of axonopathy to AD and the pathways involved in AD-related axonal abnormalities are not defined. The study of 5-HT6 promotes research on axonal abnormalities in AD. Furthermore, a previous study indicated that BBS8, a component of the primary cilia, influences axon targeting . The hedgehog signaling pathway, which is active in primary cilia, regulates axonal length in hippocampal neurons . Therefore, the study of 5-HT6 in APP/PS1 mice would contribute to exploring the relationship between primary cilia and axons.
5-HT6 is highly expressed in the central nervous system (CNS) [6, 14]. High levels of 5-HT6 expression have also been reported in the striatum, olfactory tubercle, nucleus accumbens, cortex, and hippocampus [14, 15]. 5-HT6 mRNA is exclusively expressed in neurons . In neurons, 5-HT6 is especially localized at the primary cilia [16, 17]. Primary cilia are protrusions of the cell membrane that are approximately 1–5 μm in length and anchored to the ciliary basal body (a modified mother centriole) . Cilia play critical roles in signaling in the neuronal system. For example, in the absence of the hedgehog signaling pathway ligand sonic hedgehog (shh), Patched receptors located at the cilia prevent Smoothened receptors from entering the cilia. In the presence of shh, Patched receptors move away from the cilia, allowing Smoothened receptors to accumulate in the cilia [19, 20]. Various human diseases are related to ciliary dysfunction, including cystic disease of the kidney, polydactyly, brain malformations, hydrocephalus, blindness, anosmia, obesity, and cognitive deficits . However, the functions of primary cilia in neurons and neuronal diseases, such as neurodegenerative diseases including AD, are largely unknown. Since 5-HT6 is mainly localized in the cilia, we hypothesized that 5-HT6 plays important roles in AD by regulating ciliary function.
Plasmids, chemicals, and antibodies
The mouse EGFP-tagged 5-HT6 plasmid was provided by Kirk Mykytyn (Ohio State University, USA). We constructed 13 mutagenesis plasmids. In each clone, the codons for amino acids 69, 70, 72, 106, 230, 234, 262, 265, 284, and 350 were replaced by triplet nucleotides. Mouse GFP-tagged ARL13B was amplified from cDNA from the mouse brain. 5-HT6-mcherry plasmid was constructed. PEGFP-N3 and mcherry-N3 were purchased from Addgene.
SB271046 drug (Santa Cruz Biotechnology), 5-HT6R siRNA (QIAGEN), siRNAs of 5-HT1A, 5-HT1B, 5-HT1D, 5-HT1F, 5-HT2A, 5-HT2B, 5-HT2C, 5-HT3A, 5-HT4, 5-HT5A, 5-HT5B, and 5-HT7 (Santa Cruz Biotechnology), 5-HT6R antibody (Abcam), MAP2 antibody (Abcam), AC3 antibody (Santa Cruz Biotechnology), SSTR3 antibody (Santa Cruz Biotechnology), ARL13B antibody (NeuroMab), ankyrin G antibody (Invitrogen), neurofascin 186 antibody (Abcam), γ-tubulin antibody (Sigma), IFT88 antibody (Proteintech), calbindin antibody (Swant), CamkII antibody (Millipore), Gad antibody (Abcam), Nav1.1 antibody (NeuroMab), Nav1.6 antibody (Alomone), Nav1.2 antibody (NeuroMab), Kv1.1 antibody (Alomone), and AlexaFluor fluorescence antibodies (Invitrogen) were utilized.
Cell culture and transfection
Primary culture neurons were cultured from the hippocampus of newborn C57 and APP/PS1 mice. The use of mice was in accordance with the regulations of the Peking University Animal Care and Use Committee (LSC-ZhangY-1). Fetal hippocampal samples were dissociated from the brain and digested by 0.25% trypsin (Invitrogen). Digestion was stopped by adding DMEM-F12 medium (Gibco) with 10% FBS (Gibco), after which the tissue was dispersed by a pipettor. After 2 minutes of precipitation, the supernatant was collected and centrifuged at 500 × g for 2 minutes. The cells were resuspended in DMEM-F12 medium with 10% FBS and plated on a coverslip coated with poly-D-lysine (Sigma) in 5% circulating CO2. Neurobasal medium (Gibco) containing Pen-Strep (Invitrogen), B27 (Gibco), and GlutaMAX (Thermo Fisher) was added to the medium after 4 hours. Half of the medium was replaced with fresh medium every 3 days.
Calcium phosphate transfection was used to transfer siRNA or a plasmid into hippocampal neurons. Neurons were transfected after 3–8 days. A mixed solution containing plasmids, CaCl2, H2O, and HBS (pH 7.1) was prepared and centrifuged. After adding the mixture to each well for 1–1.5 hours, HBS (pH 6.8) was added to each well to clear the precipitate. Finally, the original medium was added to every well and the cells were placed in a 37 °C incubator.
Cells were washed in PBS and fixed in 4% PFA (Sigma) at room temperature. Next, cells were permeabilized in 0.1% Triton at 4 °C, after which they were blocked in 5% donkey serum at room temperature and incubated with primary antibodies at 4 °C for 24 hours. Secondary antibodies were added to the cells for 1 hour in the dark. Nuclei were stained by DAPI (Sigma) for 15 minutes. Finally, coverslips with cells were put onto slides for imaging.
Each mouse was perfused with 0.9% saline and 4% PFA by myocardial perfusion. Brain slices were coated with OCT (Tissue-Tek) and sectioned by freezing microtomy into sections 35 μm thick. Next, the tissue samples were permeabilized in 0.3% Triton, blocked in 5% donkey serum, incubated in primary and secondary antibodies, and stained with DAPI. The brain slices were blocked in protection liquid.
HEK293T cells were transfected using Lipofectamine 2000 (Invitrogen) for 48 hours. Cells were collected using RIPA lysis buffer (R&D) and proteins were obtained by centrifugation. Proteins were incubated with protein A + G agarose beads (Beyotime) and antibodies overnight. Immunoprecipitated proteins were collected after washing three times with PBS (Sigma). The bicinchoninic acid assay (Pierce) was used to measure protein concentration. Proteins were denatured at 100 °C for 5 minutes and separated in 10% SDS-PAGE at 80 mA for 2 hours. They were then transferred to PVDF membrane (Millipore) at 100 mA for 2 hours. The membrane was blocked in 5% BSA (Sigma) in TBS (Sigma) with 0.1% Tween 20 (TBST) at room temperature for 1 hour. The antibodies were diluted at suitable concentration and added into TBST with 5% BSA. Primary antibodies were incubated overnight at 4 °C. HRP-conjugated second antibodies were added after washing three times for 10 minutes by TBST. Enhanced chemiluminescence was used to detect optical density of HRP after washing three times for 10 minutes. BioRad ChemiDox (BioRad) was used to analyze optical density of HRP.
Morris water maze test
Twelve-month-old male APPswe/PS1ΔE9 and C57 WT mice were subjected to the water maze test. Mice were kept at room temperature with food and water. Both APP/PS1 and WT mice were divided into two groups. Every group included 12 male APP/PS1 mice. The mice in the SB271046 group were given an intraperitoneal injection of SB271046 (10 mg/kg) 2 hours before water maze training and tested for 7 days. In the vehicle group, mice were injected with an equal amount of 0.9% NaCl for 7 continuous days.
The water maze test consisted of a circular pool (60 cm in diameter and 50 cm high). The pool was filled with 21–22 °C water at a depth of 30 cm. The hidden platform was 11 cm in diameter and placed 1.5 cm under the surface of the water. The pool was divided into four quadrants. The hidden platform was in one of the four quadrants. Mice were tested in the four quadrants with a fixed order. The starting location was random in the four quadrants. Mice spent four sessions of 90 s in the pool to adapt to the water environment 1 day before the experiment. The test began when a mouse was placed on the quadrants facing the wall and ended when the mouse got onto the platform. If a mouse got onto the platform within 90 s, it remained for 10 s on the platform. If the mouse failed to get onto the platform within 90 s, it was guided to the platform and stayed there for 10 s. The mice were trained and recorded in four quadrants per day for 6 continuous days to measure the swimming speed, time to the platform, and swimming distance with the hidden platform. On day 7, the hidden platform was removed. Then, we tested the time spent in the target quadrant without the platform and the crossing times at the original location of the platform on the seventh day. The position of the mouse was traced by the automatic tracking system (CSI) in real time.
Twelve-month-old male APPswe/PS1ΔE9 and C57 WT mice were used to test in the Y maze. APP/PS1 mice were divided into two groups (n = 12). One group was injected with SB271046 (10 mg/kg) and the other was injected with equal vehicle (0.9% NaCl). C57 WT mice were also divided into SB271046 and vehicle groups (n = 12). Drug or vehicle was injected 2 hours before the Y maze. The apparatus consisted of three arms which were 34 cm long, 8 cm wide, and 14.5 cm high. There are different images pasted on the end of three inner arms, containing a red circle, a black square, and a white triangle. The test began when a mouse was placed on the apparatus facing the wall. After 10 minutes, the experiment was ended. The automatic tracking system (CSI) traced the position of the mouse in real time and calculated the times when mice chose the correct arm. The correct arm was the one which was different form the last and next arms which mice chose.
Twelve-month-old male APPswe/PS1ΔE9 and C57 WT mice were used. Both APP/PS1 and C57 WT mice were divided into two groups (n = 12). One group was injected with SB271046 (10 mg/kg) and the other was injected with equal vehicle (0.9% NaCl). Drug or vehicle was injected 2 hours before training or testing for 2 continuous days. Mice were adapted on an operant box for 6 minutes on the day before the experiment. Training started for a period of 6 minutes on day 1. Mice received a 60-dB noise for a period of 30 s, followed by a 0.8-mA foot shock for a period of 5 s. This stimulus was repeated three times with an interval of 150 or 90 s. The mice were then placed in the chamber for 5 minutes on day 2. The device then recorded the duration of the freezing time of the mice.
For univariate data analysis, data were detected whether they were normally distributed or not. If data were equally distributed, a t test was used in two groups and one-way ANOVA followed by Dunnett’s test were used in more than two groups. If data were nonnormally distributed, a Mann–Whitney U test was used in two groups and the Kruskal–Wallis test followed by Dunnett’s test were used in more than two groups. For multivariable data analysis, two-way ANOVA followed by Bonferroni test or three-way ANOVA followed by Holm–Sidak test were used.
5-HT6 regulated primary cilia length
5-HT6 is mainly located at the primary cilia, where it plays an important role . Our results confirmed that 5-HT6 was localized to the primary cilia of hippocampal neurons, colocalized with two other cilia-specific proteins ARL13B and SSTR3 (Additional file 1: Figure S1A). Exogenous 5-HT6-enhanced green fluorescent protein (EGFP) was also located at the cilia, which were marked by ARL13B, γ-tubulin, and intraflagellar transport protein 88 (IFT88) (Additional file 1: Figure S1B). Therefore, we speculated that 5-HT6 may be involved in AD pathogenesis via regulating cilia function. Overexpression of 5-HT6 increased the ciliary length of WT C57 mouse hippocampal neurons (Fig. 1a, b). In addition, overexpression of 5-HT6 increased the number of cilia branches (Fig. 1a, b). The cilia became shorter after 5-HT6 was knocked down by its siRNA (Fig. 1c, d). Also, we tested the effect of other siRNA of 5-HT receptors on cilia length, including 5-HT1A, 5-HT1B, 5-HT1D, 5-HT1F, 5-HT2A, 5-HT2B, 5-HT2C, 5-HT3A, 5-HT4, 5-HT5A, 5-HT5B, 5-HT6, and 5-HT7. Both siRNAs of 5-HT6 and 5-HT3A decreased cilia length (Additional file 2: Figure S2A, B). 5-HT6 regulated cilia length in excitatory and inhibitory neurons. Overexpression of 5-HT6 increased cilia length, but siRNA of 5-HT6 decreased cilia length in excitatory neurons (Additional file 3: Figure S3A, B). Similar results were found in inhibitory neurons (Additional file 3: Figure S3C, D). Furthermore, ciliary length was measured as 5-HT6 expression was increased gradually in hippocampal neurons, confirming that enhanced 5-HT6 expression was associated with cilia elongation (Fig. 1e, f). Moreover, ciliary length increased as the concentration of 5-HT, a typical agonist of 5-HT6, was increased (Fig. 1g, h).
The primary cilia were elongated in the hippocampus of APP/PS1 mice
We first measured the length of cilia in cultured APP/PS1 and WT mouse neurons. APP/PS1 neurons had increased ciliary length. siRNA of 5-HT6 decreased cilia length in hippocampal neurons of APP/PS1 and WT mice (Fig. 1i, j). To measure the length of the primary cilia in vivo, brain slices of APP/PS1 and WT mice (6 months old) were prepared using cryosectioning and immunostaining. Adenylate cyclase 3 (AC3) and 4′,6-diamidino-2-phenylindole (DAPI) were used to label the primary cilia and indicate the nuclei, respectively. The primary cilia were elongated in the hippocampus of APP/PS1 mice (Fig. 2a–c). The number of cells with cilia, which was normalized by the number of nuclei, was similar in the APP/PS1 and WT mice (Fig. 2c). To verify elongation of the cilia of APP/PS1 mice in vitro, hippocampal neurons of newborn mice were cultured and stained for cilia markers. ADP-ribosylation factor-like protein 13B (ARL13B) was used to indicate the primary cilia . The data confirmed that the length of the primary cilia of hippocampal neurons from APP/PS1 mice increased significantly in comparison with that of WT mice (Fig. 2d, e). Thus, both in-vivo and in-vitro experiments showed that the cilia of the hippocampal neurons of APP/PS1 mice were elongated.
5-HT6 expression levels were detected by western blot analysis in different brain sections from 6-month-old APP/PS1 or WT mice (Fig. 2f). In the olfactory blub, striatum, diencephalon, and cerebellum of APP/PS1 mice, the expression level of 5-HT6 protein was not significantly different from that of WT mice. However, expression levels of 5-HT6 in APP/PS1 mice was remarkably upregulated in the hippocampus and in comparison with WT mice (Fig. 2g). Therefore, we suggested that cilia elongation in APP/PS1 mice might be induced by upregulation of 5-HT6 levels.
5-HT6 regulated cilia morphology and protein localization
To identify the key sites of 5-HT6 involved in its regulation of ciliary length, 13 reported mutations of 5-HT6 (F69L, T70I, D72A, D106A, A230F, Q234F, K262A, K265A, F284A, S350A, S350D, A230F + Q234F, and F69L + T70I + D72A) were examined. The related sequence information is presented in Additional file 4: Figure S4B. Mutations D72A, D106A, and F69L + T70I + D72A restored the ciliary length and the number of primary cilia branches when overexpressed in WT neurons (Fig. 3a–c and Additional file 4: Figure S4A). Mutations F69L, T70I, D72A, D106A, F284A, A230F + Q234F, and F69L + T70I + D72A restored the percentage of neurons with branched cilia (Fig. 3d). These results indicated that sites D72 and D106 and the combination of F69, T70, and D72 were primarily responsible for ciliary regulation by 5-HT6.
In addition, we also found that 5-HT6 mutations altered cilia protein localization. ARL13B, normally located at the primary cilia, disappeared from the primary cilia with 5-HT6 overexpression (Fig. 4a). Mutations D72A, D106A, and F69L + T70I + D72A restored the ciliary location of ARL13B (Fig. 4b). 5-HT6 affected ARL13B localization at the primary cilia. Therefore, we speculated that ARL13B was involved in 5-HT6 regulation of ciliary length. CO-IP experiments indicated that ARL13B and 5-HT6 interacted in HEK293T cells (Fig. 5a). Anti-ARL13B antibodies were able to pull down 5-HT6 protein, while anti-5-HT6 antibodies were able to pull down ARL13B protein. An ARL13B-GFP vector was constructed and transfected into hippocampal neurons with the 5-HT6 vector, after which the length of the primary cilia was measured (Fig. 5b). Cotransfection of ARL13B and 5-HT6 shortened the primary cilia, in comparison with transfection of 5-HT6 alone (Fig. 5c), indicating that ARL13B recovered cilia elongation induced by overexpression of 5-HT6. Total protein expression of ARL13B was not influenced by 5-HT6 overexpression or knock-down (Fig. 5d). 5-HT6 mutants D72A and F69L + T70I + D72A altered the activity of cyclic adenosine monophosphate (cAMP). When 5-HT concentrations in HEK293T cells were increased gradually, cells with the D72A and F69L + T70I + D72A mutations showed a reduced cAMP concentration in comparison with that of cells with WT 5-HT6 (Fig. 5e). These results showed that ARL13B and downstream cAMP activation were involved in 5-HT6 regulation of ciliary length.
5-HT6 altered axonal length and AIS morphology
Interestingly, ankyrin G (AnkG), a typical marker for the AIS , was translocated to the primary cilia of hippocampal neurons with overexpression of 5-HT6 (Fig. 4c). Mutating 5-HT6 did not restore the AIS localization of AnkG (Additional file 5: Figure S5A). To determine whether 5-HT6 affects total protein levels of AnkG, we measured their expression levels in HEK293T cells with 5-HT6 overexpression or knock-down. We found that 5-HT6 had no influence on the total protein levels of AnkG (Additional file 5: Figure S5B). Together, these results revealed that 5-HT6 influenced the localization of AnkG without affecting its total protein abundance.
The AIS is the region initiating action potential at the axon. We measured the AIS length with 5-HT6 overexpression. The AIS length was calculated from the 1/3 brightest point near the soma to the 1/3 brightest point far from the soma . Overexpression of 5-HT6 decreased the axonal length in both WT and APP/PS1 mouse neurons, while siRNA against 5-HT6 increased axonal length (Fig. 6a, b). The axonal length was short in APP/PS1 mice, compared with WT mice (Fig. 6b). In addition, 5-HT6 regulated the morphology of AIS in hippocampal neurons. AnkG and AIS marker neurofascin 186 (NF186) were colocalized at the AIS of hippocampal neurons (Fig. 6c). Since 5-HT6 altered AnkG subcellular localization, we used NF186 as the indictor for the AIS. Overexpression of 5-HT6 increased the AIS length, whereas siRNA against 5-HT6 decreased the AIS length (Fig. 6d, e). Mutations F69L, T70I, D72A, D106A, A230F, and K262A recovered the increased AIS length caused by 5-HT6 overexpression (Fig. 6f and Additional file 6: Figure S6). The distance between the AIS and cell body was influenced by 5-HT6 as well (Fig. 6g, h). Upregulation of 5-HT6 decreased the distance from the AIS to the cell body. Mutations D72A and K265A significantly recovered the decreased distance between the AIS and cell body (Fig. 6i and Additional file 7: Figure S7). In addition, 5-HT6 also affected the percentage of neurons with branched AIS (Fig. 6j, k). Overexpression of 5-HT6 increased the percentage of neurons with branched AIS. The F69L + T70I + D72A mutation of 5-HT6 reduced the percentage of neurons with branched AIS (Fig. 6l and Additional file 8: Figure S8). These results suggested that 5-HT6 influenced the morphology of the AIS in hippocampal neurons, including the AIS length, the distance from the cell body to the AIS, and the percentage of neurons with branched AIS.
In order to explore the function of changed AIS morphology in neurons, some ion channels specifically located at AIS were tested. Nav 1.2, Nav1.6, and Kv1.1 were located at the AIS. The position pattern of Nav1.2 at AIS was altered (Additional file 9: Figure S9A). Overexpression of 5-HT6 increased the AIS length and percentage of branched AIS marked by Nav1.2, and decreased distance from the AIS to the cell body (Additional file 9: Figure S9B). siRNA of 5-HT6 caused the opposite changes. There was no difference in the location of Nav1.6 with regulation of 5-HT6 (Additional file 9: Figure S9C, D). The position pattern of Kv1.1 was also unaltered (Additional file 9: Figure S9E, F).
The AIS of APP/PS1 mouse neurons showed changes similar to those caused by overexpression of 5-HT6 in WT hippocampal neurons (Fig. 7a). The AIS length of APP/PS1 neurons increased significantly in comparison with that of WT mice, while the distance between the AIS and cell body was shorter (Fig. 7b). Moreover, the percentage of hippocampal neurons with branched AIS was higher in APP/PS1 mice than in WT mice (Fig. 7b). siRNA to 5-HT6 successfully recovered the morphologies in APP/PS1 neurons almost to the WT level (Fig. 7c). Together, these results showed that the upregulated level of 5-HT6 was associated with axonal length and AIS morphology alterations in APP/PS1 mice.
5-HT6 antagonist SB271046 restored cognitive impairment in APP/PS1 mice
The typical behavior defect of APP/PS1 mice is cognitive impairment. The Morris water maze test was applied to assess cognition after APP/PS1 mice were administered SB271046, a 5-HT6 antagonist. The APP/PS1 mice received 10 mg/kg SB271046 every day for 7 days by intraperitoneal injection. The control group received 0.9% NaCl in the same manner and on the same schedule. The concentration of SB271046 in serum of WT mice increased gradually (Additional file 10: Figure S11C). APP/PS1 mice had typical cognitive impairment, embodied in the longer swimming distance to platform, longer time to platform, shorter time in target quadrant, and fewer times of crossing the original platform (Fig. 7d). The swimming speed of SB270146-treated mice was not different between APP/PS1 and WT mice (Fig. 7d). SB271046 decreased the swimming distance and time to platform, and increased time in the target quadrant and crossing times (Fig. 7d). After the water maze test, the cilia length and AIS morphology of mice were calculated. SB271046 decreased cilia length in APP/PS1 and WT mice (Fig. 7e, f). AIS morphology was recovered by SB271046 (Fig. 7g, h). SB 271046 decreased the AIS length and percentage of branched AIS, and increased the distance between the AIS and cell body (Fig. 7h). Also, siRNA of 5-HT6 increased arborization of dendrites in APP/PS1 and WT mice in cultured neurons. The arborization of dendrites of neurons in APP/PS1 mice was lower than that of WT mice in vitro (Additional file 11: Figure S10A-B). SB271046 also improved arborization of dendrites of APP/PS1 and WT mice in vivo (Additional file 11: Figure S10C, D).
In addition, the Y maze and fear conditioning were used to tested the function of 5-HT6 in cognition. SB271046 increased the percentage of correct times in APP/PS1 and WT mice (Additional file 10: Figure S11A). The percentage of correct times which APP/PS1 mice chose was lower than that of WT mice. Fear conditioning was a typical test to assess the memory of mice. SB271046 enhanced duration of freezing in APP/PS1 and WT mice (Additional file 10: Figure S11B). The duration of freezing time in APP/PS1 mice was lower than that of WT mice. Collectively, results indicated that 5-HT6 antagonist SB271046 recovered the cognitive impairment of APP/PS1 mice.
Previous studies showed that 5-HT6 is located at the primary cilia, where it plays an important role . Therefore, we suspected that 5-HT6 may influence cognition in AD patients by regulating the primary cilia function. The results of our study showed that overexpression of 5-HT6 increased the length of primary cilia. Meanwhile, 5-HT6 was upregulated in the hippocampus of APP/PS1 mice. 5-HT6, a G protein coupled receptor, not only regulated the cilia length and branch but also had an effect on the location of ARL13B and AnkG in primary cilia. 5-HT3A, which belongs to the Cys-loop superfamily of ligand-gated ion channels, also affected the cilia length. Both 5-HT6 and 5-HT3A are 5-HT receptors. This indicated that 5-HT receptors might cause the imbalance of some proteins and compromise ciliary signaling. The disorder of proteins in cilia could impair the ability of cilia to detect and respond to signaling in the local extracellular environment. As a result, ciliary function was influenced in APP/PS1 mice.
In addition, 5-HT6 regulated AIS morphology and axonal length. Compared to WT mice, the axonal length was shorter in APP/PS1 mice. The synapse, the function of which is transmitting a signal to the next neuron, is at the end of the axon. Therefore, the shorter axonal length in APP/PS1 mice might affect synaptic signal transmission in a distant position. In hippocampal neurons of APP/PS1 mice, the AIS was elongated, showed a greater degree of branching, and was closer to the cell body. There are many ion channels anchored at the AIS of neurons. The longer AIS length in APP/PS1 mice might influence the location of these ion channels, and then the initiation of action potential might be affected. Besides, the longer AIS might also affect protein transporting into the axon, which would lead to the variation of axonal length. The decreased axonal length and increased AIS length were comprehensive and complex results in APP/PS1 mice and the concrete mechanism should be explored in future experiments.
Interestingly, our finding that ciliary length was increased in APP/PS1 mice differs from that of a previous study on 3xAD transgenic mice, in which cilia length was reduced in hippocampal dentate granule cells . We speculated that three factors might underlie this difference. Different animal models of AD were used in each study. 3xAD-transgenic mice harbor APPSwe, PS1M148V, and tauP301L, while the transgenic mice used in this study harbor APPswe and PS1ΔE9. The use of different models of AD may lead to differing results. In addition, different types of cell were assessed in each study. The study using 3xAD-transgenic mice only measured the length of primary cilia in dentate granule cells, while cilia length was measured in dentate granule cells and pyramidal neurons in CA1 and CA3 in our study. Finally, female 3xAD-transgenic mice were used in their study, and male APP/PS1 mice were tested in our experiments. Different gender of the mice might also affect results. The reason for the differences of cilia length among animal models of AD merits further study.
Typical symptoms of patients with AD include severe deficits in learning and memory. Several studies have shown that 5-HT6 antagonists can improve learning and memory ability . In this study, we aimed to understand the mechanism through which 5-HT6 antagonists improve cognition in patients with AD. Numerous convincing studies have shown that there is a strong correlation between the primary cilia and cognition, indicating that 5-HT6 may affect cognition through the primary cilia. Children with primary ciliary dyskinesia, a typical ciliopathy due to structure and functional abnormalities of respiratory cilia, showed significant behavioral and social competence problems . The cognitive functioning and memory performances of patients with Bardet-Biedl Syndrome (BBS), which is also a kind of ciliopathy, was due to hippocampal dysgenesis . Joubert syndrome (JS) in patients with intellectual disability is also related to primary cilia . Primary cilia deficiency can lead to ciliopathies, including Joubert syndrome (JBTS), polycystic kidney disease, developmental disorders, cognitive impairment, and cognitive damage . Interfering with the normal function of BBS complex delays development and leads to mental retardation and learning impairment . SB399885, another antagonist of 5-HT6, could shorten primary cilia of striatal neurons after 24 hours in culture . It also reversed the age-dependent deficit of rats in the water maze test . In mice, interfering with genes required for proper cilia function delays spatial learning and recognition . Mutation of SSTR3, which is located at the primary cilia, impairs hippocampal LTP in mice [35, 36]. AC3, which was located at primary cilia, was required for novel object learning and extinction of contextual memory . Our water maze test, Y maze, and fear condition results showed that 5-HT6 antagonist SB271046 recovered cognitive impairment in APP/PS1 mice, indicating that 5-HT6 has a close relationship with the cognition of APP/PS1 mice. Meanwhile, our results showed that primary cilia length was recovered by SB271046 in behavioral experiments, which strongly suggested that 5-HT6 might affect the cognition through regulating the primary cilia. The more specific inner mechanism deserves future study.
Our data showed that 5-HT6 regulated ciliary length through ARL13B, influenced the function of neuronal primary cilia, affected the AIS morphology and axonal length in APP/PS1 mice. Besides, SB271046, 5-HT6 antagonist, could recover the cognitive impairment of APP/PS1 mice. Based on our results and those of reported studies, we suggest that 5-HT6 plays a critical role in AD development through regulating the morphology and function of neuronal primary cilia, which is possibly related to the AIS and axon alterations in AD development.
Serotonin 6 receptor
Adenylyl cyclase 3
Axon initial segment
Amyloid precursor protein/presenilin 1
ADP-ribosylation factor-like protein 13B
Ca2+/calmodulin-dependent protein kinase II
Central nervous system
Enhanced green fluorescent protein
Familial Alzheimer’s disease
Glutamic acid decarboxylase
Intraflagellar transport protein 88
Potassium voltage-gated channel subfamily A member 1
Microtubule-associated protein 2
Sodium channel, voltage-gated, type I, alpha subunit
Sodium channel, voltage-gated, type II, alpha subunit
Sodium channel, voltage gated, type VIII, alpha subunit
Somatostatin receptor 3
Price DL, Sisodia SS, Gandy SE. Amyloid beta amyloidosis in Alzheimer's disease. Curr Opin Neurol. 1995;8(4):268–74.
Tsai SJ, Liu HC, Liu TY, Wang YC, Hong CJ. Association analysis of the 5-HT6 receptor polymorphism C267T in Alzheimer's disease. Neurosci Lett. 1999;276(2):138–9.
Liu HC, Hong CJ, Liu CY, Lin KN, Tsai SJ, Liu TY, Chi CW, Wang PN. Association analysis of the 5-HT6 receptor polymorphism C267T with depression in patients with Alzheimer's disease. Psychiatry Clin Neurosci. 2001;55(4):427–9.
Upton N, Chuang TT, Hunter AJ, Virley DJ. 5-HT6 receptor antagonists as novel cognitive enhancing agents for Alzheimer's disease. Neurotherapeutics. 2008;5(3):458–69.
Ramirez MJ. 5-HT6 receptors and Alzheimer's disease. Alzheimers Res Ther. 2013;5(2):15.
Wicke K, Haupt A, Bespalov A. Investigational drugs targeting 5-HT6 receptors for the treatment of Alzheimer's disease. Expert Opin Investig Drugs. 2015;24(12):1515–28.
Wilkinson D, Windfeld K, Colding-Jorgensen E. Safety and efficacy of idalopirdine, a 5-HT6 receptor antagonist, in patients with moderate Alzheimer's disease (LADDER): a randomised, double-blind, placebo-controlled phase 2 trial. Lancet Neurol. 2014;13(11):1092–9.
Guadiana SM, Semple-Rowland S, Daroszewski D, Madorsky I, Breunig JJ, Mykytyn K, Sarkisian MR. Arborization of dendrites by developing neocortical neurons is dependent on primary cilia and type 3 adenylyl cyclase. J Neurosci. 2013;33(6):2626–38.
Jacobshagen M, Niquille M, Chaumont-Dubel S, Marin P, Dayer A. The serotonin 6 receptor controls neuronal migration during corticogenesis via a ligand-independent Cdk5-dependent mechanism. Development. 2014;141(17):3370–7.
Bayer TA, Wirths O, Majtenyi K, Hartmann T, Multhaup G, Beyreuther K, Czech C. Key factors in Alzheimer's disease: beta-amyloid precursor protein processing, metabolism and intraneuronal transport. Brain Pathol. 2001;11(1):1–11.
Wirths O, Weis J, Szczygielski J, Multhaup G, Bayer TA. Axonopathy in an APP/PS1 transgenic mouse model of Alzheimer's disease. Acta Neuropathol. 2006;111(4):312–9.
Tadenev AL, Kulaga HM, May-Simera HL, Kelley MW, Katsanis N, Reed RR. Loss of Bardet-Biedl syndrome protein-8 (BBS8) perturbs olfactory function, protein localization, and axon targeting. Proc Natl Acad Sci U S A. 2011;108(25):10320–5.
Yao PJ, Petralia RS, Ott C, Wang YX, Lippincott-Schwartz J, Mattson MP. Dendrosomatic sonic hedgehog signaling in hippocampal neurons regulates axon elongation. J Neurosci. 2015;35(49):16126–41.
Monsma Jr FJ, Shen Y, Ward RP, Hamblin MW, Sibley DR. Cloning and expression of a novel serotonin receptor with high affinity for tricyclic psychotropic drugs. Mol Pharmacol. 1993;43(3):320–7.
Ward RP, Hamblin MW, Lachowicz JE, Hoffman BJ, Sibley DR, Dorsa DM. Localization of serotonin subtype 6 receptor messenger RNA in the rat brain by in situ hybridization histochemistry. Neuroscience. 1995;64(4):1105–11.
Gerard C, Martres MP, Lefevre K, Miquel MC, Verge D, Lanfumey L, Doucet E, Hamon M, el Mestikawy S. Immuno-localization of serotonin 5-HT6 receptor-like material in the rat central nervous system. Brain Res. 1997;746(1–2):207–19.
Mitchell ES, Neumaier JF. 5-HT6 receptors: a novel target for cognitive enhancement. Pharmacol Ther. 2005;108(3):320–33.
Louvi A, Grove EA. Cilia in the CNS: the quiet organelle claims center stage. Neuron. 2011;69(6):1046–60.
Chen MH, Wilson CW, Li YJ, Law KK, Lu CS, Gacayan R, Zhang X, Hui CC, Chuang PT. Cilium-independent regulation of Gli protein function by Sufu in Hedgehog signaling is evolutionarily conserved. Genes Dev. 2009;23(16):1910–28.
Corbit KC, Aanstad P, Singla V, Norman AR, Stainier DY, Reiter JF. Vertebrate Smoothened functions at the primary cilium. Nature. 2005;437(7061):1018–21.
Fry AM, Leaper MJ, Bayliss R. The primary cilium: guardian of organ development and homeostasis. Organogenesis. 2014;10(1):62–8.
Brailov I, Bancila M, Brisorgueil MJ, Miquel MC, Hamon M, Verge D. Localization of 5-HT(6) receptors at the plasma membrane of neuronal cilia in the rat brain. Brain Res. 2000;872(1-2):271–5.
Higginbotham H, Eom TY, Mariani LE, Bachleda A, Hirt J, Gukassyan V, Cusack CL, Lai C, Caspary T, Anton ES. Arl13b in primary cilia regulates the migration and placement of interneurons in the developing cerebral cortex. Dev Cell. 2012;23(5):925–38.
Hedstrom KL, Ogawa Y, Rasband MN. AnkyrinG is required for maintenance of the axon initial segment and neuronal polarity. J Cell Biol. 2008;183(4):635–40.
Grubb MS, Burrone J. Activity-dependent relocation of the axon initial segment fine-tunes neuronal excitability. Nature. 2010;465(7301):1070–4.
Chakravarthy B, Gaudet C, Menard M, Brown L, Atkinson T, Laferla FM, Ito S, Armato U, Dal Pra I, Whitfield J. Reduction of the immunostainable length of the hippocampal dentate granule cells' primary cilia in 3xAD-transgenic mice producing human Abeta(1-42) and tau. Biochem Biophys Res Commun. 2012;427(1):218–22.
Routledge C, Bromidge SM, Moss SF, Price GW, Hirst W, Newman H, Riley G, Gager T, Stean T, Upton N, et al. Characterization of SB-271046: a potent, selective and orally active 5-HT(6) receptor antagonist. Br J Pharmacol. 2000;130(7):1606–12.
Carotenuto M, Esposito M, Di Pasquale F, De Stefano S, Santamaria F. Psychological, cognitive and maternal stress assessment in children with primary ciliary dyskinesia. World J Pediatr. 2013;9(4):312–7.
Bennouna-Greene V, Kremer S, Stoetzel C, Christmann D, Schuster C, Durand M, Verloes A, Sigaudy S, Holder-Espinasse M, Godet J, et al. Hippocampal dysgenesis and variable neuropsychiatric phenotypes in patients with Bardet-Biedl syndrome underline complex CNS impact of primary cilia. Clin Genet. 2011;80(6):523–31.
Doherty D. Joubert syndrome: insights into brain development, cilium biology, and complex disease. Semin Pediatr Neurol. 2009;16(3):143–54.
Beales PL, Elcioglu N, Woolf AS, Parker D, Flinter FA. New criteria for improved diagnosis of Bardet-Biedl syndrome: results of a population survey. J Med Genet. 1999;36(6):437–46.
Brodsky M, Lesiak AJ, Croicu A, Cohenca N, Sullivan JM, Neumaier JF. 5-HT6 receptor blockade regulates primary cilia morphology in striatal neurons. Brain Res. 2017;1660:10–9.
Hirst WD, Stean TO, Rogers DC, Sunter D, Pugh P, Moss SF, Bromidge SM, Riley G, Smith DR, Bartlett S, et al. SB-399885 is a potent, selective 5-HT6 receptor antagonist with cognitive enhancing properties in aged rat water maze and novel object recognition models. Eur J Pharmacol. 2006;553(1–3):109–19.
Amador-Arjona A, Elliott J, Miller A, Ginbey A, Pazour GJ, Enikolopov G, Roberts AJ, Terskikh AV. Primary cilia regulate proliferation of amplifying progenitors in adult hippocampus: implications for learning and memory. J Neurosci. 2011;31(27):9933–44.
Stanic D, Malmgren H, He H, Scott L, Aperia A, Hokfelt T. Developmental changes in frequency of the ciliary somatostatin receptor 3 protein. Brain Res. 2009;1249:101–12.
Einstein EB, Patterson CA, Hon BJ, Regan KA, Reddi J, Melnikoff DE, Mateer MJ, Schulz S, Johnson BN, Tallent MK. Somatostatin signaling in neuronal cilia is critical for object recognition memory. J Neurosci. 2010;30(12):4306–14.
Wang Z, Phan T, Storm DR. The type 3 adenylyl cyclase is required for novel object learning and extinction of contextual memory: role of cAMP signaling in primary cilia. J Neurosci. 2011;31(15):5557–61.
The authors thank Hongxia Lv, Xiaochen Li, Hui Li, and Yu Wu (College of Life Sciences, Peking University) for their help in fluorescence imaging.
This work was supported by the National Science Foundation of China (NSFC) Fund for Distinguished Young Scholars (81425009) and NSFC major research grants (31630028 and 91632305).
Availability of data and materials
The use of mice was in accordance with the regulations of the Peking University Animal Care and Use Committee (LSC-ZhangY-1).
Consent for publication
All authors approved the manuscript for submission.
All authors declare no actual or potential conflicts of interest, including any financial, personal or other relationships, with other people or organizations within 3 years of beginning the work submitted that could inappropriately influence (bias) their work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Showing 5-HT6 located at the primary cilia of hippocampal neurons. (A) Endogenous 5-HT6 located at the primary cilia of neurons. ARL13B and SSTR3 were markers of primary cilia. MAP2 was the neuronal marker. γ-tubulin was the marker of the centrosome from which primary cilia grow. (B) Exogenous 5-HT6, which was expressed from 5-HT6-GFP vector, located at the primary cilia of neurons. ARL13B and IFT88 were markers of primary cilia. Scale bars, 10 μm. (DOC 8716 kb)
Showing siRNAs of 5-HT6 and 5-HT3A decreased cilia length. (A) Both 5-HT6 and 5-HT3A affected cilia length in hippocampal neurons. AC3 indicated cilia. (B) siRNA of 5-HT6 and 5-HT3A reduced cilia length. 5-HT3A-siRNA, **p < 0.01; 5-HT6-siRNA, *p < 0.05 (Kruskal–Wallis test, p < 0.001). Scale bars, 10 μm. All data presented as mean ± SEM. (DOC 4600 kb)
Showing the effect of 5-HT6 on cilia length in excitatory and inhibitory neurons. (A) 5-HT6 regulated cilia length in excitatory neurons. 5-HT6 indicated cilia. CamkII was the marker of excitatory neuron. (B) Overexpression of 5-HT6 increased cilia length. siRNA of 5-HT6 reduced cilia length. 5-HT6, ***p < 0.001; 5-HT6 siRNA, *p < 0.05 (Kruskal–Wallis test, p < 0.001). (C) 5-HT6 regulated cilia length in inhibitory neurons. Gad was the maker of inhibitory neuron. 5-HT6 indicated cilia. (D) Cilia length increased by overexpression of 5-HT6 and decreased by siRNA of 5-HT6. 5-HT6, ***p < 0.001; 5-HT6 siRNA, *p < 0.05 (Kruskal–Wallis test, p < 0.001). Scale bars, 10 μm. All data presented as mean ± SEM. (DOC 4724 kb)
Showing 5-HT6 mutations influenced the number of ciliary branches and the percentage of neurons with branched cilia. (A) Mutations D72A, D106A, and F69L + T70I + D72A recovered the number of branches in cilia. Mutations F69L, D72A, D106A, F284A, A230F + Q234F, and F69L + T70I + D72A restored the percentage of neurons with branched cilia. 5-HT6-EGFP indicated primary cilia. (B) Sequence information of all overexpression, siRNA, and specific point mutations. Underline indicates sequence after the mutation; () indicates original sequence. Scale bars, 10 μm. (DOC 6514 kb)
Showing 5-HT6 influenced the location of AnkG between AIS and cilia. (A) AnkG (purple) located at the primary cilia following overexpression of 5-HT6 in hippocampal neurons. Mutations of 5-HT6 did not recover the location of AnkG in primary cilia. 5-HT6-EGFP (green) indicated primary cilia. (B) 5-HT6 had no effect on total protein level of AnkG in HEK293T cells. Kruskal–Wallis test, p = 0.77. β-actin was the internal control. Scale bars, 10 μm. All data presented as mean ± SEM. (DOC 6089 kb)
Showing 5-HT6 mutations influenced AIS length in hippocampal neurons. Mutations F69L, T70I, D72A, D106A, A230F, and K262A decreased AIS length, compared with that of overexpression of 5-HT6. 5-HT6-EGFP (green) indicated primary cilia. NF186 (purple) was the marker of the AIS. Scale bars, 10 μm. (DOC 6200 kb)
Showing 5-HT6 mutations influenced the distance between the AIS and cell body in hippocampal neurons. Mutations D72A and K265A restored the distance between the AIS and cell body. 5-HT6-EGFP (green) indicated primary cilia. NF186 (purple) was a marker of the AIS. Scale bars, 10 μm. (DOC 6784 kb)
Showing 5-HT6 mutations influenced the percentage of cells with branched AIS in hippocampal neurons. Mutations F69L + T70I + D72A decreased the percentage of neurons with branched AIS, compared with that of 5-HT6 overexpression. 5-HT6-EGFP (green) indicated primary cilia. NF186 (purple) was the marker of the AIS. Scale bars, 10 μm. (DOC 7506 kb)
Showing 5-HT6 regulated AIS morphology marker by Nav1.2. (A) 5-HT6 affected the Nav1.2 positioning pattern in the AIS. Nav1.2 was located at the AIS of neurons specifically. (B) Overexpression of 5-HT6 increased AIS length and percentage of branched AIS marked by Nav1.2, and decreased distance from AIS to cell body. siRNA of 5-HT6 had the opposite results. AIS length: 5-HT6, *p < 0.05; 5-HT6 siRNA, *p < 0.05 (Kruskal–Wallis test, p < 0.01). Distance to cell body: 5-HT6, *p < 0.05; 5-HT6 siRNA, *p < 0.05 (Kruskal–Wallis test, p < 0.001). Branched AIS: 5-HT6, *p < 0.05; 5-HT6 siRNA, *p < 0.05 (Kruskal–Wallis test, p < 0.01). (C) 5-HT6 had no effect on Nav1.6 positioning pattern in AIS. Nav1.6 was located at AIS of neuron specifically. (D) There was no difference in positioning pattern of Nav1.6 with regulation of 5-HT6. AIS length (Kruskal–Wallis test, p = 0.94), distance to cell body (Kruskal–Wallis test, p = 0.59), branched AIS (Kruskal–Wallis test, p = 0.44). (E) 5-HT6 had no effect on Kv1.1 positioning pattern in AIS. (F) Positioning pattern of kv1.1 unaltered. AIS length (Kruskal–Wallis test, p = 0.06), distance to cell body (Kruskal–Wallis test, p = 0.68), branched AIS (Kruskal–Wallis test, p = 0.38). Scale bars, 10 μm. All data presented as mean ± SEM. (DOC 6017 kb)
Showing SB271046 restored cognitive impairment of APP/PS1 mice. (A) SB271046 improved the percentage of correct times in APP/PS1 and WT mice in Y maze. APP/PS1 mice chose more wrong arms, compared with WT mice. Vehicle versus SB271046 in WT, *p < 0.05; vehicle versus SB271046 in APP/PS1, *p < 0.05; WT versus APP/PS1 in vehicle, *p < 0.05 (two-way ANOVA: interaction, F 1,34 = 0.02, p = 0.88; treatment, F 1,34 = 11.88, p < 0.01; gene, F 1,34 = 7.84, p < 0.01). (B) SB271046 increased the duration of freezing in APP/PS1 and WT mice in fear conditioning. The duration of freezing of APP/PS1 mice was longer than that of WT mice. Vehicle versus SB271046 in WT, *p < 0.05; vehicle versus SB271046 in APP/PS1, *p < 0.05; WT versus APP/PS1 in vehicle, *p < 0.05 (two-way ANOVA: interaction, F 1,33 = 1.41, p = 0.24; treatment, F 1,33 = 11.51, p < 0.01; gene, F 1,33 = 15.35, p < 0.001). (C) Concentration of SB271046 in serum in water maze test. HPLC-MS was used to measure concentration of SB271046 in serum after 2 hours of injecting drug. Concentration of SB271046 was increased gradually from 2 μM on day 1 to 12 μM ion day 7. Kruskal–Wallis test, p < 0.01. (DOC 190 kb)
Showing 5-HT6 affected the arborization of dendrites. (A) 5-HT6 influenced arborization of dendrites of excitatory neurons in vitro. CamKII was the maker of excitatory neurons. (B) siRNA of 5-HT6 increased the arborization of dendrites in cultured excitatory neuron APP/PS1 and WT mice. The arborization of dendrites was lower in neurons of APP/PS1 mice, compared with WT mice (three-way ANOVA: distance, F 18,4465 = 445.26, p < 0.001; treatment, F 1,4465 = 5.3, p < 0.05; gene, F 1,4465 = 5.69, p < 0.05; distance × treatment, F 18,4465 = 4.85, p < 0.001; distance × gene, F 18,4465 = 24.85, p < 0.001; treatment × gene, F 1,4465 = 102.84, p < 0.001; distance × treatment × gene, F 18,4465 = 6.49, p < 0.001). (C) 5-HT6 influenced arborization of dendrites of excitatory neuron in vivo. (D) SB271046 enhanced arborization of dendrites of APP/PS1 and WT mice. Arborization of dendrites was lower in neurons of APP/PS1 mice, compared with WT mice (three-way ANOVA: distance, F 14,4005 = 59.3, p < 0.001; treatment, F 1,4005 = 11.3, p < 0.001; gene, F 1,4005 = 18.91, p < 0.001; distance × treatment, F 14,4005 = 0.98, p = 0.48; distance × gene, F 14,4005 = 1.97, p < 0.05; treatment × gene, F 1,4005 = 13.74, p < 0.001; distance × treatment × gene, F 14,4005 = 0.88, p = 0.59). Scale bars, 10 μm. All data presented as mean ± SEM. (DOC 6890 kb)