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

In silico search for modifier genes associated with pancreatic and liver disease in Cystic Fibrosis

Abstract

Cystic Fibrosis is the most common lethal autosomal recessive disorder in the white population, affecting among other organs, the lung, the pancreas and the liver. Whereas Cystic Fibrosis is a monogenic disease, many studies reveal a very complex relationship between genotype and clinical phenotype. Indeed, the broad phenotypic spectrum observed in Cystic Fibrosis is far from being explained by obvious genotype-phenotype correlations and it is admitted that Cystic Fibrosis disease is the result of multiple factors, including effects of the environment as well as modifier genes. Our objective was to highlight new modifier genes with potential implications in the lung, pancreatic and liver outcomes of the disease. For this purpose we performed a system biology approach which combined, database mining, literature mining, gene expression study and network analysis as well as pathway enrichment analysis and protein-protein interactions. We found that IFI16, CCNE2 and IGFBP2 are potential modifiers in the altered lung function in Cystic Fibrosis. We also found that EPHX1, HLA-DQA1, HLA-DQB1, DSP and SLC33A1, GPNMB, NCF2, RASGRP1, LGALS3 and PTPN13, are potential modifiers in pancreas and liver, respectively. Associated pathways indicate that immune system is likely involved and that Ubiquitin C is probably a central node, linking Cystic Fibrosis to liver and pancreatic disease. We highlight here new modifier genes with potential implications in Cystic Fibrosis. Nevertheless, our in silico analysis requires functional analysis to give our results a physiological relevance.

Introduction

Cystic Fibrosis (CF) is the most common lethal autosomal recessive disorder in the white population. Its incidence is one in 2,500 with a carrier frequency of one in 25. The Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene, causing CF, was identified in 1989 [1] and located on chromosome 7q31.2, spanning a transcription unit of about 216.7 kb with 27 exons [2]. It encodes a transmembrane protein (CFTR, 1,480 amino acids) which is an ATP-binding cassette transporter functioning as a chloride (Cl-) channel [16]. The CFTR channel’s opening requires phosphorylation by cAMP-dependent protein kinases [68] and hydrolyzable MgATP [8]. The CFTR protein is located at the apical membrane of polarized epithelial cells in diverse tissues including the lungs, sweat ducts, pancreas, gastrointestinal tract, and vas deferens [9, 10].

There are currently 2,009 mutations listed in the CFTR mutation database (http://www.genet.sickkids.on.ca/cftr/app). The most common one, F508del-CFTR (a 3 bp deletion in exon 10 causing loss of the amino acid phenylalanine at position 508), encodes a cAMP-regulated Cl- channel that is retained in the endoplasmic reticulum (ER) during translation and folding and is targeted to the early proteasomal degradation [11]. Besides an altered folding, the F508del-CFTR protein exhibits an altered function with an open time period of approximately 0.1 to 0.3 seconds [12, 13].

The phenotype due to a given mutation depends on its interaction with the second mutated CFTR allele and largely on disease modifiers. Whereas CF is a monogenic disease, there is a very complex relationship between genotype and phenotype [14]. A good genotype-phenotype correlation is observed in the pancreas. Severe and mild mutations are associated with pancreatic insufficiency (PI) and pancreatic sufficiency, respectively [15]. Although CFTR genotype is clearly predictive of PI [16], the genotype-phenotype link is unclear regarding lung and liver diseases. Environmental factors such as infection, nutritional status or socioeconomic status may influence in the pulmonary phenotype in CF. Nevertheless, they can’t explain the degree of variability observed in patients exhibiting a same CFTR genotypes. This is shown by studies in CF twins indicating that other genetic factors may explain the observed variability [17, 18]. Modifier genes may play a significant role in determining the severity of CF. Furthermore, for a single mutation such as F508del, the pulmonary status may range from insignificant to severe [19, 20]. It is therefore obvious that, whereas classic genotype-phenotype studies in CF are important, they are not sufficient and have to be complemented by the search for environmental effects on the phenotype of patients, to increase the basic knowledge of the disease and to develop new therapeutic approaches.

The influence of genetic modulators on CF [21] shows the importance to search for potential candidate gene acting upon CF phenotype. For a given CF genotype, modifier genes may influence the phenotype via the same or parallel pathways [22], leading to complex studies. Therefore, in many works, the candidate gene approach is used. Candidate genes are search in the interactome of the disease gene or in a pathway that is indirectly involved in the course of the disease. Some genome-wide association meta-analysis also identified modifier loci of lung disease severity [23]. In the present work, a different methodology was used.

To search for genetic factors that may play a role in CF, we used a system biology approach. This approach combined database and literature mining, gene expression, and network analysis. Pathway enrichment analysis and protein-protein interactions (PPI) were also searched, to give our results a physiological relevance. This approach allowed us to examine functional relationships between reported genes, and led us to identify novel genes and enriched pathways that may play a role in CF, regarding its lung, pancreatic and liver affections. We found that genes involved in immunity including Ubiquitin C are likely modifiers.

Materials and methods

Data sources and gene selection from databases

The gene selection was performed using a previously used methodology [24].

Genes were first identified utilizing the Online Mendelian Inheritance in Man (OMIM) database (http://www.omim.org/) because it is considered to be the best curated resource for genotype-phenotype relationships [25, 26]. In OMIM, we used the following key words: CF, CF associated genes, PI and Liver Disease. CF and CF associated genes were used to retrieve genes directly involved in CF as well as genes potentially implicated in the disease.

The Comparative Toxicogenomics Database (CTD, http://ctdbase.org/, [27]) was also used. It curates relationships between genes and human diseases in a unique fashion which integrates gene/protein-disease relationships. In CTD, CF, Pancreatic diseases and Liver Disease key words were used. In CTD, disease-gene associations are reported as curated or inferred. We selected curated associations due to a higher confidence than inferred associations. Solely the genes that may be biomarkers of a disease or play a role in the etiology of a disease were selected.

Lastly, the Human Genome Epidemiology encyclopedia (HuGE Navigator, http://www.hugenavigator.net/HuGENavigator/startPagePubLit.do, [28]) which mines the scientific literature on human gene-disease associations was used. It is a database of population-based epidemiologic studies of human genes [29]. In HuGE Navigator we used the CF, PI and Liver Disease key words. The genes with a genetic association above 5 were selected (score formula is described in [30]).

Different key words were used among databases because some of them did not permit to retrieve any result.

Venn diagrams were drawn to find common genes from these databases.

Comparison with published gene lists from bibliography

Differentially expressed genes in CF were retrieved from works published in PubMed (http://www.ncbi.nlm.nih.gov/pubmed). Up- and down-regulated genes were retrieved from a paper [31] in which 4 independent studies [3235] and a re-analyzed list of genes [36] were compared. We used here a list of 75 up-regulated and 114 down-regulated genes which are shared by these studies and having the same direction of expression in at least two studies [3137]. These genes are given in Table 1.

thumbnail
Table 1. Retrieved differentially expressed genes in CF (adapted from [31]).

The lower part of the table shows genes shared between two or more studies when all six studies are combined. Common genes in at least three studies are underlined.

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

To search for up- and down-regulated genes in Liver Disease, we focused on gene expression in Human non alcoholic steato-hepatitis (NASH, [3841]). The rationale that liver disease in NASH and liver of patients with CF share common pathways is based on the fact that steatosis is one of the most common hepatic affection in CF (prevalence of 23–75%). It is observed in 70% of children with liver disease, independently of their nutritional status. Together with a uniform hyperechogenicity of their liver, some pseudomasses may be seen by ultrasound and are due to lobulated fatty structures. In about 60% of the cases, the steatosis is associated with hepatic enzymes elevation (aminotransferases). Because this is observed in children, it is likely not due to alcohol consumption. Therefore, beside the many liver affections observed in livers from CF patients, they present many features of NASH (for review: [42]). Liver Disease key word was not used because it did not permit to retrieve studies with gene analysis. In a first study [38], 16 genes were differentially expressed in subjects with NASH when compared with healthy controls. 12 and 4 genes were significantly under- and 4 over-expressed in NASH, respectively. In a second study [39] comparing NASH samples versus controls, 14 genes were found with more than a two-fold difference. In the third study [40], comparison of non-obese controls patients with NASH exhibited 34 differentially expressed genes (> two fold). We noticed that the expression of 19 genes among 34, were not different in the obese and non-obese controls, showing that the second study using obese controls could be used here. In a last study [41], 9724 genes were differentially expressed in the normal liver versus NASH, (21,619 genes were unchanged).

To search for disease-specific genes in Chronic Pancreatitis, a study showing 152 and 34 genes with increased and decreased expression, respectively, was used [43]. Full name and OMIM identification for these 186 genes was searched. We used another study from which we extracted specific genes in Chronic Pancreatitis [44].

In our study, some genes were withdrawn because the corresponding record was deleted by NCBI or because the term was not found in Gene records. Pseudogenes and non human genes (artifacts) were also withdrawn.

Expression analysis

Datasets restricted to humans, were collected from the Gene Expression Omnibus (GEO, [45]) and GEO2R was used for the analysis [46]. For this purpose, experimental groups were re-assigned. The following datasets were passed through the GEO2R online application: GSE40445 (Gene expression in CF-vs-non CF airway epithelial cells from nasal brushing), GSE48452 (Gene expression in control liver-vs-NASH) and GSE44314 (Gene expression in type 1 diabetes-vs-healthy controls). Type 1 diabetes was used because no PI dataset was found. Samples were assigned to the case or control groups. The GEO2R uses the limma (Linear Models for Microarray Data) package from Bioconductor and was designed to analyze complex experiments involving comparisons between many RNA targets simultaneously, with the idea to fit to a linear model to the expression data for each gene. Whereas limma provides several p-value options, we applied the adjustment to the p-values also called multiple-testing corrections, to correct the occurrence of false positive results. The Benjamini & Hochberg false discovery rate method, commonly used to adjust microarray data, was selected. After GEO2R dataset analysis, genes of interest were selected by using the criteria of a Bonferroni-corrected p-value ≤ 0.05. According to the annotation files, probes that were not corresponding to any genes and genes which were present several times within a single study were eliminated. Differentially expressed transcripts between normal and pathologic were selected and common genes between GSE40445, GSE48452 and GSE44314 were searched by drawing a Venn diagram (not shown).

Analysis of enrichment of KEGG biochemical pathways in CTD

Our target genes were mapped to biological process pathways using CTD (http://ctd.mdibl.org). Biocurators at CTD manually curated gene-disease relationships from the literature as well as pathways. The core data are integrated to construct pathway networks. The VennViewer function, which permits to compare associated datasets, was used for the genes that were found to be potential modifiers in CF, Liver Disease and Chronic Pancreatitis. The Pathway associations function was selected to compare annotated enriched pathways, as well as pathways for disease for these genes. Associated data for interacting genes and diseases, inferred KEGG and REACTOME pathways, or enriched GO terms were downloaded. Annotated pathways for genes are associations established by KEGG and REACTOME curation. The significance of a given pathway is reported as a p-value.

Gene-disease association search

Gene-disease associations were search in CTD in which associations are extracted from the published literature or are derived from the OMIM database using the mim2gene file from the NCBI Gene database. Curated associations were retrieved according to their Bonferroni-corrected p-value < 0.01.

PPI networks design

PPI networks were obtained by STRING 10.0 (Search Tool for the Retrieval of Interacting Genes and Proteins, http://string-db.org/). STRING is a global resource for searching connections (edges) between genes or proteins (nodes) currently covering 9,643,763 proteins. Each edge has a confidence value between 0 and 1, lowest and highest confidence, respectively. The gene’s names were used as an initial input and experimental data sets were selected to display interactive networks. The parameters were defined as follows: all PPI prediction methods were enabled, except for High-throughput Lab Experiments; maximum of 5 interactions by node; cut-off criterion of combined score ≥ 0.7 (interactions at high confidence or better); Homo sapiens.

Results and discussion

Study design

Fig 1 depicts the workflow of our study. We retrieved genes from OMIM, CTD and HuGE. Common genes were searched in CF, pancreatic disease and Liver Disease. Differentially expressed genes in CF were also retrieved from published results in PubMed. Finally, Datasets were collected from the GEO and re-analyzed. Because the CF, Chronic Pancreatitis and Liver Disease keywords did not give any result, type 1 diabetes-vs-healthy control and control liver-vs-NASH were used in GEO. Candidate genes from these databases were compared and analyzed to search for potential common genes and pathways. PPI were further searched to give our results a physiological relevance.

thumbnail
Fig 1. Workflow of the study.

Common genes in CF, PI and Liver Disease were retrieved from OMIM, CTD and HuGE. Candidate genes were also retrieved from published works (PubMed) and datasets were re-analyzed in GEO. Potential modifier genes from the 3 different origins were compared and analyzed to search for pathways and protein-protein interactions.

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

Gene selection from databases

OMIM is a comprehensive, daily-updated human phenotype database, containing more than 12,000 genes of all human genetic diseases. The searches in OMIM using CF, CF associated genes, PI and Liver Disease key words led to 100, 109, 152, 1192 genes, respectively (Fig 2). The intersection of CF and CF associated genes showed 42 different genes and 46 genes were shared by PI and Liver Disease. The 6 genes shared by CF, CF associated genes and PI keywords were S100A8/S100A9 (OMIM Accession ID: 123885), CFTR (OMIM Accession ID: 602421), SCNN1G (OMIM Accession ID: 600761), TGFB1 (OMIM Accession ID: 190180), SERPINA1 (OMIM Accession ID: 107400) and PKD1 (OMIM Accession ID: 601313). In the intersection of CF plus CF associated genes with Liver Disease key words, we found 24 genes (Fig 3). Finally, the number of genes in the intersection of the pooled CF and CF associated genes search and the pooled PI and Liver Disease search was only 4. These genes were CFTR (OMIM Accession ID: 602421), TGFB1 (OMIM Accession ID: 190180), SERPINA1 (OMIM Accession ID: 107400) and PKD1 (OMIM Accession ID: 601313).

thumbnail
Fig 2. Venn diagrams of the OMIM search.

42 different genes were found to be common in CF and CF associated genes (left). 46 genes were shared by PI and Liver Disease (right).

https://doi.org/10.1371/journal.pone.0173822.g002

thumbnail
Fig 3. Venn diagrams showing the common genes in CF, CF associated genes and Liver Disease (OMIM).

24 genes were found in common in CF, CF associated genes and Liver Disease.

https://doi.org/10.1371/journal.pone.0173822.g003

Using the CTD database, 11, 142 and 813 genes were retrieved using the CF, Pancreatic diseases and Liver Disease key words, respectively. Fig 4 presents a Venn diagram showing the shared genes between these key words. Only 3 genes were found to be shared among the 3 used key words: TGFB1 (NCBI Accession ID: 7040), TNFRSF1A (NCBI Accession ID: 7132), CFM1 (NCBI Accession ID: 10167). Surprisingly, among some differences with the OMIM search, the CFTR gene was not found to be present in the list of genes retrieved using the Liver Disease key word in CTD.

thumbnail
Fig 4. Venn diagrams of the CTD database search.

The 3 shared genes in CF, Pancreatic diseases and Liver Disease were TGFB1, TNFRSF1A and CFM1.

https://doi.org/10.1371/journal.pone.0173822.g004

HuGE Navigator database is an integrated, searchable, Web-based knowledge base which mines the scientific literature on human genetic associations and human genome epidemiology [29, 30]. In HuGE Navigator, common genes (score > 5) in PI, CF and Liver Disease were CFTR, SPINK1, GSTP1, PRSS1, ADRB2, GSTM1, GSTT1, HRAS, MIF, OGG1 and TNF. The corresponding OMIM Accession IDs are given in Fig 5.

thumbnail
Fig 5. Venn diagrams of the HuGE Navigator database search.

Common genes in PI, CF and Liver Disease were CFTR, SPINK1, GSTP1, PRSS1, ADRB2, GSTM1, GSTT1, HRAS, MIF, OGG1 and TNF. (score > 5).

https://doi.org/10.1371/journal.pone.0173822.g005

Differential gene expressions retrieved from bibliographic analysis

Differentially expressed genes in CF were retrieved from PubMed. First, a small-scale microarray study performed with human native nasal epithelial cells (F508del homozygous patients vs. controls) was used [31]. We also retrieved genes from other papers aimed to compare gene expression in CF vs normal cells [32]. Another work using non-CF and F508del-CFTR homozygote samples, showed significant changes in the expression in 24 genes (two-sample t-test, p < 0.00001). A three-filter comparative analysis showed that 18 genes were significantly increased and 6 genes were decreased in CF relative to and non-CF samples, respectively [33]. We also retrieved genes from a microarray study in which results from 12 CF and 11 non-CF participants were used [34], and from data in which functional CFTR was absent of the plasma membrane [35]. Finally, we used a paper in which 4 studies evaluating the effect of the F508del-CFTR mutation on airway epithelial cells gene expression were analyzed [36]. A profiled gene expression in CF and non- CF nasal and bronchial epithelium samples, using Illumina HumanRef-8 Expression BeadChips was used. It showed that 863 genes were differentially expressed between CF and non- CF bronchial epithelium and that only 15 were differentially expressed between CF and non- CF nasal epithelium [37]. This indicated that within airways, gene expression varies depending on the region that is studied. Up- and down-regulated genes were compared within these studies [3237] and common genes were searched (Table 1). We found that 75 genes were up-regulated and 114 were down-regulated in human airway cells expressing F508del-CFTR (Table 1, lower part).

Nonalcoholic fatty liver disease has a large spectrum ranging from simple steatosis to NASH, which may lead to progressive fibrosis. In the first study that we used [38], the abundance of intra-hepatic messenger RNA for a broad array of genes was measured. From this study we retrieved differentially expressed genes in NASH-vs-normal liver. Among 6,412 genes, only 16 were differentially expressed in subjects with NASH when compared with controls (Table 2). We used another study using microarrays [39]. 14 genes for NASH-vs-obese controls were found to be up-regulated (Table 2). In a third study, genes in NASH patients and controls were selected [40]. Whereas, 34 genes were differentially expressed in NASH vs non-obese controls, 19 of these genes had no significant differences in obese vs non-obese, suggesting a stronger association of these genes to NASH. Therefore, we used this list of 15 genes (Table 2). Finally, a study with normal, steatotis, NASH with fatty liver, and NASH without fatty liver samples was analyzed. Among 11,633 genes with altered expression out of 33,252 genes, 39 genes were changed in expression, between normal and NASH [21]. Thus, a list of 84 differentially expressed genes from 4 different studies was used in the present study (Table 2).

thumbnail
Table 2. Retrieved differentially expressed genes in Liver Disease (from [3841]).

When possible, names or abbreviations as well as OMIM accession number were added (in italic). The number sign # is used because α1-antitrypsin deficiency is caused by mutation in the SERPINA1 gene (OMIM: 107400). “?” is used when no result was obtained in our OMIM search for the corresponding gene name or when several results could be retrieved. The lower list is retrieved from [21] shows the Solute Carrier Family (SLC, sodium/potassium/chloride transporter family) with differential gene expression in NASH.

https://doi.org/10.1371/journal.pone.0173822.t002

For pancreatic disease, a study with normal and Chronic Pancreatitis specimens was used [43]. Comparison of the expression of 5,600 genes between the normal and Chronic Pancreatitis was performed. GenBank accession numbers of 152 genes with increased expression ([43], Footnote 1) and 34 with decreased genes levels in Chronic Pancreatitis ([43], Footnote 2) were retrieved. We analyzed those genes one by one and removed non human genes and pseudo genes. This led us to list 129 increased genes and 23 decreased genes in Chronic Pancreatitis (Table 3). Because 152 genes were simultaneously increased in pancreatic cancer, only 5 of 5,600 genes were significantly over expressed in Chronic Pancreatitis compared with normal pancreas: Mucin-6 (GenBank Accession Number: L07517), COMP (GenBank Accession Number: L32137), TPSB1 (GenBank Accession Number: M33493), Rearranged Ig-lambda light chain (GenBank Accession Number: X57809) and CRISP-3 (GenBank Accession Number: X95240). An analysis of 6,800 different genes expressed in samples of normal pancreas and Chronic Pancreatitis was studied [44]. 107 genes were predicted to be expressed within cells with Chronic Pancreatitis ([44], Table 1). Genes from both studies [43 and 44] were retrieved and non human, pseudo genes and common genes were withdrawn. Finally, we found 23 decreased genes and 229 increased genes in Chronic Pancreatitis, when compared to normal pancreas (Table 4).

thumbnail
Table 3. Retrieved differentially expressed genes in Chronic Pancreatitis (from [43]).

Some genes were withdrawn according to the reasons explained in Methods.

https://doi.org/10.1371/journal.pone.0173822.t003

thumbnail
Table 4. Combined differentially expressed genes in Chronic Pancreatitis [43, 44].

Some genes were withdrawn according to the reasons explained in Methods. OMIM ID and full names were searched for each gene. 23 genes were found to be decreased and 229 genes were found to be increased in Chronic Pancreatitis (p<0.05).

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

Differential gene expressions retrieved from microarray studies

We also acquired gene array data from the GEO database, a public available archive of individual microarrays studies [46, 47]. GSE40445 [31], is a study of gene expression in human native nasal epithelial cells from F508del-CFTR homozygous patients and non-CF controls. GSE48452 [48] is an expression profiling performed with samples of control human liver, healthy obese, steatosis and NASH). Only controls and NASH groups were re-analyzed here. For PI, we failed to find any array analysis. Therefore, we used a gene expression of type 1 diabetes-vs-classical type 1A diabetes-vs-healthy controls array (GSE44314) [49]. We solely re-analyzed classical type 1A diabetes-vs-healthy controls. 246, 211 and 145 genes were retrieved from our GSE40445, GSE48452 and GSE44314 re-analysis, respectively (Table 5). A Venn diagram was drawn (not shown) but no common genes among the genes selected from the 3 studies were found. Therefore, if there are some common genes in the 3 studies, they do not appear among the most highly regulated genes of each study.

thumbnail
Table 5. List of the genes retrieved from SE40445, GSE48452 and GSE44314.

Genes that were found in at least two groups are underlined.

https://doi.org/10.1371/journal.pone.0173822.t005

Common genes in CF, retrieved from bibliography and microarray studies

We then search for common genes between the genes which were obtained by the bibliographic search for CF genes (Table 1), Liver Disease (Table 2) and Chronic Pancreatitis (Tables 3 and 4) and the genes which were retrieved by the re-analyzing of GSE40445, GSE48452 and GSE44314. The result is summarized in Table 6.

thumbnail
Table 6. Genes involved in CF, Chronic Pancreatitis and Liver Disease.

Common genes present in our re-analyzed data of the GEO database and in our bibliographic analysis. Common genes in CF and in pancreatic and in liver affections are underlined.

https://doi.org/10.1371/journal.pone.0173822.t006

Only 4 genes were found to be common to both groups of genes retrieved from the bibliographic analysis and the GEO database analysis, regarding CF. These genes were IFRD1 (Interferon-related developmental regulator 1, OMIM Accession ID: 603502), IFI16 (Interferon gamma inducible protein 16, OMIM Accession ID: 147586), CCNE2 (Cyclin E2, OMIM Accession ID: 603775) and IGFBP2 (Insulin-like growth factor-binding protein 2, OMIM Accession ID: 146731). This low number of common genes is likely due to the fact that in microarray experiments, the recorded intensities depend on the conditions under which the measurements are made [50].

IFRD1, a histone deacetylase (HDAC)-dependent transcriptional co-regulator expressed during terminal neutrophil differentiation, was already identified of as a modifier gene in CF [51]. It was also found to be common for our CF, CF associated genes and Liver disease keywords search in OMIM.

We failed to find published works from PubMed, showing that IFI16, CCNE2 and IGFBP2 are also possibly involved in CF. Therefore, we searched in previously published works some lines of evidence to link these genes to CF.

IFI16 has a possible role in chronic inflammatory autoimmune disorders. The IFI16 protein is normally expressed in nuclei but may be mislocalized in the cytoplasm and secreted. Indeed, significant levels of extracellular IFI16 protein have been identified in the sera of patients with autoimmune diseases [52]. These data provide evidence for a function of IFI16 upon inflammation and tissue damage, as it is also observed in CF.

The human CCNE2 gene encodes a 404-amino-acid protein, related to cyclin E. It is associated with Cdk2 in a functional kinase complex [53] and is involved in the G1 progression [54]. Activated neutrophils elastase is found in high concentrations in the airways of CF patients [55] and has injurious effects on airway epithelial cells. Because the treatment of normal human bronchial epithelial by elastase results in G arrest, due to cyclin E complex inhibition ([56]), we propose that CCNE2 is possibly involved in the pathophysiology of CF.

The insulin-like growth factors (IGF) stimulates growth of multiple cell types. IGF-II access to the cell surface receptor is mediated by IGFBP-2 [57]. It was previously described that serum IGF-II levels were significantly lower in CF patients than in controls, whereas serum IGFBP-2 was significantly higher [58]. IGF-II/IGFBP-2 molar ratios were also found to be significantly lower in CF. Because chronic inflammation is an important modulator of the IGF/IGFBP system in CF and because IGF2 promotes development of pancreatic beta cells, related to some forms of diabetes mellitus, we also propose that IGFBP-2 is a modifier gene in CF.

Candidate genes in CF were selected due to their involvement in CF-related diseases or in similar diseases. Published studies have mostly examined genes involved in immune or inflammatory response (for review: [59]). These genes are: ACE, ADRβ2, ATB, CAPN10, ClCN2, DEFβ4, ENaC, FcβRII, GCLC, GSTM1, GSTM3, GSTP1, GSTT1, HFE, HLA-I, HLA-II, HLA-III, HSP70, IFN-γ, IL1-β, IL6, IL10, IL18, KCNJ11, MBL2, MIF, NOS1, NOS3, MASP-2, PPARβ, SERPINA1, SERPINA3, SFTPA1, SFTPA2, TLR4, TGFB1, TNFα, TNFα-receptor, and TNFβ [59]. Among the published modifier genes [59], we found using databases that TGFB1, SERPINA1, TNFα-receptor, GSTP1, ADRβ2, GSTM1, GSTT1 and TNFα are also involved in pancreatic and liver diseases (Figs 3 and 4). Among candidate gene modifiers of clinical phenotypes in CF, we also retrieved MBL2, EDNRA, TGFB1, IFRD1, IL8, MSRA, ADIPOR2, TCF7L2 and SERPINA1 (for review: [60]). When compared to our database search, we found that TGFB1 was present in OMIM and CTD and that IFRD1 and SERPINA1 were only present in OMIM, indicating that they are also possibly involved in pancreatic and liver diseases. IFRD1 was also found in our CF-vs-non CF re-analysis of SE40445 (Table 6).

In conclusion, we propose here that IFI16, CCNE2 and IGFBP2 are new candidates of modifiers genes in CF.

Common genes in CF and Chronic Pancreatitis and in CF and Liver Disease retrieved from bibliography and microarray studies.

Common genes in CF-vs-non CF re-analysis of SE40445 and bibliographic search for Chronic Pancreatitis and Liver Disease were EPHX1 (Epoxide Hydrolase 1 Microsomal, OMIM Accession ID: 132810) and SLC33A1 (Solute carrier family 33 (Acetyl-CoA Transporter), member 1, OMIM Accession ID: 603690), respectively.

EPHX1 gene encodes a ubiquitous enzyme (α/β hydrolases), localized in the ER. Mutations in the EPHX1 gene contribute to several diseases [61] in link with tobacco exposure and lung carcinoma [62]. Regarding CF, a low activity of EPHX1 was found to be a risk factor for chronic obstructive pulmonary disease (COPD) in Caucasian [63]. Genetic polymorphisms of GSTP1 and EPHX1 were shown to correlate with oxidative stress markers and lung function in COPD [64]. It has to be noticed that GSTP1 belongs to our list of common genes for PI, CF and Liver Disease search In HuGE Navigator. Mutations in EPHX1 were also suggested to modify the severity of respiratory disorders in CF [65].

SLC33A1 (AT1) is a transporter of the ER membrane permitting the entry of acetyl-CoA into the ER lumen. Its down-regulation induces autophagy and cell death [66]. CFTR needs N-linked glycosylation. Its altered glycosyalation leads to a misfolded protein such as F508del-CFTR which is degradated [67]. The accumulation of this incorrectly folded protein in the ER, triggers the Unfolded Protein Response (UPR; [68, 69]). In CF, excessive UPR may be alleviated by the inhibition of ATF6 which possesses an endogenous inhibitor: XBP1 [69, 45] which controls ERAD by activating the expression of SLC33A1 [70]. Therefore, we propose here that SLC33A1 is likely a strong modulator in CF.

Differentially expressed genes in the bibliographic search for CF and in the Diabetes-vs-healthy re-analysis of GSE44314 were the Major Histocompatibility Complex Class II genes DRB1 and DQA1 (OMIM Accession ID: 142857 and 146880, respectively) and the Desmoplakin gene (DSP, OMIM Accession ID: 125647) (Table 6). HLA-DQA1 and HLA-DQB1 are likely genetic markers of the Celiac Disease [71] which was suggested to be a risk factor in CF [72, 73]. Therefore, together with our present results, we propose HLA-DQA1 and HLA-DQB1 genes, as new modifiers in CF. Desmoplakin is a cytoskeletal linker protein conferring the structural integrity of tissues by linking intermediate filament cytoskeleton from cell to cell (for review: [74]). In remodeled airways, as in CF, abnormalities of the cytokeratins and desmoplakins are observed [75]. According to the importance of intermediate filaments and cell-to-cell communications in CF [76, 77], DSP gene expression is likely involved.

Differentially expressed genes in the bibliographic search for CF and in NASH-vs-healthy re-analysis of GSE48452 were: HLA-DRA (OMIM Accession ID: 142860), GPNMB (Glycoprotein NMB, OMIM Accession ID: 604368), NCF2 (Neutrophil Cytosolic Factor 2, OMIM Accession ID: 608515), RASGRP1 (RAS Guanyl Nucleotide-Releasing Protein 1, OMIM Accession ID: 603962), LGALS3 (Lectin, Galactosidase-Binding, Soluble 3, OMIM Accession ID: 153619), PTPN13 (Protein-Tyrosine Phosphatase, Nonreceptor-Type 13; OMIM Accession ID: 600267) (Table 6).

Involvement of Class II major histocompatibility complex (MHC) genes in CF were discussed above. Because HLA Class II Polymorphism is already known to be a modifier of the pulmonary phenotype in CF [78, 79] our finding is not new.

GPNMB is a transmembrane glycoprotein expressed in the heart, lung, and small intestine [80]. The GPNMB gene is a p53- and androgen-dysregulated gene with antiproliferative and antitumorigenic effects in prostate [81]. To our knowledge, its role in CF has never been studied.

Chronic granulomatous disease (CGD) is the most common phagocyte immunodeficiency due to mutations in a gene encoding NADPH oxidase in phagocytes. As in CF, CGD is characterized by recurrent infections and inflammation in the lungs. Because patients with mutated NCF2 gene present severe infections and CGD [82, 83] NCF2 gene is a good candidate for modifying CF.

RASGRP1 plays a critical role in T cell receptor (TCR) stimulation and is essential in the Natural Killer (NK) cell activation and in the immune response [84, 85]. The NK cells are of particular importance in CF as they play a central role in clearing P. aeruginosa from the lung [86, 87].

LGALS3 gene encodes a β-galactoside-binding lectin which plays a role in apoptosis and innate immunity. Because LGALS3 corrects F508del-CFTR trafficking [88], its gene is a good candidate as a modifier in CF.

The protein encoded by PTPN13 is a tyrosine phosphatase with five PDZ (postsynaptic density protein 95/discs large/zonula occludens-1) domains that bind to plasma membrane and cytoskeleton. The C terminus of CFTR binds various PDZ domain-containing proteins including EBP50, CAL and members of the GRASP family [89]. Despite it was suggested a possible role of reduced PDZ interactions in the accelerated internalization of F508del-CFTR [90], the role of CFTR-PDZ interactions remains incompletely resolved [88]. The PTPN13 gene is therefore of interest and the encoded protein may bind and regulate CFTR.

New potential modifier genes

In conclusion, we found three groups of potential modifier genes. From our bibliographic analysis and microarray studies, we propose here that IFI16 (NCBI ID: 3428), CCNE2 (NCBI ID: 9134) and IGFBP2 (NCBI ID: 3485) are new modifiers candidates in CF. Common genes in CF and Chronic Pancreatitis are EPHX1 (NCBI ID: 2052), HLA-DQA1 (NCBI ID: 3117) and -DQB1 (NCBI ID: 3119) and DSP (NCBI ID: 1832). Modifier genes in CF, in link to Liver Disease, are SLC33A1 (NCBI ID: 9197), GPNMB (NCBI ID: 10457), NCF2 (NCBI ID: 4688), RASGRP1 (NCBI ID: 10125), LGALS3 (NCBI ID: 3958) and PTPN13 (NCBI ID: 5783).

Gene-pathway association of the new potential modifier genes

To better understand the biological function of our new potential modifier genes, we search for their annotated pathways, which are associations established by KEGG [91] and REACTOME [92]. Biological pathways describe biological processes. Therefore, they can be used to integrate and visualize gene products in different conditions (healthy vs disease). KEGG and REACTOME pathway data show known molecular interactions and reaction networks. These data, integrating genes and diseases provide insights into molecular networks involving our depicted candidate genes. Indeed, pathway enrichment analysis permits to decipher biological functions associated with these genes [93]. First, the three common genes in CF (IFI16, CCNE2 and IGFBP2), which were retrieved from Table 6, were submitted to the pathway search, using VennViewer in CTD. 1, 1 and 9 pathways were retrieved for IFI16, IGFBP2 and CCNE2, respectively (Fig 6). No pathway was common between those genes. The 4 common genes in CF and Chronic Pancreatitis and the 6 common genes in CF and Liver Disease were further submitted to the pathway search (Table 7). The only common pathway between CF, Chronic Pancreatitis and Liver Disease was “Immune System” (REACTOME: 6900). The implicated genes were IFI16 for CF, HLA-DQB1 for CF plus Chronic Pancreatitis and LGALS3, NCF2 and RASGRP1 for CF plus Liver Disease (Table 7). Therefore, genes of the Immune System are likely involved in CF and in its association with Chronic Pancreatitis and Liver Disease.

thumbnail
Fig 6. Gene-pathway association of the new potential modifier genes in KEGG and REACTOME.

The 3 common genes in CF, which were retrieved from Table 6, were submitted to the pathway search, using VennViewer in CTD. 1, 1 and 9 pathways were retrieved for IFI16, IGFBP2 and CCNE2, respectively. No pathway was common between those genes.

https://doi.org/10.1371/journal.pone.0173822.g006

thumbnail
Table 7. Associated pathways with the retrieved potential modifiers in CF, in CF and Chronic Pancreatitis and in CF and Liver Disease.

Common pathways are underlined.

https://doi.org/10.1371/journal.pone.0173822.t007

We further search for the diseases that are statistically enriched among our genes. Modifiers candidates in CF (IFI16, CCNE2 and IGFBP2), the common genes in CF and Chronic Pancreatitis (EPHX1, HLA-DQA1 -DQB1, DSP) and in CF in link to Liver Disease (SLC33A1, GPNMB, NCF2, RASGRP1, LGALS3 and PTPN13) were entered as a single data set in CTD. According to CTD, a disease was considered enriched if the proportion of genes annotated to it in our test set was significantly larger than the proportion of all genes annotated to it in the genome. The results are presented Table 8.

thumbnail
Table 8. Associated diseases with the retrieved potential modifiers in CF, in CF and Chronic Pancreatitis and in CF and Liver Disease.

The list of genes was submitted to CTD search and associated diseases were retrieved. In bold: diseases which may present relevant involvement in CF physiopathology.

https://doi.org/10.1371/journal.pone.0173822.t008

Biological function and PPI networks association of the new potential modifier genes

For a full description of the proteins which are encoded by our candidate modifiers genes, we search for their potential interactions using the recent 10.0 version of STRING [94]. We used the basic interaction unit in STRING for functional association between proteins, derived from known experimental data. Given the importance of interactions in protein function, we search for groups of interacting proteins into functional sets, in physical complexes.

We first search for high confidence (≥ 0.7) interactions involving the proteins encoded by IFI16, CCNE2 and IGFBP2 (CF candidates). We added the CFTR protein to the query. As shown in Fig 7A, a network of 9 nodes and 10 edges was obtained (clustering coefficient: 0.744). Whereas IGFBP2 did not belong to the network, we observed that CFTR and IFI16 have a common partner: Ubiquitin C (UBC, OMIM: 191340). IFI16 and UBC are involved in the innate immune signaling pathways [95]. Furthermore, Ubiquitin modification is required for the proteasomal targeting of CFTR, which is ubiquitinated in the ER during assembly and during recycling from the cell surface [96]. Therefore, modifications in the IFI16 and/or UBC gene expression could lead to modifications of the maturation of CFTR, showing the importance of IFI16 as a modifier gene in CF, in link to immunity. CCNE2 binds to CDK2 in a catalytically active complex which modulates the cell cycle progression [97] and CDK inhibition corrects F508del-CFTR proteins [98]. Together with the role of UBC upon the proteasomal targeting of CFTR, the physical and functional complexes involving CCNE2, UBC and CFTR (Fig 7A) indicates that CCNE2 is likely a modifier gene in CF.

thumbnail
Fig 7. PPI networks association of the new potential modifier genes.

The recent 10.0 version of STRING (High confidence ≥ 0.7) was used to search for PPI. A. PPI between proteins encoded by IFI16, CCNE2 and IGFBP2 (CF candidates). We observed that CFTR and IFI16 have a common partner: UBC. B. Interactions involving the proteins encoded by common genes in CF and Chronic Pancreatitis (EPHX1, HLA-DQA1 -DQB1, DSP and CFTR). EPHX1, DSP and CFTR were found to be linked by UBC. C. Networks formed by the proteins encoded by the genes found in CF in link to Liver Disease (SLC33A1, GPNMB, NCF2, RASGRP1, LGALS3 and PTPN13). SLC33A1 forms a protein complex together with CFTR, with UBC as an intermediate. D. Network formed by the proteins encoded by the genes involved in CF, in CF plus Chronic Pancreatitis and in CF plus Liver Disease, which were found in a network linked to the CFTR protein. UBC was observed as a central node. UBC is likely a key component in CF physiopathology. The proteins used in the search tool are marked by a red arrow.

https://doi.org/10.1371/journal.pone.0173822.g007

Interactions involving the proteins encoded by common genes in CF and Chronic Pancreatitis (EPHX1, HLA-DQA1 -DQB1, DSP and CFTR) were searched (Fig 7B). Whereas, HLA-DQA1 and HLA-DQB1 did not belong to any complex, EPHX1, DSP and CFTR were linked by UBC. The Ubiquitin-proteasome system (UPS) is an important regulator for the intracellular trafficking of proteins and a UPS-dependent stabilization of cell-cell contacts involving DSP was shown [99]. Because a low activity of EPHX1 is a risk factor for COPD in Caucasian [63] and because an increased UPS activity is observed in COPD [100], a functional link between EPHX1 and UBC is likely involved in CF and Chronic Pancreatitis.

Fig 7C shows the networks formed by the proteins encoded by the genes found in CF in link to Liver Disease (SLC33A1, GPNMB, NCF2, RASGRP1, LGALS3 and PTPN13). GPNMB, RASGRP1, LGALS3 and PTPN13 were not found to be involved in a network, with the used search parameters. NCF2 was linked to NCF4 and to NCF1. NCF1-NCF2-NCF4 (NCF1-4, neutrophil cytosolic factors 1 to 4) from a multicomponent enzyme system (NADPH-oxidase) responsible for oxidative bursts. We failed to find any information regarding a potential involvement of the NCF1-NCF2-NCF4 enzymatic complex in CF, in PubMed. Interestingly, we observed that SLC33A1 forms a protein complex together with CFTR, with UBC as an intermediate. As mentioned above, IRE1/XBP1 controls ERAD by activating the expression of SLC33A1 [70]. UBC is also involved in ERAD which is a key element of the CFTR degradation. This reinforces our proposition that SLC33A1 is a strong modulator in CF.

Finally, a network was drawn using the proteins encoded by the genes involved in CF, in CF plus Chronic Pancreatitis and in CF plus Liver Disease, which were found in a network linked to the CFTR protein (Fig 7D). UBC was observed as a central node. It is well known that F508del-CFTR protein degradation involves Ubiquitin modification. It was previoulsy shown that the de-ubiquitinating enzyme, Ubiquitin C-terminal hydrolase-L1 (UCH-L1), is highly expressed in CF airway epithelial cells and that there is a positive correlation between UCH-L1 expression and steady state levels of Wt-CFTR protein and F508del-CFTR [101]. Although it is not sufficient to rescue F508del-CFTR, the effect of UCH-L1 upon CFTR processing shows its potential roles in CF. UBC is thus a key component in CF.

Conclusion

CF is characterized by progressive lung disease with inflammation, pancreatic exocrine insufficiency and fatty liver disease. Liver disease is found in two-thirds of sick childrens [102104]. Inflammation is indeed a hallmark of CF and is characterized by bacterial infection, neutrophils infiltrations in the lung and high levels of cytokines. Whereas, inflammatory response is deregulated and excessive, it is not known whether abnormal CFTR is responsible or a consequence. The most common hypothesis is that a defective CFTR leads to a decreased airway surface liquid which fails to clear infected secretions from the lung, triggering the excessive inflammatory response. It is therefore still controversial whether this hyper-inflammation is solely the result of the chronic infection or is a primary event due to CFTR defects [105]. Regarding liver, three types of Liver Disease are observed in CF patients: steatosis, cirrhosis and biliary fibrosis. Fatty infiltration of the liver being found in 30% at biopsy and in 60% at the autopsy of the CF patients, is the most common hepatic alteration [106]. In CF, PI phenotype is a hallmark in patients with two severe alleles such as F508del [16] and CF Related Diabetes (CFRD) is observed in more than 50% of patients. Whereas CFRD is a unique entity, it exhibits common features of both type 1 and 2 diabetes. However, the hallmark of CFRD is insulin deficiency, as in type 1 diabetes in which pancreatic beta cells are destroyed [107].

Despite residual involvement of CFTR’s sequence variation upon lung function, modifiers are of main importance. Indeed, it is clear that differences in the genotype in CF are not solely responsible for the disease variability [108] and large-scale genome-wide association studies (GWAS) were undertaken to search for genetic determinants of phenotypic variation (for review, [23]). Nevertheless, further studies are still needed. Therefore, our aim was to depict new candidate genes with modifying activity in CF, in association with liver and pancreatic diseases. Using databases, bibliographic mining and gene expression analysis we found some candidate genes specific to lung, liver and pancreas. Nevertheless, two groups of different genes could be built, depending on the way they were retrieved. Unsurprisingly, databases search led to known genes involved in CF. Therefore, we focused on genes obtained by microarray results re-analysis, compared to those from our bibliographic analysis. We highlighted 3 genes specific to CF, 4 genes specific to CF and Chronic Pancreatitis and 6 genes for CF and Liver Disease. Nevertheless, we failed to find common genes in CF, Chronic Pancreatitis and Liver Disease. Pathways associated with these retrieved potential modifiers indicated that immune system is likely involved and that Ubiquitin C is the central node linking CF to liver and pancreatic disease. A schematic representation of the results is proposed in Fig 8. Despite we propose here some new genes involved in CF, we are fully aware that some other modifier genes may be missing and that our in silico analysis requires functional analysis to give our results a physiological relevance. Indeed, our analysis may present some limitations due to the statistical combination of results from several studies which may present variable quality and heterogeneity of the used tools. This could lead to misleading results, inherent to any systematic meta-analyses which may be non-exhaustive or over interpreted. Database and literature analysis can also be limiting because of incomplete knowledge base leading to byproduct results. Because in the manual curation of the literature we have not performed an exhaustive review of all the genes and because the used studies have been performed using different methodologies and have different power, no statistical assessment was possible here. Finally, gene expression profiles may be affected by the variability between the used cell models, their intrinsic properties and their type of statistical analysis.

thumbnail
Fig 8. Schematic representation of the results.

Candidate genes are shown in each organ. The corresponding proteins presenting PPI are linked by a red line.

https://doi.org/10.1371/journal.pone.0173822.g008

In conclusion, we propose here a methodology that may be used for some other genetic disease with great variability and identify new candidate modifier genes in CF, related to lung, pancreas and liver, such as Ubiquitin C.

Author Contributions

  1. Conceptualization: PT.
  2. Data curation: PT.
  3. Formal analysis: PT.
  4. Investigation: PT.
  5. Methodology: PT.
  6. Project administration: PT.
  7. Resources: PT.
  8. Software: PT.
  9. Supervision: PT EG CF.
  10. Validation: PT EG CF.
  11. Visualization: PT.
  12. Writing – original draft: PT.
  13. Writing – review & editing: PT EG CF.

References

  1. 1. Riordan JR, Rommens JM, Kerem BS, Alon N, Rozmahel R, Grzelczak Z, et al. Identification of the cystic fibrosis gene: clonic and characterization of complementary DNA. Science 1989;245: 1066–1073. pmid:2475911
  2. 2. Tsui LC, Dorfman R. The cystic fibrosis gene: a molecular genetic perspective. Cold Spring Harb Perspect Med. 2013;3(2):a009472. pmid:23378595
  3. 3. Welsh MJ, Tsui LC, Boat TF, Beaudet AL in: Scriver C.R., Beaudet A.L., Sly W.S. and Valle D. In The Metabolic and Molecular Bases of Inherited Disease (7th edition., McGraw-Hill, New York), 1995;pp. 3799–3876.
  4. 4. Drumm ML, Pope HA, Cliff WH, Rommens JM, Marvin SA, Tsui L, et al. Correction of the cystic fibrosis defect in vitro by retrovirus-madiated gene transfer. Cell 1990;62:1227–1233. pmid:1698126
  5. 5. Rich DP, Anderson MP, Gregory RJ, Cheng SH, Paul S, Jefferson DM, et al. Expression of the cystic fibrosis transmembrane conductance regulator corrects defective chloride channel regulation in cystic fibrosis airway epithelial cells. Nature 1990;347:358–363. pmid:1699126
  6. 6. Szellas T, Nagel G. Apparent affinity of CFTR for ATP is increased by continous kinase activity. FEBS Letters 2003;535:141–146. pmid:12560093
  7. 7. Treharne KJ, Xu Z, Chen JH, Best OG, Cassidy DM, Gruenert DC et al. Inhibition of protein kinase CK2 closes the CFTR Cl− channel, but has no effect on the cystic fibrosis mutant DeltaF508-CFTR. Cell Physiol Biochem. 2009;24:347–360. pmid:19910675
  8. 8. Howell LD, Borchardt R, Kole J, Kaz AM, Randak C Cohn JA. Protein kinase A regulates ATP hydrolysis and dimerization by a CFTR domain. Biochem. J. 2004;378:151–159. pmid:14602047
  9. 9. Bradbury NA. Intracellular CFTR: localization and function. Physiol Rev. 1999;79(Suppl 1):S175–S191.
  10. 10. Bertrand CA, Frizzell RA. The role of regulated CFTR trafficking in epithelial secretion. Am J Physiol Cell Physiol. 2003;285(1):C1–C18. pmid:12777252
  11. 11. Ward CL, Omura S, Kopito RR. Degradation of CFTR by the ubiquitin-proteasome pathway. Cell. 1995;83(1):121–7. pmid:7553863
  12. 12. Haws CM, Nepomuceno IB, Krouse ME, Wakelee H, Law T, Xia Y, et al. Delta F508-CFTR channels: kinetics, activation by forskolin, and potentiation by xanthines. Am J Physiol. 1996;270(5 Pt 1):C1544–55. pmid:8967457
  13. 13. Schultz BD, Frizzell RA, Bridges RJ. Rescue of dysfunctional deltaF508-CFTR chloride channel activity by IBMX. J Membr Biol. 1999;170(1):51–66. pmid:10398760
  14. 14. Mekus F, Ballmann M, Bronsveld I, Bijman J, Veeze H, Tummler B. Categories of deltaF508 homozygous cystic fibrosis twin and sibling pairs with distinct phenotypic characteristics. Twin Res. 2000;3(4):277–293. pmid:11463149
  15. 15. Ferrari M, Cremonesi L. Genotype-phenotype correlation in cystic fibrosis patients. Ann Biol Clin (Paris) 1996;54(6):235–241.
  16. 16. Kristidis P, Bozon D, Corey M, Markiewicz D, Rommens J, Tsui LC, et al. Genetic determination of exocrine pancreatic function in cystic fibrosis. Am J Hum Genet. 1992;50:1178–84. pmid:1376016
  17. 17. Bombieri C, Seia M, Castellani C. Genotypes and phenotypes in cystic fibrosis and cystic fibrosis transmembrane regulator-related disorders. Semin Respir Crit Care Med. 2015;36(2):180–93. pmid:25826586
  18. 18. Bronsveld I, Mekus F, Bijman J, Ballmann M, de Jonge HR, Laabs U, et al. Chloride conductance and genetic background modulate the cystic fibrosis phenotype of Delta F508 homozygous twins and siblings. J Clin Invest. 2001;108:1705–15. pmid:11733566
  19. 19. Kerem E, Corey M, Kerem BS, et al. The relation between genotype and phenotype in cystic fibrosis—analysis of the most common mutation (delta F508) N Engl J Med. 1990;323(22):1517–1522. pmid:2233932
  20. 20. Hubert D, Bienvenu T, Desmazes-Dufeu N, et al. Genotype-phenotype relationships in a cohort of adult cystic fibrosis patients. Eur Respir J. 1996;9(11):2207–2214. pmid:8947061
  21. 21. Mekus F, Laabs U, Veeze H, Tummler B. Genes in the vicinity of CFTR modulate the cystic fibrosis phenotype in highly concordant or discordant F508del homozygous sib pairs. Hum Genet. 2003;112(1):1–11. pmid:12483292
  22. 22. Slavotinek A, Biesecker LG. Genetic modifiers in human development and malformation syndromes, including chaperone proteins. Hum Mol Genet. 2003;12(Spec No 1):R45–R50.
  23. 23. Corvol H, Blackman SM, Boëlle PY, Gallins PJ, Pace RG et al. Genome-wide association meta-analysis identifies five modifier loci of lung disease severity in cystic fibrosis. Nat Commun. 2015;6:8382. pmid:26417704
  24. 24. Lipner EM, Garcia BJ, Strong M. Network Analysis of Human Genes Influencing Susceptibility to Mycobacterial Infections. PLoS One. 2016;11(1):e0146585. pmid:26751573
  25. 25. Hamosh A. Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders. Nucleic Acids Research. 2004;33:D514–D517.
  26. 26. Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. OMIM.org: Online Mendelian Inheritance in Man (OMIM), an online catalog of human genes and genetic disorders. Nucleic Acids Research. 2015;43:D789–D798. pmid:25428349
  27. 27. Davis AP, Murphy CG, Johnson R, Lay JM, Lennon-Hopkins K, Saraceni-Richards C, et al. The Comparative Toxicogenomics Database: update 2013. Nucleic Acids Research. 2013;41:D1104–1114. pmid:23093600
  28. 28. Yu W, Gwinn M, Clyne M, Yesupriya A, Khoury MJ. A navigator for human genome epidemiology. Nat Genet. 2008;40(2):124–5. pmid:18227866
  29. 29. Yesupriya A, Evangelou E, Kavvoura FK, Patsopoulos NA, Clyne M, Walsh MC, et al. Reporting of human genome epidemiology (HuGE) association studies: an empirical assessment. BMC Med Res Methodol. 2008;8:31. pmid:18492284
  30. 30. Yu W, Yesupriya A, Wulf A, Qu J, Khoury MJ, Gwinn M. An open source infrastructure for managing knowledge and finding potential collaborators in a domain-specific subset of PubMed, with an example from human genome epidemiology. BMC Bioinformatics 2007;8: 436. pmid:17996092
  31. 31. Clarke LA, Sousa L, Barreto C, Amaral MD. Changes in transcriptome of native nasal epithelium expressing F508del-CFTR and intersecting data from comparable studies. Respiratory Research 2013;14:38. pmid:23537407
  32. 32. Virella-Lowell I, Herlihy J, Liu B, Lopez C, Cruz P, Muller C, et al. Effects of CFTR, interleukin-10, and Pseudomonas aeruginosa on gene expression profiles in a CF bronchial epithelial cell Line. Mol Ther. 2004;10(3):562–573. pmid:15336656
  33. 33. Zabner J, Scheetz TE, Almabrazi HG, Casavant TL, Huang J, Keshavjee S, et al. CFTR DeltaF508 mutation has minimal effect on the gene expression profile of differentiated human airway epithelia. Am J Physiol Lung Cell Mol Physiol. 2005;289(4):L545–L553. pmid:15937068
  34. 34. Wright JM, Merlo CA, Reynolds JB, Zeitlin PL, Garcia JG, Guggino WB, et al. Respiratory epithelial gene expression in patients with mild and severe cystic fibrosis lung disease. Am J Respir Cell Mol Biol. 2006;35(3):327–336. pmid:16614352
  35. 35. Verhaeghe C, Remouchamps C, Hennuy B, Vanderplasschen A, Chariot A, Tabruyn S, et al. Role of IKK and ERK pathways in intrinsic inflammation of cystic fibrosis airways. Biochem Pharmacol. 2007;73(12):1982–1994. pmid:17466952
  36. 36. Hampton T, Stanton B. A novel approach to analyze gene expression data demonstrates that the DeltaF508 mutation in CFTR downregulates the antigen presentation pathway. Am J Physiol Lung Cell Mol Physiol. 2010;298(4):L473–L482. pmid:20044437
  37. 37. Ogilvie V, Passmore M, Hyndman L, Jones L, Stevenson B, Wilson A, et al. Differential global gene expression in cystic fibrosis nasal and bronchial epithelium. Genomics 2011;98(5):327–36. pmid:21756994
  38. 38. Sreekumar R, Rosado B, Rasmussen D, Charlton M. Hepatic gene expression in histologically progressive nonalcoholic steatohepatitis. Hepatology 2003;38:244–251. pmid:12830008
  39. 39. Younossi ZM, Baranova A, Ziegler K, Del Giacco L, Schlauch K, Born TL, et al. A genomic and proteomic study of the spectrum of nonalcoholic fatty liver disease. Hepatology 2005;42:665–674. pmid:16116632
  40. 40. Younossi ZM, Gorreta F, Ong JP, Schlauch K, Giacco LD, Elariny H, et al. Hepatic gene expression in patients with obesity-related non-alcoholic steatohepatitis. Liver Int. 2005;25:760–771. pmid:15998427
  41. 41. Lake AD, Novak P, Fisher CD, Jackson JP, Hardwick RN, Dean Billheimer D et al. Analysis of Global and Absorption, Distribution, Metabolism, and Elimination Gene Expression in the Progressive Stages of Human Nonalcoholic Fatty Liver Disease. Drug Metab Dispos 2011;39(10):1954–1960. pmid:21737566
  42. 42. Flass T, Narkewicz MR. Cirrhosis and other liver disease in cystic fibrosis. J Cyst Fibros. 2013;12(2):116–124. pmid:23266093
  43. 43. Friess H, Ding J, Kleeff J, Liao Q, Berberat PO, Hammer J, et al. Identification of Disease-specific Genes in Chronic Pancreatitis Using DNA Array Technology. Annals of Surgery 2001;234(6):769–779. pmid:11729383
  44. 44. Binkley CE, Zhang L, Greenson JK, Giordano TJ, Kuick R, Misek D, et al. The Molecular Basis of Pancreatic—FibrosisCommon Stromal Gene Expression in Chronic Pancreatitis and Pancreatic Adenocarcinoma. Pancreas 2004;29(4):254–263. pmid:15502640
  45. 45. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Research 2002;30(1):207–210. pmid:11752295
  46. 46. Barrett T, Wilhite SE, Ledoux P, et al. NCBIGEO: archive for functional genomics data sets-update. Nucleic Acids Research 2013;41(1):D991–D995.
  47. 47. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, et al. NCBI GEO: archive for highthroughput functional genomic data. Nucleic Acids Research 2009;37:D885–D890. pmid:18940857
  48. 48. Ahrens M, Ammerpohl O, von Schönfels W, Kolarova J et al. DNA methylation analysis in nonalcoholic fatty liver disease suggests distinct disease-specific and remodeling signatures after bariatric surgery. Cell Metab. 2013;18(2):296–302. pmid:23931760
  49. 49. Nakata S, Imagawa A, Miyata Y, Yoshikawa A et al. Low gene expression levels of activating receptors of natural killer cells (NKG2E and CD94) in patients with fulminant type 1 diabetes. Immunol Lett. 2013;156(1–2):149–55. pmid:24177169
  50. 50. Chen JJ, Hsueh HM, Delongchamp RR, Lin CJ, Tsai CA. Reproducibility of microarray data: a further analysis of microarray quality control (MAQC) data. BMC Bioinformatics 2007;8:412–426. pmid:17961233
  51. 51. Gu Y, Harley IT, Henderson LB, Aronow BJ, Vietor I, Huber LA, et al. Identification of IFRD1 as a modifier gene for cystic fibrosis lung disease. Nature 2009;458(7241):1039–42. pmid:19242412
  52. 52. Mondini M, Costa S, Sponza S, Gugliesi F, Gariglio M, Landolfo S. The interferon-inducible HIN-200 gene family in apoptosis and inflammation: implication for autoimmunity. Autoimmunity 2010;43:226–31. pmid:20187706
  53. 53. Payton M, Coats S. Cyclin E2, the cycle continues. Int J Biochem Cell Biol. 2002;34(4):315–20. pmid:11854029
  54. 54. Lauper N, Beck ARP, Cariou S, Richman L, Hofmann K, Reith W, et al. Cyclin E2: a novel CDK2 partner in the late G1 and S phases of the mammalian cell cycle. Oncogene 1998;17(20):2637–2643. pmid:9840927
  55. 55. Nakamura H, Yoshimura K, McElvaney NG, Crystal RG. Neutrophil elastase in respiratory epithelial lining fluid of individuals with cystic fibrosis induces interleukin-8 gene expression in a human bronchial epithelial cell line. J Clin Invest. 1992;89:1478–1484. pmid:1569186
  56. 56. Fischer BM, Zheng S, Fan R, Voynow JA. Neutrophil elastase inhibition of cell cycle progression in airway epithelial cells in vitro is mediated by p27kip1. Am J Physiol Lung Cell Mol Physiol. 2007; 293(3):L762–8. pmid:17586698
  57. 57. Kelley Kevin M, Oh Y, Gargosky SE, Gucev Z, Matsumoto T, Hwa V, et al. Insulin-like growth factor-binding proteins (IGFBPs) and their regulatory dynamics. The International Journal of Biochemistry & Cell Biology 1996;28:619–637.
  58. 58. Street ME, Ziveri MA, Spaggiari C, Viani I, Volta C, Grzincich GL, et al. Inflammation is a modulator of the insulin-like growth factor (IGF)/IGF-binding protein system inducing reduced bioactivity of IGFs in cystic fibrosis. Eur J Endocrinol. 2006;154(1):47–52. pmid:16381990
  59. 59. Collaco JM, Cutting GR. Update on gene modifiers in cystic fibrosis. Curr Opin Pulm Med. 2008;14(6):559–66. pmid:18812833
  60. 60. Michael RK, Drumm M. The Influence of Genetics on Cystic Fibrosis Phenotypes. Cold Spring Harb Perspect Med. 2012;2:a009548. pmid:23209180
  61. 61. Václavíková R, Hughes DJ, Souček P. Microsomal epoxide hydrolase 1 (EPHX1): Gene, structure, function, and role in human disease. Gene 2015;571(1):1–8. pmid:26216302
  62. 62. Peluso ME, Munnia A, Srivatanakul P, Jedpiyawongse A, Sangrajrang S, Ceppi M, et al. DNA adducts and combinations of multiple lung cancer at-risk alleles in environmentally exposed and smoking subjects. Environ. Mol. Mutagen. 2013;54:375–383. pmid:23797975
  63. 63. Li H, Fu WP, Hong ZH. Microsomal epoxide hydrolase gene polymorphisms and risk of chronic obstructive pulmonary disease: a comprehensive meta-analysis. Oncol. Lett. 2013;5:1022–1030. pmid:23426996
  64. 64. Vibhuti A, Arif E, Deepak D, Singh B, Qadar Pasha MA. Genetic polymorphisms of GSTP1 and mEPHX correlate with oxidative stress markers and lung function in COPD. Biochem Biophys Res Commun. 2007;359(1):136–42. pmid:17532303
  65. 65. Korytina GF, Ianbaeva DG, Viktorova TV. Role of polymorphic variants of cytochrome P450 genes (CYP1A1, CYP2E1) and microsomal epoxide hydrolase (mEPHX) in pathogenesis of cystic fibrosis and chronic respiratory tract diseases. Mol Biol (Mosk). 2003;37(5):784–92.
  66. 66. Jonas MC, Pehar M, Puglielli L. AT-1 is the ER membrane acetyl-CoA transporter and is essential for cell viability. J. Cell Sci. 2010;123:3378–3388. pmid:20826464
  67. 67. Farinha CM, Amaral MD. Most F508del-CFTR is targeted to degradation at an early folding checkpoint and independently of calnexin. Mol Cell Biol. 2005;25(12):5242–52. pmid:15923638
  68. 68. Chakrabarti A, Chen AW, Varner JD. A review of the mammalian unfolded protein response. Biotechnol Bioeng. 2011;108(12):2777–93. pmid:21809331
  69. 69. Kerbiriou M, Le Drévo MA, Férec C, Trouvé P. Coupling cystic fibrosis to endoplasmic reticulum stress: differential role of Grp78 and ATF6. BBA Molecular Basis of Diseases 2007;1772(11–12):1236–49.
  70. 70. Pehar M, Cabell Jonas M, Hare TM, Puglielli L. SLC33A1/AT-1 Protein Regulates the Induction of Autophagy Downstream of IRE1/XBP1 Pathway. The Journal of biological chemistry 2012;287(35):29921–29930. pmid:22787145
  71. 71. Zubillaga P, Vidales MC, Zubillaga I, Ormaechea V, García-Urkía N, Vitoria JC. HLA-DQA1 and HLA-DQB1 Genetic Markers and Clinical Presentation in Celiac Disease. Journal of Pediatric Gastroenterology and Nutrition 2002;34:548–554. pmid:12050583
  72. 72. Walkowiak J, Blask-Osipa A, Lisowska A, Oralewska B, Pogorzelski A, Cichy W, et al. Cystic fibrosis is a risk factor for celiac disease. Acta Biochim Pol. 2010;57(1):115–8. pmid:20300660
  73. 73. Fluge G, Olesen HV, Gilljam M, Meyer P, Pressler T, Storrösten OT, et al. Co-morbidity of cystic fibrosis and celiac disease in Scandinavian cystic fibrosis patients. J Cyst Fibros. 2009;8(3):198–202. pmid:19303374
  74. 74. Schmidt A, Heid HW, Schäfer S, Nuber UA, Zimbelmann R, Franke WW. Desmosomes and cytoskeletal architecture in epithelial differentiation: cell type-specific plaque components and intermediate filament anchorage. Eur J Cell Biol. 1994;65(2):229–45. pmid:7720719
  75. 75. Brezillon S, Dupuit F, Hinnrasky J, Marchand V, Kahin N, Tummler B, et al. Decreased expression of CFTR protein in remodeled human nasal epithelium from non-cystic fibrosis patients. Lab. Invest. 1995;72:191–200. pmid:7531792
  76. 76. Edelman A. Cytoskeleton and CFTR. The International Journal of Biochemistry & Cell Biology 2014;52:68–72.
  77. 77. Losa D, Chanson M. The lung communication network. Cell Mol Life Sci. 2015;72(15):2793–808. pmid:26100513
  78. 78. Aron Y, Polla BS, Bienvenu T, Dall'ava J, Dusser D, Hubert D. HLA class II polymorphism in cystic fibrosis. A possible modifier of pulmonary phenotype. Am J Respir Crit Care Med. 1999; 159(5 Pt1):1464–8.
  79. 79. Corvol H, Blackman SM, Boëlle PY, Gallins PJ, Pace RG, Stonebraker JR, et al. Genome-wide association meta-analysis identifies five modifier loci of lung disease severity in cystic fibrosis. Nat Commun. 2015;6:8382. pmid:26417704
  80. 80. Safadi FF, Xu J, Smock SI, Rico MC, Owen TA, Popoff SN. Cloning and characterization of osteoactivin, a novel cDNA expressed in osteoblasts. J Cell Biochem. 2002;84:12–26.
  81. 81. Tsui KH, Chang YL, Feng TH, Chang PL, Juang HH. Glycoprotein transmembrane nmb: an androgen-downregulated gene attenuates cell invasion and tumorigenesis in prostate carcinoma cells. Prostate 2012;72(13):1431–42. pmid:22290289
  82. 82. Kuhns DB, Alvord WG, Heller T, Feld JJ, Pike KM, Marciano BE, et al. Residual NADPH oxidase and survival in chronic granulomatous disease. N Engl J Med. 2010;363(27):2600–10. pmid:21190454
  83. 83. Ben-Farhat K, Ben-Mustapha I, Ben-Ali M, Rouault K, Hamami S, Mekki N, et al. A Founder Effect of c.257 + 2T > C Mutation in NCF2 Gene Underlies Severe Chronic Granulomatous Disease in Eleven Patients. J Clin Immunol. 2016;[Epub ahead of print].
  84. 84. Ebinu JO, Bottorff DA, Chan EY, Stang SL, Dunn RJ, Stone JC. RasGRP, a Ras guanyl nucleotide- releasing protein with calcium- and diacylglycerol-binding motifs. Science 1998;280:1082–1086. pmid:9582122
  85. 85. Dower NA, Stang SL, Bottorff DA, Ebinu JO, Dickie P, Ostergaard HL, et al. RasGRP is essential for mouse thymocyte differentiation and TCR signaling. Nat Immunol. 2000;1(4):317–21. pmid:11017103
  86. 86. Lee SH, Yun S, Lee J, Kim MJ, Piao ZH, Jeong M, et al. RasGRP1 is required for human NK cell function. J Immunol. 2009;183(12):7931–8. pmid:19933860
  87. 87. Nieuwenhuis EE, Matsumoto T, Lindenbergh D, Willemsen R, Kaser A, Simons-Oosterhuis Y, et al. Cd1d-dependent regulation of bacterial colonization in the intestine of mice. J Clin Invest. 2009;119(5):1241–50. pmid:19349688
  88. 88. Trzcinska-Daneluti AM, Ly D, Huynh L, Jiang C, Fladd C, Rotin D. High-content functional screen to identify proteins that correct F508del-CFTR function. Mol Cell Proteomics 2009;8(4):780–90. pmid:19088066
  89. 89. Guggino WB, Stanton BA. New insights into cystic fibrosis: molecular switches that regulate CFTR. Nat. Rev. Mol. Cell Biol. 2006;7:426–436. pmid:16723978
  90. 90. Valentine CD, Lukacs GL, Verkman AS, Haggie PM. Reduced PDZ interactions of rescued ΔF508CFTR increases its cell surface mobility. J Biol Chem. 2012;287(52):43630–8. pmid:23115232
  91. 91. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–D484. pmid:18077471
  92. 92. Croft D. Building models using Reactome pathways as templates In Silico Systems Biology: Springer 2013;p. 273–83.
  93. 93. Khatri P, Sirota M, Butte AJ. Ten years of pathway analysis: Current approaches and outstanding challenges. PLoS Comput. Biol. 2012;8:373.
  94. 94. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52. pmid:25352553
  95. 95. Cui J, Chen Y, Wang HY, Wang RF. Mechanisms and pathways of innate immune activation and regulation in health and cancer. Hum Vaccin Immunother. 2014;10(11):3270–85. pmid:25625930
  96. 96. Sato S, Ward CL, Kopito RR. Cotranslational ubiquitination of cystic fibrosis transmembrane conductance regulator in vitro. J Biol Chem. 1998;273(13):7189–92. pmid:9516408
  97. 97. Gudas JM, Payton M, Thukral S, Chen E, Bass M, Robinson MO., et al. Coats S. Cyclin E2, a Novel G1 Cyclin That Binds Cdk2 and Is Aberrantly Expressed in Human Cancers. Mol. Cell. Biol. 1999;1:612–622.
  98. 98. Norez C, Vandebrouck C, Bertrand J, Noel S, Durieu E, Oumata N, et al. Roscovitine is a proteostasis regulator that corrects the trafficking defect of F508del-CFTR by a CDK-independent mechanism. Br J Pharmacol. 2014;171(21):4831–4849. pmid:25065395
  99. 99. Löffek S, Bruckner-Tuderman L, Magin TM. Involvement of the ubiquitin-proteasome system in the stabilization of cell-cell contacts in human keratinocytes. Exp Dermatol. 2012; 21(10):791–3. pmid:22882483
  100. 100. Ottenheijm CA, Heunks LM, Li YP, Jin B, Minnaard R, van Hees HW, et al. Activation of the ubiquitin-proteasome pathway in the diaphragm in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2006;174(9):997–1002. pmid:16917114
  101. 101. Henderson MJ, Vij N, Zeitlin PL. Ubiquitin C-terminal hydrolase-L1 protects cystic fibrosis transmembrane conductance regulator from early stages of proteasomal degradation. J Biol Chem. 2010;285(15):11314–25. pmid:20147297
  102. 102. Kelly T, Buxbaum J. Gastrointestinal Manifestations of Cystic Fibrosis Dig Dis Sci. 2015;60:1903–1913. pmid:25648641
  103. 103. Karlas T, Neuschulz M, Oltmanns A, et al. Non-invasive evaluation of cystic fibrosis related liver disease in adults with ARFI, transient elastography and different fibrosis scores. Plos One. 2012;7:8.
  104. 104. Kobelska-Dubiel N, Klincewicz B, Cichy W. Liver disease in cystic fibrosis. Prz Gastroenterol. 2014;9(3):136–141. pmid:25097709
  105. 105. Machen TE. Innate immune response in CF airway epithelia: hyperinflammatory? Am J Physiol Cell Physiol. 2006;291:C218–C230. pmid:16825601
  106. 106. Akata D, Akhan O. Liver manifestations of cystic fibrosis. European Journal of Radiology. 2007;61(1):11–17. pmid:17174503
  107. 107. Laguna TA, Nathan BM, Moran A. Managing diabetes in cystic fibrosis. Diabetes Obes Metab. 2010;12(10):858–864. pmid:20920037
  108. 108. Kerem B, Kerem E. The molecular basis for disease variability in cystic fibrosis. Eur J Hum Genet. 1996; 4:65–73. pmid:8744024