- Research article
- Open Access
Altered neuronal physiology, development, and function associated with a common chromosome 15 duplication involving CHRNA7
BMC Biology volume 19, Article number: 147 (2021)
Copy number variants (CNVs) linked to genes involved in nervous system development or function are often associated with neuropsychiatric disease. While CNVs involving deletions generally cause severe and highly penetrant patient phenotypes, CNVs leading to duplications tend instead to exhibit widely variable and less penetrant phenotypic expressivity among affected individuals. CNVs located on chromosome 15q13.3 affecting the alpha-7 nicotinic acetylcholine receptor subunit (CHRNA7) gene contribute to multiple neuropsychiatric disorders with highly variable penetrance. However, the basis of such differential penetrance remains uncharacterized. Here, we generated induced pluripotent stem cell (iPSC) models from first-degree relatives with a 15q13.3 duplication and analyzed their cellular phenotypes to uncover a basis for the dissimilar phenotypic expressivity.
The first-degree relatives studied included a boy with autism and emotional dysregulation (the affected proband-AP) and his clinically unaffected mother (UM), with comparison to unrelated control models lacking this duplication. Potential contributors to neuropsychiatric impairment were modeled in iPSC-derived cortical excitatory and inhibitory neurons. The AP-derived model uniquely exhibited disruptions of cellular physiology and neurodevelopment not observed in either the UM or unrelated controls. These included enhanced neural progenitor proliferation but impaired neuronal differentiation, maturation, and migration, and increased endoplasmic reticulum (ER) stress. Both the neuronal migration deficit and elevated ER stress could be selectively rescued by different pharmacologic agents. Neuronal gene expression was also dysregulated in the AP, including reduced expression of genes related to behavior, psychological disorders, neuritogenesis, neuronal migration, and Wnt, axonal guidance, and GABA receptor signaling. The UM model instead exhibited upregulated expression of genes in many of these same pathways, suggesting that molecular compensation could have contributed to the lack of neurodevelopmental phenotypes in this model. However, both AP- and UM-derived neurons exhibited shared alterations of neuronal function, including increased action potential firing and elevated cholinergic activity, consistent with increased homomeric CHRNA7 channel activity.
These data define both diagnosis-associated cellular phenotypes and shared functional anomalies related to CHRNA7 duplication that may contribute to variable phenotypic penetrance in individuals with 15q13.3 duplication. The capacity for pharmacological agents to rescue some neurodevelopmental anomalies associated with diagnosis suggests avenues for intervention for carriers of this duplication and other CNVs that cause related disorders.
Reciprocal copy number variants (CNVs) related to neurodevelopmental and neuropsychiatric disorders result from single-copy deletion or duplication of susceptible genomic intervals [1,2,3]. While CNVs involving deletion generally cause severe, highly penetrant patient phenotypes, CNVs involving duplication often exhibit widely variable and less penetrant phenotypic expressivity among affected individuals [4,5,6]. Among these, either chromosome 15q13.3 deletion or duplication causes clinical phenotypes including autism spectrum disorder (ASD), intellectual disability (ID), mood disorder, language delay, or schizophrenia [7,8,9]. 15q13.3 microdeletions usually cause severe cognitive deficits, behavioral abnormalities, and highly penetrant ASD [7, 10]. By comparison, 15q13.3 microduplication often causes milder phenotypes, including borderline ID, ASD, and attention deficit hyperactivity disorder (ADHD) [7, 11]. Notably, 15q13.3 microduplication is present in 1.25% of reported ADHD probands but also 0.61% of control subjects, suggesting that many individuals tolerate 15q13.3 duplication without a clinical diagnosis, due to its poor phenotypic penetrance [7, 8, 11].
15q13.3 microduplications either involve only the α-7 nicotinic acetylcholine receptor subunit gene (CHRNA7), with or without the first exon of OTUD7A, or instead duplicate multiple genes [7, 12,13,14,15]. CHRNA7 is a prototypical genetic contributor to complex psychiatric disease: this protein family forms ligand-gated ion channels, which are stimulated by choline and acetylcholine to trigger calcium, sodium, and potassium cation flux . CHRNA7 duplication causes motor delays, hypotonia, ASD, ID, schizophrenia, and epilepsy, with particular phenotypes varying by individual [7, 15,16,17]. Although the clinical significance of increased CHRNA7 receptor dosage is unclear, no other copy number variation is detected in most reported cases, indicating that CHRNA7 duplication is pathogenic, but with reduced and highly variable penetrance [7, 8, 16].
Patient-derived induced pluripotent stem cell (iPSC) models provide a powerful approach for modeling neurodevelopmental disorders. In recent years, we and others have characterized phenotypes associated with clinical diagnoses in these models, including modeling syndromic and de novo cases of ASD and monogenic or polygenic contributors to disease [18,19,20,21,22]. These studies identified diagnosis-associated phenotypes, including gene expression changes, differential regulation of developmental signaling, and altered neurogenesis and synaptogenesis [19,20,21, 23,24,25,26]. In addition to identifying these diagnosis-associated phenotypes, some studies found specific targets amenable to pharmacological rescue [18, 27, 28].
While CNVs at 15q13.3 involving CHRNA7 duplication contribute to many neuropsychiatric disorders, consequences of these CNVs have not been extensively modeled. Mice with CHRNA7 knockout, mimicking some 15q13.3 deletions, did not exhibit behavioral phenotypes [10, 11, 29]. CHRNA7 duplication has not been modeled in vivo and the size of 15q13.3 duplications precludes their modeling by genome engineering in rodents. CHRNA7 overexpression in a mouse neuroblastoma cell line altered receptor sensitivity to choline and varenicline . Two studies have characterized human iPSC-derived models involving duplications in the chromosome 15q region [28, 31]. One focused on a large duplication at 15q11-q13.1 involving 33 genes; a maternally expressed gene (UBE3A) was upregulated in this interval and this effect could be pharmacologically rescued . A second study utilized three iPSC-derived excitatory neuron progenitor cell (NPC) models each from subjects with either 15q13.3 duplication or deletion . These models exhibited increased or decreased CHRNA7 expression, respectively, while both duplication and deletion reduced CHRNA7-dependent calcium flux, indicating diminished channel activity. In the duplication models, diminished calcium flux was linked to elevated expression of several chaperone genes that control CHRNA7 protein folding and trafficking through the endoplasmic reticulum (ER) to the cell surface. Based upon these data, the authors proposed a model in which increased CHRNA7 levels impaired the efficiency of CHRNA7 trafficking through the ER, resulting in lower levels of CHRNA7 protein reaching the cell membrane, which diminished CHRNA7 channel activity .
CHRNA7 duplications are common contributors to psychopathophysiology, exhibit widely variable phenotypic penetrance, and are difficult to model in animals. While there is some evidence for altered physiology of NPCs carrying CHRNA7 duplications , potential contributors to the variable phenotypic penetrance among 15q13.3 duplication carriers are entirely unknown. How this CNV more broadly affects neurodevelopment, global gene expression, and electrophysiological function of neurons has also not been evaluated. Therefore, here we generated iPSC models from two first-degree relatives with the same 15q13.3 copy number variant, which duplicates only CHRNA7. These subjects include a boy with distinct features of autism and emotional dysregulation (the affected proband, AP) and his clinically unaffected mother (the UM). These models were compared to unrelated male and female control subjects lacking this duplication (the UC-M and UC-F). We used these models to assess the consequences of CHRNA7 duplication on the development and function of both cortical excitatory neurons (cExN) and inhibitory interneurons (cIN), as disruption of both neuronal cell types frequently contributes to neurodevelopmental disorders. This work defined a suite of phenotypic penetrance-associated neurodevelopmental anomalies in the AP, including deficits in neurite extension, neuronal migration, neuronal maturation, and ER stress, and related perturbations of neuronal gene expression. While these neurodevelopmental deficits were not present in the other model with CHRNA7 duplication (UM), both models carrying the duplication (UM and AP) shared multiple functional abnormalities, as defined by electrophysiology, including increased action potential firing and enhanced choline responsiveness.
To our knowledge, this is the first cellular modeling study to characterize cases of CHRNA7 duplication involving two first-degree relatives with differential diagnoses. A range of molecular, cellular, and functional assays was used to define alterations of neuronal gene expression, neurodevelopment, and altered functional properties; these were linked either to phenotypic penetrance or to CHRNA7 duplication, both with and without diagnosis, and some diagnosis-associated phenotypes were amenable to pharmacological rescue. The distinct neuronal anomalies related to clinical phenotype versus CHRNA7 duplication that were defined here could contribute to the differential phenotypic penetrance seen in other duplication carriers and provide models for pre-clinical testing of potential interventions for individuals with this disorder.
iPSC-derived neural progenitor cells from the affected proband exhibit increased proliferation
The selected pedigree of four individuals includes three subjects with the same 15q13.3 duplication, the mother, who has no clinical diagnosis (UM), her older son, who exhibits distinct features of autism and emotional dysregulation (the affected proband, AP), and his younger, affected brother, who exhibits mild ASD, ADHD, and mood disorder traits. The family pedigree and the clinical phenotypes of these three mutation carriers are summarized in Table 1 and Additional file 1. Chromosomal Microarray Analysis (CMA) defined one CNV in these three subjects, a ~400 kilobase gain at chromosome 15, band q13.3 (see the “Methods” section, Table 1, Additional file 1). Only one gene, CHRNA7, is in the duplicated region. The father in this pedigree does not carry this duplication (Additional file 1).
Renal epithelial cells from the UM and AP were reprogrammed to generate three clonal iPSC lines per subject and compared with single clonal iPSC lines derived from unrelated, unaffected male and female control subjects (UC-M and UC-F). All lines exhibited similar morphology, expressed the pluripotency marker POU5F1, and were karyotypically normal (Additional file 2). The development and/or function of cortical excitatory neurons (cExNs) and inhibitory cortical interneurons (cINs) is frequently disrupted to contribute to neurodevelopmental disorders [32, 33]. Therefore, we used two clonal lines each from the UM and AP, and one clonal line each from the UC-M and UC-F, to generate cExN and cIN neural progenitor cells (cExNPCs and cINPCs) (Fig. 1a). After 15 days of NPC specification, cExNPCs and cINPCs were maintained as a monolayer. CHRNA7 expression levels, assessed by RT-qPCR, were significantly increased in both AP- and UM-derived NPCs, relative to UC-M or UC-F-derived NPCs (Fig. 1b). During cExNPC and cINPC maintenance, AP-derived NPCs exhibited more rapid proliferation than UM, UC-M, or UC-F-derived NPCs, as measured by seeding equal cell numbers and quantitation after four days (Fig. 1c, e). We further assessed cell cycle kinetics by FACS analysis of propidium iodide (PI)-stained cExNPCs and cINPCs, finding that AP-derived NPCs had significantly higher percentages of both S- and M-phase cells than one or more of the unrelated control NPCs. Percentages of cells in S and 4N phases of the cell cycle are shown for cExNPCs and cINPCs, respectively (Fig. 1d, f), with summary data for all cell cycle phases in Additional file 3, and biological replicate experiments and clones used for each experiment in this manuscript summarized in Additional file 4.
Defects in neurite extension and production of VGAT-expressing punctae in AP-derived neurons
To assess the differentiation potential of the UC-M, UM, and AP lines, cExNPCs and cINPCs were combined at a 1:1 ratio to generate cortical neuroids (Fig. 2a). This co-culture accelerates neuronal differentiation and maturation and models interactions between cortical excitatory and inhibitory neurons that occur during cortical development in vivo. Upon dissociation and plating of these day 15 neuroids, AP-derived neuroids exhibited increased expression of NPC markers DLX2, TBR1, and PAX6, versus UM and UC neuroids, and also had higher frequencies of cells expressing the proliferative marker Ki67 (Fig. 2b–f). This could relate to increased AP-derived NPC proliferation (Fig. 1). We also plated differentiating neuroids without dissociation, allowing them to extend neurites. By 5 days, AP-derived neuroids exhibited a neurite extension deficit not seen in either UM- or UC-derived neuroids. This was quantified using light microscopy and further visualized by MAP2 staining (Fig. 3a–c). On differentiation day 15, neuroids were also dissociated and plated neurons assessed by immunocytochemistry (ICC) for MAP2, confirming neurite length deficits in AP-derived neurons (Fig. 3d–d’). ICC for GABA and glutamate transporters, VGAT and VGLUT, further assessed synaptic vesicle formation and transport in cINs and cExNs, respectively. While VGLUT-expressing punctae were present in cExN neurites derived from all three lines, AP-derived cINs exhibited significantly reduced formation of VGAT-expressing punctae along neurites, versus UM and UC-M models (Fig. 3e–e”). Both increased expression of NPC and proliferative markers and impaired acquisition of characteristics of differentiated and mature neurons suggested that the AP-derived neurons were relatively immature, relative to UM-, UC-F-, and UC-M-derived neurons.
Transcriptome comparisons in cortical neuroids
To identify gene expression differences distinguishing neurons from the four models, we performed RNA-seq analysis on cortical neuroids after 15 days of differentiation. Four independent biological replicates were generated and clustered by Principal Component Analysis (PCA) of processed reads (Fig. 4a). We defined differentially expressed genes (DEGs) by pairwise comparisons between sample types (FDR <0.05 and log2 fold-difference >1 cut-off values; see Methods, Additional file 5). We focused first on genes that were differentially expressed in the AP versus (vs.) the UM, as these subjects have a shared genetic background including CHRNA7 duplication, but exhibit differential phenotypes (Fig. 4b, blue). Hierarchical clustering of these DEGs, with the UC samples included for comparison, demonstrated that most exhibited increased expression in the AP, versus the UM (Fig. 4c). We used Ingenuity Pathway Analysis (IPA) to define enriched pathways and disease-related networks for AP-specific DEGs (Additional file 6). These included molecular pathways involved in axon guidance signaling, the cell cycle, Wnt signaling, GABA receptor activity, neuroinflammation, and gap junction signaling (Fig. 4d; Additional file 6), with most Wnt signaling-related genes exhibiting diminished expression in AP- versus UM-derived neurons (Fig. 4e). IPA disease network analysis also defined clusters of DEGs related to nervous system development, developmental disorders, behavior, and psychological disorders (Fig. 4f). For example, genes related to cognition and neuritogenesis exhibited predominantly reduced expression in AP-derived neurons (Fig. 4g-h). This was congruent with AP-specific neurite extension deficits (Fig. 3a).
We also curated AP-specific DEGs by comparison with both UM-, UC-M-, and UC-F-derived neurons (1052 genes; Additional file 7A). These included many genes with diminished expression in the AP relative to all other models (Additional file 7B) and encompassed many molecular pathways and disease-related GO terms similar to those identified in AP versus UM comparisons, including Wnt, axon guidance, and GABA receptor signaling, nervous system development, and psychological disorders (Additional file 7C-F; Additional file 8). For example, gene clusters related to neurite growth and neuronal migration had predominantly reduced expression in the AP versus other models (Additional file 7C-F; Additional file 8), suggesting that similar pathways exhibited impaired expression in the AP, by comparison with neurons from both related (UM) and unrelated (UC) models.
As another way to identify specific pathways that could be contributing to the AP's differential affectation, we cross-referenced all genes with differential expression between the AP and the three other models (UM, UC-M, UC-F) with genes that have genome-wide significance for de novo protein-coding variants associated with neurodevelopmental disorders [34, 35] or that were significant in a transcription-wide association study (TWAS) of autism . The intersecting genes included ANXA1, BCL11A, BCL11B, CACNA1A, COL2A1, DSCAM, ERF, FBN1, FOXP2, GATA3, KCNJ6, KCNS3, KIF11, MEIS2, NFIA, PRKAR1B, SATB1, SLC6A1, SPRY2, TBR1, TCF4, TMEM42, TRPM3, and ZEB2 (Additional file 8). To look at potential genetic interactions, physical interactions, co-expression, shared protein domains, and co-localization we entered this list into GeneMania , which defined a network of dysregulated DEGs (Additional file 7G). We also tested whether there was specifically an enrichment of protein-protein interactions (PPI) using STRING-DB . In this PPI analysis, we observed 15 edges and expected 4 edges which showed an enrichment (p = 6.62 x 10-6). Gene ontology analysis of this network revealed enrichments for molecular functions including DNA-binding transcription repressor activity (FDR = 1.0 x 10-5) and voltage-gated cation channel activity (FDR = 6.9 x 10-3), suggesting that the genes in the categories that are shown in Additional file 7G may be contributors to the AP’s clinical phenotype.
UM-specific differential gene expression in cortical neuroids compared to unrelated controls
Since both related models (AP/UM) carry the same CHRNA7 duplication, we also defined UM-specific DEGs, by comparison with unrelated controls (Additional file 9A). A DEG cluster exhibited increased expression, versus all other samples, while UM- and AP-derived neurons also shared a cluster with decreased expression versus the UC samples (Additional file 9B). UM-specific DEGs exhibited enrichment for some of the same gene ontology (GO) terms identified for AP-specific DEGs (e.g., GABA receptor and gap junction signaling, behavior/psychological disorders; Additional file 9C-D; Additional file 10). However, many DEGs and enriched GO terms and pathways obtained from UC model comparisons differed between the UM and AP. For example, Wnt signaling-related genes were not differentially expressed in the UM versus UC comparisons. Furthermore, while UM-specific DEGs also included a neuritogenesis-related network, genes in this network were more highly expressed in the UM versus the UC-M and UC-F models (Additional file 9E-F), while neuritogenesis-related genes had diminished expression in the AP, versus both UM and UC models (Fig. 4h, Additional file 7F). These DEGs included cytoskeletal genes involved in neuritogenesis, including ACTG2, ACTN1, and/or ACTN2 [39,40,41,42,43]; these genes exhibited diminished expression in the AP versus the UM and UC controls, but elevated expression in the UM relative to the unrelated controls. We used cortical neuroids derived from a second clonal line for the AP and UM models to confirm a subset of these findings by RT-qPCR, obtaining findings congruent with the RNA-seq analysis (Additional file 11A-B).
To assess the relative contributions that sex, age, or the individual genetic backgrounds of the subjects used to derive these models may have had on the identification of DEGs in this work, we also performed variancePartition analysis. This analysis indicated that the individual from whom the samples were derived was the major contributor to differential expression, while the age and sex of study subjects were minor contributors (Additional file 11C). We further characterized the contribution that sex of the model may have had on the identification of DEGs by comparing each set of DEGs (Additional file 5) with genes that were observed to be sex-differentially expressed in the human brain in a recent study , finding that fewer than 1% of the DEGs identified here were previously characterized as exhibiting sex-differential expression (Additional file 5, right columns). This is congruent with the covariance analysis and suggests that the sex of the samples was not a major contributor to the DEGs identified.
Interneuron migration is diminished in AP-derived neurons and this defect is partially reversed by a Wnt agonist
AP-derived neurons (versus the UM) exhibited diminished expression of genes that regulate neuronal migration (Fig. 5a), suggesting that they could have impaired migration. In vivo, cortical interneurons undergo tangential migration from the ventral telencephalon to the cortex. Therefore, we adapted an organoid-based neuronal migration assay (see the “Methods” section): a synapsin-GFP expressing cExN neuroid and synapsin-RFP expressing cIN neuroid were apposed, and neuronal migration from one to the other was evaluated (Fig. 5b). cINs from AP neuroids exhibited reduced migration by comparison with the other models (Fig. 5c-d). By contrast, the migration of AP-derived cExNs was only modestly impaired, while UM-derived cExNs exhibited diminished migration versus the control models, suggesting some impairment of migration in both CHRNA7 duplication carriers (Fig. 5c, e). As Wnt signaling has been linked to neuronal migration , we hypothesized that diminished expression of Wnt signaling-related genes (Fig. 4e) could contribute to the AP model’s reduced cIN migration. Assessment of the AP vs. UM DEGs with the Percayai CompBio tool likewise revealed a network of differentially expressed genes linked to both Wnt signaling and neuronal migration (Additional file 6C-D), some of which (e.g., GJA1/TGFBR2) were also identified as Wnt- and neuronal migration-related DEGs in Ingenuity Pathway Analysis (Figs. 4e, 5a). These data suggested that deficits in Wnt signaling could contribute to reduced migration of the AP-derived neurons. Indeed, treatment of AP neuroids with the Wnt agonist CHIR-99021 partially reversed their impaired cIN migration (Fig. 5c-d). However, as a number of other pathways implicated in neuronal migration were also dysregulated in the AP, including axonal guidance, actin cytoskeleton, integrin, and netrin signaling (e.g., Additional file 6), this dysregulation could also have contributed to the deficits in neuronal migration observed in the AP model.
Increased ER stress selective to the AP-derived NPCs is rescued by the ryanodine receptor antagonist JTV-519
In prior work, microduplication at chromosome 15q13.3 increased CHRNA7 expression and caused moderately increased expression of several endoplasmic reticulum (ER) chaperone and unfolded protein response (UPR)-related ER stress markers . We assessed whether our models likewise exhibited altered expression of ER chaperones in neuroids, but did not observe increased expression of these or other ER chaperone or UPR-related genes in AP-derived samples versus the UM and UC-M models (Fig. 6a; not obtained as DEGs in Additional file 5). However, the ryanodine receptor RYR3, which is most highly expressed in the brain, exhibited increased expression in the AP versus the other samples (Fig. 6b). As ryanodine receptors modulate calcium homeostasis following increased ER stress, we hypothesized that the AP model’s elevated ER stress could involve altered calcium homeostasis, rather than elevated ER chaperone or UPR pathway activities. To assess ER stress in these models, we introduced an expression construct encoding a Secreted ER Calcium Monitoring Protein (SERCaMP)-luciferase stress sensor , which monitors stress-related calcium release from the ER, into cINPCs (see the “Methods” section). Interestingly, the AP’s cINPCs selectively exhibited increased ER stress, while the UM did not differ from the UC models (Fig. 6c). As ER stress both activates the UPR and triggers calcium release through ryanodine receptors, we tested whether chemical antagonists of either process could rescue the AP model’s elevated ER stress. While neither UPR antagonist (PBA/TUDCA) exhibited rescue activity, one of two ryanodine receptor antagonists tested (JTV-519, but not Dantrolene) suppressed the AP’s elevated ER stress (Fig. 6d).
Electrophysiological characterization of cortical neurons derived from the AP, UM, and UC-M models
To assess whether these models had altered neuronal function, synapsin-GFP-labeled cExNPCs and synapsin–RFP-labeled cINPCs were differentiated as neuroid co-cultures, then were dissociated and grown for 20 days on rat cortical astrocyte feeders to promote neuronal maturation (see Methods; Additional file 12A). Whole-cell voltage and current-clamp recordings were obtained separately from cExNs and cINs from the AP, UM, and UC-M models (18–27 cells assessed/group). For a complete data summary, see Additional file 13. Visual inspection suggested that UM-derived neurons were consistently larger than AP- or UC-M-derived neurons; this was confirmed by cell soma diameter measurements (Fig. 7a; Additional file 12B). This was associated with higher UM cell capacitance (proportional to surface area) and lower input resistance (Fig. 7b). Voltage steps from −80 to 0 mV evoked a fast peak of inward current mediated by tetrodotoxin-sensitive sodium channels, followed by steady-state outward current mediated by voltage-gated potassium channels. Comparisons between models revealed that both the UM and AP had lower outward current density than the UC-M (Fig. 7c). A 2-way ANOVA analysis among all groups also indicated that cINs had a slightly higher input resistance than cExNs (Fig. 7d). Interestingly, although separate recordings were generated and analyzed for cExNs and cINs across the three models, most significant differences were commonly observed in both cExNs and cINs. Therefore, we combined cIN and cExN recording data points for results in Figs. 7 and 8.
Because CHRNA7 encodes the α-7 nicotinic acetylcholine (ACh) receptor subunit , we measured integrated current elicited by brief exposure to 500μM ACh in a subset of cells that were characterized by whole-cell current and voltage clamp. Most of these cells were also tested with 500μM choline, a relatively selective agonist for homomeric α7 receptors  with much weaker activity at heteromeric α7β2 receptors or other neuronally expressed nicotinic receptors [48, 49]. Interestingly, UC-M neurons exhibited a significantly larger response to ACh than AP cells (Fig. 7e-f, note log scale in F) and a larger ACh to choline ratio compared to both UM and AP neurons (Fig. 7g). These data suggest that homomeric α7 receptors may mediate a higher proportion of the total ACh-evoked current in both UM and AP neurons, relative to UC-M neurons.
Current clamp recordings (Fig. 8) revealed additional functional alterations shared by UM- and AP-derived neurons. UM and AP cells fired a larger numbers of action potentials on average than UC-M cells, for current injections up to 120 pA (Fig. 8a, b), and displayed a less depolarized threshold for action potential initiation (Fig. 8c). There was also a trend for lower initial spike frequency that reached significance for the UM model (Fig. 8d). In addition to these shared alterations, UM and AP neurons exhibited several functional differences: UM neurons fired a larger maximal number of spikes during an 800 msec depolarizing pulse than UC-M or AP neurons (Fig. 8e) and also exhibited a larger peak amplitude and shorter half width to the first spike than UC-M or AP neurons (Fig. 8f-g). AP neurons uniquely exhibited a substantially lower rheobase (the minimal depolarizing current required to reach threshold; Fig. 8h) and a more significant decline in average action potential peak amplitude with each succeeding spike (Fig. 8a, i), versus both UM and UC-M models. Therefore, by contrast with AP-specific neurodevelopmental phenotypes, the AP and UM models shared a number of functional anomalies not seen in the UC-M model, while also exhibiting some unique characteristics.
Spontaneous excitatory and inhibitory synaptic currents (Additional file 12C) were detected in recordings from all three models; however, only a subset of recorded cells displayed synaptic events and the frequency was very low in all but a few cells. This result demonstrates that, while both cExNs and cINs of all three genotypes are able to form functional connections, the number of fully functional contacts is likely to be low, even after 20+ days of maturation on astrocyte feeder cells.
Cellular correlates of differential phenotypic penetrance among carriers of CHRNA7 duplication
Like many CNVs involving microduplication, duplications at 15q13.3 cause multiple psychiatric phenotypes, with patients exhibiting widely varied clinical phenotypic penetrance. How duplication of CHRNA7 alone versus other genes in this interval causes these phenotypes, and the basis of their variable penetrance, were largely uncharacterized. Therefore, here we derived iPSC models from a pedigree where multiple individuals carried the same CNV, involving only CHRNA7 duplication, but had differential diagnoses. We compared clonal iPSC models from the UM, who has no clinical phenotypes, her son (AP), who has multiple clinical neurodevelopmental phenotypes, and controls lacking this duplication, using cExN and cIN derivation to examine neurodevelopment and function. This study design enabled the identification of potential contributors to differential phenotypic penetrance between models with the same CHRNA7 duplication on a partially shared genetic background. While both the AP and UM had elevated CHRNA7 expression, the AP model exhibited many neurodevelopmental phenotypes not seen in the UM, UC-M, or UC-F models. These included reduced neurite extension and length, diminished cIN neuronal migration and formation of VGAT-expressing synaptic punctae, elevated ER stress, and correspondingly reduced gene expression in related molecular pathways and processes. Intriguingly, many pathways with reduced expression in the AP instead exhibited increased expression in the UM model, versus both the AP and controls. This elevated molecular signaling could have compensated for some consequences of the genetic liability of CHRNA7 duplication, enabling the UM to overcome neurodevelopmental impairment that likely contributed to the AP’s clinical phenotypic penetrance. However, electrophysiological analyses of cExNs and cINs matured on astrocyte feeder layers revealed functional anomalies shared by AP- and UM-derived neurons versus controls, including increased action potential firing and enhanced choline responsiveness. These functional changes are consistent with elevated CHRNA7 channel activity and altered ligand sensitivity, and indicate common functional anomalies likely conferred by CHRNA7 duplication in both AP- and UM-derived neurons; these, if also manifested in vivo, are nonetheless insufficient to trigger clinical phenotypes in the UM.
This study generally illustrates that patient-derived neuronal models can define cellular and molecular anomalies associated with, and potentially contributing to, differential clinical phenotypic penetrance. Although many CNVs involving microduplication likewise exhibit highly variable phenotypic penetrance, the basis of this phenomenon is not well understood. Here, we distinguished both contributors to differential phenotypic penetrance and those common to multiple carriers of the same genetic liability, regardless of the clinical phenotype. We could also identify changes in molecular signaling that may have compensated for this genetic liability, contributing to a lack of clinical phenotype in the UM. This work indicates that patient-derived iPSC models may have general utility for linking both genetic liabilities and differing phenotypic penetrance stemming from those liabilities to specific cellular, molecular, or functional alterations that may be prior or ongoing contributors to clinical phenotype. These phenotypes then provide a basis and experimental platform for identifying molecular and pharmacological agents that can rescue diagnosis-associated phenotypes, as we did in identifying distinct pharmacological agents that could rescue both the AP model’s ER stress and impaired cIN migration.
Shared functional alterations among models with CHRNA7 duplication
Functional assessment of neurons differentiated from these models revealed a number of electrophysiological abnormalities common to both individuals with CHRNA7 duplication (AP and UM), regardless of clinical diagnosis. The CHRNA7 receptor forms either homomeric or heteromeric pentamers; these exhibit differential responsiveness, with homomeric α7 receptors having selective choline sensitivity, while heteromeric receptors are more ACh responsive . Here, both AP- and UM-derived neurons exhibited elevated choline responsiveness and diminished ACh responsiveness, versus controls. This is congruent with the possibility that CHRNA7 duplication, which increased CHRNA7 expression in the AP and UM, may increase the production of homomeric α7 channels, altering electrophysiological function and ligand sensitivity. Other functional abnormalities were also shared by UM- and AP-derived neurons, versus controls, including a reduced outward potassium (K) current and increased number of action potentials. Together, these results suggest that increased channel activity may enhance entry of calcium, sodium, and potassium ions into neurons, increasing action potential firing in both models with 15q13.3 duplication.
Altered developmental mechanisms that track with the AP model’s severe phenotype
Congruent with a prior study , we observed elevated ER stress in iPSC-derived NPCs carrying CHRNA7 duplication in the AP but not the UM model. However, the ER stress phenotypes observed in these two studies differ substantially. Nicotinic acetylcholine receptor trafficking in cells involves assembling five subunits into a pentameric receptor in the endoplasmic reticulum (ER), followed by trafficking to the plasma membrane [50, 51]. In Gillentine et al., elevated ER stress was attributed to the unfolded protein response (UPR), based upon increased expression of several ER chaperone- and UPR-related marker genes . By contrast, our mRNA expression analysis did not reveal altered expression of these or other UPR-related markers in the AP. However, altered calcium homeostasis is another aspect of ER stress : to maintain homeostasis, the ER releases calcium via several mechanisms, one of which involves ryanodine receptor signaling activation in response to altered cytoplasmic calcium, through calcium-induced calcium release (CICR) . Our analysis revealed elevated levels of both calcium receptors CACNA1A, CACNA1B, and CACNA2D2, and RYR3, a ryanodine receptor with brain-enriched expression , in the AP, suggesting potentially altered calcium homeostasis. Accordingly, a reporter assay detected elevated ER stress-linked CICR in the AP, versus the other models, and the ryanodine receptor antagonist JTV-519 selectively reduced this response to baseline conditions, while UPR antagonists did not. These results further implicate altered calcium homeostasis rather than UPR pathway activation as the major contributor to the AP model’s elevated ER stress here. Another major difference between these findings is that the NPCs in Gillentine et al. exhibited characteristics consistent with reduced CHRNA7 protein trafficking to the cell membrane and reduced CHRNA7 channel activity. By contrast, when we conducted electrophysiology on neurons from our models after maturation on astrocyte feeder layers, both CHRNA7 duplication carriers (AP and UM) displayed hallmarks of increased CHRNA7 channel activity, including increased action potential firing, and altered ligand sensitivity. Therefore, although both studies observed elevated ER stress in models with CHRNA7 duplication, here reporter assays and pharmacological rescue demonstrated that, for AP-derived models, this involved altered calcium homeostasis rather than the UPR, and that CHRNA7 duplication in the AP and UM increased, rather than diminished, channel activity. It will be interesting to assess these phenotypes in additional models with CHRNA7 duplication in future work, to determine whether effects on ER stress and channel function indeed vary between models derived from different subjects; these could relate to other aspects of particular models, such as interactions with the patient’s genetic background, specific clinical phenotypes, or degree of phenotypic penetrance.
In this study, the AP model exhibited many neurodevelopmental anomalies, including increased NPC proliferation but reduced neurite outgrowth and length, and concordant dysregulation of genes mediating axonal guidance, gap junctional connectivity, and neuritogenesis, with reduced expression of most genes in these pathways in the AP versus the UM and controls (Additional file 14 [54,55,56,57,58,59,60,61,62,63,64,65];). Altered NPC proliferation, neuronal differentiation, neurite extension, neuronal maturation, and synapse formation have also been observed in other iPSC models derived from ASD patients [18,19,20, 23, 27, 54, 66,67,68] and in clinical studies of ASD and other neuropsychiatric disorders [54,55,56,57,58]. These findings are intriguing since the AP and UM share the same CHRNA7 duplication, but only the AP exhibits these neurodevelopmental alterations, which may contribute to his pathogenesis. UM-dysregulated genes were linked to many similar GO terms, including axonal guidance, integrin and gap junctions, and behavior- and nervous system development; however, these were often upregulated in the UM versus both AP and UC models, suggesting compensatory molecular changes that may have normalized some corresponding neurodevelopmental processes in the UM, potentially contributing to an absence of clinical phenotypes. In contrast to these AP-specific neurodevelopmental deficits, electrophysiological analysis of more mature cells revealed functional anomalies that were common to both AP and UM cells, raising the possibility that maturation or the influence of astrocyte feeder cells may modify the phenotypic relationships among the UM-C, UM, and AP models.
Altered GABAergic neurotransmission is implicated in ASD affectation and may result from disrupted synaptic development, which can imbalance excitatory and inhibitory neuronal activity in the cortex [33, 69]. Most in vitro modeling studies characterize neurodevelopmental phenotypes only in excitatory, glutamatergic neurons [18, 27, 70, 71]; therefore, the extent to which GABAergic neurodevelopment is altered in ASD iPSC models remains largely unexplored. Here, we observed multiple phenotypes only in AP-derived cINs but not cExNs, suggesting that assessing both neuronal types may enable the identification of additional cell-type-specific potential contributors to disease. For example, we observed the reduced formation of VGAT-expressing punctae in AP-derived cINs, which could alter cIN function, potentially impairing synaptic activity to disrupt excitatory-inhibitory balance in the AP. We did not attempt to quantify a difference in functional connectivity owing to the low overall prevalence of spontaneous synaptic activity in all three of the models at the time point evaluated; however, this could be a promising topic for future study. In the cortex, inhibitory neurons constitute a minority population (~20%) [72, 73], are specified in sub-cortical progenitor territories, and undergo tangential migration to their cortical targets. Migration of these and other neurons is critical for normal neurodevelopment [74, 75]. Accordingly, here, the AP exhibited decreased cIN migration and diminished expression of corresponding suites of genes that regulate migration (Fig. 5a; Additional file 14). Other studies of ASD and other neuropsychiatric disorders have likewise revealed evidence of abnormalities in neurogenesis, neuronal migration, and neuronal maturation [2, 61, 76, 77]. Furthermore, Wnt signaling contributes to neuronal migration, and here many Wnt signaling genes had reduced expression in the AP model, while the Wnt agonist CHIR-99021 partially rescued this model’s migration deficit (Figs. 4e, 5; Additional file 14). Modulating Wnt signaling has been explored for ASD treatment; however, as both reduced and elevated Wnt activities are implicated in ASD-related behavioral and cognitive alterations, further studies are required to stratify how Wnt pathway biomarkers are altered in neuropsychiatric disorders including ASD [63,64,65].
This work assessed how CHRNA7 duplication affects neurodevelopment and function and identified potential contributors to variable phenotypic penetrance among 15q13.3 duplication carriers. As the study design involved comparing first-degree relatives with the same duplication but differential diagnoses, age, and sex represent potential confounds for data interpretation. ASD is more prevalent in males than females in the population, with evidence suggesting that females may require a higher load of genetic liability to exhibit clinical phenotypes . However, here the relationship of differential penetrance to sex differences is unclear, as 15q13.3 is autosomal rather than sex chromosome-linked, and current iPSC modeling approaches cannot recapitulate random X-chromosome inactivation that female somatic tissues, including the brain, undergo [79, 80]. However, in our RNA-seq analysis, very few (<3%) DEGs were either X-chromosome-linked or exhibit sex-biased expression in the human brain [81, 82]. Age is another unavoidable confound in studies involving modeling multiple subjects in the same pedigree. It is currently unclear whether iPSC lines derived from old versus young donors exhibit differences in their potential for differentiation or senescence, and this issue is controversial [83, 84]. To address these confounds, we performed variancePartition analysis, which indicated that differences between samples were predominantly driven by subject identity and genetic background, while age and sex were minimal contributors to DEG identification. Another confounding variable in many iPSC modeling studies is differences in genetic background. Here, the AP and UM have a partially shared genetic background, including the same 424kb duplication at 15q13.3. The large size of this duplication precludes CRISPR-mediated correction, which could enable the identification of duplication-linked phenotypes on an isogenic background. Finally, we could not access biological material from the father, who represents an unaffected male first-degree relative control. As 50% of the AP’s genome derives from the father, it would have been interesting to understand whether additional paternally inherited genetic liabilities could have contributed to the AP's phenotypic penetrance. However, biomaterials were not available for study from this subject. Despite these limitations, this work still identified cellular, molecular, and functional signatures that differed in clinically affected versus unaffected individuals carrying the same contributory genetic duplication on a partially shared genetic background, as well as defining functional alterations shared by both CHRNA7 duplication carriers, that distinguished them from unrelated controls lacking CHRNA7 duplication.
It is challenging to identify variables that contribute to the wide variability of clinical phenotypic penetrance among individuals carrying microduplications in the same genomic region. These anomalies usually cannot be recapitulated in animal models, confounding experimental study, while both the genomic interval, which genes are involved, and heterogeneous genetic background contributors in unrelated individuals can confound iPSC model cross-comparison. Here, using related individuals with the same duplication but variable phenotypes was informative in identifying neurodevelopmental deficits and elevated ER stress specific to AP-derived models and potentially contributing to his clinical diagnoses. Conversely, upregulation of overlapping developmental signaling pathways in the UM may have contributed to a lack of both neurodevelopmental phenotypes in this model and clinical diagnoses in this subject. Finally, despite their variable diagnoses, both models carrying CHRNA7 duplication exhibited shared functional abnormalities. Therefore, these findings highlight the potential for iPSC models to identify cellular and molecular anomalies linked to either the presence versus absence of a psychiatric disorder-linked CNV, such as 15q13.3 duplication, and to differences in phenotypic penetrance among multiple individuals carrying the same CNV.
It would be interesting to examine whether neurodevelopmental deficits, particularly in cINs, are associated with patient phenotype in models from other CHRNA7 duplication carriers, to determine whether the phenotypes identified here frequently predict clinical diagnoses in this patient population. It would also be intriguing to know whether these phenotypes are more broadly generalizable indicators of phenotypic penetrance in other CNVs involving genome microduplication, many of which exhibit similarly variable clinical phenotypic penetrance. Defining general hallmarks of and contributors to phenotypic penetrance in CNVs involving microduplication, at either 15q13.3 or other susceptible genomic intervals, would be informative in defining anomalies that contribute to these disorders, while these iPSC models provide a platform for molecular and pharmacological screening for phenotypic rescue, to develop interventions for treating affected individuals.
Clinical phenotypes of these three subjects with 15q13.3 duplication are summarized in Table 1. The AP’s pregnancy was planned and there was no history of in utero exposure to drugs, alcohol, or tobacco. There were no complications during the pregnancy, and the AP was delivered normally at full term, with a birth weight of seven pounds, six ounces. As an infant, the AP was delayed in sleeping through the night, which did not begin to occur until after 1 year of age. He has no history of delay in motor development but has a history of significant delay in language development, producing his first word after the age of 2 years. He started speech therapy at 2 years of age, resulting in rapid improvement in his language development. The AP was 12 years old when this study was initiated, with a history of ADHD, depression, and ASD. Prior to age five, he reportedly did not make eye contact and did not exhibit age-appropriate social reciprocity. While in the third grade, the AP had frequent crying episodes and was overwhelmed by homework, which brought him to psychological evaluation. At that time, he manifested both autistic features and pronounced mood lability on exam, which was manifested on several occasions by the child becoming morose and tearful when told even slightly sad stories that would not have elicited such a reaction in a typical child his age. During periods of depressive symptomatology, his mother reported that he had difficulty falling asleep, low appetite, and decreased interest in his favorite activities, such as sports. In middle school, he manifested poor concentration and difficulty with time management, standardized test-taking, and organizing tasks and activities. These issues were treated by using cognitive behavioral therapy and play therapy. He was subsequently treated with sertraline, followed by escitalopram. On these selective serotonin reuptake inhibitors, the AP’s anxiety was significantly lessened, but he then experienced a residual lack of motivation and his perseverative traits were not improved by treatment. Based on his developmental history, including language delay, impairment in social reciprocity and non-verbal communication, and repetitive thinking, it was clinically determined that a significant contribution to his overall impairment was autistic perseveration and rigidity, for which a trial of risperidone was initiated and resulted in significant clinical improvement over the ensuing years. Ultimately, he was successfully weaned from risperidone and was reasonably well adapted in high school. The AP’s mother reported a history of mild depression, anxiety, and obsessive-compulsive traits, while the AP’s 8-year-old brother had subtle autistic traits that were less pronounced than those of the AP and behavioral features of emotional dysregulation that were more pronounced than those of the AP; he met criteria for disruptive mood dysregulation disorder, ADHD, and mood disorder.
Cytogenetics Microarray (CMA) analysis was performed for research testing by the Washington University Cytogenetics and Molecular (CMP) Pathology Laboratory, using the Affymetrix CytoScan HD array. This array includes 2.6 million copy number markers, 1.9 million non-polymorphic probes, and nearly 750,000 single-nucleotide polymorphism (SNP) probes. Average intragenic marker spacing is equivalent to 1 probe per 880 basepairs. Analysis of these data by the CMP Laboratory, after alignment to hg19, defined a 424 kb gain at 15q13.3 in samples from the affected proband (AP) and unaffected mother (UM) and a 444 kb gain in the same location in the affected brother. This copy number variant was not present in the father.
The Washington University Genome Engineering and Induced Pluripotency Center (GEiC) derived multiple clonal iPSC lines from individuals in this family. Briefly, fresh urine samples were procured from the AP and UM and were used to obtain renal epithelial cells, which were reprogrammed using the CytoTune-iPSC 2.0 Sendai virus-based reprogramming kit (Thermo Fischer Scientific), as per the manufacturer’s instructions. iPSC clones were picked and three clonal cell lines were derived from each study subject. Clones number 1 and 3 from each subject were used for experimentation. These newly derived lines were compared with established male and female iPSC control lines from unrelated individuals (BJFF6 and AN1.1, denoted UC-M and UC-F in this study) provided by the Washington University Genome Engineering and Induced Pluripotency Center.
iPSC cultures and differentiation
iPSC line derivation and directed differentiation were performed by modification of our previously described methods [22, 85]. iPSC lines were grown on Matrigel (Corning) under feeder-free conditions using mTeSR1 (STEMCELL Technologies), authenticated by STR profiling, and tested for mycoplasma contamination regularly during culture. For directed differentiation to generate cortical excitatory neural progenitor cells (cExNPCs), iPSCs were dissociated into single cells with Accutase (Life Technologies) and 40,000 cells were seeded in V-bottom 96 well non-adherent plates (Corning). Embryoid bodies (EBs) were generated by centrifugation of the plate at 200xg for 5 min and were then incubated in 5% CO2 at 37°C, in cExNPC differentiation medium with 10μM Y-27632 (Tocris Biosciences). cExNPC differentiation medium includes Neurobasal-A (Life Technologies), 1X B-27 supplement without Vitamin A (Life Technologies), 10μM SB-431542 (Tocris Biosciences), and 100nM LDN-193189 (Tocris Biosciences). On day 4, EBs were picked with wide bore P1000 tips and were transferred to Poly-l-Ornithine- (20μg/ml) and laminin- (10μg/ml) coated plates. Every other day media without Y-27632 was replenished, and on day 15 Neural Rosette Selection Reagent (STEMCELL Technologies) was used to isolate cExNPCs from rosettes, as per the manufacturer’s instructions. cExNPCs were grown as a monolayer using cExNPC differentiation media up to 15 passages.
Cortical inhibitory neural progenitor cells (cINPCs) were generated by directed differentiation in media which included Neurobasal-A (Life Technologies), 1X B-27 supplement without Vitamin A (Life Technologies), 10μM SB-431542 (Tocris Biosciences), 100nM LDN-193189 (Tocris Biosciences), 1μM Purmorphamine (Calbiochem), and 2μM XAV-939 (Tocris Biosciences). Y-27632 was also included in this media until day eight. For cINPC differentiation, EBs were generated as described above for cExNPC differentiation. On day 4, EBs were transferred to non-adherent plates and were placed on an orbital shaker at 80rpm in an incubator with 5% CO2 and 37°C. cINPC media were replenished every other day, and on day 10, EBs were transferred to Matrigel- and laminin- (5μg/ml) coated plates. On day 15, cINPCs were dissociated with Accutase and were either cryopreserved and/or grown in monolayer culture on Matrigel- and laminin-coated for up to 15 passages. For analysis of both cExNPC and cINPC growth properties, equal numbers of cells for each line were seeded on Matrigel- and laminin- (5μg/ml) coated plates and the total number of cells was counted after 4 days.
For differentiation and maturation of neurons, cortical neuroids were generated by seeding both 2X104 cExNPCs and 2X104 cINPCs into each well of a V-bottom 96 well non-adherent plate in maturation media. Plates were spun at 200xg for 5 min and were incubated in 5% CO2 at 37°C in maturation media, with addition of Y-27632 for the first 4 days of culture. The composition of maturation media includes Neurobasal-A and 1X B-27 supplement without vitamin A, while DAPT (10μM; Tocris) was included in the media from day 7 to day 11, and 200μM ascorbic acid (Sigma Aldrich), 20ng/ml BDNF (Peprotech), and 200μM cAMP were included in the media from day 11 to day 15. On day 4, neuroids were transferred to non-adherent plates and were placed on an orbital shaker at 80rpm in an incubator with 5% CO2 and 37°C. On day 5, neuroids were moved to Matrigel- and laminin- (5μg/ml) coated plates and further incubated in 5% CO2 at 37°C. Media were replenished every other day until day 15.
Immunocytochemistry (ICC) and immunoblotting
ICC and immunoblotting experiments were performed as previously described . In brief, for ICC, cortical neuroids were dissociated after 15 days of maturation and cells were plated in eight-well chamber slides coated with Matrigel and laminin (5μg/ml). After 24 h, cells were washed with PBS without calcium and magnesium and were fixed in 4% paraformaldehyde for 15–20 min. See  for detailed protocol. Primary and secondary antibodies used for these experiments and for immunoblotting are provided in Additional file 4. Images were taken using a spinning-disk confocal microscope (Quorum) and an Olympus inverted microscope using MetaMorph software. ImageJ was used to process images and for quantification: 15–20 random fields were imaged from three to five biological replicate experiments, which included work with two different clones per sample type for the AP and UM models and one clone for the UC-M/F models; total numbers of both immune-positive and all DAPI stained cell nuclei quantified are shown in Additional file 4. For each experimental finding in this manuscript, the number of biological replicate experiments and the clones used for each replicate of each type of experiment is also summarized in Additional file 4.
Approximately one million cExNPCs or cINPCs between passages 4–9 were used for FACS analysis, and experiments were performed as described previously . The plot shows the median value gathered from seven biological replicates, using two different clones per sample type for the AP and UM models and one clone for the UC-M/F models, calculated by using a Kruskal-Wallis non-parametric test. P values: *P < 0.05, **P < 0.01, and ***P < 0.001.
RNA-Sequencing and RT-qPCR
After 15 days of cortical neuroid differentiation as described above, total RNA was collected from the AP, UM, UC-M, and UC-F lines, using the NucleoSpin RNA II kit (Takara) per the manufacturer’s instructions. RNA was quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific) and the Agilent Bioanalyzer 2100 was used to assess RNA integrity, with only samples with an RNA Integrity Number of >8 used for sequencing and analysis. RNA-Sequencing (RNA-Seq) library preparation and Illumina Sequencing were performed by the Genome Technology Access Center at Washington University. The Illumina Hi-Seq3000 was used to obtain single-end 50 base pair reads, with approximately 30 million unique reads per sample obtained after alignment. Four independent biological replicates per cell line were analyzed by RNA-Seq. For RT-qPCR, 1μg total RNA was reverse transcribed using iScript Reverse Transcription Supermix (Bio-Rad) and equal quantities of cDNA were used as a template for RT-qPCR using the Applied Biosystems Fast Real-Time quantitative PCR platform. GAPDH or RPL30 were used as endogenous controls for normalization. Four biological replicate experiments using one clonal line for each sample type were used for RNA-seq analysis (n = 4), while a second clonal line for the UM and AP models was used to generate RNA for RT-qPCR validation of a subset of the RNA-Seq findings. P values for RT-qPCR validation: *P < 0.05, **P < 0.01, and ***P < 0.001 were determined by unpaired t testing.
Bioinformatics and IPA analysis
RNA-Seq data analysis was performed as described in [22, 85] to obtain differentially regulated genes (DEG). Briefly, STAR version 2.5.4b was used to align the RNA-Seq reads to the human genome assembly hg38 . To derive uniquely aligned unambiguous reads, we used Subread:featureCount, version 1.6.3 with GENCODE gene annotation  and gene-level transcripts were imported into the R-bioconductor package . After excluding genes expressed at <1.0 counts per million (CPM), differentially expressed genes (DEG) were curated based upon a Log2 fold change >1 and a Benjamini and Hochberg FDR of <0.05. DEGs were used to perform hierarchical clustering analysis using ClustVis  and to perform Ingenuity Pathway Analysis (IPA) (Qiagen), as described previously . To determine the contribution of different covariates to gene expression, we also performed variancePartition analysis, including the individual sample types (UC-M, UC-F, UM, and AP), age (young and old), and sex (male and female) as variables .
To measure neurite extension, cortical neuroids were generated and cultured to promote differentiation and maturation as described above. On day 6, images were acquired using an inverted light microscope, with data collected for three or more independent biological replicate experiments encompassing work with two clonal lines per subject for the AP and UM, and one clonal line for the UC-F and UC-M. Neurite extension length from adherently plated neuroids was measured using ImageJ, as the distance between two circles drawn at the border of the plated neuroid and at the tips of neurites extending from that neuroid. Since neurites extend in all directions from each plated neuroid, each datapoint is a calculation of neurite length from the EB border to the tip of the neurites, calculated on a per-EB basis and normalized by the number of EBs quantified per sample type (e.g., UC-M). The plot shows the median value calculated from 50–80 data points per sample type, with data points gathered from seven biological replicates.
ICC for MAP2 in adherent neuroids was conducted as described above, with image acquisition using a spinning-disk (Quorum) confocal microscope and an Axiovision inverted microscope. Day 15 neuroids were also dissociated and the neurons plated on Matrigel- and laminin-coated plated and stained with MAP2.
To assess neurite length, images were acquired using a spinning-disk (Quorum) confocal microscope and an Axiovision inverted microscope. Images were processed with Imaris software (Bitplane) and neurite length was measured using the filament tracer application and normalized to the number of nuclei stained with DAPI in the neuroid, which was measured with the particle application in Imaris. The neuronal soma area was measured from MAP2 stained images using ImageJ. The data points are normalized on a per neuron basis, with neurons quantified from at least 20 different images (with at least five neurons measured per image) per biological replicate, with a total of three independent biological replicates from two clones analyzed.
To quantify VGAT- and VGLUT-expressing punctae, dissociated and plated neuroids generated as described above were immunostained for the respective antibodies, and images were taken using a spinning-disk (Quorum) confocal microscope and Axiovision inverted microscope. Punctae were measured using the synaptic counter plugin in ImageJ. Each finding was obtained in three or more independent biological replicate experiments encompassing work with two clonal lines per subject for the AP and UM, and one clonal line for the UC-M.
To study the migration of cExN and cIN neurons, we developed an approach that utilized fused co-culture of two 3D spheres consisting of cExNs and cINs. These 3D spheres were generated by transducing cExNPCs and cINPCs separately with either a lentiviral synapsin-eGFP or a synapsin-RFP expression construct, respectively. 30,000 of these transduced cExNPCs or cINPCs per well were then seeded into separate wells of a V bottom 96 well plate in 100μl of maturation media containing 10μM Y-27632. The V bottom plate was centrifuged at 200xg for 5 min at room temperature and then incubated in 5% CO2 at 37°C. On day 2, 50μl of media was replaced with fresh media without disturbing the spheres. On day 4, cExN and cIN EBs were selected with wide bore P1000 tips and moved to a U bottom plate. One cExN and one cIN sphere were placed side by side in each U bottom well. Placement of spheres in close apposition caused them to undergo fusion without further manipulation, enabling assessment of neuronal migration. On day 6, these fused spheres were moved with wide bore P1000 tips to a coverslip placed in a 3-cm plate and coated with matrigel and laminin (5μg/ml). On day 10, images were acquired using a spinning-disk (Quorum) confocal microscope and Axiovision inverted microscope and image processing was performed using ImageJ. Migration was assessed in such fused co-cultures. The ability of Wnt signaling to rescue migration deficits observed in the AP line was tested by the addition of 10μM CHIR-99021 in DMSO, with parallel treatment of control spheres with equal quantities of DMSO. Three or more independent biological replicate experiments were performed in two clonal lines for the AP and UM, and in one clonal line for the UC-F and UC-M. Each data point is the number of neurons that migrated from the cIN into cExN EB or from the cExN into the cIN EB, normalized to each fused EB (e.g., EB pair). Each data point represents the number of migrating cExNs per fused EB or the number of cINs per fused EB, with the median value calculated from ~40–80 fused EBs per sample type. The plot shows the median value calculated from 40 to 80 data points per sample type, with data points gathered from six biological replicates, using two different clones per sample type. P values were calculated by using a Kruskal-Wallis non-parametric test. *P < 0.05, **P < 0.01, and ***P < 0.001.
ER stress luciferase assay
To test the effects of endoplasmic reticulum (ER) stress on cINPCs, we used an expression construct encoding a stress sensor . This construct encodes a Guassia luciferase protein fusion, with replacement of the first 18 amino acids with the signal peptide from the mesencephalic astrocyte-derived neurotrophic factor (MANF) protein and carboxy-terminal fusion to MANF’s final 5 amino acids, which encode a stress sensor. To detect ER stress, 35,000 cINPCs were seeded on Matrigel- and laminin- (5μg/ml) coated 96 well plates in cIN differentiation media containing Y-27632. After 24 h, the cells were transfected with the stress sensor expression construct using FuGENE 6 (Promega) transfection reagent. 48 h after transfection, 50μl of supernatant was removed and assayed for luciferase activity using the BioLux Gaussia Luciferase Assay Detection System (New England Biolabs). For rescue experiments, small molecules were obtained from Sigma Aldrich or Tocris Biosciences, reconstituted in DMSO or PBS-Ca2+/Mg2+, and added in the medium after 24 h of transfection, at the final concentrations indicated (Tudca-50μM, PBA-500μM, Dantrolene sodium-1μM, and JTV-519-10μM). Luciferase levels were measured after 48 h of small molecule treatment as above. Luciferase data for rescue experiments performed in the presence of small molecules were normalized to the DMSO control, and three or more independent biological replicate experiments were performed, using two clonal lines for the AP and UM and one clonal line for the UC-F and UC-M.
cExNPCs and cINPCs were transduced with Synapsin-GFP and Synapsin-RFP expression constructs, respectively, before performing cortical neuroid maturation. At day 0 of cortical neuroid differentiation, 2X104 cExNPCs and 2X104 cINPCs were mixed in each well of a V-bottom 96-well non-adherent plate in maturation media. Maturation followed the approach described above. At day 15, neuroids were dissociated using Accutase (Life Technologies), were seeded onto a layer of rat cortical astrocytes, prepared as described previously , and were grown for another three weeks using Neurobasal-A, 1XB27 with vitamin A (Life Technologies), and supplementation with BDNF (20ng), cAMP (200μM), and ascorbic acid (200μM). iPSC-derived neuron cultures were perfused at 1 ml/min with room temperature (22°C) Tyrode’s solution (in mM): 150 NaCl, 4 KCl, 2 MgCl2, 2 CaCl2, 10 glucose, 10 HEPES, with pH adjusted to pH 7.4 with NaOH. Recording electrodes had an open-tip resistance of 2–6 MOhm when filled with (in mM): 140 K-glucuronate, 10 NaCl, 5 MgCl2, 0.2 EGTA, 5 Na-ATP, 1 Na-GTP, and 10 HEPES, pH adjusted to 7.4 with KOH. Whole-cell currents and membrane potentials were recorded with an Axopatch 200A amplifier (Molecular Devices). Voltage clamp recordings were used to determine cell capacitance and input resistance as well as peak inward sodium current and steady-state outward potassium currents during depolarizing voltage steps from a holding potential of −80 mV . Current clamp recordings and currents evoked by choline and acetylcholine (ACh) were obtained in a modified extracellular solution (in mM): 120 NaCl, 3 KCl, 10 glucose, 1 NaH2PO4, 4 NaHCO3, 5 HEPES, pH adjusted to 7.4 with NaOH, delivered from an 8-barrelled local perfusion pipette positioned very close to the recorded cell in order to minimize desensitization. Although desensitization can be further reduced by recording from outside-out patches that enable faster solution exchange, the present study employed whole-cell recordings in order to allow for current and voltage-clamp analysis of cellular physiology as well as agonist-gated currents all within the same cells. In future experiments, higher speed agonist delivery would allow for a more detailed kinetic comparison among the genotypes.
All statistical analyses were performed using IBM SPSS Statistics (v.27) or Sigma-STAT. Prior to analyses, data was screened for missing values and fit of distributions with assumptions underlying univariate analysis. This included the Shapiro-Wilk test on z-score-transformed data and qqplot investigations for normality, Levene’s test for homogeneity of variance, and boxplot and z-score (±3.29) investigation for identification of influential outliers. Means and standard errors were computed for each measure. Non-parametric Kruskal-Wallis and Mann-Whitney U tests were used to analyze data. The critical alpha value for all analyses was P <0.05, unless otherwise stated. Multiple pairwise comparisons were subjected to Bonferroni correction, where appropriate. The datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.
Availability of data and materials
All data generated and analyzed during this study are included in this published article, its additional information files, or in the publically available Gene Expression Omnibus (GEO) repository (Series GSE143908) . Data values are provided in Additional files 4, 5, or 13, as indicated. iPSC lines used in this study are available from the investigator by material transfer agreement with Washington University.
Autism spectrum disorder
Induced pluripotent stem cells
Copy number variant
α-7 nicotinic acetylcholine receptor subunit
Attention deficit and hyperactivity disorder
Unfolded protein response
Cortical excitatory neuron
Cortical inhibitory neuron
Neural progenitor cell
Genome Technology Access Center
Reverse transcription and quantitative PCR
Differentially expressed gene
Ingenuity Pathway Analysis
Principal Component Analysis
Grayton HM, Fernandes C, Rujescu D, Collier DA. Copy number variations in neurodevelopmental disorders. Prog Neurobiol. 2012;99(1):81–91. https://doi.org/10.1016/j.pneurobio.2012.07.005.
Cook EH Jr, Scherer SW. Copy-number variations associated with neuropsychiatric conditions. Nature. 2008;455(7215):919–23. https://doi.org/10.1038/nature07458.
Deshpande A, Weiss LA. Recurrent reciprocal copy number variants: roles and rules in neurodevelopmental disorders. Dev Neurobiol. 2018;78(5):519–30. https://doi.org/10.1002/dneu.22587.
Rosenfeld JA, Coe BP, Eichler EE, Cuckle H, Shaffer LG. Estimates of penetrance for recurrent pathogenic copy-number variations. Genet Med. 2013;15(6):478–81. https://doi.org/10.1038/gim.2012.164.
Newman S, Hermetz KE, Weckselblatt B, Rudd MK. Next-generation sequencing of duplication CNVs reveals that most are tandem and some create fusion genes at breakpoints. Am J Hum Genet. 2015;96(2):208–20. https://doi.org/10.1016/j.ajhg.2014.12.017.
Fernandez BA, Roberts W, Chung B, Weksberg R, Meyn S, Szatmari P, et al. Phenotypic spectrum associated with de novo and inherited deletions and duplications at 16p11.2 in individuals ascertained for diagnosis of autism spectrum disorder. J Med Genet. 2010;47(3):195–203. https://doi.org/10.1136/jmg.2009.069369.
Gillentine MA, Schaaf CP. The human clinical phenotypes of altered CHRNA7 copy number. Biochem Pharmacol. 2015;97(4):352–62. https://doi.org/10.1016/j.bcp.2015.06.012.
Szafranski P, Schaaf CP, Person RE, Gibson IB, Xia Z, Mahadevan S, et al. Structures and molecular mechanisms for common 15q13.3 microduplications involving CHRNA7: benign or pathological? Hum Mutat. 2010;31(7):840–50. https://doi.org/10.1002/humu.21284.
Lowther C, Costain G, Stavropoulos DJ, Melvin R, Silversides CK, Andrade DM, et al. Delineating the 15q13.3 microdeletion phenotype: a case series and comprehensive review of the literature. Genet Med. 2015;17(2):149–57. https://doi.org/10.1038/gim.2014.83.
Ziats MN, Goin-Kochel RP, Berry LN, Ali M, Ge J, Guffey D, et al. The complex behavioral phenotype of 15q13.3 microdeletion syndrome. Genet Med. 2016;18(11):1111–8. https://doi.org/10.1038/gim.2016.9.
Williams NM, Franke B, Mick E, Anney RJ, Freitag CM, Gill M, et al. Genome-wide analysis of copy number variants in attention deficit hyperactivity disorder: the role of rare variants and duplications at 15q13.3. Am J Psychiatry. 2012;169(2):195–204. https://doi.org/10.1176/appi.ajp.2011.11060822.
Makoff AJ, Flomen RH. Detailed analysis of 15q11-q14 sequence corrects errors and gaps in the public access sequence to fully reveal large segmental duplications at breakpoints for Prader-Willi, Angelman, and inv dup(15) syndromes. Genome Biol. 2007;8(6):R114. https://doi.org/10.1186/gb-2007-8-6-r114.
International Schizophrenia C. Rare chromosomal deletions and duplications increase risk of schizophrenia. Nature. 2008;455(7210):237–41. https://doi.org/10.1038/nature07239.
van Bon BW, Mefford HC, Menten B, Koolen DA, Sharp AJ, Nillesen WM, et al. Further delineation of the 15q13 microdeletion and duplication syndromes: a clinical spectrum varying from non-pathogenic to a severe outcome. J Med Genet. 2009;46(8):511–23. https://doi.org/10.1136/jmg.2008.063412.
Al Ageeli E, Drunat S, Delanoe C, Perrin L, Baumann C, Capri Y, et al. Duplication of the 15q11-q13 region: clinical and genetic study of 30 new cases. Eur J Med Genet. 2014;57(1):5–14. https://doi.org/10.1016/j.ejmg.2013.10.008.
Sinkus ML, Graw S, Freedman R, Ross RG, Lester HA, Leonard S. The human CHRNA7 and CHRFAM7A genes: a review of the genetics, regulation, and function. Neuropharmacology. 2015;96(Pt B):274–88.
Adam MP, Ardinger HH, Pagon RA, et al., editors. GeneReviews® [Internet]. Seattle (WA): University of Washington, Seattle; 1993-2021. Available from: https://www.ncbi.nlm.nih.gov/books/NBK1116/.
Mekahli D, Bultynck G, Parys JB, De Smedt H, Missiaen L. Endoplasmic-reticulum calcium depletion and disease. Cold Spring Harb Perspect Biol. 2011;3(6):a004317. Published 2011 Jun 1. https://doi.org/10.1101/cshperspect.a004317.
Deshpande A, Yadav S, Dao DQ, Wu ZY, Hokanson KC, Cahill MK, et al. Cellular phenotypes in human iPSC-derived neurons from a genetic model of autism spectrum disorder. Cell Rep. 2017;21(10):2678–87. https://doi.org/10.1016/j.celrep.2017.11.037.
Kathuria A, Nowosiad P, Jagasia R, Aigner S, Taylor RD, Andreae LC, et al. Stem cell-derived neurons from autistic individuals with SHANK3 mutation show morphogenetic abnormalities during early development. Mol Psychiatry. 2018;23(3):735–46. https://doi.org/10.1038/mp.2017.185.
Wang P, Lin M, Pedrosa E, Hrabovsky A, Zhang Z, Guo W, et al. CRISPR/Cas9-mediated heterozygous knockout of the autism gene CHD8 and characterization of its transcriptional networks in neurodevelopment. Mol Autism. 2015;6(1):55. https://doi.org/10.1186/s13229-015-0048-6.
Lewis EMA, Meganathan K, Baldridge D, Gontarz P, Zhang B, Bonni A, et al. Cellular and molecular characterization of multiplex autism in human induced pluripotent stem cell-derived neurons. Mol Autism. 2019;10(1):51. https://doi.org/10.1186/s13229-019-0306-0.
Griesi-Oliveira K, Acab A, Gupta AR, Sunaga DY, Chailangkarn T, Nicol X, et al. Modeling non-syndromic autism and the impact of TRPC6 disruption in human neurons. Mol Psychiatry. 2015;20(11):1350–65. https://doi.org/10.1038/mp.2014.141.
Woodbury-Smith M, Deneault E, Yuen RKC, Walker S, Zarrei M, Pellecchia G, et al. Mutations in RAB39B in individuals with intellectual disability, autism spectrum disorder, and macrocephaly. Mol Autism. 2017;8(1):59. https://doi.org/10.1186/s13229-017-0175-3.
DeRosa BA, El Hokayem J, Artimovich E, Garcia-Serje C, Phillips AW, Van Booven D, et al. Convergent pathways in idiopathic autism revealed by time course transcriptomic analysis of patient-derived neurons. Sci Rep. 2018;8(1):8423. https://doi.org/10.1038/s41598-018-26495-1.
Liu X, Campanac E, Cheung HH, Ziats MN, Canterel-Thouennon L, Raygada M, et al. Idiopathic autism: cellular and molecular phenotypes in pluripotent stem cell-derived neurons. Mol Neurobiol. 2017;54(6):4507–23. https://doi.org/10.1007/s12035-016-9961-8.
Marchetto MC, Belinson H, Tian Y, Freitas BC, Fu C, Vadodaria K, et al. Altered proliferation and networks in neural cells derived from idiopathic autistic individuals. Mol Psychiatry. 2017;22(6):820–35. https://doi.org/10.1038/mp.2016.95.
Germain ND, Chen PF, Plocik AM, Glatt-Deeley H, Brown J, Fink JJ, et al. Gene expression analysis of human induced pluripotent stem cell-derived neurons carrying copy number variants of chromosome 15q11-q13.1. Mol Autism. 2014;5:44.
Yin J, Chen W, Yang H, Xue M, Schaaf CP. Chrna7 deficient mice manifest no consistent neuropsychiatric and behavioral phenotypes. Sci Rep. 2017;7(1):39941. https://doi.org/10.1038/srep39941.
Wang Y, Xiao C, Indersmitten T, Freedman R, Leonard S, Lester HA. The duplicated alpha7 subunits assemble and form functional nicotinic receptors with the full-length alpha7. J Biol Chem. 2014;289(38):26451–63. https://doi.org/10.1074/jbc.M114.582858.
Gillentine MA, Yin J, Bajic A, Zhang P, Cummock S, Kim JJ, et al. Functional consequences of CHRNA7 copy-number alterations in induced pluripotent stem cells and neural progenitor cells. Am J Hum Genet. 2017;101(6):874–87. https://doi.org/10.1016/j.ajhg.2017.09.024.
Donovan AP, Basson MA. The neuroanatomy of autism - a developmental perspective. J Anat. 2017;230(1):4–15. https://doi.org/10.1111/joa.12542.
Zikopoulos B, Barbas H. Altered neural connectivity in excitatory and inhibitory cortical circuits in autism. Front Hum Neurosci. 2013;7:609.
Coe BP, Stessman HAF, Sulovari A, Geisheker MR, Bakken TE, Lake AM, et al. Neurodevelopmental disease genes implicated by de novo mutation and copy number variation morbidity. Nat Genet. 2019;51(1):106–16. https://doi.org/10.1038/s41588-018-0288-4.
Kaplanis J, Samocha KE, Wiel L, et al. Evidence for 28 genetic disorders discovered by combining healthcare and research data. Nature. 2020;586(7831):757–62. https://doi.org/10.1038/s41586-020-2832-5.
Pain O, Pocklington AJ, Holmans PA, Bray NJ, O'Brien HE, Hall LS, et al. Novel insight into the etiology of autism spectrum disorder gained by integrating expression data with genome-wide association statistics. Biol Psychiatry. 2019;86(4):265–73. https://doi.org/10.1016/j.biopsych.2019.04.034.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(Web Server issue):W214–20.
Snel B, Lehmann G, Bork P, Huynen MA. STRING: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene. Nucleic Acids Res. 2000;28(18):3442–4. https://doi.org/10.1093/nar/28.18.3442.
Maino B, Ciotti MT, Calissano P, Cavallaro S. Transcriptional analysis of apoptotic cerebellar granule neurons following rescue by gastric inhibitory polypeptide. Int J Mol Sci. 2014;15(4):5596–622. https://doi.org/10.3390/ijms15045596.
Miwa T, Manabe Y, Kurokawa K, Kamada S, Kanda N, Bruns G, et al. Structure, chromosome location, and expression of the human smooth muscle (enteric type) gamma-actin gene: evolution of six human actin genes. Mol Cell Biol. 1991;11(6):3296–306. https://doi.org/10.1128/mcb.11.6.3296-3306.1991.
Torii T, Miyamoto Y, Nakamura K, Maeda M, Yamauchi J, Tanoue A. Arf6 guanine-nucleotide exchange factor, cytohesin-2, interacts with actinin-1 to regulate neurite extension. Cell Signal. 2012;24(9):1872–82. https://doi.org/10.1016/j.cellsig.2012.05.016.
Munoz-Cobo JP, Sanchez-Hernandez N, Gutierrez S, El Yousfi Y, Montes M, Gallego C, et al. Transcriptional elongation regulator 1 affects transcription and splicing of genes associated with cellular morphology and cytoskeleton dynamics and is required for neurite outgrowth in neuroblastoma cells and primary neuronal cultures. Mol Neurobiol. 2017;54(10):7808–23. https://doi.org/10.1007/s12035-016-0284-6.
Kalinowska M, Chavez AE, Lutzu S, Castillo PE, Bukauskas FF, Francesconi A. Actinin-4 governs dendritic spine dynamics and promotes their remodeling by metabotropic glutamate receptors. J Biol Chem. 2015;290(26):15909–20. https://doi.org/10.1074/jbc.M115.640136.
Lopes-Ramos CM, Chen CY, Kuijjer ML, Paulson JN, Sonawane AR, Fagny M, et al. Sex differences in gene expression and regulatory networks across 29 human tissues. Cell Rep. 2020;31(12):107795. https://doi.org/10.1016/j.celrep.2020.107795.
Bocchi R, Egervari K, Carol-Perdiguer L, Viale B, Quairiaux C, De Roo M, et al. Perturbed Wnt signaling leads to neuronal migration delay, altered interhemispheric connections and impaired social behavior. Nat Commun. 2017;8(1):1158. https://doi.org/10.1038/s41467-017-01046-w.
Henderson MJ, Wires ES, Trychta KA, Richie CT, Harvey BK. SERCaMP: a carboxy-terminal protein modification that enables monitoring of ER calcium homeostasis. Mol Biol Cell. 2014;25(18):2828–39. https://doi.org/10.1091/mbc.e14-06-1141.
Alkondon M, Pereira EF, Barbosa CT, Albuquerque EX. Neuronal nicotinic acetylcholine receptor activation modulates gamma-aminobutyric acid release from CA1 neurons of rat hippocampal slices. J Pharmacol Exp Ther. 1997;283(3):1396–411.
Khiroug SS, Harkness PC, Lamb PW, Sudweeks SN, Khiroug L, Millar NS, et al. Rat nicotinic ACh receptor alpha7 and beta2 subunits co-assemble to form functional heteromeric nicotinic receptor channels. J Physiol. 2002;540(Pt 2):425–34. https://doi.org/10.1113/jphysiol.2001.013847.
Zwart R, Strotton M, Ching J, Astles PC, Sher E. Unique pharmacology of heteromeric alpha7beta2 nicotinic acetylcholine receptors expressed in Xenopus laevis oocytes. Eur J Pharmacol. 2014;726:77–86. https://doi.org/10.1016/j.ejphar.2014.01.031.
St John PA. Cellular trafficking of nicotinic acetylcholine receptors. Acta Pharmacol Sin. 2009;30(6):656–62. https://doi.org/10.1038/aps.2009.76.
Miller DT, Shen Y, Weiss LA, Korn J, Anselm I, Bridgemohan C, et al. Microdeletion/duplication at 15q13.2q13.3 among individuals with features of autism and other neuropsychiatric disorders. J Med Genet. 2009;46(4):242–8. https://doi.org/10.1136/jmg.2008.059907.
Krebs J, Agellon LB, Michalak M. Ca(2+) homeostasis and endoplasmic reticulum (ER) stress: an integrated view of calcium signaling. Biochem Biophys Res Commun. 2015;460(1):114–21. https://doi.org/10.1016/j.bbrc.2015.02.004.
Hakamata Y, Nakai J, Takeshima H, Imoto K. Primary structure and distribution of a novel ryanodine receptor/calcium release channel from rabbit brain. FEBS Lett. 1992;312(2-3):229–35. https://doi.org/10.1016/0014-5793(92)80941-9.
Jellinger K, Armstrong D, Zoghbi HY, Percy AK. Neuropathology of Rett syndrome. Acta Neuropathol. 1988;76(2):142–58. https://doi.org/10.1007/BF00688098.
McFadden K, Minshew NJ. Evidence for dysregulation of axonal growth and guidance in the etiology of ASD. Front Hum Neurosci. 2013;7:671.
Wang Z, Li P, Wu T, Zhu S, Deng L, Cui G. Axon guidance pathway genes are associated with schizophrenia risk. Exp Ther Med. 2018;16(6):4519–26. https://doi.org/10.3892/etm.2018.6781.
Lilja J, Ivaska J. Integrin activity in neuronal connectivity. J Cell Sci. 2018;131(12).
Carter MD, Shah CR, Muller CL, Crawley JN, Carneiro AM, Veenstra-VanderWeele J. Absence of preference for social novelty and increased grooming in integrin beta3 knockout mice: initial studies and future directions. Autism Res. 2011;4(1):57–67. https://doi.org/10.1002/aur.180.
Thanseem I, Nakamura K, Anitha A, Suda S, Yamada K, Iwayama Y, et al. Association of transcription factor gene LMX1B with autism. PLoS One. 2011;6(8):e23738. https://doi.org/10.1371/journal.pone.0023738.
Guglielmi L, Servettini I, Caramia M, Catacuzzeno L, Franciolini F, D'Adamo MC, et al. Update on the implication of potassium channels in autism: K(+) channelautism spectrum disorder. Front Cell Neurosci. 2015;9:34.
Reiner O, Karzbrun E, Kshirsagar A, Kaibuchi K. Regulation of neuronal migration, an emerging topic in autism spectrum disorders. J Neurochem. 2016;136(3):440–56. https://doi.org/10.1111/jnc.13403.
Blueprint Genetics [Internet]. Neuronal Migration Disorder Panel. 2021. https://blueprintgenetics.com/tests/panels/malformations/neuronal-migration-disorder-panel/.
Kwan V, Unda BK, Singh KK. Wnt signaling networks in autism spectrum disorder and intellectual disability. J Neurodev Disord. 2016;8(1):45. https://doi.org/10.1186/s11689-016-9176-3.
Kalkman HO. A review of the evidence for the canonical Wnt pathway in autism spectrum disorders. Mol Autism. 2012;3(1):10. https://doi.org/10.1186/2040-2392-3-10.
Bae SM, Hong JY. The Wnt signaling pathway and related therapeutic drugs in autism spectrum disorder. Clin Psychopharmacol Neurosci. 2018;16(2):129–35. https://doi.org/10.9758/cpn.2018.16.2.129.
Belichenko PV, Oldfors A, Hagberg B, Dahlstrom A. Rett syndrome: 3-D confocal microscopy of cortical pyramidal dendrites and afferents. Neuroreport. 1994;5(12):1509–13. https://doi.org/10.1097/00001756-199407000-00025.
Marchetto MC, Brennand KJ, Boyer LF, Gage FH. Induced pluripotent stem cells (iPSCs) and neurological disease modeling: progress and promises. Hum Mol Genet. 2011;20(R2):R109–15. https://doi.org/10.1093/hmg/ddr336.
Forrest MP, Parnell E, Penzes P. Dendritic structural plasticity and neuropsychiatric disease. Nat Rev Neurosci. 2018;19(4):215–34. https://doi.org/10.1038/nrn.2018.16.
Oblak AL, Gibbs TT, Blatt GJ. Decreased GABA(B) receptors in the cingulate cortex and fusiform gyrus in autism. J Neurochem. 2010;114(5):1414–23. https://doi.org/10.1111/j.1471-4159.2010.06858.x.
Ross PJ, Zhang WB, Mok RSF, Zaslavsky K, Deneault E, D'Abate L, et al. Synaptic dysfunction in human neurons with autism-associated deletions in PTCHD1-AS. Biol Psychiatry. 2020;87(2):139–49. https://doi.org/10.1016/j.biopsych.2019.07.014.
Thomas CA, Tejwani L, Trujillo CA, Negraes PD, Herai RH, Mesci P, et al. Modeling of TREX1-dependent autoimmune disease using human stem cells highlights L1 accumulation as a source of neuroinflammation. Cell Stem Cell. 2017;21(3):319–31 e8. https://doi.org/10.1016/j.stem.2017.07.009.
Kessaris N, Magno L, Rubin AN, Oliveira MG. Genetic programs controlling cortical interneuron fate. Curr Opin Neurobiol. 2014;26:79–87. https://doi.org/10.1016/j.conb.2013.12.012.
Han W, Sestan N. Cortical projection neurons: sprung from the same root. Neuron. 2013;80(5):1103–5. https://doi.org/10.1016/j.neuron.2013.11.016.
Lambert de Rouvroit C, Goffinet AM. Neuronal migration. Mech Dev. 2001;105(1-2):47–56. https://doi.org/10.1016/S0925-4773(01)00396-3.
Kriegstein AR, Noctor SC. Patterns of neuronal migration in the embryonic cortex. Trends Neurosci. 2004;27(7):392–9. https://doi.org/10.1016/j.tins.2004.05.001.
Avino TA, Hutsler JJ. Abnormal cell patterning at the cortical gray-white matter boundary in autism spectrum disorders. Brain Res. 2010;1360:138–46. https://doi.org/10.1016/j.brainres.2010.08.091.
Piven J, Berthier ML, Starkstein SE, Nehme E, Pearlson G, Folstein S. Magnetic resonance imaging evidence for a defect of cerebral cortical development in autism. Am J Psychiatry. 1990;147(6):734–9. https://doi.org/10.1176/ajp.147.6.734.
Robinson EB, Lichtenstein P, Anckarsater H, Happe F, Ronald A. Examining and interpreting the female protective effect against autistic behavior. Proc Natl Acad Sci U S A. 2013;110(13):5258–62. https://doi.org/10.1073/pnas.1211070110.
Sahakyan A, Kim R, Chronis C, Sabri S, Bonora G, Theunissen TW, et al. Human naive pluripotent stem cells model X chromosome dampening and X inactivation. Cell Stem Cell. 2017;20(1):87–101. https://doi.org/10.1016/j.stem.2016.10.006.
Dandulakis MG, Meganathan K, Kroll KL, Bonni A, Constantino JN. Complexities of X chromosome inactivation status in female human induced pluripotent stem cells-a brief review and scientific update for autism research. J Neurodev Disord. 2016;8(1):22. https://doi.org/10.1186/s11689-016-9155-8.
Huret JL, Ahmad M, Arsaban M, Bernheim A, Cigna J, Desangles F, et al. Atlas of genetics and cytogenetics in oncology and haematology in 2013. Nucleic Acids Res. 2013;41(Database issue):D920–4. https://doi.org/10.1093/nar/gks1082.
Tukiainen T, Villani AC, Yen A, Rivas MA, Marshall JL, Satija R, et al. Landscape of X chromosome inactivation across human tissues. Nature. 2017;550(7675):244–8. https://doi.org/10.1038/nature24265.
Strassler ET, Aalto-Setala K, Kiamehr M, Landmesser U, Krankel N. Age is relative-impact of donor age on induced pluripotent stem cell-derived cell functionality. Front Cardiovasc Med. 2018;5:4. https://doi.org/10.3389/fcvm.2018.00004.
Mahmoudi S, Brunet A. Aging and reprogramming: a two-way street. Curr Opin Cell Biol. 2012;24(6):744–56. https://doi.org/10.1016/j.ceb.2012.10.004.
Meganathan K, Lewis EMA, Gontarz P, Liu S, Stanley EG, Elefanty AG, et al. Regulatory networks specifying cortical interneurons from human embryonic stem cells reveal roles for CHD2 in interneuron development. Proc Natl Acad Sci U S A. 2017;114(52):E11180–E9. https://doi.org/10.1073/pnas.1712365115.
Metsalu T, Vilo J. ClustVis: a web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 2015;43(W1):W566–70. https://doi.org/10.1093/nar/gkv468.
Hoffman GE, Schadt EE. variancePartition: interpreting drivers of variation in complex gene expression studies. BMC Bioinformatics. 2016;17(1):483.
Meganathan K, Prakasam R, Baldridge D, Gontarz P, Zhang B, Urano F, et al. Alterations in neuronal physiology, development, and function associated with a common duplication of chromosome 15 involving CHRNA7. Gene Expression Omnibus. 2021. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=gse143908.
We thank the family for providing biomaterials for use in this study. We thank the Genome Engineering and iPSC Center (GEiC) at Washington University School of Medicine (WUSM) for deriving the iPSC lines used here. We thank the Alvin J. Siteman Cancer Center at WUSM, which is supported in part by the NCI Cancer Center Support Grant #P30 CA091842, for the use of their Flow Cytometry Core. We thank the Genome Technology Access Center in the Department of Genetics at WUSM (supported by grants #P30 CA91842 from the Siteman Cancer Center and ICTS/CTSA Grant # UL1 TR000448 from NIH/NCRR) for genomic sequencing services. We also thank WUSM Cytogenetics & Molecular Pathology for karyotyping services. In addition, we thank the reviewers for their thoughtful comments, which have improved this manuscript. A preprint of this manuscript was deposited into bioRxiv https://www.biorxiv.org/content/10.1101/2020.01.28.922187v2.
This work was supported by grants U01 HG007530 (NINDS/NIH Common Fund), R01NS114551, R01GM66815, and the M-CM Network, Jakob Gene Fund, and Washington University (WU) Center for Regenerative Medicine, WU Children’s Discovery Institute, and WU Institute for Clinical and Translational Sciences to KLK and by U54 HD087011-05S1 to JNC and KLK (Project PI, administrative supplement). This work was also partly supported by grants from the National Institutes of Health/NIDDK (DK112921, DK020579) to F. Urano.
Ethics approval and consent to participate
Subjects were consented for iPSC line generation by the Washington University Institutional Review Board of the Human Research Protection Office under human studies protocol #201409091 (Dr. John Constantino).
Consent for publication
Consent to publish data was provided by all subjects.
F. Urano received JTV-519 from the National Center for Advancing Translational Sciences for developing small molecule-based therapies for ER stress-related disorders and shares the intellectual property rights related to JTV-519 with the National Institutes of Health. F. Urano is an inventor of the patent related to ER calcium stabilizers, 10,441,574, B2 TREATMENT FOR WOLFRAM SYNDROME AND OTHER ER STRESS DISORDERS. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Characteristics of subjects in pedigree modeled here. (.jpg). (A) The study samples were derived from a pedigree with 15q13.3 duplication, with differential clinical affectation indicated by shading of subjects. The affected proband (AP) is represented in black, his affected brother, shown in gray, exhibits subtle autistic traits and has more volatile emotional dysregulation than the AP, while the unaffected mother (UM) and father are shown in white. Renal epithelial cells from the family members indicated (*) were used to derive iPSC models. (B) CNV array data for the AP and UM shows the signal intensity (log2 weighted ratio) and predicted copy number across the duplicated region of 15q13.3. The region lacking signal was not covered by the CNV array used for this analysis.
Characterization of iPSC models. (.jpg). Renal epithelial cell-derived iPSC lines from the UC-M, UM, and AP subjects (A) exhibit normal human stem cell colony morphology in bright field images (scale bar = 250 μm), (B) express the pluripotency marker OCT4/POU5F1 (scale bar = 150 μm), and (C) have a normal karyotype.
Cell cycle analysis of neural progenitor cells (NPCs). (.jpg). (A) cExNPCs and (B) cINPCs were stained with propidium iodide for DNA content and analyzed by FACS. Percentages of cells in each phase of the cell cycle were quantified for each model. Values shown are from seven independent biological replicate experiments (n = 7), using two clonal lines for the UM and AP, and one clonal line for the UC-M and UC-F. p-values were calculated by using a Kruskal-Wallis non-parametric test, as described in the Methods: *P < 0.05, **P < 0.01, ***P < 0.001.
Tables describing clones used and replicates performed, antibodies, cell counts, and data values. (.xlsx). (A) Indication of which clonal line was used for each biological replicate experiment and the number of biological replicates that were performed to procure the data for each figure panel. (B) Antibodies used for immunocytochemistry and immunoblotting, with the dilutions, suppliers, and host species indicated. (C) Number of Ki67-expressing NPCs and total number of DAPI-stained nuclei quantified to define the Ki67-expressing fraction in Fig. 2F. For Fig. 3D, total number of neurons used for quantitation of neurite length, based upon number of DAPI-expressing nuclei. Additional sheets provide data values for each figure panel involving <6 replicate experiments are also provided here or as indicated.
Differentially expressed genes. (.xlsx). DEGs were obtained by pairwise comparisons of normalized RNA-seq expression data from the four models, including calculation of the log2 fold change, FDR corrected p-values (padj), and average RPKM values across the sample types analyzed. See Methods for further information. (n = 4).
Ingenuity Pathways Analysis (IPA) of DEGs specific to the AP by comparison with the UM. (.xlsx). Significantly enriched: (A) Pathway and (B) Disease-related GO terms are shown, with the p-value and the DEGs from which each term was derived. (C) The top 1500 DEGs from the AP vs. UM comparison (based upon fdr/adj p-value) were input into the Percayai suite's CompBio pathway analysis tool to obtain networks associated with these DEGs. These networks were filtered with a concept filter utilizing the search terms "Wnt AND neuronal migration" to identify a network of interconnected themes (shown); concepts and underlying genes (entities) within each theme in the network above are shown. (D) Genes (entities) present in the Wnt AND Neuron Migration/Cadherin AND neural/interneuron themes are shown. These include two genes (GJA1/TGFBR2) that overlapped between Wnt-related and neuronal migration-related DEGs also called by Ingenuity Pathway Analysis. n = 4; data values are in Additional file 5.
Differentially expressed genes in the AP, by comparison with all three other models. (.jpg). (A) Venn diagram shows numbers of differentially expressed genes (DEGs) obtained from pairwise comparisons of the AP versus (vs.) UM, AP vs. UC-M, and AP vs. UC-F models. AP-specific DEGs, based upon comparisons to at least two of the other datasets, are shaded in blue. These AP-specific DEGs were further analyzed by: (B) Hierarchical clustering analysis, visualizing comparisons with the other three sample types, and by (C-H) Ingenuity Pathway Analysis (IPA), which identified (C) enriched pathways and (D) disease-related GO terms. In C-D, the number of DEGs enriched for each term present is represented on the x-axis, with red and blue colors indicating up- and down-regulated genes, respectively. p-values for each enriched GO term are indicated. (E-F) IPA disease terms enriched in these AP-specific DEGs include gene networks associated with (E) Cellular movement and (F) Nervous system development and function. The numbers of up-and down-regulated genes present in the networks are indicated. Within each network, red and green symbols indicate up- and down-regulated genes respectively, while color intensity indicates the relative degree of differential expression. (G) Interaction network of differentially expressed genes with known significance in neurodevelopmental disorders (NDDs). In this network, the nodes are each of the genes, the size of the node corresponds to the size of the differential expression, and the color indicates the direction of the differential expression change (red = negative, blue = positive). n = 4; data values are in Additional file 5.
Ingenuity Pathways Analysis (IPA) of DEGs specific to the AP, by comparison with two or more of the other (UM, UC-M, and/or UC-F) models. (.xlsx). Significantly enriched: (A) Pathway and (B) Disease-related GO terms are shown, with the p-value and the DEGs from which each term was derived. (C-G) All genes with differential expression between the AP and the three other models (gene list in C: comparison to UM, UC-M, UC-F) were compared with genes (D-E) that were significant in a transcription-wide association study (TWAS) of autism or that (F-G) had genome-wide significance for de novo protein-coding variants associated with neurodevelopmental disorders, as described in the results text. n = 4; data values are in Additional file 5.
Differentially expressed genes in the UM, by comparison with the UC-M and UC-F control models. (.jpg). (A) Venn diagram shows numbers of differentially expressed genes (DEGs) obtained from pairwise comparisons between the UM vs. the UC-M or UC-F models. UM-specific DEGs are shown in blue. (B-F) These UM-specific DEGs were further analyzed by (B) Hierarchical clustering analysis, with comparisons to all three other models shown, and (C) using Ingenuity Pathway Analysis (IPA), which identified UM-enriched (C) pathways and (D) disease-related GO terms. In C-D, the number of DEGs enriched for each term present is represented on the x-axis, with red and blue colors indicating up- and down-regulated genes, respectively. p-values for each enriched GO term are indicated. (E-F) IPA disease terms enriched in these AP-specific DEGs include gene networks associated with (E) Nervous system development and function and (F) Behavior. The numbers of up-and down-regulated genes present in the networks are indicated. Within each network, red and green symbols indicate up- and down-regulated genes respectively, while color intensity indicates the relative degree of differential expression. n = 4; data values are in Additional file 5.
Ingenuity Pathways Analysis (IPA) of DEGs specific to the UM, by comparison with the UC-M and UC-F models. (.xlsx). Significantly enriched: (A) Pathway and (B) Disease-related GO terms are shown, with the p-value and the DEGs from which each term was derived. n = 4; data values are in Additional file 5.
AP-specific differential gene expression defined by RNA-seq analysis of differentiated cortical neuroids was validated by RT-qPCR and variancePartition analysis. (.jpg). Genes defined as differentially expressed in the AP, by comparison with the UM, by RNA-seq analysis were selected from top AP-enriched gene networks, including axon guidance molecules, integrins, ion channels, and transcription factors, and were validated by RT-qPCR. (A) RPKM values for these DEGs were obtained using RNA-seq analysis (n = 4; data values are in Additional file 5). (B) These were compared with relative gene expression in these models, as defined by RT-qPCR (n = 3; data values are in Additional file 4). RT-qPCR analysis was performed using samples obtained from three independent biological replicate experiments (n = 3) and was performed using a second set of clonal lines derived from the AP and UM that was different than the AP and UM clonal lines used for the RNA-seq analysis. P values *P < 0.05, **P < 0.01, ***P < 0.001 were determined by an unpaired t-test. (C) variancePartition analysis indicates the percent variance that is attributable to individuals from whom the samples were procured, and the study subjects’ age and sex.
Schematic for maturation of cExNs and cINs. (.jpg). (A) Differentiation scheme used to obtain neurons for electrophysiology. cExNPCs and cINPCs, labelled respectively with Synapsin (Syn)-GFP and Syn–RFP, were differentiated in co-culture as cortical neuroids for 15 days, and then further matured by replating on a rat cortical astrocyte feeder layer. (B) MAP2 staining of dissociated cortical neuroids demonstrated that the UM-derived neurons had increased soma size, as quantified in Fig. 7A. (Scale bar = 75 μM). (C) Spontaneous fast excitatory (*) and slower inhibitory (^) postsynaptic currents recorded at -80 mV in a UM-derived cIN. The EPSC is replotted on a faster time base to the right.
Comparison of electrophysiological properties of neurons derived in vitro from the UC-M, AP, and UM iPSC models. (.xlsx). Data is shown for (A) all cells (cExNs and cINs), (B) cExNs, and (C) cINs. (A) Data columns for cells of each genotype show the median, first, and third quartiles, and count (# of cells) for each parameter for all cells recorded, by combining data from cExNs and cINs for each genotype. Columns to the right show neuronal cell type-specific data, by combining cExN data or cIN data for all three genotypes. Light red shading denotes Significant Difference by 2-WAY ANOVA on ranks with post hoc comparison by Student-Newman-Keuls method (p < 0.05) with genotype (UC-M : AP : UM) and neuronal cell type (cExN : cIN) as the two factors. (B-C) Values for (B) cExNs or (C) cINs of each genotype are shown. Tan shading denotes Significant Difference by 1-WAY ANOVA on ranks with post hoc comparison by Student-Newman-Keuls method (p < 0.05), with darker tan indicating that one of the genotypes was different from the other two and lighter tan indicating that two of the three genotypes were different from each other. Capacitance and input resistance were determined from 10 mV voltage steps from -80 mV. First spike threshold, amplitude and half-width were determined for the first action potential recorded at threshold. Maximum first frequency and maximum average firing frequency were determined for 800 msec depolarizing steps that elicited the maximum number of spikes under current clamp. Peak inward sodium current, steady-state outward potassium currents and currents evoked by ACh and choline were recorded under whole-cell voltage clamp. Values are provided.
Examples of genes in selected gene ontology terms described in text. (.pdf). Genes in selected classes are shown, with their classes, associated network and pathways, and physiological roles, and supporting references (cited in main text).
About this article
Cite this article
Meganathan, K., Prakasam, R., Baldridge, D. et al. Altered neuronal physiology, development, and function associated with a common chromosome 15 duplication involving CHRNA7. BMC Biol 19, 147 (2021). https://doi.org/10.1186/s12915-021-01080-7
- Induced pluripotent stem cells
- Cortical neurons
- Psychiatric disease
- Neurodevelopmental disorders
- Copy number variants
- Chromosome 15q13.3 duplication