Replication of genome-wide association studies (GWAS) loci for fasting plasma glucose in African-Americans

Melatonin secretion and the incidence of type 2 diabetes

IMPORTANCE

Loss-of-function mutations in the melatonin receptor are associated with insulin resistance and type 2 diabetes. Additionally, in a cross-sectional analysis of persons without diabetes, lower nocturnal melatonin secretion was associated with increased insulin resistance.

OBJECTIVE

To study the association between melatonin secretion and the risk of developing type 2 diabetes.

DESIGN, SETTING, AND PARTICIPANTS

Case-control study nested within the Nurses' Health Study cohort. Among participants without diabetes who provided urine and blood samples at baseline in 2000, we identified 370 women who developed type 2 diabetes from 2000-2012 and matched 370 controls using risk-set sampling.

MAIN OUTCOME MEASURES

Associations between melatonin secretion at baseline and incidence of type 2 diabetes were evaluated with multivariable conditional logistic regression controlling for demographic characteristics, lifestyle habits, measures of sleep quality, and biomarkers of inflammation and endothelial dysfunction.

RESULTS

The median urinary ratios of 6-sulfatoxymelatonin to creatinine were 28.2 ng/mg (5%-95% range, 5.5-84.2 ng/mg) among cases and 36.3 ng/mg (5%-95% range, 6.9-110.8 ng/mg) among controls. Women with lower ratios of 6-sulfatoxymelatonin to creatinine had increased risk of diabetes (multivariable odds ratio, 1.48 [95% CI, 1.11-1.98] per unit decrease in the estimated log ratio of 6-sulfatoxymelatonin to creatinine). Compared with women in the highest ratio category of 6-sulfatoxymelatonin to creatinine, those in the lowest category had a multivariable odds ratio of 2.17 (95% CI, 1.18-3.98) of developing type 2 diabetes. Women in the highest category of melatonin secretion had an estimated diabetes incidence rate of 4.27 cases/1000 person-years compared with 9.27 cases/1000 person-years in the lowest category.

CONCLUSIONS AND RELEVANCE

Lower melatonin secretion was independently associated with a higher risk of developing type 2 diabetes. Further research is warranted to assess if melatonin secretion is a modifiable risk factor for diabetes within the general population.

Combined genome-wide linkage and association analyses of fasting glucose level in healthy twins and families of Korea

This study was undertaken to identify genetic polymorphisms that are associated with the risk of an elevated fasting glucose (FG) level using genome-wide analyses. We explored a quantitative trait locus (QTL) for FG level in a genome-wide study from a Korean twin-family cohort (the Healthy Twin Study) using a combined linkage and family-based association analysis approach. We investigated 1,754 individuals, which included 432 families and 219 pairs of monozygotic twins. Regions of chromosomes 2q23.3-2q31.1, 15q26.1-15q26.3, 16p12.1, and 20p13-20p12.2, were found to show evidence of linkage with FG level, and several markers in these regions were found to be significantly associated with FG level using family-based or general association tests. In particular, a single-nucleotide polymorphism (rs6138953) on the PTPRA gene in the 20p13 region (combined P = 1.8 × 10(-6)) was found to be associated with FG level, and the PRKCB1 gene (in 16p12.1) to be possibly associated with FG level. In conclusion, multiple regions of chromosomes 2q23.3-2q31.1, 15q26.1-15q26.3, 16p12.1, and 20p13-20p12.2 are associated with FG level in our Korean twin-family cohort. The combined approach of genome-wide linkage and family-based association analysis is useful to identify novel or known genetic regions concerning FG level in a family cohort study.

Animal models of GWAS-identified type 2 diabetes genes

More than 65 loci, encoding up to 500 different genes, have been implicated by genome-wide association studies (GWAS) as conferring an increased risk of developing type 2 diabetes (T2D). Whilst mouse models have in the past been central to understanding the mechanisms through which more penetrant risk genes for T2D, for example, those responsible for neonatal or maturity-onset diabetes of the young, only a few of those identified by GWAS, notably TCF7L2 and ZnT8/SLC30A8, have to date been examined in mouse models. We discuss here the animal models available for the latter genes and provide perspectives for future, higher throughput approaches towards efficiently mining the information provided by human genetics.

Why personalized medicine will fail if we stay the course

Genomic science and associated technologies are providing scientists and clinicians with novel insights that are transforming the delivery of healthcare and the overall well-being of society. However, these insights inform us that historical population sampling approaches for investigating rare and common genetic variations are not representative of the complex ancestral backgrounds of today's patients. In order for personalized medicine to be meaningful and applicable to the global populations, we will need to know how common and rare genetic variants found in different parts of the world influence health and drug response. This article demonstrates the importance of increasing ethnic and racial diversity among participants in genomic research, highlights areas of opportunity for improving our understanding of genomic diversity among populations, and provides examples of successful models that help to resolve these concerns.

Transferability and fine-mapping of glucose and insulin quantitative trait loci across populations: CARe, the Candidate Gene Association Resource

AIMS/HYPOTHESIS

Hyperglycaemia disproportionately affects African-Americans (AfAs). We tested the transferability of 18 single-nucleotide polymorphisms (SNPs) associated with glycaemic traits identified in European ancestry (EuA) populations in 5,984 non-diabetic AfAs.

METHODS

We meta-analysed SNP associations with fasting glucose (FG) or insulin (FI) in AfAs from five cohorts in the Candidate Gene Association Resource. We: (1) calculated allele frequency differences, variations in linkage disequilibrium (LD), fixation indices (F(st)s) and integrated haplotype scores (iHSs); (2) tested EuA SNPs in AfAs; and (3) interrogated within ± 250 kb around each EuA SNP in AfAs.

RESULTS

Allele frequency differences ranged from 0.6% to 54%. F(st) exceeded 0.15 at 6/16 loci, indicating modest population differentiation. All iHSs were <2, suggesting no recent positive selection. For 18 SNPs, all directions of effect were the same and 95% CIs of association overlapped when comparing EuA with AfA. For 17 of 18 loci, at least one SNP was nominally associated with FG in AfAs. Four loci were significantly associated with FG (GCK, p = 5.8 × 10(-8); MTNR1B, p = 8.5 × 10(-9); and FADS1, p = 2.2 × 10(-4)) or FI (GCKR, p = 5.9 × 10(-4)). At GCK and MTNR1B the EuA and AfA SNPs represented the same signal, while at FADS1, and GCKR, the EuA and best AfA SNPs were weakly correlated (r(2) <0.2), suggesting allelic heterogeneity for association with FG at these loci.

CONCLUSIONS/INTERPRETATION

Few glycaemic SNPs showed strict evidence of transferability from EuA to AfAs. Four loci were significantly associated in both AfAs and those with EuA after accounting for varying LD across ancestral groups, with new signals emerging to aid fine-mapping.

Fine-mapping of the 6q25 locus identifies a novel SNP associated with breast cancer risk in African-American women

The rs2046210 single nucleotide polymorphism (SNP) in the 6q25.1 region was identified in a breast cancer genome-wide association study of Chinese women. The SNP has been replicated in European ancestry populations, but replication efforts have failed in African ancestry populations. We evaluated a total of 13 tagging SNPs in the linkage disequilibrium block around rs2046210 in a case-control study of breast cancer nested within the Black Women's Health Study, which included 1191 cases and 1941 controls. Replication of initial significant findings was carried out in 665 cases and 821 controls of African ancestry from the Women's Circle of Health Study (WCHS). No significant association was found for rs2046210 in univariate analysis. A new SNP, rs2046211, was significantly associated with reduced risk of breast cancer and was replicated in data from WCHS. In joint analyses that included both SNPs, the rs2046210-A allele was associated with increased risk of breast cancer [odds ratio (OR) = 1.14; 95% confidence interval (CI) = 1.02-1.28], and the rs2046211-G allele was associated with reduced risk (OR = 0.80; 95% CI = 0.67-0.95). Haplotype analysis confirmed these results and showed that the rs2046210-A allele is present in high-risk (rs2046211-C/rs2046210-A) and low-risk (rs2046211-G/rs2046210-A) haplotypes. Our results confirm the importance of 6q25.1 as a breast cancer susceptibility region. We replicated the rs2046210 association, after accounting for the haplotype background that included rs2046211 in African-American women, and we report the presence of a novel signal that is tagged by rs2046211.

Transferability and fine mapping of genome-wide associated loci for lipids in African Americans

Background

A recent, large genome-wide association study (GWAS) of European ancestry individuals has identified multiple genetic variants influencing serum lipids. Studies of the transferability of these associations to African Americans remain few, an important limitation given interethnic differences in serum lipids and the disproportionate burden of lipid-associated metabolic diseases among African Americans.

Methods

We attempted to evaluate the transferability of 95 lipid-associated loci recently identified in European ancestry individuals to 887 non-diabetic, unrelated African Americans from a population-based sample in the Washington, DC area. Additionally, we took advantage of the generally reduced linkage disequilibrium among African ancestry populations in comparison to European ancestry populations to fine-map replicated GWAS signals.

Results

We successfully replicated reported associations for 10 loci (CILP2/SF4, STARD3, LPL, CYP7A1, DOCK7/ANGPTL3, APOE, SORT1, IRS1, CETP, and UBASH3B). Through trans-ethnic fine-mapping, we were able to reduce associated regions around 75% of the loci that replicated.

Conclusions

Between this study and previous work in African Americans, 40 of the 95 loci reported in a large GWAS of European ancestry individuals also influence lipid levels in African Americans. While there is now evidence that the lipid-influencing role of a number of genetic variants is observed in both European and African ancestry populations, the still considerable lack of concordance highlights the importance of continued ancestry-specific studies to elucidate the genetic underpinnings of these traits.

Multiple loci associated with renal function in African Americans

The incidence of chronic kidney disease varies by ethnic group in the USA, with African Americans displaying a two-fold higher rate than European Americans. One of the two defining variables underlying staging of chronic kidney disease is the glomerular filtration rate. Meta-analysis in individuals of European ancestry has identified 23 genetic loci associated with the estimated glomerular filtration rate (eGFR). We conducted a follow-up study of these 23 genetic loci using a population-based sample of 1,018 unrelated admixed African Americans. We included in our follow-up study two variants in APOL1 associated with end-stage kidney disease discovered by admixture mapping in admixed African Americans. To address confounding due to admixture, we estimated local ancestry at each marker and global ancestry. We performed regression analysis stratified by local ancestry and combined the resulting regression estimates across ancestry strata using an inverse variance-weighted fixed effects model. We found that 11 of the 24 loci were significantly associated with eGFR in our sample. The effect size estimates were not significantly different between the subgroups of individuals with two copies of African ancestry vs. two copies of European ancestry for any of the 11 loci. In contrast, allele frequencies were significantly different at 10 of the 11 loci. Collectively, the 11 loci, including four secondary signals revealed by conditional analyses, explained 14.2% of the phenotypic variance in eGFR, in contrast to the 1.4% explained by the 24 loci in individuals of European ancestry. Our findings provide insight into the genetic basis of variation in renal function among admixed African Americans.

Genome-wide association study identifies novel loci association with fasting insulin and insulin resistance in African Americans

Insulin resistance (IR) is a key determinant of type 2 diabetes (T2D) and other metabolic disorders. This genome-wide association study (GWAS) was designed to shed light on the genetic basis of fasting insulin (FI) and IR in 927 non-diabetic African Americans. 5 396 838 single-nucleotide polymorphisms (SNPs) were tested for associations with FI or IR with adjustments for age, sex, body mass index, hypertension status and first two principal components. Genotyped SNPs (n = 12) with P < 5 × 10(-6) in African Americans were carried forward for de novo genotyping in 570 non-diabetic West Africans. We replicated SNPs in or near SC4MOL and TCERG1L in West Africans. The meta-analysis of 1497 African Americans and West Africans yielded genome-wide significant associations for SNPs in the SC4MOL gene: rs17046216 (P = 1.7 × 10(-8) and 2.9 × 10(-8) for FI and IR, respectively); and near the TCERG1L gene with rs7077836 as the top scoring (P = 7.5 × 10(-9) and 4.9 × 10(-10) for FI and IR, respectively). In silico replication in the MAGIC study (n = 37 037) showed weak but significant association (adjusted P-value of 0.0097) for rs34602777 in the MYO5A gene. In addition, we replicated previous GWAS findings for IR and FI in Europeans for GCKR, and for variants in four T2D loci (FTO, IRS1, KLF14 and PPARG) which exert their action via IR. In summary, variants in/near SC4MOL, and TCERG1L were associated with FI and IR in this cohort of African Americans and were replicated in West Africans. SC4MOL is under-expressed in an animal model of T2D and plays a key role in lipid biosynthesis, with implications for the regulation of energy metabolism, obesity and dyslipidemia. TCERG1L is associated with plasma adiponectin, a key modulator of obesity, inflammation, IR and diabetes.

Fasting glucose GWAS candidate region analysis across ethnic groups in the Multiethnic Study of Atherosclerosis (MESA)

Genetic variants associated with fasting glucose in European ancestry populations are increasingly well understood. However, the nature of the associations between these single nucleotide polymorphisms (SNPs) and fasting glucose in other racial and ethnic groups is unclear. We sought to examine regions previously identified to be associated with fasting glucose in Caucasian genome-wide association studies (GWAS) across multiple ethnicities in the Multiethnic Study of Atherosclerosis (MESA). Nondiabetic MESA participants with fasting glucose measured at the baseline exam and with GWAS genotyping were included; 2,349 Caucasians, 664 individuals of Chinese descent, 1,366 African Americans, and 1,171 Hispanics. Genotype data were generated from the Affymetrix 6.0 array and imputation in IMPUTE. Fasting glucose was regressed on SNP dosage data in each ethnic group adjusting for age, gender, MESA study center, and ethnic-specific principal components. SNPs from the three gene regions with the strongest associations to fasting glucose in previous Caucasian GWAS (MTNR1B / GCK / G6PC2) were examined in depth. There was limited power to replicate associations in other ethnic groups due to smaller allele frequencies and limited sample size; SNP associations may also have differed across ethnic groups due to differing linkage disequilibrium patterns with causal variants. rs10830963 in MTNR1B and rs4607517 in GCK demonstrated consistent magnitude and direction of association with fasting glucose across ethnic groups, although the associations were often not nominally significant. In conclusion, certain SNPs in MTNR1B and GCK demonstrate consistent effects across four racial and ethnic groups, narrowing the putative region for these causal variants.

C-reactive protein (CRP) promoter polymorphisms influence circulating CRP levels in a genome-wide association study of African Americans

C-reactive protein (CRP) is an acute phase reactant protein produced primarily by the liver. Circulating CRP levels are influenced by genetic and non-genetic factors, including infection and obesity. Genome-wide association studies (GWAS) provide an unbiased approach towards identifying loci influencing CRP levels. None of the six GWAS for CRP levels has been conducted in an African ancestry population. The present study aims to: (i) identify genetic variants that influence serum CRP in African Americans (AA) using a genome-wide association approach and replicate these findings in West Africans (WA), (ii) assess transferability of major signals for CRP reported in European ancestry populations (EA) to AA and (iii) use the weak linkage disequilibrium (LD) structure characteristic of African ancestry populations to fine-map the previously reported CRP locus. The discovery cohort comprised 837 unrelated AA, with the replication of significant single-nucleotide polymorphisms (SNPs) assessed in 486 WA. The association analysis was conducted with 2 366 856 genotyped and imputed SNPs under an additive genetic model with adjustment for appropriate covariates. Genome-wide and replication significances were set at P < 5 × 10(-8) and P < 0.05, respectively. Ten SNPs in (CRP pseudogene-1) CRPP1 and CRP genes were associated with serum CRP (P = 2.4 × 10(-09) to 4.3 × 10(-11)). All but one of the top-scoring SNPs associated with CRP in AA were successfully replicated in WA. CRP signals previously identified in EA samples were transferable to AAs, and we were able to fine-map this signal, reducing the region of interest from the 25 kb of LD around the locus in the HapMap CEU sample to only 8 kb in our AA sample.

Genetic determinants of glucose homeostasis

Type 2 diabetes is a complex metabolic disorder characterised by varying degrees of impairment in insulin secretion and resistance to the action of insulin. Considerable progress has been made recently in understanding the genetic determinants of diabetes. A logical next step is to describe how these variants relate to the underlying pathophysiological processes that lead to diabetes as this may provide insights into pathways to disease. These quantitative traits are, of course, of direct interest in themselves and a growing literature is now emerging on the genetic determinants of insulin secretion and insulin resistance. This review article focuses on describing the complex associations between type 2 diabetes risk variants and quantitative glycaemic traits and the relationship between variants initially discovered in association studies of these traits and risk of type 2 diabetes.

Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci

To identify genetic factors contributing to type 2 diabetes (T2D), we performed large-scale meta-analyses by using a custom ∼50,000 SNP genotyping array (the ITMAT-Broad-CARe array) with ∼2000 candidate genes in 39 multiethnic population-based studies, case-control studies, and clinical trials totaling 17,418 cases and 70,298 controls. First, meta-analysis of 25 studies comprising 14,073 cases and 57,489 controls of European descent confirmed eight established T2D loci at genome-wide significance. In silico follow-up analysis of putative association signals found in independent genome-wide association studies (including 8,130 cases and 38,987 controls) performed by the DIAGRAM consortium identified a T2D locus at genome-wide significance (GATAD2A/CILP2/PBX4; p = 5.7 × 10(-9)) and two loci exceeding study-wide significance (SREBF1, and TH/INS; p < 2.4 × 10(-6)). Second, meta-analyses of 1,986 cases and 7,695 controls from eight African-American studies identified study-wide-significant (p = 2.4 × 10(-7)) variants in HMGA2 and replicated variants in TCF7L2 (p = 5.1 × 10(-15)). Third, conditional analysis revealed multiple known and novel independent signals within five T2D-associated genes in samples of European ancestry and within HMGA2 in African-American samples. Fourth, a multiethnic meta-analysis of all 39 studies identified T2D-associated variants in BCL2 (p = 2.1 × 10(-8)). Finally, a composite genetic score of SNPs from new and established T2D signals was significantly associated with increased risk of diabetes in African-American, Hispanic, and Asian populations. In summary, large-scale meta-analysis involving a dense gene-centric approach has uncovered additional loci and variants that contribute to T2D risk and suggests substantial overlap of T2D association signals across multiple ethnic groups.

Characterization of Bglu3, a mouse fasting glucose locus, and identification of Apcs as an underlying candidate gene

Bglu3 is a quantitative trait locus for fasting glucose on distal chromosome 1 identified in an intercross between C57BL/6 (B6) and C3H/HeJ (C3H) apolipoprotein E-deficient (apoE(-/-)) mice. This locus was subsequently replicated in two separate mouse intercrosses. The objective of this study was to characterize Bglu3 through construction and analysis of a congenic strain and identify underlying candidate genes. Congenic mice were constructed by introgressing a genomic region harboring Bglu3 from C3H.apoE(-/-) into B6.apoE(-/-) mice. Mice were started with a Western diet at 6 wk of age and maintained on the diet for 12 wk. Gene expression in the liver was analyzed by microarrays. Congenic mice had significantly higher fasting glucose levels and developed more significant glucose intolerance compared with B6.apoE(-/-) mice on the Western diet. Microarray analysis revealed 336 genes to be differentially expressed in the liver of congenic mice. Further pathway analysis suggested a role for acute phase response signaling in regulating glucose intolerance. Apcs, encoding an acute phase response protein serum amyloid P (SAP), is located underneath the linkage peak of Bglu3. Multiple single nucleotide polymorphisms between B6 and C3H mice were detected within and surrounding Apcs. Apcs expression in the liver was significantly higher in congenic and C3H mice compared with B6 mice. The Western diet consumption led to a gradual rise in plasma SAP levels, which was accompanied by rising fasting glucose in both B6 and C3H apoE(-/-) mice. Expression of C3H Apcs in B6.apoE(-/-) mice aggravated glucose intolerance. Bglu3 is confirmed to be a locus affecting diabetes susceptibility, and Apcs is a probable candidate gene.

UGT1A1 is a major locus influencing bilirubin levels in African Americans

Total serum bilirubin is associated with several clinical outcomes, including cardiovascular disease, diabetes and drug metabolism. We conducted a genome-wide association study in 619 healthy unrelated African Americans in an attempt to replicate reported findings in Europeans and Asians and to identify novel loci influencing total serum bilirubin levels. We analyzed a dense panel of over two million genotyped and imputed SNPs in additive genetic models adjusting for age, sex, and the first two significant principal components from the sample covariance matrix of genotypes. Thirty-nine SNPs spanning a 78 kb region within the UGT1A1 displayed P-values <5 × 10⁻⁸. The lowest P-value was 1.7 × 10⁻²² for SNP rs887829. None of SNPs in the UGT1A1 remained statistically significant in conditional association analyses that adjusted for rs887829. In addition, SNP rs10929302 located in phenobarbital response enhancer module was significantly associated with bilirubin level with a P-value of 1.37 × 10⁻¹¹; this enhancer module is believed to have a critical role in phenobarbital treatment of hyperbilirubinemia. Interestingly, the lead SNP, rs887829, is in strong linkage disequilibrium (LD) (r²≥0.74) with rs10929302. Taking advantage of the lower LD and shorter haplotypes in African-ancestry populations, we identified rs887829 as a more refined proxy for the causative variant influencing bilirubin levels. Also, we replicated the reported association between variants in SEMA3C and bilirubin levels. In summary, UGT1A1 is a major locus influencing bilirubin levels and the results of this study promise to contribute to understanding of the etiology and treatment of hyperbilirubinaemia in African-ancestry populations.

The role of osteocalcin in the endocrine cross-talk between bone remodelling and energy metabolism

Bone remodelling, which maintains bone mass constant during adulthood, is an energy-demanding process. This, together with the observation that the adipocyte-derived hormone leptin is a major inhibitor of bone remodelling, led to the hypothesis that bone cells regulate energy metabolism through an endocrine mechanism. Studies to test this hypothesis identified osteocalcin, a hormone secreted by osteoblasts, as a positive regulator of insulin secretion, insulin resistance and energy expenditure. Remarkably, insulin signalling in osteoblasts is a positive regulator of osteocalcin production and activation via its ability to indirectly enhance bone resorption by osteoclasts. In contrast, leptin is a potent inhibitor of osteocalcin function through its effect on the sympathetic tone. Hence, osteocalcin is part of a complex signalling network between bone and the organs more classically associated with the regulation of energy homeostasis, such as the pancreas and adipose tissue. This review summarises the molecular and cellular bases of the present knowledge on osteocalcin biology and discusses the potential relevance of osteocalcin to human metabolism and pathology.