Transcriptomic analysis of the human habenula in schizophrenia.

Importance: Habenula (Hb) pathophysiology is involved in many neuropsychiatric disorders, including schizophrenia. Deep brain stimulation and pharmacological targeting of the Hb are emerging as promising therapeutic treatments. However, little is known about the cell type-specific transcriptomic organization of the human Hb or how it is altered in schizophrenia. Objective: To define the molecular neuroanatomy of the human habenula and identify transcriptomic changes in individuals with schizophrenia compared to neurotypical controls. Design Setting and Participants: This study utilized Hb-enriched postmortem human brain tissue. Single nucleus RNA-sequencing (snRNA-seq) and single molecule fluorescent in situ hybridization (smFISH) experiments were conducted to identify molecularly defined Hb cell types and map their spatial location (n=3-7 donors). Bulk RNA-sequencing and cell type deconvolution were used to investigate transcriptomic changes in Hb-enriched tissue from 35 individuals with schizophrenia and 33 neurotypical controls. Gene expression changes associated with schizophrenia in the Hb were compared to those previously identified in the dorsolateral prefrontal cortex (DLPFC), hippocampus, and caudate. Main Outcomes and Measures: Semi-supervised snRNA-seq cell type clustering. Transcript visualization and quantification of smFISH probes. Bulk RNA-seq cell type deconvolution using reference snRNA-seq data. Schizophrenia-associated gene differential expression analysis adjusting for Hb and thalamus fractions, RNA degradation-associated quality surrogate variables, and other covariates. Cross-brain region schizophrenia-associated gene expression comparison. Results: snRNA-seq identified 17 cell type clusters across 16,437 nuclei, including 3 medial and 7 lateral Hb populations. Cell types were conserved with those identified in a rodent model. smFISH for cell type marker genes validated snRNA-seq Hb cell types and depicted the spatial organization of subpopulations. Bulk RNA-seq analyses yielded 45 schizophrenia-associated differentially expressed genes (FDR < 0.05), with 32 (71%) unique to Hb-enriched tissue. Conclusions: These results identify topographically organized cell types with distinct molecular signatures in the human Hb. They further demonstrate unique transcriptomic changes in the epithalamus associated with schizophrenia, thereby providing molecular insights into the role of Hb in neuropsychiatric disorders.

Dissecting transcriptomic signatures of neuronal differentiation and maturation using iPSCs.

Human induced pluripotent stem cells (hiPSCs) are a powerful model of neural differentiation and maturation. We present a hiPSC transcriptomics resource on corticogenesis from 5 iPSC donor and 13 subclonal lines across 9 time points over 5 broad conditions: self-renewal, early neuronal differentiation, neural precursor cells (NPCs), assembled rosettes, and differentiated neuronal cells. We identify widespread changes in the expression of both individual features and global patterns of transcription. We next demonstrate that co-culturing human NPCs with rodent astrocytes results in mutually synergistic maturation, and that cell type-specific expression data can be extracted using only sequencing read alignments without cell sorting. We lastly adapt a previously generated RNA deconvolution approach to single-cell expression data to estimate the relative neuronal maturity of iPSC-derived neuronal cultures and human brain tissue. Using many public datasets, we demonstrate neuronal cultures are maturationally heterogeneous but contain subsets of neurons more mature than previously observed.

Characterizing the nuclear and cytoplasmic transcriptomes in developing and mature human cortex uncovers new insight into psychiatric disease gene regulation.

Transcriptome compartmentalization by the nuclear membrane provides both stochastic and functional buffering of transcript activity in the cytoplasm, and has recently been implicated in neurodegenerative disease processes. Although many mechanisms regulating transcript compartmentalization are also prevalent in brain development, the extent to which subcellular localization differs as the brain matures has yet to be addressed. To characterize the nuclear and cytoplasmic transcriptomes during brain development, we sequenced both RNA fractions from homogenate prenatal and adult human postmortem cortex using poly(A)+ and Ribo-Zero library preparation methods. We find that while many genes are differentially expressed by fraction and developmental expression changes are similarly detectable in nuclear and cytoplasmic RNA, the compartmented transcriptomes become more distinct as the brain matures, perhaps reflecting increased utilization of nuclear retention as a regulatory strategy in adult brain. We examined potential mechanisms of this developmental divergence including alternative splicing, RNA editing, nuclear pore composition, RNA-binding protein motif enrichment, and RNA secondary structure. Intron retention is associated with greater nuclear abundance in a subset of transcripts, as is enrichment for several splicing factor binding motifs. Finally, we examined disease association with fraction-regulated gene sets and found nuclear-enriched genes were also preferentially enriched in gene sets associated with neurodevelopmental psychiatric disorders. These results suggest that although gene-level expression is globally comparable between fractions, nuclear retention of transcripts may play an underappreciated role in developmental regulation of gene expression in brain, particularly in genes whose dysregulation is related to neuropsychiatric disorders.

Comprehensive assessment of multiple biases in small RNA sequencing reveals significant differences in the performance of widely used methods.

BACKGROUND: RNA sequencing offers advantages over other quantification methods for microRNA (miRNA), yet numerous biases make reliable quantification challenging. Previous evaluations of these biases have focused on adapter ligation bias with limited evaluation of reverse transcription bias or amplification bias. Furthermore, evaluations of the quantification of isomiRs (miRNA isoforms) or the influence of starting amount on performance have been very limited. No study had yet evaluated the quantification of isomiRs of altered length or compared the consistency of results derived from multiple moderate starting inputs. We therefore evaluated quantifications of miRNA and isomiRs using four library preparation kits, with various starting amounts, as well as quantifications following removal of duplicate reads using unique molecular identifiers (UMIs) to mitigate reverse transcription and amplification biases. RESULTS: All methods resulted in false isomiR detection; however, the adapter-free method tested was especially prone to false isomiR detection. We demonstrate that using UMIs improves accuracy and we provide a guide for input amounts to improve consistency. CONCLUSIONS: Our data show differences and limitations of current methods, thus raising concerns about the validity of quantification of miRNA and isomiRs across studies. We advocate for the use of UMIs to improve accuracy and reliability of miRNA quantifications.

Regional Heterogeneity in Gene Expression, Regulation, and Coherence in the Frontal Cortex and Hippocampus across Development and Schizophrenia.

The hippocampus formation, although prominently implicated in schizophrenia pathogenesis, has been overlooked in large-scale genomics efforts in the schizophrenic brain. We performed RNA-seq in hippocampi and dorsolateral prefrontal cortices (DLPFCs) from 551 individuals (286 with schizophrenia). We identified substantial regional differences in gene expression and found widespread developmental differences that were independent of cellular composition. We identified 48 and 245 differentially expressed genes (DEGs) associated with schizophrenia within the hippocampus and DLPFC, with little overlap between the brain regions. 124 of 163 (76.6%) of schizophrenia GWAS risk loci contained eQTLs in any region. Transcriptome-wide association studies in each region identified many novel schizophrenia risk features that were brain region-specific. Last, we identified potential molecular correlates of in vivo evidence of altered prefrontal-hippocampal functional coherence in schizophrenia. These results underscore the complexity and regional heterogeneity of the transcriptional correlates of schizophrenia and offer new insights into potentially causative biology.

BDNF-TrkB signaling in oxytocin neurons contributes to maternal behavior.

Brain-derived neurotrophic factor (Bdnf) transcription is controlled by several promoters, which drive expression of multiple transcripts encoding an identical protein. We previously reported that BDNF derived from promoters I and II is highly expressed in hypothalamus and is critical for regulating aggression in male mice. Here we report that BDNF loss from these promoters causes reduced sexual receptivity and impaired maternal care in female mice, which is concomitant with decreased oxytocin (Oxt) expression during development. We identify a novel link between BDNF signaling, oxytocin, and maternal behavior by demonstrating that ablation of TrkB selectively in OXT neurons partially recapitulates maternal care impairments observed in BDNF-deficient females. Using translating ribosome affinity purification and RNA-sequencing we define a molecular profile for OXT neurons and delineate how BDNF signaling impacts gene pathways critical for structural and functional plasticity. Our findings highlight BDNF as a modulator of sexually-dimorphic hypothalamic circuits that govern female-typical behaviors.

xMD-miRNA-seq to generate near in vivo miRNA expression estimates in colon epithelial cells.

Accurate, RNA-seq based, microRNA (miRNA) expression estimates from primary cells have recently been described. However, this in vitro data is mainly obtained from cell culture, which is known to alter cell maturity/differentiation status, significantly changing miRNA levels. What is needed is a robust method to obtain in vivo miRNA expression values directly from cells. We introduce expression microdissection miRNA small RNA sequencing (xMD-miRNA-seq), a method to isolate cells directly from formalin fixed paraffin-embedded (FFPE) tissues. xMD-miRNA-seq is a low-cost, high-throughput, immunohistochemistry-based method to capture any cell type of interest. As a proof-of-concept, we isolated colon epithelial cells from two specimens and performed low-input small RNA-seq. We generated up to 600,000 miRNA reads from the samples. Isolated epithelial cells, had abundant epithelial-enriched miRNA expression (miR-192; miR-194; miR-200b; miR-200c; miR-215; miR-375) and overall similar miRNA expression patterns to other epithelial cell populations (colonic enteroids and flow-isolated colon epithelium). xMD-derived epithelial cells were generally not contaminated by other adjacent cells of the colon as noted by t-SNE analysis. xMD-miRNA-seq allows for simple, economical, and efficient identification of cell-specific miRNA expression estimates. Further development will enhance rapid identification of cell-specific miRNA expression estimates in health and disease for nearly any cell type using archival FFPE material.

Dynamic regulation of RNA editing in human brain development and disease.

RNA editing is increasingly recognized as a molecular mechanism regulating RNA activity and recoding proteins. Here we surveyed the global landscape of RNA editing in human brain tissues and identified three unique patterns of A-to-I RNA editing rates during cortical development: stable high, stable low and increasing. RNA secondary structure and the temporal expression of adenosine deaminase acting on RNA (ADAR) contribute to cis- and trans-regulatory mechanisms of these RNA editing patterns, respectively. Interestingly, the increasing pattern was associated with neuronal maturation, correlated with mRNA abundance and potentially influenced miRNA binding energy. Gene ontology analyses implicated the increasing pattern in vesicle or organelle membrane-related genes and glutamate signaling pathways. We also found that the increasing pattern was selectively perturbed in spinal cord injury and glioblastoma. Our findings reveal global and dynamic aspects of RNA editing in brain, providing new insight into epitranscriptional regulation of sequence diversity.

Engineering targeted chromosomal amplifications in human breast epithelial cells.

Chromosomal amplifications are among the most common genetic alterations found in human cancers. However, experimental systems to study the processes that lead to specific, recurrent amplification events in human cancers are lacking. Moreover, some common amplifications, such as that at 8p11-12 in breast cancer, harbor multiple driver oncogenes, which are poorly modeled by conventional overexpression approaches. We sought to develop an experimental system to model recurrent chromosomal amplification events in human cell lines. Our strategy is to use homologous-recombination-mediated gene targeting to deliver a dominantly selectable, amplifiable marker to a specified chromosomal location. We used adeno-associated virus vectors to target human MCF-7 breast cancer cells at the ZNF703 locus, in the recurrent 8p11-12 amplicon, using the E. coli inosine monophosphate dehydrogenase (IMPDH) enzyme as a marker. We applied selective pressure using IMPDH inhibitors. Surviving clones were found to have increased copy number of ZNF703 (average 2.5-fold increase) by droplet digital PCR and FISH. Genome-wide array comparative genomic hybridization confirmed that amplifications had occurred on the short arm of chromosome 8, without changes on 8q or other chromosomes. Patterns of amplification were variable and similar to those seen in primary human breast cancers, including "sawtooth" patterns, distal copy number loss, and large continuous regions of copy number gain. This system will allow study of the cis- and trans-acting factors that are permissive for chromosomal amplification and provide a model to analyze oncogene cooperativity in amplifications harboring multiple candidate driver genes.

PIK3CA and AKT1 mutations have distinct effects on sensitivity to targeted pathway inhibitors in an isogenic luminal breast cancer model system.

PURPOSE: Activating mutations in the phosphoinositide-3-kinase (PI3K)/AKT/mTOR pathway are present in the majority of breast cancers and therefore are a major focus of drug development and clinical trials. Pathway mutations have been proposed as predictive biomarkers for efficacy of PI3K-targeted therapies. However, the precise contribution of distinct PI3K pathway mutations to drug sensitivity is unknown. EXPERIMENTAL DESIGN: We describe the creation of a physiologic human luminal breast cancer cell line model to study the phenotype of these mutations using the MCF-7 cell line. We used somatic cell gene targeting to "correct" PIK3CA E545K-mutant alleles in MCF-7 cells to wild-type sequence. The AKT1 E17K hotspot mutation was knocked in on this wild-type background. RESULTS: Loss of mutant PIK3CA dramatically reduced phosphorylation of AKT proteins and several known AKT targets, but other AKT target proteins and downstream effectors of mTOR were not affected. PIK3CA wild-type cells exhibited reduced proliferation in vitro and in vivo. Knockin of the AKT1 E17K hotspot mutation on this PIK3CA wild-type background restored pathway signaling, proliferation, and tumor growth in vivo. PIK3CA, but not AKT1 mutation, increased sensitivity to the PI3K inhibitor GDC-0941 and the allosteric AKT inhibitor MK-2206. CONCLUSIONS: AKT1 E17K is a bona fide oncogene in a human luminal breast cancer context. Distinct PI3K pathway mutations confer differential sensitivity to drugs targeting the pathway at different points and by distinct mechanisms. These findings have implications for the use of tumor genome sequencing to assign patients to targeted therapies.

Functional analysis of non-hotspot AKT1 mutants found in human breast cancers identifies novel driver mutations: implications for personalized medicine.

The phosphatidylinositol 3-kinase (PI3-kinase)-Akt-mTOR pathway is mutated at high frequency in human breast cancer, and this pathway is the focus of active drug discovery and clinical investigation. Trials of personalized cancer therapy seek to leverage knowledge of cancer gene mutations by using mutations to guide the choice of targeted therapies. At the same time, cancer genome sequencing studies are identifying low frequency variants of unknown significance in known cancer genes, as well as genes of unknown function. We have performed functional analysis of six non-hotspot AKT1 pleckstrin homology domain mutants identified in recent large-scale breast cancer sequencing studies. Three of these mutants cause constitutive activation of Akt1 in the absence of growth factors, leading to phosphorylation of downstream target proteins. Like the hotspot E17K mutation, these mutants confer constitutive membrane localization of Akt1. Finally, the same three mutants showed oncogenic activity in a cellular transformation assay. The other three mutants were inactive in all assays. These findings validate novel driver mutations in AKT1, and extend the number and type of mutations that activate the PI3-kinase pathway in human breast cancers.

Progression from mitotic catastrophe to germ cell death in Caenorhabditis elegans lis-1 mutants requires the spindle checkpoint.

Deletion of the lissencephaly disease gene LIS-1 in humans causes an extreme disorganization of the brain associated with significant reduction in cortical neurons. Here we show that deletion or RNA interference (RNAi) of Caenorhabditis elegans lis-1 results in a reduction in germline nuclei and causes a variety of cellular, developmental, and neurological defects throughout development. Our analysis of the germline defects suggests that the reduction in nuclei number stems from dysfunctional mitotic spindles resulting in cell cycle arrest and eventually programmed cell death (apoptosis). Deletion of the spindle checkpoint gene mdf-1 blocks lis-1(lf)-induced cell cycle arrest and germline apoptosis, placing the spindle checkpoint pathway upstream of the programmed cell death pathway. These results suggest that apoptosis may contribute to the cell-sparse pathology of lissencephaly.