Genetic regulators of sputum mucin concentration and their associations with COPD phenotypes.

Hyper-secretion and/or hyper-concentration of mucus is a defining feature of multiple obstructive lung diseases, including chronic obstructive pulmonary disease (COPD). Mucus itself is composed of a mixture of water, ions, salt and proteins, of which the gel-forming mucins, MUC5AC and MUC5B, are the most abundant. Recent studies have linked the concentrations of these proteins in sputum to COPD phenotypes, including chronic bronchitis (CB) and acute exacerbations (AE). We sought to determine whether common genetic variants influence sputum mucin concentrations and whether these variants are also associated with COPD phenotypes, specifically CB and AE. We performed a GWAS to identify quantitative trait loci for sputum mucin protein concentration (pQTL) in the Sub-Populations and InteRmediate Outcome Measures in COPD Study (SPIROMICS, n = 708 for total mucin, n = 215 for MUC5AC, MUC5B). Subsequently, we tested for associations of mucin pQTL with CB and AE using regression modeling (n = 822-1300). Replication analysis was conducted using data from COPDGene (n = 5740) and by examining results from the UK Biobank. We identified one genome-wide significant pQTL for MUC5AC (rs75401036) and two for MUC5B (rs140324259, rs10001928). The strongest association for MUC5B, with rs140324259 on chromosome 11, explained 14% of variation in sputum MUC5B. Despite being associated with lower MUC5B, the C allele of rs140324259 conferred increased risk of CB (odds ratio (OR) = 1.42; 95% confidence interval (CI): 1.10-1.80) as well as AE ascertained over three years of follow up (OR = 1.41; 95% CI: 1.02-1.94). Associations between rs140324259 and CB or AE did not replicate in COPDGene. However, in the UK Biobank, rs140324259 was associated with phenotypes that define CB, namely chronic mucus production and cough, again with the C allele conferring increased risk. We conclude that sputum MUC5AC and MUC5B concentrations are associated with common genetic variants, and the top locus for MUC5B may influence COPD phenotypes, in particular CB.

BPIFB1 loss alters airway mucus properties and diminishes mucociliary clearance.

Airway mucociliary clearance (MCC) is required for host defense and is often diminished in chronic lung diseases. Effective clearance depends upon coordinated actions of the airway epithelium and a mobile mucus layer. Dysregulation of the primary secreted airway mucin proteins, MUC5B and MUC5AC, is associated with a reduction in the rate of MCC; however, how other secreted proteins impact the integrity of the mucus layer and MCC remains unclear. We previously identified the gene Bpifb1/Lplunc1 as a regulator of airway MUC5B protein levels using genetic approaches. Here, we show that BPIFB1 is required for effective MCC in vivo using Bpifb1 knockout (KO) mice. Reduced MCC in Bpifb1 KO mice occurred in the absence of defects in epithelial ion transport or reduced ciliary beat frequency. Loss of BPIFB1 in vivo and in vitro altered biophysical and biochemical properties of mucus that have been previously linked to impaired MCC. Finally, we detected colocalization of BPIFB1 and MUC5B in secretory granules in mice and the protein mesh of secreted mucus in human airway epithelia cultures. Collectively, our findings demonstrate that BPIFB1 is an important component of the mucociliary apparatus in mice and a key component of the mucus protein network.NEW & NOTEWORTHY BPIFB1, also known as LPLUNC1, was found to regulate mucociliary clearance (MCC), a key aspect of host defense in the airway. Loss of this protein was also associated with altered biophysical and biochemical properties of mucus that have been previously linked to impaired MCC.

Collaborative cross strain CC011/UncJ as a novel mouse model of T2-high, severe asthma.

Among asthmatics, there is significant heterogeneity in the clinical presentation and underlying pathophysiological mechanisms, leading to the recognition of multiple disease endotypes (e.g., T2-high vs. T2-low). This heterogeneity extends to severe asthmatics, who may struggle to control symptoms even with high-dose corticosteroid treatment and other therapies. However, there are limited mouse models available to model the spectrum of severe asthma endotypes. We sought to identify a new mouse model of severe asthma by first examining responses to chronic allergen exposure among strains from the Collaborative Cross (CC) mouse genetics reference population, which contains greater genetic diversity than other inbred strain panels previously used for models of asthma. Mice from five CC strains and the often-used classical inbred strain BALB/cJ were chronically exposed to house dust mite (HDM) allergen for five weeks followed by measurements of airway inflammation. CC strain CC011/UncJ (CC011) exhibited extreme responses to HDM including high levels of airway eosinophilia, elevated lung resistance, and extensive airway wall remodeling, and even fatalities among ~ 50% of mice prior to study completion. Compared to BALB/cJ mice, CC011 mice had stronger Th2-mediated airway responses demonstrated by significantly elevated total and HDM-specific IgE and increased Th2 cytokines during tests of antigen recall, but not enhanced ILC2 activation. Airway eosinophilia in CC011 mice was completely dependent upon CD4+ T-cells. Notably, we also found that airway eosinophilia in CC011 mice was resistant to dexamethasone steroid treatment. Thus, the CC011 strain provides a new mouse model of T2-high, severe asthma driven by natural genetic variation likely acting through CD4+ T-cells. Future studies aimed at determining the genetic basis of this phenotype will provide new insights into mechanisms underlying severe asthma.

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.

Obesity reprograms the pulmonary polyunsaturated fatty acid-derived lipidome, transcriptome, and gene-oxylipin networks.

Obesity exacerbates inflammation upon lung injury; however, the mechanisms by which obesity primes pulmonary dysregulation prior to external injury are not well studied. Herein, we tested the hypothesis that obesity dysregulates pulmonary PUFA metabolism that is central to inflammation initiation and resolution. We first show that a high-fat diet (HFD) administered to C57BL/6J mice increased the relative abundance of pulmonary PUFA-containing triglycerides and the concentration of PUFA-derived oxylipins (particularly prostaglandins and hydroxyeicosatetraenoic acids), independent of an increase in total pulmonary PUFAs, prior to onset of pulmonary inflammation. Experiments with a genetic model of obesity (ob/ob) generally recapitulated the effects of the HFD on the pulmonary oxylipin signature. Subsequent pulmonary next-generation RNA sequencing identified complex and unique transcriptional regulation with the HFD. We found the HFD increased pathways related to glycerophospholipid metabolism and immunity, including a unique elevation in B cell differentiation and signaling. Furthermore, we conducted computational integration of lipidomic with transcriptomic data. These analyses identified novel HFD-driven networks between glycerophospholipid metabolism and B cell receptor signaling with specific PUFA-derived pulmonary oxylipins. Finally, we confirmed the hypothesis by demonstrating that the concentration of pulmonary oxylipins, in addition to inflammatory markers, were generally increased in mice consuming a HFD upon ozone-induced acute lung injury. Collectively, these data show that a HFD dysregulates pulmonary PUFA metabolism prior to external lung injury, which may be a mechanism by which obesity primes the lungs to respond poorly to infectious and/or inflammatory challenges.

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.

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.

Development of a large-scale computer-controlled ozone inhalation exposure system for rodents.

Objective: Complete systems for laboratory-based inhalation toxicology studies are typically not commercially available; therefore, inhalation toxicologists utilize custom-made exposure systems. Here we report on the design, construction, testing, operation and maintenance of a newly developed in vivo rodent ozone inhalation exposure system. Materials and methods: Key design requirements for the system included large-capacity exposure chambers to facilitate studies with large sample sizes, automatic and precise control of chamber ozone concentrations, as well as automated data collection on airflow and environmental conditions. The exposure system contains two Hazelton H-1000 stainless steel and glass exposure chambers, each providing capacity for up to 180 mice or 96 rats. We developed an empirically tuned proportional-integral-derivative control loop that provides stable ozone concentrations throughout the exposure period (typically 3h), after a short ramp time (∼8 min), and across a tested concentration range of 0.2-2 ppm. Specific details on the combination of analog and digital input/output system for environmental data acquisition, control and safety systems are provided, and we outline the steps involved in maintenance and calibration of the system. Results: We show that the exposure system produces consistent ozone exposures both within and across experiments, as evidenced by low coefficients of variation in chamber ozone concentration and consistent biological responses (airway inflammation) in mice, respectively. Conclusion: Thus, we have created a large and robust ozone exposure system, facilitating future studies on the health effects of ozone in rodents.

Meta-analysis of airway epithelium gene expression in asthma.

Differential gene expression in the airway epithelium of patients with asthma versus controls has been reported in several studies. However, there is no consensus on which genes are reproducibly affected in asthma. We sought to identify a consensus list of differentially expressed genes (DEGs) using a meta-analysis approach.We identified eight studies with data that met defined inclusion criteria. These studies comprised 355 cases and 193 controls and involved sampling either bronchial or nasal epithelium. We conducted study-level analyses, followed by a meta-analysis. Likewise, we applied a meta-analysis framework to the results of study-level pathway enrichment.We identified 1273 DEGs, 431 of which had not been identified in previous studies. 450 DEGs exhibited large effect sizes and were robust to study population differences in age, sex, race/ethnicity, medication use, smoking status and exacerbations. The magnitude of differential expression of these 450 genes was highly similar in bronchial and nasal airway epithelia. Meta-analysis of pathway enrichment revealed a number of consistently dysregulated biological pathways, including putative transcriptional and post-transcriptional regulators.In total, we identified a set of genes that is consistently dysregulated in asthma, that links to known and novel biological pathways, and that will inform asthma subtype identification.

Identification of microRNAs associated with allergic airway disease using a genetically diverse mouse population.

BACKGROUND: Allergic airway diseases (AADs) such as asthma are characterized in part by granulocytic airway inflammation. The gene regulatory networks that govern granulocyte recruitment are poorly understood, but evidence is accruing that microRNAs (miRNAs) play an important role. To identify miRNAs that may underlie AADs, we used two complementary approaches that leveraged the genotypic and phenotypic diversity of the Collaborative Cross (CC) mouse population. In the first approach, we sought to identify miRNA expression quantitative trait loci (eQTL) that overlap QTL for AAD-related phenotypes. Specifically, CC founder strains and incipient lines of the CC were sensitized and challenged with house dust mite allergen followed by measurement of granulocyte recruitment to the lung. Total lung RNA was isolated and miRNA was measured using arrays for CC founders and qRT-PCR for incipient CC lines. RESULTS: Among CC founders, 92 miRNAs were differentially expressed. We measured the expression of 40 of the most highly expressed of these 92 miRNAs in the incipient lines of the CC and identified 18 eQTL corresponding to 14 different miRNAs. Surprisingly, half of these eQTL were distal to the corresponding miRNAs, and even on different chromosomes. One of the largest-effect local miRNA eQTL was for miR-342-3p, for which we identified putative causal variants by bioinformatic analysis of the effects of single nucleotide polymorphisms on RNA structure. None of the miRNA eQTL co-localized with QTL for eosinophil or neutrophil recruitment. In the second approach, we constructed putative miRNA/mRNA regulatory networks and identified three miRNAs (miR-497, miR-351 and miR-31) as candidate master regulators of genes associated with neutrophil recruitment. Analysis of a dataset from human keratinocytes transfected with a miR-31 inhibitor revealed two target genes in common with miR-31 targets correlated with neutrophils, namely Oxsr1 and Nsf. CONCLUSIONS: miRNA expression in the allergically inflamed murine lung is regulated by genetic loci that are smaller in effect size compared to mRNA eQTL and often act in trans. Thus our results indicate that the genetic architecture of miRNA expression is different from mRNA expression. We identified three miRNAs, miR-497, miR-351 and miR-31, that are candidate master regulators of genes associated with neutrophil recruitment. Because miR-31 is expressed in airway epithelia and is predicted to target genes with known links to neutrophilic inflammation, we suggest that miR-31 is a potentially novel regulator of airway inflammation.

Integrative genetic analysis of allergic inflammation in the murine lung.

Airway allergen exposure induces inflammation among individuals with atopy that is characterized by altered airway gene expression, elevated levels of T helper type 2 cytokines, mucus hypersecretion, and airflow obstruction. To identify the genetic determinants of the airway allergen response, we employed a systems genetics approach. We applied a house dust mite mouse model of allergic airway disease to 151 incipient lines of the Collaborative Cross, a new mouse genetic reference population, and measured serum IgE, airway eosinophilia, and gene expression in the lung. Allergen-induced serum IgE and airway eosinophilia were not correlated. We detected quantitative trait loci (QTL) for airway eosinophilia on chromosome (Chr) 11 (71.802-87.098 megabases [Mb]) and allergen-induced IgE on Chr 4 (13.950-31.660 Mb). More than 4,500 genes expressed in the lung had gene expression QTL (eQTL), the majority of which were located near the gene itself. However, we also detected approximately 1,700 trans-eQTL, and many of these trans-eQTL clustered into two regions on Chr 2. We show that one of these loci (at 147.6 Mb) is associated with the expression of more than 100 genes, and, using bioinformatics resources, fine-map this locus to a 53 kb-long interval. We also use the gene expression and eQTL data to identify a candidate gene, Tlcd2, for the eosinophil QTL. Our results demonstrate that hallmark allergic airway disease phenotypes are associated with distinct genetic loci on Chrs 4 and 11, and that gene expression in the allergically inflamed lung is controlled by both cis and trans regulatory factors.

Genetic analysis of complex traits in the emerging Collaborative Cross.

The Collaborative Cross (CC) is a mouse recombinant inbred strain panel that is being developed as a resource for mammalian systems genetics. Here we describe an experiment that uses partially inbred CC lines to evaluate the genetic properties and utility of this emerging resource. Genome-wide analysis of the incipient strains reveals high genetic diversity, balanced allele frequencies, and dense, evenly distributed recombination sites-all ideal qualities for a systems genetics resource. We map discrete, complex, and biomolecular traits and contrast two quantitative trait locus (QTL) mapping approaches. Analysis based on inferred haplotypes improves power, reduces false discovery, and provides information to identify and prioritize candidate genes that is unique to multifounder crosses like the CC. The number of expression QTLs discovered here exceeds all previous efforts at eQTL mapping in mice, and we map local eQTL at 1-Mb resolution. We demonstrate that the genetic diversity of the CC, which derives from random mixing of eight founder strains, results in high phenotypic diversity and enhances our ability to map causative loci underlying complex disease-related traits.

Multi-omics in nasal epithelium reveals three axes of dysregulation for asthma risk in the African Diaspora populations.

Asthma has striking disparities across ancestral groups, but the molecular underpinning of these differences is poorly understood and minimally studied. A goal of the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) is to understand multi-omic signatures of asthma focusing on populations of African ancestry. RNASeq and DNA methylation data are generated from nasal epithelium including cases (current asthma, N = 253) and controls (never-asthma, N = 283) from 7 different geographic sites to identify differentially expressed genes (DEGs) and gene networks. We identify 389 DEGs; the top DEG, FN1, was downregulated in cases (q = 3.26 × 10-9) and encodes fibronectin which plays a role in wound healing. The top three gene expression modules implicate networks related to immune response (CEACAM5; p = 9.62 × 10-16 and CPA3; p = 2.39 × 10-14) and wound healing (FN1; p = 7.63 × 10-9). Multi-omic analysis identifies FKBP5, a co-chaperone of glucocorticoid receptor signaling known to be involved in drug response in asthma, where the association between nasal epithelium gene expression is likely regulated by methylation and is associated with increased use of inhaled corticosteroids. This work reveals molecular dysregulation on three axes - increased Th2 inflammation, decreased capacity for wound healing, and impaired drug response - that may play a critical role in asthma within the African Diaspora.

Mice with FVB-derived sequence on chromosome 17 succumb to disseminated virus infection due to aberrant NK cell and T cell responses.

Zoonotic arenavirus infections can result in viral hemorrhagic disease, characterized by platelet loss, petechia, and multi-organ injury. The mechanisms governing these outcomes are likely impacted by virus strain and infection dose, as well as an individual's genetic background and immune constitution. To better understand the processes leading to severe pathogenesis, we compared two strains of inbred mice, C57BL/6J (B6) and FVB/NJ (FVB), that have diametrically opposed outcomes during disseminated lymphocytic choriomeningitis virus (LCMV) infection. Infection caused minimal pathogenesis in B6 mice, whereas FVB mice developed acute hepatitis and perished due, in part, to aberrant NK cell and T cell responses. Susceptible mice showed an outgrowth of cytolytic CD4+ T cells and loss of Treg cells. B6 congenic mice with the FVB allele at a 25Mb locus on chromosome 17 recapitulated FVB pathogenesis upon infection. A locus containing a limited number of variants in immune-related genes greatly impacts survival during infection.

Regional Differences in Mucociliary Clearance in the Upper and Lower Airways.

As the nasal cavity is the portal of entry for inspired air in mammals, this region is exposed to the highest concentration of inhaled particulate matter and pathogens, which must be removed to keep the lower airways sterile. Thus, one might expect vigorous removal of these substances via mucociliary clearance (MCC) in this region. We have investigated the rate of MCC in the murine nasal cavity compared to the more distal airways (trachea). The rate of MCC in the nasal cavity (posterior nasopharynx, PNP) was ∼3-4× greater than on the tracheal wall. This appeared to be due to a more abundant population of ciliated cells in the nasal cavity (∼80%) compared to the more sparsely ciliated trachea (∼40%). Interestingly, the tracheal ventral wall exhibited a significantly lower rate of MCC than the tracheal posterior membrane. The trachealis muscle underlying the ciliated epithelium on the posterior membrane appeared to control the surface architecture and likely in part the rate of MCC in this tracheal region. In one of our mouse models (Bpifb1 KO) exhibiting a 3-fold increase in MUC5B protein in lavage fluid, MCC particle transport on the tracheal walls was severely compromised, yet normal MCC occurred on the tracheal posterior membrane. While a blanket of mucus covered the surface of both the PNP and trachea, this mucus appeared to be transported as a blanket by MCC only in the PNP. In contrast, particles appeared to be transported as discrete patches or streams of mucus in the trachea. In addition, particle transport in the PNP was fairly linear, in contrast transport of particles in the trachea often followed a more non-linear route. The thick, viscoelastic mucus blanket that covered the PNP, which exhibited ∼10-fold greater mass of mucus than did the blanket covering the surface of the trachea, could be transported over large areas completely devoid of cells (made by a breach in the epithelial layer). In contrast, particles could not be transported over even a small epithelial breach in the trachea. The thick mucus blanket in the PNP likely aids in particle transport over the non-ciliated olfactory cells in the nasal cavity and likely contributes to humidification and more efficient particle trapping in this upper airway region.

Strain-dependent genomic factors affect allergen-induced airway hyperresponsiveness in mice.

Asthma is etiologically and clinically heterogeneous, making the genomic basis of asthma difficult to identify. We exploited the strain-dependence of a murine model of allergic airway disease to identify different genomic responses in the lung. BALB/cJ and C57BL/6J mice were sensitized with the immunodominant allergen from the Dermatophagoides pteronyssinus species of house dust mite (Der p 1), without exogenous adjuvant, and the mice then underwent a single challenge with Der p 1. Allergic inflammation, serum antibody titers, mucous metaplasia, and airway hyperresponsiveness were evaluated 72 hours after airway challenge. Whole-lung gene expression analyses were conducted to identify genomic responses to allergen challenge. Der p 1-challenged BALB/cJ mice produced all the key features of allergic airway disease. In comparison, C57BL/6J mice produced exaggerated Th2-biased responses and inflammation, but exhibited an unexpected decrease in airway hyperresponsiveness compared with control mice. Lung gene expression analysis revealed genes that were shared by both strains and a set of down-regulated genes unique to C57BL/6J mice, including several G-protein-coupled receptors involved in airway smooth muscle contraction, most notably the M2 muscarinic receptor, which we show is expressed in airway smooth muscle and was decreased at the protein level after challenge with Der p 1. Murine strain-dependent genomic responses in the lung offer insights into the different biological pathways that develop after allergen challenge. This study of two different murine strains demonstrates that inflammation and airway hyperresponsiveness can be decoupled, and suggests that the down-modulation of expression of G-protein-coupled receptors involved in regulating airway smooth muscle contraction may contribute to this dissociation.

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.

Inferring the Allelic Series at QTL in Multiparental Populations.

Multiparental populations (MPPs) are experimental populations in which the genome of every individual is a mosaic of known founder haplotypes. These populations are useful for detecting quantitative trait loci (QTL) because tests of association can leverage inferred founder haplotype descent. It is difficult, however, to determine how haplotypes at a locus group into distinct functional alleles, termed the allelic series. The allelic series is important because it provides information about the number of causal variants at a QTL and their combined effects. In this study, we introduce a fully Bayesian model selection framework for inferring the allelic series. This framework accounts for sources of uncertainty found in typical MPPs, including the number and composition of functional alleles. Our prior distribution for the allelic series is based on the Chinese restaurant process, a relative of the Dirichlet process, and we leverage its connection to the coalescent to introduce additional prior information about haplotype relatedness via a phylogenetic tree. We evaluate our approach via simulation and apply it to QTL from two MPPs: the Collaborative Cross (CC) and the Drosophila Synthetic Population Resource (DSPR). We find that, although posterior inference of the exact allelic series is often uncertain, we are able to distinguish biallelic QTL from more complex multiallelic cases. Additionally, our allele-based approach improves haplotype effect estimation when the true number of functional alleles is small. Our method, Tree-Based Inference of Multiallelism via Bayesian Regression (TIMBR), provides new insight into the genetic architecture of QTL in MPPs.

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.