Skip to content
naulon

No retraction, correction or concern found when last checked on Oct 18, 2026.

Genetic variants associated with susceptibility to idiopathic pulmonary fibrosis in people of European ancestry: a genome-wide association study

Richard J Allen, Joanne Porte, Rebecca Braybrooke and 40 more

The Lancet Respiratory Medicine · 2017 · CC BY

doi:10.1016/s2213-2600(17)30387-9

Abstract

Background. Idiopathic pulmonary fibrosis (IPF) is a chronic progressive lung disease with high mortality, uncertain cause, and few treatment options. Studies have identified a significant genetic risk associated with the development of IPF; however, mechanisms by which genetic risk factors promote IPF remain unclear. We aimed to identify genetic variants associated with IPF susceptibility and provide mechanistic insight using gene and protein expression analyses.

Methods. We used a two-stage approach: a genome-wide association study in patients with IPF of European ancestry recruited from nine different centres in the UK and controls selected from UK Biobank (stage 1) matched for age, sex, and smoking status; and a follow-up of associated genetic variants in independent datasets of patients with IPF and controls from two independent US samples from the Chicago consortium and the Colorado consortium (stage 2). We investigated the effect of novel signals on gene expression in large transcriptomic and genomic data resources, and examined expression using lung tissue samples from patients with IPF and controls.

Findings. 602 patients with IPF and 3366 controls were selected for stage 1. For stage 2, 2158 patients with IPF and 5195 controls were selected. We identified a novel genome-wide significant signal of association with IPF susceptibility near A-kinase anchoring protein 13 (AKAP13; rs62025270, odds ratio [OR] 1·27 [95% CI 1·18–1·37], p=1·32 × 10−9) and confirmed previously reported signals, including in mucin 5B (MUC5B; rs35705950, OR 2·89 [2·56–3·26], p=1·12 × 10−66) and desmoplakin (DSP; rs2076295, OR 1·44 [1·35–1·54], p=7·81 × 10−28). For rs62025270, the allele A associated with increased susceptibility to IPF was also associated with increased expression of AKAP13 mRNA in lung tissue from patients who had lung resection procedures (n=1111). We showed that AKAP13 is expressed in the alveolar epithelium and lymphoid follicles from patients with IPF, and AKAP13 mRNA expression was 1·42-times higher in lung tissue from patients with IPF (n=46) than that in lung tissue from controls (n=51).

Interpretation. AKAP13 is a Rho guanine nucleotide exchange factor regulating activation of RhoA, which is known to be involved in profibrotic signalling pathways. The identification of AKAP13 as a susceptibility gene for IPF increases the prospect of successfully targeting RhoA pathway inhibitors in patients with IPF.

Funding. UK Medical Research Council, National Heart, Lung, and Blood Institute of the US National Institutes of Health, Agencia Canaria de Investigación, Innovación y Sociedad de la Información, Spain, UK National Institute for Health Research, and the British Lung Foundation.

Full text

From PubMed Central, U.S. National Library of Medicine, as current on Oct 18, 2026. Paragraph numbers are ours, so an answer can point at one.

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive, fibrotic lung disease of unknown cause with a poor prognosis. The incidence of IPF in the UK is about 4·5–7·5 per 100 000 and is increasing, with a median survival of approximately 3 years., Since 2010, two new therapies, pirfenidone and nintedanib, have been approved for the treatment of IPF, but these treatments only slow disease progression and do not halt or reverse pulmonary fibrosis., Furthermore, these drugs are not universally effective and the mechanisms of their antifibrotic effects are unknown. Therefore, a detailed understanding of the genetic risk factors for IPF and their associated molecular pathways is urgently required.

The current theory suggests that IPF is characterised by initial damage to the alveolar epithelium, which then signals to various cell types, predominantly fibroblasts and macrophages, promoting tissue damage and extracellular matrix synthesis that leads to parenchymal destruction and alveolar replacement by dense fibrotic tissue. Studies from patients with familial pulmonary fibrosis have identified telomerase-related and surfactant protein-related genes that are associated with epithelial dysfunction.,, Similarly, genome-wide association studies (GWAS) of IPF have reported various independent genetic association signals related to epithelial cell function of genome-wide significance,,,, including lung defence (such as mucin 5B [MUC5B]), telomere maintenance (such as telomerase reverse transcriptase [TERT] and CST complex subunit [STN1]) and cell–cell adhesion (such as desmoplakin [DSP] and dipeptidyl peptidase 9 [DPP9]). Telomerase mutations are associated with short telomeres, which lead to increased disease progression, whereas MUC5B and Toll interacting protein (TOLLIP) variants are associated with reduced disease progression., However, these genetic abnormalities account for only about 30% of the genetic risk associated with IPF, and the molecular mechanisms that become dysregulated remain to be established.

Evidence before this studyWe searched Web of Science between Aug 21, 2015, and June 7, 2017, with the search terms “pulmonary fibrosis” and “genome wide” with no restrictions on publication date or language. Previous genome-wide association studies (GWAS) have reported 16 independent signals associated, at genome-wide significance, with susceptibility to idiopathic pulmonary fibrosis (IPF). Although previous studies have shown that epithelial cell function, lung defence, cell–cell adhesion, and telomere maintenance might play important roles in IPF, the precise mechanisms through which the genes identified in previous GWAS promote IPF are still poorly understood.Added value of this studyTo our knowledge, this study is the largest genetic study of IPF done to date, bringing together 2760 patients with IPF and 8561 controls. We present a novel genome-wide significant genetic association signal and novel gene expression data that implicate A-kinase anchoring protein 13 (AKAP13) as an IPF susceptibility gene. AKAP13 is a Rho guanine nucleotide exchange factor (RhoGEF) that regulates activation of Rho A. RhoA is a molecule with a known role in profibrotic signalling pathways; however, AKAP13 has not previously been implicated in the pathogenesis of IPF. These studies provide genetic veracity for targeting Rho signalling in IPF.Implications of all the available evidenceThe identification of an association between AKAP13, an epithelial RhoGEF, and IPF supports the role of epithelial processes in IPF pathogenesis. Furthermore, AKAP13 has a role in a pharmacologically tractable molecular pathway that could lead to new treatments for IPF.

Evidence before this study

We searched Web of Science between Aug 21, 2015, and June 7, 2017, with the search terms “pulmonary fibrosis” and “genome wide” with no restrictions on publication date or language. Previous genome-wide association studies (GWAS) have reported 16 independent signals associated, at genome-wide significance, with susceptibility to idiopathic pulmonary fibrosis (IPF). Although previous studies have shown that epithelial cell function, lung defence, cell–cell adhesion, and telomere maintenance might play important roles in IPF, the precise mechanisms through which the genes identified in previous GWAS promote IPF are still poorly understood.

Added value of this study

To our knowledge, this study is the largest genetic study of IPF done to date, bringing together 2760 patients with IPF and 8561 controls. We present a novel genome-wide significant genetic association signal and novel gene expression data that implicate A-kinase anchoring protein 13 (AKAP13) as an IPF susceptibility gene. AKAP13 is a Rho guanine nucleotide exchange factor (RhoGEF) that regulates activation of Rho A. RhoA is a molecule with a known role in profibrotic signalling pathways; however, AKAP13 has not previously been implicated in the pathogenesis of IPF. These studies provide genetic veracity for targeting Rho signalling in IPF.

Implications of all the available evidence

The identification of an association between AKAP13, an epithelial RhoGEF, and IPF supports the role of epithelial processes in IPF pathogenesis. Furthermore, AKAP13 has a role in a pharmacologically tractable molecular pathway that could lead to new treatments for IPF.

Identifying further genetic risk factors will enable the identification of molecular pathways involved in the pathogenesis of IPF that could potentially be targeted with novel treatments. Molecular targets with supporting genetic evidence are twice as likely to be successful in clinical development as those with no genetic support. We designed and implemented the largest UK-based IPF case-control GWAS to date (stage 1) and followed up promising genetic association signals in two independent IPF case-control studies (stage 2). We aimed to identify genetic variants associated with IPF susceptibility and provide mechanistic insight into those genetic association signals using gene and protein expression analyses.

Study design

This GWAS used a two-stage approach to identify novel genome-wide significant signals associated with susceptibility to IPF. In stage 1, a GWAS was done using patients with IPF from nine different centres in the UK and controls from UK Biobank. Variants showing a significant association (p<5 × 10−6) in stage 1 were then further analysed in independent samples of patients with IPF and controls from two US samples, the Colorado consortium and the Chicago consortium (stage 2). We defined statistically significant associations with susceptibility to IPF as variants that met genome-wide significance (p<5 × 10−8) after meta-analysis of stages 1 and 2 together. Variants that became less significant in the meta-analysis than in stage 1 alone were not reported as showing an association.

Stage 1 comprised patients with IPF selected from nine different centres across the UK (appendix p 1). All diagnoses were made in accordance with accepted international criteria.,, Controls were selected from UK Biobank and were matched for age, sex, and smoking distributions, showed no evidence of having any interstitial lung disease, and had genetic data available (appendix p 1). Stringent quality control testing of all samples was done, such as removing individuals with poor call rates, heterozygosity outliers, duplicates, related individuals, ancestry outliers, and sex mismatches (appendix p 1). All individuals were of European ancestry.

Stage 2 comprised two additional independent case-control studies that have been previously described by Noth and colleagues (Chicago consortium) and Fingerlin and colleagues (Colorado consortium). For both these studies, patients with IPF were diagnosed using American Thoracic Society and European Respiratory Society guidelines., All studies had appropriate institutional review board or ethics approval.

Procedures and statistical analysis

For stage 1, patients with IPF and a third of the UK Biobank controls were genotyped using the Affymetrix Axiom UK BiLEVE array (Affymetrix, Santa Clara, CA, USA). The remainder of the controls were genotyped using the Affymetrix Axiom UK Biobank array (95% identical to the UK BiLEVE array). Genotyping and imputation procedures for studies contributing to stage 1 and 2 are described in the appendix p 1. Because the UK Biobank controls genotyped on the Axiom UK BiLEVE array had originally been selected on the basis of lung function and smoking behaviour, it was reasonable to assume that allele frequency differences between the controls genotyped on the two arrays could feasibly be driven by either technical array artefacts or genuine associations with lung function and smoking. To account for this in our analysis, and to avoid incorrectly reporting associations driven by array as associations with IPF susceptibility, we did additional sensitivity analyses (appendix p 2).

The genome-wide association analysis of IPF susceptibility was done assuming an additive genetic effect and conditioning on age, sex, and the first ten principal components to adjust for ancestry. The analysis was run using a score test because of its computational efficiency, using SNPTEST version 2.5.2. For variants with minor allele count less than 400 and score test p<5 × 10−3, the analysis was rerun using the Firth test using EPACTS version 3.2.4 (appendix p 2).

Independent variants reaching a threshold of p<5 × 10−6 in association testing in stage 1 were followed up in stage 2. Conditional analyses were used to identify additional independent signals in the same genomic region (appendix p 2).

In the Chicago consortium, analyses were done using the Firth test, adjusting for age and sex. In the Colorado consortium, analyses were done using the Firth test on a logistic regression model, adjusting for sex and the top three dimensions from a multidimensional scaling model.

For variants that were significantly associated with IPF susceptibility in stage 1 only after conditioning on another variant, the analysis was run in stage 2 conditioning on the same variant. The results from these stage 2 studies were meta-analysed using a fixed-effects model.

For variants that showed a significant (p<5 × 10−6) association with IPF risk in stage 1, available association test statistics from stages 1 and 2 were meta-analysed using a fixed-effects model. Signals that had genome-wide significance (p<5 × 10−8) when meta-analysing stages 1 and 2 were reported as significantly associated with susceptibility to IPF. A Bayesian approach was used to fine-map those signals to create 95% credible sets (a set of variants that was 95% likely to contain the causal variant; appendix p 3).

To investigate whether the variants showing an association with susceptibility to IPF were also associated with survival time in IPF, a Cox proportional hazards model was fitted in a subset of the stage 1 patients with IPF that had data for survival time. The model made adjustments for age, sex, the first ten principal components, and the recruiting study centre. Analysis was done using the Survival package in R version 2.2.3.

We searched for evidence that genetic variants associated with susceptibility to IPF were independently associated with expression of a particular gene, because an association could implicate a particular gene as the driver of the signal. Variants that were significantly associated with susceptibility to IPF and proxy variants (correlated with linkage disequilibrium r2>0·8) were investigated in three expression quantitative trait locus (eQTL) datasets: a lung eQTL database comprising individuals from three cohorts (University of British Columbia, Laval, and Groningen), a blood eQTL database, and in the Genotype-Tissue Expression project (GTEx) cohort (multiple tissues; appendix p 3).

RT-PCR gene expression analysis was done on RNA extracted from human lung tissues of 46 patients with IPF and 51 controls from the Lung Tissue Research Consortium using standard methods (appendix p 3). Initially, differential gene expression in the lung was compared between the patients with IPF and controls using two-tailed Student's t test comparing the change in CT values between groups. Linear regression analysis was used to compare A-kinase anchoring protein 13 (AKAP13) expression in controls and patients with IPF while also controlling for age, smoking (ever or never and pack-years), and sex.

Formalin-fixed paraffin-embedded human lung tissue sections were obtained from tissue distant from the tumour that was obtained during lung resection from controls and tissue from patients with IPF taken at either post-mortem examination or lung transplantation. Immunohistochemistry was done using standard methods (appendix p 3).

We performed in-silico analyses to establish whether known drugs target the proteins identified in the genome-wide and expression analyses or the proteins that interact with them (appendix p 4).

Role of the funding source

The funders had no role in the study design. GlaxoSmithKline Research and Design participated in collection of data and had access to the raw data from a subset of the UK IPF patient data. LVW and RGJ were involved in all stages of study development and delivery, had full access to all data in the study, and had final responsibility for the decision to submit for publication.

Results

After sample quality control testing and genotype imputation, 602 patients with IPF and 3366 UK Biobank controls (selected on April 18, 2016) were included in the stage 1 analysis of 13 076 821 variants (figure 1A). The patients were diagnosed between June, 1996, and July, 2013, with the exception of 52 patients who were historical cases with unknown dates of diagnosis. In total, 44 independent signals were associated with susceptibility to IPF (p<5 × 10−6) and were followed up in stage 2 (figure 2).

A total of 2158 patients with IPF and 5195 controls were included in stage 2 (figure 1B, table 1). Of the 44 variants identified in stage 1, data were available for 27 variants in one or both of the stage 2 studies (appendix pp 8–9). Of these 27 variants, three were significant (p<1·85 × 10−3, after Bonferroni correction for 27 tests) in stage 2, and had genome-wide significance (p<5 × 10−8) in the meta-analysis of stage 1 and 2 together (table 2). A novel association for the single-nucleotide polymorphism (SNP) rs62025270, which was the most strongly associated SNP in a broad signal covering AKAP13 and kelch like family member 25 (KLHL25; figure 3), had genome-wide significance in the meta-analysis of stages 1 and 2 (minor allele A; odds ratio [OR] 1·27 [95% CI 1·18–1·37], p=1·32 × 10−9). This SNP was also significant in stage 2 alone after Bonferroni adjustment for 27 tests (OR 1·22 [95% CI 1·11–1·33], p=9·96 × 10−6). The 95% credible set (a set of variants that was 95% likely to contain the causal variant) produced for this signal contained 113 variants (appendix pp 10–11). The other two genome-wide significant signals were previously reported variants in MUC5B (rs35705950) and DSP (rs2076295).,,

We explored whether genotyping array differences could have influenced the association test result for SNP rs62025270. In the control–control comparison (appendix p 6), rs62025270 showed no association with the array (p=0·35, appendix p 12). We repeated the association analysis for rs62025270 restricted to controls who were genotyped on the same array as the patients (602 patients vs 1231 controls); the results were consistent with the original analysis (stage 1 OR 1·45, [95% CI 1·20–1·74], p=1·19 × 10−4; meta-analysis OR 1·25 [95% CI 1·17–1·33], p=1·86 × 10−8). We reimputed the region with only variants genotyped on both arrays using the new Haplotype Reference Consortium imputation reference panel; the results were again consistent with the original analysis (stage 1 OR 1·48 [95% CI 1·25–1·75], p=4·24 × 10−6; meta-analysis OR 1·27 [95% CI 1·19–1·35], p=1·55 × 10−9).

Previous GWAS,,, reported 16 independent signals associated with susceptibility to IPF with genome-wide significance. In stage 1 of our study, which analysed previously unreported data only, nine of these signals had a p value less than 0·05 and a direction of effect consistent with the previously reported effects. Furthermore, seven signals met a Bonferroni adjusted significance threshold for 16 tests (p<3·13 × 10−3) with a consistent direction of effect to previous GWAS; namely variants in MUC5B, DSP, family with sequence similarity 13 member A (FAM13A), TERT, mucin 2 (MUC2), DPP9, and TOLLIP (appendix p 13). Five of the remaining seven signals that were not significant had a consistent direction of effect with previous GWAS.

SNP rs62025270[A] was associated with increased expression of AKAP13 in the lung. To identify the functional consequence of the rs62025270 polymorphism, we assessed associations between rs62025270 and nearby variants and expression of genes (eQTLs). In non-diseased whole-lung tissue from the lung eQTL databases (University of British Columbia n=339, Laval n=409, and Groningen n=363), the minor allele A of rs62025270, associated with increased susceptibility to IPF, was associated with increased expression of AKAP13 (figure 4). SNP rs62025270 is in linkage disequilibrium (r2=0·60) with rs17636666, the variant that was most significantly associated with expression of AKAP13 in lung tissue. This co-localisation of IPF susceptibility association signal and strongest AKAP13 gene expression association signal (eQTL) in lung tissue (figure 3) further suggests that altered expression of AKAP13 has a role in IPF susceptibility. The minor allele A of rs62025270 (via proxy SNP rs2554, r2=0·93) was also associated with expression of AKAP13 in whole blood (n=5311; appendix pp 7 and 14), although it was not in linkage disequilibrium (r2=0·03) with the SNP most strongly associated with AKAP13 expression in blood (rs870689). This finding suggests differences in regulation of AKAP13 expression between lung and whole blood, which have different cell-type compositions. Furthermore, the minor allele A of rs62025270 was associated with decreased expression of AKAP13 in blood, rather than increased expression as seen in lung tissue. The SNP rs62025270 was also found to be associated with expression of RNA genes RP11–158M2.3, RP11–158M2.4, RP11–158M2.5, RP11–815J21.3, and RP11–815J21.4 in a range of tissues in GTEx (appendix p 14).

Having identified that the minor allele A of rs62025270 was associated with increased susceptibility to IPF and increased expression of AKAP13 in lung tissue, we did a further analysis in lung tissue to identify the cellular location of AKAP13 in the lung and compared levels of AKAP13 mRNA from patients with IPF versus controls. Morphological assessment of histological sections from 10 control and 10 patients with IPF showed that AKAP13 protein was expressed primarily in bronchial epithelium (figure 5A) and alveolar type 1 and 2 cells (figure 5B) in control lung samples. High AKAP13 expression was observed in fibrotic regions of lung tissue from patients with IPF, and the cells expressing AKAP13 were primarily epithelial cells, macrophages, and lymphoid aggregates (figure 5C–F). High AKAP13 expression was also seen in epithelial cells from less fibrotic alveolar regions of lungs of patients with IPF, with more AKAP13-expressing cells in alveoli of patients with IPF (figure 5G) than in alveoli of controls (figure 5H). Real-time PCR analysis of samples from 46 patients with IPF and 51 controls established that concentrations of AKAP13 mRNA were 1·42-times higher in whole lung tissue homogenates from patients with IPF than in lung tissue from controls (p=0·0011; figure 5I). Linear regression analysis confirmed that the significant increase in AKAP13 mRNA expression in lung tissue of patients with IPF was maintained after controlling for age, sex, and smoking history (p=2·03 × 10−4).

We tested the association of rs62025270 with survival time in all patients from stage 1 for whom survival data were available (n=565; 360 with recorded deaths and 205 with censored survival times, with a minimum follow-up time of 1112 days before censoring). The median survival time after diagnosis of IPF was 1621 days (4·4 years; IQR 764–2815 days). rs62025270 showed no association with survival time (hazard ratio 1·01 [95% CI 0·86–1·19], p=0·878).

11 proteins that interact with AKAP13 were identified as known targets for existing drugs or compounds in development, including aspirin, dextromethophan, and GSK-690693 (appendix pp 15–17).

Discussion

We studied 2760 patients with IPF and 8561 controls and identified a novel genome-wide significant association of variant rs62025270 that implicates AKAP13 as a potential driver of IPF pathogenesis. We showed that the minor allele (A) of rs62025270 (minor allele frequency [MAF]=25%), associated with increased susceptibility to IPF, was associated with increased expression of AKAP13 in lung tissue. Morphological analysis of control and diseased lung tissue identified AKAP13 expression primarily in epithelial cells, with some additional expression in lymphoid tissue in patients with IPF. Furthermore, we showed increased expression of AKAP13 in fibrotic lung tissue from patients with IPF compared with samples from healthy controls.

The high AKAP13 expression in epithelial cells further supports the hypothesis that IPF is a disease characterised by epithelial susceptibility to injury, consistent with findings of other genetic risk factors associated with epithelial cell function. AKAP13 is a RhoA guanine nucleotide exchange factor (RhoGEF) that is crucial for murine cardiac development. AKAP13-null mice are embryonic lethal because of arrested cardiac development at embryonic day 10. Conversely, enhanced AKAP13 signalling has been shown to promote profibrotic signals in cardiac fibroblasts and leads to cardiac hypertrophy in mice. Although AKAP13 has not previously been implicated in the pathogenesis of IPF, it plays a key role in profibrotic signalling pathways. Thrombin and lysophosphatidic acid have been identified as key profibrotic mediators acting via a Gq and RhoA signalling pathway in epithelial cells to promote αvβ6 integrin-mediated transforming growth factor β (TGFβ) activation., AKAP13 coordinates thrombin and lysophosphatidic acid receptor signalling through Ga12 to RhoA in fibroblasts and astrocytes, and has also been shown to functionally interact with Gq, promoting RhoA activation.,, Therefore, AKAP13 might promote epithelial αvβ6 integrin-mediated TGFβ activation in response to epithelial injury, thereby promoting IPF. A phase 2 clinical trial assessing the safety and tolerability of a humanised monoclonal antibody directed against the αvβ6 integrin heterodimer (compound BG00011, formerly known as STX-100) in participants with IPF is underway (NCT01371305).

We also identified AKAP13 in lymphoid follicles within fibrotic regions of lung tissue from patients with IPF. Lymphoid follicles without germinal centres are well described in IPF and are thought to be composed primarily of B cells and dendritic cells.,,,, Increased concentrations of plasma B lymphocyte stimulating factor have been reported in patients with IPF compared with either controls or patients with chronic obstructive pulmonary disease, and higher concentrations were associated with severe disease. Furthermore, analysis of gene expression profiles from lymphoid follicles identified them as the source of C-X-C motif chemokine ligand 13, a biomarker of severe IPF,, suggesting that lymphoid follicles might play an important undefined role in the pathogenesis of IPF. The role of AKAP13 in B-cell function is similarly poorly understood, although it has been shown to be responsible for glucocorticoid responsiveness in lymphocytes stimulated with lysophosphatidic acid, and is involved in lymphocyte responses to osmotic stress.

Because of the known function of AKAP13, it probably has some potential downstream molecular interactions that could be targeted therapeutically, most notably RhoA and Rho kinase. 11 proteins that interact with AKAP13 were identified as known targets for existing drugs or compounds in development (appendix pp 15–17). One of the existing drugs was aspirin, which targets prostaglandin-endoperoxide synthase 1 (PTGS1) and PTGS2, and PTGS2 has a well described role in the pathogenesis of IPF.,,, Another of the drugs was dextromethorphan, a well known cough-suppressant that antagonises N-methyl-D-aspartate receptors. Notably, a novel pan AKT inhibitor, GSK-690693, was also identified as a possible asset for targeting AKAP13 binding partners, and the AKT pathway is emerging as an important pathway in fibrogenesis.

To our knowledge this is the first GWAS implicating AKAP13 at genome-wide significance in IPF susceptibility. Additionally, we confirmed seven previously reported associations with IPF susceptibility: the T allele of SNP rs35705950 in the MUC5B promoter (stage 1 p=1·86 × 10−37) was the strongest genetic risk factor for susceptibility to IPF with an OR of 4·11, then DSP (stage 1 p=4·14 × 10−12), MUC2 (stage 1 p=1·33 × 10−5), TOLLIP (stage 1 p=1·60 × 10−5), TERT (stage 1 p=8·25 × 10−5), DPP9 (stage 1 p=4·17 × 10−4), and FAM13A (stage 1 p=0·002). 17 variants were identified in stage 1 that have not been previously reported to be associated with IPF (of which 14 were low frequency or rare), but we could not investigate all of these in the stage 2 analyses because not all data were available. These could be additional true positive signals of association with susceptibility to IPF that will reach genome-wide significance in larger and more densely imputed studies.

A strength of our study is that we were able to bring together the largest sample size of patients with IPF and controls available to date and apply a robust and widely used two-stage study design,, to identify a novel genetic signal of association with IPF susceptibility. Furthermore, our stage 1 study samples were genotyped using the Affymetrix Axiom UK BiLEVE and UK Biobank arrays, which are optimised for imputation of individuals with European ancestry, and were genome-wide imputed using the most comprehensive imputation panel resource available at the time. This enabled analysis of 13 076 821 variants genome-wide with MAF of more than 0·1%.

This study also has some limitations. Although clinical investigators were careful to exclude other fibrotic lung diseases, a small number of non-IPF cases of pulmonary fibrosis might have been included. However, the advantage of collating a large dataset is that the effects of misclassification on the analysis are reduced. Furthermore, we expect that misclassification would lead to attenuation of signals rather than produce false-positive findings, which was supported by our replication of nine previously reported signals of association with IPF in our previously unreported stage 1 data. This included five signals that reached a Bonferroni-adjusted significance threshold for 16 tests (p<3·13 × 10−3), and two additional signals that reached genome-wide significance (MUC5B p=1·86 × 10−37 and DSP p=4·14 × 10−12). In addition, the diagnostic practices between the UK and USA might be different but the concordant results from the three case-control studies reassures us that these observations are robust.

A previous study of 119 patients with IPF and 50 donor lung controls from the Lung Tissue Resource Consortium reported a small decrease of AKAP13 expression in patients with IPF compared with controls using gene expression data from the Affymetrix Gene ST1·0 array; whereas our data using RT-PCR showed a moderate increase of AKAP13 expression in patients with IPF compared with controls. The discordance between these results could be due to a number of factors, including the difference in normalisation strategies between the two techniques and sample heterogeneity. The eQTL data that we present from control lung tissue shows that the minor allele A of rs62025270, associated with susceptibility to IPF, is also associated with increased expression of AKAP13 in three independent studies.,, The immunohistochemistry results provide evidence that the predominant cell types expressing AKAP13, in both normal and fibrotic lung, are epithelial cells and resident immune cells, whereas the dense fibrotic regions of lung with an expanded population of myofibroblasts (figure 5C–F) do not express AKAP13. Furthermore, eQTL data from whole blood reveals the minor allele A of rs62025270 is associated with reduced expression of AKAP13 mRNA. Therefore, these data provide strong independent evidence that increased AKAP13 expression in alveolar cells has a role in IPF, and the cellular heterogeneity within the fibrotic lung samples used for RNA analysis is likely to explain this discordance.

In summary, we report a novel genome-wide significant association for IPF susceptibility with SNP rs62025270 and present evidence that this SNP might exert its effect via expression of the nearby gene AKAP13. AKAP13 is a RhoGEF that is known to interact with a central fibrogenic pathway—G protein-coupled receptor activation of RhoA—and expression is increased in fibrotic regions of lungs from patients with IPF. Expression of AKAP13 occurs primarily in the epithelium thereby reinforcing the idea that epithelial processes are central to IPF pathogenesis. The identification of a RhoGEF that regulates a pathway known to be involved in the pathogenesis of IPF, a pathway for which drugs are in development, could potentially provide a novel target for antifibrotic therapy for IPF.

References 55

  1. 1.

    Navaratnam V, Fleming KM, West J. The rising incidence of idiopathic pulmonary fibrosis in the UK. Thorax. 2011;66:462-467.

    Not in our library10.1136/thx.2010.148031

  2. 2.

    Coward WR, Saini G, Jenkins G. The pathogenesis of idiopathic pulmonary fibrosis. Ther Adv Respir Dis. 2010;4:367-388.

    Not in our library10.1177/1753465810379801

  3. 3.

    Ley B, Collard HR, King TE. Clinical course and prediction of survival in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2011;183:431-440.

    Not in our library10.1164/rccm.201006-0894ci

  4. 4.

    King TE, Bradford WZ, Castro-Bernardini S. A phase 3 trial of pirfenidone in patients with idiopathic pulmonary fibrosis. N Engl J Med. 2014;370:2083-2092.

    Not in our library10.1056/nejmoa1402582

  5. 5.

    Richeldi L, du Bois RM, Raghu G. Efficacy and safety of nintedanib in idiopathic pulmonary fibrosis. N Engl J Med. 2014;370:2071-2082.

    Not in our library10.1056/nejmoa1402584

  6. 6.

    Thomas AQ, Lane K, Phillips J. Heterozygosity for a surfactant protein C gene mutation associated with usual interstitial pneumonitis and cellular nonspecific interstitial pneumonitis in one kindred. Am J Respir Crit Care Med. 2002;165:1322-1328.

    Not in our library10.1164/rccm.200112-123oc

  7. 7.

    van Moorsel CH, van Oosterhout MF, Barlo NP. Surfactant protein C mutations are the basis of a significant portion of adult familial pulmonary fibrosis in a Dutch cohort. Am J Respir Crit Care Med. 2010;182:1419-1425.

    Not in our library10.1164/rccm.200906-0953oc

  8. 8.

    Stuart BD, Choi J, Zaidi S. Exome sequencing links mutations in PARN and RTEL1 with familial pulmonary fibrosis and telomere shortening. Nat Genet. 2015;47:512-517.

    Not in our library10.1038/ng.3278

  9. 9.

    Mushiroda T, Wattanapokayakit S, Takahashi A. A genome-wide association study identifies an association of a common variant in TERT with susceptibility to idiopathic pulmonary fibrosis. J Med Genet. 2008;45:654-656.

    Not in our library10.1136/jmg.2008.057356

  10. 10.

    Noth I, Zhang Y, Ma SF. Genetic variants associated with idiopathic pulmonary fibrosis susceptibility and mortality: a genome-wide association study. Lancet Respir Med. 2013;1:309-317.

    Not in our library10.1016/s2213-2600(13)70045-6

  11. 11.

    Fingerlin TE, Murphy E, Zhang W. Genome-wide association study identifies multiple susceptibility loci for pulmonary fibrosis. Nat Genet. 2013;45:613-620.

    Not in our library10.1038/ng.2609

  12. 12.

    Fingerlin TE, Zhang W, Yang IV. Genome-wide imputation study identifies novel HLA locus for pulmonary fibrosis and potential role for auto-immunity in fibrotic idiopathic interstitial pneumonia. BMC Genet. 2016;17:74.

    Not in our library10.1186/s12863-016-0377-2

  13. 13.

    Mathai SK, Newton CA, Schwartz DA, Garcia CK. Pulmonary fibrosis in the era of stratified medicine. Thorax. 2016;71:1154-1160.

    Not in our library10.1136/thoraxjnl-2016-209172

  14. 14.

    Stuart BD, Lee JS, Kozlitina J. Effect of telomere length on survival in patients with idiopathic pulmonary fibrosis: an observational cohort study with independent validation. Lancet Respir Med. 2014;2:557-565.

    Not in our library10.1016/s2213-2600(14)70124-9

  15. 15.

    Peljto AL, Zhang Y, Fingerlin TE. Association between the MUC5B promoter polymorphism and survival in patients with idiopathic pulmonary fibrosis. JAMA. 2013;309:2232-2239.

    Not in our library10.1001/jama.2013.5827

  16. 16.

    Nelson MR, Tipney H, Painter JL. The support of human genetic evidence for approved drug indications. Nat Genet. 2015;47:856-860.

    Not in our library10.1038/ng.3314

  17. 17.

    Travis WD, King TE, Bateman ED. American Thoracic Society/European Respiratory Society international multidisciplinary consensus classification of the idiopathic interstitial pneumonias. Am J Respir Crit Care Med. 2002;165:277-304.

    Not in our library10.1164/ajrccm.165.2.ats01

  18. 18.

    Raghu G, Collard HR, Egan JJ. 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:788-824.

    Not in our library10.1164/rccm.2009-040gl

  19. 19.

    Raghu G, Rochwerg B, Zhang Y. An official ATS/ERS/JRS/ALAT clinical practice guideline: treatment of idiopathic pulmonary fibrosis, an update of the 2011 clinical practice guideline. Am J Respir Crit Care Med. 2015;192:e3-e19.

    Not in our library10.1164/rccm.201506-1063st

  20. 20.
  21. 21.

    Marchini J, Howie B, Myers S, McVean G, Donnelly P. A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007;39:906-913.

    Not in our library10.1038/ng2088

  22. 22.

    Kang H. Efficient and parallelisable association container toolbox (EPACTS). 2014.

    Not in our library

  23. 23.

    McCarthy S, Das S, Kretzschmar W. A reference panel of 64,976 haplotypes for genotype imputation. Nat Genet. 2016;48:1279-1283.

    Not in our library10.1038/ng.3643

  24. 24.

    Mayers CM, Wadell J, McLean K. The rho guanine nucleotide exchange factor AKAP13 (BRX) is essential for cardiac development in mice. J Biol Chem. 2010;285:12344-12354.

    Not in our library10.1074/jbc.m110.106856

  25. 25.

    Cavin S, Maric D, Diviani D. A-kinase anchoring protein-Lbc promotes pro-fibrotic signaling in cardiac fibroblasts. Biochim Biophys Acta. 2014;1843:335-345.

    Not in our library10.1016/j.bbamcr.2013.11.008

  26. 26.

    del Vescovo CD, Cotecchia S, Diviani D. A-kinase-anchoring protein-Lbc anchors IκB kinase β to support interleukin-6-mediated cardiomyocyte hypertrophy. Mol Cell Biol. 2013;33:14-27.

    Not in our library10.1128/mcb.00887-12

  27. 27.

    Jenkins RG, Su X, Su G. Ligation of protease-activated receptor 1 enhances alpha(v)beta6 integrin-dependent TGF-beta activation and promotes acute lung injury. J Clin Invest. 2006;116:1606-1614.

    Not in our library10.1172/jci27183

  28. 28.

    Xu MY, Porte J, Knox AJ. Lysophosphatidic acid induces alphavbeta6 integrin-mediated TGF-beta activation via the LPA2 receptor and the small G protein G alpha(q). Am J Pathol. 2009;174:1264-1279.

    Not in our library10.2353/ajpath.2009.080160

  29. 29.

    Diviani D, Soderling J, Scott JD. AKAP-Lbc anchors protein kinase A and nucleates Galpha 12-selective Rho-mediated stress fiber formation. J Biol Chem. 2001;276:44247-44257.

    Not in our library10.1074/jbc.m106629200

  30. 30.

    Majumdar M, Seasholtz TM, Buckmaster C, Toksoz D, Brown JH. A rho exchange factor mediates thrombin and Galpha(12)-induced cytoskeletal responses. J Biol Chem. 1999;274:26815-26821.

    Not in our library10.1074/jbc.274.38.26815

  31. 31.

    Pi M, Spurney RF, Tu Q, Hinson T, Quarles LD. Calcium-sensing receptor activation of rho involves filamin and rho-guanine nucleotide exchange factor. Endocrinology. 2002;143:3830-3838.

    Not in our library10.1210/en.2002-220240

  32. 32.

    Dutt P, Kjoller L, Giel M, Hall A, Toksoz D. Activated Galphaq family members induce Rho GTPase activation and Rho-dependent actin filament assembly. FEBS Lett. 2002;531:565-569.

    Not in our library10.1016/s0014-5793(02)03625-6

  33. 33.

    Bear MD, Li M, Liu Y, Giel-Moloney MA, Fanburg BL, Toksoz D. The lbc rho guanine nucleotide exchange factor alpha-catulin axis functions in serotonin-induced vascular smooth muscle cell mitogenesis and RhoA/ROCK activation. J Biol Chem. 2010;285:32919-32926.

    Not in our library10.1074/jbc.m109.062513

  34. 34.

    Campbell DA, Poulter LW, Janossy G, du Bois RM. Immunohistological analysis of lung tissue from patients with cryptogenic fibrosing alveolitis suggesting local expression of immune hypersensitivity. Thorax. 1985;40:405-411.

    Not in our library10.1136/thx.40.6.405

  35. 35.

    Wallace WA, Howie SE, Krajewski AS, Lamb D. The immunological architecture of B-lymphocyte aggregates in cryptogenic fibrosing alveolitis. J Pathol. 1996;178:323-329.

    Not in our library10.1002/(sici)1096-9896(199603)178:3<323::aid-path467>3.0.co;2-7

  36. 36.

    Marchal-Somme J, Uzunhan Y, Marchand-Adam S. Cutting edge: nonproliferating mature immune cells form a novel type of organized lymphoid structure in idiopathic pulmonary fibrosis. J Immunol. 2006;176:5735-5739.

    Not in our library10.4049/jimmunol.176.10.5735

  37. 37.

    Marchal-Sommé J, Uzunhan Y, Marchand-Adam S. Dendritic cells accumulate in human fibrotic interstitial lung disease. Am J Respir Crit Care Med. 2007;176:1007-1014.

    Not in our library10.1164/rccm.200609-1347oc

  38. 38.

    DePianto DJ, Chandriani S, Abbas AR. Heterogeneous gene expression signatures correspond to distinct lung pathologies and biomarkers of disease severity in idiopathic pulmonary fibrosis. Thorax. 2015;70:48-56.

    Not in our library10.1136/thoraxjnl-2013-204596

  39. 39.

    Xue J, Kass DJ, Bon J. Plasma B lymphocyte stimulator and B cell differentiation in idiopathic pulmonary fibrosis patients. J Immunol. 2013;191:2089-2095.

    Not in our library10.4049/jimmunol.1203476

  40. 40.

    Vuga LJ, Tedrow JR, Pandit KV. CXC motif chemokine 13 (CXCL13) is a prognostic biomarker of idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2014;189:966-974.

    Not in our library10.1164/rccm.201309-1592oc

All 55 references

Cited by papers we hold 11

Reads

Read by people on this page

None yet.

Checked or read by AI through naulon

None yet.

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.

Are you an author?

Is this your paper?

Claim it to add your note and be paid when agents cite it. Free.