Post transcriptional modifications of RNA are powerful mechanisms by which eukaryotes expand their genetic diversity. For instance, researchers estimate that most transcripts in humans undergo alternative splicing and alternative polyadenylation. These splicing events produce distinct RNA molecules, which in turn yield distinct protein isoforms and/or influence RNA stability, translation, nuclear export, and RNA/protein cellular localization. Due to their pervasiveness and impact, we hypothesized that alternative splicing and alternative polyadenylation in brain can contribute to a predisposition for voluntary alcohol consumption. Using the HXB/BXH recombinant inbred rat panel (a subset of the Hybrid Rat Diversity Panel), we generated over one terabyte of brain RNA sequencing data (total RNA) and identified novel splice variants (via StringTie) and alternative polyadenylation sites (via aptardi) to determine the transcriptional landscape in the brains of these animals. After establishing an analysis pipeline to ascertain high quality transcripts, we quantitated transcripts and integrated genotype data to identify candidate transcript coexpression networks and individual candidate transcripts associated with predisposition to voluntary alcohol consumption in the two-bottle choice paradigm. For genes that were previously associated with this trait (e.g., Lrap, Ift81, and P2rx4) (Saba et al., Febs. J., 282, 3556-3578, Saba et al., Genes. Brain. Behav., 20, e12698), we were able to distinguish between transcript variants to provide further information about the specific isoforms related to the trait. We also identified additional candidate transcripts associated with the trait of voluntary alcohol consumption (i.e., isoforms of Mapkapk5, Aldh1a7, and Map3k7). Consistent with our previous work, our results indicate that transcripts and networks related to inflammation and the immune system in brain can be linked to voluntary alcohol consumption. Overall, we have established a pipeline for including the quantitation of alternative splicing and alternative polyadenylation variants in the transcriptome in the analysis of the relationship between the transcriptome and complex traits.
- Publikační typ
- časopisecké články MeSH
Brown adipose tissue (BAT) has been suggested to play an important role in lipid and glucose metabolism in rodents and possibly also in humans. In the current study, we used genetic and correlation analyses in the BXH/HXB recombinant inbred (RI) strains, derived from Brown Norway (BN) and spontaneously hypertensive rats (SHR), to identify genetic determinants of BAT function. Linkage analyses revealed a quantitative trait locus (QTL) associated with interscapular BAT mass on chromosome 4 and two closely linked QTLs associated with glucose oxidation and glucose incorporation into BAT lipids on chromosome 2. Using weighted gene coexpression network analysis (WGCNA) we identified 1,147 gene coexpression modules in the BAT from BXH/HXB rats and mapped their module eigengene QTLs. Through an unsupervised analysis, we identified modules related to BAT relative mass and function. The Coral4.1 coexpression module is associated with BAT relative mass (includes Cd36 highly connected gene), and the Darkseagreen coexpression module is associated with glucose incorporation into BAT lipids (includes Hiat1, Fmo5, and Sort1 highly connected transcripts). Because multiple statistical criteria were used to identify candidate modules, significance thresholds for individual tests were not adjusted for multiple comparisons across modules. In summary, a systems genetic analysis using genomic and quantitative transcriptomic and physiological information has produced confirmation of several known genetic factors and significant insight into novel genetic components functioning in BAT and possibly contributing to traits characteristic of the metabolic syndrome.
- MeSH
- genetická predispozice k nemoci genetika MeSH
- glukosa metabolismus MeSH
- hnědá tuková tkáň metabolismus MeSH
- krysa rodu rattus MeSH
- lokus kvantitativního znaku genetika MeSH
- metabolický syndrom genetika metabolismus MeSH
- potkani inbrední BN MeSH
- potkani inbrední SHR MeSH
- zvířata MeSH
- Check Tag
- krysa rodu rattus MeSH
- mužské pohlaví MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, N.I.H., Extramural MeSH
BACKGROUND: A statistical pipeline was developed and used for determining candidate genes and candidate gene coexpression networks involved in 2 alcohol (i.e., ethanol [EtOH]) metabolism phenotypes, namely alcohol clearance and acetate area under the curve in a recombinant inbred (RI) (HXB/BXH) rat panel. The approach was also used to provide an indication of how EtOH metabolism can impact the normal function of the identified networks. METHODS: RNA was extracted from alcohol-naïve liver tissue of 30 strains of HXB/BXH RI rats. The reconstructed transcripts were quantitated, and data were used to construct gene coexpression modules and networks. A separate group of rats, comprising the same 30 strains, were injected with EtOH (2 g/kg) for measurement of blood EtOH and acetate levels. These data were used for quantitative trait loci (QTL) analysis of the rate of EtOH disappearance and circulating acetate levels. The analysis pipeline required calculation of the module eigengene values, the correction of these values with EtOH metabolism rates and acetate levels across the rat strains, and the determination of the eigengene QTLs. For a module to be considered a candidate for determining phenotype, the module eigengene values had to have significant correlation with the strain phenotypic values and the module eigengene QTLs had to overlap the phenotypic QTLs. RESULTS: Of the 658 transcript coexpression modules generated from liver RNA sequencing data, a single module satisfied all criteria for being a candidate for determining the alcohol clearance trait. This module contained 2 alcohol dehydrogenase genes, including the gene whose product was previously shown to be responsible for the majority of alcohol elimination in the rat. This module was also the only module identified as a candidate for influencing circulating acetate levels. This module was also linked to the process of generation and utilization of retinoic acid as related to the autonomous immune response. CONCLUSIONS: We propose that our analytical pipeline can successfully identify genetic regions and transcripts which predispose a particular phenotype and our analysis provides functional context for coexpression module components.
- MeSH
- ethanol aplikace a dávkování metabolismus MeSH
- játra účinky léků metabolismus MeSH
- krysa rodu rattus MeSH
- metabolická clearance účinky léků fyziologie MeSH
- multifaktoriální dědičnost účinky léků fyziologie MeSH
- pití alkoholu genetika metabolismus MeSH
- potkani inbrední BN MeSH
- potkani inbrední SHR MeSH
- potkani transgenní MeSH
- strojové učení bez učitele * MeSH
- systémová biologie metody MeSH
- zvířata MeSH
- Check Tag
- krysa rodu rattus MeSH
- mužské pohlaví MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, N.I.H., Extramural MeSH
Gene co-expression analysis has proven to be a powerful tool for ascertaining the organization of gene products into networks that are important for organ function. An organ, such as the liver, engages in a multitude of functions important for the survival of humans, rats, and other animals; these liver functions include energy metabolism, metabolism of xenobiotics, immune system function, and hormonal homeostasis. With the availability of organ-specific transcriptomes, we can now examine the role of RNA transcripts (both protein-coding and non-coding) in these functions. A systems genetic approach for identifying and characterizing liver gene networks within a recombinant inbred panel of rats was used to identify genetically regulated transcriptional networks (modules). For these modules, biological consensus was found between functional enrichment analysis and publicly available phenotypic quantitative trait loci (QTL). In particular, the biological function of two liver modules could be linked to immune response. The eigengene QTLs for these co-expression modules were located at genomic regions coincident with highly significant phenotypic QTLs; these phenotypes were related to rheumatoid arthritis, food preference, and basal corticosterone levels in rats. Our analysis illustrates that genetically and biologically driven RNA-based networks, such as the ones identified as part of this research, provide insight into the genetic influences on organ functions. These networks can pinpoint phenotypes that manifest through the interaction of many organs/tissues and can identify unannotated or under-annotated RNA transcripts that play a role in these phenotypes.
- MeSH
- genová ontologie MeSH
- imunitní systém metabolismus MeSH
- játra imunologie metabolismus MeSH
- lod skóre MeSH
- lokus kvantitativního znaku MeSH
- potkani inbrední SHR MeSH
- RNA genetika metabolismus MeSH
- sekvenční analýza RNA MeSH
- transkriptom MeSH
- vazebná nerovnováha MeSH
- zvířata MeSH
- Check Tag
- mužské pohlaví MeSH
- ženské pohlaví MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, N.I.H., Extramural MeSH