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.

Use of adenine base editing and homology-independent targeted integration strategies to correct the cystic fibrosis causing variant, W1282X.

Small molecule drugs known as modulators can treat ~90% of people with cystic fibrosis (CF), but do not work for premature termination codon variants such as W1282X (c.3846G>A). Here we evaluated two gene editing strategies, Adenine Base Editing (ABE) to correct W1282X, and Homology-Independent Targeted Integration (HITI) of a CFTR superexon comprising exons 23-27 (SE23-27) to enable expression of a CFTR mRNA without W1282X. In Flp-In-293 cells stably expressing a CFTR expression minigene bearing W1282X, ABE corrected 24% of W1282X alleles, rescued CFTR mRNA from nonsense mediated decay and restored protein expression. However, bystander editing at the adjacent adenine (c.3847A>G), caused an amino acid change (R1283G) that affects CFTR maturation and ablates ion channel activity. In primary human nasal epithelial cells homozygous for W1282X, ABE corrected 27% of alleles, but with a notably lower level of bystander editing, and CFTR channel function was restored to 16% of wild-type levels. Using the HITI approach, correct integration of a SE23-27 in intron 22 of the CFTR locus in 16HBEge W1282X cells was detected in 5.8% of alleles, resulting in 7.8% of CFTR transcripts containing the SE23-27 sequence. Analysis of a clonal line homozygous for the HITI-SE23-27 produced full-length mature protein and restored CFTR anion channel activity to 10% of wild-type levels, which could be increased three-fold upon treatment with the triple combination of CF modulators. Overall, these data demonstrate two different editing strategies can successfully correct W1282X, the second most common class I variant, with a concomitant restoration of CFTR function.

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.

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.

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.

Investigation of CFTR Function in Human Nasal Epithelial Cells Informs Personalized Medicine.

We broaden the clinical versatility of human nasal epithelial (HNE) cells. HNEs were isolated from 10 participants harboring CFTR variants: nine with rare variants (Q359R [n=2], G480S, R334W [n=5], and R560T) and one person harboring R117H;7T;TG10/5T;TG12. Cultures were differentiated at air-liquid interface. CFTR function was measured in Ussing chambers at three conditions - baseline, ivacaftor, and elexacaftor+tezacaftor+ivacaftor (ETI). Four participants initiated modulators. Q359R HNEs had 5.4% (%WT) baseline CFTR function and 25.5% with ivacaftor. With therapy, sweat [Cl-] decreased and symptoms resolved. G480S HNEs had 4.1% baseline and 32.1% CFTR function with ETI. Clinically, FEV1 increased and sweat [Cl-] decreased (119 to 46mmol/L) with ETI. In vitro cultures derived from five individuals harboring R334W showed a moderate increase in CFTR function with exposure to modulators. For one of these participants, ETI was begun in vivo; symptoms and FEV1 improved. c.1679G>C (R560T) HNEs had <4% baseline CFTR function and no modulator response. RNA analysis confirmed that c.1679G>C completely mis-splices. A symptomatic patient harboring R117H;7T;TG10/5T;TG12 exhibited reduced CFTR function (17.5%) in HNEs, facilitating mild CF diagnosis. HNEs responded to modulators (ivacaftor: 32.8%, ETI: 55.5%) and, since beginning therapy, lung function improved. While reaffirming HNE use for guiding therapeutic approaches, we inform predictions on modulator response (e.g. R334W) and closely assess variants affecting splicing (e.g. c.1679G>C). Notably, functional studies in HNEs harboring R117H;7T;TG10/5T;TG12 facilitated mild CF diagnosis, suggesting use for HNE functional studies as a clinical diagnostic test.

Two rare variants that affect the same amino acid in CFTR have distinct responses to ivacaftor.

Some residues in the cystic fibrosis transmembrane conductance regulator (CFTR) channel are the site of more than one CFTR variant that cause cystic fibrosis. Here, we investigated the function of S1159F and S1159P, two variants associated with different clinical phenotypes, which affect the same pore-lining residue in transmembrane segment 12 that are both strongly potentiated by ivacaftor when expressed in CFBE41o- bronchial epithelial cells. To study the single-channel behaviour of CFTR, we applied the patch-clamp technique to Chinese hamster ovary cells heterologously expressing CFTR variants incubated at 27°C to enhance channel residence at the plasma membrane. S1159F- and S1159P-CFTR formed Cl- channels activated by cAMP-dependent phosphorylation and gated by ATP that exhibited thermostability at 37°C. Both variants modestly reduced the single-channel conductance of CFTR. By severely attenuating channel gating, S1159F- and S1159P-CFTR reduced the open probability (Po ) of wild-type CFTR by ≥75% at ATP (1 mM); S1159F-CFTR caused the greater decrease in Po consistent with its more severe clinical phenotype. Ivacaftor (10-100 nM) doubled the Po of both CFTR variants without restoring Po values to wild-type levels, but concomitantly, ivacaftor decreased current flow through open channels. For S1159F-CFTR, the reduction of current flow was marked at high (supersaturated) ivacaftor concentrations (0.5-1 µM) and voltage-independent, identifying an additional detrimental action of elevated ivacaftor concentrations. In conclusion, S1159F and S1159P are gating variants, which also affect CFTR processing and conduction, but not stability, necessitating the use of combinations of CFTR modulators to optimally restore their channel activity. KEY POINTS: Dysfunction of the ion channel cystic fibrosis transmembrane conductance regulator (CFTR) causes the genetic disease cystic fibrosis (CF). This study investigated two rare pathogenic CFTR variants, S1159F and S1159P, which affect the same amino acid in CFTR, to understand the molecular basis of disease and response to the CFTR-targeted therapy ivacaftor. Both rare variants diminished CFTR function by modestly reducing current flow through the channel and severely inhibiting ATP-dependent channel gating with S1159F exerting the stronger adverse effect, which correlates with its association with more severe disease. Ivacaftor potentiated channel gating by both rare variants without restoring their activity to wild-type levels, but concurrently reduced current flow through open channels, particularly those of S1159F-CFTR. Our data demonstrate that S1159F and S1159P cause CFTR dysfunction by multiple mechanisms that require combinations of CFTR-targeted therapies to fully restore channel function.

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.

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

Elucidating the functional consequence of molecular defects underlying genetic diseases enables appropriate design of therapeutic options. Treatment of cystic fibrosis (CF) is an exemplar of this paradigm as the development of CFTR modulator therapies has allowed for targeted and effective treatment of individuals harboring specific genetic variants. However, the mechanism of these drugs limits effectiveness to particular classes of variants that allow production of CFTR protein. Thus, assessment of the molecular mechanism of individual variants is imperative for proper assignment of these precision therapies. This is particularly important when considering variants that affect pre-mRNA splicing, thus limiting success of the existing protein-targeted therapies. Variants affecting splicing can occur throughout exons and introns and the complexity of the process of splicing lends itself to a variety of outcomes, both at the RNA and protein levels, further complicating assessment of disease liability and modulator response. To investigate the scope of this challenge, we evaluated splicing and downstream effects of 52 naturally occurring CFTR variants (exonic = 15, intronic = 37). Expression of constructs containing select CFTR intronic sequences and complete CFTR exonic sequences in cell line models allowed for assessment of RNA and protein-level effects on an allele by allele basis. Characterization of primary nasal epithelial cells obtained from individuals harboring splice variants corroborated in vitro data. Notably, we identified exonic variants that result in complete missplicing and thus a lack of modulator response (e.g. c.2908G>A, c.523A>G), as well as intronic variants that respond to modulators due to the presence of residual normally spliced transcript (e.g. c.4242+2T>C, c.3717+40A>G). Overall, our data reveals diverse molecular outcomes amongst both exonic and intronic variants emphasizing the need to delineate RNA, protein, and functional effects of each variant in order to accurately assign precision therapies.

C FTR variants are associated with chronic bronchitis in smokers.

INTRODUCTION: Loss-of-function variants in both copies of the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause cystic fibrosis (CF); however, there is evidence that reduction in CFTR function due to the presence of one deleterious variant can have clinical consequences. Here, we hypothesise that CFTR variants in individuals with a history of smoking are associated with chronic obstructive pulmonary disease (COPD) and related phenotypes. METHODS: Whole-genome sequencing was performed through the National Heart, Lung, and Blood Institute TOPMed (TransOmics in Precision Medicine) programme in 8597 subjects from the COPDGene (Genetic Epidemiology of COPD) study, an observational study of current and former smokers. We extracted clinically annotated CFTR variants and performed single-variant and variant-set testing for COPD and related phenotypes. Replication was performed in 2118 subjects from the ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) study. RESULTS: We identified 301 coding variants within the CFTR gene boundary: 147 of these have been reported in individuals with CF, including 36 CF-causing variants. We found that CF-causing variants were associated with chronic bronchitis in variant-set testing in COPDGene (one-sided p=0.0025; OR 1.53) and in meta-analysis of COPDGene and ECLIPSE (one-sided p=0.0060; OR 1.52). Single-variant testing revealed that the F508del variant was associated with chronic bronchitis in COPDGene (one-sided p=0.015; OR 1.47). In addition, we identified 32 subjects with two or more CFTR variants on separate alleles and these subjects were enriched for COPD cases (p=0.010). CONCLUSIONS: Cigarette smokers who carry one deleterious CFTR variant have higher rates of chronic bronchitis, while presence of two CFTR variants may be associated with COPD. These results indicate that genetically mediated reduction in CFTR function contributes to COPD related phenotypes, in particular chronic bronchitis.

Accurate assignment of disease liability to genetic variants using only population data.

PURPOSE: The growing size of public variant repositories prompted us to test the accuracy of pathogenicity prediction of DNA variants using population data alone. METHODS: Under the a priori assumption that the ratio of the prevalence of variants in healthy population vs that in affected populations form 2 distinct distributions (pathogenic and benign), we used a Bayesian method to assign probability to a variant belonging to either distribution. RESULTS: The approach, termed Bayesian prevalence ratio (BayPR), accurately parsed 300 of 313 expertly curated CFTR variants: 284 of 296 pathogenic/likely pathogenic variants in 1 distribution and 16 of 17 benign/likely benign variants in another. BayPR produced an area under the receiver operating characteristic curve of 0.99 for 103 functionally confirmed missense CFTR variants, which is equal to or exceeds 10 commonly used algorithms (area under the receiver operating characteristic curve range = 0.54-0.99). Application of BayPR to expertly curated variants in 8 genes associated with 7 Mendelian conditions led to the assignment of a disease-causing probability of ≥80% to 1350 of 1374 (98.3%) pathogenic/likely pathogenic variants and of ≤20% to 22 of 23 (95.7%) benign/likely benign variants. CONCLUSION: Irrespective of the variant type or functional effect, the BayPR approach provides probabilities of pathogenicity for DNA variants responsible for Mendelian disorders using only the variant counts in affected and unaffected population samples.

Characterizing mucociliary clearance in young children with cystic fibrosis.

BACKGROUND: This research characterized mucociliary clearance (MCC) in young children with cystic fibrosis (CF). METHODS: Fourteen children (5-7 years old) with CF underwent: two baseline MCC measurements (Visits 1 and 2); one MCC measurement approximately 1 year later (Visit 3); and measurements of lung clearance index (LCI), a measure of ventilation inhomogeneity. RESULTS: Median (range) percent MCC through 60 min (MCC60) was similar on Visits 1 and 2 with 11.0 (0.9-33.7) and 12.8 (2.7-26.8), respectively (p = 0.95), and reproducible (Spearman Rho = 0.69; p = 0.007). Mucociliary clearance did not change significantly over 1 year with median percent MCC60 on Visit 3 [12.8 (3.7-17.6)] similar to Visit 2 (p = 0.58). Lower percent MCC60 on Visit 3 was significantly associated with higher LCI scores on Visit 3 (N = 14; Spearman Rho = -0.56; p = 0.04). CONCLUSIONS: Tests of MCC were reproducible and reliable over a 2-week period and stable over a 1-year period in 5-7-year-old children with CF. Lower MCC values were associated with increased ventilation inhomogeneity. These results suggest that measurements of MCC could be used in short-term clinical trials of interventions designed to modulate MCC and as a new, non-invasive test to evaluate early lung pathology in children with CF. IMPACT: This is the first study to characterize mucociliary clearance (MCC) in children with cystic fibrosis (CF) who were 5-7 years old. Measurements of mucociliary clearance were reproducible and reliable over a 2-week period and stable over a 1-year period. Variability in MCC between children was associated with differences in ventilation homogeneity, such that children with lower MCC values had increased ventilation inhomogeneity. These results suggest that measurements of MCC could be used in short-term clinical trials of interventions designed to modulate MCC and as a new, non-invasive test to evaluate early lung pathology in children with CF.

Genetic association and transcriptome integration identify contributing genes and tissues at cystic fibrosis modifier loci.

Cystic Fibrosis (CF) exhibits morbidity in several organs, including progressive lung disease in all patients and intestinal obstruction at birth (meconium ileus) in ~15%. Individuals with the same causal CFTR mutations show variable disease presentation which is partly attributed to modifier genes. With >6,500 participants from the International CF Gene Modifier Consortium, genome-wide association investigation identified a new modifier locus for meconium ileus encompassing ATP12A on chromosome 13 (min p = 3.83x10(-10)); replicated loci encompassing SLC6A14 on chromosome X and SLC26A9 on chromosome 1, (min p<2.2x10(-16), 2.81x10(-11), respectively); and replicated a suggestive locus on chromosome 7 near PRSS1 (min p = 2.55x10(-7)). PRSS1 is exclusively expressed in the exocrine pancreas and was previously associated with non-CF pancreatitis with functional characterization demonstrating impact on PRSS1 gene expression. We thus asked whether the other meconium ileus modifier loci impact gene expression and in which organ. We developed and applied a colocalization framework called the Simple Sum (SS) that integrates regulatory and genetic association information, and also contrasts colocalization evidence across tissues or genes. The associated modifier loci colocalized with expression quantitative trait loci (eQTLs) for ATP12A (p = 3.35x10(-8)), SLC6A14 (p = 1.12x10(-10)) and SLC26A9 (p = 4.48x10(-5)) in the pancreas, even though meconium ileus manifests in the intestine. The meconium ileus susceptibility locus on chromosome X appeared shifted in location from a previously identified locus for CF lung disease severity. Using the SS we integrated the lung disease association locus with eQTLs from nasal epithelia of 63 CF participants and demonstrated evidence of colocalization with airway-specific regulation of SLC6A14 (p = 2.3x10(-4)). Cystic Fibrosis is realizing the promise of personalized medicine, and identification of the contributing organ and understanding of tissue specificity for a gene modifier is essential for the next phase of personalizing therapeutic strategies.

Verifying nomenclature of DNA variants in submitted manuscripts: Guidance for journals.

Documenting variation in our genomes is important for research and clinical care. Accuracy in the description of DNA variants is therefore essential. To address this issue, the Human Variome Project convened a committee to evaluate the feasibility of requiring authors to verify that all variants submitted for publication complied with a widely accepted standard for description. After a pilot study of two journals, the committee agreed that requiring authors to verify that variants complied with Human Genome Variation Society nomenclature is a reasonable step toward standardizing the worldwide inventory of human variation.

Functional Assays Are Essential for Interpretation of Missense Variants Associated with Variable Expressivity.

Missense DNA variants have variable effects upon protein function. Consequently, interpreting their pathogenicity is challenging, especially when they are associated with disease variability. To determine the degree to which functional assays inform interpretation, we analyzed 48 CFTR missense variants associated with variable expressivity of cystic fibrosis (CF). We assessed function in a native isogenic context by evaluating CFTR mutants that were stably expressed in the genome of a human airway cell line devoid of endogenous CFTR expression. 21 of 29 variants associated with full expressivity of the CF phenotype generated <10% wild-type CFTR (WT-CFTR) function, a conservative threshold for the development of life-limiting CF lung disease, and five variants had moderately decreased function (10% to ∼25% WT-CFTR). The remaining three variants in this group unexpectedly had >25% WT-CFTR function; two were higher than 75% WT-CFTR. As expected, 14 of 19 variants associated with partial expressivity of CF had >25% WT-CFTR function; however, four had minimal to no effect on CFTR function (>75% WT-CFTR). Thus, 6 of 48 (13%) missense variants believed to be disease causing did not alter CFTR function. Functional studies substantially refined pathogenicity assignment with expert annotation and criteria from the American College of Medical Genetics and Genomics and Association for Molecular Pathology. However, four algorithms (CADD, REVEL, SIFT, and PolyPhen-2) could not differentiate between variants that caused severe, moderate, or minimal reduction in function. In the setting of variable expressivity, these results indicate that functional assays are essential for accurate interpretation of missense variants and that current prediction tools should be used with caution.

Cystic fibrosis-related diabetes onset can be predicted using biomarkers measured at birth.

PURPOSE: Cystic fibrosis (CF), caused by pathogenic variants in the CF transmembrane conductance regulator (CFTR), affects multiple organs including the exocrine pancreas, which is a causal contributor to cystic fibrosis-related diabetes (CFRD). Untreated CFRD causes increased CF-related mortality whereas early detection can improve outcomes. METHODS: Using genetic and easily accessible clinical measures available at birth, we constructed a CFRD prediction model using the Canadian CF Gene Modifier Study (CGS; n = 1,958) and validated it in the French CF Gene Modifier Study (FGMS; n = 1,003). We investigated genetic variants shown to associate with CF disease severity across multiple organs in genome-wide association studies. RESULTS: The strongest predictors included sex, CFTR severity score, and several genetic variants including one annotated to PRSS1, which encodes cationic trypsinogen. The final model defined in the CGS shows excellent agreement when validated on the FGMS, and the risk classifier shows slightly better performance at predicting CFRD risk later in life in both studies. CONCLUSION: We demonstrated clinical utility by comparing CFRD prevalence rates between the top 10% of individuals with the highest risk and the bottom 10% with the lowest risk. A web-based application was developed to provide practitioners with patient-specific CFRD risk to guide CFRD monitoring and treatment.

Cystic fibrosis genetics: from molecular understanding to clinical application.

The availability of the human genome sequence and tools for interrogating individual genomes provide an unprecedented opportunity to apply genetics to medicine. Mendelian conditions, which are caused by dysfunction of a single gene, offer powerful examples that illustrate how genetics can provide insights into disease. Cystic fibrosis, one of the more common lethal autosomal recessive Mendelian disorders, is presented here as an example. Recent progress in elucidating disease mechanism and causes of phenotypic variation, as well as in the development of treatments, demonstrates that genetics continues to play an important part in cystic fibrosis research 25 years after the discovery of the disease-causing gene.

Integrated Transcriptomic and Proteomic Analysis of Human Eccrine Sweat Glands Identifies Missing and Novel Proteins.

The eccrine sweat gland is an exocrine gland that is involved in the secretion of sweat for control of temperature. Malfunction of the sweat glands can result in disorders such as miliaria, hyperhidrosis and bromhidrosis. Understanding the transcriptome and proteome of sweat glands is important for understanding their physiology and role in diseases. However, no systematic transcriptome or proteome analysis of sweat glands has yet been reported. Here, we isolated eccrine sweat glands from human skin by microdissection and performed RNA-seq and proteome analysis. In total, ∼138,000 transcripts and ∼6,100 proteins were identified. Comparison of the RNA-seq data of eccrine sweat glands to other human tissues revealed the closest resemblance to the cortex region of kidneys. The proteome data showed enrichment of proteins involved in secretion, reabsorption, and wound healing. Importantly, protein level identification of the calcium ion channel TRPV4 suggests the importance of eccrine sweat glands in re-epithelialization of wounds and prevention of dehydration. We also identified 2 previously missing proteins from our analysis. Using a proteogenomic approach, we identified 7 peptides from 5 novel genes, which we validated using synthetic peptides. Most of the novel proteins were from short open reading frames (sORFs) suggesting that many sORFs still remain to be annotated in the human genome. This study presents the first integrated analysis of the transcriptome and proteome of the human eccrine sweat gland and would become a valuable resource for studying sweat glands in physiology and disease.

Potentially lethal cystic fibrosis gene variant in the orangutan.

A syndrome of chronic upper and lower airway disease leading to increased morbidity and mortality occurs primarily in captive orangutans. Similarities in symptoms to the inherited human respiratory disease, cystic fibrosis, led us to hypothesize that orangutan respiratory disease is a result of variants in the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene. We identified the nonsense variant, c.484A>T (p.Lys162X), in heterozygosity in an unaffected orangutan. Analysis of the pedigree of this orangutan confirmed that both his sire and deceased fetus also harbored the c.484A>T allele. An expression minigene harboring c.484A>T produced no full-length CFTR protein in HEK293 cells. Finally, the c.484A>T CFTR messenger RNA abundance was severely reduced in primary nasal epithelial cells of the orangutan indicating that c.484A>T (p.Lys162X) is potentially lethal. Genetic screening of the captive orangutan population could be used to prevent transmission of this potentially lethal variant, and thus aid in the conservation of this critically endangered species.

Capitalizing on the heterogeneous effects of CFTR nonsense and frameshift variants to inform therapeutic strategy for cystic fibrosis.

CFTR modulators have revolutionized the treatment of individuals with cystic fibrosis (CF) by improving the function of existing protein. Unfortunately, almost half of the disease-causing variants in CFTR are predicted to introduce premature termination codons (PTC) thereby causing absence of full-length CFTR protein. We hypothesized that a subset of nonsense and frameshift variants in CFTR allow expression of truncated protein that might respond to FDA-approved CFTR modulators. To address this concept, we selected 26 PTC-generating variants from four regions of CFTR and determined their consequences on CFTR mRNA, protein and function using intron-containing minigenes expressed in 3 cell lines (HEK293, MDCK and CFBE41o-) and patient-derived conditionally reprogrammed primary nasal epithelial cells. The PTC-generating variants fell into five groups based on RNA and protein effects. Group A (reduced mRNA, immature (core glycosylated) protein, function <1% (n = 5)) and Group B (normal mRNA, immature protein, function <1% (n = 10)) variants were unresponsive to modulator treatment. However, Group C (normal mRNA, mature (fully glycosylated) protein, function >1% (n = 5)), Group D (reduced mRNA, mature protein, function >1% (n = 5)) and Group E (aberrant RNA splicing, mature protein, function > 1% (n = 1)) variants responded to modulators. Increasing mRNA level by inhibition of NMD led to a significant amplification of modulator effect upon a Group D variant while response of a Group A variant was unaltered. Our work shows that PTC-generating variants should not be generalized as genetic 'nulls' as some may allow generation of protein that can be targeted to achieve clinical benefit.