Replication of LDL SWAs hits in PROSPER/PHASE as validation for future (pharmaco)genetic analyses

<p><b>Background:</b> The PHArmacogenetic study of Statins in the Elderly at risk (PHASE) is a genome wide association study in the PROspective Study of Pravastatin in the Elderly at risk for vascular disease (PROSPER) that investigates the genetic variation responsible for the individual variation in drug response to pravastatin. Statins lower LDL-cholesterol in general by 30%, however not in all subjects. Moreover, clinical response is highly variable and adverse effects occur in a minority of patients. In this report we first describe the rationale of the PROSPER/PHASE project and second show that the PROSPER/PHASE study can be used to study pharmacogenetics in the elderly.</p> <p><b>Methods:</b> The genome wide association study (GWAS) was conducted using the Illumina 660K-Quad beadchips following manufacturer's instructions. After a stringent quality control 557,192 SNPs in 5,244 subjects were available for analysis. To maximize the availability of genetic data and coverage of the genome, imputation up to 2.5 million autosomal CEPH HapMap SNPs was performed with MACH imputation software. The GWAS for LDL-cholesterol is assessed with an additive linear regression model in PROBABEL software, adjusted for age, sex, and country of origin to account for population stratification.</p> <p><b>Results:</b> Forty-two SNPs reached the GWAS significant threshold of p = 5.0e-08 in 5 genomic loci (APOE/APOC1; LDLR; FADS2/FEN1; HMGCR; PSRC1/CELSR5). The top SNP (rs445925, chromosome 19) with a p-value of p = 2.8e-30 is located within the APOC1 gene and near the APOE gene. The second top SNP (rs6511720, chromosome 19) with a p-value of p = 5.22e-15 is located within the LDLR gene. All 5 genomic loci were previously associated with LDL-cholesterol levels, no novel loci were identified. Replication in WOSCOPS and CARE confirmed our results.</p> <p><b>Conclusion:</b> With the GWAS in the PROSPER/PHASE study we confirm the previously found genetic associations with LDL-cholesterol levels. With this proof-of-principle study we show that the PROSPER/PHASE study can be used to investigate genetic associations in a similar way to population based studies. The next step of the PROSPER/PHASE study is to identify the genetic variation responsible for the variation in LDL-cholesterol lowering in response to statin treatment in collaboration with other large trials.</p&gt

Novel genetic analysis for case-control genome-wide association studies: quantification of power and genomic prediction accuracy

Genome-wide association studies (GWAS) are routinely conducted for both quantitative and binary (disease) traits. We present two analytical tools for use in the experimental design of GWAS. Firstly, we present power calculations quantifying power in a unified framework for a range of scenarios. In this context we consider the utility of quantitative scores (e.g. endophenotypes) that may be available on cases only or both cases and controls. Secondly, we consider, the accuracy of prediction of genetic risk from genome-wide SNPs and derive an expression for genomic prediction accuracy using a liability threshold model for disease traits in a case-control design. The expected values based on our derived equations for both power and prediction accuracy agree well with observed estimates from simulations

Fifteen new risk loci for coronary artery disease highlight arterial-wall-specific mechanisms

Coronary artery disease (CAD) is a leading cause of morbidity and mortality worldwide. Although 58 genomic regions have been associated with CAD thus far, most of the heritability is unexplained, indicating that additional susceptibility loci await identification. An efficient discovery strategy may be larger-scale evaluation of promising associations suggested by genome-wide association studies (GWAS). Hence, we genotyped 56,309 participants using a targeted gene array derived from earlier GWAS results and performed meta-analysis of results with 194,427 participants previously genotyped, totaling 88,192 CAD cases and 162,544 controls. We identified 25 new SNP-CAD associations (P < 5 × 10(-8), in fixed-effects meta-analysis) from 15 genomic regions, including SNPs in or near genes involved in cellular adhesion, leukocyte migration and atherosclerosis (PECAM1, rs1867624), coagulation and inflammation (PROCR, rs867186 (p.Ser219Gly)) and vascular smooth muscle cell differentiation (LMOD1, rs2820315). Correlation of these regions with cell-type-specific gene expression and plasma protein levels sheds light on potential disease mechanisms

Proinflammatory genotype is associated with the frailty phenotype in the English Longitudinal Study of Ageing

Background: Frailty is a state of increased vulnerability to poor resolution of homeostasis after a stressor event, which increases the risk of adverse outcomes including falls, disability and death. The underlying pathophysiological pathways of frailty are not known but the hypothalamic–pituitary–adrenal axis and heightened chronic systemic inflammation appear to be major contributors. Methods: We used the English Longitudinal Study of Ageing dataset of 3160 individuals over the age of 50 and assessed their frailty status according to the Fried-criteria. We selected single nucleotide polymorphisms in genes involved in the steroid hormone or inflammatory pathways and performed linear association analysis using age and sex as covariates. To support the biological plausibility of any genetic associations, we selected biomarker levels for further analyses to act as potential endophenotypes of our chosen genetic loci. Results: The strongest association with frailty was observed in the Tumor Necrosis Factor (TNF) (rs1800629, P = 0.001198, β = 0.0894) and the Protein Tyrosine Phosphatase, Receptor type, J (PTPRJ) (rs1566729, P = 0.001372, β = 0.09397) genes. Rs1800629 was significantly associated with decreased levels of high-density lipoprotein (HDL) (P = 0.00949) and cholesterol levels (P = 0.00315), whereas rs1566729 was associated with increased levels of HDL (P = 0.01943). After correcting for multiple testing none of the associations remained significant. Conclusions: We provide potential evidence for the involvement of a multifunctional proinflammatory cytokine gene (TNF) in the frailty phenotype. The implication of this gene is further supported by association with the endophenotype biomarker results

Challenges in Using Cultured Primary Rodent Hepatocytes or Cell Lines to Study Hepatic HDL Receptor SR-BI Regulation by Its Cytoplasmic Adaptor PDZK1

Background: PDZK1 is a four PDZ-domain containing cytoplasmic protein that binds to a variety of membrane proteins via their C-termini and can influence the abundance, localization and/or function of its target proteins. One of these targets in hepatocytes in vivo is the HDL receptor SR-BI. Normal hepatic expression of SR-BI protein requires PDZK1 - <5% of normal hepatic SR-BI is seen in the livers of PDZK1 knockout mice. Progress has been made in identifying features of PDZK1 required to control hepatic SR-BI in vivo using hepatic expression of wild-type and mutant forms of PDZK1 in wild-type and PDZK1 KO transgenic mice. Such in vivo studies are time consuming and expensive, and cannot readily be used to explore many features of the underlying molecular and cellular mechanisms. Methodology/Principal Findings: Here we have explored the potential to use either primary rodent hepatocytes in culture using 2D collagen gels with newly developed optimized conditions or PDZK1/SR-BI co-transfected cultured cell lines (COS, HEK293) for such studies. SR-BI and PDZK1 protein and mRNA expression levels fell rapidly in primary hepatocyte cultures, indicating this system does not adequately mimic hepatocytes in vivo for analysis of the PDZK1 dependence of SR-BI. Although PDZK1 did alter SR-BI protein expression in the cell lines, its influence was independent of SR-BI’s C-terminus, and thus is not likely to occur via the same mechanism as that which occurs in hepatocytes in vivo. Conclusions/Significance: Caution must be exercised in using primary hepatocytes or cultured cell lines when studying the mechanism underlying the regulation of hepatic SR-BI by PDZK1. It may be possible to use SR-BI and PDZK1 expression as sensitive markers for the in vivo-like state of hepatocytes to further improve primary hepatocyte cell culture conditions.National Institutes of Health (U.S.) (Grant HL052212)National Institutes of Health (U.S.) (Grant HL066105)National Institutes of Health (U.S.) (Grant ES015241)National Institutes of Health (U.S.) (Grant GM068762

Znf202 Affects High Density Lipoprotein Cholesterol Levels and Promotes Hepatosteatosis in Hyperlipidemic Mice

Background: The zinc finger protein Znf202 is a transcriptional suppressor of lipid related genes and has been linked to hypoalphalipoproteinemia. A functional role of Znf202 in lipid metabolism in vivo still remains to be established. Methodology and Principal Findings: We generated mouse Znf202 expression vectors, the functionality of which was established in several in vitro systems. Next, effects of adenoviral znf202 overexpression in vivo were determined in normo- as well as hyperlipidemic mouse models. Znf202 overexpression in mouse hepatoma cells mhAT3F2 resulted in downregulation of members of the Apoe/c1/c2 and Apoa1/c3/a4 gene cluster. The repressive activity of Znf202 was firmly confirmed in an apoE reporter assay and Znf202 responsive elements within the ApoE promoter were identified. Adenoviral Znf202 transfer to Ldlr-/- mice resulted in downregulation of apoe, apoc1, apoa1, and apoc3 within 24 h after gene transfer. Interestingly, key genes in bile flux (abcg5/8 and bsep) and in bile acid synthesis (cyp7a1) were also downregulated. At 5 days post-infection, the expression of the aforementioned genes was normalized, but mice had developed severe hepatosteatosis accompanied by hypercholesterolemia and hypoalphalipoproteinemia. A much milder phenotype was observed in wildtype mice after 5 days of hepatic Znf202 overexpression. Interestingly and similar to Ldl-/- mice, HDL-cholesterol levels in wildtype mice were lowered after hepatic Znf202 overexpression. Conclusion/Significance: Znf202 overexpression in vivo reveals an important role of this transcriptional regulator in liver lipid homeostasis, while firmly establishing the proposed key role in the control of HDL levels

Are C-Reactive Protein Associated Genetic Variants Associated with Serum Levels and Retinal Markers of Microvascular Pathology in Asian Populations from Singapore?

Introduction:C-reactive protein (CRP) levels are associated with cardiovascular disease and systemic inflammation. We assessed whether CRP-associated loci were associated with serum CRP and retinal markers of microvascular disease, in Asian populations.Methods:Genome-wide association analysis (GWAS) for serum CRP was performed in East-Asian Chinese (N = 2,434) and Malays (N = 2,542) and South-Asian Indians (N = 2,538) from Singapore. Leveraging on GWAS data, we assessed, in silico, association levels among the Singaporean datasets for 22 recently identified CRP-associated loci. At loci where directional inconsistencies were observed, quantification of inter-ethnic linkage disequilibrium (LD) difference was determined. Next, we assessed association for a variant at CRP and retinal vessel traits [central retinal artery equivalent (CRAE) and central retinal vein equivalent (CRVE)] in a total of 24,132 subjects of East-Asian, South-Asian and European ancestry.Results:Serum CRP was associated with SNPs in/near APOE, CRP, HNF1A and LEPR (p-values ≤4.7×10-8) after meta-analysis of Singaporean populations. Using a candidate-SNP approach, we further replicated SNPs at 4 additional loci that had been recently identified to be associated with serum CRP (IL6R, GCKR, IL6 and IL1F10) (p-values ≤0.009), in the Singaporean datasets. SNPs from these 8 loci explained 4.05% of variance in serum CRP. Two SNPs (rs2847281 and rs6901250) were detected to be significant (p-value ≤0.036) but with opposite effect directions in the Singaporean populations as compared to original European studies. At these loci we did not detect significant inter-population LD differences. We further did not observe a significant association between CRP variant and CRVE or CRAE levels after meta-analysis of all Singaporean and European datasets (p-value >0.058).Conclusions:Common variants associated with serum CRP, first detected in primarily European studies, are also associated with CRP levels in East-Asian and South-Asian populations. We did not find a causal link between CRP and retinal measures of microvascular disease

Population genomics of cardiometabolic traits: design of the University College London-London School of Hygiene and Tropical Medicine-Edinburgh-Bristol (UCLEB) Consortium.

Substantial advances have been made in identifying common genetic variants influencing cardiometabolic traits and disease outcomes through genome wide association studies. Nevertheless, gaps in knowledge remain and new questions have arisen regarding the population relevance, mechanisms, and applications for healthcare. Using a new high-resolution custom single nucleotide polymorphism (SNP) array (Metabochip) incorporating dense coverage of genomic regions linked to cardiometabolic disease, the University College-London School-Edinburgh-Bristol (UCLEB) consortium of highly-phenotyped population-based prospective studies, aims to: (1) fine map functionally relevant SNPs; (2) precisely estimate individual absolute and population attributable risks based on individual SNPs and their combination; (3) investigate mechanisms leading to altered risk factor profiles and CVD events; and (4) use Mendelian randomisation to undertake studies of the causal role in CVD of a range of cardiovascular biomarkers to inform public health policy and help develop new preventative therapies

Quantifying Missing Heritability at Known GWAS Loci

Recent work has shown that much of the missing heritability of complex traits can be resolved by estimates of heritability explained by all genotyped SNPs. However, it is currently unknown how much heritability is missing due to poor tagging or additional causal variants at known GWAS loci. Here, we use variance components to quantify the heritability explained by all SNPs at known GWAS loci in nine diseases from WTCCC1 and WTCCC2. After accounting for expectation, we observed all SNPs at known GWAS loci to explain 1.29 X more heritability than GWAS-associated SNPs on average (P = 3.3 X 10[superscript -5]). For some diseases, this increase was individually significant:2.07 X for Multiple Sclerosis (MS) (P = 6.5 X 10 [superscript -9]) and for Crohn's Disease (CD) (P = 1.3 X 10[superscript -3]); all analyses of autoimmune diseases excluded the well-studied MHC region. Additionally, we found that GWAS loci from other related traits also explained significant heritability. The union of all autoimmune disease loci explained 7.15 X more MS heritability than known MS SNPs (P 20,000 Rheumatoid Arthritis (RA) samples typed on ImmunoChip, with 2.37 X more heritability from all SNPs at GWAS loci (P = 2.3 X 10[superscript -6]) and more heritability from all autoimmune disease loci (P < 1 X 10[superscript -16]) compared to known RA SNPs (including those identified in this cohort). Our methods adjust for LD between SNPs, which can bias standard estimates of heritability from SNPs even if all causal variants are typed. By comparing adjusted estimates, we hypothesize that the genome-wide distribution of causal variants is enriched for low-frequency alleles, but that causal variants at known GWAS loci are skewed towards common alleles. These findings have important ramifications for fine-mapping study design and our understanding of complex disease architecture.National Institutes of Health (U.S.) (Grant R03HG006731)National Institutes of Health (U.S.) (Fellowship F32GM106584

Using LASSO regression to detect predictive aggregate effects in genetic studies

We use least absolute shrinkage and selection operator (LASSO) regression to select genetic markers and phenotypic features that are most informative with respect to a trait of interest. We compare several strategies for applying LASSO methods in risk prediction models, using the Genetic Analysis Workshop 17 exome simulation data consisting of 697 individuals with information on genotypic and phenotypic features (smoking, age, sex) in 5-fold cross-validated fashion. The cross-validated averages of the area under the receiver operating curve range from 0.45 to 0.63 for different strategies using only genotypic markers. The same values are improved to 0.69–0.87 when both genotypic and phenotypic information are used. The ability of the LASSO method to find true causal markers is limited, but the method was able to discover several common variants (e.g., FLT1) under certain conditions