Deep resequencing of CFTR in 762 F508del homozygotes reveals clusters of non-coding variants associated with cystic fibrosis disease traits

Estimation of Chloride Channel Residual Function and Assessment of Targeted Drugs Efficiency in the Presence of a Complex Allele [L467F;F508del] in the CFTR Gene

Complex alleles of the CFTR gene complicate the diagnosis of cystic fibrosis (CF), the classification of its pathogenic variants, affect the clinical picture of the disease and can affect the efficiency of targeted drugs. The total frequency of complex allele [L467F;F508del] in the Russian population of patients with CF is 0.74%, and in patients with the F508del/F508del genotype, its frequency reaches 8%. This article presents multi-faceted study of the complex allele [L467F;F508del] in a cohort of patients with genotypes [L467F;F508del]/class I (c.3532_3535dup, c.1766+2T>C, W1310X, 712-1G>T), and data for a unique patient with the genotype [L467F;F508del]/[L467F;F508del]. Using the intestinal current measurement method, it was demonstrated the absence of CFTR function for [L467F;F508del]/class I and [L467F;F508del]/[L467F;F508del] genotypes. In intestinal organoids, it was shown that [L467F;F508del] in combination with class I variants and in the homozygotes abolishes the efficacy of both two-component (ivacaftor+lumacaftor; ivacaftor+tezacaftor) and three-component (ivacaftor+tezacaftor+elexacaftor) targeted drugs. When prescribing ivacaftor+tezacaftor+elexacaftor to three patients, they did not have a clinical effect after 6-12 months.

Mutation accumulation in H. sapiens F508del CFTR countermands dN/dS type genomic analysis

Understanding the mechanisms that underlie de novo mutations (DNMs) can be essential for interpreting human evolution, including aspects such as rapidly diverging genes, conservation of non-coding regulatory elements, and somatic DNA adaptation, among others. DNM accumulation in Homo sapiens is often limited to evaluation of human trios or quads across a single generation. Moreover, human SNPs in exons, pseudogenes, or other non-coding elements can be ancient and difficult to date, including polymorphisms attributable to founder effects and identity by descent. In this report, we describe multigenerational evolution of a human coding locus devoid of natural selection, and delineate patterns and principles by which DNMs have accumulated over the past few thousand years. We apply a data set comprising cystic fibrosis transmembrane conductance regulator (CFTR) alleles from 2,393 individuals homozygous for the F508del defect. Additional polymorphism on the F508del background diversified subsequent to a single mutational event during recent human history. Because F508del CFTR is without function, SNPs observed on this haplotype are effectively attributable to factors that govern accumulating de novo mutations. We show profound enhancement of transition, synonymous, and positionally repetitive polymorphisms, indicating appearance of DNMs in a manner evolutionarily designed to protect protein coding DNA against mutational attrition while promoting diversity.

Small-molecule correctors divert CFTR-F508del from ERAD by stabilizing sequential folding states

Over 80% of people with cystic fibrosis (CF) carry the F508del mutation in the cystic fibrosis transmembrane conductance regulator (CFTR), a chloride ion channel at the apical plasma membrane (PM) of epithelial cells. F508del impairs CFTR folding causing it to be destroyed by endoplasmic reticulum associated degradation (ERAD). Small-molecule correctors, which act as pharmacological chaperones to divert CFTR-F508del from ERAD, are the primary strategy for treating CF, yet corrector development continues with only a rudimentary understanding of how ERAD targets CFTR-F508del. We conducted genome-wide CRISPR/Cas9 knockout screens to systematically identify the molecular machinery that underlies CFTR-F508del ERAD. Although the ER-resident ubiquitin ligase, RNF5 was the top E3 hit, knocking out RNF5 only modestly reduced CFTR-F508del degradation. Sublibrary screens in an RNF5 knockout background identified RNF185 as a redundant ligase and demonstrated that CFTR-F508del ERAD is robust. Gene-drug interaction experiments illustrated that correctors tezacaftor (VX-661) and elexacaftor (VX-445) stabilize sequential, RNF5-resistant folding states. We propose that binding of correctors to nascent CFTR-F508del alters its folding landscape by stabilizing folding states that are not substrates for RNF5-mediated ubiquitylation.

Finding Needles in the Haystack: Strategies for Uncovering Noncoding Regulatory Variants

Despite accumulating evidence implicating noncoding variants in human diseases, unraveling their functionality remains a significant challenge. Systematic annotations of the regulatory landscape and the growth of sequence variant data sets have fueled the development of tools and methods to identify causal noncoding variants and evaluate their regulatory effects. Here, we review the latest advances in the field and discuss potential future research avenues to gain a more in-depth understanding of noncoding regulatory variants.

Small molecule correctors divert CFTR-F508del from ERAD by stabilizing sequential folding states

Over 80% of people with cystic fibrosis (CF) carry the F508del mutation in the cystic fibrosis transmembrane conductance regulator (CFTR), a chloride ion channel at the apical plasma membrane (PM) of epithelial cells. F508del impairs CFTR folding causing it to be destroyed by endoplasmic reticulum associated degradation (ERAD). Small molecule correctors, which act as pharmacological chaperones to divert CFTR-F508del from ERAD, are the primary strategy for treating CF, yet corrector development continues with only a rudimentary understanding of how ERAD targets CFTR-F508del. We conducted genome-wide CRISPR/Cas9 knockout screens to systematically identify the molecular machinery that underlies CFTR-F508del ERAD. Although the ER-resident ubiquitin ligase, RNF5 was the top E3 hit, knocking out RNF5 only modestly reduced CFTR-F508del degradation. Sublibrary screens in an RNF5 knockout background identified RNF185 as a redundant ligase, demonstrating that CFTR-F508del ERAD is highly buffered. Gene-drug interaction experiments demonstrated that correctors tezacaftor (VX-661) and elexacaftor (VX-445) stabilize sequential, RNF5-resistant folding states. We propose that binding of correctors to nascent CFTR-F508del alters its folding landscape by stabilizing folding states that are not substrates for RNF5-mediated ubiquitylation.

Tissue-Specific Regulation of CFTR Gene Expression

More than 2000 variations are described within the CFTR (Cystic Fibrosis Transmembrane Regulator) gene and related to large clinical issues from cystic fibrosis to mono-organ diseases. Although these CFTR-associated diseases have been well documented, a large phenotype spectrum is observed and correlations between phenotypes and genotypes are still not well established. To address this issue, we present several regulatory elements that can modulate CFTR gene expression in a tissue-specific manner. Among them, cis-regulatory elements act through chromatin loopings and take part in three-dimensional structured organization. With tissue-specific transcription factors, they form chromatin modules and can regulate gene expression. Alterations of specific regulations can impact and modulate disease expressions. Understanding all those mechanisms highlights the need to expand research outside the gene to enhance our knowledge.

Targeted locus amplification reveals heterogeneity between and within CFTR genotypes and association with CFTR function in patient-derived intestinal organoids

BACKGROUND

Cystic fibrosis (CF) disease severity can be highly variable, even between people with CF (pwCF) with similar genotypes. Here we use patient-derived intestinal organoids to study the influence of genetic variation within the cystic fibrosis transmembrane conductance regulator (CFTR) gene on CFTR function.

METHODS

Organoids of F508del/class I, F508del/S1251N and pwCF with only one detected CF-causing mutation were cultured. Allele-specific CFTR variation was investigated using targeted locus amplification (TLA), CFTR function was measured using the forskolin-induced swelling assay and mRNA levels were quantified using RT-qPCR.

RESULTS

We were able to distinguish CFTR genotypes based on TLA data. Additionally, we observed heterogeneity within genotypes, which we were able to link to CFTR function for S1251N alleles.

CONCLUSIONS

Our results indicate that the paired analysis of CFTR intragenic variation and CFTR function can gain insights in the underlying CFTR defect for individuals where the disease phenotype does not match the CFTR mutations detected during diagnosis.

Personalized Selection of a CFTR Modulator for a Patient with a Complex Allele [L467F;F508del]

The presence of complex alleles in the CFTR gene can lead to difficulties in diagnosing cystic fibrosis and cause resistance to therapy with CFTR modulators. Tezacaftor/ivacaftor therapy for 8 months in a patient with the initially established F508del/F508del genotype did not lead to an improvement in her condition—there was no change in spirometry and an increase in the patient’s weight, while there was only a slight decrease in NaCl values, measured by a sweat test. The intestinal current measurements of the patient’s rectal biopsy showed no positive dynamics in the rescue of CFTR function while taking tezacaftor/ivacaftor. The assumption that the patient had an additional mutation in the cis position was confirmed by sequencing the CFTR gene, and the complex allele [L467F;F508del] was identified. Based on the rescue of CFTR function by elexacaftor/tezacaftor/ivacaftor obtained using forskolin-induced swelling on intestinal organoids, the patient was prescribed therapy with this targeted drug. The use of elexacaftor/tezacaftor/ivacaftor for 7 months resulted in a significant improvement in the patient’s clinical condition.

Evaluation of the Complex p.[Leu467Phe;Phe508del] CFTR Allele in the Intestinal Organoids Model: Implications for Therapy

In the cohort of Russian patients with cystic fibrosis, the p.[Leu467Phe;Phe508del] complex allele (legacy name [L467F;F508del]) of the CFTR gene is understudied. In this research, we present the results of frequency evaluation of the [L467F;F508del] complex allele in the Russian Federation among patients with a F508del/F508del genotype, its effect on the clinical course of cystic fibrosis, the intestinal epithelium ionic channel function, and the effectiveness of target therapy. The frequency of the [L467F;F508del] complex allele among patients with homozygous F508del was determined with multiplex ligase-dependent probe amplification followed by polymerase chain reaction and fragment analysis. The function of ionic channels, including the residual CFTR function, and the effectiveness of CFTR modulators was analyzed using intestinal current measurements on rectal biopsy samples and the forskolin-induced swelling assay on organoids. The results showed that the F508del/[L467F;F508del] genotype is present in 8.2% of all Russian patients with F508del in a homozygous state. The clinical course of the disease in patients with the F508del/[L467F;F508del] genotype is severe and does not vary from the course in the cohort with homozygous F508del, although the CFTR channel function is significantly lower. For patients with the F508del/[L467F;F508del] genotype, we can recommend targeted therapy using a combined ivacaftor + tezacaftor + elexacaftor medication.

Modulator Therapy in Cystic Fibrosis Patients with cis Variants in F508del Complex Allele: A Short-Term Observational Case Series

Previous studies reported the influence of cis variants in F508del cystic fibrosis (CF) patients in their responses to CFTR modulators. The current study is a prospective, observational study involving three patients with CF and pancreatic insufficiency, carrying a complex allele including F508del with A238V, I1027T, or L467F. We report clinical data before and after 4 weeks of treatment with tezacaftor (TEZ)/ivacaftor (IVA), elexacaftor (ELX)/TEZ/IVA, and lumacaftor (LUM)/IVA for patients with complex alleles A238V, I1027T, and L467F, respectively. The 50-year-old patient bearing F508del;A238V/D1152H showed a normal sweat test (13 mEq/L) and improvements in forced expiratory volume in the first second (FEV₁) (+7 points), body mass index (BMI) (+0.85), and respiratory CF Questionnaire-Revised (CFQ-R) domain (+22.2 points). The 12-year-old patient bearing F508del;I1027T/R709X showed an improvement in a sweat test (−40 mEq/l), FEV₁ (+9 points) and the respiratory CFQ-R domain (+16.7 points). No changes in outcomes were observed for the 6-year-old patient F508del;L467F/F508del. Our data highlight that the reported variants do not modify the phenotypic expression of F508del. Searching L467F is crucial in CF patients with F508del nonresponsive to ELX/TEZ/IVA. Further data are needed to evaluate the clinical effect of these variants after a longer follow up.

High frequency of complex CFTR alleles associated with c.1521_1523delCTT (F508del) in Russian cystic fibrosis patients

Cystic fibrosis (CF, MIM# 219,700) is an autosomal recessive disease caused by pathogenic variants within the CFTR gene. It was shown that genetic variants located in cis can affect disease severity or treatment response because of additive or epistatic effects. Studies on the prevalence of complex alleles in Russian CF patients have just begun.

Aim

To evaluate frequencies and genetic background of complex alleles carrying c.1521_1523delCTT (F508del) and c.1399C>T (L467F), c.2562T>G (T854=) or c.4389G>A (Q1463=) in cis; to determine clinical consequences of complex allele c.[1399C>T;1521_1523delCTT] ([L467;F508del]) in Russian CF patients.

Methods

Sequencing of coding regions of CFTR gene and analysis of polymorphic markers in CF patients carrying F508del variant. Comparing of clinical features in two groups patients having genotypes [L467F;F508del];[F508del] (group 1) and [F508del];[F508del] (group 2).

Results

Frequency of [L467F;F508del] allele linked to 2–2–21–6–17–13 haplotype was 4.42%, of [F508del;T854=;Q1463=] allele linked to haplotype 1–2–21–6–17–13 – 2.2% in F508del chromosomes. No differences in disease severity in patients carrying complex allele [L467F;F508del] and patients homozygous for F508del was found.

Conclusion

The frequency of complex alleles associated with F508del was at least 6.6% in Russian CF patients, which should be taken into account for the decision on optimal treatment options with CFTR modulators.

Systematic estimation of cystic fibrosis prevalence in Chinese and genetic spectrum comparison to Caucasians

Background

Cystic fibrosis (CF) is a common, life-threatening genetic disease in Caucasians but rarely reported in Chinese population. The prevalence and population-specific genetic spectrum of CF in China needs to be systematically estimated and compared with Caucasians.

Materials and methods

We reviewed 30,951 exome-sequencing samples, including 20,909 pediatric patient samples and 10,042 parent samples, from Chinese Children's Rare Disease Genetic Testing Clinical Collaboration System (CCGT). After the in-lab filtration process, 477 candidate variants of CFTR gene were left and 53 variants were manually curated as pathogenic/likely-pathogenic (P/LP). These P/LP variants were adopted to estimate CF prevalence in three methods: the carrier frequency method, the permutation-combinations method and the Bayesian framework method. Allele frequencies of the 477 CFTR variants were compared with non-Finland European (NFE) and East Asian (EAS) from gnomAD database. To investigate the haplotype structure difference of CFTR, another 2067 whole-genome-sequencing samples from CCGT and 195 NFE from 1000 genome project were analyzed by Shapeit4 software.

Result

With the 53 manually curated P/LP variants in CFTR gene, we excluded individuals identified or suspected with CF and their parents in our cohorts and estimated the Chinese CF prevalence is approximately 1/128,434. Only 21 (39.6%) of the 53 variants were included in Caucasian specific CF screening panels, resulting in significantly under-estimation of CF prevalence in our children cohort (1/143,171 vs. 1/1,387,395, P = 5e−24) and parent’s cohort (1/110,127 vs. 1/872,437, P = 7e−10). The allele frequencies of six pathogenic variants (G970D, D979A, M469V, G622D, L88X, 1898+5G->T) were significantly higher in our cohorts compared with gnomAD-NFE population (all P-value < 0.1). Haplotype analysis showed more haplotype diversity in Chinese compared to Caucasians. In addition, G970D and F508del were founder mutation of Chinese and Caucasians with two SNPs (rs213950-rs1042077) identified as related genotype in exon region.

Conclusions

Chinese population showed significantly different genetic spectrum pattern in CFTR gene compared with Caucasian population, and thus a Chinese-specific CF screening panel is needed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-022-02279-9.

Variation in CFTR-dependent ‘β-sweating’ among healthy adults

The genetic disease cystic fibrosis (CF) results when mutations in the gene for the anion channel CFTR reduce CFTR’s activity below a critical level. CFTR activity = N·P_O·γ (number of channels x open probability x channel conductance). Small molecules are now available that partially restore CFTR function with dramatic improvements in health of CF subjects. Continued evaluation of these and other compounds in development will be aided by accurate assessments of CFTR function. However, measuring CFTR activity in vivo is challenging and estimates vary widely. The most accurate known measure of CFTR activity in vivo is the ‘β/M’ ratio of sweat rates, which is produced by stimulation with a β-adrenergic agonist cocktail referenced to the same individual’s methacholine-stimulated sweat rate. The most meaningful metric of CFTR activity is to express it as a percent of normal function, so it is critical to establish β/M carefully in a population of healthy control subjects. Here, we analyze β/M from a sample of 50 healthy adults in which sweat rates to cholinergic and β-adrenergic agonists were measured repeatedly (3 times) in multiple, (~50) identified sweat glands from each individual (giving ~20,000 measurements). The results show an approximately 7-fold range, 26–187% of the WT average set to 100%. These provide a benchmark against which other measures of CFTR activity can be compared. Factors contributing to β/M variation in healthy controls are discussed.

Mutation profiling of the c.1521_1523delCTT (p.Phe508del, F508del) CFTR allele using haplotype-resolved long-read next generation sequencing

Current approaches to characterize the mutational profile of the cystic fibrosis transmembrane conductance regulator (CFTR) gene are based on targeted mutation analysis (TMA) or whole gene studies derived from short-read next generation sequencing (NGS). However, these methods lack phasing capability which, in certain scenarios, can provide clinically valuable information. In the present work, we performed near-full length CFTR using Single-Molecule Real-Time Sequencing to produce haplotype-resolved data from both homozygous and heterozygous individuals for mutation c.1521_1523delCTT (p.Phe508del, F508del). This approach utilizes target enrichment of the CFTR gene using biotinylated probes, facilitates multiplexing samples in the same sequencing run, and utilizes fully-automated bioinformatics pipelines for error correction and variant calling. We show a remarkable conservation of F508del haplotype, consistent with the single gene founder effect, as well as diverse mutational profiles in non-F508del alleles. By the same method, 105 single nucleotide polymorphisms (SNPs) exhibiting invariant linkage to F508del CFTR (which better define the founder haplotype) were identified. High level homology between F508del sequences derived from heterozygotes, and those obtained from homozygous individuals, demonstrate accuracy of this method to produce haplotype resolved sequencing. The studies provide a new diagnostic technology for detailed analysis of complex CFTR alleles linked to disease severity. This article is protected by copyright. All rights reserved.

Cystic fibrosis

Cystic fibrosis is a monogenic disease considered to affect at least 100 000 people worldwide. Mutations in CFTR, the gene encoding the epithelial ion channel that normally transports chloride and bicarbonate, lead to impaired mucus hydration and clearance. Classical cystic fibrosis is thus characterised by chronic pulmonary infection and inflammation, pancreatic exocrine insufficiency, male infertility, and might include several comorbidities such as cystic fibrosis-related diabetes or cystic fibrosis liver disease. This autosomal recessive disease is diagnosed in many regions following newborn screening, whereas in other regions, diagnosis is based on a group of recognised multiorgan clinical manifestations, raised sweat chloride concentrations, or CFTR mutations. Disease that is less easily diagnosed, and in some cases affecting only one organ, can be seen in the context of gene variants leading to residual protein function. Management strategies, including augmenting mucociliary clearance and aggressively treating infections, have gradually improved life expectancy for people with cystic fibrosis. However, restoration of CFTR function via new small molecule modulator drugs is transforming the disease for many patients. Clinical trial pipelines are actively exploring many other approaches, which will be increasingly needed as survival improves and as the population of adults with cystic fibrosis increases. Here, we present the current understanding of CFTR mutations, protein function, and disease pathophysiology, consider strengths and limitations of current management strategies, and look to the future of multidisciplinary care for those with cystic fibrosis.

CFTR Cooperative Cis-Regulatory Elements in Intestinal Cells

About 8% of the human genome is covered with candidate cis-regulatory elements (cCREs). Disruptions of CREs, described as "cis-ruptions" have been identified as being involved in various genetic diseases. Thanks to the development of chromatin conformation study techniques, several long-range cystic fibrosis transmembrane conductance regulator (CFTR) regulatory elements were identified, but the regulatory mechanisms of the CFTR gene have yet to be fully elucidated. The aim of this work is to improve our knowledge of the CFTR gene regulation, and to identity factors that could impact the CFTR gene expression, and potentially account for the variability of the clinical presentation of cystic fibrosis as well as CFTR-related disorders. Here, we apply the robust GWAS3D score to determine which of the CFTR introns could be involved in gene regulation. This approach highlights four particular CFTR introns of interest. Using reporter gene constructs in intestinal cells, we show that two new introns display strong cooperative effects in intestinal cells. Chromatin immunoprecipitation analyses further demonstrate fixation of transcription factors network. These results provide new insights into our understanding of the CFTR gene regulation and allow us to suggest a 3D CFTR locus structure in intestinal cells. A better understand of regulation mechanisms of the CFTR gene could elucidate cases of patients where the phenotype is not yet explained by the genotype. This would thus help in better diagnosis and therefore better management. These cis-acting regions may be a therapeutic challenge that could lead to the development of specific molecules capable of modulating gene expression in the future.

Genomic, transcriptomic, and protein landscape profile of CFTR and cystic fibrosis

Cystic Fibrosis (CF) is caused most often by removal of amino acid 508 (Phe508del, deltaF508) within CFTR, yet dozens of additional CFTR variants are known to give rise to CF and many variants in the genome are known to contribute to CF pathology. To address CFTR coding variants, we developed a sequence-to-structure-to-dynamic matrix for all amino acids of CFTR using 233 vertebrate species, CFTR structure within a lipid membrane, and 20 ns of molecular dynamic simulation to assess known variants from the CFTR1, CFTR2, ClinVar, TOPmed, gnomAD, and COSMIC databases. Surprisingly, we identify 18 variants of uncertain significance within CFTR from diverse populations that are heritable and a likely cause of CF that have been understudied due to nonexistence in Caucasian populations. In addition, 15 sites within the genome are known to modulate CF pathology, where we have identified one genome region (chr11:34754985-34836401) that contributes to CF through modulation of expression of a noncoding RNA in epithelial cells. These 15 sites are just the beginning of understanding comodifiers of CF, where utilization of eQTLs suggests many additional genomics of CFTR expressing cells that can be influenced by genomic background of CFTR variants. This work highlights that many additional insights of CF genetics are needed, particularly as pharmaceutical interventions increase in the coming years.

Variant classifications, databases and genotype-phenotype correlations

Because CFTR gene studies now represent one of the most frequent genetic analyses routinely performed worldwide, the number of rare CFTR variants identified in various clinical situations, regularly increases. To provide appropriate diagnosis and prognosis to CF patients as well as appropriate genetic counseling to families, the clinical impact and the phenotypic spectrum of variants identified by diagnostic techniques need to be characterized. Three complementary locus specific databases, called CFTR1, CFTR2 and CFTR-France were developed to address these issues. Besides, the growing knowledge of the CF pathophysiology and the technical evolution in molecular biology allowed to identify candidate modifier genes, regulatory loci, epigenetic profiles and trans-regulators that could help to refine genotype-phenotype correlations at the individual level. These different factors may contribute to the large phenotypic variability between patients with CF, even when they carry identical CFTR variants, regarding lung function, meconium ileus susceptibility or the risk for developing CFTR-related diabetes and liver disease. Finally, the availability of new therapies that target the CFTR protein for numbers of CF patients led to the identification of 'good' and 'poor' responders, thus raising questions of pharmacogenetics factors that may influence treatment efficiency as a novel feature of the complexity of CF patients' management. © 2020 French Society of Pediatrics. Published by Elsevier Masson SAS. All rights reserved.

Current and future diagnosis of cystic fibrosis: Performance and limitations

Cystic fibrosis (CF) is the most frequent genetic disorder in the Caucasian population benefiting from systematic newborn screening tests. It is also the most frequent indication of prenatal and preimplantation genetic diagnosis for a single gene disorder. During the past thirty years, thanks in part to the evolution of diagnostic techniques, our knowledge on CFTR genetics and pathophysiological mechanisms involved in CF have significantly improved. With the implementation of newborn screening in France and in several countries, the diagnosis now often occurs in clinically asymptomatic infants and this has modified the criteria for CF diagnosis. Recently, guidelines for CF diagnosis have been reformulated in Europe and the US, in regard to sweat chloride usual values and disease terminology. This review describes the methods and molecular approaches that are used in routine practice or are being developed to detect CFTR protein dysfunction and to identify disease-causing CFTR variants. Ultimately, an optimal use of all these functional and genetic resources may improve patient care and therapeutic decision-making. © 2020 French Society of Pediatrics. Published by Elsevier Masson SAS. All rights reserved.

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.

Cystic fibrosis gene mutations and polymorphisms in Saudi men with infertility

BACKGROUND

Some mutations of the cystic fibrosis transmembrane regulator (CFTR) gene may impair spermatogenesis or cause a congenital absence of the vas deferens that manifests as isolated male infertility.

OBJECTIVE

Assess the frequency and analyze the spectrum of CFTR gene variations in Saudi men with primary infertility.

DESIGN

Prospective, cross-sectional.

SETTING

Tertiary care specialist hospital in Jeddah.

PATIENTS AND METHODS

Genomic DNA was extracted from peripheral blood samples of Saudi men who presented with primary infertility to the outpatient andrology clinic with either azoospermia or oligoasthenoteratozoospermia. Polymerase chain reaction and direct sequencing were used to identify all variants of the CFTR gene.

MAIN OUTCOME MEASURES

Proportion of the patients with a mutant CFTR gene and the spectrum of CFTR gene variations.

SAMPLE SIZE

50 infertile Saudi men.

RESULTS

This study identified 10 CFTR gene variants in 7 (14%) subjects (100 chromosomes). The detected variants and polymorphisms were: c.1408G>A, c.4389G>A, c.2562T>G, c.869+11C>T, c.2909-92G>A, c.3469-65C>A, c.1210-6delT, c.1210-6T>A, c.2988+1G>A, and c.1210-13GT>TG.

CONCLUSION

We demonstrated that 14% of the study subjects had one or more CFTR mutations and these were compounded in most of the affected patients. The spectrum of CFTR gene mutations in these subjects was similar to the mutations reported in other studies throughout the world.

LIMITATIONS

Small sample size and the lack of a control group.

CONFLICTS OF INTEREST

None.

Rare variant association testing in the non-coding genome

The development of next-generation sequencing technologies has opened-up some new possibilities to explore the contribution of genetic variants to human diseases and in particular that of rare variants. Statistical methods have been developed to test for association with rare variants that require the definition of testing units and, in these testing units, the selection of qualifying variants to include in the test. In the coding regions of the genome, testing units are usually the different genes and qualifying variants are selected based on their functional effects on the encoded proteins. Extending these tests to the non-coding regions of the genome is challenging. Testing units are difficult to define as the non-coding genome organisation is still rather unknown. Qualifying variants are difficult to select as the functional impact of non-coding variants on gene expression is hard to predict. These difficulties could explain why very few investigators so far have analysed the non-coding parts of their whole genome sequencing data. These non-coding parts yet represent the vast majority of the genome and some studies suggest that they could play a major role in disease susceptibility. In this review, we discuss recent experimental and statistical developments to gain knowledge on the non-coding genome and how this knowledge could be used to include rare non-coding variants in association tests. We describe the few studies that have considered variants from the non-coding genome in association tests and how they managed to define testing units and select qualifying variants.

Molecular Diagnosis and Genetic Counseling of Cystic Fibrosis and Related Disorders: New Challenges

Identification of the cystic fibrosis transmembrane conductance regulator (CFTR) gene and its numerous variants opened the way to fantastic breakthroughs in diagnosis, research and treatment of cystic fibrosis (CF). The current and future challenges of molecular diagnosis of CF and CFTR-related disorders and of genetic counseling are here reviewed. Technological advances have enabled to make a diagnosis of CF with a sensitivity of 99% by using next generation sequencing in a single step. The detection of heretofore unidentified variants and ethnic-specific variants remains challenging, especially for newborn screening (NBS), CF carrier testing and genotype-guided therapy. Among the criteria for assessing the impact of variants, population genetics data are insufficiently taken into account and the penetrance of CF associated with CFTR variants remains poorly known. The huge diversity of diagnostic and genetic counseling indications for CFTR studies makes assessment of variant disease-liability critical. This is especially discussed in the perspective of wide genome analyses for NBS and CF carrier screening in the general population, as future challenges.

The influence of CFTR complex alleles on precision therapy of cystic fibrosis

CFTR is an extensively studied gene and multiple sequence variants have been identified, many of which still need to be defined as neutral or disease causing. Complex alleles are defined when at least two variants are identified on the same allele. Each pathogenic variant can affect distinct steps of the CFTR biogenesis. As CFTR modulators are being developed to alleviate specific defects, pathogenic variants need to be characterized to propose adequate treatments. Conversely, cis-variants can affect treatment response when defects are additive or if they alter the binding or efficacy of the modulator. Hence, complex alleles increase the complexity of CFTR variant classification and need to be assigned as neutral, disease causing or modulating treatment efficacy. This review was based on a symposium session presented at the 16th ECFS Basic Science Conference, Dubrovnik, Croatia, 27 to 30 March, 2019.

Forskolin-induced swelling of intestinal organoids correlates with disease severity in adults with cystic fibrosis and homozygous F508del mutations

BACKGROUND

CFTR function measurements in intestinal organoids may help to better characterise individual disease expression in F508del homozygous people. Our objective was to study correlations between CFTR function as measured with forskolin-induced swelling in rectal organoids with clinical parameters in adult patients with homozygous F508del mutations.

METHODS

Multicentre observational study. Thirty-four adults underwent rectal biopsy, pulmonary function tests (FEV₁ and FVC), chest X-ray and chest CT. Body-mass index (BMI) was assessed at study visit and exacerbation rate was determined during five years prior to study visit. Organoids were cultured and measured after stimulation with 5 µm forskolin for three hours to quantitate CFTR residual function.

FINDINGS

FIS was positively correlated with FEV₁ (r = 0.36, 95% CI 0.02-0.62, p = 0.04) and BMI (r = 0.42, 95% CI 0.09-0.66, p = 0.015). FIS was negatively correlated with PRAGMA-CF CT score for% of disease (r = -0.37, 95% CI -0.62- -0.03, p = 0.049). We found no significant correlation between FIS and chest radiography score for CF (r = -0.16, 95% CI -0.48-0.20, p = 0.44). We observed a trend between higher FIS and a lower mean number of exacerbations over the last 5 years of observation, but this was not statistically significant (Poisson regression, p = 0.089).

INTERPRETATION

FIS of intestinal organoids varied between subjects with homozygous F508del and correlated with pulmonary and nutritional parameters. These findings suggest that differences at low CFTR residual function may contribute to clinical heterogeneity in F508del homozygous patients and small changes in CFTR residual function might impact long-term disease expression.

Predictive factors for lumacaftor/ivacaftor clinical response

BACKGROUND

Ivacaftor-lumacaftor combination therapy corrects the F508 del-CFTR mutated protein which causes Cystic Fibrosis. The clinical response of the patients treated with the combination therapy is highly variable. This study aimed to determine factors involved in the individual's response to lumacaftor-ivacaftor therapy.

METHODS

Sweat test was assessed at baseline and after 6 months of ivacaftor-lumacaftor treatment in 41 homozygous F508del children and young adults. β-adrenergic peak sweat secretion, nasal potential difference (NPD) and intestinal current measurements (ICM) were performed in patients accepting these tests. Seric level of lumacaftor and ivacaftor were determined and additional CFTR variant were searched.

RESULTS

Sweat chloride concentration significantly decreased after treatment, whereas the β-adrenergic peak sweat response did not vary in 9 patients who underwent these tests. The average level of F508del-CFTR activity rescue reached up to 15% of the normal level in intestinal epithelium, as studied by ICM in 12 patients (p = .03) and 20% of normal in the nasal epithelium in NPD tests performed in 21 patients (NS). There was no significant correlation between these changes and improvements in FEV₁ at 6 months. Serum drug levels did not correlate with changes in FEV₁, BMI-Zscore or other CFTR activity biomarkers. Additional exonic variants were identified in 4 patients. The F87L-I1027T-F508del-CFTR complex allele abolished the lumacaftor corrector effect.

CONCLUSION

This observational study investigates a number of potential factors linked to the clinical response of F508del homozygous patients treated with lumacaftor-ivacaftor combination therapy. Lumacaftor and ivacaftor blood levels are not associated with the clinical response. Additional exonic variants may influence protein correction.

Functional characterization and phenotypic spectrum of three recurrent disease-causing deep intronic variants of the CFTR gene

BACKGROUND

The CFTR genotype remains incomplete in 1% of Cystic Fibrosis (CF) cases, because only one or no disease-causing variants is detected after extended analysis. This fraction is probably higher in CFTR-Related Disorders (CFTR-RD). Deep-intronic CFTR variants are putative candidates to fill this gap. However, the recurrence, phenotypic spectrum and full molecular characterization of newly reported variants are unknown.

METHODS

Minigenes and analysis of CFTR transcripts in nasal epithelial cells were used to determine the impact on CFTR splicing of intronic variants that we previously identified by next generation sequencing of the whole CFTR locus. Phenotypic data were collected in 19 patients with CF and CFTR-RD, in whom one of the deep intronic variants has been detected.

RESULTS

Three deep-intronic variants promoted the inclusion of pseudo-exons (PE) in the CFTR transcript, hindering the synthesis of a functional protein. The c.2989-313A > T variant, detected in four patients with CF or CFTR-RD from three different families, led to the inclusion of a 118 bp PE. The c.3469-1304C > G variant promoted the inclusion of a 214 bp-PE and was identified in five patients with CF from four families. Haplotype analysis confirmed that this variant was associated with one CF chromosome of African origin. The most represented variant in our cohort was the c.3874-4522A > G, detected in 10 patients with various phenotypes, from male infertility to CF with pancreatic insufficiency.

CONCLUSION

These three deep intronic CFTR variants are associated with a large phenotypic spectrum, including typical CF. They should be included in CF diagnostic testing and carrier screening 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.

Novel long-range regulatory mechanisms controlling PKD2 gene expression

Background

Cis-regulatory elements control gene expression over large distances through the formation of chromatin loops, which allow contact between enhancers and gene promoters. Alterations in cis-acting regulatory systems could be linked to human genetic diseases. Here, we analyse the spatial organization of a large region spanning the polycystic kidney disease 2 (PKD2) gene, one of the genes responsible of autosomal dominant polycystic kidney disease (ADPKD).

Results

By using chromosome conformation capture carbon copy (5C) technology in primary human renal cyst epithelial cells, we identify novel contacts of the PKD2 promoter with chromatin regions, which display characteristics of regulatory elements. In parallel, by using functional analysis with a reporter assay, we demonstrate that three DNAse I hypersensitive sites regions are involved in the regulation of PKD2 gene expression.

Conclusions

Finally, through alignment of CCCTC-binding factor (CTCF) sites, we suggest that these novel enhancer elements are brought to the PKD2 promoter by chromatin looping via the recruitment of CTCF.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4892-6) contains supplementary material, which is available to authorized users.

Current and future molecular approaches in the diagnosis of cystic fibrosis

INTRODUCTION

Cystic Fibrosis is among the first diseases to have general population genetic screening tests and one of the most common indications of prenatal and preimplantation genetic diagnosis for single gene disorders. During the past twenty years, thanks to the evolution of diagnostic techniques, our knowledge of CFTR genetics and pathophysiological mechanisms involved in cystic fibrosis has significantly improved. Areas covered: Sanger sequencing and quantitative methods greatly contributed to the identification of more than 2,000 sequence variations reported worldwide in the CFTR gene. We are now entering a new technological age with the generalization of high throughput approaches such as Next Generation Sequencing and Droplet Digital PCR technologies in diagnostics laboratories. These powerful technologies open up new perspectives for scanning the entire CFTR locus, exploring modifier factors that possibly influence the clinical evolution of patients, and for preimplantation and prenatal diagnosis. Expert commentary: Such breakthroughs would, however, require powerful bioinformatics tools and relevant functional tests of variants for analysis and interpretation of the resulting data. Ultimately, an optimal use of all those resources may improve patient care and therapeutic decision-making.

Cis variants identified in F508del complex alleles modulate CFTR channel rescue by small molecules

Molecules correcting the trafficking (correctors) and gating defects (potentiators) of the cystic fibrosis causing mutation c.1521_1523delCTT (p.Phe508del) begin to be a useful treatment for CF patients bearing p.Phe508del. This mutation has been identified in different genetic contexts, alone or in combination with variants in cis. Until now, 21 exonic variants in cis of p.Phe508del have been identified, albeit at a low frequency. The aim of this study was to evaluate their impact on the efficacy of CFTR-directed corrector/potentiator therapy (Orkambi). The analysis by minigene showed that two out of 15 cis variants tested increased exon skipping (c.609C > T and c.2770G > A). Four cis variants were studied functionally in the absence of p.Phe508del, one of which was found to be deleterious for protein maturation c.1399C > T (p.Leu467Phe). In the presence of p.Phe508del, this variant was the only to prevent the response to Orkambi treatment. This study showed that some patients carrying p.Phe508del complex alleles are predicted to poorly respond to corrector/potentiator treatments. Our results underline the importance to validate treatment efficacy in the context of complex alleles.

CFTR-France, a national relational patient database for sharing genetic and phenotypic data associated with rare CFTR variants

Most of the 2,000 variants identified in the CFTR (cystic fibrosis transmembrane regulator) gene are rare or private. Their interpretation is hampered by the lack of available data and resources, making patient care and genetic counseling challenging. We developed a patient-based database dedicated to the annotations of rare CFTR variants in the context of their cis- and trans-allelic combinations. Based on almost 30 years of experience of CFTR testing, CFTR-France (https://cftr.iurc.montp.inserm.fr/cftr) currently compiles 16,819 variant records from 4,615 individuals with cystic fibrosis (CF) or CFTR-RD (related disorders), fetuses with ultrasound bowel anomalies, newborns awaiting clinical diagnosis, and asymptomatic compound heterozygotes. For each of the 736 different variants reported in the database, patient characteristics and genetic information (other variations in cis or in trans) have been thoroughly checked by a dedicated curator. Combining updated clinical, epidemiological, in silico, or in vitro functional data helps to the interpretation of unclassified and the reassessment of misclassified variants. This comprehensive CFTR database is now an invaluable tool for diagnostic laboratories gathering information on rare variants, especially in the context of genetic counseling, prenatal and preimplantation genetic diagnosis. CFTR-France is thus highly complementary to the international database CFTR2 focused so far on the most common CF-causing alleles.

Alteration of protein function by a silent polymorphism linked to tRNA abundance

Synonymous single nucleotide polymorphisms (sSNPs) are considered neutral for protein function, as by definition they exchange only codons, not amino acids. We identified an sSNP that modifies the local translation speed of the cystic fibrosis transmembrane conductance regulator (CFTR), leading to detrimental changes to protein stability and function. This sSNP introduces a codon pairing to a low-abundance tRNA that is particularly rare in human bronchial epithelia, but not in other human tissues, suggesting tissue-specific effects of this sSNP. Up-regulation of the tRNA cognate to the mutated codon counteracts the effects of the sSNP and rescues protein conformation and function. Our results highlight the wide-ranging impact of sSNPs, which invert the programmed local speed of mRNA translation and provide direct evidence for the central role of cellular tRNA levels in mediating the actions of sSNPs in a tissue-specific manner.

Synonymous single nucleotide polymorphisms (sSNPs) occur at high frequency in the human genome and are associated with ~50 diseases in humans; the responsible molecular mechanisms remain enigmatic. Here, we investigate the impact of the common sSNP, T2562G, on cystic fibrosis transmembrane conductance regulator (CFTR). Although this sSNP, by itself, does not cause cystic fibrosis (CF), it is prevalent in patients with CFTR-related disorders. T2562G sSNP modifies the local translation speed at the Thr854 codon, leading to changes in CFTR stability and channel function. This sSNP introduces a codon pairing to a low-abundance tRNA, which is particularly rare in human bronchial epithelia, but not in other human tissues, suggesting a tissue-specific effect of this sSNP. Enhancement of the cellular concentration of the tRNA cognate to the mutant ACG codon rescues the stability and conduction defects of T2562G-CFTR. These findings reveal an unanticipated mechanism—inverting the programmed local speed of mRNA translation in a tRNA-dependent manner—for sSNP-associated diseases.