The Puerto Rican Infant Metagenomic and Epidemiologic Study of Respiratory Outcomes (PRIMERO): Design and Baseline Characteristics for a Birth Cohort Study of Early-life Viral Respiratory Illnesses and Airway Dysfunction in Puerto Rican Children.

Epidemiologic studies demonstrate an association between early-life respiratory illnesses (RIs) and the development of childhood asthma. However, it remains uncertain whether these children are predisposed to both conditions or if early-life RIs induce alterations in airway function, immune responses, or other human biology that contribute to the development of asthma. Puerto Rican children experience a disproportionate burden of early-life RIs and asthma, making them an important population for investigating this complex interplay. PRIMERO, the Puerto Rican Infant Metagenomics and Epidemiologic Study of Respiratory Outcomes , recruited pregnant women and their newborns to investigate how the airways develop in early life among infants exposed to different viral RIs, and will thus provide a critical understanding of childhood asthma development. As the first asthma birth cohort in Puerto Rico, PRIMERO will prospectively follow 2,100 term healthy infants. Collected samples include post-term maternal peripheral blood, infant cord blood, the child's peripheral blood at the year two visit, and the child's nasal airway epithelium, collected using minimally invasive nasal swabs, at birth, during RIs over the first two years of life, and at annual healthy visits until age five. Herein, we describe the study's design, population, recruitment strategy, study visits and procedures, and primary outcomes.

Data Science Curriculum in the iField.

Many disciplines, including the broad Field of Information (iField), have been offering Data Science (DS) programs. There have been significant efforts exploring an individual discipline's identity and unique contributions to the broader DS education landscape. To advance DS education in the iField, the iSchool Data Science Curriculum Committee (iDSCC) was formed and charged with building and recommending a DS education framework for iSchools. This paper reports on the research process and findings of a series of studies to address important questions: What is the iField identity in the multidisciplinary DS education landscape? What is the status of DS education in iField schools? What knowledge and skills should be included in the core curriculum for iField DS education? What are the jobs available for DS graduates from the iField? What are the differences between graduate-level and undergraduate-level DS education? Answers to these questions will not only distinguish an iField approach to DS education but also define critical components of DS curriculum. The results will inform individual DS programs in the iField to develop curriculum to support undergraduate and graduate DS education in their local context.

Multi-omic approach associates blood methylome with bronchodilator drug response in pediatric asthma.

BACKGROUND: Albuterol is the drug most widely used as asthma treatment among African Americans despite having a lower bronchodilator drug response (BDR) than other populations. Although BDR is affected by gene and environmental factors, the influence of DNA methylation is unknown. OBJECTIVE: This study aimed to identify epigenetic markers in whole blood associated with BDR, study their functional consequences by multi-omic integration, and assess their clinical applicability in admixed populations with a high asthma burden. METHODS: We studied 414 children and young adults (8-21 years old) with asthma in a discovery and replication design. We performed an epigenome-wide association study on 221 African Americans and replicated the results on 193 Latinos. Functional consequences were assessed by integrating epigenomics with genomics, transcriptomics, and environmental exposure data. Machine learning was used to develop a panel of epigenetic markers to classify treatment response. RESULTS: We identified 5 differentially methylated regions and 2 CpGs genome-wide significantly associated with BDR in African Americans located in FGL2 (cg08241295, P = 6.8 × 10-9) and DNASE2 (cg15341340, P = 7.8 × 10-8), which were regulated by genetic variation and/or associated with gene expression of nearby genes (false discovery rate < 0.05). The CpG cg15341340 was replicated in Latinos (P = 3.5 × 10-3). Moreover, a panel of 70 CpGs showed good classification for those with response and nonresponse to albuterol therapy in African American and Latino children (area under the receiver operating characteristic curve for training, 0.99; for validation, 0.70-0.71). The DNA methylation model showed similar discrimination as clinical predictors (P > .05). CONCLUSIONS: We report novel associations of epigenetic markers with BDR in pediatric asthma and demonstrate for the first time the applicability of pharmacoepigenetics in precision medicine of respiratory diseases.

Nasal airway transcriptome-wide association study of asthma reveals genetically driven mucus pathobiology.

To identify genetic determinants of airway dysfunction, we performed a transcriptome-wide association study for asthma by combining RNA-seq data from the nasal airway epithelium of 681 children, with UK Biobank genetic association data. Our airway analysis identified 95 asthma genes, 58 of which were not identified by transcriptome-wide association analyses using other asthma-relevant tissues. Among these genes were MUC5AC, an airway mucin, and FOXA3, a transcriptional driver of mucus metaplasia. Muco-ciliary epithelial cultures from genotyped donors revealed that the MUC5AC risk variant increases MUC5AC protein secretion and mucus secretory cell frequency. Airway transcriptome-wide association analyses for mucus production and chronic cough also identified MUC5AC. These cis-expression variants were associated with trans effects on expression; the MUC5AC variant was associated with upregulation of non-inflammatory mucus secretory network genes, while the FOXA3 variant was associated with upregulation of type-2 inflammation-induced mucus-metaplasia pathway genes. Our results reveal genetic mechanisms of airway mucus pathobiology.

Epigenome-wide association study of lung function in Latino children and youth with asthma.

INTRODUCTION: DNA methylation studies have associated methylation levels at different CpG sites or genomic regions with lung function. Moreover, genetic ancestry has been associated with lung function in Latinos. However, no epigenome-wide association study (EWAS) of lung function has been performed in this population. Here, we aimed to identify DNA methylation patterns associated with lung function in pediatric asthma among Latinos. RESULTS: We conducted an EWAS in whole blood from 250 Puerto Rican and 148 Mexican American children and young adults with asthma. A total of five CpGs exceeded the genome-wide significance threshold of p = 1.17 × 10-7 in the combined analyses from Puerto Ricans and Mexican Americans: cg06035600 (MAP3K6, p = 6.13 × 10-8) showed significant association with pre-bronchodilator Tiffeneau-Pinelli index, the probes cg00914963 (TBC1D16, p = 1.04 × 10-7), cg16405908 (MRGPRE, p = 2.05 × 10-8), and cg07428101 (MUC2, p = 5.02 × 10-9) were associated with post-bronchodilator forced vital capacity (FVC), and cg20515679 (KCNJ6) with post-bronchodilator Tiffeneau-Pinelli index (p = 1.13 × 10-8). However, these markers did not show significant associations in publicly available data from Europeans (p > 0.05). A methylation quantitative trait loci analysis revealed that methylation levels at these CpG sites were regulated by genetic variation in Latinos and the Biobank-based Integrative Omics Studies (BIOS) consortium. Additionally, two differentially methylated regions in REXOC and AURKC were associated with pre-bronchodilator Tiffeneau-Pinelli index (adjusted p < 0.05) in Puerto Ricans and Mexican Americans. Moreover, we replicated some of the previous differentially methylated signals associated with lung function in non-Latino populations. CONCLUSIONS: We replicated previous associations of epigenetic markers with lung function in whole blood and identified novel population-specific associations shared among Latino subgroups.

Racial/ethnic differences in eligibility for asthma biologics among pediatric populations.

BACKGROUND: Asthma is a heterogeneous disease. Clinical blood parameters differ by race/ethnicity and are used to distinguish asthma subtypes and inform therapies. Differences in subtypes may explain population-specific trends in asthma outcomes. However, these differences in racial/ethnic minority pediatric populations are unclear. OBJECTIVE: We investigated the association of blood parameters and asthma subtypes with asthma outcomes and examined population-specific eligibility for biologic therapies in minority pediatric populations. METHODS: Using data from 2 asthma case-control studies of pediatric minority populations, we performed case-control (N = 3738) and case-only (N = 2743) logistic regressions to quantify the association of blood parameters and asthma subtypes with asthma outcomes. Heterogeneity of these associations was tested using an interaction term between race/ethnicity and each exposure. Differences in therapeutic eligibility were investigated using chi-square tests. RESULTS: Race/ethnicity modified the association between total IgE and asthma exacerbations. Elevated IgE level was associated with worse asthma outcomes in Puerto Ricans. Allergic asthma was associated with worse outcomes in Mexican Americans, whereas eosinophilic asthma was associated with worse outcomes in Puerto Ricans. A lower proportion of Puerto Ricans met dosing criteria for allergic asthma-directed biologic therapy than other groups. A higher proportion of Puerto Ricans qualified for eosinophilic asthma-directed biologic therapy than African Americans. CONCLUSIONS: We found population-specific associations between blood parameters and asthma subtypes with asthma outcomes. Our findings suggest that eligibility for asthma biologic therapies differs across pediatric racial/ethnic populations. These findings call for more studies in diverse populations for equitable treatment of minority patients with asthma.

Outcomes of emergency endovascular versus open repair for abdominal aortic aneurysm rupture.

PURPOSE: Ruptured abdominal aortic aneurysm (rAAA) is one of the most common aortic emergencies in vascular surgery and is associated with high operative mortality and morbidity rates despite recent treatment advances. We evaluated operative mortality risks for the outcomes of emergency endovascular aneurysm repair (eEVAR) vs. open repair in rAAA. METHODS: Twenty patients underwent eEVAR (n = 12) or open repair (n = 8) for rAAA between 2016 and 2020. We adopted the EVAR first strategy since 2018. Primary endpoints included in-hospital mortality and 1-year survival. The outcome variables were analyzed with Fisher exact, Mann-Whitney test, and linear by linear association. The Kaplan-Meier method was used to estimate survival. RESULTS: There were 13 males (65.0%) and the median age of the study cohort was 78.0 years (range, 49-88 years). In-hospital mortality occurred in 7 patients (35.0%); 5 (50.0%) in the early period and 2 (20.0%) in the later period of this series. According to the procedure type, 4 (50.0%) and 3 (25.0%) in-hospital mortalities occurred in the open repair and eEVAR patients, respectively. In 6 patients (50.0%), eEVAR was performed on unfavorable anatomy. The 1-year survival of eEVAR vs. open repair group was 75% ± 12.5% and 50% ± 17.7%, respectively. On univariate analysis, preoperative high-risk indices, postoperative acute renal failure requiring dialysis, pulmonary complications, and prolonged mechanical ventilation were associated with higher operative mortality. CONCLUSION: The current data showed relatively superior outcomes with eEVAR vs. open repair for rAAA, even in some patients with unfavorable anatomy supporting the feasibility, efficacy, and safety of EVAR first strategy.

A genome-wide association study of severe asthma exacerbations in Latino children and adolescents.

Severe asthma exacerbations are a major cause of school absences and healthcare costs in children, particularly those in high-risk racial/ethnic groups.To identify susceptibility genes for severe asthma exacerbations in Latino children and adolescents, we conducted a meta-analysis of genome-wide association studies (GWAS) in 4010 Latino youth with asthma in four independent cohorts, including 1693 Puerto Ricans, 1019 Costa Ricans, 640 Mexicans, 256 Brazilians and 402 members of other Latino subgroups. We then conducted methylation quantitative trait locus, expression quantitative trait locus and expression quantitative trait methylation analyses to assess whether the top single nucleotide polymorphism (SNP) in the meta-analysis is linked to DNA methylation and gene expression in nasal (airway) epithelium in separate cohorts of Puerto Rican and Dutch children and adolescents.In the meta-analysis of GWAS, an SNP in FLJ22447 (rs2253681) was significantly associated with 1.55 increased odds of severe asthma exacerbation (95% CI 1.34-1.79, p=6.3×10-9). This SNP was significantly associated with DNA methylation of a CpG site (cg25024579) at the FLJ22447 locus, which was in turn associated with increased expression of KCNJ2-AS1 in nasal airway epithelium from Puerto Rican children and adolescents (ß=0.10, p=2.18×10-7).SNP rs2253681 was significantly associated with both DNA methylation of a cis-CpG in FLJ22447 and severe asthma exacerbations in Latino youth. This may be partly explained by changes in airway epithelial expression of a gene recently implicated in atopic asthma in Puerto Rican children and adolescents (KCNJ2-AS1).

A genome-wide association study of asthma hospitalizations in adults.

BACKGROUND: Little is known about the genetic determinants of severe asthma exacerbations. OBJECTIVES: We aimed to identify genetic variants associated with asthma hospitalizations. METHODS: We conducted a genome-wide association study of asthma hospitalizations in 34,167 white British adults with asthma, 1,658 of whom had at least 1 asthma-related hospitalization. This analysis was conducted by using logistic regression under an additive genetic model with adjustment for age, sex, body mass index, smoking status, and the first 5 principal components derived from genotypic data. We then analyzed data from 2 cohorts of Latino children and adolescents for replication and conducted quantitative trait locus and functional annotation analyses. RESULTS: At the chromosome 6p21.3 locus, the single-nucleotide polymorphism (SNP) rs56151658 (8 kb from the promoter of HLA-DQB1) was most significantly associated with asthma hospitalizations (for test allele A, odds ratio = 1.36 [95% CI = 1.22-1.52]; P = 3.11 × 10-8); 21 additional SNPs in this locus were associated with asthma hospitalizations at a P value less than 1 × 10-6. In the replication cohorts, multiple SNPs in strong linkage disequilibrium with rs56151658 were associated with severe asthma exacerbations at a P value of .01 or less in the same direction of association as in the discovery cohort. Three HLA genes (HLA-DQA2, HLA-DRB6, and HLA-DOB) were also shown to mediate the estimated effects of the SNPs associated with asthma hospitalizations through effects on gene expression in lung tissue. CONCLUSIONS: We identified strong candidate genes for asthma hospitalizations in adults in the region for class II HLA genes through genomic, quantitative trait locus, and summary data-based mendelian randomization analyses.

Integrative genomic analysis in African American children with asthma finds three novel loci associated with lung function.

Bronchodilator (BD) drugs are commonly prescribed for treatment and management of obstructive lung function present with diseases such as asthma. Administration of BD medication can partially or fully restore lung function as measured by pulmonary function tests. The genetics of baseline lung function measures taken before BD medication have been extensively studied, and the genetics of the BD response itself have received some attention. However, few studies have focused on the genetics of post-BD lung function. To address this gap, we analyzed lung function phenotypes in 1103 subjects from the Study of African Americans, Asthma, Genes, and Environment, a pediatric asthma case-control cohort, using an integrative genomic analysis approach that combined genotype, locus-specific genetic ancestry, and functional annotation information. We integrated genome-wide association study (GWAS) results with an admixture mapping scan of three pulmonary function tests (forced expiratory volume in 1 s [FEV1 ], forced vital capacity [FVC], and FEV1 /FVC) taken before and after albuterol BD administration on the same subjects, yielding six traits. We identified 18 GWAS loci, and five additional loci from admixture mapping, spanning several known and novel lung function candidate genes. Most loci identified via admixture mapping exhibited wide variation in minor allele frequency across genotyped global populations. Functional fine-mapping revealed an enrichment of epigenetic annotations from peripheral blood mononuclear cells, fetal lung tissue, and lung fibroblasts. Our results point to three novel potential genetic drivers of pre- and post-BD lung function: ADAMTS1, RAD54B, and EGLN3.

Long-term Clinical and Hemodynamic Outcomes of Edge-to-Edge Repair for Tricuspid Regurgitation.

BACKGROUND: The edge-to-edge technique (Alfieri stitch) has been widely adopted in mitral valve but not tricuspid valve (TV) repair. We evaluated long-term clinical and hemodynamic outcomes of tricuspid edge-to-edge repair. METHODS: We retrospectively reviewed 237 patients (mean age, 58.4 ± 13.9 years) who had undergone tricuspid edge-to-edge repair from January 2001 to March 2019 in our institution. Tricuspid annuloplasty was performed in 175 patients (73.8%) using ring (91 [38.4%]) or suture (84 [35.4%]) annuloplasty. Concomitant procedures were mitral valve replacement (133 [56.1%]), mitral valve repair (52 [21.9%]), aortic valve replacement (41 [17.3%]), and maze operation (138 [58.2%]). RESULTS: Postoperative echocardiography revealed mild or less tricuspid regurgitation (TR) in 220 patients (92.8%). Early mortality (<30 days) occurred in 9 patients (3.8%), reoperation for bleeding in 16 (6.8%), and low cardiac output syndrome in 15 (6.4%). Freedom from all-cause mortality was 87.2% at 5 years and 80.6% at 10 years. Freedom from moderate or severe TR was 97.1% at 5 years and 84.9% at 10 years. Transtricuspid pressure gradient was 3.8 ± 4.2 mm Hg at discharge and 3.2 ± 5.6 mm Hg at the last follow-up (P = .60). Freedom from significant tricuspid stenosis (transtricuspid pressure gradient ≥5 mm Hg) was 96.4% and 88.4% at 5 and 10 years, respectively. There was 1 early TV reoperation for severe TR on postoperative day 3 (0.4%). Freedom from TV reoperation was 99.1% at 10 years. CONCLUSIONS: Tricuspid edge-to-edge repair showed acceptable long-term clinical and hemodynamic results and is an effective and safe option in TV surgery.

Type 2 and interferon inflammation regulate SARS-CoV-2 entry factor expression in the airway epithelium.

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2, an emerging virus that utilizes host proteins ACE2 and TMPRSS2 as entry factors. Understanding the factors affecting the pattern and levels of expression of these genes is important for deeper understanding of SARS-CoV-2 tropism and pathogenesis. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci for both ACE2 and TMPRSS2, that vary in frequency across world populations. We find TMPRSS2 is part of a mucus secretory network, highly upregulated by type 2 (T2) inflammation through the action of interleukin-13, and that the interferon response to respiratory viruses highly upregulates ACE2 expression. IL-13 and virus infection mediated effects on ACE2 expression were also observed at the protein level in the airway epithelium. Finally, we define airway responses to common coronavirus infections in children, finding that these infections generate host responses similar to other viral species, including upregulation of IL6 and ACE2. Our results reveal possible mechanisms influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

On the cross-population generalizability of gene expression prediction models.

The genetic control of gene expression is a core component of human physiology. For the past several years, transcriptome-wide association studies have leveraged large datasets of linked genotype and RNA sequencing information to create a powerful gene-based test of association that has been used in dozens of studies. While numerous discoveries have been made, the populations in the training data are overwhelmingly of European descent, and little is known about the generalizability of these models to other populations. Here, we test for cross-population generalizability of gene expression prediction models using a dataset of African American individuals with RNA-Seq data in whole blood. We find that the default models trained in large datasets such as GTEx and DGN fare poorly in African Americans, with a notable reduction in prediction accuracy when compared to European Americans. We replicate these limitations in cross-population generalizability using the five populations in the GEUVADIS dataset. Via realistic simulations of both populations and gene expression, we show that accurate cross-population generalizability of transcriptome prediction only arises when eQTL architecture is substantially shared across populations. In contrast, models with non-identical eQTLs showed patterns similar to real-world data. Therefore, generating RNA-Seq data in diverse populations is a critical step towards multi-ethnic utility of gene expression prediction.

Single-Cell and Population Transcriptomics Reveal Pan-epithelial Remodeling in Type 2-High Asthma.

The type 2 cytokine-high asthma endotype (T2H) is characterized by IL-13-driven mucus obstruction of the airways. To further investigate this incompletely understood pathobiology, we characterize IL-13 effects on human airway epithelial cell cultures using single-cell RNA sequencing, finding that IL-13 generates a distinctive transcriptional state for each cell type. Specifically, we discover a mucus secretory program induced by IL-13 in all cell types which converts both mucus and defense secretory cells into a metaplastic state with emergent mucin production and secretion, while leading to ER stress and cell death in ciliated cells. The IL-13-remodeled epithelium secretes a pathologic, mucin-imbalanced, and innate immunity-depleted proteome that arrests mucociliary motion. Signatures of IL-13-induced cellular remodeling are mirrored by transcriptional signatures characteristic of the nasal airway epithelium within T2H versus T2-low asthmatic children. Our results reveal the epithelium-wide scope of T2H asthma and present candidate therapeutic targets for restoring normal epithelial function.

Type 2 and interferon inflammation strongly regulate SARS-CoV-2 related gene expression in the airway epithelium.

Coronavirus disease 2019 (COVID-19) outcomes vary from asymptomatic infection to death. This disparity may reflect different airway levels of the SARS-CoV-2 receptor, ACE2, and the spike protein activator, TMPRSS2. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci (eQTL) for both ACE2 and TMPRSS2, that vary in frequency across world populations. Importantly, we find TMPRSS2 is part of a mucus secretory network, highly upregulated by T2 inflammation through the action of interleukin-13, and that interferon response to respiratory viruses highly upregulates ACE2 expression. Finally, we define airway responses to coronavirus infections in children, finding that these infections upregulate IL6 while also stimulating a more pronounced cytotoxic immune response relative to other respiratory viruses. Our results reveal mechanisms likely influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

Differential asthma odds following respiratory infection in children from three minority populations.

RATIONALE: Severe early-life respiratory illnesses, particularly those caused by respiratory syncytial virus (RSV) and human rhinovirus (HRV), are strongly associated with the development of asthma in children. Puerto Rican children in particular have a strikingly high asthma burden. However, prior studies of the potential associations between early-life respiratory illnesses and asthma in Puerto Rican and other minority populations have been limited. OBJECTIVES: We sought to determine whether early-life respiratory illness was associated with asthma in Puerto Rican, Mexican American, and African American children. METHODS: Using a logistic regression analysis, we examined the association between early-life respiratory illnesses (report of upper respiratory infection (URI), pneumonia, bronchitis, and bronchiolitis/RSV) within the first two years of life and physician-diagnosed asthma after the age of two in a large cohort of Puerto Rican, Mexican American, and African American children. MEASUREMENTS AND MAIN RESULTS: While early-life respiratory illnesses were associated with greater asthma odds in Puerto Ricans, Mexican Americans, and African Americans, these associations were stronger among Puerto Rican children. Specifically, in Puerto Ricans, the odds was 6.15 (95% CI: 4.21-9.05) if the child reported at least one of the following respiratory illness: URI, pneumonia, bronchitis or bronchiolitis. The odds were also higher in Puerto Ricans when considering these conditions separately. CONCLUSIONS: We observed population-specific associations between early-life respiratory illnesses and asthma, which were especially significant and stronger in Puerto Ricans. Taken together with the known high burden of RSV in Puerto Rico, our results may help explain the high burden of asthma in Puerto Ricans.

Lung Function in African American Children with Asthma Is Associated with Novel Regulatory Variants of the KIT Ligand KITLG/SCF and Gene-By-Air-Pollution Interaction.

Baseline lung function, quantified as forced expiratory volume in the first second of exhalation (FEV1), is a standard diagnostic criterion used by clinicians to identify and classify lung diseases. Using whole-genome sequencing data from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine project, we identified a novel genetic association with FEV1 on chromosome 12 in 867 African American children with asthma (P = 1.26 × 10-8, ß = 0.302). Conditional analysis within 1 Mb of the tag signal (rs73429450) yielded one major and two other weaker independent signals within this peak. We explored statistical and functional evidence for all variants in linkage disequilibrium with the three independent signals and yielded nine variants as the most likely candidates responsible for the association with FEV1 Hi-C data and expression QTL analysis demonstrated that these variants physically interacted with KITLG (KIT ligand, also known as SCF), and their minor alleles were associated with increased expression of the KITLG gene in nasal epithelial cells. Gene-by-air-pollution interaction analysis found that the candidate variant rs58475486 interacted with past-year ambient sulfur dioxide exposure (P = 0.003, ß = 0.32). This study identified a novel protective genetic association with FEV1, possibly mediated through KITLG, in African American children with asthma. This is the first study that has identified a genetic association between lung function and KITLG, which has established a role in orchestrating allergic inflammation in asthma.

Genome-Wide Analysis Reveals Mucociliary Remodeling of the Nasal Airway Epithelium Induced by Urban PM2.5.

Air pollution particulate matter <2.5 µm (PM2.5) exposure is associated with poor respiratory outcomes. Mechanisms underlying PM2.5-induced lung pathobiology are poorly understood but likely involve cellular and molecular changes to the airway epithelium. We extracted and chemically characterized the organic and water-soluble components of air pollution PM2.5 samples, then determined the whole transcriptome response of human nasal mucociliary airway epithelial cultures to a dose series of PM2.5 extracts. We found that PM2.5 organic extract (OE), but not water-soluble extract, elicited a potent, dose-dependent transcriptomic response from the mucociliary epithelium. Exposure to a moderate OE dose modified the expression of 424 genes, including activation of aryl hydrocarbon receptor signaling and an IL-1 inflammatory program. We generated an OE-response gene network defined by eight functional enrichment groups, which exhibited high connectivity through CYP1A1, IL1A, and IL1B. This OE exposure also robustly activated a mucus secretory expression program (>100 genes), which included transcriptional drivers of mucus metaplasia (SPDEF and FOXA3). Exposure to a higher OE dose modified the expression of 1,240 genes and further exacerbated expression responses observed at the moderate dose, including the mucus secretory program. Moreover, the higher OE dose significantly increased the MUC5AC/MUC5B gel-forming mucin expression ratio and strongly downregulated ciliated cell expression programs, including key ciliating cell transcription factors (e.g., FOXJ1 and MCIDAS). Chronic OE stimulation induced mucus metaplasia-like remodeling characterized by increases in MUC5AC+ secretory cells and MUC5AC mucus secretions. This epithelial remodeling may underlie poor respiratory outcomes associated with high PM2.5 exposure.