Genome-wide association studies of the PR interval in African Americans

The PR interval on the electrocardiogram reflects atrial and atrioventricular nodal conduction time. The PR interval is heritable, provides important information about arrhythmia risk, and has been suggested to differ among human races. Genome-wide association (GWA) studies have identified common genetic determinants of the PR interval in individuals of European and Asian ancestry, but there is a general paucity of GWA studies in individuals of African ancestry. We performed GWA studies in African American individuals from four cohorts (n = 6,247) to identify genetic variants associated with PR interval duration. Genotyping was performed using the Affymetrix 6.0 microarray. Imputation was performed for 2.8 million single nucleotide polymorphisms (SNPs) using combined YRI and CEU HapMap phase II panels. We observed a strong signal (rs3922844) within the gene encoding the cardiac sodium channel (SCN5A) with genome-wide significant association (p<2.5×10⁻⁸) in two of the four cohorts and in the meta-analysis. The signal explained 2% of PR interval variability in African Americans (beta  = 5.1 msec per minor allele, 95% CI  = 4.1–6.1, p = 3×10⁻²³). This SNP was also associated with PR interval (beta = 2.4 msec per minor allele, 95% CI = 1.8–3.0, p = 3×10⁻¹⁶) in individuals of European ancestry (n = 14,042), but with a smaller effect size (p for heterogeneity <0.001) and variability explained (0.5%). Further meta-analysis of the four cohorts identified genome-wide significant associations with SNPs in SCN10A (rs6798015), MEIS1 (rs10865355), and TBX5 (rs7312625) that were highly correlated with SNPs identified in European and Asian GWA studies. African ancestry was associated with increased PR duration (13.3 msec, p = 0.009) in one but not the other three cohorts. Our findings demonstrate the relevance of common variants to African Americans at four loci previously associated with PR interval in European and Asian samples and identify an association signal at one of these loci that is more strongly associated with PR interval in African Americans than in Europeans.     

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.     

The role of inflammatory pathway genetic variation on maternal metabolic phenotypes during pregnancy

Background

Since mediators of inflammation are associated with insulin resistance, and the risk of developing diabetes mellitus and gestational diabetes, we hypothesized that genetic variation in members of the inflammatory gene pathway impact glucose levels and related phenotypes in pregnancy. We evaluated this hypothesis by testing for association between genetic variants in 31 inflammatory pathway genes in the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) cohort, a large multiethnic multicenter study designed to address the impact of glycemia less than overt diabetes on pregnancy outcome.

Results

Fasting, 1-hour, and 2-hour glucose, fasting and 1-hour C-peptide, and HbA1c levels were measured in blood samples obtained from HAPO participants during an oral glucose tolerance test at 24-32 weeks gestation. We tested for association between 458 SNPs mapping to 31 genes in the inflammatory pathway and metabolic phenotypes in 3836 European ancestry and 1713 Thai pregnant women. The strongest evidence for association was observed with TNF alpha and HbA1c (rs1052248; 0.04% increase per allele C; p-value = 4.4×10⁻⁵), RETN and fasting plasma glucose (rs1423096; 0.7 mg/dl decrease per allele A; p-value = 1.1×10⁻⁴), IL8 and 1 hr plasma glucose (rs2886920; 2.6 mg/dl decrease per allele T; p-value = 1.3×10⁻⁴), ADIPOR2 and fasting C-peptide (rs2041139; 0.55 ug/L decrease per allele A; p-value = 1.4×10⁻⁴), LEPR and 1-hour C-peptide (rs1171278; 0.62 ug/L decrease per allele T; p-value = 2.4×10⁻⁴), and IL6 and 1-hour plasma glucose (rs6954897; −2.29 mg/dl decrease per allele G, p-value = 4.3×10⁻⁴).

Conclusions

Based on the genes surveyed in this study the inflammatory pathway is unlikely to have a strong impact on maternal metabolic phenotypes in pregnancy although variation in individual members of the pathway (e.g. RETN, IL8, ADIPOR2, LEPR, IL6, and TNF alpha,) may contribute to metabolic phenotypes in pregnant women.     

A Genome-Wide Search for Linkage of Estimated Glomerular Filtration Rate (eGFR) in the Family Investigation of Nephropathy and Diabetes (FIND)

Objective

Estimated glomerular filtration rate (eGFR), a measure of kidney function, is heritable, suggesting that genes influence renal function. Genes that influence eGFR have been identified through genome-wide association studies. However, family-based linkage approaches may identify loci that explain a larger proportion of the heritability. This study used genome-wide linkage and association scans to identify quantitative trait loci (QTL) that influence eGFR.

Methods

Genome-wide linkage and sparse association scans of eGFR were performed in families ascertained by probands with advanced diabetic nephropathy (DN) from the multi-ethnic Family Investigation of Nephropathy and Diabetes (FIND) study. This study included 954 African Americans (AA), 781 American Indians (AI), 614 European Americans (EA) and 1,611 Mexican Americans (MA). A total of 3,960 FIND participants were genotyped for 6,000 single nucleotide polymorphisms (SNPs) using the Illumina Linkage IVb panel. GFR was estimated by the Modification of Diet in Renal Disease (MDRD) formula.

Results

The non-parametric linkage analysis, accounting for the effects of diabetes duration and BMI, identified the strongest evidence for linkage of eGFR on chromosome 20q11 (log of the odds [LOD] = 3.34; P = 4.4×10⁻⁵) in MA and chromosome 15q12 (LOD = 2.84; P = 1.5×10⁻⁴) in EA. In all subjects, the strongest linkage signal for eGFR was detected on chromosome 10p12 (P = 5.5×10⁻⁴) at 44 cM near marker rs1339048. A subsequent association scan in both ancestry-specific groups and the entire population identified several SNPs significantly associated with eGFR across the genome.

Conclusion

The present study describes the localization of QTL influencing eGFR on 20q11 in MA, 15q21 in EA and 10p12 in the combined ethnic groups participating in the FIND study. Identification of causal genes/variants influencing eGFR, within these linkage and association loci, will open new avenues for functional analyses and development of novel diagnostic markers for DN.     

Differential effects of MYH9 and APOL1 risk variants on FRMD3 Association with Diabetic ESRD in African Americans

Single nucleotide polymorphisms (SNPs) in MYH9 and APOL1 on chromosome 22 (c22) are powerfully associated with non-diabetic end-stage renal disease (ESRD) in African Americans (AAs). Many AAs diagnosed with type 2 diabetic nephropathy (T2DN) have non-diabetic kidney disease, potentially masking detection of DN genes. Therefore, genome-wide association analyses were performed using the Affymetrix SNP Array 6.0 in 966 AA with T2DN and 1,032 non-diabetic, non-nephropathy (NDNN) controls, with and without adjustment for c22 nephropathy risk variants. No associations were seen between FRMD3 SNPs and T2DN before adjusting for c22 variants. However, logistic regression analysis revealed seven FRMD3 SNPs significantly interacting with MYH9—a finding replicated in 640 additional AA T2DN cases and 683 NDNN controls. Contrasting all 1,592 T2DN cases with all 1,671 NDNN controls, FRMD3 SNPs appeared to interact with the MYH9 E1 haplotype (e.g., rs942280 interaction p-value = 9.3E⁻⁷ additive; odds ratio [OR] 0.67). FRMD3 alleles were associated with increased risk of T2DN only in subjects lacking two MYH9 E1 risk haplotypes (rs942280 OR = 1.28), not in MYH9 E1 risk allele homozygotes (rs942280 OR = 0.80; homogeneity p-value = 4.3E⁻⁴). Effects were weaker stratifying on APOL1. FRMD3 SNPS were associated with T2DN, not type 2 diabetes per se, comparing AAs with T2DN to those with diabetes lacking nephropathy. T2DN-associated FRMD3 SNPs were detectable in AAs only after accounting for MYH9, with differential effects for APOL1. These analyses reveal a role for FRMD3 in AA T2DN susceptibility and accounting for c22 nephropathy risk variants can assist in detecting DN susceptibility genes.     

Genome-wide association study of blood pressure extremes identifies variant near UMOD associated with hypertension

Hypertension is a heritable and major contributor to the global burden of disease. The sum of rare and common genetic variants robustly identified so far explain only 1%–2% of the population variation in BP and hypertension. This suggests the existence of more undiscovered common variants. We conducted a genome-wide association study in 1,621 hypertensive cases and 1,699 controls and follow-up validation analyses in 19,845 cases and 16,541 controls using an extreme case-control design. We identified a locus on chromosome 16 in the 5′ region of Uromodulin (UMOD; rs13333226, combined P value of 3.6×10⁻¹¹). The minor G allele is associated with a lower risk of hypertension (OR [95%CI]: 0.87 [0.84–0.91]), reduced urinary uromodulin excretion, better renal function; and each copy of the G allele is associated with a 7.7% reduction in risk of CVD events after adjusting for age, sex, BMI, and smoking status (H.R. = 0.923, 95% CI 0.860–0.991; p = 0.027). In a subset of 13,446 individuals with estimated glomerular filtration rate (eGFR) measurements, we show that rs13333226 is independently associated with hypertension (unadjusted for eGFR: 0.89 [0.83–0.96], p = 0.004; after eGFR adjustment: 0.89 [0.83–0.96], p = 0.003). In clinical functional studies, we also consistently show the minor G allele is associated with lower urinary uromodulin excretion. The exclusive expression of uromodulin in the thick portion of the ascending limb of Henle suggests a putative role of this variant in hypertension through an effect on sodium homeostasis. The newly discovered UMOD locus for hypertension has the potential to give new insights into the role of uromodulin in BP regulation and to identify novel drugable targets for reducing cardiovascular risk.     

Common variants on chromosome 9p21 are associated with normal tension glaucoma

Although intraocular pressure (IOP) is the most definitive cause of glaucoma, a subtype of open angle glaucoma (OAG) termed normal tension glaucoma (NTG), which occurs in spite of normal IOP, accounts for a large part of glaucoma cases, especially in Japan. To find common genetic variants contributing to NTG in Japanese patients, we conducted a genome-wide association study (GWAS). We performed the first screening for 531,009 autosomal SNPs with a discovery cohort of 286 cases and 557 controls, and then a second screening for the top 30 suggestive loci in an independent cohort of 183 cases and 514 controls. Our findings identified a significantly associated SNP; rs523096 [combined p-value = 7.40× 10⁻⁸, odds ratio (OR)  = 2.00 with 95% confidence interval (CI) 1.55–2.58] located 10 kbp upstream of CDKN2B on chromosome 9p21. Moreover, analysis of another independent case-control set successfully replicated the results of the screening studies (combined values of all 3 stages p = 4.96 × 10⁻¹¹, OR  = 2.13 with 95% CI 1.69–2.68). The SNPs near rs523096 were recently reported to be associated with OAG associated with elevated IOP in primary open-angle glaucoma (POAG), the predominant subtype of glaucoma in Caucasian populations. Our results revealed that the 9p21 locus is also associated with NTG in Japanese. In addition, we identified SNPs more strongly associated with NTG.     

A genome-wide association search for type 2 diabetes genes in African Americans

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10⁻⁸). SNP rs7560163 (P = 7.0×10⁻⁹, OR (95% CI) = 0.75 (0.67–0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10⁻⁵) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.     

Association of eGFR-Related Loci Identified by GWAS with Incident CKD and ESRD

Family studies suggest a genetic component to the etiology of chronic kidney disease (CKD) and end stage renal disease (ESRD). Previously, we identified 16 loci for eGFR in genome-wide association studies, but the associations of these single nucleotide polymorphisms (SNPs) for incident CKD or ESRD are unknown. We thus investigated the association of these loci with incident CKD in 26,308 individuals of European ancestry free of CKD at baseline drawn from eight population-based cohorts followed for a median of 7.2 years (including 2,122 incident CKD cases defined as eGFR <60ml/min/1.73m² at follow-up) and with ESRD in four case-control studies in subjects of European ancestry (3,775 cases, 4,577 controls). SNPs at 11 of the 16 loci (UMOD, PRKAG2, ANXA9, DAB2, SHROOM3, DACH1, STC1, SLC34A1, ALMS1/NAT8, UBE2Q2, and GCKR) were associated with incident CKD; p-values ranged from p = 4.1e-9 in UMOD to p = 0.03 in GCKR. After adjusting for baseline eGFR, six of these loci remained significantly associated with incident CKD (UMOD, PRKAG2, ANXA9, DAB2, DACH1, and STC1). SNPs in UMOD (OR = 0.92, p = 0.04) and GCKR (OR = 0.93, p = 0.03) were nominally associated with ESRD. In summary, the majority of eGFR-related loci are either associated or show a strong trend towards association with incident CKD, but have modest associations with ESRD in individuals of European descent. Additional work is required to characterize the association of genetic determinants of CKD and ESRD at different stages of disease progression.     

Association of genetic variants in complement factor H and factor H-related genes with systemic lupus erythematosus susceptibility

Systemic lupus erythematosus (SLE), a complex polygenic autoimmune disease, is associated with increased complement activation. Variants of genes encoding complement regulator factor H (CFH) and five CFH-related proteins (CFHR1-CFHR5) within the chromosome 1q32 locus linked to SLE, have been associated with multiple human diseases and may contribute to dysregulated complement activation predisposing to SLE. We assessed 60 SNPs covering the CFH-CFHRs region for association with SLE in 15,864 case-control subjects derived from four ethnic groups. Significant allelic associations with SLE were detected in European Americans (EA) and African Americans (AA), which could be attributed to an intronic CFH SNP (rs6677604, in intron 11, P _meta = 6.6×10⁻⁸, OR = 1.18) and an intergenic SNP between CFHR1 and CFHR4 (rs16840639, P _meta = 2.9×10⁻⁷, OR = 1.17) rather than to previously identified disease-associated CFH exonic SNPs, including I62V, Y402H, A474A, and D936E. In addition, allelic association of rs6677604 with SLE was subsequently confirmed in Asians (AS). Haplotype analysis revealed that the underlying causal variant, tagged by rs6677604 and rs16840639, was localized to a ∼146 kb block extending from intron 9 of CFH to downstream of CFHR1. Within this block, the deletion of CFHR3 and CFHR1 (CFHR3-1Δ), a likely causal variant measured using multiplex ligation-dependent probe amplification, was tagged by rs6677604 in EA and AS and rs16840639 in AA, respectively. Deduced from genotypic associations of tag SNPs in EA, AA, and AS, homozygous deletion of CFHR3-1Δ (P _meta = 3.2×10⁻⁷, OR = 1.47) conferred a higher risk of SLE than heterozygous deletion (P _meta = 3.5×10⁻⁴, OR = 1.14). These results suggested that the CFHR3-1Δ deletion within the SLE-associated block, but not the previously described exonic SNPs of CFH, might contribute to the development of SLE in EA, AA, and AS, providing new insights into the role of complement regulators in the pathogenesis of SLE.     

Bivariate genome-wide association analyses of femoral neck bone geometry and appendicular lean mass

Objective

Femoral neck geometric parameters (FNGPs), such as periosteal diameter (W), cross-sectional area (CSA), cortical thickness (CT), buckling ratio (BR), and section modulus (Z), are highly genetically correlated with body lean mass. However, the specific SNPs/genes shared by these phenotypes are largely unknown.

Methods

To identify the specific SNPs/genes shared between FNGPs and appendicular lean mass (ALM), we performed an initial bivariate genome-wide association study (GWAS) by scanning ∼690,000 SNPs in 1,627 unrelated Han Chinese adults (802 males and 825 females) and a follow-up replicate study in 2,286 unrelated US Caucasians.

Results

We identified 13 interesting SNPs that may be important for both FNGPs and ALM. Two SNPs, rs681900 located in the HK2 (hexokinase 2) gene and rs11859916 in the UMOD (uromodulin) gene, were bivariately associated with FNGPs and ALM (p = 7.58×10⁻⁶ for ALM-BR and p = 2.93×10⁻⁶ for ALM-W, respectively). The associations were then replicated in Caucasians, with corresponding p values of 0.024 for rs681900 and 0.047 for rs11859916. Meta-analyses yielded combined p values of 3.05×10⁻⁶ and 2.31×10⁻⁶ for rs681900 and rs11859916, respectively. Our findings are consistent with previous biological studies that implicated HK2 and UMOD in both FNGPs and ALM. Our study also identified a group of 11 contiguous SNPs, which spanned a region of ∼130 kb, were bivariately associated with FNGPs and ALM, with p values ranging from 3.06×10⁻⁷ to 4.60×10⁻⁶ for ALM-BR. The region contained two neighboring miRNA coding genes, MIR873 (MicroRNA873) and MIR876 (MicroRNA876).

Conclusion

Our study implicated HK2, UMOD, MIR873 and MIR876, as pleiotropic genes underlying variation of both FNGPs and ALM, thus suggesting their important functional roles in co-regulating both FNGPs and ALM.     

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.     

Genome-wide association meta-analysis of cortical bone mineral density unravels allelic heterogeneity at the RANKL locus and potential pleiotropic effects on bone

Previous genome-wide association (GWA) studies have identified SNPs associated with areal bone mineral density (aBMD). However, this measure is influenced by several different skeletal parameters, such as periosteal expansion, cortical bone mineral density (BMD_C) cortical thickness, trabecular number, and trabecular thickness, which may be under distinct biological and genetic control. We have carried out a GWA and replication study of BMD_C, as measured by peripheral quantitative computed tomography (pQCT), a more homogenous and valid measure of actual volumetric bone density. After initial GWA meta-analysis of two cohorts (ALSPAC n = 999, aged ∼15 years and GOOD n = 935, aged ∼19 years), we attempted to replicate the BMD_C associations that had p<1×10⁻⁵ in an independent sample of ALSPAC children (n = 2803) and in a cohort of elderly men (MrOS Sweden, n = 1052). The rs1021188 SNP (near RANKL) was associated with BMD_C in all cohorts (overall p = 2×10⁻¹⁴, n = 5739). Each minor allele was associated with a decrease in BMD_C of ∼0.14SD. There was also evidence for an interaction between this variant and sex (p = 0.01), with a stronger effect in males than females (at age 15, males −6.77mg/cm³ per C allele, p = 2×10⁻⁶; females −2.79 mg/cm³ per C allele, p = 0.004). Furthermore, in a preliminary analysis, the rs1021188 minor C allele was associated with higher circulating levels of sRANKL (p<0.005). We show this variant to be independent from the previously aBMD associated SNP (rs9594738) and possibly from a third variant in the same RANKL region, which demonstrates important allelic heterogeneity at this locus. Associations with skeletal parameters reflecting bone dimensions were either not found or were much less pronounced. This finding implicates RANKL as a locus containing variation associated with volumetric bone density and provides further insight into the mechanism by which the RANK/RANKL/OPG pathway may be involved in skeletal development.     

Genome wide association (GWA) study for early onset extreme obesity supports the role of fat mass and obesity associated gene (FTO) variants

Background

Obesity is a major health problem. Although heritability is substantial, genetic mechanisms predisposing to obesity are not very well understood. We have performed a genome wide association study (GWA) for early onset (extreme) obesity.

Methodology/Principal Findings

a) GWA (Genome-Wide Human SNP Array 5.0 comprising 440,794 single nucleotide polymorphisms) for early onset extreme obesity based on 487 extremely obese young German individuals and 442 healthy lean German controls; b) confirmatory analyses on 644 independent families with at least one obese offspring and both parents. We aimed to identify and subsequently confirm the 15 SNPs (minor allele frequency ≥10%) with the lowest p-values of the GWA by four genetic models: additive, recessive, dominant and allelic. Six single nucleotide polymorphisms (SNPs) in FTO (fat mass and obesity associated gene) within one linkage disequilibrium (LD) block including the GWA SNP rendering the lowest p-value (rs1121980; log-additive model: nominal p = 1.13×10⁻⁷, corrected p = 0.0494; odds ratio (OR)_CT 1.67, 95% confidence interval (CI) 1.22–2.27; OR_TT 2.76, 95% CI 1.88–4.03) belonged to the 15 SNPs showing the strongest evidence for association with obesity. For confirmation we genotyped 11 of these in the 644 independent families (of the six FTO SNPs we chose only two representing the LD bock). For both FTO SNPs the initial association was confirmed (both Bonferroni corrected p<0.01). However, none of the nine non-FTO SNPs revealed significant transmission disequilibrium.

Conclusions/Significance

Our GWA for extreme early onset obesity substantiates that variation in FTO strongly contributes to early onset obesity. This is a further proof of concept for GWA to detect genes relevant for highly complex phenotypes. We concurrently show that nine additional SNPs with initially low p-values in the GWA were not confirmed in our family study, thus suggesting that of the best 15 SNPs in the GWA only the FTO SNPs represent true positive findings.     

Genetic determinants for body iron store and type 2 diabetes risk in US men and women

Background

High body iron store has been associated with an increased risk of type 2 diabetes (T2D); it remains unknown whether the genetic variants related to body iron status affect T2D risk. We aimed at comprehensively investigating the associations between the genetic variants related to body iron status and the T2D risk.

Methodology/Principal Findings

Six common SNPs related to body iron status from recent genome-wide association (GWA) studies were determined in the Nurses’ Health Study (NHS; 1,467 diabetic cases and 1,754 controls) and the Health Professionals Follow-up Study (HPFS; 1,124, diabetic cases and 1,298 controls). Plasma levels of ferritin, soluble transferrin receptor (sTfR), and transferrin were measured in NHS. Significant associations were observed for loci in TPMRSS6 with sTfR (P = 3.47×10⁻⁶), TF with transferrin (P = 0.0002 to 1.72×10⁻¹⁰); and HFE with ferritin (P = 0.017 to 1.6×10⁻⁸), sTfR (P = 0.007 to 7.9×10⁻⁶), and transferrin (P = 0.006 to 0.0007). The six SNPs together explained 5.7%, 2.7%, and 13.3% of the variation in plasma levels of ferritin, sTfR, and transferrin. After adjustment for the conventional risk factors, the T allele of SNP rs855791 in the TPMRSS6 gene was significantly associated with a 19% decreased risk of T2D (OR = 0.81; 95% CI = 0.66–0.98; P = 0.03) in men. Multiple tests attenuated this significant association to null. No associations were observed in women. SNPs at HFE and TF were not associated with diabetes risk in either sex. Dietary iron intake did not modify the associations of the newly identified loci with diabetes risk.

Conclusions/Significance

The newly identified iron-related SNP rs855791 in TPMRSS6 was nominally associated with a decreased risk of T2D in men but not in women. The apparent differences by gender warrant further study.     

A two-stage meta-analysis identifies several new loci for Parkinson's disease

A previous genome-wide association (GWA) meta-analysis of 12,386 PD cases and 21,026 controls conducted by the International Parkinson''s Disease Genomics Consortium (IPDGC) discovered or confirmed 11 Parkinson''s disease (PD) loci. This first analysis of the two-stage IPDGC study focused on the set of loci that passed genome-wide significance in the first stage GWA scan. However, the second stage genotyping array, the ImmunoChip, included a larger set of 1,920 SNPs selected on the basis of the GWA analysis. Here, we analyzed this set of 1,920 SNPs, and we identified five additional PD risk loci (combined p<5×10⁻¹⁰, PARK16/1q32, STX1B/16p11, FGF20/8p22, STBD1/4q21, and GPNMB/7p15). Two of these five loci have been suggested by previous association studies (PARK16/1q32, FGF20/8p22), and this study provides further support for these findings. Using a dataset of post-mortem brain samples assayed for gene expression (n = 399) and methylation (n = 292), we identified methylation and expression changes associated with PD risk variants in PARK16/1q32, GPNMB/7p15, and STX1B/16p11 loci, hence suggesting potential molecular mechanisms and candidate genes at these risk loci.     

A genome-wide association study of nephrolithiasis in the Japanese population identifies novel susceptible Loci at 5q35.3, 7p14.3, and 13q14.1

Nephrolithiasis is a common nephrologic disorder with complex etiology. To identify the genetic factor(s) for nephrolithiasis, we conducted a three-stage genome-wide association study (GWAS) using a total of 5,892 nephrolithiasis cases and 17,809 controls of Japanese origin. Here we found three novel loci for nephrolithiasis: RGS14-SLC34A1-PFN3-F12 on 5q35.3 (rs11746443; P = 8.51×10⁻¹², odds ratio (OR) = 1.19), INMT-FAM188B-AQP1 on 7p14.3 (rs1000597; P = 2.16×10⁻¹⁴, OR = 1.22), and DGKH on 13q14.1 (rs4142110; P = 4.62×10⁻⁹, OR = 1.14). Subsequent analyses in 21,842 Japanese subjects revealed the association of SNP rs11746443 with the reduction of estimated glomerular filtration rate (eGFR) (P = 6.54×10⁻⁸), suggesting a crucial role for this variation in renal function. Our findings elucidated the significance of genetic variations for the pathogenesis of nephrolithiasis.     

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.     

Genome-Wide Association Study of Gene by Smoking Interactions in Coronary Artery Calcification

Many GWAS have identified novel loci associated with common diseases, but have focused only on main effects of individual genetic variants rather than interactions with environmental factors (GxE). Identification of GxE interactions is particularly important for coronary heart disease (CHD), a major preventable source of morbidity and mortality with strong non-genetic risk factors. Atherosclerosis is the major cause of CHD, and coronary artery calcification (CAC) is directly correlated with quantity of coronary atherosclerotic plaque. In the current study, we tested for genetic variants influencing extent of CAC via interaction with smoking (GxS), by conducting a GxS discovery GWAS in Genetic Epidemiology Network of Arteriopathy (GENOA) sibships (N = 915 European Americans) followed by replication in Framingham Heart Study (FHS) sibships (N = 1025 European Americans). Generalized estimating equations accounted for the correlation within sibships in strata-specific groups of smokers and nonsmokers, as well as GxS interaction. Primary analysis found SNPs that showed suggestive associations (p≤10⁻⁵) in GENOA GWAS, but these index SNPs did not replicate in FHS. However, secondary analysis was able to replicate candidate gene regions in FHS using other SNPs (+/−250 kb of GENOA index SNP). In smoker and nonsmoker groups, replicated genes included TCF7L2 (p = 6.0×10⁻⁵) and WWOX (p = 4.5×10⁻⁶); and TNFRSF8 (p = 7.8×10⁻⁵), respectively. For GxS interactions, replicated genes included TBC1D4 (p = 6.9×10⁻⁵) and ADAMTS9 (P = 7.1×10⁻⁵). Interestingly, these genes are involved in inflammatory pathways mediated by the NF-κB axis. Since smoking is known to induce chronic and systemic inflammation, association of these genes likely reflects roles in CAC development via inflammatory pathways. Furthermore, the NF-κB axis regulates bone remodeling, a key physiological process in CAC development. In conclusion, GxS GWAS has yielded evidence for novel loci that are associated with CAC via interaction with smoking, providing promising new targets for future population-based and functional studies of CAC development.