Potential Clinical Utility of MUC5B und TOLLIP Single Nucleotide Polymorphisms (SNPs) in the Management of Patients with IPF

doi:10.21203/rs.3.rs-133655/v1

Download PDF

Research

Potential Clinical Utility of MUC5B und TOLLIP Single Nucleotide Polymorphisms (SNPs) in the Management of Patients with IPF

https://doi.org/10.21203/rs.3.rs-133655/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 27 Feb, 2021

Read the published version in Orphanet Journal of Rare Diseases →

You are reading this older preprint version

Read the latest preprint version →

Background: Genetic variants of TOLLIP and MUC5B, both on chromosome 11, have been reported to be associated with the development and/or prognosis of idiopathic pulmonary fibrosis (IPF). This retrospective study was conducted to investigate the association of MUC5B and TOLLIP SNPs with disease outcome in IPF. 62 IPF patients and 50 healthy controls (HC) from our Institution were genotyped for SNPs within MUC5B (rs35705950) and TOLLIP (rs3750920 and rs5743890). Correlation of SNPs genotypes with survival, acute exacerbation (AE) or disease progression (defined as a decline of ≥ 5% in FVC and or >10% in DLco in one year) was investigated.

Results: The MUC5B rs35705950 minor allele (T) was more frequent in IPF subjects than in HC (35% vs 9% p<0.001). TOLLIP SNPs alleles and genotype distribution did not differ between IPF and HC and did not vary according to gender, age, BMI and lung functional impairment at baseline. The minor allele (C) in TOLLIP rs5743890 was associated with worse survival and with disease progression in all performed analyses. The MUC5B rs35705950 or the TOLLIP rs3750920 minor allele, were not associated with disease progression or acute exacerbation.

Conclusion: We confirm that the minor allele of MUC5B rs35705950 is associated with IPF. The minor allele of TOLLIP rs5743890 appears to be a predictor of worse survival and more rapid disease progression, therefore being of potential utility to stratify IPF patients at baseline.

Molecular Genetics

Epidemiology

IPF

MUC5B

TOLLIP

disease progression

Idiopathic pulmonary fibrosis is a chronic, progressing and ultimately fatal lung disease, characterized by loss of lung volume and development of respiratory insufficiency over time. Whereas lung functional and imaging parameters are currently used to monitor disease progression in the clinical routine, circulating biomarkers still have limited clinical impact due to issues with reproducibility and inter-individual variability(1). On the other hand, data from genetic investigations in interstitial lung diseases, mainly in IPF, indicate that single nucleotide polymorphisms (SNPs) associated with the development of lung fibrosis correlate with disease outcome or response to certain treatments(2–4).

MUC5B affects the rheological properties of airway mucus, mucociliary transport, and airway defense(5). The MUC5B promoter variant rs35705950 has been validated as a contributor to the expression of MUC5B in the lung(6). This polymorphism may be related to the methylation in the MUC5B promoter region and the disruption of the activities of some transcriptional factors(7). The minor allele (T) in rs35705950 SNP confers the risk of lung fibrosis(8) but patients with the TT genotype show better survival (3, 6). It has been hypothesized that persistent bronchiolar epithelial injury and the over- production of MUC5B by airway progenitor cells result in the development of honeycomb cysts (UIP Pattern) and IPF (9).

Genetic variants of TOLLIP have also been reported to be associated with the development and/or prognosis of IPF (2). The TOLLIP protein is an inhibitor of the toll-like receptors (TLRs) 2 and 4which both are active in the lung, thereby suppressing tumor necrosis factor-a (TNF-a) and IL-6 production(10). The rs5743890 minor allele is associated with a 50% reduction in TOLLIP production whereas rs3750920, located within TOLLIP exon 3, has been linked to decreased TOLLIP mRNA expression and differential responses to treatment in IPF (4). In particular, significant interaction was observed between N-acetylcysteine (NAC) therapy and rs3750920 within TOLLIP. A trend toward harm was observed in patients with a CC genotype treated with NAC, but not in those with a TT genotype, suggesting that a genotype-stratified prospective clinical trial should be conducted before any recommendation regarding the use of off-label NAC to treat IPF (4).

The aim of this study was to investigate the association of MUC5B and TOLLIP SNPs with disease course and outcome in IPF patients.

Study Subjects

We retrospectively studied 62 patients with IPF followed in our Institution from 2012 to 2019. As a comparison group, 50 healthy subjects (HC) were included. All subjects were caucasians. Diagnosis of IPF was made according to the ATS/ERS guidelines 2011 (11) and 2018 (12). Two patients had familiar IPF. The comorbidity burden was assessed using the Charlson comorbidity Index (CCI), which scores 19 comorbidities and age to provide a total score between 0 and 37(13).

Definition of Disease Progression in ILD

Disease progression from baseline was defined as a deterioration of self-reported symptoms (worsening of dyspnoea, cough) and relative decrease in FVC > 10% pred., and/or a decrease in DLco ≥ 15% pred. and/or increase in existing or appearance of new densities compatible with IPF at HRCT (14, 15). Otherwise, the patients were defined as stable. Acute exacerbation was defined according to the ATS/ERS consensus paper (16).

DNA extraction and Genotyping.

Genomic DNA was extracted from peripheral blood leukocytes by using a silica-membrane- based nucleic acid purification kit (Qiagen DNA Mini Kit, Qiagen, USA), and stored at − 80 °C before use. Three SNPs at the chr11p15.5 locus were chosen for this analysis: rs35705950 (position chr11:1241221), rs5743890 (position chr11:1325829) and rs3750920 (position chr11:1309956).

Genotyping of the MUC5B rs35705950, TOLLIP rs3750920 and rs5743890 single-nucleotide polymorphism involved a commercially available SNP genotyping assay (TaqMan SNP Genotyping Assay) and the Applied Biosystems 7500 Fast RT-PCR System (Life Technologies Corp. Carlsbad, California, USA) as reported previously (4, 17, 18).

Pulmonary Function Tests and Blood Gas Analysis

Measurements including forced vital capacity (FVC), forced expiratory volume in one second (FEV₁), total lung capacity (TLC), diffusing capacity of the lung for carbon monoxide (DLco), partial pressure of oxygen in arterial blood (PaO₂), arterial carbon dioxide tension (PCO₂), arterial oxygen saturation (SaO₂) and alveolar-arterial oxygen gradient (AaDO₂) were performed at the time of the blood sample collection. The protocol for pulmonary function measurements conformed to ATS recommendations(19). Values were expressed as percentages of predicted normal values.

Statistical Analysis

Continuous variables were evaluated for a normal distribution with the Kolmogorov-Smirnov test. Parametric data are presented as mean ± SEM. Categorical variables are presented as either a percentage of the total or numerically, as appropriate. Spearman’s or Pearson’s correlation coefficient was obtained for correlations. Haplotype frequencies for multiple loci were estimated by means of HPlus v2.5 software for biallelic markers. The measures of pair-wise linkage disequilibrium between SNPs were calculated by using the expectation-maximization (EM) algorithm with SNPAlyze (Version 5.1; Dynacom Co., Ltd., Yokohama, Japan). Chi-square test and Fisher’s exact test were used to test for deviation from Hardy-Weinberg equilibrium (HWE).

Univariate and multivariate Cox proportional hazard regression model was conducted to study the independent effect of age, gender, smoking history, ethnicity, MUC5B and TOLLIP genotypes, and pulmonary function tests on the disease course. The Kaplan-Meier method with the log-rank test was used to analyze whether SNP genotypes were associated with the disease outcome. P values of < 0.05 were considered statistically significant. All statistical analyses were performed using SPSS 27.0 (SPSS Inc., Chicago, IL, USA).

Characteristics of Study Subjects

Demographics and clinical characteristics of the studied subjects are shown in Table 1. Comorbidities are shown separately in Supplement Table 1 (Table S1). Follow up time of IPF patients was 7.5 ± 4.5 yrs. Forty-five patients (72.5%) received antifibrotic treatment, 15 (24%) received N-acetylcysteine (NAC) monotherapy and 2 no treatment over time. Sixteen patients (26%) died and 12 (19%) underwent lung transplantation during follow-up. 19 (31%) patients experienced acute exacerbation. Median survival time, defined as time to death or lung transplantation, was 5.9 ± 3.2 years.

MUC5B and TOLLIP SNPs allele and genotypes distribution

Allele and genotype distribution for each SNP in IPF and HC is shown In Table 2. MUC5B rs35705950, TOLLIP rs5743890 and rs3750920 alleles were all in HWE. The frequency of MUC5B rs35705950 minor allele (T) was significantly different between IPF patients and HC (p < 0.001). For TOLLIP rs5743890 and rs3750920 the frequency of minor alleles did not differ between IPF and HC (Table 2). Genotype distribution did not differ significantly between IPF and HC for the three SNPs studied.

Four major haplotypes with three polymorphic sites were established. The haplotype TTT occurred more frequently and GTT less frequently in IPF patients than HC (Table S2).

Measures of pair-wise linkage disequilibrium between MUC5B and TOLLIP SNPs are shown in Table S3. MUC5B rs35705950 and TOLLIP rs5743890 were in high linkage disequilibrium (D´-1, R² 0.038, p = 0.0035).

Correlations between SNPs and clinical characteristics

None of the SNP genotypes was associated with gender, BMI, smoking habits, lung function, BAL cell differential or comorbidities at baseline (data not shown). No difference in the frequency of antifibrotic or NAC treatment was observed between genotypes.

MUC5B und TOLLIP SNPs and pulmonary dysfunction at baseline and over time

Lung function tests at baseline in patients with TOLLIP rs5743890 C/T genotype (n = 14) did not differ from those in patients with T/T genotype (n = 48) (Table S4). Patients with TOLLIP rs5743890 C/T genotype tended to have a greater decline in FVC, DLco % pred. and PaO2 mmHg per year compared to those with T/T genotype (Table S4). No differences were seen in lung function tests at baseline or decline over time according to MUC5B and TOLLIP rs3750920 genotype (data not shown).

Correlation between MUC5B and TOLLIP genotype and survival

The frequency of death or lung transplantation did not vary across MUC5B and TOLLIP genotypes. Kaplan-Meier analysis for MUC5B and TOLLIP rs3750920 did not show any correlation with survival, and median survival time did not vary according to genotype (Fig. 1a and 1c). Patients with TOLLIP rs5743890 C/T genotype had worse survival compared to T/T genotype, T being the major allele (p = 0.014) (Fig. 1b, Table S4).

Correlation between MUC5B and TOLLIP genotype and disease progression

Disease progression as defined in the methods occurred in a total of 37 (62%) of 60 patients where data were available. Median time to progression was 50 months (IQR: 34–68). No association was observed between MUC5B and TOLLIP rs3750920 genotype and disease progression.

The frequency of disease progression significantly varied across TOLLIP rs5743890 genotype: all patients with C/T genotype had disease progression compared to 50% of patients carrying T/T (p = 0.001) (Tab S5) and median time to progression tended to be shorter in patients with the C/T genotype, 43.5 ± 23 months vs 63 ± 27 months compared to those with the T/T genotype (p = 0.058). Kaplan Kaplan-Meier analysis of disease progression according to TOLLIP rs5743890 genotype in IPF patients is shown in Fig. 2a.

Correlation between MUC5B and TOLLIP SNP genotype and acute exacerbation

Acute exacerbation occurred in 19/62 (31%) patients. Median time to first acute exacerbation was 47 months (IQR: 38–73). No association between MUC5B or TOLLIP genotype and frequency of acute exacerbation was seen (Table S5). Whereas Kaplan-Meier analysis did not show any significant association between AE and TOLLIP rs5743890 genotype (Log rank p = 0.326), a clear tendency for an unfavorable effect of the minor allele can be observed (Fig. 2b).

Uni and multivariate analysis for predictors of survival, disease progression and AE

We performed univariate and multivariate analysis by using Cox regression to investigate the role of several factors as predictors for survival and disease outcome. We did not find any association between MUC5B or TOLLIP rs3750920 and survival or AE (Table 3 and S6). In the uni- and multivariate analysis for identifying predictors of death or lung transplantation, TOLLIP rs5743890 (C/T genotype) was the strongest one, followed by Dlco % pred., also after adjustment for a number of covariates (HR, 3.434, 95% CI 1.360–8.671, p = 0.009) (Table 3).

The results of uni- and multivariate analysis for predictors of disease progression are shown in Table 4. In the univariate analysis, FVC % pred. and TOLLIP rs5743890 C/T genotype were the strongest predictors of disease progression (Table 3) but in the multivariate analysis only TOLLIP rs5743890 C/T genotype was identified as a predictor of disease progression after adjustment for age, gender, BMI, smoking status and lung function at baseline as covariates (HR 2.949, 95% CI 1.327–6.556, p = 0.008) (Table 3).

In this study, we confirm that the MUC5B rs35705950 minor allele (T) is associated with IPF. The TOLLIP gene variant rs5743890 appears to correlate with survival and disease progression, therefore being of potential utility to stratify IPF patients at diagnosis.

The MUC5B promoter variant rs35705950 has been recognized as the dominant risk factor for developing (IPF) (7, 17, 20, 21). A recent metanalysis encompassing 2733 IPF patients and 5044 controls from 8 different cohorts confirmed that the minor T allele is significantly and strongly associated with an increased risk of IPF compared with the G allele (odds ratio [OR] 4.85, 95% confidence interval [CI] 3.79–6.21) in an allele dose-dependent manner. The strength of this association was more remarkable in Caucasians (22). In line with these previous studies, we found that the minor allele T (G/T and T/T genotype) is more frequent in IPF patients than in healthy controls (p < 0.0001). There was no difference with regard to age, gender or functional impairment at diagnosis between carriers or non-carriers of the minor allele T in our cohort. Borie et al. already reported similar results for French patients (23). We did not observe any association of MUC5B variants with survival, disease progression or acute exacerbation. In previous studies the associations with survival were contradictory. A significant association of the minor allele T with improved survival has been described in sporadic IPF in Caucasians (significant reduction in 2-year cumulative incidence of death)(3) whereas an association with increased mortality was reported in Han patients (24). A recent study from the Netherlands on a total of 170 IPF patients did not find an association of the minor allele T with survival in patients with sporadic IPF (N = 115), in contrast to MUC5B minor allele carriers with familial IPF (N = 55), who had a better survival (carriers 53 months vs non-carriers 37 months, p = 0.01) (25). Consistently with previous studies (3, 26), Van der Vis et al. did not find any association of MUC5B minor allele with severity of lung function impairment at diagnosis both in sporadic and familiar IPF(27). In 187 Han patients with IPF, carriers of MUC5B minor allele had significantly lower FVC and DLco values compared to those with homozygote major allele genotype (p < 0.0001 for both) (24). Further, Stock et al. did not observe any association of MUC5B minor T allele with lung function impairment at diagnosis, but with a longer time to decline in FVC (HR 0.59, p = 0.052), while no association was seen with time to decline in DLco (p = 0.9)(26). We did not find a significant association of the minor allele T with lung function impairment at diagnosis or decline over time. These findings are in line with a previous study from our group comparing the MUC5B polymorphism in German and Japanese IPF patients(18).

We did not find an association between TOLLIP SNPs alleles and IPF. This is in contrast to two previous genome-wide studies where the TOLLIP rs5743890 minor allele (C) was found to be protective and associated with reduced susceptibility to IPF(2, 20). Intriguingly, in one of these studies, IPF patients carrying this protective minor allele had increased mortality in three different cohorts (on average HR 1.72 with 95% CI 1.24–2.38, p = 0.0012). In the lung tissue from these patients the expression of TOLLIP protein was reduced by 20% (2). We also observed that patients carrying this minor allele had a reduced survival compared to those with the T/T genotype, with a remarkable difference of 20 months in median survival time. This association with survival was confirmed by multivariate analysis, after adjustment for a number of covariates. Consistently, patients with the C/T genotype had a higher rate of disease progression and tended to have a greater decline in FVC over time than patients with the T/T genotype. We could also demonstrate that the TOLLIP rs5743890 C/T genotype was the strongest predictor of disease progression in a multivariate analysis (HR 3.212, p = 0.003). So far, this is the first study showing an association of the minor allele (C) in TOLLIP rs5743890 with disease progression in IPF. Taken together, these data suggest a potential role of this SNP as a genetic biomarker for risk stratification of IPF patients in a real-world setting.

With respect to the other investigated TOLLIP SNP rs3750920, we did not find any significant association of the genotypes with survival or disease progression. However, a tendency for a protective role for AE of the minor allele (T) was observed (Fig. 3). Oldham et al. found that the TOLLIP rs3750920 polymorphisms may influence the response to NAC therapy in patients with IPF. Patients with the TOLLIP rs3750920 homozygote genotype for the minor allele may benefit from NAC therapy, whereas NAC may cause harm to patients carrying the homozygote genotype for the major allele. (4) .Due to the low number of patients treated with NAC in our cohort (24%), we were not able to replicate these observations. We did not find any association of MUC5B or TOLLIP SNPs with clinical response to treatment with antifibrotic drugs.

The importance of TOLLIP in IPF is emphasized by the relationship between biological impact and potential disease outcome. Besides the well-known role in inflammation and autoimmunity by enhancing signaling pathways associated with IL-1β, TNF-a, IFN-b, IL-13, toll-like receptor (TLR) (28) and autophagy(29), a recent study revealed that TOLLIP protects bronchial epithelial cells from bleomycin-induced apoptosis(30). A study investigating TOLLIP expression in lung tissue of IPF patients revealed that the TOLLIP expressing cells were macrophages, alveolar type II, and basal cells. A mechanistic assay using basal cells showed that TOLLIP was able to reduce mitochondrial ROS accumulation following bleomycin-induced mitochondrial damage of basal cells (30). All these data together corroborate the hypothesis of a protective effect of TOLLIP in several cell populations against oxidative damage and fibrosis. This could be a possible explanation of the link between TOLLIP SNPs and worse outcome in IPF. Further validation in bigger cohorts is needed.

Despite the novel findings of this study, it has several limitations. First, the number of patients is rather low for a genetic study and a validation cohort was not included. We found for example, that alleles in MUC5B and TOLLIP rs5743890, despite the highest linkage disequilibrium, were differently distributed, suggesting a possible genetic instability of the locus or, simpler, that a larger number of subjects is needed to confirm allele distribution. Second, there was an imbalance in age and smoking habits between patients and healthy controls. Further, only two subjects had familial IPF, so that a comparison with sporadic IPF patients, as other studies did(27), was not possible in our cohort. In conclusion, we confirm that the minor allele T at MUC5B rs35705950 is associated with IPF and that the minor allele C at TOLLIP rs5743890 is associated with an increased risk of death and disease progression. Although the results are promising, further validation in multi-center and multi-ethnic studies is needed.

AaDO₂: Alveolar arterial oxygen gradient; BAL: Bronchoalveolar lavage; CCI: Charlson comorbidity index; DLco: Diffusing capacity of the lung for carbon monoxide; FVC: Forced vital capacity; GERD: Gastroesophageal reflux disease; GWAS: genome wide association study HRCT: High resolution computed tomography; ILD: interstitial lung diseases; KM: Kaplan-Meier analysis; OSAS: PaO₂: Partial pressure of oxygen in arterial blood; ROC: SNPs: Single nucleotide polymorphisms; TLC: Total lung capacity.

Ethics approval and consent to participate.

The experiments in this study comply with the current laws in Germany. The study was approved by the local IRB (20-9236-BO). Patients and healthy subjects provided written informed consent.

Consent for publication

Not applicable.

Availability of data and materials.

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors have reported that no potential conflicts of interest exist with any companies/organizations whose products or services may be discussed in this article.

Funding

This study was supported by „Arbeitsgemeinschaft zur Förderung der Pneumologie an der Ruhrlandklinik (AFPR)“ and the German IPF patients association „Lungenfibrose eV.“

Authors' contributions

FB contributed to the conception and design of the study, collecting samples, performing genotyping, statistical analysis, drafting and finalizing the manuscript. IC and MZ performed allele, genotype and haplotype analysis, contributed to the interpretation of the results and finalization of the manuscript.

EB, SO and DT contributed to collecting samples, analysing and interpreting the data. UC and CT contributed to the conception and design of the study, data interpretation and finalization of the manuscript. All authors have read and approved the final manuscript.

Acknowledgements

No one to acknowledge.

Campo I, Zorzetto M, Bonella F. Facts and promises on lung biomarkers in interstitial lung diseases. Expert Rev Respir Med. 2015:1-21.
Noth I, Zhang Y, Ma SF, Flores C, Barber M, Huang Y, et al. Genetic variants associated with idiopathic pulmonary fibrosis susceptibility and mortality: a genome-wide association study. Lancet Respir Med. 2013;1(4):309-17.
Peljto AL, Zhang Y, Fingerlin TE, Ma SF, Garcia JG, Richards TJ, et al. Association between the MUC5B promoter polymorphism and survival in patients with idiopathic pulmonary fibrosis. JAMA. 2013;309(21):2232-9.
Oldham JM, Ma SF, Martinez FJ, Anstrom KJ, Raghu G, Schwartz DA, et al. TOLLIP, MUC5B, and the Response to N-Acetylcysteine among Individuals with Idiopathic Pulmonary Fibrosis. Am J Respir Crit Care Med. 2015;192(12):1475-82.
Thornton DJ, Rousseau K, McGuckin MA. Structure and function of the polymeric mucins in airways mucus. Annu Rev Physiol. 2008;70:459-86.
Wang H, Zhuang Y, Peng H, Cao M, Li Y, Xu Q, et al. The relationship between MUC5B promoter, TERT polymorphisms and telomere lengths with radiographic extent and survival in a Chinese IPF cohort. Sci Rep. 2019;9(1):15307.
Helling BA, Gerber AN, Kadiyala V, Sasse SK, Pedersen BS, Sparks L, et al. Regulation of MUC5B Expression in Idiopathic Pulmonary Fibrosis. Am J Respir Cell Mol Biol. 2017;57(1):91-9.
Moore C, Blumhagen RZ, Yang IV, Walts A, Powers J, Walker T, et al. Resequencing Study Confirms That Host Defense and Cell Senescence Gene Variants Contribute to the Risk of Idiopathic Pulmonary Fibrosis. Am J Respir Crit Care Med. 2019;200(2):199-208.
Zhang Q, Wang Y, Qu D, Yu J, Yang J. The Possible Pathogenesis of Idiopathic Pulmonary Fibrosis considering MUC5B. Biomed Res Int. 2019;2019:9712464.
Jezierska A, Kolosova IA, Verin AD. Toll Like Receptors Signaling Pathways as a Target for Therapeutic Interventions. Curr Signal Transduct Ther. 2011;6(3):428-40.
Raghu G, Collard HR, Egan JJ, Martinez FJ, Behr J, Brown KK, et al. An official ATS/ERS/JRS/ALAT statement: idiopathic pulmonary fibrosis: evidence-based guidelines for diagnosis and management. Am J Respir Crit Care Med. 2011;183(6):788-824.
Raghu G, Remy-Jardin M, Myers JL, Richeldi L, Ryerson CJ, Lederer DJ, et al. Diagnosis of Idiopathic Pulmonary Fibrosis. An Official ATS/ERS/JRS/ALAT Clinical Practice Guideline. Am J Respir Crit Care Med. 2018;198(5):e44-e68.
Quan H, Li B, Couris CM, Fushimi K, Graham P, Hider P, et al. Updating and validating the Charlson comorbidity index and score for risk adjustment in hospital discharge abstracts using data from 6 countries. Am J Epidemiol. 2011;173(6):676-82.
Bonella F, Ohshimo S, Miaotian C, Griese M, Guzman J, Costabel U. Serum KL-6 is a predictor of outcome in pulmonary alveolar proteinosis. Orphanet J Rare Dis. 2013;8:53.
Schmidt SL, Tayob N, Han MK, Zappala C, Kervitsky D, Murray S, et al. Predicting pulmonary fibrosis disease course from past trends in pulmonary function. Chest. 2014;145(3):579-85.
Collard HR, Ryerson CJ, Corte TJ, Jenkins G, Kondoh Y, Lederer DJ, et al. Acute Exacerbation of Idiopathic Pulmonary Fibrosis. An International Working Group Report. Am J Respir Crit Care Med. 2016;194(3):265-75.
Seibold MA, Wise AL, Speer MC, Steele MP, Brown KK, Loyd JE, et al. A common MUC5B promoter polymorphism and pulmonary fibrosis. N Engl J Med. 2011;364(16):1503-12.
Horimasu Y, Ohshimo S, Bonella F, Tanaka S, Ishikawa N, Hattori N, et al. MUC5B promoter polymorphism in Japanese patients with idiopathic pulmonary fibrosis. Respirology. 2015;20(3):439-44.
Standardization of Spirometry, 1994 Update. American Thoracic Society. Am J Respir Crit Care Med. 1995;152(3):1107-36.
Fingerlin TE, Murphy E, Zhang W, Peljto AL, Brown KK, Steele MP, et al. Genome-wide association study identifies multiple susceptibility loci for pulmonary fibrosis. Nat Genet. 2013.
Hancock LA, Hennessy CE, Solomon GM, Dobrinskikh E, Estrella A, Hara N, et al. Muc5b overexpression causes mucociliary dysfunction and enhances lung fibrosis in mice. Nat Commun. 2018;9(1):5363.
Zhu QQ, Zhang XL, Zhang SM, Tang SW, Min HY, Yi L, et al. Association Between the MUC5B Promoter Polymorphism rs35705950 and Idiopathic Pulmonary Fibrosis: A Meta-analysis and Trial Sequential Analysis in Caucasian and Asian Populations. Medicine (Baltimore). 2015;94(43):e1901.
Borie R, Crestani B, Dieude P, Nunes H, Allanore Y, Kannengiesser C, et al. The MUC5B variant is associated with idiopathic pulmonary fibrosis but not with systemic sclerosis interstitial lung disease in the European Caucasian population. PLoS One. 2013;8(8):e70621.
Jiang H, Hu Y, Shang L, Li Y, Yang L, Chen Y. Association between MUC5B polymorphism and susceptibility and severity of idiopathic pulmonary fibrosis. Int J Clin Exp Pathol. 2015;8(11):14953-8.
Juge PA, Lee JS, Ebstein E, Furukawa H, Dobrinskikh E, Gazal S, et al. MUC5B Promoter Variant and Rheumatoid Arthritis with Interstitial Lung Disease. N Engl J Med. 2018;379(23):2209-19.
Stock CJ, Sato H, Fonseca C, Banya WA, Molyneaux PL, Adamali H, et al. Mucin 5B promoter polymorphism is associated with idiopathic pulmonary fibrosis but not with development of lung fibrosis in systemic sclerosis or sarcoidosis. Thorax. 2013;68(5):436-41.
van der Vis JJ, Snetselaar R, Kazemier KM, ten Klooster L, Grutters JC, van Moorsel CH. Effect of Muc5b promoter polymorphism on disease predisposition and survival in idiopathic interstitial pneumonias. Respirology. 2016;21(4):712-7.
Brissoni B, Agostini L, Kropf M, Martinon F, Swoboda V, Lippens S, et al. Intracellular trafficking of interleukin-1 receptor I requires Tollip. Curr Biol. 2006;16(22):2265-70.
Kowalski E, Geng S, Rathes A, Lu R, Li L. Toll-interacting protein differentially modulates HIF1alpha and STAT5-mediated genes in fibroblasts. J Biol Chem. 2018;293(31):12239-47.
Li X, Kim SE, Chen TY, Wang J, Yang X, Tabib T, et al. Toll interacting protein protects bronchial epithelial cells from bleomycin-induced apoptosis. FASEB J. 2020.

TABLE 1. Demographics and characteristics of the studied subjects.

Variable	IPF (n=62)	HC (n=50)	p
Age (year)	63.5±11	42±2	<0.001
Gender (M/F), n	43/8	37/13	0.228
Smoking habits (non/ex/current), n	9/40/4	33/7/10	<0.001
Pack years, median (min-max)	28 (1-80)	15 (1-40)	0.018
BMI (Kg/m²)	27±4	25.5±3	0.064
PaO₂ (mmHg)	73±12	-	-
FVC (% pred)	69±21	-	-
DLco (% pred)	43±14	-	-
BAL cell differentials -Total cells (x104·mL−1) -Macrophages (%) -Lymphocytes (%) -Neutrophils (%) -Eosinophils (%)	18±3 69±25 10±10 15±23 3.5±3	-	-
Antifibrotic treatment (yes/no)	45/17	-	-
Duration of antifibrotic treatment (months)	43±16	-	-
Time to progression (months), median (IQR)	50 (34-68)	-	-
Time to death or lung transplant (months), median (IQR)	70 (40-95)	-	-
CCI (score)	3.7±1.7
Otherwise indicated, values are expressed as mean±SD. Charlson comorbidity index

TABLE 2. Genotype and allele distribution of the SNP rs35705950 of MUC5B, rs5743890 and rs3750920 of TOLLIP for IPF and HC subjects.

TABLE 3 Univariate and multivariate Cox proportional hazard model evaluating predictors of death or lung transplantation.

Variables	β	HR	(95% CI)	p value
Univariate analysis
Age at diagnosis (continuous)	0.030	1.030	0.976-1.088	0.282
Gender (male)	0.393	1.481	0.583-3.764	0.409
Smoking history (smoker)	0.543	1.721	0.672-4.404	0.258
CCI** (>3)	-0.178	0.837	0.280-2.498	0.749
BMI (continuous)	0.068	1.070	0.996-1.149	0.064
FVC % pred (continuous)	-0.006	0.994	0.980-1.008	0.406
PaO₂ mmHg (continuous)	0.103	1.108	0.989-1.241	0.076
DLco % pred (continuous)	0.006	0.966	0.942-0.990	0.006
Antifibrotic treatment (yes)	0.549	1.332	0.648-3.136	0.074
MUC5B (T allele)	-0.615	0.541	0.069-4.250	0.559
TOLLIP rs5743890 (C/T genotype)	0.799	2.224	1.151-4.297	0.017
TOLLIP rs3750920 (T allele)	-0.267	0.773	0.325-1.837	0.560
Multivariate analysis*
TOLLIP rs5743890 (C/T genotype)	1,234	3.434	1.360-8.671	0.009
DLco (continuous)	-0.071	0.931	0.898-0.966	0.001
* Model obtained by using backward conditional stepwise regression through 7 steps, including age, gender, smoking status, BMI, FVC, DLco and PaO2 as covariates. **Charlson Comorbidity Index

TABLE 4. Univariate and multivariate Cox proportional hazard model evaluating predictors of disease progression.

Variables	β	HR	(95% CI)	p value
Univariate analysis
Age at diagnosis (continuous)	0.025	1.026	0.991-1.062	0.152
Gender (male)	0.001	1.001	0.385-2.606	0.998
Smoking history (non smoker)	-0.663	0.515	0.060-4.437	0.546
CCI** (>3)	0.873	1.011	0.260-2.406	0.842
BMI (continuous)	0.028	1.029	0.951-1.113	0.481
FVC % pred (continuous)	0.018	1.018	1.002-1.034	0.027
PaO₂ mmHg (continuous)	0.016	1.016	0.982-1.052	0.363
DLco % pred (continuous)	0.002	1.002	0.977-1.028	0.884
Antifibrotic treatment (yes)	0.669	1.602	0.952-3.004	0.068
MUC5B (T allele)	-1.378	0.252	0.034-1.862	0.177
TOLLIP rs5743890 (C/T genotype)	0.923	3.517	1.276-5.964	0.008
TOLLIP rs3750920 (T allele)	-0.136	0.873	0.290-2.626	0.809
Multivariate analysis*
TOLLIP rs5743890 (C/T genotype)	1.082	2.949	1.327 6.556	0.008
* Model obtained by using backward conditional stepwise regression (11 steps), including age, gender, smoking status, BMI and antifibrotic treatment as covariates. **Charlson Comorbidity Index

SupplementalMaterial5OJRD.docx

Download PDF

Journal Publication

published 27 Feb, 2021

Read the published version in Orphanet Journal of Rare Diseases →

Review #1 received at journal
05 Jan, 2021
Editorial decision: Minor revision
05 Jan, 2021
Reviewer #2 agreed at journal
22 Dec, 2020
Review #2 received at journal
22 Dec, 2020
Reviewers invited by journal
20 Dec, 2020
Reviewer #1 agreed at journal
20 Dec, 2020
Editor assigned by journal
19 Dec, 2020
Submission checks completed at journal
19 Dec, 2020
Editor invited by journal
19 Dec, 2020
First submitted to journal
14 Dec, 2020

You are reading this older preprint version

Read the latest preprint version →

Potential Clinical Utility of MUC5B und TOLLIP Single Nucleotide Polymorphisms (SNPs) in the Management of Patients with IPF

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Methods

Study Subjects

Definition of Disease Progression in ILD

DNA extraction and Genotyping.

Pulmonary Function Tests and Blood Gas Analysis

Statistical Analysis

Results

Characteristics of Study Subjects

MUC5B and TOLLIP SNPs allele and genotypes distribution

Correlations between SNPs and clinical characteristics

MUC5B und TOLLIP SNPs and pulmonary dysfunction at baseline and over time

Correlation between MUC5B and TOLLIP genotype and survival

Correlation between MUC5B and TOLLIP genotype and disease progression

Correlation between MUC5B and TOLLIP SNP genotype and acute exacerbation

Uni and multivariate analysis for predictors of survival, disease progression and AE

Discussion

Abbreviations

Declarations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 1