Bram P. Prins

Poor replication of candidate genes for major depressive disorder using genome-wide association data.

Data from the Genetic Association Information Network (GAIN) genome-wide association study (GWAS) in major depressive disorder (MDD) were used to explore previously reported candidate gene and single-nucleotide polymorphism (SNP) associations in MDD. A systematic literature search of candidate genes associated with MDD in case–control studies was performed before the results of the GAIN MDD study became available. Measured and imputed candidate SNPs and genes were tested in the GAIN MDD study encompassing 1738 cases and 1802 controls. Imputation was used to increase the number of SNPs from the GWAS and to improve coverage of SNPs in the candidate genes selected. Tests were carried out for individual SNPs and the entire gene using different statistical approaches, with permutation analysis as the final arbiter. In all, 78 papers reporting on 57 genes were identified, from which 92 SNPs could be mapped. In the GAIN MDD study, two SNPs were associated with MDD: C5orf20 (rs12520799; P=0.038; odds ratio (OR) AT=1.10, 95% CI 0.95–1.29; OR TT=1.21, 95% confidence interval (CI) 1.01–1.47) and NPY (rs16139; P=0.034; OR C allele=0.73, 95% CI 0.55–0.97), constituting a direct replication of previously identified SNPs. At the gene level, TNF (rs76917; OR T=1.35, 95% CI 1.13–1.63; P=0.0034) was identified as the only gene for which the association with MDD remained significant after correction for multiple testing. For SLC6A2 (norepinephrine transporter (NET)) significantly more SNPs (19 out of 100; P=0.039) than expected were associated while accounting for the linkage disequilibrium (LD) structure. Thus, we found support for involvement in MDD for only four genes. However, given the number of candidate SNPs and genes that were tested, even these significant may well be false positives. The poor replication may point to publication bias and false-positive findings in previous candidate gene studies, and may also be related to heterogeneity of the MDD phenotype as well as contextual genetic or environmental factors.

Role of common and rare variants in SCN10A: results from the Brugada syndrome QRS locus gene discovery collaborative study

AIMS Brugada syndrome (BrS) remains genetically heterogeneous and is associated with slowed cardiac conduction. We aimed to identify genetic variation in BrS cases at loci associated with QRS duration. METHODS AND RESULTS A multi-centre study sequenced seven candidate genes (SCN10A, HAND1, PLN, CASQ2, TKT, TBX3, and TBX5) in 156 Caucasian SCN5A mutation-negative BrS patients (80% male; mean age 48) with symptoms (64%) and/or a family history of sudden death (47%) or BrS (18%). Forty-nine variants were identified: 18 were rare (MAF <1%) and non-synonymous; and 11/18 (61.1%), mostly in SCN10A, were predicted as pathogenic using multiple bioinformatics tools. Allele frequencies were compared with the Exome Sequencing and UK10K Projects. SKAT methods tested rare variation in SCN10A finding no statistically significant difference between cases and controls. Co-segregation analysis was possible for four of seven probands carrying a novel pathogenic variant. Only one pedigree (I671V/G1299A in SCN10A) showed co-segregation. The SCN10A SNP V1073 was, however, associated strongly with BrS [66.9 vs. 40.1% (UK10K) OR (95% CI) = 3.02 (2.35-3.87), P = 8.07 × 10-19]. Voltage-clamp experiments for NaV1.8 were performed for SCN10A common variants V1073, A1073, and rare variants of interest: A200V and I671V. V1073, A200V and I671V, demonstrated significant reductions in peak INa compared with ancestral allele A1073 (rs6795970). CONCLUSION Rare variants in the screened QRS-associated genes (including SCN10A) are not responsible for a significant proportion of SCN5A mutation negative BrS. The common SNP SCN10A V1073 was strongly associated with BrS and demonstrated loss of NaV1.8 function, as did rare variants in isolated patients.

Genetics of coronary artery disease: Genome-wide association studies and beyond

Genome-wide association (GWA) studies on coronary artery disease (CAD) have been very successful, identifying a total of 32 susceptibility loci so far. Although these loci have provided valuable insights into the etiology of CAD, their cumulative effect explains surprisingly little of the total CAD heritability. In this review, we first highlight and describe the type of genetic variants potentially underlying the missing heritability of CAD: single nucleotide polymorphisms (SNPs) or structural variants, each of which may either be common or rare. Although finding missing heritability is important, we further argue in this review that it constitutes only a first step towards a fuller understanding of the etiology of CAD development. To close the gap between the genotype and phenotype, we propose a systems genetics approach in the post-GWA study era. This approach that integrates genetic, epigenetic, transcriptomic, proteomic, metabolic and intermediate outcome variables has potential to significantly aid the understanding of CAD etiology.

Investigating the Causal Relationship of C-Reactive Protein with 32 Complex Somatic and Psychiatric Outcomes: A Large-Scale Cross-Consortium Mendelian Randomization Study

Background C-reactive protein (CRP) is associated with immune, cardiometabolic, and psychiatric traits and diseases. Yet it is inconclusive whether these associations are causal. Methods and Findings We performed Mendelian randomization (MR) analyses using two genetic risk scores (GRSs) as instrumental variables (IVs). The first GRS consisted of four single nucleotide polymorphisms (SNPs) in the CRP gene (GRSCRP), and the second consisted of 18 SNPs that were significantly associated with CRP levels in the largest genome-wide association study (GWAS) to date (GRSGWAS). To optimize power, we used summary statistics from GWAS consortia and tested the association of these two GRSs with 32 complex somatic and psychiatric outcomes, with up to 123,865 participants per outcome from populations of European ancestry. We performed heterogeneity tests to disentangle the pleiotropic effect of IVs. A Bonferroni-corrected significance level of less than 0.0016 was considered statistically significant. An observed p-value equal to or less than 0.05 was considered nominally significant evidence for a potential causal association, yet to be confirmed. The strengths (F-statistics) of the IVs were 31.92–3,761.29 and 82.32–9,403.21 for GRSCRP and GRSGWAS, respectively. CRP GRSGWAS showed a statistically significant protective relationship of a 10% genetically elevated CRP level with the risk of schizophrenia (odds ratio [OR] 0.86 [95% CI 0.79–0.94]; p < 0.001). We validated this finding with individual-level genotype data from the schizophrenia GWAS (OR 0.96 [95% CI 0.94–0.98]; p < 1.72 × 10−6). Further, we found that a standardized CRP polygenic risk score (CRPPRS) at p-value thresholds of 1 × 10−4, 0.001, 0.01, 0.05, and 0.1 using individual-level data also showed a protective effect (OR < 1.00) against schizophrenia; the first CRPPRS (built of SNPs with p < 1 × 10−4) showed a statistically significant (p < 2.45 × 10−4) protective effect with an OR of 0.97 (95% CI 0.95–0.99). The CRP GRSGWAS showed that a 10% increase in genetically determined CRP level was significantly associated with coronary artery disease (OR 0.88 [95% CI 0.84–0.94]; p < 2.4 × 10−5) and was nominally associated with the risk of inflammatory bowel disease (OR 0.85 [95% CI 0.74–0.98]; p < 0.03), Crohn disease (OR 0.81 [95% CI 0.70–0.94]; p < 0.005), psoriatic arthritis (OR 1.36 [95% CI 1.00–1.84]; p < 0.049), knee osteoarthritis (OR 1.17 [95% CI 1.01–1.36]; p < 0.04), and bipolar disorder (OR 1.21 [95% CI 1.05–1.40]; p < 0.007) and with an increase of 0.72 (95% CI 0.11–1.34; p < 0.02) mm Hg in systolic blood pressure, 0.45 (95% CI 0.06–0.84; p < 0.02) mm Hg in diastolic blood pressure, 0.01 ml/min/1.73 m2 (95% CI 0.003–0.02; p < 0.005) in estimated glomerular filtration rate from serum creatinine, 0.01 g/dl (95% CI 0.0004–0.02; p < 0.04) in serum albumin level, and 0.03 g/dl (95% CI 0.008–0.05; p < 0.009) in serum protein level. However, after adjustment for heterogeneity, neither GRS showed a significant effect of CRP level (at p < 0.0016) on any of these outcomes, including coronary artery disease, nor on the other 20 complex outcomes studied. Our study has two potential limitations: the limited variance explained by our genetic instruments modeling CRP levels in blood and the unobserved bias introduced by the use of summary statistics in our MR analyses. Conclusions Genetically elevated CRP levels showed a significant potentially protective causal relationship with risk of schizophrenia. We observed nominal evidence at an observed p < 0.05 using either GRSCRP or GRSGWAS—with persistence after correction for heterogeneity—for a causal relationship of elevated CRP levels with psoriatic osteoarthritis, rheumatoid arthritis, knee osteoarthritis, systolic blood pressure, diastolic blood pressure, serum albumin, and bipolar disorder. These associations remain yet to be confirmed. We cannot verify any causal effect of CRP level on any of the other common somatic and neuropsychiatric outcomes investigated in the present study. This implies that interventions that lower CRP level are unlikely to result in decreased risk for the majority of common complex outcomes.

Beyond Genome-Wide Association Studies: New Strategies for Identifying Genetic Determinants of Hypertension

Genetic linkage and association methods have long been the most important tools for gene identification in humans. These approaches can either be hypothesis-based (i.e., candidate-gene studies) or hypothesis-free (i.e., genome-wide studies). The first part of this review offers an overview of the latest successes in gene finding for blood pressure (BP) and essential hypertension using these DNA sequence–based discovery techniques. We further emphasize the importance of post–genome-wide association study (post-GWAS) analysis, which aims to prioritize genetic variants for functional follow-up. Whole-genome next-generation sequencing will eventually be necessary to provide a more comprehensive picture of all DNA variants affecting BP and hypertension. The second part of this review discusses promising novel approaches that move beyond the DNA sequence and aim to discover BP genes that are differentially regulated by epigenetic mechanisms, including microRNAs, histone modification, and methylation.

Eosinophil Count Is a Common Factor for Complex Metabolic and Pulmonary Traits and Diseases: The LifeLines Cohort Study

There is ongoing debate on the association between eosinophil count and diseases, as previous studies were inconsistent. We studied the relationship of eosinophil count with 22 complex metabolic, cardiac, and pulmonary traits and diseases. From the population-based LifeLines Cohort Study (N = 167,729), 13,301 individuals were included. We focused on relationship of eosinophil count with three classes of metabolic (7 traits, 2 diseases), cardiac (6 traits, 2 diseases), and pulmonary (2 traits, 2 diseases) outcomes. Regression analyses were applied in overall, women and men, while adjusted for age, sex, BMI and smoking. A p-value of <0.00076 was considered statistically significant. 58.2% of population were women (mean±SD 51.3±11.1 years old). In overall, one-SD higher of ln-eosinophil count was associated with a 0.04 (±SE ±0.002;p = 6.0×10−6) SD higher levels in ln-BMI, 0.06 (±0.007;p = 3.1×10−12) SD in ln-TG, 0.04 (±0.003;p = 7.0×10−6) SD in TC, 0.04 (±0.004;p = 6.3×10−7) SD in LDL, 0.04 (±0.006;p = 6.0×10−6) SD in HbA1c; and with a 0.05 (±0.004;p = 1.7×10−8) SD lower levels in HDL, 0.05 (±0.007;p = 3.4×10−23) SD in FEV1, and 0.09 (±0.001;p = 6.6×10−28) SD in FEV1/FVC. A higher ln-eosinophil count was associated with 1.18 (95%CI 1.09–1.28;p = 2.0×10−5) odds ratio of obesity, 1.29 (1.19–1.39;p = 1.1×10−10) of metabolic syndrome, 1.40 (1.25–1.56;p = 2.7×10−9) of COPD and 1.81 (1.61–2.03;p = 1.0×10−23) of asthma. Similar results were found in women. We found no association between ln-eosinophil count either with blood pressure indices in overall, women and men; or with BMI, LDL, HbA1c and obesity in men. In a large population based cohort, we confirmed eosinophil count as a potential factor implicated in metabolic and pulmonary outcomes.

In Silico Post Genome-Wide Association Studies Analysis of C-Reactive Protein Loci Suggests an Important Role for Interferons

Background—Genome-wide association studies (GWASs) have successfully identified several single nucleotide polymorphisms (SNPs) associated with serum levels of C-reactive protein (CRP). An important limitation of GWASs is that the identified variants merely flag the nearby genomic region and do not necessarily provide a direct link to the biological mechanisms underlying their corresponding phenotype. Here we apply a bioinformatics-based approach to uncover the functional characteristics of the 18 SNPs that had previously been associated with CRP at a genome-wide significant level. Methods and Results—In the first phase of in silico sequencing, we explore the vicinity of GWAS SNPs to identify all linked variants. In the second phase of expression quantitative trait loci analysis, we attempt to identify all nearby genes whose expression levels are associated with the corresponding GWAS SNPs. These 2 phases generate several relevant genes that serve as input to the next phase of functional network analysis. Our in silico sequencing analysis using 1000 Genomes Project data identified 7 nonsynonymous SNPs, which are in moderate to high linkage disequilibrium (r 2>0.5) with the GWAS SNPs. Our expression quantitative trait loci analysis, which was based on one of the largest single data sets of genome-wide expression probes (n>5000) identified 23 significantly associated expression probes belonging to 15 genes (false discovery rate <0.01). The final phase of functional network analysis revealed 93 significantly enriched biological processes (false discovery rate <0.01). Conclusions—Our post-GWAS analysis of CRP GWAS SNPs confirmed the previously known overlap between CRP and lipids biology. Additionally, it suggested an important role for interferons in the metabolism of CRP.

Sequencing of SCN5A Identifies Rare and Common Variants Associated With Cardiac Conduction: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium

Background—The cardiac sodium channel SCN5A regulates atrioventricular and ventricular conduction. Genetic variants in this gene are associated with PR and QRS intervals. We sought to characterize further the contribution of rare and common coding variation in SCN5A to cardiac conduction. Methods and Results—In Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study, we performed targeted exonic sequencing of SCN5A (n=3699, European ancestry individuals) and identified 4 common (minor allele frequency >1%) and 157 rare variants. Common and rare SCN5A coding variants were examined for association with PR and QRS intervals through meta-analysis of European ancestry participants from CHARGE, National Heart, Lung, and Blood Institute’s Exome Sequencing Project (n=607), and the UK10K (n=1275) and by examining Exome Sequencing Project African ancestry participants (n=972). Rare coding SCN5A variants in aggregate were associated with PR interval in European and African ancestry participants (P=1.3×10−3). Three common variants were associated with PR and QRS interval duration among European ancestry participants and one among African ancestry participants. These included 2 well-known missense variants: rs1805124 (H558R) was associated with PR and QRS shortening in European ancestry participants (P=6.25×10−4 and P=5.2×10−3, respectively) and rs7626962 (S1102Y) was associated with PR shortening in those of African ancestry (P=2.82×10−3). Among European ancestry participants, 2 novel synonymous variants, rs1805126 and rs6599230, were associated with cardiac conduction. Our top signal, rs1805126 was associated with PR and QRS lengthening (P=3.35×10−7 and P=2.69×10−4, respectively) and rs6599230 was associated with PR shortening (P=2.67×10−5). Conclusions—By sequencing SCN5A, we identified novel common and rare coding variants associated with cardiac conduction.

An in silico Post-GWAS Analysis of C-Reactive Protein Loci Suggests an Important Role for Interferons

Background—Genome-wide association studies (GWASs) have successfully identified several single nucleotide polymorphisms (SNPs) associated with serum levels of C-reactive protein (CRP). An important limitation of GWASs is that the identified variants merely flag the nearby genomic region and do not necessarily provide a direct link to the biological mechanisms underlying their corresponding phenotype. Here we apply a bioinformatics-based approach to uncover the functional characteristics of the 18 SNPs that had previously been associated with CRP at a genome-wide significant level. Methods and Results—In the first phase of in silico sequencing, we explore the vicinity of GWAS SNPs to identify all linked variants. In the second phase of expression quantitative trait loci analysis, we attempt to identify all nearby genes whose expression levels are associated with the corresponding GWAS SNPs. These 2 phases generate several relevant genes that serve as input to the next phase of functional network analysis. Our in silico sequencing analysis using 1000 Genomes Project data identified 7 nonsynonymous SNPs, which are in moderate to high linkage disequilibrium (r 2>0.5) with the GWAS SNPs. Our expression quantitative trait loci analysis, which was based on one of the largest single data sets of genome-wide expression probes (n>5000) identified 23 significantly associated expression probes belonging to 15 genes (false discovery rate <0.01). The final phase of functional network analysis revealed 93 significantly enriched biological processes (false discovery rate <0.01). Conclusions—Our post-GWAS analysis of CRP GWAS SNPs confirmed the previously known overlap between CRP and lipids biology. Additionally, it suggested an important role for interferons in the metabolism of CRP.

Evidence for large-scale gene-by-smoking interaction effects on pulmonary function

Abstract Background: Smoking is the strongest environmental risk factor for reduced pulmonary function. The genetic component of various pulmonary traits has also been demonstrated, and at least 26 loci have been reproducibly associated with either FEV1 (forced expiratory volume in 1 second) or FEV1/FVC (FEV1/forced vital capacity). Although the main effects of smoking and genetic loci are well established, the question of potential gene-by-smoking interaction effect remains unanswered. The aim of the present study was to assess, using a genetic risk score approach, whether the effect of these 26 loci on pulmonary function is influenced by smoking. Methods: We evaluated the interaction between smoking exposure, considered as either ever vs never or pack-years, and a 26-single nucleotide polymorphisms (SNPs) genetic risk score in relation to FEV1 or FEV1/FVC in 50 047 participants of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) and SpiroMeta consortia. Results: We identified an interaction (βint = –0.036, 95% confidence interval, –0.040 to –0.032, P = 0.00057) between an unweighted 26 SNP genetic risk score and smoking status (ever/never) on the FEV1/FVC ratio. In interpreting this interaction, we showed that the genetic risk of falling below the FEV1/FVC threshold used to diagnose chronic obstructive pulmonary disease is higher among ever smokers than among never smokers. A replication analysis in two independent datasets, although not statistically significant, showed a similar trend in the interaction effect. Conclusions: This study highlights the benefit of using genetic risk scores for identifying interactions missed when studying individual SNPs and shows, for the first time, that persons with the highest genetic risk for low FEV1/FVC may be more susceptible to the deleterious effects of smoking.