Protospacer modification improves base editing of a canonical splice site variant and recovery of CFTR function in human airway epithelial cells

Canonical splice site variants affecting the 5' GT and 3' AG nucleotides of introns result in severe missplicing and account for about 10% of disease-causing genomic alterations. Treatment of such variants has proven challenging due to the unstable mRNA or protein isoforms that typically result from disruption of these sites. Here, we investigate CRISPR-Cas9-mediated adenine base editing for such variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene. We validate a CFTR expression minigene (EMG) system for testing base editing designs for two different targets. We then use the EMG system to test non-standard single-guide RNAs with either shortened or lengthened protospacers to correct the most common cystic fibrosis-causing variant in individuals of African descent (c.2988+1G>A). Varying the spacer region length allowed placement of the editing window in a more efficient context and enabled use of alternate protospacer adjacent motifs. Using these modifications, we restored clinically significant levels of CFTR function to human airway epithelial cells from two donors bearing the c.2988+1G>A variant.

CFTR bearing variant p.Phe312del exhibits function inconsistent with phenotype and negligible response to ivacaftor

The chloride channel dysfunction caused by deleterious cystic fibrosis transmembrane conductance regulator (CFTR) variants generally correlates with severity of cystic fibrosis (CF). However, 3 adults bearing the common severe variant p.Phe508del (legacy: F508del) and a deletion variant in an ivacaftor binding region of CFTR (p.Phe312del; legacy: F312del) manifested only elevated sweat chloride concentration (sw[Cl-]; 87-105 mEq/L). A database review of 25 individuals with F312del and a CF-causing variant revealed elevated sw[Cl-] (75-123 mEq/L) and variable CF features. F312del occurs at a higher-than-expected frequency in the general population, confirming that individuals with F312del and a CF-causing variant do not consistently develop overt CF features. In primary nasal cells, CFTR bearing F312del and F508del generated substantial chloride transport (66.0% ± 4.5% of WT-CFTR) but did not respond to ivacaftor. Single-channel analysis demonstrated that F312del did not affect current flow through CFTR, minimally altered gating, and ablated the ivacaftor response. When expressed stably in CF bronchial epithelial (CFBE41o-) cells, F312del-CFTR demonstrated residual function (50.9% ± 3.3% WT-CFTR) and a subtle decrease in forskolin response compared with WT-CFTR. F312del provides an exception to the established correlation between CFTR chloride transport and CF phenotype and informs our molecular understanding of ivacaftor response.

Evaluation of both exonic and intronic variants for effects on RNA splicing allows for accurate assessment of the effectiveness of precision therapies

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.

Genetic variants that impact pre-mRNA splicing are a common cause of genetic disease and have varying downstream molecular consequences. As a result, precision therapies that function at the protein level are not always effective for these variants and thus careful assessment is necessary. Here we evaluate RNA-level effects of 52 variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene and show that study of splicing and its consequences allows for more accurate assignment of precision therapies.

Increased expression of anion transporter SLC26A9 delays diabetes onset in cystic fibrosis

Diabetes is a common complication of cystic fibrosis (CF) that affects approximately 20% of adolescents and 40%-50% of adults with CF. The age at onset of CF-related diabetes (CFRD) (marked by clinical diagnosis and treatment initiation) is an important measure of the disease process. DNA variants associated with age at onset of CFRD reside in and near SLC26A9. Deep sequencing of the SLC26A9 gene in 762 individuals with CF revealed that 2 common DNA haplotypes formed by the risk variants account for the association with diabetes. Single-cell RNA sequencing (scRNA-Seq) indicated that SLC26A9 is predominantly expressed in pancreatic ductal cells and frequently coexpressed with CF transmembrane conductance regulator (CFTR) along with transcription factors that have binding sites 5' of SLC26A9. These findings were replicated upon reanalysis of scRNA-Seq data from 4 independent studies. DNA fragments derived from the 5' region of SLC26A9-bearing variants from the low-risk haplotype generated 12%-20% higher levels of expression in PANC-1 and CFPAC-1 cells compared with the high- risk haplotype. Taken together, our findings indicate that an increase in SLC26A9 expression in ductal cells of the pancreas delays the age at onset of diabetes, suggesting a CFTR-agnostic treatment for a major complication of CF.