Using a modular massively parallel reporter assay to discover context-specific regulatory grammars in type 2 diabetes.

Recent genome-wide association studies have established that most complex disease-associated loci are found in noncoding regions where defining their function is nontrivial. In this study, we leverage a modular massively parallel reporter assay (MPRA) to uncover sequence features linked to context-specific regulatory activity. We screened enhancer activity across a panel of 198-bp fragments spanning over 10k type 2 diabetes- and metabolic trait-associated variants in the 832/13 rat insulinoma cell line, a relevant model of pancreatic beta cells. We explored these fragments' context sensitivity by comparing their activities when placed up-or downstream of a reporter gene, and in combination with either a synthetic housekeeping promoter (SCP1) or a more biologically relevant promoter corresponding to the human insulin gene ( INS ). We identified clear effects of MPRA construct design on measured fragment enhancer activity. Specifically, a subset of fragments (n = 702/11,656) displayed positional bias, evenly distributed across up- and downstream preference. A separate set of fragments exhibited promoter bias (n = 698/11,656), mostly towards the cell-specific INS promoter (73.4%). To identify sequence features associated with promoter preference, we used Lasso regression with 562 genomic annotations and discovered that fragments with INS promoter-biased activity are enriched for HNF1 motifs. HNF1 family transcription factors are key regulators of glucose metabolism disrupted in maturity onset diabetes of the young (MODY), suggesting genetic convergence between rare coding variants that cause MODY and common T2D-associated regulatory variants. We designed a follow-up MPRA containing HNF1 motif-enriched fragments and observed several instances where deletion or mutation of HNF1 motifs disrupted the INS promoter-biased enhancer activity, specifically in the beta cell model but not in a skeletal muscle cell line, another diabetes-relevant cell type. Together, our study suggests that cell-specific regulatory activity is partially influenced by enhancer-promoter compatibility and indicates that careful attention should be paid when designing MPRA libraries to capture context-specific regulatory processes at disease-associated genetic signals.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation.

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

A Locus on Chromosome 15 Contributes to Acute Ozone-induced Lung Injury in Collaborative Cross Mice.

Ozone (O3)-induced respiratory toxicity varies considerably within the human population and across inbred mouse strains, indicative of gene-environment interactions (GxE). Though previous studies have identified several quantitative trait loci (QTL) and candidate genes underlying responses to O3 exposure, precise mechanisms of susceptibility remain incompletely described. We sought to update our understanding of the genetic architecture of O3 responsiveness using the Collaborative Cross (CC) recombinant inbred mouse panel. We evaluated hallmark O3-induced inflammation and injury phenotypes in 56 CC strains after exposure to filtered air or 2 ppm O3, and performed focused genetic analysis of variation in lung injury, as reflected by protein in lung lavage fluid. Strain-dependent responses to O3 were clear, and QTL mapping revealed two novel loci on Chr (Chromosomes) 10 (peak, 26.2 Mb; 80% confidence interval [CI], 24.6-43.6 Mb) and 15 (peak, 47.1 Mb; 80% CI, 40.2-54.9 Mb), the latter surpassing the 95% significance threshold. At the Chr 15 locus, C57BL/6J and CAST/EiJ founder haplotypes were associated with higher lung injury responses compared with all other CC founder haplotypes. With further statistical analysis and a weight of evidence approach, we delimited the Chr 15 QTL to an ∼2 Mb region containing 21 genes (10 protein coding) and nominated three candidate genes, namely Oxr1, Rspo2, and Angpt1. Gene and protein expression data further supported Oxr1 and Angpt1 as priority candidate genes. In summary, we have shown that O3-induced lung injury is modulated by genetic variation, identified two high priority candidate genes, and demonstrated the value of the CC for detecting GxE.

Integrative analysis reveals mouse strain-dependent responses to acute ozone exposure associated with airway macrophage transcriptional activity.

Acute ozone (O3) exposure is associated with multiple adverse cardiorespiratory outcomes, the severity of which varies across individuals in human populations and inbred mouse strains. However, molecular determinants of response, including susceptibility biomarkers that distinguish who will develop severe injury and inflammation, are not well characterized. We and others have demonstrated that airway macrophages (AMs) are an important resident immune cell type that are functionally and transcriptionally responsive to O3 inhalation. Here, we sought to explore influences of strain, exposure, and strain-by-O3 exposure interactions on AM gene expression and identify transcriptional correlates of O3-induced inflammation and injury across six mouse strains, including five Collaborative Cross (CC) strains. We exposed adult mice of both sexes to filtered air (FA) or 2 ppm O3 for 3 h and measured inflammatory and injury parameters 21 h later. Mice exposed to O3 developed airway neutrophilia and lung injury with strain-dependent severity. In AMs, we identified a common core O3 transcriptional response signature across all strains, as well as a set of genes exhibiting strain-by-O3 exposure interactions. In particular, a prominent gene expression contrast emerged between a low- (CC017/Unc) and high-responding (CC003/Unc) strain, as reflected by cellular inflammation and injury. Further inspection indicated that differences in their baseline gene expression and chromatin accessibility profiles likely contribute to their divergent post-O3 exposure transcriptional responses. Together, these results suggest that aspects of O3-induced respiratory responses are mediated through altered AM transcriptional signatures and further confirm the importance of gene-environment interactions in mediating differential responsiveness to environmental agents.

Ozone-induced changes in the murine lung extracellular vesicle small RNA landscape.

Inhalation exposure to ozone (O3 ) causes adverse respiratory health effects that result from airway inflammation, a complex response mediated in part by changes to airway cellular transcriptional programs. These programs may be regulated by microRNAs transferred between cells (e.g., epithelial cells and macrophages) via extracellular vesicles (EV miRNA). To explore this, we exposed female C57BL/6J mice to filtered air (FA), 1, or 2 ppm O3 by inhalation and collected bronchoalveolar lavage fluid (BALF) 21 h later for markers of airway inflammation, EVs, and EV miRNA. Both concentrations of O3  significantly increased markers of inflammation (neutrophils), injury (total protein), and the number of EV-sized particles in the BALF. Imagestream analysis indicated a substantial portion of particles was positive for canonical EV markers (CD81, CD51), and Siglec-F, a marker of alveolar macrophages. Using high-throughput small RNA sequencing, we identified several differentially expressed (DE) BALF EV miRNAs after 1 ppm (16 DE miRNAs) and 2 ppm (99 DE miRNAs) O3 versus FA exposure. O3 concentration-response patterns in EV miRNA expression were apparent, particularly for miR-2137, miR-126-3p, and miR-351-5p. Integrative analysis of EV miRNA expression and airway cellular mRNA expression identified EV miR-22-3p as a candidate regulator of transcriptomic responses to O3 in airway macrophages. In contrast, we did not identify candidate miRNA regulators of mRNA expression data from conducting airways (predominantly composed of epithelial cells). In summary, our data show that O3 exposure alters EV release and EV miRNA expression, suggesting that further investigation of EVs may provide insight into their effects on airway macrophage function and other mechanisms of O3 -induced respiratory inflammation.

Baseline and innate immune response characterization of a Zfp30 knockout mouse strain.

Airway neutrophilia is correlated with disease severity in a number of chronic and acute pulmonary diseases, and dysregulation of neutrophil chemotaxis can lead to host tissue damage. The gene Zfp30 was previously identified as a candidate regulator of neutrophil recruitment to the lungs and secretion of CXCL1, a potent neutrophil chemokine, in a genome-wide mapping study using the Collaborative Cross. ZFP30 is a putative transcriptional repressor with a KRAB domain capable of inducing heterochromatin formation. Using a CRISPR-mediated knockout mouse model, we investigated the role that Zfp30 plays in recruitment of neutrophils to the lung using models of allergic airway disease and acute lung injury. We found that the Zfp30 null allele did not affect CXCL1 secretion or neutrophil recruitment to the lungs in response to various innate immune stimuli. Intriguingly, despite the lack of neutrophil phenotype, we found there was a significant reduction in the proportion of live Zfp30 homozygous female mutant mice produced from heterozygous matings. This deviation from the expected Mendelian ratios implicates Zfp30 in fertility or embryonic development. Overall, our results indicate that Zfp30 is an essential gene but does not influence neutrophilic inflammation in this particular knockout model.

Transcriptional Profiling of the Murine Airway Response to Acute Ozone Exposure.

Ambient ozone (O3) exposure has serious consequences on respiratory health, including airway inflammation and injury. Decades of research have yielded thorough descriptions of these outcomes; however, less is known about the molecular processes that drive them. The aim of this study was to further describe the cellular and molecular responses to O3 exposure in murine airways, with a particular focus on transcriptional responses in 2 critical pulmonary tissue compartments: conducting airways (CA) and airway macrophages (AM). After exposing adult, female C57BL/6J mice to filtered air, 1 or 2 ppm O3, we assessed hallmark responses including airway inflammation (cell counts and cytokine secretion) and injury (epithelial permeability), followed by gene expression profiling of CA and AM by RNA-seq. As expected, we observed concentration-dependent increases in airway inflammation and injury. Conducting airways and AM both exhibited changes in gene expression to both 1 and 2 ppm O3 that were largely compartment-specific. In CA, genes associated with epithelial barrier function, detoxification processes, and cellular proliferation were altered, while O3 affected genes involved in innate immune signaling, cytokine production, and extracellular matrix remodeling in AM. Further, CA and AM also exhibited notable differences in concentration-response expression patterns for large numbers of genes. Overall, our study has described transcriptional responses to acute O3 exposure, revealing both shared and unique gene expression patterns across multiple concentrations of O3 and in 2 important O3-responsive tissues. These profiles provide broad mechanistic insight into pulmonary O3 toxicity, and reveal a variety of targets for focused follow-up studies.

Molecular classification of Crohn's disease reveals two clinically relevant subtypes.

OBJECTIVE: The clinical presentation and course of Crohn's disease (CD) is highly variable. We sought to better understand the cellular and molecular mechanisms that guide this heterogeneity, and characterise the cellular processes associated with disease phenotypes. DESIGN: We examined both gene expression and gene regulation (chromatin accessibility) in non-inflamed colon tissue from a cohort of adult patients with CD and control patients. To support the generality of our findings, we analysed previously published expression data from a large cohort of treatment-naïve paediatric CD and control ileum. RESULTS: We found that adult patients with CD clearly segregated into two classes based on colon tissue gene expression-one that largely resembled the normal colon and one where certain genes showed expression patterns normally specific to the ileum. These classes were supported by changes in gene regulatory profiles observed at the level of chromatin accessibility, reflective of a fundamental shift in underlying molecular phenotypes. Furthermore, gene expression from the ilea of a treatment-naïve cohort of paediatric patients with CD could be similarly subdivided into colon-like and ileum-like classes. Finally, expression patterns within these CD subclasses highlight large-scale differences in the immune response and aspects of cellular metabolism, and were associated with multiple clinical phenotypes describing disease behaviour, including rectal disease and need for colectomy. CONCLUSIONS: Our results strongly suggest that these molecular signatures define two clinically relevant forms of CD irrespective of tissue sampling location, patient age or treatment status.