- Open Access
Modeling neurological diseases with induced pluripotent cells reprogrammed from immortalized lymphoblastoid cell lines
Molecular Brain volume 9, Article number: 88 (2016)
Patient-specific induced pluripotent stem cells (iPSCs) facilitate understanding of the etiology of diseases, discovery of new drugs and development of novel therapeutic interventions. A frequently used starting source of cells for generating iPSCs has been dermal fibroblasts (DFs) isolated from skin biopsies. However, there are also numerous repositories containing lymphoblastoid B-cell lines (LCLs) generated from a variety of patients. To date, this rich bioresource of LCLs has been underused for generating iPSCs, and its use would greatly expand the range of targeted diseases that could be studied by using patient-specific iPSCs. However, it remains unclear whether patient’s LCL-derived iPSCs (LiPSCs) can function as a disease model. Therefore, we generated Parkinson’s disease patient-specific LiPSCs and evaluated their utility as tools for modeling neurological diseases. We established iPSCs from two LCL clones, which were derived from a healthy donor and a patient carrying PARK2 mutations, by using existing non-integrating episomal protocols. Whole genome sequencing (WGS) and comparative genomic hybridization (CGH) analyses showed that the appearance of somatic variations in the genomes of the iPSCs did not vary substantially according to the original cell types (LCLs, T-cells and fibroblasts). Furthermore, LiPSCs could be differentiated into functional neurons by using the direct neurosphere conversion method (dNS method), and they showed several Parkinson’s disease phenotypes that were similar to those of DF-iPSCs. These data indicate that the global LCL repositories can be used as a resource for generating iPSCs and disease models. Thus, LCLs are the powerful tools for generating iPSCs and modeling neurological diseases.
Modeling with induced pluripotent cells (iPSCs) is highly useful in research of neural diseases [1–4] because it is difficult to obtain patient-derived cells from the central nervous system that recapitulates the disease pathology. Dermal fibroblasts (DFs) have widely been used as a source of patient-specific iPSCs because they were used to derive the first human iPSCs . However, the isolation of DF cell lines requires invasive skin biopsies from patients. Recently, we have reported that T-cell derived iPSCs (TiPSCs) and fibroblast-derived iPSCs can be used to model neural diseases by using robust neural induction protocols . TiPSCs can be established from a small amount of peripheral blood, such that patients can provide samples through a minimally invasive procedure . Although peripheral blood cells, including T-cells, appear to be ideal sources for generating iPSCs, cloned and frozen T-cells are unstable, owing to their variable reprogramming efficiency depending on the donors and conditions of each sample (e.g., culture medium or number of passages) [8–10]. In particular, efficient wide-scale parallel reprogramming of samples requires a stable source of cells to achieve high-throughput processing and to model polygenic or sporadic diseases. Lymphoblastoid B-cell lines (LCLs) are stable peripheral B-cell lines that are transformed by infection with Epstein-Barr virus (EBV). LCLs are easy to maintain, and various types of well-characterized LCL clones established from patients are already available in worldwide repositories and are usually linked to patient clinical history, long-term genotype and phenotype data, and molecular/functional studies of various diseases [11, 12]. LCL banks have been of great importance in providing reference material for rare genetic diseases, as well as in managing large amounts of DNA for the genetic analysis of complex conditions in population- and family-based disease collections [13, 14]. A number of major facilities currently establish and manage extensive collections of cell lines (https://catalog.coriell.org/; http://ja.brc.riken.jp/; https://www.eagle-i.net/; www.ecacc.org.uk; www.alspac.bris.ac.uk; www.lgcpromochem-atcc.com; www.rutgers.edu) for the international research community.
To date, five groups have reported the successful reprogramming of LCLs from healthy donors and patients into iPSCs [15–19]. However, it is unclear whether LCL-derived iPSCs (LiPSCs) can be used for disease-specific analyses. Although immortalization of B-cells by EBV infection  may be a problematic procedure in terms of the cellular characteristics of LiPSCs, little is known about the effects of EBV infection on the properties of LiPSCs or even LCLs. Although EBV transformation commonly maintains the genome composition of the original cells , long-term LCL culture may cause several aberrations, including genomic aberrations such as chromosomal aneuploidy, down-regulation of p16/Rb, mutation of the p53 gene, modulation of apoptosis and sensitivity to various chemical agents . Therefore, detailed analyses of the genomic structure and cellular properties of LCLs and LiPSCs are necessary because unexpected genomic mutations would confound the disease-specific phenotypes in disease models generated with iPSCs. Additionally, some genomic aberrations and/or epigenetic memories derived from the source cells can affect the differentiation ability of iPSCs [23–25].
Although neurological disease models generated with LiPSCs would accelerate the progress of patient-specific iPSC studies, these concerns regarding the effects of EBV must be elucidated. Here, we established LiPSCs from both healthy donor and patient with the mutation of PARK2 (parkin) known as one of the causative genes for familial Parkinson’s disease,) to evaluate the characteristics of LiPSCs, including whole genomic sequencing and to confirm the utility of LiPSCs as tools for modeling neurological diseases.
Establishment and characterization of iPSCs derived from LCLs
Previous reports have shown that reprogramming LCLs to iPSCs by using existing non-integrating episomal protocols provided no identifiable iPSC clones, even after 35–40 days [17, 18]. Recently, R. Barrett and colleagues have developed a new episomal (OriP/EBNA1) plasmid reprogramming method for LCLs, which reduces the time required to generate LiPSCs and enhances the efficiency of reprograming LCLs into iPSCs . We modified this method for on-feeder cultures and applied it to LCLs from two donors: a healthy person, “KA” and a patient with autosomal recessive juvenile Parkinson’s disease (PARK2), “PB” (Fig. 1a). LiPSCs were generated with a similar efficiency to that of a previously described method , and the reprogramming efficiencies of LKA and LPB were 0.002 and 0.0015 %, respectively. After optimizing the protocol for generating iPSCs from LCLs, we established several LiPSC clones, including three healthy donor-derived LiPSC lines, LKA10, LKA29, and LKA36 and four PARK2 patient-derived LiPSC lines, LPB1, LPB3, LPB7 and LPB8, which had a normal karyotype (Additional file 1: Figure S1A). In addition, we used four DF-iPSC lines derived from the same healthy donor, eKA3, eKA4, KA11 and KA23, as the control reference lines  in this study.
Representative morphologies of LiPSC colonies are shown in Fig. 1b. All LiPSC lines exhibited typical PSC characteristics, including tightly packed colonies, a high cell nuclear-cytoplasmic ratio, and the production of the surface and nuclear pluripotency proteins, TRA1-60 and OCT4 and were indistinguishable from the DF-iPSC lines (Fig. 1b). LiPSCs also expressed endogenous pluripotency genes at a similar level to DF-iPSCs (Figs. 1c–e). These data indicated that LiPSCs and DF-iPSCs were indistinguishable in terms of morphology and the expression of pluripotent markers at both the protein and mRNA levels (Fig. 1b–e). In addition, we confirmed whether LiPSCs have differentiation potentials into three-germ layers by in vitro differentiation analysis via EB (Additional file 1: Figure S1B and C). Thus, all LiPSC clones had a pluripotency. Because EBNA-1 has been reported to be required for the establishment of persistent EBV infection and survival of host B-cells , we next examined the expression of EBNA1 and additional EBV-related genes (EBNA-2, BZLF-1, LMP-1 and OriP) in LiPSCs. A PCR analysis of genomic DNA showed that all the EBV-related latency elements were eventually eliminated from the established LiPSCs, suggesting the loss of EBV-associated elements as a result of the reprogramming process (Fig. 1f).
The global gene expression profiles in iPSCs at passages below 20 were evaluated with a microarray analysis to explore the detailed differences between LiPSCs, TiPSCs and DF-iPSCs resulting from the origin of the iPSCs. With the exception of the genes that were expressed at low levels in all samples, the data were normalized and subjected to principal component analysis (PCA) (Fig. 1g) and hierarchical clustering (Fig. 1h). Although the LiPSCs, TiPSCs and DF-iPSCs were relatively close to each other in the PCA analysis (Fig. 1g), hierarchical clustering placed the LiPSCs, TiPSCs and DF-iPSCs into different groups (Fig. 1h). These data also suggested that the original cell type influences the properties of hiPSCs.
Immortalization by EBV does not affect the number of de novo mutations and structural variations in LiPSCs
We performed array-based comparative genomic hybridization (aCGH) and whole genome sequence (WGS) analyses to examine the somatic structural variations (SVs) and single nucleotide variations (SNVs) in LiPSCs (Fig. 2a). A comparison of the genomes of the LiPSC clones and LCLs revealed a deletion (233,645 bp) at 19p13.3 in all the LiPSC clones examined from the healthy donor, KA (Fig. 2b–d). Although the number of reads was limited (approximately 6 % of the total reads), the presence of the reads spanning the breakpoint was confirmed not only in the LiPSC clones but also in the LCLs (Fig. 2e), thus strongly suggesting that the deletions in the LiPSC clones were already present in a subpopulation of their original cells, LCLs and were not detected by the aCGH analysis.
We then searched for the appearance of SNVs by comparing the LiPSC clones with LCLs by using the WGS analysis. In this analysis, we focused on non-synonymous SNVs in coding regions, and all the variations were further validated by direct nucleotide sequence analysis. The analysis revealed the appearance of 4–6 non-synonymous variations in the LiPSCs clones derived from KA. The analysis revealed the appearance of 9–12 non-synonymous SNVs in the PB-LiPSC clones compared with the original cell source, LCLs (Fig. 2a).
We identified a somatic mutation in SLC26A5 (rs758296903) in all the LiPSC clones from KA (LKA10, LKA29 and LKA36) compared with their original cells. A detailed examination of the short reads revealed that 8.6 % of the reads from the LCLs carried the mutation, thus indicating that the SNV was already present in a subpopulation of the original cells. Other variations observed in the TiPSCs and DF-iPSCs compared with T-cells and DFs, respectively, are shown in Fig. 2a and Additional file 1: Figure S1. These data indicate that the reprogramming and/or immortalization processes might cause some somatic mutations. However, the total number of somatic mutations observed in the genomes of iPSCs compared with their corresponding original cells did not vary among the cell sources of origin (LCLs, T-cells and fibroblasts: Fig. 2a), thus suggesting that LiPSCs may also have the same properties and functions as hiPSCs, similarly to TiPSCs and DF-iPSCs.
Differentiation of neural cells from LiPSCs through a direct neurosphere conversion method
In our previous reports, we have shown that TiPSCs are poorly differentiated into the ectodermal lineage, and it is difficult to induce TiPSCs to differentiate into neuronal cells via EB formation by using spontaneous neuronal differentiation protocols. To overcome this limitation, we developed a neurosphere (NS)-based differentiation method (direct NS conversion method: dNS method). TiPSCs were differentiated into neural cells with similar efficiency as DF-iPSCs with the dNS method . We used the dNS method to differentiate the LiPSCs into neural cells (Fig. 3a) because both TiPSCs and LiPSCs were derived from peripheral blood cells and supposedly have similar differentiation propensities. With the dNS method, LiPSCs efficiently formed similar number of NSs (Fig. 3b) compared with DF-iPSCs (Fig. 3c). Morphological analysis of NSs focusing on their circularity (Fig. 3d) and diameter (Fig. 3e) also demonstrated that there were no significant differences between LiPSC- and DF-iPSC-NSs (Figs. 3d and e). We quantified the expression of pluripotent markers, NANOG and OCT4 and neural markers, PAX6 and NESTIN, in LiPSC-derived NSs by qPCR analysis. Although the expression of NANOG in LiPSC-derived NSs was lower in all clones and DF-iPSC-derived NSs (Fig. 3d), OCT4 expression was maintained in some LiPSC-derived NS clones and DF-iPSC-derived NS clones (Fig. 3e). However, these differences in NANOG (Fig. 3h) and OCT4 expression (Fig. 3i) were not statistically significant among the cell types of origin. The expression of neural stem cell markers did not differ significantly among the various NSs we tested, regardless of the original cell type and donor (Fig. 3f, g, j and k). These results indicate that LiPSCs can differentiate into NSs via the dNS method as efficiently as DF-iPSCs.
LiPSCs can be differentiated into various types of functional neurons to a similar level as DF-iPSCs by using the dNS method
We next examined whether LiPSC-derived NSs could be induced to differentiate into neural cells as efficiently as DF-iPSCs. LiPSC-derived NSs were dissociated into single cells and cultured for neural differentiation. Thirteen days later, the differentiated cells were immunostained with antibodies against a neural marker, βIII-TUBULIN and an astrocyte marker, glial fibrillary acid protein (GFAP) (Fig. 4a). LiPSC-derived NSs differentiated into neurons and astrocytes at a similar ratio as DF-iPSCs, indicating that LiPSC-derived NSs had similar propensities to differentiate into neurons and astrocytes as the DF-iPSC-derived NSs (Fig. 4b–d). The immunocytochemical analysis revealed that these βIII-TUBULIN-positive neurons included tyrosine hydroxylase (TH)-positive dopaminergic neurons (Fig. 4e and f), gamma-aminobutyric acid (GABA)-positive GABAergic neurons (Fig. 4e and g), and vesicular glutamate transporter 1-positive glutamatergic neurons (Fig. 4e and h). Figs. 4d and i show a summary of the data from these differentiated cell types and neuronal subtypes, suggesting that LiPSCs differentiated into various types of cells as efficiently as DF-iPSCs, regardless of the differences among donors.
We recorded voltage-sensitive currents in 60-day-old LiPSC-derived neurons to confirm their electrophysiological properties. LiPSC-derived NSs were infected with a lentivirus expressing a human Synapsin I promoter-driven GFP (CSIV–hSynI-GFP-IRES2-NeoR) soon after single dissociated NS cells were plated . After neuronal maturation, voltage-dependent Na+ and K+ currents and TTX-sensitive voltage-gated membrane currents were detected (Fig. 4j and k), and action potentials were elicited in the majority of the LiPSC-derived neuronal cells analyzed by depolarizing the membrane in current clamp mode (Fig. 4l). These results suggested that LiPSCs can be differentiated into functional neurons through the dNS method.
Neurons derived from PARK2-LiPSCs exhibited impaired mitochondrial activity and phenotypes
We differentiated LiPSCs derived from a patient with parkin mutations, a familial form of Parkinson’s disease, PARK2, into neurons to determine whether hiPSCs established from LCLs could be used as a model of neurological disease. A consistent neurochemical abnormality found in Parkinson’s disease is the degeneration of dopaminergic neurons in the substantia nigra. Therefore, we exposed single dissociated hiPSCs to FGF8, sonic hedgehog (Shh), purmorphamine (PMA) and CHIR99021 during the NS formation, as described in our previous study , to generate midbrain dopaminergic neuron (mDAN)-enriched culture. After neuronal maturation for 13 days in culture, immunostaining of the neural epithelial cells demonstrated that nearly 20 % of the βIII-TUBULIN-positive cells expressed the dopaminergic neuron marker TH, with a similar differentiation ratio between the DF-iPSCs and LiPSCs (Fig. 5b and c). These data indicate that the LiPSCs and DF-iPSCs were able to differentiate into mDANs with similar efficiencies (Fig. 5a–c).
A PCR analysis of the genomic DNA confirmed that PARK2 patient-derived LiPSCs (LPBs; LPB1, LPB3 and LPB7) contained a homozygous deletion in exons 6 and 7 of the parkin gene (Fig. 5d), which encodes a component of an E3 ubiquitin ligase involved in mitochondrial homeostasis . We have previously reported that hiPSCs established from DFs of a PARK2 patient exhibited abnormal turnover of damaged mitochondria  and increased production of reactive oxygen species (ROS) in their neurons. Therefore, we treated the LiPSC-derived neurons with carbonyl cyanide m-chlorophenyl hydrazone (CCCP), which triggers the loss of mitochondrial membrane potential and results in the removal of damaged mitochondria. We visualized the area of the inner mitochondrial membrane (IMM) by using an antibody against the IMM marker Complex-III Core I (C-III Core I) to determine the extent to which the damaged mitochondria were eliminated after CCCP treatment. Compared with the untreated samples, the CCCP-treated samples elicited a dramatic decrease in the IMM area in the control neurons (LKA clones), but not in the PARK2 neurons (LPB clones) (Fig. 5e and f). Using the same CCCP-treated neuronal samples, we performed immunostaining with an anti-TH antibody to quantify the ratio of TH-positive dopaminergic neurons. In both the DF-iPSC- and LiPSC-derived neurons, we observed significant decreases in the numbers of TH-positive neurons in PARK2-derived cells (LPB clones) after the CCCP treatment, indicating that dopaminergic neurons derived from PARK2-iPSCs were more vulnerable to mitochondrial stress than control-iPSCs (Fig. 5g).
Finally, we evaluated ROS production in the neurons derived from the control- and PARK2-LiPSCs by using the CellROX® Green Reagent, which is weakly fluorescent in a reduced state and exhibits bright green photostable fluorescence after oxidation by ROS, with absorption/emission maxima of ~ 485/520 nm. The reactive CellROX fluorescence was significantly increased in the PARK2 neurons, and this phenotype prominently appeared in the TH-positive neurons, indicating that ROS production was increased in the PARK2 neurons, particularly the mDANs (Fig. 5h and i). These phenotypes were observed in both DF-iPSC- and LiPSC-derived neurons (Figs. 5e–i). These results strongly suggest that the hiPSCs derived from patient LCLs can be used as models of neurological diseases.
Patient-specific iPSC technology makes it possible to recapitulate disease phenotypes in vitro, thus significantly facilitating the elucidation of disease processes and the development of therapeutic drugs. However, invasive skin biopsies are generally performed to obtain a patient’s original cells, DFs and generate iPSCs. Several groups have succeeded in establishing iPSCs from peripheral blood cells, T-cells and LCLs, thus extending the use of iPSC technologies [15–19, 30, 31]. In particular, our previous study has demonstrated the utility of TiPSCs as tools for modeling diseases. However, it had been unclear whether LiPSCs can function as a disease model. Because LCLs are stable cell lines that are easy to culture, they are expected to be the ideal resource for efficient wide-scale parallel reprogramming to achieve high-throughput processing for modeling polygenic or sporadic diseases. Moreover, various types of LCL clones established from patients are already available and stored in global repositories . Therefore, it would be very important to obtain experimental proof of the utility of LiPSCs as disease models to increase the number of cases, perform disease-specific analyses based on statistical analysis, and expand the target-disease areas for iPSC technology.
In addition to this study, several groups have established iPSCs from peripheral blood cells and have investigated their characteristics as iPSCs [6, 7, 15–19, 30, 31]. However, little is known about the effects of the differences among original cells on genomic mutations in the derived iPSCs. In particular, LiPSCs are highly influenced by the original cells because LCLs are generated through an EBV infection-mediated immortalization process. This study clearly showed that the differences among LCLs compared with the original cells had no substantial effects on the number of somatic single nucleotide variations and structural variations in the established iPSCs compared with those obtained from T-cells or DFs (Fig. 2). These data broaden the applications of LiPSCs, and strongly support the utility of LCLs as original cells for generating iPSCs. Moreover, our results indicated that the reprogramming processes can cause some somatic mutations, regardless of the cell type of origin (Fig. 2 and Additional file 2: Figure S2). Although some of these mutations have recently been described , further studies are needed to clarify their effects.
In recent studies of human clinical genetics, progress has been made in the identification of the genes responsible for neurological disorders [33, 34] and the genes increasing the onset risk of neuropsychiatric disorders . However, there are many diseases for which the causative genes and genetic risk factors have not yet been identified. Recent whole exome sequencing studies have revealed that these sporadic cases are related to many variants, indicating the cumulative genetic mutations that trigger the onset of the disease. It is important to establish hiPSCs from a sufficient number of patients and characterize multiple clones for statistical analyses to elucidate the pathogenic mechanisms of these diseases. Well-characterized LCLs are already available in global repositories and are linked to patient clinical history, long-term genotype and phenotype data, and molecular/functional studies. Then, the clarification of phenotypes using hiPSC disease model could elucidate the pathological mechanism of such sporadic diseases, and provide a path to the identification of their related genes. Therefore, in addition to T-cells, LCLs are a suitable source for generating hiPSCs to model sporadic or common diseases.
Human iPSC generation from LCLs
Peripheral blood mononuclear cells (PBMCs) were obtained from a healthy donor, “KA” and a patient with autosomal recessive juvenile Parkinson’s disease (PARK2), “PB”, by centrifuging heparinized blood over a Ficoll-Paque PREMIUM (GE Healthcare, Chicago, USA) gradient. KA- or PB-derived PBMCs were immortalized by an EBV infection according to the protocol of SRL Medisearch Incorporation and were transformed into LCLs. LCLs were cultured in RPMI 1640 (Gibco, Massachusetts, USA) supplemented with 10 % fetal bovine serum (FBS) at 37 °C and 5 % CO2 in a humidified incubator. After several passages, the LCLs were electroporated with the Neon™ Transfection System 100 μL Kit (MPK10096; Thermo Fisher Scientific, Massachusetts, USA) using 1.0 μg of each episomal plasmid (Addgene, Cambridge, USA) expressing 6 factors: OCT4, SOX2, KLF4, l-MYC, LIN28 and p53 shRNAs (pCXLE-hOCT3/4-shp53, pCXLE-hSK and pCXLE-hUL), according to the manufacturer’s instructions. The transfected LCLs were rapidly transferred to a 6-well plate at a density of 2.0 × 106 cells/well and incubated for 24 h. At 24 h after electroporation, the medium was replaced with hiPSC medium. After an additional 24 h, the cells were passaged to a 100 mm dish containing mitomycin C-inactivated mouse SNL feeder cells at a density of 5.0 × 104–5.0 × 105 cells/dish and cultured in hiPSC medium, which was changed every other day until colonies were picked. The generated hiPSCs were maintained on mitomycin C-inactivated mouse SNL feeder cells in hiPSC medium. All hiPSC lines analyzed in this study were between passage 7 and 20.
All human primary cells were used after appropriate written informed consent was given to the commercial providers. All experimental procedures for biopsy and reprogramming were approved by the Ethics Committee of the Keio University School of Medicine (No. 20080016).
Neural differentiation in vitro
The dNS method, which we have previously reported, was used for the neural differentiation of hiPSCs . For neural induction, hiPSCs were dissociated into single cells by incubation with TrypLETM Select (Life Technologies, Massachusetts, USA) for 5 min and pipetting. The cells were cultured at a density of 10 cells/μL in a T25 flask (Nunclon, Massachusetts, USA) in MHM supplemented with B27 (Gibco, Massachusetts, USA), 20 ng/mL FGF-2 (Wako, Osaka, Japan), 10 μM Y-27632 (Wako, Osaka, Japan) and 10 ng/mL hLIF (Nacalai Tesque, Kyoto, Japan) in 4 % oxygen for 7 days. NSs were repeatedly passaged by dissociation into single cells, and then cultured at a density of 50 cells/μL in the same manner as the primary sphere formation. NSs were used at passage 3 for analysis. For terminal differentiation, the dissociated NSs were plated on PO- (Sigma-Aldrich, Missouri, USA) and fibronectin- (Sigma-Aldrich, Missouri, USA) coated coverslips and cultured in MHM containing B27 (Gibco, Massachusetts, USA), 10 ng/mL brain-derived neurotrophic factor (BDNF; R&D Systems, Minnesota, USA), 10 ng/mL glial cell-differentiated neurotrophic factor (GDNF; R&D Systems, Minnesota, USA), 200 μM ascorbic acid (Sigma-Aldrich, Missouri, USA) and 1 mM dibutyryl-cAMP (Sigma-Aldrich, Missouri, USA) for 10–60 days.
Immunocytochemical analysis of hiPSCs and neurons derived from hiPSCs
Cells were fixed in phosphate-buffered saline (PBS) containing 4 % paraformaldehyde (PFA) for 30 min at room temperature. Thereafter, all cells were blocked with 5 % FBS and Triton X-100 and incubated with the primary antibodies described in Additional file 3: Table S1. The cells were then rinsed with PBS, and incubated with species-specific Alexa Fluor 488-, Alexa Fluor 555- or Alexa Fluor 647-conjugated secondary antibodies (1:500; Invitrogen, Massachusetts, USA); this was followed by Hoechst 33258 (0.5 μg/mL; Sigma-Aldrich, Missouri, USA) to counterstain the nuclei. The images were obtained using a universal fluorescence microscope (Axioplan2; Carl Zeiss AG, Oberkochen, Germany) or a confocal laser scanning microscope (LSM-710; Carl Zeiss AG, Oberkochen, Germany).
CGH array analysis
Genomic DNA was extracted using a DNeasy Blood & Tissue Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s protocol. Five hundred micrograms of genomic DNA was subjected to the Human CGH array (4X180K; Agilent Technologies, California, USA) according to the manufacturer’s protocol. All combinations (original cells and iPSCs) were analyzed. The data were analyzed using the Agilent Genomic Workbench. Rearrangements involving immunoglobulin gene regions and T-cell receptor gene regions that are rearranged in T-cell and B-cell lineages, respectively, were not taken into consideration in the analysis of the structural variations.
Whole genome sequence analysis
One microgram of genomic DNA was subjected to whole genome sequence analysis. Sequencing libraries were constructed using a TruSeq DNA PCR-free library preparation kit (Illumina, California, USA). The sequences were analyzed using HiSeq2500 (Illumina, California, USA) in the rapid mode (150 bp, paired end). After alignment to the reference genome (hg19) using BWA  and removal of multiply aligned reads and duplicate reads, basecalls were performed with SAMtools . Somatic variants were detected as previously described [38, 39]. We confirmed that all of the candidate somatic mutations consisting of single nucleotide substitutions were predicted to result in changes in the amino acid sequence by direct nucleotide sequence analysis. The SNVs inside T-cell receptor regions and immunoglobulin regions were excluded. If necessary, whole genome amplification using REPLI-g (QIAGEN, Hilden, Germany) was applied to the direct nucleotide sequence analysis.
Total RNA was extracted with an RNAeasy Kit (QIAGEN, Hilden, Germany) and the RNA quality was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, California, USA). Total RNA (100 ng) was reverse-transcribed, labeled with biotin using a 3’IVT Express Kit (Affymetrix, California, USA) and hybridized to a GeneChip® Human Genome U133 plus 2.0 Array (Affymetrix, California, USA). The arrays were washed and stained using a GeneChip Fluidics Station 450 (Affymetrix, California, USA) and then scanned with the GeneChip Scanner 3000 7G System (Affymetrix, California, USA) according to the manufacturer’s instructions. The raw probe intensity files were MAS5-normalized and log (base2) transformed by using GeneSpring GX 13.1 software (Agilent Technologies, California, USA). The gene set was filtered based on the expression levels to remove the genes that were not expressed in all samples. Principal component analysis (PCA) was performed using the normalized data. For the hierarchical clustering, the normalized data were calculated on the basis of Euclidean correlations with average linkages.
We used NCBI GEO microarray data GSE76832 for T-cell, DF and iPSCs established from them (T-cell(KA), DF(KA), TKA4, TKA9, eKA3 and KA11) .
Total RNA was isolated using an RNeasy Mini Kit (QIAGEN, Hilden, Germany) and 1 μg of RNA was used to generate cDNAs by using a reverse transcription (RT) system (Promega, Wisconsin, USA). RT-PCR was performed as previously described . Values were normalized to ACTB. Quantitative RT-PCR (qPCR) was performed on an ABI PRISM Sequence Detection System 7900HT (Applied BioSystems, Massachusetts, USA) by using SYBR premix ExTaq Tli RNaseH Plus (Takara Bio, Shiga, Japan). The primers are described in Additional file 4: Table S2.
PCR amplification of genomic DNA
Genomic DNA was purified from HDFs, LCLs and hiPSCs using a DNeasy Kit (QIAGEN, Hilden, Germany). The PCR conditions were described previously (Additional file 4: Table S2) .
For the electrophysiology experiments, the culture medium was replaced with a physiological solution (118 mM NaCl, 2.5 mM KCl, 26 mM NaHCO3, 1 mM NaH2PO4, 10 mM glucose, 4 mM MgCl2, and 4 mM CaCl2). Tetrodotoxin (TTX, 1 μM) was bath-applied. The electrodes (5–8 MΩ) were filled with whole-cell pipette solution (120 mM potassium acetate, 20 mM KCl, 0.1 mM CaCl2, 5 mM MgCl2, 0.2 mM EGTA, 5 mM ATP, and 10 mM HEPES, pH 7.3) [42, 43]. The whole-cell recordings of GFP-expressing neurons were configured using an EPC-7 amplifier (HEKA Elektronik Dr. Schulze GmbH, Lambrecht/Pfalz, Germany) and a Digidata 1200 acquisition board (Axon Instruments, California, USA). The membrane potential was clamped at −60 mV. Membrane resistance (Rm), series resistance (Rs) and membrane capacitance (Cm) were monitored. Only recordings with Rm > 100 MΩ and Rs < 20 MΩ were included in the analysis.
Carbonyl cyanide m-chlorophenylhydrazone treatment
Neurons were cultured with 30 μM carbonyl cyanide m-chlorophenyl hydrazone (CCCP; Sigma-Aldrich, Missouri, USA) or DMSO for 48 h. The cells were then fixed, stained for βIII-TUBULIN and Complex-III Core I (CIII-Core I), and counterstained with Hoechst. The cytoplasmic area was extracted to quantify the IMM area of neurons, as shown in Fig. 5e. The IMM area of the neurons was quantified from the digitized values using IN CELL Analyzer 6000 (GE Healthcare, Chicago, USA).
Oxidative stress analysis
The ROS levels were determined by measuring the CellROX fluorescence using the CellROX® Green Reagent to detect oxidative stress (Life Technologies, Massachusetts, USA). Briefly, the neurons were incubated with the CellROX Reagent for 30 min at 37 °C, after which they were washed with PBS and then fixed with 4 % PFA for 30 min at room temperature. Then, the cells were incubated with the primary antibody against βIII-TUBULIN (1:1000; Sigma-Aldrich, Missouri, USA) overnight at 4 °C, washed with PBS and incubated with an Alexa Fluor 555-conjugated secondary antibody (1:500; Invitrogen, Massachusetts, USA) for 1 h at room temperature. The fluorescence of the βIII-TUBULIN-positive neurons was measured by using an IN Cell Analyzer 6000 (GE Healthcare, Chicago, USA).
Comparative genomic hybridization
- dNS method:
Direct neurosphere conversion method
Induced pluripotent stem cells
Lymphoblastoid B-cell lines
Midbrain dopaminergic neuron
Reactive oxygen species
Single nucleotide variation
T-cell derived iPSC
Whole genome sequencing
Imaizumi Y, Okano H. Modeling human neurological disorders with induced pluripotent stem cells. J Neurochem. 2014;129(3):388–99.
Marchetto MC, Gage FH. Modeling brain disease in a dish: really? Cell Stem Cell. 2012;10(6):642–5.
Mattis VB, Svendsen CN. Induced pluripotent stem cells: a new revolution for clinical neurology? Lancet Neurol. 2011;10(4):383–94.
Okano H, Yamanaka S. iPS cell technologies: significance and applications to CNS regeneration and disease. Mol Brain. 2014;7:22.
Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, Yamanaka S. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131(5):861–72.
Matsumoto T, Fujimori K, Andoh-Noda T, Ando T, Kuzumaki N, Toyoshima M, Tada H, Imaizumi K, Ishikawa M, Yamaguchi R, et al. Functional Neurons Generated from T Cell-Derived Induced Pluripotent Stem Cells for Neurological Disease Modeling. Stem Cell Rep. 2016;6(3):422–35.
Seki T, Yuasa S, Fukuda K. Generation of induced pluripotent stem cells from a small amount of human peripheral blood using a combination of activated T cells and Sendai virus. Nat Protoc. 2012;7(4):718–28.
Klebanoff CA, Gattinoni L, Restifo NP. CD8+ T-cell memory in tumor immunology and immunotherapy. Immunol Rev. 2006;211:214–24.
Nishimura T, Kaneko S, Kawana-Tachikawa A, Tajima Y, Goto H, Zhu D, Nakayama-Hosoya K, Iriguchi S, Uemura Y, Shimizu T, et al. Generation of rejuvenated antigen-specific T cells by reprogramming to pluripotency and redifferentiation. Cell Stem Cell. 2013;12(1):114–26.
Wherry EJ. T cell exhaustion. Nat Immunol. 2011;12(6):492–9.
Sie L, Loong S, Tan EK. Utility of lymphoblastoid cell lines. J Neurosci Res. 2009;87(9):1953–9.
Wheeler HE, Dolan ME. Lymphoblastoid cell lines in pharmacogenomic discovery and clinical translation. Pharmacogenomics. 2012;13(1):55–70.
Beck JC, Beiswanger CM, John EM, Satariano E, West D. Successful transformation of cryopreserved lymphocytes: a resource for epidemiological studies. Cancer Epidemiol Biomark Prev. 2001;10(5):551–4.
Worthington J, Ollier WE, Leach MK, Smith I, Hay EM, Thomson W, Pepper L, Carthy D, Farhan A, Martin S, et al. The Arthritis and Rheumatism Council’s National Repository of Family Material: pedigrees from the first 100 rheumatoid arthritis families containing affected sibling pairs. Br J Rheumatol. 1994;33(10):970–6.
Barrett R, Ornelas L, Yeager N, Mandefro B, Sahabian A, Lenaeus L, Targan SR, Svendsen CN, Sareen D. Reliable generation of induced pluripotent stem cells from human lymphoblastoid cell lines. Stem Cells Transl Med. 2014;3(12):1429–34.
Bueno C, Sardina JL, Di Stefano B, Romero-Moya D, Munoz-Lopez A, Ariza L, Chillon MC, Balanzategui A, Castano J, Herreros A, et al. Reprogramming human B cells into induced pluripotent stem cells and its enhancement by C/EBPalpha. Leukemia. 2016;30(3):674–82.
Choi SM, Liu H, Chaudhari P, Kim Y, Cheng L, Feng J, Sharkis S, Ye Z, Jang YY. Reprogramming of EBV-immortalized B-lymphocyte cell lines into induced pluripotent stem cells. Blood. 2011;118(7):1801–5.
Rajesh D, Dickerson SJ, Yu J, Brown ME, Thomson JA, Seay NJ. Human lymphoblastoid B-cell lines reprogrammed to EBV-free induced pluripotent stem cells. Blood. 2011;118(7):1797–800.
Thomas SM, Kagan C, Pavlovic BJ, Burnett J, Patterson K, Pritchard JK, Gilad Y. Reprogramming LCLs to iPSCs Results in Recovery of Donor-Specific Gene Expression Signature. PLoS Genet. 2015;11(5):e1005216.
Amoli MM, Carthy D, Platt H, Ollier WE. EBV Immortalization of human B lymphocytes separated from small volumes of cryo-preserved whole blood. Int J Epidemiol. 2008;37 Suppl 1:i41–5.
Nilsson K. Human B-lymphoid cell lines. Hum Cell. 1992;5(1):25–41.
Sugimoto M, Tahara H, Ide T, Furuichi Y. Steps involved in immortalization and tumorigenesis in human B-lymphoblastoid cell lines transformed by Epstein-Barr virus. Cancer Res. 2004;64(10):3361–4.
Kim K, Doi A, Wen B, Ng K, Zhao R, Cahan P, Kim J, Aryee MJ, Ji H, Ehrlich LI, et al. Epigenetic memory in induced pluripotent stem cells. Nature. 2010;467(7313):285–90.
Kim K, Zhao R, Doi A, Ng K, Unternaehrer J, Cahan P, Huo H, Loh YH, Aryee MJ, Lensch MW, et al. Donor cell type can influence the epigenome and differentiation potential of human induced pluripotent stem cells. Nat Biotechnol. 2011;29(12):1117–9.
Polo JM, Liu S, Figueroa ME, Kulalert W, Eminli S, Tan KY, Apostolou E, Stadtfeld M, Li Y, Shioda T, et al. Cell type of origin influences the molecular and functional properties of mouse induced pluripotent stem cells. Nat Biotechnol. 2010;28(8):848–55.
Altmann M, Pich D, Ruiss R, Wang J, Sugden B, Hammerschmidt W. Transcriptional activation by EBV nuclear antigen 1 is essential for the expression of EBV’s transforming genes. Proc Natl Acad Sci U S A. 2006;103(38):14188–93.
Zhou Z, Kohda K, Ibata K, Kohyama J, Akamatsu W, Yuzaki M, Okano HJ, Sasaki E, Okano H. Reprogramming non-human primate somatic cells into functional neuronal cells by defined factors. Mol Brain. 2014;7:24.
Shimura H, Hattori N, Kubo S, Mizuno Y, Asakawa S, Minoshima S, Shimizu N, Iwai K, Chiba T, Tanaka K, et al. Familial Parkinson disease gene product, parkin, is a ubiquitin-protein ligase. Nat Genet. 2000;25(3):302–5.
Imaizumi Y, Okada Y, Akamatsu W, Koike M, Kuzumaki N, Hayakawa H, Nihira T, Kobayashi T, Ohyama M, Sato S, et al. Mitochondrial dysfunction associated with increased oxidative stress and alpha-synuclein accumulation in PARK2 iPSC-derived neurons and postmortem brain tissue. Mol Brain. 2012;5:35.
Loh YH, Hartung O, Li H, Guo C, Sahalie JM, Manos PD, Urbach A, Heffner GC, Grskovic M, Vigneault F, et al. Reprogramming of T cells from human peripheral blood. Cell Stem Cell. 2010;7(1):15–9.
Seki T, Yuasa S, Oda M, Egashira T, Yae K, Kusumoto D, Nakata H, Tohyama S, Hashimoto H, Kodaira M, et al. Generation of induced pluripotent stem cells from human terminally differentiated circulating T cells. Cell Stem Cell. 2010;7(1):11–4.
Okamura K, Toyoda M, Hata K, Nakabayashi K, Umezawa A. Whole-exome sequencing of fibroblast and its iPS cell lines derived from a patient diagnosed with xeroderma pigmentosum. Genom Data. 2015;6:4–6.
Kalinderi K, Bostantjopoulou S, Fidani L. The genetic background of Parkinson’s disease: current progress and future prospects. Acta Neurol Scand. 2016. doi:10.1111/ane.12563.
Li HF, Wu ZY. Genotype-phenotype correlations of amyotrophic lateral sclerosis. Trans Neurodegeneration. 2016;5:3.
Gatt JM, Burton KL, Williams LM, Schofield PR. Specific and common genes implicated across major mental disorders: a review of meta-analysis studies. J Psychiatr Res. 2015;60:1–13.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics (Oxford, England). 2009;25(14):1754–60.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The Sequence Alignment/Map format and SAMtools. Bioinformatics (Oxford, England). 2009;25(16):2078–9.
Isojima T, Doi K, Mitsui J, Oda Y, Tokuhiro E, Yasoda A, Yorifuji T, Horikawa R, Yoshimura J, Ishiura H, et al. A recurrent de novo FAM111A mutation causes Kenny-Caffey syndrome type 2. J Bone Miner Res Off J Am Soc Bone Miner Res. 2014;29(4):992–8.
Koh K, Ishiura H, Miwa M, Doi K, Yoshimura J, Mitsui J, Goto J, Morishita S, Tsuji S, Takiyama Y. Exome sequencing shows a novel de novo mutation in ATL1. Neuropsychiatry Clin Neurosci. 2014;2(1):1–4.
Okada Y, Matsumoto A, Shimazaki T, Enoki R, Koizumi A, Ishii S, Itoyama Y, Sobue G, Okano H. Spatiotemporal recapitulation of central nervous system development by murine embryonic stem cell-derived neural stem/progenitor cells. Stem Cells (Dayton, Ohio). 2008;26(12):3086–98.
Kitada T, Asakawa S, Hattori N, Matsumine H, Yamamura Y, Minoshima S, Yokochi M, Mizuno Y, Shimizu N. Mutations in the parkin gene cause autosomal recessive juvenile parkinsonism. Nature. 1998;392(6676):605–8.
Tada H, Kuroki Y, Funabashi T, Kamiya Y, Goto T, Suyama K, Sano A, Mitsushima D, Etgen AM, Takahashi T. Phasic synaptic incorporation of GluR2-lacking AMPA receptors at gonadotropin-releasing hormone neurons is involved in the generation of the luteinizing hormone surge in female rats. Neuroscience. 2013;248:664–9.
Tada H, Okano HJ, Takagi H, Shibata S, Yao I, Matsumoto M, Saiga T, Nakayama KI, Kashima H, Takahashi T, et al. Fbxo45, a novel ubiquitin ligase, regulates synaptic activity. J Biol Chem. 2010;285(6):3840–9.
Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14(2):178–92.
The authors thank all of the members of H.O.’s laboratory for their encouragement and support.
This work was supported by funding from the Project for the Realization of Regenerative Medicine and Support for Core Institutes for iPS Cell Research from the Ministry of Education, Culture, Sports, Science and Technology of Japan (MEXT) to H.O., the Research Center Network for Realization Research Centers/Projects of Regenerative Medicine (the Program for Intractable Disease Research utilizing disease-specific iPS Cells) from the Japan Science and Technology Agency (JST) and Japan Agency for Medical Research and Development (AMED) to W.A. and H.O., the New Energy and Industrial Technology Development Organization (NEDO) to W.A. and H.O., the Japan Society for the Promotion of Science (JSPS) to W.A., and a Grant-in-Aid for the Global COE Program from MEXT to Keio University.
Availability of data and materials
The data discussed in this publication have been deposited in the NCBI Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) database and are accessible through GEO Series accession numbers GSE76832  and GSE82159.
KF, TTe, TM, RH, ST, WA and HO conceived and designed the experiments. KF, TTe, HI, JM, KD, JY, HT, TM, MI, TTa, ST and WA performed the experiments and analyzed data. KF, TTe, HI, JM, KD, JY, ST, WA and HO wrote and edited the manuscript. RH, NH, SM and TTa contributed reagents, materials and analysis tools. All authors read and approved the final manuscript.
H.O. is a paid Scientific Advisory Board Member of SanBio Co., Ltd. T.M. is employed by Ajinomoto Co., Inc. M.I. is employed by Sumitomo Dainippon Pharma Co., Ltd.1
Consent for publication
Written informed consent to publish was obtained from the patient and the healthy donor.
Ethics approval and consent to participate
All human primary cells were used after appropriate written informed consent was given to the commercial providers. All experimental procedures for biopsy and reprogramming were approved by the Ethics Committee of the Keio University School of Medicine (No. 20080016).
Additional file 1: Figure S1.
In vitro differentiation analysis and karyotyping of LiPSC. (A) Representative karyotype of the established iPSC line from LCL-PB. (B) Quantitative analysis of immunocytochemistry for markers of three-germ layers; βIII-TUBULIN (ectoderm), αSMA (mesoderm) and AFP (endoderm) based on their fluorescence intensities (n = 3 independent experiments; means ± SEM; n.s., not significant; ANOVA). (C) Representative images of immunocytochemistry for the in vitro three-germ layer markers (βIII-TUBULIN, αSMA and AFP). Scale bars = 50 μm. (PDF 500 kb)
Additional file 2: Figure S2.
Somatic mutations in TiPSCs caused by the reprogramming process. a A substantial number of nonsynonymous variants were suggested to disappear in TKA4, all of which are located on chromosome 9, and Sanger sequencing showed an apparent mosaic loss of heterozygosity. An example of the Sanger sequence analysis of apparently lost nonsynonymous variant is shown. A small peak indicating a variant was observed in TKA4, suggesting a mosaic loss of heterozygosity. b The CGH analysis of TKA4 compared with the T-cells showed there were no copy number alterations in chromosome 9. However, a detailed inspection of the ratios of variants/reference bases revealed that the ratios were significantly different from 0.5 in the long arm of chromosome 9 of TKA4, supporting mosaic loss of heterozygosity in the region. Because the CGH analysis did not support copy number alterations (upper panel), these findings are more likely to represent somatic parental disomy involving the long arm of chromosome 9. A similar observation has recently been described . (PDF 580 kb)
Additional file 3: Table S1.
List of antibodies. (PDF 400 kb)
Additional file 4: Table S2.
List of primers. (PDF 397 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Fujimori, K., Tezuka, T., Ishiura, H. et al. Modeling neurological diseases with induced pluripotent cells reprogrammed from immortalized lymphoblastoid cell lines. Mol Brain 9, 88 (2016). https://doi.org/10.1186/s13041-016-0267-6
- Lymphoblastoid B-cell line
- Disease modeling
- Neurological disorder
- Induced pluripotent stem cells
- Genomic mutation in reprogramming process