ABSTRACT
BACKGROUND: Cystic fibrosis (CF) is caused by deleterious variants in each CFTR gene. We investigated the utility of whole-gene CFTR sequencing when fewer than two pathogenic or likely pathogenic (P/LP) variants were detected by conventional testing (sequencing of exons and flanking introns) of CFTR. METHODS: Individuals with features of CF and a CF-diagnostic sweat chloride concentration with zero or one P/LP variants identified by conventional testing enrolled in the CF Mutation Analysis Program (MAP) underwent whole-gene CFTR sequencing. Replication was performed on individuals enrolled in the CF Genome Project (CFGP), followed by phenotype review and interrogation of other genes. RESULTS: Whole-gene sequencing identified a second P/LP variant in 20/43 MAP enrollees (47 %) and 10/22 CFGP enrollees (45 %) who had one P/LP variant after conventional testing. No P/LP variants were detected when conventional testing was negative (MAP: n = 43; CFGP: n = 13). Genome-wide analysis was unable to find an alternative etiology in CFGP participants with fewer than two P/LP CFTR variants and CF could not be confirmed in 91 % following phenotype re-review. CONCLUSIONS: Whole-gene CFTR analysis is beneficial in individuals with one previously-identified P/LP variant and a CF-diagnostic sweat chloride. Negative conventional CFTR testing indicates that the phenotype should be re-evaluated.
Subject(s)
Cystic Fibrosis Transmembrane Conductance Regulator , Cystic Fibrosis , Humans , Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Cystic Fibrosis/genetics , Female , Male , Phenotype , Mutation , Child , DNA Mutational Analysis , Adolescent , Adult , Genetic Testing/methods , Child, PreschoolABSTRACT
Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.
Subject(s)
Gene Frequency , Genotype , Polymorphism, Single Nucleotide , Software , Humans , Cohort Studies , Linkage Disequilibrium , Genome-Wide Association Study/methods , Genome, Human , Quality Control , Machine Learning , Whole Genome Sequencing/standards , Whole Genome Sequencing/methodsABSTRACT
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.
ABSTRACT
Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.
Subject(s)
Cystic Fibrosis , Humans , Cystic Fibrosis/genetics , Genome-Wide Association Study/methods , Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Patient Acuity , Lung , Microtubule-Associated Proteins/geneticsABSTRACT
Whole-genome sequencing (WGS) is the gold standard for fully characterizing genetic variation but is still prohibitively expensive for large samples. To reduce costs, many studies sequence only a subset of individuals or genomic regions, and genotype imputation is used to infer genotypes for the remaining individuals or regions without sequencing data. However, not all variants can be well imputed, and the current state-of-the-art imputation quality metric, denoted as standard Rsq, is poorly calibrated for lower-frequency variants. Here, we propose MagicalRsq, a machine-learning-based method that integrates variant-level imputation and population genetics statistics, to provide a better calibrated imputation quality metric. Leveraging WGS data from the Cystic Fibrosis Genome Project (CFGP), and whole-exome sequence data from UK BioBank (UKB), we performed comprehensive experiments to evaluate the performance of MagicalRsq compared to standard Rsq for partially sequenced studies. We found that MagicalRsq aligns better with true R2 than standard Rsq in almost every situation evaluated, for both European and African ancestry samples. For example, when applying models trained from 1,992 CFGP sequenced samples to an independent 3,103 samples with no sequencing but TOPMed imputation from array genotypes, MagicalRsq, compared to standard Rsq, achieved net gains of 1.4 million rare, 117k low-frequency, and 18k common variants, where net gains were gained numbers of correctly distinguished variants by MagicalRsq over standard Rsq. MagicalRsq can serve as an improved post-imputation quality metric and will benefit downstream analysis by better distinguishing well-imputed variants from those poorly imputed. MagicalRsq is freely available on GitHub.
Subject(s)
Genome-Wide Association Study , Polymorphism, Single Nucleotide , Humans , Genome-Wide Association Study/methods , Polymorphism, Single Nucleotide/genetics , Calibration , Genotype , Machine LearningABSTRACT
Individuals with cystic fibrosis (CF) develop complications of the gastrointestinal tract influenced by genetic variants outside of CFTR. Cystic fibrosis-related diabetes (CFRD) is a distinct form of diabetes with a variable age of onset that occurs frequently in individuals with CF, while meconium ileus (MI) is a severe neonatal intestinal obstruction affecting â¼20% of newborns with CF. CFRD and MI are slightly correlated traits with previous evidence of overlap in their genetic architectures. To better understand the genetic commonality between CFRD and MI, we used whole-genome-sequencing data from the CF Genome Project to perform genome-wide association. These analyses revealed variants at 11 loci (6 not previously identified) that associated with MI and at 12 loci (5 not previously identified) that associated with CFRD. Of these, variants at SLC26A9, CEBPB, and PRSS1 associated with both traits; variants at SLC26A9 and CEBPB increased risk for both traits, while variants at PRSS1, the higher-risk alleles for CFRD, conferred lower risk for MI. Furthermore, common and rare variants within the SLC26A9 locus associated with MI only or CFRD only. As expected, different loci modify risk of CFRD and MI; however, a subset exhibit pleiotropic effects indicating etiologic and mechanistic overlap between these two otherwise distinct complications of CF.
Subject(s)
Cystic Fibrosis , Diabetes Mellitus , Infant, Newborn, Diseases , Intestinal Obstruction , Cystic Fibrosis/complications , Cystic Fibrosis/genetics , Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Diabetes Mellitus/genetics , Genome-Wide Association Study , Humans , Infant, Newborn , Intestinal Obstruction/complications , Intestinal Obstruction/geneticsABSTRACT
Cystic fibrosis (CF) is a severe genetic disorder that can cause multiple comorbidities affecting the lungs, the pancreas, the luminal digestive system and beyond. In our previous genome-wide association studies (GWAS), we genotyped approximately 8,000 CF samples using a mixture of different genotyping platforms. More recently, the Cystic Fibrosis Genome Project (CFGP) performed deep (approximately 30×) whole genome sequencing (WGS) of 5,095 samples to better understand the genetic mechanisms underlying clinical heterogeneity among patients with CF. For mixtures of GWAS array and WGS data, genotype imputation has proven effective in increasing effective sample size. Therefore, we first performed imputation for the approximately 8,000 CF samples with GWAS array genotype using the Trans-Omics for Precision Medicine (TOPMed) freeze 8 reference panel. Our results demonstrate that TOPMed can provide high-quality imputation for patients with CF, boosting genomic coverage from approximately 0.3-4.2 million genotyped markers to approximately 11-43 million well-imputed markers, and significantly improving polygenic risk score (PRS) prediction accuracy. Furthermore, we built a CF-specific CFGP reference panel based on WGS data of patients with CF. We demonstrate that despite having approximately 3% the sample size of TOPMed, our CFGP reference panel can still outperform TOPMed when imputing some CF disease-causing variants, likely owing to allele and haplotype differences between patients with CF and general populations. We anticipate our imputed data for 4,656 samples without WGS data will benefit our subsequent genetic association studies, and the CFGP reference panel built from CF WGS samples will benefit other investigators studying CF.
ABSTRACT
Chronic Pseudomonas aeruginosa (Pa) infection is associated with increased morbidity and mortality in people with cystic fibrosis (CF). There is no gold standard definition of chronic Pa infection in CF. We compared chronic Pa definitions using encounter-based versus annualized data in the Early Pseudomonas Infection Control (EPIC) Observational study cohort, and subsequently compared annualized chronic Pa definitions across a range of U.S. cohorts spanning decades of CF care. We found that an annualized chronic Pa definition requiring at least 1 Pa+ culture in 3 of 4 consecutive years ("Green 3/4") resulted in chronic Pa metrics similar to established encounter-based modified Leeds criteria definitions, including a similar age at and proportion who fulfilled chronic Pa criteria, and a similar proportion with sustained Pa infection after meeting the chronic Pa definition. The Green 3/4 chronic Pa definition will be valuable for longitudinal analyses in cohorts with limited culture frequency.
Subject(s)
Cystic Fibrosis/microbiology , Pseudomonas Infections/diagnosis , Terminology as Topic , Child , Child, Preschool , Chronic Disease , Cohort Studies , Humans , Infant , Pseudomonas aeruginosa , Registries , Time FactorsABSTRACT
BACKGROUND: Cystic fibrosis (CF) is a recessive condition caused by variants in each CF transmembrane conductance regulator (CFTR) allele. Clinically affected individuals without two identified causal variants typically have no further interrogation of CFTR beyond examination of coding regions, but the development of variant-specific CFTR-targeted treatments necessitates complete understanding of CFTR genotype. METHODS: Whole genome sequences were analyzed on 5,058 individuals with CF. We focused on the full CFTR gene sequence and identified disease-causing variants in three phases: screening for known and structural variants; discovery of novel loss-of-function variants; and investigation of remaining variants. RESULTS: All variants identified in the first two phases and coding region variants found in the third phase were interpreted according to CFTR2 or ACMG criteria (n = 371; 16 [4.3%] previously unreported). Full gene sequencing enabled delineation of 18 structural variants (large insertions or deletions), of which two were novel. Additional CFTR variants of uncertain effect were found in 76 F508del homozygotes and in 21 individuals with other combinations of CF-causing variants. Both causative variants were identified in 98.1% (n = 4,960) of subjects, an increase of 2.3 percentage points from the 95.8% (n = 4,847) who had a registry- or chart-reported disease-causing CFTR genotype. Of the remaining 98 individuals, 78 carried one variant that has been associated with CF (CF-causing [n = 70] or resulting in varying clinical consequences n = 8]). CONCLUSIONS: Complete CFTR gene sequencing in 5,058 individuals with CF identified at least one DNA variant in 99.6% of the cohort that is targetable by current molecular or emerging gene-based therapeutic technologies.
Subject(s)
Cystic Fibrosis , Alleles , Cohort Studies , Cystic Fibrosis/diagnosis , Cystic Fibrosis/genetics , Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Genotype , Humans , MutationABSTRACT
BACKGROUND: Chronic infection with Pseudomonas aeruginosa (P. aeruginosa) is a leading cause of death in patients with cystic fibrosis (CF). Immunobiology of P. aeruginosa infection is complex and not well understood. Chronically infected CF patients generate high levels of antibodies to P. aeruginosa, but this response does not lead to clinical improvement. Therefore, additional studies aimed at identification and understanding of the host factors that influence naturally occurring immune responses to P. aeruginosa are needed. In this investigation, we evaluated the contribution of immunoglobulin GM (γ marker) and KM (κ marker) allotypes to the antibody responses to P. aeruginosa lipopolysaccharide (LPS) O1, O6, O11, and alginate antigens and the broadly-conserved surface polysaccharide expressed by many microbial pathogens, poly-N-acetyl-D-glucosamine (PNAG), in 58 chronically infected CF patients. METHODS: IgG1 markers GM 3 and 17 and IgG2 markers GM 23- and 23+ were determined by a pre-designed TaqMan® genotyping assay. The κ chain determinants KM 1 and 3 were characterized by PCR-RFLP. Antibodies to the LPS O antigens, alginate, and PNAG were measured by an ELISA. RESULTS: Several significant associations were noted with KM alleles. Particular KM 1/3 genotypes were individually and epistatically (with GM 3/17) associated with the level of IgG antibodies to O1, O11, alginate, and PNAG antigens. CONCLUSIONS: Immunoglobulin GM and KM genotypes influence the magnitude of humoral immunity to LPS O, alginate, and PNAG antigens. These results, if confirmed in a larger study population, will be helpful in devising novel immunotherapeutic approaches against P. aeruginosa.
Subject(s)
Cystic Fibrosis/complications , Immunoglobulin Gm Allotypes/immunology , Immunoglobulin Km Allotypes/immunology , Pseudomonas Infections/immunology , Antibody Formation , Antigens, Bacterial/immunology , Female , Genotype , Humans , Immunoglobulin Gm Allotypes/genetics , Immunoglobulin Km Allotypes/genetics , Male , Persistent Infection , Young AdultABSTRACT
BACKGROUND: The CFTR modulator ivacaftor has been variably effective in treating individuals with cystic fibrosis (CF) who harbor CFTR gating variants such as G551D, as well as other classes of CFTR variants when used with other modulators. Because CFTR genotype does not fully explain this variability, defining genetic modifiers of response to modulator therapy is of particular interest to the field of individualized CF drug therapy. Previous studies have proposed that a variant in SLC26A9 (rs7512462) is associated with lung disease severity and with response to treatment with ivacaftor in individuals with CF who carry G551D or gating variants. METHODS: Given the implications for CF treatment, we re-examined the reported associations in three cohorts; patients enrolled in the Twin and Siblings study at Johns Hopkins University, the CF modifier study at the University of North Carolina at Chapel Hill, and the prospective G551D Observational (GOAL) study. The GOAL study was specifically designed to measure lung function response to ivacaftor. RESULTS: We find no association between SLC26A9 (rs7512462) genotype and lung disease severity (n = 272) or change in lung function at one-, three-, and six-month intervals following ivacaftor treatment(n = 141) in individuals with CF who carry at least one G551D variant. CONCLUSIONS: Our inability to replicate this association indicates that rs7512462 genotype should not be used in treatment decisions.
Subject(s)
Antiporters/genetics , Cystic Fibrosis/physiopathology , Sulfate Transporters/genetics , Aminophenols/therapeutic use , Cystic Fibrosis/drug therapy , Cystic Fibrosis/genetics , Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Female , Genotype , Humans , Male , Polymorphism, Single Nucleotide , Prospective Studies , Quinolones/therapeutic use , Respiratory Function Tests , Severity of Illness Index , Young AdultABSTRACT
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.
Subject(s)
Cystic Fibrosis/genetics , Gene Expression/genetics , Genes, Modifier/genetics , Quantitative Trait Loci/genetics , Cohort Studies , Female , Gene Expression Regulation/genetics , Genome-Wide Association Study/methods , Genomics/methods , Humans , MaleABSTRACT
CONTEXT: Individuals with cystic fibrosis (CF) develop a distinct form of diabetes characterized by ß-cell dysfunction and islet amyloid accumulation similar to type 2 diabetes (T2D), but generally have normal insulin sensitivity. CF-related diabetes (CFRD) risk is determined by both CFTR, the gene responsible for CF, and other genetic variants. OBJECTIVE: To identify genetic modifiers of CFRD and determine the genetic overlap with other types of diabetes. DESIGN AND PATIENTS: A genome-wide association study was conducted for CFRD onset on 5740 individuals with CF. Weighted polygenic risk scores (PRSs) for type 1 diabetes (T1D), T2D, and diabetes endophenotypes were tested for association with CFRD. RESULTS: Genome-wide significance was obtained for variants at a novel locus (PTMA) and 2 known CFRD genetic modifiers (TCF7L2 and SLC26A9). PTMA and SLC26A9 variants were CF-specific; TCF7L2 variants also associated with T2D. CFRD was strongly associated with PRSs for T2D, insulin secretion, postchallenge glucose concentration, and fasting plasma glucose, and less strongly with T1D PRSs. CFRD was inconsistently associated with PRSs for insulin sensitivity and was not associated with a PRS for islet autoimmunity. A CFRD PRS comprising variants selected from these PRSs (with a false discovery rate < 0.1) and the genome-wide significant variants was associated with CFRD in a replication population. CONCLUSIONS: CFRD and T2D have more etiologic and mechanistic overlap than previously known, aligning along pathways involving ß-cell function rather than insulin sensitivity. Two CFRD risk loci are unrelated to T2D and may affect multiple aspects of CF. An 18-variant PRS stratifies risk of CFRD in an independent population.
Subject(s)
Cystic Fibrosis/complications , Cystic Fibrosis/genetics , Diabetes Mellitus, Type 2/genetics , Diabetes Mellitus/etiology , Genes, Modifier , Adolescent , Adult , Child , Cohort Studies , Cystic Fibrosis/epidemiology , Diabetes Mellitus/epidemiology , Diabetes Mellitus/genetics , Diabetes Mellitus, Type 2/epidemiology , Female , France/epidemiology , Genetic Predisposition to Disease , Genome-Wide Association Study , Humans , Male , Quantitative Trait, Heritable , Risk Factors , Young AdultABSTRACT
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.
Subject(s)
Cystic Fibrosis/physiopathology , Forced Expiratory Volume/physiology , Hypertension, Portal/physiopathology , Liver Cirrhosis/physiopathology , Lung/physiopathology , Adolescent , Adult , Case-Control Studies , Child , Cystic Fibrosis/complications , Cystic Fibrosis/diagnosis , Female , Humans , Hypertension, Portal/diagnosis , Hypertension, Portal/etiology , Liver/pathology , Liver Cirrhosis/diagnosis , Liver Cirrhosis/etiology , Male , Retrospective Studies , Severity of Illness Index , Spirometry , Young AdultABSTRACT
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.
Subject(s)
Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Cystic Fibrosis/genetics , Genetic Variation , Lung Diseases/genetics , RNA/genetics , Adolescent , Adult , Cohort Studies , Cystic Fibrosis/complications , Cystic Fibrosis/pathology , Disease Progression , Female , Gene Expression Profiling , Gene Expression Regulation , Genome-Wide Association Study , Genomics , Humans , Lung Diseases/etiology , Lung Diseases/pathology , Male , Nasal Mucosa/pathology , Prognosis , RNA/analysis , Risk Assessment , Severity of Illness Index , Young AdultABSTRACT
[This corrects the article DOI: 10.1038/hgv.2016.20.].
ABSTRACT
Extensive phenotypic variability is commonly observed in individuals with Mendelian disorders, even among those with identical genotypes in the disease-causing gene. To determine whether variants within and surrounding CFTR contribute to phenotypic variability in cystic fibrosis (CF), we performed deep sequencing of CFTR in 762 patients homozygous for the common CF-causing variant, F508del. In phase 1, ~200 kb encompassing CFTR and extending 10 kb 5' and 5 kb 3' of the gene was sequenced in 486 F508del homozygotes selected from the extremes of sweat chloride concentration. In phase 2, a 510 kb region, which included the entire topologically associated domain of CFTR, was sequenced in 276 F508del homozygotes drawn from extremes of lung function. An additional 163 individuals who carried F508del and a different CF-causing variant were sequenced to inform haplotype construction. Region-based burden testing of both common and rare variants revealed seven regions of significance (α=0.01), five of which overlapped known regulatory elements or chromatin interactions. Notably, the -80 kb locus known to interact with the CFTR promoter was associated with variation in both CF traits. Haplotype analysis revealed a single rare recombination event (1.9% frequency) in intron 15 of CFTR bearing the F508del variant. Otherwise, the majority of F508del chromosomes were markedly similar, consistent with a single origin of the F508del allele. Together, these high-resolution variant analyses of the CFTR locus suggest a role for non-coding regulatory motifs in trait variation among individuals carrying the common CF allele.
ABSTRACT
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.
ABSTRACT
RATIONALE: Expanding the use of cystic fibrosis transmembrane conductance regulator (CFTR) potentiators and correctors for the treatment of cystic fibrosis (CF) requires precise and accurate biomarkers. Sweat chloride concentration provides an in vivo assessment of CFTR function, but it is unknown the degree to which CFTR mutations account for sweat chloride variation. OBJECTIVES: To estimate potential sources of variation for sweat chloride measurements, including demographic factors, testing variability, recording biases, and CFTR genotype itself. METHODS: A total of 2,639 sweat chloride measurements were obtained in 1,761 twins/siblings from the CF Twin-Sibling Study, French CF Modifier Gene Study, and Canadian Consortium for Genetic Studies. Variance component estimation was performed by nested mixed modeling. MEASUREMENTS AND MAIN RESULTS: Across the tested CF population as a whole, CFTR gene mutations were found to be the primary determinant of sweat chloride variability (56.1% of variation) with contributions from variation over time (e.g., factors related to testing on different days; 13.8%), environmental factors (e.g., climate, family diet; 13.5%), other residual factors (e.g., test variability; 9.9%), and unique individual factors (e.g., modifier genes, unique exposures; 6.8%) (likelihood ratio test, P < 0.001). Twin analysis suggested that modifier genes did not play a significant role because the heritability estimate was negligible (H2 = 0; 95% confidence interval, 0.0-0.35). For an individual with CF, variation in sweat chloride was primarily caused by variation over time (58.1%) with the remainder attributable to residual/random factors (41.9%). CONCLUSIONS: Variation in the CFTR gene is the predominant cause of sweat chloride variation; most of the non-CFTR variation is caused by testing variability and unique environmental factors. If test precision and accuracy can be improved, sweat chloride measurement could be a valuable biomarker for assessing response to therapies directed at mutant CFTR.
Subject(s)
Chlorides/metabolism , Cystic Fibrosis/metabolism , Sweat/metabolism , Adult , Biomarkers/metabolism , Cystic Fibrosis Transmembrane Conductance Regulator/metabolism , Female , Humans , Male , Young AdultABSTRACT
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.