- Open Access
The gene expression profiles of induced pluripotent stem cells from individuals with childhood cerebral adrenoleukodystrophy are consistent with proposed mechanisms of pathogenesis
Stem Cell Research & Therapy volume 3, Article number: 39 (2012)
X-linked adrenoleukodystrophy (X-ALD) is a complex disorder with variable expressivity that affects the nervous, adrenocortical and male reproductive systems. Although ABCD1 mutations are known to provide the genetic basis for X-ALD, its pathogenesis is not fully elucidated. While elevated very long chain fatty acid (VLCFA) levels in blood and reduced VLCFA catabolic activity in cultured fibroblasts are biomarkers used to identify ABCD1 mutation carriers, the roles peroxisomal lipid metabolism play in disease etiology are unknown.
Primary skin fibroblasts from two male patients with the childhood cerebral form of the disease (CCALD) caused by ABCD1 frameshift or missense mutations and three healthy donors were transduced with retroviral vectors expressing the OCT4, SOX2, KLF4 and c-MYC factors. Candidate induced pluripotent stem cells (iPSCs) were subject to global gene expression, DNA methylation, DNA copy number variation, and genotyping analysis and tested for pluripotency through in vitro differentiation and teratoma formation. Saturated VLCFA (sVLCFA) and plasmalogen levels in primary fibroblasts and iPSCs from healthy donors as well as CCALD patients were determined through mass spectroscopy.
Skin fibroblasts from CCALD patients and healthy donors were reprogrammed into validated iPSCs. Unlike fibroblasts, CCALD patient iPSCs show differentially expressed genes (DEGs) relevant to both peroxisome abundance and neuroinflammation. Also, in contrast to fibroblasts, iPSCs from patients showed no significant difference in sVLCFA levels relative to those from controls. In all cell types, the plasmalogen levels tested did not correlate with ABCD1 mutation status.
Normal ABCD1 gene function is not required for reprogramming skin fibroblasts into iPSCs or maintaining pluripotency. Relative to DEGs found in fibroblasts, DEGs uncovered in comparisons of CCALD patient and control iPSCs are more consistent with major hypotheses regarding disease pathogenesis. These DEGs were independent of differences in sVLCFA levels, which did not vary according to ABCD1 mutation status. The highlighted genes provide new leads for pathogenic mechanisms that can be explored in animal models and human tissue specimens. We suggest that these iPSC resources will have applications that include assisting efforts to identify genetic and environmental modifiers and screening for therapeutic interventions tailored towards affected cell populations and patient genotypes.
X-linked adrenoleukodystrophy (X-ALD) is a complex disorder caused by mutations in the ABCD1 gene that encodes an integral peroxisome membrane protein belonging to the ATP-binding cassette transporter superfamily [1–4]. X-ALD primarily affects the nervous system, adrenal cortex and testes with highly variable clinical presentations that are influenced by modifier genes and the environment [2, 3]. Males with ABCD1 mutations develop childhood cerebral ALD (CCALD) about 33% of the time and adult onset adrenomyeloneuropathy (AMN) about 45% of the time [2, 3]. CCALD patients typically show symptoms between five and nine years of age with rapid cerebral demyelination and adrenocortical atrophy. Within a few years of onset, they suffer dementia and progressive neurological deficits that eventually lead to death. In contrast, AMN patients show a later onset of disease (20 to 40 years of age) and present with adrenal insufficiency, a distal axonopathy in the spinal cord and peripheral neuropathy that results in progressive spastic paraparesis with debilitating end stage disease [2, 3]. Approximately 10% of hemizygotes develop primary adrenocortical insufficiency (Addison's disease) with no evidence of nervous system dysfunction [2, 3]. Disease prognosis is challenging since mutations do not correlate with clinical phenotypes  and male siblings with the same ABCD1 mutation, including monozygotic male twins [6, 7], can have dramatically different clinical presentations . Although hemizygotes typically show the most severe clinical manifestations of disease, about half of female ABCD1 mutation carriers develop AMN-like symptoms later in life [2, 3].
The molecular mechanisms underlying the inflammatory brain demyelination found in CCALD patients are not fully understood. It has been hypothesized to be related to the accumulation of saturated very long chain fatty acids (sVLCFAs) in specific central nervous system (CNS) cell types (for example, oligodendrocytes and microglial cells) and/or lipid classes (for example, ganglioside, phosphatidylcholine and cholesterol ester fractions) [9–11]. Other hypotheses have focused on the roles of oxidative stress [12–16], myelin sheath integrity , oligodendrocyte apoptosis and microglial cell activation [1, 2], and CNS cell membrane receptors .
Here, we report the generation and genomic characterization of CCALD patient-specific induced pluripotent stem cell (iPSC) model systems that can provide a platform to investigate cell autonomous processes relevant to X-ALD pathogenesis. The gene expression and biochemical profiles of these patient-specific iPSCs provide a novel perspective that supports the leading hypotheses regarding disease pathogenesis. Furthermore, our resources provide a first step required for the development and interpretation of patient-specific model systems that investigate genetic modifiers, environmental risk factors and non-cell autonomous processes relevant to the etiology of X-ALD.
Materials and methods
Cell culture conditions
Primary dermal fibroblast cultures from CCALD patients and controls were obtained from the Peroxisomal Disease Laboratory at the Kennedy Krieger Institute and Coriell Institute Cell Repositories, respectively. All cells described herein were cultured at 37°C with 5% CO2. Human primary dermal fibroblasts and mitomycin inactivated mouse embryonic fibroblasts (MEFs) were cultured in fibroblast media (DMEM supplemented with 10% FBS, L-glutamine, penicillin/streptomycin, vitamin solution, and essential and nonessential amino acids (Life technologies, Foster City, CA, USA)), as previously described . iPSCs were cultured on a layer of mitomycin C-inactivated MEF feeder cells in iPSC medium (DMEM:F12 medium supplemented with 20% KSR, L-glutamine, penicillin/streptomycin, nonessential amino acids, β-mecaptoethanol and bFGF (Life technologies, Foster City, CA, USA)) [20, 21].
Five different pMX retroviral vectors designed to deliver green fluorescent protein (GFP) and human OCT4, SOX2, KLF4 and C-MYC cDNA sequences were obtained from Addgene (Cambridge, MA, USA) [20, 21]. Primary human fibroblasts were twice transduced with a mixture of all five retroviruses as described [20, 21]. Transduction efficiency was evaluated by GFP expression. After four days, cells were re-plated onto MEF feeders and cultured in hESC medium containing 1 mM valproic acid . By four weeks, candidate iPSC colonies were manually picked and clonally expanded [20, 21]. A full list of the analyses conducted on each of the candidate iPSCs is described below and provided in Additional file 1.
Protein pluripotency biomarker analysis
Alkaline phosphatase staining was performed using the leukocyte alkaline phosphatase kit (Sigma-Aldrich, St. Louis, MO, USA). For immunostaining, cells were fixed in 4% paraformaldehyde for 20 minutes, permeabilized with 1% Triton X-100 for 5 minutes except for surface marker staining, and blocked in 1% BSA in 1 × PBS for 1 hour at room temperature. Primary antibody staining was performed at 4°C overnight with antibodies against OCT4 and NANOG (goat polyclonal IgG against human, R&D Systems, Minneapolis, MN, USA), SOX2 (goat polyclonal IgG against mouse, rat and human) and SSEA4 (mouse monoclonal IgG3 against human, Santa Cruz Biotechnology, Santa Cruz, CA, USA), TRA-1-60 (mouse monoclonal IgM against human, Millipore, Billerica, MA, USA), TuJ1 (rabbit monoclonal IgG1 against mammalian, Covance Research Products, Princeton, NJ, USA), α-SMA (mouse monoclonal IgG2a against human, mouse and so on, Sigma-Aldrich, St. Louis, MO, USA), and AFP (mouse monoclonal IgG2a against human, pig and canine, Life technologies, Foster City, CA, USA). Secondary antibody staining was performed at room temperature for one hour with appropriate fluorescence conjugated secondary antibodies from Life technologies, Foster City, CA, USA and Jackson ImmunoResearch, West Grove, PA, USA. Nuclei were visualized by staining with 100 ng/ml DAPI (Life technologies, Foster City, CA, USA).
Gene expression profiling
Total RNA samples (100 ng per sample) were converted into biotin-labeled cRNA targets (Affymetrix GeneChip® IVT Labeling Kit, Santa Clara, CA, USA), processed and analyzed on Affymetrix Human Genome 133A 2.0 or 133 Plus 2.0 GeneChips, as previously described . Using WebArray software (San Diego, CA, USA), we applied the RMA algorithm to generate log2-transformed gene expression values and linear model statistical analysis (limma) to identify differentially expressed genes (DEGs) with false discovery rates (FDRs) calculated using the spacings LOESS histogram (SPLOSH) method [23, 24] (Additional file 2). We performed hierarchical clustering analysis using Partek Genomics Suite software (St. Louis, MO, USA). We conducted GeneOntology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses using WebGestalt software (Nashville, TN, USA) [25–27]. We used the DAVID (Database for Annotation, Visualization and Integrated Discovery) v6.7 bioinformatics resource  for the annotation of gene functions. Scaled gene expression scores and .cel files are available at the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) repository under Series Accession Number GSE34308.
DNA methylation profiling
Genomic DNA was extracted from cultured cells as described [29, 30] and analyzed on 450 K Infinium Methylation BeadChips (Illumina, San Diego, CA, USA), which interrogate the methylation status of over 485,000 CpG sites distributed across the human genome. The resulting data were analyzed using GenomeStudio software (that is, unmethylated to fully methylated) for each locus (Additional file 3). Bisulfite DNA sequencing was conducted as previously described [29, 30].
Genotyping and copy number analysis
Genomic DNA from cultured cells was also analyzed using Human CytoSNP-12 Infinium HD BeadChips (Illumina) that interrogate the genotypes of 299,140 human single nucleotide polymorphisms (SNPs). Data filtering and analysis were performed in GenomeStudio (Illumina). Copy number analysis was performed using CNVPartition version 2.4.4 with a confidence threshold set at 50 and a minimum of 10 SNP probes per CNV region, as previously described . In multiple samples, we performed the global genotyping analysis two independent times and only assigned a copy number change (CNC) if both analyses were in agreement (Additional file 4). Dideoxysequencing of ABCD1 exons 1, 8 and 9 was performed as previously described [32, 33].
In vitro differentiation and teratoma assays
iPSCs were detached from culture dishes with collagenase IV, maintained in suspension to induce embryoid body (EB) formation and subjected to an in vitro differentiation procedure, as described [34, 35]. For teratoma analysis, iPSCs from a confluent 10 cm2 plate were harvested and subcutaneously injected to the dorsal flanks of immunodeficient (SCID) mice (Jackson Laboratory, Bar Harbor, ME, USA), as described . Nine weeks after injection, teratomas were excised, fixed in 10% formalin, sectioned and stained with hematoxylin and eosin.
We used liquid chromatography-tandem mass spectrometry (LC-MS/MS) to measure C26:0-lysophosphorylcholine and phosphatidylethanolamine (PE) plasmalogen levels in cell lysates processed by methanol extraction as described in reference . Herein, C26:0-lysophosphorylcholine measurements were used to evaluate VLCFA levels. The tetradeuterated analog of 1-O-hexadecyl-2-lysn-sn-3 phosphorylcholine was used to quantify PE plasmalogens. PE plasmalogens were identified based on the fragmentation patterns reported in reference .
Derivation of candidate iPSCs from CCALD patient fibroblasts
Primary skin fibroblast cultures from three healthy donors and two CCALD patients (Table 1) were infected with retroviruses designed to express the human OCT4, SOX2, KLF4 and c-MYC genes (Methods). We observed iPSC-like colonies for two weeks and clonally expanded TRA-1-60 positive colonies for four weeks, consistent with prior reports of reprogramming skin fibroblasts from healthy human donors . All candidate iPSC colonies maintained the expected morphological features and expressed protein biomarkers of pluripotency (Figure 1A, Additional file 1).
Genotypes and DNA copy number profiles of iPSCs
We confirmed that the patient iPSCs had the expected mutant ABCD1 genotypes and that control iPSCs lacked these specific ABCD1 mutations by dideoxysequencing. As determined by BeadArray analysis, the genotypes of over 290,000 SNPs in iPSCs and original fibroblasts were > 99.9% concordant. Based on the same genotyping data, we did not detect copy number changes (CNCs, that is, insertions or deletions greater than 10-kb in length) in patient CCALD1-1, CCALD1-2 and CCALD2-1 iPSCs or Control1-3, Control1-4 and Control2-1 iPSCs (Additional file 4). Consistent with prior reports of reprogrammed human cells [39–41], we detected CNCs in 8/14 (57%) iPSCs analyzed (Additional file 4). These iPSCs had one (CCALD1-3, CCALD1-4 and Control2-3), two (Control1-2 and Control3-1), three (Control2-2 and Control2-4) or five (Control1-1) separate genomic regions affected by a CNC (Additional file 4).
Gene expression profiles of CCALD and control donor cells
We validated the robust expression of previously reported iPSC signature genes in control and CCALD donor-derived iPSCs and skin fibroblasts  based on a subset of the data generated from global expression profiling of over 18,000 transcripts (Figure 2A). Unsupervised hierarchical clustering analysis based on the expression of preselected pluripotency biomarkers or the most variable transcripts (that is, coefficient of variation (CV) > 0.10 across all samples)) produced two distinct clusters consisting of skin fibroblasts and the iPSCs (Figure 1 and Additional file 5).
DNA methylation profiles of CCALD and control donor cells
We performed global DNA methylation analysis interrogating over 485,000 CpG sites of all starting fibroblast cultures and reprogrammed iPSCs (Table 1) (Methods). Hierarchical clustering analysis demonstrated that the iPSCs and fibroblasts have distinct DNA methylation profiles that were independent of ABCD1 mutation status (Figure 2B).
In vitro differentiation of CCALD patient iPSCs and in vivo teratoma formation
All control and CCALD patient-derived iPSCs subjected to in vitro embryoid body assays produced cells derived from all three germ layers (Methods) (Figure 1B). Furthermore, all three iPSCs tested (one control and two CCALD patient donors) formed teratomas when injected into immune-deficient mice (Methods). In all these teratomas, histological analysis demonstrated the presence of tissues representative of all three germ layers (Figure 1C).
Comparative lipid profiling of cells from healthy donors and those with CCALD
Consistent with prior reports , cultured CCALD patient skin fibroblasts grown in fibroblast growth media showed 4.3-fold elevated VLCFA levels (based on C26:0-lysophosphorylcholine measurements, see Materials and methods), but similar PE plasmalogen levels relative to fibroblasts from healthy donors (Figure 3A, C)). Similarly, we found 3.8-fold elevated VLCFA levels, but comparable PE plasmalogen levels, in patients relative to control fibroblasts grown under the same conditions in iPSC media (Figure 3A, C). In contrast, no significant differences were found for either VLCFA or PE plasmalogen levels in patient and control iPSCs cultured under the same conditions in iPSC media (Figure 3). Nevertheless, all fibroblast cultures had approximately two-fold elevated PE plasmalogen levels when grown in iPSC relative to fibroblast growth media. The fact that iPSCs cannot be maintained in the undifferentiated state in fibroblast media precluded lipid analyses under these conditions.
Differential gene expression among patient and control cells
Among the two patient and three control fibroblasts used for reprogramming, no differentially expressed genes (DEGs) ( > 1.2-fold change, FDR < 0.1) were found. This was not unexpected due to the limited number of samples analyzed. To enhance our ability to detect potential DEGs, we conducted a larger-scale gene expression analysis ( > 47,000 transcripts) using cultured skin fibroblasts from five healthy control donors and five CCALD patients. In these studies, we identified 127 DEGs (44 with lower and 83 with higher expression in patient relative to control fibroblasts) (Additional File 6).
Based on GeneOntology (GO) analysis, we found a total of 13 functional categories enriched (≥ 4 genes, B-H corrected P < 0.05) for DEGs with higher expression in patient relative to control fibroblasts (Additional file 7). Of these, the most specific GO category was nuclear lumen genes. KEGG analysis did not show any enriched categories (≥ 4 genes, B-H corrected P < 0.05) for DEGs with higher expression in patient relative to control fibroblasts. In contrast, DEGs with lower expression in patient relative to healthy control fibroblasts were enriched for one KEGG category ('Pathways in Cancer' as represented by the SMAD3, LAMA4, PML and DAPK1 genes), but no GO categories.
Given the possible gaps in public databases of gene functions relevant to peroxisome biology and X-ALD pathogenesis, we used the DAVID Bioinformatics resource to annotate the function of DEGs and manually searched for genes relevant to peroxisome biology, lipid metabolism, oxidative stress and neuroinflammation (Additional file 8). The only peroxisomal gene was AGPS, which is involved in plasmalogen (a type of ether phospholipid) biosynthesis, and it had higher expression in patient fibroblasts. Nevertheless, we note that the PE plasmalogen levels in patient and healthy control fibroblasts were similar under both growth conditions (Figure 3C). Two gene encoding enzymes involved in the Lands cycle, deacylation/reacylation reactions responsible for glycerophospholipid remodeling , were present with the higher expression of LYPLA1 (deacylation), but lower expression of MBOAT7 (reacylation), in patients relative to control fibroblasts. Two DEGs were involved in sphingosine metabolism with SMPD1, which converts sphingomyelin to ceramide, showing lower expression in patients relative to control fibroblasts and SGPL1, which degrades sphingosine-1-phosphate, showing higher expression in patients relative to control fibroblasts. This is of interest given the differences in the sphingolipid composition of white matter from the brains of CCALD and healthy individuals . No genes involved in classic oxidative stress responses were found in the list. The DEGs strongly or primarily related to neuroinflammation included CBLB, RAB27A and TCF3, which showed higher expression in patient fibroblasts, as well as CD81, MASP1 and GRK5, which showed higher expression in control fibroblasts.
Based on our manual curation, we found that the identity of approximately 40% of the DEGs (46 total) were consistent with the expression profiles of cultured fibroblasts related to the site of skin biopsy (Additional file 6). All these genes showed the highest variability in expression based on biopsy sites, as described in reference . We also note that the expression profiles of 4/6 DEGs described above as being involved in neuroinflammation, (RAB27A, CD81, MASP1 and GRK5), are also influenced by the biopsy site. Although all the fibroblasts in our study were obtained from the upper limbs, the control and patient donor cells were collected and expanded at different laboratories (Methods), which could influence their gene expression signatures.
We identified 75 DEGs (22 with higher and 53 with lower expression in patient relative to control cells) based on the gene expression profiles of five CCALD iPSCs from two CCALD donors and nine control iPSCs from three healthy donors (Additional file 9). There was no overlap with the Affymetrix probe IDs of the DEGs uncovered in the cultured skin fibroblasts from the five healthy controls and five CCALD patient donors discussed above. Different Affymetrix probe IDs interrogated the CEP57 gene indicated it was a DEG in both systems, but in opposing directions (that is, higher in patient relative to control fibroblasts, but higher in control relative to patient iPSCs). Based on GO analysis, we found a total of 14 functional categories enriched (≥ 4 genes, B-H corrected P < 0.05) for DEGs with higher expression in patient relative to control cells (Additional file 10). These included blood vessel morphogenesis, regulation of cellular protein metabolic process and carboxylic acid metabolic process. In contrast, GO analysis identified no enriched categories for DEGs with higher expression in healthy control cells. KEGG analysis did not identify any enriched pathways (≥ 4 genes, B-H corrected P < 0.05) for DEGs with higher expression in either the patient or control cells.
Although GO and KEGG analysis did not highlight biological processes proposed to be relevant to disease, inspection of the DEG functions based on the DAVID Bioinformatics resource (Additional file 11) uncovered genes associated with major hypotheses pertinent to X-ALD pathogenesis. Among the relevant genes with reduced expression in CCALD patient relative to healthy donor-derived iPSCs were PEX11B (confirmed by qPCR in Figure 4) and CD200. The former plays a pivotal role in peroxisome proliferation and maintenance [46, 47]. Decreased CD200 expression is associated with the activation and accumulation of macrophages, including brain microglia, and causes inflammatory responses in other systems [48, 49]. DEGs with higher expression in patient relative to control iPSCs were also related to hypotheses relevant to X-ALD pathogenesis and lipid metabolism. ULK1 is the mammalian homolog of the yeast Atg1 gene, which plays a critical role in the autophagy-mediated turnover of peroxisomes in yeast . PLA2G2A is involved in phospholipid turnover . NAAA, THBS1 and BSG all have functions related to neuroinflammation [52–61]. SLC7A8 is a transporter of thyroid hormones, which can induce peroxisomal biogenesis and β-oxidation as well as the ABCD2 expression, whose induction can correct biochemical functions of X-ALD patient fibroblasts .
Robust differences in DNA methylation frequently found between fibroblasts and iPSCs are not associated with ABCD1 mutation status
In our global DNA methylation analysis (Methods), the starting 5 fibroblasts and 14 iPSCs showed over 62,000 loci where there was a 0.25 unit difference in average β-values and B-H corrected P < 0.05 (Additional file 12). To focus on the most robust differentially methylated loci (DML), we identified 744 sites that were hypomethylated (β < 0.2) in all samples of one group and hypermethylated (β > 0.8) in all samples in the remaining group (fibroblasts and iPSCs). We identified 266 distinct genes proximal to at least one DML with approximately half of them (49%) being poorly expressed ( < 60th percentile) in all samples regardless of its methylation status, consistent with our prior results involving the co-analysis of DNA methylation and gene expression profiles in normal and cancerous B-cell populations [29, 30]. Nevertheless, 60 genes (23%) were at least modestly higher expressed when the proximal CpG locus was in the hypomethylated relative to hypermethylated state ( > 1.5-fold change, FDR < 0.1 and expressed > 60th percentile in the hypomethylated state).
In contrast, we only found five robust DMLs in comparisons of two patient and three control fibroblasts and one robust DML in comparisons of five patient and nine control iPSCs (Additional file 13). The DML present in the iPSC analysis was also present in the fibroblast analysis, with higher methylation being found in all CCALD patients relative to control donor cells regardless of ABCD1 mutation status. This shared DML was proximal to the PRDM15 gene, whose expression was not interrogated in our global GeneChip gene expression assays. The remaining four DMLs in fibroblasts were proximal to the PAX3, CCDC140, UTRN and BAIAP2 genes. All three of the genes interrogated by our GeneChip expression assays (PAX3, UTRN and BAIAP) were poorly expressed in all fibroblasts regardless of ABCD1 mutation status.
Local gene expression is not substantially affected by CNCs found in iPSCs
To begin to address the influence that CNCs present in iPSC have on their transcriptome, we focused on the expression profiles of genes residing in the affected genomic regions. A total of 11 amplified segments containing 22 unique genes were found in 8 iPSCs. Only six of these unique genes showed elevated expression in the amplified relative to the diploid samples (that is, a > 1.2-fold increase relative to the mean expression scores of diploid samples and falling outside their range). This included the ID1 gene in CCALD1-3, WWC1 gene in CCALD1-4, and IQCA1, CXCR7, SQLE and KIAA0196 genes in Control1-1.
Three iPSCs (Control2-2, Control2-4 and Control3-1) showed evidence of having at least one genomic deletion, with evidence in each case that one allele was retained. Collectively, five unique genes were present in the four deleted genomic regions in these iPSCs. There was no evidence of reduced expression in the samples with reduced copy number (that is, > 1.2-fold decrease relative to the mean expression scores of diploid samples and falling outside their range).
Amplified or deleted segments show no differences in DNA methylation status
A total of 745 DNA methylation assays interrogated loci located within amplified regions present in control or patient iPSCs. In all cases, the DNA methylation status of such genomic regions was similar regardless of whether it was in the diploid or amplified state. In fact, we observed no evidence of a block of DNA methylation change associated with a CNC (that is, three or more contiguous assays wherein the β-value of the amplified segment was greater than 0.2 units outside the range of the diploid samples).
Next, we accessed the methylation status of genomic regions subject to a loss of copy number in iPSCs. A total of 79 DNA methylation assays interrogate loci with the genomic regions of heterozygous deletion. The affected samples included Control2-iPS2 (chr13:q14.2), Control2-iPS4 (chr2:q33.3) and Control3-iPS1 (chr3:p14.2 and chr5:p15.2). Again, we observed no evidence of a block of DNA methylation change associated with a CNC (that is, three or more contiguous assays wherein the β-value of the deleted segment fell more than 0.2 units outside the range of the diploid samples).
X-ALD is a complex peroxisomal disorder with variable expressivity. Although its primary genetic basis has been known for some time [1–4], the exact nature of X-ALD pathogenesis and its genetic and environmental modifiers have not been elucidated. Here, we generated iPSC resources for the longer-term purpose of developing novel tissue culture models for elucidating the pathogenesis of X-ALD and screening for more effective drug therapies. In keeping with prior reports , skin fibroblasts from ABCD1 mutation carriers can be reprogrammed to form iPSCs with the hallmark molecular properties of pluripotency, including the expression of appropriate gene and protein biomarkers and changes in DNA methylation levels, as summarized in Additional file 1. Patient iPSCs could be differentiated into embryoid bodies and differentiated in vitro into representative cell types of all three germ layers. Most importantly, patient iPSCs formed teratomas with evidence of cell types from all three germ layers.
Consistent with prior reports [64–66], we identified de novo CNCs over 10-kb in length in approximately half of our iPSCs. These CNCs were only found in healthy donor iPSCs and not those generated from patient cells (Additional file 4). The observed genomic deletions always provided evidence of affecting only one allele and genomic amplifications always involved a limited increase in copy number (three to four copies maximum). Due to the fact that we conducted global expression and DNA methylation analyses on these samples, we could investigate the effects that these CNCs have on the expression of genes located within affected genomic segments. In almost all circumstances, their expression levels were within the range of diploid samples. Although multiple factors likely contribute to these observations, we favor the explanation that this primarily reflects the effects of selection whereby CNCs are only tolerated in iPSCs if they involve genomic regions that do not influence the initiation of reprogramming or maintenance of pluripotency.
As a result of our genomic characterization of these cell resources, we acquired global gene expression data from patient and control fibroblasts. Many DEGs were previously reported to be associated with the site of biopsy . This is reasonable given that the patient and control fibroblasts were acquired from different institutions even though all biopsies involved the upper limbs of donors. We sought to determine if there was enrichment for functional categories or biological processes in the DEGs, keeping in mind the limitations of using cultured cells to study complex diseases involving interactions between multiple organ systems. Only very broad functional categories or KEGG pathways were highlighted in these analyses, with none of them showing a direct relation to disease.
Since there are likely to be gaps in public databases of processes relevant to peroxisome biology and X-ALD pathogenesis, we conducted a manual inspection of gene annotations provided by the DAVID bioinformatics resource and found multiple DEGs involved in immune related processes, but only two (CBLB and RAB27A) of these genes were not associated with the site of biopsy. CBLB plays a critical role in antigen-induced immune tolerance and Cblb-deficient mice immunized with myelin basic protein are more susceptible to experimental autoimmune encephalomyelitis (EAE), a model for multiple sclerosis [67, 68]. RAB27A mutations can lead to an uncontrolled T lymphocyte and macrophage activation syndrome in humans, with some individuals showing possible leukocyte brain infiltration . In one Saudi Arabian kindred, RAB27 mutations were associated with immunodeficiency and progressive demyelination of brain white matter .
The DEGs found in patient and control iPSCs did not overlap with those found in fibroblasts and instead were consistent with several leading hypotheses regarding X-ALD pathogenesis. This suggests that the reprogramming process can minimize the confounding influence the site of skin biopsy has on the gene expression profiles of cultured fibroblasts. In particular, we highlight the reduced expression of PEX11B, a major controller of peroxisome proliferation and neuroinflammatory genes, in patient relative to control iPSCs. Pex11B null mice show several pathologic features, including neuronal migration defects, enhanced neuronal apoptosis, developmental delay, hypotonia and neonatal lethality . Despite these severe phenotypes, Pex11B null mice displays only mild defects in peroxisomal fatty acid beta-oxidation and ether lipid biosynthesis . Intriguingly, the deletion of a single Pex11B allele leads to a slightly increased number of peroxisomes, increased levels of oxidative stress in brain tissue, and neuronal cell death in mice . In addition, the ULK1, whose yeast homolog plays a critical role in the autophagy-mediated peroxisome turnover , showed higher expression in CCALD patient relative to control iPSCs. It is tempting to speculate that if nervous and immune system cells of ABCD1 mutation carriers indeed have aberrant PEX11B and ULK1 activity, this could lead to differences in peroxisome abundance and/or activity that increase reactive oxygen species (ROS) levels and promote X-ALD pathogenesis. It remains to be determined if ABCD1 mutation carriers have abnormal peroxisome abundance in their pertinent nervous and immune system cells and tissues.
In a similar vein, the increased NAAA, THBS1, BSG (aka CD147 aka EMMPRIN) and NOTCH1 gene expression in patients relative to control iPSCs is supportive of hypotheses regarding a predisposition to neuroinflammation that is a prelude to devastating autoimmune responses. NAAA hydrolyzes palmitoylethanolamide (PEA), a naturally occurring lipid amide that, when administered as a drug, inhibits inflammatory responses . In principle, increasing leukocyte NAAA levels could reduce PEA levels and promote inflammation. In fact, a chemical inhibitor of NAAA function attenuates inflammation and tissue damage and improves recovery of motor function in mice with spinal cord trauma . Intriguingly, CD200 has been proposed to play a role in the immune privileged status of the CNS when CD200-mediated immune suppression occurs via neuron-microglial as well as glial-glial interactions in inflammatory conditions . THBS1 is linked to neuroinflammatory processes involving astrocyte and microglia through its role in processing and activating the TGF-β ligand  and is also implicated in responses to oxidative stress . Likewise, Notch1 is involved in microglial associated inflammation . Also of relevance are emerging reports that BSG acts a master regulator of matrix metalloproteinases (MMP) implicated in most diseases involving neuroinflammation and thus has been proposed to play a role in the immune-privileged status of the CNS .
Although we highlight the possible implications of the gene expression profiles observed in patient iPSCs, we note alternative hypotheses regarding their origins and biological significance. While the iPSCs described in this study have the hallmark properties of pluripotency, their gene expression profiles could reflect subtle ABCD1 mutation status-dependent differences in their predisposition to differentiate into specific cell types and lineages. Comparisons of the gene expression profiles of mature cell types derived from patient- and healthy donor-derived iPSCs will be especially informative. The persistence or elimination of groups of DEGs reflective of biological processes and pathways could provide a means of assessing the tissue specificity of disease and enhance the ability to discern biologically informative gene expression signatures from noise resulting from confounding variables, such as tissue culture conditions.
Although ABCD1 mutation carriers show elevated sVLCFA levels in their blood and urine and reduced sVLCFA catabolic activity in their cultured fibroblasts, the role of sVLCFA in disease pathogenesis is still under discussion. The significance of decreased plasmalogen levels in the patients' brain white matter also is unclear [44, 78]. As expected, CCALD patient fibroblasts had elevated VLCFA levels, but similar PE plasmalogen levels, relative to those from healthy donors (Figure 3A, C). Likewise, iPSCs from CCALD patient and healthy control donors also showed similar PE plasmalogen levels (Figure 3D). The fact that all patient and control iPSCs tested had low VLCFA levels, based on C26:0-lysophosphorylcholine measurements, (Figure 3B) is puzzling, yet consistent with prior reports . VLCFA levels are determined by their rate of synthesis (that is, catalyzed by fatty acid elongating enzyme ELOVL1 ), degradation (affected by ABCD1 and ABCD2 function ) and uptake of these fatty acids from the culture medium. As such, one hypothesis is that the rate of VLCFA synthesis is lower in iPSCs relative to fibroblasts under the culture conditions evaluated. Consistent with this hypothesis and prior reports , the expression of the ELOVL1 gene was significantly lower (2.6-fold, FDR < 2.2 × 10-9) in iPSCs relative to fibroblasts. Nevertheless, we note that ELOVL1 was not differentially expressed in patient relative to control fibroblasts or iPSCs. An alternate hypothesis that the ABCD2 gene is compensating for the impaired ABCD1 function in patient iPSCs; however, ABCD2 was not differentially expressed in patient relative to control fibroblasts or iPSCs. This does not preclude the possibilities that ABCD2 activity is being increased on the protein level or that another gene is playing a major role in significantly lowering VLCFA levels in CCALD iPSCs. We also note a prior hypothesis that the rapid growth rate of iPSCs could reduce their VLCFA levels, independent of their ABCD1 mutation status . Fibroblasts have altered morphology and slowed growth in iPSC media relative to fibroblast media, which according to the growth rate hypothesis could contribute to their reduced VLCFA levels. Given that iPSCs can rapidly differentiate in fibroblast media, iPSC growth media provides an imperfect, but necessary, compromise in direct comparisons between cultured fibroblasts and iPSCs. We note the potential contribution of MEF feeder cells to iPSC lipid profiles and the advantages of using feeder-free media in future experiments.
Future applications and directions
The impending implementation of newborn screening for X-ALD based on blood lipid profiles will increase the demand for model systems to screen for more effective therapeutic interventions [36, 81]. Early detection would provide physicians with a window of opportunity to treat presymptomatic patients prior to the development of CCALD, and may also prevent or delay AMN onset. Therapeutic interventions, such as Lorenzo's Oil, help prevent the onset of cerebral disease in some individuals, but are not effective for the majority of CCALD patients and, likewise, there are no effective options for AMN [75–77]. A compelling attribute of iPSC model systems is that they represent the exact ABCD1 mutations found in the patient population and thus provide an opportunity to test therapeutic agents tailored to a patient's genotype in cell populations most affected by disease. Examples of genotype-dependent therapeutic strategies include nonsense suppressor drugs  and molecular chaperones  for individuals with nonsense and missense mutations, respectively.
The fact that CCALD iPSCs show gene expression profiles similar to those derived from healthy controls may reflect the fact that X-ALD clinical symptoms do not manifest at birth but, instead, occur in early childhood or later in life. Given that ABCD1 mutant mice show clinical aspects of X-ALD with increasing age , it is possible that later passage CCALD iPSCs and their derivatives may manifest gene expression profiles and/or functional properties more consistent with disease pathogenesis and progression. In this regard, a comparison of the properties of iPSCs and their derivatives previously obtained from other CCALD and AMN patients  as a function of in vitro passage number could be informative. Despite the promise of iPSC approaches, it will remain a significant challenge to generate and optimize in vitro model systems for X-ALD and other complex disorders that involve multiple organ systems as well as unknown gene-environment interactions and genetic modifiers.
We have reprogrammed skin fibroblasts from CCALD patients and control donor primary fibroblasts into iPSCs that show all the fundamental hallmark molecular and cellular properties of pluripotency. The DEGs found in comparisons of patient- and healthy donor-derived iPSCs are consistent with emerging hypotheses regarding the role of peroxisomes, oxidative stress and neuroinflammation in the pathogenesis of X-ALD. Given that the sVLFCA profiles of the patient iPSCs are similar to those of control iPSCs, these gene expression profiles are not dependent on the presence of the hallmark lipid biomarkers of disease. Additional gene expression and functional analyses involving differentiated cell types derived from CCALD and control iPSCs could be especially informative given our preliminary results. This would include cell types related to the CNS (for example, oligodendrocytes, astrocytes and microglia), adrenocortical and male reproductive (for example, Leydig and Sertoli cells) aspects of disease. Furthermore, investigations involving patient tissue samples and animal models are required in order to determine if the observed fibroblast and iPSC gene expression profiles are reflective of pathogenic mechanisms or are simply specific to our cultured cells.
We thank P. Watkins, G. Raymond, N. Braverman and R. Pelikan for thoughtful discussion and/or assistance in data analysis. We thank D. Weisenberger and D. Van Den Berg at the USC Epigenome Center for conducting the Illumina BeadArray DNA methylation and SNP genotyping assays and for advice in data analysis. This study was funded by the National Institutes of Health (GM072447 and GM072447-S1) and a California Institute of Regenerative Medicine Pre-doctoral Training Grant (WYY).
childhood cerebral adrenoleukodystrophy
copy number change
central nervous system
coefficient of variation
Database for Annotation: Visualization and Integrated Discovery
differentially expressed genes
Dulbecco's Modified Eagle Medium
differentially methylated loci
experimental autoimmune encephalomyelitis
fetal bovine serum
false discovery rates
green fluorescent protein
induced pluripotent stem cell
Kyoto Encyclopedia of Genes and Genomes
liquid chromatography-tandem mass spectrometry
mouse embryonic fibroblasts
National Center for Biotechnology Information
reactive oxygen species
single nucleotide polymorphisms
saturated very long chain fatty acid(s)
Kemp S, Wanders R: Biochemical aspects of X-linked adrenoleukodystrophy. Brain Pathol. 2010, 20: 831-837.
Berger J, Pujol A, Aubourg P, Forss-Petter S: Current and future pharmacological treatment strategies in X-linked adrenoleukodystrophy. Brain Pathol. 2010, 20: 845-856. 10.1111/j.1750-3639.2010.00393.x.
Berger J, Gartner J: X-linked adrenoleukodystrophy: clinical, biochemical and pathogenetic aspects. Biochim Biophys Acta. 2006, 1763: 1721-1732. 10.1016/j.bbamcr.2006.07.010.
Moser HW, Mahmood A, Raymond GV: X-linked adrenoleukodystrophy. Nat Clin Pract Neurol. 2007, 3: 140-151.
Kok F, Neumann S, Sarde CO, Zheng S, Wu KH, Wei HM, Bergin J, Watkins PA, Gould S, Sack G, Moser H, Mandel JL, Smith KD: Mutational analysis of patients with X-linked adrenoleukodystrophy. Hum Mutat. 1995, 6: 104-115. 10.1002/humu.1380060203.
Korenke GC, Fuchs S, Krasemann E, Doerr HG, Wilichowski E, Hunneman DH, Hanefeld F: Cerebral adrenoleukodystrophy (ALD) in only one of monozygotic twins with an identical ALD genotype. Ann Neurol. 1996, 40: 254-257. 10.1002/ana.410400221.
Wilichowski E, Ohlenbusch A, Korenke GC, Hunneman DH, Hanefeld F: Identical mitochondrial DNA in monozygotic twins with discordant adrenoleukodystrophy phenotype. Ann Neurol. 1998, 43: 835-836.
Mosser J, Douar AM, Sarde CO, Kioschis P, Feil R, Moser H, Poustka AM, Mandel JL, Aubourg P: Putative X-linked adrenoleukodystrophy gene shares unexpected homology with ABC transporters. Nature. 1993, 361: 726-730. 10.1038/361726a0.
Theda C, Moser AB, Powers JM, Moser HW: Phospholipids in X-linked adrenoleukodystrophy white matter: fatty acid abnormalities before the onset of demyelination. J Neurol Sci. 1992, 110: 195-204. 10.1016/0022-510X(92)90028-J.
Garashi M, Belchis D, Suzuki K: Brain gangliosides in adrenoleukodystrophy. J Neurochem. 1976, 27: 327-328. 10.1111/j.1471-4159.1976.tb01593.x.
Asheuer M, Bieche I, Laurendeau I, Moser A, Hainque B, Vidaud M, Aubourg P: Decreased expression of ABCD4 and BG1 genes early in the pathogenesis of X-linked adrenoleukodystrophy. Hum Mol Genet. 2005, 14: 1293-1303. 10.1093/hmg/ddi140.
Singh I, Pujol A: Pathomechanisms underlying X-adrenoleukodystrophy: a three-hit hypothesis. Brain Pathol. 2010, 20: 838-844.
Vargas CR, Wajner M, Sirtori LR, Goulart L, Chiochetta M, Coelho D, Latini A, Llesuy S, Bello-Klein A, Giugliani R, Deon M, Mello CF: Evidence that oxidative stress is increased in patients with X-linked adrenoleukodystrophy. Biochim Biophys Acta. 2004, 1688: 26-32. 10.1016/j.bbadis.2003.10.004.
Powers JM, Pei Z, Heinzer AK, Deering R, Moser AB, Moser HW, Watkins PA, Smith KD: Adreno-leukodystrophy: oxidative stress of mice and men. J Neuropathol Exp Neurol. 2005, 64: 1067-1079. 10.1097/01.jnen.0000190064.28559.a4.
Gilg AG, Singh AK, Singh I: Inducible nitric oxide synthase in the central nervous system of patients with X-adrenoleukodystrophy. J Neuropathol Exp Neurol. 2000, 59: 1063-1069.
Fourcade S, Lopez-Erauskin J, Galino J, Duval C, Naudi A, Jove M, Kemp S, Villarroya F, Ferrer I, Pamplona R, Portero-Otin M, Pujol A: Early oxidative damage underlying neurodegeneration in X-adrenoleukodystrophy. Hum Mol Genet. 2008, 17: 1762-1773. 10.1093/hmg/ddn085.
Ho JK, Moser H, Kishimoto Y, Hamilton JA: Interactions of a very long chain fatty acid with model membranes and serum albumin. Implications for the pathogenesis of adrenoleukodystrophy. J Clin Invest. 1995, 96: 1455-1463. 10.1172/JCI118182.
Singh H, Brogan M, Johnson D, Poulos A: Peroxisomal beta-oxidation of branched chain fatty acids in human skin fibroblasts. J Lipid Res. 1992, 33: 1597-1605.
Karaman MW, Houck ML, Chemnick LG, Nagpal S, Chawannakul D, Sudano D, Pike BL, Ho VV, Ryder OA, Hacia JG: Comparative analysis of gene-expression patterns in human and african great ape cultured fibroblasts. Genome Res. 2003, 13: 1619-1630. 10.1101/gr.1289803.
Park IH, Lerou PH, Zhao R, Huo H, Daley GQ: Generation of human-induced pluripotent stem cells. Nat Protoc. 2008, 3: 1180-1186. 10.1038/nprot.2008.92.
Takahashi K, Okita K, Nakagawa M, Yamanaka S: Induction of pluripotent stem cells from fibroblast cultures. Nat Protoc. 2007, 2: 3081-3089. 10.1038/nprot.2007.418.
Huangfu D, Maehr R, Guo W, Eijkelenboom A, Snitow M, Chen AE, Melton DA: Induction of pluripotent stem cells by defined factors is greatly improved by small-molecule compounds. Nat Biotechnol. 2008, 26: 795-797. 10.1038/nbt1418.
Xia X, McClelland M, Wang Y: WebArray: an online platform for microarray data analysis. BMC Bioinformatics. 2005, 6: 306-10.1186/1471-2105-6-306.
Wang Y, McClelland M, Xia XQ: Analyzing microarray data using WebArray. Cold Spring Harb Protoc. 2009, 8: 5260-
Zhang B, Kirov S, Snoddy J: WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 2005, 33: W741-748. 10.1093/nar/gki475.
Tazi KA, Quioc JJ, Abdel-Razek W, Tellier Z, Guichard C, Ogier-Denis E, Lebrec D, Moreau R: Protein array technology to investigate cytokine production by monocytes from patients with advanced alcoholic cirrhosis: an ex vivo pilot study. Hepatol Res. 2009, 39: 706-715. 10.1111/j.1872-034X.2009.00498.x.
Kirov SA, Zhang B, Snoddy JR: Association analysis for large-scale gene set data. Methods Mol Biol. 2007, 408: 19-33. 10.1007/978-1-59745-547-3_2.
Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.
Pike BL, Greiner TC, Wang X, Weisenburger DD, Hsu YH, Renaud G, Wolfsberg TG, Kim M, Weisenberger DJ, Siegmund KD, Ye W, Groshen S, Mehrian-Shai R, Delabie J, Chan WC, Laird PW, Hacia JG: DNA methylation profiles in diffuse large B-cell lymphoma and their relationship to gene expression status. Leukemia. 2008, 22: 1035-1043. 10.1038/leu.2008.18.
Wang XM, Greiner TC, Bibikova M, Pike BL, Siegmund KD, Sinha UK, Muschen M, Jaeger EB, Weisenburger DD, Chan WC, Shibata D, Fan JB, Hacia JG: Identification and functional relevance of de novo DNA methylation in cancerous B-cell populations. J Cell Biochem. 2010, 109: 818-827.
Laurent LC, Ulitsky I, Slavin I, Tran H, Schork A, Morey R, Lynch C, Harness JV, Lee S, Barrero MJ, Ku S, Martynova M, Semechkin R, Galat V, Gottesfeld J, Izpisua Belmonte JC, Murry C, Keirstead HS, Park HS, Schmidt U, Laslett AL, Muller FJ, Nievergelt CM, Shamir R, Loring JF: Dynamic changes in the copy number of pluripotency and cell proliferation genes in human ESCs and iPSCs during reprogramming and time in culture. Cell Stem Cell. 8: 106-118.
Boehm CD, Cutting GR, Lachtermacher MB, Moser HW, Chong SS: Accurate DNA-based diagnostic and carrier testing for X-linked adrenoleukodystrophy. Mol Genet Metab. 1999, 66: 128-136. 10.1006/mgme.1998.2779.
Yik WY, Steinberg SJ, Moser AB, Moser HW, Hacia JG: Identification of novel mutations and sequence variation in the Zellweger syndrome spectrum of peroxisome biogenesis disorders. Hum Mutat. 2009, 30: E467-480. 10.1002/humu.20932.
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: 861-872. 10.1016/j.cell.2007.11.019.
Itskovitz-Eldor J, Schuldiner M, Karsenti D, Eden A, Yanuka O, Amit M, Soreq H, Benvenisty N: Differentiation of human embryonic stem cells into embryoid bodies compromising the three embryonic germ layers. Mol Med. 2000, 6: 88-95.
Hubbard WC, Moser AB, Liu AC, Jones RO, Steinberg SJ, Lorey F, Panny SR, Vogt RF, Macaya D, Turgeon CT, Tortorelli S, Raymond GV: Newborn screening for X-linked adrenoleukodystrophy (X-ALD): validation of a combined liquid chromatography-tandem mass spectrometric (LC-MS/MS) method. Mol Genet Metab. 2009, 97: 212-220. 10.1016/j.ymgme.2009.03.010.
Zemski Berry KA, Murphy RC: Electrospray ionization tandem mass spectrometry of glycerophosphoethanolamine plasmalogen phospholipids. J Am Soc Mass Spectrom. 2004, 15: 1499-1508. 10.1016/j.jasms.2004.07.009.
Lowry WE, Richter L, Yachechko R, Pyle AD, Tchieu J, Sridharan R, Clark AT, Plath K: Generation of human induced pluripotent stem cells from dermal fibroblasts. Proc Natl Acad Sci USA. 2008, 105: 2883-2888. 10.1073/pnas.0711983105.
Hussein SM, Batada NN, Vuoristo S, Ching RW, Autio R, Narva E, Ng S, Sourour M, Hamalainen R, Olsson C, Lundin K, Mikkola M, Trokovic R, Peitz M, Brustle O, Bazett-Jones DP, Alitalo K, Lahesmaa R, Nagy A, Otonkoski T: Copy number variation and selection during reprogramming to pluripotency. Nature. 2011, 471: 58-62. 10.1038/nature09871.
Taapken SM, Nisler BS, Newton MA, Sampsell-Barron TL, Leonhard KA, McIntire EM, Montgomery KD: Karotypic abnormalities in human induced pluripotent stem cells and embryonic stem cells. Nat Biotechnol. 2011, 29: 313-314. 10.1038/nbt.1835.
Martins-Taylor K, Nisler BS, Taapken SM, Compton T, Crandall L, Montgomery KD, Lalande M, Xu RH: Recurrent copy number variations in human induced pluripotent stem cells. Nat Biotechnol. 2011, 29: 488-491. 10.1038/nbt.1890.
Corzo D, Gibson W, Johnson K, Mitchell G, LePage G, Cox GF, Casey R, Zeiss C, Tyson H, Cutting GR, Raymond GV, Smith KD, Watkins PA, Moser AB, Moser HW, Steinberg SJ: Contiguous deletion of the X-linked adrenoleukodystrophy gene (ABCD1) and DXS1357E: a novel neonatal phenotype similar to peroxisomal biogenesis disorders. Am J Hum Genet. 2002, 70: 1520-1531. 10.1086/340849.
Shimizu T: Lipid mediators in health and disease: enzymes and receptors as therapeutic targets for the regulation of immunity and inflammation. Annu Rev Pharmacol Toxicol. 2009, 49: 123-150. 10.1146/annurev.pharmtox.011008.145616.
Wilson R, Sargent JR: Lipid and fatty acid composition of brain tissue from adrenoleukodystrophy patients. J Neurochem. 1993, 61: 290-297. 10.1111/j.1471-4159.1993.tb03567.x.
Rinn JL, Bondre C, Gladstone HB, Brown PO, Chang HY: Anatomic demarcation by positional variation in fibroblast gene expression programs. PLoS Genet. 2006, 2: e119-10.1371/journal.pgen.0020119.
Koch J, Pranjic K, Huber A, Ellinger A, Hartig A, Kragler F, Brocard C: PEX11 family members are membrane elongation factors that coordinate peroxisome proliferation and maintenance. J Cell Sci. 2010, 123: 3389-3400. 10.1242/jcs.064907.
Schrader M, Reuber BE, Morrell JC, Jimenez-Sanchez G, Obie C, Stroh TA, Valle D, Schroer TA, Gould SJ: Expression of PEX11beta mediates peroxisome proliferation in the absence of extracellular stimuli. J Biol Chem. 1998, 273: 29607-29614. 10.1074/jbc.273.45.29607.
Hoek RM, Ruuls SR, Murphy CA, Wright GJ, Goddard R, Zurawski SM, Blom B, Homola ME, Streit WJ, Brown MH, Barclay AN, Sedgwick JD: Down-regulation of the macrophage lineage through interaction with OX2 (CD200). Science. 2000, 290: 1768-1771.
Liu Y, Bando Y, Vargas-Lowy D, Elyaman W, Khoury SJ, Huang T, Reif K, Chitnis T: CD200R1 agonist attenuates mechanisms of chronic disease in a murine model of multiple sclerosis. J Neurosci. 2010, 30: 2025-2038. 10.1523/JNEUROSCI.4272-09.2010.
Chen HY, White E: Role of autophagy in cancer prevention. Cancer Prev Res (Phila). 2011, 4: 973-983. 10.1158/1940-6207.CAPR-10-0387.
Nakanishi M, Rosenberg DW: Roles of cPLA2alpha and arachidonic acid in cancer. Biochim Biophys Acta. 2006, 1761: 1335-1343. 10.1016/j.bbalip.2006.09.005.
Li Y, Tong X, Rumala C, Clemons K, Wang S: Thrombospondin1 deficiency reduces obesity-associated inflammation and improves insulin sensitivity in a diet-induced obese mouse model. PLoS One. 2011, 6: e26656-10.1371/journal.pone.0026656.
Martin-Manso G, Galli S, Ridnour LA, Tsokos M, Wink DA, Roberts DD: Thrombospondin 1 promotes tumor macrophage recruitment and enhances tumor cell cytotoxicity of differentiated U937 cells. Cancer Res. 2008, 68: 7090-7099. 10.1158/0008-5472.CAN-08-0643.
Tournier JN, Rossi Paccani S, Quesnel-Hellmann A, Baldari CT: Anthrax toxins: a weapon to systematically dismantle the host immune defenses. Mol Aspects Med. 2009, 30: 456-466. 10.1016/j.mam.2009.06.002.
Arora K, Gwinn WM, Bower MA, Watson A, Okwumabua I, MacDonald HR, Bukrinsky MI, Constant SL: Extracellular cyclophilins contribute to the regulation of inflammatory responses. J Immunol. 2005, 175: 517-522.
Solorzano C, Zhu C, Battista N, Astarita G, Lodola A, Rivara S, Mor M, Russo R, Maccarrone M, Antonietti F, Duranti A, Tontini A, Cuzzocrea S, Tarzia G, Piomelli D: Selective N-acylethanolamine-hydrolyzing acid amidase inhibition reveals a key role for endogenous palmitoylethanolamide in inflammation. Proc Natl Acad Sci USA. 2009, 106: 20966-20971. 10.1073/pnas.0907417106.
Buyse M, Sitaraman SV, Liu X, Bado A, Merlin D: Luminal leptin enhances CD147/MCT-1-mediated uptake of butyrate in the human intestinal cell line Caco2-BBE. J Biol Chem. 2002, 277: 28182-28190. 10.1074/jbc.M203281200.
Maccarrone M, Battista N, Centonze D: The endocannabinoid pathway in Huntington's disease: a comparison with other neurodegenerative diseases. Prog Neurobiol. 2007, 81: 349-379. 10.1016/j.pneurobio.2006.11.006.
Tsuboi K, Sun YX, Okamoto Y, Araki N, Tonai T, Ueda N: Molecular characterization of N-acylethanolamine-hydrolyzing acid amidase, a novel member of the choloylglycine hydrolase family with structural and functional similarity to acid ceramidase. J Biol Chem. 2005, 280: 11082-11092. 10.1074/jbc.M413473200.
Isenberg JS, Jia Y, Fukuyama J, Switzer CH, Wink DA, Roberts DD: Thrombospondin-1 inhibits nitric oxide signaling via CD36 by inhibiting myristic acid uptake. J Biol Chem. 2007, 282: 15404-15415. 10.1074/jbc.M701638200.
Stolz A, Hammond L, Lou H, Takikawa H, Ronk M, Shively JE: cDNA cloning and expression of the human hepatic bile acid-binding protein. A member of the monomeric reductase gene family. J Biol Chem. 1993, 268: 10448-10457.
Fourcade S, Savary S, Gondcaille C, Berger J, Netik A, Cadepond F, El Etr M, Molzer B, Bugaut M: Thyroid hormone induction of the adrenoleukodystrophy-related gene (ABCD2). Mol Pharmacol. 2003, 63: 1296-1303. 10.1124/mol.63.6.1296.
Jang J, Kang HC, Kim HS, Kim JY, Huh YJ, Kim DS, Yoo JE, Lee JA, Lim B, Lee J, Yoon TM, Park IH, Hwang DY, Daley GQ, Kim DW: Induced pluripotent stem cell models from X-linked adrenoleukodystrophy patients. Ann Neurol. 2011, 70: 402-409. 10.1002/ana.22486.
Laurent LC, Ulitsky I, Slavin I, Tran H, Schork A, Morey R, Lynch C, Harness JV, Lee S, Barrero MJ, Ku S, Martynova M, Semechkin R, Galat V, Gottesfeld J, Izpisua Belmonte JC, Murry C, Keirstead HS, Park HS, Schmidt U, Laslett AL, Muller FJ, Nievergelt CM, Shamir R, Loring JF: Dynamic changes in the copy number of pluripotency and cell proliferation genes in human ESCs and iPSCs during reprogramming and time in culture. Cell Stem Cell. 2011, 8: 106-118. 10.1016/j.stem.2010.12.003.
Mayshar Y, Ben-David U, Lavon N, Biancotti JC, Yakir B, Clark AT, Plath K, Lowry WE, Benvenisty N: Identification and classification of chromosomal aberrations in human induced pluripotent stem cells. Cell Stem Cell. 2010, 7: 521-531. 10.1016/j.stem.2010.07.017.
Gore A, Li Z, Fung HL, Young JE, Agarwal S, Antosiewicz-Bourget J, Canto I, Giorgetti A, Israel MA, Kiskinis E, Lee JH, Loh YH, Manos PD, Montserrat N, Panopoulos AD, Ruiz S, Wilbert ML, Yu J, Kirkness EF, Izpisua Belmonte JC, Rossi DJ, Thomson JA, Eggan K, Daley GQ, Goldstein LS, Zhang K: Somatic coding mutations in human induced pluripotent stem cells. Nature. 2011, 471: 63-67. 10.1038/nature09805.
Chiang YJ, Kole HK, Brown K, Naramura M, Fukuhara S, Hu RJ, Jang IK, Gutkind JS, Shevach E, Gu H: Cbl-b regulates the CD28 dependence of T-cell activation. Nature. 2000, 403: 216-220. 10.1038/35003235.
Gruber T, Hermann-Kleiter N, Hinterleitner R, Fresser F, Schneider R, Gastl G, Penninger JM, Baier G: PKC-theta modulates the strength of T cell responses by targeting Cbl-b for ubiquitination and degradation. Sci Signal. 2009, 2: ra30-10.1126/scisignal.2000046.
Menasche G, Pastural E, Feldmann J, Certain S, Ersoy F, Dupuis S, Wulffraat N, Bianchi D, Fischer A, Le Deist F, de Saint Basile G: Mutations in RAB27A cause Griscelli syndrome associated with haemophagocytic syndrome. Nat Genet. 2000, 25: 173-176. 10.1038/76024.
de Saint Basile G: Inherited Hemophagocytic Syndromes. Primary Immunodeficiency Diseases: A Molecular and Genetic Approach. Edited by: Ochs HD, Smith CIE, Puck JM. 2007, New York: Oxford University Press, 578-588.
Li X, Baumgart E, Morrell JC, Jimenez-Sanchez G, Valle D, Gould SJ: PEX11 beta deficiency is lethal and impairs neuronal migration but does not abrogate peroxisome function. Mol Cell Biol. 2002, 22: 4358-4365. 10.1128/MCB.22.12.4358-4365.2002.
Ahlemeyer B, Gottwald M, Baumgart-Vogt E: Deletion of a single allele of the Pex11beta gene is sufficient to cause oxidative stress, delayed differentiation and neuronal death in mouse brain. Dis Model Mech. 2012, 5: 125-140. 10.1242/dmm.007708.
Koning N, Swaab DF, Hoek RM, Huitinga I: Distribution of the immune inhibitory molecules CD200 and CD200R in the normal central nervous system and multiple sclerosis lesions suggests neuron-glia and glia-glia interactions. J Neuropathol Exp Neurol. 2009, 68: 159-167. 10.1097/NEN.0b013e3181964113.
Lanz TV, Ding Z, Ho PP, Luo J, Agrawal AN, Srinagesh H, Axtell R, Zhang H, Platten M, Wyss-Coray T, Steinman L: Angiotensin II sustains brain inflammation in mice via TGF-beta. J Clin Invest. 2010, 120: 2782-2794. 10.1172/JCI41709.
Ning M, Sarracino DA, Kho AT, Guo S, Lee SR, Krastins B, Buonanno FS, Vizcaino JA, Orchard S, McMullin D, Wang X, Lo EH: Proteomic temporal profile of human brain endothelium after oxidative stress. Stroke. 2011, 42: 37-43. 10.1161/STROKEAHA.110.585703.
Wei Z, Chigurupati S, Arumugam TV, Jo DG, Li H, Chan SL: Notch activation enhances the microglia-mediated inflammatory response associated with focal cerebral ischemia. Stroke. 2011, 42: 2589-2594. 10.1161/STROKEAHA.111.614834.
Agrawal SM, Yong VW: The many faces of EMMPRIN - roles in neuroinflammation. Biochim Biophys Acta. 2011, 1812: 213-219. 10.1016/j.bbadis.2010.07.018.
Khan M, Singh J, Singh I: Plasmalogen deficiency in cerebral adrenoleukodystrophy and its modulation by lovastatin. J Neurochem. 2008, 106: 1766-1779.
Ofman R, Dijkstra IM, van Roermund CW, Burger N, Turkenburg M, van Cruchten A, van Engen CE, Wanders RJ, Kemp S: The role of ELOVL1 in very long-chain fatty acid homeostasis and X-linked adrenoleukodystrophy. EMBO Mol Med. 2010, 2: 90-97. 10.1002/emmm.201000061.
Pujol A, Ferrer I, Camps C, Metzger E, Hindelang C, Callizot N, Ruiz M, Pampols T, Giros M, Mandel JL: Functional overlap between ABCD1 (ALD) and ABCD2 (ALDR) transporters: a therapeutic target for X-adrenoleukodystrophy. Hum Mol Genet. 2004, 13: 2997-3006. 10.1093/hmg/ddh323.
Raymond GV, Jones RO, Moser AB: Newborn screening for adrenoleukodystrophy: implications for therapy. Mol Diagn Ther. 2007, 11: 381-384. 10.1007/BF03256261.
Moser HW: Peroxisomal disorders. Clin Biochem. 1991, 24: 343-351. 10.1016/0009-9120(91)80009-R.
Pujol A, Hindelang C, Callizot N, Bartsch U, Schachner M, Mandel JL: Late onset neurological phenotype of the X-ALD gene inactivation in mice: a mouse model for adrenomyeloneuropathy. Hum Mol Genet. 2002, 11: 499-505. 10.1093/hmg/11.5.499.
The authors declare that they have no competing interests.
XMW and WYY derived the iPSCs from patient and healthy donor fibroblast cultures and conducted immunostaining analysis of protein pluripotency biomarkers. XMW conducted all the global gene expression, genetic and epigenetic analyses. PZ assisted XMW in the teratoma analysis and also provided technical advice in cellular reprogramming and maintaining iPSCs. WL provided technical advice in cellular reprogramming and maintaining iPSCs. DS supervised the histological analysis of teratomas and advised XMW in the data analysis. SJS carried out the biochemical analyses in this project. JGH and XMW were involved in the overall design and conception of the project, statistical analysis of all data sets and wrote the manuscript with the help of all the other authors. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Log-transformed gene expression scores from childhood cerebral adrenoleukodystrophy (CCALD) patient and healthy donor control fibroblasts and induced pluripotent stem cells (iPSCs). A complete list of log-transformed gene expression scores from CCALD patient and healthy donor control fibroblasts and iPSCs is provided. (XLSX 7 MB)
Additional file 3: DNA methylation assay beta-values from childhood cerebral adrenoleukodystrophy (CCALD) patient and healthy donor control fibroblasts and induced pluripotent stem cells (iPSCs). A partial list of DNA methylation assay beta-values from CCALD patient and healthy donor control fibroblasts and iPSCs is provided. (XLSX 18 MB)
Additional file 5: Hierarchical clustering analysis gene expression data from childhood cerebral adrenoleukodystrophy (CCALD) patient and control fibroblasts and induced pluripotent stem cells (iPSCs) based on pluripotency genes. The analysis was based on gene expression data from 30 pluripotency genes reported in reference  and performed with average-linkage and Euclidean distance. (PPTX 261 KB)
Additional file 6: Differential gene expression between childhood cerebral adrenoleukodystrophy (CCALD) patient and healthy donor control fibroblasts. A complete list of all gene expression scores and differentially expressed genes (DEGs) between CCALD patient and healthy donor control fibroblasts is provided. (XLSX 13 MB)
Additional file 8: The Database for Annotation, Visualization and Integrated Discovery (DAVID) annotation of differentially expressed genes (DEGs) found in childhood cerebral adrenoleukodystrophy (CCALD) patient and control donor fibroblasts. DAVID annotation of DEGs found in CCALD patient and control donor fibroblasts is provided. (XLSX 111 KB)
Additional file 9: Differential gene expression between childhood cerebral adrenoleukodystrophy (CCALD) patient and healthy donor control iPSCs. A complete list of differentially expressed genes (DEGs) between CCALD patient and healthy donor control induced pluripotent stem cells (iPSCs) is provided. (XLSX 37 KB)
Additional file 10: Gene Ontology (GO) analysis of differentially expressed genes (DEGs) between childhood cerebral adrenoleukodystrophy (CCALD) patient and healthy donor control induced pluripotent stem cells (iPSCs). GO analysis of DEGs between CCALD patient and healthy donor control iPSCs is provided. (XLSX 77 KB)
Additional file 11: Database for Annotation, Visualization and Integrated Discovery (DAVID) annotation of differentially expressed genes (DEGs) found in childhood cerebral adrenoleukodystrophy (CCALD) patient and control donor induced pluripotent stem cells (iPSCs). DAVID annotation of DEGs found in CCALD patient and control donor iPSCs is provided. (XLSX 78 KB)
About this article
Cite this article
Wang, X., Yik, W.Y., Zhang, P. et al. The gene expression profiles of induced pluripotent stem cells from individuals with childhood cerebral adrenoleukodystrophy are consistent with proposed mechanisms of pathogenesis. Stem Cell Res Ther 3, 39 (2012). https://doi.org/10.1186/scrt130
- Very Long Chain Fatty Acid
- DAVID Bioinformatics Resource
- Differentially Methylated Locus
- Plasmalogen Level