The MUC5B-associated variant rs35705950 resides within an enhancer subject to lineage- and disease-dependent epigenetic remodeling

Epigenetics as a versatile regulator of fibrosis

Fibrosis, a process caused by excessive deposition of extracellular matrix (ECM), is a common cause and outcome of organ failure and even death. Researchers have made many efforts to understand the mechanism of fibrogenesis and to develop therapeutic strategies; yet, the outcome remains unsatisfactory. In recent years, advances in epigenetics, including chromatin remodeling, histone modification, DNA methylation, and noncoding RNA (ncRNA), have provided more insights into the fibrotic process and have suggested the possibility of novel therapy for organ fibrosis. In this review, we summarize the current research on the epigenetic mechanisms involved in organ fibrosis and their possible clinical applications.

Genomic Fingerprint Associated with Familial Idiopathic Pulmonary Fibrosis: A Review

Idiopathic pulmonary fibrosis (IPF) is a severe interstitial lung disease; although the recent introduction of two anti-fibrosis drugs, pirfenidone and Nidanib, have resulted in a significant reduction in lung function decline, IPF is still not curable. Approximately 2-20% of patients with IPF have a family history of the disease, which is considered the strongest risk factor for idiopathic interstitial pneumonia. However, the genetic predispositions of familial IPF (f-IPF), a particular type of IPF, remain largely unknown. Genetics affect the susceptibility and progression of f-IPF. Genomic markers are increasingly being recognized for their contribution to disease prognosis and drug therapy outcomes. Existing data suggest that genomics may help identify individuals at risk for f-IPF, accurately classify patients, elucidate key pathways involved in disease pathogenesis, and ultimately develop more effective targeted therapies. Since several genetic variants associated with the disease have been found in f-IPF, this review systematically summarizes the latest progress in the gene spectrum of the f-IPF population and the underlying mechanisms of f-IPF. The genetic susceptibility variation related to the disease phenotype is also illustrated. This review aims to improve the understanding of the IPF pathogenesis and facilitate his early detection.

Epithelial Endoplasmic Reticulum Stress Enhances the Risk of Muc5b-associated Lung Fibrosis

The gain-of-function minor allele of the MUC5B (mucin 5B, oligomeric mucus/gel-forming) promoter (rs35705950) is the strongest risk factor for idiopathic pulmonary fibrosis (IPF), a devastating fibrotic lung disease that leads to progressive respiratory failure in adults. We have previously demonstrated that Muc5b overexpression in mice worsens lung fibrosis after bleomycin exposure and have hypothesized that excess Muc5b promotes endoplasmic reticulum (ER) stress and apoptosis, stimulating fibrotic lung injury. Here, we report that ER stress pathway members ATF4 (activating transcription factor 4) and ATF6 coexpress with MUC5B in epithelia of the distal IPF airway and honeycomb cyst and that this is more pronounced in carriers of the gain-of-function MUC5B promoter variant. Similarly, in mice exposed to bleomycin, Muc5b expression is temporally associated with markers of ER stress. Using bulk and single-cell RNA sequencing in bleomycin-exposed mice, we found that pathologic ER stress-associated transcripts Atf4 and Ddit3 (DNA damage inducible transcript 3) were elevated in alveolar epithelia of SFTPC-Muc5b transgenic (SFTPC-Muc5bTg) mice relative to wild-type (WT) mice. Activation of the ER stress response inhibits protein translation for most genes by phosphorylation of Eif2α (eukaryotic translation initiation factor 2 alpha), which prevents guanine exchange by Eif2B and facilitates translation of Atf4. The integrated stress response inhibitor (ISRIB) facilitates interaction of phosphorylated Eif2α with Eif2B, overcoming translation inhibition associated with ER stress and reducing Atf4. We found that a single dose of ISRIB diminished Atf4 translation in SFTPC-Muc5bTg mice after bleomycin injury. Moreover, ISRIB resolved the exaggerated fibrotic response of SFTPC-Muc5bTg mice to bleomycin. In summary, we demonstrate that MUC5B and Muc5b expression is associated with pathologic ER stress and that restoration of normal translation with a single dose of ISRIB promotes lung repair in bleomycin-injured Muc5b-overexpressing mice.

MUC5B rs35705950 minor allele associates with older age and better survival in idiopathic pulmonary fibrosis

BACKGROUND AND OBJECTIVE

The minor T-allele of the MUC5B promoter polymorphism rs35705950 is strongly associated with idiopathic pulmonary fibrosis (IPF). However, conflicting results have been reported on the relationship between the MUC5B minor allele and survival and it is unknown whether a specific subgroup of IPF patients might benefit from MUC5B minor allele carriage. We investigated the association between MUC5B rs35705950, survival and patient characteristics in a real-world population of European IPF patients.

METHODS

In this retrospective study, 1751 patients with IPF from 8 European centres were included. MUC5B rs35705950 genotype, demographics, clinical characteristics at diagnosis and survival data were analysed.

RESULTS

In a multi-variate Cox proportional hazard model the MUC5B minor allele was a significant independent predictor of survival when adjusted for age, sex, high resolution computed tomography pattern, smoking behaviour and pulmonary function tests in IPF. MUC5B minor allele carriers were significantly older at diagnosis (p = 0.001). The percentage of MUC5B minor allele carriers increased significantly with age from 44% in patients aged <56 year, to 63% in patients aged >75. In IPF patients aged <56, the MUC5B minor allele was not associated with survival. In IPF patients aged ≥56, survival was significantly better for MUC5B minor allele carriers (45 months [CI: 42-49]) compared to non-carriers (29 months [CI: 26-33]; p = 4 × 10^-12 ).

CONCLUSION

MUC5B minor allele carriage associates with a better median transplant-free survival of 16 months in the European IPF population aged over 56 years. MUC5B genotype status might aid disease prognostication in clinical management of IPF patients.

Genetics in Idiopathic Pulmonary Fibrosis: A Clinical Perspective

BACKGROUND

Unraveling the genetic background in a significant proportion of patients with both sporadic and familial IPF provided new insights into the pathogenic pathways of pulmonary fibrosis.

AIM

The aim of the present study is to overview the clinical significance of genetics in IPF.

PERSPECTIVE

It is fascinating to realize the so-far underestimated but dynamically increasing impact that genetics has on aspects related to the pathophysiology, accurate and early diagnosis, and treatment and prevention of this devastating disease. Genetics in IPF have contributed as no other in unchaining the disease from the dogma of a "a sporadic entity of the elderly, limited to the lungs" and allowed all scientists, but mostly clinicians, all over the world to consider its many aspects and "faces" in all age groups, including its co-existence with several extra pulmonary conditions from cutaneous albinism to bone-marrow and liver failure.

CONCLUSION

By providing additional evidence for unsuspected characteristics such as immunodeficiency, impaired mucus, and surfactant and telomere maintenance that very often co-exist through the interaction of common and rare genetic variants in the same patient, genetics have created a generous and pluralistic yet unifying platform that could lead to the understanding of the injurious and pro-fibrotic effects of many seemingly unrelated extrinsic and intrinsic offending factors. The same platform constantly instructs us about our limitations as well as about the heritability, the knowledge and the wisdom that is still missing.

A MUC5B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (N_cases = 4,325; N_controls = 507,640; OR = 0.89 [0.82-0.97]; P = 6.86 × 10^-3) and joint meta-analyses with the HGI (N_cases = 13,320; N_controls = 1,508,841; OR, 0.90 [0.86-0.95]; P = 8.99 × 10^-5). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (N_cases = 19,168/N_controls = 492,854; OR, 0.98 [0.95-1.01]; P = 0.06) but was nominally significant (P < 0.05) in the joint meta-analysis with the HGI (N_cases = 44,820; N_controls = 1,775,827; OR, 0.97 [0.95-1.00]; P = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

Colocalization of Gene Expression and DNA Methylation with Genetic Risk Variants Supports Functional Roles of MUC5B and DSP in Idiopathic Pulmonary Fibrosis

Rationale: Common genetic variants have been associated with idiopathic pulmonary fibrosis (IPF). Objectives: To determine functional relevance of the 10 IPF-associated common genetic variants we previously identified. Methods: We performed expression quantitative trait loci (eQTL) and methylation quantitative trait loci (mQTL) mapping, followed by co-localization of eQTL and mQTL with genetic association signals and functional validation by luciferase reporter assays. Illumina multi-ethnic genotyping arrays, mRNA sequencing, and Illumina 850k methylation arrays were performed on lung tissue of participants with IPF (234 RNA and 345 DNA samples) and non-diseased controls (188 RNA and 202 DNA samples). Measurements and Main Results: Focusing on genetic variants within 10 IPF-associated genetic loci, we identified 27 eQTLs in controls and 24 eQTLs in cases (false-discovery-rate-adjusted P < 0.05). Among these signals, we identified associations of lead variants rs35705950 with expression of MUC5B and rs2076295 with expression of DSP in both cases and controls. mQTL analysis identified CpGs in gene bodies of MUC5B (cg17589883) and DSP (cg08964675) associated with the lead variants in these two loci. We also demonstrated strong co-localization of eQTL/mQTL and genetic signal in MUC5B (rs35705950) and DSP (rs2076295). Functional validation of the mQTL in MUC5B using luciferase reporter assays demonstrates that the CpG resides within a putative internal repressor element. Conclusions: We have established a relationship of the common IPF genetic risk variants rs35705950 and rs2076295 with respective changes in MUC5B and DSP expression and methylation. These results provide additional evidence that both MUC5B and DSP are involved in the etiology of IPF.

Dysregulated Cell-Cell Communication Characterizes Pulmonary Fibrosis

Idiopathic pulmonary fibrosis (IPF) is a progressive disease of older adults characterized by fibrotic replacement of functional gas exchange units in the lung. The strongest risk factor for IPF is a genetic variantin the promoter region of the gel-forming mucin, MUC5B. To better understand how the MUC5B variant influences development of fibrosis, we used the NicheNet R package and leveraged publicly available single-cell RNA sequencing data to identify and evaluate how epithelia participating in gas exchange are influenced by ligands expressed in control, MUC5B variant, and fibrotic environments. We observed that loss of type-I alveolar epithelia (AECI) characterizes the single-cell RNA transcriptome in fibrotic lung and validated the pattern of AECI loss using single nuclear RNA sequencing. Examining AECI transcriptomes, we found enrichment of transcriptional signatures for IL6 and AREG, which we have previously shown to mediate aberrant epithelial fluidization in IPF and murine bleomycin models. Moreover, we found that the protease ADAM17, which is upstream of IL6 trans-signaling, was enriched in control MUC5B variant donors. We used immunofluorescence to validate a role for enhanced expression of ADAM17 among MUC5B variants, suggesting involvement in IPF pathogenesis and maintenance.

Integrated genomics approaches identify transcriptional mediators and epigenetic responses to Afghan desert particulate matter in small airway epithelial cells

Military Deployment to Southwest Asia and Afghanistan and exposure to toxic airborne particulates have been associated with an increased risk of developing respiratory disease, collectively termed deployment-related respiratory diseases (DRRDs). Our knowledge about how particulates mediate respiratory disease is limited, precluding the appropriate recognition or management. Central to this limitation is the lack of understanding of how exposures translate into dysregulated cell identity with dysregulated transcriptional programs. The small airway epithelium is involved in both the pathobiology of DRRD and fine particulate matter deposition. To characterize small airway epithelial cell epigenetic and transcriptional responses to Afghan desert particulate matter (APM) and investigate the functional interactions of transcription factors that mediate these responses, we applied two genomics assays, the assay for transposase accessible chromatin with sequencing (ATAC-seq) and Precision Run-on sequencing (PRO-seq). We identified activity changes in a series of transcriptional pathways as candidate regulators of susceptibility to subsequent insults, including signal-dependent pathways, such as loss of cytochrome P450 or P53/P63, and lineage-determining transcription factors, such as GRHL2 loss or TEAD3 activation. We further demonstrated that TEAD3 activation was unique to APM exposure despite similar inflammatory responses when compared with wood smoke particle exposure and that P53/P63 program loss was uniquely positioned at the intersection of signal-dependent and lineage-determining transcriptional programs. Our results establish the utility of an integrated genomics approach in characterizing responses to exposures and identifying genomic targets for the advanced investigation of the pathogenesis of DRRD.

COVID-19 Survivor Patients Carrying the Rs35705950 Risk Allele in MUC5B Have Higher Plasma Levels of Mucin 5B

BACKGROUND

Genetic susceptibility to infectious diseases is partly due to the variation in the human genome, and COVID-19 is not the exception. This study aimed to identify whether risk alleles of known genes linked with emphysema (SERPINA1) and pulmonary fibrosis (MUC5B) are associated with severe COVID-19, and whether plasma mucin 5B differs according to patients' outcomes.

MATERIALS AND METHODS

We included 1258 Mexican subjects diagnosed with COVID-19. We genotyped rs2892474 and rs17580 of the SERPINA1 gene and rs35705950 of MUC5B. Based on the rs35705950 genotypes, mucin 5B plasma protein levels were quantified.

RESULTS

Homozygous for the risk alleles of the three polymorphisms were found in less than 5% of the study population, but no statistically significant difference in the genotype or allele association analysis. At the protein level, non-survivors carrying one or two copies of the risk allele rs35705950 in MUC5B (GT + TT) had lower levels of mucin 5B compared to the survivors (0.0 vs. 0.17 ng/mL, p = 0.0013).

CONCLUSION

The polymorphisms rs28929474 and rs17580 of SERPINA1 and rs35705950 of MUC5B are not associated with the risk of severe COVID-19 in the Mexican population. COVID-19 survivor patients bearing one or two copies of the rs35705950 risk allele have higher plasma levels of mucin 5B.

Macrophage programming is regulated by a cooperative interaction between fatty acid binding protein 5 and peroxisome proliferator-activated receptor γ

Resolution of inflammation is an active process that is tightly regulated to achieve repair and tissue homeostasis. In the absence of resolution, persistent inflammation underlies the pathogenesis of chronic lung disease such as chronic obstructive pulmonary disease (COPD) with recurrent exacerbations. Over the course of inflammation, macrophage programming transitions from pro-inflammatory to pro-resolving, which is in part regulated by the nuclear receptor Peroxisome Proliferator-Activated Receptor γ (PPARγ). Our previous work demonstrated an association between Fatty Acid Binding Protein 5 (FABP5) expression and PPARγ activity in peripheral blood mononuclear cells of healthy and COPD patients. However, a role for FABP5 in macrophage programming has not been examined. Here, using a combination of in vitro and in vivo approaches, we demonstrate that FABP5 is necessary for PPARγ activation. In turn, PPARγ acts directly to increase FABP5 expression in primary human alveolar macrophages. We further illustrate that lack of FABP5 expression promotes a pro-inflammatory macrophage programming with increased secretion of pro-inflammatory cytokines and increased chromatin accessibility for pro-inflammatory transcription factors (e.g., NF-κB and MAPK). And finally, real-time cell metabolic analysis using the Seahorse technology shows an inhibition of oxidative phosphorylation in FABP5-deficient macrophages. Taken together, our data indicate that FABP5 and PPARγ reciprocally regulate each other's expression and function, consistent with a novel positive feedback loop between the two factors that mediates macrophage pro-resolving programming. Our studies highlight the importance of defining targets and regulatory mechanisms that control the resolution of inflammation and may serve to inform novel interventional strategies directed towards COPD.

The MUC5B Promoter Polymorphism Associates With Severe COVID-19 in the European Population

Background: Diversity in response on exposure to severe acute respiratory syndrome coronavirus 2 may be related to the innate immune response in the elderly. The mucin MUC5B is an important component of the innate immune response and expression levels are associated with the MUC5B promoter polymorphism, rs35705950. The high expressing T-allele is a risk allele for the non-infectious aging lung disease idiopathic pulmonary fibrosis (IPF). We investigated if MUC5B rs35705950 associates with severe COVID-19.

Methods: In this retrospective candidate gene case-control study we recruited 108 Dutch patients (69% male, median age 66 years, 77% white) requiring hospitalization for COVID-19 (22% ICU stay, 24% died). For validation, genotypes were obtained from the UK-Biobank (n = 436, 57% male, median age 70 years, 27% died), for replication data from the severe COVID-19 GWAS group from Italy (n = 835) and Spain (n = 775) was used, each with a control cohort (n = 356,735, n = 1,255, n = 950, respectively). MUC5B association analysis was performed including adjustment for age and sex.

Results: The rs35705950 T-allele frequency was significantly lower in Dutch white patients (n = 83) than in controls (0.04 vs. 0.10; p = 0.02). This was validated in the UK biobank cohort (0.08 vs. 0.11; p = 0.001). While age and sex differed significantly between cases and control, comparable results were obtained with age and sex as confounding variables in a multivariate analysis. The association was replicated in the Italian (p = 0.04), and Spanish (p = 0.03) case-control cohorts. Meta-analysis showed a negative association for the T-allele with COVID-19 (OR = 0.75 (CI: 0.67–0.85); p = 6.63 × 10⁻⁶).

Conclusions: This study shows that carriage of the T-allele of MUC5B rs35705950 confers protection from development of severe COVID-19. Because the T-allele is a known risk allele for IPF, this study provides further evidence for the existence of trade-offs between optimal mucin expression levels in the aging lung.

Deconvolution of multiplexed transcriptional responses to wood smoke particles defines rapid aryl hydrocarbon receptor signaling dynamics

The heterogeneity of respirable particulates and compounds complicates our understanding of transcriptional responses to air pollution. Here, we address this by applying precision nuclear run-on sequencing and the assay for transposase-accessible chromatin sequencing to measure nascent transcription and chromatin accessibility in airway epithelial cells after wood smoke particle (WSP) exposure. We used transcription factor enrichment analysis to identify temporally distinct roles for ternary response factor-serum response factor complexes, the aryl hydrocarbon receptor (AHR), and NFκB in regulating transcriptional changes induced by WSP. Transcription of canonical targets of the AHR, such as CYP1A1 and AHRR, was robustly increased after just 30 min of WSP exposure, and we discovered novel AHR-regulated pathways and targets including the DNA methyltransferase, DNMT3L. Transcription of these genes and associated enhancers rapidly returned to near baseline by 120 min after exposure. The kinetics of AHR- and NFκB-regulated responses to WSP were distinguishable based on the timing of both transcriptional responses and chromatin remodeling, with induction of several cytokines implicated in maintaining NFκB-mediated responses through 120 min of exposure. In aggregate, our data establish a direct and primary role for AHR in mediating airway epithelial responses to WSP and identify crosstalk between AHR and NFκB signaling in controlling proinflammatory gene expression. This work also defines an integrated genomics-based strategy for deconvoluting multiplexed transcriptional responses to heterogeneous environmental exposures.