Genetic modifiers of body mass index in individuals with cystic fibrosis.

To identify modifier loci underlying variation in body mass index (BMI) in persons with cystic fibrosis (pwCF), we performed a genome-wide association study (GWAS). Utilizing longitudinal height and weight data, along with demographic information and covariates from 4,393 pwCF, we calculated AvgBMIz representing the average of per-quarter BMI Z scores. The GWAS incorporated 9.8M single nucleotide polymorphisms (SNPs) with a minor allele frequency (MAF) > 0.005 extracted from whole-genome sequencing (WGS) of each study subject. We observed genome-wide significant association with a variant in FTO (FaT mass and Obesity-associated gene; rs28567725; p value = 1.21e-08; MAF = 0.41, ß = 0.106; n = 4,393 individuals) and a variant within ADAMTS5 (A Disintegrin And Metalloproteinase with ThromboSpondin motifs 5; rs162500; p value = 2.11e-10; MAF = 0.005, ß = -0.768; n = 4,085 pancreatic-insufficient individuals). Notably, BMI-associated variants in ADAMTS5 occur on a haplotype that is much more common in African (AFR, MAF = 0.183) than European (EUR, MAF = 0.006) populations (1000 Genomes project). A polygenic risk score (PRS) calculated using 924 SNPs (excluding 17 in FTO) showed significant association with AvgBMIz (p value = 2.2e-16; r2 = 0.03). Association between variants in FTO and the PRS correlation reveals similarities in the genetic architecture of BMI in CF and the general population. Inclusion of Black individuals in whom the single-gene disorder CF is much less common but genomic diversity is greater facilitated detection of association with variants that are in LD with functional SNPs in ADAMTS5. Our results illustrate the importance of population diversity, particularly when attempting to identify variants that manifest only under certain physiologic conditions.

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.

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.

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.

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.

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.

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.

Silencing XIST on the future active X: Searching human and bovine preimplantation embryos for the repressor.

X inactivation is the means of equalizing the dosage of X chromosomal genes in male and female eutherian mammals, so that only one X is active in each cell. The XIST locus (in cis) on each additional X chromosome initiates the transcriptional silence of that chromosome, making it an inactive X. How the active X in both males and females is protected from inactivation by its own XIST locus is not well understood in any mammal. Previous studies of autosomal duplications suggest that gene(s) on the short arm of human chromosome 19 repress XIST on the active X. Here, we examine the time of transcription of some candidate genes in preimplantation embryos using single-cell RNA sequencing data from human embryos and qRT-PCR from bovine embryos. The candidate genes assayed are those transcribed from 19p13.3-13.2, which are widely expressed and can remodel chromatin. Our results confirm that XIST is expressed at low levels from the future active X in embryos of both sexes; they also show that the XIST locus is repressed in both sexes when pluripotency factors are being upregulated, during the 4-8 cell and morula stages in human and bovine embryos - well before the early blastocyst (E5) when XIST on the inactive X in females starts to be upregulated. Our data suggest a role for DNMT1, UHRF1, SAFB and SAFB2 in XIST repression; they also exclude XACT and other 19p candidate genes and provide the transcriptional timing for some genes not previously assayed in human or bovine preimplantation embryos.

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.

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.

CFTR bearing variant p.Phe312del exhibits function inconsistent with phenotype and negligible response to ivacaftor.

The chloride channel dysfunction caused by deleterious cystic fibrosis transmembrane conductance regulator (CFTR) variants generally correlates with severity of cystic fibrosis (CF). However, 3 adults bearing the common severe variant p.Phe508del (legacy: F508del) and a deletion variant in an ivacaftor binding region of CFTR (p.Phe312del; legacy: F312del) manifested only elevated sweat chloride concentration (sw[Cl-]; 87-105 mEq/L). A database review of 25 individuals with F312del and a CF-causing variant revealed elevated sw[Cl-] (75-123 mEq/L) and variable CF features. F312del occurs at a higher-than-expected frequency in the general population, confirming that individuals with F312del and a CF-causing variant do not consistently develop overt CF features. In primary nasal cells, CFTR bearing F312del and F508del generated substantial chloride transport (66.0% ± 4.5% of WT-CFTR) but did not respond to ivacaftor. Single-channel analysis demonstrated that F312del did not affect current flow through CFTR, minimally altered gating, and ablated the ivacaftor response. When expressed stably in CF bronchial epithelial (CFBE41o-) cells, F312del-CFTR demonstrated residual function (50.9% ± 3.3% WT-CFTR) and a subtle decrease in forskolin response compared with WT-CFTR. F312del provides an exception to the established correlation between CFTR chloride transport and CF phenotype and informs our molecular understanding of ivacaftor response.

Increased expression of anion transporter SLC26A9 delays diabetes onset in cystic fibrosis.

Diabetes is a common complication of cystic fibrosis (CF) that affects approximately 20% of adolescents and 40%-50% of adults with CF. The age at onset of CF-related diabetes (CFRD) (marked by clinical diagnosis and treatment initiation) is an important measure of the disease process. DNA variants associated with age at onset of CFRD reside in and near SLC26A9. Deep sequencing of the SLC26A9 gene in 762 individuals with CF revealed that 2 common DNA haplotypes formed by the risk variants account for the association with diabetes. Single-cell RNA sequencing (scRNA-Seq) indicated that SLC26A9 is predominantly expressed in pancreatic ductal cells and frequently coexpressed with CF transmembrane conductance regulator (CFTR) along with transcription factors that have binding sites 5' of SLC26A9. These findings were replicated upon reanalysis of scRNA-Seq data from 4 independent studies. DNA fragments derived from the 5' region of SLC26A9-bearing variants from the low-risk haplotype generated 12%-20% higher levels of expression in PANC-1 and CFPAC-1 cells compared with the high- risk haplotype. Taken together, our findings indicate that an increase in SLC26A9 expression in ductal cells of the pancreas delays the age at onset of diabetes, suggesting a CFTR-agnostic treatment for a major complication of CF.

Cystic fibrosis transmembrane conductance regulator function, not TAS2R38 gene haplotypes, predict sinus surgery in children and young adults with cystic fibrosis.

BACKGROUND: Chronic rhinosinusitis symptomatology begins in early childhood individuals with cystic fibrosis (CF). Cystic fibrosis transmembrane conductance regulator (CFTR) function contributes to sinus development and disease. Genetic variants of the bitter taste receptor TAS2R38 have been suggested to contribute to sinus disease severity in individuals without CF. Our objective was to explore whether functional TAS2R38 haplotypes and CFTR function are associated with sinus disease or the need for sinus surgery in individuals with CF. METHODS: We conducted a retrospective study using prospectively collected data from the CF Twin-Sibling Study. The function of CFTR was assessed via chloride conductance. Genotyping of the TAS2R38 gene identified patients who were homozygous for the functional haplotype, heterozygous, or homozygous for nonfunctional haplotypes. Clustered multivariate logistic regression was performed, controlling for sex and family relationship. RESULTS: A total of 1291 patients were evaluated. Patients with ≤1% CFTR function were 1.56 times more likely to require sinus surgery than those with >1% CFTR function (p = 0.049). CFTR function did not correlate significantly with the presence of sinus disease (p = 0.30). In addition, there were no statistically significant differences in diagnosis of sinus disease or need for sinus surgery between patients with functional and nonfunctional TAS2R38 haplotypes. CONCLUSION: CFTR function correlates with need for sinus surgery, whereas TAS2R38 function does not appear to contribute to sinus disease or the need for sinus surgery in patients with CF.

Genetic Modifiers of Cystic Fibrosis-Related Diabetes Have Extensive Overlap With Type 2 Diabetes and Related Traits.

CONTEXT: Individuals with cystic fibrosis (CF) develop a distinct form of diabetes characterized by ß-cell dysfunction and islet amyloid accumulation similar to type 2 diabetes (T2D), but generally have normal insulin sensitivity. CF-related diabetes (CFRD) risk is determined by both CFTR, the gene responsible for CF, and other genetic variants. OBJECTIVE: To identify genetic modifiers of CFRD and determine the genetic overlap with other types of diabetes. DESIGN AND PATIENTS: A genome-wide association study was conducted for CFRD onset on 5740 individuals with CF. Weighted polygenic risk scores (PRSs) for type 1 diabetes (T1D), T2D, and diabetes endophenotypes were tested for association with CFRD. RESULTS: Genome-wide significance was obtained for variants at a novel locus (PTMA) and 2 known CFRD genetic modifiers (TCF7L2 and SLC26A9). PTMA and SLC26A9 variants were CF-specific; TCF7L2 variants also associated with T2D. CFRD was strongly associated with PRSs for T2D, insulin secretion, postchallenge glucose concentration, and fasting plasma glucose, and less strongly with T1D PRSs. CFRD was inconsistently associated with PRSs for insulin sensitivity and was not associated with a PRS for islet autoimmunity. A CFRD PRS comprising variants selected from these PRSs (with a false discovery rate < 0.1) and the genome-wide significant variants was associated with CFRD in a replication population. CONCLUSIONS: CFRD and T2D have more etiologic and mechanistic overlap than previously known, aligning along pathways involving ß-cell function rather than insulin sensitivity. Two CFRD risk loci are unrelated to T2D and may affect multiple aspects of CF. An 18-variant PRS stratifies risk of CFRD in an independent population.