Metasubtract: an R‐package to analytically produce leave‐one‐out meta‐analysis GWAS summary statistics

SUMMARY: statistics from a meta-analysis of genome-wide association studies (meta-GWAS) can be used for many follow-up analyses. One valuable application is the creation of polygenic scores. However, if polygenic scores are calculated in a validation cohort that was part of the meta-GWAS consortium, this cohort is not independent and analyses will therefore yield inflated results. The R package 'MetaSubtract' was developed to subtract the results of the validation cohort from meta-GWAS summary statistics analytically. The statistical formulas for a meta-analysis were inverted to compute corrected summary statistics of a meta-GWAS leaving one (or more) cohort(s) out. These formulas have been implemented in MetaSubtract for different meta-analyses methods (fixed effects inverse variance or square root sample size weighted z-score) accounting for no, single or double genomic control correction. Results obtained by MetaSubtract correlate very well to those calculated using the traditional way, i.e. by performing a meta-analysis leaving out the validation cohort. In conclusion, MetaSubtract allows researchers to compute meta-GWAS summary statistics that are independent of the GWAS results of the validation cohort without requiring access to the cohort level GWAS results of the corresponding meta-GWAS consortium. AVAILABILITY AND IMPLEMENTATION: https://cran.r-project.org/web/packages/MetaSubtract. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online

Genetic pre-screening for glaucoma in population-based epidemiology:protocol for a double-blind prospective screening study within Lifelines (EyeLife)

BACKGROUND: Early detection of glaucoma is paramount to maintain patients' eyesight, however glaucomatous vision loss tends to begin in the periphery with up to 50% of patients unaware they are affected. Because glaucomatous vision loss is permanent, screening appears attractive, but currently is not cost-effective. Therefore we aim to investigate the utility of genetic pre-screening for glaucoma in a population-based setting, called EyeLife. METHODS: EyeLife adopts a double blind prospective design with contrasting groups. Selected participants (n = 1600) from the Lifelines cohort are 55 years of age or older, and of either the highest or lowest 20% of the genetic risk distribution for glaucoma. We obtained a highly curated list of genetic variants from the literature to obtain each participants' genetic risk for glaucoma. Participants will undergo comprehensive ophthalmic screening. The primary outcome is the relative risk of glaucoma given a high genetic risk compared to a low genetic risk. DISCUSSION: If genetic pre-screening is successful, it will increase the yield of a glaucoma screening program by focusing on high-risk individuals. This, in turn, may improve long-term visual health of middle-aged and elderly people. TRIAL REGISTRATION: Ethics approval was obtained on January 31, 2019, and the study was retrospectively registered with the Netherlands Trial Register ( NL8718 ) on the 17th of June, 2020

Ethnic differences in prevalence of Dupuytren disease can partly be explained by known genetic risk variants

Dupuytren disease (DD), a fibroproliferative disorder of the palmar fascia that causes flexion contractures in the fingers, is prevalent in people of North-Western European descent and less so in other ethnicities. DD is a complex disorder, influenced by genetic risk variants. We aimed to study if the marked differences in prevalences in DD between ethnic (sub)groups could be explained by differences in allele frequencies of the 26 known genetic risk variants of DD. Therefore, genetic risk scores (GRS) composed of the 26 DD risk variants were calculated for the 26 populations from the 1000 Genomes database and correlated to observed DD prevalences from literature. For comparison, GRSs were generated for 10,000 sets of 26 random SNPs and also correlated to the observed DD prevalences to determine the significance of the observed correlation. To determine whether differences in allele frequencies between ethnicities were caused by natural selection, fixation indices (Fst) were calculated from the 26 SNPs and from the sets of 26 random SNPs for comparison. Observed prevalences could be determined from literature for 10 populations. Their correlation with the GRS composed of DD SNPs proved to be 0.60 (p = 0.0003). The Fsts between British and other populations were low for European, ad mixed American, and South-Asian populations, and moderate for East-Asians. African populations were significantly different from expected values determined from the random sets. In conclusion, the 26 known genetic risk variants associated with DD explain for a substantial part (R-2 = 0.36) the differing DD prevalences observed between ethnicities

Novel genes for QTc interval. How much heritability is explained, and how much is left to find?

The corrected QT (QTc) interval is a complex quantitative trait, believed to be influenced by several genetic and environmental factors. It is a strong prognostic indicator of cardiovascular mortality in patients with and without cardiac disease. More than 700 mutations have been described in 12 genes (LQT1-LQT12) involved in congenital long QT syndrome. However, the heritability (genetic contribution) of QTc interval in the general population cannot be adequately explained by these long QT syndrome genes. In order to further investigate the genetic architecture underlying QTc interval in the general population, genome-wide association studies, in which up to one million single nucleotide polymorphisms are assayed in thousands of individuals, are now being employed and have already led to the discovery of variants in seven novel loci and five loci that are known to cause congenital long or short QT syndrome. Here we show that a combined risk score using 11 of these loci explains about 10% of the heritability of QTc. Additional discovery of both common and rare variants will yield further etiological insight and accelerate clinical applications

The associations of CNR1 SNPs and haplotypes with vulnerability and treatment response phenotypes in Han Chinese with major depressive disorder:A case-control association study

BACKGROUND: Understanding how genetic polymorphisms are associated with the pathophysiology of major depressive disorder (MDD) may aid in diagnosis and the development of personalized treatment strategies. CNR1 is the gene coding Cannabinoid type 1 receptor which is highly involved in emotional processing and in regulating neurotransmitter releases. We aimed to investigate the associations of CNR1 single‐nucleotide polymorphisms (SNPs) with MDD susceptibility and treatment response. METHODS: The study reported data on 181 Han Chinese with MDD and 80 healthy controls. The associations of CNR1 genetic polymorphisms with MDD susceptibility and treatment response were examined, wherein the MDD patients were subgrouped further by responding to antidepressant treatment, compared with healthy controls separately. RESULTS: The CNR1 SNPs rs806367 and rs6454674 and haplotype C‐T‐T‐C of rs806366, rs806367, rs806368, and rs806370 were associated with increased susceptibility for MDD and antidepressant treatment resistance, but the association was not detected in other SNPs or the haplotype block of rs806368 and rs806370. CONCLUSION: The CNR1 is a promising candidate for the genetic association study of MDD. Larger and well‐characterized samples are required to confirm the genetic association of CNR1 with MDD because of the limitations such as relatively small sample size and lack of information for correcting confounding factors

GWASinspector:comprehensive quality control of genome-wide association study results

Quality control (QC) of genome wide association study (GWAS) result files has become increasingly difficult due to advances in genomic technology. The main challenges include continuous increases in the number of polymorphic genetic variants contained in recent GWASs and reference panels, the rising number of cohorts participating in a GWAS consortium, and inclusion of new variant types. Here, we present GWASinspector, a flexible R package for comprehensive QC of GWAS results. This package is compatible with recent imputation reference panels, handles insertion/deletion and multi-allelic variants, provides extensive QC reports and efficiently processes big data files. Reference panels covering three human genome builds (NCBI36, GRCh37 and GRCh38) are available. GWASinspector has a user friendly design and allows easy set-up of the QC pipeline through a configuration file. In addition to checking and reporting on individual files, it can be used in preparation of a meta-analysis by testing for systemic differences between studies and generating cleaned, harmonized GWAS files. Comparison with existing GWAS QC tools shows that the main advantages of GWASinspector are its ability to more effectively deal with insertion/deletion and multi-allelic variants and its relatively low memory use

Mitochondrial Genome Study Identifies Association Between Primary Open-Angle Glaucoma and Variants in MT-CYB, MT-ND4 Genes and Haplogroups

Background and purpose: Primary open-angle glaucoma (POAG) is an optic neuropathy characterized by death of retinal ganglion cells and atrophy of the optic nerve head. The susceptibility of the optic nerve to damage has been shown to be mediated by mitochondrial dysfunction. In this study, we aimed to determine a possible association between mitochondrial SNPs or haplogroups and POAG. Methods: Mitochondrial DNA single nucleotide polymorphisms (mtSNPs) were genotyped using the Illumina Infinium Global Screening Array-24 (GSA) 700K array set. Genetic analyses were performed in a POAG case-control study involving the cohorts, Groningen Longitudinal Glaucoma Study-Lifelines Cohort Study and Amsterdam Glaucoma Study, including 721 patients and 1951 controls in total. We excluded samples not passing quality control for nuclear genotypes and samples with low call rate for mitochondrial variation. The mitochondrial variants were analyzed both as SNPs and haplogroups. These were determined with the bioinformatics software HaploGrep, and logistic regression analysis was used for the association, as well as for SNPs. Results: Meta-analysis of the results from both cohorts revealed a significant association between POAG and the allele A of rs2853496 [odds ratio (OR) = 0.64; p = 0.006] within the MT-ND4 gene, and for the T allele of rs35788393 (OR = 0.75; p = 0.041) located in the MT-CYB gene. In the mitochondrial haplogroup analysis, the most significant p-value was reached by haplogroup K (p = 1.2 × 10(−05)), which increases the risk of POAG with an OR of 5.8 (95% CI 2.7–13.1). Conclusion: We identified an association between POAG and polymorphisms in the mitochondrial genes MT-ND4 (rs2853496) and MT-CYB (rs35788393), and with haplogroup K. The present study provides further evidence that mitochondrial genome variations are implicated in POAG. Further genetic and functional studies are required to substantiate the association between mitochondrial gene polymorphisms and POAG and to define the pathophysiological mechanisms of mitochondrial dysfunction in glaucoma

Killer Cell Immunoglobulin-Like Receptor Haplotype B Modulates Susceptibility to EBV-Associated Classic Hodgkin Lymphoma

Tumor cells of classic Hodgkin lymphoma (cHL) are derived from antigen presenting B cells that are infected by Epstein Barr virus (EBV) in ~30% of patients. Polymorphic Killer cell immunoglobulin-like receptors (KIRs) expressed on NK cells interact with human leukocyte antigen (HLA) class I and play a key role in immune surveillance against virally infected cells and tumor cells. We investigated the effect of KIR types on cHL susceptibility overall (n=211) and in EBV-stratified subgroups using the Dutch GoNL cohort as controls (n=498). The frequency of the KIR haplotype B subgroup was significantly different between EBV+ and EBV− cHL patients (62% vs. 77%, p=0.04) and this difference was more pronounced in nodular sclerosis (NS) cHL (49% vs. 79%, p=0.0003). The frequency of KIR haplotype B subgroup was significantly lower in EBV+ NS cHL compared to controls (49% vs. 67%, p=0.01). Analyses of known KIR – HLA interaction pairs revealed lower carrier frequencies of KIR2DS2 – HLA-C1 (29% vs. 46%, p=0.03) and KIR2DL2 – HLA-C1 (29% vs. 45%, p=0.04) in EBV+ NS cHL patients compared to controls. Carriers of the KIR haplotype B subgroup are less likely to develop EBV+ NS cHL, probably because of a more efficient control over EBV-infected B cells