Phosphatase and tensin homolog (PTEN) is a major negative regulator of the phosphatidylinositol-3-kinase (PI3K)/Akt/mechanistic target of rapamycin (mTOR) pathway. Loss-of-function mutations in PTEN have been found in a subset of patients with macrocephaly and autism spectrum disorder (ASD). PTEN loss in neurons leads to somal hypertrophy, aberrant migration, dendritic overgrowth, increased spine density, and hyperactivity of neuronal circuits. These neuronal overgrowth phenotypes are present on Pten knock-out (KO) and reconstitution with autism-associated point mutations. The mechanism underlying dendritic overgrowth in Pten deficient neurons is unclear. In this study, we examined how Pten loss impacts microtubule (MT) dynamics in both sexes using retroviral infection and transfection strategies to manipulate PTEN expression and tag the plus-end MT binding protein, end-binding protein 3 (EB3). We found Pten KO neurons sprout more new processes over time compared with wild-type (WT) neurons. We also found an increase in MT polymerization rate in Pten KO dendritic growth cones. Reducing MT polymerization rate to the WT level was sufficient to reduce dendritic overgrowth in Pten KO neurons in vitro and in vivo Finally, we found that rescue of dendritic overgrowth via inhibition of MT polymerization was sufficient to improve the performance of Pten KO mice in a spatial memory task. Taken together, our data suggests that one factor underlying PTEN loss dependent dendritic overgrowth is increased MT polymerization. This opens the possibility for an intersectional approach targeting MT polymerization and mTOR with low doses of inhibitors to achieve therapeutic gains with minimal side effects in pathologies associated with loss of neuronal PTEN function.SIGNIFICANCE STATEMENT Loss of Pten function because of genetic deletion or expression of mutations associated with autism spectrum disorder (ASD), results in overgrowth of neurons including increased total dendritic length and branching. We have discovered that this overgrowth is accompanied by increased rate of microtubule (MT) polymerization. The increased polymerization rate is insensitive to acute inhibition of mechanistic target of rapamycin (mTOR)C1 or protein synthesis. Direct pharmacological inhibition of MT polymerization can slow the polymerization rate in Pten knock-out (KO) neurons to rates seen in wild-type (WT) neurons. Correction of the MT polymerization rate rescues increased total dendritic arborization and spatial memory. Our studies suggest that phosphatase and tensin homolog (PTEN) inhibits dendritic growth through parallel regulation of protein synthesis and cytoskeletal polymerization.

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWAS) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N=52) and nonsmokers (N=171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and using a two-stage multiple testing approach with eigenMT and Bonferroni corrections. We found >2 million significant meQTL variants (padj<0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects; five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTLs for 958 unique eGenes (padj<0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTIN1 and ITIH4 colocalized across all data types (GWAS + meQTL + eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Cigarette smoking during pregnancy is a major public health concern. While there are well-described consequences in early child development, there is very little known about the effects of maternal smoking on human cortical biology during prenatal life. We therefore performed a genome-wide differential gene expression analysis using RNA sequencing (RNA-seq) on prenatal (N = 33; 16 smoking-exposed) as well as adult (N = 207; 57 active smokers) human postmortem prefrontal cortices. Smoking exposure during the prenatal period was directly associated with differential expression of 14 genes; in contrast, during adulthood, despite a much larger sample size, only two genes showed significant differential expression (FDR < 10%). Moreover, 1,315 genes showed significantly different exposure effects between maternal smoking during pregnancy and direct exposure in adulthood (FDR < 10%)-these differences were largely driven by prenatal differences that were enriched for pathways previously implicated in addiction and synaptic function. Furthermore, prenatal and age-dependent differentially expressed genes were enriched for genes implicated in non-syndromic autism spectrum disorder (ASD) and were differentially expressed as a set between patients with ASD and controls in postmortem cortical regions. These results underscore the enhanced sensitivity to the biological effect of smoking exposure in the developing brain and offer insight into how maternal smoking during pregnancy affects gene expression in the prenatal human cortex. They also begin to address the relationship between in utero exposure to smoking and the heightened risks for the subsequent development of neuropsychiatric disorders.

Epigenome-wide association studies of human disease and other quantitative traits are becoming increasingly common. A series of papers reporting age-related changes in DNA methylation profiles in peripheral blood have already been published. However, blood is a heterogeneous collection of different cell types, each with a very different DNA methylation profile.

DNA methylation levels change with age. Recent studies have identified biomarkers of chronological age based on DNA methylation levels. It is not yet known whether DNA methylation age captures aspects of biological age.

Retroviruses expressing a fluorescent protein, Cas9, and a small guide RNA are used to mimic nonsense PTEN mutations from autism patients in developing mouse neurons. We compare the cellular phenotype elicited by CRISPR-Cas9 to those elicited using shRNA or Cre/Lox technologies and find that knockdown or knockout (KO) produced a corresponding moderate or severe neuronal hypertrophy in all cells. In contrast, the Cas9 approach produced missense and nonsense Pten mutations, resulting in a mix of KO-equivalent hypertrophic and wild type-like phenotypes. Importantly, despite this mixed phenotype, the neuronal hypertrophy resulting from Pten loss was evident on average in the population of manipulated cells. Having reproduced the known Pten KO phenotype using the CRISPR-Cas9 system we design viruses to target a gene that has recently been associated with autism, KATNAL2. Katnal2 deletion in the mouse results in decreased dendritic arborization of developing neurons. We conclude that retroviral implementation of the CRISPR-Cas9 system is an efficient system for cellular phenotype discovery in wild-type animals.

The primate-specific brain voltage-gated potassium channel isoform Kv11.1-3.1 has been identified as a novel therapeutic target for the treatment of schizophrenia. While this ether-a-go-go related K(+)channel has shown clinical relevance, drug discovery efforts have been hampered due to low and inconsistent activity in cell-based assays. This poor activity is hypothesized to result from poor trafficking via the lack of an intact channel-stabilizing Per-Ant-Sim (PAS) domain. Here we characterize Kv11.1-3.1 cellular localization and show decreased channel expression and cell surface trafficking relative to the PAS-domain containing major isoform, Kv11.1-1A. Using small molecule inhibition of proteasome degradation, cellular expression and plasma membrane trafficking are rescued. These findings implicate the importance of the unfolded-protein response and endoplasmic reticulum associated degradation pathways in the expression and regulation of this schizophrenia risk factor. Utilizing this identified phenomenon, an electrophysiological and high throughput in-vitro fluorescent assay platform has been developed for drug discovery in order to explore a potentially new class of cognitive therapeutics.

Neurogenesis is a highly-regulated process occurring in the dentate gyrus that has been linked to learning, memory, and antidepressant efficacy. MicroRNAs (miRNAs) have been previously shown to play an important role in the regulation of neuronal development and neurogenesis in the dentate gyrus via modulation of gene expression. However, this mode of regulation is both incompletely described in the literature thus far and highly multifactorial. In this study, we designed sensors and detected relative levels of expression of 10 different miRNAs and found miR-338-3p was most highly expressed in the dentate gyrus. Comparison of miR-338-3p expression with neuronal markers of maturity indicates miR-338-3p is expressed most highly in the mature neuron. We also designed a viral "sponge" to knock down in vivo expression of miR-338-3p. When miR-338-3p is knocked down, neurons sprout multiple primary dendrites that branch off of the soma in a disorganized manner, cellular proliferation is upregulated, and neoplasms form spontaneously in vivo. Additionally, miR-338-3p overexpression in glioblastoma cell lines slows their proliferation in vitro. Further, low miR-338-3p expression is associated with increased mortality and disease progression in patients with glioblastoma. These data identify miR-338-3p as a clinically relevant tumor suppressor in glioblastoma.

Genome-wide association studies have identified 108 schizophrenia risk loci, but biological mechanisms for individual loci are largely unknown. Using developmental, genetic and illness-based RNA sequencing expression analysis in human brain, we characterized the human brain transcriptome around these loci and found enrichment for developmentally regulated genes with novel examples of shifting isoform usage across pre- and postnatal life. We found widespread expression quantitative trait loci (eQTLs), including many with transcript specificity and previously unannotated sequence that were independently replicated. We leveraged this general eQTL database to show that 48.1% of risk variants for schizophrenia associate with nearby expression. We lastly found 237 genes significantly differentially expressed between patients and controls, which replicated in an independent dataset, implicated synaptic processes, and were strongly regulated in early development. These findings together offer genetics- and diagnosis-related targets for better modeling of schizophrenia risk. This resource is publicly available at http://eqtl.brainseq.org/phase1 .

Significance analysis plays a major role in identifying and ranking genes, transcription factor binding sites, DNA methylation regions, and other high-throughput features associated with illness. We propose a new approach, called gene set bagging, for measuring the probability that a gene set replicates in future studies. Gene set bagging involves resampling the original high-throughput data, performing gene-set analysis on the resampled data, and confirming that biological categories replicate in the bagged samples.

The dentate gyrus (DG) is a region of the adult rodent brain that undergoes continuous neurogenesis. Seizures and loss or dysfunction of GABAergic synapses onto adult-born dentate granule cells (GCs) alter their dendritic growth and migration, resulting in dysmorphic and hyperexcitable GCs. Additionally, transplants of fetal GABAergic interneurons in the DG of mice with temporal lobe epilepsy (TLE) result in seizure suppression, but it is unknown whether increasing interneurons with these transplants restores GABAergic innervation to adult-born GCs. Here, we address this question by birth-dating GCs with retrovirus at different times up to 12 weeks after pilocarpine-induced TLE in adult mice. Channelrhodopsin 2 (ChR2)-enhanced yellow fluorescent protein (EYFP)-expressing medial-ganglionic eminence (MGE)-derived GABAergic interneurons from embryonic day (E)13.5 mouse embryos were transplanted into the DG of the TLE mice and GCs with transplant-derived inhibitory post-synaptic currents (IPSCs) were identified by patch-clamp electrophysiology and optogenetic interrogation. Putative synaptic sites between GCs and GABAergic transplants were also confirmed by intracellular biocytin staining, immunohistochemistry, and confocal imaging. 3D reconstructions of dendritic arbors and quantitative morphometric analyses were carried out in >150 adult-born GCs. GABAergic inputs from transplanted interneurons correlated with markedly shorter GC dendrites, compared to GCs that were not innervated by the transplants. Moreover, these effects were confined to distal dendritic branches and a short time window of six to eight weeks. The effects were independent of seizures as they were also observed in naïve mice with MGE transplants. These findings are consistent with the hypothesis that increased inhibitory currents over a smaller dendritic arbor in adult-born GCs may reduce their excitability and lead to seizure suppression.

DNA methylation (DNAm) is a critical regulator of both development and cellular identity and shows unique patterns in neurons. To better characterize maturational changes in DNAm patterns in these cells, we profile the DNAm landscape at single-base resolution across the first two decades of human neocortical development in NeuN+ neurons using whole-genome bisulfite sequencing and compare them to non-neurons (primarily glia) and prenatal homogenate cortex.

We present recount3, a resource consisting of over 750,000 publicly available human and mouse RNA sequencing (RNA-seq) samples uniformly processed by our new Monorail analysis pipeline. To facilitate access to the data, we provide the recount3 and snapcount R/Bioconductor packages as well as complementary web resources. Using these tools, data can be downloaded as study-level summaries or queried for specific exon-exon junctions, genes, samples, or other features. Monorail can be used to process local and/or private data, allowing results to be directly compared to any study in recount3. Taken together, our tools help biologists maximize the utility of publicly available RNA-seq data, especially to improve their understanding of newly collected data. recount3 is available from http://rna.recount.bio .

Single-cell gene expression technologies are powerful tools to study cell types in the human brain, but efforts have largely focused on cortical brain regions. We therefore created a single-nucleus RNA-sequencing resource of 70,615 high-quality nuclei to generate a molecular taxonomy of cell types across five human brain regions that serve as key nodes of the human brain reward circuitry: nucleus accumbens, amygdala, subgenual anterior cingulate cortex, hippocampus, and dorsolateral prefrontal cortex. We first identified novel subpopulations of interneurons and medium spiny neurons (MSNs) in the nucleus accumbens and further characterized robust GABAergic inhibitory cell populations in the amygdala. Joint analyses across the 107 reported cell classes revealed cell-type substructure and unique patterns of transcriptomic dynamics. We identified discrete subpopulations of D1- and D2-expressing MSNs in the nucleus accumbens to which we mapped cell-type-specific enrichment for genetic risk associated with both psychiatric disease and addiction.

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

DNA methylation (DNAm) is an epigenetic regulator of gene expression and a hallmark of gene-environment interaction. Using whole-genome bisulfite sequencing, we have surveyed DNAm in 344 samples of human postmortem brain tissue from neurotypical subjects and individuals with schizophrenia. We identify genetic influence on local methylation levels throughout the genome, both at CpG sites and CpH sites, with 86% of SNPs and 55% of CpGs being part of methylation quantitative trait loci (meQTLs). These associations can further be clustered into regions that are differentially methylated by a given SNP, highlighting the genes and regions with which these loci are epigenetically associated. These findings can be used to better characterize schizophrenia GWAS-identified variants as epigenetic risk variants. Regions differentially methylated by schizophrenia risk-SNPs explain much of the heritability associated with risk loci, despite covering only a fraction of the genomic space. We provide a comprehensive, single base resolution view of association between genetic variation and genomic methylation, and implicate schizophrenia GWAS-associated variants as influencing the epigenetic plasticity of the brain.

DNA methylation (DNAm) is important in brain development and is potentially important in schizophrenia. We characterized DNAm in prefrontal cortex from 335 non-psychiatric controls across the lifespan and 191 patients with schizophrenia and identified widespread changes in the transition from prenatal to postnatal life. These DNAm changes manifest in the transcriptome, correlate strongly with a shifting cellular landscape and overlap regions of genetic risk for schizophrenia. A quarter of published genome-wide association studies (GWAS)-suggestive loci (4,208 of 15,930, P < 10(-100)) manifest as significant methylation quantitative trait loci (meQTLs), including 59.6% of GWAS-positive schizophrenia loci. We identified 2,104 CpGs that differ between schizophrenia patients and controls that were enriched for genes related to development and neurodifferentiation. The schizophrenia-associated CpGs strongly correlate with changes related to the prenatal-postnatal transition and show slight enrichment for GWAS risk loci while not corresponding to CpGs differentiating adolescence from later adult life. These data implicate an epigenetic component to the developmental origins of this disorder.

regionReport is an R package for generating detailed interactive reports from region-level genomic analyses as well as feature-level RNA-seq. The report includes quality-control checks, an overview of the results, an interactive table of the genomic regions or features of interest and reproducibility information. regionReport provides specialised reports for exploring DESeq2, edgeR, or derfinder differential expression analyses results. regionReport is also flexible and can easily be expanded with report templates for other analysis pipelines.

Using a functional approach to investigate the epigenetics of type 2 diabetes (T2D), we combine three lines of evidence-diet-induced epigenetic dysregulation in mouse, epigenetic conservation in humans, and T2D clinical risk evidence-to identify genes implicated in T2D pathogenesis through epigenetic mechanisms related to obesity. Beginning with dietary manipulation of genetically homogeneous mice, we identify differentially DNA-methylated genomic regions. We then replicate these results in adipose samples from lean and obese patients pre- and post-Roux-en-Y gastric bypass, identifying regions where both the location and direction of methylation change are conserved. These regions overlap with 27 genetic T2D risk loci, only one of which was deemed significant by GWAS alone. Functional analysis of genes associated with these regions revealed four genes with roles in insulin resistance, demonstrating the potential general utility of this approach for complementing conventional human genetic studies by integrating cross-species epigenomics and clinical genetic risk.

Long noncoding RNAs (lncRNAs) have emerged as key coordinators of biological and cellular processes. Characterizing lncRNA expression across cells and tissues is key to understanding their role in determining phenotypes, including human diseases. We present here FC-R2, a comprehensive expression atlas across a broadly defined human transcriptome, inclusive of over 109,000 coding and noncoding genes, as described in the FANTOM CAGE-Associated Transcriptome (FANTOM-CAT) study. This atlas greatly extends the gene annotation used in the original recount2 resource. We demonstrate the utility of the FC-R2 atlas by reproducing key findings from published large studies and by generating new results across normal and diseased human samples. In particular, we (a) identify tissue-specific transcription profiles for distinct classes of coding and noncoding genes, (b) perform differential expression analysis across thirteen cancer types, identifying novel noncoding genes potentially involved in tumor pathogenesis and progression, and (c) confirm the prognostic value for several enhancer lncRNAs expression in cancer. Our resource is instrumental for the systematic molecular characterization of lncRNA by the FANTOM6 Consortium. In conclusion, comprised of over 70,000 samples, the FC-R2 atlas will empower other researchers to investigate functions and biological roles of both known coding genes and novel lncRNAs.