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.

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.

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.

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.

Recent advances in developing therapeutics for cystic fibrosis.

Despite hope that a cure was imminent when the causative gene was cloned nearly 30 years ago, cystic fibrosis (CF [MIM: 219700]) remains a life-shortening disease affecting more than 70 000 individuals worldwide. However, within the last 6 years the Food and Drug Administration's approval of Ivacaftor, the first drug that corrects the defective cystic fibrosis transmembrane conductance regulator protein [CFTR (MIM: 602421)] in patients with the G551D mutation, marks a watershed in the development of novel therapeutics for this devastating disease. Here we review recent progress in diverse research areas, which all focus on curing CF at the genetic, biochemical or physiological level. In the near future it seems probable that development of mutation-specific therapies will be the focus, since it is unlikely that any one approach will be efficient in correcting the more than 2000 disease-associated variants. We discuss the new drugs and combinations of drugs that either enhance delivery of misfolded CFTR protein to the cell membrane, where it functions as an ion channel, or that activate channel opening. Next we consider approaches to correct the causative genetic lesion at the DNA or RNA level, through repressing stop mutations and nonsense-mediated decay, modulating splice mutations, fixing errors by gene editing or using novel routes to gene replacement. Finally, we explore how modifier genes, loci elsewhere in the genome that modify CF disease severity, may be used to restore a normal phenotype. Progress in all of these areas has been dramatic, generating enthusiasm that CF may soon become a broadly treatable disease.

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.

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.

SLCO5A1 and synaptic assembly genes contribute to impulsivity in juvenile myoclonic epilepsy.

Elevated impulsivity is a key component of attention-deficit hyperactivity disorder (ADHD), bipolar disorder and juvenile myoclonic epilepsy (JME). We performed a genome-wide association, colocalization, polygenic risk score, and pathway analysis of impulsivity in JME (n = 381). Results were followed up with functional characterisation using a drosophila model. We identified genome-wide associated SNPs at 8q13.3 (P = 7.5 × 10-9) and 10p11.21 (P = 3.6 × 10-8). The 8q13.3 locus colocalizes with SLCO5A1 expression quantitative trait loci in cerebral cortex (P = 9.5 × 10-3). SLCO5A1 codes for an organic anion transporter and upregulates synapse assembly/organisation genes. Pathway analysis demonstrates 12.7-fold enrichment for presynaptic membrane assembly genes (P = 0.0005) and 14.3-fold enrichment for presynaptic organisation genes (P = 0.0005) including NLGN1 and PTPRD. RNAi knockdown of Oatp30B, the Drosophila polypeptide with the highest homology to SLCO5A1, causes over-reactive startling behaviour (P = 8.7 × 10-3) and increased seizure-like events (P = 6.8 × 10-7). Polygenic risk score for ADHD genetically correlates with impulsivity scores in JME (P = 1.60 × 10-3). SLCO5A1 loss-of-function represents an impulsivity and seizure mechanism. Synaptic assembly genes may inform the aetiology of impulsivity in health and disease.

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-7 and p = 1.4 × 10-4, 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.

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.

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.

Sex-specific disease modifiers in juvenile myoclonic epilepsy.

Juvenile myoclonic epilepsy (JME) is a common idiopathic generalised epilepsy with variable seizure prognosis and sex differences in disease presentation. Here, we investigate the combined epidemiology of sex, seizure types and precipitants, and their influence on prognosis in JME, through cross-sectional data collected by The Biology of Juvenile Myoclonic Epilepsy (BIOJUME) consortium. 765 individuals met strict inclusion criteria for JME (female:male, 1.8:1). 59% of females and 50% of males reported triggered seizures, and in females only, this was associated with experiencing absence seizures (OR = 2.0, p < 0.001). Absence seizures significantly predicted drug resistance in both males (OR = 3.0, p = 0.001) and females (OR = 3.0, p < 0.001) in univariate analysis. In multivariable analysis in females, catamenial seizures (OR = 14.7, p = 0.001), absence seizures (OR = 6.0, p < 0.001) and stress-precipitated seizures (OR = 5.3, p = 0.02) were associated with drug resistance, while a photoparoxysmal response predicted seizure freedom (OR = 0.47, p = 0.03). Females with both absence seizures and stress-related precipitants constitute the prognostic subgroup in JME with the highest prevalence of drug resistance (49%) compared to females with neither (15%) and males (29%), highlighting the unmet need for effective, targeted interventions for this subgroup. We propose a new prognostic stratification for JME and suggest a role for circuit-based risk of seizure control as an avenue for further investigation.

Trait impulsivity in Juvenile Myoclonic Epilepsy.

OBJECTIVE: Impulsivity is a multidimensional construct that can predispose to psychopathology. Meta-analysis demonstrates an association between response impulsivity and Juvenile Myoclonic Epilepsy (JME), a common genetic generalized epilepsy. Here, we test the hypotheses that trait impulsivity is (i) elevated in JME compared to controls; (ii) moderated by specific seizure characteristics; and (iii) associated with psychiatric adverse effects of antiepileptic drugs (AEDs). METHODS: 322 participants with JME and 126 age and gender-matched controls completed the Barratt's Impulsiveness Scale (BIS-brief) alongside information on seizure history and AED use. We compared group BIS-brief scores and assessed associations of JME BIS-brief scores with seizure characteristics and AED adverse effects. RESULTS: The mean BIS-brief score in JME was 18.1 ± 4.4 compared with 16.2 ± 4.1 in controls (P = 0.0007). Elevated impulsivity was associated with male gender (P = 0.027), frequent absence seizures (P = 0.0004) and lack of morning predominance of myoclonus (P = 0.008). High impulsivity significantly increased the odds of a psychiatric adverse event on levetiracetam (P = 0.036), but not any other psychiatric or somatic adverse effects. INTERPRETATION: Trait impulsivity is elevated in JME and comparable to scores in personality and neurotic disorders. Increased seizure frequency and absence of circadian seizure pattern moderate BIS score, suggesting disruption of both cortico-striatal and thalamocortical networks as a shared mechanism between seizures and impulsivity in JME. These findings warrant consideration of impulsivity as a distinct target of intervention, and as a stratifying factor for AED treatment in JME, and perhaps other types of epilepsy. The role of impulsivity in treatment adherence and psychosocial outcome requires further investigation.

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.

A microRNA-328 binding site in PAX6 is associated with centrotemporal spikes of rolandic epilepsy.

OBJECTIVE: Rolandic epilepsy is a common genetic focal epilepsy of childhood characterized by centrotemporal sharp waves on electroencephalogram. In previous genome-wide analysis, we had reported linkage of centrotemporal sharp waves to chromosome 11p13, and fine mapping with 44 SNPs identified the ELP4-PAX6 locus in two independent US and Canadian case-control samples. Here, we aimed to find a causative variant for centrotemporal sharp waves using a larger sample and higher resolution genotyping array. METHODS: We fine-mapped the ELP4-PAX6 locus in 186 individuals from rolandic epilepsy families and 1000 population controls of European origin using the Illumina HumanCoreExome-12 v1.0 BeadChip. Controls were matched to cases on ethnicity using principal component analysis. We used generalized estimating equations to assess association, followed up with a bioinformatics survey and literature search to evaluate functional significance. RESULTS: Homozygosity at the T allele of SNP rs662702 in the 3' untranslated region of PAX6 conferred increased risk of CTS: Odds ratio = 12.29 (95% CI: 3.20-47.22), P = 2.6 × 10(-4) and is seen in 3.9% of cases but only 0.3% of controls. INTERPRETATION: The minor T allele of SNP rs662702 disrupts regulation by microRNA-328, which is known to result in increased PAX6 expression in vitro. This study provides, for the first time, evidence of a noncoding genomic variant contributing to the etiology of a common human epilepsy via a posttranscriptional regulatory mechanism.

GLP-1 receptor activation indirectly reduces hepatic lipid accumulation but does not attenuate development of atherosclerosis in diabetic male ApoE(-/-) mice.

Glucagon-like peptide-1 receptor (GLP-1R) agonists reduce lipid accumulation in peripheral tissues, attenuating atherosclerosis and hepatic steatosis in preclinical studies. We examined whether GLP-1R activation decreases atherosclerosis progression in high-fat diet-fed male ApoE(-/-) mice after administration of streptozotocin and treatment with the long-acting GLP-1R agonist taspoglutide administered once monthly vs. metformin in the drinking water for 12 wk. Taspoglutide did not reduce plaque area or lipid content in the aortic arch or abdominal aorta, and no significant change in aortic macrophage accumulation was detected after taspoglutide or metformin. In contrast, hepatic triglyceride levels were significantly reduced in livers from taspoglutide-treated mice. Both peripheral and intracerebroventricular administration of exendin-4 rapidly decreased plasma triglyceride levels in fasted mice, and taspoglutide therapy in ApoE(-/-) mice modulated the expression of hepatic genes controlling fatty acid uptake and oxidation. We were unable to detect expression of the entire Glp1r coding sequence in macrophages isolated from ApoE(-/-), C57BL/6, and IL10(-/-) mice. Similarly, Glp1r mRNA transcripts were not detected in RNA from isolated murine hepatocytes. Using Western blotting and tissue extracts from Glp1r(+/+) and Glp1r(-/-) mice, and cells transfected with a tagged murine GLP-1R cDNA, we could not validate the sensitivity and specificity of three different GLP-1R antisera commonly used for the detection of GLP-1R protein. Taken together, these findings illustrate divergent actions of GLP-1R agonists on atherosclerosis progression and accumulation of ectopic lipid in ApoE(-/-) mice and highlight the importance of indirect GLP-1R actions for the control of hepatic lipid accumulation.

Assessment of the effects of pH, formulation and deformulation on the conformation of interferon alpha-2 by NMR.

This article reports the results of our investigation of the effects of pH and various formulations on the conformation of interferon (IFN) alpha-2a and IFN alpha-2b using the NMR fingerprinting assay. Samples of (15)N-IFN alpha-2 were produced and their activity was inferred by comparing their NMR spectra with those recorded for the corresponding European Directorate for the Quality of Medicines (EDQM) reference standards. The proteins were then mixed with appropriate excipients to reproduce formulations used in innovator products of Roferon-A and Intron-A and deformulated via cation-exchange chromatography. The conformation of IFN alpha-2 was monitored by two-dimensional (2D)-NMR spectroscopy at various pHs, after formulation and deformulation procedures. Our results show that the process does not alter the conformation of IFN alpha-2 and that the optimal pH for deformulation is 4.0 +/- 0.5. Variation in pH below 3.0 causes the protein to unfold, whereas above pH 4.5, the three-dimensional (3D) fold is maintained, but the NMR spectra indicate a propensity to oligomerize. This behaviour is reversible upon readjusting the pH to 3.5-4.5. Here, we demonstrate the applicability of NMR to assess the structure of protein therapeutics. The proposed method can assist in validating analytical methods that require deformulation of IFN-based products.