AAV vector-derived elements integrate into Cas9-generated double-strand breaks and disrupt gene transcription.

We previously developed an adeno-associated virus (AAV) Cas9 gene therapy for Angelman syndrome that integrated into the genome and prematurely terminated Ube3a-ATS. Here, we assessed the performance of 3 additional AAV vectors containing S. aureus Cas9 in vitro and in vivo, and 25 vectors containing N. meningitidis Cas9 in vitro, all targeting single sites within Ube3a-ATS. We found that none of these single-target gRNA vectors were as effective as multi-target gRNA vectors at reducing Ube3a-ATS expression in neurons. We also developed an anchored multiplex PCR sequencing method and analysis pipeline to quantify the relative frequency of all possible editing events at target sites, including AAV integration and unresolved double-strand breaks. We found that integration of AAV was the most frequent editing event (67%-89% of all edits) at three different single target sites, surpassing insertions and deletions (indels). None of the most frequently observed indels were capable of blocking transcription when incorporated into a Ube3a-ATS minigene reporter, whereas two vector derived elements-the poly(A) and reverse promoter-reduced downstream transcription by up to 50%. Our findings suggest that the probability that a gene trapping AAV integration event occurs is influenced by which vector-derived element(s) are integrated and by the number of target sites.

Stimulating Wnt signaling reveals context-dependent genetic effects on gene regulation in primary human neural progenitors.

Gene regulatory effects have been difficult to detect at many non-coding loci associated with brain-related traits, likely because some genetic variants have distinct functions in specific contexts. To explore context-dependent gene regulation, we measured chromatin accessibility and gene expression after activation of the canonical Wnt pathway in primary human neural progenitors (n = 82 donors). We found that TCF/LEF motifs and brain-structure-associated and neuropsychiatric-disorder-associated variants were enriched within Wnt-responsive regulatory elements. Genetically influenced regulatory elements were enriched in genomic regions under positive selection along the human lineage. Wnt pathway stimulation increased detection of genetically influenced regulatory elements/genes by 66%/53% and enabled identification of 397 regulatory elements primed to regulate gene expression. Stimulation also increased identification of shared genetic effects on molecular and complex brain traits by up to 70%, suggesting that genetic variant function during neurodevelopmental patterning can lead to differences in adult brain and behavioral traits.

Autism-linked UBE3A gain-of-function mutation causes interneuron and behavioral phenotypes when inherited maternally or paternally in mice.

The E3 ubiquitin ligase Ube3a is biallelically expressed in neural progenitors and glial cells, suggesting that UBE3A gain-of-function mutations might cause neurodevelopmental disorders irrespective of parent of origin. Here, we engineered a mouse line that harbors an autism-linked UBE3AT485A (T503A in mouse) gain-of-function mutation and evaluated phenotypes in animals that inherited the mutant allele paternally, maternally, or from both parents. We find that paternally and maternally expressed UBE3AT503A results in elevated UBE3A activity in neural progenitors and glial cells. Expression of UBE3AT503A from the maternal allele, but not the paternal one, leads to a persistent elevation of UBE3A activity in neurons. Mutant mice display behavioral phenotypes that differ by parent of origin. Expression of UBE3AT503A, irrespective of its parent of origin, promotes transient embryonic expansion of Zcchc12 lineage interneurons. Phenotypes of Ube3aT503A mice are distinct from Angelman syndrome model mice. Our study has clinical implications for a growing number of disease-linked UBE3A gain-of-function mutations.

Standards of Fluid Biomarker Collection and Pre-analytical Processes in Humans and Mice: Recommendations by the Ataxia Global Initiative Working Group on Biomarkers.

The Ataxia Global Initiative (AGI) aims to serve as a platform to facilitate clinical trial readiness for the hereditary ataxias. Clinical trials for these diseases have been hampered by the lack of objective measures to study disease onset, progression, and treatment efficacy. While these issues are not unique to the genetic ataxias, the relative rarity of these diseases makes the need for such measures even more pressing to achieve statistical power in clinical trials. In this report, we have described the efforts of the AGI fluid biomarker working group (WG) in developing uniform protocols for biomarker sampling and storage, both for human and preclinical studies in mice. By reducing collection variability, we anticipate reduced noise in downstream biomarker analysis that will improve statistical power and minimize the necessary sample size. The emphasis has been on defining and standardizing the sampling and pre-analytical work-up of minimal set of biological samples, specifically blood plasma and serum, keeping in mind the need for harmonization of collection and storage that can be achieved with relatively limited cost and resources. An optional package is detailed for those centers that have the resources and commitment for additional biofluids/sample processing and storage. Finally, we have delineated similar standardized protocols for mice that will be important for preclinical studies in the field.

Wnt activity reveals context-specific genetic effects on gene regulation in neural progenitors.

Gene regulatory effects in bulk-post mortem brain tissues are undetected at many non-coding brain trait-associated loci. We hypothesized that context-specific genetic variant function during stimulation of a developmental signaling pathway would explain additional regulatory mechanisms. We measured chromatin accessibility and gene expression following activation of the canonical Wnt pathway in primary human neural progenitors from 82 donors. TCF/LEF motifs, brain structure-, and neuropsychiatric disorder-associated variants were enriched within Wnt-responsive regulatory elements (REs). Genetically influenced REs were enriched in genomic regions under positive selection along the human lineage. Stimulation of the Wnt pathway increased the detection of genetically influenced REs/genes by 66.2%/52.7%, and led to the identification of 397 REs primed for effects on gene expression. Context-specific molecular quantitative trait loci increased brain-trait colocalizations by up to 70%, suggesting that genetic variant effects during early neurodevelopmental patterning lead to differences in adult brain and behavioral traits.

MicroRNA-eQTLs in the developing human neocortex link miR-4707-3p expression to brain size.

Expression quantitative trait loci (eQTL) data have proven important for linking non-coding loci to protein-coding genes. But eQTL studies rarely measure microRNAs (miRNAs), small non-coding RNAs known to play a role in human brain development and neurogenesis. Here, we performed small-RNA sequencing across 212 mid-gestation human neocortical tissue samples, measured 907 expressed miRNAs, discovering 111 of which were novel, and identified 85 local-miRNA-eQTLs. Colocalization of miRNA-eQTLs with GWAS summary statistics yielded one robust colocalization of miR-4707-3p expression with educational attainment and brain size phenotypes, where the miRNA expression increasing allele was associated with decreased brain size. Exogenous expression of miR-4707-3p in primary human neural progenitor cells decreased expression of predicted targets and increased cell proliferation, indicating miR-4707-3p modulates progenitor gene regulation and cell fate decisions. Integrating miRNA-eQTLs with existing GWAS yielded evidence of a miRNA that may influence human brain size and function via modulation of neocortical brain development.

Cellular Genome-wide Association Study Identifies Common Genetic Variation Influencing Lithium-Induced Neural Progenitor Proliferation.

BACKGROUND: Bipolar disorder is a highly heritable neuropsychiatric condition affecting more than 1% of the human population. Lithium salts are commonly prescribed as a mood stabilizer for individuals with bipolar disorder. Lithium is clinically effective in approximately half of treated individuals, and their genetic backgrounds are known to influence treatment outcomes. While the mechanism of lithium's therapeutic action is unclear, it stimulates adult neural progenitor cell proliferation, similar to some antidepressant drugs. METHODS: To identify common genetic variants that modulate lithium-induced proliferation, we conducted an EdU incorporation assay in a library of 80 genotyped human neural progenitor cells treated with lithium. These data were used to perform a genome-wide association study to identify common genetic variants that influence lithium-induced neural progenitor cell proliferation. We manipulated the expression of a putatively causal gene using CRISPRi/a (clustered regularly interspaced short palindromic repeats interference/activation) constructs to experimentally verify lithium-induced proliferation effects. RESULTS: We identified a locus on chr3p21.1 associated with lithium-induced proliferation. This locus is also associated with bipolar disorder risk, schizophrenia risk, and interindividual differences in intelligence. We identified a single gene, GNL3, whose expression temporally increased in an allele-specific fashion following lithium treatment. Experimentally increasing the expression of GNL3 led to increased proliferation under baseline conditions, while experimentally decreasing GNL3 expression suppressed lithium-induced proliferation. CONCLUSIONS: Our experiments reveal that common genetic variation modulates lithium-induced neural progenitor proliferation and that GNL3 expression is necessary for the full proliferation-stimulating effects of lithium. These results suggest that performing genome-wide associations in genetically diverse human cell lines is a useful approach to discover context-specific pharmacogenomic effects.

Multi-omic profiling reveals the ataxia protein sacsin is required for integrin trafficking and synaptic organization.

Autosomal recessive spastic ataxia of Charlevoix-Saguenay (ARSACS) is a childhood-onset cerebellar ataxia caused by mutations in SACS, which encodes the protein sacsin. Cellular ARSACS phenotypes include mitochondrial dysfunction, intermediate filament disorganization, and progressive death of cerebellar Purkinje neurons. It is unclear why the loss of sacsin causes these deficits or why they manifest as cerebellar ataxia. Here, we perform multi-omic profiling in sacsin knockout (KO) cells and identify alterations in microtubule dynamics and mislocalization of focal adhesion (FA) proteins, including multiple integrins. Deficits in FA structure, signaling, and function can be rescued by targeting PTEN, a negative regulator of FA signaling. ARSACS mice possess mislocalization of ITGA1 in Purkinje neurons and synaptic disorganization in the deep cerebellar nucleus (DCN). The sacsin interactome reveals that sacsin regulates interactions between cytoskeletal and synaptic adhesion proteins. Our findings suggest that disrupted trafficking of synaptic adhesion proteins is a causal molecular deficit in ARSACS.

ToxCast chemical library Wnt screen identifies diethanolamine as an activator of neural progenitor proliferation.

Numerous autism spectrum disorder (ASD) risk genes are associated with Wnt signaling, suggesting that brain development may be especially sensitive to genetic perturbation of this pathway. Additionally, valproic acid, which modulates Wnt signaling, increases risk for ASD when taken during pregnancy. We previously found that an autism-linked gain-of-function UBE3A T485A mutant construct hyperactivated canonical Wnt signaling, providing a genetic means to elevate Wnt signaling above baseline levels. To identify environmental use chemicals that enhance or suppress Wnt signaling, we screened the ToxCast Phase I and II libraries in cells expressing this autism-linked UBE3A T485A gain-of-function mutant construct. Using structural comparisons, we identify classes of chemicals that stimulated Wnt signaling, including ethanolamines, as well as chemicals that inhibited Wnt signaling, such as agricultural pesticides, and synthetic hormone analogs. To prioritize chemicals for follow-up, we leveraged predicted human exposure data, and identified diethanolamine (DEA) as a chemical that stimulates Wnt signaling in UBE3A T485A -transfected cells, and has a high potential for prenatal exposure in humans. DEA enhanced proliferation in primary human neural progenitor cell lines (phNPC), but did not affect expression of canonical Wnt target genes in NPCs or primary mouse neuron cultures. Instead, we found DEA increased expression of the H3K9 methylation sensitive gene CALB1, consistent with competitive inhibition of the methyl donor enzymatic pathways.

Cell-type-specific effects of genetic variation on chromatin accessibility during human neuronal differentiation.

Common genetic risk for neuropsychiatric disorders is enriched in regulatory elements active during cortical neurogenesis. However, it remains poorly understood as to how these variants influence gene regulation. To model the functional impact of common genetic variation on the noncoding genome during human cortical development, we performed the assay for transposase accessible chromatin using sequencing (ATAC-seq) and analyzed chromatin accessibility quantitative trait loci (QTL) in cultured human neural progenitor cells and their differentiated neuronal progeny from 87 donors. We identified significant genetic effects on 988/1,839 neuron/progenitor regulatory elements, with highly cell-type and temporally specific effects. A subset (roughly 30%) of chromatin accessibility-QTL were also associated with changes in gene expression. Motif-disrupting alleles of transcriptional activators generally led to decreases in chromatin accessibility, whereas motif-disrupting alleles of repressors led to increases in chromatin accessibility. By integrating cell-type-specific chromatin accessibility-QTL and brain-relevant genome-wide association data, we were able to fine-map and identify regulatory mechanisms underlying noncoding neuropsychiatric disorder risk loci.

Cas9 gene therapy for Angelman syndrome traps Ube3a-ATS long non-coding RNA.

Angelman syndrome (AS) is a severe neurodevelopmental disorder caused by a mutation or deletion of the maternally inherited UBE3A allele. In neurons, the paternally inherited UBE3A allele is silenced in cis by a long non-coding RNA called UBE3A-ATS. Here, as part of a systematic screen, we found that Cas9 can be used to activate ('unsilence') paternal Ube3a in cultured mouse and human neurons when targeted to Snord115 genes, which are small nucleolar RNAs that are clustered in the 3' region of Ube3a-ATS. A short Cas9 variant and guide RNA that target about 75 Snord115 genes were packaged into an adeno-associated virus and administered to a mouse model of AS during the embryonic and early postnatal stages, when the therapeutic benefit of restoring Ube3a is predicted to be greatest1,2. This early treatment unsilenced paternal Ube3a throughout the brain for at least 17 months and rescued anatomical and behavioural phenotypes in AS mice. Genomic integration of the adeno-associated virus vector into Cas9 target sites caused premature termination of Ube3a-ATS at the vector-derived polyA cassette, or when integrated in the reverse orientation, by transcriptional collision with the vector-derived Cas9 transcript. Our study shows that targeted genomic integration of a gene therapy vector can restore the function of paternally inherited UBE3A throughout life, providing a path towards a disease-modifying treatment for a syndromic neurodevelopmental disorder.

High-throughput screening and classification of chemicals and their effects on neuronal gene expression using RASL-seq.

We previously used RNA-seq to identify chemicals whose effects on neuronal gene expression mimicked transcriptional signatures of autism, aging, and neurodegeneration. However, this approach was costly and time consuming, which limited our study to testing a single chemical concentration on mixed sex cortical neuron cultures. Here, we adapted a targeted transcriptomic method (RASL-seq, similar to TempO-seq) to interrogate changes in expression of a set of 56 signature genes in response to a library of 350 chemicals and chemical mixtures at four concentrations in male and female mouse neuronal cultures. This enabled us to replicate and expand our previous classifications, and show that transcriptional responses were largely equivalent between sexes. Overall, we found that RASL-seq can be used to accelerate the pace at which chemicals and mixtures that transcriptionally mimic autism and other neuropsychiatric diseases can be identified, and provides a cost-effective way to quantify gene expression with a panel of marker genes.

The autism-linked UBE3A T485A mutant E3 ubiquitin ligase activates the Wnt/ß-catenin pathway by inhibiting the proteasome.

UBE3A is a HECT domain E3 ubiquitin ligase whose dysfunction is linked to autism, Angelman syndrome, and cancer. Recently, we characterized a de novo autism-linked UBE3A mutant (UBE3AT485A) that disrupts phosphorylation control of UBE3A activity. Through quantitative proteomics and reporter assays, we found that the UBE3AT485A protein ubiquitinates multiple proteasome subunits, reduces proteasome subunit abundance and activity, stabilizes nuclear ß-catenin, and stimulates canonical Wnt signaling more effectively than wild-type UBE3A. We also found that UBE3AT485A activates Wnt signaling to a greater extent in cells with low levels of ongoing Wnt signaling, suggesting that cells with low basal Wnt activity are particularly vulnerable to UBE3AT485A mutation. Ligase-dead UBE3A did not stimulate Wnt pathway activation. Overexpression of several proteasome subunits reversed the effect of UBE3AT485A on Wnt signaling. We also observed that subunits that interact with UBE3A and affect Wnt signaling are located along one side of the 19S regulatory particle, indicating a previously unrecognized spatial organization to the proteasome. Altogether, our findings indicate that UBE3A regulates Wnt signaling in a cell context-dependent manner and that an autism-linked mutation exacerbates these signaling effects. Our study has broad implications for human disorders associated with UBE3A gain or loss of function and suggests that dysfunctional UBE3A might affect additional proteins and pathways that are sensitive to proteasome activity.

Evolutionary patterns of metazoan microRNAs reveal targeting principles in the let-7 and miR-10 families.

MicroRNAs (miRNAs) regulate gene output by targeting degenerate elements in mRNAs and have undergone drastic expansions in higher metazoan genomes. The evolutionary advantage of maintaining copies of highly similar miRNAs is not well understood, nor is it clear what unique functions, if any, miRNA family members possess. Here, we study evolutionary patterns of metazoan miRNAs, focusing on the targeting preferences of the let-7 and miR-10 families. These studies reveal hotspots for sequence evolution with implications for targeting and secondary structure. High-throughput screening for functional targets reveals that each miRNA represses sites with distinct features and regulates a large number of genes with cooperative function in regulatory networks. Unexpectedly, given the high degree of similarity, single-nucleotide changes grant miRNA family members with distinct targeting preferences. Together, our data suggest complex functional relationships among miRNA duplications, novel expression patterns, sequence change, and the acquisition of new targets.

Differential expression of conserved and novel microRNAs during tail regeneration in the lizard Anolis carolinensis.

BACKGROUND: Lizards are evolutionarily the most closely related vertebrates to humans that can lose and regrow an entire appendage. Regeneration in lizards involves differential expression of hundreds of genes that regulate wound healing, musculoskeletal development, hormonal response, and embryonic morphogenesis. While microRNAs are able to regulate large groups of genes, their role in lizard regeneration has not been investigated. RESULTS: MicroRNA sequencing of green anole lizard (Anolis carolinensis) regenerating tail and associated tissues revealed 350 putative novel and 196 known microRNA precursors. Eleven microRNAs were differentially expressed between the regenerating tail tip and base during maximum outgrowth (25 days post autotomy), including miR-133a, miR-133b, and miR-206, which have been reported to regulate regeneration and stem cell proliferation in other model systems. Three putative novel differentially expressed microRNAs were identified in the regenerating tail tip. CONCLUSIONS: Differentially expressed microRNAs were identified in the regenerating lizard tail, including known regulators of stem cell proliferation. The identification of 3 putative novel microRNAs suggests that regulatory networks, either conserved in vertebrates and previously uncharacterized or specific to lizards, are involved in regeneration. These findings suggest that differential regulation of microRNAs may play a role in coordinating the timing and expression of hundreds of genes involved in regeneration.

A human 3'UTR clone collection to study post-transcriptional gene regulation.

BACKGROUND: 3'untranslated regions (3'UTRs) are poorly understood portions of eukaryotic mRNAs essential for post-transcriptional gene regulation. Sequence elements in 3'UTRs can be target sites for regulatory molecules such as RNA binding proteins and microRNAs (miRNAs), and these interactions can exert significant control on gene networks. However, many such interactions remain uncharacterized due to a lack of high-throughput (HT) tools to study 3'UTR biology. HT cloning efforts such as the human ORFeome exemplify the potential benefits of genomic repositories for studying human disease, especially in relation to the discovery of biomarkers and targets for therapeutic agents. Currently there are no publicly available human 3'UTR libraries. To address this we have prepared the first version of the human 3'UTRome (h3'UTRome v1) library. The h3'UTRome is produced to a single high quality standard using the same recombinational cloning technology used for the human ORFeome, enabling universal operating methods and high throughput experimentation. The library is thoroughly sequenced and annotated with simple online access to information, and made publically available through gene repositories at low cost to all scientists with minimal restriction. RESULTS: The first release of the h3'UTRome library comprises 1,461 human 3'UTRs cloned into Gateway® entry vectors, ready for downstream analyses. It contains 3'UTRs for 985 transcription factors, 156 kinases, 171 RNA binding proteins, and 186 other genes involved in gene regulation and in disease. We demonstrate the feasibility of the h3'UTRome library by screening a panel of 87 3'UTRs for targeting by two miRNAs: let-7c, which is implicated in tumorigenesis, and miR-221, which is implicated in atherosclerosis and heart disease. The panel is enriched with genes involved in the RAS signaling pathway, putative novel targets for the two miRNAs, as well as genes implicated in tumorigenesis and heart disease. CONCLUSIONS: The h3'UTRome v1 library is a modular resource that can be utilized for high-throughput screens to identify regulatory interactions between trans-acting factors and 3'UTRs, Importantly, the library can be customized based on the specifications of the researcher, allowing the systematic study of human 3'UTR biology.

Detection of miRNA Targets in High-throughput Using the 3'LIFE Assay.

Luminescent Identification of Functional Elements in 3'UTRs (3'LIFE) allows the rapid identification of targets of specific miRNAs within an array of hundreds of queried 3'UTRs. Target identification is based on the dual-luciferase assay, which detects binding at the mRNA level by measuring translational output, giving a functional readout of miRNA targeting. 3'LIFE uses non-proprietary buffers and reagents, and publically available reporter libraries, making genome-wide screens feasible and cost-effective. 3'LIFE can be performed either in a standard lab setting or scaled up using liquid handling robots and other high-throughput instrumentation. We illustrate the approach using a dataset of human 3'UTRs cloned in 96-well plates, and two test miRNAs, let-7c and miR-10b. We demonstrate how to perform DNA preparation, transfection, cell culture and luciferase assays in 96-well format, and provide tools for data analysis. In conclusion 3'LIFE is highly reproducible, rapid, systematic, and identifies high confidence targets.

3'LIFE: a functional assay to detect miRNA targets in high-throughput.

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene output at the post-transcriptional level by targeting degenerate elements primarily in 3'untranslated regions (3'UTRs) of mRNAs. Individual miRNAs can regulate networks of hundreds of genes, yet for the majority of miRNAs few, if any, targets are known. Misexpression of miRNAs is also a major contributor to cancer progression, thus there is a critical need to validate miRNA targets in high-throughput to understand miRNAs' contribution to tumorigenesis. Here we introduce a novel high-throughput assay to detect miRNA targets in 3'UTRs, called Luminescent Identification of Functional Elements in 3'UTRs (3'LIFE). We demonstrate the feasibility of 3'LIFE using a data set of 275 human 3'UTRs and two cancer-relevant miRNAs, let-7c and miR-10b, and compare our results to alternative methods to detect miRNA targets throughout the genome. We identify a large number of novel gene targets for these miRNAs, with only 32% of hits being bioinformatically predicted and 27% directed by non-canonical interactions. Functional analysis of target genes reveals consistent roles for each miRNA as either a tumor suppressor (let-7c) or oncogenic miRNA (miR-10b), and preferentially target multiple genes within regulatory networks, suggesting 3'LIFE is a rapid and sensitive method to detect miRNA targets in high-throughput.