Theratyping cystic fibrosis in vitro in ALI-culture and organoid models generated from patient-derived nasal epithelial Conditionally Reprogrammed Stem Cells

QUESTION

Cystic Fibrosis (CF) is due to pathogenic variants in the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene. Recent improvement enabled pharmacologic therapy aiming at restoring mutated CFTR expression and function. CFTR "modulators" have revolutionised the CF therapeutic landscape, particularly the last approved Trikafta. This drug-combination is indicated by FDA and very recently by EMA for genotypes carrying at least one copy of CFTR with F508del pathogenic variant. However, several genotypes, are not eligible for Trikafta treatment, yet.

MATERIALS/PATIENTS AND METHODS

We exploited an innovative cellular approach allowing highly efficient in vitro-expansion of airway epithelial stem cells (AESC) through conditional reprogramming (CRC) from nasal brushing of CF patients. This approach, coupled to development of AESC-derived personalised disease models, as organoids and air liquid interface (ALI) cultures, revealed highly suitable for CFTR pharmacological-testing.

RESULTS AND ANSWER TO THE QUESTION

We fully validated the experimental models and implemented the CFTR functional assays and biochemical CFTR protein characterisation, that allowed to evaluate the efficacy of clinically available modulators in restoring CFTR maturation and function of each patient-derived "avatar" (theratyping). F508del homozygous genotypes, used as controls, confirmed the higher clinical activity of Trikafta in comparison with older modulators. Trikafta showed its efficacy also on three rare genotypes previously not eligible for modulators-treatment, opening the way to clinical translation. Finally, encouraging results for innovative drug combinations were also obtained.

A new complex allele of the CFTR gene partially explains the variable phenotype of the L997F mutation

PURPOSE

To evaluate the role of complex alleles, with two or more mutations in cis position, of the cystic fibrosis transmembrane conductance regulator (CFTR) gene in the definition of the genotype-phenotype relationship in cystic fibrosis (CF), and to evaluate the functional significance of the highly controversial L997F CFTR mutation.

METHODS

We evaluated the diagnosis of CF or CFTR-related disorders in 12 unrelated subjects with highly variable phenotypes. According to a first CFTR mutational analysis, subjects appeared to be compound heterozygotes for a classic mutation and the L997F mutation. A further CFTR mutational analysis was conducted by means of a protocol of extended sequencing, particularly suited to the detection of complex alleles.

RESULTS

We detected a new [R117L; L997F] CFTR complex allele in the four subjects with the highest sweat test values and CF. The eight subjects without the complex allele showed the most varied biochemical and clinical outcome and were diagnosed as having mild CF, CFTR-related disorders, or even no disease.

CONCLUSIONS

The new complex allele partially explains the variable phenotype in CF subjects with the L997F mutation. CFTR complex alleles are likely to have a role in the definition of the genotype-phenotype relationship in CF. Whenever apparently identical CFTR-mutated genotypes are found in subjects with divergent phenotypes, an extensive mutational search is mandatory.

A Genotypic-Oriented View of CFTR Genetics Highlights Specific Mutational Patterns Underlying Clinical Macrocategories of Cystic Fibrosis

Cystic fibrosis (CF) is a monogenic disease caused by mutations of the cystic fibrosis transmembrane conductance regulator (CFTR) gene. The genotype-phenotype relationship in this disease is still unclear, and diagnostic, prognostic and therapeutic challenges persist. We enrolled 610 patients with different forms of CF and studied them from a clinical, biochemical, microbiological and genetic point of view. Overall, there were 125 different mutated alleles (11 with novel mutations and 10 with complex mutations) and 225 genotypes. A strong correlation between mutational patterns at the genotypic level and phenotypic macrocategories emerged. This specificity appears to largely depend on rare and individual mutations, as well as on the varying prevalence of common alleles in different clinical macrocategories. However, 19 genotypes appeared to underlie different clinical forms of the disease. The dissection of the pathway from the CFTR mutated genotype to the clinical phenotype allowed to identify at least two components of the variability usually found in the genotype-phenotype relationship. One component seems to depend on the genetic variation of CFTR, the other component on the cumulative effect of variations in other genes and cellular pathways independent from CFTR. The experimental dissection of the overall biological CFTR pathway appears to be a powerful approach for a better comprehension of the genotype-phenotype relationship. However, a change from an allele-oriented to a genotypic-oriented view of CFTR genetics is mandatory, as well as a better assessment of sources of variability within the CFTR pathway.

DAT gene polymorphisms (rs28363170, rs393795) and levodopa-induced dyskinesias in Parkinson's disease

L-dopa-induced dyskinesias (LID) is a common motor side effect of levodopa therapy of Parkinson's disease (PD). The identified predictors may only partially account for the risk of developing LID and genetic factors may contribute to this variability. The present study is aimed to investigate whether polymorphisms in the dopamine transporter gene (DAT) are associated with the risk of developing LID. Genotyping of the 40-bp VNTR (rs28363170) and rs393795 (A/C) polymorphisms of the DAT gene was performed in a well-characterized cohort of 181 Italian PD patients in treatment with L-DOPA for 3 years or more. The results of our study show that there is no difference in dyskinesias prevalence among carriers of the two DAT gene polymorphisms. However, the combination of the two genotypes 10R/10R (rs28363170) and A carrier (rs393795) of the DAT gene reduces the risk of LID occurrence during long-term therapy with l-DOPA with respect to the PD subjects who did not carry these alleles (OR = 0.31; 95% CI, 0.09-0.88). Also based on a logistic regression analysis, the 10R/10R and the A carrier allele of the rs393795 polymorphisms of the DAT gene, could reduce the susceptibility to develop LID during levodopa therapy adjusted by demographical and clinical variables (OR = 0.19; 95% CI, 0.05-0.69). Additional studies further investigating the rs28363170 and rs393795 polymorphisms with LID in PD are needed to clarify their role in different ethnicities.

DNA Methylation Patterns Correlate with the Expression of SCNN1A, SCNN1B, and SCNN1G (Epithelial Sodium Channel, ENaC) Genes

The interplay between the cystic fibrosis transmembrane conductance regulator (CFTR) and the epithelial sodium channel (ENaC) in respiratory epithelia has a crucial role in the pathogenesis of cystic fibrosis (CF). The comprehension of the mechanisms of transcriptional regulation of ENaC genes is pivotal to better detail the pathogenic mechanism and the genotype-phenotype relationship in CF, as well as to realize therapeutic approaches based on the transcriptional downregulation of ENaC genes. Since we aimed to study the epigenetic transcriptional control of ENaC genes, an assessment of their expression and DNA methylation patterns in different human cell lines, nasal brushing samples, and leucocytes was performed. The mRNA expression of CFTR and ENaC subunits α, β and γ (respectively SCNN1A, SCNN1B, and SCNN1G genes) was studied by real time PCR. DNA methylation of 5'-flanking region of SCNN1A, SCNN1B, and SCNN1G genes was studied by HpaII/PCR. The levels of expression and DNA methylation of ENaC genes in the different cell lines, brushing samples, and leukocytes were very variable. The DNA regions studied of each ENaC gene showed different methylation patterns. A general inverse correlation between expression and DNA methylation was evidenced. Leukocytes showed very low expression of all the 3 ENaC genes corresponding to a DNA methylated pattern. The SCNN1A gene resulted to be the most expressed in some cell lines that, accordingly, showed a completely demethylated pattern. Coherently, a heavy and moderate methylated pattern of, respectively, SCNN1B and SCNN1G genes corresponded to low levels of expression. As exceptions, we found that dexamethasone treatment appeared to stimulate the expression of all the 3 ENaC genes, without an evident modulation of the DNA methylation pattern, and that in nasal brushing a considerable expression of all the 3 ENaC genes were found despite an apparent methylated pattern. At least part of the expression modulation of ENaC genes seems to depend on the DNA methylation patterns of specific DNA regions. This points to epigenetics as a controlling mechanism of ENaC function and as a possible therapeutic approach for CF.

A template for mutational data analysis of the CFTR gene

BACKGROUND

Automated DNA sequencing produces large amounts of data that need to be analyzed by appropriate software. Personalization of software can be a difficult and time-consuming task, especially if a large number of mutations have to be analyzed.

METHODS

The Applied BioSystems SeqScape software, based on the KB basecaller algorithm, is a versatile tool that can be used for mutational analysis and for data quality assessment of sequences belonging to any gene of interest. Using this software we analyzed over 1400 sequences of CFTR exons and adjacent intronic zones, representing over 500,000 bases.

RESULTS

We present an up to date specific template and a linked set of instructions for automated labeling of all point mutations and polymorphisms of the CFTR gene, whose mutations cause cystic fibrosis (the most common genetic disease among Caucasian individuals). We also describe our refined software settings for mutational analysis, in order to keep to a minimum the need of manual validation.

CONCLUSIONS

The use of our template greatly simplifies the mutational analysis of the CFTR gene, reducing human intervention. In our opinion, it might not only be useful to researchers that already perform CFTR mutational analysis by sequencing methods but it should also improve the approach in those laboratories that already use ABI PRISM instrumentation for a limited mutational analysis of the CFTR gene. Similar mutational templates can also be used for other disease causing genes, thus improving molecular genetics protocols.

CFTR expression and activity from the human CFTR locus in BAC vectors, with regulatory regions, isolated by a single-step procedure

We have assembled two BAC vectors containing a single fragment spanning the entire CFTR locus and including the upstream and downstream regions. The two vectors differ in size of the upstream region, and were recovered in Escherichia coli, with intact BAC DNAs prepared for structural and functional analyses. Sequence analysis allowed precise mapping of the inserts. We show that the CFTR gene was wild type and is categorized as the most frequent haplotype in Caucasian populations, identified by the following polymorphisms: (GATT)₇ in intron 6a; (TG)₁₁T₇ in intron 8; V470 at position 470. CFTR expression and activity were analyzed in model cells by RT-PCR, quantitative real-time PCR, western blotting, indirect immunofluorescence and electrophysiological methods, which show the presence of an active CFTR Cl ⁻ channel. Finally, and supporting the hypothesis that CFTR functions as a receptor for Pseudomonas aeruginosa, we show that CFTR-expressing cells internalized more bacteria than parental cells that do not express CFTR. Overall, these data demonstrate that the BAC vectors contain a functional CFTR fragment and have unique features, including derivation from a single fragment, availability of a detailed genomic map and the possibility to use standard extraction procedures for BAC DNA preparations.

Downregulation of epithelial sodium channel (ENaC) activity in cystic fibrosis cells by epigenetic targeting

The pathogenic mechanism of cystic fibrosis (CF) includes the functional interaction of the cystic fibrosis transmembrane conductance regulator (CFTR) protein with the epithelial sodium channel (ENaC). The reduction of ENaC activity may constitute a therapeutic option for CF. This hypothesis was evaluated using drugs that target the protease-dependent activation of the ENaC channel and the transcriptional activity of its coding genes. To this aim we used: camostat, a protease inhibitor; S-adenosyl methionine (SAM), showed to induce DNA hypermethylation; curcumin, known to produce chromatin condensation. SAM and camostat are drugs already clinically used in other pathologies, while curcumin is a common dietary compound. The experimental systems used were CF and non-CF immortalized human bronchial epithelial cell lines as well as human bronchial primary epithelial cells. ENaC activity and SCNN1A, SCNN1B and SCNN1G gene expression were analyzed, in addition to SCNN1B promoter methylation. In both immortalized and primary cells, the inhibition of extracellular peptidases and the epigenetic manipulations reduced ENaC activity. Notably, the reduction in primary cells was much more effective. The SCNN1B appeared to be the best target to reduce ENaC activity, in respect to SCNN1A and SCNN1G. Indeed, SAM treatment resulted to be effective in inducing hypermethylation of SCNN1B gene promoter and in lowering its expression. Importantly, CFTR expression was unaffected, or even upregulated, after treatments. These results open the possibility of CF patients’ treatment by epigenetic targeting.

Two novel and correlated CF-causing insertions in the (TG)mTn tract of the CFTR gene

Two novel and related pathogenic variants of the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene were structurally and functionally characterized. These alterations have not been previously described in literature. Two patients with diagnosis of Cystic Fibrosis (CF) based on the presence of one mutated allele, p.Phe508del, pathological sweat test and clinical symptoms were studied. To complete the genotypes of both patients, an extensive genetic and functional analysis of the CFTR gene was performed. Extensive genetic characterization confirmed the presence of p.Phe508del pathogenic variant and revealed, in both patients, the presence of an insertion of part of intron 10 in intron 9 of the CFTR gene, within the (TG)m repeat, with a variable poly-T stretch. The molecular lesions resulted to be very similar in both patients, with only a difference in the number of T in the poly-T stretch. The functional characterization at RNA level revealed a complete anomalous splicing, without exon 10, from the allele with the insertion of both patients. Consequently, the alleles with the insertions are expected not to contribute to the formation of a functional CFTR protein. Molecular and functional features of these alterations are compatible with the definition of novel CF-causing variants of the CFTR gene. This also allowed the completion of the genetic characterization of both patients.

Quantitative Evaluation of CFTR Gene Expression: A Comparison between Relative Quantification by Real-Time PCR and Absolute Quantification by Droplet Digital PCR

In the precision medicine era of cystic fibrosis (CF), therapeutic interventions, by the so-called modulators, target the cystic fibrosis transmembrane conductance regulator (CFTR) protein. The levels of targetable CFTR proteins are a main variable in the success of patient-specific therapy. In turn, the CFTR protein level depends, at least in part, on the level of CFTR mRNA. Many mechanisms can modulate the CFTR mRNA level, for example, transcriptional rate, stability of the mRNA, epigenetics, and pathogenic variants that can affect mRNA production and degradation. Independently from the causes of variable CFTR mRNA levels, their exact quantitative assessment is of great importance in CF. Methods with high analytical sensitivity, precision, and accuracy are mandatory for the quantitative evaluation aimed at the amelioration of the diagnostic, prognostic, and therapeutic aspects. This paper compares, for the first time, two CFTR gene expression quantification methods: a well-established method for the relative quantification of CFTR mRNA using a real-time PCR and an innovative method for its absolute quantification using a droplet digital PCR. No comprehensive methods for absolute CFTR quantification via droplet digital PCR have been published so far. The accurate quantification of CFTR expression at the mRNA level is a critical step for the personalized therapeutic approaches of CF.