Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A family-based genome-wide association study of chronic rhinosinusitis with nasal polyps implicates several genes in the disease pathogenesis

  • Anton Bohman ,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing

    anton.bohman@gmail.com

    Affiliation Department of Otorhinolaryngology, Uppsala University Hospital, Uppsala, Sweden

  • Julius Juodakis,

    Roles Data curation, Formal analysis, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Obstetrics and Gynecology, Institute of Clinical Sciences, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

  • Martin Oscarsson,

    Roles Conceptualization, Investigation, Resources, Writing – review & editing

    Affiliation Department of Otorhinolaryngology, Skaraborg Hospital, Skövde, Sweden

  • Jonas Bacelis,

    Roles Methodology, Supervision, Writing – review & editing

    Affiliation Department of Obstetrics and Gynecology, Institute of Clinical Sciences, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

  • Mats Bende,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Validation, Writing – review & editing

    Affiliation Department of Otorhinolaryngology, Skaraborg Hospital, Skövde, Sweden

  • Åsa Torinsson Naluai

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

    Affiliation Department of Microbiology and Immunology, Institute of Biomedicine, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

Abstract

Background

The pathogenesis of chronic rhinosinusitis with nasal polyps is largely unknown. Previous studies have given valuable information about genetic variants associated with this disease but much is still unexplained. Our goal was to identify genetic markers and genes associated with susceptibility to chronic rhinosinusitis with nasal polyps using a family-based genome-wide association study.

Methods

427 patients (293 males and 134 females) with CRSwNP and 393 controls (175 males and 218 females) were recruited from several Swedish hospitals. SNP association values were generated using DFAM (implemented in PLINK) and Efficient Mixed Model Association eXpedited (EMMAX). Analyses of pathway enrichment, gene expression levels and expression quantitative trait loci were then performed in turn.

Results

None of the analysed SNPs reached genome wide significant association of 5.0 x 10−8. Pathway analyses using our top 1000 markers with the most significant association p-values resulted in 138 target genes. A comparison between our target genes and gene expression data from the NCBI Gene Expression Omnibus database showed significant overlap for 36 of these genes. Comparisons with data from expression quantitative trait loci showed the most skewed allelic distributions in cases with chronic rhinosinusitis with nasal polyps compared with controls for the genes HLCS, HLA-DRA, BICD2, VSIR and SLC5A1.

Conclusion

Our study indicates that HLCS, HLA-DRA, BICD2, VSIR and SLC5A1 could be involved in the pathogenesis of chronic rhinosinusitis with nasal polyps. HLA-DRA has been associated with chronic rhinosinusitis with nasal polyps in previous studies and HLCS, BICD2, VSIR and SLC5A1 may be new targets for future research.

Introduction

Chronic rhinosinusitis (CRS) as defined by the European Position Paper on Rhinosinusitis and Nasal Polyps (EPOS) [1] is classified into chronic rhinosinusitis with nasal polyps (CRSwNP) and chronic rhinosinusitis without nasal polyps (CRSsNP). CRSwNP is a disease characterized by benign outgrowths from the middle meatus of the nasal cavity and chronic sinonasal inflammation. CRSwNP is a common chronic disease and depending on the geographical area, 2–4% of the population is afflicted [24]. The disease causes individual suffering and a decreased quality of life [5,6]. Risk factors include asthma, male sex and increasing age. The disease often requires a combination of surgical and medical treatment. However, CRSwNP often recurs even after therapy.

The aetiology of the disease is unknown. Several environmental factors have been suggested and previous studies have also shown an increased prevalence among relatives [7,8] and a higher rate of positive family history of CRSwNP among those affected [911], confirming a genetic susceptibility to the disease. Compared to the general population, having an afflicted family member increases the risk of disease five times [7]. Genetic studies on patients with CRSwNP could help to explain the pathogenesis of the disease and over time identify new drug targets leading to a more effective, individually tailored, therapy.

Genetic association can be explored using candidate gene studies or genome-wide association studies (GWAS). Candidate gene studies usually investigate a small number of single-nucleotide polymorphisms (SNPs) or other types of genetic variation, in order to find or reject associations between the genetic variants and the disease in question. These studies rely on previous knowledge and hypotheses regarding which SNPs to suspect and investigate. In comparison, a GWAS investigates hundreds of thousands of SNPs across the whole genome and is therefore not reliant on previous knowledge or hypotheses regarding the pathogenesis of the investigated disease or trait. A large number of GWASs have been performed for various complex diseases such as diabetes and asthma which has led to the finding of novel genetic pathways [12]. There is currently no published GWAS performed only on patients with CRSwNP but there is a pooling-based GWAS done on patients with CRS (both CRSsNP and CRSwNP) [13] as well as several studies of candidate genes [14]. These studies have implicated several genes and pathways such as the cystic fibrosis transmembrane conductance regulator gene (CFTR) [15,16] and, among others, genes involved in immunity [1720], inflammation [21,22], tissue remodelling [20] and arachidonic acid metabolism [23]

Even though GWA studies have been successful in detecting genetic variants associated with many common diseases, the inability to explain most of the estimated heritability makes linkage analysis an alternative to detect possible rare variants. To this date, no published linkage analysis have been performed on subjects with CRSwNP. However, one such study has been performed on 8 subjects with CRS (not specified whether any of them had CRSwNP) which found a linkage signal on chromosome 7q [24]. A combination of a GWAS and linkage analysis such as a family-based GWAS could be a more potent way of identifying both common and rare variants associated with CRSwNP [25].

Aim

The aim of this study is to identify SNPs and genes associated with CRSwNP susceptibility using a family-based approach.

Materials and methods

Material

367 patients with CRSwNP (250 men and 117 women, mean age 52.3 years) from three Swedish ear, nose and throat clinics were recruited. These subjects were all known patients at their respective clinics, all of them fulfilled the EPOS criteria for CRSwNP [1] and had at least intermittently been on either intranasal or systemic corticosteroids, most of them had undergone at least one operation for the condition. In order to increase power level for genome-wide analysis, patients with associated diseases such as asthma or aspirin-exacerbated respiratory disease (AERD) were not excluded. A total of 453 first-degree relatives (218 men, 235 women, mean age 49.4 years) were also recruited.

The study was carried out in accordance with the Declaration of Helsinki and was approved by the Ethics Committee at the University of Gothenburg, Sweden. Written consent was obtained from all participants.

Investigation

Nasal endoscopy was performed on all participants using a 2.7 mm rigid endoscope (KARL STORZ) and the participants were subsequently phenotyped as either having CRSwNP or being free from this disease. Additional data about asthma and corticosteroid medication (used to counter symptoms from either the upper or lower airways) was obtained via a structured interview. Peripheral blood samples were collected from each individual.

DNA was extracted from whole blood using an in-house protocol at KBiosciences (LGC Genomics, Hoddeston UK). The HiSeq Illumina platform was used for genotyping. 144 of the samples were run on Illumina Omni Express bead chips and the remaining 676 on the Illumina Core Exome array.

Quality control

Individuals were removed if they showed >2% missing calls (all samples passed), heterozygosity >3 SDs above or below mean, >100 Mendelian errors). Markers with >2% missing calls, >3 Mendelian errors or those with minor allele frequency = 0 (i.e. monomorphic loci) were excluded.

SNPs with low linkage disequilibrium among each other (LD) (r2<0.2) were selected for population structure analysis. As a reference, samples from the 1000 Genomes project, Phase 3 were retrieved. Principal component analysis was then performed using GCTA, a genome-wide complex trait analysis software [26], and the first three principal components were investigated. As a measure of non-European admixture in each sample, we calculated the Euclidean distance E from that sample to the mass centre of the CEU reference group. Individuals with E > 5 SDE were removed (2 samples).

Only autosomal markers shared in both genotyping platforms were retained. Finally, principal components analysis was performed to check for batch effects. Visual inspection of sample scores along the first three principal components showed no differences between batches.

Association testing

Two methods were used to generate SNP association values. First, we used DFAM, implemented in PLINK, which combines the transmission disequilibrium test (TDT), the sibling TDT and an allelic test for unrelated cases and controls in a single Cochran-Mantel-Haenszel test for each marker [27]. The second method was EMMAX (Efficient Mixed Model Association eXpedited), implemented in Golden Helix SNP & Variation Suite v8.3.4 [28]. This method involves computing an empirical relatedness matrix of the samples, and using this relatedness as a covariate in linear regression for each marker [29]. We performed this test using additive, dominant and recessive models, and the smallest p-value from the three models was assigned to each SNP. A commonly used level of significance in conventional GWA-studies on unrelated subjects is 5x10-8 [30] but many variants that do not reach this level of significance can still be true disease-influencing variation. The decision was therefore made to perform post-GWAS analyses for the 1000 SNPs with the lowest p-values for each method rather than adhering to a strict threshold for the level of significance.

Pathway enrichment analysis

The top 1000 markers with the most significant association p-values were combined into intervals of SNPs in high LD (defined by pairwise r2>0.25). This process was performed separately for DFAM and EMMAX association results.

INRICH software was then used to detect possible pathway enrichment within these regions [31]. To calculate the significance of such overlaps, the process was repeated 50,000 times with random genomic regions, however matched in size and SNP density. INRICH analysis was performed separately using DFAM or EMMAX results and Gene Ontology (GO) or Kyoto Encyclopedia of Genes and Genomes (KEGG)-based gene-sets. The 20 gene-sets with the highest enrichment p-values were retrieved from each of these setups. INRICH produced a list of genes which were located close to the top GWAS ‘hits’ in the genome and that share functional annotations. All genes retrieved in this way from the four INRICH analyses (DFAM+GO, EMMAX+GO, DFAM+KEGG, EMMAX+KEGG) were combined together, creating a list of target NP genes implicated in this study.

Gene expression data

Publicly-available gene expression data, collected by Plager et al. [32], was retrieved from NCBI Gene Expression Omnibus (GEO) database ([33]; series accession number GSE23552). Per authors’ recommendation, two samples (aCRSm1 and aCRSm2) were excluded, leaving 20 case samples (all from patients with CRSwNP) and 17 control samples from either allergic rhinitis patients or healthy individuals. Expression levels between the case and control groups were compared using the GEO2R interface. All genes corresponding to probes with significant difference in expression levels (Benjamini-Hochberg FDR < 0.05) comprise the differentially-expressed gene set.

Expression quantitative trait loci (eQTL) analysis

Two datasets were used for eQTL analysis: Blood eQTL from Westra et al. [34] and Multiple Tissue Human Expression Resource (MuTHER) project [35]. These datasets are produced by microarray genotyping and expression profiling of selected tissue samples. SNP variations are then associated with gene expression patterns, resulting in a list of regulatory SNPs for each gene.

In MuTHER project, the regulatory effects of each SNP were determined in adipose, skin tissues and lymphoblastoid cell lines (LCL). For each SNP-gene pair we have retained either LCL or skin data, corresponding to the tissue with more significant regulatory effect (i.e. lower p-value in skin samples meant that LCL data was discarded for that SNP-gene pair). We also excluded all SNPs with p-values > 0.05 or absolute effect size (regression coefficient β) of < 0.01. FDR of 0.5 was used as a cut-off for the Blood dataset, with no additional limits on effect size.

To check for directed eQTL enrichment, all eQTL SNPs for each gene of interest were extracted and classified according to the direction of their regulatory effect (up-regulating or down-regulating). The frequency of the allele bearing the reported regulatory effect was then determined in our GWAS cases and controls using PLINK [27]. The marker was then assigned to a bin depending on whether the regulatory allele shows higher frequency in cases or in controls. In this way, a 2x2 contingency table was constructed for each gene, where all SNPs fall into one of four bins (up-regulating + less frequent in cases; up-regulating + more frequent in cases; down-regulating + less frequent in cases; down-regulating + more frequent in cases). Fisher’s test was used to test whether the regulatory effect and frequency difference are dependent.

However, Fisher’s test does not account for the effect of LD between SNPs. Therefore, the empirical significance was calculated using an iterative procedure. Genes were ordered according to the number of eQTL SNPs remaining after all filters; genes found in the differentially-expressed NP set (as described in the previous section) were removed; for each gene of interest with n SNPs, 500 genes with the same number n of SNPs are retrieved; if less than 500 genes have the required number of SNPs, genes with n+1 (then n+2, n+3…) SNPs are also retrieved, and n SNPs are randomly selected for analysis in those genes. Each gene is analysed in the same manner as the target gene. Resulting empirical distribution of p-values is used to determine the empirical significance for the gene of interest.

The workflow of the analysis is shown in Fig 1.

thumbnail
Fig 1. Workflow of analysis.

A. Initially, SNP association p-values are produced by an association test. Based on these values, top 1000 SNPs are selected and annotated to nearby genes. B. INRICH software is then used to detect over-represented gene-sets (empty circles denote all genes that were not detected by GWAS). Genes from top 20 such sets are retrieved, and only the ones overlapping with GWAS hits are analysed further. C. Using publicly available expression data from NP samples, this gene list is filtered to retain only differentially expressed genes. D. For each of these remaining targets, known eQTLs are assigned into bins based on effect direction (up- or down- regulating) and frequency distribution in our genotyping data (higher frequency in cases or in controls). Fisher's exact test is used to evaluate the observed distribution. 500 genes with equal or higher count of eQTLs are analysed in the same way, and statistic values generated from this control set are then compared with the target gene statistic to estimate the empirical significance.

https://doi.org/10.1371/journal.pone.0185244.g001

Results

Genetic analysis

Six samples were removed due to heterozygosity > 3 SDs above or below mean, three samples were removed due to > 100 Mendelian errors and one sample due to a mismatch between genotyped and indicated sex. The final dataset after quality-control consists of 782 individuals and 233 409 SNPs. Additional data about the subjects who passed quality control are presented in Table 1. Of the 406 individuals with nasal polyps, 350 were index patients with CRSwNP, 22 of the non-index patients knew they had polyps beforehand and 34 had newly discovered polyps.

thumbnail
Table 1. Descriptive statistics of the participants who passed quality control.

https://doi.org/10.1371/journal.pone.0185244.t001

The Manhattan plot from the DFAM analysis is provided in Fig 2 and the Manhattan plot from the EMMAX analysis in Fig 3. Table 2 lists the 30 SNPs with the strongest associations from the DFAM analysis and Table 3 lists the 30 SNPs with the strongest association values from the EMMAX analysis. None of tested SNPs reached a significance of 5x10-8, the top ranking SNP from the DFAM analysis was rs4629180 with a p-value of 1.47x10-6, the top ranking SNP from the EMMAX analysis was rs2491026 with a p-value of 0.00014.

From the pathway enrichment analysis we extracted the top 20 gene-sets from each of the four INRICH analyses resulting in a combined list of 138 target CRSwNP genes from this study (S1 Table).

Out of our 138 target genes from the INRICH analysis, 36 genes showed a significant difference in mRNA expression levels between nasal polyp tissue and normal tissue (from Plager et al. [32]) (Table 4)

thumbnail
Table 4. Overlap between genes from INRICH analysis and expression data (from Plager et al. [32]).

https://doi.org/10.1371/journal.pone.0185244.t004

Finally, the eQTL analysis for the Blood eQTL dataset showed significantly skewed distributions of eQTLs in cases with CRSwNP compared to controls for HLCS (empirical p-value 0.014) HLA-DRA (empirical p-value 0.02) and BICD2 (empirical p 0.046) (Table 5). The same analysis performed with MuTHER eQTL dataset showed significantly skewed eQTL distribution in cases for VSIR (empirical p-value 0.006), HLCS (empirical p-value 0.014) and BICD2 (empirical p-value 0.016). SLC5A1 also had a skewed distribution with an empirical p-value of 0.052, these results are provided in Table 6.

Discussion

None of the SNPs in this study reached the suggested genome-wide significance level of 5x10-8. However, by using pathway enrichment and post-GWAS analyses, we identified five interesting genes that could be involved in the pathogenesis of CRSwNP in this study: HLCS, HLA-DRA, BICD2, VSIR and SLC5A1. Of these, only HLA-DRA has been presented in previous studies on subjects with CRSwNP [36].

The present study is the first GWAS performed only on subjects with CRSwNP. Using family relationship data and non-transmitted genetic variation for control, as in the TDT, is both a strength and a potential weakness. In most chromosomal regions, we expect the related individuals to be more similar compared with completely unrelated controls. Therefore, when there are differences between related individuals these are more likely to be due to the disease than to general differences in a population (population stratification). In practice this is a strength because it would be expected to lead to less false positive results. However, due to the increased risk of CRSwNP among relatives [7] and the increased prevalence of polyps with higher age [4], some of the relatives who we have defined as not having polyps could develop CRSwNP later in life and therefore be falsely classified as controls in this study. This could possibly have led to missed markers and genes of potential importance. However, most of the relatives in this study are middle-aged or older (mean age 49.4 years) and the prevalence of nasal polyps among them is 13% which makes it unlikely that more than a few percent of the relatives are falsely classified as phenotype negative. One could also argue that there could be subjects with asymptomatic polyps among the 55 relatives who had nasal polyps during endoscopy. However, 21 of them knew they had polyps beforehand and only 34 were unaware of this. The heritability of CRSwNP makes it much more likely that first-degree relatives of patients with CRSwNP would inherit variants associated with CRSwNP than inherit variants associated with asymptomatic nasal polyps. In the event that there are many relatives with polyps but without the predisposition to develop CRS it would indeed influence the GWAS findings, however excluding relatives with polyps would most likely influence the result in a more negative way and severely limit the study.

The HLCS gene was the most significant gene from the eQTL analysis in the Blood dataset and the second most significant in the MuTHER dataset. HLCS is under-expressed in nasal polyp tissue with significantly increased frequency of down-regulating alleles among cases in the eQTL analysis from both the MuTHER and the Blood datasets. This gene encodes the enzyme Holocarboxylase synthetase, which is important for biotin metabolism. Holocarboxylase synthetase itself has not been implicated in CRSwNP but a study from 2013 showed that the enzyme catalyses biotinylation of heat shock protein 72 thereby inducing the expression of the gene RANTES (regulated on activation normal T-expressed and presumably secreted) [37]. RANTES is implicated in multiple studies of CRSwNP; for example a study by Chao et al. found a positive correlation between plasma RANTES protein levels and severity of disease among patients with CRSwNP[38]. RANTES protein has also been detected in nasal polyps using immunological staining [39]. Although it seems counter-intuitive that HLCS is under-expressed when RANTES protein levels are higher in polyp tissue compared with controls, the up-regulation of RANTES might be a counter reaction to initially too low levels during an early disease phase. Further studies are needed to increase the knowledge about the role of HLCS in the pathogenesis of CRSwNP.

HLA-DRA is over-expressed in polyp tissue and has a significantly skewed distribution of eQTLs from the Blood dataset where cases have an increased number of down-regulating alleles. HLA-DRA is one of the major histocompatibility complex, class II genes. Polymorphisms in this gene have been associated with the presence of nasal polyps in asthmatic patients [36]. Additionally, polymorphisms in other HLA class II genes have been linked to CRSwNP [18,40]. Using HLA typing on a series of 29 patients with nasal polyps, with or without asthma, Moloney and Oliver found a significant increase in the haplotype AI/B8 in patients with both nasal polyps and asthma[41].

BICD2 is over-expressed in nasal polyps and up-regulating alleles more common in controls and down regulating slightly more common in cases. The gene product is bicaudal D homolog 2, which has been shown to induce microtubule movement [42]. It is also linked to dominant congenital spinal muscular atrophy [43].

The relationship between gene-expression and eQTLs is reversed for HLA-DRA and BICD2 where cases have an increased number of down-regulating eQTLs even though the gene is over-expressed in polyp tissue. Over-production in the diseased state could possibly be the result of the body compensating for an under-production in the pre-disease state, which hypothetically could have contributed to the development of the disease.

VSIR (V-set immunoregulatory receptor) is over-expressed in polyp tissue and up-regulating alleles from the MuTHER dataset are more common in our cases with CRSwNP compared with unaffected individuals. The gene codes for the protein V-type immunoglobulin domain-containing suppressor of T-cell activation, a member of the Ig superfamily. An experimental study has suggested that it could facilitate tumour invasiveness by regulating cell surface membrane-type 1 matrix metalloproteinase [44]. Lines et al. published an article showing that it also acts a negative checkpoint regulator that suppresses T cell activation [45]. It has not been implicated in CRSwNP in previous studies.

Even though the empirical p-value is slightly higher than 0.05 (empirical p-value 0.052) this study also implicates SLC5A1 as borderline significant. SLC5A1 is under-expressed in polyp tissue and also has an increased frequency of down-regulating alleles in our CRSwNP cases. The gene-product, solute carrier family 5 (sodium/glucose cotransporter) member 1 (SGLT1), is part of a family of sodium-dependent glucose transporters. Once again, this gene has not been associated with CRSwNP but one article suggests a positive substrate cross-regulation of SGLT1 and CFTR [46]. The CFTR gene is highly associated with cystic fibrosis, which often has CRSwNP as one of its clinical features [47]. Furthermore, a study from Varon et al. found an association between CRSwNP (without any other clinical features of cystic fibrosis) and mutations in the CFTR locus [48].

In order to reach a power level necessary for genome-wide analysis we decided to include all patients with CRSwNP regardless of any other diseases which could be associated to this condition. Another reason for this is that the only largely accepted subgrouping of CRS is the division into either CRSsNP or CRSwNP. This is in all likelihood a gross simplification of the pathophysiological and genetic mechanisms behind these conditions and CRS as defined by EPOS is probably a result of a large number of different sub-diseases, each with their own genetic and/or environmental background of which little is known at the present. In this study we did not exclude participants based on potential subgroups since there is still uncertainty surrounding a division into subgroups other than CRSsNP and CRSwNP and we chose to focus on the phenotype CRSwNP itself. Similarly, we chose not to record allergy history due to the controversy surrounding allergy as a possible association to CRSwNP [49]. However, two possible genetic subgroups of CRSwNP; CRSwNP with concomitant asthma and aspirin-exacerbated respiratory disease (AERD) warrant attention in our minds.

Since 51.5% of the participants with nasal polyps also had asthma it is possible that some of the associations could be due to an association with asthma. However, this is less likely since 17.3% of the subjects in the healthy group without polyps also had asthma. The high number of participants with both CRSwNP and asthma could have diluted the association analysis if their SNP-profiles differ significantly from subjects with CRSwNP but without asthma and potentially have made us miss markers of importance in the association analysis. This situation is likely to be countered by the gene enrichment and pathway analyses. HLCS, BICD2, and SLC5A1 have not been connected to asthma, VSIR was implicated with lung function decline in non-asthmatic patients in a genome-wide study published in 2012 but this association could not be confirmed by replication [50]. HLA-DRA is associated to asthma [51] but also, as mentioned above, associated to the presence of nasal polyps in a cohort of asthmatic patients [36].

AERD is another condition linked to CRSwNP, a meta-analysis published in 2015 found a large variation in the prevalence of AERD among patients with CRSwNP among the included studies, the overall prevalence was 9.7% [52]. One of the included studies is from the same geographical region as 96% of our test subjects and found that the prevalence of AERD among subjects with CRSwNP was 6/82 but the numbers were thought to be too small for any meaningful statistical analysis [53]. Although some of our participants probably suffer from AERD, the overall number is in all likelihood too small to influence the association analysis and post-GWAS analyses significantly. HLCS, BICD2, SLC5A1 and VSIR have not been associated with AERD in previous studies. HLA-DRA has not been linked to AERD but other HLA class II genes have [54], however, the same study that linked HLA-DRA to the presence of nasal polyps among asthmatic patients found two HLA-DRA polymorphisms to be potential markers for nasal polyp development in aspirin-tolerant asthma compared to the AERD subgroup [36].

With these caveats in mind, this study is the first of its kind. It is currently the only GWAS performed on CRSwNP and the only study that explores linkage and family-based genome-wide association with regards to this condition. Despite the issue of accurate phenotyping discussed above, this study suggests four novel genes as potential targets of interest for future research as well as once again implicate HLA-DRA.

Conclusion

This study suggests that HLA-DRA as well as four additional genes; HLCS, VSIR, BICD2 and SLC5A1, which have not been previously identified as associated with chronic rhinosinusitis with nasal polyps, could be important for the development of this disease.

Supporting information

S1 Table. List of target genes from the top 20 gene-sets in the INRICH analysis.

https://doi.org/10.1371/journal.pone.0185244.s001

(DOCX)

Acknowledgments

We thank the Gothenburg Genomic Core Facility, the Science for Life Laboratory SNP&SEQ Technology Platform in Uppsala and the Biobanking and Molecular Resource Infrastructure (BBMRI.se) for support. Furthermore, the authors are grateful to Christel Larson for organizing the meetings with all of the test subjects. This study would not be possible without the generous participation of all the families and patients who contributed to the study.

References

  1. 1. Fokkens WJ, Lund VJ, Mullol J, Bachert C, Alobid I, Baroody F, et al. (2012) European Position Paper on Rhinosinusitis and Nasal Polyps 2012. Rhinol Suppl: 3 p preceding table of contents, 1–298.
  2. 2. Klossek JM, Neukirch F, Pribil C, Jankowski R, Serrano E, Chanal I, et al. (2005) Prevalence of nasal polyposis in France: a cross-sectional, case-control study. Allergy 60: 233–237. pmid:15647046
  3. 3. Hedman J, Kaprio J, Poussa T, Nieminen MM (1999) Prevalence of asthma, aspirin intolerance, nasal polyposis and chronic obstructive pulmonary disease in a population-based study. Int J Epidemiol 28: 717–722. pmid:10480701
  4. 4. Johansson L, Akerlund A, Holmberg K, Melen I, Bende M (2003) Prevalence of nasal polyps in adults: the Skovde population-based study. Ann Otol Rhinol Laryngol 112: 625–629. pmid:12903683
  5. 5. Sahlstrand-Johnson P, Ohlsson B, Von Buchwald C, Jannert M, Ahlner-Elmqvist M (2011) A multi-centre study on quality of life and absenteeism in patients with CRS referred for endoscopic surgery. Rhinology 49: 420–428. pmid:21991567
  6. 6. Hopkins C, Gillett S, Slack R, Lund VJ, Browne JP (2009) Psychometric validity of the 22-item Sinonasal Outcome Test. Clin Otolaryngol 34: 447–454. pmid:19793277
  7. 7. Bohman A, Oscarsson M, Holmberg K, Johansson L, Millqvist E, Nasic S, et al. (2015) Heredity of nasal polyps. Rhinology 53: 25–28. pmid:25756074
  8. 8. Oakley GM, Curtin K, Orb Q, Schaefer C, Orlandi RR, Alt JA (2015) Familial risk of chronic rhinosinusitis with and without nasal polyposis: genetics or environment. Int Forum Allergy Rhinol.
  9. 9. Toledano Munoz A, Herraiz Puchol C, Navas Molinero C, Garcia Simal M, Navarro Cunchillos M, Galindo Campillo AN (2008) [Epidemiological study in patients with nasal polyposis]. Acta Otorrinolaringol Esp 59: 438–443. pmid:19080774
  10. 10. Alexiou A, Sourtzi P, Dimakopoulou K, Manolis E, Velonakis E (2011) Nasal polyps: heredity, allergies, and environmental and occupational exposure. J Otolaryngol Head Neck Surg 40: 58–63. pmid:21303603
  11. 11. Cohen NA, Widelitz JS, Chiu AG, Palmer JN, Kennedy DW (2006) Familial aggregation of sinonasal polyps correlates with severity of disease. Otolaryngol Head Neck Surg 134: 601–604. pmid:16564380
  12. 12. Burdett T, Hall P, Hastings E, Hindorff L, Junkins H, Klemm A, et al. (2015) The NHGRI-EBI Catalog of published genome-wide association studies.
  13. 13. Bosse Y, Bacot F, Montpetit A, Rung J, Qu HQ, Engert JC, et al. (2009) Identification of susceptibility genes for complex diseases using pooling-based genome-wide association scans. Hum Genet 125: 305–318. pmid:19184112
  14. 14. Hsu J, Avila PC, Kern RC, Hayes MG, Schleimer RP, Pinto JM (2013) Genetics of chronic rhinosinusitis: state of the field and directions forward. J Allergy Clin Immunol 131: 977–993, 993 e971-975. pmid:23540616
  15. 15. Kuchynkova Z, Macek M Jr., Holcat M, Macek M (1995) [Detection of the G551D mutation in a patient with nasal polyps]. Cas Lek Cesk 134: 212–213. pmid:7538905
  16. 16. Meth MJ, Serota M, Rosenthal DW, Santiago MT, Cavuoto MA (2009) High Frequency of CF Transmembrane Conductance Regulator (CFTR) Mutations in a Population with Persistent Asthma and/or Chronic Rhinosinusitis. Journal of Allergy and Clinical Immunology 123: S159.
  17. 17. Keles B, Cora T, Acar H, Arbag H, Inan Z, Ozturk K, et al. (2008) Evaluation of HLA-A, -B, -Cw, and -DRB1 alleles frequency in Turkish patients with nasal polyposis. Otolaryngology–Head and Neck Surgery 139: 580–585. pmid:18922348
  18. 18. Molnar-Gabor E, Endreffy E, Rozsasi A (2000) HLA-DRB1, -DQA1, and -DQB1 genotypes in patients with nasal polyposis. Laryngoscope 110: 422–425. pmid:10718431
  19. 19. Schubert MS, Hutcheson PS, Graff RJ, Santiago L, Slavin RG (2004) HLA-DQB1*03 in allergic fungal sinusitis and other chronic hypertrophic rhinosinusitis disorders. Journal of Allergy and Clinical Immunology 114: 1376–1383. pmid:15577839
  20. 20. Zielinska-Blizniewska H, Sitarek P, Milonski J, Dziki L, Przybylowska K, Olszewski J, et al. (2012) Association of the −33C/G OSF-2 and the 140A/G LF gene polymorphisms with the risk of chronic rhinosinusitis with nasal polyps in a Polish population. Molecular Biology Reports 39: 5449–5457. pmid:22173721
  21. 21. Palikhe NS, Kim S-H, Cho B-Y, Choi G-S, Kim J-H, Ye Y-M, et al. (2010) IL-13 Gene Polymorphisms are Associated With Rhinosinusitis and Eosinophilic Inflammation in Aspirin Intolerant Asthma. Allergy, Asthma & Immunology Research 2: 134–140.
  22. 22. Bernstein JM, Anon JB, Rontal M, Conroy J, Wang C, Sucheston L (2009) Genetic polymorphisms in chronic hyperplastic sinusitis with nasal polyposis. The Laryngoscope 119: 1258–1264. pmid:19405090
  23. 23. Sitarek P, Zielinska-Blizniewska H, Dziki L, Milonski J, Przybylowska K, Mucha B, et al. (2012) Association of the -14C/G MET and the -765G/C COX-2 gene polymorphisms with the risk of chronic rhinosinusitis with nasal polyps in a Polish population. DNA Cell Biol 31: 1258–1266. pmid:22416915
  24. 24. Pinto JM, Hayes MG, Schneider D, Naclerio RM, Ober C (2008) A genomewide screen for chronic rhinosinusitis genes identifies a locus on chromosome 7q. Laryngoscope 118: 2067–2072. pmid:18622306
  25. 25. Ott J, Kamatani Y, Lathrop M (2011) Family-based designs for genome-wide association studies. Nat Rev Genet 12: 465–474. pmid:21629274
  26. 26. Yang JA, Lee SH, Goddard ME, Visscher PM (2011) GCTA: A Tool for Genome-wide Complex Trait Analysis. American Journal of Human Genetics 88: 76–82. pmid:21167468
  27. 27. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. (2007) PLINK: A tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics 81: 559–575. pmid:17701901
  28. 28. Bozeman M SNP & Variation Suite™. 8.4 ed: Golden Helix, Inc.
  29. 29. Kang HM, Sul JH, Service SK, Zaitlen NA, Kong SY, Freimer NB, et al. (2010) Variance component model to account for sample structure in genome-wide association studies. Nature Genetics 42: 348–U110. pmid:20208533
  30. 30. Clarke GM, Anderson CA, Pettersson FH, Cardon LR, Morris AP, Zondervan KT (2011) Basic statistical analysis in genetic case-control studies. Nature protocols 6: 121–133. pmid:21293453
  31. 31. Lee PH, O'Dushlaine C, Thomas B, Purcell SM (2012) INRICH: interval-based enrichment analysis for genome-wide association studies. Bioinformatics 28: 1797–1799. pmid:22513993
  32. 32. Plager DA, Kahl JC, Asmann YW, Nilson AE, Pallanch JF, Friedman O, et al. (2010) Gene Transcription Changes in Asthmatic Chronic Rhinosinusitis with Nasal Polyps and Comparison to Those in Atopic Dermatitis. Plos One 5.
  33. 33. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. (2013) NCBI GEO: archive for functional genomics data sets-update. Nucleic Acids Research 41: D991–D995. pmid:23193258
  34. 34. Westra HJ, Peters MJ, Esko T, Yaghootkar H, Schurmann C, Kettunen J, et al. (2013) Systematic identification of trans eQTLs as putative drivers of known disease associations. Nature Genetics 45: 1238–U1195. pmid:24013639
  35. 35. Grundberg E, Small KS, Hedman AK, Nica AC, Buil A, Keildson S, et al. (2012) Mapping cis- and trans-regulatory effects across multiple tissues in twins. Nature Genetics 44: 1084-+. pmid:22941192
  36. 36. Kim JH, Park BL, Cheong HS, Pasaje CF, Bae JS, Park JS, et al. (2012) HLA-DRA polymorphisms associated with risk of nasal polyposis in asthmatic patients. Am J Rhinol Allergy 26: 12–17. pmid:22391069
  37. 37. Xue J, Zhou J, Zempleni J (2013) Holocarboxylase synthetase catalyzes biotinylation of heat shock protein 72, thereby inducing RANTES expression in HEK-293 cells. Am J Physiol Cell Physiol 305: C1240–1245. pmid:24133061
  38. 38. Chao PZ, Chou CM, Chen CH (2012) Plasma RANTES and eotaxin levels are correlated with the severity of chronic rhinosinusitis. Eur Arch Otorhinolaryngol 269: 2343–2348. pmid:22271279
  39. 39. Beck LA, Stellato C, Beall LD, Schall TJ, Leopold D, Bickel CA, et al. (1996) Detection of the chemokine RANTES and endothelial adhesion molecules in nasal polyps. J Allergy Clin Immunol 98: 766–780. pmid:8876553
  40. 40. Fajardo-Dolci G, Solorio-Abreu J, Romero-Alvarez JC, Zavaleta-Villa B, Cerezo-Camacho O, Jimenez-Lucio R, et al. (2006) DQA1 and DQB1 association and nasal polyposis. Otolaryngol Head Neck Surg 135: 243–247. pmid:16890076
  41. 41. Moloney JR, Oliver RT (1980) HLA antigens, nasal polyps and asthma. Clin Otolaryngol Allied Sci 5: 183–189. pmid:7408219
  42. 42. Hoogenraad CC, Wulf P, Schiefermeier N, Stepanova T, Galjart N, Small JV, et al. (2003) Bicaudal D induces selective dynein-mediated microtubule minus end-directed transport. EMBO J 22: 6004–6015. pmid:14609947
  43. 43. Oates EC, Rossor AM, Hafezparast M, Gonzalez M, Speziani F, MacArthur DG, et al. (2013) Mutations in BICD2 cause dominant congenital spinal muscular atrophy and hereditary spastic paraplegia. Am J Hum Genet 92: 965–973. pmid:23664120
  44. 44. Sakr MA, Takino T, Domoto T, Nakano H, Wong RW, Sasaki M, et al. (2010) GI24 enhances tumor invasiveness by regulating cell surface membrane-type 1 matrix metalloproteinase. Cancer Sci 101: 2368–2374. pmid:20666777
  45. 45. Lines JL, Pantazi E, Mak J, Sempere LF, Wang L, O'Connell S, et al. (2014) VISTA is an immune checkpoint molecule for human T cells. Cancer Res 74: 1924–1932. pmid:24691993
  46. 46. Mailleau C, Capeau J, Brahimi-Horn MC (1998) Interrelationship between the Na+/glucose cotransporter and CFTR in Caco-2 cells: relevance to cystic fibrosis. J Cell Physiol 176: 472–481. pmid:9699500
  47. 47. Ramsey B, Richardson MA (1992) Impact of sinusitis in cystic fibrosis. J Allergy Clin Immunol 90: 547–552. pmid:1527348
  48. 48. Varon R, Magdorf K, Staab D, Wahn HU, Krawczak M, Sperling K, et al. (1995) Recurrent nasal polyps as a monosymptomatic form of cystic fibrosis associated with a novel in-frame deletion (591del18) in the CFTR gene. Hum Mol Genet 4: 1463–1464. pmid:7581390
  49. 49. Wilson KF, McMains KC, Orlandi RR (2014) The association between allergy and chronic rhinosinusitis with and without nasal polyps: an evidence-based review with recommendations. International Forum of Allergy & Rhinology 4: 93–103.
  50. 50. Imboden M, Bouzigon E, Curjuric I, Ramasamy A, Kumar A, Hancock DB, et al. (2012) Genome-wide association study of lung function decline in adults with and without asthma. Journal of Allergy and Clinical Immunology 129: 1218–1228. pmid:22424883
  51. 51. Ramasamy A, Kuokkanen M, Vedantam S, Gajdos ZK, Couto Alves A, Lyon HN, et al. (2012) Genome-Wide Association Studies of Asthma in Population-Based Cohorts Confirm Known and Suggested Loci and Identify an Additional Association near HLA. PLOS ONE 7: e44008. pmid:23028483
  52. 52. Rajan JP, Wineinger NE, Stevenson DD, White AA (2015) Prevalence of aspirin-exacerbated respiratory disease among asthmatic patients: A meta-analysis of the literature. Journal of Allergy and Clinical Immunology 135: 676–681.e671. pmid:25282015
  53. 53. Johansson L, Bramerson A, Holmberg K, Melen I, Akerlund A, Bende M (2004) Clinical relevance of nasal polyps in individuals recruited from a general population-based study. Acta Otolaryngol 124: 77–81. pmid:14977082
  54. 54. Esmaeilzadeh H, Nabavi M, Amirzargar AA, Aryan Z, Arshi S, Bemanian MH, et al. (2015) HLA-DRB and HLA-DQ genetic variability in patients with aspirin-exacerbated respiratory disease. Am J Rhinol Allergy 29: e63–69. pmid:25975240