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.

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.

Dual-trait selection for ethanol consumption and withdrawal: genetic and transcriptional network effects.

BACKGROUND: Data from C57BL/6J (B6) × DBA/2J (D2) F2 intercrosses (B6xD2 F2 ), standard and recombinant inbred strains, and heterogeneous stock mice indicate that a reciprocal (or inverse) genetic relationship exists between alcohol consumption and withdrawal severity. Furthermore, some genetic studies have detected reciprocal quantitative trait loci (QTLs) for these traits. We used a novel mouse model developed by simultaneous selection for both high alcohol consumption/low withdrawal and low alcohol consumption/high withdrawal and analyzed the gene expression and genome-wide genotypic differences. METHODS: Randomly chosen third selected generation (S3 ) mice (N = 24/sex/line), bred from a B6xD2 F2 , were genotyped using the Mouse Universal Genotyping Array, which provided 2,760 informative markers. QTL analysis used a marker-by-marker strategy with the threshold for a significant log of the odds (LOD) set at 10. Gene expression in the ventral striatum was measured using the Illumina Mouse 8.2 array. Differential gene expression and the weighted gene co-expression network analysis (WGCNA) were implemented. RESULTS: Significant QTLs for consumption/withdrawal were detected on chromosomes (Chr) 2, 4, 9, and 12. A suggestive QTL mapped to Chr 6. Some of the QTLs overlapped with known QTLs mapped for 1 of the traits individually. One thousand seven hundred and forty-five transcripts were detected as being differentially expressed between the lines; there was some overlap with known withdrawal genes (e.g., Mpdz) located within QTL regions. WGCNA revealed several modules of co-expressed genes showing significant effects in both differential expression and intramodular connectivity; a module richly annotated with kinase-related annotations was most affected. CONCLUSIONS: Marked effects of selection on expression and network structure were detected. QTLs overlapping with differentially expressed genes on Chr 2 (distal) and 4 suggest that these are cis-eQTLs (Chr 2: Kif3b, Kcnq2; Chr 4: Mpdz, Snapc3). Other QTLs identified were on Chr 2 (proximal), 9, and 12. Network results point to involvement of kinase-related mechanisms and outline the need for further efforts such as interrogation of noncoding RNAs.

Selection for drinking in the dark alters brain gene coexpression networks.

BACKGROUND: Heterogeneous stock (HS/NPT) mice have been used to create lines selectively bred in replicate for elevated drinking in the dark (DID). Both selected lines routinely reach a blood ethanol (EtOH) concentration (BEC) of 1.00 mg/ml or greater at the end of the 4-hour period of access in Day 2. The mechanisms through which genetic differences influence DID are currently unclear. Therefore, the current study examines the transcriptome, the first stage at which genetic variability affects neurobiology. Rather than focusing solely on differential expression (DE), we also examine changes in the ways that gene transcripts collectively interact with each other, as revealed by changes in coexpression patterns. METHODS: Naïve mice (N = 48/group) were genotyped using the Mouse Universal Genotyping Array, which provided 3,683 informative markers. Quantitative trait locus (QTL) analysis used a marker-by-marker strategy with the threshold for a significant logarithm of odds (LOD) set at 10.6. Gene expression in the ventral striatum was measured using the Illumina Mouse 8.2 array. Differential gene expression and the weighted gene coexpression network analysis (WGCNA) were implemented largely as described elsewhere. RESULTS: Significant QTLs for elevated BECs after DID were detected on chromosomes 4, 14, and 16; the latter 2 were associated with gene-poor regions. None of the QTLs overlapped with known QTLs for EtOH preference drinking. Ninety-four transcripts were detected as being differentially expressed in both selected lines versus HS controls; there was no overlap with known preference genes. The WGCNA revealed 2 modules as showing significant effects of both selections on intramodular connectivity. A number of genes known to be associated with EtOH phenotypes (e.g., Gabrg1, Glra2, Grik1, Npy2r, and Nts) showed significant changes in connectivity. CONCLUSIONS: We found marked and consistent effects of selection on coexpression patterns; DE changes were more modest and less concordant. The QTLs and differentially expressed genes detected here are distinct from the preference phenotype. This is consistent with behavioral data and suggests that the DID and preference phenotypes are markedly different genetically.

Genetic diversity and striatal gene networks: focus on the heterogeneous stock-collaborative cross (HS-CC) mouse.

BACKGROUND: The current study focused on the extent genetic diversity within a species (Mus musculus) affects gene co-expression network structure. To examine this issue, we have created a new mouse resource, a heterogeneous stock (HS) formed from the same eight inbred strains that have been used to create the collaborative cross (CC). The eight inbred strains capture > 90% of the genetic diversity available within the species. For contrast with the HS-CC, a C57BL/6J (B6) × DBA/2J (D2) F2 intercross and the HS4, derived from crossing the B6, D2, BALB/cJ and LP/J strains, were used. Brain (striatum) gene expression data were obtained using the Illumina Mouse WG 6.1 array, and the data sets were interrogated using a weighted gene co-expression network analysis (WGCNA). RESULTS: Genes reliably detected as expressed were similar in all three data sets as was the variability of expression. As measured by the WGCNA, the modular structure of the transcriptome networks was also preserved both on the basis of module assignment and from the perspective of the topological overlap maps. Details of the HS-CC gene modules are provided; essentially identical results were obtained for the HS4 and F2 modules. Gene ontology annotation of the modules revealed a significant overrepresentation in some modules for neuronal processes, e.g., central nervous system development. Integration with known protein-protein interactions data indicated significant enrichment among co-expressed genes. We also noted significant overlap with markers of central nervous system cell types (neurons, oligodendrocytes and astrocytes). Using the Allen Brain Atlas, we found evidence of spatial co-localization within the striatum for several modules. Finally, for some modules it was possible to detect an enrichment of transcription binding sites. The binding site for Wt1, which is associated with neurodegeneration, was the most significantly overrepresented. CONCLUSIONS: Despite the marked differences in genetic diversity, the transcriptome structure was remarkably similar for the F2, HS4 and HS-CC. These data suggest that it should be possible to integrate network data from simple and complex crosses. A careful examination of the HS-CC transcriptome revealed the expected structure for striatal gene expression. Importantly, we demonstrate the integration of anatomical and network expression data.

Electrolytic lesions of the medial nucleus accumbens shell selectively decrease ethanol consumption without altering preference in a limited access procedure in C57BL/6J mice.

The central extended amygdala (cExtA) is a limbic region proposed to play a key role in drug and alcohol addiction and to contain the medial nucleus accumbens shell (MNAc shell). The aim of this study was to examine the involvement of the MNAc shell in ethanol and sucrose consumption in a limited and free access procedure in the C57BL/6J (B6) mouse. Separate groups of mice received bilateral electrolytic lesions of the MNAc shell or sham surgery, and following recovery from surgery, were allowed to voluntarily consume ethanol (15% v/v) in a 2 h limited access 2-bottle-choice procedure. Following 1 week of limited access ethanol consumption, mice were given 1 week of limited access sucrose consumption. A separate group of lesioned and sham mice were given free access (24 h) to ethanol in a 2-bottle choice procedure and were run in parallel to the mice receiving limited access consumption. Electrolytic lesions of the MNAc shell decreased ethanol (but not sucrose) consumption in a limited access procedure, but did not alter free access ethanol consumption. These results suggest that the MNAc shell is a component of the underlying neural circuitry contributing to limited access alcohol consumption in the B6 mouse.

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.

Alignment of the transcriptome with individual variation in animals selectively bred for High Drinking-In-the-Dark (HDID).

Among animals at risk for excessive ethanol consumption such as the HDID selected mouse lines, there is considerable individual variation in the amount of ethanol consumed and the associated blood ethanol concentrations (BECs). For the HDID lines, this variation occurs even though the residual genetic variation associated with the DID phenotype has been largely exhausted and thus is most likely associated with epigenetic factors. Here we focus on the question of whether the genes associated with individual variation in HDID-1 mice are different from those associated with selection (risk) (Iancu et al., 2013). Thirty-three HDID-1 mice were phenotyped for their BECs at the end of a standard DID trial, were sacrificed 3 weeks later, and RNA-Seq was used to analyze the striatal transcriptome. The data obtained illustrate that there is considerable overlap of the risk and variation gene sets, both focused on the fine-tuning of synaptic plasticity.

Diversity Outbred Mice at 21: Maintaining Allelic Variation in the Face of Selection.

Multi-parent populations (MPPs) capture and maintain the genetic diversity from multiple inbred founder strains to provide a resource for high-resolution genetic mapping through the accumulation of recombination events over many generations. Breeding designs that maintain a large effective population size with randomized assignment of breeders at each generation can minimize the impact of selection, inbreeding, and genetic drift on allele frequencies. Small deviations from expected allele frequencies will have little effect on the power and precision of genetic analysis, but a major distortion could result in reduced power and loss of important functional alleles. We detected strong transmission ratio distortion in the Diversity Outbred (DO) mouse population on chromosome 2, caused by meiotic drive favoring transmission of the WSB/EiJ allele at the R2d2 locus. The distorted region harbors thousands of polymorphisms derived from the seven non-WSB founder strains and many of these would be lost if the sweep was allowed to continue. To ensure the utility of the DO population to study genetic variation on chromosome 2, we performed an artificial selection against WSB/EiJ alleles at the R2d2 locus. Here, we report that we have purged the WSB/EiJ allele from the drive locus while preserving WSB/EiJ alleles in the flanking regions. We observed minimal disruption to allele frequencies across the rest of the autosomal genome. However, there was a shift in haplotype frequencies of the mitochondrial genome and an increase in the rate of an unusual sex chromosome aneuploidy. The DO population has been restored to genome-wide utility for genetic analysis, but our experience underscores that vigilant monitoring of similar genetic resource populations is needed to ensure their long-term utility.

Cosplicing network analysis of mammalian brain RNA-Seq data utilizing WGCNA and Mantel correlations.

Across species and tissues and especially in the mammalian brain, production of gene isoforms is widespread. While gene expression coordination has been previously described as a scale-free coexpression network, the properties of transcriptome-wide isoform production coordination have been less studied. Here we evaluate the system-level properties of cosplicing in mouse, macaque, and human brain gene expression data using a novel network inference procedure. Genes are represented as vectors/lists of exon counts and distance measures sensitive to exon inclusion rates quantifies differences across samples. For all gene pairs, distance matrices are correlated across samples, resulting in cosplicing or cotranscriptional network matrices. We show that networks including cosplicing information are scale-free and distinct from coexpression. In the networks capturing cosplicing we find a set of novel hubs with unique characteristics distinguishing them from coexpression hubs: heavy representation in neurobiological functional pathways, strong overlap with markers of neurons and neuroglia, long coding lengths, and high number of both exons and annotated transcripts. Further, the cosplicing hubs are enriched in genes associated with autism spectrum disorders. Cosplicing hub homologs across eukaryotes show dramatically increasing intronic lengths but stable coding region lengths. Shared transcription factor binding sites increase coexpression but not cosplicing; the reverse is true for splicing-factor binding sites. Genes with protein-protein interactions have strong coexpression and cosplicing. Additional factors affecting the networks include shared microRNA binding sites, spatial colocalization within the striatum, and sharing a chromosomal folding domain. Cosplicing network patterns remain relatively stable across species.

Differential network analysis reveals genetic effects on catalepsy modules.

We performed short-term bi-directional selective breeding for haloperidol-induced catalepsy, starting from three mouse populations of increasingly complex genetic structure: an F2 intercross, a heterogeneous stock (HS) formed by crossing four inbred strains (HS4) and a heterogeneous stock (HS-CC) formed from the inbred strain founders of the Collaborative Cross (CC). All three selections were successful, with large differences in haloperidol response emerging within three generations. Using a custom differential network analysis procedure, we found that gene coexpression patterns changed significantly; importantly, a number of these changes were concordant across genetic backgrounds. In contrast, absolute gene-expression changes were modest and not concordant across genetic backgrounds, in spite of the large and similar phenotypic differences. By inferring strain contributions from the parental lines, we are able to identify significant differences in allelic content between the selected lines concurrent with large changes in transcript connectivity. Importantly, this observation implies that genetic polymorphisms can affect transcript and module connectivity without large changes in absolute expression levels. We conclude that, in this case, selective breeding acts at the subnetwork level, with the same modules but not the same transcripts affected across the three selections.

Intra-amygdalar okadaic acid enhances conditioned taste aversion learning and CREB phosphorylation in rats.

Protein phosphatases (PPs) regulate many substrates implicated in learning and memory. Conditioned taste aversion (CTA) learning, in which animals associate a novel taste paired with a toxin and subsequently avoid the taste, is dependent on several serine/threonine phosphatase substrates and the PP1-binding protein spinophilin. In order to examine the effects of PP1/2A blockade on CTA acquisition and extinction, rats received bilateral infusions of okadaic acid (OA) (100nM, 1microl/hemisphere) or vehicle (0.15M NaCl) into the amygdala either 5min prior to, or 5min after, a single pairing of sodium saccharin (0.125%, 10-min access) and LiCl or NaCl (0.15M, 3ml/kg i.p.). Two-bottle, 24-h preference tests were conducted for 13days to measure CTA expression and extinction. Rats conditioned with saccharin and LiCl showed a decreased preference for saccharin, and OA administered before (but not after) the pairing of saccharin and LiCl resulted in a significantly stronger CTA that did not extinguish over 13days. The enhancement of the CTA was not due to aversive effects of OA, because rats given OA and a pairing of saccharin and NaCl did not acquire a CTA. Finally, OA administration increased levels of phosphorylated CREB immunoreactivity following a CTA trial. Together, these results suggest a critical role for PP1/2A during normal CTA learning. Because CTA learning was enhanced only when OA was given prior to conditioning, phosphatase activity may be a constraint on learning during the taste-toxin interval but not during acquisition and consolidation processes that occur after toxin administration.

Estradiol increases Pet-1 and serotonin transporter mRNA in the midbrain raphe nuclei of ovariectomized rats.

Previous research has shown that estradiol increases the anorexia associated with serotonin (5-HT) neurotransmission. To examine further the putative relationship between estradiol and 5-HT, we investigated whether estradiol increases the expression of Pet-1 and the 5-HT transporter (5-HTT), two genes implicated in the development and regulation of the 5-HT system. Ovariectomized (OVX) rats (n=5-6/group) were treated with 0, 2, or 10 microg estradiol benzoate (EB) in sesame oil on 2 consecutive days. Food intake and body weight were recorded 2 days later when EB-treated rats typically display signs of behavioral estrus (e.g., reduced feeding). Following the collection of behavioral data, rats were perfused, brains were removed, and coronal sections were cut through the midbrain raphe nuclei. Pet-1 and 5-HTT mRNA levels were quantified throughout the dorsal and median raphe nuclei (DRN and MRN) by conducting in situ hybridization on free-floating tissue sections using (35)S-labeled cDNA probes. As expected, EB treatment decreased food intake and body weight on the day that modeled estrus. At this same time, EB treatment increased Pet-1 and 5-HTT mRNA levels within the DRN and MRN. We conclude that a physiologically relevant regimen of estradiol treatment in OVX rats increases Pet-1 and 5-HTT mRNA levels in the midbrain raphe nuclei at a time when the anorexigenic effect of estradiol is apparent. Further studies are required to determine whether the increased expression of Pet-1 and 5-HTT mRNA plays a causal role in the anorexigenic effect of estradiol.

Laser-assisted microdissection.

When analyzing alcohol's effects on the brain, researchers often want to look at small clusters of cells that can be studied in isolation from the surrounding brain tissue rather than at the entire brain or larger brain areas. This implies that relatively small numbers of cells have to be retrieved from the brain and studied in culture or subjected to biochemical analyses. The challenge then becomes how to isolate small numbers of cells from a specific brain region without including unwanted cells. One approach to solving this problem is to use a technology known as laser-assisted microdissection (LMD). This article reviews some of the principles of LMD and its use in alcohol research.

Strategies to study the neuroscience of alcoholism: introduction.

Alcohol use and abuse are widespread in the U.S. population. Moreover, for each drinker, alcohol consumption, particularly at excessive levels, has a vast range of effects on the body. Accordingly, research programs aimed at understanding alcohol's effects on the individual as well as on society are similarly varied and widespread. Much of this research focuses on alcohol's impact on the brain and individual nerve cells (i.e., neurons). A detailed survey of the strategies used to investigate the neural mechanisms associated with alcohol use and abuse would easily fill multiple volumes. Instead, this Special Section provides brief reviews of topics largely associated with two areas of research: 1) What strategies can researchers use to image the acute and chronic effects of alcohol on brain function? 2) How can investigators detect the genes, gene products, and gene networks associated with alcohol-related traits?