Recent advances in 3D culture systems have led to the generation of brain organoids that resemble different human brain regions; however, a 3D organoid model of the midbrain containing functional midbrain dopaminergic (mDA) neurons has not been reported. We developed a method to differentiate human pluripotent stem cells into a large multicellular organoid-like structure that contains distinct layers of neuronal cells expressing characteristic markers of human midbrain. Importantly, we detected electrically active and functionally mature mDA neurons and dopamine production in our 3D midbrain-like organoids (MLOs). In contrast to human mDA neurons generated using 2D methods or MLOs generated from mouse embryonic stem cells, our human MLOs produced neuromelanin-like granules that were structurally similar to those isolated from human substantia nigra tissues. Thus our MLOs bearing features of the human midbrain may provide a tractable in vitro system to study the human midbrain and its related diseases.

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.

A major goal of drug abuse research is to identify and understand drug-induced changes in brain function that are common to many or all drugs of abuse. As these may underlie drug dependence and addiction, the purpose of the present study was to examine if different drugs of abuse effect changes in gene expression that converge in common molecular pathways. Microarray analysis was employed to assay brain gene expression in postmortem anterior prefrontal cortex (aPFC) from 42 human cocaine, cannabis and/or phencyclidine abuse cases and 30 control cases, which were characterized by toxicology and drug abuse history. Common transcriptional changes were demonstrated for a majority of drug abuse cases (N = 34), representing a number of consistently changed functional classes: Calmodulin-related transcripts (CALM1, CALM2, CAMK2B) were decreased, while transcripts related to cholesterol biosynthesis and trafficking (FDFT1, APOL2, SCARB1), and Golgi/endoplasmic reticulum (ER) functions (SEMA3B, GCC1) were all increased. Quantitative PCR validated decreases in calmodulin 2 (CALM2) mRNA and increases in apolipoprotein L, 2 (APOL2) and semaphorin 3B (SEMA3B) mRNA for individual cases. A comparison between control cases with and without cardiovascular disease and elevated body mass index indicated that these changes were not due to general cellular and metabolic stress, but appeared specific to the use of drugs. Therefore, humans who abused cocaine, cannabis and/or phencyclidine share a decrease in transcription of calmodulin-related genes and increased transcription related to lipid/cholesterol and Golgi/ER function. These changes represent common molecular features of drug abuse, which may underlie changes in synaptic function and plasticity that could have important ramifications for decision-making capabilities in drug abusers.

Cytochrome P450 (P450) is a super-family of drug metabolizing enzymes. P450 enzymes have dual function; they can metabolize drugs to pharmacologically inactive metabolites facilitating their excretion or biotransform them to pharmacologically active metabolites which may have longer half-life than the parent drug. The variable pharmacological response to psychoactive drugs typically seen in population groups is often not accountable by considering dissimilarities in hepatic metabolism. Metabolism in brain specific nuclei may play a role in pharmacological modulation of drugs acting on the CNS and help explain some of the diverse response to these drugs seen in patient population. P450 enzymes are also present in brain where drug metabolism can take place and modify therapeutic action of drugs at the site of action. We have earlier demonstrated an intrinsic difference in the biotransformation of alprazolam (ALP) in brain and liver, relatively more alpha-hydroxy alprazolam (alpha-OHALP) is formed in brain as compared to liver. In the present study we show that recombinant CYP3A43 metabolizes ALP to both alpha-OHALP and 4-hydroxy alprazolam (4-OHALP) while CYP3A4 metabolizes ALP predominantly to its inactive metabolite, 4-OHALP. The expression of CYP3A43 mRNA in human brain samples correlates with formation of relatively higher levels of alpha-OH ALP indicating that individuals who express higher levels of CYP3A43 in the brain would generate larger amounts of alpha-OHALP. Further, the expression of CYP3A43 was relatively higher in brain as compared to liver across different ethnic populations. Since CYP3A enzymes play a prominent role in the metabolism of drugs, the higher expression of CYP3A43 would generate metabolite profile of drugs differentially in human brain and thus impact the pharmacodynamics of psychoactive drugs at the site of action.

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.

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.

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.

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.

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.

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 .

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), a crucial epigenetic mark, plays a key role in gene regulation, mammalian development, and various human diseases. Single-cell technologies enable the profiling of DNAm states at cytosines within the DNA sequence of individual cells, but they often suffer from limited coverage of CpG sites. In this study, we introduce scMeFormer, a transformer-based deep learning model designed to impute DNAm states for each CpG site in single cells. Through comprehensive evaluations, we demonstrate the superior performance of scMeFormer compared to alternative models across four single-nucleus DNAm datasets generated by distinct technologies. Remarkably, scMeFormer exhibits high-fidelity imputation, even when dealing with significantly reduced coverage, as low as 10% of the original CpG sites. Furthermore, we applied scMeFormer to a single-nucleus DNAm dataset generated from the prefrontal cortex of four schizophrenia patients and four neurotypical controls. This enabled the identification of thousands of differentially methylated regions associated with schizophrenia that would have remained undetectable without imputation and added granularity to our understanding of epigenetic alterations in schizophrenia within specific cell types. Our study highlights the power of deep learning in imputing DNAm states in single cells, and we expect scMeFormer to be a valuable tool for single-cell DNAm studies.

Cellular deconvolution of bulk RNA-sequencing (RNA-seq) data using single cell or nuclei RNA-seq (sc/snRNA-seq) reference data is an important strategy for estimating cell type composition in heterogeneous tissues, such as human brain. Several deconvolution methods have been developed and they have been previously benchmarked against simulated data, pseudobulked sc/snRNA-seq data, or cell type proportions derived from immunohistochemistry reference data. A major limitation preventing the improvement of deconvolution algorithms has been the lack of highly integrated datasets with orthogonal measurements of gene expression and estimates of cell type proportions on the same tissue block. The performance of existing deconvolution algorithms has not yet been explored across different RNA extraction methods (e.g. cytosolic, nuclear, or whole cell RNA), different library preparation types (e.g. mRNA enrichment vs. ribosomal RNA depletion), or with matched single cell reference datasets.

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.

The human prefrontal cortex (PFC), a mastermind of the brain, is one of the last brain regions to mature. To investigate the role of epigenetics in the development of PFC, we examined DNA methylation in ∼14,500 genes at ∼27,000 CpG loci focused on 5' promoter regions in 108 subjects range in age from fetal to elderly. DNA methylation in the PFC shows unique temporal patterns across life. The fastest changes occur during the prenatal period, slow down markedly after birth and continue to slow further with aging. At the genome level, the transition from fetal to postnatal life is typified by a reversal of direction, from demethylation prenatally to increased methylation postnatally. DNA methylation is strongly associated with genotypic variants and correlates with expression of a subset of genes, including genes involved in brain development and in de novo DNA methylation. Our results indicate that promoter DNA methylation in the human PFC is a highly dynamic process modified by genetic variance and regulating gene transcription. Additional discovery is made possible with a stand-alone application, BrainCloudMethyl.

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.