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.

We developed a variant-annotation method that combines sequence-based machine-learning classification with a context-dependent algorithm for selecting splice variants. Our approach is distinctive in that it compares the splice potential of a sequence bearing a variant with the splice potential of the reference sequence. After training, classification accurately identified 168 of 180 (93.3%) canonical splice sites of five genes. The combined method, CryptSplice, identified and correctly predicted the effect of 18 of 21 (86%) known splice-altering variants in CFTR, a well-studied gene whose loss-of-function variants cause cystic fibrosis (CF). Among 1,423 unannotated CFTR disease-associated variants, the method identified 32 potential exonic cryptic splice variants, two of which were experimentally evaluated and confirmed. After complete CFTR sequencing, the method found three cryptic intronic splice variants (one known and two experimentally verified) that completed the molecular diagnosis of CF in 6 of 14 individuals. CryptSplice interrogation of sequence data from six individuals with X-linked dyskeratosis congenita caused by an unknown disease-causing variant in DKC1 identified two splice-altering variants that were experimentally verified. To assess the extent to which disease-associated variants might activate cryptic splicing, we selected 458 pathogenic variants and 348 variants of uncertain significance (VUSs) classified as high confidence from ClinVar. Splice-site activation was predicted for 129 (28%) of the pathogenic variants and 75 (22%) of the VUSs. Our findings suggest that cryptic splice-site activation is more common than previously thought and should be routinely considered for all variants within the transcribed regions of genes.

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.

Elevated sweat chloride levels, failure to thrive (FTT), and lung disease are characteristic features of cystic fibrosis (CF, OMIM #219700). Here we describe variants in CA12 encoding carbonic anhydrase XII in two pedigrees exhibiting CF-like phenotypes. Exome sequencing of a white American adult diagnosed with CF due to elevated sweat chloride, recurrent hyponatremia, infantile FTT and lung disease identified deleterious variants in each CA12 gene: c.908-1 G>A in a splice acceptor and a novel frameshift insertion c.859_860insACCT. In an unrelated consanguineous Omani family, two children with elevated sweat chloride, infantile FTT, and recurrent hyponatremia were homozygous for a novel missense variant (p.His121Gln). Deleterious CFTR variants were absent in both pedigrees. CA XII protein was localized apically in human bronchiolar epithelia and basolaterally in the reabsorptive duct of human sweat glands. Respiratory epithelial cell RNA from the adult proband revealed only aberrant CA12 transcripts and in vitro analysis showed greatly reduced CA XII protein. Studies of ion transport across respiratory epithelial cells in vivo and in culture revealed intact CFTR-mediated chloride transport in the adult proband. CA XII protein bearing either p.His121Gln or a previously identified p.Glu143Lys missense variant localized to the basolateral membranes of polarized Madin-Darby canine kidney (MDCK) cells, but enzyme activity was severely diminished when assayed at physiologic concentrations of extracellular chloride. Our findings indicate that loss of CA XII function should be considered in individuals without CFTR mutations who exhibit CF-like features in the sweat gland and lung.

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.

Elucidating the functional consequence of molecular defects underlying genetic diseases enables appropriate design of therapeutic options. Treatment of cystic fibrosis (CF) is an exemplar of this paradigm as the development of CFTR modulator therapies has allowed for targeted and effective treatment of individuals harboring specific genetic variants. However, the mechanism of these drugs limits effectiveness to particular classes of variants that allow production of CFTR protein. Thus, assessment of the molecular mechanism of individual variants is imperative for proper assignment of these precision therapies. This is particularly important when considering variants that affect pre-mRNA splicing, thus limiting success of the existing protein-targeted therapies. Variants affecting splicing can occur throughout exons and introns and the complexity of the process of splicing lends itself to a variety of outcomes, both at the RNA and protein levels, further complicating assessment of disease liability and modulator response. To investigate the scope of this challenge, we evaluated splicing and downstream effects of 52 naturally occurring CFTR variants (exonic = 15, intronic = 37). Expression of constructs containing select CFTR intronic sequences and complete CFTR exonic sequences in cell line models allowed for assessment of RNA and protein-level effects on an allele by allele basis. Characterization of primary nasal epithelial cells obtained from individuals harboring splice variants corroborated in vitro data. Notably, we identified exonic variants that result in complete missplicing and thus a lack of modulator response (e.g. c.2908G>A, c.523A>G), as well as intronic variants that respond to modulators due to the presence of residual normally spliced transcript (e.g. c.4242+2T>C, c.3717+40A>G). Overall, our data reveals diverse molecular outcomes amongst both exonic and intronic variants emphasizing the need to delineate RNA, protein, and functional effects of each variant in order to accurately assign precision therapies.