High Risk Population Isolate Reveals Low Frequency Variants Predisposing to Intracranial Aneurysms

3% of the population develops saccular intracranial aneurysms (sIAs), a complex trait, with a sporadic and a familial form. Subarachnoid hemorrhage from sIA (sIA-SAH) is a devastating form of stroke. Certain rare genetic variants are enriched in the Finns, a population isolate with a small founder population and bottleneck events. As the sIA-SAH incidence in Finland is >2× increased, such variants may associate with sIA in the Finnish population. We tested 9.4 million variants for association in 760 Finnish sIA patients (enriched for familial sIA), and in 2,513 matched controls with case-control status and with the number of sIAs. The most promising loci (p<5E-6) were replicated in 858 Finnish sIA patients and 4,048 controls. The frequencies and effect sizes of the replicated variants were compared to a continental European population using 717 Dutch cases and 3,004 controls. We discovered four new high-risk loci with low frequency lead variants. Three were associated with the case-control status: 2q23.3 (MAF 2.1%, OR 1.89, p 1.42×10-9); 5q31.3 (MAF 2.7%, OR 1.66, p 3.17×10-8); 6q24.2 (MAF 2.6%, OR 1.87, p 1.87×10-11) and one with the number of sIAs: 7p22.1 (MAF 3.3%, RR 1.59, p 6.08×-9). Two of the associations (5q31.3, 6q24.2) replicated in the Dutch sample. The 7p22.1 locus was strongly differentiated; the lead variant was more frequent in Finland (4.6%) than in the Netherlands (0.3%). Additionally, we replicated a previously inconclusive locus on 2q33.1 in all samples tested (OR 1.27, p 1.87×10-12). The five loci explain 2.1% of the sIA heritability in Finland, and may relate to, but not explain, the increased incidence of sIA-SAH in Finland. This study illustrates the utility of population isolates, familial enrichment, dense genotype imputation and alternate phenotyping in search for variants associated with complex diseases.Public Library of Science open acces

Genetic associations of protein-coding variants in human disease

Genome-wide association studies (GWAS) have identified thousands of genetic variants linked to the risk of human disease. However, GWAS have so far remained largely underpowered in relation to identifying associations in the rare and low-frequency allelic spectrum and have lacked the resolution to trace causal mechanisms to underlying genes(1). Here we combined whole-exome sequencing in 392,814 UK Biobank participants with imputed genotypes from 260,405 FinnGen participants (653,219 total individuals) to conduct association meta-analyses for 744 disease endpoints across the protein-coding allelic frequency spectrum, bridging the gap between common and rare variant studies. We identified 975 associations, with more than one-third being previously unreported. We demonstrate population-level relevance for mutations previously ascribed to causing single-gene disorders, map GWAS associations to likely causal genes, explain disease mechanisms, and systematically relate disease associations to levels of 117 biomarkers and clinical-stage drug targets. Combining sequencing and genotyping in two population biobanks enabled us to benefit from increased power to detect and explain disease associations, validate findings through replication and propose medical actionability for rare genetic variants. Our study provides a compendium of protein-coding variant associations for future insights into disease biology and drug discovery. A meta-analysis combining whole-exome sequencing data from UK Biobank participants and imputed genotypes from FinnGen participants enables identification of genetic associations with human disease in the rare and low-frequency allelic spectrumPeer reviewe

TAFFEL: Independent Enrichment Analysis of gene sets

<p>Abstract</p> <p>Background</p> <p>A major challenge in genomic research is identifying significant biological processes and generating new hypotheses from large gene sets. Gene sets often consist of multiple separate biological pathways, controlled by distinct regulatory mechanisms. Many of these pathways and the associated regulatory mechanisms might be obscured by a large number of other significant processes and thus not identified as significant by standard gene set enrichment analysis tools.</p> <p>Results</p> <p>We present a novel method called Independent Enrichment Analysis (IEA) and software TAFFEL that eases the task by clustering genes to subgroups using Gene Ontology categories and transcription regulators. IEA indicates transcriptional regulators putatively controlling biological functions in studied condition.</p> <p>Conclusions</p> <p>We demonstrate that the developed method and TAFFEL tool give new insight to the analysis of differentially expressed genes and can generate novel hypotheses. Our comparison to other popular methods showed that the IEA method implemented in TAFFEL can find important biological phenomena, which are not reported by other methods.</p

Efficient and accurate frailty model approach for genome-wide survival association analysis in large-scale biobanks

With decades of electronic health records linked to genetic data, large biobanks provide unprecedented opportunities for systematically understanding the genetics of the natural history of complex diseases. Genome-wide survival association analysis can identify genetic variants associated with ages of onset, disease progression and lifespan. We propose an efficient and accurate frailty model approach for genome-wide survival association analysis of censored time-to-event (TTE) phenotypes by accounting for both population structure and relatedness. Our method utilizes state-of-the-art optimization strategies to reduce the computational cost. The saddlepoint approximation is used to allow for analysis of heavily censored phenotypes (>90%) and low frequency variants (down to minor allele count 20). We demonstrate the performance of our method through extensive simulation studies and analysis of five TTE phenotypes, including lifespan, with heavy censoring rates (90.9% to 99.8%) on similar to 400,000 UK Biobank participants with white British ancestry and similar to 180,000 individuals in FinnGen. We further analyzed 871 TTE phenotypes in the UK Biobank and presented the genome-wide scale phenome-wide association results with the PheWeb browser.Peer reviewe

Novel patients with NHLRC2 variants expand the phenotypic spectrum of FINCA disease

PurposeFINCA disease (Fibrosis, Neurodegeneration and Cerebral Angiomatosis, OMIM 618278) is an infantile-onset neurodevelopmental and multiorgan disease. Since our initial report in 2018, additional patients have been described. FINCA is the first human disease caused by recessive variants in the highly conserved NHLRC2 gene. Our previous studies have shown that Nhlrc2-null mouse embryos die during gastrulation, indicating the essential role of the protein in embryonic development. Defect in NHLRC2 leads to cerebral neurodegeneration and severe pulmonary, hepatic and cardiac fibrosis. Despite having a structure suggestive of an enzymatic role and the clinical importance of NHLRC2 in multiple organs, the specific physiological role of the protein is unknown.MethodsThe clinical histories of five novel FINCA patients diagnosed with whole exome sequencing were reviewed. Segregation analysis of the biallelic, potentially pathogenic NHLRC2 variants was performed using Sanger sequencing. Studies on neuropathology and NHLRC2 expression in different brain regions were performed on autopsy samples of three previously described deceased FINCA patients.ResultsOne patient was homozygous for the pathogenic variant c.442G &gt; T, while the other four were compound heterozygous for this variant and two other pathogenic NHLRC2 gene variants. All five patients presented with multiorgan dysfunction with neurodevelopmental delay, recurrent infections and macrocytic anemia as key features. Interstitial lung disease was pronounced in infancy but often stabilized. Autopsy samples revealed widespread, albeit at a lower intensity than the control, NHLRC2 expression in the brain.ConclusionThis report expands on the characteristic clinical features of FINCA disease. Presentation is typically in infancy, and although patients can live to late adulthood, the key clinical and histopathological features are fibrosis, infection susceptibility/immunodeficiency/intellectual disability, neurodevelopmental disorder/neurodegeneration and chronic anemia/cerebral angiomatosis (hence the acronym FINCA) that enable an early diagnosis confirmed by genetic investigations

Migraine polygenic risk score associates with efficacy of migraine-specific drugs

Objective To assess whether the polygenic risk score (PRS) for migraine is associated with acute and/or prophylactic migraine treatment response. Methods We interviewed 2,219 unrelated patients at the Danish Headache Center using a semistructured interview to diagnose migraine and assess acute and prophylactic drug response. All patients were genotyped. A PRS was calculated with the linkage disequilibrium pred algorithm using summary statistics from the most recent migraine genome-wide association study comprising similar to 375,000 cases and controls. The PRS was scaled to a unit corresponding to a twofold increase in migraine risk, using 929 unrelated Danish controls as reference. The association of the PRS with treatment response was assessed by logistic regression, and the predictive power of the model by area under the curve using a case-control design with treatment response as outcome. Results A twofold increase in migraine risk associates with positive response to migraine-specific acute treatment (odds ratio [OR] = 1.25 [95% confidence interval (CI) = 1.05-1.49]). The association between migraine risk and migraine-specific acute treatment was replicated in an independent cohort consisting of 5,616 triptan users with prescription history (OR = 3.20 [95% CI = 1.26-8.14]). No association was found for acute treatment with non-migraine-specific weak analgesics and prophylactic treatment response. Conclusions The migraine PRS can significantly identify subgroups of patients with a higher-than-average likelihood of a positive response to triptans, which provides a first step toward genetics-based precision medicine in migraine.Peer reviewe

Meta-analysis fine-mapping is often miscalibrated at single-variant resolution

Funding Information: We acknowledge all the participants and researchers of the 23 biobanks that have contributed to the GBMI. Biobank-specific acknowledgments are included in the Data S3 . We thank H. Huang, A.R. Martin, B.M. Neale, Y. Okada, K. Tsuo, J.C. Ulirsch, Y. Wang, and all the members of Finucane and Daly labs for their helpful feedback. M.K. was supported by a Nakajima Foundation Fellowship and the Masason Foundation . H.K.F. was funded by NIH grant DP5 OD024582 . Publisher Copyright: © 2022 The Author(s)Meta-analysis is pervasively used to combine multiple genome-wide association studies (GWASs). Fine-mapping of meta-analysis studies is typically performed as in a single-cohort study. Here, we first demonstrate that heterogeneity (e.g., of sample size, phenotyping, imputation) hurts calibration of meta-analysis fine-mapping. We propose a summary statistics-based quality-control (QC) method, suspicious loci analysis of meta-analysis summary statistics (SLALOM), that identifies suspicious loci for meta-analysis fine-mapping by detecting outliers in association statistics. We validate SLALOM in simulations and the GWAS Catalog. Applying SLALOM to 14 meta-analyses from the Global Biobank Meta-analysis Initiative (GBMI), we find that 67% of loci show suspicious patterns that call into question fine-mapping accuracy. These predicted suspicious loci are significantly depleted for having nonsynonymous variants as lead variant (2.7×; Fisher's exact p = 7.3 × 10−4). We find limited evidence of fine-mapping improvement in the GBMI meta-analyses compared with individual biobanks. We urge extreme caution when interpreting fine-mapping results from meta-analysis of heterogeneous cohorts.Peer reviewe

Quantifying the Impact of Rare and Ultra-rare Coding Variation across the Phenotypic Spectrum

There is a limited understanding about the impact of rare protein-truncating variants across multiple phenotypes. We explore the impact of this class of variants on 13 quantitative traits and 10 diseases using whole-exome sequencing data from 100,296 individuals. Protein-truncating variants in genes intolerant to this class of mutations increased risk of autism, schizophrenia, bipolar disorder, intellectual disability, and ADHD. In individuals without these disorders, there was an association with shorter height, lower education, increased hospitalization, and reduced age at enrollment. Gene sets implicated from GWASs did not show a significant protein-truncating variants burden beyond what was captured by established Mendelian genes. In conclusion, we provide a thorough investigation of the impact of rare deleterious coding variants on complex traits, suggesting widespread pleiotropic risk.Peer reviewe

A novel variant in SMG9 causes intellectual disability, confirming a role for nonsense-mediated decay components in neurocognitive development

Biallelic loss-of-function variants in the SMG9 gene, encoding a regulatory subunit of the mRNA nonsense-mediated decay (NMD) machinery, are reported to cause heart and brain malformation syndrome. Here we report five patients from three unrelated families with intellectual disability (ID) and a novel pathogenic SMG9 c.551 T > C p.(Val184Ala) homozygous missense variant, identified using exome sequencing. Sanger sequencing confirmed recessive segregation in each family. SMG9 c.551T > C p.(Val184Ala) is most likely an autozygous variant identical by descent. Characteristic clinical findings in patients were mild to moderate ID, intention tremor, pyramidal signs, dyspraxia, and ocular manifestations. We used RNA sequencing of patients and age- and sex-matched healthy controls to assess the effect of the variant. RNA sequencing revealed that the SMG9 c.551T > C variant did not affect the splicing or expression level of SMG9 gene products, and allele-specific expression analysis did not provide evidence that the nonsense mRNA-induced NMD was affected. Differential gene expression analysis identified prevalent upregulation of genes in patients, including the genes SMOX, OSBP2, GPX3, and ZNF155. These findings suggest that normal SMG9 function may be involved in transcriptional regulation without affecting nonsense mRNA-induced NMD. In conclusion, we demonstrate that the SMG9 c.551T > C missense variant causes a neurodevelopmental disorder and impacts gene expression. NMD components have roles beyond aberrant mRNA degradation that are crucial for neurocognitive development.Peer reviewe