Multiple apical plasma membrane constituents are associated with susceptibility to meconium ileus in individuals with cystic fibrosis

Disorder within order: Identification of the disordered loop of STAS domain as the inhibitory domain in SLC26A9 chloride channel

The chloride transporter-channel SLC26A9 is mediated by a reciprocal regulatory mechanism through the interaction between its cytoplasmic sulfate transporter and anti-sigma (STAS) domain and the R domain of cystic fibrosis (CF) transmembrane regulator. In vertebrate Slc26a9s, the STAS domain structures are interrupted by a disordered loop which is conserved in mammals but is variable in nonmammals. Despite the numerous studies involving the STAS domains in SLC26 proteins, the role of the disordered loop region has not been identified. Deletion of the entire Slc26a9-STAS domain results in loss of Cl- channel function. Surprisingly, we found that partial or full deletion of the STAS-disordered loop substantially increases the SLC26A9 chloride transport-channel activity. Bioinformatics analysis reveals that the disordered loops there are three subregions: a K/R-rich region, a "middle" region, and an ordered S/T-rich motif. In this study, the role of this STAS-disordered loop is investigated by using serial deletions and the ordered S/T-rich motif is examined by serial alanine substitution. Substitutions of alanine for serine or threonine in the 620 to 628 S/T-rich motif decrease SLC26A9 chloride channel activity. These experiments parse the functional roles of SLC26A9-STAS-disordered loop and its subdivisions modifying overall SLC26A9 activities. Recently, SLC26A9 has emerged as one of the potential substitutes for abnormal CF transmembrane regulator in CF. Our findings suggest that deletion of variable loop of human SLC26A9 may provide a new gene therapy strategy in the treatment of CF disease.

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

The solute carrier family 26 member 9 modifies rapidly progressing cystic fibrosis associated with homozygous F508del CFTR mutation

BACKGROUND AND AIMS

Cystic fibrosis (CF) is an autosomal recessive disease caused by mutations to the CF transmembrane conductance regulator (CFTR). Symptoms and severity of the disease can be quite variable suggesting modifier genes play an important role.

MATERIALS AND METHODS

Exome sequencing was performed on six individuals carrying homozygous deltaF508 for CFTR genotype but present with rapidly progressing CF (RPCF). Data was analyzed using an unbiased genome-wide genetic burden test against 3076 controls. Single cell RNA sequencing data from LungMAP was utilized to evaluate unique and co-expression of candidate genes, and structural modeling to evaluate the deleterious effects of identified candidate variants.

RESULTS

We have identified solute carrier family 26 member 9 (SLC26A9) as a modifier gene to be associated with RPCF. Two rare missense SLC26A9 variants were discovered in three of six individuals deemed to have RPCF: c.229G > A; p.G77S (present in two patients), and c.1885C > T; p.P629S. Co-expression of SLC26A9 and CFTR mRNA is limited across different lung cell types, with the highest level of co-expression seen in human (6.3 %) and mouse (9.0 %) alveolar type 2 (AT2) cells. Structural modeling suggests deleterious effects of these mutations as they are in critical protein domains which might affect the anion transport capability of SLC26A9.

CONCLUSION

The enrichment of rare and potentially deleterious SLC26A9 mutations in patients with RPCF suggests SLC26A9 may act as an alternative anion transporter in CF and is a modifier gene associated with this lung phenotype.

Diagnosis and Management of Simple and Complicated Meconium Ileus in Cystic Fibrosis, a Systematic Review

The early management of neonates with meconium ileus (MI) and cystic fibrosis (CF) is highly variable across countries and is not standardized. We conducted a systematic review according to the Preferred Reporting Items for Systematic Reviews and Meta-analyses statement. The protocol was registered in PROSPERO (CRD42024522838). Studies from three providers of academic search engines were checked for inclusion criteria, using the following search terms: meconium ileus AND cystic fibrosis OR mucoviscidosis. Regarding the patient population studied, the inclusion criteria were defined using our predefined PICOT framework: studies on neonates with simple or complicated meconium which were confirmed to have cystic fibrosis and were conservatively managed or surgically treated. Results: A total of 566 publications from the last 10 years were verified by the authors of this review to find the most recent and relevant data, and only 8 met the inclusion criteria. Prenatally diagnosed meconium pseudocysts, bowel dilation, and ascites on ultrasound are predictors of neonatal surgery and risk factor for negative 12-month clinical outcomes in MI-CF newborns. For simple MI, conservative treatment with hypertonic solutions enemas can be effective in more than 25% of cases. If repeated enemas fail to disimpact the bowels, the Bishop-Koop stoma is a safe option. No comprehensive research has been conducted so far to determine the ideal surgical protocol for complicated MI. We only found three studies that reported the types of stomas performed and another study comparing the outcomes of patients depending on the surgical management; the conclusions are contradictory especially since the number of cases analyzed in each study was small. Between 18% and 38% of patients with complicated MI will require reoperation for various complications and the mortality rate varies between 0% and 8%. Conclusion: This study reveals a lack of strong data to support management decisions, unequivocally shows that the care of infants with MI is not standardized, and suggests a great need for international collaborative studies.

Pancreatic enzymes digest obstructive meconium from cystic fibrosis pig intestines

Introduction

Meconium ileus (MI) is a life-threatening obstruction of the intestines affecting ∼15% of newborns with cystic fibrosis (CF). Current medical treatments for MI often fail, requiring surgical intervention. MI typically occurs in newborns with pancreatic insufficiency from CF. Meconium contains mucin glycoprotein, a potential substrate for pancreatic enzymes or mucolytics. Our study aim was to determine whether pancreatic enzymes in combination with mucolytic treatments dissolve obstructive meconium using the CF pig model.

Methods

We collected meconium from CF pigs at birth and submerged it in solutions with and without pancreatic enzymes, including normal saline, 7% hypertonic saline, and the reducing agents N-acetylcysteine (NAC) and dithiothreitol (DTT). We digested meconium at 37 °C with agitation, and measured meconium pigment release by spectrophotometry and residual meconium solids by filtration.

Results and discussion

In CF pigs, meconium appeared as a solid pigmented mass obstructing the ileum. Meconium microscopically contained mucus glycoprotein, cellular debris, and bile pigments. Meconium fragments released pigments with maximal absorption at 405 nm after submersion in saline over approximately 8 h. Pancreatic enzymes significantly increased pigment release and decreased residual meconium solids. DTT did not improve meconium digestion and the acidic reducing agent NAC worsened digestion. Pancreatic enzymes digested CF meconium best at neutral pH in isotonic saline. We conclude that pancreatic enzymes digest obstructive meconium from CF pigs, while hydrating or reducing agents alone were less effective. This work suggests a potential role for pancreatic enzymes in relieving obstruction due to MI in newborns with CF.

Convergent gene losses and pseudogenizations in multiple lineages of stomachless fishes

The regressive evolution of independent lineages often results in convergent phenotypes. Several teleost groups display secondary loss of the stomach, and four gastric genes, atp4a, atp4b, pgc, and pga2 have been co-deleted in agastric (stomachless) fish. Analyses of genotypic convergence among agastric fishes showed that four genes, slc26a9, kcne2, cldn18a, and vsig1, were co-deleted or pseudogenized in most agastric fishes of the four major groups. kcne2 and vsig1 were also deleted or pseudogenized in the agastric monotreme echidna and platypus, respectively. In the stomachs of sticklebacks, these genes are expressed in gastric gland cells or surface epithelial cells. An ohnolog of cldn18 was retained in some agastric teleosts but exhibited an increased non-synonymous substitution when compared with gastric species. These results revealed novel convergent gene losses at multiple loci among the four major groups of agastric fish, as well as a single gene loss in the echidna and platypus.

SLC26 Anion Transporters

Solute carrier family 26 (SLC26) is a family of functionally diverse anion transporters found in all kingdoms of life. Anions transported by SLC26 proteins include chloride, bicarbonate, and sulfate, but also small organic dicarboxylates such as fumarate and oxalate. The human genome encodes ten functional homologs, several of which are causally associated with severe human diseases, highlighting their physiological importance. Here, we review novel insights into the structure and function of SLC26 proteins and summarize the physiological relevance of human members.

SLC26A9 in airways and intestine: secretion or absorption?

SLC26A9 is one out of 11 proteins that belong to the SLC26A family of anion transporters. Apart from expression in the gastrointestinal tract, SLC26A9 is also found in the respiratory system, in male tissues and in the skin. SLC26A9 has gained attention because of its modifier role in the gastrointestinal manifestation of cystic fibrosis (CF). SLC26A9 appears to have an impact on the extent of intestinal obstruction caused by meconium ileus. SLC26A9 supports duodenal bicarbonate secretion, but was assumed to provide a basal Cl- secretory pathway in airways. However, recent results show that basal airway Cl- secretion is due to cystic fibrosis conductance regulator (CFTR), while SLC26A9 may rather secrete HCO₃-, thereby maintaining proper airway surface liquid (ASL) pH. Moreover, SLC26A9 does not secrete but probably supports reabsorption of fluid particularly in the alveolar space, which explains early death by neonatal distress in Slc26a9-knockout animals. While the novel SLC26A9 inhibitor S9-A13 helped to unmask the role of SLC26A9 in the airways, it also provided evidence for an additional role in acid secretion by gastric parietal cells. Here we discuss recent data on the function of SLC26A9 in airways and gut, and how S9-A13 may be useful in unraveling the physiological role of SLC26A9.

Transgenic ferret models define pulmonary ionocyte diversity and function

Speciation leads to adaptive changes in organ cellular physiology and creates challenges for studying rare cell-type functions that diverge between humans and mice. Rare cystic fibrosis transmembrane conductance regulator (CFTR)-rich pulmonary ionocytes exist throughout the cartilaginous airways of humans1,2, but limited presence and divergent biology in the proximal trachea of mice has prevented the use of traditional transgenic models to elucidate ionocyte functions in the airway. Here we describe the creation and use of conditional genetic ferret models to dissect pulmonary ionocyte biology and function by enabling ionocyte lineage tracing (FOXI1-CreERT2::ROSA-TG), ionocyte ablation (FOXI1-KO) and ionocyte-specific deletion of CFTR (FOXI1-CreERT2::CFTRL/L). By comparing these models with cystic fibrosis ferrets3,4, we demonstrate that ionocytes control airway surface liquid absorption, secretion, pH and mucus viscosity-leading to reduced airway surface liquid volume and impaired mucociliary clearance in cystic fibrosis, FOXI1-KO and FOXI1-CreERT2::CFTRL/L ferrets. These processes are regulated by CFTR-dependent ionocyte transport of Cl- and HCO3-. Single-cell transcriptomics and in vivo lineage tracing revealed three subtypes of pulmonary ionocytes and a FOXI1-lineage common rare cell progenitor for ionocytes, tuft cells and neuroendocrine cells during airway development. Thus, rare pulmonary ionocytes perform critical CFTR-dependent functions in the proximal airway that are hallmark features of cystic fibrosis airway disease. These studies provide a road map for using conditional genetics in the first non-rodent mammal to address gene function, cell biology and disease processes that have greater evolutionary conservation between humans and ferrets.

VUS: Variant of uncertain significance or very unclear situation?

The advancements in population screening, including newborn screening, enables the identification of disease-causing variants and timely initiation of treatment. However, screening may also identify mild variants, non-disease variants, and variants of uncertain significance (VUS). The identification of a VUS poses a challenge in terms of diagnostic uncertainty and confusion. X-linked adrenoleukodystrophy (ALD) serves as an illustrative example of this complex issue. ALD is a monogenic neurometabolic disease with a complex clinical presentation and a lack of predictive tests for clinical severity. Despite the success of ALD newborn screening, a significant proportion (62%) of missense variants identified through newborn screening exhibit uncertainty regarding their pathogenicity. Resolving this issue requires ongoing efforts to accurately classify variants and refine screening protocols. While it is undisputable that ALD newborn screening greatly benefits boys with the disease, the identification of VUS underscores the need for continuous research and collaboration in improving screening practices.

Pathogenic Relationships in Cystic Fibrosis and Renal Diseases: CFTR, SLC26A9 and Anoctamins

The Cl--transporting proteins CFTR, SLC26A9, and anoctamin (ANO1; ANO6) appear to have more in common than initially suspected, as they all participate in the pathogenic process and clinical outcomes of airway and renal diseases. In the present review, we will therefore concentrate on recent findings concerning electrolyte transport in the airways and kidneys, and the role of CFTR, SLC26A9, and the anoctamins ANO1 and ANO6. Special emphasis will be placed on cystic fibrosis and asthma, as well as renal alkalosis and polycystic kidney disease. In essence, we will summarize recent evidence indicating that CFTR is the only relevant secretory Cl- channel in airways under basal (nonstimulated) conditions and after stimulation by secretagogues. Information is provided on the expressions of ANO1 and ANO6, which are important for the correct expression and function of CFTR. In addition, there is evidence that the Cl- transporter SLC26A9 expressed in the airways may have a reabsorptive rather than a Cl--secretory function. In the renal collecting ducts, bicarbonate secretion occurs through a synergistic action of CFTR and the Cl-/HCO3- transporter SLC26A4 (pendrin), which is probably supported by ANO1. Finally, in autosomal dominant polycystic kidney disease (ADPKD), the secretory function of CFTR in renal cyst formation may have been overestimated, whereas ANO1 and ANO6 have now been shown to be crucial in ADPKD and therefore represent new pharmacological targets for the treatment of polycystic kidney disease.

Tissue-Specific Regulation of CFTR Gene Expression

More than 2000 variations are described within the CFTR (Cystic Fibrosis Transmembrane Regulator) gene and related to large clinical issues from cystic fibrosis to mono-organ diseases. Although these CFTR-associated diseases have been well documented, a large phenotype spectrum is observed and correlations between phenotypes and genotypes are still not well established. To address this issue, we present several regulatory elements that can modulate CFTR gene expression in a tissue-specific manner. Among them, cis-regulatory elements act through chromatin loopings and take part in three-dimensional structured organization. With tissue-specific transcription factors, they form chromatin modules and can regulate gene expression. Alterations of specific regulations can impact and modulate disease expressions. Understanding all those mechanisms highlights the need to expand research outside the gene to enhance our knowledge.

Validating organoid-derived human intestinal monolayers for personalized therapy in cystic fibrosis

Highly effective drugs modulating the defective protein encoded by the CFTR gene have revolutionized cystic fibrosis (CF) therapy. Preclinical drug-testing on human nasal epithelial (HNE) cell cultures and 3-dimensional human intestinal organoids (3D HIO) are used to address patient-specific variation in drug response and to optimize individual treatment for people with CF. This study is the first to report comparable CFTR functional responses to CFTR modulator treatment among patients with different classes of CFTR gene variants using the three methods of 2D HIO, 3D HIO, and HNE. Furthermore, 2D HIO showed good correlation to clinical outcome markers. A larger measurable CFTR functional range and access to the apical membrane were identified as advantages of 2D HIO over HNE and 3D HIO, respectively. Our study thus expands the utility of 2D intestinal monolayers as a preclinical drug testing tool for CF.

Toward a Systematic Assessment of Sex Differences in Cystic Fibrosis

(1) Background: Cystic fibrosis (CF) is a disease with well-documented clinical differences between female and male patients. However, this gender gap is very poorly studied at the molecular level. (2) Methods: Expression differences in whole blood transcriptomics between female and male CF patients are analyzed in order to determine the pathways related to sex-biased genes and assess their potential influence on sex-specific effects in CF patients. (3) Results: We identify sex-biased genes in female and male CF patients and provide explanations for some sex-specific differences at the molecular level. (4) Conclusion: Genes in key pathways associated with CF are differentially expressed between sexes, and thus may account for the gender gap in morbidity and mortality in CF.

High-quality read-based phasing of cystic fibrosis cohort informs genetic understanding of disease modification

Phasing of heterozygous alleles is critical for interpretation of cis-effects of disease-relevant variation. We sequenced 477 individuals with cystic fibrosis (CF) using linked-read sequencing, which display an average phase block N50 of 4.39 Mb. We use these samples to construct a graph representation of CFTR haplotypes, demonstrating its utility for understanding complex CF alleles. These are visualized in a Web app, CFTbaRcodes, that enables interactive exploration of CFTR haplotypes present in this cohort. We perform fine-mapping and phasing of the chr7q35 trypsinogen locus associated with CF meconium ileus, an intestinal obstruction at birth associated with more severe CF outcomes and pancreatic disease. A 20-kb deletion polymorphism and a PRSS2 missense variant p.Thr8Ile (rs62473563) are shown to independently contribute to meconium ileus risk (p = 0.0028, p = 0.011, respectively) and are PRSS2 pancreas eQTLs (p = 9.5 × 10⁻⁷ and p = 1.4 × 10⁻⁴, respectively), suggesting the mechanism by which these polymorphisms contribute to CF. The phase information from linked reads provides a putative causal explanation for variation at a CF-relevant locus, which also has implications for the genetic basis of non-CF pancreatitis, to which this locus has been reported to contribute.

Endocrine Complications of Cystic Fibrosis

Endocrine comorbidities have become increasingly important medical considerations as improving cystic fibrosis (CF) care increases life expectancy. Although the underlying pathophysiology of CF-related diabetes remains elusive, the use of novel technologies and therapeutics seeks to improve both CF-related outcomes and quality of life. Improvements in the overall health of those with CF have tempered concerns about pubertal delay and short stature; however, other comorbidities such as hypogonadism and bone disease are increasingly recognized. Following the introduction of highly effective modulator therapies there are many lessons to be learned about their long-term impact on endocrine comorbidities.

Hypothesis-driven genome-wide association studies provide novel insights into genetics of reading disabilities

Reading Disability (RD) is often characterized by difficulties in the phonology of the language. While the molecular mechanisms underlying it are largely undetermined, loci are being revealed by genome-wide association studies (GWAS). In a previous GWAS for word reading (Price, 2020), we observed that top single-nucleotide polymorphisms (SNPs) were located near to or in genes involved in neuronal migration/axon guidance (NM/AG) or loci implicated in autism spectrum disorder (ASD). A prominent theory of RD etiology posits that it involves disturbed neuronal migration, while potential links between RD-ASD have not been extensively investigated. To improve power to identify associated loci, we up-weighted variants involved in NM/AG or ASD, separately, and performed a new Hypothesis-Driven (HD)-GWAS. The approach was applied to a Toronto RD sample and a meta-analysis of the GenLang Consortium. For the Toronto sample (n = 624), no SNPs reached significance; however, by gene-set analysis, the joint contribution of ASD-related genes passed the threshold (p~1.45 × 10-2, threshold = 2.5 × 10-2). For the GenLang Cohort (n = 26,558), SNPs in DOCK7 and CDH4 showed significant association for the NM/AG hypothesis (sFDR q = 1.02 × 10-2). To make the GenLang dataset more similar to Toronto, we repeated the analysis restricting to samples selected for reading/language deficits (n = 4152). In this GenLang selected subset, we found significant association for a locus intergenic between BTG3-C21orf91 for both hypotheses (sFDR q < 9.00 × 10-4). This study contributes candidate loci to the genetics of word reading. Data also suggest that, although different variants may be involved, alleles implicated in ASD risk may be found in the same genes as those implicated in word reading. This finding is limited to the Toronto sample suggesting that ascertainment influences genetic associations.

Modifier Factors of Cystic Fibrosis Phenotypes: A Focus on Modifier Genes

Although cystic fibrosis (CF) is recognized as a monogenic disease, due to variants within the CFTR (Cystic Fibrosis Transmembrane Regulator) gene, an extreme clinical heterogeneity is described among people with CF (pwCF). Apart from the exocrine pancreatic status, most studies agree that there is little association between CFTR variants and disease phenotypes. Environmental factors have been shown to contribute to this heterogeneity, accounting for almost 50% of the variability of the lung function of pwCF. Nevertheless, pwCF with similar CFTR variants and sharing the same environment (such as in siblings) may have highly variable clinical manifestations not explained by CFTR variants, and only partly explained by environmental factors. It is recognized that genetic variants located outside the CFTR locus, named "modifier genes", influence the clinical expression of the disease. This short review discusses the latest studies that have described modifier factors associated with the various CF phenotypes as well as the response to the recent CFTR modulator therapies.

The SLC26A9 inhibitor S9-A13 provides no evidence for a role of SLC26A9 in airway chloride secretion but suggests a contribution to regulation of ASL pH and gastric proton secretion

The solute carrier 26 family member A9 (SLC26A9) is an epithelial anion transporter that is assumed to contribute to airway chloride secretion and surface hydration. Whether SLC26A9 or CFTR is responsible for airway Cl^- transport under basal conditions is still unclear, due to the lack of a specific inhibitor for SLC26A9. In the present study, we report a novel potent and specific inhibitor for SLC26A9, identified by screening of a drug-like molecule library and subsequent chemical modifications. The most potent compound S9-A13 inhibited SLC26A9 with an IC₅₀ of 90.9 ± 13.4 nM. S9-A13 did not inhibit other members of the SLC26 family and had no effects on Cl^- channels such as CFTR, TMEM16A, or VRAC. S9-A13 inhibited SLC26A9 Cl^- currents in cells that lack expression of CFTR. It also inhibited proton secretion by HGT-1 human gastric cells. In contrast, S9-A13 had minimal effects on ion transport in human airway epithelia and mouse trachea, despite clear expression of SLC26A9 in the apical membrane of ciliated cells. In both tissues, basal and stimulated Cl^- secretion was due to CFTR, while acidification of airway surface liquid by S9-A13 suggests a role of SLC26A9 for airway bicarbonate secretion.

Flux coupling, not specificity, shapes the transport and phylogeny of SLC6 glycine transporters

ATB[Formula: see text] (SLC6A14) is a member of the amino acid transporter branch of the SLC6 family along with GlyT1 (SLC6A9) and GlyT2 (SLC6A5), two glycine-specific transporters coupled to 2:1 and 3:1 Na[Formula: see text]:Cl[Formula: see text], respectively. In contrast, ATB[Formula: see text] exhibits broad substrate specificity for all neutral and cationic amino acids, and its ionic coupling remains unsettled. Using the reversal potential slope method, we demonstrate a 3:1:1 Na[Formula: see text]:Cl[Formula: see text]:Gly stoichiometry for ATB[Formula: see text] that is consistent with its 2.1 e/Gly charge coupling. Like GlyT2, ATB[Formula: see text] behaves as a unidirectional transporter with virtually no glycine efflux at negative potentials after uptake, except by heteroexchange as remarkably shown by leucine activation of NMDARs in Xenopus oocytes coexpressing both membrane proteins. Analysis and computational modeling of the charge movement of ATB[Formula: see text] reveal a higher affinity for sodium in the absence of substrate than GlyT2 and a gating mechanism that locks Na[Formula: see text] into the apo-transporter at depolarized potentials. A 3:1 Na[Formula: see text]:Cl[Formula: see text] stoichiometry justifies the concentrative transport properties of ATB[Formula: see text] and explains its trophic role in tumor growth, while rationalizing its phylogenetic proximity to GlyT2 despite their extreme divergence in specificity.

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.

Linear Mixed-Effect Models Through the Lens of Hardy-Weinberg Disequilibrium

For genetic association studies with related individuals, the linear mixed-effect model is the most commonly used method. In this report, we show that contrary to the popular belief, this standard method can be sensitive to departure from Hardy-Weinberg equilibrium (i.e., Hardy-Weinberg disequilibrium) at the causal SNPs in two ways. First, when the trait heritability is treated as a nuisance parameter, although the association test has correct type I error control, the resulting heritability estimate can be biased, often upward, in the presence of Hardy-Weinberg disequilibrium. Second, if the true heritability is used in the linear mixed-effect model, then the corresponding association test can be biased in the presence of Hardy-Weinberg disequilibrium. We provide some analytical insights along with supporting empirical results from simulation and application studies.

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.

SLC6A14 Impacts Cystic Fibrosis Lung Disease Severity via mTOR and Epithelial Repair Modulation

Cystic fibrosis (CF), due to pathogenic variants in CFTR gene, is associated with chronic infection/inflammation responsible for airway epithelium alteration and lung function decline. Modifier genes induce phenotype variability between people with CF (pwCF) carrying the same CFTR variants. Among these, the gene encoding for the amino acid transporter SLC6A14 has been associated with lung disease severity and age of primary airway infection by the bacteria Pseudomonas aeruginosa. In this study, we investigated whether the single nucleotide polymorphism (SNP) rs3788766, located within SLC6A14 promoter, is associated with lung disease severity in a large French cohort of pwCF. We also studied the consequences of this SNP on SLC6A14 promoter activity using a luciferase reporter and the role of SLC6A14 in the mechanistic target of rapamycin kinase (mTOR) signaling pathway and airway epithelial repair. We confirm that SLC6A14 rs3788766 SNP is associated with lung disease severity in pwCF (p = 0.020; n = 3,257, pancreatic insufficient, aged 6-40 years old), with the minor allele G being deleterious. In bronchial epithelial cell lines deficient for CFTR, SLC6A14 promoter activity is reduced in the presence of the rs3788766 G allele. SLC6A14 inhibition with a specific pharmacological blocker reduced 3H-arginine transport, mTOR phosphorylation, and bronchial epithelial repair rates in wound healing assays. To conclude, our study highlights that SLC6A14 genotype might affect lung disease severity of people with cystic fibrosis via mTOR and epithelial repair mechanism modulation in the lung.

Factors Predisposing the Response to Lumacaftor/Ivacaftor in People with Cystic Fibrosis

Lumacaftor/ivacaftor (LUMA-IVA) therapy is prescribed to people with cystic fibrosis (pwCF) homozygous for the Phe508del-CFTR variant to restore CFTR protein function. There is, however, large inter-individual variability in treatment response. Here, we seek to identify clinical and/or genetic factors that may modulate the response to this CFTR modulator therapy. A total of 765 pwCF older than 12 years under LUMA-IVA therapy and with lung function and nutritional measurements available before and after treatment initiation were included. Response to treatment was determined by the change in lung function and nutritional status, from baseline and over the first two years after initiation, and it was assessed by weighted generalized estimating equation models. Gains in lung function and nutritional status were observed after 6 months of treatment (on average 2.11 ± 7.81% for percent predicted FEV₁ and 0.44 ± 0.77 kg/m² for BMI) and sustained over the 2 years. We observed that the more severe patients gained the most in lung function and nutritional status. While females started with a nutritional status more impaired than males, they had a larger response and regained BMI Z-score values similar to men after 2 years of treatment. We observed no association between variants in solute carrier (SLC) genes and the respiratory function response to LUMA-IVA, but the SLC6A14 rs12839137 variant was associated with the nutritional response. Further investigations, including other genomic regions, will be needed to fully explore the inter-individual variability of the response to LUMA-IVA.

A flexible summary statistics-based colocalization method with application to the mucin cystic fibrosis lung disease modifier locus

Mucus obstruction is a central feature in the cystic fibrosis (CF) airways. A genome-wide association study (GWAS) of lung disease by the CF Gene Modifier Consortium (CFGMC) identified a significant locus containing two mucin genes, MUC20 and MUC4. Expression quantitative trait locus (eQTL) analysis using human nasal epithelia (HNE) from 94 CF-affected Canadians in the CFGMC demonstrated MUC4 eQTLs that mirrored the lung association pattern in the region, suggesting that MUC4 expression may mediate CF lung disease. Complications arose, however, with colocalization testing using existing methods: the locus is complex and the associated SNPs span a 0.2 Mb region with high linkage disequilibrium (LD) and evidence of allelic heterogeneity. We previously developed the Simple Sum (SS), a powerful colocalization test in regions with allelic heterogeneity, but SS assumed eQTLs to be present to achieve type I error control. Here we propose a two-stage SS (SS2) colocalization test that avoids a priori eQTL assumptions, accounts for multiple hypothesis testing and the composite null hypothesis, and enables meta-analysis. We compare SS2 to published approaches through simulation and demonstrate type I error control for all settings with the greatest power in the presence of high LD and allelic heterogeneity. Applying SS2 to the MUC20/MUC4 CF lung disease locus with eQTLs from CF HNE revealed significant colocalization with MUC4 (p = 1.31 × 10-5) rather than with MUC20. The SS2 is a powerful method to inform the responsible gene(s) at a locus and guide future functional studies. SS2 has been implemented in the application LocusFocus.

SLC26A9 as a Potential Modifier and Therapeutic Target in Cystic Fibrosis Lung Disease

SLC26A9 belongs to the solute carrier family 26 (SLC26), which comprises membrane proteins involved in ion transport mechanisms. On the basis of different preliminary findings, including the phenotype of SlC26A9-deficient mice and its possible role as a gene modifier of the human phenotype and treatment response, SLC26A9 has emerged as one of the most interesting alternative targets for the treatment of cystic fibrosis (CF). However, despite relevant clues, some open issues and controversies remain. The lack of specific pharmacological modulators, the elusive expression reported in the airways, and its complex relationships with CFTR and the CF phenotype prevent us from conclusively understanding the contribution of SLC26A9 in human lung physiology and its real potential as a therapeutic target in CF. In this review, we summarized the various studies dealing with SLC26A9 expression, molecular structure, and function as an anion channel or transporter; its interaction and functional relationships with CFTR; and its role as a gene modifier and tried to reconcile them in order to highlight the current understanding and the gap in knowledge regarding the contribution of SLC26A9 to human lung physiology and CF disease and treatment.

Hydrokinetic pancreatic function and insulin secretion are moduled by Cl- uniporter Slc26a9 in mice

AIM

Slc26a9 is a member of the Slc26 multifunctional anion transporter family. Polymorphisms in Slc26a9 are associated with an increased incidence of meconium ileus and diabetes in cystic fibrosis patients. We investigated the expression of Slc26a9 in the murine pancreatic ducts, islets and parenchyma, and elucidated its role in pancreatic ductal electrolyte and fluid secretion and endocrine function.

METHODS

Pancreatic Slc26a9 and CFTR mRNA expression, fluid and bicarbonate secretion were assessed in slc26a9^-/- mice and their age- and sex-matched wild-type (wt) littermates. Glucose and insulin tolerance tests were performed.

RESULTS

Compared with stomach, the mRNA expression of Slc26a9 was low in pancreatic parenchyma, 20-fold higher in microdissected pancreatic ducts than parenchyma, and very low in islets. CFTR mRNA was ~10 fold higher than Slc26a9 mRNA expression in each pancreatic cell type. Significantly reduced pancreatic fluid secretory rates and impaired glucose tolerance were observed in female slc26a9^-/- mice, whereas alterations in male mice did not reach statistical significance. No significant difference was observed in peripheral insulin resistance in slc26a9^-/- compared to sex- and aged-matched wt controls. In contrast, isolated slc26a9^-/- islets in short term culture displayed no difference in insulin content, but a significantly reduced glucose-stimulated insulin secretion compared to age- and sex-matched wt islets, suggesting that the impaired glucose tolerance in the absence of Slc26a9 expression these is a pancreatic defect.

CONCLUSIONS

Deletion of Slc26a9 is associated with a reduction in pancreatic fluid secretion and impaired glucose tolerance in female mice. The results underline the importance of Slc26a9 in pancreatic physiology.

SLC26A9 is selected for endoplasmic reticulum associated degradation (ERAD) via Hsp70-dependent targeting of the soluble STAS domain

SLC26A9, a member of the solute carrier protein family, transports chloride ions across various epithelia. SLC26A9 also associates with other ion channels and transporters linked to human health, and in some cases these heterotypic interactions are essential to support the biogenesis of both proteins. Therefore, understanding how this complex membrane protein is initially folded might provide new therapeutic strategies to overcome deficits in the function of SLC26A9 partners, one of which is associated with Cystic Fibrosis. To this end, we developed a novel yeast expression system for SLC26A9. This facile system has been used extensively with other ion channels and transporters to screen for factors that oversee protein folding checkpoints. As commonly observed for other channels and transporters, we first noted that a substantial fraction of SLC26A9 is targeted for endoplasmic reticulum associated degradation (ERAD), which destroys folding-compromised proteins in the early secretory pathway. We next discovered that ERAD selection requires the Hsp70 chaperone, which can play a vital role in ERAD substrate selection. We then created SLC26A9 mutants and found that the transmembrane-rich domain of SLC26A9 was quite stable, whereas the soluble cytosolic STAS domain was responsible for Hsp70-dependent ERAD. To support data obtained in the yeast model, we were able to recapitulate Hsp70-facilitated ERAD of the STAS domain in human tissue culture cells. These results indicate that a critical barrier to nascent membrane protein folding can reside within a specific soluble domain, one that is monitored by components associated with the ERAD machinery.

Genetic Modifying Factors of Cystic Fibrosis Phenotype: A Challenge for Modern Medicine

Cystic fibrosis (CF) is a monogenic autosomal recessive disease caused by cystic fibrosis transmembrane conductance regulator (CFTR) gene mutations. CF is characterized by a high phenotypic variability present even in patients with the same genotype. This is due to the intervention of modifier genes that interact with both the CFTR gene and environmental factors. The purpose of this review is to highlight the role of non-CFTR genetic factors (modifier genes) that contribute to phenotypic variability in CF. We analyzed literature data starting with candidate gene studies and continuing with extensive studies, such as genome-wide association studies (GWAS) and whole exome sequencing (WES). The results of both types of studies revealed that the number of modifier genes in CF patients is impressive. Their identification offers a new perspective on the pathophysiological mechanisms of the disease, paving the way for the understanding of other genetic disorders. In conclusion, in the future, genetic analysis, such as GWAS and WES, should be performed routinely. A challenge for future research is to integrate their results in the process of developing new classes of drugs, with a goal to improve the prognosis, increase life expectancy, and enhance quality of life among CF patients.

Synergy in Cystic Fibrosis Therapies: Targeting SLC26A9

SLC26A9, a constitutively active Cl^- transporter, has gained interest over the past years as a relevant disease modifier in several respiratory disorders including Cystic Fibrosis (CF), asthma, and non-CF bronchiectasis. SLC26A9 contributes to epithelial Cl^- secretion, thus preventing mucus obstruction under inflammatory conditions. Additionally, SLC26A9 was identified as a CF gene modifier, and its polymorphisms were shown to correlate with the response to drugs modulating CFTR, the defective protein in CF. Here, we aimed to investigate the relationship between SLC26A9 and CFTR, and its role in CF pathogenesis. Our data show that SLC26A9 expression contributes to enhanced CFTR expression and function. While knocking-down SLC26A9 in human bronchial cells leads to lower wt- and F508del-CFTR expression, function, and response to CFTR correctors, the opposite occurs upon its overexpression, highlighting SLC26A9 relevance for CF. Accordingly, F508del-CFTR rescue by the most efficient correctors available is further enhanced by increasing SLC26A9 expression. Interestingly, SLC26A9 overexpression does not increase the PM expression of non-F508del CFTR traffic mutants, namely those unresponsive to corrector drugs. Altogether, our data indicate that SLC26A9 stabilizes CFTR at the ER level and that the efficacy of CFTR modulator drugs may be further enhanced by increasing its expression.

Separating the contributions of SLC26A9 and CFTR to anion secretion in primary human bronchial epithelia

Aberrant anion secretion across the bronchial epithelium is associated with airway disease, most notably in cystic fibrosis. Although the cystic fibrosis transmembrane conductance regulator (CFTR) is recognized as the primary source of airway anion secretion, alternative anion transport mechanisms play a contributing role. An alternative anion transporter of growing interest is SLC26A9, a constitutively active chloride channel that has been shown to interact with CFTR and may also contribute to bicarbonate secretion. Interest in SLC26A9 has been fueled by genome-wide association studies that suggest it is a significant modifier of CF disease severity. Despite this growing evidence that SLC26A9 plays an important role in the airway, its presence and function in bronchial epithelia remain poorly understood, in part, because its activity is difficult to separate from the activity of CFTR. Here, we present results using primary human bronchial epithelia (HBE) from multiple patient sources to confirm that SLC26A9 mRNA is present in HBE and that its constitutive channel activity is unaffected by knockdown of CFTR. Furthermore, SLC26A9 and CFTR show differential responses to common inhibitors of anion secretion. Finally, we assess the impact of bicarbonate on the activity of SLC26A9 and CFTR. These results confirm that SLC26A9 is the primary source of constitutive anion secretion across HBE, and should inform future studies focused on activation of SLC26A9 as an alternative anion channel in CF. These results should provide a strong foundation to investigate how single-nucleotide polymorphisms in SLC26A9 modulate airway disease.

Small-molecule drugs for cystic fibrosis: Where are we now?

The cystic fibrosis (CF) lung disease is due to the lack/dysfunction of the CF Transmembrane Conductance Regulator (CFTR), a chloride channel expressed by epithelial cells as the main regulator of ion and fluid homeostasis. More than 2000 genetic variation in the CFTR gene are known, among which those with identified pathomechanism have been divided into six VI mutation classes. A major advancement in the pharmacotherapy of CF has been the development of small-molecule drugs hitting the root of the disease, i.e. the altered ion and fluid transport through the airway epithelium. These drugs, called CFTR modulators, have been advanced to the clinics to treat nearly 90% of CF patients, including the CFTR potentiator ivacaftor, approved for residual function mutations (Classes III and IV), and combinations of correctors (lumacaftor, tezacaftor, elexacaftor) and ivacaftor for patients bearing at least one the F508del mutation, the most frequent mutation belonging to class II. To cover the 10% of CF patients without etiological therapies, other novel small-molecule CFTR modulators are in evaluation of their effectiveness in all the CFTR mutation classes: read-through agents for Class I, correctors, potentiators and amplifiers from different companies for Class II-V, stabilizers for Class VI. In alternative, other solute carriers, such as SLC26A9 and SLC6A14, are the focus of intensive investigation. Finally, other molecular targets are being evaluated for patients with no approved CFTR modulator therapy or as means of enhancing CFTR modulatory therapy, including small molecules forming ion channels, inhibitors of the ENaC sodium channel and potentiators of the calcium-activated chloride channel TMEM16A. This paper aims to give an up-to-date overview of old and novel CFTR modulators as well as of novel strategies based on small-molecule drugs. Further investigations in in-vivo and cell-based models as well as carrying out large prospective studies will be required to determine if novel CFTR modulators, stabilizers, amplifiers, and the ENaC inhibitors or TMEM16A potentiators will further improve the clinical outcomes in CF management.

New generation ENaC inhibitors detach cystic fibrosis airway mucus bundles via sodium/hydrogen exchanger inhibition

Cystic fibrosis (CF) is a recessive inherited disease caused by mutations affecting anion transport by the epithelial ion channel cystic fibrosis transmembrane conductance regulator (CFTR). The disease is characterized by mucus accumulation in the airways and intestine, but the major cause of mortality in CF is airway mucus accumulation, leading to bacterial colonization, inflammation and respiratory failure. Several drug targets are under evaluation to alleviate airway mucus obstruction in CF and one of these targets is the epithelial sodium channel ENaC. To explore effects of ENaC inhibitors on mucus properties, we used two model systems to investigate mucus characteristics, mucus attachment in mouse ileum and mucus bundle transport in piglet airways. We quantified mucus attachment in explants from CFTR null (CF) mice and tracheobronchial explants from newborn CFTR null (CF) piglets to evaluate effects of ENaC or sodium/hydrogen exchanger (NHE) inhibitors on mucus attachment. ENaC inhibitors detached mucus in the CF mouse ileum, although the ileum lacks ENaC expression. This effect was mimicked by two NHE inhibitors. Airway mucus bundles were immobile in untreated newborn CF piglets but were detached by the therapeutic drug candidate AZD5634 (patent WO, 2015140527). These results suggest that the ENaC inhibitor AZD5634 causes detachment of CF mucus in the ileum and airway via NHE inhibition and that drug design should focus on NHE instead of ENaC inhibition.

Pharmacological Modulation of Ion Channels for the Treatment of Cystic Fibrosis

Cystic fibrosis (CF) is a life-shortening monogenic disease caused by mutations in the gene encoding the CF transmembrane conductance regulator (CFTR) protein, an anion channel that transports chloride and bicarbonate across epithelia. Despite clinical progress in delaying disease progression with symptomatic therapies, these individuals still develop various chronic complications in lungs and other organs, which significantly restricts their life expectancy and quality of life. The development of high-throughput assays to screen drug-like compound libraries have enabled the discovery of highly effective CFTR modulator therapies. These novel therapies target the primary defect underlying CF and are now approved for clinical use for individuals with specific CF genotypes. However, the clinically approved modulators only partially reverse CFTR dysfunction and there is still a considerable number of individuals with CF carrying rare CFTR mutations who remain without any effective CFTR modulator therapy. Accordingly, additional efforts have been pursued to identify novel and more potent CFTR modulators that may benefit a larger CF population. The use of ex vivo individual-derived specimens has also become a powerful tool to evaluate novel drugs and predict their effectiveness in a personalized medicine approach. In addition to CFTR modulators, pro-drugs aiming at modulating alternative ion channels/transporters are under development to compensate for the lack of CFTR function. These therapies may restore normal mucociliary clearance through a mutation-agnostic approach (ie, independent of CFTR mutation) and include inhibitors of the epithelial sodium channel (ENaC), modulators of the calcium-activated channel transmembrane 16A (TMEM16, or anoctamin 1) or of the solute carrier family 26A member 9 (SLC26A9), and anionophores. The present review focuses on recent progress and challenges for the development of ion channel/transporter-modulating drugs for the treatment of CF.

The X factor: A robust and powerful approach to X-chromosome-inclusive whole-genome association studies

The X‐chromosome is often excluded from genome‐wide association studies because of analytical challenges. Some of the problems, such as the random, skewed, or no X‐inactivation model uncertainty, have been investigated. Other considerations have received little to no attention, such as the value in considering nonadditive and gene–sex interaction effects, and the inferential consequence of choosing different baseline alleles (i.e., the reference vs. the alternative allele). Here we propose a unified and flexible regression‐based association test for X‐chromosomal variants. We provide theoretical justifications for its robustness in the presence of various model uncertainties, as well as for its improved power when compared with the existing approaches under certain scenarios. For completeness, we also revisit the autosomes and show that the proposed framework leads to a more robust approach than the standard method. Finally, we provide supporting evidence by revisiting several published association studies. Supporting Information for this article are available online.

Score tests for scale effects, with application to genomic analysis

Tests for variance or scale effects due to covariates are used in many areas and recently, in genomic and genetic association studies. We study score tests based on location-scale models with arbitrary error distributions that allow incorporation of additional adjustment covariates. Tests based on Gaussian and Laplacian double generalized linear models are examined in some detail. Numerical properties of the tests under Gaussian and other error distributions are examined. Our results show that the use of model-based asymptotic distributions with score tests for scale effects does not control type 1 error well in many settings of practical relevance. We consider simple statistics based on permutation distribution approximations, which correspond to well-known statistics derived by another approach. They are shown to give good type 1 error control under different error distributions and under covariate distribution imbalance. The methods are illustrated through a differential gene expression analysis involving breast cancer tumor samples.

Cystic Fibrosis Lung Disease Modifiers and Their Relevance in the New Era of Precision Medicine

Our understanding of cystic fibrosis (CF) has grown exponentially since the discovery of the cystic fibrosis transmembrane conductance regulator (CFTR) gene in 1989. With evolving genetic and genomic tools, we have come to better understand the role of CFTR genotypes in the pathophysiology of the disease. This, in turn, has paved the way for the development of modulator therapies targeted at mutations in the CFTR, which are arguably one of the greatest advances in the treatment of CF. These modulator therapies, however, do not target all the mutations in CFTR that are seen in patients with CF and, furthermore, a variation in response is seen in patients with the same genotype who are taking modulator therapies. There is growing evidence to support the role of non-CFTR modifiers, both genetic and environmental, in determining the variation seen in CF morbidity and mortality and also in the response to existing therapies. This review focusses on key findings from studies using candidate gene and genome-wide approaches to identify CF modifier genes of lung disease in cystic fibrosis and considers the interaction between modifiers and the response to modulator therapies. As the use of modulator therapies expands and we gain data around outcomes, it will be of great interest to investigate this interaction further. Going forward, it will also be crucial to better understand the relative influence of genomic versus environmental factors. With this understanding, we can truly begin to deliver personalised care by better profiling the likely disease phenotype for each patient and their response to treatment.

A generalized robust allele-based genetic association test

The allele-based association test, comparing allele frequency difference between case and control groups, is locally most powerful. However, application of the classical allelic test is limited in practice, because the method is sensitive to the Hardy-Weinberg equilibrium (HWE) assumption, not applicable to continuous traits, and not easy to account for covariate effect or sample correlation. To develop a generalized robust allelic test, we propose a new allele-based regression model with individual allele as the response variable. We show that the score test statistic derived from this robust and unifying regression framework contains a correction factor that explicitly adjusts for potential departure from HWE and encompasses the classical allelic test as a special case. When the trait of interest is continuous, the corresponding allelic test evaluates a weighted difference between individual-level allele frequency estimate and sample estimate where the weight is proportional to an individual's trait value, and the test remains valid under Y-dependent sampling. Finally, the proposed allele-based method can analyze multiple (continuous or binary) phenotypes simultaneously and multiallelic genetic markers, while accounting for covariate effect, sample correlation, and population heterogeneity. To support our analytical findings, we provide empirical evidence from both simulation and application studies.

SLC26A9 SNP rs7512462 is not associated with lung disease severity or lung function response to ivacaftor in cystic fibrosis patients with G551D-CFTR

BACKGROUND

The CFTR modulator ivacaftor has been variably effective in treating individuals with cystic fibrosis (CF) who harbor CFTR gating variants such as G551D, as well as other classes of CFTR variants when used with other modulators. Because CFTR genotype does not fully explain this variability, defining genetic modifiers of response to modulator therapy is of particular interest to the field of individualized CF drug therapy. Previous studies have proposed that a variant in SLC26A9 (rs7512462) is associated with lung disease severity and with response to treatment with ivacaftor in individuals with CF who carry G551D or gating variants.

METHODS

Given the implications for CF treatment, we re-examined the reported associations in three cohorts; patients enrolled in the Twin and Siblings study at Johns Hopkins University, the CF modifier study at the University of North Carolina at Chapel Hill, and the prospective G551D Observational (GOAL) study. The GOAL study was specifically designed to measure lung function response to ivacaftor.

RESULTS

We find no association between SLC26A9 (rs7512462) genotype and lung disease severity (n = 272) or change in lung function at one-, three-, and six-month intervals following ivacaftor treatment(n = 141) in individuals with CF who carry at least one G551D variant.

CONCLUSIONS

Our inability to replicate this association indicates that rs7512462 genotype should not be used in treatment decisions.

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.

Bayesian model averaging for the X-chromosome inactivation dilemma in genetic association study

X-chromosome is often excluded from the so called "whole-genome" association studies due to the differences it exhibits between males and females. One particular analytical challenge is the unknown status of X-inactivation, where one of the two X-chromosome variants in females may be randomly selected to be silenced. In the absence of biological evidence in favor of one specific model, we consider a Bayesian model averaging framework that offers a principled way to account for the inherent model uncertainty, providing model averaging-based posterior density intervals and Bayes factors. We examine the inferential properties of the proposed methods via extensive simulation studies, and we apply the methods to a genetic association study of an intestinal disease occurring in about 20% of cystic fibrosis patients. Compared with the results previously reported assuming the presence of inactivation, we show that the proposed Bayesian methods provide more feature-rich quantities that are useful in practice.

Variant classifications, databases and genotype-phenotype correlations

Because CFTR gene studies now represent one of the most frequent genetic analyses routinely performed worldwide, the number of rare CFTR variants identified in various clinical situations, regularly increases. To provide appropriate diagnosis and prognosis to CF patients as well as appropriate genetic counseling to families, the clinical impact and the phenotypic spectrum of variants identified by diagnostic techniques need to be characterized. Three complementary locus specific databases, called CFTR1, CFTR2 and CFTR-France were developed to address these issues. Besides, the growing knowledge of the CF pathophysiology and the technical evolution in molecular biology allowed to identify candidate modifier genes, regulatory loci, epigenetic profiles and trans-regulators that could help to refine genotype-phenotype correlations at the individual level. These different factors may contribute to the large phenotypic variability between patients with CF, even when they carry identical CFTR variants, regarding lung function, meconium ileus susceptibility or the risk for developing CFTR-related diabetes and liver disease. Finally, the availability of new therapies that target the CFTR protein for numbers of CF patients led to the identification of 'good' and 'poor' responders, thus raising questions of pharmacogenetics factors that may influence treatment efficiency as a novel feature of the complexity of CF patients' management. © 2020 French Society of Pediatrics. Published by Elsevier Masson SAS. All rights reserved.

Transport of L-Arginine Related Cardiovascular Risk Markers

L-arginine and its derivatives, asymmetric and symmetric dimethylarginine (ADMA and SDMA) and L-homoarginine, have emerged as cardiovascular biomarkers linked to cardiovascular outcomes and various metabolic and functional pathways such as NO-mediated endothelial function. Cellular uptake and efflux of L-arginine and its derivatives are facilitated by transport proteins. In this respect the cationic amino acid transporters CAT1 and CAT2 (SLC7A1 and SLC7A2) and the system y⁺L amino acid transporters (SLC7A6 and SLC7A7) have been most extensively investigated, so far, but the number of transporters shown to mediate the transport of L-arginine and its derivatives is constantly increasing. In the present review we assess the growing body of evidence regarding the function, expression, and clinical relevance of these transporters and their possible relation to cardiovascular diseases.

The GENDULF algorithm: mining transcriptomics to uncover modifier genes for monogenic diseases

Modifier genes are believed to account for the clinical variability observed in many Mendelian disorders, but their identification remains challenging due to the limited availability of genomics data from large patient cohorts. Here, we present GENDULF (GENetic moDULators identiFication), one of the first methods to facilitate prediction of disease modifiers using healthy and diseased tissue gene expression data. GENDULF is designed for monogenic diseases in which the mechanism is loss of function leading to reduced expression of the mutated gene. When applied to cystic fibrosis, GENDULF successfully identifies multiple, previously established disease modifiers, including EHF, SLC6A14, and CLCA1. It is then utilized in spinal muscular atrophy (SMA) and predicts U2AF1 as a modifier whose low expression correlates with higher SMN2 pre-mRNA exon 7 retention. Indeed, knockdown of U2AF1 in SMA patient-derived cells leads to increased full-length SMN2 transcript and SMN protein expression. Taking advantage of the increasing availability of transcriptomic data, GENDULF is a novel addition to existing strategies for prediction of genetic disease modifiers, providing insights into disease pathogenesis and uncovering novel therapeutic targets.

How to determine the mechanism of action of CFTR modulator compounds: A gateway to theranostics

The greatest challenge of 21st century biology is to fully understand mechanisms of disease to drive new approaches and medical innovation. Parallel to this is the huge biomedical endeavour of treating people through personalized medicine. Until now all CFTR modulator drugs that have entered clinical trials have been genotype-dependent. An emerging alternative is personalized/precision medicine in CF, i.e., to determine whether rare CFTR mutations respond to existing (or novel) CFTR modulator drugs by pre-assessing them directly on patient's tissues ex vivo, an approach also now termed theranostics. To administer the right drug to the right person it is essential to understand how drugs work, i.e., to know their mechanism of action (MoA), so as to predict their applicability, not just in certain mutations but also possibly in other diseases that share the same defect/defective pathway. Moreover, an understanding the MoA of a drug before it is tested in clinical trials is the logical path to drug discovery and can increase its chance for success and hence also approval. In conclusion, the most powerful approach to determine the MoA of a compound is to understand the underlying biology. Novel large datasets of intervenients in most biological processes, namely those emerging from the post-genomic era tools, are available and should be used to help in this task.

Amino Acid Transporter SLC6A14 (ATB0,+) - A Target in Combined Anti-cancer Therapy

Cancer cells are characterized by quick growth and proliferation, demanding constant supply of various nutrients. Several plasma membrane transporters delivering such compounds are upregulated in cancer. Solute carrier family 6 member 14 (SLC6A14), known as amino acid transporter B^0,+ (ATB^0,+) transports all amino acids with exception of the acidic ones: aspartate and glutamate. Its malfunctioning is correlated with several pathological states and it is upregulated in solid tumors. The high expression of SLC6A14 is prognostic and unfavorable in pancreatic cancer, while in breast cancer it is expressed in estrogen receptor positive cells. As many plasma membrane transporters it resides in endoplasmic reticulum (ER) membrane after translation before further trafficking through Golgi to the cell surface. Transporter exit from ER is strictly controlled. The proper folding of SLC6A14 was shown to be controlled from the cytoplasmic side by heat shock proteins, further exit from ER and formation of coatomer II (COPII) coated vesicles depends on specific interaction with COPII cargo-recognizing subunit SEC24C, phosphorylated by kinase AKT. Inhibition of heat shock proteins, known to be upregulated in cancer, directs SLC6A14 to degradation. Targeting proteins regulating SLC6A14 trafficking is proposed as an additional pharmacological treatment of cancer.