Validating organoid-derived human intestinal monolayers for personalized therapy in cystic fibrosis.

Highly effective drugs modulating the defective protein encoded by the CFTR gene have revolutionized cystic fibrosis (CF) therapy. Preclinical drug-testing on human nasal epithelial (HNE) cell cultures and 3-dimensional human intestinal organoids (3D HIO) are used to address patient-specific variation in drug response and to optimize individual treatment for people with CF. This study is the first to report comparable CFTR functional responses to CFTR modulator treatment among patients with different classes of CFTR gene variants using the three methods of 2D HIO, 3D HIO, and HNE. Furthermore, 2D HIO showed good correlation to clinical outcome markers. A larger measurable CFTR functional range and access to the apical membrane were identified as advantages of 2D HIO over HNE and 3D HIO, respectively. Our study thus expands the utility of 2D intestinal monolayers as a preclinical drug testing tool for CF.

Caution advised in the use of CFTR modulator treatment for individuals harboring specific CFTR variants.

In December 2020, the U.S. Food and Drug Administration (FDA) expanded the list of CFTR variants approved for treatment with CFTR modulators drugs from 39 to 183. Clinicians should be aware that individuals harboring certain variants approved for treatment may not respond to or benefit from this therapy. After review, the expert panel leading the CFTR2 project identified four categories of variants that may not result in a clinical response to modulator treatment: 15 variants assigned as non CF-causing; 45 variants of unknown significance; six variants known or suspected to cause mis-splicing as their primary defect rather than an amino acid substitution; and eight variants known to occur together in cis with another deleterious variant not expected to lead to CFTR protein (nonsense or frameshift). The potential risks and benefits of CFTR modulator therapy should be considered carefully for individuals harboring these variants.

Genetic evidence supports the development of SLC26A9 targeting therapies for the treatment of lung disease.

Over 400 variants in the cystic fibrosis (CF) transmembrane conductance regulator (CFTR) are CF-causing. CFTR modulators target variants to improve lung function, but marked variability in response exists and current therapies do not address all CF-causing variants highlighting unmet needs. Alternative epithelial ion channel/transporters such as SLC26A9 could compensate for CFTR dysfunction, providing therapeutic targets that may benefit all individuals with CF. We investigate the relationship between rs7512462, a marker of SLC26A9 activity, and lung function pre- and post-treatment with CFTR modulators in Canadian and US CF cohorts, in the general population, and in those with chronic obstructive pulmonary disease (COPD). Rs7512462 CC genotype is associated with greater lung function in CF individuals with minimal function variants (for which there are currently no approved therapies; p = 0.008); and for gating (p = 0.033) and p.Phe508del/ p.Phe508del (p = 0.006) genotypes upon treatment with CFTR modulators. In parallel, human nasal epithelia with CC and p.Phe508del/p.Phe508del after Ussing chamber analysis of a combination of approved and experimental modulator treatments show greater CFTR function (p = 0.0022). Beyond CF, rs7512462 is associated with peak expiratory flow in a meta-analysis of the UK Biobank and Spirometa Consortium (p = 2.74 × 10-44) and provides p = 0.0891 in an analysis of COPD case-control status in the UK Biobank defined by spirometry. These findings support SLC26A9 as a therapeutic target to improve lung function for all people with CF and in individuals with other obstructive lung diseases.

Expression of cystic fibrosis lung disease modifier genes in human airway models.

BACKGROUND: Variation in respiratory response to cystic fibrosis (CF) small molecule therapies is due in part to the contribution of CF lung disease modifier genes. Cultured human bronchial epithelia (HBE) is the gold standard respiratory model for assessing CF therapeutic efficacy but it is hard to access. Cultured human nasal epithelia (HNE) is proposed as a more accessible surrogate model but it is unknown whether the expression profile of the modifier genes are comparable between HNE and HBE which we assess here. METHODS: RNA-sequencing was conducted on paired cultured and fresh HNE and HBE (n = 71 samples) collected from 21 individuals with CF. Genome-wide gene expression was first compared between cultured and fresh cells and then between cultured HNE and HBE based on an equivalence testing procedure we implemented. The co-expression relationships of CFTR and CF lung disease modifier genes were compared between cultured HNE and HBE to determine equivalent interactions. RESULTS: The culturing process had little impact on the expression level of CF lung disease modifier genes. Over 90% of expressed genes showed significant equivalent expression level across cultured HNE and HBE (expression fold-change<2, FDR<0.1), including CFTR and CF lung disease modifier genes. The difference in co-expression relationships among these genes was not significant (p-value=0.99), suggesting their functional interactions are likely to be consistent in the two models. CONCLUSIONS: Cultured HNE recapitulates the expression profile of CF lung disease modifier genes in cultured HBE, suggesting the biological processes involving these genes are likely to be consistent across the two models.

Positive epistasis between disease-causing missense mutations and silent polymorphism with effect on mRNA translation velocity.

Epistasis refers to the dependence of a mutation on other mutation(s) and the genetic context in general. In the context of human disorders, epistasis complicates the spectrum of disease symptoms and has been proposed as a major contributor to variations in disease outcome. The nonadditive relationship between mutations and the lack of complete understanding of the underlying physiological effects limit our ability to predict phenotypic outcome. Here, we report positive epistasis between intragenic mutations in the cystic fibrosis transmembrane conductance regulator (CFTR)-the gene responsible for cystic fibrosis (CF) pathology. We identified a synonymous single-nucleotide polymorphism (sSNP) that is invariant for the CFTR amino acid sequence but inverts translation speed at the affected codon. This sSNP in cis exhibits positive epistatic effects on some CF disease-causing missense mutations. Individually, both mutations alter CFTR structure and function, yet when combined, they lead to enhanced protein expression and activity. The most robust effect was observed when the sSNP was present in combination with missense mutations that, along with the primary amino acid change, also alter the speed of translation at the affected codon. Functional studies revealed that synergistic alteration in ribosomal velocity is the underlying mechanism; alteration of translation speed likely increases the time window for establishing crucial domain-domain interactions that are otherwise perturbed by each individual mutation.

LocusFocus: Web-based colocalization for the annotation and functional follow-up of GWAS.

Genome-wide association studies (GWAS) have primarily identified trait-associated loci in the non-coding genome. Colocalization analyses of SNP associations from GWAS with expression quantitative trait loci (eQTL) evidence enable the generation of hypotheses about responsible mechanism, genes and tissues of origin to guide functional characterization. Here, we present a web-based colocalization browsing and testing tool named LocusFocus (https://locusfocus.research.sickkids.ca). LocusFocus formally tests colocalization using our established Simple Sum method to identify the most relevant genes and tissues for a particular GWAS locus in the presence of high linkage disequilibrium and/or allelic heterogeneity. We demonstrate the utility of LocusFocus, following up on a genome-wide significant locus from a GWAS of meconium ileus (an intestinal obstruction in cystic fibrosis). Using LocusFocus for colocalization analysis with eQTL data suggests variation in ATP12A gene expression in the pancreas rather than intestine is responsible for the GWAS locus. LocusFocus has no operating system dependencies and may be installed in a local web server. LocusFocus is available under the MIT license, with full documentation and source code accessible on GitHub at https://github.com/naim-panjwani/LocusFocus.

Genetic Modifiers of Cystic Fibrosis-Related Diabetes Have Extensive Overlap With Type 2 Diabetes and Related Traits.

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.

Correlating Cystic Fibrosis Transmembrane Conductance Regulator Function with Clinical Features to Inform Precision Treatment of Cystic Fibrosis.

Rationale: The advent of precision treatment for cystic fibrosis using small-molecule therapeutics has created a need to estimate potential clinical improvements attributable to increases in cystic fibrosis transmembrane conductance regulator (CFTR) function. Objectives: To derive CFTR function of a variety of CFTR genotypes and correlate with key clinical features (sweat chloride concentration, pancreatic exocrine status, and lung function) to develop benchmarks for assessing response to CFTR modulators. Methods: CFTR function assigned to 226 unique CFTR genotypes was correlated with the clinical data of 54,671 individuals enrolled in the Clinical and Functional Translation of CFTR (CFTR2) project. Cross-sectional FEV1% predicted measurements were plotted by age at which measurement was obtained. Shifts in sweat chloride concentration and lung function reported in CFTR modulator trials were compared with function-phenotype correlations to assess potential efficacy of therapies. Measurements and Main Results: CFTR genotype function exhibited a logarithmic relationship with each clinical feature. Modest increases in CFTR function related to differing genotypes were associated with clinically relevant improvements in cross-sectional FEV1% predicted over a range of ages (6-82 yr). Therapeutic responses to modulators corresponded closely to predictions from the CFTR2-derived relationship between CFTR genotype function and phenotype. Conclusions: Increasing CFTR function in individuals with severe disease will have a proportionally greater effect on outcomes than similar increases in CFTR function in individuals with mild disease and should reverse a substantial fraction of the disease process. This study provides reference standards for clinical outcomes that may be achieved by increasing CFTR function.

Genetic association and transcriptome integration identify contributing genes and tissues at cystic fibrosis modifier loci.

Cystic Fibrosis (CF) exhibits morbidity in several organs, including progressive lung disease in all patients and intestinal obstruction at birth (meconium ileus) in ~15%. Individuals with the same causal CFTR mutations show variable disease presentation which is partly attributed to modifier genes. With >6,500 participants from the International CF Gene Modifier Consortium, genome-wide association investigation identified a new modifier locus for meconium ileus encompassing ATP12A on chromosome 13 (min p = 3.83x10(-10)); replicated loci encompassing SLC6A14 on chromosome X and SLC26A9 on chromosome 1, (min p<2.2x10(-16), 2.81x10(-11), respectively); and replicated a suggestive locus on chromosome 7 near PRSS1 (min p = 2.55x10(-7)). PRSS1 is exclusively expressed in the exocrine pancreas and was previously associated with non-CF pancreatitis with functional characterization demonstrating impact on PRSS1 gene expression. We thus asked whether the other meconium ileus modifier loci impact gene expression and in which organ. We developed and applied a colocalization framework called the Simple Sum (SS) that integrates regulatory and genetic association information, and also contrasts colocalization evidence across tissues or genes. The associated modifier loci colocalized with expression quantitative trait loci (eQTLs) for ATP12A (p = 3.35x10(-8)), SLC6A14 (p = 1.12x10(-10)) and SLC26A9 (p = 4.48x10(-5)) in the pancreas, even though meconium ileus manifests in the intestine. The meconium ileus susceptibility locus on chromosome X appeared shifted in location from a previously identified locus for CF lung disease severity. Using the SS we integrated the lung disease association locus with eQTLs from nasal epithelia of 63 CF participants and demonstrated evidence of colocalization with airway-specific regulation of SLC6A14 (p = 2.3x10(-4)). Cystic Fibrosis is realizing the promise of personalized medicine, and identification of the contributing organ and understanding of tissue specificity for a gene modifier is essential for the next phase of personalizing therapeutic strategies.

Screening for Regulatory Variants in 460 kb Encompassing the CFTR Locus in Cystic Fibrosis Patients.

It is estimated that up to 5% of cystic fibrosis transmembrane conductance regulator (CFTR) pathogenic alleles are unidentified. Some of these errors may lie in noncoding regions of the locus and affect gene expression. To identify regulatory element variants in the CFTR locus, SureSelect targeted enrichment of 460 kb encompassing the gene was optimized to deep sequence genomic DNA from 80 CF patients with an unequivocal clinical diagnosis but only one or no CFTR-coding region pathogenic variants. Bioinformatics tools were used to identify sequence variants and predict their impact, which were then assayed in transient reporter gene luciferase assays. The effect of five variants in the CFTR promoter and four in an intestinal enhancer of the gene were assayed in relevant cell lines. The initial analysis of sequence data revealed previously known CF-causing variants, validating the robustness of the SureSelect design, and showed that 85 of 160 CF alleles were undefined. Of a total 1737 variants revealed across the extended 460-kb CFTR locus, 51 map to known CFTR cis-regulatory elements, and many of these are predicted to alter transcription factor occupancy. Four promoter variants and all those in the intestinal enhancer significantly repress reporter gene activity. These data suggest that CFTR regulatory elements may harbor novel CF disease-causing variants that warrant further investigation, both for genetic screening protocols and functional assays.

SLC6A14, an amino acid transporter, modifies the primary CF defect in fluid secretion.

The severity of intestinal disease associated with Cystic Fibrosis (CF) is variable in the patient population and this variability is partially conferred by the influence of modifier genes. Genome-wide association studies have identified SLC6A14, an electrogenic amino acid transporter, as a genetic modifier of CF-associated meconium ileus. The purpose of the current work was to determine the biological role of Slc6a14, by disrupting its expression in CF mice bearing the major mutation, F508del. We found that disruption of Slc6a14 worsened the intestinal fluid secretion defect, characteristic of these mice. In vitro studies of mouse intestinal organoids revealed that exacerbation of the primary defect was associated with reduced arginine uptake across the apical membrane, with aberrant nitric oxide and cyclic GMP-mediated regulation of the major CF-causing mutant protein. Together, these studies highlight the role of this apical transporter in modifying cellular nitric oxide levels, residual function of the major CF mutant and potentially, its promise as a therapeutic target.

Improving imputation in disease-relevant regions: lessons from cystic fibrosis.

Does genotype imputation with public reference panels identify variants contributing to disease? Genotype imputation using the 1000 Genomes Project (1KG; 2504 individuals) displayed poor coverage at the causal cystic fibrosis (CF) transmembrane conductance regulator (CFTR) locus for the International CF Gene Modifier Consortium. Imputation with the larger Haplotype Reference Consortium (HRC; 32,470 individuals) displayed improved coverage but low sensitivity of variants clinically relevant for CF. A hybrid reference that combined whole genome sequencing (WGS) from 101 CF individuals with the 1KG imputed a greater number of single-nucleotide variants (SNVs) that would be analyzed in a genetic association study (r2 ≥ 0.3 and MAF ≥ 0.5%) than imputation with the HRC, while the HRC excelled in the lower frequency spectrum. Using the 1KG or HRC as reference panels missed the most common CF-causing variants or displayed low imputation accuracy. Designs that incorporate population-specific WGS can improve imputation accuracy at disease-specific loci, while imputation using public data sets can omit disease-relevant genotypes.

Phenotypic profiling of CFTR modulators in patient-derived respiratory epithelia.

Pulmonary disease is the major cause of morbidity and mortality in patients with cystic fibrosis, a disease caused by mutations in the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene. Heterogeneity in CFTR genotype-phenotype relationships in affected individuals plus the escalation of drug discovery targeting specific mutations highlights the need to develop robust in vitro platforms with which to stratify therapeutic options using relevant tissue. Toward this goal, we adapted a fluorescence plate reader assay of apical CFTR-mediated chloride conductance to enable profiling of a panel of modulators on primary nasal epithelial cultures derived from patients bearing different CFTR mutations. This platform faithfully recapitulated patient-specific responses previously observed in the "gold-standard" but relatively low-throughput Ussing chamber. Moreover, using this approach, we identified a novel strategy with which to augment the response to an approved drug in specific patients. In proof of concept studies, we also validated the use of this platform in measuring drug responses in lung cultures differentiated from cystic fibrosis iPS cells. Taken together, we show that this medium throughput assay of CFTR activity has the potential to stratify cystic fibrosis patient-specific responses to approved drugs and investigational compounds in vitro in primary and iPS cell-derived airway cultures.

Mammographic density defined by higher than conventional brightness thresholds better predicts breast cancer risk.

Background: Mammographic density defined by the conventional pixel brightness threshold, and adjusted for age and body mass index (BMI), is a well-established risk factor for breast cancer. We asked if higher thresholds better separate women with and without breast cancer. Methods: We studied Australian women, 354 with breast cancer over-sampled for early-onset and family history, and 944 unaffected controls frequency-matched for age at mammogram. We measured mammographic dense area and percent density using the CUMULUS software at the conventional threshold, which we call Cumulus , and at two increasingly higher thresholds, which we call Altocumulus and Cirrocumulus , respectively. All measures were Box-Cox transformed and adjusted for age and BMI. We estimated the odds per adjusted standard deviation (OPERA) using logistic regression and the area under the receiver operating characteristic curve (AUC). Results: Altocumulus and Cirrocumulus were correlated with Cumulus (r ∼ 0.8 and 0.6 , respectively) . For dense area, the OPERA was 1.62, 1.74 and 1.73 for Cumulus, Altocumulus and Cirrocumulus , respectively (all P < 0.001). After adjusting for Altocumulus and Cirrocumulus , Cumulus was not significant ( P > 0.6). The OPERAs for percent density were less but gave similar findings. The mean of the standardized adjusted Altocumulus and Cirrocumulus dense area measures was the best predictor; OPERA = 1.87 [95% confidence interval (CI): 1.64-2.14] and AUC = 0.68 (0.65-0.71). Conclusions: The areas of higher mammographically dense regions are associated with almost 30% stronger breast cancer risk gradient, explain the risk association of the conventional measure and might be more aetiologically important. This has substantial implications for clinical translation and molecular, genetic and epidemiological research.

SBDS-Deficient Cells Have an Altered Homeostatic Equilibrium due to Translational Inefficiency Which Explains their Reduced Fitness and Provides a Logical Framework for Intervention.

Ribosomopathies are a family of inherited disorders caused by mutations in genes necessary for ribosomal function. Shwachman-Diamond Bodian Syndrome (SDS) is an autosomal recessive disease caused, in most patients, by mutations of the SBDS gene. SBDS is a protein required for the maturation of 60S ribosomes. SDS patients present exocrine pancreatic insufficiency, neutropenia, chronic infections, and skeletal abnormalities. Later in life, patients are prone to myelodisplastic syndrome and acute myeloid leukemia (AML). It is unknown why patients develop AML and which cellular alterations are directly due to the loss of the SBDS protein. Here we derived mouse embryonic fibroblast lines from an SbdsR126T/R126T mouse model. After their immortalization, we reconstituted them by adding wild type Sbds. We then performed a comprehensive analysis of cellular functions including colony formation, translational and transcriptional RNA-seq, stress and drug sensitivity. We show that: 1. Mutant Sbds causes a reduction in cellular clonogenic capability and oncogene-induced transformation. 2. Mutant Sbds causes a marked increase in immature 60S subunits, limited impact on mRNA specific initiation of translation, but reduced global protein synthesis capability. 3. Chronic loss of SBDS activity leads to a rewiring of gene expression with reduced ribosomal capability, but increased lysosomal and catabolic activity. 4. Consistently with the gene signature, we found that SBDS loss causes a reduction in ATP and lactate levels, and increased susceptibility to DNA damage. Combining our data, we conclude that a cell-specific fragile phenotype occurs when SBDS protein drops below a threshold level, and propose a new interpretation of the disease.

Sources of Variation in Sweat Chloride Measurements in Cystic Fibrosis.

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.

Cystic fibrosis gene modifier SLC26A9 modulates airway response to CFTR-directed therapeutics.

Cystic fibrosis is realizing the promise of personalized medicine. Recent advances in drug development that target the causal CFTR directly result in lung function improvement, but variability in response is demanding better prediction of outcomes to improve management decisions. The genetic modifier SLC26A9 contributes to disease severity in the CF pancreas and intestine at birth and here we assess its relationship with disease severity and therapeutic response in the airways. SLC26A9 association with lung disease was assessed in individuals from the Canadian and French CF Gene Modifier consortia with CFTR-gating mutations and in those homozygous for the common Phe508del mutation. Variability in response to a CFTR-directed therapy attributed to SLC26A9 genotype was assessed in Canadian patients with gating mutations. A primary airway model system determined if SLC26A9 shows modification of Phe508del CFTR function upon treatment with a CFTR corrector. In those with gating mutations that retain cell surface-localized CFTR we show that SLC26A9 modifies lung function while this is not the case in individuals homozygous for Phe508del where cell surface expression is lacking. Treatment response to ivacaftor, which aims to improve CFTR-channel opening probability in patients with gating mutations, shows substantial variability in response, 28% of which can be explained by rs7512462 in SLC26A9 (P = 0.0006). When homozygous Phe508del primary bronchial cells are treated to restore surface CFTR, SLC26A9 likewise modifies treatment response (P = 0.02). Our findings indicate that SLC26A9 airway modification requires CFTR at the cell surface, and that a common variant in SLC26A9 may predict response to CFTR-directed therapeutics.

Prevalence of meconium ileus marks the severity of mutations of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene.

RATIONALE: Meconium ileus (MI) is a perinatal complication in cystic fibrosis (CF), which is only minimally influenced by environmental factors. We derived and examined MI prevalence (MIP) scores to assess CFTR phenotype-phenotype correlation for severe mutations. METHOD: MIP scores were established using a Canadian CF population (n = 2,492) as estimates of the proportion of patients with MI among all patients carrying the same CFTR mutation, focusing on patients with p.F508del as the second allele. Comparisons were made to the registries from the US CF Foundation (n = 43,432), Italy (Veneto/Trentino/Alto Adige regions) (n = 1,788), and Germany (n = 3,596). RESULTS: The prevalence of MI varied among the different registries (13-21%). MI was predominantly prevalent in patients with pancreatic insufficiency carrying "severe" CFTR mutations. In this severe spectrum MIP scores further distinguished between mutation types, for example, G542X (0.31) with a high, F508del (0.22) with a moderate, and G551D (0.08) with a low MIP score. Higher MIP scores were associated with more severe clinical phenotypes, such as a lower forced expiratory volume in 1 second (P = 0.01) and body mass index z score (P = 0.04). CONCLUSIONS: MIP scores can be used to rank CFTR mutations according to their clinical severity and provide a means to expand delineation of CF phenotypes.Genet Med 18 4, 333-340.