Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10⁻⁶ were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10⁻⁸ for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10⁻⁸ for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10⁻⁸; RREB1: p = 5.7×10⁻⁸). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.     

A Genome-Wide Association Study Identifies Potential Susceptibility Loci for Hirschsprung Disease

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder characterized by the absence of intramural nervous plexuses along variable lengths of the hindgut. Although RET is a well-established risk factor, a recent genome-wide association study (GWAS) of HSCR has identified NRG1 as an additional susceptibility locus. To discover additional risk loci, we performed a GWAS of 123 sporadic HSCR patients and 432 unaffected controls using a large-scale platform with coverage of over 1 million polymorphic markers. The result was that our study replicated the findings of RET-CSGALNACT2-RASGEF1A genomic region (_rawP = 5.69×10⁻¹⁹ before a Bonferroni correction; _corrP = 4.31×10⁻¹³ after a Bonferroni correction) and NRG1 as susceptibility loci. In addition, this study identified SLC6A20 (_adjP = 2.71×10⁻⁶), RORA (_adjP = 1.26×10⁻⁵), and ABCC9 (_adjP = 1.86×10⁻⁵) as new potential susceptibility loci under adjusting the already known loci on the RET-CSGALNACT2-RASGEF1A and NRG1 regions, although none of the SNPs in these genes passed the Bonferroni correction. In further subgroup analysis, the RET-CSGALNACT2-RASGEF1A genomic region was observed to have different significance levels among subgroups: short-segment (S-HSCR, _corrP = 1.71×10⁻⁵), long-segment (L-HSCR, _corrP = 6.66×10⁻⁴), and total colonic aganglionosis (TCA, _corrP>0.05). This differential pattern in the significance level suggests that other genomic loci or mechanisms may affect the length of aganglionosis in HSCR subgroups during enteric nervous system (ENS) development. Although functional evaluations are needed, our findings might facilitate improved understanding of the mechanisms of HSCR pathogenesis.     

Insight in Genome-Wide Association of Metabolite Quantitative Traits by Exome Sequence Analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.     

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10⁻⁸–1.2×10⁻⁴³). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10⁻⁴). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10⁻³, n = 22,044), increased triglycerides (p = 2.6×10⁻¹⁴, n = 93,440), increased waist-to-hip ratio (p = 1.8×10⁻⁵, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10⁻³, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10⁻¹³, n = 96,748) and decreased BMI (p = 1.4×10⁻⁴, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.     

Modeling of Environmental Effects in Genome-Wide Association Studies Identifies SLC2A2 and HP as Novel Loci Influencing Serum Cholesterol Levels

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N _SE = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N _NS = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p _{SE,unadjusted} = 3.6×10⁻⁵, p _SE,adjusted = 2.2×10⁻⁶, p _{NS,unadjusted} = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p _{SE,unadjusted} = 1.1×10⁻³, p _SE,adjusted = 3.8×10⁻⁴, p _{NS,unadjusted} = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.     

Insights into the Genetic Architecture of Early Stage Age-Related Macular Degeneration: A Genome-Wide Association Study Meta-Analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.     

Obesity-Related Genomic Loci Are Associated with Type 2 Diabetes in a Han Chinese Population

Background and Aims

Obesity is a well-known risk factor for type 2 diabetes. Genome-wide association studies have identified a number of genetic loci associated with obesity. The aim of this study is to examine the contribution of obesity-related genomic loci to type 2 diabetes in a Chinese population.

Methods

We successfully genotyped 18 obesity-related single nucleotide polymorphisms among 5338 type 2 diabetic patients and 4663 controls. Both individual and joint effects of these single nucleotide polymorphisms on type 2 diabetes and quantitative glycemic traits (assessing β-cell function and insulin resistance) were analyzed using logistic and linear regression models, respectively.

Results

Two single nucleotide polymorphisms near MC4R and GNPDA2 genes were significantly associated with type 2 diabetes before adjusting for body mass index and waist circumference (OR (95% CI) = 1.14 (1.06, 1.22) for the A allele of rs12970134, P = 4.75×10⁻⁴; OR (95% CI) = 1.10 (1.03, 1.17) for the G allele of rs10938397, P = 4.54×10⁻³). When body mass index and waist circumference were further adjusted, the association of MC4R with type 2 diabetes remained significant (P = 1.81×10⁻²) and that of GNPDA2 was attenuated (P = 1.26×10⁻¹), suggesting the effect of the locus including GNPDA2 on type 2 diabetes may be mediated through obesity. Single nucleotide polymorphism rs2260000 within BAT2 was significantly associated with type 2 diabetes after adjusting for body mass index and waist circumference (P = 1.04×10⁻²). In addition, four single nucleotide polymorphisms (near or within SEC16B, BDNF, MAF and PRL genes) showed significant associations with quantitative glycemic traits in controls even after adjusting for body mass index and waist circumference (all P values<0.05).

Conclusions

This study indicates that obesity-related genomic loci were associated with type 2 diabetes and glycemic traits in the Han Chinese population.     

Genome-Wide Association Study of Breast Cancer in the Japanese Population

Breast cancer is the most common malignancy among women in worldwide including Japan. Several studies have identified common genetic variants to be associated with the risk of breast cancer. Due to the complex linkage disequilibrium structure and various environmental exposures in different populations, it is essential to identify variants associated with breast cancer in each population, which subsequently facilitate the better understanding of mammary carcinogenesis. In this study, we conducted a genome-wide association study (GWAS) as well as whole-genome imputation with 2,642 cases and 2,099 unaffected female controls. We further examined 13 suggestive loci (P<1.0×10⁻⁵) using an independent sample set of 2,885 cases and 3,395 controls and successfully validated two previously-reported loci, rs2981578 (combined P-value of 1.31×10⁻¹², OR = 1.23; 95% CI = 1.16–.30) on chromosome 10q26 (FGFR2), rs3803662 (combined P-value of 2.79×10⁻¹¹, OR = 1.21; 95% CI = 1.15–.28) and rs12922061 (combined P-value of 3.97×10⁻¹⁰, OR = 1.23; 95% CI = 1.15–.31) on chromosome 16q12 (TOX3-LOC643714). Weighted genetic risk score on the basis of three significantly associated variants and two previously reported breast cancer associated loci in East Asian population revealed that individuals who carry the most risk alleles in category 5 have 2.2 times higher risk of developing breast cancer in the Japanese population than those who carry the least risk alleles in reference category 1. Although we could not identify additional loci associated with breast cancer, our study utilized one of the largest sample sizes reported to date, and provided genetic status that represent the Japanese population. Further local and international collaborative study is essential to identify additional genetic variants that could lead to a better, accurate prediction for breast cancer.     

Genetic Determinants of Long-Term Changes in Blood Lipid Concentrations: 10-Year Follow-Up of the GLACIER Study

Recent genome-wide meta-analyses identified 157 loci associated with cross-sectional lipid traits. Here we tested whether these loci associate (singly and in trait-specific genetic risk scores [GRS]) with longitudinal changes in total cholesterol (TC) and triglyceride (TG) levels in a population-based prospective cohort from Northern Sweden (the GLACIER Study). We sought replication in a southern Swedish cohort (the MDC Study; N = 2,943). GLACIER Study participants (N = 6,064) were genotyped with the MetaboChip array. Up to 3,495 participants had 10-yr follow-up data available in the GLACIER Study. The TC- and TG-specific GRSs were strongly associated with change in lipid levels (β = 0.02 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻¹¹ for TC; β = 0.02 mmol/l per effect allele per decade follow-up, P = 5.0×10⁻⁵ for TG). In individual SNP analysis, one TC locus, apolipoprotein E (APOE) rs4420638 (β = 0.12 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻⁵), and two TG loci, tribbles pseudokinase 1 (TRIB1) rs2954029 (β = 0.09 mmol/l per effect allele per decade follow-up, P = 5.1×10⁻⁴) and apolipoprotein A-I (APOA1) rs6589564 (β = 0.31 mmol/l per effect allele per decade follow-up, P = 1.4×10⁻⁸), remained significantly associated with longitudinal changes for the respective traits after correction for multiple testing. An additional 12 loci were nominally associated with TC or TG changes. In replication analyses, the APOE rs4420638, TRIB1 rs2954029, and APOA1 rs6589564 associations were confirmed (P≤0.001). In summary, trait-specific GRSs are robustly associated with 10-yr changes in lipid levels and three individual SNPs were strongly associated with 10-yr changes in lipid levels.     

Identification, replication, and fine-mapping of Loci associated with adult height in individuals of african ancestry

Adult height is a classic polygenic trait of high heritability (h ² ∼0.8). More than 180 single nucleotide polymorphisms (SNPs), identified mostly in populations of European descent, are associated with height. These variants convey modest effects and explain ∼10% of the variance in height. Discovery efforts in other populations, while limited, have revealed loci for height not previously implicated in individuals of European ancestry. Here, we performed a meta-analysis of genome-wide association (GWA) results for adult height in 20,427 individuals of African ancestry with replication in up to 16,436 African Americans. We found two novel height loci (Xp22-rs12393627, P = 3.4×10⁻¹² and 2p14-rs4315565, P = 1.2×10⁻⁸). As a group, height associations discovered in European-ancestry samples replicate in individuals of African ancestry (P = 1.7×10⁻⁴ for overall replication). Fine-mapping of the European height loci in African-ancestry individuals showed an enrichment of SNPs that are associated with expression of nearby genes when compared to the index European height SNPs (P<0.01). Our results highlight the utility of genetic studies in non-European populations to understand the etiology of complex human diseases and traits.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.     

Targeting 160 Candidate Genes for Blood Pressure Regulation with a Genome-Wide Genotyping Array

The outcome of Genome-Wide Association Studies (GWAS) has challenged the field of blood pressure (BP) genetics as previous candidate genes have not been among the top loci in these scans. We used Affymetrix 500K genotyping data of KORA S3 cohort (n = 1,644; Southern-Germany) to address (i) SNP coverage in 160 BP candidate genes; (ii) the evidence for associations with BP traits in genome-wide and replication data, and haplotype analysis. In total, 160 gene regions (genic region±10 kb) covered 2,411 SNPs across 11.4 Mb. Marker densities in genes varied from 0 (n = 11) to 0.6 SNPs/kb. On average 52.5% of the HAPMAP SNPs per gene were captured. No evidence for association with BP was obtained for 1,449 tested SNPs. Considerable associations (P<10⁻³) were detected for the genes, where >50% of HAPMAP SNPs were tagged. In general, genes with higher marker density (>0.2 SNPs/kb) revealed a better chance to reach close to significance associations. Although, none of the detected P-values remained significant after Bonferroni correction (P<0.05/2319, P<2.15×10⁻⁵), the strength of some detected associations was close to this level: rs10889553 (LEPR) and systolic BP (SBP) (P = 4.5×10⁻⁵) as well as rs10954174 (LEP) and diastolic BP (DBP) (P = 5.20×10⁻⁵). In total, 12 markers in 7 genes (ADRA2A, LEP, LEPR, PTGER3, SLC2A1, SLC4A2, SLC8A1) revealed considerable association (P<10⁻³) either with SBP, DBP, and/or hypertension (HYP). None of these were confirmed in replication samples (KORA S4, HYPEST, BRIGHT). However, supportive evidence for the association of rs10889553 (LEPR) and rs11195419 (ADRA2A) with BP was obtained in meta-analysis across samples stratified either by body mass index, smoking or alcohol consumption. Haplotype analysis highlighted LEPR and PTGER3. In conclusion, the lack of associations in BP candidate genes may be attributed to inadequate marker coverage on the genome-wide arrays, small phenotypic effects of the loci and/or complex interaction with life-style and metabolic parameters.     

A genome-wide association study identifies novel alleles associated with hair color and skin pigmentation

We conducted a multi-stage genome-wide association study of natural hair color in more than 10,000 men and women of European ancestry from the United States and Australia. An initial analysis of 528,173 single nucleotide polymorphisms (SNPs) genotyped on 2,287 women identified IRF4 and SLC24A4 as loci highly associated with hair color, along with three other regions encompassing known pigmentation genes. We confirmed these associations in 7,028 individuals from three additional studies. Across these four studies, SLC24A4 rs12896399 and IRF4 rs12203592 showed strong associations with hair color, with p=6.0×10⁻⁶² and p=7.46×10⁻¹²⁷, respectively. The IRF4 SNP was also associated with skin color (p=6.2×10⁻¹⁴), eye color (p=6.1×10⁻¹³), and skin tanning response to sunlight (p=3.9×10⁻⁸⁹). A multivariable analysis pooling data from the initial GWAS and an additional 1,440 individuals suggested that the association between rs12203592 and hair color was independent of rs1540771, a SNP between the IRF4 and EXOC2 genes previously found to be associated with hair color. After adjustment for rs12203592, the association between rs1540771 and hair color was not significant (p=0.52). One variant in the MATP gene was associated with hair color. A variant in the HERC2 gene upstream of the OCA2 gene showed the strongest and independent association with hair color compared with other SNPs in this region, including three previously reported SNPs. The signals detected in a region around the MC1R gene were explained by MC1R red hair color alleles. Our results suggest that the IRF4 and SLC24A4 loci are associated with human hair color and skin pigmentation.     

Meta-Analysis of Genome-Wide Association Studies in African Americans Provides Insights into the Genetic Architecture of Type 2 Diabetes

Type 2 diabetes (T2D) is more prevalent in African Americans than in Europeans. However, little is known about the genetic risk in African Americans despite the recent identification of more than 70 T2D loci primarily by genome-wide association studies (GWAS) in individuals of European ancestry. In order to investigate the genetic architecture of T2D in African Americans, the MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium examined 17 GWAS on T2D comprising 8,284 cases and 15,543 controls in African Americans in stage 1 analysis. Single nucleotide polymorphisms (SNPs) association analysis was conducted in each study under the additive model after adjustment for age, sex, study site, and principal components. Meta-analysis of approximately 2.6 million genotyped and imputed SNPs in all studies was conducted using an inverse variance-weighted fixed effect model. Replications were performed to follow up 21 loci in up to 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry. We identified three known loci (TCF7L2, HMGA2 and KCNQ1) and two novel loci (HLA-B and INS-IGF2) at genome-wide significance (4.15×10⁻⁹⁴<P<5×10⁻⁸, odds ratio (OR) = 1.09 to 1.36). Fine-mapping revealed that 88 of 158 previously identified T2D or glucose homeostasis loci demonstrated nominal to highly significant association (2.2×10⁻²³ < locus-wide P<0.05). These novel and previously identified loci yielded a sibling relative risk of 1.19, explaining 17.5% of the phenotypic variance of T2D on the liability scale in African Americans. Overall, this study identified two novel susceptibility loci for T2D in African Americans. A substantial number of previously reported loci are transferable to African Americans after accounting for linkage disequilibrium, enabling fine mapping of causal variants in trans-ethnic meta-analysis studies.     

Association of Genetic Variants Influencing Lipid Levels with Coronary Artery Disease in Japanese Individuals

Background/Objective

In Japanese populations, we performed a replication study of genetic loci previously identified in European-descent populations as being associated with lipid levels and risk of coronary artery disease (CAD).

Methods

We genotyped 48 single nucleotide polymorphisms (SNPs) from 22 candidate loci that had previously been identified by genome-wide association (GWA) meta-analyses for low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and/or triglycerides in Europeans. We selected 22 loci with 2 parallel tracks from 95 reported loci: 16 significant loci (p<1×10⁻³⁰ in Europeans) and 6 other loci including those with suggestive evidence of lipid associations in 1292 GWA-scanned Japanese samples. Genotyping was done in 4990 general population samples, and 1347 CAD cases and 1337 controls. For 9 SNPs, we further examined CAD associations in an additional panel of 3052 CAD cases and 6335 controls.

Principal Findings

Significant lipid associations (one-tailed p<0.05) were replicated for 18 of 22 loci in Japanese samples, with significant inter-ethnic heterogeneity at 4 loci–APOB, APOE-C1, CETP, and APOA5–and allelic heterogeneity. The strongest association was detected at APOE rs7412 for LDL-C (p = 1.3×10⁻⁴¹), CETP rs3764261 for HDL-C (p = 5.2×10⁻²⁴), and APOA5 rs662799 for triglycerides (p = 5.8×10⁻⁵⁴). CAD association was replicated and/or verified for 4 loci: SORT1 rs611917 (p = 1.7×10⁻⁸), APOA5 rs662799 (p = 0.0014), LDLR rs1433099 (p = 2.1×10⁻⁷), and APOE rs7412 (p = 6.1×10⁻¹³).

Conclusions

Our results confirm that most of the tested lipid loci are associated with lipid traits in the Japanese, further indicating that in genetic susceptibility to lipid levels and CAD, the related metabolic pathways are largely common across the populations, while causal variants at individual loci can be population-specific.     

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci

Crohn''s disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD–susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2–4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.     

HLA-DRB1 Shared Epitope-Dependent DR-DQ Haplotypes Are Associated with Both Anti-CCP–Positive and –Negative Rheumatoid Arthritis in Chinese Han

The association between Human Leukocyte Antigen (HLA) class II and rheumatoid arthritis (RA) has been extensively studied, but few reported DR-DQ haplotype. Here we investigated the association of HLA-DRB1, DQA1, DQB1, and DR-DQ haplotypes with RA susceptibility and with anti-CCP antibodies in 281 RA patients and 297 control in Han population. High-resolution genotyping were performed. The HLA-DRB1 shared epitope (SE)-encoding allele *0405 displayed the most significant RA association (P = 1.35×10⁻⁶). The grouped DRB1 SE alleles showed great association with RA (P = 3.88×10⁻¹³). The DRB1 DRRAA alleles displayed significant protective effects (P = 0.021). The SE-dependent DR-DQ haplotype SE-DQ3/4/5 remained strong association with both anti-CCP -positive (P = 3.71×10⁻¹³) and -negative RA (P = 3.89×10⁻⁵). Our study revealed that SE alleles and its haplotypes SE-DQ3/4/5 were highly associated with RA susceptibility in Han population. The SE-DQ3/4/5 haplotypes were associated with both anti-CCP positive RA and -negative RA.     

New Variants Including ARG1 Polymorphisms Associated with C-Reactive Protein Levels Identified by Genome-Wide Association and Pathway Analysis

C-reactive protein (CRP) is a general marker of systemic inflammation and cardiovascular disease (CVD). The genetic contribution to differences in CRP levels remains to be explained, especially in non-European populations. Thus, the aim of this study was to identify genetic loci associated with CRP levels in Korean population. We performed genome-wide association studies (GWAS) using SNPs from 8,529 Korean individuals (7,626 for stage 1 and 903 for stage 2). We also performed pathway analysis. We identified a new genetic locus associated with CRP levels upstream of ARG1 gene (top significant SNP: rs9375813, P_meta = 2.85×10⁻⁸), which encodes a key enzyme of the urea cycle counteract the effects of nitric oxide, in addition to known CRP (rs7553007, P_meta = 1.72×10⁻¹⁶) and HNF1A loci (rs2259816, P_meta = 2.90×10⁻¹⁰). When we evaluated the associations between the CRP-related SNPs with cardiovascular disease phenotypes, rs9375813 (ARG1) showed a marginal association with hypertension (P = 0.0440). To identify more variants and pathways, we performed pathway analysis and identified six candidate pathways comprised of genes related to inflammatory processes and CVDs (CRP, HNF1A, PCSK6, CD36, and ABCA1). In addition to the previously reported loci (CRP, HNF1A, and IL6) in diverse ethnic groups, we identified novel variants in the ARG1 locus associated with CRP levels in Korean population and a number of interesting genes related to inflammatory processes and CVD through pathway analysis.     

European bone mineral density loci are also associated with BMD in East-Asian populations

Most genome-wide association (GWA) studies have focused on populations of European ancestry with limited assessment of the influence of the sequence variants on populations of other ethnicities. To determine whether markers that we have recently shown to associate with Bone Mineral Density (BMD) in Europeans also associate with BMD in East-Asians we analysed 50 markers from 23 genomic loci in samples from Korea (n = 1,397) and two Chinese Hong Kong sample sets (n = 3,869 and n = 785). Through this effort we identified fourteen loci that associated with BMD in East-Asian samples using a false discovery rate (FDR) of 0.05; 1p36 (ZBTB40, P = 4.3×10⁻⁹), 1p31 (GPR177, P = 0.00012), 3p22 (CTNNB1, P = 0.00013), 4q22 (MEPE, P = 0.0026), 5q14 (MEF2C, P = 1.3×10⁻⁵), 6q25 (ESR1, P = 0.0011), 7p14 (STARD3NL, P = 0.00025), 7q21 (FLJ42280, P = 0.00017), 8q24 (TNFRSF11B, P = 3.4×10⁻⁵), 11p15 (SOX6, P = 0.00033), 11q13 (LRP5, P = 0.0033), 13q14 (TNFSF11, P = 7.5×10⁻⁵), 16q24 (FOXL1, P = 0.0010) and 17q21 (SOST, P = 0.015). Our study marks an early effort towards the challenge of cataloguing bone density variants shared by many ethnicities by testing BMD variants that have been established in Europeans, in East-Asians.     

Novel Association of HK1 with Glycated Hemoglobin in a Non-Diabetic Population: A Genome-Wide Evaluation of 14,618 Participants in the Women's Genome Health Study

Type 2 diabetes is a leading cause of morbidity and mortality. While genetic variants have been found to influence the risk of type 2 diabetes mellitus, relatively few studies have focused on genes associated with glycated hemoglobin, an index of the mean blood glucose concentration of the preceding 8–12 weeks. Epidemiologic studies and randomized clinical trials have documented the relationship between glycated hemoglobin levels and the development of long-term complications in diabetes; moreover, higher glycated hemoglobin levels in the subdiabetic range have been shown to predict type 2 diabetes risk and cardiovascular disease. To examine the common genetic determinants of glycated hemoglobin levels, we performed a genome-wide association study that evaluated 337,343 SNPs in 14,618 apparently healthy Caucasian women. The results show that glycated hemoglobin levels are associated with genetic variation at the GCK (rs730497; P = 2.8×10⁻¹²), SLC30A8 (rs13266634; P = 9.8×10⁻⁸), G6PC2 (rs1402837; P = 6.8×10⁻¹⁰), and HK1 (rs7072268; P = 6.4×10⁻⁹) loci. While associations at the GCK, SLC30A8, and G6PC2 loci are confirmatory, the findings at HK1 are novel. We were able to replicate this novel association in an independent validation sample of 455 additional non-diabetic men and women. HK1 encodes the enzyme hexokinase, the first step in glycolysis and a likely candidate for the control of glucose metabolism. This observed genetic association between glycated hemoglobin levels and HK1 polymorphisms paves the way for further studies of the role of HK1 in hemoglobin glycation, glucose metabolism, and diabetes.