Genetic variation in severe cystic fibrosis liver disease is associated with novel mechanisms for disease pathogenesis.

BACKGROUND AND AIMS: It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. APPROACH AND RESULTS: Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. CONCLUSION: These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD.

Mining GWAS and eQTL data for CF lung disease modifiers by gene expression imputation.

Genome wide association studies (GWAS) have identified several genomic loci with candidate modifiers of cystic fibrosis (CF) lung disease, but only a small proportion of the expected genetic contribution is accounted for at these loci. We leveraged expression data from CF cohorts, and Genotype-Tissue Expression (GTEx) reference data sets from multiple human tissues to generate predictive models, which were used to impute transcriptional regulation from genetic variance in our GWAS population. The imputed gene expression was tested for association with CF lung disease severity. By comparing and combining results from alternative approaches, we identified 379 candidate modifier genes. We delved into 52 modifier candidates that showed consensus between approaches, and 28 of them were near known GWAS loci. A number of these genes are implicated in the pathophysiology of CF lung disease (e.g., immunity, infection, inflammation, HLA pathways, glycosylation, and mucociliary clearance) and the CFTR protein biology (e.g., cytoskeleton, microtubule, mitochondrial function, lipid metabolism, endoplasmic reticulum/Golgi, and ubiquitination). Gene set enrichment results are consistent with current knowledge of CF lung disease pathogenesis. HLA Class II genes on chr6, and CEP72, EXOC3, and TPPP near the GWAS peak on chr5 are most consistently associated with CF lung disease severity across the tissues tested. The results help to prioritize genes in the GWAS regions, predict direction of gene expression regulation, and identify new candidate modifiers throughout the genome for potential therapeutic development.

Analysis of a large cohort of cystic fibrosis patients with severe liver disease indicates lung function decline does not significantly differ from that of the general cystic fibrosis population.

Previous reports of lung function in cystic fibrosis (CF) patients with liver disease have shown worse, similar, or even better forced expiratory volume in 1 second (FEV1), compared to CF patients without liver disease. Varying definitions of CF liver disease likely contribute to these inconsistent relationships reported between CF lung function and liver disease. We retrospectively evaluated spirometric data in 179 subjects (62% male; 58% Phe508del homozygous) with severe CF liver disease (CFLD; defined by presence of portal hypertension due to cirrhosis). FEV1 values were referenced to both a normal population (FEV1% predicted) and CF population (CF-specific FEV1 percentile). We utilized a linear mixed model with repeated measures to assess changes in lung function (before and after diagnosis of CFLD), relative to both the normal and CF populations. At diagnosis of CFLD, the mean FEV1 was 81% predicted, or at the 53rd percentile referenced to CF patients without CFLD. There was a significant difference in post-CFLD slope compared to pre-CFLD slope (post-pre) using FEV1% predicted (-1.94, p-value < 0.0001). However, there was insignificant evidence of this difference using the CF-specific FEV1 percentile measure (-0.99, p-value = 0.1268). Although FEV1% predicted values declined in patients following CFLD diagnosis, there was not significant evidence of lung function decline in CF-specific FEV1 percentiles. Thus, the observed study cohort indicates diagnosis of severe CFLD was not associated with worsened CF lung disease when compared to a large CF reference population.

Airway Mucosal Host Defense Is Key to Genomic Regulation of Cystic Fibrosis Lung Disease Severity.

RATIONALE: The severity of cystic fibrosis (CF) lung disease varies widely, even for Phe508del homozygotes. Heritability studies show that more than 50% of the variability reflects non-cystic fibrosis transmembrane conductance regulator (CFTR) genetic variation; however, the full extent of the pertinent genetic variation is not known. OBJECTIVES: We sought to identify novel CF disease-modifying mechanisms using an integrated approach based on analyzing "in vivo" CF airway epithelial gene expression complemented with genome-wide association study (GWAS) data. METHODS: Nasal mucosal RNA from 134 patients with CF was used for RNA sequencing. We tested for associations of transcriptomic (gene expression) data with a quantitative phenotype of CF lung disease severity. Pathway analysis of CF GWAS data (n = 5,659 patients) was performed to identify novel pathways and assess the concordance of genomic and transcriptomic data. Association of gene expression with previously identified CF GWAS risk alleles was also tested. MEASUREMENTS AND MAIN RESULTS: Significant evidence of heritable gene expression was identified. Gene expression pathways relevant to airway mucosal host defense were significantly associated with CF lung disease severity, including viral infection, inflammation/inflammatory signaling, lipid metabolism, apoptosis, ion transport, Phe508del CFTR processing, and innate immune responses, including HLA (human leukocyte antigen) genes. Ion transport and CFTR processing pathways, as well as HLA genes, were identified across differential gene expression and GWAS signals. CONCLUSIONS: Transcriptomic analyses of CF airway epithelia, coupled to genomic (GWAS) analyses, highlight the role of heritable host defense variation in determining the pathophysiology of CF lung disease. The identification of these pathways provides opportunities to pursue targeted interventions to improve CF lung health.

Corrigendum: Novel variation at chr11p13 associated with cystic fibrosis lung disease severity.

[This corrects the article DOI: 10.1038/hgv.2016.20.].

Novel variation at chr11p13 associated with cystic fibrosis lung disease severity.

Published genome-wide association studies (GWASs) identified an intergenic region with regulatory features on chr11p13 associated with cystic fibrosis (CF) lung disease severity. Targeted resequencing in n=377, followed by imputation to n=6,365 CF subjects, was used to identify unrecognized genetic variants (including indels and microsatellite repeats) associated with phenotype. Highly significant associations were in strong linkage disequilibrium and were seen only in Phe508del homozygous CF subjects, indicating a CFTR genotype-specific mechanism.

Features of Severe Liver Disease With Portal Hypertension in Patients With Cystic Fibrosis.

BACKGROUND & AIMS: Liver disease is the third leading cause of death in patients with cystic fibrosis (CF), but features of patients with CF, severe liver disease, and portal hypertension have not been characterized fully. METHODS: We performed a retrospective analysis of data from 561 patients with CF (63% male, 99% with pancreatic insufficiency), liver disease (hepatic parenchymal abnormalities consistent with cirrhosis, confirmed by imaging), and portal hypertension (esophageal varices, portosystemic collaterals, or splenomegaly), with no alternate causes of liver disease. All patients were enrolled in the Genetic Modifier Study of Severe CF Liver Disease at 76 international centers, from January 1999 through July 2013. RESULTS: Male patients were diagnosed with liver disease at a younger age than female patients (10 vs 11 y; P = .01). Splenomegaly was observed in 99% of patients and varices in 71%. Levels of liver enzymes were near normal in most patients. Thrombocytopenia affected 70% of patients and was more severe in patients with varices (88 × 10(9)/L vs 145 × 10(9)/L; P < .0001). Ninety-one patients received liver transplants (16%), at a median age of 13.9 years. Compared with patients who did not receive liver transplants, patients who received liver transplants had lower platelet counts (78 × 10(9)/L vs 113 × 10(9)/L; P < .0001), higher international normalized ratios (P < .0001), and lower levels of albumin (P = .0002). The aminotransferase to platelet ratio index (APRi) and fibrosis index based on 4 factor (FIB-4) values were higher than the diagnostic thresholds for CF liver disease in 96% and in 90% of patients, respectively. Patients who received liver transplants or who had varices had higher APRi and FIB-4 values than patients who did not. CONCLUSIONS: In patients with CF, severe liver disease develops early in childhood (approximately 10 years of age), and is more common in boys than in girls. Patients with varices and those who receive liver transplants have more abnormal platelet counts and APRi and FIB-4 scores.

Genome-wide association meta-analysis identifies five modifier loci of lung disease severity in cystic fibrosis.

The identification of small molecules that target specific CFTR variants has ushered in a new era of treatment for cystic fibrosis (CF), yet optimal, individualized treatment of CF will require identification and targeting of disease modifiers. Here we use genome-wide association analysis to identify genetic modifiers of CF lung disease, the primary cause of mortality. Meta-analysis of 6,365 CF patients identifies five loci that display significant association with variation in lung disease. Regions on chr3q29 (MUC4/MUC20; P=3.3 × 10(-11)), chr5p15.3 (SLC9A3; P=6.8 × 10(-12)), chr6p21.3 (HLA Class II; P=1.2 × 10(-8)) and chrXq22-q23 (AGTR2/SLC6A14; P=1.8 × 10(-9)) contain genes of high biological relevance to CF pathophysiology. The fifth locus, on chr11p12-p13 (EHF/APIP; P=1.9 × 10(-10)), was previously shown to be associated with lung disease. These results provide new insights into potential targets for modulating lung disease severity in CF.

Gene expression in transformed lymphocytes reveals variation in endomembrane and HLA pathways modifying cystic fibrosis pulmonary phenotypes.

Variation in cystic fibrosis (CF) phenotypes, including lung disease severity, age of onset of persistent Pseudomonas aeruginosa (P. aeruginosa) lung infection, and presence of meconium ileus (MI), has been partially explained by genome-wide association studies (GWASs). It is not expected that GWASs alone are sufficiently powered to uncover all heritable traits associated with CF phenotypic diversity. Therefore, we utilized gene expression association from lymphoblastoid cells lines from 754 p.Phe508del CF-affected homozygous individuals to identify genes and pathways. LPAR6, a G protein coupled receptor, associated with lung disease severity (false discovery rate q value = 0.0006). Additional pathway analyses, utilizing a stringent permutation-based approach, identified unique signals for all three phenotypes. Pathways associated with lung disease severity were annotated in three broad categories: (1) endomembrane function, containing p.Phe508del processing genes, providing evidence of the importance of p.Phe508del processing to explain lung phenotype variation; (2) HLA class I genes, extending previous GWAS findings in the HLA region; and (3) endoplasmic reticulum stress response genes. Expression pathways associated with lung disease were concordant for some endosome and HLA pathways, with pathways identified using GWAS associations from 1,978 CF-affected individuals. Pathways associated with age of onset of persistent P. aeruginosa infection were enriched for HLA class II genes, and those associated with MI were related to oxidative phosphorylation. Formal testing demonstrated that genes showing differential expression associated with lung disease severity were enriched for heritable genetic variation and expression quantitative traits. Gene expression provided a powerful tool to identify unrecognized heritable variation, complementing ongoing GWASs in this rare disease.

Meconium ileus in cystic fibrosis is not linked to central repetitive region length variation in MUC1, MUC2, and MUC5AC.

BACKGROUND: Mucins are excellent candidates for contributing to the presence of meconium ileus (MI) in cystic fibrosis (CF) due to their extensive genetic variation and known function in intestinal physiology. The length of variants in mucin central repetitive regions has not been explored as "risk" factors for MI in CF. METHODS: We investigated the length polymorphisms in the central repetitive regions of MUC1, MUC2, and MUC5AC by Southern blot and tested for association with MI in CF subjects. RESULTS: No significant associations were found for the allele sizes of any of the genes with respect to the prevalence of MI (p values=0.33, 0.16, and 0.71 for MUC1, MUC2, and MUC5AC, respectively). CONCLUSIONS: The genetic length variants in the central repetitive region of three MUC genes studied are not associated with MI in subjects with CF.

Genome reference and sequence variation in the large repetitive central exon of human MUC5AC.

Despite modern sequencing efforts, the difficulty in assembly of highly repetitive sequences has prevented resolution of human genome gaps, including some in the coding regions of genes with important biological functions. One such gene, MUC5AC, encodes a large, secreted mucin, which is one of the two major secreted mucins in human airways. The MUC5AC region contains a gap in the human genome reference (hg19) across the large, highly repetitive, and complex central exon. This exon is predicted to contain imperfect tandem repeat sequences and multiple conserved cysteine-rich (CysD) domains. To resolve the MUC5AC genomic gap, we used high-fidelity long PCR followed by single molecule real-time (SMRT) sequencing. This technology yielded long sequence reads and robust coverage that allowed for de novo sequence assembly spanning the entire repetitive region. Furthermore, we used SMRT sequencing of PCR amplicons covering the central exon to identify genetic variation in four individuals. The results demonstrated the presence of segmental duplications of CysD domains, insertions/deletions (indels) of tandem repeats, and single nucleotide variants. Additional studies demonstrated that one of the identified tandem repeat insertions is tagged by nonexonic single nucleotide polymorphisms. Taken together, these data illustrate the successful utility of SMRT sequencing long reads for de novo assembly of large repetitive sequences to fill the gaps in the human genome. Characterization of the MUC5AC gene and the sequence variation in the central exon will facilitate genetic and functional studies for this critical airway mucin.

Genetic modifiers of cystic fibrosis-related diabetes.

Diabetes is a common age-dependent complication of cystic fibrosis (CF) that is strongly influenced by modifier genes. We conducted a genome-wide association study in 3,059 individuals with CF (644 with CF-related diabetes [CFRD]) and identified single nucleotide polymorphisms (SNPs) within and 5' to the SLC26A9 gene that associated with CFRD (hazard ratio [HR] 1.38; P = 3.6 × 10(-8)). Replication was demonstrated in 694 individuals (124 with CFRD) (HR, 1.47; P = 0.007), with combined analysis significant at P = 9.8 × 10(-10). SLC26A9 is an epithelial chloride/bicarbonate channel that can interact with the CF transmembrane regulator (CFTR), the protein mutated in CF. We also hypothesized that common SNPs associated with type 2 diabetes also might affect risk for CFRD. A previous association of CFRD with SNPs in TCF7L2 was replicated in this study (P = 0.004; combined analysis P = 3.8 × 10(-6)), and type 2 diabetes SNPs at or near CDKAL1, CDKN2A/B, and IGF2BP2 were associated with CFRD (P < 0.004). These five loci accounted for 8.3% of the phenotypic variance in CFRD onset and had a combined population-attributable risk of 68%. Diabetes is a highly prevalent complication of CF, for which susceptibility is determined in part by variants at SLC26A9 (which mediates processes proximate to the CF disease-causing gene) and at four susceptibility loci for type 2 diabetes in the general population.

Multiple apical plasma membrane constituents are associated with susceptibility to meconium ileus in individuals with cystic fibrosis.

Variants associated with meconium ileus in cystic fibrosis were identified in 3,763 affected individuals by genome-wide association study (GWAS). Five SNPs at two loci near SLC6A14 at Xq23-24 (minimum P = 1.28 × 10(-12) at rs3788766) and SLC26A9 at 1q32.1 (minimum P = 9.88 × 10(-9) at rs4077468) accounted for ~5% of phenotypic variability and were replicated in an independent sample of affected individuals (n = 2,372; P = 0.001 and 0.0001, respectively). By incorporating the knowledge that disease-causing mutations in CFTR alter electrolyte and fluid flux across surface epithelium into a hypothesis-driven GWAS (GWAS-HD), we identified associations with the same SNPs in SLC6A14 and SLC26A9 and established evidence for the involvement of SNPs in a third solute carrier gene, SLC9A3. In addition, GWAS-HD provided evidence of association between meconium ileus and multiple genes encoding constituents of the apical plasma membrane where CFTR resides (P = 0.0002; testing of 155 apical membrane genes jointly and in replication, P = 0.022). These findings suggest that modulating activities of apical membrane constituents could complement current therapeutic paradigms for cystic fibrosis.

Mucin variable number tandem repeat polymorphisms and severity of cystic fibrosis lung disease: significant association with MUC5AC.

Variability in cystic fibrosis (CF) lung disease is partially due to non-CFTR genetic modifiers. Mucin genes are very polymorphic, and mucins play a key role in the pathogenesis of CF lung disease; therefore, mucin genes are strong candidates as genetic modifiers. DNA from CF patients recruited for extremes of lung phenotype was analyzed by Southern blot or PCR to define variable number tandem repeat (VNTR) length polymorphisms for MUC1, MUC2, MUC5AC, and MUC7. VNTR length polymorphisms were tested for association with lung disease severity and for linkage disequilibrium (LD) with flanking single nucleotide polymorphisms (SNPs). No strong associations were found for MUC1, MUC2, or MUC7. A significant association was found between the overall distribution of MUC5AC VNTR length and CF lung disease severity (p = 0.025; n = 468 patients); plus, there was robust association of the specific 6.4 kb HinfI VNTR fragment with severity of lung disease (p = 6.2×10(-4) after Bonferroni correction). There was strong LD between MUC5AC VNTR length modes and flanking SNPs. The severity-associated 6.4 kb VNTR allele of MUC5AC was confirmed to be genetically distinct from the 6.3 kb allele, as it showed significantly stronger association with nearby SNPs. These data provide detailed respiratory mucin gene VNTR allele distributions in CF patients. Our data also show a novel link between the MUC5AC 6.4 kb VNTR allele and severity of CF lung disease. The LD pattern with surrounding SNPs suggests that the 6.4 kb allele contains, or is linked to, important functional genetic variation.

Genome-wide association and linkage identify modifier loci of lung disease severity in cystic fibrosis at 11p13 and 20q13.2.

A combined genome-wide association and linkage study was used to identify loci causing variation in cystic fibrosis lung disease severity. We identified a significant association (P = 3.34 × 10(-8)) near EHF and APIP (chr11p13) in p.Phe508del homozygotes (n = 1,978). The association replicated in p.Phe508del homozygotes (P = 0.006) from a separate family based study (n = 557), with P = 1.49 × 10(-9) for the three-study joint meta-analysis. Linkage analysis of 486 sibling pairs from the family based study identified a significant quantitative trait locus on chromosome 20q13.2 (log(10) odds = 5.03). Our findings provide insight into the causes of variation in lung disease severity in cystic fibrosis and suggest new therapeutic targets for this life-limiting disorder.

A novel lung disease phenotype adjusted for mortality attrition for cystic fibrosis genetic modifier studies.

Genetic studies of lung disease in cystic fibrosis (CF) are hampered by the lack of a severity measure that accounts for chronic disease progression and mortality attrition. Further, combining analyses across studies requires common phenotypes that are robust to study design and patient ascertainment. Using data from the North American Cystic Fibrosis Modifier Consortium (Canadian Consortium for CF Genetic Studies, Johns Hopkins University CF Twin and Sibling Study, and University of North Carolina/Case Western Reserve University Gene Modifier Study), the authors calculated age-specific CF percentile values of FEV1 which were adjusted for CF age-specific mortality data. The phenotype was computed for 2,061 patients representing the Canadian CF population, 1,137 extreme phenotype patients in the UNC/Case Western study, and 1,323 patients from multiple CF sib families in the CF Twin and Sibling Study. Despite differences in ascertainment and median age, our phenotype score was distributed in all three samples in a manner consistent with ascertainment differences, reflecting the lung disease severity of each individual in the underlying population. The new phenotype score was highly correlated with the previously recommended complex phenotype, but the new phenotype is more robust for shorter follow-up and for extreme ages. A disease progression and mortality-adjusted phenotype reduces the need for stratification or additional covariates, increasing statistical power, and avoiding possible distortions. This approach will facilitate large-scale genetic and environmental epidemiological studies which will provide targeted therapeutic pathways for the clinical benefit of patients with CF.

Adenosine deaminase 1 and concentrative nucleoside transporters 2 and 3 regulate adenosine on the apical surface of human airway epithelia: implications for inflammatory lung diseases.

Adenosine is a multifaceted signaling molecule mediating key aspects of innate and immune lung defenses. However, abnormally high airway adenosine levels exacerbate inflammatory lung diseases. This study identifies the mechanisms regulating adenosine elimination from the apical surface of human airway epithelia. Experiments conducted on polarized primary cultures of nasal and bronchial epithelial cells showed that extracellular adenosine is eliminated by surface metabolism and cellular uptake. The conversion of adenosine to inosine was completely inhibited by the adenosine deaminase 1 (ADA1) inhibitor erythro-9-(2-hydroxy-3-nonyl)adenine (EHNA). The reaction exhibited Km and Vmax values of 24 microM and 0.14 nmol x min(-1) x cm(-2). ADA1 (not ADA2) mRNA was detected in human airway epithelia. The adenosine/mannitol permeability coefficient ratio (18/1) indicated a minor contribution of paracellular absorption. Adenosine uptake was Na+-dependent and was inhibited by the concentrative nucleoside transporter (CNT) blocker phloridzin but not by the equilibrative nucleoside transporter (ENT) blocker dipyridamole. Apparent Km and Vmax values were 17 microM and 7.2 nmol x min(-1) x cm(-2), and transport selectivity was adenosine = inosine = uridine > guanosine = cytidine > thymidine. CNT3 mRNA was detected throughout the airways, while CNT2 was restricted to nasal epithelia. Inhibition of adenosine elimination by EHNA or phloridzin raised apical adenosine levels by >3-fold and stimulated IL-13 and MCP-1 secretion by 6-fold. These responses were reproduced by the adenosine receptor agonist 5'-(N-ethylcarboxamido)adenosine (NECA) and blocked by the adenosine receptor antagonist, 8-(p-sulfophenyl) theophylline (8-SPT). This study shows that adenosine elimination on human airway epithelia is mediated by ADA1, CNT2, and CNT3, which constitute important regulators of adenosine-mediated inflammation.

Glycocalyx restricts adenoviral vector access to apical receptors expressed on respiratory epithelium in vitro and in vivo: role for tethered mucins as barriers to lumenal infection.

Inefficient adenoviral vector (AdV)-mediated gene transfer to the ciliated respiratory epithelium has hindered gene transfer strategies for the treatment of cystic fibrosis lung disease. In part, the inefficiency is due to an absence of the coxsackie B and adenovirus type 2 and 5 receptor (CAR) from the apical membranes of polarized epithelia. In this study, using an in vitro model of human ciliated airway epithelium, we show that providing a glycosylphosphatidylinositol (GPI)-linked AdV receptor (GPI-CAR) at the apical surface did not significantly improve AdV gene transfer efficiency because the lumenal surface glycocalyx limited the access of AdV to apical GPI-CAR. The highly glycosylated tethered mucins were considered to be significant glycocalyx components that restricted AdV access because proteolytic digestion and inhibitors of O-linked glycosylation enhanced AdV gene transfer. To determine whether these in vitro observations are relevant to the in vivo situation, we generated transgenic mice expressing GPI-CAR at the surface of the airway epithelium, crossbred these mice with mice that were genetically devoid of tethered mucin type 1 (Muc1), and tested the efficiency of gene transfer to murine airways expressing apical GPI-human CAR (GPI-hCAR) in the presence and absence of Muc1. We determined that AdV gene transfer to the murine airway epithelium was inefficient even in GPI-hCAR transgenic mice but that the gene transfer efficiency improved in the absence of Muc1. However, the inability to achieve a high gene transfer efficiency, even in mice with a deletion of Muc1, suggested that other glycocalyx components, possibly other tethered mucin types, also provide a significant barrier to AdV interacting with the airway lumenal surface.