Synthesis-View: visualization and interpretation of SNP association results for multi-cohort, multi-phenotype data and meta-analysis

<p>Abstract</p> <p>Background</p> <p>Initial genome-wide association study (GWAS) discoveries are being further explored through the use of large cohorts across multiple and diverse populations involving meta-analyses within large consortia and networks. Many of the additional studies characterize less than 100 single nucleotide polymorphisms (SNPs), often include multiple and correlated phenotypic measurements, and can include data from multiple-sites, multiple-studies, as well as multiple race/ethnicities. New approaches for visualizing resultant data are necessary in order to fully interpret results and obtain a broad view of the trends between DNA variation and phenotypes, as well as provide information on specific SNP and phenotype relationships.</p> <p>Results</p> <p>The Synthesis-View software tool was designed to visually synthesize the results of the aforementioned types of studies. Presented herein are multiple examples of the ways Synthesis-View can be used to report results from association studies of DNA variation and phenotypes, including the visual integration of p-values or other metrics of significance, allele frequencies, sample sizes, effect size, and direction of effect.</p> <p>Conclusions</p> <p>To truly allow a user to visually integrate multiple pieces of information typical of a genetic association study, innovative views are needed to integrate multiple pieces of information. As a result, we have created "Synthesis-View" software for the visualization of genotype-phenotype association data in multiple cohorts. Synthesis-View is freely available for non-commercial research institutions, for full details see <url>https://chgr.mc.vanderbilt.edu/synthesisview</url>.</p

The promoter polymorphism -232C/G of the PCK1 gene is associated with type 2 diabetes in a UK-resident South Asian population

Background: The PCK1 gene, encoding cytosolic phosphoenolpyruvate carboxykinase (PEPCK-C), has previously been implicated as a candidate gene for type 2 diabetes (T2D) susceptibility. Rodent models demonstrate that over-expression of Pck1 can result in T2D development and a single nucleotide polymorphism (SNP) in the promoter region of human PCK1 (-232C/G) has exhibited significant association with the disease in several cohorts. Within the UK-resident South Asian population, T2D is 4 to 6 times more common than in indigenous white Caucasians. Despite this, few studies have reported on the genetic susceptibility to T2D in this ethnic group and none of these has investigated the possible effect of PCK1 variants. We therefore aimed to investigate the association between common variants of the PCK1 gene and T2D in a UK-resident South Asian population of Punjabi ancestry, originating predominantly from the Mirpur area of Azad Kashmir, Pakistan. \ud \ud Methods: We used TaqMan assays to genotype five tagSNPs covering the PCK1 gene, including the -232C/G variant, in 903 subjects with T2D and 471 normoglycaemic controls. \ud \ud Results: Of the variants studied, only the minor allele (G) of the -232C/G SNP demonstrated a significant association with T2D, displaying an OR of 1.21 (95% CI: 1.03 - 1.42, p = 0.019). \ud \ud Conclusion: This study is the first to investigate the association between variants of the PCK1 gene and T2D in South Asians. Our results suggest that the -232C/G promoter polymorphism confers susceptibility to T2D in this ethnic group. \ud \ud Trial registration: UKADS Trial Registration: ISRCTN38297969

Calibrating the Performance of SNP Arrays for Whole-Genome Association Studies

To facilitate whole-genome association studies (WGAS), several high-density SNP genotyping arrays have been developed. Genetic coverage and statistical power are the primary benchmark metrics in evaluating the performance of SNP arrays. Ideally, such evaluations would be done on a SNP set and a cohort of individuals that are both independently sampled from the original SNPs and individuals used in developing the arrays. Without utilization of an independent test set, previous estimates of genetic coverage and statistical power may be subject to an overfitting bias. Additionally, the SNP arrays' statistical power in WGAS has not been systematically assessed on real traits. One robust setting for doing so is to evaluate statistical power on thousands of traits measured from a single set of individuals. In this study, 359 newly sampled Americans of European descent were genotyped using both Affymetrix 500K (Affx500K) and Illumina 650Y (Ilmn650K) SNP arrays. From these data, we were able to obtain estimates of genetic coverage, which are robust to overfitting, by constructing an independent test set from among these genotypes and individuals. Furthermore, we collected liver tissue RNA from the participants and profiled these samples on a comprehensive gene expression microarray. The RNA levels were used as a large-scale set of quantitative traits to calibrate the relative statistical power of the commercial arrays. Our genetic coverage estimates are lower than previous reports, providing evidence that previous estimates may be inflated due to overfitting. The Ilmn650K platform showed reasonable power (50% or greater) to detect SNPs associated with quantitative traits when the signal-to-noise ratio (SNR) is greater than or equal to 0.5 and the causal SNP's minor allele frequency (MAF) is greater than or equal to 20% (N = 359). In testing each of the more than 40,000 gene expression traits for association to each of the SNPs on the Ilmn650K and Affx500K arrays, we found that the Ilmn650K yielded 15% times more discoveries than the Affx500K at the same false discovery rate (FDR) level

Association of PCSK1 rs6234 with Obesity and Related Traits in a Chinese Han Population

Background: Common variants in PCSK1 have been reported to be associated with obesity in populations of European origin. We aimed to replicate this association in Chinese. Methodology/Principal Findings: Two PCSK1 variants rs6234 and rs6235 (in strong LD with each other, r 2 = 0.98) were genotyped in a population-based cohort of 3,210 Chinese Hans. The rs6234 was used for further association analyses with obesity and related traits. We found no significant association of rs6234 with obesity, overweight, BMI, waist circumference, or body fat percentage (P.0.05) in all participants. However, the rs6234 G-allele showed a significant association with increased risk of combined phenotype of obesity and overweight (OR 1.21[1.03–1.43], P = 0.0193) and a trend toward association with obesity (OR 1.25[0.98–1.61], P = 0.08) in men, but not in women (P$0.29). Consistently, the rs6234 G-allele showed significant association with increased BMI (P = 0.0043), waist circumference (P = 0.008) and body fat percentage (P = 0.0131) only in men, not in women (P$0.24). Interestingly, the rs6234 G-allele was significantly associated with increased HOMA-B (P = 0.0059) and decreased HOMA-S (P = 0.0349) in all participants. Conclusion/Significance: In this study, we found modest evidence for association of the PCSK1 rs6234 with BMI and overweight in men only but not in women, which suggested that PCSK1 rs6234 might not be an important contributor t

Genome-wide association study identifies loci associated with liability to alcohol and drug dependence that is associated with variability in reward-related ventral striatum activity in African- and European-Americans.

Genetic influences on alcohol and drug dependence partially overlap, however, specific loci underlying this overlap remain unclear. We conducted a genome-wide association study (GWAS) of a phenotype representing alcohol or illicit drug dependence (ANYDEP) among 7291 European-Americans (EA; 2927 cases) and 3132 African-Americans (AA: 1315 cases) participating in the family-based Collaborative Study on the Genetics of Alcoholism. ANYDEP was heritable (h 2 in EA = 0.60, AA = 0.37). The AA GWAS identified three regions with genome-wide significant (GWS; P &lt; 5E-08) single nucleotide polymorphisms (SNPs) on chromosomes 3 (rs34066662, rs58801820) and 13 (rs75168521, rs78886294), and an insertion-deletion on chromosome 5 (chr5:141988181). No polymorphisms reached GWS in the EA. One GWS region (chromosome 1: rs1890881) emerged from a trans-ancestral meta-analysis (EA + AA) of ANYDEP, and was attributable to alcohol dependence in both samples. Four genes (AA: CRKL, DZIP3, SBK3; EA: P2RX6) and four sets of genes were significantly enriched within biological pathways for hemostasis and signal transduction. GWS signals did not replicate in two independent samples but there was weak evidence for association between rs1890881 and alcohol intake in the UK Biobank. Among 118 AA and 481 EA individuals from the Duke Neurogenetics Study, rs75168521 and rs1890881 genotypes were associated with variability in reward-related ventral striatum activation. This study identified novel loci for substance dependence and provides preliminary evidence that these variants are also associated with individual differences in neural reward reactivity. Gene discovery efforts in non-European samples with distinct patterns of substance use may lead to the identification of novel ancestry-specific genetic markers of risk

Expression of the Multiple Sclerosis-Associated MHC Class II Allele HLA-DRB1*1501 Is Regulated by Vitamin D

Multiple sclerosis (MS) is a complex trait in which allelic variation in the MHC class II region exerts the single strongest effect on genetic risk. Epidemiological data in MS provide strong evidence that environmental factors act at a population level to influence the unusual geographical distribution of this disease. Growing evidence implicates sunlight or vitamin D as a key environmental factor in aetiology. We hypothesised that this environmental candidate might interact with inherited factors and sought responsive regulatory elements in the MHC class II region. Sequence analysis localised a single MHC vitamin D response element (VDRE) to the promoter region of HLA-DRB1. Sequencing of this promoter in greater than 1,000 chromosomes from HLA-DRB1 homozygotes showed absolute conservation of this putative VDRE on HLA-DRB1*15 haplotypes. In contrast, there was striking variation among non–MS-associated haplotypes. Electrophoretic mobility shift assays showed specific recruitment of vitamin D receptor to the VDRE in the HLA-DRB1*15 promoter, confirmed by chromatin immunoprecipitation experiments using lymphoblastoid cells homozygous for HLA-DRB1*15. Transient transfection using a luciferase reporter assay showed a functional role for this VDRE. B cells transiently transfected with the HLA-DRB1*15 gene promoter showed increased expression on stimulation with 1,25-dihydroxyvitamin D3 (P = 0.002) that was lost both on deletion of the VDRE or with the homologous “VDRE” sequence found in non–MS-associated HLA-DRB1 haplotypes. Flow cytometric analysis showed a specific increase in the cell surface expression of HLA-DRB1 upon addition of vitamin D only in HLA-DRB1*15 bearing lymphoblastoid cells. This study further implicates vitamin D as a strong environmental candidate in MS by demonstrating direct functional interaction with the major locus determining genetic susceptibility. These findings support a connection between the main epidemiological and genetic features of this disease with major practical implications for studies of disease mechanism and prevention

A Genome-Wide Association Study of the Metabolic Syndrome in Indian Asian Men

We conducted a two-stage genome-wide association study to identify common genetic variation altering risk of the metabolic syndrome and related phenotypes in Indian Asian men, who have a high prevalence of these conditions. In Stage 1, approximately 317,000 single nucleotide polymorphisms were genotyped in 2700 individuals, from which 1500 SNPs were selected to be genotyped in a further 2300 individuals. Selection for inclusion in Stage 1 was based on four metabolic syndrome component traits: HDL-cholesterol, plasma glucose and Type 2 diabetes, abdominal obesity measured by waist to hip ratio, and diastolic blood pressure. Association was tested with these four traits and a composite metabolic syndrome phenotype. Four SNPs reaching significance level p<5×10−7 and with posterior probability of association >0.8 were found in genes CETP and LPL, associated with HDL-cholesterol. These associations have already been reported in Indian Asians and in Europeans. Five additional loci harboured SNPs significant at p<10−6 and posterior probability >0.5 for HDL-cholesterol, type 2 diabetes or diastolic blood pressure. Our results suggest that the primary genetic determinants of metabolic syndrome are the same in Indian Asians as in other populations, despite the higher prevalence. Further, we found little evidence of a common genetic basis for metabolic syndrome traits in our sample of Indian Asian men

A Replication Study of GWAS-Derived Lipid Genes in Asian Indians: The Chromosomal Region 11q23.3 Harbors Loci Contributing to Triglycerides

Recent genome-wide association scans (GWAS) and meta-analysis studies on European populations have identified many genes previously implicated in lipid regulation. Validation of these loci on different global populations is important in determining their clinical relevance, particularly for development of novel drug targets for treating and preventing diabetic dyslipidemia and coronary artery disease (CAD). In an attempt to replicate GWAS findings on a non-European sample, we examined the role of six of these loci (CELSR2-PSRC1-SORT1 rs599839; CDKN2A-2B rs1333049; BUD13-ZNF259 rs964184; ZNF259 rs12286037; CETP rs3764261; APOE-C1-C4-C2 rs4420638) in our Asian Indian cohort from the Sikh Diabetes Study (SDS) comprising 3,781 individuals (2,902 from Punjab and 879 from the US). Two of the six SNPs examined showed convincing replication in these populations of Asian Indian origin. Our study confirmed a strong association of CETP rs3764261 with high-density lipoprotein cholesterol (HDL-C) (p = 2.03×10−26). Our results also showed significant associations of two GWAS SNPs (rs964184 and rs12286037) from BUD13-ZNF259 near the APOA5-A4-C3-A1 genes with triglyceride (TG) levels in this Asian Indian cohort (rs964184: p = 1.74×10−17; rs12286037: p = 1.58×10−2). We further explored 45 SNPs in a ∼195 kb region within the chromosomal region 11q23.3 (encompassing the BUD13-ZNF259, APOA5-A4-C3-A1, and SIK3 genes) in 8,530 Asian Indians from the London Life Sciences Population (LOLIPOP) (UK) and SDS cohorts. Five more SNPs revealed significant associations with TG in both cohorts individually as well as in a joint meta-analysis. However, the strongest signal for TG remained with BUD13-ZNF259 (rs964184: p = 1.06×10−39). Future targeted deep sequencing and functional studies should enhance our understanding of the clinical relevance of these genes in dyslipidemia and hypertriglyceridemia (HTG) and, consequently, diabetes and CAD

Genetic Markers of Obesity Risk: Stronger Associations with Body Composition in Overweight Compared to Normal-Weight Children

Genetic factors are important determinants of overweight. We examined whether there are differential effect sizes depending on children's body composition. We analysed data of n = 4,837 children recorded in the Avon Longitudinal Study of Parents and Children (ALSPAC), applying quantile regression with sex- and age-specific standard deviation scores (SDS) of body mass index (BMI) or with body fat mass index and fat-free mass index at 9 years as outcome variables and an "obesity-risk-allele score" based on eight genetic variants known to be associated with childhood BMI as the explanatory variable. The quantile regression coefficients increased with increasing child's BMI-SDS and fat mass index percentiles, indicating larger effects of the genetic factors at higher percentiles. While the associations with BMI-SDS were of similar size in medium and high BMI quantiles (40th percentile and above), effect sizes with fat mass index increased over the whole fat mass index distribution. For example, the fat mass index of a normal-weight (50th percentile) child was increased by 0.13 kg/m(2) (95% confidence interval (CI): 0.09, 0.16) per additional allele, compared to 0.24 kg/m(2) per allele (95% CI: 0.15, 0.32) in children at the 90th percentile. The genetic associations with fat-free mass index were weaker and the quantile regression effects less pronounced than those on fat mass index. Genetic risk factors for childhood overweight appear to have greater effects on fatter children. Interaction of known genetic factors with environmental or unknown genetic factors might provide a potential explanation of these findings