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/).

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.

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.

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.

Leveraging Open Electronic Health Record Data and Environmental Exposures Data to Derive Insights Into Rare Pulmonary Disease.

Research on rare diseases has received increasing attention, in part due to the realized profitability of orphan drugs. Biomedical informatics holds promise in accelerating translational research on rare disease, yet challenges remain, including the lack of diagnostic codes for rare diseases and privacy concerns that prevent research access to electronic health records when few patients exist. The Integrated Clinical and Environmental Exposures Service (ICEES) provides regulatory-compliant open access to electronic health record data that have been integrated with environmental exposures data, as well as analytic tools to explore the integrated data. We describe a proof-of-concept application of ICEES to examine demographics, clinical characteristics, environmental exposures, and health outcomes among a cohort of patients enriched for phenotypes associated with cystic fibrosis (CF), idiopathic bronchiectasis (IB), and primary ciliary dyskinesia (PCD). We then focus on a subset of patients with CF, leveraging the availability of a diagnostic code for CF and serving as a benchmark for our development work. We use ICEES to examine select demographics, co-diagnoses, and environmental exposures that may contribute to poor health outcomes among patients with CF, defined as emergency department or inpatient visits for respiratory issues. We replicate current understanding of the pathogenesis and clinical manifestations of CF by identifying co-diagnoses of asthma, chronic nasal congestion, cough, middle ear disease, and pneumonia as factors that differentiate patients with poor health outcomes from those with better health outcomes. We conclude by discussing our preliminary findings in relation to other published work, the strengths and limitations of our approach, and our future directions.

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.

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.

Leveraging TOPMed imputation server and constructing a cohort-specific imputation reference panel to enhance genotype imputation among cystic fibrosis patients.

Cystic fibrosis (CF) is a severe genetic disorder that can cause multiple comorbidities affecting the lungs, the pancreas, the luminal digestive system and beyond. In our previous genome-wide association studies (GWAS), we genotyped approximately 8,000 CF samples using a mixture of different genotyping platforms. More recently, the Cystic Fibrosis Genome Project (CFGP) performed deep (approximately 30×) whole genome sequencing (WGS) of 5,095 samples to better understand the genetic mechanisms underlying clinical heterogeneity among patients with CF. For mixtures of GWAS array and WGS data, genotype imputation has proven effective in increasing effective sample size. Therefore, we first performed imputation for the approximately 8,000 CF samples with GWAS array genotype using the Trans-Omics for Precision Medicine (TOPMed) freeze 8 reference panel. Our results demonstrate that TOPMed can provide high-quality imputation for patients with CF, boosting genomic coverage from approximately 0.3-4.2 million genotyped markers to approximately 11-43 million well-imputed markers, and significantly improving polygenic risk score (PRS) prediction accuracy. Furthermore, we built a CF-specific CFGP reference panel based on WGS data of patients with CF. We demonstrate that despite having approximately 3% the sample size of TOPMed, our CFGP reference panel can still outperform TOPMed when imputing some CF disease-causing variants, likely owing to allele and haplotype differences between patients with CF and general populations. We anticipate our imputed data for 4,656 samples without WGS data will benefit our subsequent genetic association studies, and the CFGP reference panel built from CF WGS samples will benefit other investigators studying CF.

Expression of a Truncated Form of ODAD1 Associated with an Unusually Mild Primary Ciliary Dyskinesia Phenotype.

Primary ciliary dyskinesia (PCD) is a rare lung disease caused by mutations that impair the function of motile cilia, resulting in chronic upper and lower respiratory disease, reduced fertility, and a high prevalence of situs abnormalities. The disease is genetically and phenotypically heterogeneous, with causative mutations in > 50 genes identified, and clinical phenotypes ranging from mild to severe. Absence of ODAD1 (CCDC114), a component of the outer dynein arm docking complex, results in a failure to assemble outer dynein arms (ODAs), mostly immotile cilia, and a typical PCD phenotype. We identified a female (now 34 years old) with an unusually mild clinical phenotype who has a homozygous non-canonical splice mutation (c.1502+5G>A) in ODAD1. To investigate the mechanism for the unusual phenotype, we performed molecular and functional studies of cultured nasal epithelial cells. We demonstrate that this splice mutation results in the expression of a truncated protein that is attached to the axoneme, indicating that the mutant protein retains partial function. This allows for the assembly of some ODAs and a significant level of ciliary activity that may result in the atypically mild clinical phenotype. The results also suggest that partial restoration of ciliary function by therapeutic agents could lead to significant improvement of disease symptoms.

Hereditary Mucin Deficiency Caused by Biallelic Loss of Function of MUC5B.

Rationale: Mucin homeostasis is fundamental to airway health. Upregulation of airway mucus glycoprotein MUC5B is observed in diverse common lung diseases and represents a potential therapeutic target. In mice, Muc5b is required for mucociliary clearance and for controlling inflammation after microbial exposure. The consequences of its loss in humans are unclear. Objectives: The goal of this study was to identify and characterize a family with congenital absence of MUC5B protein. Methods: We performed whole-genome sequencing in an adult proband with unexplained bronchiectasis, impaired pulmonary function, and repeated Staphylococcus aureus infection. Deep phenotyping over a 12-year period included assessments of pulmonary radioaerosol mucociliary clearance. Genotyping with reverse phenotyping was organized for eight family members. Extensive experiments, including immunofluorescence staining and mass spectrometry for mucins, were performed across accessible sample types. Measurements and Main Results: The proband, and her symptomatic sibling who also had extensive sinus disease with nasal polyps, were homozygous for a novel splicing variant in the MUC5B gene (NM_002458.2: c.1938 + 1G>A). MUC5B was absent from saliva, sputum, and nasal samples. Mucociliary clearance was impaired in the proband, and large numbers of apoptotic macrophages were present in sputum. Three siblings heterozygous for the familial MUC5B variant were asymptomatic but had a shared pattern of mild lung function impairments. Conclusions: Congenital absence of MUC5B defines a new category of genetic respiratory disease. The human phenotype is highly concordant with that of the Muc5b-/- murine model. Further study of individuals with decreased MUC5B production could provide unique mechanistic insights into airway mucus biology.

Otolaryngology Manifestations of Primary Ciliary Dyskinesia: A Multicenter Study.

OBJECTIVE: This project aims to prospectively and objectively assess otolaryngological manifestations and quality of life of children with primary ciliary dyskinesia (PCD) and compare these findings with healthy pediatric controls. STUDY DESIGN: Cross-sectional. SETTING: Two high-volume pediatric PCD specialty centers. METHODS: Standardized clinical assessment; Sino-Nasal Outcome Test 22 (SNOT-22); Hearing Environment and Reflection Quality of Life (HEAR-QL); Reflux Symptom Index (RSI); standardized physical examination of the sinonasal, laryngeal, and otological systems; and investigations including pure-tone audiograms (PTAs) and sinonasal cultures were collected. RESULTS: Forty-seven children with PCD and 25 control participants were recruited. Children with PCD had more upper airway symptoms than healthy children. They had significantly higher scores in both SNOT-22 and RSI, indicating worse sinonasal and reflux symptoms, with worse quality of life on the HEAR-QL index compared to healthy children (P < .05). Fifty-two percent of children with PCD-related hearing loss were not aware of their hearing deficit that was present on audiological assessment, and only 23% of children who had ventilation tubes had chronic otorrhea, most of which was easily controlled with ototopic drops. Furthermore, although all children with PCD had chronic rhinosinusitis, only 36% of them were using topical nasal treatment. The most common bacteria cultured from the middle meatus were Staphylococcus aureus in 11 of 47 (23%), followed by Streptococcus pneumoniae in 10 of 47 (21%). CONCLUSION: This multisite cohort highlights the importance of otolaryngology involvement in the management of children with PCD. More rigorous otolaryngological management may lead to reductions in overall morbidity and improve quality of life for children with PCD.

Comparing encounter-based and annualized chronic pseudomonas infection definitions in cystic fibrosis.

Chronic Pseudomonas aeruginosa (Pa) infection is associated with increased morbidity and mortality in people with cystic fibrosis (CF). There is no gold standard definition of chronic Pa infection in CF. We compared chronic Pa definitions using encounter-based versus annualized data in the Early Pseudomonas Infection Control (EPIC) Observational study cohort, and subsequently compared annualized chronic Pa definitions across a range of U.S. cohorts spanning decades of CF care. We found that an annualized chronic Pa definition requiring at least 1 Pa+ culture in 3 of 4 consecutive years ("Green 3/4") resulted in chronic Pa metrics similar to established encounter-based modified Leeds criteria definitions, including a similar age at and proportion who fulfilled chronic Pa criteria, and a similar proportion with sustained Pa infection after meeting the chronic Pa definition. The Green 3/4 chronic Pa definition will be valuable for longitudinal analyses in cohorts with limited culture frequency.

Complete CFTR gene sequencing in 5,058 individuals with cystic fibrosis informs variant-specific treatment.

BACKGROUND: Cystic fibrosis (CF) is a recessive condition caused by variants in each CF transmembrane conductance regulator (CFTR) allele. Clinically affected individuals without two identified causal variants typically have no further interrogation of CFTR beyond examination of coding regions, but the development of variant-specific CFTR-targeted treatments necessitates complete understanding of CFTR genotype. METHODS: Whole genome sequences were analyzed on 5,058 individuals with CF. We focused on the full CFTR gene sequence and identified disease-causing variants in three phases: screening for known and structural variants; discovery of novel loss-of-function variants; and investigation of remaining variants. RESULTS: All variants identified in the first two phases and coding region variants found in the third phase were interpreted according to CFTR2 or ACMG criteria (n = 371; 16 [4.3%] previously unreported). Full gene sequencing enabled delineation of 18 structural variants (large insertions or deletions), of which two were novel. Additional CFTR variants of uncertain effect were found in 76 F508del homozygotes and in 21 individuals with other combinations of CF-causing variants. Both causative variants were identified in 98.1% (n = 4,960) of subjects, an increase of 2.3 percentage points from the 95.8% (n = 4,847) who had a registry- or chart-reported disease-causing CFTR genotype. Of the remaining 98 individuals, 78 carried one variant that has been associated with CF (CF-causing [n = 70] or resulting in varying clinical consequences n = 8]). CONCLUSIONS: Complete CFTR gene sequencing in 5,058 individuals with CF identified at least one DNA variant in 99.6% of the cohort that is targetable by current molecular or emerging gene-based therapeutic technologies.