Identification of missing variants by combining multiple analytic pipelines

Calibrated rare variant genetic risk scores for complex disease prediction using large exome sequence repositories

Rare variants are collectively numerous and may underlie a considerable proportion of complex disease risk. However, identifying genuine rare variant associations is challenging due to small effect sizes, presence of technical artefacts, and heterogeneity in population structure. We hypothesize that rare variant burden over a large number of genes can be combined into a predictive rare variant genetic risk score (RVGRS). We propose a method (RV-EXCALIBER) that leverages summary-level data from a large public exome sequencing database (gnomAD) as controls and robustly calibrates rare variant burden to account for the aforementioned biases. A calibrated RVGRS strongly associates with coronary artery disease (CAD) in European and South Asian populations by capturing the aggregate effect of rare variants through a polygenic model of inheritance. The RVGRS identifies 1.5% of the population with substantial risk of early CAD and confers risk even when adjusting for known Mendelian CAD genes, clinical risk factors, and a common variant genetic risk score.

Whole-Exome Sequencing and C9orf72 Analysis in Primary Progressive Aphasia

Primary progressive aphasia (PPA) is mainly considered a sporadic disease and few studies have systematically analyzed its genetic basis. We here report the analyses of C9orf72 genotyping and whole-exome sequencing data in a consecutive and well-characterized cohort of 50 patients with PPA. We identified three pathogenic GRN variants, one of them unreported, and two cases with C9orf72 expansions. In addition, one likely pathogenic variant was found in the SQSTM1 gene. Overall, we found 12%of patients carrying pathogenic or likely pathogenic variants. These results support the genetic role in the pathophysiology of a proportion of patients with PPA.

Impact of variant-level batch effects on identification of genetic risk factors in large sequencing studies

Genetic studies have shifted to sequencing-based rare variants discovery after decades of success in identifying common disease variants by Genome-Wide Association Studies using Single Nucleotide Polymorphism chips. Sequencing-based studies require large sample sizes for statistical power and therefore often inadvertently introduce batch effects because samples are typically collected, processed, and sequenced at multiple centers. Conventionally, batch effects are first detected and visualized using Principal Components Analysis and then controlled by including batch covariates in the disease association models. For sequencing-based genetic studies, because all variants included in the association analyses have passed sequencing-related quality control measures, this conventional approach treats every variant as equal and ignores the substantial differences still remaining in variant qualities and characteristics such as genotype quality scores, alternative allele fractions (fraction of reads supporting alternative allele at a variant position) and sequencing depths. In the Alzheimer’s Disease Sequencing Project (ADSP) exome dataset of 9,904 cases and controls, we discovered hidden variant-level differences between sample batches of three sequencing centers and two exome capture kits. Although sequencing centers were included as a covariate in our association models, we observed differences at the variant level in genotype quality and alternative allele fraction between samples processed by different exome capture kits that significantly impacted both the confidence of variant detection and the identification of disease-associated variants. Furthermore, we found that a subset of top disease-risk variants came exclusively from samples processed by one exome capture kit that was more effective at capturing the alternative alleles compared to the other kit. Our findings highlight the importance of additional variant-level quality control for large sequencing-based genetic studies. More importantly, we demonstrate that automatically filtering out variants with batch differences may lead to false negatives if the batch discordances come largely from quality differences and if the batch-specific variants have better quality.

Genomic variant-identification methods may alter Mycobacterium tuberculosis transmission inferences

Pathogen genomic data are increasingly used to characterize global and local transmission patterns of important human pathogens and to inform public health interventions. Yet, there is no current consensus on how to measure genomic variation. To test the effect of the variant-identification approach on transmission inferences for Mycobacterium tuberculosis, we conducted an experiment in which five genomic epidemiology groups applied variant-identification pipelines to the same outbreak sequence data. We compared the variants identified by each group in addition to transmission and phylogenetic inferences made with each variant set. To measure the performance of commonly used variant-identification tools, we simulated an outbreak. We compared the performance of three mapping algorithms, five variant callers and two variant filters in recovering true outbreak variants. Finally, we investigated the effect of applying increasingly stringent filters on transmission inferences and phylogenies. We found that variant-calling approaches used by different groups do not recover consistent sets of variants, which can lead to conflicting transmission inferences. Further, performance in recovering true variation varied widely across approaches. While no single variant-identification approach outperforms others in both recovering true genome-wide and outbreak-level variation, variant-identification algorithms calibrated upon real sequence data or that incorporate local reassembly outperform others in recovering true pairwise differences between isolates. The choice of variant filters contributed to extensive differences across pipelines, and applying increasingly stringent filters rapidly eroded the accuracy of transmission inferences and quality of phylogenies reconstructed from outbreak variation. Commonly used approaches to identify M. tuberculosis genomic variation have variable performance, particularly when predicting potential transmission links from pairwise genetic distances. Phylogenetic reconstruction may be improved by less stringent variant filtering. Approaches that improve variant identification in repetitive, hypervariable regions, such as long-read assemblies, may improve transmission inference.

Alzheimer Disease Pathology-Associated Polymorphism in a Complex Variable Number of Tandem Repeat Region Within the MUC6 Gene, Near the AP2A2 Gene

We found evidence of late-onset Alzheimer disease (LOAD)-associated genetic polymorphism within an exon of Mucin 6 (MUC6) and immediately downstream from another gene: Adaptor Related Protein Complex 2 Subunit Alpha 2 (AP2A2). PCR analyses on genomic DNA samples confirmed that the size of the MUC6 variable number tandem repeat (VNTR) region was highly polymorphic. In a cohort of autopsied subjects with quantitative digital pathology data (n = 119), the size of the polymorphic region was associated with the severity of pTau pathology in neocortex. In a separate replication cohort of autopsied subjects (n = 173), more pTau pathology was again observed in subjects with longer VNTR regions (p = 0.031). Unlike MUC6, AP2A2 is highly expressed in human brain. AP2A2 expression was lower in a subset analysis of brain samples from persons with longer versus shorter VNTR regions (p = 0.014 normalizing with AP2B1 expression). Double-label immunofluorescence studies showed that AP2A2 protein often colocalized with neurofibrillary tangles in LOAD but was not colocalized with pTau proteinopathy in progressive supranuclear palsy, or with TDP-43 proteinopathy. In summary, polymorphism in a repeat-rich region near AP2A2 was associated with neocortical pTau proteinopathy (because of the unique repeats, prior genome-wide association studies were probably unable to detect this association), and AP2A2 was often colocalized with neurofibrillary tangles in LOAD.

The MUC6/AP2A2 Locus and Its Relevance to Alzheimer's Disease: A Review

We recently reported evidence of Alzheimer's disease (AD)-linked genetic variation within the mucin 6 (MUC6) gene on chromosome 11p, nearby the adaptor-related protein complex 2 subunit alpha 2 (AP2A2) gene. This locus has interesting features related to human genomics and clinical research. MUC6 gene variants have been reported to potentially influence viral-including herpesvirus-immunity and the gut microbiome. Within the MUC6 gene is a unique variable number of tandem repeat (VNTR) region. We discovered an association between MUC6 VNTR repeat expansion and AD pathologic severity, particularly tau proteinopathy. Here, we review the relevant literature. The AD-linked VNTR polymorphism may also influence AP2A2 gene expression. AP2A2 encodes a polypeptide component of the adaptor protein complex, AP-2, which is involved in clathrin-coated vesicle function and was previously implicated in AD pathogenesis. To provide background information, we describe some key knowledge gaps in AD genetics research. The "missing/hidden heritability problem" of AD is highlighted. Extensive portions of the human genome, including the MUC6 VNTR, have not been thoroughly evaluated due to limitations of existing high-throughput sequencing technology. We present and discuss additional data, along with cautionary considerations, relevant to the hypothesis that MUC6 repeat expansion influences AD pathogenesis.

Rare genetic variants in interleukin-37 link this anti-inflammatory cytokine to the pathogenesis and treatment of gout

OBJECTIVE

Gout is characterised by severe interleukin (IL)-1-mediated joint inflammation induced by monosodium urate crystals. Since IL-37 is a pivotal anti-inflammatory cytokine suppressing the activity of IL-1, we conducted genetic and functional studies aimed at elucidating the role of IL-37 in the pathogenesis and treatment of gout.

METHODS

Variant identification was performed by DNA sequencing of all coding bases of IL37 using molecular inversion probe-based resequencing (discovery cohort: gout n=675, controls n=520) and TaqMan genotyping (validation cohort: gout n=2202, controls n=2295). Predictive modelling of the effects of rare variants on protein structure was followed by in vitro experiments evaluating the impact on protein function. Treatment with recombinant IL-37 was evaluated in vitro and in vivo in a mouse model of gout.

RESULTS

We identified four rare variants in IL37 in six of the discovery gout patients; p.(A144P), p.(G174Dfs*16), p.(C181*) and p.(N182S), whereas none emerged in healthy controls (Fisher's exact p-value=0.043). All variants clustered in the functional domain of IL-37 in exon 5 (p-value=5.71×10-5). Predictive modelling and functional studies confirmed loss of anti-inflammatory functions and we substantiated the therapeutic potential of recombinant IL-37 in the treatment of gouty inflammation. Furthermore, the carrier status of p.(N182S)(rs752113534) was associated with increased risk (OR=1.81, p-value=0.031) of developing gout in hyperuricaemic individuals of Polynesian ancestry.

CONCLUSION

Here, we provide genetic as well as mechanistic evidence for the role of IL-37 in the pathogenesis of gout, and highlight the therapeutic potential of recombinant IL-37 for the treatment of gouty arthritis.