HAPRAP: a haplotype-based iterative method for statistical fine mapping using GWAS summary statistics.

MOTIVATION: Fine mapping is a widely used approach for identifying the causal variant(s) at disease-associated loci. Standard methods (e.g. multiple regression) require individual level genotypes. Recent fine mapping methods using summary-level data require the pairwise correlation coefficients (r(2)) of the variants. However, haplotypes rather than pairwise r(2), are the true biological representation of linkage disequilibrium (LD) among multiple loci. In this paper, we present an empirical iterative method, HAPlotype Regional Association analysis Program (HAPRAP), that enables fine mapping using summary statistics and haplotype information from an individual-level reference panel. RESULTS: Simulations with individual-level genotypes show that the results of HAPRAP and multiple regression are highly consistent. In simulation with summary-level data, we demonstrate that HAPRAP is less sensitive to poor LD estimates. In a parametric simulation using Genetic Investigation of ANthropometric Traits (GIANT) height data, HAPRAP performs well with a small training sample size (N<2000) while other methods become suboptimal. Moreover, HAPRAP's performance is not affected substantially by SNPs with low minor allele frequencies. We applied the method to existing quantitative trait and binary outcome meta-analyses (human height, QTc interval and gallbladder disease); all previous reported association signals were replicated and two additional variants were independently associated with human height. Due to the growing availability of summary level data, the value of HAPRAP is likely to increase markedly for future analyses (e.g. functional prediction and identification of instruments for Mendelian randomization). AVAILABILITY: The HAPRAP package and documentation are available online: http://apps.biocompute.org.uk/haprap

Genetic studies of accelerometer-based sleep measures yield new insights into human sleep behaviour

Sleep is an essential human function but its regulation is poorly understood. Using accelerometer data from 85,670 UK Biobank participants, we perform a genome-wide association study of 8 derived sleep traits representing sleep quality, quantity and timing, and validate our findings in 5,819 individuals. We identify 47 genetic associations at P < 5 × 10 −8 , of which 20 reach a stricter threshold of P < 8 × 10 −10 . These include 26 novel associations with measures of sleep quality and 10 with nocturnal sleep duration. The majority of identified variants associate with a single sleep trait, except for variants previously associated with restless legs syndrome. For sleep duration we identify a missense variant (p.Tyr727Cys) in PDE11A as the likely causal variant. As a group, sleep quality loci are enriched for serotonin processing genes. Although accelerometer-derived measures of sleep are imperfect and may be affected by restless legs syndrome, these findings provide new biological insights into sleep compared to previous efforts based on self-report sleep measures. © 2019, The Author(s).</p

Exome-wide analysis of rare coding variation identifies novel associations with COPD and airflow limitation in MOCS3, IFIT3 and SERPINA12

Background Several regions of the genome have shown to be associated with COPD in genome-wide association studies of common variants. Objective To determine rare and potentially functional single nucleotide polymorphisms (SNPs) associated with the risk of COPD and severity of airflow limitation. Methods 3226 current or former smokers of European ancestry with lung function measures indicative of Global Initiative for Chronic Obstructive Lung Disease (GOLD) 2 COPD or worse were genotyped using an exome array. An analysis of risk of COPD was carried out using ever smoking controls (n=4784). Associations with %predicted FEV1 were tested in cases. We followed-up signals of interest (p<10−5) in independent samples from a subset of the UK Biobank population and also undertook a more powerful discovery study by meta-analysing the exome array data and UK Biobank data for variants represented on both arrays. Results Among the associated variants were two in regions previously unreported for COPD; a low frequency non-synonymous SNP in MOCS3 (rs7269297, pdiscovery=3.08×10−6, preplication=0.019) and a rare SNP in IFIT3, which emerged in the meta-analysis (rs140549288, pmeta=8.56×10−6). In the meta-analysis of % predicted FEV1 in cases, the strongest association was shown for a splice variant in a previously unreported region, SERPINA12 (rs140198372, pmeta=5.72×10−6). We also confirmed previously reported associations with COPD risk at MMP12, HHIP, GPR126 and CHRNA5. No associations in novel regions reached a stringent exome-wide significance threshold (p<3.7×10−7). Conclusions This study identified several associations with the risk of COPD and severity of airflow limitation, including novel regions MOCS3, IFIT3 and SERPINA12, which warrant further study

A gene-centric study of common carotid artery remodelling

BACKGROUND: Expansive remodelling is the process of compensatory arterial enlargement in response to atherosclerotic stimuli. The genetic determinants of this process are poorly characterized. METHODS: Genetic association analyses of inter-adventitial common carotid artery diameter (ICCAD) in the IMPROVE study (n = 3427) using the Illumina 200k Metabochip was performed. Single nucleotide polymorphisms (SNPs) that met array-wide significance were taken forward for analysis in three further studies (n = 5704), and tested for association with Abdominal Aortic Aneurysm (AAA). RESULTS: rs3768445 on Chromosome 1q24.3, in a cluster of protein coding genes (DNM3, PIGC, C1orf105) was associated with larger ICCAD in the IMPROVE study. For each copy of the rare allele carried, ICCAD was on average 0.13 mm greater (95% CI 0.08-0.18 mm, P = 8.2 × 10(-8)). A proxy SNP (rs4916251, R(2) = 0.99) did not, however, show association with ICCAD in three follow-up studies (P for replication = 0.29). There was evidence of interaction between carotid intima-media thickness (CIMT) and rs4916251 on ICCAD in two of the cohorts studies suggesting that it plays a role in the remodelling response to atherosclerosis. In meta-analysis of 5 case-control studies pooling data from 5007 cases and 43,630 controls, rs4916251 was associated with presence of AAA 1.10, 95% CI 1.03-1.17, p = 2.8 × 10(-3), I(2) = 18.8, Q = 0.30). A proxy SNP, rs4916251 was also associated with increased expression of PIGC in aortic tissue, suggesting that this may the mechanism by which this locus affects vascular remodelling. CONCLUSIONS: Common variation at 1q24.3 is associated with expansive vascular remodelling and risk of AAA. These findings support a hypothesis that pathways involved in systemic vascular remodelling play a role in AAA development

A gene-centric study of common carotid artery remodelling

BACKGROUND: Expansive remodelling is the process of compensatory arterial enlargement in response to atherosclerotic stimuli. The genetic determinants of this process are poorly characterized. METHODS: Genetic association analyses of inter-adventitial common carotid artery diameter (ICCAD) in the IMPROVE study (n = 3427) using the Illumina 200k Metabochip was performed. Single nucleotide polymorphisms (SNPs) that met array-wide significance were taken forward for analysis in three further studies (n = 5704), and tested for association with Abdominal Aortic Aneurysm (AAA). RESULTS: rs3768445 on Chromosome 1q24.3, in a cluster of protein coding genes (DNM3, PIGC, C1orf105) was associated with larger ICCAD in the IMPROVE study. For each copy of the rare allele carried, ICCAD was on average 0.13 mm greater (95% CI 0.08-0.18 mm, P = 8.2 × 10(-8)). A proxy SNP (rs4916251, R(2) = 0.99) did not, however, show association with ICCAD in three follow-up studies (P for replication = 0.29). There was evidence of interaction between carotid intima-media thickness (CIMT) and rs4916251 on ICCAD in two of the cohorts studies suggesting that it plays a role in the remodelling response to atherosclerosis. In meta-analysis of 5 case-control studies pooling data from 5007 cases and 43,630 controls, rs4916251 was associated with presence of AAA 1.10, 95% CI 1.03-1.17, p = 2.8 × 10(-3), I(2) = 18.8, Q = 0.30). A proxy SNP, rs4916251 was also associated with increased expression of PIGC in aortic tissue, suggesting that this may the mechanism by which this locus affects vascular remodelling. CONCLUSIONS: Common variation at 1q24.3 is associated with expansive vascular remodelling and risk of AAA. These findings support a hypothesis that pathways involved in systemic vascular remodelling play a role in AAA development

Mendelian randomization of blood lipids for coronary heart disease

AIMS: To investigate the causal role of high-density lipoprotein cholesterol (HDL-C) and triglycerides in coronary heart disease (CHD) using multiple instrumental variables for Mendelian randomization. METHODS AND RESULTS: We developed weighted allele scores based on single nucleotide polymorphisms (SNPs) with established associations with HDL-C, triglycerides, and low-density lipoprotein cholesterol (LDL-C). For each trait, we constructed two scores. The first was unrestricted, including all independent SNPs associated with the lipid trait identified from a prior meta-analysis (threshold P < 2 × 10(-6)); and the second a restricted score, filtered to remove any SNPs also associated with either of the other two lipid traits at P ≤ 0.01. Mendelian randomization meta-analyses were conducted in 17 studies including 62,199 participants and 12,099 CHD events. Both the unrestricted and restricted allele scores for LDL-C (42 and 19 SNPs, respectively) associated with CHD. For HDL-C, the unrestricted allele score (48 SNPs) was associated with CHD (OR: 0.53; 95% CI: 0.40, 0.70), per 1 mmol/L higher HDL-C, but neither the restricted allele score (19 SNPs; OR: 0.91; 95% CI: 0.42, 1.98) nor the unrestricted HDL-C allele score adjusted for triglycerides, LDL-C, or statin use (OR: 0.81; 95% CI: 0.44, 1.46) showed a robust association. For triglycerides, the unrestricted allele score (67 SNPs) and the restricted allele score (27 SNPs) were both associated with CHD (OR: 1.62; 95% CI: 1.24, 2.11 and 1.61; 95% CI: 1.00, 2.59, respectively) per 1-log unit increment. However, the unrestricted triglyceride score adjusted for HDL-C, LDL-C, and statin use gave an OR for CHD of 1.01 (95% CI: 0.59, 1.75). CONCLUSION: The genetic findings support a causal effect of triglycerides on CHD risk, but a causal role for HDL-C, though possible, remains less certain

Mendelian randomization of blood lipids for coronary heart disease

AIMS: To investigate the causal role of high-density lipoprotein cholesterol (HDL-C) and triglycerides in coronary heart disease (CHD) using multiple instrumental variables for Mendelian randomization. METHODS AND RESULTS: We developed weighted allele scores based on single nucleotide polymorphisms (SNPs) with established associations with HDL-C, triglycerides, and low-density lipoprotein cholesterol (LDL-C). For each trait, we constructed two scores. The first was unrestricted, including all independent SNPs associated with the lipid trait identified from a prior meta-analysis (threshold P < 2 × 10(-6)); and the second a restricted score, filtered to remove any SNPs also associated with either of the other two lipid traits at P ≤ 0.01. Mendelian randomization meta-analyses were conducted in 17 studies including 62,199 participants and 12,099 CHD events. Both the unrestricted and restricted allele scores for LDL-C (42 and 19 SNPs, respectively) associated with CHD. For HDL-C, the unrestricted allele score (48 SNPs) was associated with CHD (OR: 0.53; 95% CI: 0.40, 0.70), per 1 mmol/L higher HDL-C, but neither the restricted allele score (19 SNPs; OR: 0.91; 95% CI: 0.42, 1.98) nor the unrestricted HDL-C allele score adjusted for triglycerides, LDL-C, or statin use (OR: 0.81; 95% CI: 0.44, 1.46) showed a robust association. For triglycerides, the unrestricted allele score (67 SNPs) and the restricted allele score (27 SNPs) were both associated with CHD (OR: 1.62; 95% CI: 1.24, 2.11 and 1.61; 95% CI: 1.00, 2.59, respectively) per 1-log unit increment. However, the unrestricted triglyceride score adjusted for HDL-C, LDL-C, and statin use gave an OR for CHD of 1.01 (95% CI: 0.59, 1.75). CONCLUSION: The genetic findings support a causal effect of triglycerides on CHD risk, but a causal role for HDL-C, though possible, remains less certain

Cystatin C and Cardiovascular Disease: A Mendelian Randomization Study.

BACKGROUND: Epidemiological studies show that high circulating cystatin C is associated with risk of cardiovascular disease (CVD), independent of creatinine-based renal function measurements. It is unclear whether this relationship is causal, arises from residual confounding, and/or is a consequence of reverse causation. OBJECTIVES: The aim of this study was to use Mendelian randomization to investigate whether cystatin C is causally related to CVD in the general population. METHODS: We incorporated participant data from 16 prospective cohorts (n = 76,481) with 37,126 measures of cystatin C and added genetic data from 43 studies (n = 252,216) with 63,292 CVD events. We used the common variant rs911119 in CST3 as an instrumental variable to investigate the causal role of cystatin C in CVD, including coronary heart disease, ischemic stroke, and heart failure. RESULTS: Cystatin C concentrations were associated with CVD risk after adjusting for age, sex, and traditional risk factors (relative risk: 1.82 per doubling of cystatin C; 95% confidence interval [CI]: 1.56 to 2.13; p = 2.12 × 10(-14)). The minor allele of rs911119 was associated with decreased serum cystatin C (6.13% per allele; 95% CI: 5.75 to 6.50; p = 5.95 × 10(-211)), explaining 2.8% of the observed variation in cystatin C. Mendelian randomization analysis did not provide evidence for a causal role of cystatin C, with a causal relative risk for CVD of 1.00 per doubling cystatin C (95% CI: 0.82 to 1.22; p = 0.994), which was statistically different from the observational estimate (p = 1.6 × 10(-5)). A causal effect of cystatin C was not detected for any individual component of CVD. CONCLUSIONS: Mendelian randomization analyses did not support a causal role of cystatin C in the etiology of CVD. As such, therapeutics targeted at lowering circulating cystatin C are unlikely to be effective in preventing CVD

Cystatin C and Cardiovascular Disease: A Mendelian Randomization Study.

BACKGROUND: Epidemiological studies show that high circulating cystatin C is associated with risk of cardiovascular disease (CVD), independent of creatinine-based renal function measurements. It is unclear whether this relationship is causal, arises from residual confounding, and/or is a consequence of reverse causation. OBJECTIVES: The aim of this study was to use Mendelian randomization to investigate whether cystatin C is causally related to CVD in the general population. METHODS: We incorporated participant data from 16 prospective cohorts (n = 76,481) with 37,126 measures of cystatin C and added genetic data from 43 studies (n = 252,216) with 63,292 CVD events. We used the common variant rs911119 in CST3 as an instrumental variable to investigate the causal role of cystatin C in CVD, including coronary heart disease, ischemic stroke, and heart failure. RESULTS: Cystatin C concentrations were associated with CVD risk after adjusting for age, sex, and traditional risk factors (relative risk: 1.82 per doubling of cystatin C; 95% confidence interval [CI]: 1.56 to 2.13; p = 2.12 × 10(-14)). The minor allele of rs911119 was associated with decreased serum cystatin C (6.13% per allele; 95% CI: 5.75 to 6.50; p = 5.95 × 10(-211)), explaining 2.8% of the observed variation in cystatin C. Mendelian randomization analysis did not provide evidence for a causal role of cystatin C, with a causal relative risk for CVD of 1.00 per doubling cystatin C (95% CI: 0.82 to 1.22; p = 0.994), which was statistically different from the observational estimate (p = 1.6 × 10(-5)). A causal effect of cystatin C was not detected for any individual component of CVD. CONCLUSIONS: Mendelian randomization analyses did not support a causal role of cystatin C in the etiology of CVD. As such, therapeutics targeted at lowering circulating cystatin C are unlikely to be effective in preventing CVD

Blood pressure loci identified with a gene-centric array

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10−7 study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r2 = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10−7 at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies