Integrative QTL analysis of gene expression and chromatin accessibility identifies multi-tissue patterns of genetic regulation.

Gene transcription profiles across tissues are largely defined by the activity of regulatory elements, most of which correspond to regions of accessible chromatin. Regulatory element activity is in turn modulated by genetic variation, resulting in variable transcription rates across individuals. The interplay of these factors, however, is poorly understood. Here we characterize expression and chromatin state dynamics across three tissues-liver, lung, and kidney-in 47 strains of the Collaborative Cross (CC) mouse population, examining the regulation of these dynamics by expression quantitative trait loci (eQTL) and chromatin QTL (cQTL). QTL whose allelic effects were consistent across tissues were detected for 1,101 genes and 133 chromatin regions. Also detected were eQTL and cQTL whose allelic effects differed across tissues, including local-eQTL for Pik3c2g detected in all three tissues but with distinct allelic effects. Leveraging overlapping measurements of gene expression and chromatin accessibility on the same mice from multiple tissues, we used mediation analysis to identify chromatin and gene expression intermediates of eQTL effects. Based on QTL and mediation analyses over multiple tissues, we propose a causal model for the distal genetic regulation of Akr1e1, a gene involved in glycogen metabolism, through the zinc finger transcription factor Zfp985 and chromatin intermediates. This analysis demonstrates the complexity of transcriptional and chromatin dynamics and their regulation over multiple tissues, as well as the value of the CC and related genetic resource populations for identifying specific regulatory mechanisms within cells and tissues.

Tissue- and strain-specific effects of a genotoxic carcinogen 1,3-butadiene on chromatin and transcription.

Epigenetic effects of environmental chemicals are under intense investigation to fill existing knowledge gaps between environmental/occupational exposures and adverse health outcomes. Chromatin accessibility is one prominent mechanism of epigenetic control of transcription, and understanding of the chemical effects on both could inform the causal role of epigenetic alterations in disease mechanisms. In this study, we hypothesized that baseline variability in chromatin organization and transcription profiles among various tissues and mouse strains influence the outcome of exposure to the DNA damaging chemical 1,3-butadiene. To test this hypothesis, we evaluated DNA damage along with comprehensive quantification of RNA transcripts (RNA-seq), identification of accessible chromatin (ATAC-seq), and characterization of regions with histone modifications associated with active transcription (ChIP-seq for acetylation at histone 3 lysine 27, H3K27ac). We collected these data in the lung, liver, and kidney of mice from two genetically divergent strains, C57BL/6J and CAST/EiJ, that were exposed to clean air or to 1,3-butadiene (~600 ppm) for 2 weeks. We found that tissue effects dominate differences in both gene expression and chromatin states, followed by strain effects. At baseline, xenobiotic metabolism was consistently more active in CAST/EiJ, while immune system pathways were more active in C57BL/6J across tissues. Surprisingly, even though all three tissues in both strains harbored butadiene-induced DNA damage, little transcriptional effect of butadiene was observed in liver and kidney. Toxicologically relevant effects of butadiene in the lung were on the pathways of xenobiotic metabolism and inflammation. We also found that variability in chromatin accessibility across individuals (i.e., strains) only partially explains the variability in transcription. This study showed that variation in the basal states of epigenome and transcriptome may be useful indicators for individuals or tissues susceptible to genotoxic environmental chemicals.

Colonic epithelial miR-31 associates with the development of Crohn's phenotypes.

BACKGROUND: Crohn's disease (CD) is highly heterogeneous, due in large part to variability in cellular processes that underlie the natural history of CD, thereby confounding effective therapy. There is a critical need to advance understanding of the cellular mechanisms that drive CD heterogeneity. METHODS: We performed small RNA sequencing of adult colon tissue from CD and NIBD controls. Colonic epithelial cells and immune cells were isolated from colonic tissues, and microRNA-31 (miR-31) expression was measured. miR-31 expression was measured in colonoid cultures generated from controls and patients with CD. We performed small RNA-sequencing of formalin-fixed paraffin-embedded colon and ileum biopsies from treatment-naive pediatric patients with CD and controls and collected data on disease features and outcomes. RESULTS: Small RNA-sequencing and microRNA profiling in the colon revealed 2 distinct molecular subtypes, each with different clinical associations. Notably, we found that miR-31 expression was a driver of these 2 subtypes and, further, that miR-31 expression was particularly pronounced in epithelial cells. Colonoids revealed that miR-31 expression differences are preserved in this ex vivo system. In adult patients, low colonic miR-31 expression levels at the time of surgery were associated with worse disease outcome as measured by need for an end ileostomy and recurrence of disease in the neoterminal ileum. In pediatric patients, lower miR-31 expression at the time of diagnosis was associated with future development of fibrostenotic ileal CD requiring surgeryCONCLUSIONS. These findings represent an important step forward in designing more effective clinical trials and developing personalized CD therapies. FUNDING: This work was supported by CCF Career Development Award (SZS), R01-ES024983 from NIEHS (SZS and TSF), 1R01DK104828-01A1 from NIDDK (SZS and TSF), P01-DK094779-01A1 from NIDDK (SZS), P30-DK034987 from NIDDK (SZS), 1-16-ACE-47 ADA Pathway Award (PS), UNC Nutrition Obesity Research Center Pilot & Feasibility Grant P30DK056350 (PS), CCF PRO-KIIDS NETWORK (SZS and PS), UNC CGIBD T32 Training Grant from NIDDK (JBB), T32 Training Grant (5T32GM007092-42) from NIGMS (MH), and SHARE from the Helmsley Trust (SZS). The UNC Translational Pathology Laboratory is supported, in part, by grants from the National Cancer Institute (3P30CA016086) and the UNC University Cancer Research Fund (UCRF) (PS).