Significant transcriptional changes in 15q duplication but not Angelman syndrome deletion stem cell-derived neurons
© The Author(s). 2018
Received: 7 June 2017
Accepted: 15 January 2018
Published: 27 January 2018
The inability to analyze gene expression in living neurons from Angelman (AS) and Duplication 15q (Dup15q) syndrome subjects has limited our understanding of these disorders at the molecular level.
Here, we use dental pulp stem cells (DPSC) from AS deletion, 15q Duplication, and neurotypical control subjects for whole transcriptome analysis. We identified 20 genes unique to AS neurons, 120 genes unique to 15q duplication, and 3 shared transcripts that were differentially expressed in DPSC neurons vs controls.
Copy number correlated with gene expression for most genes across the 15q11.2-q13.1 critical region. Two thirds of the genes differentially expressed in 15q duplication neurons were downregulated compared to controls including several transcription factors, while in AS differential expression was restricted primarily to the 15q region. Here, we show significant downregulation of the transcription factors FOXO1 and HAND2 in neurons from 15q duplication, but not AS deletion subjects suggesting that disruptions in transcriptional regulation may be a driving factor in the autism phenotype in Dup15q syndrome. Downstream analysis revealed downregulation of the ASD associated genes EHPB2 and RORA, both genes with FOXO1 binding sites. Genes upregulated in either Dup15q cortex or idiopathic ASD cortex both overlapped significantly with the most upregulated genes in Dup15q DPSC-derived neurons.
Finding a significant increase in both HERC2 and UBE3A in Dup15q neurons and significant decrease in these two genes in AS deletion neurons may explain differences between AS deletion class and UBE3A specific classes of AS mutation where HERC2 is expressed at normal levels. Also, we identified an enrichment for FOXO1-regulated transcripts in Dup15q neurons including ASD-associated genes EHPB2 and RORA indicating a possible connection between this syndromic form of ASD and idiopathic cases.
The study of neurogenetic syndromes in the laboratory setting is inherently complicated by a lack of easy access to living neurons for gene expression and electrophysiological studies. Although a variety of stem cell types has been used in the laboratory, including embryonic stem cells, blood stem cells, and most notably induced pluripotent stem cells (iPSCs) derived from skin biopsies , there are still technical issues complicating the use of these limited resources. As an alternative to iPSC, dental pulp stem cells (DPSC) are neural crest-derived multipotent cells that can be easily obtained from deep within the pulp of shed (deciduous) teeth that either fall out on their own or are extracted in the pediatric clinic setting when ready for exfoliation. DPSC have been differentiated into a variety of cell types including osteoblasts, hepatocytes, insulin-producing cells, and even neurons . Furthermore, neurons derived from DPSC display action potentials (Na+ currents) [3, 4] and are positive for neuronal markers like MAP2 and TUJ1 . In our hands, DPSC-derived neurons have not been able to survive in vitro long enough to become mature neurons unless they are first immortalized. We are confident, however, that these immature neuronal cultures are an excellent resource for gene expression studies in the disease state, especially for neurodevelopmental disorders where the phenotypes appear at or near birth. In addition, we recently showed that DPSC maintain an epigenetic landscape more similar to embryonic stem cells than iPSC making them ideal for gene expression studies . DPSC can also be easily made into neurons without the use of viral vectors. Finally, the collection of DPSC from individuals with these syndromes is quite easy and non-invasive. Sample collection can even take place through shipments of exfoliated teeth by parents from remote locations where the subjects reside.
Recently, our laboratory performed an extensive molecular characterization of both DPSC and DPSC-derived neuronal cultures . We established that the process of differentiation from DPSC to neurons is primarily driven by the downregulation of the REST transcriptional repressor, turning on a number of genes involved in neuronal differentiation and function producing cells that show a strong induction of the neuron-specific MAP2 marker and residual expression of glial marker GFAP indicating that these cultures are representative of the mixed neuron/glia characteristics of the human brain .
While it is apparent that the key components of the Angelman syndrome phenotype are caused by loss of function of the maternally expressed UBE3A gene in the brain [12, 13], it is still unclear how loss of UBE3A can directly lead to the observed neurological defects including synaptic plasticity problems, ataxia, seizures, and intellectual disability . Significant effort has surrounded the search for direct targets of the UBE3A ubiquitin ligase in model systems from flies to mice [15–20]. However, little progress has been made in identifying a single substrate or even pathway that can explain the plasticity, motor, and learning/memory deficits in these animal models. One aspect of UBE3A function that has largely been ignored is the observation that UBE3A can act as a transcriptional regulator in addition to its role as a ubiquitin ligase [21, 22].
Almost nothing is known about the underlying pathogenesis of Duplication 15q syndrome, caused by interstitial and isodicentric duplications of 15q11.2-q13.1 that may involve several genes in the duplicated region to produce the autism spectrum disorder (ASD) and seizure phenotypes found in patients with this disorder [7, 10, 23]. Some studies have indicated that the maternally expressed UBE3A gene may be a major contributor to the autism phenotype in these individuals since in both interstitial and isodicentric duplications, there is a clear bias towards maternally derived or inherited duplications where ASD is present [24–26]. Some evidence also suggests that the seizure phenotype in these individuals may be related to over-expression of a GABA (A) receptor subunit gene cluster also located in the critical BP2-BP3 duplicated region . However, it has been difficult to definitively establish that the GABA receptor genes are indeed elevated in 15q Duplication post-mortem brain tissue [27, 28] or blood . The contribution of other genes, both within the duplicated region and elsewhere in the genome, to the 15q Duplication syndrome phenotype is almost entirely speculative at this time since no comprehensive analysis of gene expression from Duplication 15q patient neurons is currently available.
In this study, we interrogated genome-wide gene expression changes in DPSC and DPSC neuronal cultures from individuals with both AS deletions and 15q duplications in an effort to establish which gene networks may be involved in the phenotypes of these two syndromes, but also to directly compare the consequences on gene expression for duplication/deletion of the same chromosomal region that ties these two disorders at the genomic level. Whole genome gene expression analysis of patient-derived DPSC neuronal cultures from these disorders has revealed several key regulators of gene expression and neurologically relevant genes which may substantially contribute to the 15q duplication syndrome phenotype, while the number of gene expression changes in the Angelman syndrome set remained quite small. Perhaps most surprising was the extensive change observed for transcription factor encoded genes in the 15q Duplication syndrome subjects, suggesting global transcriptional changes mediated by transcription factors FOXO1 and/or HAND2 may be driving the pathogenesis of ASD in this syndrome.
Generation of dental pulp stem cell cultures
Control teeth were obtained through the Department of Pediatric Dentistry at the University of Tennessee Health Science Center (UTHSC). Teeth from AS deletion and Dup15q subjects were shipped from various locations by the parents of children with these conditions. The UTHSC Institutional Review Board (IRB) approved this study, and informed consent was obtained in accordance with IRB policy. Immediately following extraction of a loose tooth in the clinic, it was placed in transportation media (DMEM/F12 50/50 mix with HEPES (Fisher Scientific, Waltham, MA) and 100 U/mL penicillin, 100 μg/mL streptomycin). DPSC were isolated and cultured as previously described with slight modifications . Briefly, the pulp was minced and digested in a solution of 3 mg/mL Collagenase type I and 4 mg/mL Dispase II for 1 h at 37 °C. Cells were seeded in poly-D-Lysine 12-well dishes and maintained under standard conditions (37 °C, 5% CO2) in DMEM/F12 1:1, 10% fetal bovine serum (FBS) (Fisher Scientific, Waltham, MA), 10% newborn calf serum (NCS) (Fisher Scientific, Waltham, MA) and 100 U/mL penicillin, and 100 μg/mL streptomycin (Pen/Strep). Sub-confluent cultures were passaged regularly with 0.1 μM HyQTase (HyClone).
Fluorescent in situ hybridization breakpoint analysis
Fluorescent in situ hybridization was used to determine the approximate breakpoints for either deletion or duplication DPSC lines. We used the pepsin method of 0.1 mg pepsin (Sigma-Aldrich, St. Louis, MO) in 40 ml of 0.01N HCl in order to permeabilize the nuclear membrane for access to the DNA. Three different probes localized in chromosome 15q were used: NIPA1 (BP1-BP2) and TRPM1 (BP4-BP5) in the red spectrum and SNRPN (BP2-BP3) in the green spectrum (SURE FISH, Agilent, Foster City, CA). Hybridization was carried out at RT overnight. Fluorescent signal analysis was performed according to published protocols  using a Leica DM6000 upright fluorescent microscope running the Leica application Suite Imaging software (Leica Microsystems, Wetzlar, Germany). At least 15 preparations were counted per cell line.
Early passage (2–3) primary DPSC were converted to neurons according to a previously published protocol . 20,000 cells/cm2 were seeded in poly-D-lysine-coated 12-well plates or T-25 flasks in DMEM/F12 (1:1), 2.5% FCS, 100 U/mL penicillin, and 100 μg/mL streptomycin, and cultured for 24 h. Epigenetic reprogramming was performed by exposing the DPSC to 10 μM 5-azacytidine (Acros Scientific, Geel, Belgium) in DMEM/F12 containing 2.5% FCS (Fisher Scientific, Waltham, MA) and 10 ng/mL bFGF for 48 h. Neural differentiation was induced by exposing the cells to 250 μM IBMX, 50 μM forskolin, 200 nM TPA (Sigma Aldrich, St. Louis, MO), 1 mM db-cAMP, 10 ng/mL bFGF, 10 ng/mL NGF (Invitrogen, Carlsbad, CA), 30 ng/mL NT-3 (Peprotech, Rocky Hill, NJ), and 1% insulin-transferrin-sodium selenite premix (ITS) in DMEM/F12 for 3 days. At the end of the neural induction treatment, the cells were washed with 1X PBS. Neuronal maturation was performed by maintaining the cells in Neurobasal A media (Invitrogen, Carlsbad, CA) supplemented with 1 mM dbcAMP, 1% N2, 1% B27, and 30 ng/mL NT-3 and 1X Glutamax for 3 weeks.
RNA sequencing of DPSC and DPSC-neurons
Total RNA was extracted from ~ 3-week cultured neurons and DPSC using the miRNeasy Mini kit from Qiagen (Foster City, CA). Total RNA was analyzed on an Agilent Bioanalyzer 6000 pico chip and determined to have a RIN of > 8.0. 5 ng of total RNA and was then used to prepare dscDNA using the NuGEN Ovation V2 kit. Five hundred nanograms of this dscDNA then used to prepare libraries for sequencing using the AB Library Builder™ Fragment library Kit for 5500 Genetic Analysis Systems on a Library Builder system. Libraries were amplified 7 cycles before 5500 Wildfire primers added by 5 cycles of fusion primer amplification as described in the 5500 Wildfire manual (Applied Biosystems, Foster City, CA). Before sequencing on the ABI system, small aliquots of this material pooled and sequenced on an Ion Torrent PGM 314 chip after additional amplification with PGM fusion primers to estimate sample abundance. Barcode quantification data from the PGM was used to pool the samples with Wildfire primers equally. Following this final pooling, the library pools were examined and quantified on an Agilent high sensitivity DNA chip, immobilized on flow cells for the 5500 Wildfire instrument, and sequenced (50 bp reads) on this instrument. On average, a total of 30 million reads were produced per sample for downstream analysis. Complete sequencing reads for this study are available on the Gene Expression Omnibus (GEO) website using query [GEO: GSE67122].
ABI files were converted from XSQ format to separate quality and colorspace fasta files for analysis using the TopHat with Bowtie pipeline . Short transcript fragments were mapped to the reference human genome (GRCh38; accession number GRCh38.77) using TopHat v2.0.13  with Bowtie v2.1.0  and Samtools v0.1.19. Mapped reads were indexed to known genes using HTseq v0.6.1 . Differentially expressed transcripts were identified using the DESeq2 bioconductor package . Heatmaps were generated using the R software package (v 3.1.1) . Venn diagrams were generated online using Gene List Venn Diagram software, additional descriptive analysis was performed using Gene Set Enrichment Analysis (GSEA) using a rank-ordered gene list based on fold change that was used to identify pathways and transcription factor binding motifs enriched in particular data sets . Statistical significance for GSEA analysis was set at q ≤ 0.25 as per the developer’s suggestion. Gene sets were obtained from MSigDB version 4.0. Gene sets that were found to be significantly up- or downregulated (q ≤ 0.05) were also analyzed as a gene list with default stringency settings using The Database for Annotation, Visualization and Integrated Discovery (DAVID v6.7).
Quantitative real-time PCR analysis
Total RNA was extracted from differentiated neurons using a Direct-zol RNA miniprep plus kit (Zymogen R0270) according to the manufacturer’s instructions and quantified by spectrophotometry (NanoDrop Technologies). Five hundred nanograms of total RNA was used as input for each cDNA synthesis reaction using the Transcriptor First Strand cDNA Synthesis Kit (Roche, 04379012001). cDNA reactions were diluted 1:10 and three technical replicates were performed on each sample for each gene tested using the default cycling parameters of the Light Cycler 480 system (Roche). For primers and probes, see Additional file 1: Table S9. The crossing point (Cp) value was calculated using the absolute quantitation algorithm (Roche) for each sample. Delta Cp was calculated by subtracting the GAPDH control gene Cp from the Cp of the test gene. Statistical significance was determined by unpaired t tests corrected for multiple comparisons.
Quantitative western blot analysis
Infrared (IR) Western blots and quantification from DPSC neurons were performed as previously described in Urraca et al.  using six unrelated DPSC lines from neurotypical controls and idic15q. These six independent DPSC lines were not run for RNAseq analysis as they were collected after the RNAseq runs. We used this sample to validate the RNAseq findings in an independent cohort. The α-RORA antibody (Aviva Biosystems # OAAN02076) was used at a dilution of 1:500 in conjunction with Li-Cor α-Rabbit 800 secondary antibody at 1:2500. Loading control was α-GAPDH with α-Goat 680 at 1:2500. Each individual sample was normalized to the GAPDH signal as a loading control, and quantification was performed using all 12 samples in three independent Western blots by Student’s t test (control vs idic15).
For the analysis of neural differentiation markers, we used three control and three idic15 lines per group. All primary antibodies used were from AbCam (Cambridge, MA). Primary concentrations were 1:1000 for GABAA (ab98968), NLG1 (ab26305), PSD95 (ab12093), and UBE3A (ab126765) and 1:2000 for MAP2 (ab5392) and GAPDH (ab157156). Secondary IR antibodies and quantification were performed as above for RORA analysis.
Comparison to postmortem Dup15q and ASD brain tissue
The most significant differentially expressed genes between control and Dup15q DPSC neurons were directly compared to a combined list of temporal and frontal cortical gene expression from a study using postmortem brain from Dup15q and ASD individuals (see ). This analysis was also repeated after removing any bias from genes in the 15q11.2-q13.1 duplicated region. We performed PANTHER over-representation tests (release 20150430) on the significantly overlapping intersections between the 1500 genes that went up most significantly in DPSC neurons in Dup15q vs controls and the genes that went up significantly in postmortem Dup15q cortex or postmortem idiopathic ASD cortex at an FDR of 5% . These comparisons were performed using default settings for every class of annotations, filtering only those terms significant at a Bonferroni-corrected p ≤ 0.05 for either of the two overlaps. The background was the intersection of genes expressed in Parikshak et al. and genes assigned a p value by DESeq in the current study.
Identification of breakpoints in DPSC lines
FISH results to determine breakpoints in DPSC lines
All Control DPSC showed the normal number of signals for each probe indicating no duplications or deletions in the region between BP1-BP5 in Fig. 1. AS sample TP-055 did not show a signal for probe NIPA1 indicating that this deletion is a BP1-BP3 deletion. AS samples TP-059 and TP-078 appear to be BP2-BP3 deletions. Sample TP-041 appears to be a BP2-BP3 interstitial duplication of 15q, while samples TP-044 and TP-058 are isodicentric duplications that include four copies of the entire BP1-BP5 region. See Additional file 2: Fig. S1 for example FISH images.
mRNA-seq analysis of the 15q duplication/deletion critical region
Whole genome mRNA sequencing (mRNA-seq) was performed on nine subjects for both DPSC and DPSC neuronal cultures with a total of 18 samples (DPSC: 3 control, 3 AS deletion, 3 15q duplication and DPSC-neurons: 3 control, 3 AS deletion, 3 15q duplication). DPSC lines were differentiated into neurons using previously published protocols , and the DPSC neuronal cultures used for mRNA-seq studies were allowed to mature for at least 3 weeks prior to RNA extraction. Previous studies indicate that longer neuronal maturation times (up to 10 weeks) result in significant cell death and consequently poor RNA yield and that massive gene expression changes have already occurred in these cultures by ~ 3 weeks including a switch from MAP2 negative to MAP2 positive status . Based on these studies and the potential high RNA yield, we chose to use ≥ 3-week old DPSC neuron cultures for the present study. Tooth samples consisted of three neurotypical controls obtained locally, three AS deletion cases, and three 15q duplication cases (one interstitial BP2-BP3 duplication and two idic(15) cases which included BP4-BP5) collected by the parent or guardian and shipped to the laboratory.
Next, we looked carefully at gene expression for genes located within the duplicated or deleted region as compared to neurotypical controls. A heatmap of the average gene expression across the three lines per group for protein coding genes across the proximal 15q region from the centromere to BP5 revealed that gene expression is elevated at multiple loci in the Dup15q samples as compared to neurotypical controls in both DPSC and DPSC-derived neurons (Fig. 2b). The AS deletion lines did not differ significantly in gene expression as compared to neurotypical controls for either DPSC or DPSC-derived neurons with a few notable exceptions: UBE3A expression levels followed copy number in both DPSC and neurons going up in 15q duplication and down in AS deletion; the GABA A receptor gene cluster consisting of GABRB3, GABRG3, and GABRA5 each showed significantly elevated gene expression in 15q duplication neural cultures but were also slightly elevated in AS deletion samples; and the HERC2 gene located at the distal breakpoint for the AS deletion (BP3) also showed significantly elevated gene expression in 15q duplication neural cultures and was significantly decreased in AS deletion cultures (Fig. 2c).
Jonckheere trend test for correlation between copy number and gene expression
Gene expression in AS and 15q duplication compared to controls
Global protein coding mRNA gene expression was compared between AS deletion and typical controls as well as 15q duplication and neurotypical controls for both DPSC and DPSC-derived neurons using DEseq to assess fold change and significantly differential whole transcript level expression . Only 13 genes were found to be significantly different in AS deletion subjects and 7 genes in 15q duplication vs controls for the undifferentiated DPSC at an adjusted p value ≤ 0.05. None of these 20 genes was shared between the two disease states. After differentiation into neurons, 123 genes were significantly different (adj. p value ≤ 0.05) in the 15q duplication samples vs control and 23 genes in the AS deletion samples (adj. p value ≤ 0.05). Two of these genes, DPT and AKRIC1, were elevated in both 15q duplication and AS deletion samples (Fc = + 2.8) while MED12L gene expression decreased (Fc = − 2.8) in both groups. None of these genes was located in the PWS/AS critical region. Boxplots of all significantly different genes discussed in the text that are not located in the 15q11.2-q13.1 critical region can be found in Additional file 2: Fig. S2. However, when gene expression was compared between 15q duplication and AS deletion samples (both DPSC and DPSC-derived neurons), there were significant differences in gene expression for several genes in the 15q region. In DPSC NIPA2, UBE3A, HERC2, and ATP10A, all showed significant increase in gene expression (adj. p value ≤ 0.05) in 15q duplication cells vs AS deletion cells with a fold change ranging from 2.3–3.7. In DPSC-derived neurons, six genes (NIPA2, UBE3A, HERC2, NDNL2, CYFIP1, TUBGCP5) showed significantly higher expression in 15q duplication neurons (adj. p value ≤ 0.05) with a fold change ranging from 2.3–3.7 depending on the gene. Perhaps, the most interesting genes differentially regulated in the 15q duplication neurons were the CCL7 gene, a chemokine associated with autoimmune response (adj. p value = 2 × 10−3; logFc = + 2.0), and the two transcription factor genes FOXO1, also known as FOXO1A (adj-pvalue = 3X10−3; Fc = − 3.4), a transcription factor that functions in the AKT signaling cascade and HAND2 (adj. p value = 2.4 × 10−5; Fc = − 6.5). HAND2 plays a role in the development of the sympathetic nervous system and could be involved in several Dup15q phenotypes [41–43]. In the AS deletion neurons, the most interesting single genes that changed were SHISH2 (adj. p value = 3.5 × 10−6; Fc = − 4.9), an antagonist to both WNT and FGF signaling pathways, and TRHDE (adj. p value = 3 × 10−4; Fc = + 5.7) an extracellular peptidase that cleaves neuropeptide thyrotropin-releasing hormone and CRY2 (adj. p value = 2 × 10−2; Fc = − 2.0). CRY2 is a circadian clock gene, and circadian rhythm defects have been found in both fly and mouse models of AS [44, 45].
We tested the null hypothesis that the gene expression trends observed were random or independent of the disease state (i.e., control vs Dup15q vs AS). In both cases, chi-squared testing revealed that the downregulation of gene expression is significant in 15q duplication neuron cultures. Of the 123 genes differentially expressed in the 15q duplication disease state vs control, 44 genes were upregulated and 79 genes were downregulated. A chi-squared fit test revealed that the number of up- and downregulated genes was not equally represented (p value = .6 × 10−3). These results indicate that there are more significantly downregulated than upregulated genes in the 15q duplication neuronal cultures that did not occur by chance alone.
Network analysis of differentially regulated Dup15q genes
Since there was a larger set of genes differentially regulated between Dup15q and control neurons (120) than AS and controls (20), we focused on more detailed studies of the relationships among the proteins encoded by genes up- or downregulated in the Dup15q vs control neurons (Fig. 3a). Using STRING analysis, we identified a core network of 23 proteins predicted to interact physically or through experimental evidence with the UBA52 protein providing the central link to several other proteins mostly downregulated in Dup15q neurons vs controls including the transcription factor FOXO1 and the E2 ubiquitin ligase UBE2A (Fig. 3b). A second cluster of three proteins that were significantly upregulated Fc ≥ + 2.8 fold was also identified that included the secreted immune system factors CCL7, MMP3, and MMP1 (Fig. 3b). These cytokines may play a role in the differentiation process from DPSC to DSPC-neurons, but subsequent analysis of peptide cytokine levels in the media of mature DPSC neuronal cultures failed to recapitulate the gene expression findings (data not shown).
Gene set enrichment analysis (GSEA) of differentially regulated genes
Perhaps, as expected for the AS deletion samples, since very few genes were differentially regulated, there were only a few pathway enrichment results that approached the significance level (FDR q-value ≤ 25%), considered a loose FDR cutoff value. Also of note, there was a significant upregulation for genes that encode steroid biosynthetic pathways (KEGG: “Steroid Biosynthesis” FDR q-value ≤ 0.07% and “Steroid Hormone Biosynthesis” FDR q-value≤ 1.8%; GO: “Steroid Metabolic Process” FDR q-value≤ 24%) and downregulation for genes involved in neural activity like “Chloride Channels” FDR q-value≤ 24% and “Ligand Gated Channels” FDR q-value≤ 23%.
In the Dup15q sample analysis, we identified pathways that were both significantly upregulated and downregulated by this duplication in DPSC-neurons. Consistent with the finding that the cytokine CCL7 showed a > 2 fold increase in gene expression in Dup15q neurons vs controls (adj. p value = 2.6 × 10−3), we found a significant enrichment for cytokine/chemokine activity in the Dup15q samples (FDR q-value ≤ 6% for GO: “Chemokine Activity,” “Chemokine Receptor Binding,” “Cytokine Activity,” and KEGG: Cytokine Cytokine Receptor Interaction FDR q-value ≤ 21%). There was also an increase in Ras/Rho GTPase activity genes (GO: “RAS GTPase activity” FDR q-value ≤ 25%, “Rho GTPase Activity” FDR q-value ≤ 11%, and “GTPase Activator Activity” FDR q-value ≤ 12%). Although not as pronounced as the 2.16 net enrichment score in “Steroid Biosynthesis” identified in AS deletion neurons, there was still a significant 1.51 net enrichment score for these transcripts in the Dup15q neurons as well.
Transcriptional downregulation between Dup15q neurons and controls proved much more interesting. First, there was a general downregulation for genes involved in chromatin organization and biogenesis (GO: FDR q-value ≤ 24%). Although they did not reach statistical significance, there were other downregulated genes involved in this process as well in the nominal p value range of ≤ 1% (notably GO: “Chromatin Assembly and Disassembly” and “Establishment or Maintenance of Chromatin Architecture” nominal p values ≤ 0.05). Second, there was a significant enrichment for genes potentially regulated by transcription factors OCT1 and CEBP in the downregulated set (FDR q-value < 2.7% and < 12% respectively). In addition, there was a nominal enrichment for genes with FOXO1 regulatory binding sites (p value = 7.0 × 10−3), a transcription factor encoding gene which was significantly downregulated in the Dup15q neurons (adj. p value = 3 × 10−3; Fc = − 3.5). The complete list of GSEA enrichment results is available in Additional file 1: Tables S1–S9.
Gene expression in 15q duplication DPSC derived neurons correlates to gene expression in Dup15q and ASD postmortem cortex samples
Changes in dup15q DPSC neurons overlap changes seen in postmortem cortex in both dup15q and idiopathic ASD
Fold enrichment (P value)
↑ in Dupl5q cortex
↓ in Dupl5q cortex
↑ In idiopathic ASD cortex
↓ In idiopathic ASD cortex
↑ 1500 in DPSC neurons
↓ 1500 in DPSC neurons
All genes except 15qll-13 genes
↑ 1500 in DPSC neurons
↓ 1500 in DPSC neurons
DPSC are an excellent in vitro model to study gene expression and even protein expression changes in neurogenetic disorders. They are easily accessible through noninvasive methods and can be shipped from long distances directly to the laboratory . As with any in vitro model, it is not perfect. In the near future, it will be crucial to develop protocols to differentiate these cells into specific neural subtypes coupled with extensive molecular and morphological characterization. That said, it was recently shown that gene expression profiling in neurons derived from iPSC is quite similar to neurons from DPSC and that genes involved in neural development and neurological disorders appear to be conserved in terms of regulation . In addition, the epigenetic landscape of DPSC is more similar to embryonic stem cells than other stem cell types used for neurogenetics research, including iPSC . DPSC-derived neurons in culture are heterogeneous with a mixed of neurons and glia cells. After maturation, there is an increase in neuronal markers, also demonstrated here in Fig. 6, and a decrease in glia markers (4,6); however, the proportion of these cells can vary from culture to culture. In this study, we took a new approach to the evaluation of gene expression in neurons from individuals with two distinct but genomically related neurodevelopmental disorders: Angelman syndrome (AS) deletion (missing one copy of the region), interstitial duplication of 15q11.2-q13.1 (a single extra copy of the region), and the isodicentric duplication of 15q (two extra copies of the region). Here, we perform the first direct comparison of gene expression in 3-week-old live neurons from subjects with duplication vs deletion of the 15q region, both of which can include autism spectrum disorder as part of the phenotype [26, 47]. Previous studies indicated that gene expression in 15q Duplication syndrome lymphoblasts shared a general expression signature with both Fragile X syndrome and idiopathic autism cases perhaps linked by the CYFIPI gene, a gene we also found elevated in DPSC neurons from idic(15) subjects . In another study, gene expression analysis using post mortem brain tissue indicated that at least for genes within the duplicated region (most notably GABRB3 and UBE3A) that expression correlates with copy number . Here, we were able to confirm that copy number correlates significantly with gene expression in the PWS/AS critical region for both Angelman deletion and 15q duplication syndromes (Table 2). The upward trend between copy number and gene expression was significant for two important genes in neurons: UBE3A, a gene encoding the HECT domain ubiquitin protein ligase responsible for the key features of the AS phenotype , and a related HECT domain protein HERC2, a protein which also has been shown to physically interact with and functionally enhance the activity of UBE3A [49, 50].
One limitation in interpreting these studies is the observation that gene expression of UBE3A is still detected in 3-week-old DPSC neuronal cultures that came from maternal AS deletion subjects. There may be several explanations for this phenomena: (1) we were not able to use for gene expression studies DPSC neuronal cultures matured past ~ 5 weeks for gene expression studies because as these cultures mature, the number live neurons decreases; (2) DPSC are neural crest-derived cells from a lineage that may not show imprinted gene expression for UBE3A under the current culturing conditions, although this will have a stronger effect on interpretation of the AS deletion results than the Dup15q results since UBE3A duplications are primarily maternal in origin and therefore UBE3A will be expressed on the duplicated chromosome; (3) the proportion of glial cells, which biallelically express UBE3A in humans, to neurons is high in cultured DPSC neurons, masking the mono-allelic expression of UBE3A in neurons . All of these possibilities will require further investigations using more mature DPSC neuron cultures, differing culture conditions, and possibly a shift to the use of induced dental pulp stem cells, if needed, as in a recent ASD study .
The coordinated gene expression for UBE3A and HERC2, located within the duplicated/deleted region on 15q, may have a significant impact on the phenotypes in both AS and 15q duplication syndromes. In the case of AS, the majority of individuals with this syndrome have a deletion (70%), and among the deletion cases, the BP2-BP3 deletion is more frequent (50%) but still includes both UBE3A and HERC2 . The other ~ 15% of individuals have imprinting center defects, mutations in maternal UBE3A, or paternal uniparental disomy. Individuals with a UBE3A mutations have a milder AS phenotype with severe developmental delay, absence of speech, inappropriate laughter and EEG abnormalities as the most consistent findings and ataxia, epilepsy and microcephaly can present as mild or even absent [53, 54], suggesting that loss of UBE3A expression in neurons alone is not enough to recapitulate all aspects of the AS phenotype. In fact, regulation of UBE3A by HERC2 in deletion cases may enhance the phenotype substantially and could certainly play a role in affecting UBE3A function not only in non-imprinted tissues but also in tissues where a preference for maternal-specific expression has been found like heart and liver  or even in parts of the brain where there is biallelic expression for UBE3A in GFAP positive cells lining the lateral ventricles where neural stem cells develop . In fact, mutations in HERC2 can correlate with an Angelman-like phenotype clearly implicating the physical interaction between HERC2 and UBE3A in some aspect of the syndrome . Given that HERC2 mutations can also cause phenotypes on the autism spectrum in some cohorts , we propose that the elevated levels of gene expression for both UBE3A and HERC2 found here in neurons from 15q Duplication syndrome subjects may have a more profound impact on the autism phenotype in this disorder than previously considered .
With regard to the 15q Duplication syndrome, this study provides important clues to the underlying molecular etiology of the disorder in patient-derived neurons. Gross downregulation of developmental transcription factors such as FOXO1 (Fc = − 3.4), HAND2 (Fc = − 7.3), and SOX11 (Fc = − 4.7) in 15q duplication neurons may have significant pleiotropic effects on other genes regulated by these transcription factors. We identified an enrichment for the downregulation of genes with FOXO1 binding sites (nominal p value = 7X10− 3), for example, that implicates genes regulated by FOXO1 in the syndrome. These genes include proteins involved in chromatin reorganization as well as some specific genes involved in regulating action potentials and dendrite formation. Indeed, there is specific downregulation of both EPHB2 and RORA, two potential FOXO1 gene targets, at the transcript level (Fig. 6a). The validation of downregulation of RORA protein in an independent set of idic15q DPSC neuron cultures vs controls lends further support to the idea that these transcription factors, especially FOXO1, may be the branch point for global downregulation of autism-associated nuclear receptor encoding genes like RORA in Dup15q syndrome.
Although the identification of a transcription factor like FOXO1 that may regulate several downstream ASD-related transcripts is an important finding in this study, perhaps more important was the finding that RORA protein levels are significantly decreased in idic15q vs control neurons. RORA transcript levels were previously found to be significantly decreased in blood from ASD vs control individuals . Moreover, RORA, which is a nuclear hormone receptor, has been shown in neuronal cell lines to transcriptionally regulate ASD-associated genes . RORA is therefore an excellent candidate for future molecular studies of autism in Dup15q syndrome, a major portion of the phenotype .
Moreover, these transcription factors have now been found to be dysregulated in other ASD samples. This prompted us to investigate the correlation between our gene expression data in 15q duplication neurons and a recent study looking at gene expression changes in postmortem brain samples from both Dup15q syndrome and idiopathic ASD cases. We found a clear overlap in 1/3 of upregulated genes in all three comparisons (Table 3). Contrasting DPSC-derived neurons to postmortem cortical samples should be taken with caution as our in vitro model does not consist exclusively of cortical neurons. However, these comparisons are instructive for two reasons: firstly, it shows that DPSC-derived neurons in Dup15q syndrome recapitulate the molecular changes seen in postmortem Dup15q and idiopathic autism brain, even and especially outside the 15q region and that this may therefore be a model with a tractable molecular handle to learn about neurobiology in idiopathic ASD more generally; secondly, it demonstrates that at least some components of the convergent molecular signature in autism brains are inherent to the local neurobiology and are not downstream neurological responses to growing up autistic in a neurotypical world.
It is perhaps no surprise that gene expression in AS neurons did not differ significantly from control neurons for many genes since the primary gene expected to contribute to pathogenesis in AS encodes an E3 ubiquitin ligase that targets other proteins for degradation or other cellular functions at the protein level. Several potential protein targets contributing to pathogenesis have been identified over the years, but no targets have yet emerged that can clearly explain all of the neurological phenotypes in AS on their own [16, 19, 44, 50, 60, 61]. It has been known for some time, however, that UBE3A can act as a transcriptional co-activator of steroid hormone receptors , and our own work in Drosophila indicates that Dube3a expression can regulate gene expression for GTP cyclohydrolase I leading to increased or decreased dopamine in the brains of flies . Here, we show that patient-derived DPSC neuronal cultures do show aberrant gene expression for 23 genes compared to controls. Of these, the one gene that does seem consistent with other UBE3A studies is the CRY2 gene (adj. p value = 2 × 10−2; Fc = − 2.0). CRY2 is a circadian clock gene, and circadian rhythm defects have been found in both fly and mouse models of AS [44, 45]. It remains to be seen if this regulation of CRY2 is contributing to these circadian defects, however, and the limited number of transcriptional changes found in the AS neurons seems to imply that transcriptional regulation of genes other than UBE3A is not a key contributor to the pathogenesis of AS in neurons.
In conclusion, this study demonstrates that DPSC-derived neurons express several neuronal differentiation markers (MAP2, GABAA, NLGN1, and PSD95) and may prove a useful tool for the study of neurogenetic disease. Here, we demonstrated that global gene expression in immature neurons is not significantly altered in AS deletion DPSC neurons but that changes in gene expression found in the 15q duplication DPSC neurons indicate that this particular duplication results in downregulation, rather than up, for genes located outside of the Dup15q critical region through transcription factors like HAND2, OCT1, and FOXO1 and subsequently the ASD associated transcriptional activator RORA. Finally, we found more changes in gene expression in Dup15q than AS cell lines, implying that the mechanism of AS may be primarily through the degradation of UBE3A protein substrates as previously proposed in other studies.
The authors would like to thank the families who donated teeth to this study and the support of the Angelman Syndrome Foundation and the Duplication 15q Alliance in advertising this study to potential subjects. We also acknowledge the assistance of the UTHSC Molecular Resource Center in the production of next-generation sequencing data and the UTHSC CITG for the use of server space for the analysis.
These experiments were funded in part by NIH R21NS075709-02 to L.T.R. and NIH R01DK107535 to D.B.
Availability of data and materials
Complete sequencing reads for this study are available on the Gene Expression Omnibus (GEO) website using query [GEO: GSE67122].
NA, DB, and LTR conceived of the study and analyzed the data. RM, MD, and LTR collected samples. KH, AKV, SG, CV, and JR grew cultures and performed experiments. SS performed FISH studies. TGB, QTT, and DB executed RNAseq and bioinformatics analysis. NA and LTR wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The UTHSC Institutional Review Board (IRB) approved this study, and informed consent was obtained for both participation and publication of results in accordance with IRB policy.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis 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.
- Maury Y, Gauthier M, Peschanski M, Martinat C. Human pluripotent stem cells for disease modelling and drug screening. BioEssays. 2012;34:61–71.View ArticlePubMedGoogle Scholar
- Sedgley CM, Botero TM. Dental stem cells and their sources. Dent Clin N Am. 2012;56:549–61.View ArticlePubMedGoogle Scholar
- Chen J, Lin M, Foxe JJ, Pedrosa E, Hrabovsky A, Carroll R, Zheng D, Lachman HM. Transcriptome comparison of human neurons generated using induced pluripotent stem cells derived from dental pulp and skin fibroblasts. PLoS One. 2013;8:e75682.View ArticlePubMedPubMed CentralGoogle Scholar
- Kiraly M, Porcsalmy B, Pataki A, Kadar K, Jelitai M, Molnar B, Hermann P, Gera I, Grimm WD, Ganss B, et al. Simultaneous PKC and cAMP activation induces differentiation of human dental pulp stem cells into functionally active neurons. Neurochem Int. 2009;55:323–32.View ArticlePubMedGoogle Scholar
- Urraca N, Memon R, El-Iyachi I, Goorha S, Valdez C, Tran QT, Scroggs R, Miranda-Carboni GA, Donaldson M, Bridges D, Reiter LT. Characterization of neurons from immortalized dental pulp stem cells for the study of neurogenetic disorders. Stem Cell Res. 2015;15:722–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Dunaway K, Goorha S, Matelski L, Urraca N, Lein PJ, Korf I, Reiter LT, LaSalle JM. Dental pulp stem cells model early life and imprinted DNA methylation patterns. Stem Cells. 2017;35:981–8.View ArticlePubMedGoogle Scholar
- Hogart A, Wu D, Lasalle JM, Schanen NC. The comorbidity of autism with the genomic disorders of chromosome 15q11.2-q13. Neurobiol Dis. 2010;38:181–91.View ArticlePubMedGoogle Scholar
- Moreno-De-Luca A, Ledbetter DH, Martin CL. Genetic [corrected] insights into the causes and classification of [corrected] cerebral palsies. Lancet Neurol. 2012;11:283–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Pinto D, Delaby E, Merico D, Barbosa M, Merikangas A, Klei L, Thiruvahindrapuram B, Xu X, Ziman R, Wang Z, et al. Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am J Hum Genet. 2014;94:677–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Finucane BM, Lusk L, Arkilo D, Chamberlain S, Devinsky O, Dindot S, Jeste SS, LaSalle JM, Reiter LT, Schanen NC, et al. 15q Duplication Syndrome and Related Disorders. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Stephens K, Amemiya A, editors. GeneReviews® [Internet]. Seattle: University of Washington, Seattle; 1993–2018. 2016.Google Scholar
- Fang P, Lev-Lehman E, Tsai TF, Matsuura T, Benton CS, Sutcliffe JS, Christian SL, Kubota T, Halley DJ, Meijers-Heijboer H, et al. The spectrum of mutations in UBE3A causing Angelman syndrome. Hum Mol Genet. 1999;8:129–35.View ArticlePubMedGoogle Scholar
- Kishino T, Lalande M, Wagstaff J. UBE3A/E6-AP mutations cause Angelman syndrome. Nat Genet. 1997;15:70–3.View ArticlePubMedGoogle Scholar
- Matsuura T, Sutcliffe JS, Fang P, Galjaard RJ, Jiang YH, Benton CS, Rommens JM, Beaudet AL. De novo truncating mutations in E6-AP ubiquitin-protein ligase gene (UBE3A) in Angelman syndrome. Nat Genet. 1997;15:74–7.View ArticlePubMedGoogle Scholar
- Dagli A, Buiting K, Williams CA. Molecular and clinical aspects of Angelman syndrome. Mol Syndromol. 2012;2:100–12.PubMedGoogle Scholar
- Ferdousy F, Bodeen W, Summers K, Doherty O, Wright O, Elsisi N, Hilliard G, O'Donnell JM, Reiter LT. Drosophila Ube3a regulates monoamine synthesis by increasing GTP cyclohydrolase I activity via a non-ubiquitin ligase mechanism. Neurobiol Dis. 2011;41:669–77.View ArticlePubMedGoogle Scholar
- Greer PL, Hanayama R, Bloodgood BL, Mardinly AR, Lipton DM, Flavell SW, Kim TK, Griffith EC, Waldon Z, Maehr R, et al. The Angelman syndrome protein Ube3A regulates synapse development by ubiquitinating arc. Cell. 2010;140:704–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Jiang YH, Armstrong D, Albrecht U, Atkins CM, Noebels JL, Eichele G, Sweatt JD, Beaudet AL. Mutation of the Angelman ubiquitin ligase in mice causes increased cytoplasmic p53 and deficits of contextual learning and long-term potentiation. Neuron. 1998;21:799–811.View ArticlePubMedGoogle Scholar
- Mishra A, Godavarthi SK, Jana NR. UBE3A/E6-AP regulates cell proliferation by promoting proteasomal degradation of p27. Neurobiol Dis. 2009;36:26–34.View ArticlePubMedGoogle Scholar
- Reiter LT, Seagroves TN, Bowers M, Bier E. Expression of the Rho-GEF Pbl/ECT2 is regulated by the UBE3A E3 ubiquitin ligase. Hum Mol Genet. 2006;15:2825–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Weeber EJ, Jiang YH, Elgersma Y, Varga AW, Carrasquillo Y, Brown SE, Christian JM, Mirnikjoo B, Silva A, Beaudet AL, Sweatt JD. Derangements of hippocampal calcium/calmodulin-dependent protein kinase II in a mouse model for Angelman mental retardation syndrome. The Journal of neuroscience : the official journal of the Society for Neuroscience. 2003;23:2634–44.Google Scholar
- Huibregtse JM, Scheffner M, Howley PM. Cloning and expression of the cDNA for E6-AP, a protein that mediates the interaction of the human papillomavirus E6 oncoprotein with p53. Mol Cell Biol. 1993;13:775–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramamoorthy S, Nawaz Z. E6-associated protein (E6-AP) is a dual function coactivator of steroid hormone receptors. Nucl Recept Signal. 2008;6:e006.PubMedPubMed CentralGoogle Scholar
- Battaglia A. The inv dup (15) or idic (15) syndrome (Tetrasomy 15q). Orphanet journal of rare diseases. 2008;3:30.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolton PF, Dennis NR, Browne CE, Thomas NS, Veltman MW, Thompson RJ, Jacobs P. The phenotypic manifestations of interstitial duplications of proximal 15q with special reference to the autistic spectrum disorders. Am J Med Genet. 2001;105:675–85.View ArticlePubMedGoogle Scholar
- Cook EH Jr, Lindgren V, Leventhal BL, Courchesne R, Lincoln A, Shulman C, Lord C, Courchesne E. Autism or atypical autism in maternally but not paternally derived proximal 15q duplication. Am J Hum Genet. 1997;60:928–34.PubMedPubMed CentralGoogle Scholar
- Urraca N, Cleary J, Brewer V, Pivnick EK, McVicar K, Thibert RL, Schanen NC, Esmer C, Lamport D, Reiter LT. The interstitial duplication 15q11.2-q13 syndrome includes autism, mild facial anomalies and a characteristic EEG signature. Autism research : official journal of the International Society for Autism Research. 2013;6:268–79.View ArticleGoogle Scholar
- Hogart A, Leung KN, Wang NJ, Wu DJ, Driscoll J, Vallero RO, Schanen NC, LaSalle JM. Chromosome 15q11-13 duplication syndrome brain reveals epigenetic alterations in gene expression not predicted from copy number. J Med Genet. 2009;46:86–93.View ArticlePubMedGoogle Scholar
- Hogart A, Nagarajan RP, Patzel KA, Yasui DH, Lasalle JM. 15q11-13 GABAA receptor genes are normally biallelically expressed in brain yet are subject to epigenetic dysregulation in autism-spectrum disorders. Hum Mol Genet. 2007;16:691–703.View ArticlePubMedPubMed CentralGoogle Scholar
- Nishimura Y, Martin CL, Lopez AV, Spence SJ, Alvarez-Retuerto AI, Sigman M, Steindler C, Pellegrini S, Schanen NC, Warren ST, Geschwind DH. Genome-wide expression profiling of lymphoblastoid cell lines distinguishes different forms of autism and reveals shared pathways. Hum Mol Genet. 2007;16:1682–98.View ArticlePubMedGoogle Scholar
- Goorha S, Reiter LT. Culturing and neuronal differentiation of human dental pulp stem cells. Curr Protoc Hum Genet. 2017;92:21 26 21–10.Google Scholar
- Van Hummelen P, Lowe XR, Wyrobek AJ. Simultaneous detection of structural and numerical chromosome abnormalities in sperm of healthy men by multicolor fluorescence in situ hybridization. Hum Genet. 1996;98:608–15.View ArticlePubMedGoogle Scholar
- Pollier J, Rombauts S, Goossens A. Analysis of RNA-Seq data with TopHat and Cufflinks for genome-wide expression analysis of jasmonate-treated plants and plant cultures. Methods Mol Biol. 2013;1011:305–15.View ArticlePubMedGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.View ArticlePubMedPubMed CentralGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.View ArticlePubMedPubMed CentralGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2014;2014:1–4.Google Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Cold Spring Harbor: BioRxIV; 2014.Google Scholar
- R Core Team: R: a language and environment for statistical computing.; 2013.Google Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Parikshak NN, Swarup V, Belgard TG, Irimia M, Ramaswami G, Gandal MJ, Hartl C, Leppa V, Ubieta LT, Huang J, et al. Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism. Nature. 2016;540:423–7.View ArticlePubMedGoogle Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.View ArticlePubMedPubMed CentralGoogle Scholar
- Lucas ME, Muller F, Rudiger R, Henion PD, Rohrer H. The bHLH transcription factor hand2 is essential for noradrenergic differentiation of sympathetic neurons. Development. 2006;133:4015–24.View ArticlePubMedGoogle Scholar
- Morikawa Y, D'Autreaux F, Gershon MD, Cserjesi P. Hand2 determines the noradrenergic phenotype in the mouse sympathetic nervous system. Dev Biol. 2007;307:114–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmidt M, Lin S, Pape M, Ernsberger U, Stanke M, Kobayashi K, Howard MJ, Rohrer H. The bHLH transcription factor Hand2 is essential for the maintenance of noradrenergic properties in differentiated sympathetic neurons. Dev Biol. 2009;329:191–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Gossan NC, Zhang F, Guo B, Jin D, Yoshitane H, Yao A, Glossop N, Zhang YQ, Fukada Y, Meng QJ. The E3 ubiquitin ligase UBE3A is an integral component of the molecular circadian clock through regulating the BMAL1 transcription factor. Nucleic Acids Res. 2014;42:5765–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu Y, Bolduc FV, Bell K, Tully T, Fang Y, Sehgal A, Fischer JA. A Drosophila model for Angelman syndrome. Proc Natl Acad Sci U S A. 2008;105:12399–404.View ArticlePubMedPubMed CentralGoogle Scholar
- Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, Mill J, Cantor RM, Blencowe BJ, Geschwind DH. Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011;474:380–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Peters SU, Horowitz L, Barbieri-Welge R, Taylor JL, Hundley RJ. Longitudinal follow-up of autism spectrum features and sensory behaviors in Angelman syndrome by deletion class. J Child Psychol Psychiatry. 2012;53:152–9.View ArticlePubMedGoogle Scholar
- Scoles HA, Urraca N, Chadwick SW, Reiter LT, Lasalle JM. Increased copy number for methylated maternal 15q duplications leads to changes in gene and protein expression in human cortical samples. Molecular autism. 2011;2:19.View ArticlePubMedPubMed CentralGoogle Scholar
- Kuhnle S, Kogel U, Glockzin S, Marquardt A, Ciechanover A, Matentzoglu K, Scheffner M. Physical and functional interaction of the HECT ubiquitin-protein ligases E6AP and HERC2. J Biol Chem. 2011;286:19410–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Martinez-Noel G, Galligan JT, Sowa ME, Arndt V, Overton TM, Harper JW, Howley PM. Identification and proteomic analysis of distinct UBE3A/E6AP protein complexes. Mol Cell Biol. 2012;32:3095–106.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamasaki K, Joh K, Ohta T, Masuzaki H, Ishimaru T, Mukai T, Niikawa N, Ogawa M, Wagstaff J, Kishino T. Neurons but not glial cells show reciprocal imprinting of sense and antisense transcripts of Ube3a. Hum Mol Genet. 2003;12:837–47.View ArticlePubMedGoogle Scholar
- Griesi-Oliveira K, Acab A, Gupta AR, Sunaga DY, Chailangkarn T, Nicol X, Nunez Y, Walker MF, Murdoch JD, Sanders SJ, et al. Modeling non-syndromic autism and the impact of TRPC6 disruption in human neurons. Mol Psychiatry. 2014;20(11):1350–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Moncla A, Malzac P, Livet MO, Voelckel MA, Mancini J, Delaroziere JC, Philip N, Mattei JF. Angelman syndrome resulting from UBE3A mutations in 14 patients from eight families: clinical manifestations and genetic counselling. J Med Genet. 1999;36:554–60.PubMedPubMed CentralGoogle Scholar
- Rapakko K, Kokkonen H, Leisti J. UBE3A gene mutations in Finnish Angelman syndrome patients detected by conformation sensitive gel electrophoresis. Am J Med Genet A. 2004;126A:248–52.View ArticlePubMedGoogle Scholar
- Gustin RM, Bichell TJ, Bubser M, Daily J, Filonova I, Mrelashvili D, Deutch AY, Colbran RJ, Weeber EJ, Haas KF. Tissue-specific variation of Ube3a protein expression in rodents and in a mouse model of Angelman syndrome. Neurobiol Dis. 2010;39:283–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Dindot SV, Antalffy BA, Bhattacharjee MB, Beaudet AL. The Angelman syndrome ubiquitin ligase localizes to the synapse and nucleus, and maternal deficiency results in abnormal dendritic spine morphology. Hum Mol Genet. 2008;17:111–8.View ArticlePubMedGoogle Scholar
- Harlalka GV, Baple EL, Cross H, Kuhnle S, Cubillos-Rojas M, Matentzoglu K, Patton MA, Wagner K, Coblentz R, Ford DL, et al. Mutation of HERC2 causes developmental delay with Angelman-like features. J Med Genet. 2013;50:65–73.View ArticlePubMedGoogle Scholar
- Puffenberger EG, Jinks RN, Wang H, Xin B, Fiorentini C, Sherman EA, Degrazio D, Shaw C, Sougnez C, Cibulskis K, et al. A homozygous missense mutation in HERC2 associated with global developmental delay and autism spectrum disorder. Hum Mutat. 2012;33:1639–46.View ArticlePubMedGoogle Scholar
- Sarachana T, Hu VW. Genome-wide identification of transcriptional targets of RORA reveals direct regulation of multiple genes associated with autism spectrum disorder. Mol Autism. 2013;4:14.View ArticlePubMedPubMed CentralGoogle Scholar
- Jensen L, Farook MF, Reiter LT. Proteomic profiling in Drosophila reveals potential Dube3a regulation of the actin cytoskeleton and neuronal homeostasis. PLoS One. 2013;8:e61952.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee SY, Ramirez J, Franco M, Lectez B, Gonzalez M, Barrio R, Mayor U. Ube3a, the E3 ubiquitin ligase causing Angelman syndrome and linked to autism, regulates protein homeostasis through the proteasomal shuttle Rpn10. Cell Mol Life Sci. 2014;71:2747–58.View ArticlePubMedGoogle Scholar