Effects of repeated alcohol abstinence on within-subject prefrontal cortical gene expression in rhesus macaques.

Male rhesus monkeys (n = 24) had a biopsy of prefrontal cortical area 46 prior to chronic ethanol self-administration (n = 17) or caloric control (n = 7). Fourteen months of daily self-administration (water vs. 4% alcohol, 22 h access/day termed "open-access") was followed by two cycles of prolonged abstinence (5 weeks) each followed by 3 months of open-access alcohol and a final abstinence followed by necropsy. At necropsy, a biopsy of Area 46, contralateral to the original biopsy, was obtained. Gene expression data (RNA-Seq) were collected comparing biopsy/necropsy samples. Monkeys were categorized by drinking status during the final post-abstinent drinking phase as light (LD), binge (BD), heavy (HD) and very heavy (VHD drinkers). Comparing pre-ethanol to post-abstinent biopsies, four animals that converted from HD to VHD status had significant ontology enrichments in downregulated genes (necropsy minus biopsy n = 286) that included immune response (FDR < 9 × 10-7) and plasma membrane changes (FDR < 1 × 10-7). Genes in the immune response category included IL16 and 18, CCR1, B2M, TLR3, 6 and 7, SP2 and CX3CR1. Upregulated genes (N = 388) were particularly enriched in genes associated with the negative regulation of MAP kinase activity (FDR < 3 × 10-5), including DUSP 1, 4, 5, 6 and 18, SPRY 2, 3, and 4, SPRED2, BMP4 and RGS2. Overall, these data illustrate the power of the NHP model and the within-subject design of genomic changes due to alcohol and suggest new targets for treating severe escalated drinking following repeated alcohol abstinence attempts.

Effect of chronic ethanol consumption in rhesus macaques on the nucleus accumbens core transcriptome.

The nucleus accumbens core (NAcc) has been repeatedly demonstrated to be a key component of the circuitry associated with excessive ethanol consumption. Previous studies have illustrated that in a nonhuman primate (NHP) model of chronic ethanol consumption, there is significant epigenetic remodeling of the NAcc. In the current study, RNA-Seq was used to examine genome-wide gene expression in eight each of control, low/binge (LD*), and high/very high (HD*) rhesus macaque drinkers. Using an FDR < 0.05, zero genes were significantly differentially expressed (DE) between LD* and controls, six genes between HD* and LD*, and 734 genes between HD* and controls. Focusing on HD* versus control DE genes, the upregulated genes (N = 366) were enriched in genes with annotations associated with signal recognition particle (SRP)-dependent co-translational protein targeting to membrane (FDR < 3 × 10-59 ), structural constituent of ribosome (FDR < 3 × 10-47 ), and ribosomal subunit (FDR < 5 × 10-48 ). Downregulated genes (N = 363) were enriched in annotations associated with behavior (FDR < 2 × 10-4 ), membrane organization (FDR < 1 × 10-4 ), inorganic cation transmembrane transporter activity (FDR < 2 × 10-3 ), synapse part (FDR < 4 × 10-10 ), glutamatergic synapse (FDR < 1 × 10-6 ), and GABAergic synapse (FDR < 6 × 10-4 ). Ingenuity Pathway Analysis (IPA) revealed that EIF2 signaling and mTOR pathways were significantly upregulated in HD* animals (FDR < 3 × 10-33 and <2 × 10-16 , respectively). Overall, the data supported our working hypothesis; excessive consumption would be associated with transcriptional differences in GABA/glutamate-related genes.

Identifying RNA Biomarkers and Molecular Pathways Involved in Multiple Subtypes of Uveitis.

PURPOSE: Uveitis is a heterogeneous collection of diseases. We tested the hypothesis that despite the diversity of uveitides, there could be common mechanisms shared by multiple subtypes, and that evidence of these common mechanisms may be detected as gene expression profiles in whole blood. DESIGN: Cohort study. METHODS: Ninety subjects with uveitis including axial spondyloarthritis (n = 17), sarcoidosis (n = 13), inflammatory bowel disease (n = 12), tubulointerstitial nephritis with uveitis (n = 10), or idiopathic uveitis (n = 38) as well as 18 healthy controls were enrolled, predominantly at Oregon Health & Science University. RNA-Seq data generated from peripheral, whole blood identified 19,859 unique transcripts. We analyzed gene expression pathways via Kyoto Encyclopedia of Genes and Genomes and Gene Ontology (GO). We validated our list of upregulated genes by comparison to a previously published study on peripheral blood gene expression among 50 subjects with diverse forms of uveitis. RESULTS: Both the Kyoto Encyclopedia of Genes and Genomes and GO analysis identified multiple shared pathways or GO terms with a P value of <.0001. Almost all pathways related to the immune response and/or response to an infection. A total of 119 individual transcripts were upregulated by at least 1.5-fold and false discovery rate <.05, and 61 were downregulated by similar criteria. Comparing mRNA from our study with a false discovery rate <.05 and the prior report, we identified 10 common gene transcripts: ICAM1, IL15RA, IL15, IRF1, IL10RB, GSK3A, TYK2, MEF2A, MEF2B, and MEF2D. CONCLUSIONS: Many forms of uveitis share overlapping mechanisms. These data support the concept that a single therapeutic approach could benefit diverse forms of this disease.

Revising the Diagnosis of Idiopathic Uveitis by Peripheral Blood Transcriptomics.

PURPOSE: To test the hypothesis that idiopathic uveitis can be categorized into subtypes based on gene expression from blood. DESIGN: Case control study. METHODS: We applied RNA-Seq to peripheral blood from patients with uveitis associated with 1 of 4 systemic diseases, including axial spondyloarthritis (n = 17), sarcoidosis (n = 13), inflammatory bowel disease (n = 12), tubulo-interstitial nephritis with uveitis (n = 10), or idiopathic uveitis (n = 38) as well as 18 healthy control subjects evaluated predominantly at Oregon Health and Science University. A high-dimensional negative binomial regression model implemented in the edgeR R package compared each disease group with the control subjects. The 20 most distinctive genes for each diagnosis were extracted. Of 80 genes, there were 75 unique genes. A classification algorithm was developed by fitting a gradient boosting tree with 5-fold cross-validation. Messenger RNA from subjects with idiopathic uveitis were analyzed to see if any fit clinically and by gene expression pattern with one of the diagnosable entities. RESULTS: For uveitis associated with a diagnosable systemic disease, gene expression profiling achieved an overall accuracy of 85% (balanced average of sensitivity plus specificity, P < .001). Although most patients with idiopathic uveitis presumably have none of these 4 associated systemic diseases, gene expression profiles helped to reclassify 11 of 38 subjects. CONCLUSIONS: Peripheral blood gene expression profiling is a potential adjunct in accurate differential diagnosis of the cause of uveitis. Validation of these results and characterization of the gene expression profile from additional discrete diagnoses could enhance the value of these observations.

Phenotypic and gene expression features associated with variation in chronic ethanol consumption in heterogeneous stock collaborative cross mice.

Of the more than 100 studies that have examined relationships between excessive ethanol consumption and the brain transcriptome, few rodent studies have examined chronic consumption. Heterogeneous stock collaborative cross mice freely consumed ethanol vs. water for 3 months. Transcriptional differences were examined for the central nucleus of the amygdala, a brain region known to impact ethanol preference. Early preference was modestly predictive of final preference and there was significant escalation of preference in females only. Genes significantly correlated with female preference were enriched in annotations for the primary cilium and extracellular matrix. A single module in the gene co-expression network was enriched in genes with an astrocyte annotation. The key hub node was the master regulator, orthodenticle homeobox 2 (Otx2). These data support an important role for the extracellular matrix, primary cilium and astrocytes in ethanol preference and consumption differences among individual female mice of a genetically diverse population.

Accelerating prediction of chemical shift of protein structures on GPUs: Using OpenACC.

Experimental chemical shifts (CS) from solution and solid state magic-angle-spinning nuclear magnetic resonance (NMR) spectra provide atomic level information for each amino acid within a protein or protein complex. However, structure determination of large complexes and assemblies based on NMR data alone remains challenging due to the complexity of the calculations. Here, we present a hardware accelerated strategy for the estimation of NMR chemical-shifts of large macromolecular complexes based on the previously published PPM_One software. The original code was not viable for computing large complexes, with our largest dataset taking approximately 14 hours to complete. Our results show that serial code refactoring and parallel acceleration brought down the time taken of the software running on an NVIDIA Volta 100 (V100) Graphic Processing Unit (GPU) to 46.71 seconds for our largest dataset of 11.3 million atoms. We use OpenACC, a directive-based programming model for porting the application to a heterogeneous system consisting of x86 processors and NVIDIA GPUs. Finally, we demonstrate the feasibility of our approach in systems of increasing complexity ranging from 100K to 11.3M atoms.

RNA-Seq of human whole blood: Evaluation of globin RNA depletion on Ribo-Zero library method.

Peripheral blood is a highly accessible biofluid providing a rich source of information about human physiology and health status. However, for studies of the blood transcriptome with RNA sequencing (RNA-Seq) techniques, high levels of hemoglobin mRNAs (hgbRNA) present in blood can occupy valuable sequencing space, impacting detection and quantification of non-hgbRNAs. In this study, we evaluated two methods for preparing ribosomal RNA (rRNA)-depleted sequencing libraries for RNA-Seq of whole blood, one of which is also designed to deplete hgbRNAs. Two experiments were performed: one evaluating library performance across 6 human blood samples and the other examining library reproducibility and performance in a two-subject subset. We find that addition of hgbRNA depletion to the rRNA-depletion protocol for library preparation from blood RNA effectively reduces highly abundant hgbRNA reads; however, it does not result in a statistically significant increase in differentially expressed genes in our patient-control study. Bioinformatic removal of globin gene counts in non-hgbRNA depleted libraries provides improvement in overall performance of these libraries. We conclude that use of a standard ribosomal RNA depletion method for library preparation coupled with bioinformatic removal of globin gene counts is sufficient for reproducible and sensitive measurement of both coding and noncoding RNAs in the blood transcriptome.

Chronic Voluntary Ethanol Drinking in Cynomolgus Macaques Elicits Gene Expression Changes in Prefrontal Cortical Area 46.

BACKGROUND: Genome-wide profiling to examine brain transcriptional features associated with excessive ethanol (EtOH) consumption has been applied to a variety of species including rodents, nonhuman primates (NHPs), and humans. However, these data were obtained from cross-sectional samples which are particularly vulnerable to individual variation when obtained from small outbred populations typical of human and NHP studies. In the current study, a novel within-subject design was used to examine the effects of voluntary EtOH consumption on prefrontal cortex (PFC) gene expression in a NHP model. METHODS: Two cohorts of cynomolgus macaques (n = 23) underwent a schedule-induced polydipsia procedure to establish EtOH self-administration followed by 6 months of daily open access to EtOH (4% w/v) and water. Individual daily EtOH intakes ranged from an average of 0.7 to 3.7 g/kg/d. Dorsal lateral PFC area 46 (A46) brain biopsies were collected in EtOH-naïve and control monkeys; contralateral A46 biopsies were collected from the same monkeys following the 6 months of fluid consumption. Gene expression changes were assessed using RNA-Seq paired analysis, which allowed for correction of individual baseline differences in gene expression. RESULTS: A total of 675 genes were significantly down-regulated following EtOH consumption; these were functionally enriched for immune response, cell adhesion, plasma membrane, and extracellular matrix. A total of 567 genes that were up-regulated following EtOH consumption were enriched in microRNA target sites and included target sites associated with Toll-like receptor pathways. The differentially expressed genes were also significantly enriched in transcription factor binding sites. CONCLUSIONS: The data presented here are the first to use a longitudinal biopsy strategy to examine how chronic EtOH consumption affects gene expression in the primate PFC. Prominent effects were seen in both cell adhesion and neuroimmune pathways; the latter contained both pro- and antiinflammatory genes. The data also indicate that changes in miRNAs and transcription factors may be important epigenetic regulators of EtOH consumption.

RNA-seq from archival FFPE breast cancer samples: molecular pathway fidelity and novel discovery.

BACKGROUND: Formalin-fixed, paraffin-embedded (FFPE) tissues for RNA-seq have advantages over fresh frozen tissue including abundance and availability, connection to rich clinical data, and association with patient outcomes. However, FFPE-derived RNA is highly degraded and chemically modified, which impacts its utility as a faithful source for biological inquiry. METHODS: True archival FFPE breast cancer cases (n = 58), stored at room temperature for 2-23 years, were utilized to identify key steps in tissue selection, RNA isolation, and library choice. Gene expression fidelity was evaluated by comparing FFPE data to public data obtained from fresh tissues, and by employing single-gene, gene set and transcription network-based regulon analyses. RESULTS: We report a single 10 µm section of breast tissue yields sufficient RNA for RNA-seq, and a relationship between RNA quality and block age that was not linear. We find single-gene analysis is limiting with FFPE tissues, while targeted gene set approaches effectively distinguish ER+ from ER- breast cancers. Novel utilization of regulon analysis identified the transcription factor KDM4B to associate with ER+ disease, with KDM4B regulon activity and gene expression having prognostic significance in an independent cohort of ER+ cases. CONCLUSION: Our results, which outline a robust FFPE-RNA-seq pipeline for broad use, support utilizing FFPE tissues to address key questions in the breast cancer field, including the delineation between indolent and life-threatening disease, biological stratification and molecular mechanisms of treatment resistance.

Genomic landscape of neutrophilic leukemias of ambiguous diagnosis.

Chronic neutrophilic leukemia (CNL), atypical chronic myeloid leukemia (aCML), and myelodysplastic/myeloproliferative neoplasms, unclassifiable (MDS/MPN-U) are a group of rare and heterogeneous myeloid disorders. There is strong morphologic resemblance among these distinct diagnostic entities as well as a lack of specific molecular markers and limited understanding of disease pathogenesis, which has made diagnosis challenging in certain cases. The treatment has remained empirical, resulting in dismal outcomes. We, therefore, performed whole-exome and RNA sequencing of these rare hematologic malignancies and present the most complete survey of the genomic landscape of these diseases to date. We observed a diversity of combinatorial mutational patterns that generally do not cluster within any one diagnosis. Gene expression analysis reveals enrichment, but not cosegregation, of clinical and genetic disease features with transcriptional clusters. In conclusion, these groups of diseases represent a continuum of related diseases rather than discrete diagnostic entities.

Clinical resistance to crenolanib in acute myeloid leukemia due to diverse molecular mechanisms.

FLT3 mutations are prevalent in AML patients and confer poor prognosis. Crenolanib, a potent type I pan-FLT3 inhibitor, is effective against both internal tandem duplications and resistance-conferring tyrosine kinase domain mutations. While crenolanib monotherapy has demonstrated clinical benefit in heavily pretreated relapsed/refractory AML patients, responses are transient and relapse eventually occurs. Here, to investigate the mechanisms of crenolanib resistance, we perform whole exome sequencing of AML patient samples before and after crenolanib treatment. Unlike other FLT3 inhibitors, crenolanib does not induce FLT3 secondary mutations, and mutations of the FLT3 gatekeeper residue are infrequent. Instead, mutations of NRAS and IDH2 arise, mostly as FLT3-independent subclones, while TET2 and IDH1 predominantly co-occur with FLT3-mutant clones and are enriched in crenolanib poor-responders. The remaining patients exhibit post-crenolanib expansion of mutations associated with epigenetic regulators, transcription factors, and cohesion factors, suggesting diverse genetic/epigenetic mechanisms of crenolanib resistance. Drug combinations in experimental models restore crenolanib sensitivity.

Regional Differences and Similarities in the Brain Transcriptome for Mice Selected for Ethanol Preference From HS-CC Founders.

The high genetic complexity found in heterogeneous stock (HS-CC) mice, together with selective breeding, can be used to detect new pathways and mechanisms associated with ethanol preference and excessive ethanol consumption. We predicted that these pathways would provide new targets for therapeutic manipulation. Previously (Colville et al., 2017), we observed that preference selection strongly affected the accumbens shell (SH) genes associated with synaptic function and in particular genes associated with synaptic tethering. Here we expand our analyses to include substantially larger sample sizes and samples from two additional components of the "addiction circuit," the central nucleus of the amygdala (CeA) and the prelimbic cortex (PL). At the level of differential expression (DE), the majority of affected genes are region-specific; only in the CeA did the DE genes show a significant enrichment in GO annotation categories, e.g., neuron part. In all three brain regions the differentially variable genes were significantly enriched in a single network module characterized by genes associated with cell-to-cell signaling. The data point to glutamate plasticity as being a key feature of selection for ethanol preference. In this context the expression of Dlg2 which encodes for PSD-93 appears to have a key role. It was also observed that the expression of the clustered protocadherins was strongly associated with preference selection.

Cross-species molecular dissection across alcohol behavioral domains.

This review summarizes the proceedings of a symposium presented at the "Alcoholism and Stress: A Framework for Future Treatment Strategies" conference held in Volterra, Italy on May 9-12, 2017. Psychiatric diseases, including alcohol-use disorders (AUDs), are influenced through complex interactions of genes, neurobiological pathways, and environmental influences. A better understanding of the common neurobiological mechanisms underlying an AUD necessitates an integrative approach, involving a systematic assessment of diverse species and phenotype measures. As part of the World Congress on Stress and Alcoholism, this symposium provided a detailed account of current strategies to identify mechanisms underlying the development and progression of AUDs. Dr. Sean Farris discussed the integration and organization of transcriptome and postmortem human brain data to identify brain regional- and cell type-specific differences related to excessive alcohol consumption that are conserved across species. Dr. Brien Riley presented the results of a genome-wide association study of DSM-IV alcohol dependence; although replication of genetic associations with alcohol phenotypes in humans remains challenging, model organism studies show that COL6A3, KLF12, and RYR3 affect behavioral responses to ethanol, and provide substantial evidence for their role in human alcohol-related traits. Dr. Rob Williams expanded upon the systematic characterization of extensive genetic-genomic resources for quantifying and clarifying phenotypes across species that are relevant to precision medicine in human disease. The symposium concluded with Dr. Robert Hitzemann's description of transcriptome studies in a mouse model selectively bred for high alcohol ("binge-like") consumption and a non-human primate model of long-term alcohol consumption. Together, the different components of this session provided an overview of systems-based approaches that are pioneering the experimental prioritization and validation of novel genes and gene networks linked with a range of behavioral phenotypes associated with stress and AUDs.

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.

On the relationships in rhesus macaques between chronic ethanol consumption and the brain transcriptome.

This is the first description of the relationship between chronic ethanol self-administration and the brain transcriptome in a non-human primate (rhesus macaque). Thirty-one male animals self-administered ethanol on a daily basis for over 12 months. Gene transcription was quantified with RNA-Seq in the central nucleus of the amygdala (CeA) and cortical Area 32. We constructed coexpression and cosplicing networks, and we identified areas of preservation and areas of differentiation between regions and network types. Correlations between intake and transcription included largely distinct gene sets and annotation categories across brain regions and between expression and splicing; positive and negative correlations were also associated with distinct annotation groups. Membrane, synaptic and splicing annotation categories were over-represented in the modules (gene clusters) enriched in positive correlations (CeA); our cosplicing analysis further identified the genes affected only at the exon inclusion level. In the CeA coexpression network, we identified Rab6b, Cdk18 and Igsf21 among the intake-correlated hubs, while in the Area 32, we identified a distinct hub set that included Ppp3r1 and Myeov2. Overall, the data illustrate that excessive ethanol self-administration is associated with broad expression and splicing mechanisms that involve membrane and synapse genes.

Anti-Müllerian hormone is produced heterogeneously in primate preantral follicles and is a potential biomarker for follicle growth and oocyte maturation in vitro.

PURPOSE: The main goals of this study were to investigate the expression of anti-Müllerian hormone (AMH) and its receptor (AMHR2) during follicular development in primates, and to evaluate the potential of AMH as a biomarker for follicle growth and oocyte maturation in vitro. METHODS: The mRNA and protein expression of AMH and AMHR2 were determined using isolated follicles and ovarian sections from rhesus macaques (n = 4) by real-time PCR and immunohistochemistry, respectively. Isolated secondary follicles were cultured individually. Follicle growth and media AMH concentrations were assessed by ELISA. The mRNA expression profiles, obtained from RNA sequencing, of in vitro- and in vivo-developed antral follicles were compared. Secondary follicles from additional animals (n = 35) were cultured. Follicle growth, oocyte maturation, and media AMH concentrations were evaluated for forecasting follicular development in vitro by AMH levels. RESULTS: AMH immunostaining was heterogeneous in the population of preantral follicles that were also stained for AMHR2. The mRNA expression profiles were comparable between in vivo- and in vitro-developed follicles. AMH levels produced by growing follicles were higher than those of nongrowing follicles in culture. With a cutoff value of 1.40 ng/ml, 85 % of nongrowing follicles could be identified while eliminating only 5 % of growing follicles. Growing follicles that generated metaphase II-stage oocytes secreted greater amounts of AMH than did those yielding immature germinal vesicle-stage oocytes. CONCLUSIONS: AMH, co-expressed with AMHR2, was produced heterogeneously by preantral follicles in macaques with levels correlated positively with follicle growth and oocyte maturation. AMH may serve as a biomarker for primate follicular development in vitro.

Splicing landscape of the eight collaborative cross founder strains.

BACKGROUND: The Collaborative Cross (CC) is a large panel of genetically diverse recombinant inbred mouse strains specifically designed to provide a systems genetics resource for the study of complex traits. In part, the utility of the CC stems from the extensive genome-wide annotations of founder strain sequence and structural variation. Still missing, however, are transcriptome-specific annotations of the CC founder strains that could further enhance the utility of this resource. RESULTS: We provide a comprehensive survey of the splicing landscape of the 8 CC founder strains by leveraging the high level of alternative splicing within the brain. Using deep transcriptome sequencing, we found that a majority of the splicing landscape is conserved among the 8 strains, with ~65% of junctions being shared by at least 2 strains. We, however, found a large number of potential strain-specific splicing events as well, with an average of ~3000 and ~500 with ≥3 and ≥10 sequence read coverage, respectively, within each strain. To better understand strain-specific splicing within the CC founder strains, we defined criteria for and identified high-confidence strain-specific splicing events. These splicing events were defined as exon-exon junctions 1) found within only one strain, 2) with a read coverage ≥10, and 3) defined by a canonical splice site. With these criteria, a total of 1509 high-confidence strain-specific splicing events were identified, with the majority found within two of the wild-derived strains, CAST and PWK. Strikingly, the overwhelming majority, 94%, of these strain-specific splicing events are not yet annotated. Strain-specific splicing was also located within genomic regions recently reported to be over- and under-represented within CC populations. CONCLUSIONS: Phenotypic characterization of CC populations is increasing; thus these results will not only aid in further elucidating the transcriptomic architecture of the individual CC founder strains, but they will also help in guiding the utilization of the CC populations in the study of complex traits. This report is also the first to establish guidelines in defining and identifying strain-specific splicing across different mouse strains.

Introduction to sequencing the brain transcriptome.

High-throughput next-generation sequencing is now entering its second decade. However, it was not until 2008 that the first report of sequencing the brain transcriptome appeared (Mortazavi, Williams, Mccue, Schaeffer, & Wold, 2008). These authors compared short-read RNA-Seq data for mouse whole brain with microarray results for the same sample and noted both the advantages and disadvantages of the RNA-Seq approach. While RNA-Seq provided exon level resolution, the majority of the reads were provided by a small proportion of highly expressed genes and the data analysis was exceedingly complex. Over the past 6 years, there have been substantial improvements in both RNA-Seq technology and data analysis. This volume contains 11 chapters that detail various aspects of sequencing the brain transcriptome. Some of the chapters are very methods driven, while others focus on the use of RNA-Seq to study such diverse areas as development, schizophrenia, and drug abuse. This chapter briefly reviews the transition from microarrays to RNA-Seq as the preferred method for analyzing the brain transcriptome. Compared with microarrays, RNA-Seq has a greater dynamic range, detects both coding and noncoding RNAs, is superior for gene network construction, detects alternative spliced transcripts, and can be used to extract genotype information, e.g., nonsynonymous coding single nucleotide polymorphisms. RNA-Seq embraces the complexity of the brain transcriptome and provides a mechanism to understand the underlying regulatory code; the potential to inform the brain-behavior-disease relationships is substantial.

The genetics of gene expression in complex mouse crosses as a tool to study the molecular underpinnings of behavior traits.

Complex Mus musculus crosses provide increased resolution to examine the relationships between gene expression and behavior. While the advantages are clear, there are numerous analytical and technological concerns that arise from the increased genetic complexity that must be considered. Each of these issues is discussed, providing an initial framework for complex cross study design and planning.

Utilizing RNA-Seq data for de novo coexpression network inference.

MOTIVATION: RNA-Seq experiments have shown great potential for transcriptome profiling. While sequencing increases the level of biological detail, integrative data analysis is also important. One avenue is the construction of coexpression networks. Because the capacity of RNA-Seq data for network construction has not been previously evaluated, we constructed a coexpression network using striatal samples, derived its network properties and compared it with microarray-based networks. RESULTS: The RNA-Seq coexpression network displayed scale-free, hierarchical network structure. We detected transcripts groups (modules) with correlated profiles; modules overlap distinct ontology categories. Neuroanatomical data from the Allen Brain Atlas reveal several modules with spatial colocalization. The network was compared with microarray-derived networks; correlations from RNA-Seq data were higher, likely because greater sensitivity and dynamic range. Higher correlations result in higher network connectivity, heterogeneity and centrality. For transcripts present across platforms, network structure appeared largely preserved. From this study, we present the first RNA-Seq data de novo network inference.