Caution advised in the use of CFTR modulator treatment for individuals harboring specific CFTR variants.

In December 2020, the U.S. Food and Drug Administration (FDA) expanded the list of CFTR variants approved for treatment with CFTR modulators drugs from 39 to 183. Clinicians should be aware that individuals harboring certain variants approved for treatment may not respond to or benefit from this therapy. After review, the expert panel leading the CFTR2 project identified four categories of variants that may not result in a clinical response to modulator treatment: 15 variants assigned as non CF-causing; 45 variants of unknown significance; six variants known or suspected to cause mis-splicing as their primary defect rather than an amino acid substitution; and eight variants known to occur together in cis with another deleterious variant not expected to lead to CFTR protein (nonsense or frameshift). The potential risks and benefits of CFTR modulator therapy should be considered carefully for individuals harboring these variants.

Genetic evidence supports the development of SLC26A9 targeting therapies for the treatment of lung disease.

Over 400 variants in the cystic fibrosis (CF) transmembrane conductance regulator (CFTR) are CF-causing. CFTR modulators target variants to improve lung function, but marked variability in response exists and current therapies do not address all CF-causing variants highlighting unmet needs. Alternative epithelial ion channel/transporters such as SLC26A9 could compensate for CFTR dysfunction, providing therapeutic targets that may benefit all individuals with CF. We investigate the relationship between rs7512462, a marker of SLC26A9 activity, and lung function pre- and post-treatment with CFTR modulators in Canadian and US CF cohorts, in the general population, and in those with chronic obstructive pulmonary disease (COPD). Rs7512462 CC genotype is associated with greater lung function in CF individuals with minimal function variants (for which there are currently no approved therapies; p = 0.008); and for gating (p = 0.033) and p.Phe508del/ p.Phe508del (p = 0.006) genotypes upon treatment with CFTR modulators. In parallel, human nasal epithelia with CC and p.Phe508del/p.Phe508del after Ussing chamber analysis of a combination of approved and experimental modulator treatments show greater CFTR function (p = 0.0022). Beyond CF, rs7512462 is associated with peak expiratory flow in a meta-analysis of the UK Biobank and Spirometa Consortium (p = 2.74 × 10-44) and provides p = 0.0891 in an analysis of COPD case-control status in the UK Biobank defined by spirometry. These findings support SLC26A9 as a therapeutic target to improve lung function for all people with CF and in individuals with other obstructive lung diseases.

Expression of cystic fibrosis lung disease modifier genes in human airway models.

BACKGROUND: Variation in respiratory response to cystic fibrosis (CF) small molecule therapies is due in part to the contribution of CF lung disease modifier genes. Cultured human bronchial epithelia (HBE) is the gold standard respiratory model for assessing CF therapeutic efficacy but it is hard to access. Cultured human nasal epithelia (HNE) is proposed as a more accessible surrogate model but it is unknown whether the expression profile of the modifier genes are comparable between HNE and HBE which we assess here. METHODS: RNA-sequencing was conducted on paired cultured and fresh HNE and HBE (n = 71 samples) collected from 21 individuals with CF. Genome-wide gene expression was first compared between cultured and fresh cells and then between cultured HNE and HBE based on an equivalence testing procedure we implemented. The co-expression relationships of CFTR and CF lung disease modifier genes were compared between cultured HNE and HBE to determine equivalent interactions. RESULTS: The culturing process had little impact on the expression level of CF lung disease modifier genes. Over 90% of expressed genes showed significant equivalent expression level across cultured HNE and HBE (expression fold-change<2, FDR<0.1), including CFTR and CF lung disease modifier genes. The difference in co-expression relationships among these genes was not significant (p-value=0.99), suggesting their functional interactions are likely to be consistent in the two models. CONCLUSIONS: Cultured HNE recapitulates the expression profile of CF lung disease modifier genes in cultured HBE, suggesting the biological processes involving these genes are likely to be consistent across the two models.

Positive epistasis between disease-causing missense mutations and silent polymorphism with effect on mRNA translation velocity.

Epistasis refers to the dependence of a mutation on other mutation(s) and the genetic context in general. In the context of human disorders, epistasis complicates the spectrum of disease symptoms and has been proposed as a major contributor to variations in disease outcome. The nonadditive relationship between mutations and the lack of complete understanding of the underlying physiological effects limit our ability to predict phenotypic outcome. Here, we report positive epistasis between intragenic mutations in the cystic fibrosis transmembrane conductance regulator (CFTR)-the gene responsible for cystic fibrosis (CF) pathology. We identified a synonymous single-nucleotide polymorphism (sSNP) that is invariant for the CFTR amino acid sequence but inverts translation speed at the affected codon. This sSNP in cis exhibits positive epistatic effects on some CF disease-causing missense mutations. Individually, both mutations alter CFTR structure and function, yet when combined, they lead to enhanced protein expression and activity. The most robust effect was observed when the sSNP was present in combination with missense mutations that, along with the primary amino acid change, also alter the speed of translation at the affected codon. Functional studies revealed that synergistic alteration in ribosomal velocity is the underlying mechanism; alteration of translation speed likely increases the time window for establishing crucial domain-domain interactions that are otherwise perturbed by each individual mutation.

LocusFocus: Web-based colocalization for the annotation and functional follow-up of GWAS.

Genome-wide association studies (GWAS) have primarily identified trait-associated loci in the non-coding genome. Colocalization analyses of SNP associations from GWAS with expression quantitative trait loci (eQTL) evidence enable the generation of hypotheses about responsible mechanism, genes and tissues of origin to guide functional characterization. Here, we present a web-based colocalization browsing and testing tool named LocusFocus (https://locusfocus.research.sickkids.ca). LocusFocus formally tests colocalization using our established Simple Sum method to identify the most relevant genes and tissues for a particular GWAS locus in the presence of high linkage disequilibrium and/or allelic heterogeneity. We demonstrate the utility of LocusFocus, following up on a genome-wide significant locus from a GWAS of meconium ileus (an intestinal obstruction in cystic fibrosis). Using LocusFocus for colocalization analysis with eQTL data suggests variation in ATP12A gene expression in the pancreas rather than intestine is responsible for the GWAS locus. LocusFocus has no operating system dependencies and may be installed in a local web server. LocusFocus is available under the MIT license, with full documentation and source code accessible on GitHub at https://github.com/naim-panjwani/LocusFocus.

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.

Correlating Cystic Fibrosis Transmembrane Conductance Regulator Function with Clinical Features to Inform Precision Treatment of Cystic Fibrosis.

Rationale: The advent of precision treatment for cystic fibrosis using small-molecule therapeutics has created a need to estimate potential clinical improvements attributable to increases in cystic fibrosis transmembrane conductance regulator (CFTR) function. Objectives: To derive CFTR function of a variety of CFTR genotypes and correlate with key clinical features (sweat chloride concentration, pancreatic exocrine status, and lung function) to develop benchmarks for assessing response to CFTR modulators. Methods: CFTR function assigned to 226 unique CFTR genotypes was correlated with the clinical data of 54,671 individuals enrolled in the Clinical and Functional Translation of CFTR (CFTR2) project. Cross-sectional FEV1% predicted measurements were plotted by age at which measurement was obtained. Shifts in sweat chloride concentration and lung function reported in CFTR modulator trials were compared with function-phenotype correlations to assess potential efficacy of therapies. Measurements and Main Results: CFTR genotype function exhibited a logarithmic relationship with each clinical feature. Modest increases in CFTR function related to differing genotypes were associated with clinically relevant improvements in cross-sectional FEV1% predicted over a range of ages (6-82 yr). Therapeutic responses to modulators corresponded closely to predictions from the CFTR2-derived relationship between CFTR genotype function and phenotype. Conclusions: Increasing CFTR function in individuals with severe disease will have a proportionally greater effect on outcomes than similar increases in CFTR function in individuals with mild disease and should reverse a substantial fraction of the disease process. This study provides reference standards for clinical outcomes that may be achieved by increasing CFTR function.

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.

Improving imputation in disease-relevant regions: lessons from cystic fibrosis.

Does genotype imputation with public reference panels identify variants contributing to disease? Genotype imputation using the 1000 Genomes Project (1KG; 2504 individuals) displayed poor coverage at the causal cystic fibrosis (CF) transmembrane conductance regulator (CFTR) locus for the International CF Gene Modifier Consortium. Imputation with the larger Haplotype Reference Consortium (HRC; 32,470 individuals) displayed improved coverage but low sensitivity of variants clinically relevant for CF. A hybrid reference that combined whole genome sequencing (WGS) from 101 CF individuals with the 1KG imputed a greater number of single-nucleotide variants (SNVs) that would be analyzed in a genetic association study (r2 ≥ 0.3 and MAF ≥ 0.5%) than imputation with the HRC, while the HRC excelled in the lower frequency spectrum. Using the 1KG or HRC as reference panels missed the most common CF-causing variants or displayed low imputation accuracy. Designs that incorporate population-specific WGS can improve imputation accuracy at disease-specific loci, while imputation using public data sets can omit disease-relevant genotypes.

Mammographic density defined by higher than conventional brightness thresholds better predicts breast cancer risk.

Background: Mammographic density defined by the conventional pixel brightness threshold, and adjusted for age and body mass index (BMI), is a well-established risk factor for breast cancer. We asked if higher thresholds better separate women with and without breast cancer. Methods: We studied Australian women, 354 with breast cancer over-sampled for early-onset and family history, and 944 unaffected controls frequency-matched for age at mammogram. We measured mammographic dense area and percent density using the CUMULUS software at the conventional threshold, which we call Cumulus , and at two increasingly higher thresholds, which we call Altocumulus and Cirrocumulus , respectively. All measures were Box-Cox transformed and adjusted for age and BMI. We estimated the odds per adjusted standard deviation (OPERA) using logistic regression and the area under the receiver operating characteristic curve (AUC). Results: Altocumulus and Cirrocumulus were correlated with Cumulus (r ∼ 0.8 and 0.6 , respectively) . For dense area, the OPERA was 1.62, 1.74 and 1.73 for Cumulus, Altocumulus and Cirrocumulus , respectively (all P < 0.001). After adjusting for Altocumulus and Cirrocumulus , Cumulus was not significant ( P > 0.6). The OPERAs for percent density were less but gave similar findings. The mean of the standardized adjusted Altocumulus and Cirrocumulus dense area measures was the best predictor; OPERA = 1.87 [95% confidence interval (CI): 1.64-2.14] and AUC = 0.68 (0.65-0.71). Conclusions: The areas of higher mammographically dense regions are associated with almost 30% stronger breast cancer risk gradient, explain the risk association of the conventional measure and might be more aetiologically important. This has substantial implications for clinical translation and molecular, genetic and epidemiological research.

Sources of Variation in Sweat Chloride Measurements in Cystic Fibrosis.

RATIONALE: Expanding the use of cystic fibrosis transmembrane conductance regulator (CFTR) potentiators and correctors for the treatment of cystic fibrosis (CF) requires precise and accurate biomarkers. Sweat chloride concentration provides an in vivo assessment of CFTR function, but it is unknown the degree to which CFTR mutations account for sweat chloride variation. OBJECTIVES: To estimate potential sources of variation for sweat chloride measurements, including demographic factors, testing variability, recording biases, and CFTR genotype itself. METHODS: A total of 2,639 sweat chloride measurements were obtained in 1,761 twins/siblings from the CF Twin-Sibling Study, French CF Modifier Gene Study, and Canadian Consortium for Genetic Studies. Variance component estimation was performed by nested mixed modeling. MEASUREMENTS AND MAIN RESULTS: Across the tested CF population as a whole, CFTR gene mutations were found to be the primary determinant of sweat chloride variability (56.1% of variation) with contributions from variation over time (e.g., factors related to testing on different days; 13.8%), environmental factors (e.g., climate, family diet; 13.5%), other residual factors (e.g., test variability; 9.9%), and unique individual factors (e.g., modifier genes, unique exposures; 6.8%) (likelihood ratio test, P < 0.001). Twin analysis suggested that modifier genes did not play a significant role because the heritability estimate was negligible (H2 = 0; 95% confidence interval, 0.0-0.35). For an individual with CF, variation in sweat chloride was primarily caused by variation over time (58.1%) with the remainder attributable to residual/random factors (41.9%). CONCLUSIONS: Variation in the CFTR gene is the predominant cause of sweat chloride variation; most of the non-CFTR variation is caused by testing variability and unique environmental factors. If test precision and accuracy can be improved, sweat chloride measurement could be a valuable biomarker for assessing response to therapies directed at mutant CFTR.

Cystic fibrosis gene modifier SLC26A9 modulates airway response to CFTR-directed therapeutics.

Cystic fibrosis is realizing the promise of personalized medicine. Recent advances in drug development that target the causal CFTR directly result in lung function improvement, but variability in response is demanding better prediction of outcomes to improve management decisions. The genetic modifier SLC26A9 contributes to disease severity in the CF pancreas and intestine at birth and here we assess its relationship with disease severity and therapeutic response in the airways. SLC26A9 association with lung disease was assessed in individuals from the Canadian and French CF Gene Modifier consortia with CFTR-gating mutations and in those homozygous for the common Phe508del mutation. Variability in response to a CFTR-directed therapy attributed to SLC26A9 genotype was assessed in Canadian patients with gating mutations. A primary airway model system determined if SLC26A9 shows modification of Phe508del CFTR function upon treatment with a CFTR corrector. In those with gating mutations that retain cell surface-localized CFTR we show that SLC26A9 modifies lung function while this is not the case in individuals homozygous for Phe508del where cell surface expression is lacking. Treatment response to ivacaftor, which aims to improve CFTR-channel opening probability in patients with gating mutations, shows substantial variability in response, 28% of which can be explained by rs7512462 in SLC26A9 (P = 0.0006). When homozygous Phe508del primary bronchial cells are treated to restore surface CFTR, SLC26A9 likewise modifies treatment response (P = 0.02). Our findings indicate that SLC26A9 airway modification requires CFTR at the cell surface, and that a common variant in SLC26A9 may predict response to CFTR-directed therapeutics.

Prevalence of meconium ileus marks the severity of mutations of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene.

RATIONALE: Meconium ileus (MI) is a perinatal complication in cystic fibrosis (CF), which is only minimally influenced by environmental factors. We derived and examined MI prevalence (MIP) scores to assess CFTR phenotype-phenotype correlation for severe mutations. METHOD: MIP scores were established using a Canadian CF population (n = 2,492) as estimates of the proportion of patients with MI among all patients carrying the same CFTR mutation, focusing on patients with p.F508del as the second allele. Comparisons were made to the registries from the US CF Foundation (n = 43,432), Italy (Veneto/Trentino/Alto Adige regions) (n = 1,788), and Germany (n = 3,596). RESULTS: The prevalence of MI varied among the different registries (13-21%). MI was predominantly prevalent in patients with pancreatic insufficiency carrying "severe" CFTR mutations. In this severe spectrum MIP scores further distinguished between mutation types, for example, G542X (0.31) with a high, F508del (0.22) with a moderate, and G551D (0.08) with a low MIP score. Higher MIP scores were associated with more severe clinical phenotypes, such as a lower forced expiratory volume in 1 second (P = 0.01) and body mass index z score (P = 0.04). CONCLUSIONS: MIP scores can be used to rank CFTR mutations according to their clinical severity and provide a means to expand delineation of CF phenotypes.Genet Med 18 4, 333-340.

Genome-wide association meta-analysis identifies five modifier loci of lung disease severity in cystic fibrosis.

The identification of small molecules that target specific CFTR variants has ushered in a new era of treatment for cystic fibrosis (CF), yet optimal, individualized treatment of CF will require identification and targeting of disease modifiers. Here we use genome-wide association analysis to identify genetic modifiers of CF lung disease, the primary cause of mortality. Meta-analysis of 6,365 CF patients identifies five loci that display significant association with variation in lung disease. Regions on chr3q29 (MUC4/MUC20; P=3.3 × 10(-11)), chr5p15.3 (SLC9A3; P=6.8 × 10(-12)), chr6p21.3 (HLA Class II; P=1.2 × 10(-8)) and chrXq22-q23 (AGTR2/SLC6A14; P=1.8 × 10(-9)) contain genes of high biological relevance to CF pathophysiology. The fifth locus, on chr11p12-p13 (EHF/APIP; P=1.9 × 10(-10)), was previously shown to be associated with lung disease. These results provide new insights into potential targets for modulating lung disease severity in CF.

A Joint Location-Scale Test Improves Power to Detect Associated SNPs, Gene Sets, and Pathways.

Gene-based, pathway, and other multivariate association methods are motivated by the possibility of GxG and GxE interactions; however, accounting for such interactions is limited by the challenges associated with adequate modeling information. Here we propose an easy-to-implement joint location-scale (JLS) association testing framework for single-variant and multivariate analysis that accounts for interactions without explicitly modeling them. We apply the JLS method to a gene-set analysis of cystic fibrosis (CF) lung disease, which is influenced by multiple environmental and genetic factors. We identify and replicate an association between the constituents of the apical plasma membrane and CF lung disease (p = 0.0099 and p = 0.0180, respectively) and highlight a role for the SLC9A3-SLC9A3R1/2-EZR complex in contributing to CF lung disease. Many association studies could benefit from re-analysis with the JLS method that leverages complex genetic architecture for SNP, gene, and pathway identification. Analytical verification, simulation, and additional proof-of-principle applications support our approach.

Deficiency of the ribosome biogenesis gene Sbds in hematopoietic stem and progenitor cells causes neutropenia in mice by attenuating lineage progression in myelocytes.

Shwachman-Diamond syndrome is a congenital bone marrow failure disorder characterized by debilitating neutropenia. The disease is associated with loss-of-function mutations in the SBDS gene, implicated in ribosome biogenesis, but the cellular and molecular events driving cell specific phenotypes in ribosomopathies remain poorly defined. Here, we established what is to our knowledge the first mammalian model of neutropenia in Shwachman-Diamond syndrome through targeted downregulation of Sbds in hematopoietic stem and progenitor cells expressing the myeloid transcription factor CCAAT/enhancer binding protein α (Cebpa). Sbds deficiency in the myeloid lineage specifically affected myelocytes and their downstream progeny while, unexpectedly, it was well tolerated by rapidly cycling hematopoietic progenitor cells. Molecular insights provided by massive parallel sequencing supported cellular observations of impaired cell cycle exit and formation of secondary granules associated with the defect of myeloid lineage progression in myelocytes. Mechanistically, Sbds deficiency activated the p53 tumor suppressor pathway and induced apoptosis in these cells. Collectively, the data reveal a previously unanticipated, selective dependency of myelocytes and downstream progeny, but not rapidly cycling progenitors, on this ubiquitous ribosome biogenesis protein, thus providing a cellular basis for the understanding of myeloid lineage biased defects in Shwachman-Diamond syndrome.

In Vivo Senescence in the Sbds-Deficient Murine Pancreas: Cell-Type Specific Consequences of Translation Insufficiency.

Genetic models of ribosome dysfunction show selective organ failure, highlighting a gap in our understanding of cell-type specific responses to translation insufficiency. Translation defects underlie a growing list of inherited and acquired cancer-predisposition syndromes referred to as ribosomopathies. We sought to identify molecular mechanisms underlying organ failure in a recessive ribosomopathy, with particular emphasis on the pancreas, an organ with a high and reiterative requirement for protein synthesis. Biallelic loss of function mutations in SBDS are associated with the ribosomopathy Shwachman-Diamond syndrome, which is typified by pancreatic dysfunction, bone marrow failure, skeletal abnormalities and neurological phenotypes. Targeted disruption of Sbds in the murine pancreas resulted in p53 stabilization early in the postnatal period, specifically in acinar cells. Decreased Myc expression was observed and atrophy of the adult SDS pancreas could be explained by the senescence of acinar cells, characterized by induction of Tgfß, p15(Ink4b) and components of the senescence-associated secretory program. This is the first report of senescence, a tumour suppression mechanism, in association with SDS or in response to a ribosomopathy. Genetic ablation of p53 largely resolved digestive enzyme synthesis and acinar compartment hypoplasia, but resulted in decreased cell size, a hallmark of decreased translation capacity. Moreover, p53 ablation resulted in expression of acinar dedifferentiation markers and extensive apoptosis. Our findings indicate a protective role for p53 and senescence in response to Sbds ablation in the pancreas. In contrast to the pancreas, the Tgfß molecular signature was not detected in fetal bone marrow, liver or brain of mouse models with constitutive Sbds ablation. Nevertheless, as observed with the adult pancreas phenotype, disease phenotypes of embryonic tissues, including marked neuronal cell death due to apoptosis, were determined to be p53-dependent. Our findings therefore point to cell/tissue-specific responses to p53-activation that include distinction between apoptosis and senescence pathways, in the context of translation disruption.

Variants in Solute Carrier SLC26A9 Modify Prenatal Exocrine Pancreatic Damage in Cystic Fibrosis.

OBJECTIVES: To test the hypothesis that multiple constituents of the apical plasma membrane residing alongside the causal cystic fibrosis (CF) transmembrane conductance regulator protein, including known CF modifiers SLC26A9, SLC6A14, and SLC9A3, would be associated with prenatal exocrine pancreatic damage as measured by newborn screened (NBS) immunoreactive trypsinogen (IRT) levels. STUDY DESIGN: NBS IRT measures and genome-wide genotype data were available on 111 subjects from Colorado, 37 subjects from Wisconsin, and 80 subjects from France. Multiple linear regression was used to determine whether any of 8 single nucleotide polymorphisms (SNPs) in SLC26A9, SLC6A14, and SLC9A3 were associated with IRT and whether other constituents of the apical plasma membrane contributed to IRT. RESULTS: In the Colorado sample, 3 SLC26A9 SNPs were associated with NBS IRT (min P=1.16×10(-3); rs7512462), but no SLC6A14 or SLC9A3 SNPs were associated (P>.05). The rs7512462 association replicated in the Wisconsin sample (P=.03) but not in the French sample (P=.76). Furthermore, rs7512462 was the top-ranked apical membrane constituent in the combined Colorado and Wisconsin sample. CONCLUSIONS: NBS IRT is a biomarker of prenatal exocrine pancreatic disease in patients with CF, and a SNP in SLC26A9 accounts for significant IRT variability. This work suggests SLC26A9 as a potential therapeutic target to ameliorate exocrine pancreatic disease.