Most cited article - PubMed ID 27401171
Uncovering the liver's role in immunity through RNA co-expression networks
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.
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.
- Keywords
- Alcohol Metabolism, HXB/BXH Recombinant Inbred Rat Panel, Liver, Quantitative Trait Locus Mapping, RNA Sequencing, Weighted Gene Coexpression Network Analysis,
- MeSH
- Ethanol administration & dosage metabolism MeSH
- Liver drug effects metabolism MeSH
- Rats MeSH
- Metabolic Clearance Rate drug effects physiology MeSH
- Multifactorial Inheritance drug effects physiology MeSH
- Alcohol Drinking genetics metabolism MeSH
- Rats, Inbred BN MeSH
- Rats, Inbred SHR MeSH
- Rats, Transgenic MeSH
- Unsupervised Machine Learning * MeSH
- Systems Biology methods MeSH
- Animals MeSH
- Check Tag
- Rats MeSH
- Male MeSH
- Animals MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
- Research Support, N.I.H., Extramural MeSH
- Names of Substances
- Ethanol 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.
- Keywords
- brown adipose tissue, coexpression modules, quantitative trait locus, recombinant inbred strains, spontaneously hypertensive rat,
- MeSH
- Genetic Predisposition to Disease genetics MeSH
- Glucose metabolism MeSH
- Adipose Tissue, Brown metabolism MeSH
- Rats MeSH
- Quantitative Trait Loci genetics MeSH
- Metabolic Syndrome genetics metabolism MeSH
- Rats, Inbred BN MeSH
- Rats, Inbred SHR MeSH
- Animals MeSH
- Check Tag
- Rats MeSH
- Male MeSH
- Animals MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
- Research Support, N.I.H., Extramural MeSH
- Names of Substances
- Glucose MeSH
