2. Generation of Isogenic FXS Full-Mutation hESC Subclones
In previous studies, several FXS-hESC lines representing variability in FX populations were compared to several control lines to explore the role of
FMR1 in FX
[23][24]. Other studies use CRISPR/Cas9 for removal of the CGG repeats or the hypermethylation of the CGG expansion mutation
[19][25], but such genetic intervention is commonly associated with off-target effects. To overcome this limitation, in the current study, isogenic subclones were isolated, freed from any genetic intervention, from the Lis_FX6 hESC line that was previously derived in our lab, following PGD and present CGG repeats from the normal-through the full-mutation range. Sixty-one single-cell colonies were manually isolated and screened for their CGG repeat number in the 5’-untranslated region of
FMR1 (
Figure 1A) and by the AmplideX PCR/CE
FMR1 Reagents (Asuragen, TX, USA) commonly used for clinical diagnosis of FXS (
Figure 1B). Our results show that researchers successfully isolated an FXS subclone that is in the full mutation range (>200 CGG repeats; clone 8A), a subclone that is in the pre-mutation range (50–200 CGG repeats; clone 7L) and another subclone that is free of the mutation (<50 CGG repeats; clone 7B) that will serve in all further experiments as an isogenic control (WT;
Figure 1C). These subclones were analyzed for CGG repeat length three times during the study period and the results demonstrate their stability over extended cell divisions. Importantly, all lines and subclones were shown to be pluripotent with a normal karyotype (see
Scheme 1) and maintained their repeats number along the differentiation process. Induced differentiation into neurons by overexpression of
NEUROG-1 already generated neurons with typical morphology within 7 days (
Figure 1E). These iNs already completely inactivated the
FMR1 gene, only in the FX full-mutation subclone, at 7 days following neuronal induction, at both the RNA and protein level (
Figure 1E,F), mimicking
FMR1 silencing in FXS developing fetuses. Researchers then characterized the expression profile of the FX-iNs and their isogenic counterparts. Our results show that by day 39, both FX- and control-iNs completely inactivated the pluripotent gene
OCT4, and positively expressed the neural genes
PAX6,
Tuj1 and
MAP2. They also expressed
BRN2 and
FOXG1, which characterize excitatory cortical neurons. The pre- and post-synaptic genes
SYN1 and
PSD-95 were also expressed in both FX and WT iNs at day 39, demonstrating their maturity. Expression analysis of the glutamate receptors showed that two of four subunits of the AMPA receptor (
GRIA1, GRIA3) were found to be expressed in both FX and WT iNs at the same level. Interestingly,
GRIA4 expression was significantly higher in FXS neurons compared to WT (
p < 0.05). Both FX and WT iNs also express the ionotropic glutamate receptors kainate 1 (
GRIK1), the NMDA receptor subunit 1 (
GRIN1), the metabotropic glutamate receptor 5 (mGluR5), and the vesicular glutamate transporter-2 (vGlut2). The GABA receptor subunit α-2 (
GABRA2) was also expressed in these neurons, although with lower expression levels compared with glutamatergic genes (
Figure 1D). Collectively, the expression of this series of neuronal markers shows that FX and WT iNs are excitatory glutamatergic neurons.
Figure 1. Generation of Fragile X (FX) human embryonic stem cell (hESC) isogenic clones: (A) Analysis of the FX-hESC isogenic clones’ CGG repeat size by PCR-based repeat-number assay of hESC. H9 (WiCell, University of Wisconsin) served as control; (B) AmplideX PCR/CE FMR1 Reagents (Asuragen) commonly used for clinical diagnosis of FXS; (C) FX mutation status of isogenic iN clones; (D) neuronal characterization of iNs. Expression of neuronal genes in FX- and control-iNs. Quantitative RT-PCR analysis was performed on RNA extracted at 39 days in vitro (DIV). The housekeeping gene GUSB served as an internal control. Three independent experiments were performed, and data are represented as mean ± SEM. T-test: * p < 0.05 (E) FMRP expression in FX and isogenic control hESC 7 & 32 DIV. β-III tubulin (Tuj1; green) and FMRP (red). Scale bar: 200 µm; (F) RNA expression of FMR1 analyzed by qRT-PCR at 39 DIV (left). Data are represented as mean ± SEM. Three independent experiments were performed. T-test: *** p < 0.001. FMRP expression using western blot analysis (right). Protein was extracted at either day 7 or at 32 DIV. Human FMRP antibody is shown, and actin served as a positive control.
Scheme 1. Characterization of hESC subclones generated from Lis_FX6: (A) Polymorphic marker analyses of CA repeats (marker), known to identify each unique hESC line, are listed. Lis_FX6 original hESC line; #—subclone derived from the original hESC line; (B) expression of the OCT4/ SSEA4/ TRA-160 pluripotent genes (green) in the two isogenic clones. Nuclear staining: DAPI (blue). Scale bar: 100 µm; (C) karyotype by SKY analysis of clone 7B at passage 48 and clone 8A at passage 46.
3. The Effect of FMRP Expression on Neuronal Function and Network Development
The formation and maintenance of synaptic networks are necessary for neuronal functions through continuous changes in the neurons’ fine structure. They therefore analyzed the morphology of the FXS-iNs at four time points following neuronal induction (7, 12, 23 and 32 DIV). The results show that while iNs mature in vitro, the neurites elongate gradually from day 7 to 32 (Figure 2). Temporarily increased neurite growth was detected in FX-iNs at day 12 DIV (p < 0.0005), compared with WT, but their length was not different from WT later during development, indicating that their delayed elongation catches up later (Figure 2B). In addition, neurite branching also increased as the iNs matured, in both FX and WT neurons. However, at day 32 of differentiation, the FX-iNs had significantly fewer branch points than the WT (15.5 and 17.5 branch points/cell body cluster, respectively; p < 0.005), probably reflecting their retarded maturity (Figure 2C). In addition, the soma perimeter of both FX and WT iNs increased during development, with FX-iNs presenting a significantly smaller soma at all four time points (Figure 2D). Collectively these results indicate that the lack of FMRP alters the maturation of neurons, resulting in a smaller neuronal soma and less-branched neurites.
Figure 2. High-throughput image analysis of Fragile X (FX) induced neurons (iNs): (A) Representative images of WT and FX-iNs following Incucyte live-cell analysis, at four time points following neuronal induction (7, 12, 23, 32 DIV). Soma and neurites are presented in different colors. Scale bar: 200 µm; (B) total neurite length in mm/cell body cluster; (C) neurite branch points/cell body cluster. Two–three independent experiments were performed including 2500–4000 neurons at each time point. Data are represented as mean ± SEM. One-way ANOVA: ** p < 0.005, *** p < 0.0005; (D) representative images of iNs’ soma perimeter labeled with Tuj1 (green), at 32 DIV. Scale bar: 50 µm (left). Quantification of soma perimeter (µm) was analyzed by the cellSens Dimension software (right). Two (23 & 32 DIV) or three (7 & 12 DIV) independent experiments were performed on each day for both neurite length and neurite branch points including 2000–4000 neurons at each time point. Data are represented as mean ± SEM. One-way ANOVA: * p < 0.01, **** p < 0.0001.
Next, researchers explored the functional correlates of these morphological features at the single-cell level. Basic functional properties of iNs were analyzed using the patch-clamp technique at >40 DIV. As indicated in Figure 3, the basic properties of FX-iNs were not different from WT; both can discharge repetitive action potentials (APs) each followed by a large afterhyperpolarization (AHP; Figure 3A). In addition, the resting membrane potential and the APs’ amplitude were not different between FX and WT neurons. AHPs are slightly, but not significantly, larger in the FX-iNs compared to WT neurons. Interestingly, the firing threshold of WT neurons was significantly lower than that of the FX-iNs (−27 mV compared to −24 mV, p < 0.01; Figure 3A). FX-iNs, similar to WT, responded to both glutamate and GABA applied by pressure from an adjacent micropipette (Figure 3B). More importantly, FX neuronal networks were able to generate slightly more spontaneous excitatory postsynaptic currents (EPSC) than WT neurons (see Figure 3C). The mean amplitude of the spontaneous synaptic currents was significantly larger in the FX-iNs (−36.4 ± 6.16 pA in controls compared to −61.2 ± 10.63 pA, p < 0.02), which can indicate a higher number of synapses in the FX neurons. The frequency of EPSC was slightly, but not significantly, higher in the FX neurons compared to the controls. The averaged decay time of individual EPSC was similar in both WT and FX (Figure 3C). These results indicate that even though the somatic perimeter is larger in the control than in the FX (Figure 2D), the synaptic currents, likely to reside on dendrites, are not different between the two groups. Tentatively, there was no apparent difference in capacitance between neurons in the two groups (not shown). Consequently, the APs’ initiation difference (Figure 3A) probably compensates for the higher excitability associated with an increase in synaptic activity in the FX neurons.
Figure 3. Patch-clamp recordings of Fragile X (FX)-compared to control-induced neurons (iNs): (A) Top: Intrinsic properties of action potentials (APs) of FX-iNs compared to WT. Sample illustration of an IV curve, depicting changes in membrane voltage in response to hyperpolarizing and depolarizing current pulses. The right is an expanded AP taken from the battery on the left. Bottom: summary of recorded neuron properties: resting membrane potential (MP); firing threshold (** p < 0.01); AP amplitude; after-hyperpolarization (AHP). WT: n = 30, and FX-iNs: n = 20; (B) iNs Glutamate and GABA responses. Responses to pressure application of glutamate onto FX-iNs at resting membrane potential (−60 mV) in voltage-clamp (left). Bottom: comparison of FX vs. WT-iNs’ response to glutamate at >40 DIV (WT, n = 15; FX-iNs, n = 10). Right: responses to GABA at resting potential (−60 mV) and at a depolarizing potential (−30 mV). Bottom: Comparison of FX vs. WT-iNs’ response to GABA at >40 DIV (WT, n = 7; FX-iNs, n = 8); (C) top: Spontaneous excitatory synaptic currents in FX-iNs compared to WT. Sample illustration of a trace recorded from FX neurons shown at low (left) and high (right) resolutions. Bottom: Comparison of the mean amplitude (size in pA) of synaptic current produced by FX- and control-iNs; p < 0.02. Frequency of spontaneous EPSC in FX- and control-neurons. The decay time of individual EPSCs. Comparisons were made between voltage-clamped (n = 33) WT neurons and FX-iNs (n = 21), held at −60 mV.
MEA is a complementary approach to patch-clamp and is a gold standard for characterizing the extracellular potential of neurons and assessing the neuronal network activity in vitro. Researchers used high-throughput MEAs which were recorded in 24-well-plates (12 electrodes/well) to compare the neural network activity of FX and WT neurons plated at the same density (105 neurons/well) and conditions (Figure 4A; for details, see in Methods: MEA recordings and data analysis). Blockage of APs by Tetrodotoxin (TTX; 0.5 µM) indicates that APs in both FXS and WT neurons results from the activation of Na+ currents (Figure 4B). Neurons were recorded at eight different time points during differentiation (16–44 DIV), and spontaneous APs were detected using a voltage-threshold-based algorithm (Figure 4C). Approximately 60% of electrodes were active in both lines, indicating sustained network activity in both FX and WT neurons. The increase in the number of active electrodes with age in culture may represent neuronal maturation within the duration of the experiment, reaching an optimal activity of >65% of all electrodes on >30 DIV. Electrode recording from FX neurons showed a nearly two-fold higher activity (73%) compared to WT (42%) (>23 DIV; p < 0.0001), and the results were confirmed with the high number of recorded neurons (Figure 4C).
Figure 4. Spontaneous activity of induced neurons (iNs) recorded with MEA: (A) High-throughput MEA for recording spontaneous activity (24-well plate; 12 electrodes/well) with representative images of WT and FX iNs. Note the similar distribution of iNs at 10 days post induction captured from the center of a single MEA well. Scale bar: 100 µm. (B) representative traces of WT-iNs at 51 DIV (left) with spikes totally suppressed by the application of the Na+ channel blocker Tetrodotoxin (TTX; 0.5 µM; right). The red (top) and blue (bottom) lines are the thresholds (4.5 times the standard deviation of the noise); (C) percentage of active electrodes at four time points post-neuronal-induction (16–19, 23–26, 30–33, 37–44 DIV). The same batch of neurons was recorded for all time points. n = 22–36 wells were sampled for each line at each time point. Data are represented as mean ± SEM. One-way ANOVA: **** p < 0.0001.
For further analysis, they included only electrodes that discharged ≥15 spikes/minute. Collectively, the FX-iNs were hyperexcitable compared to the control as seen by the significantly higher number of spikes per active electrode in FX neurons (at 30–44 DIV; Figure 5A,B). Inter-spike interval (ISI) analysis demonstrated no differences in ISI between the two lines (Figure 5C,D), probably due to the variable spike frequencies in different neurons constructing both FX and WT neuronal networks.
Figure 5. Neuronal hyperexcitability in FX-iNs: (A) Number of spikes per active electrode. Active electrodes were defined by presenting ≥15 spikes/minute; (B) representative raster plot of spikes from more than 20 electrodes at 37 DIV FX and WT-iNs, showing the pattern and density of spikes during 60 s recording. In all graphs, the voltage-threshold-based algorithm runs over 20000 Hz high-pass-filtered traces, and is adjusted to 4.5 times the standard deviation of the noise; (C) inter-spike interval (ISI). Average ISI below 0.5 s and (D) above 0.5 s. Three independent experiments were performed. The same batch of neurons was recorded for all time points. n = 20–48 wells were sampled for each line at each time point. Data are represented as mean ± SEM. One-way ANOVA: * p < 0.05, ** p < 0.01.
Since neurons in the network tended to fire in bursts, researchers analyzed and compared the burst-firing activity between the FXS and WT neurons (Figure 6). The average rate of burst discharges was significantly higher in FXS-iNs compared to WT (64 vs. 36 bursts/ 3 min, respectively, at 37–44 DIV, p < 0.0001; Figure 6A,B). These results are correlated with the inter-burst interval (IBI) being decreased in mature FXS neurons (>37 DIV). Interestingly, at an earlier stage of differentiation (16–19 DIV), the results were the opposite, with an increased IBI observed in FX neurons compared to WT, probably reflecting their delayed development and maturation (Figure 6C). In accordance with the above, the averaged number of spikes in a burst was also significantly higher in the FXS-iNs compared to WT (>30 DIV; p < 0.0001; Figure 6D). The duration of bursts, however, was not different between FX and controls (Figure 6E). These results indicate that FX neuronal networks display a less mature electrical activity pattern compared to WT.
Figure 6. Higher spontaneous burst activity in FX neurons: (A) representatives burst traces of WT (blue) and FX (red) at 37 DIV during 10 s recording (top) and higher resolution of the first 2 s from each trace (bottom); (B) number of bursts; (C) intra-burst interval (IBI); (D) number of spikes in a burst; (E) burst duration. Bursts were defined as >3 spikes with <300 ms inter-spike interval. Three independent experiments were performed at each time point (16–19, 23–26, 30–33 and 37–44 DIV). n = 50–235 electrodes were sampled for each line in each time point. Data are represented as mean ± SEM. One-way ANOVA: * p < 0.05, ** p < 0.01, **** p < 0.0001.
To further compare the functional connectivity between FXS and WT networks, a cross-correlation (CC) analysis was performed (Figure 7). CC measures the distance in time between two spikes generated by different neurons in the same well. If the spikes of one neuron tend to occur at the same time as the spikes of another neuron, the CC score will be close to 1 (red square in the heat map), indicating synchronous network activity. However, if the spikes generated by different neurons constructing the network occur at different time points from the network, activity is considered non-synchronized, and the CC score will be closer to zero (blue square in the heat map of Figure 7B). They analyzed the CC of both FXS and control iNs at four time points during development (23, 30, 37 and 44 DIV). The results show that although both FX and WT neurons mature in culture and their network activity becomes more synchronous (higher CC), in FX, the neuronal network is significantly less synchronous than that of control (0.17 vs. 0.25, respectively, at 44 DIV; p < 0.004; Figure 7C).
Figure 7. Synchrony of neuronal network activity: (A) Representative synchronous (WT, left) and asynchronous (FX, right) electrical activity between two electrodes, (AU—arbitrary unit). The signals shown are spike locations that have been convoluted with Gaussian (sigma = 0.5); (B) heat maps summarizing cross-correlation (CC) among all active electrodes (≥15 spikes/min). The panel ranges from high synchrony/high CC among neurons (red), to low synchrony/ low CC (blue). Diagonal red patterns represent autocorrelation within the recorded neurons; (C) averaged values of CC electrodes. The threshold was determined to be >0.1. The same batch of neurons was recorded for all time points. Two independent experiments were performed at each time point (23, 30,37 and 44 DIV). A total of 35–207 pairs of electrodes were analyzed for each line in each time point. Data are represented as mean ± SD. T-test: * p < 0.05; ** p < 0.005.