Introduction to sequencing the brain transcriptome

Transcriptomics : Approaches to Quantifying Gene Expression and Their Application to Studying the Human Brain

To date, the field of transcriptomics has been characterized by rapid methods development and technological advancement, with new technologies continuously rendering older ones obsolete.This chapter traces the evolution of approaches to quantifying gene expression and provides an overall view of the current state of the field of transcriptomics, its applications to the study of the human brain, and its place in the broader emerging multiomics landscape.

Decoding the Role of Astrocytes in the Entorhinal Cortex in Alzheimer’s Disease Using High-Dimensional Single-Nucleus RNA Sequencing Data and Next-Generation Knowledge Discovery Methodologies: Focus on Drugs and Natural Product Remedies for Dementia

Introduction: Alzheimer’s disease (AD) is a major cause of the development of cognitive decline and dementia. AD and associated dementias (ADRD) are the major contributors to the enormous burden of morbidity and mortality worldwide. To date, there are no robust therapies to alleviate or cure this debilitating disease. Most drug treatments focus on restoring the normal function of neurons and the cells that cause inflammation, such as microglia in the brain. However, the role of astrocytes, the brain’s housekeeping cells, in the development of AD and the initiation of dementia is still not well understood.

Objective: To decipher the role of astrocytes in the entorhinal cortex of AD patients using single nuclear RNA sequencing (snRNASeq) datasets from the Single Cell RNA-seq Database for Alzheimer’s Disease (scREAD). The datasets were originally derived from astrocytes, isolated from the entorhinal cortex of AD brain and healthy brain to decipher disease-specific signaling pathways as well as drugs and natural products that reverse AD-specific signatures in astrocytes.

Methods: We used snRNASeq datasets from the scREAD database originally derived from astrocytes isolated from the entorhinal cortex of AD and healthy brains from the Gene Expression Omnibus (GEO) (GSE138852 and GSE147528) and analyzed them using next-generation knowledge discovery (NGKD) platforms. scREAD is a user-friendly open-source interface available at https://bmbls.bmi.osumc.edu/scread/that enables more discovery-oriented strategies. snRNASeq data and metadata can also be visualized and downloaded via an interactive web application at adsn.ddnetbio.com. Differentially expressed genes (DEGs) for each snRNASeq dataset were analyzed using iPathwayGuide to compare and derive disease-specific pathways, gene ontologies, and in silico predictions of drugs and natural products that regulate AD -specific signatures in astrocytes. In addition, DEGs were analyzed using the L1000FWD and L1000CDS2 signature search programming interfaces (APIs) to identify additional drugs and natural products that mimic or reverse AD-specific gene signatures in astrocytes.

Results: We found that PI3K/AKT signaling, Wnt signaling, neuroactive ligand-receptor interaction pathways, neurodegeneration pathways, etc. were significantly impaired in astrocytes from the entorhinal cortex of AD patients. Biological processes such as glutamate receptor signaling pathway, regulation of synapse organization, cell-cell adhesion via plasma membrane adhesion molecules, and chylomicrons were negatively enriched in the astrocytes from the entorhinal cortex of AD patients. Gene sets involved in cellular components such as postsynaptic membrane, synaptic membrane, postsynapse, and synapse part were negatively enriched (p < 0.01). Moreover, molecular functions such as glutamate receptor activity, neurotransmitter receptor activity, and extracellular ligand-gated ion channels were negatively regulated in the astrocytes of the entorhinal cortex of AD patients (p < 0.01). Moreover, the application of NGKD platforms revealed that antirheumatic drugs, vitamin-E, emetine, narciclasine, cephaeline, trichostatin A, withaferin A, dasatinib, etc. can potentially reverse gene signatures associated with AD.

Conclusions: The present study highlights an innovative approach to use NGKD platforms to find unique disease-associated signaling pathways and specific synthetic drugs and natural products that can potentially reverse AD and ADRD-associated gene signatures.

Transcriptomics of Parental Care in the Hypothalamic-Septal Region of Female Zebra Finch Brain

(1) Background: The objective of this study was to uncover genomic causes of parental care. Since birds do not lactate and, therefore, do not show the gene expressional changes required for lactation, we investigate gene expression associated with parenting in caring and non-caring females in an avian species, the small passerine bird zebra finch (Taeniopygia guttata). Here, we compare expression patterns in the hypothalamic–septal region since, previously, we showed that this area is activated in parenting females. (2) Methods: Transcriptome sequencing was first applied in a dissected part of the zebra finch brain related to taking care of the nestlings as compared to a control group of social pairs without nestlings. (3) Results: We found genes differentially expressed between caring and non-caring females. When introducing a log2fold change threshold of 1.5, 13 annotated genes were significantly upregulated in breeding pairs, while 39 annotated genes were downregulated. Significant enrichments of dopamine and acetylcholine biosynthetic processes were identified among upregulated pathways, while pro-opiomelanocortin and thyroid hormone pathways were downregulated, suggesting the importance of these systems in parental care. Network analysis further suggested neuro-immunological changes in mothers. (4) Conclusions: The results confirm the roles of several hypothesized major pathways in parental care, whereas novel pathways are also proposed.

A modified approach for high-quality RNA extraction of spore-forming Bacillus subtilis at varied physiological stages

BACKGROUND

High quality RNA is required for the molecular study. Sample preparation of the spore-forming, Gram-positive bacteria like Bacillus sp., remains challenging although several methods have been proposed. Those techniques were simply developed using cell samples at certain growth stages despite some molecular studies like transcriptomic analyses require RNA samples from different physiological stages.

METHODS AND RESULTS

We developed the rapid, simple yet effective cell-lysis technique with limit use of harsh reagents by modifying the kit-based protocols. Appropriate lysozyme loading (20 mg/mL), incubation time (30 min), and temperature (37 °C) enabled cell lysis and enhanced RNA extraction from both vegetative cells and endospores of Bacillus subtilis TL7-3. High RNA Integrity Numbers and ratios of A₂₆₀/A₂₈₀ and A₂₆₀/A₂₃₀ of all RNA products collected during the batch cultivation confirmed that invert mixing with absolute ethanol prevented RNA damage during protein denaturation. With the process modification of the major steps in cell lysis and RNA extraction compared with the kit-based protocols that are typically used in laboratory work, interestingly, our modified protocol, simple-yet-effective, yielded higher concentration, purity, and integrity of RNA products from all cell samples collected at different physiological stages. While the kit-based protocols either failed to provide high RNA concentration or RNA purity and integrity for all cell samples particularly during the late-log, stationary, or sporulation.

CONCLUSIONS

Therefore, we can claim the significance of this modified protocol to be applicable for RNA extraction to those spore-forming Gram-positive bacteria not limited to B. subtilis growing at varied physiological stages.

RNA-Seq Analysis of Genetic and Transcriptome Network Effects of Dual-Trait Selection for Ethanol Preference and Withdrawal Using SOT and NOT Genetic Models

BACKGROUND

Genetic factors significantly affect alcohol consumption and vulnerability to withdrawal. Furthermore, some genetic models showing predisposition to severe withdrawal are also predisposed to low ethanol (EtOH) consumption and vice versa, even when tested independently in naïve animals.

METHODS

Beginning with a C57BL/6J × DBA/2J F2 intercross founder population, animals were simultaneously selectively bred for both high alcohol consumption and low acute withdrawal (SOT line), or vice versa (NOT line). Using randomly chosen fourth selected generation (S4) mice (N = 18-22/sex/line), RNA-Seq was employed to assess genome-wide gene expression in ventral striatum. The MegaMUGA array was used to detect genome-wide genotypic differences. Differential gene expression and the weighted gene co-expression network analysis were implemented as described elsewhere (Genes Brain Behav 16, 2017, 462).

RESULTS

The new selection of the SOT and NOT lines was similar to that reported previously (Alcohol Clin Exp Res 38, 2014, 2915). One thousand eight hundred and sixteen transcripts were detected as differentially expressed between the lines. For genes more highly expressed in the SOT line, there was enrichment in genes associated with cell adhesion, synapse organization, and postsynaptic membrane. The genes with a cell adhesion annotation included 23 protocadherins, Mpdz and Dlg2. Genes with a postsynaptic membrane annotation included Gabrb3, Gphn, Grid1, Grin2b, Grin2c, and Grm3. The genes more highly expressed in the NOT line were enriched in a network module (red) with annotations associated with mitochondrial function. Several of these genes were module hub nodes, and these included Nedd8, Guk1, Elof1, Ndufa8, and Atp6v1f.

CONCLUSIONS

Marked effects of selection on gene expression were detected. The NOT line was characterized by higher expression of hub nodes associated with mitochondrial function. Genes more highly expressed in the SOT aligned with previous findings, for example, Colville and colleagues (Genes Brain Behav 16, 2017, 462) that both high EtOH preference and consumption are associated with effects on cell adhesion and glutamate synaptic plasticity.

Regional Analysis of the Brain Transcriptome in Mice Bred for High and Low Methamphetamine Consumption

Transcriptome profiling can broadly characterize drug effects and risk for addiction in the absence of drug exposure. Modern large-scale molecular methods, including RNA-sequencing (RNA-Seq), have been extensively applied to alcohol-related disease traits, but rarely to risk for methamphetamine (MA) addiction. We used RNA-Seq data from selectively bred mice with high or low risk for voluntary MA intake to construct coexpression and cosplicing networks for differential risk. Three brain reward circuitry regions were explored, the nucleus accumbens (NAc), prefrontal cortex (PFC), and ventral midbrain (VMB). With respect to differential gene expression and wiring, the VMB was more strongly affected than either the PFC or NAc. Coexpression network connectivity was higher in the low MA drinking line than in the high MA drinking line in the VMB, oppositely affected in the NAc, and little impacted in the PFC. Gene modules protected from the effects of selection may help to eliminate certain mechanisms from significant involvement in risk for MA intake. One such module was enriched in genes with dopamine-associated annotations. Overall, the data suggest that mitochondrial function and glutamate-mediated synaptic plasticity have key roles in the outcomes of selective breeding for high versus low levels of MA intake.

Gender-Specific Effects of Selection for Drinking in the Dark on the Network Roles of Coding and Noncoding RNAs

BACKGROUND

Transcriptional differences between heterogeneous stock mice and high drinking-in-the-dark selected mouse lines have previously been described based on microarray technology coupled with network-based analysis. The network changes were reproducible in 2 independent selections and largely confined to 2 distinct network modules; in contrast, differential expression appeared more specific to each selected line. This study extends these results by utilizing RNA-Seq technology, allowing evaluation of the relationship between genetic risk and transcription of noncoding RNA (ncRNA); we additionally evaluate sex-specific transcriptional effects of selection.

METHODS

Naïve mice (N = 24/group and sex) were utilized for gene expression analysis in the ventral striatum; the transcriptome was sequenced with the Illumina HiSeq platform. Differential gene expression and the weighted gene co-expression network analysis were implemented largely as described elsewhere, resulting in the identification of genes that change expression level or (co)variance structure.

RESULTS

Across both sexes, we detect selection effects on the extracellular matrix and synaptic signaling, although the identity of individual genes varies. A majority of nc RNAs cluster in a single module of relatively low density in both the male and female network. The most strongly differentially expressed transcript in both sexes was Gm22513, a small nuclear RNA with unknown function. Associated with selection, we also found a number of network hubs that change edge strength and connectivity. At the individual gene level, there are many sex-specific effects; however, at the annotation level, results are more concordant.

CONCLUSIONS

In addition to demonstrating sex-specific effects of selection on the transcriptome, the data point to the involvement of extracellular matrix genes as being associated with the binge drinking phenotype.

Brain Transcriptome Databases: A User's Guide

Transcriptional programs instruct the generation and maintenance of diverse subtypes of neural cells, establishment of distinct brain regions, formation and function of neural circuits, and ultimately behavior. Spatiotemporal and cell type-specific analyses of the transcriptome, the sum total of all RNA transcripts in a cell or an organ, can provide insights into the role of genes in brain development and function, and their potential contribution to disorders of the brain. In the previous decade, advances in sequencing technology and funding from the National Institutes of Health and private foundations for large-scale genomics projects have led to a growing collection of brain transcriptome databases. These valuable resources provide rich and high-quality datasets with spatiotemporal, cell type-specific, and single-cell precision. Most importantly, many of these databases are publicly available via user-friendly web interface, making the information accessible to individual scientists without the need for advanced computational expertise. Here, we highlight key publicly available brain transcriptome databases, summarize the tissue sources and methods used to generate the data, and discuss their utility for neuroscience research.

Considerations for optimal use of postmortem human brains for molecular psychiatry: lessons from schizophrenia

Schizophrenia is a disabling disease impacting millions of people around the world, for which there is no known cure. Current antipsychotic treatments for schizophrenia mainly target psychotic symptoms, do little to ameliorate social or cognitive deficits, have side-effects that cause weight gain, and diabetes and 30% of people do not respond. Thus, better therapeutics for schizophrenia aimed at the route biologic changes are needed and discovering the underlying neurobiology is key to this quest. Postmortem brain studies provide the most direct and detailed way to determine the pathophysiology of schizophrenia. This chapter outlines steps that can be taken to ensure the best-quality molecular data from postmortem brain tissue are obtained. In this chapter, we also discuss targeted and high-throughput methods for examining gene and protein expression and some of the strengths and limitations of each method. We briefly consider why gene and protein expression changes may not always concur within brain tissue. We conclude that postmortem brain research that investigates gene and protein expression in well-characterized and matched brain cohorts provides an important foundation to be considered when interpreting data obtained from studies of living schizophrenia patients.

Coding and small non-coding transcriptional landscape of tuberous sclerosis complex cortical tubers: implications for pathophysiology and treatment

Tuberous Sclerosis Complex (TSC) is a rare genetic disorder that results from a mutation in the TSC1 or TSC2 genes leading to constitutive activation of the mechanistic target of rapamycin complex 1 (mTORC1). TSC is associated with autism, intellectual disability and severe epilepsy. Cortical tubers are believed to represent the neuropathological substrates of these disabling manifestations in TSC. In the presented study we used high-throughput RNA sequencing in combination with systems-based computational approaches to investigate the complexity of the TSC molecular network. Overall we detected 438 differentially expressed genes and 991 differentially expressed small non-coding RNAs in cortical tubers compared to autopsy control brain tissue. We observed increased expression of genes associated with inflammatory, innate and adaptive immune responses. In contrast, we observed a down-regulation of genes associated with neurogenesis and glutamate receptor signaling. MicroRNAs represented the largest class of over-expressed small non-coding RNA species in tubers. In particular, our analysis revealed that the miR-34 family (including miR-34a, miR-34b and miR-34c) was significantly over-expressed. Functional studies demonstrated the ability of miR-34b to modulate neurite outgrowth in mouse primary hippocampal neuronal cultures. This study provides new insights into the TSC transcriptomic network along with the identification of potential new treatment targets.

Gene expression elucidates functional impact of polygenic risk for schizophrenia

Over 100 genetic loci harbor schizophrenia-associated variants, yet how these variants confer liability is uncertain. The CommonMind Consortium sequenced RNA from dorsolateral prefrontal cortex of people with schizophrenia (N = 258) and control subjects (N = 279), creating a resource of gene expression and its genetic regulation. Using this resource, ∼20% of schizophrenia loci have variants that could contribute to altered gene expression and liability. In five loci, only a single gene was involved: FURIN, TSNARE1, CNTN4, CLCN3 or SNAP91. Altering expression of FURIN, TSNARE1 or CNTN4 changed neurodevelopment in zebrafish; knockdown of FURIN in human neural progenitor cells yielded abnormal migration. Of 693 genes showing significant case-versus-control differential expression, their fold changes were ≤ 1.33, and an independent cohort yielded similar results. Gene co-expression implicates a network relevant for schizophrenia. Our findings show that schizophrenia is polygenic and highlight the utility of this resource for mechanistic interpretations of genetic liability for brain diseases.

A RNA-Seq Analysis of the Rat Supraoptic Nucleus Transcriptome: Effects of Salt Loading on Gene Expression

Magnocellular neurons (MCNs) in the hypothalamo-neurohypophysial system (HNS) are highly specialized to release large amounts of arginine vasopressin (Avp) or oxytocin (Oxt) into the blood stream and play critical roles in the regulation of body fluid homeostasis. The MCNs are osmosensory neurons and are excited by exposure to hypertonic solutions and inhibited by hypotonic solutions. The MCNs respond to systemic hypertonic and hypotonic stimulation with large changes in the expression of their Avp and Oxt genes, and microarray studies have shown that these osmotic perturbations also cause large changes in global gene expression in the HNS. In this paper, we examine gene expression in the rat supraoptic nucleus (SON) under normosmotic and chronic salt-loading SL) conditions by the first time using "new-generation", RNA sequencing (RNA-Seq) methods. We reliably detect 9,709 genes as present in the SON by RNA-Seq, and 552 of these genes were changed in expression as a result of chronic SL. These genes reflect diverse functions, and 42 of these are involved in either transcriptional or translational processes. In addition, we compare the SON transcriptomes resolved by RNA-Seq methods with the SON transcriptomes determined by Affymetrix microarray methods in rats under the same osmotic conditions, and find that there are 6,466 genes present in the SON that are represented in both data sets, although 1,040 of the expressed genes were found only in the microarray data, and 2,762 of the expressed genes are selectively found in the RNA-Seq data and not the microarray data. These data provide the research community a comprehensive view of the transcriptome in the SON under normosmotic conditions and the changes in specific gene expression evoked by salt loading.

RNA Sequencing in Schizophrenia

Schizophrenia (SCZ) is a serious psychiatric disorder that affects 1% of general population and places a heavy burden worldwide. The underlying genetic mechanism of SCZ remains unknown, but studies indicate that the disease is associated with a global gene expression disturbance across many genes. Next-generation sequencing, particularly of RNA sequencing (RNA-Seq), provides a powerful genome-scale technology to investigate the pathological processes of SCZ. RNA-Seq has been used to analyze the gene expressions and identify the novel splice isoforms and rare transcripts associated with SCZ. This paper provides an overview on the genetics of SCZ, the advantages of RNA-Seq for transcriptome analysis, the accomplishments of RNA-Seq in SCZ cohorts, and the applications of induced pluripotent stem cells and RNA-Seq in SCZ research.