Canadian COVID-19 host genetics cohort replicates known severity associations

The HostSeq initiative recruited 10,059 Canadians infected with SARS-CoV-2 between March 2020 and March 2023, obtained clinical information on their disease experience and whole genome sequenced (WGS) their DNA. We analyzed the WGS data for genetic contributors to severe COVID-19 (considering 3,499 hospitalized cases and 4,975 non-hospitalized after quality control). We investigated the evidence for replication of loci reported by the International Host Genetics Initiative (HGI); analyzed the X chromosome; conducted rare variant gene-based analysis and polygenic risk score testing. Population stratification was adjusted for using meta-analysis across ancestry groups. We replicated two loci identified by the HGI for COVID-19 severity: the LZTFL1/SLC6A20 locus on chromosome 3 and the FOXP4 locus on chromosome 6 (the latter with a variant significant at P < 5E-8). We found novel significant associations with MRAS and WDR89 in gene-based analyses, and constructed a polygenic risk score that explained 1.01% of the variance in severe COVID-19. This study provides independent evidence confirming the robustness of previously identified COVID-19 severity loci by the HGI and identifies novel genes for further investigation.

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.

Gene copy number variation and pediatric mental health/neurodevelopment in a general population

We assessed the relationship of gene copy number variation (CNV) in mental health/neurodevelopmental traits and diagnoses, physical health and cognition in a community sample of 7100 unrelated children and youth of European or East Asian ancestry (Spit for Science). Clinically significant or susceptibility CNVs were present in 3.9% of participants and were associated with elevated scores on a continuous measure of attention-deficit/hyperactivity disorder (ADHD) traits (P = 5.0 × 10-3), longer response inhibition (a cognitive deficit found in several mental health and neurodevelopmental disorders; P = 1.0 × 10-2) and increased prevalence of mental health diagnoses (P = 1.9 × 10-6, odds ratio: 3.09), specifically ADHD, autism spectrum disorder anxiety and learning problems/learning disorder (P's < 0.01). There was an increased burden of rare deletions in gene-sets related to brain function or expression in brain associated with more ADHD traits. With the current mental health crisis, our data established a baseline for delineating genetic contributors in pediatric-onset conditions.

Language and reading impairments are associated with increased prevalence of non-right-handedness

Handedness has been studied for association with language-related disorders because of its link with language hemispheric dominance. No clear pattern has emerged, possibly because of small samples, publication bias, and heterogeneous criteria across studies. Non-right-handedness (NRH) frequency was assessed in N = 2503 cases with reading and/or language impairment and N = 4316 sex-matched controls identified from 10 distinct cohorts (age range 6-19 years old; European ethnicity) using a priori set criteria. A meta-analysis (Ncases  = 1994) showed elevated NRH % in individuals with language/reading impairment compared with controls (OR = 1.21, CI = 1.06-1.39, p = .01). The association between reading/language impairments and NRH could result from shared pathways underlying brain lateralization, handedness, and cognitive functions.

Variation in prognosis and treatment outcome in juvenile myoclonic epilepsy: a Biology of Juvenile Myoclonic Epilepsy Consortium proposal for a practical definition and stratified medicine classifications

Reliable definitions, classifications and prognostic models are the cornerstones of stratified medicine, but none of the current classifications systems in epilepsy address prognostic or outcome issues. Although heterogeneity is widely acknowledged within epilepsy syndromes, the significance of variation in electroclinical features, comorbidities and treatment response, as they relate to diagnostic and prognostic purposes, has not been explored. In this paper, we aim to provide an evidence-based definition of juvenile myoclonic epilepsy showing that with a predefined and limited set of mandatory features, variation in juvenile myoclonic epilepsy phenotype can be exploited for prognostic purposes. Our study is based on clinical data collected by the Biology of Juvenile Myoclonic Epilepsy Consortium augmented by literature data. We review prognosis research on mortality and seizure remission, predictors of antiseizure medication resistance and selected adverse drug events to valproate, levetiracetam and lamotrigine. Based on our analysis, a simplified set of diagnostic criteria for juvenile myoclonic epilepsy includes the following: (i) myoclonic jerks as mandatory seizure type; (ii) a circadian timing for myoclonia not mandatory for the diagnosis of juvenile myoclonic epilepsy; (iii) age of onset ranging from 6 to 40 years; (iv) generalized EEG abnormalities; and (v) intelligence conforming to population distribution. We find sufficient evidence to propose a predictive model of antiseizure medication resistance that emphasises (i) absence seizures as the strongest stratifying factor with regard to antiseizure medication resistance or seizure freedom for both sexes and (ii) sex as a major stratifying factor, revealing elevated odds of antiseizure medication resistance that correlates to self-report of catamenial and stress-related factors including sleep deprivation. In women, there are reduced odds of antiseizure medication resistance associated with EEG-measured or self-reported photosensitivity. In conclusion, by applying a simplified set of criteria to define phenotypic variations of juvenile myoclonic epilepsy, our paper proposes an evidence-based definition and prognostic stratification of juvenile myoclonic epilepsy. Further studies in existing data sets of individual patient data would be helpful to replicate our findings, and prospective studies in inception cohorts will contribute to validate them in real-world practice for juvenile myoclonic epilepsy management.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

HostSeq: a Canadian whole genome sequencing and clinical data resource

HostSeq was launched in April 2020 as a national initiative to integrate whole genome sequencing data from 10,000 Canadians infected with SARS-CoV-2 with clinical information related to their disease experience. The mandate of HostSeq is to support the Canadian and international research communities in their efforts to understand the risk factors for disease and associated health outcomes and support the development of interventions such as vaccines and therapeutics. HostSeq is a collaboration among 13 independent epidemiological studies of SARS-CoV-2 across five provinces in Canada. Aggregated data collected by HostSeq are made available to the public through two data portals: a phenotype portal showing summaries of major variables and their distributions, and a variant search portal enabling queries in a genomic region. Individual-level data is available to the global research community for health research through a Data Access Agreement and Data Access Compliance Office approval. Here we provide an overview of the collective project design along with summary level information for HostSeq. We highlight several statistical considerations for researchers using the HostSeq platform regarding data aggregation, sampling mechanism, covariate adjustment, and X chromosome analysis. In addition to serving as a rich data source, the diversity of study designs, sample sizes, and research objectives among the participating studies provides unique opportunities for the research community.

RoPE: A robust profile likelihood method for differential gene expression analysis

Variation in RNA-Seq data creates modeling challenges for differential gene expression (DE) analysis. Statistical approaches address conventional small sample sizes and implement empirical Bayes or non-parametric tests, but frequently produce different conclusions. Increasing sample sizes enable proposal of alternative DE paradigms. Here we develop RoPE, which uses a data-driven adjustment for variation and a robust profile likelihood ratio DE test. Simulation studies show RoPE can have improved performance over existing tools as sample size increases and has the most reliable control of error rates. Application of RoPE demonstrates that an active Pseudomonas aeruginosa infection downregulates the SLC9A3 Cystic Fibrosis modifier gene.

Identification of brain cell types underlying genetic association with word reading and correlated traits

Neuroimaging studies implicate multiple cortical regions in reading ability/disability. However, the neural cell types integral to the reading process are unknown. To contribute to this gap in knowledge, we integrated genetic results from genome-wide association studies for word reading (n = 5054) with gene expression datasets from adult/fetal human brain. Linkage disequilibrium score regression (LDSC) suggested that variants associated with word reading were enriched in genes expressed in adult excitatory neurons, specifically layer 5 and 6 FEZF2 expressing neurons and intratelencephalic (IT) neurons, which express the marker genes LINC00507, THEMIS, or RORB. Inhibitory neurons (VIP, SST, and PVALB) were also found. This finding was interesting as neurometabolite studies previously implicated excitatory-inhibitory imbalances in the etiology of reading disabilities (RD). We also tested traits that shared genetic etiology with word reading (previously determined by polygenic risk scores): attention-deficit/hyperactivity disorder (ADHD), educational attainment, and cognitive ability. For ADHD, we identified enrichment in L4 IT adult excitatory neurons. For educational attainment and cognitive ability, we confirmed previous studies identifying multiple subclasses of adult cortical excitatory and inhibitory neurons, as well as astrocytes and oligodendrocytes. For educational attainment and cognitive ability, we also identified enrichment in multiple fetal cortical excitatory and inhibitory neurons, intermediate progenitor cells, and radial glial cells. In summary, this study supports a role of excitatory and inhibitory neurons in reading and excitatory neurons in ADHD and contributes new information on fetal cell types enriched in educational attainment and cognitive ability, thereby improving our understanding of the neurobiological basis of reading/correlated traits.

Large scale genome-wide association analyses identify novel genetic loci and mechanisms in hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is an important cause of morbidity and mortality with both monogenic and polygenic components. We here report results from the largest HCM genome-wide association study (GWAS) and multi-trait analysis (MTAG) including 5,900 HCM cases, 68,359 controls, and 36,083 UK Biobank (UKB) participants with cardiac magnetic resonance (CMR) imaging. We identified a total of 70 loci (50 novel) associated with HCM, and 62 loci (32 novel) associated with relevant left ventricular (LV) structural or functional traits. Amongst the common variant HCM loci, we identify a novel HCM disease gene, SVIL, which encodes the actin-binding protein supervillin, showing that rare truncating SVIL variants cause HCM. Mendelian randomization analyses support a causal role of increased LV contractility in both obstructive and non-obstructive forms of HCM, suggesting common disease mechanisms and anticipating shared response to therapy. Taken together, the findings significantly increase our understanding of the genetic basis and molecular mechanisms of HCM, with potential implications for disease management.

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.

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.

Genomic architecture of autism from comprehensive whole-genome sequence annotation

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

Genome-wide analyses of individual differences in quantitatively assessed reading- and language-related skills in up to 34,000 people

The use of spoken and written language is a fundamental human capacity. Individual differences in reading- and language-related skills are influenced by genetic variation, with twin-based heritability estimates of 30 to 80% depending on the trait. The genetic architecture is complex, heterogeneous, and multifactorial, but investigations of contributions of single-nucleotide polymorphisms (SNPs) were thus far underpowered. We present a multicohort genome-wide association study (GWAS) of five traits assessed individually using psychometric measures (word reading, nonword reading, spelling, phoneme awareness, and nonword repetition) in samples of 13,633 to 33,959 participants aged 5 to 26 y. We identified genome-wide significant association with word reading (rs11208009, P = 1.098 × 10-8) at a locus that has not been associated with intelligence or educational attainment. All five reading-/language-related traits showed robust SNP heritability, accounting for 13 to 26% of trait variability. Genomic structural equation modeling revealed a shared genetic factor explaining most of the variation in word/nonword reading, spelling, and phoneme awareness, which only partially overlapped with genetic variation contributing to nonword repetition, intelligence, and educational attainment. A multivariate GWAS of word/nonword reading, spelling, and phoneme awareness maximized power for follow-up investigation. Genetic correlation analysis with neuroimaging traits identified an association with the surface area of the banks of the left superior temporal sulcus, a brain region linked to the processing of spoken and written language. Heritability was enriched for genomic elements regulating gene expression in the fetal brain and in chromosomal regions that are depleted of Neanderthal variants. Together, these results provide avenues for deciphering the biological underpinnings of uniquely human traits.

Heterogeneity of resting-state EEG features in juvenile myoclonic epilepsy and controls

Abnormal EEG features are a hallmark of epilepsy, and abnormal frequency and network features are apparent in EEGs from people with idiopathic generalized epilepsy in both ictal and interictal states. Here, we characterize differences in the resting-state EEG of individuals with juvenile myoclonic epilepsy and assess factors influencing the heterogeneity of EEG features. We collected EEG data from 147 participants with juvenile myoclonic epilepsy through the Biology of Juvenile Myoclonic Epilepsy study. Ninety-five control EEGs were acquired from two independent studies [Chowdhury et al. (2014) and EU-AIMS Longitudinal European Autism Project]. We extracted frequency and functional network-based features from 10 to 20 s epochs of resting-state EEG, including relative power spectral density, peak alpha frequency, network topology measures and brain network ictogenicity: a computational measure of the propensity of networks to generate seizure dynamics. We tested for differences between epilepsy and control EEGs using univariate, multivariable and receiver operating curve analysis. In addition, we explored the heterogeneity of EEG features within and between cohorts by testing for associations with potentially influential factors such as age, sex, epoch length and time, as well as testing for associations with clinical phenotypes including anti-seizure medication, and seizure characteristics in the epilepsy cohort. P-values were corrected for multiple comparisons. Univariate analysis showed significant differences in power spectral density in delta (2-5 Hz) (P = 0.0007, hedges' g = 0.55) and low-alpha (6-9 Hz) (P = 2.9 × 10-8, g = 0.80) frequency bands, peak alpha frequency (P = 0.000007, g = 0.66), functional network mean degree (P = 0.0006, g = 0.48) and brain network ictogenicity (P = 0.00006, g = 0.56) between epilepsy and controls. Since age (P = 0.009) and epoch length (P = 1.7 × 10-8) differed between the two groups and were potential confounders, we controlled for these covariates in multivariable analysis where disparities in EEG features between epilepsy and controls remained. Receiver operating curve analysis showed low-alpha power spectral density was optimal at distinguishing epilepsy from controls, with an area under the curve of 0.72. Lower average normalized clustering coefficient and shorter average normalized path length were associated with poorer seizure control in epilepsy patients. To conclude, individuals with juvenile myoclonic epilepsy have increased power of neural oscillatory activity at low-alpha frequencies, and increased brain network ictogenicity compared with controls, supporting evidence from studies in other epilepsies with considerable external validity. In addition, the impact of confounders on different frequency-based and network-based EEG features observed in this study highlights the need for careful consideration and control of these factors in future EEG research in idiopathic generalized epilepsy particularly for their use as biomarkers.

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.

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⁻⁵) 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.

Leveraging TOPMed Imputation Server and Constructing a Cohort-Specific Imputation Reference Panel to Enhance Genotype Imputation among Cystic Fibrosis Patients

Cystic fibrosis (CF) is a severe genetic disorder that can cause multiple comorbidities affecting the lungs, the pancreas, the luminal digestive system and beyond. In our previous genome-wide association studies (GWAS), we genotyped approximately 8,000 CF samples using a mixture of different genotyping platforms. More recently, the Cystic Fibrosis Genome Project (CFGP) performed deep (approximately 30×) whole genome sequencing (WGS) of 5,095 samples to better understand the genetic mechanisms underlying clinical heterogeneity among patients with CF. For mixtures of GWAS array and WGS data, genotype imputation has proven effective in increasing effective sample size. Therefore, we first performed imputation for the approximately 8,000 CF samples with GWAS array genotype using the Trans-Omics for Precision Medicine (TOPMed) freeze 8 reference panel. Our results demonstrate that TOPMed can provide high-quality imputation for patients with CF, boosting genomic coverage from approximately 0.3–4.2 million genotyped markers to approximately 11–43 million well-imputed markers, and significantly improving polygenic risk score (PRS) prediction accuracy. Furthermore, we built a CF-specific CFGP reference panel based on WGS data of patients with CF. We demonstrate that despite having approximately 3% the sample size of TOPMed, our CFGP reference panel can still outperform TOPMed when imputing some CF disease-causing variants, likely owing to allele and haplotype differences between patients with CF and general populations. We anticipate our imputed data for 4,656 samples without WGS data will benefit our subsequent genetic association studies, and the CFGP reference panel built from CF WGS samples will benefit other investigators studying CF.

Rare loss-of-function variants in type I IFN immunity genes are not associated with severe COVID-19

A recent report found that rare predicted loss-of-function (pLOF) variants across 13 candidate genes in TLR3- and IRF7-dependent type I IFN pathways explain up to 3.5% of severe COVID-19 cases. We performed whole-exome or whole-genome sequencing of 1,864 COVID-19 cases (713 with severe and 1,151 with mild disease) and 15,033 ancestry-matched population controls across 4 independent COVID-19 biobanks. We tested whether rare pLOF variants in these 13 genes were associated with severe COVID-19. We identified only 1 rare pLOF mutation across these genes among 713 cases with severe COVID-19 and observed no enrichment of pLOFs in severe cases compared to population controls or mild COVID-19 cases. We found no evidence of association of rare LOF variants in the 13 candidate genes with severe COVID-19 outcomes.

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.

Genome-wide association study of pediatric obsessive-compulsive traits: shared genetic risk between traits and disorder

Using a novel trait-based measure, we examined genetic variants associated with obsessive-compulsive (OC) traits and tested whether OC traits and obsessive-compulsive disorder (OCD) shared genetic risk. We conducted a genome-wide association analysis (GWAS) of OC traits using the Toronto Obsessive-Compulsive Scale (TOCS) in 5018 unrelated Caucasian children and adolescents from the community (Spit for Science sample). We tested the hypothesis that genetic variants associated with OC traits from the community would be associated with clinical OCD using a meta-analysis of all currently available OCD cases. Shared genetic risk was examined between OC traits and OCD in the respective samples using polygenic risk score and genetic correlation analyses. A locus tagged by rs7856850 in an intron of PTPRD (protein tyrosine phosphatase δ) was significantly associated with OC traits at the genome-wide significance level (p = 2.48 × 10-8). rs7856850 was also associated with OCD in a meta-analysis of OCD case/control genome-wide datasets (p = 0.0069). The direction of effect was the same as in the community sample. Polygenic risk scores from OC traits were significantly associated with OCD in case/control datasets and vice versa (p's < 0.01). OC traits were highly, but not significantly, genetically correlated with OCD (rg = 0.71, p = 0.062). We report the first validated genome-wide significant variant for OC traits in PTPRD, downstream of the most significant locus in a previous OCD GWAS. OC traits measured in the community sample shared genetic risk with OCD case/control status. Our results demonstrate the feasibility and power of using trait-based approaches in community samples for genetic discovery.

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.

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

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

Statistical power in COVID-19 case-control host genomic study design

The identification of genetic variation that directly impacts infection susceptibility to SARS-CoV-2 and disease severity of COVID-19 is an important step towards risk stratification, personalized treatment plans, therapeutic, and vaccine development and deployment. Given the importance of study design in infectious disease genetic epidemiology, we use simulation and draw on current estimates of exposure, infectivity, and test accuracy of COVID-19 to demonstrate the feasibility of detecting host genetic factors associated with susceptibility and severity in published COVID-19 study designs. We demonstrate that limited phenotypic data and exposure/infection information in the early stages of the pandemic significantly impact the ability to detect most genetic variants with moderate effect sizes, especially when studying susceptibility to SARS-CoV-2 infection. Our insights can aid in the interpretation of genetic findings emerging in the literature and guide the design of future host genetic studies.

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.

A Distributed Whole Genome Sequencing Benchmark Study

Population sequencing often requires collaboration across a distributed network of sequencing centers for the timely processing of thousands of samples. In such massive efforts, it is important that participating scientists can be confident that the accuracy of the sequence data produced is not affected by which center generates the data. A study was conducted across three established sequencing centers, located in Montreal, Toronto, and Vancouver, constituting Canada's Genomics Enterprise (www.cgen.ca). Whole genome sequencing was performed at each center, on three genomic DNA replicates from three well-characterized cell lines. Secondary analysis pipelines employed by each site were applied to sequence data from each of the sites, resulting in three datasets for each of four variables (cell line, replicate, sequencing center, and analysis pipeline), for a total of 81 datasets. These datasets were each assessed according to multiple quality metrics including concordance with benchmark variant truth sets to assess consistent quality across all three conditions for each variable. Three-way concordance analysis of variants across conditions for each variable was performed. Our results showed that the variant concordance between datasets differing only by sequencing center was similar to the concordance for datasets differing only by replicate, using the same analysis pipeline. We also showed that the statistically significant differences between datasets result from the analysis pipeline used, which can be unified and updated as new approaches become available. We conclude that genome sequencing projects can rely on the quality and reproducibility of aggregate data generated across a network of distributed sites.

Mining GWAS and eQTL data for CF lung disease modifiers by gene expression imputation

Genome wide association studies (GWAS) have identified several genomic loci with candidate modifiers of cystic fibrosis (CF) lung disease, but only a small proportion of the expected genetic contribution is accounted for at these loci. We leveraged expression data from CF cohorts, and Genotype-Tissue Expression (GTEx) reference data sets from multiple human tissues to generate predictive models, which were used to impute transcriptional regulation from genetic variance in our GWAS population. The imputed gene expression was tested for association with CF lung disease severity. By comparing and combining results from alternative approaches, we identified 379 candidate modifier genes. We delved into 52 modifier candidates that showed consensus between approaches, and 28 of them were near known GWAS loci. A number of these genes are implicated in the pathophysiology of CF lung disease (e.g., immunity, infection, inflammation, HLA pathways, glycosylation, and mucociliary clearance) and the CFTR protein biology (e.g., cytoskeleton, microtubule, mitochondrial function, lipid metabolism, endoplasmic reticulum/Golgi, and ubiquitination). Gene set enrichment results are consistent with current knowledge of CF lung disease pathogenesis. HLA Class II genes on chr6, and CEP72, EXOC3, and TPPP near the GWAS peak on chr5 are most consistently associated with CF lung disease severity across the tissues tested. The results help to prioritize genes in the GWAS regions, predict direction of gene expression regulation, and identify new candidate modifiers throughout the genome for potential therapeutic development.

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

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

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

CONTEXT

Individuals with cystic fibrosis (CF) develop a distinct form of diabetes characterized by β-cell dysfunction and islet amyloid accumulation similar to type 2 diabetes (T2D), but generally have normal insulin sensitivity. CF-related diabetes (CFRD) risk is determined by both CFTR, the gene responsible for CF, and other genetic variants.

OBJECTIVE

To identify genetic modifiers of CFRD and determine the genetic overlap with other types of diabetes.

DESIGN AND PATIENTS

A genome-wide association study was conducted for CFRD onset on 5740 individuals with CF. Weighted polygenic risk scores (PRSs) for type 1 diabetes (T1D), T2D, and diabetes endophenotypes were tested for association with CFRD.

RESULTS

Genome-wide significance was obtained for variants at a novel locus (PTMA) and 2 known CFRD genetic modifiers (TCF7L2 and SLC26A9). PTMA and SLC26A9 variants were CF-specific; TCF7L2 variants also associated with T2D. CFRD was strongly associated with PRSs for T2D, insulin secretion, postchallenge glucose concentration, and fasting plasma glucose, and less strongly with T1D PRSs. CFRD was inconsistently associated with PRSs for insulin sensitivity and was not associated with a PRS for islet autoimmunity. A CFRD PRS comprising variants selected from these PRSs (with a false discovery rate < 0.1) and the genome-wide significant variants was associated with CFRD in a replication population.

CONCLUSIONS

CFRD and T2D have more etiologic and mechanistic overlap than previously known, aligning along pathways involving β-cell function rather than insulin sensitivity. Two CFRD risk loci are unrelated to T2D and may affect multiple aspects of CF. An 18-variant PRS stratifies risk of CFRD in an independent population.

Genome-wide association study of word reading: Overlap with risk genes for neurodevelopmental disorders

Reading disabilities (RD) are the most common neurocognitive disorder, affecting 5% to 17% of children in North America. These children often have comorbid neurodevelopmental/psychiatric disorders, such as attention deficit/hyperactivity disorder (ADHD). The genetics of RD and their overlap with other disorders is incompletely understood. To contribute to this, we performed a genome-wide association study (GWAS) for word reading. Then, using summary statistics from neurodevelopmental/psychiatric disorders, we computed polygenic risk scores (PRS) and used them to predict reading ability in our samples. This enabled us to test the shared aetiology between RD and other disorders. The GWAS consisted of 5.3 million single nucleotide polymorphisms (SNPs) and two samples; a family-based sample recruited for reading difficulties in Toronto (n = 624) and a population-based sample recruited in Philadelphia [Philadelphia Neurodevelopmental Cohort (PNC)] (n = 4430). The Toronto sample SNP-based analysis identified suggestive SNPs (P ~ 5 × 10^-7 ) in the ARHGAP23 gene, which is implicated in neuronal migration/axon pathfinding. The PNC gene-based analysis identified significant associations (P < 2.72 × 10^-6 ) for LINC00935 and CCNT1, located in the region of the KANSL2/CCNT1/LINC00935/SNORA2B/SNORA34/MIR4701/ADCY6 genes on chromosome 12q, with near significant SNP-based analysis. PRS identified significant overlap between word reading and intelligence (R²  = 0.18, P = 7.25  × 10^-181 ), word reading and educational attainment (R²  = 0.07, P = 4.91 × 10^-48 ) and word reading and ADHD (R²  = 0.02, P = 8.70 × 10^-6 ; threshold for significance = 7.14 × 10^-3 ). Overlap was also found between RD and autism spectrum disorder (ASD) as top-ranked genes were previously implicated in autism by rare and copy number variant analyses. These findings support shared risk between word reading, cognitive measures, educational outcomes and neurodevelopmental disorders, including ASD.

VikNGS: a C++ variant integration kit for next generation sequencing association analysis

SUMMARY

Integration of next generation sequencing data (NGS) across different research studies can improve the power of genetic association testing by increasing sample size and can obviate the need for sequencing controls. If differential genotype uncertainty across studies is not accounted for, combining datasets can produce spurious association results. We developed the Variant Integration Kit for NGS (VikNGS), a fast cross-platform software package, to enable aggregation of several datasets for rare and common variant genetic association analysis of quantitative and binary traits with covariate adjustment. VikNGS also includes a graphical user interface, power simulation functionality and data visualization tools.

AVAILABILITY AND IMPLEMENTATION

The VikNGS package can be downloaded at http://www.tcag.ca/tools/index.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Immunomodulatory function of the cystic fibrosis modifier gene BPIFA1

BACKGROUND

Cystic fibrosis (CF) is characterized by a progressive decline in lung function due to airway obstruction, infection, and inflammation. CF patients are particularly susceptible to respiratory infection by a variety of pathogens, and the inflammatory response in CF is dysregulated and prolonged. BPI fold containing family A, member 1 (BPIFA1) and BPIFB1 are proteins expressed in the upper airways that may have innate immune activity. We previously identified polymorphisms in the BPIFA1/BPIFB1 region associated with CF lung disease severity.

METHODS

We evaluated whether the BPIFA1/BPIFB1 associations with lung disease severity replicated in individuals with CF participating in the International CF Gene Modifier Consortium (n = 6,365). Furthermore, we investigated mechanisms by which the BPIFA1 and BPIFB1 proteins may modify lung disease in CF.

RESULTS

The association of the G allele of rs1078761 with reduced lung function was replicated in an independent cohort of CF patients (p = 0.001, n = 2,921) and in a meta-analysis of the full consortium (p = 2.39x10-5, n = 6,365). Furthermore, we found that rs1078761G which is associated with reduced lung function was also associated with reduced BPIFA1, but not BPIFB1, protein levels in saliva from CF patients. Functional assays indicated that BPIFA1 and BPIFB1 do not have an anti-bacterial role against P. aeruginosa but may have an immunomodulatory function in CF airway epithelial cells. Gene expression profiling using RNAseq identified Rho GTPase signaling pathways to be altered in CF airway epithelial cells in response to treatment with recombinant BPIFA1 and BPIFB1 proteins.

CONCLUSIONS

BPIFA1 and BPIFB1 have immunomodulatory activity and genetic variation associated with low levels of these proteins may increase CF lung disease severity.

Impact of DNA source on genetic variant detection from human whole-genome sequencing data

BACKGROUND

Whole blood is currently the most common DNA source for whole-genome sequencing (WGS), but for studies requiring non-invasive collection, self-collection, greater sample stability or additional tissue references, saliva or buccal samples may be preferred. However, the relative quality of sequencing data and accuracy of genetic variant detection from blood-derived, saliva-derived and buccal-derived DNA need to be thoroughly investigated.

METHODS

Matched blood, saliva and buccal samples from four unrelated individuals were used to compare sequencing metrics and variant-detection accuracy among these DNA sources.

RESULTS

We observed significant differences among DNA sources for sequencing quality metrics such as percentage of reads aligned and mean read depth (p<0.05). Differences were negligible in the accuracy of detecting short insertions and deletions; however, the false positive rate for single nucleotide variation detection was slightly higher in some saliva and buccal samples. The sensitivity of copy number variant (CNV) detection was up to 25% higher in blood samples, depending on CNV size and type, and appeared to be worse in saliva and buccal samples with high bacterial concentration. We also show that methylation-based enrichment for eukaryotic DNA in saliva and buccal samples increased alignment rates but also reduced read-depth uniformity, hampering CNV detection.

CONCLUSION

For WGS, we recommend using DNA extracted from blood rather than saliva or buccal swabs; if saliva or buccal samples are used, we recommend against using methylation-based eukaryotic DNA enrichment. All data used in this study are available for further open-science investigation.

Airway Mucosal Host Defense Is Key to Genomic Regulation of Cystic Fibrosis Lung Disease Severity

RATIONALE

The severity of cystic fibrosis (CF) lung disease varies widely, even for Phe508del homozygotes. Heritability studies show that more than 50% of the variability reflects non-cystic fibrosis transmembrane conductance regulator (CFTR) genetic variation; however, the full extent of the pertinent genetic variation is not known.

OBJECTIVES

We sought to identify novel CF disease-modifying mechanisms using an integrated approach based on analyzing "in vivo" CF airway epithelial gene expression complemented with genome-wide association study (GWAS) data.

METHODS

Nasal mucosal RNA from 134 patients with CF was used for RNA sequencing. We tested for associations of transcriptomic (gene expression) data with a quantitative phenotype of CF lung disease severity. Pathway analysis of CF GWAS data (n = 5,659 patients) was performed to identify novel pathways and assess the concordance of genomic and transcriptomic data. Association of gene expression with previously identified CF GWAS risk alleles was also tested.

MEASUREMENTS AND MAIN RESULTS

Significant evidence of heritable gene expression was identified. Gene expression pathways relevant to airway mucosal host defense were significantly associated with CF lung disease severity, including viral infection, inflammation/inflammatory signaling, lipid metabolism, apoptosis, ion transport, Phe508del CFTR processing, and innate immune responses, including HLA (human leukocyte antigen) genes. Ion transport and CFTR processing pathways, as well as HLA genes, were identified across differential gene expression and GWAS signals.

CONCLUSIONS

Transcriptomic analyses of CF airway epithelia, coupled to genomic (GWAS) analyses, highlight the role of heritable host defense variation in determining the pathophysiology of CF lung disease. The identification of these pathways provides opportunities to pursue targeted interventions to improve CF lung health.

The evidential statistical paradigm in genetics

Concerns over reproducibility in research has reinvigorated the discourse on P-values as measures of statistical evidence. In a position statement by the American Statistical Association board of directors, they warn of P-value misuse and refer to the availability of alternatives. Despite the common practice of comparing P-values across different hypothesis tests in genetics, it is well-appreciated that P-values must be interpreted alongside the sample size and experimental design used for their computation. Here, we discuss the evidential statistical paradigm (EP), an alternative to Bayesian and Frequentist paradigms, that has been implemented in human genetics studies. Using applications in Cystic Fibrosis genetic association analyses, and describing recent theoretical developments, we review how to measure statistical evidence using the EP in the presence of covariates, model misspecification, and for composite hypotheses. Novel graphical displays are presented, and software for their computation is highlighted. The implications of multiple hypothesis testing for the EP are delineated in the analyses, demonstrating a view more consistent with scientific reasoning; the EP provides a theoretical justification for replication that is a requirement in genetic association studies. As genetic studies grow in size and complexity, a fresh look at measures of statistical evidence that are sensible amid the analysis of big data are required.

Genetic association analysis with pedigrees: Direct inference using the composite likelihood ratio

The likelihood function represents statistical evidence given data and a model. The evidential paradigm (EP), an alternative to Bayesian and Frequentist paradigms, provides considerable theory demonstrating evidence strength for different parameter values via the ratio of likelihoods at different parameter values; thus, enabling inference directly from the likelihood function. The likelihood function, however, can be difficult to compute; for example, in genetic association studies with a binary outcome in large pedigrees. Composite likelihood (CL) is an alternative when the real likelihood is intractable. We show CLs have the two large sample properties of the EP for reliable evidence interpretation: (1) CL supports the true value over a false value by an arbitrarily large factor; and (2) the probability of favouring a false value over the true value is small and bounded. Using simulation, and in a genetic association analysis of reading disability (RD) in large rolandic epilepsy pedigrees, we show that the CL approach yields valid statistical inference and identifies RD associated variants. When compared to analyses using generalized estimating equations, results show a similar prioritization of SNPs, although the CL approach provides additional complementary information, and more intuitive solutions to the multiple hypothesis testing problem.

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.

SLC26A9 Gene Is Associated With Lung Function Response to Ivacaftor in Patients With Cystic Fibrosis

Ivacaftor is a drug used to treat cystic fibrosis (CF) patients carrying specific gating CFTR mutations. Interpatient variability in the lung response has been shown to be partly explained by rs7512462 in the Solute Carrier Family 26 Member 9 (SLC26A9) gene. In an independent and larger cohort, we aimed to evaluate whether SLC26A9 variants contribute to the variability of the lung phenotype and if they influence the lung response to ivacaftor. We genotyped the French CF Gene Modifier Study cohort (n = 4,840) to investigate whether SLC26A9 variants were involved in the lung phenotype heterogeneity. Their influence in the response to ivacaftor was tested in the 30 treated patients who met the inclusion criteria: older than 6 years of age, percent-predicted forced expiratory volume measured in 1 s (FEV_1pp) in the 3 months before treatment initiation ranging between 40 and 90%. Response to treatment was determined by the change in FEV_1pp from baseline, averaged in 15–75 days, and the 1st-year post-treatment. We observed that SLC26A9 variants were not associated with lung function variability in untreated patients and that gain of lung function in patients treated with ivacaftor was similar to clinical trials. We confirmed that rs7512462 was associated with variability in ivacaftor-lung response, with a significant reduction in lung function improvement for patients with the C allele. Other SLC26A9 SNPs also contributed to the ivacaftor-response. Interindividual variability in lung response to ivacaftor is associated with SLC26A9 variants in French CF patients. Pharmacogenomics and personalized medicine will soon be part of CF patient care.

AGTR2 absence or antagonism prevents cystic fibrosis pulmonary manifestations

BACKGROUND

Pulmonary disease remains the primary cause of morbidity and mortality for individuals with cystic fibrosis (CF). Variants at a locus on the X-chromosome containing the type 2 angiotensin II receptor gene (AGTR2) were identified by a large GWAS as significantly associating with lung function in CF patients. We hypothesized that manipulating the angiotensin-signaling pathway may yield clinical benefit in CF.

METHODS

Genetic subset analysis was conducted on a local CF cohort to extend the GWAS findings. Next, we evaluated pulmonary function in CF mice with a deleted AGTR2 gene, and in those who were given subcutaneous injections of PD123,319, a selective AGTR2 antagonist for 12 weeks beginning at weaning.

RESULTS

The genetic subset analysis replicated the initial GWAS identified association, and confirmed the association of this locus with additional lung function parameters. Studies in genetically modified mice established that absence of the AGTR2 gene normalized pulmonary function indices in two independent CF mouse models. Further, we determined that pharmacologic antagonism of AGTR2 improved overall pulmonary function in CF mice to near wild-type levels.

CONCLUSIONS

These results identify that reduced AGTR2 signaling is beneficial to CF lung function, and suggest the potential of manipulating the angiotensin-signaling pathway for treatment and/or prevention of CF pulmonary disease. Importantly, the beneficial effects were not CF gene mutation dependent, and were able to be reproduced with pharmacologic antagonism. As there are clinically approved drugs available to target the renin-angiotensin signaling system, these findings may be quickly translated to human clinical trials.

Testing Calibration of Cox Survival Models at Extremes of Event Risk

Risk prediction models can translate genetic association findings for clinical decision-making. Most models are evaluated on their ability to discriminate, and the calibration of risk-prediction models is largely overlooked in applications. Models that demonstrate good discrimination in training datasets, if not properly calibrated to produce unbiased estimates of risk, can perform poorly in new patient populations. Poorly calibrated models arise due to missing covariates, such as genetic interactions that may be unknown or not measured. We demonstrate that models omitting interactions can lead to increased bias in predicted risk for patients at the tails of the risk distribution; i.e., those patients who are most likely to be affected by clinical decision making. We propose a new calibration test for Cox risk-prediction models that aggregates martingale residuals for subjects from extreme high and low risk groups with a test statistic maximum chosen by varying which risk groups are included in the extremes. To estimate the empirical significance of our test statistic, we simulate from a Gaussian distribution using the covariance matrix for the grouped sums of martingale residuals. Simulation shows the new test maintains control of type 1 error with improved power over a conventional goodness-of-fit test when risk prediction deviates at the tails of the risk distribution. We apply our method in the development of a prediction model for risk of cystic fibrosis-related diabetes. Our study highlights the importance of assessing calibration and discrimination in predictive modeling, and provides a complementary tool in the assessment of risk model calibration.

Reference percentiles of FEV1 for the Canadian cystic fibrosis population: comparisons across time and countries

BACKGROUND

Forced expiratory volume in 1 s (FEV1) indicates lung health in cystic fibrosis (CF). FEV1 is commonly communicated as a per cent predicted of a healthy individual sharing the same age, sex, race and height. CF-specific reference equations are complementary and calibrate a patient's FEV1 to that of their CF peers.

OBJECTIVES

(1) To derive Canadian CF-specific FEV1 reference percentiles (FEV1%iles), (2) characterize how they have changed over time and (3) compare the Canadian FEV1%iles to those for USA and European CF populations.

METHOD

CF FEV1%iles are calculated using the Canadian CF Registry and quantile regression.

RESULTS

The Canadian FEV1%iles demonstrated better lung function in more recent time periods within Canada, especially below the 50% percentile and in males. When compared to USA and European FEV1%iles for the same time period, Canadian FEV1%iles were higher.

CONCLUSION

CF-specific FEV1%iles can provide useful information about changes in lung health. An online calculator (available at cfpercentile.

RESEARCH

sickkids.ca) makes these FEV1%iles accessible.

The Personal Genome Project Canada: findings from whole genome sequences of the inaugural 56 participants

BACKGROUND

The Personal Genome Project Canada is a comprehensive public data resource that integrates whole genome sequencing data and health information. We describe genomic variation identified in the initial recruitment cohort of 56 volunteers.

METHODS

Volunteers were screened for eligibility and provided informed consent for open data sharing. Using blood DNA, we performed whole genome sequencing and identified all possible classes of DNA variants. A genetic counsellor explained the implication of the results to each participant.

RESULTS

Whole genome sequencing of the first 56 participants identified 207 662 805 sequence variants and 27 494 copy number variations. We analyzed a prioritized disease-associated data set (n = 1606 variants) according to standardized guidelines, and interpreted 19 variants in 14 participants (25%) as having obvious health implications. Six of these variants (e.g., in BRCA1 or mosaic loss of an X chromosome) were pathogenic or likely pathogenic. Seven were risk factors for cancer, cardiovascular or neurobehavioural conditions. Four other variants - associated with cancer, cardiac or neurodegenerative phenotypes - remained of uncertain significance because of discrepancies among databases. We also identified a large structural chromosome aberration and a likely pathogenic mitochondrial variant. There were 172 recessive disease alleles (e.g., 5 individuals carried mutations for cystic fibrosis). Pharmacogenomics analyses revealed another 3.9 potentially relevant genotypes per individual.

INTERPRETATION

Our analyses identified a spectrum of genetic variants with potential health impact in 25% of participants. When also considering recessive alleles and variants with potential pharmacologic relevance, all 56 participants had medically relevant findings. Although access is mostly limited to research, whole genome sequencing can provide specific and novel information with the potential of major impact for health care.

Endoplasmic Reticulum Stress and Chemokine Production in Cystic Fibrosis Airway Cells: Regulation by STAT3 Modulation

Endoplasmic reticulum (ER) stress has been recognized to play an important role in chronic inflammatory diseases such as cystic fibrosis (CF), and targeting ER stress may be useful for alleviating damaging neutrophilic inflammation in CF airways. Cellular models were used in conjunction with data from a recent CF genome-wide association study (GWAS) meta-analysis to determine modulators of ER stress-mediated inflammation. Surprisingly, cells undergoing ER stress during inflammatory stimulation showed reduced interleukin 8 (IL-8) and CXCL1 secretion (P < .001). Neutralization of CXCL1 and IL-8 reduced neutrophil chemotaxis >50% to supernatants from IL-1β-stimulated CF airway epithelial cells (P < .01). The clinical importance of these chemokines was validated by association of CXCL1 and IL8 polymorphisms with changes in lung disease severity in patients with CF (n = 6365; IL8, P = .001; CXCL1, P = .001), confirming that targeting these chemokine pathways could help improve lung disease. We determined that production of these chemokines was partially controlled by ER stress in a signal transducer and activator of transcription 3 (STAT3)-dependent manner, whereby ER stress inhibited STAT3 activation. Our findings support a role for CXCL1 and IL-8 in CF lung disease severity and identify STAT3 as a modulating pathway. Targeting these pathways may help improve health outcomes in CF.

Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder

We are performing whole-genome sequencing of families with autism spectrum disorder (ASD) to build a resource (MSSNG) for subcategorizing the phenotypes and underlying genetic factors involved. Here we report sequencing of 5,205 samples from families with ASD, accompanied by clinical information, creating a database accessible on a cloud platform and through a controlled-access internet portal. We found an average of 73.8 de novo single nucleotide variants and 12.6 de novo insertions and deletions or copy number variations per ASD subject. We identified 18 new candidate ASD-risk genes and found that participants bearing mutations in susceptibility genes had significantly lower adaptive ability (P = 6 × 10^-4). In 294 of 2,620 (11.2%) of ASD cases, a molecular basis could be determined and 7.2% of these carried copy number variations and/or chromosomal abnormalities, emphasizing the importance of detecting all forms of genetic variation as diagnostic and therapeutic targets in ASD.

Sources of Variation in Sweat Chloride Measurements in Cystic Fibrosis

RATIONALE

Expanding the use of cystic fibrosis transmembrane conductance regulator (CFTR) potentiators and correctors for the treatment of cystic fibrosis (CF) requires precise and accurate biomarkers. Sweat chloride concentration provides an in vivo assessment of CFTR function, but it is unknown the degree to which CFTR mutations account for sweat chloride variation.

OBJECTIVES

To estimate potential sources of variation for sweat chloride measurements, including demographic factors, testing variability, recording biases, and CFTR genotype itself.

METHODS

A total of 2,639 sweat chloride measurements were obtained in 1,761 twins/siblings from the CF Twin-Sibling Study, French CF Modifier Gene Study, and Canadian Consortium for Genetic Studies. Variance component estimation was performed by nested mixed modeling.

MEASUREMENTS AND MAIN RESULTS

Across the tested CF population as a whole, CFTR gene mutations were found to be the primary determinant of sweat chloride variability (56.1% of variation) with contributions from variation over time (e.g., factors related to testing on different days; 13.8%), environmental factors (e.g., climate, family diet; 13.5%), other residual factors (e.g., test variability; 9.9%), and unique individual factors (e.g., modifier genes, unique exposures; 6.8%) (likelihood ratio test, P < 0.001). Twin analysis suggested that modifier genes did not play a significant role because the heritability estimate was negligible (H2 = 0; 95% confidence interval, 0.0-0.35). For an individual with CF, variation in sweat chloride was primarily caused by variation over time (58.1%) with the remainder attributable to residual/random factors (41.9%).

CONCLUSIONS

Variation in the CFTR gene is the predominant cause of sweat chloride variation; most of the non-CFTR variation is caused by testing variability and unique environmental factors. If test precision and accuracy can be improved, sweat chloride measurement could be a valuable biomarker for assessing response to therapies directed at mutant CFTR.

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

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