The clinical utility of sequencing the entirety of CFTR.

BACKGROUND: Cystic fibrosis (CF) is caused by deleterious variants in each CFTR gene. We investigated the utility of whole-gene CFTR sequencing when fewer than two pathogenic or likely pathogenic (P/LP) variants were detected by conventional testing (sequencing of exons and flanking introns) of CFTR. METHODS: Individuals with features of CF and a CF-diagnostic sweat chloride concentration with zero or one P/LP variants identified by conventional testing enrolled in the CF Mutation Analysis Program (MAP) underwent whole-gene CFTR sequencing. Replication was performed on individuals enrolled in the CF Genome Project (CFGP), followed by phenotype review and interrogation of other genes. RESULTS: Whole-gene sequencing identified a second P/LP variant in 20/43 MAP enrollees (47 %) and 10/22 CFGP enrollees (45 %) who had one P/LP variant after conventional testing. No P/LP variants were detected when conventional testing was negative (MAP: n = 43; CFGP: n = 13). Genome-wide analysis was unable to find an alternative etiology in CFGP participants with fewer than two P/LP CFTR variants and CF could not be confirmed in 91 % following phenotype re-review. CONCLUSIONS: Whole-gene CFTR analysis is beneficial in individuals with one previously-identified P/LP variant and a CF-diagnostic sweat chloride. Negative conventional CFTR testing indicates that the phenotype should be re-evaluated.

Recessively Inherited Deficiency of Secreted WFDC2 (HE4) Causes Nasal Polyposis and Bronchiectasis.

RATIONALE: Bronchiectasis is a pathological dilatation of the bronchi in the respiratory airways associated with environmental or genetic causes (e.g., cystic fibrosis, primary ciliary dyskinesia and primary immunodeficiency disorders), but most cases remain idiopathic. OBJECTIVES: To identify novel genetic defects in unsolved cases of bronchiectasis presenting with severe rhinosinusitis, nasal polyposis, and pulmonary Pseudomonas aeruginosa infection. METHODS: DNA was analyzed by next-generation or targeted Sanger sequencing. RNA was analyzed by quantitative PCR and single-cell RNA sequencing. Patient-derived, cells, cell cultures and secretions (mucus, saliva, seminal fluid) were analyzed by Western blotting and immunofluorescence microscopy, and mucociliary activity was measured. Blood serum was analyzed by electrochemiluminescence immunoassay. Protein structure and proteomic analyses were used to assess the impact of a disease-causing founder variant. MEASUREMENTS AND MAIN RESULTS: We identified bi-allelic pathogenic variants in WFDC2 in 11 individuals from 10 unrelated families originating from the United States, Europe, Asia, and Africa. Expression of WFDC2 was detected predominantly in secretory cells of control airway epithelium and also in submucosal glands. We demonstrate that WFDC2 is below the limit of detection in blood serum and hardly detectable in samples of saliva, seminal fluid, and airway surface liquid from WFDC2-deficient individuals. Computer simulations and deglycosylation assays indicate that the disease-causing founder variant p.Cys49Arg structurally hampers glycosylation and thus secretion of mature WFDC2. CONCLUSIONS: WFDC2 dysfunction defines a novel molecular etiology of bronchiectasis characterized by the deficiency of a secreted component of the airways. A commercially available blood test combined with genetic testing allows its diagnosis. This article is open access and distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/).

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Genetic variation in severe cystic fibrosis liver disease is associated with novel mechanisms for disease pathogenesis.

BACKGROUND AND AIMS: It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. APPROACH AND RESULTS: Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. CONCLUSION: These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD.

Safety and efficacy of the epithelial sodium channel blocker idrevloride in people with primary ciliary dyskinesia (CLEAN-PCD): a multinational, phase 2, randomised, double-blind, placebo-controlled crossover trial.

BACKGROUND: Mucociliary clearance is dysfunctional in people with primary ciliary dyskinesia, resulting in the accumulation of dehydrated mucus in the airways that is difficult to clear. We undertook a study to assess the benefit on lung function of treatment with a nebulised epithelial sodium channel (ENaC) blocker, idrevloride, with or without hypertonic saline, in people with primary ciliary dyskinesia. METHODS: The CLEAN-PCD trial was a phase 2, randomised, double-blind, placebo-controlled crossover trial conducted at 32 tertiary adult and paediatric care centres and university hospitals in Canada, Denmark, Germany, Italy, the Netherlands, Poland, the UK, and the USA. People with a confirmed diagnosis of primary ciliary dyskinesia, aged 12 years or older, with a percentage of predicted FEV1 (ppFEV1) in the range of 40% to <90%, were randomly assigned in a 2:2:1:1 ratio (block size=6), stratified by ppFEV1 at screening, to one of four sequences: (1) idrevloride in hypertonic saline in treatment period 1 then hypertonic saline in treatment period 2; (2) hypertonic saline in treatment period 1 then idrevloride in hypertonic saline in treatment period 2; (3) idrevloride in treatment period 1 then placebo in treatment period 2; and (4) placebo in treatment period 1 then idrevloride in treatment period 2. The idrevloride dose was 85 µg and hypertonic saline was 4·2% NaCl. 3 mL of each study treatment was nebulised twice daily for 28 days in treatment periods 1 and 2; the two 28-day treatment periods were separated by a 28-day washout period. The primary endpoint was absolute change from baseline in ppFEV1 after 28 days. Safety assessments and reports of adverse events were made at clinic visits during each treatment period and by a follow-up telephone call 28 days after the last dose of study drug. Additionally, adverse events could be reported at a follow-up telephone call 3 days after the start of dosing and as they arose. Participants who received at least one dose of study drug were included in the safety analyses (safety set), and those who also had spirometry data were included in the efficacy analyses (full analysis set). The completed study is registered (EudraCT 2015-004917-26; ClinicalTrials.govNCT02871778). FINDINGS: Between Sep 14, 2016, and May 31, 2018, 216 patients were screened and 123 were randomly assigned to one of four crossover sequences. Across the two treatment periods, treatment with idrevloride in hypertonic saline was initiated in 80 patients and completed in 78 patients (all 78 had data available and were included in the analysis); hypertonic saline initiated in 81 patients and completed in 76 patients (75 had data available and were included in the analysis); idrevloride initiated in 37 patients and completed in 35 patients (34 had data available and were included in the analysis); and placebo initiated in 36 patients and completed in 34 patients (all 34 had data available and were included in the analysis). Greater absolute increases in ppFEV1 from baseline to 28 days of treatment were seen with idrevloride in hypertonic saline (least-squares mean absolute change from baseline 1·0 percentage points, 95% CI -0·4 to 2·4) than with hypertonic saline alone (least-squares mean absolute change from baseline of -0·5 percentage points, -2·0 to 0·9; difference 1·5 percentage points, 95% CI <0·1 to 3·0; p=0·044). There was no significant difference in ppFEV1 for the parallel comparison of idrevloride in hypertonic saline compared with placebo or the crossover comparison of idrevloride with placebo. Adverse events were similar across treatments (57 to 65% of patients). Cough occurred in a greater proportion of participants during treatments that contained idrevloride or hypertonic saline compared with placebo, and oropharyngeal pain occurred in a greater proportion of participants during idrevloride treatments than during treatment with hypertonic saline alone or placebo, whereas chest discomfort was more common during treatments that included hypertonic saline. INTERPRETATION: In this phase 2 crossover study, idrevloride in hypertonic saline was safe and associated with improved lung function over a 28-day period in people with primary ciliary dyskinesia compared with hypertonic saline alone. Larger, longer clinical studies are warranted to explore the potential benefits of idrevloride in combination with hypertonic saline in people with primary ciliary dyskinesia. FUNDING: Parion Sciences, under agreement with Vertex Pharmaceuticals.

Situs Ambiguus Is Associated With Adverse Clinical Outcomes in Children With Primary Ciliary Dyskinesia.

BACKGROUND: Primary ciliary dyskinesia (PCD) is a rare disorder of motile cilia associated with situs abnormalities. At least 12% of patients with PCD have situs ambiguus (SA), including organ laterality defects falling outside normal arrangement (situs solitus [SS]) or mirror image inversion (situs inversus totalis [SIT]). RESEARCH QUESTION: Do patients with PCD and SA achieve worse clinical outcomes compared with those with SS or SIT? STUDY DESIGN AND METHODS: This cross-sectional, multicenter study evaluated participants aged 21 years or younger with PCD. Participants were classified as having SA, including heterotaxy, or not having SA (SS or SIT). Markers of disease severity were compared between situs groups, adjusting for age at enrollment and severe CCDC39 or CCDC40 genotype, using generalized linear models and logistic and Poisson regression. RESULTS: In 397 participants with PCD (mean age, 8.4 years; range, 0.1-21), 42 patients were classified as having SA, including 16 patients (38%) with complex cardiovascular malformations or atrial isomerism, 13 patients (31%) with simple CVM, and 13 patients (31%) without cardiovascular malformations. Of these, 15 patients (36%) underwent cardiac surgery, 24 patients (57%) showed an anatomic spleen abnormality, and seven patients (17%) showed both. The remaining 355 participants did not have SA, including 152 with SIT and 203 with SS. Overall, 70 participants (17%) harbored the severe CCDC39 or CCDC40 genotype. Compared with participants without SA, those with SA showed lower median BMI z scores (P = .03), lower FVC z scores (P = .01), and more hospitalizations and IV antibiotic courses for acute respiratory infections during the 5 years before enrollment (P < .01). Participants with cardiovascular malformations requiring surgery or with anatomic spleen abnormalities showed lower median BMI z scores and more hospitalizations and IV therapies for respiratory illnesses compared with participants without SA. INTERPRETATION: Children with PCD and SA achieve worse nutritional and pulmonary outcomes with more hospitalizations for acute respiratory illnesses than those with SS or SIT combined. Poor nutrition and increased hospitalizations for respiratory infections in participants with SA and PCD are associated with cardiovascular malformations requiring cardiac surgery, splenic anomalies, or both. TRIAL REGISTRY: ClinicalTrials.gov; Nos.: NCT02389049 and NCT00323167; URL: www. CLINICALTRIALS: gov.

First reports of primary ciliary dyskinesia caused by a shared DNAH11 allele in Canadian Inuit.

BACKGROUND: Primary ciliary dyskinesia (PCD) is typically an autosomal recessive disease characterized by recurrent infections of the lower respiratory tract, frequent and severe otitis media, chronic rhinosinusitis, neonatal respiratory distress, and organ laterality defects. While severe lower respiratory tract infections and bronchiectasis are common in Inuit, PCD has not been recognized in this population. METHODS: We report a case series of seven Inuit patients with PCD identified by genetic testing in three Canadian PCD centers. RESULTS: Patients ranged from 4 to 59 years of age (at time of last evaluation) and originated in the Qikiqtaaluk region (Baffin Island, n = 5), Nunavut, or Nunavik (northern Quebec, n = 2), Canada. They had typical features of PCD, including neonatal respiratory distress (five patients), situs inversus totalis (four patients), bronchiectasis (four patients), chronic atelectasis (six patients), and chronic otitis media (six patients). Most had chronic rhinitis. Genetic evaluation demonstrated that all had homozygous pathogenic variants in DNAH11 at NM_001277115.1:c.4095+2C>A. CONCLUSIONS: The discovery of this homozygous DNAH11 variant in widely disparate parts of the Nunangat (Inuit homelands) suggests this is a founder mutation that may be widespread in Inuit. Thus, PCD may be an important cause of chronic lung, sinus, and middle ear disease in this population. Inuit with chronic lung disease, including bronchiectasis or laterality defects, should undergo genetic testing for PCD. Consideration of including PCD genetic analysis in routine newborn screening should be considered in Inuit regions.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients.

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

Airway Disease in Children with Primary Ciliary Dyskinesia: Impact of Ciliary Ultrastructure Defect and Genotype.

Rationale: Primary ciliary dyskinesia (PCD) is characterized by impaired mucociliary clearance, recurrent respiratory infections, progressive airway damage, and obstructive lung disease. Although the association of ciliary ultrastructure defect/genotype with the severity of airflow obstruction has been well characterized, their association with airway abnormalities on chest computed tomography (CT) has been minimally evaluated. Objectives: We sought to delineate the association of ciliary defect class/genotype with chest CT scores in children with PCD. Methods: Cross-sectional analysis of children with PCD (N = 146) enrolled in a prospective multicenter observational study, stratified by defect type: outer dynein arm (ODA), ODA/inner dynein arm (IDA), IDA/microtubular disorganization (MTD), and normal/near normal ultrastructure with associated genotypes. CTs were scored using the MERAGMA-PCD (Melbourne-Rotterdam Annotated Grid Morphometric Analysis for PCD), evaluating airway abnormalities in a hierarchical order: atelectasis, bronchiectasis, bronchial wall thickening, and mucus plugging/tree-in-bud opacities. The volume fraction of each component was expressed as the percentage of total lung volume. The percentage of disease was computed as the sum of all components. Regression analyses were used to describe the association between clinical predictors and CT scores. Results: Acceptable chest CTs were obtained in 141 children (71 male): 57 ODA, 20 ODA/IDA, 40 IDA/MTD, and 24 normal/near normal. The mean (standard deviation) age was 8.5 (4.6) years, forced expiratory volume in 1 second (FEV1) percent predicted was 82.4 (19.5), and %Disease was 4.6 (3.5). Children with IDA/MTD defects had a higher %Disease compared with children with ODA defects (2.71% higher [95% confidence interval (CI), 1.37-4.06; P < 0.001]), driven by higher %Mucus plugging (2.35% higher [1.43-3.26; P < 0.001]). Increasing age, lower body mass index, and lower FEV1 were associated with a higher %Disease (0.23%; 95% CI, 0.11-0.35; P < 0.001 and 0.03%; 95% CI, 0.01-0.04; P = 0.008 and 0.05%; 95% CI, 0.01-0.08; P = 0.011, respectively). Conclusions: Children with IDA/MTD defects had significantly greater airway disease on CT, primarily mucus plugging, compared with children with ODA defects.

Laterality Defects in Primary Ciliary Dyskinesia: Relationship to Ultrastructural Defect or Genotype.

Rationale: The association between organ laterality abnormalities and ciliary ultrastructural defect or genotype in primary ciliary dyskinesia is poorly understood. Objectives: To determine if there is an association between presence and/or type of laterality abnormality and ciliary ultrastructural defect or genotype. Methods: Participants with primary ciliary dyskinesia in a multicenter, prospective study were grouped based on ciliary ultrastructural defect or genotype. In a retrospective analysis of these data, the association of ciliary ultrastructural defect or genotype and likelihood of a laterality abnormality was evaluated by logistic regression adjusted for presence of two loss-of-function versus one or more not-loss-of-function variants. Results: Of 559 participants, 286 (51.2%), 215 (38.5%), and 58 (10.4%) were identified as having situs solitus, situs inversustotalis, and situs ambiguus, respectively; heterotaxy, defined as situs ambiguus with complex cardiovascular defects, was present in 14 (2.5%). Compared with the group with inner dynein arm defects with microtubular disorganization, laterality defects were more likely in the outer dynein arm defects group (odds ratio [OR], 2.07; 95% confidence interval [CI], 1.21-3.54; P < 0.01) and less likely in the normal/near normal ultrastructure group (OR, 0.04; 95% CI, 0.013-0.151; P < 0.01). Heterotaxy was present in 11 of 242 (4.5%) in the outer dynein arm defects group but 0 of 96 in the inner dynein arm defects with microtubular disorganization group (P = 0.038). Conclusion: In primary ciliary dyskinesia, risk of a laterality abnormality differs by ciliary ultrastructural defect. Pathophysiologic mechanisms underlying these differences require further exploration.

Airway Inflammation in Children with Primary Ciliary Dyskinesia.

Rationale: The role of airway inflammation in disease pathogenesis in children with primary ciliary dyskinesia (PCD) is poorly understood. Objectives: We investigated relationships between sputum inflammation measurements, age, lung function, bronchiectasis, airway infection, and ultrastructural defects in children with PCD. Methods: Spontaneously expectorated sputum was collected from clinically stable children and adolescents with PCD ages 6 years and older participating in a multicenter, observational study. Sputum protease and inflammatory cytokine concentrations were correlated with age, lung function, and chest computed tomography measures of structural lung disease, whereas differences in concentrations were compared between ultrastructural defect categories and between those with and without detectable bacterial infection. Results: Sputum from 77 children with PCD (39 females [51%]; mean [standard deviation] age, 13.9 [4.9] yr; mean [standard deviation] forced expiratory volume in 1 second [FEV1]% predicted, 80.8 [20.5]) was analyzed. Sputum inflammatory marker measurements, including neutrophil elastase activity, IL-1ß (interleukin-1ß), IL-8, and TNF-α (tumor necrosis factor α) concentrations, correlated positively with age, percentage of bronchiectasis, and percentage of total structural lung disease on computed tomography, and negatively with lung function. Correlations between neutrophil elastase concentrations and FEV1% predicted and percentage of bronchiectasis were -0.32 (95% confidence interval, -0.51 to -0.10) and 0.46 (0.14 to 0.69), respectively. Sputum neutrophil elastase, IL-1ß, and TNF-α concentrations were higher in those with detectable bacterial pathogens. Participants with absent inner dynein arm and microtubular disorganization had similar inflammatory profiles compared with participants with outer dynein arm defects. Conclusions: In this multicenter pediatric PCD cohort, elevated concentrations of sputum proteases and cytokines were associated with impaired lung function and structural damage as determined by chest computed tomography, suggesting that sputum inflammatory measurements could serve as biomarkers in PCD.

MagicalRsq: Machine-learning-based genotype imputation quality calibration.

Whole-genome sequencing (WGS) is the gold standard for fully characterizing genetic variation but is still prohibitively expensive for large samples. To reduce costs, many studies sequence only a subset of individuals or genomic regions, and genotype imputation is used to infer genotypes for the remaining individuals or regions without sequencing data. However, not all variants can be well imputed, and the current state-of-the-art imputation quality metric, denoted as standard Rsq, is poorly calibrated for lower-frequency variants. Here, we propose MagicalRsq, a machine-learning-based method that integrates variant-level imputation and population genetics statistics, to provide a better calibrated imputation quality metric. Leveraging WGS data from the Cystic Fibrosis Genome Project (CFGP), and whole-exome sequence data from UK BioBank (UKB), we performed comprehensive experiments to evaluate the performance of MagicalRsq compared to standard Rsq for partially sequenced studies. We found that MagicalRsq aligns better with true R2 than standard Rsq in almost every situation evaluated, for both European and African ancestry samples. For example, when applying models trained from 1,992 CFGP sequenced samples to an independent 3,103 samples with no sequencing but TOPMed imputation from array genotypes, MagicalRsq, compared to standard Rsq, achieved net gains of 1.4 million rare, 117k low-frequency, and 18k common variants, where net gains were gained numbers of correctly distinguished variants by MagicalRsq over standard Rsq. MagicalRsq can serve as an improved post-imputation quality metric and will benefit downstream analysis by better distinguishing well-imputed variants from those poorly imputed. MagicalRsq is freely available on GitHub.

Pleiotropic modifiers of age-related diabetes and neonatal intestinal obstruction in cystic fibrosis.

Individuals with cystic fibrosis (CF) develop complications of the gastrointestinal tract influenced by genetic variants outside of CFTR. Cystic fibrosis-related diabetes (CFRD) is a distinct form of diabetes with a variable age of onset that occurs frequently in individuals with CF, while meconium ileus (MI) is a severe neonatal intestinal obstruction affecting ∼20% of newborns with CF. CFRD and MI are slightly correlated traits with previous evidence of overlap in their genetic architectures. To better understand the genetic commonality between CFRD and MI, we used whole-genome-sequencing data from the CF Genome Project to perform genome-wide association. These analyses revealed variants at 11 loci (6 not previously identified) that associated with MI and at 12 loci (5 not previously identified) that associated with CFRD. Of these, variants at SLC26A9, CEBPB, and PRSS1 associated with both traits; variants at SLC26A9 and CEBPB increased risk for both traits, while variants at PRSS1, the higher-risk alleles for CFRD, conferred lower risk for MI. Furthermore, common and rare variants within the SLC26A9 locus associated with MI only or CFRD only. As expected, different loci modify risk of CFRD and MI; however, a subset exhibit pleiotropic effects indicating etiologic and mechanistic overlap between these two otherwise distinct complications of CF.

Small-molecule eRF3a degraders rescue CFTR nonsense mutations by promoting premature termination codon readthrough.

The vast majority of people with cystic fibrosis (CF) are now eligible for CF transmembrane regulator (CFTR) modulator therapy. The remaining individuals with CF harbor premature termination codons (PTCs) or rare CFTR variants with limited treatment options. Although the clinical modulator response can be reliably predicted using primary airway epithelial cells, primary cells carrying rare CFTR variants are scarce. To overcome this obstacle, cell lines can be created by overexpression of mouse Bmi-1 and human TERT (hTERT). Using this approach, we developed 2 non-CF and 6 CF airway epithelial cell lines, 3 of which were homozygous for the W1282X PTC variant. The Bmi-1/hTERT cell lines recapitulated primary cell morphology and ion transport function. The 2 F508del-CFTR cell lines responded robustly to CFTR modulators, which was mirrored in the parent primary cells and in the cell donors' clinical response. Cereblon E3 ligase modulators targeting eukaryotic release factor 3a (eRF3a) rescued W1282X-CFTR function to approximately 20% of WT levels and, when paired with G418, rescued G542X-CFTR function to approximately 50% of WT levels. Intriguingly, eRF3a degraders also diminished epithelial sodium channel (ENaC) function. These studies demonstrate that Bmi-1/hTERT cell lines faithfully mirrored primary cell responses to CFTR modulators and illustrate a therapeutic approach to rescue CFTR nonsense mutations.

Accounting for population structure in genetic studies of cystic fibrosis.

CFTR F508del (c.1521_1523delCTT, p.Phe508delPhe) is the most common pathogenic allele underlying cystic fibrosis (CF), and its frequency varies in a geographic cline across Europe. We hypothesized that genetic variation associated with this cline is overrepresented in a large cohort (N > 5,000) of persons with CF who underwent whole-genome sequencing and that this pattern could result in spurious associations between variants correlated with both the F508del genotype and CF-related outcomes. Using principal-component (PC) analyses, we showed that variation in the CFTR region disproportionately contributes to a PC explaining a relatively high proportion of genetic variance. Variation near CFTR was correlated with population structure among persons with CF, and this correlation was driven by a subset of the sample inferred to have European ancestry. We performed genome-wide association studies comparing persons with CF with one versus two copies of the F508del allele; this allowed us to identify genetic variation associated with the F508del allele and to determine that standard PC-adjustment strategies eliminated the significant association signals. Our results suggest that PC adjustment can adequately prevent spurious associations between genetic variants and CF-related traits and are therefore effective tools to control for population structure even when population structure is confounded with disease severity and a common pathogenic variant.

Association of Neonatal Hospital Length of Stay with Lung Function in Primary Ciliary Dyskinesia.

Rationale: Primary ciliary dyskinesia (PCD), an inherited lung disease, is characterized by abnormal ciliary function leading to progressive bronchiectasis. There is wide variability in respiratory disease severity at birth and later in life. Objectives: To evaluate the association between neonatal hospital length of stay (neonatal-LOS) and supplemental oxygen duration (SuppO2) with lung function in pediatric PCD. We hypothesized that longer neonatal-LOS and SuppO2 are associated with worse lung function (i.e., forced expiratory volume in 1 second percent predicted [FEV1pp]). Methods: We performed a secondary analysis of the Genetic Disorders of Mucociliary Clearance Consortium prospective longitudinal multicenter cohort study. Participants enrolled, during 2006-2011, were <19 years old with a confirmed PCD diagnosis and followed annually for 5 years. The exposure variables were neonatal-LOS and SuppO2, counted in days since birth. The outcome, FEV1pp, was measured annually by spirometry. The associations of neonatal-LOS and SuppO2 with FEV1pp were evaluated with a linear mixed-effects model with repeated measures and random intercepts, adjusted for age and ciliary ultrastructural defects. Results: Included were 123 participants (male, 47%; mean enrollment age, 8.3 yr [range, 0 to 18 yr]) with 578 visits (median follow-up, 5 yr). The median neonatal-LOS was 9 d (range, 1 to 90 d), and median SuppO2 was 5 d (range, 0 to 180 d). Neonatal-LOS was associated with worse lung function (-0.27 FEV1pp/d [95% confidence interval, -0.53 to -0.01]; P = 0.04). SuppO2 was not associated with lung function. Conclusions: Neonatal-LOS is associated with worse lung function in pediatric PCD, independent of age and ultrastructural defects. Future research on the mechanisms of neonatal respiratory distress and its management may help us understand the variability of lung health outcomes in PCD.

The role of SPAG1 in the assembly of axonemal dyneins in human airway epithelia.

Mutations in SPAG1, a dynein axonemal assembly factor (DNAAF) that facilitates the assembly of dynein arms in the cytoplasm before their transport into the cilium, result in primary ciliary dyskinesia (PCD), a genetically heterogenous disorder characterized by chronic oto-sino-pulmonary disease, infertility and laterality defects. To further elucidate the role of SPAG1 in dynein assembly, we examined its expression, interactions and ciliary defects in control and PCD human airway epithelia. Immunoprecipitations showed that SPAG1 interacts with multiple DNAAFs, dynein chains and canonical components of the R2TP complex. Protein levels of dynein heavy chains (DHCs) and interactions between DHCs and dynein intermediate chains (DICs) were reduced in SPAG1 mutants. We also identified a previously uncharacterized 60 kDa SPAG1 isoform, through examination of PCD subjects with an atypical ultrastructural defect for SPAG1 variants, that can partially compensate for the absence of full-length SPAG1 to assemble a reduced number of outer dynein arms. In summary, our data show that SPAG1 is necessary for axonemal dynein arm assembly by scaffolding R2TP-like complexes composed of several DNAAFs that facilitate the folding and/or binding of the DHCs to the DIC complex.