Treatment of individuals with cystic fibrosis (CF) has been transformed by small molecule therapies that target select pathogenic variants in the CF transmembrane conductance regulator (CFTR). To expand treatment eligibility, we stably expressed 43 rare missense CFTR variants associated with moderate CF from a single site in the genome of human CF bronchial epithelial (CFBE41o-) cells. The magnitude of drug response was highly correlated with residual CFTR function for the potentiator ivacaftor, the corrector lumacaftor, and ivacaftor-lumacaftor combination therapy. Response of a second set of 16 variants expressed stably in Fischer rat thyroid (FRT) cells showed nearly identical correlations. Subsets of variants were identified that demonstrated statistically significantly higher responses to specific treatments. Furthermore, nearly all variants studied in CFBE cells (40 of 43) and FRT cells (13 of 16) demonstrated greater response to ivacaftor-lumacaftor combination therapy than either modulator alone. Together, these variants represent 87% of individuals in the CFTR2 database with at least 1 missense variant. Thus, our results indicate that most individuals with CF carrying missense variants are (a) likely to respond modestly to currently available modulator therapy, while a small fraction will have pronounced responses, and (b) likely to derive the greatest benefit from combination therapy.

Missense DNA variants have variable effects upon protein function. Consequently, interpreting their pathogenicity is challenging, especially when they are associated with disease variability. To determine the degree to which functional assays inform interpretation, we analyzed 48 CFTR missense variants associated with variable expressivity of cystic fibrosis (CF). We assessed function in a native isogenic context by evaluating CFTR mutants that were stably expressed in the genome of a human airway cell line devoid of endogenous CFTR expression. 21 of 29 variants associated with full expressivity of the CF phenotype generated <10% wild-type CFTR (WT-CFTR) function, a conservative threshold for the development of life-limiting CF lung disease, and five variants had moderately decreased function (10% to ∼25% WT-CFTR). The remaining three variants in this group unexpectedly had >25% WT-CFTR function; two were higher than 75% WT-CFTR. As expected, 14 of 19 variants associated with partial expressivity of CF had >25% WT-CFTR function; however, four had minimal to no effect on CFTR function (>75% WT-CFTR). Thus, 6 of 48 (13%) missense variants believed to be disease causing did not alter CFTR function. Functional studies substantially refined pathogenicity assignment with expert annotation and criteria from the American College of Medical Genetics and Genomics and Association for Molecular Pathology. However, four algorithms (CADD, REVEL, SIFT, and PolyPhen-2) could not differentiate between variants that caused severe, moderate, or minimal reduction in function. In the setting of variable expressivity, these results indicate that functional assays are essential for accurate interpretation of missense variants and that current prediction tools should be used with caution.

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. Crosssectional FEV 1 % 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 FEV 1 % 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.

CFTR modulators have revolutionized the treatment of individuals with cystic fibrosis (CF) by improving the function of existing protein. Unfortunately, almost half of the disease-causing variants in CFTR are predicted to introduce premature termination codons (PTC) thereby causing absence of full-length CFTR protein. We hypothesized that a subset of nonsense and frameshift variants in CFTR allow expression of truncated protein that might respond to FDA-approved CFTR modulators. To address this concept, we selected 26 PTC-generating variants from four regions of CFTR and determined their consequences on CFTR mRNA, protein and function using intron-containing minigenes expressed in 3 cell lines (HEK293, MDCK and CFBE41o-) and patient-derived conditionally reprogrammed primary nasal epithelial cells. The PTC-generating variants fell into five groups based on RNA and protein effects. Group A (reduced mRNA, immature (core glycosylated) protein, function <1% (n = 5)) and Group B (normal mRNA, immature protein, function <1% (n = 10)) variants were unresponsive to modulator treatment. However, Group C (normal mRNA, mature (fully glycosylated) protein, function >1% (n = 5)), Group D (reduced mRNA, mature protein, function >1% (n = 5)) and Group E (aberrant RNA splicing, mature protein, function > 1% (n = 1)) variants responded to modulators. Increasing mRNA level by inhibition of NMD led to a significant amplification of modulator effect upon a Group D variant while response of a Group A variant was unaltered. Our work shows that PTC-generating variants should not be generalized as genetic ‘nulls’ as some may allow generation of protein that can be targeted to achieve clinical benefit.