Cord blood DNA methylation of immune and lipid metabolism genes is associated with maternal triglycerides and child adiposity.

OBJECTIVE: Fetal exposures may impact offspring epigenetic signatures and adiposity. The authors hypothesized that maternal metabolic traits associate with cord blood DNA methylation, which, in turn, associates with child adiposity. METHODS: Fasting serum was obtained in 588 pregnant women (27-34 weeks' gestation), and insulin, glucose, high-density lipoprotein cholesterol, triglycerides, and free fatty acids were measured. Cord blood DNA methylation and child adiposity were measured at birth, 4-6 months, and 4-6 years. The association of maternal metabolic traits with DNA methylation (429,246 CpGs) for differentially methylated probes (DMPs) and regions (DMRs) was tested. The association of the first principal component of each DMR with child adiposity was tested, and mediation analysis was performed. RESULTS: Maternal triglycerides were associated with the most DMPs and DMRs of all traits tested (261 and 198, respectively, false discovery rate < 0.05). DMRs were near genes involved in immune function and lipid metabolism. Triglyceride-associated CpGs were associated with child adiposity at 4-6 months (32 CpGs) and 4-6 years (2 CpGs). One, near CD226, was observed at both timepoints, mediating 10% and 22% of the relationship between maternal triglycerides and child adiposity at 4-6 months and 4-6 years, respectively. CONCLUSIONS: DNA methylation may play a role in the association of maternal triglycerides and child adiposity.

Lack of Methylation Changes in GJB2 and RB1 Non-coding Regions of Cochlear Implant Patients with Sensorineural Hearing Loss.

Objective: Recent advances in epigenetic studies continue to reveal novel mechanisms of gene regulation and control, however little is known on the role of epigenetics in sensorineural hearing loss (SNHL) in humans. We aimed to investigate the methylation patterns of two regions, one in RB1 and another in GJB2 in Filipino patients with SNHL compared to hearing control individuals. Methods: We investigated an RB1 promoter region that was previously identified as differentially methylated in children with SNHL and lead exposure. Additionally, we investigated a sequence in an enhancer-like region within GJB2 that contains four CpGs in close proximity. Bisulfite conversion was performed on salivary DNA samples from 15 children with SNHL and 45 unrelated ethnically-matched individuals. We then performed methylation-specific real-time PCR analysis (qMSP) using TaqMan® probes to determine percentage methylation of the two regions. Results: Using qMSP, both our cases and controls had zero methylation at the targeted GJB2 and RB1 regions. Conclusion: Our study showed no changes in methylation at the selected CpG regions in RB1 and GJB2 in the two comparison groups with or without SNHL. This may be due to a lack of environmental exposures to these target regions. Other epigenetic marks may be present around these regions as well as those of other HL-associated genes.

Spatially distinct molecular patterns of gene expression in idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a heterogeneous disease that is pathologically characterized by areas of normal-appearing lung parenchyma, active fibrosis (transition zones including fibroblastic foci) and dense fibrosis. Defining transcriptional differences between these pathologically heterogeneous regions of the IPF lung is critical to understanding the distribution and extent of fibrotic lung disease and identifying potential therapeutic targets. Application of a spatial transcriptomics platform would provide more detailed spatial resolution of transcriptional signals compared to previous single cell or bulk RNA-Seq studies. METHODS: We performed spatial transcriptomics using GeoMx Nanostring Digital Spatial Profiling on formalin-fixed paraffin-embedded (FFPE) tissue from 32 IPF and 12 control subjects and identified 231 regions of interest (ROIs). We compared normal-appearing lung parenchyma and airways between IPF and controls with histologically normal lung tissue, as well as histologically distinct regions within IPF (normal-appearing lung parenchyma, transition zones containing fibroblastic foci, areas of dense fibrosis, and honeycomb epithelium metaplasia). RESULTS: We identified 254 differentially expressed genes (DEGs) between IPF and controls in histologically normal-appearing regions of lung parenchyma; pathway analysis identified disease processes such as EIF2 signaling (important for cap-dependent mRNA translation), epithelial adherens junction signaling, HIF1α signaling, and integrin signaling. Within IPF, we identified 173 DEGs between transition and normal-appearing lung parenchyma and 198 DEGs between dense fibrosis and normal lung parenchyma; pathways dysregulated in both transition and dense fibrotic areas include EIF2 signaling pathway activation (upstream of endoplasmic reticulum (ER) stress proteins ATF4 and CHOP) and wound healing signaling pathway deactivation. Through cell deconvolution of transcriptome data and immunofluorescence staining, we confirmed loss of alveolar parenchymal signals (AGER, SFTPB, SFTPC), gain of secretory cell markers (SCGB3A2, MUC5B) as well as dysregulation of the upstream regulator ATF4, in histologically normal-appearing tissue in IPF. CONCLUSIONS: Our findings demonstrate that histologically normal-appearing regions from the IPF lung are transcriptionally distinct when compared to similar lung tissue from controls with histologically normal lung tissue, and that transition zones and areas of dense fibrosis within the IPF lung demonstrate activation of ER stress and deactivation of wound healing pathways.

Compartment-specific protein interactions in beryllium lung disease.

The study provides insights into proteins that may be relevant in BeS and CBD. It provides a framework to investigate the global changes in lung compartment-specific inflammatory cells to better understand the potential interplay of proteins in CBD. https://bit.ly/3PLNTXC.

Analysis of DNA methylation at birth and in childhood reveals changes associated with season of birth and latitude.

BACKGROUND: Seasonal variations in environmental exposures at birth or during gestation are associated with numerous adult traits and health outcomes later in life. Whether DNA methylation (DNAm) plays a role in the molecular mechanisms underlying the associations between birth season and lifelong phenotypes remains unclear. METHODS: We carried out epigenome-wide meta-analyses within the Pregnancy And Childhood Epigenetic Consortium to identify associations of DNAm with birth season, both at differentially methylated probes (DMPs) and regions (DMRs). Associations were examined at two time points: at birth (21 cohorts, N = 9358) and in children aged 1-11 years (12 cohorts, N = 3610). We conducted meta-analyses to assess the impact of latitude on birth season-specific associations at both time points. RESULTS: We identified associations between birth season and DNAm (False Discovery Rate-adjusted p values < 0.05) at two CpGs at birth (winter-born) and four in the childhood (summer-born) analyses when compared to children born in autumn. Furthermore, we identified twenty-six differentially methylated regions (DMR) at birth (winter-born: 8, spring-born: 15, summer-born: 3) and thirty-two in childhood (winter-born: 12, spring and summer: 10 each) meta-analyses with few overlapping DMRs between the birth seasons or the two time points. The DMRs were associated with genes of known functions in tumorigenesis, psychiatric/neurological disorders, inflammation, or immunity, amongst others. Latitude-stratified meta-analyses [higher (≥ 50°N), lower (< 50°N, northern hemisphere only)] revealed differences in associations between birth season and DNAm by birth latitude. DMR analysis implicated genes with previously reported links to schizophrenia (LAX1), skin disorders (PSORS1C, LTB4R), and airway inflammation including asthma (LTB4R), present only at birth in the higher latitudes (≥ 50°N). CONCLUSIONS: In this large epigenome-wide meta-analysis study, we provide evidence for (i) associations between DNAm and season of birth that are unique for the seasons of the year (temporal effect) and (ii) latitude-dependent variations in the seasonal associations (spatial effect). DNAm could play a role in the molecular mechanisms underlying the effect of birth season on adult health outcomes.

Integrative genetic and genomic networks identify microRNA associated with COPD and ILD.

Chronic obstructive pulmonary disease (COPD) and interstitial lung disease (ILD) are clinically and molecularly heterogeneous diseases. We utilized clustering and integrative network analyses to elucidate roles for microRNAs (miRNAs) and miRNA isoforms (isomiRs) in COPD and ILD pathogenesis. Short RNA sequencing was performed on 351 lung tissue samples of COPD (n = 145), ILD (n = 144) and controls (n = 64). Five distinct subclusters of samples were identified including 1 COPD-predominant cluster and 2 ILD-predominant clusters which associated with different clinical measurements of disease severity. Utilizing 262 samples with gene expression and SNP microarrays, we built disease-specific genetic and expression networks to predict key miRNA regulators of gene expression. Members of miR-449/34 family, known to promote airway differentiation by repressing the Notch pathway, were among the top connected miRNAs in both COPD and ILD networks. Genes associated with miR-449/34 members in the disease networks were enriched among genes that increase in expression with airway differentiation at an air-liquid interface. A highly expressed isomiR containing a novel seed sequence was identified at the miR-34c-5p locus. 47% of the anticorrelated predicted targets for this isomiR were distinct from the canonical seed sequence for miR-34c-5p. Overexpression of the canonical miR-34c-5p and the miR-34c-5p isomiR with an alternative seed sequence down-regulated NOTCH1 and NOTCH4. However, only overexpression of the isomiR down-regulated genes involved in Ras signaling such as CRKL and GRB2. Overall, these findings elucidate molecular heterogeneity inherent across COPD and ILD patients and further suggest roles for miR-34c in regulating disease-associated gene-expression.

Adipocyte hypertrophy in mesenchymal stem cells from infants of mothers with obesity.

OBJECTIVE: Fat content of adipocytes derived from infant umbilical cord mesenchymal stem cells (MSCs) predicts adiposity in children through 4 to 6 years of age. This study tested the hypothesis that MSCs from infants born to mothers with obesity (Ob-MSCs) exhibit adipocyte hypertrophy and perturbations in genes regulating adipogenesis compared with MSCs from infants of mothers with normal weight (NW-MSCs). METHODS: Adipogenesis was induced in MSCs embedded in three-dimensional hydrogel structures, and cell size and number were measured by three-dimensional imaging. Proliferation and protein markers of proliferation and adipogenesis in undifferentiated and adipocyte differentiating cells were measured. RNA sequencing was performed to determine pathways linked to adipogenesis phenotype. RESULTS: In undifferentiated MSCs, greater zinc finger protein (Zfp)423 protein content was observed in Ob- versus NW-MSCs. Adipocytes from Ob-MSCs were larger but fewer than adipocytes from NW-MSCs. RNA sequencing analysis showed that Zfp423 protein correlated with mRNA expression of genes enriched for cell cycle, MSC lineage specification, inflammation, and metabolism pathways. MSC proliferation was not different before differentiation but declined faster in Ob-MSCs upon adipogenic induction. CONCLUSIONS: Ob-MSCs have an intrinsic propensity for adipocyte hypertrophy and reduced hyperplasia during adipogenesis in vitro, perhaps linked to greater Zfp423 content and changes in cell cycle pathway gene expression.

Prenatal exposures to per- and polyfluoroalkyl substances and epigenetic aging in umbilical cord blood: The Healthy Start study.

BACKGROUND: Per- and polyfluoroalkyl substances (PFAS) are ubiquitous, environmentally persistent chemicals, and prenatal exposures have been associated with adverse child health outcomes. Prenatal PFAS exposure may lead to epigenetic age acceleration (EAA), defined as the discrepancy between an individual's chronologic and epigenetic or biological age. OBJECTIVES: We estimated associations of maternal serum PFAS concentrations with EAA in umbilical cord blood DNA methylation using linear regression, and a multivariable exposure-response function of the PFAS mixture using Bayesian kernel machine regression. METHODS: Five PFAS were quantified in maternal serum (median: 27 weeks of gestation) among 577 mother-infant dyads from a prospective cohort. Cord blood DNA methylation data were assessed with the Illumina HumanMethylation450 array. EAA was calculated as the residuals from regressing gestational age on epigenetic age, calculated using a cord-blood specific epigenetic clock. Linear regression tested for associations between each maternal PFAS concentration with EAA. Bayesian kernel machine regression with hierarchical selection estimated an exposure-response function for the PFAS mixture. RESULTS: In single pollutant models we observed an inverse relationship between perfluorodecanoate (PFDA) and EAA (-0.148 weeks per log-unit increase, 95% CI: -0.283, -0.013). Mixture analysis with hierarchical selection between perfluoroalkyl carboxylates and sulfonates indicated the carboxylates had the highest group posterior inclusion probability (PIP), or relative importance. Within this group, PFDA had the highest conditional PIP. Univariate predictor-response functions indicated PFDA and perfluorononanoate were inversely associated with EAA, while perfluorohexane sulfonate had a positive association with EAA. CONCLUSIONS: Maternal mid-pregnancy serum concentrations of PFDA were negatively associated with EAA in cord blood, suggesting a pathway by which prenatal PFAS exposures may affect infant development. No significant associations were observed with other PFAS. Mixture models suggested opposite directions of association between perfluoroalkyl sulfonates and carboxylates. Future studies are needed to determine the importance of neonatal EAA for later child health outcomes.

Multi-Omic Signatures of Sarcoidosis and Progression in Bronchoalveolar Lavage Cells.

Introduction: Sarcoidosis is a heterogeneous, granulomatous disease that can prove difficult to diagnose, with no accurate biomarkers of disease progression. Therefore, we profiled and integrated the DNA methylome, mRNAs, and microRNAs to identify molecular changes associated with sarcoidosis and disease progression that might illuminate underlying mechanisms of disease and potential genomic biomarkers. Methods: Bronchoalveolar lavage cells from 64 sarcoidosis subjects and 16 healthy controls were used. DNA methylation was profiled on Illumina HumanMethylationEPIC arrays, mRNA by RNA-sequencing, and miRNAs by small RNA-sequencing. Linear models were fit to test for effect of diagnosis and phenotype, adjusting for age, sex, and smoking. We built a supervised multi-omics model using a subset of features from each dataset. Results: We identified 46,812 CpGs, 1,842 mRNAs, and 5 miRNAs associated with sarcoidosis versus controls and 1 mRNA, SEPP1 - a protein that supplies selenium to cells, associated with disease progression. Our integrated model emphasized the prominence of the PI3K/AKT1 pathway in sarcoidosis, which is important in T cell and mTOR function. Novel immune related genes and miRNAs including LYST, RGS14, SLFN12L, and hsa-miR-199b-5p, distinguished sarcoidosis from controls. Our integrated model also demonstrated differential expression/methylation of IL20RB, ABCC11, SFSWAP, AGBL4, miR-146a-3p, and miR-378b between non-progressive and progressive sarcoidosis. Conclusions: Leveraging the DNA methylome, transcriptome, and miRNA-sequencing in sarcoidosis BAL cells, we detected widespread molecular changes associated with disease, many which are involved in immune response. These molecules may serve as diagnostic/prognostic biomarkers and/or drug targets, although future testing will be required for confirmation.

Sex-based differences in placental DNA methylation profiles related to gestational age: an NIH ECHO meta-analysis.

The placenta undergoes many changes throughout gestation to support the evolving needs of the foetus. There is also a growing appreciation that male and female foetuses develop differently in utero, with unique epigenetic changes in placental tissue. Here, we report meta-analysed sex-specific associations between gestational age and placental DNA methylation from four cohorts in the National Institutes of Health (NIH) Environmental influences on Child Health Outcomes (ECHO) Programme (355 females/419 males, gestational ages 23-42 weeks). We identified 407 cytosine-guanine dinucleotides (CpGs) in females and 794 in males where placental methylation levels were associated with gestational age. After cell-type adjustment, 55 CpGs in females and 826 in males were significant. These were enriched for biological processes critical to the immune system in females and transmembrane transport in males. Our findings are distinct between the sexes: in females, associations with gestational age are largely explained by differences in placental cellular composition, whereas in males, gestational age is directly associated with numerous alterations in methylation levels.

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.

Cluster analysis of transcriptomic datasets to identify endotypes of idiopathic pulmonary fibrosis.

BACKGROUND: Considerable clinical heterogeneity in idiopathic pulmonary fibrosis (IPF) suggests the existence of multiple disease endotypes. Identifying these endotypes would improve our understanding of the pathogenesis of IPF and could allow for a biomarker-driven personalised medicine approach. We aimed to identify clinically distinct groups of patients with IPF that could represent distinct disease endotypes. METHODS: We co-normalised, pooled and clustered three publicly available blood transcriptomic datasets (total 220 IPF cases). We compared clinical traits across clusters and used gene enrichment analysis to identify biological pathways and processes that were over-represented among the genes that were differentially expressed across clusters. A gene-based classifier was developed and validated using three additional independent datasets (total 194 IPF cases). FINDINGS: We identified three clusters of patients with IPF with statistically significant differences in lung function (p=0.009) and mortality (p=0.009) between groups. Gene enrichment analysis implicated mitochondrial homeostasis, apoptosis, cell cycle and innate and adaptive immunity in the pathogenesis underlying these groups. We developed and validated a 13-gene cluster classifier that predicted mortality in IPF (high-risk clusters vs low-risk cluster: HR 4.25, 95% CI 2.14 to 8.46, p=3.7×10-5). INTERPRETATION: We have identified blood gene expression signatures capable of discerning groups of patients with IPF with significant differences in survival. These clusters could be representative of distinct pathophysiological states, which would support the theory of multiple endotypes of IPF. Although more work must be done to confirm the existence of these endotypes, our classifier could be a useful tool in patient stratification and outcome prediction in IPF.

Incidence and Progression of Fibrotic Lung Disease in an At-Risk Cohort.

Rationale: Relatives of patients with familial interstitial pneumonia (FIP) are at increased risk for pulmonary fibrosis and develop preclinical pulmonary fibrosis (PrePF). Objectives: We defined the incidence and progression of new-onset PrePF and its relationship to survival among first-degree relatives of families with FIP. Methods: This is a cohort study of family members with FIP who were initially screened with a health questionnaire and chest high-resolution computed tomography (HRCT) scan, and approximately 4 years later, the evaluation was repeated. A total of 493 asymptomatic first-degree relatives of patients with FIP were evaluated at baseline, and 296 (60%) of the original subjects participated in the subsequent evaluation. Measurements and Main Results: The median interval between HRCTs was 3.9 years (interquartile range, 3.5-4.4 yr). A total of 252 subjects who agreed to repeat evaluation were originally determined not to have PrePF at baseline; 16 developed PrePF. A conservative estimate of the annual incidence of PrePF is 1,023 per 100,000 person-years (95% confidence interval, 511-1,831 per 100,000 person-years). Of 44 subjects with PrePF at baseline, 38.4% subjects had worsening dyspnea compared with 15.4% of those without PrePF (P = 0.002). Usual interstitial pneumonia by HRCT (P < 0.0002) and baseline quantitative fibrosis score (P < 0.001) are also associated with worsening dyspnea. PrePF at the initial screen is associated with decreased survival (P < 0.001). Conclusions: The incidence of PrePF in this at-risk population is at least 100-fold higher than that reported for sporadic idiopathic pulmonary fibrosis (IPF). Although PrePF and IPF represent distinct entities, our study demonstrates that PrePF, like IPF, is progressive and associated with decreased survival.

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.

Multiomic Signatures of Chronic Beryllium Disease Bronchoalveolar Lavage Cells Relate to T-Cell Function and Innate Immunity.

Chronic beryllium disease (CBD) is a Th1 granulomatous lung disease preceded by sensitization to beryllium (BeS). We profiled the methylome, transcriptome, and selected proteins in the lung to identify molecular signatures and networks associated with BeS and CBD. BAL cell DNA and RNA were profiled using microarrays from CBD (n = 30), BeS (n = 30), and control subjects (n = 12). BAL fluid proteins were measured using Olink Immune Response Panel proteins from CBD (n = 22) and BeS (n = 22) subjects. Linear models identified features associated with CBD, adjusting for covariation and batch effects. Multiomic integration methods identified correlated features between datasets. We identified 1,546 differentially expressed genes in CBD versus control subjects and 204 in CBD versus BeS. Of the 101 shared transcripts, 24 have significant cis relationships between gene expression and DNA methylation, assessed using expression quantitative trait methylation analysis, including genes not previously identified in CBD. A multiomic model of top DNA methylation and gene expression features demonstrated that the first component separated CBD from other samples and the second component separated control subjects from remaining samples. The top features on component one were enriched for T-lymphocyte function, and the top features on component two were enriched for innate immune signaling. We identified six differentially abundant proteins in CBD versus BeS, with two (SIT1 and SH2D1A) selected as important RNA features in the multiomic model. Our integrated analysis of DNA methylation, gene expression, and proteins in the lung identified multiomic signatures of CBD that differentiated it from BeS and control subjects.

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.

Interleukin-6-dependent epithelial fluidization initiates fibrotic lung remodeling.

Chronic disease results from the failure of tissues to maintain homeostasis. In the lung, coordinated repair of the epithelium is essential for preserving homeostasis. In animal models and human lung disease, airway epithelial cells mobilize in response to lung injury, resulting in the formation of airway-like cysts with persistent loss of functional cell types and parenchymal architecture. Using live-cell imaging of human lung epithelial cultures and mouse precision-cut lung slices, we demonstrated that distal airway epithelia are aberrantly fluidized both after injury and in fibrotic lung disease. Through transcriptomic profiling and pharmacologic stimulation of epithelial cultures, we identified interleukin-6 (IL-6) signaling as a driver of tissue fluidization. This signaling cascade occurred independently of canonical Janus kinase (JAK)-signal transducer and activator of transcription (STAT) signaling but instead was dependent on a downstream SRC family kinase (SFK)-yes-associated protein (YAP) axis. Airway epithelial-fibroblast cocultures revealed that the fibrotic mesenchyme acts as a source of IL-6 family cytokines, which drive airway fluidization. Inhibition of the IL-6-SFK-YAP cascade was sufficient to prevent fluidization in both in vitro and ex vivo models. Last, we demonstrated a reduction in fibrotic lung remodeling in mice through genetic or pharmacologic targeting of IL-6-related signaling. Together, our findings illustrate the critical role of airway epithelial fluidization in coordinating the balance between homeostatic lung repair and fibrotic airspace remodeling.